arXiv:1904.08285v2 [math.DS] 24 Apr 2019 { ulof hull 1 n ope ojgt pair conjugate complex a and euse we ievleof eigenvalue nw sthe as known P hc r titypstv,adnraie uhthat such normalised and positive, strictly are which ons h vrg itnebtennihorn onsi points neighbouring between distance average The points. ihcaatrsi polynomial characteristic with osbecodntso u oto ons(eaiet one to (relative points control our of coordinates possible | dens( v ,b c b, a, β , i : otr h yblcsqecsit iig,w hoeinter choose we tilings, into sequences symbolic the turn To h leri integer algebraic The osdrtepiiiesubstitution primitive the Consider ocniu,w okwt h rank-3 the with work we continue, To = r h eaiefeuniso h etr ntesmoi s symbolic the in letters the of frequencies relative the are and IFATO FAMDLSTWT OPE WINDOWS COMPLEX WITH SET MODEL A OF DIFFRACTION Λ | v } 1 = ) α ih ̺ e [ see ; ore rnfr fteRuyfatl,vaarpdyconve rapidly a via fractals, Rauzy the of transform Fourier on ircinmauevaaFuirmti oyl,whic cocycle, matrix Fourier Here a via structure. measure complicated diffraction windo a point whose with sets fractals model Rauzy regular are as smallest spectrum space described point the be pure is can sets has factor point system inflation dynamical whose corresponding line The real the of tiling Abstract. hslast ye on set, point typed a to leads this , β u o h ubrwt oiieiaiaypart. imaginary positive with number the for 2 | ihcnrlpit nterlf npit.Frec seque each For endpoints. left their on points control with , sapoetro ak1 htis, that 1, rank of projector a is / h s 4 ¯ M lsi number plastic u x .]frdtis t usiuinmti reads matrix substitution Its details. for 4.4] Ex. , = | hr h orsodn etadrgtegnetr read eigenvectors right and left corresponding the where , = 23 1 h elkonpatcnme usiuingvsrs oat a to rise gives substitution number plastic well-known The 3+ (3 + 3 β β 23 β β 7 + 7 + = saui,adtesals io–iaaahvn(V numbe (PV) Pisot–Vijayaraghavan smallest the and unit, a is opr [ compare ; IHE AK N W GRIMM UWE AND BAAKE MICHAEL β β α, 9+ (9 2 x 2 (1 n dens( and ) α 3 ,β β, , − fagbaccnuae,where conjugates, algebraic of √ x M 69) − ̺ 2 4 and ) Z = 1 .Telte sirdcbe ihoera root, real one with irreducible, is latter The 1. 3 .5 n x .7.I sas h Perron–Frobenius the also is It 2.17]. Ex. and 50 p. , endby defined 18 -module (9 + Λ    1 3 Λ P 0 1 1 0 0 1 1 0 0 i = 1 = ) 2 − = Λ | √ a v v P L i ∪ 69) i ˙ a dens( = Λ (2 =    n im( and 7→ b 1 3 h Z , ∪ 1 ˙ [ b | ≈ Λ β fte,te osdrda ) The 0). as considered then them, of disacoe oml o the for formula closed a admits h v Λ − gn nnt product. infinite rging 7→ ecluaetersligpure resulting the calculate we , = ] c i io–iaaahvnnumber. Pisot–Vijayaraghavan for ) n si w-iesoa internal two-dimensional in ws rmtetretpso control of types three the from , 1 β = P . c n h soitdcontrol associated the and , Λ 32472 2 h = ) aso aua egh namely length, natural of vals 7→ β , 1 h s¯ is u ,β β, , i qecsdfie by defined equences | 2 α | ab { ∈ s v R − | , i 2 + 4 = 2 | ntetrayalphabet ternary the on ,β β, 2 v ,b c b, a, = .S,teetisof entries the So, 1. = i i Z . β rayinflation ernary hc opie all comprises which , c ntesymbolic the in nce − − } 1 β . 1) = − T β 3 , β 2 − 2 ow get we so , .Below, 1. ̺ and , r, 2 MICHAEL BAAKE AND UWE GRIMM

