arxiv: v3 [cond-mat.str-el] 13 Jan 2014

Transcription

1 Quasiarticle interaction function in a 2D Fermi liquid near an antiferromagnetic critical oint Andrey V. Chubuov Deartment of Physics, University of Wisconsin-Madison, Madison, Wisconsin 53706, USA arxiv: v3 [cond-mat.str-el] 13 Jan 2014 Peter Wölfle Institute for Condensed Matter Theory and Institute for Nanotechnology, Karlsruhe Institute of Technology, Karlsruhe, Germany (Dated: March 8, 2018) 1

2 Abstract We resent the exression for the quasiarticle vertex function Γ ω (K F,P F ) (roortional to the Landau interaction function) in a 2D Fermi liquid (FL) near an instability towards antiferromagnetism. This function is relevant in many ways in the context of metallic quantum criticality. Previous studies have found that near a quantum critical oint, the system enters into a regime in which the fermionic self-energy is large near hot sots on the Fermi surface (oints on the Fermi surface connected by the antiferromagnetic ordering vector q π = (π,π)) and has much stronger deendence on frequency than on momentum. We show that in this regime, which we termed a critical FL, the conventional random-hase-aroximation- (RPA) tye aroach breas down, and to roerly calculate the vertex function one has to sum u an infinite series of terms which were exlicitly excluded in the conventional treatment. Besides, we show that, to roerly describe the sin comonent of Γ ω (K F,P F ) even in an ordinary FL, one has to add Aslamazov-Larin (AL) terms to the RPA vertex. We show that the total Γ ω (K F,P F ) is larger in a critical FL than in an ordinary FL, roughly by an extra ower of magnetic correlation length ξ, which diverges at the quantum critical oint. However, the enhancement of Γ ω (K F,P F ) is highly non-uniform: It holds only when, for one of the two momentum variables, the distance from a hot sot along the Fermi surface is much larger than for the other one. This fact renders our case different from quantum criticality at small momentum, where the enhancement of Γ ω (K F,P F ) was found to be homogeneous. We show that the charge and sin comonents of the total vertex function satisfy the universal relations following from the Ward identities related to the conservation of the article number and the total sin. We show that in a critical FL, the Ward identity involves Γ ω (K F,P F ) taen between articles on the FS. We find that the charge and sin comonents of Γ ω (K F,P F ) are identical to leading order in the magnetic correlation length. We use our results for Γ ω (K F,P F ) and for the quasiarticle residue to derive the Landau arameters F l=0 c = F l=0 s, the density of states, and the uniform (q = 0) charge and sin suscetibilities χ l=0 c = χ l=0 s. We show that the density of states N F diverges as logξ, however F l=0 c,s also diverge as logξ, such that the total χ (l=0) c,s N F /(1 + F l=0 c ) remain finite at ξ =. We show that at wea couling these suscetibilities are arametrically smaller than for free fermions. PACS numbers: 2

3 I. INTRODUCTION Fermi liquid (FL) theory is arguably the most successful low-energy theory of interacting fermions. It states that, as long as the system can be adiabatically transformed from free to interacting fermions, its low-energy roerties are determined by fermionic quasiarticles, which are qualitatively similar to free fermions 1 4. Fundamental characteristics of fermionic quasiarticles, such as the Fermi velocity vf, the effective mass m = F /vf, the residue 1/Z, and the velocities of collective two-article excitations, lie zero-sound waves or sin waves, are exressed via the fully renormalized, antisymmetrized four-fermion interaction Γ αβ,γδ (K,P;K +Q,P Q), taen in the limit of small momentum and frequency transfer Q = (q,ω) (we use the short-hand notation K = (,ω ), the sin indices follow the order of K) In Galilean-invariant systems, the effective mass of quasiarticles and the thermodynamic roerties, lie secific heat and magnetic suscetibility, are determined by the interaction betweenfermionsk andp rightonthefermisurface(fs)[k = K F = ( F,0)andP = P F = ( F,0)] and are exressed in terms of Γ ω αβ,γδ (K F,P F ;K F,P F ) Γ ω αβ,γδ (K F,P F ), which is the limit q = 0 and ω q 0 of Γ αβ,γδ (K,P;K +Q,P Q). This function is roortional to the quasiarticle interaction function introduced henomenologically by Landau 1. Its counterart Γ q αβ,γδ (K F,P), defined by setting ω q = 0 first and then taing q 0, determines the quasiarticle scattering roerties. The quasiarticle residue 1/Z is not determined by the roerties right on the FS, but nevertheless is exressed via an integral involving Γ ω αβ,αβ (K F,P), in which one of the two momenta is on the FS and the other is generally away from the FS (Ref. 2 ). In lattice systems, the thermodynamic roerties of fermions are not determined by Γ ω taen right on the FS, but Ward identities, associated with the conservation of the total number of fermions 2,5 and total sin 6, still allow one to exress the fermionic residue 1/Z, the effective mass m, and the effective magnetic g factor via integrals involving sin and charge comonents of Γ ω α,β,γδ (K F,P) (Ref.[7]). The subject of this aer is the analysis of the form of the fully renormalized Γ ω (K F,P) in a FL near a quantum-critical oint (QCP). We secifically consider a q π = (π,π) commensurate sin-density-wave (SDW) QCP in a 2D metal with a FS lie that of the high T c curates (Fig. 1). Previous wors 8 11 have demonstrated that, unless certain certain vertex 3

4 FIG. 1: The Fermi surface with hot sots (labeled as h.s.). Conjugate airs of hot sots are connected by Q = (π,π), and θ is the angle between Fermi velocities at F,hs and F,hs +Q. From Ref. 8. renormalizations are strong 12, the system near a SDW QCP enters into a regime in which the fermionic self-energy Σ(, ω) develos a much stronger deendence on frequency than on F and Z = 1+ Σ( F,ω)/ ω gets large in hot regions, where shifting F by q π does not tae a fermion far away from the FS We will be calling this regime a critical FL (CFL), following the notation in Ref. 13. Because the self-energy redominantly deends on frequency, the effective mass aroximately scales as Z, and we will see that the renormalization of Z comes from fermions in the vicinity of the FS and involves Γ ω αβ,γδ (K F,P F ). We emhasize that our calculation of the Z-factor refers to the renormalization induced by antiferromagnetic sin fluctuations. These renormalization is usually on to of renormalizations generated at higher energies. The latter are thought to be included in the Hamiltonian we use as our starting oint. The common sense aroach to construct a microscoic theory near a SDW tye QCP is to relace the original four-fermion interaction by the effective interaction mediated by soft collective bosonic fluctuations in the sin channel 14 (the sin-fermion model). This relacement can only be justified in the random hase aroximation (RPA) (see Sec. III below), but the outcome is hysically lausible and we adot the sin-fermion model as the microscoic low-energy model for a CFL. The model describes fermions with the FS as 4

