arxiv:cond-mat/ v3 [cond-mat.stat-mech] 18 Feb 2002

Transcription

1 Close-packed dimers on nonorientable surfaces arxiv:cond-mat/00035v3 cond-mat.stat-mec 8 Feb 2002 Wentao T. Lu and F. Y. Wu Department of Pysics orteastern University, Boston, Massacusetts 025 Abstract Te problem of enumerating dimers on an M net embedded on non-orientable surfaces is considered. We solve bot te Möbius strip and Klein bottle problems for all M and wit te aid of imaginary dimer weigts. Te use of imaginary weigts simplifies te analysis, and as a result we obtain new compact solutions in te form of double products. Te compact expressions also permit us to establis a general reciprocity teorem. keywords: close-packed dimers, non-orientable surfaces, reciprocity teorem Introduction A seminal development in modern lattice statistics is te solution of enumerating close-packed dimers, or perfect matcings, on a finite M simple-quartic net obtained by Kasteleyn and by Temperley and Fiser 2, 3 more tan 40 years ago. In teir solutions te simple-quartic net is assumed to possess free or periodic boundary conditions. In view of te connection wit te conformal field teory 4, were te boundary conditions play a crucial role, tere as been considerable renewed interest to consider lattice models on non-orientable surfaces 5, 6, 7, 8, 9. For close-packed dimers te present autors 6 ave obtained te generating function for te Möbius strip and te Klein bottle for even M and. Independently, Tesler 0 as solved te problem of perfect matcings on Möbius strips and deduced solutions in terms of a q-analogue of te Fibonacci numbers for all {M, }. It turns out tat te explicit expression of te generating function depends crucially on weter M and being even or odd, and te analysis differs considerably wen eiter M or is odd. In tis paper we consider te general {M, } problems for bot te Möbius strip and te Klein bottle by introducing imaginary dimer weigts. Te use of imaginary weigts simplifies te analysis and as a result we obtain compact expressions of te solutions witout te recourse of Fibonacci numbers. Te compact expressions of te solutions also permit us to establis a reciprocity teorem on te enumeration of dimers. 2 Summary of results For te convenience of references, we first summarize our main results. Details of derivation will be presented in subsequent sections. Consider an M simple-quartic net consisting M sites arranged in an array of M rows and columns. Te net forms a Möbius strip if tere is a twisted boundary condition in te orizontal direction as sown in Fig., and a Klein bottle if, in addition to te twisted boundary condition, tere is also a periodic boundary condition in te vertical direction. Let te dimer weigts be z and z v, respectively, in te orizontal and vertical directions. We are interested in te close-form evaluation of te dimer generating function Z M, z,z v = z n znv v

2 were n,n v are, respectively, te number of orizontal and vertical dimers, and te summation is taken over all close-packed dimer coverings of te net. Our results are as follows: For bot M and even, we ave 6 Z Mob M, z,z v = Z Kln M, z,z v = M/2 /2 m= n= M/2 /2 m= n= 4z 2 4n π sin zv 2 mπ cos2, 2 M + 4z 2 4n π sin2 + 4z 2 2m π v sin2 2 M were te superscripts refer to te type of te nonorientable surface under consideration. For M even and odd, we ave Z Mob M,z,z v = Re Z Kln M,z,z v = Re and for M odd and even, we ave Z Mob M, z,z v = z /2 Z Kln M, z,z v = z /2 M/2 i 2i M 2 +m+ z sin m= n= M/2 i 2i M 2 +m+ z sin m= n= /2 n= /2 n= M+/2 m= M+/2 m= For M and bot odd, te generating function is zero. 3 Te Möbius strip 4z 2 4n π sin2 2 4n π 2 4n π 2 + 2z v cos + 2z v sin, 3 mπ, 4 M + 2m π, M 5 + 4zv 2 mπ cos2, 6 M + 4z 2 4n π sin2 + 4z 2 2m π v sin2 2 M. 7 It is well-known tat tere is a one-one correspondence between dimer coverings and terms in a Pfaffian defined by te dimer weigts. However, since terms in te Pfaffian generally do not possess te same sign, te evaluation of te Pfaffian does not necessarily produce te desired dimer generating function. Te crux of matter is to attac signs, or more generally factors, to te dimer weigts so tat all terms in te Pfaffian ave te same sign, and te task is reduced to tat of evaluating a Pfaffian. Tese tasks were first acieved by Kasteleyn for te simple-quartic dimer lattice wit free and periodic boundary conditions, wo sowed ow to attac signs to dimer weigts and ow to evaluate te Pfaffian. Soon after te publication of, T. T. Wu pointed out tat te structure of te Pfaffian, and ence its evaluation, can be simplified if a factor i is associated to dimer weigts in one spatial direction. Indeed, te Wu prescription requires only uniformly directed lattice edges wit te association of a factor i to dimer weigts in te direction in wic te number of lattice sites is odd. If te number of lattice sites is even in bot directions, ten te factor i can be associated to dimers in eiter direction. For {M, } = {even,odd} for example, one replaces z by iz. 2

