arxiv:cond-mat/ v1 [cond-mat.str-el] 25 Aug 2000

Transcription

1 LA-UR-00-XXXX February, 008 submitted to Physical Review Letters June st Generalized Jordan-Wigner Transformations C.D. Batista and G. Ortiz Theoretical Division, Los Alamos National Laboratory, Los Alamos, NM (Received February, 008) arxiv:cond-mat/ v cond-mat.str-el] 5 Aug 000 We introduce a new spin-fermion mapping, for arbitrary spin generating the U() group algebra, that constitutes a natural generalization of the Jordan-Wigner transformation for =. The mapping, valid for regular lattices in any spatial dimension d, serves to unravel hidden symmetries in one representation that are manifest in the other. We illustrate the power of the transformation by finding exact solutions to lattice models previously unsolved by standard techniques. We also present a proof of the existence of the Haldane gap in = bilinear nearest-neighbors Heisenberg spin chains and discuss the relevance of the mapping to models of strongly correlated electrons. Moreover, we present a general spin-anyon mapping for the case d. Introduction. Theories of magnetism in normal matter are direct manifestations of quantum mechanics. The vast maority of phenomena occurring in magnets, such as different types of magnetically ordered states, are often simply described by interacting quantum spins ]. Although quantum spins do not behave as either pure boson or fermion operators, different representations have invoked such particle statistics, for example, the Holstein- Primakoff, chwinger, and Dyson-Maleev boson representations for arbitrary spin, or the Jordan-Wigner and Maorana fermion representations for spin = magnets. ometimes it is extremely convenient to reformulate a difficult strongly correlated problem in a way that it becomes more manageable; in some cases there is an exact dualism. This is the idea behind the bosonization techniques and the different algebra representations of a physical problem. The point is that these different representations help us understand various aspects of the same problem by transforming intricate interaction terms into simpler ones. Often, fundamental symmetries which are hidden in one representation are manifest in the other and, moreover, problems which seem untractable can even be exactly solved after the mapping. The simplest and perhaps most popular example is the equivalence between the Heisenberg-Ising = XXZ chain and a model of interacting spinless fermions through the Jordan-Wigner transformation ]. The Jordan-Wigner transformation involves the = irreducible representation of the Lie group U(). Here we generalize this spin-fermion mapping to any irreducible representation of dimension. The three generators µ (µ = x, y, z) of the Lie algebra for each lattice site satisfy the commutation relations 3], ] = z, z, ] ± = ± ±, {, } ( = ( ) ( z ) ). () We start by analyzing the one-dimensional = case. Then, we will show a generalization to arbitrary spin and spatial dimension d. = Mapping. We introduce the following composite operators f = c c, f = c c, (3) written in terms of the Hubbard operators c σ = c σ ( n σ ) and c σ = ( n σ )c σ (σ =, ), which form a subalgebra of the so-called double graded algebra pl(, ) 4]. A bar in a subindex means the negative of that number (e.g., σ = σ).] For spins on a lattice we fermionize the spins and reproduce the correct spin algebra with the following transformation = ( c K K c ), = (K c c K ), z = n n, whose inverse manifests the nonlocal character of the mapping f = expiπ k<( z k ) ], f = exp iπ k<( z k) ], c = z f, c = f z, c = z f, c = f z, µ, ν k] = iδk ǫ µνλ λ, () with ǫ the totally antisymmetric tensor. The algebra generated by the (linear and Hermitian operators) µ is the enveloping algebra of the group U(). Equivalently, in terms of the ladder operators ± = x ± iy where the string operators K = expiπ k< n k] = k< σ ( n σ), and the number operators n k = n k n k. These f-operators have the remarkable property that { } f, f = {, }, (4)

