Jean-Noël Fuchs
email: fuchs@lptmc.jussieu.fr

Rémy Mosseri
email: remy.mosseri@upmc.fr

Julien Vidal
email: vidal@lptmc.jussieu.fr

Landau level broadening, hyperuniformity, and discrete scale invariance

We study the energy spectrum of a two-dimensional electron in the presence of both a perpendicular magnetic field and a potential. In the limit where the potential is small compared to the Landau level spacing, we show that the broadening of Landau levels is simply expressed in terms of the structure factor of the potential. For potentials that are either periodic or random, we recover known results. Interestingly, for potentials with a dense Fourier spectrum made of Bragg peaks (as found, e.g., in quasicrystals), we find an algebraic broadening with the magnetic field characterized by the hyperuniformity exponent of the potential. Furthermore, if the potential is self-similar such that its structure factor has a discrete scale invariance, the broadening displays log-periodic oscillations together with an algebraic envelope.

I. INTRODUCTION

In the presence of a magnetic field, the energy spectrum of noninteracting electrons in two dimensions is known to consist of Landau levels. These discrete energy levels are responsible for many remarkable phenomena, among which is the celebrated integer quantum Hall effect [START_REF] Klitzing | New Method for High-Accuracy Determination of the Fine-Structure Constant Based on Quantized Hall Resistance[END_REF][START_REF] Klitzing | The quantized Hall effect[END_REF]. Each Landau level has a macroscopic degeneracy that is proportional to the strength of the magnetic field. This degeneracy is expected to be lifted by a generic perturbation leading to a broadening of Landau levels that may have important physical consequences. For instance, plateaus observed in the Hall resistance are directly related to the broadening induced by disorder, as realized early by Ando et al. [START_REF] Ando | Theory of quantum transport in a two-dimensional electron system under magnetic fields. I. Characteristics of level broadening and transport under strong fields[END_REF][START_REF] Ando | Theory of quantum transport in a two-dimensional electron system under magnetic fields. II. Single-site approximation under strong fields[END_REF][START_REF] Ando | Theory of quantum transport in a two-dimensional electron system under magnetic fields. III. Many-site approximation[END_REF]. Most studies on Landau level broadening focused on disordered systems (see Ref. [START_REF] Huckestein | Scaling theory of the integer quantum Hall effect[END_REF] for a review), but the role played by periodic potentials has also attracted much attention following the original work of Rauh [START_REF] Rauh | Degeneracy of landau levels in crystals[END_REF][START_REF] Rauh | On the broadening of landau levels in crystals[END_REF]. Based on a free-electron picture, Rauh's approach also allows one to qualitatively understand the Landau level broadening in the small-field limit of the Hofstadter butterfly for periodic lattices [START_REF] Hofstadter | Energy levels and wave functions of Bloch electrons in rational and irrational magnetic fields[END_REF][START_REF] Claro | Magnetic subband structure of electrons in hexagonal lattices[END_REF], although a quantitative analysis requires a semiclassical treatment [START_REF] Wilkinson | Critical properties of electron eigenstates in incommensurate systems[END_REF]. Recently, the Hofstadter butterfly of some quasiperiodic systems has been investigated, unveiling an unusual broadening of Landau levels [START_REF] Fuchs | Landau levels in quasicrystals[END_REF] different from the one expected for periodic or disordered systems, hence suggesting a nontrivial mechanism for potentials with a dense set of Bragg peaks.

The goal of the present paper is to provide a general framework to compute the broadening of Landau levels in the presence of an arbitrary potential. Our main result, given in Eq. [START_REF] Wilkinson | Critical properties of electron eigenstates in incommensurate systems[END_REF], relates the variance of the lowest Landau level (LLL) to the structure factor of the perturbing potential (an extension to higher-energy Landau levels is straightforward). This simple expression reproduces the aforementioned results for disordered and periodic cases, but it also allows us to investigate more subtle potentials (see Fig. 1 for a summary of the results). In particular, we find that when the Fourier spectrum of the potential is dense and made of Bragg peaks (as in quasicrystals), the variance of the LLL increases algebraically with the magnetic field [see Eq. [START_REF] Elser | Indexing problems in quasicrystal diffraction[END_REF]] with an exponent characterizing the hyperuniformity of the potential. This notion of hyperuniformity is commonly used to describe sets of points with an unusually large suppression of density fluctuations at long wavelengths [START_REF] Torquato | Hyperuniform states of matter[END_REF]. We also show that if the potential has a discrete scale invariance [START_REF] Sornette | Discrete-scale invariance and complex dimensions[END_REF][START_REF] Akkermans | Statistical mechanics and quantum fields on fractals[END_REF], then the variance displays log-periodic oscillations together with a power-law envelope. To illustrate these results, we consider three examples of quasiperiodic potentials, for which we compute exactly the hyperuniformity exponent and the period of these oscillations, when it exists.

Variance of the LLL Magnetic field 0 FIG. 1. Sketch of the LLL variance w 2 as a function of the magnetic field B for different perturbing potentials. Red: Disordered w 2 ∼ B (see Ref. [START_REF] Ando | Theory of quantum transport in a two-dimensional electron system under magnetic fields. I. Characteristics of level broadening and transport under strong fields[END_REF]). Green: Hyperuniform with discrete scale invariance w 2 ∼ B 2+α 2 + log-periodic oscillations, where α is the hyperuniformity exponent (this work). Blue: Periodic w 2 ∼ e -#/B (see Ref. [START_REF] Rauh | On the broadening of landau levels in crystals[END_REF]).

II. LANDAU LEVELS PERTURBED BY A POTENTIAL

To begin with, let us recall a few well-known results. The Hamiltonian describing a particle of mass m and charge e in a magnetic field B = ∇ × A is given by

H 0 = (p -eA) 2 2m . (1)
Here, we consider a two-dimensional system with a magnetic field perpendicular to the plane. Such a field can be described by the symmetric gauge where the vector potential reads A = B(-y/2, x/2, 0). The spectrum of H 0 consists in equidistant energy levels known as Landau levels,

E n = hω c (n + 1/2), ∀n ∈ N, (2)
where ω c = |eB|/m is the cyclotron frequency. Each Landau level has a degeneracy proportional to the sample area A and the magnetic field. In the following, we set h = e = 1.

Our aim is to study the behavior of the Landau levels in the presence of a time-independent potential. Thus, we consider the following general Hamiltonian:

H = H 0 + V (x, y), ( 3
)
and we assume that the magnitude of the potential is small compared to the Landau level spacing ω c |V |. In this regime, we can neglect the coupling between different Landau levels and use degenerate perturbation theory to compute the degeneracy splitting of a single level. Without loss of generality, we assume V (x, y) 0 and B > 0 in the following.