5 shown in Fig. 1 and with sin-sin interaction H sin-sin = q V αβ,γδ (q), c,α c,β c q,γc +q,δ (1.1) where V αβ,γδ (q) = V(q) σ αδ σ βγ (1.2) and the summation over sin indices is assumed. The interaction otential V(q) is eaed at q = q π and near the ea can be aroximated by V(q) = g/[(q q π ) 2 +ξ 2 ], (1.3) where ξ is the sin correlation length and g is the effective sin-fermion couling. The low-energy model is valid when interactions do not tae the fermions far away from the FS. This requires g to be smaller than the fermionic bandwidth, which for the FS in Fig. 1 is comarabletothefermienergye F. Weassumethattherelationg < E F holds. Wewillee ξ large but finite and consider the system s behavior at energies below ω sf (v F ξ 1 ) 2 /g, where the system is still in the FL regime, At higher frequencies, which we do not consider here, the system crosses over into a quantum-critical regime, where it dislays a non-fl behavior with Σ(,ω) ω a with a = 1/2 at the tree level (Refs. 8,9 and 15 ). It is worth reeating that the bare quantities v F E F,g arethose ofquasiarticles alreadyrenormalized by rocesses (e.g., the Kondo effect) occurring at higher energies and not considered here. For an ordinary FL with short-range interaction U(q), the condition that the interaction is weaer thane F imlies that thewea couling aroximationis valid. In that situation, the leading term in the vertex function Γ ω αβ,γδ (K F,P F ) is just the antisymmetrized interaction U(0)δ αγ δ βδ U( F F )δ αδ δ βγ. It seems natural, at first glance, to aly the same strategy near a QCP, i.e., identify Γ ω αβ,γδ (K F,P F ) in a CFL with the effective interaction V αβ,γδ (K F P F ), i..e., identify Γ ω αβ,γδ(k F,P F ) = V( F F ) σ αδ σ βγ = V(K F P F ) ( 3 2 δ αγδ βδ 1 ) 2 σ αγ σ βδ (1.4) where V(K F P F ) V( F F ) is given by (1.3). Additional antisymmetrization of V αβ,γδ (K F P F )isnotrequiredherebecausetheeffectivesin-mediatedinteractionisalready obtained from an antisymmetrized original four-fermion interaction (see Sec. III). The argument for the identification of Γ ω αβ,γδ (K F,P F ) with V αβ,γδ (K F P F ) in a CFL is seemingly quite general (and alicable also to the case g E F ) because within the 5

6 conventional FL strategy Γ ω αβ,γδ (K F,P F ) includes all renormalizations from fermions away from the FS, and V αβ,γδ (K F P F ) is also assumed to include all relevant renormalizations from outside of the FS. However, if we identify Γ ω αβ,γδ (K F,P F ) with V αβ,γδ (K F P F ) and use thewardidentitybetweenthechargecomonentofγ ω α,β,α,β (K F,P F )andthequasiarticlez factor (Eq. (2.7) below), we findahuge discreancy between Z along thefs, comuted this way, and Z comuted using the direct erturbative exansion in owers of g. Namely, for a fermion at a hot sot, Z comuted in the erturbative exansion 8,9 scales as g/(v F ξ 1 ) (modulo logarithmical corrections 8,9,16,17 ), while Z comuted using the FL formula with V(K F P F ) scales as (g/(v F ξ 1 )) 1/2. There is an even stronger discreancy with the sin Ward identity 6, which relates the quasiarticle Z with the sin comonent of Γ ω α,β,α,β (K F,P F ). Namely, if we use (1/2)V(K F P F ) σ αγ σ βδ for the sin comonent of Γ ω α,β,α,β (K F,P), as in Eq. (1.4), we find that the quasiarticle Z becomes smaller than 1, which is obviously incorrect (we recall that we define Z as the refactor of the ω term in the quasiarticle Green s function). As will be shown here, the resolution of the above roblem requires two amendments to the standard microscoic derivation of Fermi liquid theory 2. The first originates from the fact that the effective interaction Eq. (1.3) develos a dynamical character when the QCP is aroached. Physically the dominating dynamical roerties are generated by the Landau daming of sin fluctuations (see Eq. (3.8)). Once the interaction becomes dynamic, the standard argument, showing that quasiarticle contributions do not enter the vertex function Γ ω αβ,γδ (K F,P F ), does not hold any more. Instead, it is necessary to sum u a ladder series of forbidden diagrams involving an irreducible dynamical vertex V eff and the quasiarticlequasiholeroagator. Thequasiarticle-induced contributionstoγ ω leadtoanenhancement of Γ ω over V eff, and the enhancement becomes singular at the SDW QCP. The second amendment is the judicial choice of the irreducible vertex V eff. A lausible guideline is the concet of conserving aroximation roosed by Baym and Kadanoff 18. It amounts to deriving the irreducible vertex by functionally differentiating the self energy with resect to the single article Green s function. If we tae the self-energy in one loo aroximation (see Fig.9) and differentiate it, we then obtain, observing that the sin fluctuation roagator (the wavy line) is comosed of an RPA bubble series, that the irreducible vertex V eff is given by the sum of two different contributions: a single sin-fluctuation roagator and a certain combination of two sin fluctuation roagators, traditionally called 6