Minkowski embedding [4, Ex. 3.5] of L is a lattice in R C R3, with dens( ) = 2/√23. L × ≃ L A canonical choice for the basis matrix of and its dual, , is L L∗ 2 1 1 β β Im(α)β− Im(α)β Im(α) 2 2 2 = 1 Re(α) Re(α ) and ∗ = 2Im(α)β Im(α)β Im(α) , B   B √23  − −  0 Im(α) Im(α2) β 1+ 3 β2 3 β    − 2 − 2  2 2  2 2  4+9β 6β 2 6+2β 9β where Re(α) = β/2, Re(α ) = 1 β /2, Im(α) = − and Im(α ) = − in − − 2√23 2√23 our setting. We note in passing that duality can also be defined via the quadratic form (x,y) := xy +Re σ(x)σ(y) , where σ : Q(β) Q(α) denotes the algebraic conjugation map −→ defined by β α and its unique extension to a field isomorphism. 7→
We can extract the Fourier module from the first line of the dual basis matrix ∗ as 2 B ⊛ 2Im(α) 2 5 6β + 4β (1) L = 1, β , β Z = − L, β√23 h i 23 which is also the dynamical spectrum (in additive notation) of the tiling dynamical system as well as the (topologically equivalent) model set dynamical system defined by the control point sets; see [5] for background. Our system has pure point diffraction spectrum, because the defining point set is a regular model set (see below).

The Bragg peaks can be indexed by three integer Miller indices (n0,n1,n2), where we parameterise the wave number, in line with (1), as 5 6β + 4β2 (2) k = k(n ,n ,n ) = − n + n β + n β2 . 1 2 3 23 0 1 2 The original version of the Minkowski embedding ofZ[β] uses complex
numbers, which is natural from an algebraic perspective. However, for Fourier analysis, we better work with 2 real numbers, via C R , as initiated above for and ∗. Given the conjugation map σ, ≃ B B T the ⋆-map: Q(β) R2 is conveniently defined by x Re(σ(x)), Im(σ(x)) . If k L⊛ is −→ 7→ ∈ parameterised as in (2), this means
1 (18n 4n + 6n ) + (6n 9n + 2n )β (4n 6n + 9n )β2 k⋆ = 46 0 − 1 2 0 − 1 2 − 0 − 1 2 . 1 2n + (3n 2n )β 3n β2 2√23 2 1 − 0 − 2
! Next, we need the displacement (or digit) matrix T of our inflation
rule in direct space; see [1] for definitions and properties. With the interval lengths as chosen above, it reads ∅ ∅ 0 { } T = 0 ∅ 1 , { } { } ∅ 0 ∅  { }  which is also reflected in the point set iteration 

Λ′ = βΛ , Λ′ = βΛ (βΛ + 1) , Λ′ = βΛ . a c b a ∪ c c b It is known that we obtain a regular model set [9], with three topologically regular windows that have disjoint interiors. They are compact sets Wa,Wb,Wc C that satisfy ⊂ (3) W = αW , W = αW (αW + 1) , W = αW , a c b a ∪ c c b DIFFRACTION OF A MODEL SET WITH COMPLEX WINDOWS 3

Figure 1. Illustration of the three windows Wa (blue), Wb (red) and Wc (green) of the plastic number inflation rule. For a more detailed version of Wb and further properties, see [4, Fig. 7.3]. which defines a contractive iterated function system (IFS), so that the fixed point (and hence our window triple) is unique. Due to the nature of Λ and Λi as regular model sets, the volumes (areas) of the windows, compare [4, Ex. 7.6 and Fig. 7.3], satisfy the relations dens( ) vol(W )= v dens(Λ), which results in the values L i i Im(α) β2 1, β, 1 = 1 5 6β + 4β2, 6 2β + 9β2, 4 + 9β 6β2 . − 2√23 − − − − The three windows are complex
Rauzy fractals with a complicated topological
structure [8]; see Figure 1 for an illustration and [7, Sec. 7.4] and references therein for background on Rauzy fractals. From the general diffraction result for regular model sets, see [4] and C references therein, the weighted Dirac comb ω = i a,b,c hi δΛi with hi has diffraction ∈{ } ∈ P 2 (4) γω = I(k) δk with I(k) = hi Ai(k) ,

k L⊛ i

X∈ X

c c c because the Bombieri–Taylor property holds for primitive inflation rules [3]. Here, setting W = Wa Wb Wc, the amplitudes are given by ∪ ∪

