Arxiv:0801.3565V2 [Cond-Mat.Str-El] 25 Sep 2008 ∗ Als.As H U3 Esnegmdlcnb Di- Be Can Model Is Heisenberg It That SU(3) Showing the Chain, Also Spin SU(3) Gapless

Home , Spin model

arXiv:0801.3565v2 [cond-mat.str-el] 25 Sep 2008 ∗ als.As h U3 esnegmdlcnb di- be can model is Heisenberg it that SU(3) showing the chain, Also spin SU(3) gapless. systems the gener- component include who multiple which to [8], of Bethe-Ansatz Sutherland SU(3) family the by with alized the shown models first side some as other include symmetry, the chains AdS/CFT On the spin in integrable chains 7]. implication spin [6, their by SU(N) correspondence boosted ferromagnetic been on has recent interest In the viewpoints. different years from studied been also has particular a ef- with low-energy chain. spin identified the tube SU(3) where be spin anisotropic [5] can the field or Hamiltonian magnetic 4], fective [3, a trimer chain in SU(3) tetrahedron models the spin color involve a each in examples of state Other particles of with conserved. number language), is the SU(3) SU(3) that the avail- the condition in three to the (colors, the related by states be spanned atomic would space able internal spin an the in proposed case rotation recently this been In has lattices [2]. optical ultra- in SU(3) of atoms of systems cold in realization chains or- possible spin the A (AFM) to antiferromagnetic coupled freedom. is of spin degrees electron met- bital the transition where in [1] laboratories oxides in als effec- SU(4) realized be particular be can In systems can only experimentally. also not implemented represent but tively has models problems SU(2) these theoretical than because challenging larger is This symmetry a arisen. with matter densed lcrncades [email protected] address: Electronic rmatertclpito iw h U3 pnmodel spin SU(3) the view, of point theoretical a From con- of models in interest renewed a years, recent In eatmnod ´sc eoia autdd inisUni Ciencias F´ısica de Teórica, Facultad de Departamento iatmnod iiadl’nvri` iBlgaadINF and Bologna dell’Università di Fisica di Dipartimento ASnmes 46.,6.0a,1.5H,89.70.Cf 11.25.Hf, 64.60.ae, 64.60.F, numbers: PACS nrycnomlfil hoyi h SU(3) exc the first is the theory i of field dimensions model conformal scaling the energy the that and confirm charge central results its The (DMRG). Group malization a-lnkIsiu ¨rQatnpi.Hans-Kopfermann für Quantenoptik. Max-Planck-Institut eaayeteatfroantcS()Hiebr hi by chain Heisenberg SU(3) antiferromagnetic the analyze We .INTRODUCTION I. MGSmlto fteSU the of Simulation DMRG erth ¨rTertsh hskI Universität Dortmu I, Physik für Theoretische Lehrstuhl toHh tae4 42 otud Germany Dortmund, 44221 4, Straße Otto-Hahn .Eclsi .Ortolani F. Ercolessi, E. Dtd coe 1 2018) 31, October (Dated: 1 .Aguado M. .Asorey M. .Pasini S. esZmn-oio-itnmodel. Wess-Zumino-Novikov-Witten (3) hc sas nw ob rtcl[,1] ntrsof SU(3) universality the the terms same the and of In to that model belong 10]. class, LS should [9, chain the spin critical (CFT) SU(3) be Theory to Field known Conformal also model, bilinear is (LS) which SU(3)-symmetric Lai-Sutherland the particular chain, spin-1 a biquadratic to related rectly oarewt hs fteSU(3) the of out those turn with which the agree excitations compute first to We the irre- of the SU(3). dimensions as scaling of the numbers (IR) of quantum representations states same excited ducible the the have that numbers chain show the quantum can spin to the we provide state, due each to particular, for program In correspondence our its of model. as ability confirm well Lai-Sutherland which as the III) model to the (Section of results criticality II), new the (Section our framework present mathematical of we its description and short model Renor- a Matrix the After Density (DMRG). the Group of malization means by chain spin SU(3) 12]. [11, model (WZNW) ftevcu nageetentropy. entanglement means con- vacuum by charge further the central are the of de- of results computation theory The the by field chain. firmed spin effective our energy of low scriptions the to corresponds hr h pnvralsaeepesdi em fthe representation: of fundamental terms the in in expressed SU(3) are of variables generators spin the where ∗ nti ae epeetanmrclaayi fthe of analysis numerical a present we paper this In ecnie h olwn esnegmodel Heisenberg following the consider We F esnegModel Heisenberg AFM esddd aaoa 00 aaoa Spain Zaragoza, 50009 Zaragoza, de versidad ,VaInro4,417Blga Italy Bologna, 40127 46, Irnerio Via N, tdsae hwta h neliglow underlying the that show states ited Sr .D878Grhn,Germany Garching, D-85748 1. -Str. en fteDniyMti Renor- Matrix Density the of means rtcladtecmuainof computation the and critical s I H U3 MODEL SU(3) THE II. H = nd, J 1 X i Wess-Zumino-Novikov-Witten =1 L S i · 1 S i ZWmdlwhich model WZNW +1 (1) 2

Sa 1 a ¯ ¯ i = 2 λi , with a = 1, .., 8 and λa the eight Gell-Mann λ λ (Sz,Qz) xλ,λ matrices. The sign of J selects respectively an antifer- (1) (1) (0,0) 0 romagnetic spin chain (J > 0) or a ferromagnetic (FM) (±1/2, 1/2) (3) (1) 1/3 one (J < 0). In the following we shall concentrate only  (0, −1) on the AFM case, which has been partially considered (±1/2, −1/2) (3¯) (1) 1/3 also in ref. [13, 14].  (0, 1) In terms of the following ladder operators, T ± = λ1 iλ2, (±1/2, ±3/2) V ± = λ4 iλ5 and U ± = λ6 iλ7, the Hamiltonian± (1) (3) (3¯) 8 (0, 0) (3 times) 2/3 ± becomes ± ± < ( 1, 0) : Table I: Quantum numbers and scaling dimensions for some L of the primary fields Φλ,λ¯ of the SU(3)1 WZNW model. J 1 − − − H = T +T + V +V + U +U + h.c. 2 4 i i+1 i i+1 i i+1 Xi=1 1 3 3 1 8 8 + λi λi+1 + λi λi+1 . (2) algebra su(3) is through the Young Tableau (YT) which 2 2 can be labelled by two positive integer numbers (p, q). Once p and q are known, one can easily compute the di- This makes easier to identify two operators, Sz and Qz, given by the sums of the two diagonal Gell-Mann matri- mension d of the representation and the quantum num- ces bers associated to the isospin Sz and the hypercharge Qz according to [16, 17]: 1 √3 S = λ3 Q = λ8, (3) 1 z 2 i z 2 i d = (p + 1)(q + 1)(p + q + 2) (6) Xi Xi 2 that commute with the Hamiltonian and correspond to and conserved quantities (isospin and hypercharge). The corresponding quantum numbers label the different eigen- Sz = I, I +1, ..., I 1, I − − − (7) states of (1). 3 The Lai-Sutherland model is defined as the bilinear Qz = 2 Y biquadratic spin-1 chain 1 2 where I = 2 (r + s) and Y = (r s) 3 (p q), with 0 r p , 0 s q . In particular,− − the cases− (1, 0) L and≤ (0,≤1) give≤ respectively≤ the fundamental (3) and the ′ 2 ¯ H = J S˜i S˜i+1 + (S˜i S˜i+1) , (4) anti-fundamental (3) IR, while the singlet representation · · Xi=1 h i (1) corresponds to (0, 0). It has been proved [18] that, in analogy with the SU(2) and characterized by an SU(3) symmetry. The model (4) case, the ground state (GS) of the AFM SU(3) Hamil- and the SU(3) spin chain can be mapped one onto the tonian is a singlet and, since it is made of particles u, other by means of the following identity [15] d and s in equal number, it can be obtained in finite chains having only a number of sites which is a multiple 8 of three, L = 3M. As for the excited states, we expect 2 1 1 a a S˜iS˜i + (S˜iS˜i ) 1= + λ λ . (5) them to be in correspondence with the tower of conformal +1 +1 − 3 2 i i+1 Xa=1 states of the corresponding SU(3) WZNW model. The primary states of this theory are a finite number and are We have already mentioned in the introduction that the given [16] by fields Φλ,λ¯ , whose holomorphic (antiholo- LS model is known to be gapless and to belong to the morphic) part transforms according to a representation same universality class of the SU(3) level-1 Wess-Zumino- λ = (p, q) (λ¯ = (p′, q)) with the values of p, q (and simi- Novikov-Witten model with central charge c = 2. Due to larly of p′, q′) satisfying the condition: p + q k, k being the correspondence between the two models, the SU(3)1 the level. The conformal dimension of the primary≤ field WZNW model has to be the low energy effective criti- is then x ¯ = x + x ′ ′ with cal field theory also for the SU(3) spin chain. We shall λ,λ (p,q) (p ,q ) numerically show that the SU(3) Heisenberg chain is crit- 1 x = (p2 + q2 + pq +3p +3q), (8) ical, and from the energy state obtained from the DMRG, (p,q) 3(k + N) we shall compute the central charge and the scaling dimensions of (1) and compare them to the values predicted and a similar expression for x(p′,q′). For future reference, for the SU(3)1 WZNW model. the values of xλ,λ¯ for some primary fields in the case of The states of the spin chain can be organized accord- k = 1 are reported in Table I. ing to the irreducible representations of the affine (Kac- To end up this section, we notice that in a finite chain Moody) Lie algebra associated to SU(3). Let us recall of length L not all quantum numbers, i.e. states, may be [16] that a useful way of representing the IR’s of the Lie realized. For examples, working with periodic boundary 3 conditions and with an even number of sites, the singlet -0.49 S =0.0 Q =0.0 L=6M (1) (1) (ground) state, with x = 0, appears only for z z S =0.0 Q =1.0 L=6M+2 × z z chains with L = 6M (with M a positive and integer S =0.0 Q =-1.0 L=6M+4 ¯ z z number), while the (3) (1) (or the (1) (3)) states are S =0.5 Q =1.5 L=6M × × z z present only if L =6M +4 (or L =6M + 2), both with -0.5 Sz=0.0 Qz=2.0 L=6M+2 S =0.0 Q =-2.0 L=6M+4 x =1/3. z z Sz=1.0 Qz=0.0 L=6M

III. NUMERICAL ANALYSIS -0.51 (L)/L j E The SU(3) version of the DMRG we have used implements the following Hamiltonian

L -0.52 J 1 + − + − H = k0 Ti Ti + k1 Vi Vi (9) 2 4 +1 +1 E Xi=1 00 − 1 + k U +U + h.c.)+ z λ3λ3 + z λ8λ8 2 i i+1 2 0 i i+1 1 i i+1 -0.53 0 0.001 0.002 0.003 0.004 π 2 where kj and zj are input parameters. The model (9) /6L reproduces the AFM (FM) case when all the kj ’s and the zj ’s are equal to 1 (-1). By tuning the input parameters k and z , we can study all the possible anisotropic Figure 1: Plot of the ground state E00 and of the excited j j states for chains of different lengths (from L = 12 up to L = version of the SU(3) Heisenberg model. A very impor- 1 2 52). The solid lines have a slope respectively of 3 and 3 and tant feature of this DMRG is that it implements both have been drawn as a guide for the eye. For sake of clarity, the quantum numbers Sz and Qz given in (3). This im- not all the degeneracies have been reported. plementation considerably reduces the computation time and, on the other hand, once Sz and Qz are fixed from input, each run of the DMRG yields exclusively the ener- If the effective field theory describing our spin chain is gies of the states within those quantum-number sectors. the conformal SU(3)1 WZNW model, the central charge This is very useful when one needs to classify the exci- must be c = 2. tations according to the values of the isospin and of the However, it is possible to have a direct numerical hypercharge. derivation of c from the asymptotic behavior of the von By setting k1 = k2 = z1 = 0, we restrict to the SU(2) Neumann entropy Sn = Trn(ρn log ρn) of the reduced − 2 sector of SU(3). This has been used as a check to the density matrix ρn = Tri>nρ of a subchain with n spins program; the DMRG in this case reproduces perfectly all of a critical system of length L, as a function of n and L, the energy states of the SU(2) Heisenberg model. where ρ is the density matrix associated to the ground We study now the isotropic AFM chain with periodic state of the chain. Indeed, one has [19, 20]: boundary conditions by means of an infinite size DMRG with up to m = 2200 states in order to reduce the uncer- c L π Sn = log2 sin n + A. (12) tainty on the energies to the order of magnitude of the 3 π L truncation error. The data for the ground state and the first excited states are plotted in Fig. 1. As usual c is the central charge while A is a non-universal Let us first concentrate on the ground state, which, in constant. The DMRG computes the density matrix for agreement with theoretical predictions, it is found only a block of length n in a chain of length L, so that Sn 2 when L = 6M. The plot of E00 as a function of 1/L becomes quite simple to calculate. Fig. 2 shows the be- shows a good linear behavior; this justifies the fitting of havior of the von Neumann entropy as a function of the our data by the CFT equations for the GS: block length n in (a) and as a function of the quantity L y = log2 π /3 (obtained from (12) by setting n = L/2) E00 πcv = e∞ , (10) in (b), for values of the DMRG states equal to 1000 and L − 6L2 2000. The data confirm the linear behavior expected where e∞ and the product cv are kept as fitting param- from Eq. (12). The linear regression Sn = cy + A yields the value for the constant A = (1.774 0.002) and for the eters. We obtain: e∞ = 0.518288 and cv = 2.04419. − central charge c = (1.995 0.001). Thus± the theoretical In order to derive the value of v we need an independent ± derivation of c. The central charge for a SU(N) level-k prediction of Eq. (11) is confirmed with very high accu- WZNW model is given by [16] racy. From Fig. 2(b) it is also evident that the values obtained when keeping only 1000 states in the DMRG k (N 2 1) run are much less precise. This is the reason why we c = − . (11) k + N have then performed all calculations while keeping 2000 4

(a) L=18 (b) L x1 x2 L=24 5 ± ± L=30 6M+2 0.3414 0.0001 0.6291 0.0003 4.2 6M+4 0.3406 ± 0.0002 0.6238 ± 0.0003 L=36 4.5 L=42 (N) 0.3410 ± 0.0002 0.6265 ± 0.0003 4 4 6M - 0.6503 ± 0.0003 n L/2

S 3.8 3.5

S Table II: Results of the numerical analysis. The velocity and ± 3 the central charge are respectively: v = (1.0247 0.0005) and 3.6 c = (1.995 ± 0.001). The scaling dimensions of the ﬁrst (x1) 2.5 2000 states and the second (x2) excited states for each L are calculated 3.4 1000 states as described in the text. The mean value (N) between L = 2 0 10 20 30 40 0 0.5 1 1.5 6M +2 and L = 6M +4 for x1 and x2 is provided (see also Fig. n π log2(L/ )/3 1). For L = 6M only the ﬁrst excitation above the ground state has been considered.

Figure 2: Analysis of the von Neumann entropy. The figure on the left (a) shows Sn for m = 1000 DMRG states as a where xj is the scaling dimension of the j th excitation function of the block length n for different chains of length L; − the figure on the right (b) is the plot of the half chain entropy for a given chain of length L; E00 is given by Eq. (10) (n = L/2) for m = 1000 and m = 2000 DMRG states. The where c and v have been derived before and are reported linear fit on these data provides the values for the central in the caption of Table II. The numerical coefficients charge c and the constant A. for the scaling dimensions that one can obtain from the DMRG data of Fig. 1 are listed in Table II. As expected, the values of the allowed conformal di- states. Finally, the value of c can be substituted into the mensions are very close to the values of 1/3 and 2/3 pre- product cv derived from the GS to recover the velocity of dicted by a SU(3)1 WZNW model. the excited modes: v = (1, 0247 0.0005), which is close to the expected value [8] π/3. ± Before proceeding with the analysis of the excited IV. CONCLUSIONS states, let us check the asymptotic value of the energy density e∞. The theoretical prediction for the ground state of the S = 1 bilinear biquadratic Heisenberg Hamil- We have provided strong numerical evidence of the tonian (see Ref. [21]) is criticality of the AFM SU(3) spin chain. Also, we have confirmed that the conformal field theory describing the π EGS = ln 3 +1, (13) chain is effectively the SU(3)1 WZNW model, by compu- − − 3√3 tating the central charge and scaling dimensions of the which already takes into account the factor -1 of the lowest excited states of the model, which turn out to be l.h.s. of equation (5). Starting from the correspondence organized according to the IR of SU(3)1 Kac-Moody al- between our SU(3) chain and the biquadratic one (5), gebra. we can compare the value of e∞ we obtained with the There are many interesting generalizations of the above one predicted by equation (13): EGS = 0.703212. The models which deserve further study. In particular, a sim- match is exact to the third decimal ( 0−.703243), if one ilar ferromagnetic spin chain is connected with the non- 2 also recalls that the Hamiltonian has a− factor 1/4 (due to linear CP sigma mode at θ = π and might be useful the definition of the spin variables in terms of the SU(3) to clarify some controversial problems of the model. An- generators) so that e∞ needs to be multiplied by a factor other interesting problem is to consider larger SU(N) two, and summed to the factor 1/3 of equation (5). This symmetry groups. In two-dimensional chains, the vac- is a further numerical proof of the equivalence between uum state is of Néel type for N 4 and of Spin-Peierls type for N 5 [22]. The analysis≤ by means of DMRG the the Lai-Sutherland and the SU(3) spin model. ≥ Let us study now the excited states. From Fig. 1, one technique might shed some light on the transition mech- immediately sees that the slope of excited states depends anism. on the length of the chain. In particular, for L =6M the first excitation scales with a slope which is unmistakably different from the slope of the L =6M +2 or L =6M +4 Acknowledgments data. For small values of L the data corresponding to the same Sz but with opposite Qz are split by a finite size We would like to thank G. Morandi, C. Degli Esposti correction, while for increasing values of L they tend to Boschi, M. Roncaglia and L. Campos Venuti for inter- overlap and scale to the same asymptotic value. esting and helpful discussions. One of the authors S.P. For the excited states CFT predicts that: would also like to thank S. Rachel, R. Thomale and A. 2πv Läuchli for very constructive discussions. The work of Ej E = xj (14) − 00 L M.A. was partially supported by a cooperation grant 5

INFN-CICYT, the Spanish CICYT grant FPA2006-2315 and DGIID-DGA (grant2006-E24/2).

[1] Y. Tokura, N. Nagaosa, Science 288 (2000) 462 [13] C. Itoi and M. Kato, Phys. Rev. B 55 (1997) 8295 [2] M. Greiter, S. Rachel and D. Schuricht, Phys. Rev. B 75, [14] M. Fueheringer, S. Rachel, R. Thomale, M. Greiter and (2007) 060401; D. Schuricht and M. Greiter, Europhys. P. Schmitteckert, arXiv:08062563 Lett. 71 (2005) 987 [15] A. Schmitt, K-H Mütter and M. Karbach, J. Phys A 29 [3] S. Chen, C. Wu, S. C. Zhang and Y. Wang, Phys. Rev. (1996) 3951 B 72 (2005) 214428 [16] P. Di Francesco, P. Mathieu, D. Sénéchal, Conformal [4] S. Chen, Y. Wang, W. Q. Ning, C. Wu and H. Q. Lin, Field Theory, Springer (1997) Phys. Rev. B 74 (2006) 174424 [17] S. Chaturvedi, N. Mukunda, J. Math. Phys. 43 (2002) [5] R. Citro, E. Orignac, N. Andrei, C. Itoi, S. Qin, J. Phys. 5262 : Condens. Matter 12 (2000) 3041 [18] T. Hakobyan, Nucl. Phys. B 699 (2004) 575 [6] J. A. Minahan and K. Zarembo, JHEP 0303 (2003) 013 [19] C. Holzhey, F. Larsen and F. Wilczek, Nucl. Phys. B 424 [7] N Beisert, M Staudacher, Nucl. Phys. B 670(2003) 439 (1994) 443 [8] B. Sutherland, Phys. Rev. B 12 (1975) 3795 [20] P. Calabrese and J. Cardy, J. Stat. Mech. 0406 (2004) [9] J. K. Lai, J. Math. Phys. 15 (1974) 1675 002 [10] K. Chang, I. Affleck, G. W. Hayden and Z. G. Soos, J. [21] G. V. Uimin, JETP Lett 12 (1970) 225 Phys. : Condens. Matter 1, (1989) 153 [22] K. Harada, N. Kawashima and M. Troyer, Phys. Rev. [11] I. Affleck, Nucl. Phys. B 265 (1986) 409 Lett. 90 (2003) 117203 [12] I. Affleck, Nucl. Phys. B 305 (1988) 582