arxiv:cond-mat/ v1 [cond-mat.str-el] 17 Feb 2003

Transcription

1 Charge-ordered ferromagetic phase i magaites Tra Mih-Tie Istitute of Physics, NCST, P.O. Box 429, Boho, Haoi, Vietam, Electrophysics Departmet, Natioal Chiao Tug Uiversity, Hsichu 300, Taiwa. arxiv:cod-mat/ v1 [cod-mat.str-el] 17 Feb 2003 A mechaism for charge-ordered ferromagetic phase i magaites is proposed. The mechaism is based o the double exchage i the presece of diagoal disorder. It is modeled by a combiatio of the Isig double-exchage ad the Falicov-Kimball model. Withi the dyamical mea-field theory the charge ad spi correlatio fuctio are explicitely calculated. It is show that the system exhibits two successive phase trasitios. The first oe is the ferromagetic phase trasitio, ad the secod oe is a charge orderig. As a result a charge-ordered ferromagetic phase is stabilized at low temperature. PACS umbers: a, d, m There has bee much recet iterest i the properties of doped magaese oxides R 1 x A x MO 3 (R=rare earth, A=Ca,Sr). 1,2 These materials preset a very rich phase diagram ivolvig phases with spi, charge ad orbital order. The physically relevat electros i magaites are those from the M 3d levels, which are split by the cubic crystal field ito triply degeerate t 2g levels ad higher-eergydoubly degeerate e g levels. Electros from the e g levels are able to hop betwee M sites ad form acoductio bad. Electrosfrom the t 2g levelsare localized. The itierat electros ad local spis are correlated by the double-exchage (DE) mechaism. 3,4 The mai feature of the DE is a cooperative effect where the motio of a itierat electro favors the ferromagetic (FM) orderig of local spis ad, vice versa, the presece of the FM order facilitates the motio of the itierat electro. The DE model qualitatively describes some of the magetic properties of magaites, 2,5 ad provides a well-established startig poit toward comprehesive uderstadig of the phase diagram of magaites. Recetly experimets have show that beside the FM order a charge order ca exist i the magaites. 6,7 The charge order exists i regios with o et magetizatio ad, surprisigly, ca also occur i FM regios. 7 Dopig of A 2+ ios creates M 4+ holes i a M 3+ backgroud. The presece of two valece states M 3+ ad M 4+ may lead the compouds to a charge-ordered (CO) state for appropriate dopig. However, the DE model aloe caot explai the CO state which coexists i the FM phase. I priciple, the earest-eighbor Coulomb repulsio may stabilize a CO state. However, a large earesteighbor repulsio likely destabilizes the homogeeous FM state ad may produce a checkerboard charge order i three directios. 2 Aother possible mechaism for the CO phase stabilizatio is the couplig of itierat electros to the Jah-Teller distortios. However, the electro Jah-Teller phoo couplig ca oly stabilize a CO-FM state where the CO phase trasitio occurs before the FM trasitio. 8 At half fillig experimets have oly observed a charge order below the FM trasitio temperature. 2,7 Thereforethe Jah-Tellercouplig is ulikely resposible for the appearace of the CO-FM state at least at half fillig. I this paper we preset a possible alterative explaatio for the CO- FM state i the magaites. The key idea is a iterplay of the DE ad radomess of the A-site substitutio. The radomess is ievitably itroduced by experimets. The importace of the radomess has bee discussed both experimetally ad theoretically. 1,2 The radomess ca substatially decrease the critical temperature of the FM trasitio. 9,10,11 Here we will icorporate the radomess of A-site substitutio ito the DE model. For simplicity, we adopt the radomess by A-site substitutio as a radom local potetial of the itierat electros, although the radomess may cause other effects, for istace, radomess of the hopig or exchage itegral. 12 It is well kow that the diagoal disorder with biary distributio ca be modeled by the Falicov-Kimball (FK) model. 13 Although the FK model is simple, it cotais a rich variety of phases. I particular, it illustrates the disorder-order phase trasitio drive by electro iteractio. 14,15 Icorporatig the diagoal disorder of the FK type ito the DE model, oe may expect that a disorder-order phase trasitio could preset. Whe the phase trasitio occurs, a CO-FM phase may be stabilized at low temperature. I order to detect the phase trasitio we study the charge ad spi respose of system by usig the dyamical mea-field theory (DMFT). 16 The DMFT has extesively bee used for ivestigatig strogly correlated electro systems. 16 Withi the DMFT we explicitely calculate the charge ad spi correlatio fuctio. We fid that the system stabilizes a CO-FM state at low temperature. The system which we study is described by the followig Hamiltoia H = t d <ij>,σc iσ c jσ µ iσ 2J H Sis z z i + iσ i E w w i +U iσ w i, (1) iσ i where c iσ (c iσ ) is the creatio (aihilatio) operator of a itierat electro with spi σ at lattice site i; t/ d is the hopig parameter of the itierat electros. Here we have rescaled the hopig parameter with the dimesio

