arxiv:cond-mat/ v1 [cond-mat.stat-mech] 7 Sep 2004

Transcription

1 Critica behavior of weaky interacting bosons: A functiona renormaization group approach arxiv:cond-mat/967v cond-mat.stat-mech 7 Sep Nis Hassemann, Sascha Ledowski, and Peter opietz Institut für Theoretische Physik, Universität Frankfurt, Robert-Mayer-Strasse 8, 65 Frankfurt, Germany (Dated: October, 7) We present a detaied investigation of the momentum-dependent sef-energy Σ(k) at zero frequency of weaky interacting bosons at the critica temperature T c of Bose-Einstein condensation in dimensions 3 D <. Appying the functiona renormaization group, we cacuate the universa scaing function for the sef-energy at zero frequency but at a wave vectors within an approximation which truncates the fow equations of the irreducibe vertices at the four-point eve. The sef-energy interpoates between the critica regime k k c and the short-waveength regime k k c, k c is the crossover scae. In the critica regime, the sef-energy correcty approaches the asymptotic behavior Σ(k) k η, and in the short-waveength regime the behavior is Σ(k) k (D 3) in D > 3. In D = 3, we recover the ogarithmic divergence Σ(k) n(k/k c) encountered in perturbation theory. Our approach yieds the crossover scae k c as we as a reasonabe estimate for the critica exponent η in D = 3. From our scaing function we find for the interaction-induced shift in T c in three dimensions, T c/t c =.3an /3, a is the s-wave scattering ength and n is the density, in exceent agreement with other approaches. We aso discuss the fow of margina parameters in D = 3 and extend our truncation scheme of the renormaization group equations by incuding the six- and eight-point vertex, which yieds an improved estimate for the anomaous dimension η.53. We further cacuate the constant im k Σ(k)/k η and find good agreement with recent Monte-Caro data. PACS numbers: 3.75.Hh, 5.3.Jp, 5.7.Jk I. INTRODUCTION The physics of weaky interacting bosons has seen a reviva thanks to the improved experimenta technique of atom traps which aow for a detaied study of Bose- Einstein condensation in a controed environment. Besides the theoretica effort invested into studies of harmonicay trapped bosons (for a review, see Ref. ), this aso generated renewed interest into the behavior of homogeneous Bose gases and ed to some new insights ). Many aspects of weaky interacting Bose gases have been understood for quite some time now, such as the universaity cass O() of the condensation transition aong with an accurate knowedge of critica exponents. However, the universaity of weaky interacting bosons at the critica point is not imited to critica exponents since a weaky interacting Bose gas behaves universa at a ength scaes arger than the therma de-brogie waveength λ th 3,. The origin of this extended universaity is that at arge ength scaes quantum fuctuations become unimportant so that this regime of the Bose gas is competey described by a cassica φ -mode 5. In this work we present a detaied study of the momentumdependence of the sef-energy at the critica point of Bose- Einstein condensation at zero frequency, using the functiona renormaization group formaism in the form introduced by Wetterich 6 and by Morris 7. Some resuts of this manuscript were aready presented in a brief form 8. Here we give a detaied account of the cacuation and further incude an extensive treatment of margina terms. In the imit of weak interactions, parametrized by the s-wave scattering ength a, the sef-energy is universa not ony in the imit of sma wave vectors, k, it is proportiona to k η with a finite anomaous dimension η. It remains universa, in the sense that it can be written in a scaing form independent of a, up to momenta which ony need to be sma compared to λ th. At some crossover scae k c the zero frequency sef-energy eaves the anomaous scaing regime and enters the perturbative regime, the sef-energy correction to the dispersion becomes negigibe compared with the bare dispersion ǫ k = ρ k, with ρ = /m. (.) Here m is the bare mass. At finite temperatures and D <, a perturbative cacuation of the sef-energy is utravioet (UV) divergent. This divergence can be addressed by using the therma de-brogie ength as an UV cutoff. In D = 3, however, the perturbative regime k c < k < π/λ th remains non-trivia even in presence of a UV cutoff, since additiona ogarithmic infra-red (IR) divergences appear. As discussed in detai by Baym et a. 5, the IR divergence cannot be treated in an ad-hoc manner by introducing an IR cutoff, since doing so introduces an artificia additiona scae which directy enters quantities which shoud be universa, such as the interaction induced shift of the critica temperature. The IR divergence can be removed by a re-summation using a variety of standard many-body techniques, e. g. bubbeor adder-summation and/or sef-consistent approaches 5. However, these methods are uncontroed in the critica regime one is faced with a strong-couping probem. Renormaization group (RG) techniques are expected to perform better and severa authors have ap-

2 pied RG techniques to investigate the IR behavior of weaky interacting bosons 9,,, though no attempt was made to cacuate the momentum dependence of the sef-energy. Note that standard fied theoretica RG is confined to the critica regime k k c. In fact, even the scae k c cannot be obtained within such an approach. To interpoate between the critica and the short waveength regime, functiona RG techniques 6, 7 are a natura choice, since they track the fow of compete vertex functions rather than just a sma number of couping parameters. Note that the functiona RG approach incudes a priori aso terms which are irreevant according to their scaing behavior, but which are important for a correct description of the perturbative k c k regime. In this work we present in detai a functiona RG approach to this probem and cacuate the universa sef-energy to eading order in the sma parameter an /3 for a wave vectors (n is the boson density). Since the critica behavior at T c is cassica, we wi focus ony on the sef-energy Σ(k,iω n = ) at vanishing Matsubara frequency. We begin with the genera form of the functiona RG fow equations for the irreducibe vertices up to the fourpoint vertex in Sec. II, we sti retain a frequencies. Our approach is based on the sharp-cutoff version of the functiona RG 6, 7,, 3. In Sec. III, we turn to the effective cassica fied theory and rewrite the probem in a notation appropriate to the cassica imit. In Sec. III A the reevant and margina parameters are cassified in dimensions 3 < D < and their functiona fow equations are stated in Sec. IIIB. Sec. IV contains the centra part of our work, the cacuation of the sef-energy. We first re-examine perturbation theory and its divergence in D = 3 before we turn to the cacuation of the sef-energy within the functiona RG formaism. This formaism aows us to cacuate the scaing function σ(x) characterizing the zero frequency sef-energy at the critica point. We define σ(x) by σ(x) = (ρ k c ) Σ(k c x) Σ(), (.) such that the crossover occurs at x =. Here Σ(k) = Σ(k,iω n = ) is the exact zero frequency sef-energy at the critica temperature. Standard fied theoretica RG can ony describe the asymptotic imit x, σ(x) x η. Our approach is based on a truncation of the exact hierarchy of functiona RG fow equations at the four-point vertex, i. e. we ignore six-point and higher order vertices. Whie in this approximation the fow of margina terms is not consistenty described, the resuting approximation for the four-point vertex does in fact incude margina and infinitey many irreevant terms. The resuting fow equation for the sef-energy can then be soved and we derive the compete momentum dependence for the sefenergy in 3 D <, see Eqs. (.5,.5) beow, which constitute the centra resut of this work. A numerica evauation of Eqs. (.5,.5) in D = 3 is aso presented. We use σ(x) to cacuate the shift of the critica temperature T c of the condensation transition in Sec. V. In D = 3 we obtain T c /T c =.3 an /3 to owest order in an /3, in good agreement with recent numerica investigations, 5 and other anaytica resuts 6, 7 (for a recent review on this topic see ). In Sec. VI, we improve upon the truncation of the fow equation and account for the couping parameters which become margina in D = 3. Three additiona parameters must be taken into account, two of which are associated with the inear momentum dependence of the fourpoint vertex and one which describes the momentumindependent part of the six-point vertex. We discuss different truncation schemes of the fow equations and show that the incusion of margina terms provide an improvement for the fixed point vaue of the anomaous dimension. In Sec. VII we cacuate σ(x) in D = 3 incuding the margina terms of the four-point vertex but ignoring irreevant terms. We demonstrate how this truncation fais in the argex regime, it predicts incorrecty σ(x) x which in turn woud predict T c /T c an /3 n(an /3 ). On the other hand, in the critica regime irreevant terms ony ead to a renormaization of margina and reevant ones and σ(x) is we described by a theory were irreevant terms are not incuded. We use this approach to express the prefactor A 3 of the anomaous scaing term, σ(x) A 3 x η, as a function of η for D = 3. Our resut is in good agreement with recent Monte Caro resuts. Finay, in Sec. VIII we summarize and concude this work. II. FUNCTIONAL RG FLOW EQUATIONS FOR BOSONS Our starting point is a standard effective action describing free bosons with a two-partice interaction, which is oca at the bare eve. The action with an UV cutoff Λ (to be specified ater) is of the form S Λ { ψ,ψ} = S Λ { ψ,ψ}+s int Λ { ψ,ψ}, (.) the non-interacting part is given by SΛ { ψ,ψ} = Θ(Λ k ) iω n +ǫ k µ+σ(,i) ψ ψ, (.) with Θ(x > ) = and Θ(x < ) =. The interaction part is given by SΛ int { ψ,ψ} = Θ(Λ k )Σ Λ () Σ(,i) ψ ψ + (!) δ +,+ Γ () Λ (, ;, ) ψ ψ ψ ψ +..., (.3) the eipsis denotes three-body and higher order interactions,whichweignoreatthebareeve. Hereψ isa

