Exact Ground State and Elementary Excitations of the Spin Tetrahedron Chain
Total Page:16
File Type:pdf, Size:1020Kb
PHYSICAL REVIEW B 74, 174424 ͑2006͒ Exact ground state and elementary excitations of the spin tetrahedron chain Shu Chen,1,2 Yupeng Wang,2 W. Q. Ning,1,3 Congjun Wu,4 and H. Q. Lin1 1Department of Physics and Institute of Theoretical Physics, The Chinese University of Hongkong, Hongkong, People’s Republic of China 2Beijing National Laboratory for Condensed Matter Physics, Institute of Physics, Chinese Academy of Sciences, Beijing 100080, People’s Republic of China 3Department of Physics, Fudan University, Shanghai 200433, People’s Republic of China 4Kavli Institute for Theoretical Physics, University of California, Santa Barbara, California 93106, USA ͑Received 14 February 2006; revised manuscript received 24 August 2006; published 21 November 2006͒ We study the antiferromagnetic spin exchange models with S=1/2 and S=1 on a one-dimensional tetrahe- dron chain by both analytical and numerical approaches. The system is shown to be effectively mapped to a decoupled spin chain in the regime of strong rung coupling, and a spin sawtooth lattice in the regime of weak rung coupling with spin 2S on the top row and spin S on the lower row. The ground state for the homogeneous tetrahedron chain is found to fall into the regime of strong rung coupling. As a result, the elementary excitation for the spin-1/2 system is gapless whereas the excitation for the spin-1 system has a finite spin gap. With the aid of the exact diagonalization method, we determine the phase diagram numerically and find the existence of an additional phase in the intermediate regime. This phase is doubly degenerate and is characterized by an alternating distribution of rung singlet and rung spin 2S. We also show that the SU͑3͒ exchange model on the same lattice has completely different kind of ground state from that of its SU͑2͒ correspondence and calculate its ground state and elementary excitation analytically. DOI: 10.1103/PhysRevB.74.174424 PACS number͑s͒: 75.10.Jm ͑ ͒ I. INTRODUCTION 1D tetrahedron chain as shown in Fig. 1, with site spins Si The study of quantum antiferromagnetic spin models with residing in four corners of each tetrahedron. The basic unit, strong frustration has attracted great attention over the past i.e., a tetrahedron is composed of four spins with equal anti- decades. Early investigations of the frustrated quantum mag- ferromagnetic exchanges between each pair of spins. The nets were partly motivated by the work of Anderson to tetrahedron chain can be also viewed as a 1D pyrochlore search the resonating-valence-bond ͑RVB͒ ground state in strip in a 2D pyrochlore lattice, in which only two of four such systems.1 In the frustrated quantum magnets, the mag- corners of each tetrahedron are shared by a neighboring tet- netic ordering is generally suppressed by the frustration. rahedron. Generally, a three-dimensional ͑3D͒ pyrochlore Some well-studied frustrated magnetic systems include, for lattice is a network of corner-sharing tetrahedra and a 2D example, the kagome lattice and pyrochlore lattice, in which pyrochlore model, named also as a checkerboard-lattice the interplay between frustration and quantum fluctuation model, is obtained by a projection of the 3D lattice on a leads to rich varieties of phenomena. Recently, the frustrated plane. Different from the 1D pyrochlore strip considered in magnets are believed to be prominent candidates of realizing this article, for both the 3D and 2D pyrochlore lattices, each spin liquid states with exotic ground state and deconfined corner of the tetrahedron is shared by a neighboring tetrahe- fractional excitations.2–4 While the mechanism of deconfin- dron. As one of the most frustrated antiferromagnets, the ment in two-dimensional ͑2D͒ magnetic systems is less clear, model of spin pyrochlore lattice has been investigated by a the deconfinement of spinons in quasi-one-dimensional mag- variety of techniques including the semiclassical large-S netic system, which is closely related to the phenomenon of spin charge separation, is well investigated.5–8 Since the 1980’s, a large number of low-dimensional frustrated mag- nets have been synthesized experimentally.9 Generally speaking, strongly geometrical frustration in these systems allows the simple dimerized state to be the ground state of the low-dimensional frustrated spin system and opens a spin gap. So far, a variety of quasi-one-dimensional frustrated models have been studied theoretically.10–12 Additionally, im- portant progress has been made in trapping cold atoms under a highly controllable way very recently, and thus it stimulates intensive investigation on how to simulate the magnetic sys- tems using cold atoms. A number of schemes have been pro- posed to implement a variety of quantum spin models in optical lattices.13 FIG. 1. ͑Color online͒ Schematic pictures of ͑a͒ the homoge- In this paper, we investigate both the ground- and excited- neous spin tetrahedron chain, ͑b͒ the spin tetrahedron chain corre- state properties of the spin models on a one-dimensional sponding to model ͑2͒. 1098-0121/2006/74͑17͒/174424͑8͒ 174424-1 ©2006 The American Physical Society CHEN et al. PHYSICAL REVIEW B 74, 174424 ͑2006͒ limit, large-N expansion of the SU͑N͒ model, the contractor For the tetrahedron chain model ͑2͒, it is convenient to renormalization method based on the cluster expansion, and reformulate the Hamiltonian as follows: the bosonization method on the anisotropic limit.14–17 In ˆ ˆ ͑ˆ ˆ ͒ ˆ spite of the intensive research, even the ground state proper- HA = J͚ Si,1Si+1,1 + J͚ Si,1 + Si+1,1 · Ti ties of the 3D pyrochlore lattice are not well understood. For i i the 2D pyrochlore lattice, the numerical results based on ex- JЌ 2 + ͚ Tˆ − NJЌS͑S +1͒, ͑3͒ act diagonizations have shown that the ground state has 2 i plaquette order.18 However, for the 1D pyrochlore strip, we i can determine its ground state and elementary excitation in ˆ 2 ͑ ͒ where T =Ti Ti +1 with Ti =0,...,2S. In the strong cou- an exact manner. With the aid of numerical diagonalization i pling limit JЌ →ϱ, a pair of spins on each rung would form of the corresponding spin lattice systems with small sizes, ͑ ͒ a singlet spin dimer with Ti =0. This implies that the spins we also investigate the quantum phase transitions of the in the horizontal direction along the chain is effectively de- ground state due to the change of exchange strengths along coupled with the spins on the vertical rungs. Therefore, the perpendicular rungs. ground state of HA is a product of the ground state of spin The spin model on a 1D pyrochlore strip as shown in Fig. chain and rung singlets. Explicitly, it is represented as fol- 1 is described by the Hamiltonian: lows: ˆ ˆ ͑ ͒ ͉ ͘ ͉ ͘ ͓ ͔ ͑ ͒ H = J͚ Si · Sj, 1 GS = BA ͟ Si,2,Si,3 , 4 ͗ij͘ i where ͗ij͘ denotes sum over all the nearest neighbors along where ͉BA͘ denotes the Bethe-ansatz ground state wave ˆ ͚ ˆ ˆ the tetrahedron chain and Si represents the spin operator re- functions of the 1D Heisenberg chain Hchain=J iSi,1Si+1,1 ͓ ͔ ͓͑↑͔ ͓↓͔ ͓↓͔ ͓↑͔ ͒ ͱ siding in site i. In this work, we study both the spin-1/2 and and Si,2,Si,3 = i,2 i,3− i,2 i,3 / 2 is the dimer sin- spin-1 models on the pyrochlore strip. Our results show that glet across the ith vertical rung. The corresponding ground the excitation spectrum for the spin-1/2 system is gapless state energy is and the elementary excitation of the spin-1 system has a finite spin gap. This model can be extended to the cases A BA 3 E = E ͑N͒ − NJЌ, ͑5͒ where the strengths of bonds among the tetrahedra are not g g 4 homogeneous. Here we only consider the inhomogeneous BA͑ ͒ where Eg N is the ground state energy of the N-site case as shown in Fig. 1͑b͒, where we use JЌ to represent the Heisenberg spin-1/2 chain. From the Bethe-ansatz solution exchange strength along the vertical rung of each tetrahe- of Heisenberg model, we know the exact ground state energy dron. For convenience, we rewrite the Hamiltonian of the EBA͑N͒ N J spin tetrahedron chain corresponding to Fig. 1͑b͒ as follows: g / =−0.4431 at the infinite length limit. In fact, utilizing the Raleigh-Ritz variational ͓ˆ ͑ˆ ˆ ͒ ͑ˆ ˆ ˆ ͒ˆ ͔ principle,11,19–21 we can exactly prove the state given by Eq. HA = J͚ Si,1 Si,2 + Si,3 + Si,1 + Si,2 + Si,3 Si+1,1 i ͑4͒ is the ground state of Hamiltonian ͑2͒ as long as JЌ ജ2J. To see it clearly, we can rewrite the Hamiltonian ͑2͒ ˆ ˆ ͑ ͒ + JЌ͚ Si,2 · Si,3. 2 with JЌ =2J as the sum of a Heisenberg chain and 2N pro- i jection operators, which reads ͑ ͒ ͑ ͒ It is obvious that model 2 reduces to model 1 when JЌ 3J =J, i.e., the homogeneous tetrahedron chain ͑1͒ is a special H = J͚ Sˆ Sˆ + ͚ P3/2͑Sˆ ,Sˆ ,Sˆ ͒ A i,1 i+1,1 2 i,1 i,2 i,3 case of model ͑2͒. Since the model can be represented as a i i sum of local Hamiltonian on each tetrahedron, a classical 3J ground state is obtained whenever the total spin in the tetra- + ͚ ͓P3/2͑Sˆ ,Sˆ ,Sˆ ͒ −1͔, ͑6͒ 2 i,2 i,3 i+1,1 hedron is zero for a homogeneous model.