2 2 d of the system. S z i is the z compoet of local spi at lattice site i, ad for simplicity, it takes two values ±1. s z i = ( i i )/2, iσ = c iσ c iσ, w i is a classical variable that assumes the value 1(0) if site i is occupied (ot occupied) by A io. U is the disorder stregth ad is mapped oto the differece i the local potetial which splits eergetically favor of M 3+ ad M 4+ ios. The expectatio value x = i w i /N, (N is the umber of lattice sites), correspods to the cocetratio of ufavorable M 4+ sites. The chemical potetial µ cotrols the carrier dopig, while E w cotrols the fractio of the sites havig the additioal local potetial. We shall use the coditio + x = 1, where = iσ iσ /N is the electro dopig. This coditio determies E w for each dopig. The third term of Hamiltoia (1) is the Hud couplig of itierat ad local electros. For simplicity we oly take ito accout the Isig part of the Hud couplig. This simplificatio does ot allow ay spi-flip processes, which ca be importat at low temperature where spi-wave excitatios may gover the thermodyamics of the system. However, i the DE processes the spi of itierat electro ferromagetically aligs with the local spi, hece, the Isig part of the Hud iteractio plays a domiat role. The DMFT calculatios for thedemodelwithclassicallocalspisshowthatthesimplificatio of the Hud couplig does ot chage the self eergy of the sigle-particle Gree fuctio. 5 Moreover, withi the DMFT the umerical results for quatum local spis do ot show a sigificat differece from the oes for classical local spis. 17 Thus, oe expects that withi the DMFT the simplificatio of the Hud iteractio does ot result i a serious backwardess. J H is the stregth of the Hud couplig, ad i the followig we will take the limit J H. The first three terms of Hamiltoia (1) costitute a simplified DE model. This simplified model captures the most essetial igrediet of the DE processes. The last two terms of Hamiltoia (1) describe a biary radomess of the A-site substitutio. They together with the hopig term form the FK model. 13 It is well kow that withi the FK model the U term iduces a disorder-orderphase trasitio. 14,15 At low temperature a checkerboard orderig phase is stabilized. Hece, the model (1) may display a iterplay of the FM ad CO phase. We solvemodel(1) by the DMFT. The DMFT is based o the ifiite-dimesio limit. I the ifiite-dimesio limit the self eergy is pure local ad has o mometum depedece. The Gree fuctio of itierat electros satisfies the Dyso equatio G σ (k,iω ) = 1 iω ε(k)+µ Σ σ (iω ), (2) where ω = πt(2 + 1), ε(k) = 2t d j=1 cos(k j), ad Σ σ (iω ) is the self eergy. I the ifiite-dimesio limit the bare desity of states of itierat electros becomes ρ(ε) = exp( ε 2 /t 2 )/ πt ad we take t as the uit of eergy (t = 1). The self eergy is determied by solvig a effective sigle-site problem. The effective actio of this problem is S eff = dτdτ c σ(τ) ( τ δ(τ τ )+Λ σ (τ τ ) ) c σ (τ ) σ + dτc σ(τ) ( µ+uw J H σs z) c σ (τ)+e w w, σ where Λ σ (τ) describes the effective medium. This effective sigle-site problem ca exactly be solved. Ideed, the dyamics of the localized spi S z ad impurity w ivolved i the effective actio are idepedet, hece, we could idepedetly take the trace over S z ad w i calculatig the partitio fuctio. This is similar to the DMFT solvig of the FK model. 18 We obtai the local Gree fuctio i the limit J H G σ (iω ) = W 0σ Z σ (iω ) + W 1σ Z σ (iω ) U, (3) where Z σ (iω ) = iω +µ Λ σ (iω ), ad { W ασ = exp α =0,1 σ l [ βe w (α α)+ ( Zσ (iω ) α U )] } 1 Z σ (iω ) αu with α = 0,1. I takig the limit J H i derivig Eq. (3) we must first reormalize the chemical potetial µ µ+j H. The self eergy is determied by the Dyso equatio for the effective sigle-site problem Σ σ (iω ) = Z σ (iω ) G 1 σ (iω ). (4) Withi the DMFT, the local Gree fuctio must coicide with the sigle-site Gree fuctio of the origial lattice, i.e., G σ (iω ) = 1 G σ (k,iω ). (5) N Eqs. (2)-(5) form the complete set of equatios, which self cosistetly determie the self eergy ad Gree fuctio. We are iterested i calculatig the charge(c) ad spi (s) correlatio fuctio k χ c(s) (i,j) = (δ i ±δ i )(δ j ±δ j ) (6) i the homogeeous paramagetic (PM) phase (δ iσ = iσ iσ ). I order to calculate the charge ad spi respose of system oe has to itroduce a exteral field ito the Hamiltoia. The charge ad spi correlatio fuctio ca be obtaied by differetiatig the Gree fuctio respected to the exteral field, ad the takig the zero limit of the field. 18 Followig the stadard techiques, 18 oe ca express the correlatio fuctios i the terms of charge (c) ad spi (s) susceptibility χ c(s) (q,iω ) i mometum space χ c(s) (q) = T 2 χ c(s) (q,iω ). (7)