3 To see tat te Wu procedure is correct, one considers a standard dimer covering C 0 in wic te lattice is covered only by parallel dimers wit real weigts. Ten, te two terms in te Pfaffian corresponding to C 0 and any oter dimer covering C ave te same sign, since te superposition polygon produced by C 0 and C contains an even number of arrows pointing in one direction as well as a factor i 4n+2 =, were n is a nonnegative integer. amely, te superposition polygon is clockwise-odd. Tis use of imaginary dimer weigts is te starting point of our analysis. 3. M = even Möbius strip For M even and = even or odd, we write for definiteness M = 2M, = were M, are positive integers. We orient lattice edges as sown in Fig. a. For te time being consider more generally tat te orizontal dimers connecting te first and -t columns ave weigts z and te associated generating function Z Mob M, z,z v ;z = 2M m=0 z m T m 8 were T m T m z,z v is a multinomial in z and z v wit strictly positive coefficients. Te desired result is ten obtained by setting z = z. ote tat T 0 is precisely te generating function wit free boundary conditions. A D A E B C B D C B C C D 2 A D E B A a b Figure : Edge orientations of Möbius strips wit twisted boundary conditions in te orizontal direction. A, B, C, D, E denote repeated sites. a. M = 4, = 5. b. M = 5, = 4. Attac a factor i to all orizontal dimer weigts z. Tis leads us to consider te antisymmetric matrix Az = iz F F T I 2M + z v I F 2M F T 2M + zk + K T J 2M 9 were I is te identity matrix, F T is te transpose of F, and F 2M, K and J 2M are matrices of te order given by te subscripts, F = , K = , J 2M = ow, te Pfaffian of Az gives te correct generating function T 0 in te case of z = 0. For general z we ave te following result: 3

