Self-dual property of the Potts model in one dimension arxiv:cond-mat/9805301v2 [cond-mat.stat-mech] 26 May 1998 F. Y. Wu Department of Physics Northeastern University, Boston, Massachusetts 02115 Abstract A duality relation is derived for the Potts model in one dimension from a graphical consideration. It is shown that the partition function is selfdual with the nearest-neighbor interaction and the external field appearing as dual parameters. Zeroes of the partition function are analyzed. Particularly, we show that the duality relation implies a circle theorem in the complex temperature plane for the one-dimensional Ising model. 05.50.+q Typeset using REVTEX 1
I. INTRODUCTION It is now well-known that the Potts model [1,2] possesses a duality relation in two dimensions. In fact, it is using this duality relation that Potts [1] first determined the transition temperature for the square lattice. The purpose of this note is to point out a self-dual relation of the Potts model in one dimension, a curious and somewhat surprising result which appears to have escaped heretofore attention. For the Ising model this duality implies a new circle theorem. Consider N Potts spins placed on a ring interacting with nearest-neighbor interactions K and an external field L for one specific spin state. Number the sites by i = 1,2,,N and denote the spin state at site i by σ i = 1,2,,q. The partition function is with σ N+1 = σ 1 and q N Z N (q;k,l) = T(σ i,σ i+1 ) (1) σ i =1 i=1 T(σ i,σ i+1 ) = exp[kδ Kr (σ i,σ i+1 )]exp[lδ Kr (σ i,1)], (2) where δ Kr is the Kronecker delta function. In ensuing discussions we shall make use of the fact that, due to symmetry, the partition function (1) is the same if the external field is applied to any spin state, namely, if δ Kr (σ i,1) in (2) is replaced by δ Kr (σ i,α) for any α = 2,3,,q. Regarding T(σ,σ )astheelements ofaq q matrixt, thenthepartitionfunction(1)can be evaluated by the standard technique of the transfer matrix. The characteristic equation for T is e K+L λ e L... e L 1 e K λ... 1 det..... = 0. (3). 1 1... e K λ This equation can be factorized [4] leading to eigenvalues e K 1, which is (q 2)-fold degenerate, and λ ±, which are the two roots of the quadratic equation λ 2 (e K+L +e K +q 2)λ+e L (e K 1)(e K +q 1) = 0. (4) This gives rise to the following explicit expression for the partition function, Explicitly, we have λ ± = 1 2 Z N (q;k,l) = λ N + +λn +(q 2)(eK 1) N. (5) [ ) 1/2 ] e K+L +e K +q 2± ((e K+L e K q +2) 2 +4(q 1)e L. (6) Note that terms involving branch cuts are cancelled in expanding (5), and as a consequence Z N (q;k,l) is indeed a polynomial in e K and e L. 2
II. THE DUALITY RELATION Our main result is the self-dual relation for the partition function (5) [(e K 1)(e L 1)] N/2 Z N (q;k,l) = [(e K 1)(e L 1)] N/2 Z N (q;k,l ), (7) where the dual variables K and L are related to L and K, respectively, by (e K 1)(e L 1) = q, (e L 1)(e K 1) = q. (8) The validity of (7) can be explicitly verified using (5) for q = 2 and for specific values of N = 2,3,. But its validity for general q and N is not very obvious by looking at the solution. To establish (7) for general q and N, we construct a lattice L of N + 1 sites as shown in Fig. 1, by introducing an extra (ghost) spin interacting with all N sites with an equal interaction L. Note that the lattice L is planar and self-dual. Namely, if one places dual spins, one in each face of L including the face exterior to L, and connects dual spins with edges crossing each of the edges of L, one arrives at a lattice which is precisely L. Let Z L (q;k,l) be the partition function of the Potts model on L. Following the remark after (2), one has the identity Z L (q;k,l) = q Z N (q;k,l). (9) It is then sufficient to show that (7) holds for Z L (q;k,l). Generally, the Potts model possesses a duality relation for any planar lattice, or graph, with arbitrary edge-dependent interactions [3]. Let Z(K ij ) be the partition function of the Potts model on a planar graph of M sites with edge-dependent interactions K ij, and Z (D) (K ij) the partition function of the dual model. Then, the duality relation given in [3] can be written as q M/2 Z(K ij ) = q MD/2 Z (D) (Kij), (10) e K ij 1 e K ij 1 edges edges where M D is the number of sites of the dual graph and (e K ij 1)(e K ij 1) = q. (11) Now the lattice L is self-dual, namely, M = M D = N + 1 and Z = Z (D) = Z L. The application of (10) to L now leads to an expression which is precisely (7) with Z L in place of Z N. The duality relation (7) now follows after introducing (9). III. PARTITION FUNCTION ZEROES The zeroes of the partition function can be computed from λ N + +λ N +(q 2)(e K 1) N = 0, (12) 3
at least numerically, for any given N. The analysis is particularly simple if the three terms in (12) coalesce into two. For the q = 2 Ising model, for example, one finds zeroes located at λ + = e i(2n+1)π/n λ, n = 0,1,2,,N 1. (13) This leads to the Yang-Lee circle [5] in the complex e L plane for K 0. For the zero-field Potts model, L = 0, one has λ = e K 1 and finds from (12) the zeroes at e K +q 1 = (1 q) 1/N (e K 1), n = 0,1,2,,N 1. (14) In the limit of N, the partition function (5) is dominated by the eigenvalue with the largest magnitude. It follows that the zeroes lie continuously on the loci λ + = λ, in λ + e K 1 λ ± = e K 1, in λ ± λ, (15) generalizing a previous result [7] for real q and L = 0. Here, the loci (14) and (15) apply to real and complex K,L,q. Particularly, for L = e K = 0, Z N (q;k,l) gives the ground state entropy of the antiferromagnetic Potts model and the chromatic polynomial (in q) of a ring. In this case we have λ + = q 1, λ = e K 1 = 1, and we find the zeroes on the circle q 1 = 1 (16) in the complex q plane. This result has previously been reported by Shrock and Tsai [8]. For the Ising model the duality relation (7) provides some further implications on the partition function zeroes. To conform with usual notations, we consider N Ising spins on a ring with nearest-neighbor interactions K I and an external magnetic field L I 0. The Ising partition function is Z Ising (K I,L I ) = N [e K Iσ i σ i+1 e L Iσ i ]. (17) σ i =±1 i=1 Using the identity σ i σ i+1 = 2δ Kr (σ i,σ i+1 ) 1 and σ i = 2δ Kr (σ i,1) 1, we relate Z Ising to the 2-state Potts partition function via Z Ising (K I,L I ) = e N(L I+K I ) Z N (2;2K I,2L I ). (18) Thus, after substituting (18) into (7) and some re-arrangement, we obtain the following self-dual relation for the Ising model on L, where the dual variables are Z Ising (K I,L I ) (sinh2k I sinh2l I ) = Z Ising (KI,L I ) N/4 (sinh2ki sinh2l N/4, (19) I) e 2K I = tanhl I, e 2L I = tanhk I. (20) We note the resemblance of (19) with the duality relation of the Ising model formulated by Syozi [6]. 4
FromtheYang-Leecircle theorem [5] we know that, fork I 0, zeroes of Z Ising (K I,L I ) lie on the unit circle in the complex field e 2L I plane. Now, since we have KI 0, the duality relation (19) implies that zeroes of Z Ising (K I,L I ) lie also on the unit circle in the complex temperature tanhk I plane, a result which holds for any L I 0 and which is not readily recognized otherwise. A consequence of this is that zeroes of Z Ising (K I,L I ) lie on the pure imaginary axis in the e 2K I plane, generalizing a result previously known only for L I = 0 [7]. Our circle theorem contrasts with the square lattice result, valid only in the thermodynamic limit and in zero field, which states that in the complex tanhk I plane zeroes of the partition function Z SqIsing (K I ) lie on two circles of radius 2 and centered at ±1 [9]. Note added: After the completion of this manuscript, it has been pointed out to me that results similar to those described here have been reported previously [10]. The duality relation(7), which is derived in this manuscript from a graphical consideration, was obtained in [10] algebraically by considering the transfer matrix. I would like to thank A. Glumac for calling my attention to [10]. IV. ACKNOWLEDGEMENT Work has been supported in part by National Science Foundation Grant DMR-9614170. I would like to thank Professor R. Shrock for calling my attention to Refs. [7,8]. FIGURE CAPTION Fig. 1. The lattice L with N +1 sites. 5
REFERENCES [1] R. B. Potts, Proc. Camb. Phil. Soc. 48, 106 (1952). [2] F. Y. Wu, Rev. Mod. Phys. 54, 235 (1982). [3] F. Y. Wu and Y. K. Wang, J. Math. Phys. 17, 439 (1976). [4] The characteristic equation (3) is factorized by subtracting the second row from the third and each of the subsequent rows, followed by adding the third and all subsequent columns to the second column. [5] C. N. Yang and T. D. Lee, Phys. Rev. 87, 404 (1952); T. D. Lee and C. N. Yang, Phys. Rev. 87, 410 (1952). [6] I. Syozi, in Phase Transitions and Critical Phenomena, Vol. 1, Edited by C. Domb and M. S. Green (Academic Press, New York, 1972). [7] V. Matveev and R. Shrock, Phys. Lett. A204, 353 (1995). [8] R. Shrock and S.-H. Tsai, Phys. Rev. E55, 5165 (1997). [9] M. E. Fisher, in Lecture Notes in theoretical Physics, edited by W. E. Brittin(University of Colorado Press, Boulder, 1965), Vol. 7c. [10] A. Glumac and K. Uzelac, J. Phys. A 27, 7709 (1994). 6
N 1 L K L 2 K 3 4 Fig. 1