3 3 compex bosonicfied. We use the notation = (k,iω n ), = (βv) k,ω n, and δ = βvδ kk δ ωnω n, β is the inverse temperature, V is the voume, and ω n = πnt are bosonic Matsubara frequencies. In Eq. (.) we have incuded the exact sef-energy at vanishing momenta and frequencies as a counterterm in the definition of the free action. Throughout this work, we sha work at temperatures T T c such that the U() symmetry is not broken. The generating functiona of the one-partice irreducibe n-point vertices Γ (n) Λ (,..., n, n,..., ) of the theory with IR cutoff Λ can be expanded in terms of the fieds φ = ψ as foows Γ Λ { φ,φ} = n= (n!) n i,j= i j δ n,n+...+ Γ (n) Λ ({ i, j }) φ φ n φ n φ. (.) The functiona RG fow equations of the first few irreducibe vertices (up to the six-point vertex) for nonreativistic fermionic and bosonic many-body systems can be found in Ref. 3. We summarize beow the fow equations reevant for bosons up to the four-point vertex. The fow equations for the six-point vertex and some terms of the fow of the eight-point vertex are given in Appendix A. A. Free energy For competeness, we ist here the fow equation for the free energy Γ () Λ, athough beow we sha not discuss it further. The fow equation is given by Λ Γ () Λ = V δ(λ k ) iω n ǫ k +µ Σ(,i) n iω n ǫ k +µ Σ(,i) Γ () Λ (), (.5) Γ () Λ () is the two-point vertex, defined by Γ () Λ () = Σ Λ() Σ(,i). (.6) B. Two-point vertex The main interest of the work is the fow equation of thetwopointvertex,whichwewireateinsec.ivtothe zero energy scaing function of the sef-energy, Eq. (.). The fow equation for the two-point vertex is given by Λ Γ () Λ () = Ġ Λ ( )Γ () Λ (, ;,), (.7) δ(λ k ) Ġ Λ () = iω n ǫ k +µ Σ(,i) Γ () Λ () (.8) is the cutoff dependent singe scae propagator, with support ony at k = Λ. As is evident from Eq. (.7), the fow of the two point vertex depends on the properties of the four-point vertex at finite wave vectors and we thus need to derive the momentum-dependent fow of the fourpoint vertex. C. Four-point vertex The fow equation for the irreducibe four-point vertex is responsibe for the cross-over from the weak couping regime at sma scaes to the critica regime at arge scaes. The fow is given by Λ Γ () Λ (, ;, ) = ĠΛ ()G Λ ( )Γ () Λ (, ;,)Γ () Ġ Λ ()Γ (6) Λ (,,;,, ) Λ (,,, ) ĠΛ ()G Λ ( )+G Λ ()ĠΛ( ) Γ () Λ (, ;, )Γ () = + Λ (,;, ) ĠΛ ()G Λ ( )+G Λ ()ĠΛ( ) Γ () Λ (, ;, )Γ () Λ (,;, ) = +, (.9) = + Γ (6) Λ is the six-point vertex and Θ(Λ < k < Λ ) G Λ () = iω n ǫ k +µ Σ(,i) Γ () Λ () (.) is the cutoff reguarised propagator. The notation Θ(Λ < k < Λ ) is shorthand for Θ( k Λ) Θ( k Λ ). This equation is shown graphicay in Fig.. For the most part of this work, we sha in fact ignore the contribution

4 6 ZS BCS of the cassica sector by finite Matsubara frequencies atogether,sincethisonyeadstocorrectionsint c oforder a. It is convenient to write the singe-partice Green s function at zero Matsubara frequency in the foowing scaing form G Λ (k,i) = Z ρ Λ G (kλ), (3.3) FIG. : Diagrammatic representation of the fow equation for the four-point vertex, see Eq. (.9). We adapt the notation BCS, ZS and ZS for the diagrams from the usua fermionic anguage 8, even though this notation does not impy a physica correspondence. from the six-point vertex to this fow and work soey with the two- and four-point vertex. To cacuate fixedpoint properties, we wi however aso incude the sixpoint vertex in Sec. VI. III. EFFECTIVE CLASSICAL FIELD THEORY To discuss the cassica critica behavior, it is sufficient to retain ony the zero Matsubara frequency part of a vertices. In principe, a fieds with non-zero Matsubara frequencies can be integrated out using the fow equations given above, which eads to finite renormaizations of the parameters appearing in the effective cassica theory 9. The effective UV cutoff of the cassica theory is determined by the therma de Brogie waveength, Λ = π λ th, λ th = h πmt = ZS π mt. (3.) In D = 3 the initia vaue of the four-point vertex can be parameterized in terms of the two-body scattering ength a, Γ () Λ (,;,) = 8π a m = 6πρ a. (3.) A. Functiona RG fow equations for the rescaed vertices We assume that the finite renormaizations due to nonzero Matsubara frequencies are impicity taken into account via the initia conditions at scae Λ = Λ of the effective cassica theory. As discussed in detai in Ref. 5, to cacuate the inear shift of the critica temperature (see Sec. V), one can in fact ignore the renormaization G (q) = Θ( < q < e ) R (q), (3.) were the minus sign is introduced in Eq. (3.3) to arrive at the usua definition of the cassica Green s function. Here = n(λ/λ ) is the ogarithmic fow parameter, and the inverse dimensioness propagator is R (q) = Z q + Γ () (q)+ Z ρ ΛΣ(,i) µ, (3.5) with the dimensioness irreducibe two-point vertex Γ () (q) = Z ρ Λ Γ() Λ (Λq,i) = Z ρ Λ Σ Λ(Λq,i) Σ(,i). (3.6) The cassica wave-function normaization factor Z is given by () Γ (q) Z = q. (3.7) q = Simiar to Eq. (3.6), we define the cassica dimensioness higher order vertices for n by Γ (n) (q,...,q n;q n,...,q ) = ( D T) n Λ D(n ) n (Z /ρ ) n Γ (n) Λ ({k i = Λq i,ω ni = }). (3.8) For ater convenience we have incuded the numerica factor D = Ω D /(π) D in the definition of the vertices, Ω D = π D/ /Γ(D/) is the surface area of the D-dimensiona unit sphere. The rescaed vertices satisfy functiona fow equations of the form Γ(n) ({q i }) = n D(n ) nη n i= q i qi Γ (n) ({q i }) + Γ (n) ({q i }), (3.9) η = nz (3.) is the fowing anomaous dimension. In particuar, the rescaed two-point vertex satisfies Γ() (q) = η q q Γ () (q)+ Γ () (q), (3.)