4 Teorem: Te dimer generating function for te simple-quartic net wit a twisted boundary condition in te orizontal direction is ZM,z Mob,z v ;z = ipfaiz + + ipfa iz 2 = Re ipfaiz. Proof: It is clear tat te term in 8 corresponding to te configuration C 0 m = 0 as te correct sign. For any oter dimer configuration C, te superposition of C 0 and C forms superposition polygons containing z edges. We ave te following facts wic can be readily verified: i Te sign of a superposition polygon remains uncanged under deformation of its border wic leaves n z, te number of z edges it contains, invariant. ii Deformations of te border of a superposition polygon can cange n z only by multiples of 2, and te sign of te superposition polygon reverses wenever n z canges by 2. iii Superposition polygons aving 0 or z edges ave te sign +. As a result, we obtain PfAz Az were denotes te determinant of and = X 0 + zx z 2 X 2 z 3 X 3 2 X α = T α + z 4 T α+4 + z 8 T α+8 +, α = 0,,2,3. 3 Te teorem is now a consequence of te fact Z Mob M, z,z v ;z = X 0 + z X + z 2 X 2 + z 3 X 3. Remarks: Te teorem olds also for te Klein bottle wic, in addition to a twisted boundary condition in te orizontal direction, as a periodic boundary condition in te vertical direction see below. For te Möbius strip we ave X = X 3 = 0 wen = even. It now remains to evaluate PfA±iz. To evaluate PfA±iz = A±iz we make use of te fact tat, since F 2M F T 2M commutes wit J 2M, te 2M 2M matrix Az can be diagonalized in te 2M-subspace 6. Introducing te 2M 2M matrix U wose elements are we find U m,m = i m 2 mm 2M + sin π 2M + U 2 mm m,m = i m 2M + sin π, m,m =,2,...,2M, 4 2M + U F 2M F T 2MU m,m = 2icos φ m δ m,m U J 2M U m,m = i M+m δ m,m, m,m =,2,...,2M, 5 were φ m = mπ/2m +. Tus, we can replace te 2M 2M matrices in 9 by teir respective eigenvalues, and express Az as a product of te replaced determinants, namely, Az = i 2M 2M 4 m= A m z 6

5 were we ave taken out a common factor i from eac element of te matrix A m z = 2z v cos φ m I + z F F T + M+m zk + K T. 7 Te matrix A m z can be evaluated for general z in terms of a q-analogue of Fibonacci numbers, but for our purposes wen z = ±iz, te matrix can be diagonalized directly. Define te matrix T = F + i M+m. K =......, i M+m we can rewrite 7 wen z = iz as A m iz = 2z v cos φ m I + z T T. 9 ow T and T commute so tey can be diagonalized simultaneously wit respective eigenvalues e iθn and e iθn, were θ n = M+m+ 4n π/2. 20 Tus, we obtain A m iz = 2z v cos n= and as a result M Aiz = 2z v cos m= n= mπ 2M + + 2i M+m+ z sin mπ 2M + + 2i M+m+ z sin 4n π, n π 22 2 were we ave made use of te fact tat cos φ 2M+ m = cos φ m, 2M+ m = m, and i 2M = M. We tus obtain after taking te square root of 22 M PfAiz = 2z v cos m= n= mπ 2M + + 2i M+m+ z sin 4n π Te substitution of 23 into now yields 4. For = even te Pfaffian 23 is real and 4 reduces to 2. Tere is no suc simplification for = odd. 3.2 M = odd Möbius strip For M = odd and = even or odd, we write, for definiteness, M = 2M +, =. Since te number of rows M is odd, we now attac a factor i to dimers in te vertical direction. Again, we assign weigts z to orizontal dimers connecting te first and -t columns and consider te generating function 2M+ ZM, Mob z,z v ;z = z m T m 24 A similar expression of θ n given in 7 of Ref. 6 contains a typo were M + m+ in te exponent sould read M + m. Tis does not alter te results of Ref m=0

6 defined similar to 8. It is readily verified tat, wit lattice edge orientations sown in Fig. b, all terms in T 0 ave te same sign. It follows tat we can use te teorem of te preceding subsection were PfA is te Pfaffian of te antisymmetric matrix Az = z F F T I 2M+ iz v I F 2M+ F T 2M+ + zg H 2M+, 25 wit G = K K T and H 2M+ =. 26 Apply to 25 te unitary transformation 4 wit 2M replaced by 2M +. Te transformation diagonalizes F 2M+ F2M+ T as in 5 and, in addition, produces Tus, we obtain were φ m = mπ/2m + 2. Writing U H 2M+ U = M J 2M+. 27 Āz I U AzI U = z F F T + 2z v cos φ m I I 2M+ zg J 2M+ 28 B m z F F T + 2z v cos φ m I, 29 ten te matrix Āz assumes a form sown below in te case of 2M + = 5, Āz = B zg B 3 + zg zg B 4 B 2 zg B 5 zg. 30 Intercanging rows and columns, tis matrix can be canged into a block-diagonal form aving te same determinant, B zg zg B 5 B 2 zg. 3 zg B 4 B 3 + zg For general M we define te 2 2 matrix A m 2 z = B m M+m+ zg and use te result B M+ M+m+ zg, 32 B 2M+ m + zg = 2 + z 2 z + z