III. VARIANCE OF THE LLL

For simplicity, in the following we focus on the LLL corresponding to n = 0 for which the nonperturbed wave functions, in the thermodynamic limit, can be chosen as

ϕ l (z) = x, y|l = 1 2π l 2 B l! 2 l z l e -|z| 2 /4 , (4)
where

z = (x + iy)/l B , l = 0, 1, . . . , N φ -1 is the angular momentum, N φ = A/(2π l 2 B )
1 is the degeneracy of the LLL, and l B = 1/ √ B is the magnetic length. To characterize the broadening of the LLL due to the potential, we consider its variance defined by

w 2 = 1 N φ N φ -1 p=0 ε 2 p - ⎛ ⎝ 1 N φ N φ -1 p=0 ε p ⎞ ⎠ 2 , ( 5
)
where ε p 's are eigenenergies of H projected onto the LLL. This variance can be recast as

w 2 = 1 N φ N φ -1 l=0 N φ -1 l =0 | l|V |l | 2 - ⎛ ⎝ 1 N φ N φ -1 l=0 l|V |l ⎞ ⎠ 2 , (6)
so that one does not need to compute explicitly the ε p 's. Setting r = (x, y) = r(cos θ, sin θ ), a matrix element of the perturbation potential in the LLL basis {|l , l = 0, . . . ,

N φ -1} reads l|V |l = dq (2π ) 2 V (q) 2π √ l! l ! 2 l+l × ∞ 0 dr l B r l B 1+l+l e -r 2 2l 2 B 2π 0 dθ e iq•r e iθ (l-l ) , ( 7
)
where we introduced the Fourier transform of the potential

V (q) = dr e -iq•r V (r). (8)
In the large-N φ (thermodynamical) limit, one then gets

∞ l=0 l|V |l = V (0) 2π l 2 B , ( 9
) ∞ l=0 ∞ l =0 | l|V |l | 2 = 1 2π l 2 B dq (2π ) 2 | V (q)| 2 e -|q| 2 l 2 B /2 . ( 10
)
Finally, one obtains the following expression for the variance:

w 2 = dq (2π ) 2 S(q)e -|q| 2 l 2 B /2 , ( 11
)
where we introduced the structure factor

S(q) = | V (q)| 2 A (1 -δ q,0 ). (12)
Note that the term proportional to δ q,0 comes from the second term of Eq. ( 5) and is irrelevant only if V (0) = 0. The variance is therefore essentially equal to the integral of the structure factor over a disk of radius l -1 B , which is the main result of this paper. Before discussing the most interesting case of a potential with a dense Fourier spectrum, let us first show that this expression allows one to recover known results for simple potentials.

IV. PERIODIC POTENTIAL

For a periodic potential of strength V 0 with a single spatial frequency a -1 ,

V (x, y) = V 0 [cos(2π x/a) + cos(2π y/a)], (13)
Eq. ( 11) leads to

w 2 = V 2 0 e -2π 2 Ba 2 , ( 14
)
in agreement with the expression found by Rauh [START_REF] Rauh | On the broadening of landau levels in crystals[END_REF] (see Appendix A for details). The generalization to Fourier spectra with a finite set of frequencies is straightforward, even if the potential is no longer periodic. In the zero-field limit, the LLL broadening is exponential and controlled by the smallest frequency. The case of a dense set of frequencies is more subtle.

V. RANDOM POTENTIAL

Landau level broadening due to an uncorrelated random potential has been widely studied in the literature [START_REF] Huckestein | Scaling theory of the integer quantum Hall effect[END_REF]. For the simple case of a random potential with zero mean and whitenoise correlations,

V (r) = 0, (15)
V (r)V (r ) = (V 0 a) 2 δ(r -r ), (16)
where the overline denotes the average over disorder realizations, Eq. [START_REF] Wilkinson | Critical properties of electron eigenstates in incommensurate systems[END_REF] gives

w 2 = V 2 0 a 2 2π l 2 B = V 2 0 Ba 2 2π , ( 17
)
in agreement with the result of Ando [START_REF] Ando | Theory of quantum transport in a two-dimensional electron system under magnetic fields. III. Many-site approximation[END_REF] (see Appendix B for details). This result is very different from the one obtained for a potential with a finite number of frequencies discussed above.

For stealthy hyperuniform disorder [START_REF] Torquato | Hyperuniform states of matter[END_REF], the structure factor is identically zero in a disk of radius q 0 > 0 around the origin. A reasonable approximation is to assume that S(q) ∝ (|q|q 0 ) leading to a LLL broadening,

w 2 ∝ B e -q 2 0 2B , ( 18
)
intermediate between that of a periodic potential, Eq. ( 14), and that of uncorrelated random disorder, Eq. ( 17).

VI. POTENTIAL WITH A DENSE FOURIER SPECTRUM

The most interesting situation comes from potentials with a dense Fourier spectrum made of Bragg peaks (as found, e.g, in quasicrystals). To this end, let us consider a general potential

V (r) = V 0 a 2 N j=1 δ(r -r j ), (19)
built on a set of N scattering points located at position r j with a typical density a -2 . The random potential discussed above belongs to this family. Before proceeding further, let us stress that the exponential term in Eq. ( 11) acts as a smooth cutoff that eliminates wave vectors |q| l -1 B . To analyze the behavior of w 2 , we shall instead consider a sharp cutoff regularization by introducing the integrated intensity function

Z (k) = |q|<k dq S(q), ( 20
)
so that one has

w 2 ∼ Z (k ∼ 1/l B ). ( 21
)
This approximation clearly misses exponentially small terms so that, for the periodic case discussed previously [see Eq. ( 14)], it leads to w 2 = 0. Let us remember that Eq. ( 21) is valid in the perturbative regime where mV 0 B. In the following, we further focus on the case in which B 1/a 2 since, for many potentials, Z has a simple behavior in the k ∼ 1/l B 1/a limit. The integrated intensity function (also known as the spectral measure [START_REF] Luck | Cantor spectra and scaling of gap widths in deterministic aperiodic systems[END_REF]) is commonly used to analyze sets of points with a nontrivial structure factor [START_REF] Oguz | Hyperuniformity of quasicrystals[END_REF]. In one-dimensional quasicrystals, Z is conjectured to have a power-law envelope for k 1/a [START_REF] Luck | Cantor spectra and scaling of gap widths in deterministic aperiodic systems[END_REF][START_REF] Oguz | Hyperuniformity of quasicrystals[END_REF][START_REF] Oguz | Hyperuniformity and anti-hyperuniformity in one-dimensional substitution tilings[END_REF]. As we shall see, this is also the case

