arXiv:cond-mat/0203258v1 [cond-mat.str-el] 12 Mar 2002 AS 71.10.-w,71.27.+a PACS: pnpolmi oi tt hsc.Iscoeconnection close Its high- of physics. phenomena the state with solid in problem open T days. present in relevant problem this operators oe r eodcas h aoia uniainof quantization canonical The class. second are and model Faddeev-Jackiw the methods con- using Dirac first-order was by of point family Lagrangians particular starting strained a Our of construction decoupling. the any without ables, X t fve st rt ahitga o - oe by model t-J for integral path a Hubbard write the to directly point using is Our view Hubbard-operators. of for formulation integral operators eetto fthe of resentation oko lv oo representation Kadanoff boson slave of work parameter. inrlso h ubr operators Hubbard the of commuta- quantities. rules complicate physical a)the tion of mainly: are calculation difficulties the These in difficulties eral Hubbard-operators using written is when ilb eeati h rsn ae)i h large- the is paper) present (which the techniques in nonperturbative relevant the be of will one b), point lv atce saebsnadfrin odcul the decouple to fermion) and original boson (slave particles slave small observables. be the to of taken evaluation be perturbative can parameters a model in the of any b) e ie(e eo)ad1 and below) (see site per pansion. − c ag- xaso ae nteHbadoeao ahinteg path Hubbard-operator the on based expansion Large-N -operators. autdd inisEatsIgne´ayArmnuaand Agrimensura Ingenier´ıa y Exactas Ciencias de Facultad h oeo lcrnccreain sa motn and important an is correlations electronic of role The nti ae edvlpalarge- a develop we paper this In high- of context the in models popular most the of One nasre fpapers of series a In h large- The ohnl h on ) oua ehdi ouse to is method popular a a), point the handle To uecnutvt sthe is superconductivity J oe ae ntept nerlrpeetto for representation integral path the on based model X 7 N 8 2 , , oeao nuulbsn n fermions and bosons usual in -operator 9 3 stenme feetoi ere ffreedom of degrees electronic of number the is . lhuhthe Although . ucinlter ntrso the of terms in theory functional ai ue,i hc h rpgtr n etxaewritten for are calculated, vertex have we and formulation propagators large-N the our which in rules, matic nerlfruainfrHbadoeaos u ag- e large-N Our Hubbard-operators. for formulation integral eas aecluae h pnsi n hrecag cor charge-charge and spin-spin the calculated have also We ehv oprdordarmtcnqeadrslswt h e the with results and technique diagram our compared have We N ntepeetwr ehv eeoe ag- xaso for expansion large-N a developed have we work present the In 12 xaso a eeoe nteframe- the in developed was expansion esoe httecntan fthe of constrains the that showed We . .INTRODUCTION I. t − J 10 /N oe si em fHubbard of terms in is model , 11 t T a ecniee sasmall a as considered be can X − c ehv tde h path the studied have we t superconductivity − oeaosa ed vari- fields as -operators J J oe.Antrlrep- natural A model. 5 2 N , t plcto othe to application its oe squadratic, is model 6 Xoeaos and, (X-operators) n,uigBaym- using and, xaso o the for expansion din osasadAndr´es Greco and Foussats Adriana 3 hr r sev- are there , 5-00Rosario-Argentina. 250-2000 1 Otbr2,2018) 29, (October 4 makes N For . ex- X - J 1 nttt eFıiaRsro(N-OIE) Av.Pellegrin (UNR-CONICET). F´ısica Rosario de Instituto ,terenormalized the 0, = nscinVw ieormi icsin n conclu- and discussions main our give we Finally, sions. V expansion. large-N section of in order correlation leading charge-charge to and functions spin-spin calculate we IV, adoperators bard our problem, follows: the as attack make works To method them X-operators, nontrivial. highly between using model rules written this commutation is complicate when the quadratic for is point model starting the are expansion. which large-N [11], our in obtained sults epeetepii acltosfrtecase the method. for large-N ex- calculations the our large-N explicit develop we our present III, We understanding section for In useful pansion. be will our which for point starting the is expansion. 11] large-N [ in defined formulation diino hmclpotential chemical a of addition nerlcnb apdt hto h lv-oo rep- slave-boson gauge the radial of the that path- to in mapped corresponding resentation be the can that integral showed and Lagrangians La- effective and constructed. were functional path-integral grangian using correlation rules by the Next, commutation techniques, the Hubbard-operators. to the leads of theory constrained this hr the where )W osrc arninwitni em fHub- of terms in written Lagrangian a construct We a) re- the although main introduction, the in the mentioned we of As summary a give we section this In nscinI egv umr fteformalism the of summary a give we II section In nRf[1] efudapriua aiyo constrained of family particular a found we 11], Ref.[ In t I RLMNRE N DEFINITIONS AND PRELIMINARIES II. t − − pninfruaincnan diagram- contains formulation xpansion eainfntost edn re 1 order leading to functions relation J ntr fHbadoeaos Using operators. Hubbard of term in J oe ( model L the model itn nsi h literature. the in ones xisting = t − t X − i J U J O a ifiieHbadmdl.I section In model). Hubbard -infinite oe Hamiltonian model iαβ (1 oe ae ntepath the on based model /N ( X a ersnainand representation ral oo propagator. boson ) ) X ˙ i αβ 13 µ − ed as: reads , , 5 H h path-integral The . ( X H (1) ) /N ( X . ,wt the with ), J of 0 = t − i J 1 1 1 H(X)= t Xσ0X0σ + J (Xσσ¯ Xσσ¯ (sdetM ) 2 =1/ . (8) ij i j 2 ij i j AB ( X00)2 i,j,σX Xij;σ − i XσσXσ¯σ¯ ) µ Xσσ. (2) is just the superdeterminant of the Dirac matrix12 (or − i j − i Xi,σ equivalently the sympletic matrix of the Faddeev-Jackiw method) formed with the set of all the second class con- In (1) the coefficients aiαβ(X) are unknown and must strains of the theory. 00 1 be determined. The power two on X in (sdetMAB) 2 appears because The indices α and β can take both 0 value (empty we are working with two spin projections. This point will state) or spin index σ = (up and down state, respec- be very important for our large-N extension in the next ± ′ tively). The five Hubbard Xˆ-operators Xσσ and X00 are section. boson-like and the four Hubbard Xˆ-operators Xσ0 and Note the first constrain (4) (the first δ-function in (7)) 0σ 2 X are fermion-like . In (2), tij and Jij are hoping and is the completeness condition which means that ”double exchange parameters, respectively, between sites i and j. occupancy” at each site is forbidden. The constrains (5) b) We introduce (via Lagrangian multipliers λαβ) a set (the second δ-function in (7)) is obtained in our formal- of bosonic constrains Ωαβ which must be also determined. ism as a consequence of imposing that the commutation c) After these definitions and proposals, we impose rules of X-operators must be fulfilled. that the Dirac brackets12 between X-variables obtained From the above equations our theory corresponds to a using our constrained theory verify the correct commuta- configuration in which the bosons are totally constrained tion rules for Hubbard operators. With this prescription, and the dynamics is carried out only by the fermions. we obtain a set of differential equations for the coefficients As it was shown in Ref[11], the path integral formalism aiαβ(X) and constrains Ωa. A particular solution for the corresponding to this dynamical situation is mapped in coefficients is the slave boson representation. Equation (7) is our main formula to develop the large- i σ0 i 0σ N expansion. ai0σ = 00 Xi , aiσ0 = 00 Xi (3) 2Xi 2Xi
and for the constrains, III. THE LARGE-N EXPANSION Ω00 = X00 + Xσσ 1=0 , (4) i i i − In this section, we present the large-N expansion in the Xσ framework of the path integral representation for Hub- bard operators. We present our explicit results for the ′ σ0 0σ ′ ′ X X J = 0case(U-infinite Hubbard model) which will be very Ωσσ = Xσσ i i =0 . (5) i i − X00 useful to show how our method works and to compare our i results with previously ones found in the literature. The boson-like Lagrangian coefficients are all zero. Our starting point is the partition function (7) written The constrains are second class because none of them in the Euclidean form: ′ commutes with all the other ones. σ0 0σ ′ X X Consequently, the dynamics is given by the Lagrangian Z = Xαβ δ[X00 + Xσσ 1] δ[Xσσ i i ] D i i i − i − X00 Z Xσ i 0σ σ0 σ0 0σ i (X˙ i X + X˙ i X ) 1 L(X, X˙ )= i i H(X) (6) (sdetM ) 2 exp ( dτ L (X, X˙ )) . (9) −2 X00 − AB i E Xi,σ i × − Z ˙ and the constrains (4) and (5). The Euclidean Lagrangian LE(X, X) in (9) is:
Now, in order to write the path integral representa- 0σ σ0 14 σ0 0σ tion we use the method of Faddeev-Senjanovich . This 1 (X˙ i X + X˙ i X ) L (X, X˙ )= i i + H(X) . (10) method is currently used in quantum field theory to ob- E 2 X00 Xi,σ i tain the path integral of constrained systems. The partition function reads as: First, we make the following steps in the path integral ′ σ0 0σ (9): ′ X X ′ Z = Xαβ δ[X00 + Xσσ 1] δ[Xσσ i i ] a) we integrate over the boson variables Xσσ using the Z D i i i − i − X00 Xσ i second δ-function in (9), 1 b) the spin index σ = , is extended to a new index p 2 (sdetMAB) exp (i dt L(X, X˙ )) . (7) ± × i Z running from 1 to N. In order to get a finite theory in the N-infinite limit, we re-scale the hoping tij to tij /N where the superdeterminant in the Hamiltonian,
2 c) the completeness condition (X00 + Xpp = N/2) N i p i 1 + ˙+ 2 Leff = f˙ipf + f fip (1 δRi + δR ) can be exponentiated, as usual, by usingP the Lagrangian −2 ip ip − i i,p multipliers λi, X d) in order to prepare the path integral for the large-N N + + 2 expansion, we write the boson fields in terms of static + tij rof fjp µ f fip(1 δRi + δR ) ip − ip − i mean-field values and dynamic fluctuations i,j,pX Xi,p + 2 00 + N r0 δλi δRi + fipfip(1 δRi + δRi ) δλi Xi = Nr0(1 + δRi) , (11) − Xi Xi,p † 2 1 δRi + δR ip, (18) − Zip − i Z λi = λ0 + δλi , (12) Xip
e) finally, we make the following change of variables where we have changed µ to µ λ0 and dropped con- stant and linear terms in the fields.− + 1 p0 Looking at the effective Lagrangian (18), the Feynman fip = Xi , (13) √Nro rules can be obtained as usual. The bilinear parts give rise to the propagators and the remaining pieces are rep- resented by vertices. Besides, we assume the equation 1 0p fip = Xi , (14) (18) written in the momentum space once the Fourier √Nro transformation was performed. By following the steps a-e, we find the effective La- In summary, the Feynman rules are: grangian for the J = 0 case i) Propagators: We associate with the two component a 1 + ˙+ 1 δX = (δR, δλ) boson field, the propagator Leff = f˙ipf + f fip −2 ip ip (1 + δR ) Xi,p i 1 0 Nr0 D (q,ωn)= (19) N (0)ab 1 0 + + 1 Nr0 + tij rof fjp (µ λ ) f fip ip − − 0 ip (1 + δR ) i,j,pX Xi,p i which is represented by a dashed line in Fig.1 connecting 1 two generic components a and b. + N r δλ δR + f +f δλ 0 i i ip ip (1 + δR ) i The quantities q and ωn are the momentum and the Xi Xi,p i Bose Matsubara frequency of the bosonic field, respec- + Lghost. (15) tively. We associate with the N-component fermion field fp, In (15), Lghost is defined via the exponentiation of the the propagator superdenterminant using ghost fields. After the exten- sion to large-N, the superdeterminant becomes: δpp′ G(0)pp′ (k,νn)= (20) −iνn (εk µ) 1 N − − 2 (Nr0) (sdetMAB)i = N . (16) which is represented by a solid line in Fig.1 connecting −1 two generic components p and p’. 1+δRi In (20) the electron bare dispersion was defined as εk = 2tro(coskx + cosky), where t is the hopping between The numerator of (16) together with the Jacobian of the − transformations (11-14) contribute to the path-integral nearest neighbors sites on the square lattice. normalization factor. The denominator of (16) can be The quantities k and νn are the momentum and the seen as a N N diagonal matrix and the integral repre- Fermionic Matsubara frequency of the fermionic field, re- sentation is given× in terms of complex boson ghost fields spectively. We associate with the N-components ghost field p′ , p. Therefore, Z the propagator Z
† 1 Lghost( )= ip . (17) pp′ = δpp′ (21) Z − Zip 1+ δR Z D − Xip i which is represented by a dotted line in Fig.1 connecting To implement the 1/N expansion, the nonpolynomial two generic components p and p’. Leff should be developed in powers of δR. Up to order 1/N, it is sufficient to retain terms up to δR2. Then, ii)Vertices: The expressions of the different three-leg and four-leg vertices are
′ i ′ Λpp = ( 1) (ν + ν′)+ µ ; 1 δpp (22) a − 2
3 representing the interaction between two fermions and N 1 ΠRR ( q,ωn)= [2 nF (εk µ)(εk+q εk) −Ns 4 − − one boson (Fig.1a); Xk
′ i ′ 1 ′ 2 [nF (εk+q µ) nF (εk µ)] pp (ν + ν ) µ pp + (ε + ε ) − − − , (26) Λ = ( 1) − 2 − − 2 δ , (23) k+q k ab − 1 0 iωn + εk+q εk − 2 − − representing the interaction between two fermions and N 1 two bosons (Fig.1b); ΠλR(q,ωn)= (εk+q + εk) −Ns 2 Xk a Γ ′ = ( 1)(δpp′ , 0), (24) pp [nF (εk q µ) nF (εk µ)] − + − − − (27) × iωn + εk q εk representing the interaction between two ghosts and one − + − boson (Fig.1c); and and,
ab 1 0 ′ N [n (ε µ) n (ε µ)] Γpp′ = ( 1) − δpp F k+q F k − 0 0 Πλλ(q,ωn)= − − − . (28) −Ns iωn + εk+q εk (25) Xk − − From the Dyson equation (D )−1 = (D )−1 Π representing the interaction between two bosons and two ab (0)ab ab the dressed components D of the boson propagator− can ghosts (Fig.1d). Each vertex conserves the momentum ab be found. and energy. Our expressions (26)-(28) agree with the previously ir- In the N-infinite limit (order O(1)), we only have free reducible boson self energies calculated in the context of fermions with renormalized band due to correlations, the large-N slave boson approach5. ε = 2tr (cosk + cosk ). For a given value of µ, r k 0 x y 0 The ghost fields interact only with the boson fields (which− from (11) and the completeness condition is equal (Fig.1c and Fig.1d). The only role of the ghost fields, in to δ/2, where δ is the hole doping away from half-filling) order 1/N, is to cancel infinities in the boson self energy must be determined self-consistently. This result is in Π due to the frequency dependence of our vertices (22) agreement with previous calculation13,8. ab and (23). For the evaluation of quantities to higher orders The bare boson propagator (19), which is of order than 1/N, we must include and check the contribution of O(1/N), is renormalized by a series of electronic bubbles the ghost fields in each particular calculation. Moreover which contribute also in O(1/N). Then, in O(1/N), the it should be noted that in this case we must consider boson propagator is dressed by electronic interactions. terms beyond δR2 in developing the nonpolinomial parts By looking at the diagrammatic, it can be seen that the of L . Although this seems to be complicated, all the irreducible boson self-energy Π is given by the sum of eff ab contributions can be controlled by the small number 1/N the contributions corresponding to the one-loop diagrams via the vertices and propagators. At this point we want (in order 1/N) seen in Fig.1e. to remark that so far not many calculations go beyond the order 1/N. In order 1/N, diagramatic rules are well (d) (a) (b) (c) defined in the present paper. p p a p p a a a In summary, we have developed a diagrammatic tech-
p' p' b p' p' b nique for a large-N expansion in the same spirit of the (e) large-N expansion in quantum field theory. It means, via the order of the propagators and vertices, we can deter- = + + + mine the order of the diagram contribution. Note that our Green’s functions are calculated in terms of the origi- FIG. 1. Types of vertex, up to order O(1/N), between two nal Hubbard operators. For example, from (13) and (14) fermions and one boson (a), two fermions and two bosons we can see that our fermions fip are always proportional (b), two ghosts and one boson (c) and, two bosons and two to the Fermi-like X-operator Xop and not only in the ghosts (d). (e) represents the sum of all one-loop diagrams leading order like in the slave-boson formulation13. contributing to the irreducible boson-self energy.
The presence of the last two diagrams in Fig.1e involv- IV. CHARGE-CHARGE AND SPIN-SPIN ing ghost fields is very important. It is possible to show CORRELATION FUNCTION that these two diagrams exactly cancel the infinities, due to the frequency dependence of our vertices, of the two In this section we calculate numerically the charge- first diagrams in Fig.1e. Therefore, using our Feynman charge and spin-spin correlation functions on the square rules, the final (and finite) result for each component of lattice for the nearest neighbor hopping t. We choose the boson-self energy Πab is: t =1.0 as the energy unit and the temperature T =0.0K.
4 We can define retarded density-density D˜ and spin- ˜ 15 2 spin S Green’s functions as 1.5 1 (a) Q=(π/3,π/3) 1 0.5 D˜ =
5 X-operator8,9 or slave boson path integral approach13. 8 R. Zeyher and M. Kuli´c, Phys.Reb. B 53, 2850(1996). The first one cannot be defined in a diagrammatic form 9 R. Zeyher and A. Greco, Eur. Phys. J.B 6, 473 (1998). 10 while the second one can be developed in terms of Feyn- A. Foussats, A. Greco and O. S. Zandron, Annals of Physics man diagrams but needs the introduction of spinons and (NY) 275, 238(1999); A. Foussats, A. Greco and O. S. Zandron, Annals of Physics(NY) 279, 263(2000). holons. In the formulation presented in this paper we 11 have diagrammatic rules: propagators and vertices in A. Foussats, A. Greco, C. Repetto, O. P. Zandron and O. terms of the original X-operators. S. Zandron, Journal of Physics A 33, 5849(2000). 12 60 One important point is that in our diagram technique L. Faddeev and R. Jackiw, Phys. Rev. Lett , 1692 (1988); P. Dirac, Lectures on Quantum Mechanics (New the fermion operator fσ is proportional to the fermion X- operator X0σ for all orders in the large N expansion. York:Yeshiva University Pressi, 1964); K. Sundermayer, This is very important for the interpretation of the ex- Constrained Dynamics (Berlin: Spinger, 1982); J. Gov- aerts, Int. J. Mod. Physics A 5, 3625(1990). perimental spectral susceptibilities. In the slave boson 13 Z. Wang, Int. Journal of Modern Physiscs B 6, 155(1992). approach the spinon is only proportional to the fermion 14 1 0σ 13 L. Faddeev, Theoret. Mat. Phys. , 1(1970); P. Sen- X operator in the leading order of 1/N expansion . janovich, Ann. Phys. (NY) 100, 227(1976). It is also important to note that our diagram tech- 15 L. Gehlhoff and R. Zeyher; Phys. Rev. B 52, 4635(1995). nique contains ghost fields that cancel infinities to order 16 T. Tohyama, P. Horsch and S. Maekawa, Phys. Rev. Lett. O(1/N) in the boson propagator. In the slave boson ap- 74, 980(1995). proach in the radial gauge, the ghost fields necessary to 17 Y. Izyumov and Y. Skrabin; Statistical Mechanics of Mag- 2 18 treat the Jacobian, cancel infinities in O(1/N ) . netically Ordered Systems, Consultansts Bureau, New Using our large-N expansion, we have calculated the York (1988). boson propagator and the charge-charge and spin-spin 18 N. Read and D.Newns, J. Phys. C 16, 3273(1983). correlation functions in O(1). In sections III and IV we have shown that these quantities agree with those calculated by slave boson13 and X-operator canonical approach15. These results, which are far nontrivial, show the consistency of our formulation. We showed that in the N-infinite limit we have well de- fined electronic excitations with renormalized band due to correlations. This means (and we think) that the large- N expansion, in the form presented in this paper, is a better approximation to describe high-Tc cuprates in op- timal doping and overdoped than in underdoped close to antiferromagnetism. We believe that our formulation will be useful (and simpler than previous ones) to calculate quantities in the next order in O(1/N), such as superconductivity and electronic self-energy corrections. We are also planning to extend the method to the case of finite J. Acknowledgments The authors thank to O. Zandron for valuable discus- sions. A.G. wishes to acknowledge many interesting dis- cussions with A. Dobry, O. Foj´on and R.Zeyher.
1 P.W. Anderson, The Theory of superconductivity in High- T c Cuprates (Princeton University Press, Princeton, New Jersey, USA,1997). 2 J. Hubbard, Proc. R. Soc. London Ser. A 276, 238 (1963). 3 A. Izyumov, Physics-Uspekhi 40, 445 (1997). 4 J. C. Le Guillou and E. Ragoucy, Phys. Rev. B 52 , 2403 (1995). 5 G. Kotliar and J. Liu, Phys. Rev. Lett. 61 ,1784 (1988). 6 M. Grilli and G. Kotliar, Phys.Rev. lett. 64, 1170 (1990). 7 G. Baym and L. Kadanoff, Phys. Rev. B 124, 287 (1961).
6