7 Azlamazov-Larin (AL) terms (see Fig. 11b). We will show below that the AL terms lay a decisive role for the sin art of Γ ω, in fact fliing the sign and renormalizing the magnitude, such that the critical contributions to Γ ω diverge equally strong in the charge and sin channels. Both amendments have been earlier included in the analysis of a FL near a nematic instability 13,19, where one encounters a CFL driven by nematic fluctuations concentrated near q = 0. In that case the RPA result for Γ ω was found to be (i) corrected by AL terms 19 and (ii) renormalized 13 by the divergent factor of Z. The renormalizations due to AL terms also lay substantial role near a metal-insulator transition in 2D disordered systems 20. In this aer, we analyze the role of both amendments for Γ ω αβ,γδ (K,P) near a SDW transition with q π = (π,π). We show that Γ ω (K,P) again differs substantially from the sin-fermion interaction V(K P). We first demonstrate that AL corrections to the sin art of Γ ω αβ,γδ in Eq. (1.4) have to be included even in the wea couling regime, when the quasiarticle Z is close to 1. We show that the AL contribution to Γ ω αβ,γδ (K,P) for K and P near the Fermi surface is recisely 2V(K P) σ αγ σ βδ. Adding the AL contributions to the V αβ,γδ in (1.4) we obtain the irreducible vertex V eff, which reresents the new bare vertex function V eff (K,P) = 3 2 V(K P)(δ αγδ βδ + σ αγ σ βδ ) (1.5) where V(K P) is given by Eq. (1.3) for K and P on the FS (i.e., K = K F = ( F,0) and P = P F = ( F,0)) and by Eq. (3.8) for K and P slightly away from the FS. We show that V eff (K,P) satisfies the charge and sin Ward identities in the wea couling regime, as it should. We then identify the series of additional contributions to Γ ω (K F P F ), which are small in an ordinary FL at wea couling but become O(1) in a CFL. Emloying the sin structure of the vertex Γ ω αβ,γδ (K,P) = Γ c(k,p)δ αγ δ β,δ + Γ s (K,P) σ αγ σ β,δ, we obtain and solve two integral equations for Γ a (K,P), a = c,s, in the charge and sin sectors. These equations are simlified for K K F and P P F and may be solved by a suitable ansatz once we introduce f, a = Γ a(k F,P F )/Va eff (K F P F ), where scalar variables and are deviations along the FS from the corresonding hot sots and V eff a (K F P F ) = (3/2)V( F F ) are the sin and charge comonents of V eff αβ,γδ (K F P F ) in (1.5), defined analogously to Γ a. We show that f c, = fs, = f,, i.e., the δ αγ δ βδ + σ αγ σ βδ structure of the vertex survives. We comare the quasiarticle residue 1/Z obtained in (i) the direct diagrammatic calcu- 7

8 lation and (ii) using sin and charge Ward identities with the total Γ ω αβ,γδ (K,P), and show that they agree with each other, as they should. In other words, the total vertex function Γ ω αβ,γδ (K,P), which we obtain, satisfies sin and charge Ward identities. We again emhasize that the contributions to Γ ω discussed here and in the following come in addition to the usual contribution originating from the bare interaction renormalized by the incoherent art of the Green s function. The former contributions are generated by the critical fluctuations and are found to dominate near the QCP. Our results are similar, but not identical, to the ones in the nematic case 13. The ey difference is that in our SDW case the function f a, is not a constant, lie it was near a nematic transition, but rather deends strongly on the ositions of the FS momenta and relative to the corresonding hot sots. Secifically, if and are comarable, f, = O(1), i.e., Γ ω (K F,P F ) is roughly the same as the sin-fermion interaction V eff (K F P F ). If, however, one deviation is arametrically larger than the other, f, becomes large and, roughly, acquires an extra factor of Z, lie in the case of a nematic QCP. The hysical consequence of this momentum-sensitive enhancement of f, is the ultimate connection between the FL descrition of fermions in the hot and cold regions on the FS. Namely, Z for a fermion in a hot region, where the FL becomes critical near the QCP, is determined by Γ ω (K F,P F ) in which the characteristic momenta F are located at the boundary to a cold region, where Z remains O(1) even at the QCP. This connection between hot and cold fermions is not easily seen in erturbation theory where Z for a hot fermion is determined solely by fermions in hot regions, at least at one-loo order. We also analyze the interlay between the contributions to f, from rocesses with even and odd numbers of sin-fermion scattering events. For the rocesses with an odd number of scatterings, F and F differ by aroximately q π, and for an even number of scattering events, F and F are close to each other. A similar searation into vertices with small and large F F has been erformed in Ref. 15 in the context of the calculation of the conductivity near a QCP. In our case, the corresonding contributions to f, are f π, and f 0,. The full f, is the sum of the two: f, = f +, = fπ, +f0,. We comute fπ, and f0, searately and find that their difference f, = fπ, f0, is also a highly non-trivial function of and. Our results for f +, = f, and f, are summarized in Figs. 18 and 20. We use the result for f, to determine the density of states (DOS), the Landau function, and the uniform sin and charge suscetibilities. We show that the DOS N F diverges as 8

9 logξ uon the aroach to the SDW QCP. We introduce the Landau function F(K F,P F ) by straightforward extension of theformula relating F(K F,P F ) andγ ω (K F,P F ) inthe isotroic case. Because the quasiarticle Z deends on the osition on the FS, the Landau function F(K F,P F ) deends searately on and rather than on their difference. In this situation, one cannot use FL formulas relating artial comonents of F(K F,P F ) to charge and sin suscetibilities and has to obtain the suscetibilities by exlicitly summing u bubble diagrams with self-energy and vertex corrections. We demonstrate how to do this and ay secial attention to the difference between contributions coming from the infinitesimal vicinity of the FS and from states at small but still finite distances from the FS. We show that the charge and sin suscetibilities (χ c and χ s, resectively) are identical and for both higher-loo terms form a geometrical series, lie in an isotroic FL. We argue that, in this situation, one can effectively describe χ c,s by using a FL-lie formula in which F c (K F,P F ), averaged over both momenta, lays the role of the l = 0 Landau interaction comonent. The aer is organized as follows. In the next section we briefly review basic facts about FL theory, viewed from a microscoic ersective, introduce the fully renormalized, antisymmetrized vertex function, and resent the relation between Γ ω and Z. In Sec.III we derive the vertex function between low-energy fermions near a SDW QCP, but before the system enters the CFL regime. We first resent, in Sec. III A, a conventional, RPA-tye, common-sense derivation of Γ ω. We show that the common-sense Γ ω (K,P) coincides with the effective dynamical four-fermion interaction V(K P) mediated by soft collective excitations in the sin channel. In Sec. IIIB we argue that the RPA analysis is incomlete near a magnetic transition, even in the ordinary FL regime. We show that AL terms are as imortant as RPA-tye terms and the Γ ω, which satisfies both sin and charge Ward identitiesintheordinaryflregime, isnotgivenbythedirectsin-mediatedinteractionv(k P), butbythemodified, effectiveinteractionbetweennear-criticalfermions, V eff αβ,γδ (K P),which is the sum of the direct sin-fluctuation exchange and AL terms. In Sec.IV we argue that Landau daming can be neglected only if one is interested in the behavior of interacting electrons above a certain frequency ω L, but must be et if one is interested in roerties of fermions at the smallest frequencies. We discuss the Landau daming induced crossover between the ordinary (non-critical) FL and the CFL, evaluate Z in a direct loo exansion and show that there is a discreancy between Z obtained this way and Z obtained using the FL formula, if the effective interaction V eff αβ;γδ (K F P) is used for Γ ω αβ;γδ (K F,P). 9