7 to obtain Az = Āz = 2 + z 2 z + z 2 It terefore remains to evaluate A m 2 z. M m= A m 2 z. 34 Te matrix A m 2 z can be diagonalized for z = ±iz. To proceed, it is convenient to multiply from te rigt by a 2 2 matrix wose determinant is to obtain Ā m iz A m 0 2 iz I 0 = 2z v cos φ m I 2 + z F F T 0 + i M+m+ z 0 G 0 0 = 2z v cos φ m I 2 + z Q 2 Q 2, 35 were F Q 2 = i M+m K i M+m K F. 36 ow Q 2 commutes wit Q 2 so tey can be diagonalized simultaneously wit te respective eigenvalues e i θ n and e i θ n, were θ n = 2n /2,n =,2, 2. Tus we obtain A m 2 iz = Ām 2 iz = and, taking te square root of 34 wit z replaced by iz, PfAiz = iz /2 /2 M n= m= 4z 2 sin 2 θn + 4zv 2 cos 2 φm, 37 n= 4z 2 2 2n π sin + 4zv 2 cos 2 mπ M + 2 Terefore te generating function vanises identically for = odd. For = even we substitute 38 into and make use of te fact tat te two sets sin 2 2n π/2 and sin 2 4n π/2, n =, 2,, /2 are identical. Tis leads to te result 6. 4 Te Klein bottle A Klein bottle is constructed by inserting extra vertical edges wit weigt z v to te boundary of te Möbius strips of Figs. a and b, so tat tere is a periodic boundary condition in te vertical direction. Te extra edges are oriented upward as te oter vertical edges. Te consideration of te Klein bottle parallels to tat of te Möbius strip. Again, we need to consider te cases of even and odd M separately. 4. M = even Klein bottle For a 2M Klein bottle wit orizontal edges connecting te first and -t row aving weigts z, we ave as in 8 te generating function Z Kln M,z,z v ;z = 7 2M m=0 z m T m. 39

8 Te dimer weigts now generate te antisymmetric matrix A Kln z = Az z v I K 2M K T 2M 40 were Az is given by 9. Following te same analysis as in te case of te Möbius strip, since i - iii still old we find te desired dimer generating function again given by te teorem or, explicitly, ZM,z Kln,z v ;z = ipfa Kln iz + + ipfa Kln iz. 4 2 To evaluate PfA Kln z, we note tat te 2M 2M matrices F 2M K 2M F T 2M + KT 2M and J 2M commute, and can be diagonalized simultaneously by te 2M 2M matrix W wose elements are W mm = W mm = 4M e i2m 2m π/4m i m+m +M e i2m 2m π/4m, 4M e i2m 2m π/4m + i m+m +M e i2m 2m π/4m. 42 We find W F 2M K 2M F T 2M + K T 2MW m,m = 2isin α m δ m,m W J 2M W m,m = i M+m δ m,m, 43 were α m = 2m π/2m. Diagonalizing Az in te 2M-subspace, we obtain were A Kln z = i 2M 2M m= A m z 44 A m z = 2z v sin α m I + z F F T + M+m zk + K T. 45 Tis expression is te same as 7 for te Möbius strip, except wit te substitution of cos φ m by sin α m. Tus, following te same analysis, we obtain M PfAiz = 2z v sin m= n= 2m π 2M + 2i M+m+ z sin 4n π Te substitution of 46 into 4 now gives te result 5. For = even, te Pfaffian 46 is real and 5 reduces to M = odd Klein bottle For a 2M + Klein bottle, te inserted vertical edges ave dimer weigts iz v. Te consideration ten parallels tat of te preceeding sections. Particularly, te desired dimer generating function is also given by 4, but now wit A Kln z = Az + iz v I K 2M+ K T 2M