for two-dimensional quasicrystals. Assuming Z (k) ∼ k→0 k 2+α ,
Eq. ( 21) leads to

w 2 ∼ B 2+α 2 for mV 0 B 1/a 2 , ( 22
)
establishing a relation between the broadening of the LLL and the so-called hyperuniformity exponent α that characterizes the potential. For α > 0, the potential is hyperuniform [START_REF] Oguz | Hyperuniformity of quasicrystals[END_REF], whereas α < 0 refers to hypo-uniformity (or antihyperuniformity [START_REF] Oguz | Hyperuniformity and anti-hyperuniformity in one-dimensional substitution tilings[END_REF]). The special case α = 0 corresponds to a potential with a constant S, such as the random potential considered previously [see Eq. ( 17)]. Interestingly, if Z further manifests a discrete scale invariance, i.e., if there exists λ > 1 such that

Z (k/λ) = Z (k)/λ 2+α , ( 23
)
then one has

Z (k) = k 2+α F (ln k/ ln λ), ( 24
)
where

F (x + 1) = F (x) (see

the examples below and

Refs. [START_REF] Sornette | Discrete-scale invariance and complex dimensions[END_REF][START_REF] Akkermans | Statistical mechanics and quantum fields on fractals[END_REF] for a review). As a result, the LLL variance w 2 displays log-periodic oscillations together with a power-law envelope in the small-B limit.

VII. EXAMPLES OF QUASIPERIODIC POTENTIALS

For illustration, let us consider some potentials of the form given in Eq. ( 19) where the points correspond to vertices of two-dimensional quasiperiodic tilings (see Appendix C). For each tiling considered below, we computed exactly the structure factor S, the hyperuniformity exponent α characterizing the power-law behavior of Z (k) ∼ k→0 k 2+α , and the discrete scale invariance factor λ defined in Eq. ( 23) when it exists. Numerical results displayed in Fig. 2 have been obtained by integrating more and more Bragg peaks of smaller and smaller intensities. In each case, we checked that the results were converged in the range considered. Units are taken such that V 0 = A/(Na 2 ) and a = 1, where a is the edge length of the hypercubic lattice that is used to build the tiling in the standard cut-and-project method [START_REF] De Bruijn | Algebraic theory of Penrose's non-periodic tilings of the plane. I, Indag[END_REF][START_REF] Duneau | Quasiperiodic Patterns[END_REF][START_REF] Kalugin | Al 0.86 Mn 0.14 : A six-dimensional crystal[END_REF][START_REF] Elser | Indexing problems in quasicrystal diffraction[END_REF].

Let us first consider the twofold-symmetric Rauzy tiling [START_REF] Vidal | Generalized quasiperiodic Rauzy tilings[END_REF]. The hyperuniformity exponent is α = 4, but Z has no discrete scale invariance (see Appendix D). By contrast, for the eightfold-symmetric Ammann-Beenker tiling [START_REF] Beenker | Algebraic theory of non periodic tilings of the plane by two simple building blocks: A square and a rhombus[END_REF][START_REF] Grünbaum | Tilings and Patterns[END_REF][START_REF] Duneau | Approximants of quasiperiodic structures generated by the inflation mapping[END_REF], the hyperuniformity exponent is α = 2, and Z has a discrete scale invariance with λ = 1 + √ 2 (see Appendix F). For the fivefold-symmetric Penrose tilings, the hyperuniformity exponent is α = 6, and Z has a discrete scale invariance with λ = τ 2 , where τ

= 1+ √ 5 2
is the golden ratio (see Appendix G).

VIII. SUBSTITUTION TILINGS AND DISCRETE SCALE INVARIANCE

As recently conjectured by Oguz et al. [START_REF] Oguz | Hyperuniformity and anti-hyperuniformity in one-dimensional substitution tilings[END_REF], the behavior of Z in one-dimensional substitution tilings is determined by the eigenvalues of the substitution matrix. More precisely, for nonperiodic binary substitutions associated with a 2 × 2 substitution matrix with eigenvalues when k tends to zero, so that

λ 1 > |λ 2 | > 0, one has Z (k/λ 1 ) = Z (k)(λ 2 /λ 1 ) 2 (25)
Z (k) = k 1+α F (ln k/ ln λ 1 ), α = 1 -2 ln |λ 2 | ln λ 1 , (26)
with F (x + 1) = F (x). In two dimensions, it is very likely that the existence of substitution rules with inflation/deflation also implies the existence of discrete scale invariance for Z. This is clearly the case for the Ammann-Beenker and Penrose tilings, which, contrary to the Rauzy tiling, can be built by inflation/deflation. However, we have not found a simple expression for the hyperuniformity exponent [such as the one given Eq. ( 26)] for two-dimensional binary substitutions.

IX. OUTLOOK

In this paper, we obtained a simple relation between the Landau level broadening and the integrated intensity function Z of the perturbing potential. For potentials with a dense Fourier spectrum made of Bragg peaks, this relation implies that the variance of the LLL is driven by the hyperuniformity exponent α [see Eq. ( 22)]. In the absence of a complete classification of the possible behavior of Z, a first step to go beyond would consist in analyzing two-dimensional potentials with a singular continuous Fourier spectrum for which one expects more complex behavior of Z as observed in one dimension [START_REF] Luck | Cantor spectra and scaling of gap widths in deterministic aperiodic systems[END_REF][START_REF] Oguz | Hyperuniformity and anti-hyperuniformity in one-dimensional substitution tilings[END_REF]. For instance, one may find noninteger exponents α or even nonalgebraic decay. Another important issue would be to consider the influence of Landau level mixing, which is known to have dramatic effects on the localization properties of the eigenstates [START_REF] Haldane | Landau Level Mixing and Levitation of Extended States in Two Dimensions[END_REF], and hence on integer quantum Hall physics (see, e.g., Ref. [START_REF] Kivelson | Global phase diagram in the quantum Hall effect[END_REF]). It would also be important to bridge the gap between the perturbed free-electron results and the one obtained numerically in tight-binding models [START_REF] Fuchs | Landau levels in quasicrystals[END_REF]. A possible route would be to develop the analog of Wilkinson semiclassical treatment [START_REF] Wilkinson | Critical properties of electron eigenstates in incommensurate systems[END_REF] for nonperiodic potential.

