arxiv: v2 [cond-mat.str-el] 28 Jul 2015

Transcription

1 UV/IR mixing in on-ermi Liquids Ipsita Mandal 1 and Sung-Sik Lee 1, 1 Perimeter Institute for Theoretical Physics, 1 Caroline St.., Waterloo O L Y5, Canada Department of Physics & Astronomy, McMaster University, 180 Main St. W., Hamilton O L8S 4M1, Canada Dated: September, 018 arxiv: v [cond-mat.str-el] 8 Jul 015 We study low-energy effective field theories for non-ermi liquids with ermi surfaces of general dimensions and co-dimensions. When the dimension of ermi surface is greater than one, low-energy particle-hole excitations remain strongly coupled with each other across the entire ermi surface. In this case, even the observables that are local in the momentum space such as the Green s functions become dependent on the size of the ermi surface in singular ways, resulting in a UV/IR mixing. By tuning the dimension and co-dimension of the ermi surface independently, we find perturbative non-ermi liquid fixed points controlled by both UV/IR mixing and interactions. I. ITRODUCTIO There have been intensive efforts to understand unconventional metallic states that lie outside the framework of Laudau ermi liquid theory Among the goals are to construct minimal field theories that capture universal low-energy physics, understand the dynamics in controlled ways, and eventually come up with a systematic classification for non-ermi liquids. on-ermi liquids can arise when a gapless boson is coupled with a ermi surface. One of the important criteria that determines the universal properties of non-ermi liquids is the momentum carried by the critical boson. If the critical boson carries zero momentum, fermions lose coherence across the entire ermi surface. The examples include the Ising-nematic critical point 11,1,0 4 and the ermi surface coupled with an emergent gauge field 5 8. When the critical boson carries a finite momentum at the spin density wave SDW or charge density wave CDW critical points 1,14,15,19, the electrons on hot spots or hot lines play a special role because they remain strongly coupled with the critical boson in the low-energy limit. Another important criterion that characterizes different types of non-ermi liquids is the geometry of ermi surface. Although non-ermi liquids do not have a finite jump in the electron occupation number, ermi surface can be still well-defined through weaker non-analyticities such as power-law singularities of the electron spectral function 9. The ermi surface, identified from a nonanalyticity of the spectral function in a non-ermi liquid state, is inherited from that of the underlying ermi liquid before the coupling with a gapless boson is turned on. The kinematic constraints imposed by the parent ermi surface geometry are important in determining the nature of the resulting non-ermi liquid. In this paper, our goal is to understand how the nature of non-ermi liquids depends on the ermi surface geometry for those cases the critical boson carries zero momentum. When the shape of the ermi surface is globally convex in momentum space, there is no special point on the ermi surface 10. Then ermi surface geometry is classified in terms of the dimension and co-dimension of ermi surface. Throughout the paper, we will use m for the dimension of ermi surface and d for the space dimension. The co-dimension of ermi surface is then d m. d controls the strength of quantum fluctuations, and m controls the extensiveness of gapless modes. Although d and m are discrete in reality, we will treat them as continuously tunable parameters to find controlled examples from which physical dimensions can be approached. Regarding d, theories below upper critical dimensions flow to interacting non-ermi liquids at low energies, as systems above upper critical dimensions are expected to be described by ermi liquids. Concerning non-ermi liquids below upper critical dimensions, theories with m = 1 are fundamentally different from those with m > 1. This is due to an emergent locality in momentum space that is present for m = 1, but not for m > The locality has to do with the fact that observables that are local in momentum space, such as Green s functions, can be extracted from local patches in momentum space without having to refer to global properties of ermi surface 41. By exploiting the locality in momentum space, controlled non-ermi liquid fixed points are found in patch descriptions for m = 1 8,10,16,18,19. In contrast to the case with m = 1, non-ermi liquids with m > 1 are less well-understood. The naive scaling based on the patch description breaks down as the size of ermi surface k qualitatively modifies the scaling through the Landau damping. This is due to a UV/IR mixing, low-energy physics is affected by gapless modes on the entire ermi surface in a way that their effects cannot be incorporated within the patch description through renormalization of local properties of the ermi surface. In this sense, k becomes a naked scale for m > 1. Let us elaborate on the origin of the UV/IR mixing. In a renormalizable relativistic quantum field theory, the UV cut-off enters the low-energy effective theory only through the renormalized parameters that can be defined at a momentum scale far below the UV cut-off. or example, in QED, one can extract any observable at a momentum k 1 << Λ solely in terms of the renormalized mass and charge measured at another momentum k << Λ to the leading order in k 1 /Λ, Λ is a large momentum cut-off. This insensitivity of the long-distance physics to

2 the short-distance physics does not necessarily hold in the presence of ermi surface. When m > 1, low-energy observables defined at a momentum near ermi surface such as Green s functions cannot be described solely in terms of the effective couplings defined near that momentum. When a critical boson, that is coupled to ermi surface, has a large number of flavors or velocity much higher than the ermi velocity, one can find a wide range of energy scale over which the effect of k can be ignored 4. However, the UV/IR mixing becomes eventually important in the low-energy limit as long as the the number of flavors and the velocity ratio are finite either through the Landau damping or through a superconducting instability 4. In this paper, we provide a controlled analysis that shows how interactions and the UV/IR mixing interplay to determine low-energy scalings in non-ermi liquids with general m > 1. The paper is organized as follows: In Sec. II, we introduce a theory which describes the Ising-nematic quantum critical point for a system with an m-dimensional ermi surface embedded in d spatial dimensions. or general m, we identify the upper critical dimension d c m at which one-loop quantum corrections exhibit logarithmic divergences. Using ɛ = d c m d as an expansion parameter, one can perturbatively access the non-ermi liquid states that arise in d < d c m. Sec. III is devoted to the description of the dimensional regularization scheme and the beta function that describes the flow to the non- ermi liquid fixed point. In Sec. IV, the RG equations for the renormalized Green s functions are derived. Based on the one-loop results, the dynamical critical exponent, anomalous dimensions and two-point functions are also computed. Sec. V discusses the physical manifestations of UV/IR mixing in physical dimensions. In Sec. VI, we demonstrate that the expansion is controlled in the small ɛ limit with fixed, is the number of fermion flavors. In particular, we show that the one-loop critical exponents are not modified by the two-loop and possibly by all higher-loop corrections for m > 1 due to the UV/IR mixing. We finish with a summary and some outlook in Sec. VII. Details on the computation of the eynman diagrams upto two-loop order can be found in the appendix. II. MODEL We first consider an m-dimensional ermi surface, which is coupled with a critical boson whose momentum is centered at Q = 0 in d = m + 1 space dimensions. One way of characterizing non-ermi liquids is through scaling behaviors of the fermion and boson Green s functions. or this purpose, it is convenient to focus on the point say K at which the fermion Green s function is defined. At low energies, fermions are mainly scattered along the tangential directions by the critical boson. In the presence of the inversion symmetry, fermions a b IG. 1. a A compact ermi surface can be divided into two halves centered at ±K. or each half, a separate fermionic field is introduced. b The compact ermi surface is approximated by two sheets of non-compact ermi surfaces with a momentum regularization that suppresses modes far away from ±K. near K are most strongly coupled with fermions near its anti-podal point K, whose tangent space coincides with that of K. With this in mind, we divide the closed ermi surface into two halves centred at momenta K and K respectively, and introduce separate fermionic fields ψ +,j and ψ,j representing the corresponding halves, as shown in ig.. In this coordinate system, the action is written as S = [ dk ψ s,j k ik 0 + sk 1 + L k + HL k ]ψ s,j k s=±,j + 1 [ ] dk k0 + k1 + L k φ k φk + 1 dk dq e s φq ψ s,j k + q ψ s,jk. 1 s=±,j Here k is the d + 1-dimensional energy-momentum vector with dk dd+1 k. ψ π d+1 +,j k 0, k i ψ,j k 0, k i represents the fermion field with flavor j = 1,,..,, frequency k 0 and momentum Ki + k i Ki + k i with 1 i d. k 1 and L k k, k,..., k d represent the momentum components perpendicular and parallel to the ermi surface at ±K, respectively. The momentum is rescaled such that the absolute value of the ermi velocity, and the quadratic curvature of the ermi surface at ±K, are equal to one. Because the ermi surface is locally parabolic, it is natural to set the scaling dimension of k 1 L k to be 1 1/. HL k = d c i1,..,in n= i 1,..,i n= n k k i1..k in denotes cubic and higher order terms in L k, k is a scale of dimension 1. The range of L k in dk is set by the size of the ermi surface, which is of order k 1/ in this coordinate