5 with Γ () (q) = q Ġ (q ) Γ () (q,q ;q,q), (3.) q = d D q Ω D, (3.3) and Ġ (q) = δ( q ) R (q) 5. (3.) The functiona fow equation for the dimensioness fourpoint vertex foows from Eq. (3.9). It invoves the inhomogeneity Γ () (q,q ;q,q ) = Ġ (q) Γ (6) (q,q,q;q,q,q ) Ġ (q) G (q ) Γ () (q,q ;q,q) Γ () (q,q ;q,q ) q q q =q +q q Ġ (q) G (q )+ G (q)ġ(q ) Γ() (q,q ;q,q ) Γ () (q,q;q,q ) q q =q q +q Ġ (q) G (q )+ G (q)ġ(q ) Γ() (q,q ;q,q ) Γ () (q,q;q,q ). (3.5) q =q q +q q The fow equations for the six-point and eight-point vertices aso have the form Eq. (3.9) Γ (6) (8) and Γ invove various combinations of the two-, six-, four-, eightand ten-point vertices. A terms entering the inhomogeneity of the six-point vertex and the terms needed for a cacuation up to second order in the reevant and margina parameters of the inhomogeneity of the eightpoint can be found in Appendix A, see Eqs. (A,A3). with scaing dimension ǫ = D. The exact fow equations of these parameters are r = ( η )r + u = ( D η )u + Γ () (), (3.8) Γ () (,;,). (3.9) B. Cassification of couping parameters Athough we are utimatey interested in the fow of vertex functions, it is usefu to first consider the fow of margina and reevant couping parameters, since irreevant parameters become oca functions of the reevant and margina ones at the fixed point. To propery organize the fow of irreevant terms it is thus necessary to know the fow of the reevant and margina parameters. We first investigate here the reevant and margina terms for 3 < D <, ignoring additiona margina terms in D = 3.. Reevant couping parameters In 3 D < there are two reevant couping parameters, r = () Z Γ () = im q ρ Λ Σ Λ(Λq,i) Σ(,i), (3.6) with scaing dimension +, and u = Γ () (,;,), (3.7). Margina couping parameter For 3 < D < the ony margina parameter is the wave-function renormaization Z. The exact fow equation is Z = η Z, (3.) the fowing anomaous dimension is and η = = (q) q Γ () R () q = q Γ() (q,ˆq ;ˆq,q) ˆq q =, (3.)... ˆq = d D q Ω D δ( q )... (3.) denotes the integra over the unit sphere. The surface area Ω D of the unit sphere was defined beow Eq. (3.8). In D = 3, additiona margina couping parameters appear which we ignore in this section. We wi discuss them in Sec. VI.

6 6 3. The Wison-Fisher fixed point cose to D = It is instructive to take a coser ook at the perturbative one-oop RG fow equations for the reevant couping parameters. If we retain from the two-point vertex ony the reevant and margina part, we may approximate u u * R (q) r +q, (3.3) so that G (q) Θ( < q < e ) r +q, (3.) Ġ (q) δ( q ) r +q. (3.5) If in addition we retain ony the reevant part of the four-point vertex, we may repace it by a momentumindependent constant, Γ () (q,q () ;q,q ) Γ (,;,) u. (3.6) Within this approximation we obtain from Eq. (3.), and Eq. (3.5) reduces to Γ () Γ () () u +r, (3.7) (,;,) 5 (+r ). (3.8) Because the anomaous dimension is reated to the momentum dependence of the four-point vertex which is irreevant in D > 3, we may set η within this approximation. Furthermore, on the critica trajectory r = O(u ) can aso be ignored in Eq. (3.8) as ong as u remains sma compared with unity. Then we obtain, with ǫ = D, u u = ǫu 5 u, (3.9) r = r +u. (3.3) Eq. (3.9) is easiy soved. The soution can be written in the form u = u e ǫ( c) + = e ǫ(c ) +, (3.3) u = 5 ǫ (3.3) is the vaue of u at the RG fixed point, and the ogarithmic crossover scae is c = ( ) ǫ n u ( ) u ǫ n u. (3.33) u FIG. : Typica fow of the couping u as a function of for u u and D = 3. Throughout this work we assume that u u, corresponding to weak bare interactions between the bosons. Note that the right-hand side of Eq. (3.3) is expressed in terms of the Fermi function. If we think of ǫ as the inverse temperature and c as the chemica potentia, the quaitative behavior of u /u is cear: within a narrow (on the scae c ) interva of width /ǫ centered at = c the ratio u /u raises from the sma vaue e ǫ(c ) to a vaue cose to unity, see Fig.. Given the soution u in Eq. (3.3), the RG equation (3.3) for r is easiy soved, r = ǫ 5 +e r + ǫ 5 ǫ 5 c e d e ǫ( c) +. (3.3) To obtain a fixed point at, the intia vaue r has to be fine tuned such that r = ǫ 5 + ǫ 5 Then we may write r = r +u e d e ǫ( c) +. (3.35) d e ( ) e ǫ( c) +, (3.36) r = u = ǫ 5. (3.37) In Sec. IV we show that, up to a numerica factor of the order of unity, the momentum scae k c = Λ e c associated with the ogarithmic scae factor c, i. e. /ǫ ( ) /ǫ u u k c = Λ Λ, (3.38) u can be identified with the crossover scae the critica k η -form of the energy dispersion begins to emerge. u