10 Sections V-VII are the central sections of the aer. In Sec. V we argue that Γ ω αβ,γδ (K,P) in a CFL must differ from V eff αβ,γδ (K P) and identify the diagrammatic series for the fully renormalized Γ ω αβ,γδ (K,P) for K K F and P P F, in which V eff αβ,γδ (K P) is the first term. We show that the exression for Z, obtained using the Ward identities with this fully renormalized Γ ω, coincides with Z obtained in a direct loo exansion. In Sec.VI we resent our solution of the integral equation for f, = Γ ω (K F,P F )/V eff (K F P F ) and discuss the momentum-selective enhancement of f, and the inter-connection between hot and cold regions on the FS. We also discuss in this section the interlay between f π, and f0,. In Sec. VII we resent the full exression for the vertex function and discuss the quasiarticle interaction function (the Landau function), the density of states, and the uniform sin and charge suscetibilities. Sec. VIII resents our conclusions. II. FERMI LIQUID, BASIC FORMULAS Fermi liquid theory was first develoed as a henomenological theory for isotroic fermionic systems, obeying Galilean invariance and the conservation laws for total charge and sin, on the ostulate that interactions do not change the relation between the fermionic density and the volume of the FS. Subsequently, the FL theory was alied to conduction electrons in a metal 21, on which we naturally focus here, and was also re-formulated in the microscoic (diagrammatic) language, using the notion of the fully renormalized, antisymmetrized vertex function Γ αβ;γδ (K,P,Q), taen in the limit of small momentum transfer q and small frequency transfer ω q (we will continue to use 3D notation Q = (q,ω q ), etc.). The universal relation between fermionic density and the volume of the FS has been shown diagrammatically, in order-by-order erturbative calculations, and is commonly nown as the Luttinger theorem [2]. We will wor with Matsubara fermionic Green s functions in the limit when the temerature T 0. We slit the fermionic Green s function into quasiarticle and incoherent arts as G(,ω) = G q +G inc. (2.1) The quasiarticle art of the Green s function has the form G q = 1 Z [iω v F,( F )] 1, (2.2) 10

11 where Z is the inverse quasiarticle residue, and vf, = F/m, with m /m = Z, is the renormalized Fermi velocity. We omitted the quasiarticle daming, considering that G q will only be used sufficiently close to the Fermi surface. In isotroic systems, Z and vf are constants, in lattice systems both generally deend on the location of F along the FS. G inc accounts for the incoherent art of the fermionic roagator. In the following we will mae use of the fact that near a QCP the relevant fermionic states are located close to the Fermi surface, and are described by the quasiarticle art of the Green s function. The effects of the incoherent art are assumed to be included in the effective interaction to be described later. FIG. 2: Diagrammatic reresentation of the vertex function Γ α,β,γ,δ (K,P,Q). The shaded vertices are fully renormalized 4-fermion interactions The vertex function Γ αβ,γδ (K,P) is the fully renormalized and anti-symmetrized interaction between quasiarticles. Grahically, Γ αβ,γδ (K,P,Q) is the combination of the two terms shown in Fig. 2. The second one is obtained from the first one by interchanging the two outgoing fermions and changing the overall sign as required by the Pauli rincile. Of articular interest for the FL theory is the limit of Q = (q,ω q ) when q is strictly zero and ω q tends to zero, nown as the ω-limit : Γ ω αβ,γδ(k,p,q) = lim ωq 0Γ αβ,γδ (,0;,0,ω q ;, ω q ) = Γ ω αβ,γδ(k,p) According to FL theory, the function Γ ω αβ,γδ (K F,P F ), taen between articles at the FS (K F = ( F,0)), is roortional to the quasiarticle interaction function (the Landau function) F αβ,γδ (K F,P F ). For isotroic, Galilean-invariant systems Z = Z and m = m are indeendent of the location of F on the FS, and the relation between Γ ω αβ,γδ (K F,P F ) and F αβ,γδ (K F,P F ) is 2 F αβ,γδ(k F,P F ) = 2N F Z 2 Γ ω αβ,γδ(k F,P F ) (2.3) 11

12 where N F m is the fully renormalized DOS er sin at the Fermi level. Although Eq. (2.3) loos lie a simle roortionality, it is actually a highly non-linear integral relation as m isexressed viaaarticularartialcomonentoff αβ,γδ andz isexressed viaanintegral over Γ ω αβ,γδ (K F,P) (see Eq. (2.7) below). Still, as long as Z and m /m are some constants, the functional form of F αβ,γδ (K F,P F ) is the same as that of Γ ω αβ,γδ (K F,P F ). The other limit (termed the q limit ), ω q = 0 and q 0 defines Γ q αβ,γδ (K,P). The latter is related to the quantum mechanical quasiarticle scattering amlitude. Partial comonents of Γ q and Γ ω are simly related 2. For non-isotroic systems, the definition of the Landau function is more nuanced as in general there is no straightforward relation between its artial amlitudes and observables. For ractical uroses, it is convenient to introduce F αβ,γδ (K F,P F ) via a relation similar to (2.3) but with N F Z 2 relaced by N F /(Z Z ), where now N F < m > is the total DOS (m is the momentum deendent effective mass, and <... > is the average over the FS). We will return to this issue in Sec. VII. That the Landaufunction is exressed via Γ ω (as oosed to, e.g., Γ q ) has a clear hysical meaning. We will see in the next Section that in ordinary FL theory Γ ω (K F,P F ) includes all ossible renormalizations of the interactions between fermions at the FS, which come from virtual rocesses with intermediate fermions away from the FS. The rocesses in which all intermediate fermions in the immediate vicinity of the FS are exlicitly excluded from Γ ω (K F,P F ) (but these are resent in Γ(K F,P F,Q) at an arbitrary ratio of ω q and v F q, and, in articular, in Γ q (K F,P F )). In other words, Γ ω (K F,P F ) reresents the fully irreducible interaction between fermions at the FS. Similarly, in Landau FL theory, the Landau function F has the meaning of the effective interaction between quasiarticles on the FS, which absorbs all contributions from virtual fermion excitations outside of the FS. Obviously, F and Γ ω (K F,P F ) have to be exressed via each other. The remaining renormalizations, coming from fermions in the immediate vicinity of the FS, are all catured within FL theory which relates the observables with the artial comonents of F. For SU(2) sin-invariant systems, the vertex Γ ω α,β,γδ (K F,P F ) and the Landau function F α,β,γδ (K F,P F ) can be decouled into sin and charge comonents as Γ ω α,β,γδ (K F,P F ) = Γ c (K F,P F )δ αγ δ β,δ +Γ s (K F,P F )σ αγ σ β,δ F αβ,γδ (K F,P F ) = F c (K F,P F )δ αγ δ β,δ +F s (K F,P F )σ αγ σ β,δ (2.4) In isotroic systems Γ c (K F,P F ) (F c (K F,P F )) and Γ s (K F,P F ) (F s (K F,P F )) deend on the 12

