15 June 2000
Physics Letters B 483 Ž2000. 303–308
The field theory of the q
™4
q
Potts model
G. Delfino a , John Cardy b,c a
Laboratoire de Physique Theorique et Hautes Energies, UniÕersite´ Pierre et Marie Curie, ´ Tour 16 1e r etage, 4 place Jussieu, 75252 Paris cedex 05, France ´ b Theoretical Physics, UniÕersity of Oxford, 1 Keble Road, Oxford OX1 3NP, UK c All Souls College, Oxford, UK Received 16 February 2000; accepted 28 April 2000 Editor: L. Alvarez-Gaume´
Abstract
™
The q-state Potts model in two dimensions exhibits a first-order transition for q ) 4. As q 4q the correlation length at this transition diverges. We argue that this limit defines a massive integrable quantum field theory whose lowest excitations are kinks connecting 4 q 1 degenerate ground states. We construct the S-matrix of this theory and the two-particle form factors, and hence estimate a number of universal amplitude ratios. These are in very good agreement with the results of extrapolated series in qy1r2 as well as Monte Carlo results for q s 5. q 2000 Elsevier Science B.V. All rights reserved.
The q-state Potts model, defined by the lattice Hamiltonian w1x H s yJ
Ý
d sŽ x ., sŽ y. ,
² x , y:
where the spin variable sŽ x . assumes q different values Žcolours., continues to be of fundamental importance in the description of a large variety of critical phenomena, ranging from ferromagnetism to percolation and adsorbed monolayers w2x. It was shown by Baxter that in two dimensions the ferromagnetic model undergoes a phase transition which is continuous for q F 4 and first order otherwise w3x. The Coulomb gas w4x and conformal field theory ŽCFT. w5x provided later a complete description of the second order phase transition line, which turned out to correspond to CFT’s with central charge c F 1, the value c s 1 corresponding to the end point q s 4. More recently, integrable field theory has lead to new results for the scaling limit of the off-critical model for q F 4 w6–8x.
Concerning the first order transition for q ) 4, several exact lattice results – internal energy w3x, magnetisation w9x, correlation length w10x – are known, but progress through field theoretic methods is generally prevented by the absence of a scaling limit. When q approaches 4 from above, however, the correlation length at Tc diverges and a continuum description in terms of a massiÕe quantum field theory should be possible. The identification and solution of the quantum field theory describing the limit q 4q at Tc , as well as the determination of the universal critical quantities, are the subject of this note. The correlation length at criticality of the lattice model is known to behave as w10x
™
j ; aep
2
™
'
r qy4
Ž 1.
q
as q 4 , where a is proportional to the lattice spacing, so that the continuum limit corresponds to taking the limits a 0, q 4q in such a way that j
0370-2693r00r$ - see front matter q 2000 Elsevier Science B.V. All rights reserved. PII: S 0 3 7 0 - 2 6 9 3 Ž 0 0 . 0 0 5 4 8 - 7
™ ™
G. Delfino, J. Cardyr Physics Letters B 483 (2000) 303–308
304
remains finite. The presence of an essential singularity rather than a power law divergence in Eq. Ž1. is characteristic of a perturbing operator which is only marginally relevant Žscaling dimension 2. – except that in the original description of the Potts model q is a parameter which should be invariant along RG trajectories and therefore such an interpretation needs to be treated with some care. In fact we know that at q s 4 there is such a marginal field c : it was shown by Nienhuis et al. w11x to correspond to the fugacity for vacancies in the lattice model. When c - 0 the transition is second order, with logarithmic modifications to scaling arising from the marginal irrelevance of c , but when c ) 0, the transition is first order with a correlation length at the transition which diverges as c 0q like 1r bc
™
j ; ae , Ž 2. where b is a constant. The main point is that the scaling limits in which a 0 with j fixed are identical in these two cases. This may be seen, for example, from the general structure of the RG equations near q s 4 and c s 0. Based on the analysis of Nienhuis et al. w11x, Cardy, Nauenberg and Scalapino w12x argued that these have, to lowest order, the general form
™
d crdl s bc 2 q a Ž q y 4 . q O Ž c 3 , Ž q y 4. c , Ž q y 4.
2
.,
dtrdl s Ž yt q c t c . t q O Ž c 2 t , Ž q y 4 . t ,t 2 . ,
Ž 3. Ž 4.
2
dhrdl s Ž y h q c h c . h q O Ž c h, Ž q y 4 . h,th . ,
Ž 5. where t is the deviation from the critical temperature and h is a symmetry-breaking field. a, b, c t and c h are all constants, certain combinations of which are universal w12x. An important feature of Ž4., Ž5. is that, to lowest nontrivial order, they do not involve q. This is because, as may be seen from the first equation, q y 4 is effectively O Ž c 2 .. Integrating Ž3. up to a value l˜ such that c Ž l˜. s O Ž1. then gives ˜ results for the correlation length j ; e l in agreement with Ž1., Ž2. in the two cases q ) 4 and q s 4, c ) 0. The various thermodynamic quantities are then found by integrating the other equations up to l˜; ln j . To the order stated, the results will be identical in the two cases, when expressed in terms of j . Therefore the scaling limits are identical.
Another way of understanding this is through the mapping of the lattice Potts model to a height model and thence to a Coulomb gas or sine-Gordon theory w4x. From this point of view, q is merely a parameter identified with a certain function of the coupling constant conventionally called b . q s 4 corresponds to b 2 s 8p at which point the operator corresponding to c becomes marginal. Within the standard RG picture of the sine-Gordon model, both b and c have non-trivial marginal flows. However the scaling limit, corresponding to the massive sine-Gordon theory at b 2 s 8p , is unique, and therefore describes both the cases q 4q and c 0q at q s 4. Having made this observation, we may take over the results of Ref. w13x in which it was pointed out that the scaling limit of the massive q s 4 theory is integrable, and in which the scattering theory and form factors were determined. In our case, the construction of the scattering theory goes as follows. Along the first order phase transition line, q ordered ground states are degenerate with the disordered ground state. The field theory describing the scaling limit q 4q has 4 ordered vacuum states V i , i s 1, . . . ,4. Invariance under colour permutations implies that, in the order parameter space, they lie at the vertices of a tetrahedron having the disordered vacuum V 0 at its center. The elementary excitations of the scattering theory are stable kinks K 0 i Ž u ., K i0 Ž u . interpolating between the center of the tetrahedron and the i-th vertex, and vice versa. We denote by u the rapidity variable parameterising the on-shell momenta as Ž p 0 , p 1 . s Ž mcosh u ,msinh u .. The mass of the kinks m ; jy1 measures the deviation from the conformal point q s 4. The space-time trajectory of a kink on the plane draws a domain wall separating a coloured phase from the disordered one. The space of asymptotic states is made of multi-kink sequences in which adjacent vacuum indices belonging to different kinks have to coincide. For example, up to possible bound states, the lightest excitation interpolating between two ordered vacua is K i0 Ž u 1 . K 0 j Ž u 2 .. The factorisation of multi-kink processes reduces the scattering problem to the determination of the two-kink amplitudes. Colour permutation symmetry allows only for the four elementary processes depicted in Fig. 1. The four amplitudes can be determined as a solution of the requirements of unitarity
™
™
™
G. Delfino, J. Cardyr Physics Letters B 483 (2000) 303–308
305
tance expansion, we will approximate the Žconnected. two-point correlator of a scalar operator F Ž x . as Fig. 1. The two-kink scattering amplitudes A 0 , A1 , B0 , B1 Ž i/ j ..
² V 0
Žcrucially simplified by the absence of particle production., crossing and factorisation. The scattering amplitudes of Fig. 1 are given in w13x and read A0 Ž u . s
eyi gu 2 ip y u ip y u
2 yi gu
A1 Ž u . s
e
B0 Ž u . s e igu B1 Ž u . s e
Ž 6.
S0 Ž u . s
ž
u
2p 2p
< F0Fi Ž u 1 y u 2 . < 2 ey< x < E 2 ,
Ž 11 . in the disordered phase, and
,
Hu )u 1
Ž 8.
S0 Ž u . ,
1 2
is1 u 1) u 2
Ž 7.
Ž 9.
where u is the rapidity difference of the two kinks, g s p1 ln2, and
G
du 1 du 2
ÝH
² V i
u
S Žu . , ip y u 0 ip q u S Žu . , ip y u 0
2
igu
S0 Ž u . ,
,
du 1 du 2 2
2p 2p
< Fi0F Ž u 1 y u 2 . < 2 ey< x < E 2 ,
Ž 12 .
in the i-th ordered phase. Here E2 s mŽcosh u 1 q cosh u 2 . is the energy of the two-kink asymptotic state and we introduced the two-kink form factors F0Fi Ž u 1 y u 2 . s ² V 0
Ž 13 .
u
/ž / /ž / G y
q
2 ip 2 ip u u G 12 y G 2 ip 2 ip
ž
° ~ s yexp iH ¢
`
0
x
dx
y
e
2
x cosh x 2
¶ xu sin •. pß
Fi0F Ž u 1 y u 2 . s ² V i
Ž 14 .
Ž 10 .
The absence of poles in the physical strip Im u g Ž0,ip . ensures that there are no bound states and that the four amplitudes above completely determine the scattering theory. This S-matrix shares evident analytic similarities with that of the SUŽ2.-invariant Thirring model. As a matter of fact, the latter is a realisation of the same perturbed CFT on a different particle basis Žan SUŽ2. doublet rather than our kinks.. The possibility for a single perturbed CFT to be invariant under different symmetry groups and to describe different universality classes is discussed, for example, in Ref. w14x. Making contact with the thermodynamics requires the computation of correlation functions. In our Smatrix framework these are obtained as spectral series summing over all multi-kink intermediate states. Neglecting terms of order ey4 m < x < in this large dis-
The operators of interest for us are the spin sj Ž x . s d sŽ x ., j y 1rq, and the energy ´ Ž x . s Ý y d sŽ x ., sŽ y. , whose scaling dimensions around the q s 4 fixed point are Xs s 1r8 and X´ s 1r2 w4x. Some consequences for the physics of the coexisting phases follow immediately from the structure of the scattering theory. The ‘true’ correlation length j is determined by the large distance decay of the spin-spin correlator as ² sj Ž x . sj Ž0.: ; ey< x < r j . Then it follows from Ž11., Ž12. that
j o s j d s 1r2 m
Ž 15 .
Žhere and below the subscript o Ž d . denotes quantities computed in the ordered Ždisordered. phase.. Numerical simulations w15x and large q expansions w16x suggest that the phase independence of j holds true for all q ) 4. Since the interfacial tension between two coexisting phases is given by the total mass of the lightest excitation interpolating between them, we also have so d s m and so d s so or2. The latter result is known to hold for all for q ) 4 w17x.
G. Delfino, J. Cardyr Physics Letters B 483 (2000) 303–308
306
The relation of the other interesting thermodynamic quantities Žspontaneous magnetisation M, latent heat L, susceptibility x , specific heat C, second moment correlation length j 2nd . with the connected correlators of sj and ´ , and their behaviour as q 4q are
™
M s ² V j < sj < V j : , B jy1 r8 ,
Fq Ž u .
½
`
s exp 2
dx
H0
2
x s d x ² sj Ž x . sj Ž 0 . : , Go , d j
7r4
,
Ž 18 .
sin2 Ž 2 ip y u .
x 2p
° =exp~yH ¢
`
2
x x 2 ey 2 sin Ž ip y u . 2p x cosh x sinh x 2
dx
C s d 2 x ² ´ Ž x . ´ Ž 0. : , A o , d j ,
H
1 4x
H
2
Ž 19 . For the ordered case, the two-point correlators of sj entering x and j 2nd are computed on the vacuum < V j :. Since ´ Ž x . is odd under the duality transformation exchanging the low- and high-temperature phases, we have Ad s Ao .
Ž 20 .
The critical amplitudes are normalisation dependent but can be combined into a series of universal ratios characterising the scaling limit. We can evaluate the critical amplitudes by integrating the two-particle approximations Ž11., Ž12. of the correlators. What we need to know are the two-kink form factors of the operators ´ and sj . Once again the result is contained in Ref. w13x and reads g ´ Fi0,0 i
Ž u . s . iL
e " 2 Žp qi u .
u y ip
F0 Ž u . ,
Ž 21 .
The results we obtain for the universal amplitude ratios are given in Table 1 and compared with those following from the combination of the exact w3,9,10x and series w18x lattice results for the amplitudes. The accuracy exhibited by the two-particle approximation does not come as a surprise since it is known as a common feature of this kind of computations within integrable field theory. In the present case the accuracy is enhanced by the low scaling dimensions of the spin and energy operators which lead to mild singularities for their correlators and then to a small contribution of short distances to the integrals. We estimate that the errors on our values for the amplitude ratios do not exceed order 0.1%. The scaling limit we discussed so far corresponds to q 4q. At q s 5, however, the correlation length is still some 2500 times the lattice spacing and this suggests that our results for q 4q could still provide the basis for an approximate description. For a generic value of q ) 4, the model has q q 1 degenerate ground states at the transition point and the elementary excitations are 2 q kinks going from
™
Table 1 Universal amplitude ratios for the q-state Potts model at T sTc , q 4q. The field theoretical results are obtained within the two-particle approximation.
g
s. M
4d i j y 1 6Fq Ž ip .
e
2
™
™
sj Fi0,0 iŽu . "
Žp qi u .
u cosh
¶ •. ß Ž 24 .
d 2 x < x < 2² sj Ž x . sj Ž 0 . : , Ž f o , d j . .
² ´ :d s y² ´ :o ,
,
u
0
2 j 2nd s
5
Ž 23 . F0 Ž u . s yisinh
Ž 17 .
eyx
x sinh2 x
Ž 16 .
L s ² ´ :d y ² ´ :o , L jy1 r2 ,
H
with Fy Ž u . s Fq Ž u q 2 ip .,
F " Ž u . F0 Ž u . ,
2
Ž 22 .
fd fo r fd Gd r Go Ad r L 2 Gd r B 2
Field theory
Lattice
0.6744 0.9340 1.1406 0.1047 0.06607
0.673Ž8. 0.935Ž5. 1.19Ž5. 0.105Ž3. 0.0656Ž15.
G. Delfino, J. Cardyr Physics Letters B 483 (2000) 303–308
the disordered vacuum to the q ordered vacua, and vice versa. If the correlation length is sufficiently large, an approximate scaling should still hold and then it makes sense to keep for the physical quantities the parameterisations Ž16–19., namely a power of the correlation length times an amplitude. It is easy to see, however, that in the present case we have to allow for a q-dependence of the amplitudes. Consider in fact the correlator
307
Table 2 Values obtained combining the exact results of Refs. w9,10x and the large q expansions of Ref. w18x. q
j Žlattice units. Žd. j 2nd rj Ž o. Žd. j 2nd r j 2nd R1 s Ž q y1.r q x d r x o R 2 s x o rŽ Mj . 2
4q
5
10
0.673Ž8. 0.935Ž5. 0.89Ž4. 0.0550Ž6.
2512.24 0.671Ž3. 0.934Ž7. 0.810Ž5. 0.0589Ž2.
10.559 0.6587Ž1. 0.9579Ž2. 0.80399Ž1. 0.05784Ž1.
Ga j Ž x . s ² Va < sj Ž x . sj Ž 0 . < Va : ,
a s 0,i ;i , j s 1, . . . ,q .
Ž 25 .
Its two-particle approximation in the disordered phase is
xT
q
G0 j Ž x . ,
Let us conclude this note by considering the ‘transverse’ susceptibility x T w8x obtained integrating Gi j Ž x . rather Gj j Ž x .. From Ž27., in the two-particle approximation
xo
Ý < F0si < 2 eyE < x < j
2
is1
q s qy1
< F0sjj < 2 eyE 2 < x < ,
Ž 26 .
where integration over momenta is understood. The last equality follows from colour symmetry and Ý j sj s 0 which imply F0sij,i0 s Ž qd i j y 1 . r Ž q y 1 . F0sj,j j0 .
Ž 27 .
When integrating the correlator to obtain the amplitude of, say, the susceptibility in the disordered phase, the form factors computed at q s 4 should give a good approximation as long as j is large. The explicit factor qrŽ q y 1. dictated by the number of intermediate states and symmetry, however, has to be taken into account and is expected to determine the main deviation from the q 4q value of the amplitude. Following the same reasoning, the susceptibility amplitude in the ordered phase should be basically constant in q since there is only one intermediate state. More generally, we are led to expect that the ratios listed in Table 2 are approximately constant in q for j large enough so that the continuum description is accurate. Their values determined from the results of the large q expansion and reported in the Table seem to confirm our picture. Our field theory results for R 1 and R 2 as q 4q are 0.855 and 0.0579, respectively.
™
™
1 ,
Ž q y 1.
2
.
Ž 28 .
This result is basically a consequence of the nature of the elementary excitations and is expected to hold as a good approximation for all q ) 4 at Tc , as long as j 4 a. We are not aware of lattice results on x T for comparison. In summary, we have shown that the limit q 4q in the Potts model defines an integrable massive field theory, whose S-matrix and form factors may be computed exactly. The results for integrated correlation functions are in excellent agreement with lattice-based numerical results. This shows how methods of continuum field theory are not restricted to the description of second-order transitions only.
™
References w1x R.B. Potts, Proc. Cambridge Philos. Soc. 48 Ž1952. 106. w2x F.Y. Wu, Rev. Mod. Phys. 54 Ž1982. 235. w3x R.J. Baxter, Exactly solved models of statistical mechanics, Academic Press, London, 1982. w4x B. Nienhuis, J. Stat. Phys. 34 Ž1984. 781. w5x Vl.S. Dotsenko, V.A. Fateev, Nucl. Phys. B 240 Ž1984. 312. w6x L. Chim, A.B. Zamolodchikov, Int. J. Mod. Phys. A 7 Ž1992. 5317. w7x G. Delfino, J.L. Cardy, Nucl. Phys. B 519 Ž1998. 551. w8x G. Delfino, G.T. Barkema, J.L. Cardy, Nucl. Phys. B 565 Ž2000. 519. w9x R.J. Baxter, J. Phys. A 15 Ž1982. 3329. w10x E. Buffenoir, S. Wallon, J. Phys. A 26 Ž1993. 3045. w11x B. Nienhuis, A.N. Berker, E.K. Riedel, M. Shick, Phys. Rev. Lett. 43 Ž1979. 737.
308
G. Delfino, J. Cardyr Physics Letters B 483 (2000) 303–308
w12x J.L. Cardy, M. Nauenberg, D.J. Scalapino, Phys. Rev. B 22 Ž1980. 2560. w13x G. Delfino, Nucl. Phys. B 554 Ž1999. 537. w14x G. Delfino, hep-thr9911192. w15x W. Janke, S. Kappler, Nucl. Phys. Proc. Suppl. B 42 Ž1995. 770.
w16x H. Arisue, hep-latr9909166. w17x C. Borgs, W. Janke, J. Physique I ŽFrance. 2 Ž1992. 2011. w18x H. Arisue, hep-latr0001010. H. Arisue, K. Tabata, Nucl. Phys. B 546 Ž1999. 558.