3 3 The charge ad spi susceptibility ca be obtaied by differetiatio of the Gree fuctio. 18 We obtai 2+ ( Σ(iω) ) W α χ c(s) G,W 1 α Γ c(s) α (q) α=0,1 (q,iω ) = [χ 0 (q,iω )] 1 ( Σ(iω ) ), (8) G(iω ) W 0.10 =0.5 =0.7 =0.9 where χ 0 (q,iω ) = k G(k+q,iω )G(k,iω ). The matrix Γ c(s) (q) satisfies the followig equatio B c(s) (q) Γ c(s) (q) = Q c(s) (q), (9) T c 0.05 where B c(s) (q) ad Q c(s) (q) have the followig elemets B c(s) αβ (q) = δ αβ + Q c(s) α (q) = A c(s) α (iω )η(q,iω )G(iω ) ( Σ(iω ) W β )G,W 1 β 1 G 2 (iω ) ( Σ(iω ) G(iω ) )W +η(q,iω )G(iω ), (10) A c(s) α (iω ) ( G 2 (iω ) ( ) Σ(iω ) )W 1 G(iω ) 1 G 2 (iω ) ( Σ(iω ) G(iω ) )W +η(q,iω )G(iω ) (11) with α, β = 0, 1. I derivig Eqs. (9)-(11) we have used the stadard coversio 18 [χ 0 (q,iω )] 1 = [G σ (iω )] 2 + η(q,iω )[G σ (iω )] 1, ad itroduced quatity A c(s) α (iω ) = δw α /δz (iω )±δw α /δz (iω ). I the ifiite dimesio limit all of the wave vector depedece of χ 0 (q,iω ) ad η(q,iω ) iclude i the term X(q) = d j=1 cosq j/d. Hece, the spi ad charge correlatio fuctio oly deped o mometum via X(q). Each of the derivatives appearig i Eqs. (8)-(11) ca directly be calculated from the DMFT solutio of Eqs. (2)- (5). I such the way B c(s) (q) ad Q c(s) (q) are calculable oce the self-cosistet equatios of the DMFT are solved. Equatio (9) reveals that Γ c(s) (q) will diverge at a temperature where the determiat of B c(s) (q) vaishes, while Q c(s) (q) remais fiite. This results i a T c FIG. 1: The critical temperature T c as a fuctio of dopig for U = 0.5 (squares), U = 1 (circles). The filled (ope) symbols are T c of the FM (CO) phase trasitio. The solid lie is T c of the FM trasitio without disorder (U = 0) U FIG. 2: The critical temperature T c as a fuctio of U for various dopig. The filled (ope) symbols are T c of the FM (CO) phase trasitio. uphysical chage of the sig of correspodig correlatio fuctio χ c(s) (q) so that the assumptio of the homogeeous PM phase fails for lower temperature. By a similar way oe could also calculate the spi correlatio fuctio χ S (q) of local spis. After some calculatios we obtai χ S (q) = γ 0 (q)+γ 1 (q), (12) Bαβ s (q)γ β(q) = 2W α T. (13) β=0,1 From Eqs. (9), (12)-(13) oe ca see that the spi correlatio fuctio of itierat electros ad local spis will diverge at the same temperature where the determiat of Bs (q) vaishes. This meas that the spi of itierat electros parallel aligs with local spi, ad thus is a importat feature of the DE. We calculate the charge ad spi correlatio fuctio (7) by solvig the DMFT set of self-cosistet Eqs. (2)-(5). We are oly iterested i the FM ad checkerboard CO phase stability. Hece, we oly calculate the spi correlatio fuctio at X q = 1 ad the charge correlatio fuctio at X q = 1. It is foud that the spi correlatio fuctio χ s (X q = 1) always diverges at a critical temperature. This is the sigal of the FM phase trasitio. The charge correlatio fuctio χ c (X q = 1) oly diverges for U 0. This meas that without the disorder the system always is homogeeous. I Figs. 1 ad 2 we preset the critical temperature T c of the FM ad CO phase trasitio as a fuctio of dopig ad disorder stregth. At half fillig = 0.5 both critical temperatures reach their maximal value. The T c of the FM trasitio alwaysdecreases with icreasig disorder stregth. This meas that the disorder substatially decrease T c of the FM trasitio. 9,10,11 At the same time, with icreasig U, T c of the CO phase trasitio first icreases, reaches its maximal value, ad the decreases. The behaviorof T c of the CO phase trasitio is similar to the oe i the FK model. 18 At very strog disorder (U 1) the two critical temperatures