13 angle between K F and P F. Partial comonents of F c determine, e.g., the effective mass m, the secific heat, and the velocities of zero-sound collective modes, and artial comonents of F s determine the sin suscetibility and the roerties of sin-wave excitations. In nonisotroic systems, the relations are a bit more involved, and to obtain m one generally needs to extend the Landau function to the case when one of the momenta is away from the FS [Refs. 2,22]. In this aer we will be esecially interested in the relation between Γ ω αβ,γδ (K,P) and the inverse quasiarticle residue Z = G 1 (,ω)/ (iω m ) ω=0 for on the FS. The quasiarticle Z is not determined within the Landau FL (even in the isotroic case), as it generally does not come from fermions in the immediate vicinity of the FS. Still, there exist two exact relations between Z and charge and sin comonents of Γ ω (K F,P). These relations follow from the Ward identities associated with the conservation laws for the total number of articles (in other words, the charge) and the total sin. These relations do not require Galilean invariance and hold even when Z = Z deends on the osition of F along the FS. The Ward identity associated with the article number conservation identifies G 1 αγ(,ω )/ δµ = δ αγ G 1 (,ω )/ (iω ), where δµ is a small time-deendent and satially homogeneous variation of the chemical otential, with the trile charge vertex Λ c αγ (Ω,K) = Λ c(ω,k)δ αγ, where K = (,ω ) and Ω is set to be infinitesimally small 2,5 : G 1 (,ω) (iω) = Λ c (Ω,K) Ω 0 (2.5) The trile vertex Λ c (Ω,K) is in turn exressed in terms of the fully renormalized charge comonent of Γ ω as shown diagrammatically in Fig. 3. In analytical form Λ c (Ω,K)δ αγ = δ αγ + δ βδ Γ ω αβ,γδ (K,P){ G 2 (P) } d 3 P ω (2π) 3 β,δ = δ αγ (1+2 Γ c (K,P) { G 2 (P) } ) d 3 P ω (2π) 3 (2.6) where d 3 P = d 2 dω. For K = K F, G 1 (,ω)/ (iω) ω=0 = Z, and we have Z = Λ c (Ω,K F ) = 1+2 Γ c (K F,P) { G 2 (P) } d 3 P (2.7) ω (2π) 3 Similarly, the Ward identity associated with the conservation of the total sin identifies G 1 αγ(,ω)/ δ(µ z H ) ω=0 = σ z αγ G 1 (,ω)/ (iω) ω=0 with a trile charge vertex 13

14 FIG. 3: The diagrammatic reresentation of the trile vertex Λ c,s (Ω,K) (K = (,ω )), involved in the Ward identities. For the charge Ward identity (associated with the conservation of the total number of fermions, i.e., the total charge), the bare vertex (the unshaded triangle) is δ α,γ, where α and γ are sin comonents of incoming and outgoing fermion. For the sin Ward identity (associated with the conservation of the total sin), the bare vertex is σα,γ z. Λ s αγ(ω,k) = Λ s (Ω,K)σ z αγ, where now δ(µ z H ) = δ(µ BH z ) is the Zeeman energy shift induced by a small infinitely slowly time-deendent and satially homogeneous variation of a magnetic field H z, which for definiteness we direct along the z axis. Then the same artial derivative G 1 (,ω)/ (iω) is 6 G 1 (,ω) ω=0 = Λ s (Ω,K) Ω 0 (2.8) (iω) The trile vertex Λ s (Ω,K) is exressed in terms of the fully renormalized sin comonent of Γ ω Λ s (Ω,K)σαγ z = σαγ z + σβδγ z ω αβ,γδ(k,p) { G 2 (P) } d 3 P ω (2π) 3 βδ ( = σαγ z 1+2 Γ s (K,P) { G 2 (P) } ) d 3 P ω (2π) 3 For K = K F, this reduces to Z = 1+2 Γ s (K F,P) { G 2 (P) } ω (2.9) d 3 P (2π) 3 (2.10) The fact that the left hand side of (2.7) and (2.10) are identical imlies that the sin and charge comonents of the fully renormalized Γ ω αβ,γδ = Γ cδ αγ δ βδ +Γ s σ αγ σ βδ are related by Γ c (K F,P) { G 2 (P) } d 3 P ω (2π) = Γ 3 s (K F,P) { G 2 (P) } d 3 P (2.11) ω (2π) 3 To the best of our nowledge, this relation has not been exlicitly resented in the literature, although Eqs. (2.7) and (2.10) have been resented in Refs. [2,5] and [6], resectively 14

