Bethe s Ansatz for random tiling models Les Houches, 2008 Bernard Nienhuis 1. Quasi-crystals 2. Random tilings 3. The square triangle tiling 4. Bethe Ansatz equations 5. Solution to the BAE 6. Field theoretic aspects 7. Rectangle triangle tilings with other symmetries 8. Scetch of nested B.A. for octagonal case.
Quasi-Crystallography The crystalline phase is invariant for discrete translations and rotations. Rotational symmetries are exclusively two, three, four or six-fold. Why? R 2 a R a a R 3 a R 4 a Let f be some density function. Translations: f( x) = f( x + a) Rotations: f( x) = f(r x) Then also: f( x) = f ( x + R a + R 4 a ) But a and R a + R 4 a have same direction, but irrational length ratio. If f is continuous it is constant. 1
3D is only slightly more complicated: a need not be orthogonal to rotation axis. b = a + R a + R 2 a + R 3 a + R 4 a points along rotation axis. Compare now: 5(R a + R 4 a) 2 b and 5 a b. Yet many materials exist with icosahedral symmetry. (which has 5-fold axes). Mostly binary and ternary (Al) alloys. There are also solids which are periodic in one direction and have five-, eight-, ten- or twelve-fold symmetry in the orthogonal plane. Rather than periodic these materials could be quasi-periodic. 2
Periodic functions have discrete Fourier spectrum. But the reverse is not true: take sin(x) + sin(πx) 2 1 0-1 -2 0 10 20 30 40 50 This function is almost periodic, but not periodic. The function sin(x)+sin(πx) is obviously the restriction of sin(x)+sin(y) with the line y = πx. In general, quasi-periodic functions are the restriction of a higher dimensional periodic function. 3
Let W ( r) = A m exp i m 1,m 2,,m D D j=1 with D > d (dimension of space). m j a j r Necessarily { a j j = 1,, D} are linearly dependent. But we choose them independent under the integers. Note however, that { D j=1 m j a j m Z D } lie dense in R d. View the vectors a j as d-dimensional restrictions of D-dimensional vectors, say: p j = a j + b j, with r b k Then, W ( r) is simply a restriction of a D- dimensional periodic function. The basis vectors a j can have any symmetry, crystallographic and otherwise. 4
If the D-dimensional function is zero except on periodically repeated objects of (D d) dimensions, the restriction to real space consists of δ-peaks. (Then both the function and its Fourier Transform consist of δ-functions.) In this (d = 1, D = 2) example the realspace function has δ-peaks between (irregularly) alternating long and short intervals. If D > 3 it is possible to find irrational subspaces, which have special symmetries. 5
In this example objects are placed on the sites of a 4-dimensional hypercubic lattice, and intersected with a plane. 6
But what holds these materials in check? Is the quasi-periodic arrangement the minimum of some interparticle potential? In periodic crystals each atom has the same environment. In a quasi-periodic arrangements, there are very many different local environments. Are they all potential minima? 7
It is at least as plausible, that these many local environments are almost equally likely, and that the rotational symmetry is the result of thermal averaging. Effectively, in this scenario, the minimization of potential energy is replaced by minimization of free energy. Asymmetric local configuration occur in all symmetric equivalent orientations. 8
Take this to the extreme, and make the potentials strictly short-range, and many local configurations strictly of equal energy. In this case minimization of free energy is effectively replaced by maximization of entropy. This is for instance achieved by binary hard disk packings: This naturally leads to 9
Random tilings are complete covering of space by copies of one or more rigid geometric objects. A well-known example is the dimer model: space is covered by dominos (1 2 rectangles, or 1 1 2 bricks) Is is much more interesting if the building blocks permit (or even force) non-crystallographic structures. 10
For general binary ball packings: if the potential is strictly short range, e.g. hard disks, individual configurations need not be symmetric or quasiperiodic. But the ensemble of all tight packings shows a discrete rotational symmetry, and may come out to be quasiperiodic. The tight binary disk packing is equivalent to a so-called Random Tiling with (a subset of) three kinds of triangles: The angles of the isosceles triangles may be such that they support a non-crystallographic rotational symmetry. Both the rotational symmetry and the quasiperiodicity may be subject to symmetry breaking Other well-studied tilings in 2D are rhombus tilings and in 3D rhombohedron tilings. 11
The straight intersection required by strict quasi-periodicity: is in a random tiling replaced by an intersection with a fluctuating line: A natural description would be a free field theory, compactified to the D-dimensional lattice. 12
a D = 2 + 2 example: Square-Triangle tiling The specific problem: tiling of the plane by squares and equilateral triangles. A configuration of the square triangle random tiling 13
14
strong dependence on relative densities: predominantly squares 15
even mixture 16
predominantly triangles Relation between total numbers of tiles in the different orientations: S 1 S 2 + S 2 S 3 + S 3 S 1 = 1 4 T 1T 2 expressed in areal densities: σ 1 σ 2 + σ 2 σ 3 + σ 3 σ 1 = 4 3 τ 1τ 2 17
Towards a solution: Discrete statistical problem: Yet no discrete underlying lattice tile edges can be in one of 6 orientations: (say angle with horizontal is π/12, 3π/12, 5π/12 etc.) Deform the tiles by rotating each edge to nearest multiple of π/3 (dotted direction). 18
The above configuration is then deformed into the following: edges need to be marked to distinguish original orientations 19
Now we have a lattice model: We mark edges with colored lines. Five configurations of each triangle: horizontal edges: diagonal edges: Continuity of the colored lines guarantees the tiles to fit without holes or overlaps. The number of red and green lines cutting each horizontal is conserved. 20
transfer matrix: Z N (α, β): Partition sum of bounded cyclinder, T α,β = Z 1 (α, β): weight of single layer of tiles α, β, γ : state of row of horizontal edges Z N (α, β) = γ Z N 1 (α, γ)t γ,β β γ α toroidal partition sum: Z N = α Z N (α, α) = Trace T N 21
Geometrical considerations: In the transfer matrix the numbers of red and green marked horizontal edges, say n r, and n g are conserved. How are these numbers related to the tile densities? We had for the tile numbers: S 1 S 2 + S 2 S 3 + S 3 S 1 = 1 4 T 1T 2 and their areal densities: σ 1 σ 2 + σ 2 σ 3 + σ 3 σ 1 = 4 3 τ 1τ 2 In a symmetric phase, σ 1 = σ 2 = σ 3 = σ 3 τ 1 = τ 2 = τ 2 Therefore σ = τ. In the fully symmetric phase half the area is covered by triangles and half by squares. 22
Most symmetric phase dominated by squares is e.g. σ 1 > σ 2 = σ 3 and τ 1 = τ 2. On average rectangular areas l (l + 1): S 1 = l(l + 1) S 2 = S 3 = 1 2 T 1 = T 2 = 2(2l + 1) τ = 3(2l + 1) 2l 2 + 2l + 2 + (2l + 1) 3 n g L = n r L = 1 l + 2 This density regime has square symmetry 23
Most symmetric phase with predominantly triangles: σ 1 = σ 2 = σ 3 and τ 1 > τ 2. On average hexagonal areas with sides l: T 1 = 6l 2 T 2 = 2 S 1 = S 2 = S 3 = l τ = 3l2 + 1 (1 + l 3) 2 n r L = l + 1 3l 2 + 3l + 1 n g L = 3l2 + l 3l 2 + 3l + 1 This density regime has hexagonal symmetry 24
Diagonalization of transfer matrix by Bethe s Ansatz (M. Widom) reference state (pseudo vacuum): eigenvalue: 1 single red: plane wave solution ψ(x) = e ipx u x Λ ψ(x) = (T ψ)(x) = t ψ(x + 1) Λ = tu t: weight of two triangles 25
single green: ψ(y) = v y Λ ψ(y) = (T ψ)(y) = t ψ(y 1) Λ = t/v one red and one green (isolated): ψ(x, y) = u x v y Λψ(x, y) = t 2 ψ(x + 1, y 1) Λ = t 2 u/v 26
red followed by green: Involves both x < y and y < x. Amplitude may depend on the order: ψ(x, y) = A r,g u x v y with x < y ψ(y, x) = A g,r u x v y with x > y (arguments ordered) The equation from above diagram: Λψ(x 1, x + 1) = tψ(x 2, x) + tψ(x, x + 2) with Λ = t 2 u/v Substitution of the Ansatz results in: t A r,g = A g,r ( u 2 + v 2) 27
Now consider three particles (say two red and one green), and try a product of three plane waves (wave numbers u 1, u 2 and v ). Equations when particles are isolated yield: Λ = t 3 u 1 u 2 v 1 Collisions (interchanges) between red and green particles force the amplitude to depend on the red-green sequence. We allow the red particles to interchange momentum: ψ(x 1, x 2, y) = A 1,2,v u x 1 1 ux 2 2 vy + A 2,1,v u x 1 2 ux 2 1 vy ψ(x 1, y, x 2 ) = A 1,v,2 u x 1 1 vy u x 2 2 + A 2,v,1 u x 1 2 vy u x 2 1 ψ(y, x 1, x 2 ) = A v,1,2 v y u x 1 1 ux 2 2 + A v,2,1 v y u x 1 2 ux 2 1 (arguments ordered) 28
When the green particle collides with one of the red, we find: together with A 1,2,v = t 1 A 1,v,2 ( u2 2 + v 2 ) etc. But when three particles collide at once: Eigenvalue equation: Λψ(x 2, x, x + 2) = t ψ(x 3, x 1, x + 1) + t 2 ψ(x 1, x + 1, x + 3) This results in and A 1,2,v = A 2,1,v A v,1,2 = A v,2,1 29
The natural generalization: a linear combination of products of plane waves ψ ( ) c x = σ A c σ j ( w c j σ j ) xj c : the color sequence, x : position sequence (x j < x j+1 ), (both of the non-blank edges). wj r = u j and w g k = v k are the wave numbers of the red and green particles respectively. σ : permutes these wave numbers, but does not mix red and green ones. The eigenvalue must be Λ = j t u j k t v 1 k and the amplitudes A satisfy: A,r,g,,j,k, = t 1 ( u 2 j + v 2 k A,r,r,,j,k, = A,r,r,,k,j, A,g,g,,j,k, = A,g,g,,k,j, ) A,g,r,,k,j, 30
Assume periodic boundary conditions: ) ( ) c, = ψ ψ (,c,2l this is satisfied if: 0, A,c (,k w c) 2L c, k = A k, L: horizontal system size in lattice edges. The left and right hand side of this equation can also be related by successive interchanges of k with its predecessor. c the Bethe Ansatz equations (BAE) u 2L j = ( ) n n g r+1 t n g v 2L k Λ = = ( ) n g+1 t n r n r j=1 t u j k=1 n r j=1 n g k=1 u 2 j + v 2 k u 2 j + v 2 k t v 1 k Convenient to introduce ξ = u 2 and η = v 2. 31
In these new variables the BAE read ξ L j = ( ) n r+1 t n g n g k=1 ξ j η k η L k = ( )L+n g+n r +1 t n r And the eigenvalue: Λ = n r j=1 t ξ 1/2 j Introduce the functions: F r (z) = log z F g (z) = log z n g k=1 n r j=1 Then the BAE read: n g k=1 n r j=1 η k ξ j t ( η k ) 1/2 1 L log(z η k) + n g L log t 1 L log(z ξ j) + n r L log t LF r (ξ j ) = (0 or πi) mod 2πi LF g (η j ) = (0 or πi) mod 2πi 32
Expect solutions along curves in the complex plane, given by 1.5 1 Re F r (ξ j ) = 0 Re F g (η k ) = 0 0.5 η ξ 0-0.5-1 -1.5-1.5-1 -0.5 0 0.5 1 1.5 33
At successive roots: F r (ξ j+1 ) F r (ξ j ) = 2πi L Therefore in the thermodynamic limit f r (z) F r(z) z = 2πi L (root density) With this knowledge we can transform the sums into integrals: f r (z) = 1 z 1 L n g k=1 1 z η k = 1 z 1 2πi η1 η n g f g (η) η z dη (double sign change: z η and η 1 η ng ) likewise f g (z) = 1 z + 1 2πi ξn r ξ 1 f r (ξ) ξ z dξ So, for the analyticity of f r and f g : 34
f r analytic except at η, f g analytic except at ξ. η ξ What is the change of these functions at the cuts? f r (η + ε) f r (η ε) f r (z) = 1 z 1 2πi η1 η n g f r (z + ε) f r (z ε) = 1 2πi = f g (z) f g (η) η z dη η z =ε f g (η) η z dη 35
Solution of the integral equations (P. Kalugin) So f r jumps by f g across the (η) cut. analytic continuation through the cut: f r f r f g likewise f g taken through the ξ-cut changes: f g f r + f g Now inspect the function f(z) = a r f r (z) + a g f g (z). Through the η-cut f g is analytic, and f r changes by f g. This effectively multiplies the vector ( ) 1 0 by the matrix Γ g =. 1 1 ( ar a g ) At the ξ-cut the corresponding ( ) monodromy 1 1 transformation is Γ r =. 0 1 36
Together Γ g and Γ r generate the entire group SL 2 ( Z). (2 2 integer matrices with determinant 1). However, when the triangles cover half the area (the point of 12-fold symmetry): 1.5 1 0.5 0 η ξ -0.5-1 -1.5-1.5-1 -0.5 0 0.5 1 1.5 37
And also when τ > σ (hexagonal symmetry) 1.5 1 0.5 0 η ξ -0.5-1 -1.5-1.5-1 -0.5 0 0.5 1 1.5 The cuts appear to have a common end point in the thermodynamic limit. 38
When that is the case there is no way to run through one cut and return to the point of origin without going through the other cut. Effectively there is only one monodromy operator: ( ) 1 1 Γ g Γ r = 1 0 Note that: therefore and (Γ g Γ r ) 2 = (Γ g Γ r ) 3 = ( 0 1 1 1 ( 1 0 0 1 (Γ g Γ r ) 6 = 1. ) ) But can we believe that the cuts close on the basis of numerical evidence? 39
We have three parameters to play with: n r, n g, and t. The restriction that the cuts close is one complex equation; So it is reasonable to assume that it can be satisfied. Let, therefore, in the thermodynamic limit η 1 b ξ nr b η ng b ξ 1 b Then, because (Γ g Γ r ) 6 = 1, f r (z) and f g (z) must be single-valued functions of w Let w 6 = z/b 1 1 z/b z = w6 + 1 w 6 /b + 1/b f(w)dw = f r (z)dz In this way, poles of f r (z) translate into poles of f(w) with the same residues. 40
The BAE: f r (z) = 1 z 1 L n g k=1 1 z η k f g (z) = 1 z 1 L n r j=1 1 z ξ j Therefore f r (z)dz and f g (z)dz have poles in the origin, (residu 1) and at (residu n g /L 1 and n r /L 1 respectively). What does that imply for f(w)dw? First define: b = i b e iγ ω = e πi/6 θ = e γi/3 41
η ξ z 6 7 5 3 4 0 0 0 8 2 w 0 1 12 0 11 0 10 9 Then the poles of f(w)dw are: m z w m f(w)dw/dz residue r m 1 0 ω f r (z) 1 2 ω 2 θ f r (z) n g /L 1 3 0 i f r (z) f g (z) 0 4 ω 4 θ f g (z) 1 n r /L 5 0 ω 1 f g (z) 1 6 θ f r (z) f g (z) 2 (n r + n g )/L 7 0 ω f r (z) 1 8 ω 2 θ f r (z) 1 n g /L 9 0 i f r (z) + f g (z) 0 10 ω 4 θ f g (z) n r /L 1 11 0 ω 1 f g (z) 1 12 θ f r (z) + f g (z) (n r + n g )/L 2 42
f(w) is determined completely by its poles: f(w)dw = 12 m=1 r m dw w w m. This is the solution of f(w) in terms of b, n r /L and n g /L. The locus of the roots was determined by Therefore also Re F r (ξ j ) = 0 Re F g (η k ) = 0 0 = Re df r (z) = Re f r (z) dz 0 = Re df g (z) = Re f g (z) dz Since in z = b two cuts join from different directions, f r (z = b) = f g (z = b) = 0. Likewise f(w = 0) = 0, i.e. 12 m=1 r m w m = 0 43
When the values of r m and w m are substituted, this (complex) equation reads (1 n g /L)ω 1 θ 1 + (1 n r /L)ωθ 1 = 1 The solution: n r /L = 1 2 3 cos π γ 3 n g /L = 1 2 3 cos π + γ 3 is consistent with hexagonal symmetry. (from the geometrical considerations) with γ related to the triangle density: sin γ 3 = 2τ 1 (1 τ) 3 + 2τ i.e. τ = 2 3 + 4 cos(γ/3)(2 3 + 3) cos(γ/3) + 1 For γ = 0, τ = σ = 1/2, consistent with the maximal, i.e. twelve-fold, rotational symmetry. 44
We now know f r (z), f g (z), but we still have not calculated the eigenvalue. Λ = n r j=1 t ξ 1/2 j n g k=1 t ( η k ) 1/2 So we need to know t, ξ and η. Remember: F r (z) = log z n g k=1 1 L log z η k t F g (z) = log z n r j=1 1 L log z ξ j t we know that Re F r (b) = Re F g (b) = 0. Re F r (z 0) is much like k log(t/η k )/L, except for the logarithmic singularity. Likewise in the limit z, F r has a finite contribution n g /L log t, modulo the logarithmic divergence. 45
Therefore we define F 1 (z) = F r (z) log z F 2 (z) = F r (z) (1 n g L ) log z Note: (1) Since b is on the locus of the roots, F r (z = b) is imaginary. (2) The limits F 1 (z 0) and F 2 (z ) are regularized, so that Re (F 1 (0) F 1 (b)) = log b + 1 L n g k=1 log t η k Re (F 2 ( ) F 2 (b)) = log b + n g L log t b These can also be calculated directly from integrals of the known functions f r and f g (most conveniently in the w-plane). 0 ( F 1 (0) F 1 (b) = dz f r (z) 1 ) b z w1 ( = dw f(w) 1 ) dz 0 z dw 46
This is a simple integral, since and likewise 1 z f(w) = 12 m=1 12 dz dw = m=1 r m w w m ( 1) m+1 w w m F 2 is subjected to the same treatment: F 2 ( ) F 2 (b) = = w2 0 dw ( b dz ( f(w) (1 n g) z f r (z) 1 n g z dz dw ) ) Analogously F 3 and F 4 are defined from F g : F 3 (z) = F g (z) log z F 4 (z) = F g (z) (1 n r L ) log z 47
Altogether we have four linear equations in the four (yet) unknowns: log(t), log b, j log ξ j and k log η k. From these the eigenvalue Λ can be found: Λ = n r j=1 t ξ j 1/2 n g k=1 t η k 1/2 which is the partition sum per row. Let T row be the number of triangles per row: T row L = 8 4 3 cos γ 3 (from geometric relations with n g and n r.) Then the entropy per vertex is equal to S v = 1 L log Λ T row 2L log t resulting in: 48
+2 cos +2 cos ( π ( π 6 + γ 3 6 γ 3 S v = log 108 cos 2 γ + ) ) log log [ ( π tan 12 + γ 6 [ ( π tan 12 γ 6 ) ) tan tan ( π ( π 4 + γ 6 4 γ 6 )] )] 0.14 0.12 S v 0.1 0.08 0.06 0.04 0.02 0 0 0.2 0.4 0.6 0.8 1 τ For τ < 1 2 solution only numeric. 49
The finite size correction of the entropy gives the central charge: S a = S + cπ 6L 2 When we calculate this as function of the triangle density τ: 2 c 1.5 1 0.5 0 0 0.2 0.4 0.6 0.8 1 τ Does c really depend continuously on τ? Can it be a conformal theory? 50
The spectrum of the transfer matrix scales as a scale-invariant theory: log Λ j = 2π j Λ 0 L all correlation lengths are thus proportional to the circumference of the cylinder. In a conformal theory there are many gaps (Lambda j Λ k ) that scale with an integer amplitude jk. How does that come out in this system? Gaps jk that are integer for τ >.5 turn out complex for τ <.5, but continue to satisfy j Z What does that mean? 51
The resolution of this question turns out to be the following: We propose that the two families of particles even though they interact with each other, effectively behave as free fermions. The only effect of the interaction is that they stretch each others space in an anisotropic way. Though this effect is anisotropic, the combination of the effects of both families has again square symmetry. Each of the families has central charge 1, but the contribution to the effective finite size correction is modified by a geometric factor as a result of the deformation. 52
Solvable random tilings: β a a b = 2 sin β 2 b Tilings of rectangles and isosceles triangles have been solved, i.e. the free energy has been calculated in the thermodynamic limit. Special values of β support special symmetries, in particular if β = 2π n : 53
β = 2π 4 : octogonal de Gier & BN: Phys.Rev.Lett. 76 (1996), 2918 J.Stat.Phys. 87 (1997), 415 β = 2π 5 : decagonal de Gier & BN: J.Phys.A. 31, (1998) 2141 β = 2π 6 : duodecagonal Kalugin: J.Phys.A 27 (1994) 3599 de Gier & BN: Phys.Rev.E 55 (1997), 3926 54
Sketch of the solution of the octogonal rectangle-triangle tiling. The first step is to deform the tiles to fit them to a square lattice The short edges are rotated to the nearest axial direction and the long edges to the nearest diagonal direction. The rotated edges are colored to code for the original orientation. The resulting lattice model can be attacked by standard techniques such as the transfer matrix. 55
The perfect tiling condition (i.e. tiles fit together without holes or overlap) translates into continuity of the red and blue lines. Every horizontal cuts an equal number of red (blue) lines. Conserved quantities under the action of a transfer matrix. 56
The eigenstates of the transfer matrix satisfy a nested Bethe Ansatz. i.e. piecewise plane wave solution, with amplitudes changing at collisions. As a result the eigenvalues Λ follow from Bethe Ansatz equations. ( u ± j ) L m k=1 ( wk u ± j ) = 1 n + ) n i=1 ( wk u + i j=1 ( wk u j ) = 1 Λ = m n w + k k=1 i=1 1 u + i n j=1 1 u j 57
The Bethe Ansatz equations ( u ± j ) L m k=1 ( wk u ± j ) = 1 n + ( wk u + i ) n ( wk u j ) = 1 i=1 j=1 can be solved as follows: introduce the functions F (z) L log (z) + k log (w k z) G(z) i log ( z u + ) i j log ( z u j ) Then each u ± j is a solution of F ( u ± ) j = 0 mod 2πi and each w k of G (w k ) = 0 mod 2πi Hypothesis: The equations for consecutive roots differ by 2πi. Then the functions f(z) = df (z)/dz and g(z) = dg(z)/dz measure the density of roots. 58
In the thermodynamic limit the sums can be turned into integrals, so that the Bethe Ansatz equations are now integral equations with integration contours running over the locus of the roots: f(z) = L z m 2πi g(w) z w dw g(z) = n πi f(w) z w dw where we have chosen n + = n = n. The loci of the roots are thus branch cuts in the functions f(z) and g(z) with a linear monodromy group acting on the analytic continuation af(z) + b g(z). 59
The entropy per unit area is given by: S a = 2 + 2 4 cos 2 γ/4 + cos + cos ( π ( π 4 + γ 2 4 γ 2 ) ) [ log tan log tan log 4 cos γ ( π ( π 8 + γ 4 8 γ 4 An analytic solution is available only for triangle density τ > 1/2: ) )] 0.14 0.12 S a 0.10 0.08 0.06 0.04 0.02 0 with the triangle density 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 τ = ( 2 + 1 ) 2 cos γ/2 1 + cos γ/2 τ 60
Leaves many questions: 1. Correlation functions 2. diffraction patterns 3. Nature of the phases: critical? 4. Scale- vs Conformal invariance. 5. What makes it solvable? 6. Other solvable tilings. in particular: rect.tri. 2π n 7. are the cases n {4, 5, 6} exceptional like E 6,7,8 8. What properties are universal? (i.e. dependent on symmetry alone) 9. convexity of the entropy: related to solvability? 10. 3D random binary/ternary packings 61
references For the Bethe Ansatz: R.J.Baxter Exactly solved models in Statistical Mechanics, Chapter 9, Academic Press, London (1982) M. Gaudin La fonction d onde de Bethe, Masson (1983) For some background on Quasicrystals see e.g. Quasi-crystals the state of the art, ed. P.J. Steinhardt and D.P. Di Vincenzo. (World Scientific, 1991) Chapter 15, by C.L. Henley is devoted to random tiling models. More specific on a one of the models discussed: Random square-triangle tilings, a model of 12-fold quasicrystals, M. Oxborrow, C.L. Henley, Phys.Rev. B 48, 6966 (1993) 62
Recent results discussed in the lectures: Bethe-Ansatz solution of the square-triangle random tiling model, M. Widom, Phys. Rev. Lett. 70 2094 (1993) The square-triangle random tiling model in the thermodynamic limit, P.A. Kalugin, J.Phys. A. 27, 3599 (1994) Exact solution of an Octogonal random tiling model, J. C. De Gier and B. Nienhuis, Phys. Rev. Lett. 76, 2918 (1996) The exact solution of an Octagonal rectangletriangle random tiling, J. C. De Gier and B. Nienhuis, J. Stat. Phys. 87, 415 (1997) Most of the material of these lectures is covered in: B. Nienhuis, Exact solution of random tiling models, Phys. Rep. 301, p271-292 (1998). ith small corrections published in Quantum Field theory: Perspective and prospective, Nato Science Series C530 ( les Houches 1998) eds. C. DeWit-Morette and J-B. Zuber, Kluwers 1999. 63