2 LA-UR-00-XXXX February, 008 submitted to Physical Review Letters June st which suggests an analogy between spin operators and constrained fermions. Half-odd integer spin chains have a qualitatively different excitation spectrum than integer spin chains. The Lieb, chultz, Mattis and Affleck theorem 5] establishes that the half-odd integer antiferromagnetic (AF) bilinear nearest-neighbors (NN) Heisenberg chain is gapless. The same model with integer spins is conectured to display a Haldane gap 6]. To understand the origin of the Haldane gap we analyze the form of the d = XXZ Hamiltonian using the above representation (an overall omitted constant J > 0 determines the energy scale) H xxz = z z ( x x ) y y = H z H xx. (5) It is easy to show that the constrained fermion version of this Hamiltonian is a ( = ) t-j z model 7] plus a particle non-conserving term which breaks the U() symmetry H xxz = σ ( n n )( n n ) ( c σ c σ c σ c σ H.c. ). (6) In the isotropic = limit, H xxz can be written in a compact way H Heisenberg = (Ψ ) Ψ (Ψ ) Ψ, (7) where is an irreducible matrix representation of = (3 3 matrices) while Ψ is the ( 3) vector ] Ψ n =, ( c c )K, n. (8) The charge spectrum of the ( = ) t-j z model is gapless but the spin spectrum is gapped due to the explicitly broken U() symmetry (Luther-Emery liquid) 7]. Therefore, the spectrum of the = Hamiltonian associated to the t-j z model (which has only spin excitations) is gapless. Hence the term which explicitly breaks U() must be responsible for the opening of the Haldane gap. We can prove this by considering the perturbative effect that the interaction η σ ( c σ c σ H.c.) has on the t-j z Hamiltonian. To linear order in η (> 0), Eq. 7 maps onto the ( = ) XYZ model with J x = η, J y = η, and J z =. From exact solution of this model 8], it is seen that the system is critical only when η = 0 while for η 0 a gap to all excitations opens. = Integrable Models 9]. To illustrate further the power of our spin-fermion mapping we now present exact solutions of d = models that have not been discovered by traditional techniques. These models correspond to the family of bilinear-biquadratic Hamiltonians H ( ) = = H z Hxx H z σ { H z, } Hxx c σ c σ, (9) that can be mapped onto a ( = ) t-j z Hamiltonian, whose quantum phase diagram has recently been exactly solved 7]. Another well-studied class of bilinear-biquadratic U() invariant Hamiltonians is 0] H ( ) = ( ), (0) for -. The pure Heisenberg ( =0) and Valence Bond olid models ( = 3 ) belong to the Haldane gapped phase, which extends over the whole interval except at the boundaries = ± that are quantum critical points. The case =- is known to be Bethe ansatz soluble with a unique ground state and gapless. For = we can map H () onto the supersymmetric ( = ) t-j Hamiltonian plus a NN repulsive interaction H () = σ ( c σ c σ H.c. ) s s ( n 3 4 n n ), () where s represents a = operator. This model is Bethe-ansatz soluble with a gapless phase ] and is known as the Lai-utherland solution ]. We now discuss the importance of our generalized Jordan-Wigner transformation in unraveling hidden symmetries of an arbitrary spin Hamiltonian. In Eq. 0, for example, the = U() symmetry is manifest. However, both the = U() and global U() gauge symmetries are hidden. On the other hand, in the transformed Hamiltonian, Eq., these two symmetries are manifested explicitly through rotational invariance and charge conservation. The generators of these symmetries are related through the mapping already introduced. To illustrate this, we consider the U() symmetry case. Here the generator of the transformation is Q = n which maps onto Q = (z ) in the spin representation. The total group symmetry of the Hamiltonian is U(3). Generalized Transformation. A general transformation for arbitrary spin and spatial dimension is the following Half-odd integer spin (σ F = {,...,}):