15 In a generic FL the integration over P in (2.11) is not confined to the FS, and this is the reason why the fermionic Z factor is considered as an inut for Landau FL theory rather than an integral art of it (This is resumably also the reason why the relation (2.11) has not been discussed in the ast.) Lie we said in the Introduction, we will see that tyical P P F in Eq. (2.7) get rogressively smaller as the system aroaches a QCP, and in a CFL regime near a QCP the integrals for Z in (2.7) and (2.10) are redominantly determined by P P F, for which G 2 (P) ω can be aroximated by the quasiarticle art G 2 q (P). In this ω limit, Z becomes anintegral artof FL theory, andeqn. (2.11) establishes thefundamental relation between charge and sin comonents of thevertex function Γ ω between the articles on the FS. The quasiarticle residue and the vertex function Γ ω can both be comuted in direct erturbation theory, and we will erform such calculations near a SDW QCP. Equations (2.7), (2.10), and (2.11) must be satisfied at any order of erturbation theory and can be viewed as consistency checs for erturbative calculations. III. QUASIPARTICLE VERTEX FUNCTION IN AN ORDINARY FERMI LIQ- UID In this Section we resent the derivation of the quasiarticle vertex function Γ ω α,β,γδ (K F,P F ) using the rules alicable to an ordinary FL. This conventional vertex function will serve as a bare quasiarticle function in our subsequent analysis of the vertex function in a CFL. Γ (1) (,) = αβ,γδ α γ β δ - α δ β γ FIG. 4: The vertex function Γ ω to first order in the interaction. The dashed line is the interaction otential U(q) (q = 0 in the first term and q = in the second. 15

16 Consider a system of fermions with some short-range (screened Coulomb) interaction U(q). To first order in U(q), the vertex function is simly the anti-symmetrized interaction (see Fig. 4). Γ ω αβ,γδ (K,P) = U(0)δ αγδ βδ U( )δ αδ δ βγ (3.1) We use the same overall sign convention as in Ref. 13. The first-order result can be re-exressed in terms of sin and charge comonents as Γ ω αβ,γδ(k,p) = ( U(0) U( ) 2 ) δ αγ δ βδ U( ) 2 σ αγ σ βδ (3.2) Γ (2) (,) = x αβ,γδ α β α β α β α β α β γ δ γ δ γ δ γ δ γ δ Γ (2) (,) = x Γ (2) (,) = x αβ,γδ α β α β α β α β α β γ δ γ δ γ δ γ δ γ δ FIG. 5: The diagrams that contribute to Γ ω to second order in the interaction The generic rule how to comute Γ ω αβ,γδ (K,P) beyond first order is that one has to sum u all diagrams excet the ones which contain a article-hole bubble with zero momentum transfer and vanishingly small but finite frequency transfer. The diagrams that contribute to Γ ω to second order in the interaction are shown in Fig. 5, and an examle of a forbidden second-order diagram is shown in Fig. 6. The argument why the forbidden diagrams 16

17 l α γ δ 1 + γ 1 γ 1 l δ β FIG. 6: An examle of a forbidden diagram for Γ ω. This diagram has an internal article-hole bubble with zero momentum transfer, and vanishes for a static interaction (see the text). should not be included is that the corresonding bubble contains dldω l G q (l,ω l )G q (l,ω l +ω q ). (3.3) The integral over Matsubara frequency ω l has to be done first because that integration extends over infinite range. Integrating, we find that it vanishes because the oles in the two roagators are located in the same half-lane of comlex frequency. Γ (3) (,) = αβ,γδ α β γ δ - α β δ - 2x α γ β δ γ - α β δ γ FIG. 7: RPA-tye diagrams for Γ ω to third order in the interaction A. vertex function in RPA Returning to terms of second order in U (Fig. 5) we see that they contain searately article-hole and article-article olarization bubbles. At next (third) order the two get mixed, and there is no controllable way to roceed. A commonly used aroximation for 17

18 a system near a SDW instability is to entirely neglect the renormalization in the articlearticle channel, because by itself such a renormalization does not lead to magnetic order, aroximate the interaction by a constant U, and sum the RPA-tye series of diagrams which contain article-hole olarization bubbles Π( q ) = dldω l G q (l,ω l )G q (l+q,ω l ) at q = F F. The corresonding diagrams, shown in Fig. 7 to third order in U, form a ladder series and can be summed exlicitly. The result is the familiar RPA-tye exression 14 Γ ω,rpa αβ,γδ (K F,P F ) = = U 1 UΠ( F F ) δ αγδ βδ U 2(1 UΠ( F F )) σ αδ σ βγ U 1 U 2 Π 2 ( F F ) δ αδδ βγ U 2(1+UΠ( F F )) δ αδδ βγ (3.4) Using the last line, one may slit Γ ω,rpa αβ,γδ (K F,P F ) into terms containing σ matrices and δ functions ( Γ ω αβ,γδ )sin = U 2(1 UΠ( F F )) σ αδ σ βγ ( Γ ω αβ,γδ )charge = U 2(1+UΠ( F F )) δ αδδ βγ Note, however, that this slitting is not the same as in Eq. (2.4) because the combinations of sin indices in the δδ and σσ terms differ from those in (2.4). We will return to the same notations as in (2.4) later in this section. AStoner-tyemagneticinstabilityoccurswhenU > 0andUΠ( )) = 1foraarticular (later assumed to be q π ), which (for constant U) is determined by the structure of the fermionic disersion. Once ( ) Γ ω αβ,γδ gets enhanced, it is natural to neglect the sin charge comonent of the vertex and aroximate the full anti-symmetrized vertex by its sin comonent, i.e., set Γ ω,rpa αβ,γδ (K F,P F ) = U 2(1 UΠ( F F )) σ αδ σ β,γ (3.5) This Γ ω,rpa αβ,γδ (K F,P F ) can be considered as an effective interaction between fermions at the FS, mediated by collective magnetic excitations. To mae this more transarent, one can exand near q π = (π,π) and extend Eq. (3.5) to fermions not necessarily on the FS (i.e., to non-zero frequencies ω and ω ). At a finite frequency, the article-hole olarization oerator contains a dynamical term which describes Landau daming of a collective boson by interaction with the article-hole continuum. Keeing this term, we find 1 UΠ(K P) ξ 2 + q π 2 +γ ω ω (3.6) 18

19 where γ is the Landau daming coefficient and ξ is the magnetic correlation length, which diverges at the SDW transition. Substituting into (3.5), we obtain where Γ ω,rpa αβ,γδ (K,P) = V(K P) σ αδ σ β,γ (3.7) V(K P) = g ξ 2 + q π 2 +γ ω ω (3.8) can be viewed as the effective four-fermion dynamical interaction mediated by sin fluctuations. The effective couling g in (3.8) is of order U/a 2, where a is the interatomic sacing. The daming coefficient γ also scales with g: γ g/vf 2, because the daming comes from the UΠ term. FIG. 8: Grahical reresentation of the quasiarticle vertex function Γ ω,rpa αβ,γδ (K F,P F ) in the RPA scheme. Note that the combination of sin indices in the two σ matrices is the same as in the anti-symmetrized comonent of the vertex. We resent this Γ ω,rpa αβ,γδ (K,P) grahically in Fig. 8. Note that the combination of sin indices in the two σ matrices in (3.7) is the same as in the anti-symmetrized comonent of the vertex. Re-exressing σ αδ σ βγ via δ and σ matrices involving combinations (αγ) and (βδ), as in (2.4), we find σ αδ σ βγ = δ αδ δ βγ +2δ αγ δ βδ = 3 2 δ αγδ βδ 1 2 σ αγ σ βδ (3.9) and hence Γ ω,rpa αβ,γδ (K,P) = V(K P) ( 3 2 δ αγδ βδ 1 2 σ αγ σ βδ ) (3.10) Placing K and P on the FS, we obtain (K = K F = ( F,0), P = P F = ( F,0)) ( Γ ω,rpa αβ,γδ (K g 3 F,P F ) = ξ 2 + F F q π 2 2 δ αγδ βδ 1 ) 2 σ αγ σ βδ (3.11) 19