4 4 F m B m A FIG. 3: Upper pael: the temperature depedece of the spi magetizatio m A (m B) of sublattice A(B) (the dot ad dashed lie), ad of the charge-order parameter (the solid lie). Lower pael: the temperature depedece of the free eergy F. The solid, log dashed, short dashed, ad dot-dashed lie are the free eergy i the CO-FM, homogeeous FM, CO-PM, ad homogeeous PM phase, respectively (U = 1, µ = U/2). approach to a same value. Oe also otices that T c of the CO phase trasitio always is smaller tha the FM trasitio temperature. Thus, oe may expect that the CO state is stabilized i the FM phase at low temperature. However, this CO phase stability is respected to the homogeeous PM phase, ad for safety we also study a ihomogeeous phase. We divide the lattice ito two peetratig sublattices A ad B. This lattice divisio allows us to study the checkerboard CO phase. By usig the stadard techique 16 the matrix Gree fuctio ca be writte i the followig form Ĝ 1 σ (k,iω ) = T ( iω +µ Σ A σ(iω ) ε(k) ε(k) iω +µ Σ B σ (iω ) ), where Σσ A(B) (iω ) is the self eergy of the Gree fuctio of sublattice A(B). The self eergies are determied by solvig the effective problem of sigle site of the sublattices. 16 We fid that at low temperature a checkerboard CO-FM state is stabilized. We plot the magetizatio m A(B) = 2 i A(B) sz i /N of sublattice A (B) as a fuctio of temperature i Fig. 3 (upper pael). I this figure we also plot the temperature depedece of the charge-order parameter = ( i A,σ iσ j B,σ jσ )/N. It shows that below a critical temperature the magetizatios of both sublattices exist. They equal to each other util aother critical temperature, where the charge-order parameter exists. At low temperature the system is i the checkerboard CO-FM state. I this phase the charge order coexists i the FM state, as experimetally observed. 7 I the way the system exhibits two successive phase trasitios. Iitially the system goes to the homogeeous FM phase, ad after that to the checkerboard CO phase, that the CO-FM phase is stabilized. We also calculate the free eergyofthesystem. Thefreeeergycaolybeexpressed i terms of local quatities. 16 We plot the temperature depedece of the free eergy F i Fig 3 (lower pael). It shows that the CO-FM state has lowest free eergy, hece the state must be stabilized at low temperature. I coclusios, we have proposed a mechaism for the CO-FM phase which has recetly bee observed. The mechaism is based o a combiatio of diagoal disorder ad a simple DE model with local Isig spis. Employig the DMFT we have calculated the charge ad spi correlatio fuctio. It is foud that the FM ad CO state are stabilized at low temperature. As a result the checkerboard charge order ca occur i the FM state. However, the magaites are too complicated a system to be completely described by this simple model. I particular, the phase with ihomogeeous percolatio of FM ad CO regios is beyod the scope of this paper. This work was supported by the Natioal Program of Basic Research o Natural Sciece of Vietam, Project The writig was completed at the Natioal Chiao Tug Uiversity, ad was supported by Project NSC M of ROC. 1 M. B. Salamo ad M. Jaime, Rev. Mod. Phys. 73, E. Dagotto, T. Hotta, ad A. Moreo, Phys. Rep. 344, 1 3 C. Zeer, Phys. Rev. 82, 403 (1951). 4 P. W. Aderso ad H. Hasegawa, Phys. Rev. 100, 675 (1955). 5 N. Furukawa, J. Phys. Soc. Jp. 63, 3214 (1994); ibid. 64, 2754 (1995); ibid. 65, 1174 (1996). 6 M. Uehara, S. Mori, C. H. Che, ad S.-W. Cheog, Nature 399, 560 (1999). 7 J. C. Loudo, N. D. Mathur, ad P. A. Midgley, Nature 420, 797 (2002). 8 S. Yuoki, T. Hotta, ad E. Dagotto, Phys. Rev. Lett (2000). 9 F. Zhog, J. Dog, ad Z. D. Wag, Phys. Rev. B 58, (1998). 10 B. M. Letfulov ad J. K. Freericks, Phys. Rev. B 64, M. Ausleder ad E. Koga, Phys. Rev. B 65, J. Burgy et al., Phys. Rev. Lett. 87, L. M. Falicov ad J. C. Kimball, Phys. Rev. Lett. 22, 997 (1969). 14 T. Keedy ad E. H. Lieb, Physica A 138, 320 (1986); ibid. 140, 240 (1986).

5 5 15 U. Bradt ad R. Schmidt, Z. Phys. B 63, 45 (1986); ibid. 67, 43 (1987). 16 A. Georges, G. Kotliar, W. Krauth, ad M. J. Rozeberg, Rev. Mod. Phys. 68, 13 (1996). 17 K. Nagai, T. Momoi, ad K. Kubo, J. Phys. Soc. Jp. 69, 1837 (2000). 18 U. Bradt ad C. Mielsch, Z. Phys. B 75, 365 (1989); ibid. 79, 295 (1990), ibid. 82, 37 (1991).