Finally, let us mention that the magnetic-field dependence of the Landau level broadening induced by disorder has already been measured in graphene [START_REF] Orlita | Approaching the Dirac Point in High-Mobility Multilayer Epitaxial Graphene[END_REF]. Combining such an experimental device with a nontrivial superlattice potential would allow us to measure the behaviors discussed in the present work.

=

V 2 0 4 A δ q x ,2π/a δ q y ,0 + δ q x ,-2π/a δ q y ,0 + δ q x ,0 δ q y ,2π/a + δ q x ,0 δ q y ,-2π/a , (A4) 2 [δ(q x -2π/a)δ(q y ) + δ(q x + 2π/a)δ(q y ) + δ(q x )δ(q y -2π/a) + δ(q x )δ(q y + 2π/a)], (A5)

= V 2 0 4 (2π )
where we used the Kronecker delta in order to compute the modulus square of the Fourier transform (otherwise, the square of a Dirac δ is ill-defined). Then Eq. ( 11) involves the integral of Dirac δ functions, which straightforwardly leads to Eq. ( 14).

APPENDIX B: VARIANCE FOR A RANDOM POTENTIAL

The Fourier transform of the uncorrelated random potential defined by Eqs. ( 15) and ( 16) is given by

V (q) = dr e -ir•q V (r) = 0 and | V (q)| 2 = dr e -ir•q dr e ir •q V (r)V (r ) = (V 0 a) 2 A.
The structure factor [START_REF] Fuchs | Landau levels in quasicrystals[END_REF] becomes S(q) = (V 0 a) 2 (1 -δ q,0 ). In Eq. ( 11), the Kronecker delta does not contribute, and the Gaussian integral gives Eq. ( 17):

w 2 = (V 0 a) 2 (2π ) 2 dq e -|q| 2 l 2 B /2 = (V 0 a) 2 (2π ) 2 2π l 2 B .

APPENDIX C: FOURIER TRANSFORM OF A CUT-AND-PROJECT QUASICRYSTAL

The cut-and-project (CP) method [START_REF] De Bruijn | Algebraic theory of Penrose's non-periodic tilings of the plane. I, Indag[END_REF][START_REF] Duneau | Quasiperiodic Patterns[END_REF][START_REF] Kalugin | Al 0.86 Mn 0.14 : A six-dimensional crystal[END_REF][START_REF] Elser | Indexing problems in quasicrystal diffraction[END_REF] consists in selecting points of a D-dimensional periodic lattice if their projection onto the (Dd )-dimensional perpendicular space E ⊥ belongs to a so-called acceptance window. The tiling is then obtained by projecting these selected points onto the complementary d-dimensional parallel space E .

Any vector v in hyperspace can be decomposed uniquely in terms of its projection onto parallel and perpendicular spaces as

v = v + v ⊥ . (C1)
As explained in the early papers introducing the CP method [START_REF] Duneau | Quasiperiodic Patterns[END_REF][START_REF] Kalugin | Al 0.86 Mn 0.14 : A six-dimensional crystal[END_REF][START_REF] Elser | Indexing problems in quasicrystal diffraction[END_REF], the Fourier transform of quasiperiodic tilings can be computed from the higher-dimensional space from which it stems. The main idea is that since points of the tiling are selected from a periodic tiling via an acceptance window, computing the Fourier transform of the tiling essentially amounts to computing the Fourier transform of this acceptance window.

For a tiling with N sites (vertices) at position R j and obtained by the CP method, the microscopic density is

n(r ) = N j=1 δ(r -R j ), ( C2
)
and its Fourier transform is

n(q ) = N j=1 e -i q •R j , ( C3
)
where the sum runs over all sites of the d-dimensional tiling considered. The convention that we use is that upper-case letters (such as R j ) refer to discrete points, and lower-case letters (such as r ) refer to a continuum of points.

Let R be a point of the D-dimensional hypercubic lattice, and let K be a vector of its reciprocal lattice such that K • R = 2π × integer. These vectors can be decomposed onto the parallel and perpendicular spaces such that their scalar product reads

K • R = K • R + K ⊥ • R ⊥ = 2π × integer. (C4)
Equation ( C3) is nonzero iff q = K , in which case it becomes

n(K ) = N j=1 e i K ⊥ •R ⊥ j . ( C5
)
For a quasicrystal built along an irrational plane (parallel space), the points in perpendicular space densely and uniformly fill the acceptance window such that

n(K ) = N A ⊥ dr ⊥ A ⊥ e i K ⊥ •r ⊥ , ( C6
)
where the integral is over the acceptance window in perpendicular space, and A ⊥ is its (Dd )-dimensional volume. Now, for any vector q in parallel space, the Fourier transform of the density Eq. (C3) reads

n(q ) = K δ q ,K N A ⊥ dr ⊥ A ⊥ e i K ⊥ •r ⊥ , (C7)
where the sum is performed over all vectors K of the reciprocal lattice of the hypercubic lattice. As we are considering a quasicrystal, for any K there is a unique K and therefore K ⊥ is well-defined. If {a * j ; j = 1, . . . , D} is a basis of vectors in reciprocal space, then K = j n j a * j , where n j are integers. Its parallel and perpendicular components are also functions of the same integers:

K = K (n 1 , . . . , n D ) + K ⊥ (n 1 , . . . , n D ). (C8)
Therefore, the sum over K in Eq. (C7) is actually a sum over D integers n 1 , . . . , n D , clearly showing that the Fourier transform is a pure point of rank D > d.

Let us define the structure factor in the thermodynamic limit as

S(q ) = | n(q )| 2 N (1 -δ q ,0 ). ( C9
)
For two-dimensional potentials (d = 2) of the form given by Eq. ( 19), this definition differs from the one given in Eq. ( 12) by a factor V 2 0 a 4 N/A, which disappears upon choosing units such that V 0 = A/(Na 2 ).

As explained in Ref. [START_REF] Oguz | Hyperuniformity of quasicrystals[END_REF], for a spectrum made of a dense set of Bragg peaks (discontinuous S), the integrated intensity function

Z (k) = |q |<k S(q ) dq (C10)
provides a reliable characterization of the point distribution.

Here, the integral is performed over a disk of radius k. This function is also known as the spectral measure in Ref. [START_REF] Luck | Cantor spectra and scaling of gap widths in deterministic aperiodic systems[END_REF].

For d-dimensional tilings built by the CP method, this quantity can be recast in the following form:

Z (k) = (2π ) d A |K |<k S(K ). ( C11
)
APPENDIX D: THE RAUZY TILING

Fourier transform