20 B. The role of AL diagrams At first glance, Eq. (3.11) is a natural choice for the quasiarticle vertex function in a situation when scattering by sin fluctuations is much more relevant than scattering by charge fluctuations. Uon closer exection, however, we see that this interaction does not satisfy the condition set by Ward identities. Indeed, according to (3.10), the sin and charge comonents of Γ ω have the same deendence on K P through V(K P), but the overall factors differ in sign and magnitude. The Ward identity (2.11), on the other hand, requires that the two must have the same refactors, if they both scale as V(K P). Clearly then, the exression for Γ ω in (3.11) is incomlete and one has to include further contributions for Γ ω which do not fit into the RPA scheme. V(q,Ω) Σ(,ω) =,ωω -q,ω-ω,ωω,ωω FIG. 9: One loo diagram for the fermionic self-energy On hysical grounds, it is natural to use the RPA-tye sin-mediated interaction as a building bloc for constructing further contributions to Γ ω, i.e., exress all non-rpa contributions to Γ ω in terms of RPA-renormalized, sin-mediated interaction Γ ω,rpa rather than the bare interaction U. In this nomenclature, the RPA interaction is the first-order term (one wavy line) and all other terms contain more than one interaction line. At first glance, including higher-order terms cannot resolve the issue osed by Ward identities, as higher-order terms contain higher owers of g, while Ward identities are valid for any couling and hence must hold indeendently at each order in g. However, we show below that in some higher-order terms, extra owers of g come in the combination g/γ, where γ is the rate of Landau daming of collective excitations. The latter is by itself of order g, i.e., the ratio g/γ is of order one. Because of this, certain higher-order terms are actually of the same order in g as Γ ω,rpa. We follow the same strategy as before and do not include terms containing article-hole 20

21 (a) (b) (c) (d) FIG. 10: Thefour two-loo diagrams for thevertex functionγ ω in anordinaryfl. Thedifference with Fig. 5 is that now the interaction line (the wavy line) is the RPA otential V(K P) rather than the bare interaction otential U. To avoid double counting, we didn t include terms which contain Π(K P) as these terms are already included into the RPA renormalization from U to V. Lie before, we did not include the diagram which contains a article-hole bubble with zero momentum and finite frequency transfer. bubbles with strictly zero momentum and vanishing frequency. We also do not include terms which contain article-hole bubbles with transferred 4-momentum K P, as such terms are already incororated into the RPA scheme. The four remaining second-order diagrams are shown in Fig. 10. The first diagram [Fig.(10a)] is of order g 2 without γ in the denominator, and is not a cometitor to the first-order in g term. The second diagram [Fig.(10b)] contains Γ ω,rpa (K P) multilied by the roduct of two Green s functions and one interaction line. Such a term reresents a correction to one of the two trile vertices in the effective interaction. It does contain additional g in the combination g/γ, and, moreover, theinternal integrationyields logξ (Refs. 8,9 ). At thesametime, thesummationover internal sin indices in this diagram shows that the sin structure of the interaction term given by 21

22 diagram in Fig. (10b) is the same σ αδ σ βγ as in (3.7), i.e., this term just renormalizes the overall couling g. The resence of such a term is not unexected because, as we said, there is no hysical small arameter which would distinguish the RPA series. Still, the fact that the diagram in Fig. (10b) reserves the sin structure of Γ ω,rpa imlies that this articular vertex renormalization is irrelevant to the issue of Ward identities. In a more formal way, this vertex renormalization can be made small by extending the theory to N 1 fermionic flavors 8,9,17. The last two diagrams [Figs. (10c) and (10d)] serve our urose in the sense that, on one hand, an additional g aears in the combination g/γ (see Eq.(3.16) below), and, on the other hand, the sin structure of the effective interaction given by any of these diagrams does not match that of Γ ω,rpa. Secifically, the diagram in Fig. (10c) yields a sin structure in the form ( σ αs σ tδ )( σ sγ σ βt ) = 4δ αδ δ βγ +δ αγ δ βδ = 2 σ αγ σ βδ +3δ αγ δ βδ (3.12) s,t and the diagram in Fig. (10d) yields a sin structure of the form ( σ αs σ βt )( σ sγ σ tδ ) = 5δ αγ δ βδ 4δ αδ δ βγ = 2 σ αγ σ βδ +3δ αγ δ βδ (3.13) s,t The diagrams in Figs. (10c) and (10d) have the same structure as the diagrams considered by Azlamazov and Larin(AL) in the theory of the suerconducting fluctuation contribution to the conductivity above T c in layered suerconductors 23,24, and by analogy we call these two diagrams for Γ ω the AL terms. As mentioned in the introduction, these AL terms are generated naturally within the conserving aroximation scheme of Baym and Kadanoff 18, i.e. they are exected to hel conserve article number and sin, and this is indeed what we find. The internal art of each of the two AL diagrams contains the combination of two interactions V 2 (Q) and two Green functions, which for the diagram in Fig. (10c) are G(K + Q)G(P + Q) and for diagram in Fig. (10d) are G(K +Q)G(P Q). For external momenta on the Fermi surface, ω q = ω q. We assume and then verify that tyical q transverse to the FSare small and linearize the disersion ǫ q around F as ǫ q v F, q. In this aroximation, we have G(P Q) = G(P +Q), i.e., the internal arts of the two diagrams in Figs. (10c) and (10d) differ by a minus sign. Multilying Eq. (3.13) by ( 1) and adding to (3.12) we find that the sin structure of the effective interaction from the 22