3 LA-UR-00-XXXX February, 008 submitted to Physical Review Letters June st = η c K η σ c σ c σ, σ F σ = η K c η σ c σ c σ, z = σ F σ σ F ( σ) n σ, c σ = K L σ ( ) σ P, where P = τ z τ, L τ F σ σ = ητ. τ= Integer spin (σ F = {,...,,,..., }): = η 0 ( c K K c ) η σ c σ c σ, The string operators K introduce nonlinear and nonlocal interactions between the constrained fermions. For d lattices (K = K, K i, K ] = 0) they are the so-called kink operators K = expiπ n k ], while for d 3] k< K = expi k n k = a(k,) n k ], with (4) σ F α n kσ = P α k. (5) Here, a(k,) is the angle between the spatial vector k and a fixed direction on the lattice, and a(,) is defined to be zero. We comment that the d kink operators constitute a particular case of Eq. 4 with a(k, ) = π when k < and equals zero otherwise. For d >, the string operators generalize 9] along the lines introduced in Ref. 4]. σ F σ, = η 0 (K c c K ) z = σ F σ, σ n σ, σ F {( ) σ P if σ > 0, η σ c σ c σ, Ð ÚÓÖ Þ ¼ Ð ÚÓÖ c σ = K L σ where ( ) σ P if σ < 0, P = τ z τ τ F σ, L σ = ητ, τ=0 and η σ = ( σ)( σ ) (see Fig. ). The total number of flavors is N f =, and the = case simply reduces to the traditional Jordan-Wigner transformation. These mappings enforce the condition on the Casimir operator = ( ). The generalized constrained fields c σ = c σ ( n τ ), c σ = ( n τ ) c σ τ F α τ F α () form a subalgebra of the generalized Hubbard double graded algebra, where the unconstrained operators c σ, c σ satisfy the standard fermion anticommutation relations (α =, depending upon the spin character of the representation). These generalized constrained operators (only single occupancy is allowed) anticommute for different sites } { cσ, c kσ } = { c σ, c kσ = 0, } ( n τ ) if σ = σ { c, σ, c kσ = δ τ F k α (3) τ σ c σ c σ if σ σ, and their number operators satisfy n σ n σ = δ σσ n σ. ÁÒØ Ö À Ð ¹Ó ÁÒØ Ö FIG.. Constrained fermion states per site for integer and half-odd integer spin. In both cases there are flavors and the corresponding values of z are shown in the middle column. One degree of freedom is assigned to the fermion vacuum (circle) whose relative position depends upon the spin character. There is always the freedom to perform rotations in spin space to get equivalent representations to the one presented above. However, for bilinear isotropic NN Heisenberg (spin U() rotationally invariant) Hamiltonians in the large- limit there is a fundamental difference between effective integer and half-odd integer spin cases. In the latter case a new local U() gauge symmetry emerges that is explicitly broken in the integer case. For d lattices, this is precisely what distinguishes Haldane gap systems 6] from half-odd integer spin chains that are critical. We mention that other fermionic representations are feasible. In particular, for half-odd integer cases where = N f ( Nf ) i=0 i = Nf (e.g., = 3 with N f = ) a simple transformation in terms of standard unconstrained fermions is possible 9]. 3