The two-dimensional (generalized) Rauzy tiling has been introduced in Ref. [START_REF] Vidal | Generalized quasiperiodic Rauzy tilings[END_REF]. This is a codimension-1 tiling built from the cubic lattice Z 3 (edge length a = 1) with a onedimensional perpendicular space oriented along the direction e ⊥ = (θ 2 , θ, 1), where θ is the real (Pisot-Vijayaraghavan) root of the cubic equation x 3 = x 2 + x + 1. Contrary to the Ammann-Beenker and the Penrose tilings discussed in the next Appendixes, the Rauzy tiling cannot be built by substitution rules.

For such a codimension-1 quasicrystal, the acceptance window is a segment of length A ⊥ defined as the projection of h = (1, 1, 1) onto the perpendicular space. This acceptance window also corresponds to the projection of the unit cube onto the perpendicular space. In this case, Eq. (C6) gives

| n(K )| = N sinc K ⊥ • h ⊥ 2 . (D1)

Structure factor

The structure factor is defined in Eq. (C9). Our goal is to analyze the behavior of S(q ) in the limit where |q | tends to zero. By definition, one has S(0) = 0, but its behavior for small |q | is nontrivial since S(q ) = 0 only when q coincides with the parallel component K of a reciprocallattice vector K of the cubic lattice. Thus, we are interested in computing the behavior of S when |K | goes to 0 for K = 0:

S(K ) = N sin K ⊥ •h ⊥ 2 K ⊥ •h ⊥ 2 2 , (D2) = N sin K •h 2 K ⊥ •h ⊥ 2 2 , ( D3
) |K |→0 N K • h K ⊥ • h ⊥ 2 , ( D4
)
where we used the fact that K is a reciprocal-lattice vector and h is a direct-lattice vector. The difficulty comes from the fact that, when |K | goes to 0, |K ⊥ | diverges. So, the goal is to find the relation between these two components. One way to investigate this issue is to follow the approach proposed in Ref. [START_REF] Oguz | Hyperuniformity of quasicrystals[END_REF] for the Fibonacci chain (see Appendix E). In the Z 3 canonical basis, any reciprocal-lattice vector K has coordinates 2π (l, m, n), where l, m, and n are integers. To analyze the behavior of K • h and K ⊥ • h ⊥ , let us consider the matrix

M = ⎛ ⎝ 1 1 1 1 0 0 0 1 0 ⎞ ⎠ , ( D 5 )
which satisfies

M 3 = M 2 + M + 1.
The eigenvalues of M are the Tribonacci constant θ 1.8393 and two complex conjugate eigenvalues e ±iφ / √ θ with φ 2.1762. The eigenvector associated with θ corresponds to the perpendicular direction e ⊥ . Since θ is a Pisot-Vijayaraghavan number, the action of M p onto any vector v, such that v • e ⊥ = 0, drives this vector toward the direction e ⊥ in the large-p limit.

Hence, to analyze the behavior of S in the limit where K tends to zero [see Eq. (D5)], let us consider K (p) = M p K. More precisely, we are interested in computing K (p) • h and K (p) ⊥ • h ⊥ . Keeping in mind that h ⊥ = P ⊥ (1, 1, 1) and h = (1 -P ⊥ )(1, 1, 1) (where P ⊥ is the projector onto the perpendicular space), one can easily compute these quantities. After some algebra, one gets

K (p) • h = (θ -1)θ -p+1 2 {sin(p φ)[(1 + θ 2 )m -θ (l + n)] + √ θ {(θ m -l ) sin[(1 + p)φ] + (θ n -m) sin[(1 -p)φ]}} sin(φ)[(θ -1)θ + 1] ,
(D6)

K (p) ⊥ • h ⊥ = 1 sin(φ)[(θ -1)θ + 1][2θ 3/2 cos(φ) -θ 3 -1] {-θ p (θ 4 + θ 2 + 1) sin(φ)(θ l -2 √ θm cos(φ) + n) + θ -p-1 2 [θ 3/2 sin[(p + 1)φ][(θ -1)l + θ 2 (n -m)] + (θ -1)θ sin(p φ)[l -(θ + 1)m + θ n] + √ θ sin[(p -1)φ][-l + m + (θ -1)θ n] + θ 3 (l -θ m) sin[(p + 2)φ] + (m -θ n) sin[(p -2)φ]]}. (D7)
Thus, in the large-p limit, one finds that K (p) • h vanishes as θ -p/2 , K (p) ⊥ • h ⊥ diverges as θ p , and S(K (p) ) behaves as θ -3p . As a result, one finds that

S(K ) ∼ |K |→0 |K | 6 (D8)
for all (l, m, n). However, we emphasize that, contrary to the Fibonacci chain Ref. [START_REF] Oguz | Hyperuniformity of quasicrystals[END_REF] (see also Appendix E), this power-law behavior is modulated by a bounded oscillating nonperiodic function, as can be seen in Eqs. (D6) and (D7).

Integrated intensity function

The integrated intensity function is defined in Eq. (C11). The sum over all vector K with a norm smaller than k can be decomposed into a sum over all triplets (l, m, n) and their iterated under M p . As a result, one has

Z (k) = 4π 2 A (l,m,n) ∞ p=p (l,m,n) S(K (p) ), ( D9
)
where p (l,m,n) is the smallest integer fulfilling the constraint

|K (p) | < k.
As already discussed in the previous section, in the large-p limit,

S(K (p) ) |K (p) | 6 f (l, m, n, p), ( D10
)
|K (p) | θ -p/2 g(l, m, n, p), (D11)
where, for a given triplet (l, m, n), f and g are bounded oscillating functions of p [see Eqs. (D6) and (D7)]). Thus, S is bounded both above and below,

Z -(k) Z (k) Z + (k), (D12)
where 3 , as can be seen from Eqs. (D10) and (D11), since dividing k by √ θ simply amounts to changing p (l,m,n) into p (l,m,n) + 1 in Eq. (D13). Such a relation reflects a discrete scale invariance [START_REF] Sornette | Discrete-scale invariance and complex dimensions[END_REF] (see also the next Appendix) for Z ± and implies a power-law envelope

Z ± (k) = 4π 2 A (l,m,n) c ± (l, m, n) ∞ p=p (l,m,n) θ -3p , (D13) c + (l, m, n) = max p f (l, m, n, p) 6 g(l, m, n, p), (D14) c -(l, m, n) = min p f (l, m, n, p) 6 g(l, m, n, p). (D15) Interestingly, Z ± (k/ √ θ ) = Z ± (k)/θ
Z (k) ∼ k→0 k 6 . (D16)
Note that, despite the fact that Z is defined as an integral of S, they are both characterized by a power law with the same exponent. This is a consequence of the fact that S is discontinuous (discrete) and dense.