3 system. The Yukawa coupling has been also expanded around ±K, and only the momentum independent part is kept. or the Ising-nematic quantum critical point, which we consider in this paper, we have e + = e. a IG.. or m > 1, any two points on the ermi surface have at least m 1 common tangent vectors. The action is reminiscent of patch theories that have been used to describe non-ermi liquids with onedimensional ermi surface. However, there is also an important difference between the theory with general m > 1 and the one with m = 1. or m > 1, any two points on the ermi surface have at least m 1 common tangent vectors, and the whole ermi surface remains strongly coupled at low energies. ig. illustrates this point for a spherical ermi surface embedded in a three-dimensional momentum space. Suppose that two fermions at general momenta k 1, k are scattered to k 1 + q, k q by exchanging a boson with small momentum q. Because q can be tangential to the ermi surface both at k 1 and k for m > 1, the fermions can stay near the ermi surface before and after scattering. Therefore any two fermions on the ermi surface remain strongly coupled in the low energy limit even though the processes with large momentum exchanges are suppressed. or m = 1, such low-energy scattering is not present except for the two anti-podal points. Because of the coupling that is global in the momentum space, the theory with m > 1 is ill-defined in the k limit unlike the case with m = 1. In other words, low-energy IR observables, such as the fermion and boson Green s functions near k = 0, can not be defined until global properties, such as the size and shape, of the ermi surface are specified at large momenta UV for m > 1. In ermi liquids, this UV/IR mixing is encoded in the Landau parameters which are non-local in the momentum space. It is our goal to understand consequences of the UV/IR mixing in non-ermi liquid states with m > 1. The scale k in HL k provides a large momentum cut-off along the directions parallel to the ermi surface. Although irrelevant by power counting, it is crucial to include the higher order terms to keep the information that the ermi surface is compact. In principle, one has to keep an infinite set of independent parameters c i1,..,i n that encodes the precise shape of the ermi surface away b from K. Here we consider a simplified UV regularization which keeps the essential physics of the higher-order terms, but is simple enough to be amenable to an analytic treatment. Specifically, we consider a regularized kinetic term ] {L dk ψ s,j [ik } k 0 + sk d m + L k k ψ s,j k exp. k s,j Here we keep the dispersion parabolic, but the exponential factor effectively makes the size of the ermi surface finite by damping the propagation of fermions with L k > k 1/, as is illustrated in ig. 1b. In order to control the Yukawa coupling and the strength of UV/IR mixing independently, we tune both the dimension 4 44 and the co-dimension of the ermi surface 18,19,45. To keep the analyticity of the theory in momentum space locality in real space with general co-dimensions, we introduce a spinor 18,19 Ψ T j ψ k = +,j k, ψ,j k to write an action that describes an m-dimensional ermi surface embedded in the d- dimensional momentum space: S = j i e µx/ dk Ψ ] { L j k [iγ } k K + iγ d m δ k Ψ j k exp µ k dk L k φ k φk j dkdq φq Ψ j k + q γ d m Ψ j k. Here, K k 0, k 1,..., k d m 1 includes the frequency and the first d m 1 components of the d-dimensional momentum vector, L k k d,..., k d and δ k = k d m + L k. In the d-dimensional momentum space, k 1,.., k d m L k represents the d m m directions perpendicular tangential to the ermi surface. Γ γ 0, γ 1,..., γ d m 1 represents the gamma matrices associated with K. Since we are interested in codimension 1 d m, we consider only gamma matrices with γ 0 = σ y, γ d m = σ x and Ψ Ψ γ 0. The leading terms in the quadratic action in Eq. are invariant under the scale transformations: K = K b, k d m = k d m b, L k = L k b, Ψ j k = b d+4 m 4 Ψ jk, φk = b d+4 m 4 φ k. 4 In the quadratic action of the boson, only L k φ kφk is kept, because K + kd m is irrelevant under the scaling k 0, k 1,.., k d m have dimension 1 and k d,.., k d have dimension 1/. In the presence of the -dimensional rotational symmetry, all components of k d m,..., k d should be equivalent. The reason why k d m is treated differently from L k = k d,.., k d is because bosons that are strongly coupled to the fermions

4 4 III. DIMESIOAL REGULARIZATIO To access perturbative non-ermi liquids, we fix m and tune d towards a critical dimension, at which quantum corrections depend logarithmically on Λ within the range Λ << k. The Yukawa coupling is dimensionless at d cm = 4 + m. 6 IG.. As the high-energy modes away from the ermi surface are integrated out, the ratio k /Λ grows. k is the size of the ermi surface and Λ is the energy cut-off perpendicular to the ermi surface. or m > 1, the Green s functions are singular in the k /Λ limit, which results in the UV/IR mixing. However, it turns out that this is not the actual upper critical dimension at which the quantum corrections diverge logarithmically in Λ. The shift of the upper critical dimension is a sign that k enters the low-energy physics in a way that is singular in the large k limit, resulting in UV/IR mixing. around ±K have momentum L k >> k d m. Therefore we ignore the dependence on k d m in the boson kinetic term in the effective theory that describes the regions around ±K. The scaling dimension of the Yukawa coupling is a x = 4 + m d. 5 Here, e is dimensionless and µ is a mass scale. We also define a dimensionless parameter for the ermi momentum, k = k /µ. The spinor has the energy dispersion d m 1 with two bands, E k = ± i=1 ki + δ k, which gives an m-dimensionsal ermi surface embedded in the d- dimensional momentum space, defined by the d m equations: k i = 0 for i = {1,..., d m 1} and k d m = L k. Besides k and e, the theory implicitly has a UV cutoff for K and k d m, which we denote as Λ. It is natural to choose Λ = µ, and the theory has two important dimensionless parameters : e, k = k /Λ. If k is the typical energy at which we probe the system, the limit of interest is k << Λ << k. This is because Λ sets the largest energy equivalently, momentum perpendicular to the ermi surface fermions can have, as k sets the size of the ermi surface. We will consider the renormalization group flow generated by changing Λ and requiring that low-energy observables are independent of it. This is equivalent to a coarse-graining procedure of integrating out high-energy modes away from ermi surface. Because the zero energy modes are not integrated out, k /Λ keeps increasing in the coarse graining procedure. We treat k as a dimensionful coupling constant that flows to infinity in the low-energy limit. Physically, this describes the fact that the size of the ermi surface measured in the unit of the thickness of the thin shell around the ermi surface diverges in the low-energy limit. This is illustrated in ig.. b c IG. 4. The one-loop diagrams for the boson self-energy a, the fermion self-energy b, and the vertex correction c. Solid lines represent the bare fermion propagator, as wiggly lines in b and c represent the dressed boson propagator which includes the one-loop self-energy in a. In order to identify the actual upper critical dimension, we consider the one-loop quantum corrections. Since the bare boson propagator is independent of k 0,.., k d m, the loop integrations involving it are ill-defined, unless one resums a series of diagrams that provides a non-trivial dispersion along those directions. This amounts to rearranging the perturbative expansion such that the oneloop boson self-energy is included at the zero -th order. The dressed boson propagator, which includes the one-

5 5 loop self-energy ig. 4a, is given by D 1 k = 1 L k + β d e µ x µ k m 1 K d m L k, 7 to the leading order in k/k, for K / L k, δ k / L k << k. Here β d = Γ d d+m 1 π d 1 cos{ πd } Γ d m Γd is a parameter of the theory that depends on the shape of the ermi surface. See Appendix A for the derivation of the one-loop self-energies. Since the boson propagator depends on e, the higher-loop diagrams are no longer suppressed by e, but by a fractional power of e 18. Moreover, the boson self-energy diverges in the k limit for m > 1. This is due to the fact that the Landau damping gets stronger for a system with a larger ermi surface, as the boson can decay into particle-hole excitations that encompass the entire ermi surface for m > 1. This is in contrast with the case for m = 1, a low-energy boson with a given momentum can decay into particle-hole excitations only near the isolated patches whose tangent vectors are parallel to that momentum. or m = 1, k drops out in Eq. 7, which indicates the absence of UV/IR mixing. Eq. 7 is valid when there exists at least one direction that is tangential to the ermi surface, and it should not be extended to the cases with m < 1 for which the conventional quantum field theories work well. rom now on, we will focus on the case with m > 1. The apparent lack of rotational symmetry in the space of k d m,.., k d in Eq. 7 is because the expression is valid only for the boson whose momentum is almost tangential to the ermi surface at points ±K as is discussed below Eq. 4. or the boson propagator with general momentum, L k in Eq. 7 should be replaced by kd m k d in the presence of m + 1-dimensional rotational symmetry. This is because for any given boson momentum k, one can always find a point on the ermi surface k is tangential to the ermi surface. If one chooses a coordinate system k d m = 0, the boson self-energy takes the exactly same form as in Eq. 7. Since we can do this for any k, the boson propagator with general momentum should be independent of the direction in the space of k d m,, k d. In the following, we will use the expression in Eq. 7 because we are only interested in describing the ermi surface near ±K. The bosons that are strongly coupled to that region have momentum k d m << L k, and kd m k d is reduced to L k. By using the dressed boson propagator, we have computed the one-loop fermion self-energy Σ 1 q ig. 4b in Appendix A. This blows up logarithmically in Λ at a new critical dimension d d c d ć d c m = m + m + 1, m IG. 5. The plots for d c and d c as functions of m. which is smaller than d c for general 1 < m <. ig. 5 shows the plots of d c and d c as functions of m. ow we consider the space dimension d = d c m ɛ. In the dimensional regularization scheme, the logarithmic divergence in Λ turns into a pole in 1/ɛ: Σ 1 q = e/ k m 1 m 6 u 1 ɛ + finite terms iγ Q, 9 to the leading order in q/k, u 1 = 1 π m 4π m 1 sin{π/} βd Γ m+4 Γm/Γ m m+5 Γ. m The one-loop vertex correction in ig. 4c vanishes due to a Ward identity 18. It is noted that one can tune the dimension of ermi surface from m = 1 to m = while keeping ɛ small, thus providing a controlled description for any m between 1 and. This is possible because we are tuning m and d independently. or a given m, we tune d such that ɛ = d c m d is small. To remove the UV divergences in the ɛ 0 limit, we add counterterms using the minimal subtraction scheme. Adding the counterterms to the original action, we obtain the renormalized action which gives the finite quantum effective action:

6 S ren = ] {L dk B ΨBj k B [iγ } k,b K B + iγ d m δ kb Ψ Bj k B exp k j,b + 1 dk B L k φ B k B φ B k B + i e B dk B dq B φ B q B Ψ Bj k B + q B γ d m Ψ Bj k B, 10 j 6 K = Z Z 1 K B, k d m = k B,d m, L k = L k,b, Ψk = Z 1/ Ψ Ψ B k B, φk = Z 1/ φ φ B k B, d m/ e B = Z 1/ Z µ x/ e, k = µ k, 11 Z 1 with d m d m Z Z Z Ψ = Z, Z φ = Z. 1 Z 1 Z 1 The subscript B denotes the bare quantitites. To the one-loop order, we have Z n = 1 + Zn,1 ɛ with Z 1,1 = e/ u 1 k m 1 m 6 Z,1 = 0, Z,1 = 0., 1 The one-loop beta functions, that dictate the flow of k and e with the increasing logarithmic length scale l, are given by with d k = dl k, 14 [ ] de ɛ dl = m 1 m + e u 1 ẽ 4 m + 1 e, 15 ẽ e/ k m 1 m k increases under the RG flow because the size of the ermi surface, measured in the unit of the floating energy scale µ exp l, increases at low energies. The first term in the beta function of e indicates that e remains strictly relevant at d = d c m for 1 < m <. However, the form of the loop correction the second term implies that the higher order corrections are controlled not by e, but by an effective coupling ẽ. This can be also checked for higherorder diagrams. The beta function for ẽ, which no longer contains k, is given by dẽ m + 1 ɛ = ẽ m + 1 u 1 ẽ, 17 dl to order ẽ. Eq. 17 shows that the effective coupling flows to an IR stable fixed point at ẽ = ɛ u 1 + Oɛ. 18 or small ɛ, the interacting fixed point is perturbatively accessible despite the fact that the scaling dimension of the bare coupling e stays positive in the ɛ 0 limit for 1 < m <. Although e grows at low energies, the increase of the bare coupling is compensated by the Landau damping, which also increases with the effective size of the ermi surface. The competition between the interaction and the Landau damping makes the effective coupling marginal at the new critical dimension d c. It is interesting that k drops from the effective coupling not only for m = 1 but also for m =. or the latter case, the k dependence in the Landau damping cancels out the k dependence from the phase space of intermediate states in ig. 4b. However, the UV/IR mixing is present for all m > 1 because the Landau damping diverges in the large k limit. IV. REORMALIZATIO GROUP EQUATIOS The renormalized Green s functions, defined by φk 1..φk nφ Ψk nφ +1..Ψk nφ +n ψ Ψk nφ +n ψ +1.. Ψk nφ +n ψ = G n ψ,n ψ,n φ {k i }; ẽ, k, µ δ d+1 nφ +n ψ i=1 k i n ψ +n φ j=n φ +n ψ +1 k j, satisfy the RG equations

7 { n ψ +n φ i=1 z K i Ki + k i,d m + L k i Lki k i,d m d k dl +n φ d c ɛ + 4 m + η φ + d c ɛ + 1 m 4 + d c ɛ mz 1 dẽ k dl ẽ + n ψ } 7 d c ɛ + 4 m + η ψ 4 G n ψ,n ψ,n φ {k i }; ẽ, k, µ = Here the dynamical critical exponent z and the anomalous dimensions are given by z = m + 1ɛ, η ψ = η φ = ɛ, 0 to the one-loop order at the fixed point. It is remarkable that the exponents are insensitive to the details of the ermi surface such as β d despite the fact that patch scaling is violated by k. This vindicates our use of the exponential cut-off scheme in Eq. which captures the compactness of the ermi surface in a minimal way without including the details of the shape. The finite anomalous dimensions are the result of the dynamical balance between the two strongly relevant couplings, e and k. This is opposite to the case finite anomalous dimensions result from a balance between two irrelevant couplings 19. rom this, one can write down the general scaling form of the two-point functions at the IR fixed point as: G 0,0, 1 K 1/z = φ f D L L, k d m, L k, k k k k G 1,1,0 1 K 1/z = δ k f G, δ k, L k, 1 ψ δ k k k φ = 1 z 1 ɛ ηφ = 1 + Oɛ, ψ = 1 z 1 ɛ ηψ = 1 + Oɛ. To the one-loop order, the universal scaling functions are given by [ ] f D X, Y, Z = 1 + β d ẽ X Z m 1 1, [ f G X, Y, Z = i C Γ ˆK ] 1 X + γ d m, { in the Y, Z 0 limit with fixed X. Here C = µ ɛ 1 γ ɛ 6 }, and γ is the Euler-Mascheroni constant. It is noted that f D has a singular dependence on Z in the small Z limit. The absence of the sliding symmetry makes the fermion Green s function depend on δ k and L k separately in general. V. PHYSICAL RELEVACE O THE DIMESIOAL EXPASIO The motivation for the dimensional expansion is to understand the stark difference in the behaviors of non-ermi liquids in two and three dimensions. on- ermi liquids in two dimensions can at most have onedimensional ermi surface, for which k does not play an important role in the low-energy scaling. On the other hands, k enters as an important scale in three dimensions due to UV/IR mixing. Our goal is to understand how this transition occurs in a systematic way by tuning the dimension and co-dimension continuously. Although systems with non-integer dimensions are unphysical by themselves, they provide an insight on how dimension and co-dimension contribute to different behaviors of metals in the physical dimensions. In physical dimension with d =, m =, k drops out in the effective coupling. However, the non-trivial UV/IR mixing still manifests itself in the dispersions of fermion and boson. ear d =, fermion has the dispersion k 0 k x + L k, as the boson has the scaling k 0 L k upto small corrections. Our theory provides a scaling consistent both for boson and fermion by including k as a dimensionful parameter of the theory. Boson and fermion can have different effective dynamical critical exponents at the scale-invariant fixed point, because the difference of the dynamical critical exponents is compensated by k. This is in contrast to the case with m = 1, the dispersions of the boson and fermion obey the same scaling behavior. UV/IR mixing also plays an important role in suppressing higher-loop quantum corrections for m > 1. In the following section, we will examine the effect of k in higher-loop diagrams. VI. EXPASIO PARAMETER The present work is an extension of the early work 18 which provides a controlled expansion for non-ermi liquids supporting one-dimensional ermi surface m = 1. or m = 1, it was explicitly shown that two and threeloop diagrams are suppressed by positive powers of ẽ, which is of the order ɛ. There also exists a general argument outlining why higher-loop diagrams are systematically suppressed by higher powers of ẽ. The expansion here is different from an expansion in powers of 1/, and it does not suffer from the proliferation of planar diagrams in the 1/ expansion 9,1. Due to the addition of the extra co-dimensions, the density of states is suppressed at low energies. The weaker IR singularity allows one to have a controlled expansion for a sufficiently small ɛ. In the present paper, we have generalized this to the cases the dimension of ermi surface is greater than one m > 1. or m > 1, the suppression of higher-

8 or m = 1, k -dependence drops out from every. Since ẽ Oɛ within the perturbative window, one always deals with the limit, λ cross << 1, for m = 1. 7 The situation is different for m > 1. Unlike the case with m = 1, the tangential momentum carried by the boson depends on both Λ and k. or a fixed value of ẽ Oɛ, one is always in the limit of λ cross >> 1, for m > 1, 8 8 IG. 6. A two-dimensional slice of an m-dimensional ermi surface. The typical momentum carried by a boson is proportional to α 1/ Λ d m ẽ 1 m 1 Λ 1/. or m > 1, k Λ this momentum is much larger than Λ 1/ in the low-energy limit. As a result, the momentum transferred from a boson takes a fermion near the ermi surface outside the thin sell of the UV cut-off. This leads to a suppression of the virtual particle-hole excitations by powers of Λ/k for m > 1. loop diagrams by positive powers of ẽ is unchanged. The difference for m > 1 is the presence of an additional scale k. In order to estimate the magnitudes of higher-loop corrections, we first discuss an interplay between k and Λ that plays an important role for m > 1. Suppose k = K, k d m, L k denotes the momentum that flows through a boson propagator within a two-loop or higherloop diagram. When K is order of Λ, the typical momentum carried by a boson along the tangential direction of the ermi surface is given by L k α Λ d m, 4 α = β d e µ x µ k m 1. 5 This can be seen from the form of the boson propagator, which is given in Eq. 7. If α Λ d m 1/ >> Λ 1/, the momentum imparted from the boson to fermion is much larger than Λ 1/ as is illustrated in ig. 6. In this case, the typical energy of virtual particle-hole excitations within the loop is much larger than Λ. As a result, the loop contributions are suppressed by a power of Λ/k at low energies. On the other hand, there is no such suppression if α Λ d m 1/ << Λ 1/. The crossover is controlled by the dimensionless quantity, λ cross ẽ k Λ m 1, 6 which determines whether α Λ d m 1/ α Λ d m 1/ << Λ 1/. >> Λ 1/ or at sufficiently low energies. This is because k has a positive scaling dimension, and k /Λ flows to in the low-energy limit. The crossover occurs at the energy scale Λ ẽ m 1 k. It is noted that there exists a large energy window for small ɛ and m 1, before the theory enters into the low-energy limit controlled by λ cross >> 1. There can be non-trivial quantum corrections from higher-loop diagrams in this intermediate energy scale. We postpone the detailed study of the intermediate scale to a future work. Here we focus on the low-energy limit with λ cross >> 1. In this limit, higherloop diagrams are suppressed by k as was shown in an earlier work 46, for the special case of m = in three dimensions. or general m > 1 with λ cross >> 1, the two-loop self-energies for boson and fermion are given by Π q k Σ a q ẽ m 1 ẽ m m 1 Λ k Q m L q Π 1q, 9 m 1 i γ d m δ q, 0 m 1 Σ b q ẽ m Λ i Γ Q k 1 to the leading order in Λ/k. See Appendix B for the derivation of the results. Here Π q is the two-loop boson self-energy contributed from ig. 7a. Σ a q and Σ b q are the fermion self-energy from ig. 8a, which are proportional to γ d m δ q and Γ Q respectively. Other diagrams in igs. 7b-e and 8b-c do not contribute 18. The coefficients in the expressions for Σ a,b q vanish at d m = 1. The vertex correction is related to the fermion self-energy through the Ward identity. Compared to the one-loop self-energies, the two-loop corrections are suppressed not only by ẽ but also by powers of Λ/k. Because of the suppression by 1/k, there is no logarithmic or higher-order divergence at the critical dimension. As a result, the critical exponents are not modified by the two-loop diagrams in the k limit. It is noted that the suppression by Λ/k originates

9 9 from the large Landau damping which suppresses quantum fluctuations at low energies. Since the suppression is not specific to the two-loop diagrams, we expect that all higher-loop diagrams are also suppressed by ẽ and 1/k in the low-energy limit. We have checked this explicitly for some three-loop diagrams. VII. COCLUSIO To summarize, we have extracted the scaling behaviour of non-ermi liquids with a ermi surface of general dimensions and co-dimensions based on a dimensional regularization scheme. or m > 1, the low-energy physics becomes sensitive to the size of ermi surface k, which results in UV/IR mixing. As a result, the upper critical dimension is shifted from the one predicted by the power-counting, and the perturbative expansion is controlled by a combination of the Yukawa coupling and the ermi momentum k. By tuning the dimension below the upper critical dimension, we have shown that there exists a stable non-ermi liquid fixed point both interaction and UV/IR mixing play crucial roles. We have also shown that the critical exponents at the low-energy fixed point are not modified by the two-loop diagrams, due to the UV/IR mixing for m > 1. This is likely to be the case for all higher-loop diagrams as well. So far we have not considered the four-fermion interaction V, which has the tree-level scaling dimension d+1+m/. However, scatterings in the pairing channel are enhanced by the volume of the ermi surface k m/. As a result, the effective coupling that dictates the potential instability, driven by the four-fermion interactions, is given by Ṽ = V km/, which has an enhanced scaling dimension d m. Ṽ is marginal at the tree-level for co-dimension d m = 1. or a co-dimension d m > 1, there is no perturbative pairing instability for a sufficiently small ɛ = d c d. When d m 1 ɛ with d c d ɛ, the interaction plays an important role for the pairing instability 47,48. ACKOWLEDGMETS We thank Ganapathy Baskaran, Liam itzpatrick, Shamit Kachru, Max Metlitski, Sri Raghu, and Subir Sachdev for stimulating discussions. The research was supported by SERC, ERA and the Templeton oundation. Research at the Perimeter Institute is supported in part by the Government of Canada through Industry Canada, and by the Province of Ontario through the Ministry of Research and Information. 1 T. Holstein, R. E. orton and P. Pincus, Phys. Rev. B 8, M. Y. Reizer, Phys. Rev. B 40, P. A. Lee and. agaosa, Phys. Rev. B 46, B. I. Halperin, P. A. Lee and. Read, Phys. Rev. B 47, J. Polchinski, ucl. Phys. B 4, B. L. Altshuler, L. B. Ioffe and A. J. Millis, Phys. Rev. B 50, Y. B. Kim, A. urusaki, X.-G. Wen and P. A. Lee, Phys. Rev. B 50, C. ayak and. Wilczek, ucl. Phys. B 417, ; ucl. Phys. B 40, S.-S. Lee, Phys. Rev. B 80, S. Sur and S.-S. Lee, Phys. Rev. B 90, M. A. Metlitski and S. Sachdev, Phys. Rev. B 8, M. A. Metlitski and S. Sachdev, Phys. Rev. B 8, M. J. Lawler, D. G. Barci, V. ernandez, E. radkin, and L. Oxman, Phys.Rev.B 7, ; M. J. Lawler and E. radkin, Phys.Rev.B 75, A. Abanov, A. V. Chubukov and J. Schmalian, Advances in Physics 5, A. Abanov, A. V. Chubukov, Phys. Rev. Lett. 84, ; ibid. 9, D.. Mross, J. McGreevy, H. Liu, and T.Senthil, Phys. Rev. B 8, H.-C. Jiang, M. S. Block, R. V. Mishmash, J. R. Garrison, D.. Sheng, O. I. Motrunich and M. P. A. isher, ature 49, D. Dalidovich and S.-S. Lee, Phys. Rev. B 88, S. Sur and S.-S. Lee, Phys. Rev. B 91, V. Oganesyan, S. A. Kivelson, and E. radkin, Phys. Rev. B 64, W. Metzner, D. Rohe, and S. Andergassen, Phys. Rev. Lett. 91, ; L. Dell Anna and W. Metzner, Phys. Rev. B 7, ; Phys. Rev. Lett. 98, H.-Y. Kee, E. H. Kim, and C.-H. Chung, Phys. Rev. B 68, J. Rech, C. Pèpin, and A. V. Chubukov, Phys. Rev. B 74, P. Wölfle and A. Rosch, J. Low Temp. Phys. 147, D.L. Maslov and A.V. Chubukov, Phys. Rev. B 81, J. Quintanilla and A. J. Schofield, Phys. Rev. B 74, H. Yamase and H. Kohno, J. Phys. Soc. Jpn. 69, H. Yamase, V. Oganesyan, and W. Metzner, Phys. Rev. B 7, C. J. Halboth and W. Metzner, Phys. Rev. Lett. 85, P.Jakubczyk, P. Strack, A. A. Katanin, and W. Metzner, Phys. Rev. B 77, M. Zacharias, P. Wölfle, and M. Garst, Phys. Rev. B 80,

10 10 E.-A. Kim, M.J. Lawler, P. Oreto, S. Sachdev, E. radkin, and S. A. Kivelson, Phys. Rev. B 77, Y. Huh and S. Sachdev, Phys. Rev. B 78, O. I. Motrunich, Phys. Rev. B 7, S.-S. Lee and P. A. Lee, Phys. Rev. Lett. 95, P. A. Lee,. agaosa and X.-G. Wen, Rev. Mod. Phys. 78, ; references there-in. 8 O. I. Motrunich and M. P. A. isher, Phys. Rev. B 75, T. Senthil, Phys. Rev. B 78, S.-S. Lee, Phys. Rev. B 78, If there is a superconducting instability, the locality in momentum space breaks down even for m = 1. Here we assume d m > 1, for which there is no perturbative pairing instability. 4 A. L. itzpatrick, S. Kachru, J. Kaplan and S. Raghu, Phys. Rev. B 88, ; ibid. 89, ; R. Mahajan, D.M. Ramirez, S. Kachru, S. Raghu, arxiv: G. Torroba and H. Wang, arxiv S. Chakravarty, R. E. orton, and O.. Syljuåsen, Phys. Rev. Lett., 74, T. Senthil and R. Shankar, Phys. Rev. Lett. 10, T. Schäefer and K. Schwenzer, Phys. Rev.D 70, M. A. Metlitski, D.. Mross, S. Sachdev, T. Senthil, Phys.Rev.B 91, I. Mandal and S.-S. Lee, in progress. Appendix A: Computation of the eynman Diagrams at One-loop 1. One-loop boson self-energy In this section, we compute the one-loop boson self-energy : Π 1 q = ie µ x dk Tr [γ d m G 0 k + qγ d m G 0 k], A1 the bare fermion propagator is given by G 0 k = 1 i k d m, we obtain Γ K+γ d m δ k K +δ k { } exp L k µ k. Performing the integration over K + Q + K [ K K + Q K K + Q ] exp L k +L k+q Π 1 q = e µ x dlk dk [ π d ], K K + Q δq + q d k d + K + Q + K µ k we have chosen a coordinate system such that L q = q d, 0, 0,..., 0, without loss of generality. Because of the rotational symmetry in L q, Π 1 q depends only on its magnitude. rom the expression for the integration over k d, [ ] exp k d k d + L q + L q dkd k I δq + k d L q + K + Q + K with and J dz π exp written as: π = exp L q k z µ k = L q e µ x J m 1 exp Π 1 q = L q 8 k π d m k µ k 8π 1 K + Q + K L q 8 k L q, k y, u = dz π exp z q d m L q k, A z + u + y, A for the integration over the remaining components of L k, the self-energy is { } K K + Q dk K K + Q 1 K + Q + K K + Q + K L q, k q d m L q k. A4

11 It is difficult to obtain the exact expression of Π 1 q for general values of k and q. Here we focus on the limit of physical importance k is much larger than all other scales, including external momentum and the UV cut-off of K. In this limit, we use y, u 1 y for y, u 1 to simplify the expression for the self-energy to 11 Π 1 q = e µ x µ m 1 k L q π I 1 d m, Q, A5 I 1 d m, Q = { } dk K K + Q π d m K K + Q 1. A6 Using the eynman parametrization 1 A n1 B = Γn 1 + n 1 n Γn 1 Γn 0 with n 1 = n = 1/, A = K + Q and B = K we rewrite Eq. A6 as I 1 d m, Q = The integrations over K and t give 1 1 dt π π d m 0 t 1 t I 1 d m, Q = d m 1 π d m for 0 < d m <, and the self-energy is obtained to be Γ d m t n1 1 1 t n 1 dt, A7 n1+n [ t A + 1 t B ] { dk Γ d Q d m K K + Q t K + Q + 1 t K 1 }. Γd m + 1 cos{ πd }, A8 Π 1 q = β d e µ x Q d m µ k m 1 L q, A9 β d is defined after Eq. 7. The expression in Eq. A9 is valid for in a slim region of q satisfying boson self-energy in this region. When Q L q k 1 or Q L q k 1 and q d m Q L q k 1 or q d m L q k 1. On the other hand, it is no longer valid L q 1, when k k is large but finite. ow let us consider the q d m L q 1, we need to use the expression k y, u 1 dz y + u π exp 1 z = π y + u A10 for u 1 or y 1, and the self-energy becomes with T K, Q, q d m = Π 1 q = e µ x J m 1 k π d m 8π t 1 + t, t 1 = t = K < Q K > Q dk T K, Q, q d m, dk T K, Q, q d m, A11 A1 A1 { } K K+Q K K+Q 1 K+Q + K. We can estimate t K+Q + K +qd m 1 and t as t 1 Q d Q Q + qd m, t Q q d m d m g, A14 q d m

12 1 gt = dv vd m t 1+ 4v. Hence, near L q = 0, the boson self-energy becomes Π 1 q e µ x J m 1 k f Q, qd m, A15 [ Q d f Q, q d m Q + qd m ] + C Q Q q d m d m g, A16 q d m with C > 0, a constant. The full one-loop boson self-energy crosses over from Eq. A9 to Eq. A15 around Q L q q, d m k L q 1. However, the range of q in which Eq. A15 is needed becomes vanishingly small in the large k k limit. As a result, we can use Eq. A9 for the dressed boson propagator to the leading order in 1/k as will be shown explicitly in the next section.. One-loop fermion self-energy Here we compute the one-loop fermion self-energy: Σ 1 q = ie µ x dk γ d m G 0 q k γ d m D 1 k. We will first compute Eq. A17 using D 1 k dressed with the boson-self energy in Eq. A9. Integrating over k d m, we obtain Σ 1 q = i e µ x dk Γ K Q π d m I, K Q I = = dω m d L k π m 1 0 π L k m 1 L k Π 1 k π m m 1 Γm/ sin{m + 1π/} { β d e µ x µ k m 1 K d m } m A17. A18 Here we have assumed that 0 < m + 1 <, for which the L k integral converges without the exponential damping factor in the fermion propagator. Therefore, we obtain i e / µ x / π m/ I d m, Q Σ 1 q = 6 m 1 Γm/ sin{m + 1π/} β m/ d µ k A19 m 1 m/6 to the leading order in 1/k, I d m, Q = Using the eynman parametrization A7, we obtain I d m, Q = Γ 1 + β Γβ π π d m 1 dk Γ K Q π d m K d m m/ K Q. 0 dt 1 t β 1 t dk Γβ d m 1 Γ d m β Γ d Γ Q = Γβ Γd m β + 1 π 4π d m Q β d m m 6. inally, the fermion self-energy is obtained to be Γ K Q [ t K Q + 1 tk ] 1 +β β d m 1 A0, A1 i e / µ x / Γ Q Σ 1 q = 6 π m 1/ 4π d m m 1 sin{m + 1π/} β m/ µ k m 1 m/6 Q d m 6 d Γ d m 6 Γ d m β Γ d Γm/ Γβ Γd m β + 1. A

13 1 The above expression has been obtained by using the boson self-energy in Eq. A9, which is valid in the large k limit. ow we explicitly check that our use of Eq. A9 in Eq. A17 is valid to the leading order in 1/k. Suppose we use the exact expression of the one-loop boson self-energy, which is valid for all k, to compute the fermion self-energy: Σ exact 1 q = i e µ x dk exp L q k k L k Πexact 1 k γ d m δ q k Γ Q K Q K + δ q k. A Π exact 1 deviates from Π 1 in Eq. A9 for L k < δ Λ k Λ k, Λ is the UV cut-off for K. In this region, we can safely ignore the exponential damping factor and the contribution of L k in δ q k. In this case, we have Σ exact 1 q Σ 1 q i e µ x i 1 = i e µ x dk dkd m Lk =δ L k =0 π d+1 Lk =δ L k =0 dω m dk dkd m π d+1 [ ] γ d m δ q k Γ Q K 1 dl k Q K + δq k L k Πexact 1 k 1 L k Π 1k [ ] γ d m δ q k Γ Q K Q K + δq k i 1 i, A4 Lk =0 d L k L k m 1 Lk =δ L k + e µ x J m 1 k f K, kd m, i d L k L k m 1 = L k =0 L k Π. A5 1k Integrating over L k, we obtain i 1 = 1 1, m ; m + δ ; e µ x J m 1 k f K, kd m i = 1 1, m + 1 ; m + 4 δ ; β d e µ x k m 1 K d m δ m m e µ x J m 1 k f K, kd m, A6 δ, A7 m + 1 β d e µ x k m 1 K d m 1 is the hypergeometric function. In the limit k with fixed Λ, we have i1 I, i I k m m+ 0, I is defined in Eq. A18. Therefore, Σ exact q goes to Σq in the large k limit, and Eq. A is valid to the leading order in 1/k. Appendix B: Computation of the eynman Diagrams at Two-loop IG. 7. The diagrams for two-loop boson self-energy.

14 14 IG. 8. The diagrams for two-loop fermion self-energy. IG. 9. The diagrams for two-loop vertex corrections. All the two-loop diagrams are shown in igs. 7,8 and 9. The black circles in igs. 7 d-e, 8c and 9i-j denote the one-loop counterterm for the fermion self-energy, i A 1 1 ΨΓ QΨ. Among the self-energy diagrams, only igs. 7a and 8a contribute 18. The vertex correction can be obtained from the fermion self-energy correction through the Ward identity. As we will see through explicit computations in the following sections, the two-loop diagrams are suppressed not only by ẽ but also by Λ/k in the low-energy limit for m > 1. B1 1. Two-loop contribution to boson self-energy We compute the two-loop boson self-energy shown in ig. 7 a: Π q = e4 µ x dl dp D 1 l Tr{γ d m G 0 p γ d m G 0 p + l γ d m G 0 p + l + q γ d m G 0 p + q}. Taking the trace, we obtain Π q = e4 µ x dl dp D 1 l B 1 D 1 exp L p + L p+q + L p+l + L p+l+q k B, B B 1 = [ δ p+l δ p+q+l P + L P + L + Q ] [ δ p+q δ p P + Q P ] [ P + L P + Q ] [ P + L + Q P] + [ P + L P] [ P + L + Q P + Q ] [ δ p+l P + L + Q + δ p+l+q P + L ] [ δ p+q P + δ p P + Q ], D 1 = [ δ p + P ] [ δ p+q + P + Q ] [ δ p+l + P + L ] [ δ p+l+q + P + L + Q ]. B4

15 15 Shifting the variables as p d m p d m L p, l d m l d m p d m L p+l, we can substitute δ p p d m, δ p+q p d m + L p L q + δ q, δ l+p l d m, δ p+l+q l d m + L p+l L q + δ q. Integration over p d m and l d m gives us Π q = e4 µ x dll dl dl p dp π d D 1 l B D exp L p + L p+q + L p+l + L p+l+q k B = P + L + P + L + Q P + Q + P { [ P + L P + L + Q P + L P + L + Q ] [ P + Q P P + Q P ] } [ P + L P + Q ] [ P + L + Q P ] + [ P + L P ] [ P + L + Q P + Q ] L p+l L q + δ q Lp L q + δ q, B5 [ P + L + Q P + L P + L P + L + Q] [ P + Q P P P + Q ], B6 [ ] D = 4 P P + Q P + L P + Q + L Lp+l L q + δ q + P + L + P + Q + L [ ] Lp L q + δ q + P + P + Q. B7 Without loss of generality, we can choose the coordinate system such that L q = q d, 0, 0,..., 0 with q d > 0. After making a further change of variables as L L P, P P Q, L q p d + δ q p d, B8 and integrating over p d neglecting the corresponding exponential damping part, we obtain: Π q e4 µ x dll dl du p dp π d π d 1 D 1L l, L P B L, P, Q D l, P, q exp u p k e4 µ x m 1 k dll dl dp 1 π π d π d m D 1L l, L P B L, P, Q D l, P, q, u k = k d m+,..., k d, B L, P, Q = B 4 L, P, Q DL, P, Q, D l, P, q = 8 L q D 4 L, P, Q { D L, P, Q + 4L q L l }, B 4 L, P, Q = L Q/ L + Q/ L + Q /4 P Q/ P + Q/ P + Q /4 [ L Q/ P + Q/ ] [ L + Q/ P Q/ ] + [ L Q/ P Q/ ] [ L + Q/ P + Q/ ] L + Q/ P + Q/ [ L Q/ P Q/ ] + L + Q/ P Q/ [ L Q/ P + Q/ ] + L Q/ P + Q/ [L + Q/ P Q/ ] L Q/ P Q/ [ L + Q/ P + Q/ ], D 4 L, P, Q = L Q/ L + Q/ P Q/ P + Q/, DL, P, Q = L Q/ + L + Q/ + P Q/ + P + Q/. ote that we can ignore the exponential damping part for L l. B9 B10 B11 B1 B1

16 or m > 1 and λ cross >> 1, the angular integrals along the ermi surface directions give a factor proportional to 16 π 0 dθ DL, P, Q sin m θ D L, P, Q + 4 L l L q cos θ π L l L q B14 in the limit DL,P,Q L l L q << 1, which is valid when L q >> Λ 1 λ cross. This follows from the fact that the main contribution to the integral over L l comes from L l α 1 L P d m >> Λ in the large λ cross limit. Using L l m 1 d L l L l + α L P d m = π sin d m m m π L P we perform the L l -integral to obtain α m for 0 < m <, B15 Π q e4 µ x π 48 m 1 k 1 π dl dp B 4 L, P, Q π d m D 4 L, P, Q L q α m d m m L P π sin m π. B16 The total power of e comes out to be m+, and we find Π q ẽ m e k m 1 π m 1 48 L k q sin dl dp B 4 L, P, Q m π π d m ẽ m Q m d m m m 1 D 4 L, P, Q L P L k q Π 1q, B17 to leading order in Q and ɛ. Therefore, the two-loop diagram is suppressed not only by a higher power of ẽ but also m 1 by Q k compared to the one-loop diagram.. Two-loop contribution to fermion self-energy The two-loop fermion self-energy in ig. 8a is given by Σ q = ie4 µ x dp dl D 1 p D 1 l γ d m G 0 p + q γ d m G 0 p + l + q γ d m G 0 l + q γ d m. B18 Using the gamma matrix algebra, we find that the self-energy can be divided into two parts: Σ q = Σ a q + Σ b q, B19 Σ a,b q = i e4 µ x dp dl D 1 p D 1 l [P + Q + δp+q] [P + L + Q + δp+l+q ] [L + Q + δl+q ], C a,b B0 with C a = γ d m [δ p+q δ p+l+q δ l+q δ l+q { Γ P + Q} { Γ P + L + Q} δ p+q { Γ P + L + Q} { Γ L + Q} ] δ p+l+q { Γ P + Q} { Γ L + Q}, B1 C b = [ Γ P + Q ] [ Γ P + L + Q ] [ Γ L + Q ] δ p+q δ l+q [Γ P + L + Q] δ p+l+q δ l+q [ Γ P + Q ] δ p+q δ p+l+q [ Γ L + Q ]. B Shifting the variables as p d m p d m δ q L p L q L p, l d m l d m δ q L l L q L l,

17 and integrating over p d m and l d m, we obtain Σ a q = i e4 µ x 4 Σ b q = i e4 µ x 4 dpdl dl p dl l π d m π m dpdl dl p dl l π d m π m 17 γ d m δ q L l L p C a L, P, Q D 1 p D 1 l δ q L l L p + CL,, P, Q B CL, P, Q Cb L, P, Q D 1 p D 1 l δ q L l L p + CL, P, Q, B4 CL, P, Q = P + Q + L + Q + P + L + Q, [Γ P + Q] [Γ P + L + Q] [Γ P + L + Q] [Γ L + Q] [Γ P + Q] [Γ L + Q] C a L, P, Q = 1 +, P + Q P + L + Q P + L + Q L + Q P + Q L + Q [Γ P + Q] [Γ P + L + Q] [Γ L + Q] [Γ L + Q] [Γ L + P + Q] [Γ P + Q] C b L, P, Q = +. P + Q P + L + Q L + Q L + Q L + P + Q P + Q B5 We can extract the UV divergent parts by expanding the integrand for small δ q in Eq. B, and by setting δ q = 0 in Eq. B4. In this limit, we have: Σ a q = i e4 µ x 4 Σ b q = i e4 µ x 4 dp dl dlp dl l γ d m δ q Ca L, P, Q D 1 p D 1 l { CL, P, Q 4 L l L p } π d { CL,, B6 P, Q + 4 L l L p } dpdl dl p dl l CL, P, Q Cb L, P, Q D 1 p D 1 l π d m π m 4 L l L p + CL,. B7 P, Q or m > 1, the angular integrals along the ermi surface directions give a factor proportional to π C L, P, Q 4 L l L p cos θ π Γ m 1, for dθ { 0 CL, P, Q + 4 L l L p cos θ } C sinm θ Γ m π Γ m 1 L Γ π m C l L m p L l L p, for CL,P,Q L l L p >> 1, CL,P,Q L l L p << 1, for Σ a ; and π Γ m 1 π CL, P, Q sin m CL,P,Q θ, for dθ 0 C L, P, Q + 4 L l L p cos θ C Γ m L l L p >> 1, π L l L p π C Γ m 1 L Γ l L m 1, for CL,P,Q L p l L p << 1, B9 B8 CL,P,Q for Σ b. or L l L p << 1, the leading and the second leading order terms in CL,P,Q L l L p are kept. Although the second leading term is not important for Σ b, it plays an important role for Σ a at m = because the leading term vanishes at m =. Let us first estimate Σ a. It is convenient to perform the integrations for L l and L p in the four regions separately, { 0 < L l < C L p, 0 < L p < C α 1 L d m { } C 0 < L l < L p, C < L α 1 L d m p <, { } C L p < L C l <, 0 < L p <, α 1 L d m { C L p < L l <, C α 1 L d m } < L p <, }. B0

18 Adding the result of integrations in each region, we obtain Σ a q i e4 µ x γ d m δ q π 4 dpdl π d C a L, P, 0 [ π Γ m 1 sec π 18 Γ m 1 α 4 m d m 4 m L P + C m 1 L, P, 0 { m m + 4 m + 1 log A Γ m 1 A m 1 8 α L P d m π m + 1 Γ m m + log A m m/ m + m m 1 A m 5 m m + m m 1 A m + m m 1 + A m }] 5 m m 5 m, B1 AL, P, α = CL, P, 0 L P d m α 18. B Here Q dependent terms are dropped because they are sub-leading compared to the one that depends on δ q. ote m and A m 1 m log A that A m 1 m m/ m become terms that include log A and log A at m =. rom similar integrations, Σ b q is evaluated to be Σ b q i e4 µ x π dpdl Cb L, P, Q 7 π d α m d m m L P sin. B m π or d = d c ɛ, we find that Σ a q i γ d m δ q π 4 dpdl π d C a L, P, 0 [ ẽ m 1 m 1 k π Γ m 1 sec π 18 Γ m 1 L P d m 4 m C m 1 L, P, 0 { m m + 4 m + 1 log A Γ m k m 1 L P d m π m + 1 Γ m + + m m and 1 5 m Σ b q ẽ m k m 1 A α=1 ẽ m 1 k β d m i π 7 m + m 1 5 m dpdl Cb L, P, Q π d d m m L P A α=1 ẽ m 1 k A α=1 m 1 ẽ k m β d m 1 m m/ m β d m m sin m π m 1 + log A A α=1 ẽ m 1 k m β } ] d + m 5 m, / m B4 m 1 ẽ m Λ i Γ Q k to leading order in δ q, Q and ɛ. The two-loop fermion self-energy is suppressed by ẽ and k at low energies for m > 1. Due to the Ward identity, the two-loop vertex corrections shown in ig. 9 are also suppressed. It is noted that Eqs. B4 and B5 are finite in the ɛ 0 limit. This is because the ermi energy k enters as a dimensionful parameter that further suppresses the two-loop contributions compared to the one-loop diagrams. Therefore the critical exponents do not receive quantum corrections from the two-loop diagrams in the low-energy limit. Λ B5