7 7 IV. THE SELF-ENERGY AT THE CRITICAL POINT OF BOSE-EINSTEIN CONDENSATION A. Second order perturbation theory To begin with, et us attempt to cacuate the sefenergy by means of straightforward second order perturbation theory5, which yieds in a continuum mode with UV cutoff Λ Σ(k) Σ() = T 3 d D p Γ (π) Dχ (p) (.) Θ(Λ p+k ) ǫ p+k Θ(Λ p ) ǫ p, Γ = Γ () Λ (,;,) is the bare vertex, and d D p Θ(Λ p )Θ(Λ p +p ) χ (p) = T (π) D. ǫ p ǫ p +p (.) With ǫ p = ρ p we obtain for Λ, χ (p) = D D T ρ (.3) p D D = dxx( x) (D )/ dy yd (y +) πγ( D = 3 D ) π (D ) Γ( D ) sin( π. (.) (D ) Note that D is finite for < D < ; in particuar, 3 = π /. Substituting Eq. (.3) into Eq. (.) and taking the imit Λ, we can scae out the k-dependence and obtain for D < 3 <, using the definitions u = D TΛ ǫ ρ Γ and k c = Λ (u /u ) /ǫ, ρ Σ(k) Σ() = B Dk ǫ c k ǫ, (.5) after an integration by parts the coefficient B D can be cast into the foowing form, I D = Ω D Ω D B D = d dx π 3u (D 3) D I D, (.6) dϑ(sinϑ) D ( +xcosϑ +xcosϑ+x dxx D 6 ). (.7) In Fig. 3 we show a numerica evauation of B D for 3 < D <. With I 3 = and 3 = π / we obtain for D 3 B D FIG. 3: Graph of the coefficient B D defined in Eqs. (.6) and (.7) in dimensions 3 < D. Perturbation theory is vaid as ong as the correction(.5) issmacomparedwith ρ k, i.e.fork k c B /ǫ D. Dueto the IR divergence of B D for D 3, cose to three dimensions perturbation theory breaks down even outside the critica regime. In a IR reguarized approach with a cutoff ξ, the divergence in Eq. (.5) is repaced by a n(kξ) behavior 5. Beow we sha see that the functiona RG approach shows a n(k/k c ) behavior, as expected for a theory the divergence due to density fuctuations is screened within a non-perturbative treatment. B. Functiona RG cacuation of the sef-energy We now derive the behavior of the sef-energy microscopicay using the functiona RG equations for the effective cassica fied theory.. The u -truncation for the Γ () D inhomogeneity To begin with, we need the RG fow of the fourpoint vertex. In the simpest approximation, we expand Γ () (q,q ;q,q ) in powers of u, which shoud be accurate as ong as the renormaized couping remains sma. To eading order we simpy ignore the six-point vertex and repace the four-point vertices on the right-hand side of Eq. (3.5) by their reevant part u. In this approximation, Γ () (q,q ;q,q ) u χ ( q +q ) + χ ( q q )+ χ ( q q ), (.9) B D b 3 D 3, b 3 = 3π u. (.8) 6 the generaized susceptibiity is given by

8 8 χ (q) = Ġ( q ) G ( q +q ) q Θ( < ˆq +q < e ) +r ˆq +q +r ˆq Ω π D = dϑ(sinϑ) D +r Ω D Θ( < q +qcosϑ < e ) +r +q, (.) +qcosϑ in the second and third ine we have used the approximations (3.) and (3.5), taking into account ony the reevant and margina part of the two-point vertex.. Generaized susceptibiity It is instructive to examine the behavior of the generaized susceptibiity χ (q). For e the asymptotic behavior for sma and arge q can be easiy obtained anayticay, χ (q) { (+r ) +O(q) for q e, (+r )q for q e. (.) The important point is that the eading correction for sma q is inear in q: expanding Eq. (.) in powers of q, we find and χ (q) = χ ()+ χ ()q +O(q ), (.) χ () = χ () = (+r ) (+r ), (.3) S (+r ) S, (.) S = δ(ˆq ˆq) ˆq = Ω D Ω D, (.5) S = Θ(ˆq ˆq)ˆq ˆq ˆq = S D. (.6) In D = 3 we obtain with S = / and S = / ( χ () = (+r ) ) +r, D = 3. (.7) For ater reference, et us give the exact function χ (q) in D = 3, the anguar integration in Eq. (.) can easiy be performed anayticay. The resut can be written as χ (q) = Θ(q q ) (+r )q n +r +q +qq +r +q, (.8) +qq { if q > q =, (.9) q/ if q < { if e > q q = e q q if. (.) e < q In particuar, for q min{,e } we have (+q) χ (q) = (+r )q n +r. (.) +r Thus, for sma momenta our approximation (.9) yieds in D = 3, Γ () (q,q ;q,q ) 5 (+r ) ( + u (+r ) ) +r q +q + q q + q q +O(q i ). (.) The term inear in q generates margina parameters, as mentioned in Sec.IIIB, when we iterate the RG. The inear term exists for a D, but the corresponding couping parameters are irreevant in D > 3. Because these couping parameters are not consistenty taken into account in the u -truncation given in Eq. (.9), we cannot expect that this truncation gives numericay accurate resuts in the critica regime cose to three dimensions. On the other hand, in the short-waveength regime (and aso for D 3 η), our u -truncation (.9) is sufficient. In this section we sha therefore proceed with this approximation, which produces we-defined resuts even in D = 3 provided the anomaous exponent η is cacuated sef-consistenty by soving an integra equation, see Eq. (.9) beow. In Sec. VI we sha improve on this approximation by expicity incuding a margina couping parameters in the critica regime in D = Expicit expression for the four-point vertex () With Eq. (.9) as an approximation for Γ, we are now in a position to cacuate the fow of the four-point vertex. Let us rewrite Eq. (.9) in the foowing way, Γ () (q,q ;q,q ) 5u (+r ) + Γ (mi) (q,q ;q,q ), (.3) u

9 9 Γ (mi) (q,q ;q,q ) = u χ(mi) ( q +q ) with + χ (mi) ( q q )+ χ (mi) ( q q ), (.) χ (mi) (q) = χ (q) χ () (.5) describing the momentum-dependent part. The superscript (mi) indicates that these terms are margina or irreevant. In this approximation the four-point vertex takes the form Γ () (q,q ;q,q ) u + Γ (mi) (q,q ;q,q ), (.6) with the irreevant (and in D = 3 aso margina) parts given by Γ (mi) (q,q ;q,q ) = = = d e ǫ( ) dτη τ Γ(mi) (e ( ) q,e ( ) q ;e ( ) q,e ( ) q ) d e ǫ( ) dτη τ u χ(mi) (e ( ) q +q )+ χ (mi) (e ( ) q q )+ χ(mi) (e ( ) q q ) dte ǫt t dτητ u t χ(mi) t (e t q +q )+ χ (mi) t (e t q q )+ χ(mi) t (e t q q ). (.7) C. Anomaous dimension To obtain the anomaous dimension η = Γ () (q)/ q () q = we cacuate Γ (q) via Eq. (3.). With the approximation Eqs. (.6,.7) for the four-point vertex we find Γ () Γ() (q) = (q,ˆq ;ˆq,q) R () ˆq u 3 dte ǫt t dτητ u t +r (+r ) (mi) χ t (e t ( ˆq +q ) ˆq. (.8) From Eq. (3.) we then obtain an integra equation for the fowing anomaous dimension η = dt(,t)u t e t dτητ, (.9) with (,t) = 3eǫt (+r ) q χ t (e t ( ˆq q= +q )) ˆq 3 = (D )e (D 3)t χ D(+r ) t(e t ) +e (D )t χ t (e t ), (.3) χ (x) = d χ (x)/dx and χ (x) = d χ (x)/dx. Note that on the right-hand side we have repaced χ (mi) t (q) χ t(q), because the constant part χ t () does not contribute to the derivative. The integra equation (.9) together with the one-oop fow equations for the reevant couping parameters, r = ( η )r + u +r, (.3).. -. r u η / c FIG. : Sef-consistent RG fows of u, r and η at the critica point in D = 3 (dashed ines), obtained numericay from the couped intego-differentia Eqs. (.9,.3,.3). For comparison, the Fermi-function approximation Eq.(.33) for the fows is aso shown (fu ines), with fixed point vaues matching the numerica soution. u = (ǫ η )u 5 (+r ), (.3) form a system of three couped integro-differentia equations for the three unknown functions r, u, and η. A numerica soution of these equations for D = 3 is shown in Fig.. To eading order in ǫ (i. e. cose to four dimensions) Eqs. (.9,.3,.3) decoupe and one recovers the fixed point vaues u = ǫ/5, r = ǫ/5 from Eq. (3.37). To this order it is consistent to ignore η and r on u

10 the right-hand side of Eq. (.3), so that it reduces to Eq. (3.9), with the Fermi function soution for u given in Eq. (3.3). From the numerica soution of Eqs. (.9.3) for D = 3 we find that even if r and η are taken into account, to a good approximation the quaitative behavior of u is sti we described by a Fermi function if we adjust the fixed point vaue u to the one obtained from the numerica soution. Furthermore, Fig. shows that aso η and r roughy foow the same functiona form and we therefore approximate η η r r u u e ǫ(c ) +, (.33) with the fixed point vaues taken from the numerica soution and η = im η is the anomaous dimension at the fixed point. Athough the functiona form (.33) for η is quaitativey incorrect for c ( η u, as Eq. (.33) predicts η u ) we may use Eq. (.33) as a zeroth approximation on the right-hand side of the integra equation (.9). The point is that for c the right-hand side of Eq. (.9) is not sensitive to the precise vaue of η in this regime, and the first iteration of the integra equation yieds an accurate resut even for sma. Let us now cacuate the fixed point vaues. We obtain from Eq. (.9) η = u dt(,t)e ηt, (.3) we use the notation η = η. The function (,t) is obtained from (,t) defined in Eq.(.3) by repacing the functions χ t (q) on the right-hand side by χ (q) im χ (q) In particuar, in D = 3 we have χ (q) = Θ( ˆq +q ) +r ˆq +q +r ˆq. (.35) ( (+q) +r ) Θ( q)n (+r )q +r ( (+q) +r ) +Θ(q )n ( q). (.36) +r Note that for arge t the kerne (,t) vanishes as e (D 3)t, so that in D = 3 it is crucia to retain the η-dependence on the right-hand side of Eq. (.3). This is cosey reated to the appearance of margina terms of the four-point vertex in D = 3. The sef-consistent fixed point vaues can be obtained by using the fact that, at the fixed point, we can express u and r as a function of η, r = 3 u = 9 η 3η, (.37a) ( η)( η) ( 3η). (.37b) Inserting these expressions for u and r into Eq. (.3), we obtain a sef-consistent equation for η which can be soved numericay. The numerica soution gives r.3, u.3, η.. (.38a) (.38b) (.38c) The above vaue for η is approximatey three times as arge as the generay accepted vaue η.38 obtained by severa different techniques, 3. However, given the simpicity of our truncation, it is quite satisfactory that our estimate (.38c) for η has the correct order of magnitude. In any case, our simpe u -truncation certainy gives a much better estimate for η than the sefconsistent two-oop cacuation empoyed by Baym et a. 5, which predicts a vaue of.5 for η. In Sec.VI we sha further improve on the u -truncation empoyed here by expicity taking the RG fow of margina couping parameters into account. For competeness, we show how to obtain the standard resut for η to second order in an ǫ-expansion. We can rewrite Eq. (.3) as 3u η = D(+r ) 3u = D(+r ) dλλ ǫ+η (D ) χ (λ)+λ χ (λ) ( η) dλλ ǫ+η χ (λ)+ χ (). For sma ǫ = D we obtain to eading order (.39) η = 3u 6 ( χ() χ())+ χ (), (.) the function χ(q) is obtained from χ (q) by simpy setting r. Evauating the integras in D = (see aso Ref. ), χ() =, χ() = 3 3 π, χ () = π, (.a) (.b) (.c) and using u = ǫ/5, we finay obtain η = ǫ /5+O(ǫ 3 ), in agreement with the fied theoretica resut. D. The scaing function Let us now cacuate the dimensioness scaing function σ (x) which we define by σ (x) = e ( c) Z Γ () (e ( c) x). (.) () Sincethetwo-pointvertex Γ (q) and the inhomogeneity Γ () (q)dependonyonq = q, weshainthissubsection

11 use scaar arguments for these functions to simpify the notation. For the definition (.) for σ (x) yieds the universa scaing function defined in Eq. (.), im σ (x) σ(x) = (ρ kc) Σ(k c x) Σ(). (.3) We first perform some exact manipuations. The functiona fow equation for Γ () (q) given by Eq. (3.) can be transformed into an integra equation Γ () (q) = e + + dτητ Γ () = (e q) dte t dτητ () t Γ t (e t q) = e dτητ Γ () = (e q) d e + dτητ () Γ (e ( ) q).(.) To describe a critica system, we choose the initia vaue Γ () = (q) = r to be momentum independent such that for the reevant couping r = e dτητ r + d e + dτητ () Γ () (.5) has a finite imit r = im r, just has we have done in Eq. (3.3). This is guaranteed if the initia vaue r is chosen such that Defining r = Γ (mi) (q) = d e + dτητ () Γ (). (.6) () () Γ (q) Γ (), (.7) the superscript mi indicates that these function contain ony margina and irreevant parameters, we obtain σ (x) = e ( c) Z + r d e ( c)+ dτητ (mi) Γ c x). (.8) (e The first term on the right hand side vanishes for because by construction im r = r is finite on the critica surface. Taking the imit we thus obtain (after renaming ) σ(x) = de ( c)+ dτητ (mi) Γ (e c x). (.9). Truncation of the fow equations So far, no approximation has been made. We now approximatethefunction (q) on the right-hand side Γ (mi) of Eq. (.9) by the eading term in the expansion in powers of u, see Eq. (.8), Γ (mi) 3 (q) (+r ) 3 = (+r ) dte ǫt t dτητ u t χ t (e t ˆq +q ) χ t (e t ) ˆq d e ǫ( ) dτη τ u χ (e ( ) ˆq +q ) χ (e ( ) ) ˆq. (.5) Aquantitiesontheright-handsideofEq.(.9)arenow known, so that we may cacuate the crossover function by performing the four-dimensiona integration (two anguar integrations and two integrations over the scae parameters and ) numericay. There are no divergences even in D = 3 provided the integra equation (.9) for η is soved sef-consistenty.. Resuts for the scaing function To obtain an anaytic approximation for the scaing function σ(x), we adopt again the successfu strategy used in the soution of the integra equation (.9) for η: we substitute the Fermi function ansatz (.33) for the fowing anomaous dimension η τ on the right-hand side of Eqs. (.9) and (.5). After some transformations of the integration variabes we then obtain σ(x) = 3u x η with xe c dy y 3+ǫ x ǫ +y ǫ η/ǫf(x,y;η, c), z ǫ (.5) F(x,y;η, c ) = dz x ǫ +(y/z) ǫ η/ǫ +r c+n y zx χ c+n y(z) χ x c+n y( zˆq +yˆq ) ˆq, (.5) x ˆq is an arbitrary unit vector. Note that by assumption x e c so that we may repace xe c in the ower imit of the y-integra in Eq. (.5). Consider first the regime x. Using the fact that for c we may approximate r r and χ (q) im χ (q) χ (q) (see Eq. (.35)), we may repace in this regime χ c+n y(q) χ (q) and r x c+n y r in Eq. (.5). zx Then we obtain y 3+ǫ σ(x) 3u x η dy x ǫ +y ǫ η/ǫf (x,y;η). (.53) Using D-dimensiona spherica coordinates, the function F (x,y;η) im c F(x,y;η, c ) can be written as π Ω D F (x,y;η) = dz dϑ z ǫ (sinϑ) D (+r )Ω D x ǫ +(y/z) ǫ η/ǫ χ (z) χ ( z +zycosϑ+y ). (.5)

12 From Eqs. (.53) and (.5) it is now straightforward to obtain the asymptotic behavior of σ(x) for sma x, B 3 n x + B 3 A D = 3u σ(x) A D x η, (.55) dyy 3+η F (,y;η) 3u Ω π D = dyy 3+η dzz ǫ+η dϑ (+r ) Ω D (sinϑ) D χ (z) χ ( z +zycosϑ+y ). (.56) In D = 3 we find numericay A 3.7. (.57) On the other hand, for arge x the fow parameter in χ c+n y(q) and r x c+n y is typicay sma compared with zx c. To take the effect of the fowing r approximatey into account, we repace in this regime r c+n y r and zx χ c+n y(q) χ (q), χ x (q) is obtained from χ (q) by repacingr r. Using Eq.(3.35) onecan showthat onthe criticasurfaceandforsmau, the initiavaueof r is approximatey given by r u /( ǫ). With this approximation we obtain in the regime (D 3) nx, σ(x) B D x (D 3), (.58) 3u Ω π D B D = dyy 3+ǫ dzz ǫ dϑ (+r ) Ω D (sinϑ) D χ (z) χ ( z +zycosϑ+y ), (.59) and for nx (D 3) (which incudes the imit D 3), B 3 = 3u (+r ) σ(x) B 3nx+B 3, (.6) dz χ (z) 3π u, (.6) and on the right-hand side we have set r. Numericay we find B Note that the physica sef-energy Σ(k) Σ() = ρ k c σ(k/k c) (.6) is independent of u in the regime k k c: keeping in mind thatk c = Λ (u /u ) /ǫ, it isobviousthat thefactor of u in Eqs. (.59) and (.6) combines with a suitabe power of k c so that the sef-energy is proportiona to u. A numerica evauation of σ(x) is shown if Fig. 5. To summarize, using the functiona RG, we have derived in this section the momentum-dependent sefenergy of weaky interacting bosons, covering the entire range from the critica regime up to momenta of n σ(x) A 3 x -η - c - c / c / c n x = n (k/k c ) FIG. 5: Numerica evauation of the scaing function σ(x) of the sef-energy in the imit a and D = 3. Dashed ines are the asymptotic forms predicted by Eqs. (.55,.6). the order of the inverse therma de Brogie waveength. Whie in the critica regime k k c it is simper to obtain the asymptotic ong-waveength behavior of the sefenergy within the fied theoretica RG, the functiona RG approach adopted in this work can describe the entire crossover from the critica regime k k c to the shortwaveength regime k c k π/λ th. Note that the behavior of the sef-energy in the short-waveength regime is determined by irreevant terms which are simpy discarded in the fied theoretica RG. V. INTERACTION INDUCED T c-shift Baym et a. 5 have shown that to owest order in the scattering ength a, the interaction induced shift for the critica temperature T c in D = 3 can be cacuated from the k-dependence of the zero frequency sef-energy Σ(k, ). In a straightforward generaization of their resut to arbitrary dimensions, we can write the contribution from cassica modes to the shift in T c as T c = T c Λ/k c dxx D 3 Ω D ( ) D kc πdζ(d/) Λ (ρ kc) Σ(k c x) Σ() x +(ρ kc) Σ(k c x) Σ(), (5.) Using k c /Λ = e c = (u /u ) /ǫ, and substituting our scaing function σ(x) as defined in Eq. (.), this can aso be written as T c T c = Ω D πdζ(d/) ( u u Λ/k c dxx D 3 )D D σ(x) x +σ(x). (5.)

13 3 The shift in T c is dominated by cassica fuctuations as ong as the UV cutoff Λ /k c can be removed in Eq. (5.). eeping in mind that according to Eq. (.58) σ(x) B D x (D 3) for arge x, we see that the integra <,i.e.d < /3. For D /3 the vaue of the integra in Eq. (5.) depends on the UV cutoff, such that the owest order contribution to the shift in T c is proportiona to u, with ogarithmic corrections ( u nu ) in D = /3. Hence for D /3 the shift in the critica temperature cannot be obtained from an effective cassica fied theory, whie for D < /3, the cutoff ony eads to higher order corrections (in D = 3 one finds u and u n(u ) corrections 5). Writing (5.)isUVconvergentasongas D D we obtain for D < /3, ( Ω D J D = πdζ(d/) T c = J D u D D, (5.3) T c )D D dxx D 3 σ(x) u x +σ(x). (5.) Note that J D impicity depends on u via σ(x). However, for u the coefficient J D approaches a finite imit independent of u as ong as D < /3. In D = 3 one usuay writes eeping in mind that T c T c = c an /3. (5.5) u = 6π ζ(3/) /3 an /3 for D = 3, (5.6) we have in this case c = 6π ζ(3/) /3 J 3. Within our approximation we find c =.3, in rather good agreement with the most accurate numerica investigations which give c =.3±. and c =.9±.5 5. Variationa perturbation techniques give simiar resuts, c =.3 ±. 6, and c =.7 ±. 7. More resuts from other approaches are summarized in a recent review. There is aso experimenta data from a He-VycorsystemconfirmingtheinearscaingofT c with an /3 5. The reported vaue c.66 seems however rather arge which may be attributed to uncertainties regarding the precise vaue of a in the experiment, see the discussion in Ref. 5. VI. IMPROVED DESCRIPTION OF THE CRITICAL REGIME IN D = 3: INCLUDING MARGINAL TERMS Within the u -truncation empoyed so far, we have negected the six-point vertex and approximated the inhomogeneity in the fow equation of the four-point vertex by Γ () (q,q ;q,q ) u χ ( q +q ) + χ ( q q )+ χ ( q q ), (6.) see Eq. (.9). This approximation amounts to a truncation of the exact hierarchyof fow equations on the right-hand side of Eq. (3.5) we simpy repace Γ () (q,q ;q,q ) u. (6.) A forma justification for this procedure can ony be given for D > 3: for weak interactions, irreevant couping parameters can be expanded in powers of reevant and margina ones,. For 3 < D < the ony part of the four-point vertex which is not irreevant is its constant part u. However, in D = 3 there are two additiona margina parameters reated to the momentum dependence of the four-point vertex. Hence, in this case Eq. (6.) is not consistent. For sma momenta the fowing four-point vertex has the form Γ () (q,q ;q,q ) = u +a ( q q + q q )+b q +q + Γ (i) (q,q ;q,q ), (6.3) (i) the irreevant part Γ (q,q ;q,q ) vanishes at east quadraticay if a momenta become sma. Beow, we wi denote the diagrams entering the renormaization of the four-point vertex using the usua fermionic anguage (i. e. BCS, ZS and ZS, see Fig. ). We emphasize that this terminoogy simpy describes the topoogy of the diagrams, not the physica phenomena arising from the diagrams. Note that the four-point vertex must be symmetric under the exchange q q and q q, so that both zero sound channes are characterized by the same parameter a. The parameter b characterizes the BCS channe. The exact fow equations for the couping parameters a and b are a = (3 D η )a + Γ a, (6.) b = (3 D η )b + Γ b, (6.5) the coefficients Γ a and Γ b are defined via the expansion of the function Γ () (q,q ;q,q ) in Eq. (3.5) for sma momenta, Γ () (q,q ;q (),q ) = Γ (,;,) + Γ a ( q q + q q ) + Γ b q +q +O(q i ). (6.6) Even if the initia vertex Γ() = (q,q ;q,q ) is momentum-independent (corresponding to the initia conditions a = b = ), finite vaues of these parameters are generated as we iterate the RG. Because in D = 3 these are margina, they cannot be ignored and are possiby the source of ogarithmic corrections. Besides a and b, there is a third margina parameter in D = 3, the momentum-independent part of the sixpoint vertex, v = Γ (6) (,,;,,), (6.7)

14 see Eq. (3.9). It satisfies the exact fow equation v = (6 D 3η )v + Γ (6) (,,;,,). (6.8) Again, the RG fow generates a finite v even if initiay v =. At the Wison-Fisher fixed point in D = 3 the renormaized v = im v is of order unity. To arrive at the RG equations, we adopt the foowing (n) approximation scheme: We expand Γ to some fixed power in the margina parameters a,b,v, and the reevant parameter u whie keeping higher orders in the reevant parameter r. The truncation is based on the hope that, athough a, b, v and u are not parametricay sma, these parameters sti remain numericay sma. The fixed point vaues for these parameters which we obtain within this scheme in Sec. VID are indeed much smaer than unity, which gives some a posteriori justification to our approach. To impement the scheme, it is usefu to transform Eq. (3.9) into an integra equation, ({q i }) = e ( D)n+D n (n) dτητ Γ = ({e q i }) Γ (n) + dte ( D)n+Dt n dτητ t Γ(n) t ({e t q i }). (6.9) We use these retarded expressions in the cacuations of (m) the inhomogeneities Γ. Since the right hand side in Eq. (6.9) depends on vertices which aso have the form (6.9), one quicky arrives at rather compex expressions which we need to truncate. We wi ony incude terms up to second order in a,b,v, and u. However, even with this approximation the cacuation becomes rather invoved. To keep the cacuations tractabe, we truncate further: To evauate the fow of parameters associated with the n-point vertex, we ignore irreevant contributions arising from vertices of order m with m > n +. Thus, when we cacuate the fow of parameters associated with the n-point vertex and encounter a vertex Γ (m) in the cacuation with m > n+, we ony keep its reevant and margina part instead of empoying Eq. (6.9). Effectivey, we negect in this way infinitey many irreevant parameters reated to higher order vertices. The expectation is, that these irreevant contributions are numericay sma compared to the contribtutions originating from irreevant terms of ower order vertices. In concrete terms, our truncation is as foows: When we cacuate the fow of v, we wi consider irreevant contributions from Γ (8) () but ignore contributions from Γ (in this case, there are in fact no second order terms arising (n) from inhomogeneities Γ with n so that we reay have a second order terms incuded in our fow of v ). Simiary, to cacuate the fow of the four-point vertex parameters a, b, and u, we ignore contributions arising (8) from Γ. This eads to the absence of v terms in the (8) RG equations for a, b, and u, the Γ terms ony enter indirecty via the fow of v. To cacuate the fow of r we (n) ignore irreevant contributions from Γ for n 6 which again ony enter via their contributions to the fow of a, b, and u. A. Truncating the fow equation for the six-point vertex We begin with the renormaization of v, the momentum independent part of Γ (6), using the fow equations for the six-point vertex from Appendix A, Eq.(A). Note that, working consistenty to second order, the ast two terms in Eq. (A), which are at east of third order in the margina and reevant parameters, shoud not be incuded. However, there is a contribution proportiona to v arising from the eight-point vertex which must be incuded. To see this, et us expicity write Γ (6) up to second order in the margina and reevant parameters (and zeroth order in η), Γ (6) (q,q,q 3 ;q 3,q,q ) = 3v +r {3S (,),3 G (q +q ˆq)(u +a ( q ˆq + q ˆq )+b q +q )ˆq +3S (, ),3 S (,),3 G (q 3 q 3 +ˆq)(u +a ( q 3 q 3 )+b q 3 +ˆq )ˆq +(q i q i ) } + +r dt e (D 8)t Γ(8) t (e t q,...,e t q 3,e tˆq;e tˆq,e t q 3,...,e t q ) ˆq, (6.) we used Γ (8) (q,...,q ;q,...,q ) = dt e t dτητ+(8 3D)t Γ (8) t (e t q,...,e t q ;e t q,...,e t q ), (6.) and from Eq. (A3) we have Γ (8) (q,...,q ;q,...,q ) = { v 6S (,,3), S (,,3 ), χ (q +q +q 3 q ) } +8S (,),(3,) S (, ),(3, ) χ (q +q q q ) + terms at east cubic in v, a, b, and u. (6.)

15 5 The definition of the symmetrization operators S (,),3 used in Eq. (6.) and simiar ones used in Eq. (6.) can be found in Eqs.(Ab,Aa,Ab) of Appendix A. eeping ony the zeroth orderin a momentum expansion, we have Γ (6) ({q i }) = β v u 5β v a 9β v b (6.3) dt e (D 8)t v +r t 6 χ t ()+8 χ t (e t ), β = χ () = (+r ). (6.) B. Improved truncation of the fow equation for the four-point vertex Incuding the terms which become margina in D = 3 we need to cacuate Γ () (q,q ;q,q ) + + Γ(6) +r (,BCS) Γ (q,q,ˆq;ˆq,q,q ) ˆq (q,q ;q,q ) (,ZS) Γ (q,q ;q,q ) + Γ (,ZS ) (q,q ;q,q ), (6.5), instead of Eqs. (.3,.), we now have to cacuate the zeroth and first order terms in a momentum expansionofeq.(6.5)inpowersofboththereevantand margina parameters. We discuss the terms in Eq. (6.5) beow, beginning with the contribution from the six-point vertex.. Contributions from Γ (6) To cacuate Γ(6) (q,q,ˆq;ˆq,q,q ), we write, assuming Γ ˆq (6) = =, Γ (6) (q,q,ˆq;ˆq,q,q ) = dt e 3 t dτητ+(6 D)t Γ (6) t (e t q,e t q,e tˆq;e tˆq,e t q,e t q ). (6.6) Using Eqs. (6.6,6.,6.), we can cacuate the contribution of Γ (6) () to Γ. Incuding from the eight-point vertex ony the q i = contribution and expanding the other contributions up to inear order in the externa momenta we have, Γ(6) ( ) (λq,λq,λˆq,λˆq,λq,λq ) v {u (β + χ (λ))+a β +8 χ (λ)(+λ)+8 Φ (λ)+ ϕ (λ) ˆq ) +b (8β + χ (λ)(+λ)+ Φ (λ)+ ϕ (λ) +λ q +q (β u +(β +β )a +β b ) } +λ( q q + q q )(β u +(β +β +β )a +(β +β )b ) +r dt e (D 8)t v t 6 χ t ()+8 χ t (e t ), (6.7) β = χ () = S (+r ) S (+r ) 3, (6.8) S β = (+r ), (6.9) and S and S are defined in Eqs. (.5) and (.6). We further introduced the functions and Φ (λ) = G (ˆq+ˆq λ) ˆq+ˆq λ +r, (6.) ˆq ϕ (λ) = (+r ) ˆq+ˆq λ. (6.) ˆq Since to zeroth order in λ Eq. (6.7) determines the fowofv viaeq. (6.8), the integraon the righthand side of Eq. (6.6) gives simpy v to order λ. We therefore spit the anguar average of the six-point vertex into two contributions, Γ(6) (q,q,ˆq;ˆq,q,q ) ˆq v + Γ(6i) (q,q,ˆq;ˆq,q,q ), (6.) ˆq

16 6 Γ (6i) contains the irreevant parts of the six-point vertex. To cacuate the anguar average of the irreevant part, wecannegecttheη dependence ofthe integraand aso approximate the retarded dependence of the parameters on t by an dependence, since the retardation is exponentiay damped. For arge we therefore can write the contributions of the irreevant part to Eq. (6.6) as foows Γ(6i) (q,q,ˆq;ˆq,q,q ) ˆq v dλ ( ) {6u ( χ (λ) β )+a χ (λ)(+λ)+ Φ (λ) 8β + ϕ (λ) ϕ () λ ) +λ q +q (β u +(β +β )a +β b ) }. (6.3) +b ( χ (λ)(+λ)+ Φ (λ) β + ϕ (λ) ϕ () +λ( q q + q q )(β u +(β +β +β )a +(β +β )b ) We now turn to the cacuation ofthe BCS, ZS and ZS diagrams in Eq. (6.5). eeping ony the reevant and margina terms in the expansion of the four-point vertex, see Eq. (6.3), yieds a contributions to second order in the margina and reevant parameters arising from these diagrams and the irreevant part of the four-point vertex does not contribute at this order. Γ (,ZS) (q,q ;q,q ) u +a q q χ ( q q ) (a +b )u +a q q +r Θ( < q q +ˆq < e ) q q +ˆq +r ˆq q + ˆq+q + ˆq q + ˆq+q ˆq (a +b Θ( < q q +ˆq < e ) ) +r q q +ˆq +r. Contribution of the BCS, ZS and ZS diagrams The contribution from the BCS-channe to Eq. (6.5) is Γ (,BCS) (q,q ;q,q ) u +b q +q χ ( q +q ) a u +b q +q +r Θ( < q +q ˆq < e ) q +q ˆq +r ˆq q + ˆq q + ˆq q + ˆq q ˆq a Θ( < q +q ˆq < e ) +r q +q ˆq +r ) ( ˆq q + ˆq q )( ˆq q + ˆq q, (6.) )ˆq and the contribution from the zero-sound channe is ˆq q ˆq q + ˆq+q ˆq+q ˆq Θ( < q q +ˆq < e ) a b +r q q +ˆq +r ˆq q ˆq+q + ˆq+q ˆq q. (6.5) ˆq The contribution from the other zero-sound channe (ZS in Fig. ) is obtained by repacing q q on the right-hand side of Eq. (6.5). Expanding to inear order in the momenta, we arrive at Γ (,BCS) u in, (q,q ;q,q ) β +a +u a q +q β u +(β +β )a +(β +β )u a +β u b +β a b. (6.6) Simiary, we obtain for the contribution from the zerosound-channe, Γ (,ZS) in, (q,q ;q,q ) β u +a +b q q β u +(β +β +β )a +(β +β )b +(β +β +β )u a +(β +β )u b +(β +β +β )a b. (6.7) The contribution from ZS is again obtained by repacing q q.

17 7 C. Truncated fow equation for the two-point vertex Tocacuate the fowofr, weneed to cacuate see Eq. (3.8). We write it as Γ () (), (+r ) Γ () () = Γ() (,ˆq,ˆq,) ˆq e dτητ u + dte t dτητ t Γ() t (,e tˆq,e tˆq,) ˆq (6.8) with the second order expression of Γ () given by Eq. (6.5). Thus, we have Γ () (,λˆq,λˆq,) χ (λ) (u +λb +a ) +(u +b +a λ) +(u +a +λ) φ (λ) a (u +b λ+a )+a (u +b +a λ)+b (u +a +λ) ψ (λ) a +b (u +b +a )β u +(b +a ) ϕ (λ) + +r Γ(6) (,λˆq,ˆq ;ˆq,λˆq,) ˆq, (6.9) with ψ (λ) = G (ˆq+ˆq λ) ˆq+ˆq λ. (6.3) +r ˆq D. Fow equations of margina and reevant parameters eeping ony terms up to inear order in λ, which define the fow of u, a and b, one finds Γ () () u +a +b +r, (6.3) which is just the contribution of the reevant and margina parts of Γ() (,ˆq,ˆq,). Higher orders in λ ˆq correspond to the contribution of the irreevant part of the four-point vertex. With Γ() (,ˆq,ˆq,) ˆq = u +a +b + Γ(i) (,ˆq,ˆq,) (6.3) ˆq we can write the contributions of the irreevant parts (we here ignore the contributions of the irreevant part of the six-point vertex beyond those which are impicity contained in the renormaization of u ) as Γ(i) (,ˆq,ˆq,) ˆq dλ + Γ (,ZS) Γ (,BCS) in, Γ(,BCS) λ (,λˆq;λˆq,) (,λˆq;λˆq,)+ Γ (,ZS ) (,λˆq;λˆq,) (,ZS) Γ in, (, λˆq; λˆq, ) (,λˆq;λˆq,) Γ (,ZS ) in, (, λˆq; λˆq, ) ˆq, (6.33) which is identica to the second and higher order in λ contributions contained in Eq. (6.9). It is now straightforward to write down the fow equations for the margina and reevant couping parameters. Thefowofv isdeterminedbyeq.(6.8). With Eq.(6.3) one finds v = (6 D 3η )v 3β v (u +5a +3b ) v +r dλ λ D 9 8 χ ()+9 χ (λ),(6.3) we approximated the retarded dependence on v t and χ t in Eq. (6.3) by v and χ (since the retardation is exponentiay damped) and further took for the ower bound of the integra the imit. Eq. (6.3) obviousy has a v = fixed point soution, which is stabe if 3β (u +5a +3b ) > 6 D 3η. Coecting a contributions to Eq. (6.5) from Eqs. (6.3, 6.6, 6.7), and decomposing Γ () (q,q ;q,q ) into a momentum independent contribution and inear contributions according to Eq. (6.6), we find, using Eqs. (3.9, 6., 6.5), the fow equations for the parameters characterizing the