4 LA-UR-00-XXXX February, 008 submitted to Physical Review Letters June st d Lattices and pin-anyon Mapping. The generalization of these transformations to higher dimensions gives new exact mappings between spin theories and constrained fermion systems in the presence of gauge fields. To illustrate this we write the = Hamiltonian H () in the fermion representation for d = H () = ( ) c e νσ eiaν() c σ H.c. s s eν σ,ν,ν and ( n 3 4 n n eν ), (6) A ν () = k a(k,) a(k, e ν )] n k, (7) where e ν (ν =, ) are basis vectors of the Bravais lattice connecting NN and s represent sites of the corresponding d lattice. We note that the field A ν () is associated with the change in particle statistics. It is well-known 3,4] that the same transmutation of particle statistics can be achieved via a path-integral formulation for H () where an Abelian lattice Chern-imons term is included. In this formulation a constraint (Gauss s law) requiring that the gauge flux through a plaquette be proportional to the total fermion density on the site, n, is enforced. This suggests that our spin-fermion mapping can be generalized to an spin-anyon transformation with a hardcore condition for the anyon fields 9]. In fact, one can formally take our generalized Jordan-Wigner transformation and replace the string operators K by the statistical operators K (θ) = expiθ k a(k,) n k] with 0 θ. With this choice, the c operators satisfy equal-time anyon commutation relations θ = (0) corresponds to constrained fermions(bosons)] 9]. imilar ideas apply for d lattices. One immediately sees the relevance of these transformations for the theories of magnetism and hightemperature superconductivity: A class of = Hamiltonians that can be mapped onto a lattice-gauge (Chern- imons) = t-j theory and vice versa, for example, a = t-j model, H t J = t ( c σ c ) eνσ H.c. J s s eν σ,ν,ν µ n, (8) can be exactly mapped onto a lattice-gauge bilinearbiquadratic = theory H t J = µ ( z ) J Hν z 4t 8 J eiaν() e ν,ν 4t } ] {H z J ν, eiaν() eν ( e ν ) H.c.. (9) By means of a semiclassical approximation it has been shown 5] that the ground state of H () is on the boundary between AF ( < ) and orthogonal nematic (non-uniform, > ) phases 5,0]. These two states are the result of the competition between the quadratic and quartic spin-exchange interactions. In terms of the equivalent t-j gauge theory this translates into a competition between antiferromagnetism and delocalization. Qualitatively, the string-path of the particle moving in an AF background gives rise to a linear confining potential since the number of frustrated magnetic links is proportional to the length of the path. This observation suggests that the inhomogeneous phases observed in the striped high-t c compounds can be driven by the competition between magnetism and delocalization. ummary. We introduced a general spin-fermion mapping for arbitrary spin and spatial dimension that naturally generalizes the Jordan-Wigner transformation for =. Mathematically, we established a one-to-one mapping of elements of a Lie algebra onto elements of a fermionic algebra with a hard-core constraint. everal generalizations, like a spin-anyon mapping, and important consequences result from these transformations 9]. Incidentally, we note that there are extremely powerful numerical techniques (cluster algorithms 6]) to study quantum spin systems, and our mapping allows one to extend these methods to study the equivalent fermionic problems. We thank J.E. Gubernatis for a careful reading of the manuscript. This work was sponsored by the U DOE under contract W-7405-ENG-36. ] D.C. Mattis, The Theory of Magnetism I (pringer, Berlin, 988). ] P. Jordan and E. Wigner, Z. Phys. 47, 63 (98). 3] For a d-dimensional lattice with N sites, the operator µ is defined in terms of a Kronecker product as µ = l l µ }{{} th l, where l is the ( ) ( ) unit matrix and µ is a spin- operator. Thus µ admits a matrix representation of dimension ( ) N ( ) N. 4] A.M. Tsvelik, Quantum Field Theory in Condensed Matter Physics (Cambridge, Cambridge, 995). 5] E.H. Lieb, T.D. chultz, and D.C. Mattis, Ann. Phys. 6, 407 (96); I. Affleck and E.H. Lieb, Lett. Math. Phys., 57 (986). 6] F.D.M. Haldane, Phys. Lett. A 93, 464 (983). 7] C.D. Batista and G. Ortiz, cond-mat/ ] R.J. Baxter, Exactly olved Models in tatistical Mechanics (Academic Press, London, 990). 4

5 LA-UR-00-XXXX February, 008 submitted to Physical Review Letters June st 9] More examples and details will be published elsewhere. 0] G. Fáth and J. ólyom, Phys. Rev. B 5, 360 (995). ] P. chlottmann, Phys. Rev. B 36, 577 (987). ] C.K. Lai, J. Math. Phys. 5, 675 (974); B. utherland, Phys. Rev. B, 3795 (975). 3] E. Fradkin, Phys. Rev. Lett. 63, 3 (989); D. Eliezer and G.W. emenoff, Phys. Lett. B 86, 8 (99). 4] L. Huerta and J. Zanelli, Phys. Rev. Lett. 7, 36 (993). 5] N. Papanicolaou, Nucl. Phys. B 305, 367 (988). 6] J.E. Gubernatis and N. Kawashima, in Monte Carlo and Molecular Dynamics of Condensed Matter ystems (ocietà Italiana de Fisica, Bologna, 995) edited by K. Binder and G. Ciccotti. 5