Creation an estruction of a spin gap in weakly couple quarter-fille laers B. Eegger,, H.G. Evertz, an R.M. Noack Institut für Theoretische Physik, Technische Universität Graz, A-8 Graz, Austria Institut für Theoretische Physik, Universität Frankfurt, D-68 Frankfurt, Germany Fachbereich Physik, Philipps Universität Marburg, D-5 Marburg, Germany arxiv:con-mat/55 v Oct 5 We investigate weakly couple quarter-fille laers with moel parameters relevant for NaV O 5 using ensity-matrix renormalization group calculations on an extene Hubbar moel couple to the lattice. NaV O 5 exhibits super-antiferroelectric charge orer with a zigzag pattern on each laer. We show that this orer, with a perioicity of four laers, causes a spin imerization along the laer an a corresponing spin gap of the same magnitue as that observe experimentally. The spin gap is estroye again at large charge orer ue to a restructuring of the spins. An analysis of an effective spin moel preicts a re-creation of the gap by inter-laer singlets when the charge orer increases further. PACS numbers: 7..F, 7.8.-k, 6..+m The iscovery of a phase transition at T C K [] in NaV O 5, below which charge orer an a spin gap appear [], has precipitate intensive theoretical investigation. NaV O 5 consists of well-separate planes which contain weakly couple quarter-fille vanaium laers [] (Fig.. They can be escribe by an extene Hubbar moel (EHM []. A zigzag charge orer as observe in NaV O 5 [] is alreay create in this moel by the nearest-neighbor Coulomb repulsion V for an isolate laer, but only at overly large values of V [5, 6, 7, 8, 9, ]. In a recent DMRG stuy [8], we showe that the inclusion of a strong effective Holstein coupling to the lattice, which was foun in LDA calculations [], reuces the require Coulomb repulsion to a realistic value. The DMRG calculations then yiele goo agreement between theoretical an experimental results for the amount of charge orer, the extent of lattice istortion in the c-irection, the effective spin coupling J eff in the b-irection, an the charge gap. However, on an isolate laer [7, 8] the spin gap vanishes in the thermoynamic limit. The occurrence of the spin gap in NaV O 5 appears to be intimately connecte to the coupling of laers, as inicate, e.g., by the splitting of magnon branches [,, ]. An intriguing scenario has been put forwar by Mostovoy an Khomskii []. It is base on the fact that the experimentally observe unit cell of NaV O 5 in the orere phase is a b c [5, 6]. The charge orer then has a perioicity of four parallel laers in the a-irection, an is thus super-antiferroelectric (SAF [7]. The corresponing polarization of the electrons on the rungs is illustrate by wie-narrow ellipses in Fig. (top. Each laer exhibits antiferroelectric charge orer along the laer (generally calle zigzag charge orer as well as an antiferroelectric orer to the next nearest laer. The laers correspon to effective antiferromagnetic spin-/ chains in the b-irection. For each laer, the two closest sites on the neighboring in-plane laers alternate between low charge occupation [inicate by ashe loops t a,v t b,v t ab,v b b a = ab ab J J > J J J FIG. : Quarter-fille couple laers in NaV O 5. The mapping of a rung (top onto an effective spin site (bottom is illustrate for a charge-isorere an a SAF-orere regime. The shifte perioic bounary conitions in the a- irection are inicate by numbers which ientify ientical sites. The propose singlet formation for NaV O 5 is epicte by ashe ellipses in the effective spin moel (b, bottom. in Fig. ] an large occupation. Large charge occupation shoul effectively weaken the electron hopping along the laer, both through purely electronic interactions as well as by pushing away the neighboring oxygen atoms an thus reuce the spin coupling. With SAF charge orer, the effective spin chains are therefore imerize, which was propose to lea to the formation of spin singlets [Fig., bottom] an to the observe spin gap. This scenario is ifficult to evaluate quantitatively. It has been investigate using exact iagonalization [, 9, 8] on small clusters an by means of an approximate xy-moel together with a mean-fiel approximation [7]. It has also
been argue that the mechanism propose by Mostovoy an Khomskii cannot work for NaV O 5 [9]. In the present paper, we use the Density Matrix Renormalization Group (DMRG to stuy the extene Hubbar moel with coupling to the lattice [8] for large systems of couple laers up to length, with perioicity of four laers in the a-irection. We show that the scenario propose by Mostovoy an Khomskii inee works. In the thermoynamic limit, it prouces a imerization an a spin gap of similar size to that seen in the experiment. The spin gap appears concurrently with charge orering an is foun only in the case of the -laer SAF orer; it oes not appear for charge orering with - laer perioicity. Surprisingly, the spin gap closes again at large Coulomb repulsion V. We explain this behavior by consiering an effective spin moel for which DMRG calculations on larger systems with lengths of up to rungs are possible. Moel. We use the extene Hubbar moel with coupling to the lattice as introuce in Ref., H = H EHM + H l + H e-l, where H EHM = ( t ij c iσ c jσ + h.c. +U i ij,σ n i n i + ij V ij n i n j, ( with hopping matrix elements from first-principles calculations [,, ], t a =.5 ev, t b =.5 t a, t ab =.7 t a an. t a, an a uniform on-site repulsion U = 8. t a. The lattice an the corresponing moel parameters are illustrate in Fig.. Since the nearest-neighbor Coulomb repulsion V is ifficult to estimate, we investigate the functional epenence on this parameter. Recent LDA calculations [] have shown that there is a very large Holstein-like coupling to out-of-plane movements of oxygen atoms, which can be moele by incluing the terms H l = κ z i i an H e-l = C i z in i, with κ =.5 t a an C =.5 t a, an z i in units of.5 Å, in the Hamiltonian. We treat these movements aiabatically like in Refs. an 8. The optimal configuration for the istortions z i is then a zigzag-pattern [8] with amplitue z opt, for which we take the results from an isolate laer [8]. This simplification is reasonable since we showe in Ref. 8 that istortions remain unchange in the case of two couple laers with perioic bounary conitions (pbc in the a-irection. Then, H l becomes an irrelevant constant an H e-l reuces to a zigzag alternation of the local chemical potentials. In orer to investigate the proposal by Mostovoy an Khomskii, a system with a perioicity of four laers is necessary. Within the DMRG, calculations on an extene Hubbar moel of four couple laers woul be, however, too restricte in length. We therefore employ a system of only two couple laers of length up to, but with shifte perioic bounary conitions, as illus- (J eff...5..5 V=., =.57, t ab =. t a V=., =.57, t ab =.7 t a V=.9, =, t ab =.7 t a quaratic fit (L=8- linear fit (L=- (c Hubbar moel Hubbar moel.5. /L =.5, J ab / b =/5 =., J ab / b =/5 =.9, J ab / b =/5 quaratic fit (L > quaratic fit (L=8- ( eff. spin moel eff. spin moel.5. /L FIG. : Examples of finite-size extrapolations for (a,b the spin gap an (c, the imerization of (a,c the Hubbar moel an of (b, the effective spin moel. Linear an quaratic fits in are illustrate for charge-isorere states L ( =, circles an SAF-orere states ( >. trate in Fig.. They ensure the proper perioicity an the same SAF structure of neighboring charges in the orere phase as in a four-laer system. We apply open bounary conitions (obc in the b-irection. Results. We calculate the spin gap using [7]...5..5 (J b (L = E (L, N, S z = E (L, N, S z =, where E is the groun state energy, L is the length of each laer, N = L/ is the number of electrons, an S z is the total spin in z-irection. The results are extrapolate using linear an quaratic fits in /L, incluing L = (8,,, 6,, as illustrate in Fig.. Due to broken translational invariance in the b-irection (obc, we can efine a spin imerization as = ( i Ŝz i,l (Ŝz i+,l Ŝz i,l, l=i,ii i where i counts the rungs in the b-irection. We restrict i to four rungs in the mile of each laer (l = I, II. Here Si,l z is the sum of the Sz -spin at the i-th rung. The imerization is an inicator for bon alternation an spin singlets along the laer, as illustrate in Fig. (bottom. Translational invariance in the b-irection is aitionally broken by the lattice istortions, an we then fin the contributions to to be positive on both laers. A state consisting of consecutive singlets in the a-irection [Fig., bottom] woul give the maximum value for, i.e., max = /. The charge orer parameter is given by [] m CO = N n e iq(ri Rj ( n i n j n, ij
J eff (mev (t a 5 5...8. J eff z opt quaratic fit linear fit quaratic fit linear fit (c -...6.8...6.8 V (t a.5 an z opt (.5A FIG. : Results for the EHM with coupling to the lattice, as a function of V. Extrapolate values of the charge orer parameter, the spin gap, an (c the imerization are illustrate for two couple laers (t ab =.7 t a. The effective magnetic exchange constant J eff [8] an the istortion z opt [8] of an isolate laer are also shown in. The three regimes iscusse in the text are highlighte. where Q = (π, π an N is the total number of sites on the laer. We calculate m CO for each laer separately an then average. In the thermoynamic limit, we obtain the results shown in Fig. for the charge orer parameter, the spin gap, an the imerization, as a function of V. For comparison, we also show the lattice istortion z opt an the effective spin coupling J eff, which were etermine for an isolate laer in Ref. 8. Here correspons to charge occupations ( ± / of two inequivalent sites an becomes finite concurrently with z opt. We can ientify three ifferent regimes: (i For small nearest-neighbor Coulomb repulsion V.95 t a, the lattice istortion an the charge orer are both zero. The estimation of the spin gap is ifficult, as illustrate by the relatively large isagreement between the linear an the quaratic extrapolations in that regime. From calculations on a system with pure pbc in the a-irection [8] an from consierations of an effective spin moel (see below, we expect a vanishing gap in the thermoynamic limit as long as =. The imerization, which is expecte to be zero in this regime, shows slightly negative values, probably ue to finite-size effects. (ii At a critical V C.95 t a [8], the lattice istortions become finite an cause charge orering in a SAF pattern for our choice of bounary conitions. Linear an quaratic extrapolations to L = match quite well an we fin that the spin gap an the intra-laer imerization are finite. This is in contrast to the non-saf charge-orere system of two couple laers with purely perioic bounary conitions stuie in Ref. 8, which shows neither a finite spin gap nor finite imerization in the thermoynamic limit. The concurrent appearance of intra-laer imerization an a spin gap only when SAF charge orering occurs strongly inicates that here charge orering inee causes singlet formation along the laer. (iii At larger V.6 t a, the spin gap vanishes again, to our initial surprise. Concurrently, the imerization becomes smaller, but oes not completely isappear. This behavior an the negative extrapolations for the spin gap are consistent with the effective spin system (see below. Regime (ii can be associate with NaV O 5 because it provies a escription for the concurrent appearance of lattice istortions, zigzag charge orer, an the opening of a spin gap in the low temperature phase. At the effective V =. t a previously etermine [8], the spin gap is of the same magnitue as in NaV O 5 ( exp S mev =.9 t a []. Calculations for ifferent interlaer hopping strengths t ab show that the maximal size of the gap increases with t ab (t ab = /.7/. t a,max /.5/.5 t a, but that concurrently the transition from regime (ii to (iii appears at smaller V (for t ab =. t a at V. t a, corresponing to smaller. In summary, finite SAF charge orering an inter-laer hopping t ab are both necessary to prouce imerization an a spin gap, while very large charge orer estroys the spin gap again. Effective spin moel. In orer to better unerstan the behavior of the spin gap in the EHM, we consier an effective spin moel for NaV O 5 which was introuce in Ref.. This moel allows us to investigate much larger systems up to length L = an thus to extrapolate results to the thermoynamic limit more accurately. By replacing the quarter-fille rungs (one electron on each rung by a single effective spin site, we obtain the the Heisenberg spin moel whose lattice structure is epicte in Fig. (bottom with two ifferent magnetic exchange constants, Jb along the effective chains, an Jab, inter-chain. The chains in this effective spin moel correspon to the laers of the original system. In the SAF charge-orere state with weakly couple effective chains [Fig. ], the magnetic exchange Jab ifferentiates to three ifferent values epening on the positions of the electrons on the interacting rungs []: large J = Jab ( +, meium J = Jab (, an small J = Jab (, as illustrate in Fig.. The bon alternation along the b-irection is taken to be J / = J b ( ± δ, where J b = Jb (, an we set δ( = = an δ(.. []. Assuming a linear epenence yiels δ(.76. Together with Jab / b = /5 [] (or ab / b = /5,
(J b...5 obc, eff. spin moel, J ab / b =/5 obc, eff. spin moel, J ab / b =/5 pbc, eff. spin moel, J ab / b =/5 pbc, eff. spin moel, J ab / b =/5 eff. spin moel, J ab / b = Hubbar (t ab =.7...6.8 FIG. : Effective spin moel. Quaratically extrapolate values (L = for the spin gap an the imerization, as a function of the charge orer parameter. Calculations with obc (black ots an pbc (open circles in the b-irection are shown. For comparison the spin gap for the EHM is shown, in units of J eff [Fig. ]. where J ab an b are the exchange constants at =, the parameters of the effective spin moel are completely etermine an the behavior can be investigate as a function of the strength of the charge orer parameter,. Like in the Hubbar system, we apply the DMRG to two couple spin chains (corresponing to two Hubbar laers with shifte pbc in the a-irection an obc as well as pbc in the b-irection. Due to the small interchain couplings, the finite-size effects in the a-irection are small [], justifying calculations on only two couple effective chains. A quaratic extrapolation [Fig.,(] of the finite system sizes (L = for obc, L = 6 for pbc to the thermoynamic limit provies the - epenence of imerization an spin gap. The imerization an the spin gap increase with up to a maximum (Fig. ; then, both ecrease even though the bon alternation still increases. This behavior matches well with the results from the EHM. In contrast, uncouple chains (D spin-/-heisenberg chains show a continuing increase of imerization an spin gap. The estruction of the spin gap is likely ue to a restructuring of the singlets, which now occur with a larger probability on the J -bons. For small charge orering, all inter-chain exchanges are of comparable size an are small in units of J b. The singlet formation along the b-irection is then only marginally influence. On the other han, a large value of rastically changes the inter-chain exchanges, e.g., for =.9 an J ab / b = /5, they become J /J b., J /J b., an J /J b.. The large inter-chain exchange J estroys the spin gap an causes a restructuring of the singlets. For pbc along the laer [open circles in Fig. ] a re-creation of the spin gap is observe after passing a critical charge orering m CO. In contrast, obc-calculations (black ots show a vanishing gap for all m CO, although both calculations match excellently for < m CO [Fig. ]. This strongly suggests a singlet formation on the J -bons [5] with increasing J. In the case of obc, unpaire spins are then present at the ens of the laers [Fig. ], preventing a spin gap. This explanation is supporte by a consierable increase of inter-chain spin correlations on the J -bons for > m CO. For obc, after passing the critical charge orer m CO, the spin gap extrapolates to values that are slightly negative [Fig. ]. If the extrapolation is restricte to smaller system sizes [otte line in Fig. ], the values become more negative. This is in agreement with the negative extrapolations for in the Hubbar moel in the large -regime [Fig. ]. The magnitue of these negative values gives a minimal estimate of the error in the extrapolation. Conclusions. We have stuie couple quarter-fille laers in an extene Hubbar moel with lattice coupling an with moel parameters appropriate for NaV O 5. In a super-antiferroelectric charge-orere regime, we fin spin imerization along the laer which generates a spin gap that matches the experimental result. This provies strong support for the mechanism of Mostovoy an Khomskii []. The spin gap will likely be further enhance by in-plane atomic movements ue to the super-antiferroelectric orer [, 9,,, ]. The super-antiferroelectric charge-orere state provies a possibility to unerstan the concurrent appearance of charge orer, lattice istortions, an spin gap in NaV O 5. It is in agreement with experimental observations [5] an explains the puzzling properties of the low-temperature phase. Our calculations match the effective spin moel propose by Gros an Valenti [] an show a vanishing spin gap at very large charge orering which is cause by an increase inter-laer magnetic exchange an the corresponing restructuring of the spins. This work was supporte by the Austrian science fun FWF, project P55. We thank M. Aichhorn, C. Gros, T.C. Lang, an F. Michel for useful conversations. [] M. Isobe an Y. Uea, J. Phys. Soc. Jpn. 65, 78 (996. [] For a review, see P. Lemmens et al., Phys. Rep. 75, (. [] H. Smolinski et al., Phys. Rev. Lett. 8, 56 (998. [] B. Grenier et al., Phys. Rev. Lett. 86, 5966 (. [5] H. Seo an H. Fukuyama, J. Phys. Soc. Jpn. 67, 6 (998; J. Phys. Chem. Sol. 6, 95 (999. [6] D. Sa an C. Gros, Eur. Phys. J. B 8, (. [7] M. Vojta et al., Phys. Rev. B 6, 55 (; Phys. Rev. B 6, R87 (999.
5 [8] B. Eegger et al., Phys. Rev. B 7, 85 (5. [9] C. Gros an G.Y. Chitov, Europhys. Lett. 69, 7 (5. [] M. Aichhorn et al., Phys. Rev. B 69, 58 (. [] J. Spitaler et al., Physica Scripta T9, 59 (; Phys. Rev. B 7, 57 (. [] T. Yosihama et al., J. Phys. Soc. Jpn. 67, 7 (998. [] C. Gros an R. Valentí, Phys. Rev. Lett. 8, 976 (999. [] M.V. Mostovoy an D. I. Khomskii, Soli State Commun., 59 (; M. V. Mostovoy et al., Phys. Rev. B 65, 6 (. [5] S. van Smaalen et al., Phys. Rev. B 65, 6 (. [6] H. Nakao et al., Phys. Rev. Lett. 85, 9 (; S. Grenier et al., Phys. Rev. B 65, 8 (; H. Sawa et al., J. Phys. Soc. Jpn. 7, 85 (; K. Ohwaa et al., Phys. Rev. Lett. 9, 6 (5. [7] G. Y. Chitov an C. Gros, Phys. Rev. B 69, (; C. Gros an G. Y. Chitov, Europhys. Lett. 69, 7 (5. [8] J. Riera an D. Poilblanc, Phys. Rev. B 59, 667, (999. [9] R.T. Clay an S. Mazumar, Phys. Rev. Lett. 9, 76 (5 [] P. Horsch an F. Mack, Eur. Phys. J. B 5, 67 (998. [] We checke this with the DMRG by using effective spin systems with couple chains an pbc in the a-irection for smaller systems (L at =./.5/.8. The results for the imerization an the spin gap match well with our calculations on two couple laers. [] P. Thalmeier an P. Fule, Europhys. Lett., (998. [] A. Bernert et al., Eur. Phys. J. B, 55 (. [] A. Bernert et al., Phys. Rev. B 66, 658 (.