9 To evaluate PfA Kln z, one again applies in te 2M + subspace te unitary transformation wic diagonalizes F 2M+ K 2M+ F T 2M+ + KT 2M+. Define V mm = V mm = 2M + e im2m π/2m+, 2M + e im 2m π/2m+, m,m =,2,...,2M Using te result V F 2M+ K 2M+ F T 2M+ + K T 2M+V m,m = 2isin ᾱ m δ m,m V H 2M+ V mm = e iᾱm δ m,2m+2 m, 49 were ᾱ m = 2m π/2m+,m =,2,,2M+, ten te matrix Āz = I V AzI V assumes te form in te case of 2M + = 5, Āz = B B 2 ze iᾱ 2 G B 3 + zg ze iᾱ 2 G B 4 ze iᾱ G B 5 ze iᾱ G. 50 Here B m = z Q + 2z v sin ᾱ m. Again, intercanging rows and columns, one canges Āz into te block-diagonal form B ze iᾱ G ze iᾱ G B 5 B 2 ze iᾱ 2 G. 5 ze iᾱ 2 G B 4 B 3 + zg Explicitly, for general M, te m-t block is A m 2 z = B m ze iᾱm G ze iᾱm G B 2M+ m = z F F T I 2 + 2z v sin ᾱ m I 0 0 e iᾱ m zg 0 e iᾱm 0 We proceed as in 35 by multiplying a 2 2 matrix wose determinant is from te rigt, and obtain Ā m 2 iz A m 0 2 iz I 0 = 2z v sin ᾱ m I 2 + z F F T 0 0 e iᾱ m + iz 0 G e iᾱm 0 = 2z v sin ᾱ m I 2 + z Q 2 Q 2,

10 were F Q 2 = ie iᾱm K ie iᾱm K F. 54 ow Q 2 commutes wit Q 2, and tey can be diagonalized simultaneously wit te respective eigenvalues e i β n and e i β n were β n = 2n /2,n =,2,,2. Tus we obtain Ām 2 iz = Using 53 and 33, we get 4z 2 2 2n π sin + 4z 2 2 2m π v sin 2 2M + n= PfAiz = + /2 + iz /2 M 4z 2 2 2n π sin + 4z 2 2 2m π v sin M + n= m= Tus, te generating function 4 vanises identically for = odd. For = even, we replace as before sin 2 2n π/ by sin 2 4n π/. Te substitution of 56 into 4 now leads to te result 7. 5 A reciprocity teorem Using te explicit expression of dimer enumerations on a simple-quartic lattice wit free boundaries, Stanley 2 as sown tat te enumeration expression satisfies a certain reciprocity relation, a relation rederived recently by Propp 3 from a combinatorial approac. Here, we sow tat our solutions of dimer enumerations lead to an extension of te reciprocity relation to enumerations on cylindrical, toroidal, and nonorientable surfaces. We first consider solutions 2, 3, 6 and 7 for = even. Writing 55 T Mob M, = Z Mob M,,, T Kln M, = ZM, Kln, 57 and using te identity 4 n k=0 x 2 2xcos α + 2kπ n + = x 2n 2x n cosnα + 58 repeatedly, we can rewrite our solutions for general z and z v in te form of M+/2 ZM,z Mob,z v = z M/2 x m + x m = m= 3 M 2 z M/2 v /2 n= M+/2 ZM,z Kln,z v = z M/2 t m + t m = m= 3 M 2 z M/2 v /2 n= y M+ n + M yn M y n + yn 59 y M n + y n M + + M, 60 0