dens(Λi) ⋆ dens(Λ) ⋆ Ai(k) = 1Wi (k ) = 1Wi (k ), vol(Wi) vol(W ) 4 MICHAEL BAAKE AND UWE GRIMM

I18(k)

k 0) , 1) , 1 , 1 , (1 3) , (1 5) 5 , , 15) 6 , , (3 1) 8) , 7) 20 , (3 , , 2 , 13) , 11 10 , , (1 20) (11 , 17 (6 , (6 18) , 3) 2) 2) 26 14) , , , , , 4 2 3 17) (10 24 11) , , , , , , 18 4) , (15 , (2 (1 (2 19) 23 14 5 , , 14) , (14 , , (10 25 (8 5) (3 , 19 , (13 9) 7) , , 7 I(k) , , 12) 9 , , (14 12 16) , (11 (4 , 0) 16 (5 , , (7 0 21 , , (9 (0 (12

Figure 2. Bragg peaks of the plastic number inflation rule (red, bottom;

some with their Miller index triple) as obtained from the cocycle approach, in comparison with a finite-sizec approximation by exponential sums (blue, top).

where 1K denotes the characteristic function of the set K. Though this formula looks nice,

it is difficult to calculate 1Wi directly, due to the fractal nature of the windows. Let us next explain an alternative approach that harvests the inflation nature of our point sets.

2 We start by reformulating the window IFS (3) in R . Recall thatc multiplication with α C ∈

Re(α) Im(α) R2 2 1 is matrix multiplication with Q = Im(α)− Re(α) in , where det(Q)= α = β− . As the c | | windows are interior-disjoint, Eq. (3), turns into equations in L1(R2) for the characteristic

functions of the windows and their images. Setting f (y) := 1 (y) and using the relation [ i Wi 2πi t y T 1 (y) = det(A) e h | i 1 (A y), AK+t | | K which holds for any real, invertible 2 2-matrix A and compact set K R2, one finds × ⊂

fa fa 0 0 1 1 T 2πiy1 (5) fb (y) = β− B(y) fb (Ry), with R = Q and B(y) = 1 0 e  . f f 0 1 0  c  c         Here, B(y) is the internal Fourier matrix that emerges from the (inverse) Fourier transform 2πi x⋆ y R2 of the ⋆-image of T , that is, Bij(y)= x Tij e h | i, where y . ∈ ∈ N (n) n 1 Now, for n , we can define a cocycleP via B (y) = B(y)B(Ry) B(R − y), with (1) ∈(n) n · · · n (n) B = B and B (0) = M . Next, consider the matrix function C(y) := limn β− B (y), →∞ which is well defined because the limit exists for every y R2. In fact, C is a continuous ∈ DIFFRACTION OF A MODEL SET WITH COMPLEX WINDOWS 5

m=66 log(Im(k)) m=60 m=54 m=48

m=42

(1,2,2)

Figure 3. Example of a Bragg peak and its successive approximation, for wave numbers k [1.2672395, 1.2672405] and finite systems defined by the ∈ point sets from the inflation words ̺m(a) with m 42, 48, 54, 60, 66 . The ∈ { } vertical black line denotes the location of the Bragg peak (1, 2, 2). Note that, to be able to compare different system sizes, the logarithm of the intensity is plotted.

matrix function, and one has 2 β2 β 1 β2 β 3+ β + 7β2 − − − C(0) = P = v u = β2 β 1+ β β2 β2 1 . | ih | 23  − − −  β 1 β2 β 1+ β β2  − − −   dens( Λ) Moreover, one can show that C(y) = c(y) u with c(0) = v and f(y) = dens( ) c(y) . | ih | | i | i | i L | i In particular, one has c(y) = C(y) v , which makes the functions fi accessible. | i | i For k L⊛ and i a, b, c , our amplitudes are ∈ ∈ { } 3+ β + 7β2 A (k) = dens(Λ) c (k⋆) = c (k⋆), i i 23 i and the corresponding intensities follow from Eq. (4). For any index triple, the intensity at the corresponding wave number k L⊛ is approximated by truncating the infinite product ∈ representation for C(k⋆) and calculating the amplitudes as explained above. Here, by using 10 about 50–100 terms, one obtains the peaks with a relative precision better than 10− . Since the diffraction measure γω can also be approximated (in the vague topology) by the absolute squares of exponential sums of finite approximants (divided by the system size), we illustrate the diffraction of the uniformc Dirac comb (all hi = 1) in Figure 2 in comparison to an approximation for a finite system with 114 tiles, obtained by m = 18 inflation steps from an initial tile of type a. Here, as well as in Figure 3, the function Im(k) for the ap- proximation refers to the (absolutely continuous) diffraction intensity for the corresponding finite system. The apparent mismatch emerges from the extremely slow convergence of the finite-size approximation, which should not come as a surprise due to the complicated fractal 6 MICHAEL BAAKE AND UWE GRIMM

2 Figure 4. Inverse Fourier transform fb(y) of the window Wb for y [ 4, 4] . ∈ − The left panel shows the absolute value, which takes values between 0 (blue) 2 and fb(0) = vol(Wb) = v = β β (red). The right panel shows the corre- 2 − sponding argument (with red corresponding to phase 0). nature of the windows. To expand on the latter, we compare the neighbourhood of the peak with Miller indices (1, 2, 2) at k = (1 + 8β + 10β2)/23 1.267240014 with approximations for ≈

various system sizes in Figure 3. While the computation of the spectrum requires that k L⊛, the recursion equationsc ∈ (5) for the (inverse) Fourier transforms of the three windows can be used for any y R2. ∈ Thus, the (inverse) Fourier transform of the windows can be computed efficiently, despite the complex nature of the windows. As an example, the inverse Fourier transform f = 1 of b Wb the largest window is shown in Figure 4. In general, it is difficult to compute the Fourier transform of fractals; see, for instance, [6] for a related result on self-similar fractals of zero Lebesgue measure. The cocycle method explained for this example works in full generality for all primitive PV inflation rules (with a unit inflation factor) that lead to regular model sets. In fact, it can even be applied to any primitive PV inflation rule, as well as to S-adic inflations with the same PV inflation multiplier and compatible tile sizes. When the spectrum is mixed, the cocycle method reproduces the pure point part; further details will be given in a forthcoming publication.

Acknowledgements

It is our pleasure to thank the MFO at Oberwolfach for hospitality, where this manuscript was completed. Our work was supported by the German Research Foundation (DFG), within the CRC 1283 at Bielefeld University, and by EPSRC through grant EP/S010335/1. DIFFRACTION OF A MODEL SET WITH COMPLEX WINDOWS 7

References

[1] Baake M, Frank N P, Grimm U and Robinson E A, Geometric properties of a binary non-Pisot inflation and absence of absolutely continuous diffraction, Studia Math. 247 (2019) 109–154; arXiv:1706.03976. [2] Baake M and G¨ahler F, Pair correlations of aperiodic inflation rules via renormalisation: Some interesting examples, Topology & Appl. 205 (2016) 4–27; arXiv:1511.00885. [3] Baake M, G¨ahler F and Ma˜nibo N, Renormalisation of pair correlation measures for primitive inflation rules and absence of absolutely continuous diffraction, preprint arXiv:1805.09650. [4] Baake M and Grimm U, Aperiodic Order. Vol. 1: A Mathematical Invitation, Cambridge Univer- sity Press, Cambridge (2013). [5] Baake M and Lenz D, Spectral notions of aperiodic order, Discr. Cont. Dynam. Syst. S 10 (2018) 161–190; arXiv:1601.06629. [6] Dettmann C P and Frankel N E, Structure factor of deterministic fractals with rotations, Fractals 1 (1993) 253–261. [7] Pytheas Fogg N, Substitutions in Dynamics, Arithmetics and Combinatorics, LNM 1794, Springer, Berlin (2002). [8] Siegel A and Thuswaldner J M, Topological Properties of Rauzy Fractals, M´em. Soc. Math. France 118, Soci´et´eMath´ematiques de France, Paris (2009). [9] Sing B, Pisot Substitutions and Beyond, PhD thesis, Bielefeld University (2007); available elec- tronically at https://pub.uni-bielefeld.de/publication/2302336.

Fakultat¨ fur¨ Mathematik, Universitat¨ Bielefeld, Postfach 100131, 33501 Bielefeld, Germany E-mail address: [email protected]

School of Mathematics and Statistics, The Open University, Walton Hall, Milton Keynes MK7 6AA, UK E-mail address: [email protected]