APPENDIX E: INTEGRATED INTENSITY FUNCTION OF THE FIBONACCI CHAIN

The Fibonacci chain is a one-dimensional tiling built from the square lattice Z 2 (edge length a = 1). The integrated intensity function Z of the Fibonacci chain has been widely discussed in Ref. [START_REF] Oguz | Hyperuniformity of quasicrystals[END_REF]. However, one important property has been missed. As a codimension-1 system, the Fourier transform of the Fibonacci chain can be easily computed. The structure factor is

S(K ) = N sin K .h 2 K ⊥ .h ⊥ 2 2 , ( E1
)
where h and h ⊥ are the projections of the vector

h = (1, 1) onto e = 1 √ 1+τ 2 (-1, τ ) and e ⊥ = 1 √ 1+τ 2 (τ, 1), where τ = 1+ √ 5 2
is the golden ratio. The integrated intensity function is then

Z (k) = 2π A |K |<k S(K ), ( E2
)
where A is the total length of the chain. As for the Rauzy tiling, let us consider the matrix

M = 1 1 1 0 , ( E 3 )
which satisfies M 2 = M + 1. Eigenspaces of M correspond to the perpendicular and parallel directions with eigenvalues τ and -1/τ , respectively. The small-k behavior of Z is obtained by analyzing sequences K (p) = M p K (see Ref. [START_REF] Oguz | Hyperuniformity of quasicrystals[END_REF]). One then gets, in the large-p limit,

S(K (p) ) |K (p) | 4 f (l, m), (E4)
|K (p) | τ -p g(l, m). (E5)
However, contrary to the Rauzy tiling, f and g do not depend on p. Thus, following the same line of reasoning as above, one straightforwardly gets the discrete scaling relation

Z (k/τ ) = Z (k)/τ 4 . (E6)
The solution of this equation can be written as

Z (k) = k 4 F (ln k/ ln τ ), (E7)
where F (x + 1) = F (x) (for a review on discrete scale invariance, see Ref. [START_REF] Sornette | Discrete-scale invariance and complex dimensions[END_REF]). As a result, Z has a power-law envelope together with log-periodic oscillations (see Fig. 3 for illustration). This is in stark contrast with the Rauzy tiling where only Z + and Z -obey such a discrete scale invariance but not Z itself. Practically, to compute Z, we first select a set of K points in the reciprocal lattice of Z 2 inside a given ball of radius K max around the origin. For each of these points, we consider the sequence of points K (p) with p = 0, . . . , p max , and we compute S for each corresponding K (p) (avoiding possible redundancy). Z is then obtained by summing over these Bragg peaks according to Eq. (E2). We check the convergence of the results displayed in Fig. 3 by increasing K max and p max .

APPENDIX F: THE OCTAGONAL TILING

Fourier transform

The octagonal (Ammann-Beenker) tiling [START_REF] Beenker | Algebraic theory of non periodic tilings of the plane by two simple building blocks: A square and a rhombus[END_REF][START_REF] Grünbaum | Tilings and Patterns[END_REF][START_REF] Duneau | Approximants of quasiperiodic structures generated by the inflation mapping[END_REF] is a codimension-2 tiling built from the four-dimensional hypercubic lattice Z 4 (edge length a = 1). Perpendicular and parallel spaces are spanned by the eigenvectors of the matrix,

M = ⎛ ⎜ ⎝ 1 0 -1 1 0 1 -1 -1 -1 -1 1 0 1 -1 0 1 ⎞ ⎟ ⎠, (F1)
associated with eigenvalues λ ± = 1 ± √ 2. This matrix satisfies M 2 = 2M + 1, and its eigenvalues are twofolddegenerate. Here, we choose the following orthonormal eigenbasis:

e ,1 = - 1 2 , 1 2 , 0, 1 √ 2 , e ,2 = 1 2 , 1 2 , 1 √ 2 , 0 , (F2) e ⊥,1 = 1 2 , - 1 2 , 0, 1 √ 2 , e ⊥,2 = - 1 2 , - 1 2 , 1 √ 2 , 0 ,
where the perpendicular (parallel) space is associated with λ + (λ -). The acceptance window is an octagon corresponding to the projection of the four-dimensional unit hypercube onto the perpendicular space. In this case, Eq. (C6) gives

| n(K )| = N λ + cos λ + K ⊥,1 -K ⊥,2 2 K ⊥,2 (K ⊥,1 +K ⊥,2 ) - cos λ + K ⊥,1 +K ⊥,2 2 K ⊥,2 (K ⊥,1 -K ⊥,2 ) + cos λ + K ⊥,2 -K ⊥,1 2 K ⊥,1 (K ⊥,2 +K ⊥,1 ) - cos λ + K ⊥,2 +K ⊥,1 2 K ⊥,1 (K ⊥,2 -K ⊥,1 ) , (F3)
for all reciprocal-lattice vectors K with components K , j = K • e , j and K ⊥, j = K • e ⊥, j . These expressions coincide with the one given in Ref. [START_REF] Janssen | Aperiodic Crystals[END_REF].

Structure factor

We are interested in computing the behavior of S when |K | goes to 0 for K = 0. To that end, we note that for a vector K = 2π (s, t, u, v) [where (s, t, u, v) ∈ Z 4 ], one has

λ + K ⊥,1 = K ,1 λ + + 2π (s -t + v), (F4) λ + K ⊥,2 = K ,2 λ + + 2π (-s -t + u), (F5) K ⊥,1 + K ,1 = 2π v, (F6) K ⊥,2 + K ,2 = 2π u. (F7)
A close inspection of Eq. (F3) shows that one has to distinguish three different cases.

a. Symmetry axes: S(K )∼|K | 4

As can be seen in Eq. (F3), the denominator vanishes if one of the components K ⊥,i = 0 or when K ⊥,1 = ±K ⊥,2 . When K ⊥ belongs to these four symmetry axes, the Fourier transform can be recast in a simple form. For simplicity, let us focus on the case in which K ⊥,2 = 0 (the other cases being treated similarly), for which

| n(K )| = N λ + K 2 ⊥,1 2 cos K ⊥,1 2 -2 cos λ + K ⊥,1 2
+K ⊥,1 sin λ + K ⊥,1 2 . ( F8
)
Using Eqs. (F4)-(F6), one then obtains

