Charge dynamics in the Mott insulating phase of the ionic Hubbard model

  • 格式:pdf
  • 大小:114.56 KB
  • 文档页数:4

a r X i v :c o n d -m a t /0308120v 2 [c o n d -m a t .s t r -e l ] 18 S e p 2003Charge dynamics in the Mott insulating phase of the ionic Hubbard modelA.A.AligiaCentro At´o mico Bariloche and Instituto Balseiro,Comisi´o n Nacional de Energ´ıa At´o mica,8400Bariloche,ArgentinaI extend to charge and bond operators the transformation that maps the ionic Hubbard model at half filling onto an effective spin ing the transformed operators I calculate the amplitude of the charge density wave in different dimensions D.In 1D,the charge-charge correlations at large distance d decay as d −3ln −3/2d ,in spite of the finite charge gap,due to remaining charge-spin coupling.Bond-bond correlations decay as (−1)d d −1ln −3/2d as in the usual Hubbard model.I.INTRODUCTIONThe ionic Hubbard model (IHM)has been proposed in the 80’s for the description of the neutral-ionic tran-sition in mixed-stack charge-transfer organic crystals.1,2In the 90’s the interest on the model increased due to its potential application to ferroelectric perovskites.3–9The model in any bipartite lattice can be written as:H =H 0+H t ;H 0=∆II.THE CANONICAL TRANSFORMATIONIn this Section,I use the canonical transformation plus projection onto the low-energy subspace that maps Hinto a spin Hamiltonian˜H,to transform charge and bond operators in one direction(δ=±1)n i=n i↑+n i↓,b i= δσδ(c†i+δσc iσ+H.c.),(2)and discuss some symmetry properties.For our purposes, it is enough to work up to second order in t/(U−∆). Thus˜H=P e−S He S P=P(H+[H,S]+1E m−E n,(4) where E j is the energy of the eigenstate|j of H0.From Eqs.(3)and(4),proceeding in a similar way as done below,one obtains the known result:2˜H=4t2U4).(5)The most important correction to˜H in higher order is a next-nearest-neighbor exchange of order t4,which does not affect the physics for t≪U−∆.2,14The effective Hamiltonian˜H to all orders in t is in-variant under a nearest-neighbor translation Tδ,while H is not.2This is simply a consequence of the fact that to all orders˜H is a purely spin Hamiltonian,since the charge degrees of freedom are frozen(n i≡1for all i). Then˜H is invariant under the electron-hole transforma-tion T eh:c†iσ→σc i¯σ,which leaves invariant all spin operators.Since the original Hamiltonian H is invariant under the product TδT eh,then˜H should also be invari-ant under Tδ.However other transformed operators,like ˜n i below,are not invariant under Tδ.The transformed charge operator is:˜n i=P e−S n i e S P∼=P(n i+[n i,S]+1U−(−1)i∆2c†iσ′c i+δσ′c†i+δσc iσ+ t2−2S i·S j.(8) And replacing this in Eq.(7):˜n i=1−(−1)i2U∆t2U+(−1)i∆c†i+δσc iσc†iσ′c i+δσ′+tU2−∆2δδS i·S i+δ.(10)The second member is a local measure of the asymmetry between the spin“bonds”involving site i in the directionofδ.Here there is no essential difference with the resultfor the ordinary Hubbard model.III.OBSER V ABLES AND CORRELATIONFUNCTIONSFor t≪U−∆,using Eq.(9)the amplitude of the charge density wave is given by:A=| n i−n i+δ H|=| ˜n i−˜n i+δ ˜H|=2a δ(1−4 S i·S i+δ ˜H),(11)where the subscript in the expectation values indicates the Hamiltonian with which they are calculated anda=2U∆t22ϕ(x))for the slowest decaying part of S z l,the operator product expansion S z l S z l+1≈cos(2√8 (S+0S−δ+S−0S+δ)(S+d S−d+δ+S−d S+d+δ)− S z0S zδS z d S z d+δ .(17)Thefirst term in this equation turns out to be the dom-inant one,and therefore I explain it in more detail.Per-forming a Jordan-Wigner transformation from spin op-erators to fermions with annihilation operators a j,goingto the continuum limit using a j=i j L(x)+(−i)j R(x),with x=ja,and then bosonizing one gets:δ(S+j S−j+δ+H.c.)→ δ((a†j a j+δ+H.c.)→ δ[iδL†(x)L(x+aδ)+(−i)δR†(x)R(x+aδ)+(−1)j(−i)δL†(x)R(x+aδ)+(−1)j iδR†(x)L(x+aδ)+H.c.]→2i[L†∂L∂x−(−1)j L†∂R∂x]+H.c.=2i[2L†∂L∂x+(−1)j∂∂x)2+(∂θ∂xcos(√2ϕ(x))::cos(√|x−y|ln3/2|x−y|,(19)together with Eqs.(14)to(19),Ifinally obtain(exceptfor some factor of the order of one):C d≈48U2∆2t42ϕ(x)),oneobtains:b i b i+d ≈64U2t21J.Hubbard and J.B.Torrance,Phys.Rev.Lett.47,1750(1981).2N.Nagaosa and J.Takimoto,J.Phys.Soc.Jpn.55,2735(1986).3T.Egami,S.Ishihara,and M.Tachiki,Science261,1307(1993);Phys.Rev.B49,8944(1994).4R.Resta and S.Sorella,Phys.Rev.Lett.74,4738(1995).5G.Ortiz,P.Ordej´o n,R.M.Martin,and G.Chiappe,Phys.Rev.B54,13515(1996);references therein.6R.Resta and S.Sorella,Phys.Rev.Lett.82,370(1999).7N.Gidopoulos,S.Sorella,and E.Tosatti,Eur.Phys.J.B14,217(2000).8M.Fabrizio,A.O.Gogolin,and A.A.Nersesyan,Phys.Rev.Lett83,2014(1999).9M.E.Torio,A.A.Aligia,and H.A.Ceccatto,Phys.Rev.B64,121105(R)(2001).10M.Nakamura,J.Phys.Soc.Jpn.68,3123(1999);Phys.Rev.B61,16377(2000).11M.Tsuchiizu and A.Furusaki,Phys.Rev.Lett88,056402(2002).12A.A.Aligia,Europhys.Lett.45,411(1999).13T.Wilkens and R.M.Martin,Phys.Rev.B63,235108(2001).14A.P.Kampf,M.Sekania,G.I.Japaridze,and P.Brune,J.Phys.C(in press).15A.A.Aligia,K.Hallberg,C.D.Batista,and G.Ortiz,Phys.Rev B61,7883(2000).16G.I.Japaridze and A.P.Kampf,Phys.Rev B59,12822(1999).17A.A.Aligia and L.Arrachea,Phys.Rev.B60,15332(1999).18L.Arrachea,A.A.Aligia and E.Gagliano,Phys.Rev.Lett.76,4396(1996);references therein.19S.R.Manmana,V.Meden,R.M.Noack,and K.Sch¨o nhammer,cond-mat/030774120C.D.Batista and A.A.Aligia,Phys.Rev.B47,8929(1993).21L.F.Feiner,Phys.Rev.B48,16857(1993).22M.E.Simon,A.A.Aligia,and E.R.Gagliano,Phys.Rev.B56,5637(1997).23J.Eroles,C.D.Batista,and A.A.Aligia,Phys.Rev.B59,14092(1999).24H.Eskes,A.M.Oles,M.B.J.Meinders and W.Stephan,Phys.Rev.B5017980(1994);references therein.25F.Lema and A.A.Aligia,Phys.Rev.B55,14092(1997);Physica C307,307(1998).26See for example J.des Cloizeaux and J.J.Pearson,Phys.Rev.128,2131(1962).27M.Calandra and S.Sorella,Phys.Rev.B57,11446(1998).28C.Kittel,Quantum Theory of Solids,(John Wiley&Sons,New York,1987).29A.O.Gogolin, A.A.Nersesyan,and A.M.Tsvelick,Bosonization and strongly correlated systems(UniversityPress,Cambridge,1998).30T.Giamarchi and H.J.Schulz,PRB39,4620(1989).。