23 diagrams (10c) and (10d) is 4 σ αγ σ βδ (3.14) This is also a sin-sin interaction, but it contains a different combination of sin indices comared to that in the RPA interaction in Eq. (3.1). Let us now comute the internal art of this diagram. Because we assumed that tyical q aresmall(i.e., muchsmallerthanthefermimomentum), wearoximatethefullg(k) = G(,ω ) by its quasiarticle art. We first tae and to be right at hot sots searated by q π = (π,π) (we recall that a hot sot h is a F oint on the FS in Fig. 1 for which F + q π is also located on the FS). For and at the hot sots, V(K P) = gξ 2. Using 8 v F,+q = v x q x + v y q y, v F,+q+q π = v x q x v y q y, and the fact that, by symmetry, the quasiarticle weight Z h must be the same at each of the hot sots, we exress I c = I d = G(K +Q)G(P +Q)V 2 (Q)d 3 Q/(2π) 3 with Q = (q+q π,ω q ) as I c = g 2 dqx dq y dω q 1 1 (2π) 3 Z 2 h iω q (vx q x +vy q y) iω q (vx q x vy q y) 1 (3.15) (q 2 +ξ 2 +γ ω q ) 2 Because tyical q x and q y in the two Green s functions are of order ω q and tyical q in the bosonic roarator (the interaction V(Q)) are of order ξ 2 +γ ω q, i.e., much larger, one can integrate over momenta in the two fermionic roagators (this integral is ultra-violet convergent) and set q = 0 in the bosonic roagator. The 2D integration over dq x dq y yields π 2 /(Z 2 h v xv y) indeendent of ω q. Integrating over ω q in the bosonic roagator we then obtain I c = gξ 2g 1 γ 4πZ 2 h vxv y The bosonic daming rate has been calculated before 8 and we cite just the result: (3.16) 2 γ = g πz 2 h vxv y 4 = g π(z h vf )2 sinθ (3.17) where θ is the angle between the velocities at the hot sots (sinθ = 2vx v y /(v F )2 ). Substituting into (3.16) we obtain I c = 1 8 gξ2 (3.18) We see that I c has the same structure gξ 2 as the first-order term. To roerly comare the overall factors, we now recall that we need to substitute the two terms into the right hand 23

24 side of the diagrammatic exression for the trile vertex in Fig.3, i.e., comare the diagrams shown in Fig.11. For the interaction terms given by Figs. (10c) and (10d), the internal fermion loo yields an additional 1 and, besides, one needs to sum over all hot regions F which are searated from the external F by (±π,±π). The constraint is that the velocity at should not be antiarallel to that at F, otherwise the integral in (3.15) over one of momentum comonents of q would vanish. A simle exerimentation shows that there are four allowed hot regions of F, i.e., the additional factor for the second-order diagrams of Figs. (10c) and (10d) is 4. FIG. 11: a) The contributions from direct sin-fluctuation exchange and from the two AL diagrams to the vertex function in an ordinary FL. b) The full irreducible vertex function in an ordinary FL (see Eq. (3.20)). Incororating this extra factor into I c, we find that the nominally second-order contributions to Γ ω from the AL diagrams [Figs. (10c) and (10d)] add an extra term to the vertex function for K and P at the hot sots of the form Combining this with Γ ω,rpa = gξ 2 σ αδ σ βγ, we obtain Γ ω,rpa+al αβ,γδ Γ ω,al = 2gξ 2 σ αγ σ βδ (3.19) = Γ ω,rpa αβ,γδ +Γ ω,al αβ,γδ = gξ2 ( σ αδ σ βγ +2 σ αγ σ βδ ) = 3 2 gξ2 (δ αγ δ βδ + σ αγ σ βδ ) (3.20) Weseethatnowsinandchargecomonentsareofequalsignandmagnitude, arecondition for the Ward identities to be satisfied. 24

25 The analysis of Figs. (10c) and (10d) can be extended to K and P away from hot sots. Performing the integration in the same way as before we find, to leading order in ξ, I c = (1/8)V(K P)andΓ ω,al = 2V(K P) σ αγ σ βδ. ThetotalΓ ω,rpa+al = Γ ω,rpa +Γ ω,al is given by where Γ ω,rpa+al αβ,γδ (K,P) = V eff (K P) (3.21) V eff (K P) = V(K P)( σ αδ σ βγ +2 σ αγ σ βδ ) = 3 2 V(K P)(δ αγδ βδ + σ αγ σ βδ ) (3.22) This vertex function has equal sin and charge comonents and obviously satisfies the relation (2.11) imosed by the Ward identities. We emhasize, however, that the equivalence between the comonents Γ c (K,P) and Γ s (K,P) holds only as long as both have the same deendence on K P (given, in our case, by V(K P)). This is true only when the system is sufficiently close to a magnetic transition, such that the charge comonent of the RPA interaction can be neglected, together with non-rpa contributions. In a generic FL, Γ c (K,P) and Γ s (K,P) have different functional forms and do not have to be equal, although Eq. (2.11) must, indeed, hold. The imortance of AL terms for the roer descrition of a FL has been emhasized earlierinvariouscontexts: intheanalysisoftheeffectoffluctuationsnearasuerconducting 23,24 or a ferromagnetic transition 19, or also at the disorder driven metal-insulator transition in disordered systems 20. That these terms must generally be included into the irreducible vertex function within the conserving aroximation scheme can also be seen by exlicitly differentiating the diagrammatic exression for G 1 (,ω) = ω v F ( F )+Σ(,ω). The argument that the sin-mediated interaction is a building bloc for the diagrammatic exansion imlies that the fermionic self-energy can be exressed via V(K P). The one-loo self-energy diagram is shown in Fig. 9. A variation of this self-energy caused by an external erturbation (either a wealy time-deendent comonent of the chemical otential or a wealy time-deendent magnetic field, both homogeneous in sace) generates two arts one comes from the variation of the internal fermionic line in Fig. 9 and the other is generated by the variation of a Green s function within the RPA interaction line 24. By varying the internal fermionic line, we reroduce Fig. 3 with Γ ω = Γ ω,rpa. To see the effect of the variation of the interaction line, we recall that Γ ω,rpa is obtained by summing u article-hole bubbles; using Eq. (3.5) for Γ ω,rpa, and varying Π(K P) we see that 25

26 contributions involving two fluctuation roagators are thereby generated. We show this rocedure grahically in Fig. 12. Varying each one of the two Green s functions maing u Π(K P), we obtain two additional terms for Γ ω, which are exactly the two AL diagrams. FIG. 12: The relation between the derivative of the self-energy with resect to a time-deendent correction to the chemical otential δµ(t) and the direct and AL terms in the article-hole irreducible vertex function. The cross shows which Green s function is varied. The AL terms are generated by varying the bosonic roagator V(q) which contains article-hole bubbles with transferred momentum q q π. 1. The accuracy of the calculation of Γ ω in an ordinary FL and the crossover to a critical FL The consideration above is essentially that of the conserving aroximation scheme 18,25 in the sense that the guiding rincile which forced us to loo at additional terms for Γ ω beyond Γ ω,rpa is the set of charge and sin Ward identities associated with the conservation of the total number of articles and total sin. In terms of a erturbative exansion, in the derivation of (3.21), we collected terms of order g and neglected the term of order g 2 and also the term of order glogξ, which had the same sin structure as the RPA interaction. Consider first the g 2 term [Fig. (10a)]. A simle analysis shows that in 2D the extra 26