11 were x m = G zv z cos mπ z, y n = G sin M + z v 2n π zv, t m = G sin 2 z 2m π M 6 and Gy = y + y 2 +. Tus, te following reciprocity relations are obtained by inspection: were T Mob M, = T Mob M, = ǫ M T Mob M 2, T Kln M, = T Kln M, = ǫ M T Kln M,, 62 Tere are no reciprocity relations for = odd. ǫ M = M, = 2 mod 4 = +, oterwise. 63 We ave carried out similar analyses for dimer enumerations on a simple-quartic net embedded on a cylinder and a torus, using te solutions given in 5,, and ave discovered universal rules of associating reciprocity relations to specific boundary conditions. Generally, tere are 3 different boundary conditions, or matcings, tat can be imposed between 2 opposite boundaries of an M net. Te conditions can be twisted suc as tose sown in te orizontal direction in Fig., periodic suc as on a torus, or free wic means free standing. To establis te convention we sall refer to te boundary condition between te first and te -t columns as te boundary condition in te direction, and similar tat between te first and te M-t rows as in te M direction. Ten, our findings togeter wit tose of Ref. 3 lead to te following teorem applicable to all cases: Reciprocity teorem: Let TM, be te number of close-packed dimer configurations perfect matcings on an M simple-quartic lattice wit free, periodic, or twisted boundary conditions in eiter direction. Te case of twisted boundary conditions in bot directions is excluded. If te twisted boundary condition, if occurring, is in te M direction, we restrict to M = even. Ten, we ave. TM, = ǫ M T 2 M, if te boundary condition in te M direction is free. 2. TM, = ǫ M T M, if te boundary condition in te M direction is periodic or twisted. 6 Summary and Discussions We ave evaluated te dimer generating function for an M simple-quartic net embedded on a Möbius strip and a Klein bottle for all M,. Te results are given by 2-7. Our results can also be written in terms of te q-analogue of te Fibonacci numbers F n q defined by qs s 2 = n=0 Using te first line of 59, for example, and te identity F n q + F n 2 q = x n + x n F n qs n. 64 q x x, 65

12 one can verify tat our results 2 and 6 for te Möbius strip are te same as tose given by Tesler 0. Details of te proof wic also lead to some new product identities involving te Fibonacci numbers will be given elsewere 6. We ave also deduced a reciprocity teorem for te enumeration TM, of dimers on an M lattice under arbitrary including free, periodic, and twisted boundary conditions. Finally, we point out tat te results 2-7 can be put in a compact expression valid for all cases as M+/2 Z M, z,z v = z M/2 Re i 2i M 2 +m+ 4n π sin + 2 X m, 66 2 m= n= were x is te integral part of x, and X m = = zv z zv z cos sin mπ M + 2m π M for te Moebius strip for te Klein bottle. 67 Acknowledgement Work as been supported in part by ational Science Foundation Grant DMR References P. W. Kasteleyn, Pysica 27, H.. V. Temperley and M. E. Fiser, Pil. Mag. 6, M. E. Fiser, Pys. Rev. 24, H. W. J. Blöte, J. L. Cardy, and M. P. igtingale, Pys. Rev. Lett. 56, R. Srock, Pys. Lett. A 26, W. T. Lu and F. Y. Wu, Pys. Lett. A 259, W.-J. Tzeng and F. Y. Wu, Appl. Mat. Lett. 3, W. T. Lu and F. Y. Wu, Pys. Rev. E 63, K. Kaneda, Y. Okabe, Pys. Rev. Lett. 86, G. Tesler, J. Comb. Teo. B 78, T. T. Wu, J. Mat. Pys. 3, R. Stanley, Disc. Appl. Mat. 2, J. Propp, arxiv:mat.co/ I. S. Gradsteyn and I. M. Ryzik, Table of integrals, series, and products, Academic Press, ew York 994,.394 or, specifically,.396., and

13 5 B. M. McCoy and T. T. Wu, Te Two-dimensional Ising Model, Harvard University Press, Cambridge, W. T. Lu and F. Y. Wu, under preparation. 3