S(K ) |K |→0 N 4λ 4 + K ,1 K ⊥,1 2 . ( F9
)
As was done previously, to analyze the behavior of the structure factor for small |K |, we consider K (p) = M p K. By construction, in the large-p limit, the parallel components of K (p) tend to zero as λ p -and its perpendicular components diverge as λ p + . As a result, S(K (p) ) behaves as λ -4p + so that, in this case,

S(K ) ∼ |K |→0 |K | 4 . (F10)
This result actually holds for all K ⊥ belonging to the four symmetry axes discussed above.

b. Generic cases: S(K )∼|K | 8 or S(K )∼|K | 12

When K ⊥ does not belong to the four symmetry axes defined as K ⊥,1 = 0, K ⊥,2 = 0, and K ⊥,1 = ±K ⊥,2 , one can again use Eqs. (F4)-(F7) to express the structure factor as

S(K ) |K |→0 N 4λ 4 + (K ,1 K ⊥,2 +K ,2 K ⊥,1 )(K ,1 K ⊥,1 -K ,2 K ⊥,2 ) K ⊥,1 K ⊥,2 (K ⊥,1 +K ⊥,2 )(K ⊥,1 -K ⊥,2 ) 2 , (F11)
which leads to

S(K ) ∼ |K |→0 |K | 8 . (F12)
However, as can be seen in Eq. (F11), this leading contribution may vanish for some special K. In this case, S is given by the subleading contribution, which gives

S(K ) ∼ |K |→0 |K | 12 .
(F13)

Integrated intensity function

To compute the integrated intensity function defined in Eq. (C11), we decompose the sum over all vector K with a norm smaller than k as a sum over all quadruplets (s, t, u, v) and their iterated vectors K (p) = M p K . One can then write However, since we are interested in the small-k (large-p) limit, one only keeps the dominant terms in Eq. (F14) that come from the symmetry axes and gives

Z (k) = 4π 2 A (s,t,u,v) ∞ p=p (s,t,u,v) S(K (p) ), (F14)
S(K (p) ) |K (p) | 4 f (s, t, u, v), (F15) |K (p) | λ -p + g(s, t, u, v). (F16)
We emphasize that, as for the Fibonacci chain, f and g are functions that do not depend on p, so that one straightforwardly gets the following discrete scaling relation:

Z (k/λ + ) = Z (k)/λ 4 + . (F17)
The solution of this equation can be written as

Z (k) = k 4 F (ln k/ ln λ + ), (F18)
where F (x + 1) = F (x). As a result, Z has a power-law envelope together with log-periodic oscillations.

APPENDIX G: THE PENROSE TILING

The Penrose rhombus tiling [START_REF] De Bruijn | Algebraic theory of Penrose's non-periodic tilings of the plane. I, Indag[END_REF][START_REF] Penrose | Pentaplexity: A class of non-periodic tilings of the plane[END_REF] can be built by CP from the five-dimensional hypercubic lattice Z 5 (edge length a = 1) along a well-known procedure (see, e.g., Ref. [START_REF] Duneau | Lectures on Quasicrystals[END_REF] for details). For our purpose, let us consider the following orthogonal (non-normalized) basis:

e ,1 = 2 5 (1, c 2 , c 4 , c 4 , c 2 ), e ,2 = 2 5 (0, s 2 , s 4 , -s 4 , -s 2 ), e ⊥,1 = 2 5 (1, c 4 , c 2 , c 2 , c 4 ), e ⊥,2 = 2 5 (0, s 4 , -s 2 , s 2 , -s 4 ), e = 1 10 (1, 1, 1, 1, 1), (G1)
which defines the three subspaces E , E ⊥ , and . Here, we introduced the notation c n = cos(2π n/5) and s n = sin(2π n/5). A point in Z 5 is selected whenever it projects onto the perpendicular space E ⊥ + inside a three-dimensional acceptance window which is the projection of the five-dimensional unit hypercube onto this subspace. Remarkably, selected points only fill five planes perpendicular to . Thus, the selection step only amounts to considering discrete sections of the acceptance window. Among all possible choices, the fivefoldsymmetric canonical Penrose tilings (known as star and sun [START_REF] Grünbaum | Tilings and Patterns[END_REF]) considered here correspond to the following sections: One point that is the symmetry center of the tiling, two regular pentagons of side 2 √ 2/5 cos(3π/10), and two regular pentagons of side 2 τ √ 2/5 cos(3π/10), where τ = 1+ √ 5 2 is the golden ratio.

Fourier transform

The Fourier transform of the tiling's vertices is obtained as a weighted sum of the Fourier transform of the four regular pentagons. For any vector R of the five-dimensional hypercubic lattice, the five-dimensional reciprocal lattice vectors K satisfy

K • R = K • R + K ⊥ • R ⊥ + K • R = 2π × integer. (G2)
Then, after some algebra, Eq. (C6) leads to

|ñ(K )| = N 8(s 2 -s 4 ) 5K ⊥,2 cos(c 2 K ⊥,1 + s 2 K ⊥,2 -3K ) + cos[K ⊥,1 /2 + (s 4 + s 2 )K ⊥,2 + K ] -cos(K ⊥,1 -3K ) -cos(2c 4 K ⊥,1 -K ) 5K ⊥,1 -(6s 2 + 2s 4 )K ⊥,2 - cos(c 2 K ⊥,1 -s 2 K ⊥,2 -3K ) + cos[K ⊥,1 /2 -(s 4 + s 2 )K ⊥,2 + K ] -cos(K ⊥,1 -3K ) -cos(2c 4 K ⊥,1 -K ) 5K ⊥,1 + (6s 2 + 2s 4 )K ⊥,2 + cos(c 2 K ⊥,1 -s 2 K ⊥,2 -3K ) + cos[K ⊥,1 /2 -(s 4 + s 2 )K ⊥,2 + K ] -cos(c 4 K ⊥,1 -s 4 K ⊥,2 -3K ) -cos[(1 + c 2 )K ⊥,1 + s 2 K ⊥,2 -K ] 5K ⊥,1 -(6s 4 -2s 2 )K ⊥,2 - cos(c 2 K ⊥,1 + s 2 K ⊥,2 -3K ) + cos[K ⊥,1 /2 + (s 4 + s 2 )K ⊥,2 + K ] -cos(c 4 K ⊥,1 + s 4 K ⊥,2 -3K ) -cos[(1 + c 2 )K ⊥,1 -s 2 K ⊥,2 -K ] 5K ⊥,1 + (6s 4 -2s 2 )K ⊥,2 (G3)
for all reciprocal-lattice vectors K with components K , j = K • e , j , K ⊥, j = K • e ⊥, j , and K = K • e .

Structure factor

We are interested in computing the behavior of S when |K | goes to 0 for K = 0. To that end, we note that for a vector K = 2π (s, t, u, v, w) [where (s, t, u, v, w) ∈ Z 5 ], one has c 2 K ⊥,1 + s 2 K ⊥,2 -3K = -c 4 K ,1 + s 4 K ,2 -5K + 2π v, (G4) Keeping in mind that 5K = π (s + t + u + v + w), one can finally rewrite Eq. (G3) as a function of K , j and K ⊥, j only. In the limit where |K | → 0, one then gets generically

c 2 K ⊥,
S(K ) |K |→0 N ( √ 5 -2) 2 K 2 ⊥,1 + K 2 ⊥,2 (K 2 ,1 -K 2 ,2 )K ⊥,2 -2K ,1 K ,2 K ⊥,1 K ⊥,2 5K 4 ⊥,1 -10K 2 ⊥,1 K 2 ⊥,2 + K 4 ⊥,2 2
.

(G14)

However, when one of the denominators in Eq. (G3) vanishes, one gets different expressions that are easily obtained along the same lines.

To analyze the behavior of the structure factor for small |K |, we consider the matrix M = ⎛ ⎜ ⎜ ⎜ ⎝ 0 1 0 0 1 1 0 1 0 0 0 1 0 1 0 0 0 1 0 1 1 0 0 1 0

⎞ ⎟ ⎟ ⎟ ⎠ , ( G15
)
whose eigenspaces are E , E ⊥ , and E with eigenvalues λ = 1/τ , λ ⊥ = -τ , and λ = 2, respectively. By construction, in the large-p limit, the parallel components of K (p) = M p K vanish as λ p , whereas its perpendicular components diverge as λ p ⊥ . As a result, in the large-p limit, S(K (p) ) behaves as λ 8p so that

S(K ) ∼ |K |→0 |K | 8 .
(G16)

Integrated intensity function

As for the octagonal tiling, we decompose the sum over K in the integrated intensity function defined in Eq. (C11) as a sum over all quintuplets (s, t, u, v, w) and their iteration under K (p) = M p K . One can then write

Z (k) = 4π 2 A
(s,t,u,v,w) ∞ p=p (s,t,u,v,w) S(K (p) ), (G17) where p (s,t,u,v,w) is the smallest integer fulfilling the constraint |K (p) | < k. In the small-k (large-p) limit, one can check that S(K (p) ) |K (p) | 8 f (s, t, u, v, w) (G18)

for any quintuplet (s, t, u, v, w) ∈ Z 5 . However, contrary to the octagonal tiling, the scaling of K (p) with p depends on the quintuplet. For quintuplets that do not annihilate the denominator in Eq. (G3), one gets The solution of this equation can be written as

|K (p) | τ -p g(s, t, u, v, w), ( G19
Z (k) = k 8 F (ln k/ ln τ 2 ), (G21)
where F (x + 1) = F (x). As a result, Z has a power-law envelope together with log-periodic oscillations.

FIG. 2 .

 2 FIG. 2. Integrated density function Z of the Rauzy (left), Ammann-Beenker (center), and Penrose (right) tilings (see insets for illustrations). Top: Log-log plot (blue dots) together with the algebraic envelope k 2+α (green line). Bottom: Z (k)/k 2+α as a function of ln k/ ln λ (ln k) showing 1-periodic (nonperiodic) oscillations for tilings with (without) discrete scale invariance.

3 FIG. 3 .

 33 FIG. 3. Integrated density function Z of the Fibonacci chain (see the inset). Top: Log-log plot (blue dots) together with the power-law envelope k 4 (green line). Bottom: Z (k)/k 4 as a function of ln k/ ln τ showing periodic oscillations (with period 1).

  where p (s,t,u,v) is the smallest integer fulfilling the constraint |K (p) | < k. As discussed above, the behavior of S in the largep (small-|K |) limit strongly depends on K [see Eqs. (F10)-(F13)]. This is in stark contrast with the Rauzy tiling and the Fibonacci chain, where there is the same power-law scaling for all K [see Eqs. (D10)-(E4)].

  )or, in other words, |K (p) /K (p+1) | = τ . Importantly, when one of the denominators in Eq. (G3) vanishes, one gets a weaker relation since one only has |K (p) /K (p+2) | = τ 2 . As a direct consequence, one gets the following discrete scaling relation:Z (k/τ 2 ) = Z (k)/τ16 .(G20)

  1s 2 K ⊥,2 -3K = -c 4 K ,1s 4 K ,2 -5K + 2π u, (G5) K ⊥,1 /2 + (s 4 + s 2 )K ⊥,2 + K = -K ,1 /2 + (s 4s 2 )K ,2 + 5K -2π (u + w), (G6) K ⊥,1 /2 -(s 4 + s 2 )K ⊥,2 + K = -K ,1 /2 -(s 4s 2 )K ,2 + 5K -2π (v + t ), (G7) K ⊥,1 -3K = -K ,1 -5K + 2π s, (G8) 2c 4 K ⊥,1 -K = -2c 2 K ,1 -5K + 2π (t + w), (G9) c 4 K ⊥,1 + s 4 K ⊥,2 -3K = -c 2 K ,1s 2 K ,2 -5K + 2πt, (G10) c 4 K ⊥,1s 4 K ⊥,2 -3K = -c 2 K ,1 + s 2 K ,2 -5K + 2π w, (G11) (1 + c 2 )K ⊥,1 + s 2 K ⊥,2 -K = -(1 + c 4 )K ,1 + s 4 K ⊥,2 -5K + 2π (s + v), (G12) (1 + c 2 )K ⊥,1s 2 K ⊥,2 -K = -(1 + c 4 )K ,1s 4 K ⊥,2 -5K + 2π (s + u).

	(G13)

ACKNOWLEDGMENTS

We thank M. Duneau, T. Fernique, and J.-M. Luck for fruitful discussions.

APPENDIX A: VARIANCE FOR A PERIODIC POTENTIAL

The Fourier transform of the periodic potential ( 13) is given by

2π/a δ q y ,0 + δ q x ,-2π/a δ q y ,0 + δ q x ,0 δ q y ,2π/a + δ q x ,0 δ q y ,-2π/a , (A2)

where we used the fact that Aδ q,0 = (2π ) 2 δ(q) in the thermodynamic limit. The structure factor [START_REF] Fuchs | Landau levels in quasicrystals[END_REF] becomes

2π/a δ q y ,0 + δ q x ,-2π/a δ q y ,0 + δ q x ,0 δ q y ,2π/a + δ q x ,0 δ q y ,-2π/a 2 , (A3)