arxiv: v1 [cond-mat.str-el] 15 Jan 2018

Transcription

1 Fluctuation effects at the onset of 2k F density wave order with one pair of hot spots in two-dimensional metals Jáchym Sýkora, Tobias Holder, 2 and Walter Metzner Max Planck Institute for Solid State Research, D-7569 Stuttgart, Germany 2 Department of Condensed Matter Physics, arxiv:8.4885v [cond-mat.str-el] 5 Jan 28 Weizmann Institute of Science, Rehovot, Israel 76 (Dated: January 6, 28) We analyze quantum fluctuation effects at the onset of charge or spin density wave order in two-dimensional metals with an incommensurate 2k F wave vector connecting a single pair of hot spots on the Fermi surface. We compute the momentum and frequency dependence of the fermion self-energy near the hot spots to leading order in a fluctuation expansion (one loop). Non-Fermi liquid behavior with anomalous frequency scaling and a vanishing quasi-particle weight is obtained. The momentum dependence yields a divergent renormalization of the Fermi velocity and a flattening of the Fermi surface near the hot spots. Going beyond the leading order calculation we find that the one-loop result is not self-consistent. We show that any momentumindependent self-energy with a non-fermi liquid frequency exponent wipes out the peak of the polarization function at the 2k F wave vector, and thus destroys the mechanism favoring 2k F density waves over those with generic wave vectors. However, a 2k F density wave QCP might survive in presence of a sufficiently flat renormalized Fermi surface. PACS numbers: 7..Hf,75.3.Fv,64.7.Tg I. INTRODUCTION Order parameter fluctuations at or near quantum critical points in metals interact strongly with fermionic excitations. They can destroy Landau quasi-particles, thus undermining the basis for Fermi liquid theory. Vice versa, the order parameter fluctuations are themselves strongly affected by the gapless fermionic degrees of freedom in metallic systems. An effective order parameter theory obtained by integrating out the fermions is thus

2 2 complicated by singular interactions. The assumption of a regular order parameter theory as originally proposed by Hertz 2 and Millis 3 is therefore not justified, at least not a priori. Instead, one is compelled to treat fermions and their order parameter fluctuations on equal footing. Quantum critical points in metals fall into a variety of distinct universality classes. The geometry of the ordering wave vector Q plays a crucial role. Magnetic order with Q = includes Heisenberg and Ising ferromagnets, while examples for charge order with Q = are given by the director nematic in continuum systems 4 and the Ising nematic on a square lattice. 5,6 Concerning charge and magnetic order with Q, that is, charge and spin density waves, commensurate and incommensurate wave vectors need to be distinguished, and whether Q is a nesting vector of the Fermi surface or not. The most thoroughly analyzed quantum critical point with Q is the one at the onset of commensurate antiferromagnetic order with a non-nested wave vector Q = (π, π) in two dimensions. 7,8 Quantum critical behavior at the onset of charge order with a non-nested incommensurate wave vector has been studied in the context of cuprate superconductors. 9 A special situation arises when Q is a nesting vector connecting Fermi points with antiparallel Fermi velocities. Spin and charge correlations at such wave vectors feature a well-known singularity generated by an enhanced phase space for low-energy particle-hole excitations. In isotropic systems, the singularity is located at wave vectors Q with Q = 2k F, where k F is the radius of the Fermi surface. In inversion symmetric crystalline metals with a valence band dispersion ɛ k, wave vectors Q obeying the condition ɛ (Q+G)/2 = ɛ F are nesting vectors. Here ɛ F is the Fermi energy, and G is a reciprocal lattice vector. We refer to such nesting vectors also as 2k F vectors. 2k F singularities are particularly pronounced in metals with reduced dimensionality. Charge and spin susceptibilities in low-dimensional systems frequently exhibit a peak at 2k F vectors, such that these wave vectors are favored for charge- and spin-density wave instabilities. For example, the ground state of the two-dimensional Hubbard model undergoes spin-density wave instabilities at 2k F -vectors upon increasing the coupling strength, at least within mean-field (Hartree-Fock) theory., Also d-wave bond charge order generated by antiferromagnetic fluctuations in models for cuprate superconductors 8,2 occurs naturally at 2k F vectors. 3 Quantum critical behavior at the onset of density wave order with a 2k F -vector in two-

3 3 FIG. : Fermi surface with a single pair of hot spots ±k H connected by ordering wave vector Q. dimensional metals was first analyzed by Altshuler et al. 4 They derived non-fermi liquid power-laws for the fermion self-energy and the susceptibility for the case where Q is half a reciprocal lattice vector. In the special case where Q = (π, π) is a 2k F vector, additional umklapp processes have to be taken into account. 5,6 For incommensurate 2k F vectors, Altshuler et al. found very strong infrared divergences in two-loop contributions to the susceptibility. They therefore concluded that the order-parameter fluctuations destroy the quantum critical point and replace it by a first order transition. Recently, two of us have analyzed the influence of quantum critical fluctuations at the onset of charge- or spin-density wave order with an incommensurate 2k F wave vector in twodimensional metals on single-particle excitations. 7 To this end, the fermion self-energy was computed to first order (one loop) in the fluctuation propagator, which was also computed in a one-loop approximation. A breakdown of Fermi liquid behavior was obtained at hot spots of the Fermi surface, that is, Fermi surface points that are connected by the ordering wave vector Q. Two qualitatively distinct cases had to be distinguished. In the first case, the ordering wave vector Q connects only a single pair of hot spots and points typically in axial or diagonal direction (see Fig. ). The frequency dependence of the one-loop selfenergy at the hot spots obeys a power-law with exponent 2 in this 3 case.7 In the second case, Q connects two pairs of hot spots, and the imaginary part of the real frequency one-loop self-energy features a linear frequency dependence at the hot spots. 7 In this article we extend the analysis of the case of a single hot-spot pair in several

4 4 directions. First, we compute the momentum dependence of the one-loop self-energy in the vicinity of the hot spots. In particular, we thereby determine the non-fermi liquid regime in energy-momentum space near the hot spots. Second, we analyze to what extent the one-loop susceptibility and fluctuation propagator is modified by the singular fermion self-energy. We find that this feedback effect is substantial, giving rise to various possible scenarios. The article is structured as follows. In Sec. II we discuss the asymptotic behavior of the RPA susceptibility and effective interaction at the QCP. In Sec. III the momentum and frequency dependence of the fermion self-energy near the hot spots as obtained to leading order in the effective interaction (one-loop) is presented. In Sec. IV the issue of self-consistency and stability of the QCP is investigated. Final conclusions follow in Sec. V. The appendices contain a derivation of effective interactions for two familiar microscopic models (Appendix A), and the evaluation of the loop integrals (Appendix B and C). II. RPA SUSCEPTIBILITY AND EFFECTIVE INTERACTION ɛ k. We consider a one-band system of interacting fermions with a bare dispersion relation Our calculations are based on the standard quantum many-body formalism with an imaginary frequency representation of dynamical quantities. 8 The bare fermion propagator is thus given by G (k, ik ) = ik ξ k, () where k is the frequency variable, and ξ k = ɛ k µ is the single-particle energy relative to the chemical potential µ. We assume that, in mean-field theory, the system undergoes a charge or spin density-wave instability with an incommensurate and nested (2k F ) modulation vector Q at a QCP which can be reached by tuning a suitable parameter such as density or interaction strength. Approaching the QCP from the normal metallic regime, the instability is signalled by a diverging RPA susceptibility χ(q, iq ) = χ (q, iq ) + gχ (q, iq ), (2) where g < is the coupling constant parametrizing the bare interaction in the instability

5 5 channel, and χ is the bare susceptibility χ (q, iq ) = N Π (q, iq ) d 2 k dk = N (2π) 2 2π f k q/2 2 G (k, ik ) G (k q, ik iq ). (3) N is the spin-multiplicity (N = 2 for spin- 2 fermions), and f p is a form factor related to the internal structure of the density-wave order parameter. For an order parameter with s- wave symmetry, f p is symmetric under rotations and reflections. In the following we assume f p = for definiteness. A generalization to form factors with other symmetries such as d-wave is straightforward. Note that Eq. (2) holds also for the spin susceptibility in the normal (symmetric) phase of spin-rotation invariant systems, where all components of the spin susceptibility are equal. For the Hubbard model, the RPA spin susceptibility is given by Eq. (2) with g = U/2. Throughout this article we consider the case where the ordering wave vector Q connects a single pair of hot spots k H and k H on the Fermi surface. To parametrize momenta near the ordering momentum Q we introduce coordinates normal and tangential to the line formed by the 2k F -vectors in q-space, which we denote by q r and q t, respectively. For momenta near the ordering momentum Q and low frequencies, the bare susceptibility can then be expanded as 4,7 χ (q, iq ) = χ (Q, ) a h(e q, q ) + b e q c q2 t 4m, (4) where a, b, and c are positive constants, and h(q, q ) = e q + iq + e q iq = 2 e 2 q + q 2 + e q. (5) The energy-momentum relation e q is given by e q = v F q r + q2 t 4m, (6) where v F is the Fermi velocity at ±k H and m parametrizes the Fermi surface curvature at these points (mv F is the radius of curvature); e q /v F is the oriented distance of q from the 2k F -line. In previous works 4,7 the last term in the expansion (4) was discarded. It plays however a significant role as it defines the maximum of χ (q, ) at Q with respect to variations along the 2k F -line. The prefactor of the square-root term is determined by the Fermi velocity and curvature near ±k H as a = N m/(4πv F ), while the other constants b

6 6 and c receive contributions from everywhere. For fermions with a parabolic dispersion in the continuum and a constant form factor, b and c vanish. 9 At the QCP one has gχ (Q, ) =, so that the RPA susceptibility assumes the singular form g χ (Q, ) χ(q, iq ) =. (7) a h(e q, q ) b e q + c q2 t 4m To deal with the critical order parameter fluctuations, the perturbation expansion has to be organized in powers of a dynamical effective interaction. This arises naturally as a boson propagator by decoupling the bare interaction in the instability channel via a Hubbard- Stratonovich transformation. 2,3 Alternatively it can be obtained by an RPA resummation of particle-hole bubbles or ladders. In the simplest case of a charge-density wave instability in a spinless fermion system, the RPA effective interaction can be written as D(q, iq ) = g + gχ (q, iq ). (8) As an example we derive the effective interaction for the specific case of spinless lattice fermions with a nearest neighbor interaction in Appendix A. For a charge-density wave instability in a spin- fermion system, the effective interaction has the diagonal spin structure 2 D 2 2 (q, iq ) = δ δ 2 2 D(q, iq ), (9) where, 2 (, 2) are the spin-indices of the ingoing (outgoing) fermions. For a spindensity wave, the effective interaction acquires a non-diagonal spin structure. In a spinrotation invariant system of spin- 2 fermions, it can be written as D 2 2 (q, iq ) = τ τ 2 2 D(q, iq ), () where τ = (τ x, τ y, τ z ) is the vector formed by the three Pauli matrices τ x, τ y, τ z. We derive the effective interaction for the spin density wave instability in the Hubbard model as an example in Appendix A. Expanding the bare susceptibility as in Eq. (4), the effective interaction at the QCP assumes the asymptotic form D(q, iq ) =. () a h(e q, q ) b e q + c q2 t 4m It thus features the same singularity as the RPA susceptibility. In case that the coupling g has a (regular) momentum dependence for q near Q, the coefficients b and c are not determined by χ (q, ) only, but receive additional contributions from the expansion of g (q) around Q.

7 7 III. FERMION SELF-ENERGY To leading order in the effective interaction, the fermion self-energy is given by the oneloop expression Σ(k, iω) = M d 2 q dq (2π) 2 2π D(q, iq ) G (k q, iω iq ), (2) with M = for a charge density and M = 3 for a spin density instability. We evaluate Σ(k, iω) for low frequencies ω and momenta k near a hot spot k H. The dominant contributions come from momentum transfers q near Q, such that k q is situated near the antipodal hot spot k H. We assume that the Fermi surface is convex at ±k H. Using normal and tangential coordinates for k near k H, we expand the dispersion relation to leading order as ξ k q = v F (k r q r )+ 2m (k t q t ) 2. It is convenient to perform the momentum integration by using q t and e q as integration variables. A. Imaginary part In Appendix B we show that the imaginary part of the one-loop self-energy has the following asymptotic form for low frequencies and momenta near k H : ImΣ(k, iω) = M πn sgn(ω) ω 2/3 I( ξ k, b k t ), (3) with b = b/a and the dimensionless scaling variables ξ k = ξ k ( ω / b) 2/3, kt = k t m( ω / b) /3. (4) The dimensionless scaling function I is given by I( ξ k, k min(, ξk ) dẽ q ẽ t ) = q 3/2 ξ ln + ( k /ẽ q s=± + γ s k t / ẽ q + ξ ) 2, (5) k /ẽ q where γ = c/b. The integral in Eq. (5) can be evaluated analytically for large and small ξ k and k t (see Appendix B). Analytic results for the self energy can thus be derived in various limiting cases. In particular, at the hot spot k H we obtain ImΣ(k H, iω) = 4Msgn(ω) 3N ω b + c 2/3, (6)

8 8 with c = c/a. This result is valid also for k k H as long as b ξ k 3/2 kt 3 ω and b ω. m 3/2 Eq. (6) is consistent with the result ImΣ(k H, ω + i) = 2M 3N ω/ b 2/3 for the real frequency self-energy derived for the special case c = in Ref. 7. For k t = and ξ k < with ω b ξ k 3/2, we obtain 2M ImΣ(k, iω) = N b b ω + c ξk. (7) This result remains valid for k t, as long as ω min( b, c) ξ 3/2, we find k b kt 3 m 3/2 ω. For k t = and ξ k > with ImΣ(k, iω) = 4MArcsinh b/ c πn b b + c ω ξk, (8) kt 3 which is also applicable for k t with b ω. Note the asymmetry of the prefactors in m 3/2 Eqs. (7) and (8) for ξ k < and ξ k >, respectively. The frequency range where Eq. (8) is valid is restricted not only by b ξ 3/2 k but also by c ξ 3/2 k, and thus shrinks to zero for c =. The above results were derived for a convex Fermi surface. For a concave Fermi surface, Eq. (7) holds for ξ k > and Eq. (8) for ξ k <. Finally, for ξ k = and k t with ω min( b, c) kt 3 m 3/2, we obtain ImΣ(k, iω) = 2M m N b c ω k t. (9) This result is also applicable for ξ k, as long as b ξ k 3/2 ω. Eqs. (6) and (9) are also valid if the Fermi surface is concave instead of convex at ±k H. At the hot spots the imaginary part of the self-energy vanishes sublinearly for vanishing ω. Hence, there are no stable quasi-particle excitations at these points, which implies a breakdown of Fermi liquid theory. At momenta away from the hot spots, ImΣ(k, iω) is linear in ω as in a Fermi liquid. We can therefore define a momentum dependent Z-factor ( ) Σ(k, iω) Z k = i ω. (2) ω= In a Fermi liquid, Z k is finite on the Fermi surface (k = k F ) and describes the spectral weight of quasi-particle excitations. For the one-loop self-energy at the 2k F QCP derived above, Z k is well-defined and finite for k ±k H, but it vanishes as k t upon approaching the hot spot along the Fermi surface, and as ξ k upon approaching k H perpendicularly to the Fermi surface.

9 9 B. Real part There is no non-fermi liquid contribution to the frequency dependence of the real part of the self-energy at the hot spots. The symmetry of ImΣ(k H, ω + i) = 2M 3N ω/ b 2/3 under a sign change of ω implies that there is no contribution of order ω 2/3 to ReΣ(k H, iω). However, there are singular (real) contributions to the momentum dependence of Σ(k, ), which renormalize the dispersion and the Fermi surface shape near the hot spots. parametrize the momentum dependence by the variables ξ k, which describes the distance from the Fermi surface, and k t. For momenta close to a hot spot k H in perpendicular direction to the Fermi surface (that is, k t = ), one obtains (see Appendix B) δσ(k, ) = Σ(k, ) Σ(k H, ) = M [ ] N + Θ(ξ k ) + K(γ) ξ k ln Λ ξ k + O(ξ k), (2) where Λ is an arbitrary fixed energy scale, γ = c/b, and K(γ) = 2 ( γ π + γ arctan ) γ. (22) Hence, Σ(k, ) leads to a logarithmically diverging renormalization of the bare dispersion ξ k. The prefactor in Eq. (2) is asymmetric, that is, it depends on the sign of ξ k. The function K(γ) ranges from for γ = to for γ, and it crosses the value /2 at γ = γ c 5.. For ξ k <, the self-energy makes the effective dispersion steeper, with a logarithmically diverging slope at the hot spots. For ξ k >, it becomes steeper only as long as γ < γ c. By contrast, for γ > γ c the logarithmic self-energy correction reduces the effective dispersion outside the bare Fermi surface (where ξ k > ), and even leads to a sign change upon approaching the hot spots, such that the Fermi surface at the hot spots is destroyed. However, the logarithmic divergence indicates a breakdown of perturbation theory for momenta close to the hot spots, such that the above results are questionable in the region where the logarithm is large. Moving away from the hot spot along the Fermi surface (that is, keeping ξ k = ), the self-energy varies as We δσ(k, ) = Σ(k, ) Σ(k H, ) = M N 8γ 3/2 π( + γ) 2 k 2 t 2m ln 2mΛ k 2 t. (23) This result has been extracted from a numerical evaluation of the self-energy integral. We refrained from performing a lengthy analytic derivation. Again we encounter a logarithmic

10 divergence. Due to the negative prefactor, this self-energy correction shifts the Fermi surface into the region where the bare dispersion is positive, which leads to a flattening of the Fermi surface near the hot spots. Summing the contributions from Eqs. (2) and (23), one obtains the total self-energy correction δσ(k, ) for k r and k t in the form δσ(k, ) = A sgn(ξ k) r v F k r ln k 2 t Λ v F k r + Asgn(ξ k) t 2m ln 2mΛ k 2 t, (24) where A α r A α t = M N = A α r M N [ ] + Θ(α) + K(γ), (25) 8γ 3/2 π( + γ) 2. (26) We have checked numerically that this formula describes the correct leading behavior for small k r and k t not only in radial and tangential direction, but in any direction in the momentum plane near the hot spots. Note that A α t receives a contribution from Eq. (2) with ξ k >, since the term in Eq. (2) depends on ξ k, not only on k r (see the derivation in Appendix B). While A r is always positive, the signs of A + r, A + t and A t depend on γ. A + r is positive for γ < γ c and negative for γ > γ c with γ c 5., as already discussed above. A + t is positive for γ < γ c and negative for γ > γ c with γ c.5. Hence, for γ > γ c the logarithmic divergence of the prefactor of the quadratic tangential momentum dependence flips the Fermi surface shape from convex to concave in the immediate vicinity of the hot spots. However, this effect is also most likely an artifact of applying a perturbative expression outside its range of validity. The self-energy correction renormalizes the bare dispersion ξ k to become an effective dispersion, which can be written as where ξ k = ξ k + δσ(k, ) = Z sgn(ξ k) r v F k r + Z sgn(ξ k) t 2m, (27) k 2 t Z α r Z α t = + A α Λ r ln v F k r, (28) = + A α t ln 2mΛ. kt 2 (29)

11 Logarithmic divergences are frequently a perturbative manifestation of power-law behavior. In particular, for k t =, where sgn(ξ k ) = sgn(k r ), Eqs. (28) might be just the leading order term corresponding to a power law of the form Z α r = ( ) A α Λ r. (3) v F k r Expanding Eq. (3) in powers of the exponent A α r one obtains Eq. (28). This ad hoc resummation can also be derived more systematically from a flow equation for the renormalized Fermi velocity v α F = Zα r v F. Inserting the perturbative expression Eq. (28) and taking a derivative with respect to k r yields vα F k r = v F A α r k r. Making this flow equation selfconsistent by replacing the bare Fermi velocity v F on the right hand side by the renormalized one, we obtain v F α k r = vα F A α r k r. (3) With the renormalization condition v α F (k r) = v F for v F k r = Λ, this flow equation has the unique solution v α F (k r ) = ( ) A α Λ r v F, (32) v F k r where α = sgn(k r ). The renormalized Fermi velocity thus obeys a power-law behavior as a function of k r. For k r < it diverges upon approaching k r =. For k r > it also diverges if γ < γ c, but it vanishes if γ > γ c. For k r = and k t (that is, ξ k > ) we can resum Eq. (29) in the same way to obtain a renormalized mass For γ < γ c, the exponent A + t 2 m(k t ) = ( 2mΛ k 2 t ) A + t 2m. (33) is positive such that m(k t ) vanishes upon approaching the hot spot. For γ > γ c the exponent is negative and m(k t ) diverges. In contrast to the behavior suggested by the perturbative logarithm there is no sign change of m(k t ) inducing a (rather unplausible) change of the Fermi surface shape from convex to concave. Note that A α t for all γ. < A α r Hence, the ratio of the renormalizations of tangential and radial momentum dependences always vanishes upon approaching the hot spot, indicating a flattening of the Fermi surface for any γ.

12 2 IV. SELF-CONSISTENCY AND STABILITY OF QCP The RPA calculation yields a susceptibility with maxima for certain momenta Q on the 2k F -line. These momenta determine the modulation of the density wave as obtained in mean-field theory. In the preceding section we have shown that fluctuations at the QCP give rise to a singular fermion self-energy, which leads to pronounced non-fermi liquid behavior at the hot spots. Since the square root singularity of the RPA susceptibility and effective interaction at Q is determined by states in the hot spot regions, one has to analyze to what extent this singularity is affected by the fermion self-energy. To this end, we compute the polarization function d 2 k dk Π(q, iq ) = (2π) 2 2π G(k, ik ) G(k q, ik iq ), (34) with fermion propagators G dressed by the self-energy. Since contributions with momenta near the hot spots ±k H are expected to dominate, we first approximate Σ(k, iω) by its asymptotic low-energy behavior at k H, as given by Eq. (6). Under this assumption an analytic evaluation is possible. The influence of the momentum dependence of Σ(k, iω) will be discussed later. We actually calculate the polarization function for a generalized form of the self-energy defined as where 2 3 Σ(k, iω) = icsgn(ω) ω α, (35) α < and C is a positive constant. This provides a broader view of the singularity structure and allows us to discuss the fate of the QCP for a possibly renormalized exponent α > 2 3. The Dyson equation G = G Σ yields a dressed propagator given by G(k, iω) = icsgn(ω) ω α ξ k. (36) The linear frequency term stemming from the bare propagator is subleading and has thus been discarded. In Appendix C we compute the leading momentum dependence of the static polarization function δπ(q, ) = Π(q, ) Π(Q, ) for momenta near Q. The result is a singular function of e q. For 2 3 < α < we obtain δπ(q, ) = a sgn(eq) e q α 2 + regular terms, (37)

13 3 where a < and a + > are two constants given by Eq. (C). For α and C = we recover the momentum dependence of the bare polarization function m δπ (q, ) = Θ(e q ) e q + regular terms. (38) 2πv F For α = 2/3 we find δπ(q, ) = 3 m 32πv F C e 3/2 q ln CΛ 2/3 where Λ is an ultraviolet cutoff. Inserting the prefactor C = δπ(q, ) = 37/4 N /2 64M (b + c)e 3/2 q ln e q CΛ 2/3 e q + regular terms, (39) 4M 3N( b+ c) 2/3 from Eq. (6) yields + regular terms. (4) The square-root singularity of δπ (q, ) for e q > is replaced by a power-law with a larger exponent between 2 and for 2 3 logarithmic correction for α = 2 3. The value α = 2 3 obtained for the one-loop fermion selfenergy thus defines a marginal case. For α < 2 3 < α <, and by a linear behavior with a multiplicative the singular contribution to δπ(q, ) would be super-linear and thus subleading compared to the regular linear contributions. For 2 3 α <, the fermion self-energy in the propagators generates a leading singular momentum dependence of the polarization function also for e q <, with the same functional form as for e q >. The prefactor of the singular term happens to be negative for e q <. Since it dominates over the leading regular contribution of the form be q, it spoils the minimum of Π(q, ) at q = Q and shifts it elsewhere, away from the 2k F -line. Correspondingly, the maximum of the susceptibility is also shifted away from the 2k F -line, which is inconsistent with the original assumption on the nature of the QCP. To assess the influence of the momentum dependence of the fermion self-energy, we have also computed δπ(q, ) with a momentum-dependent self-energy as obtained in Sec. III, performing a numerical evaluation of the loop-integral, where we used the perturbative, not resummed, expressions for Σ(k, ). In Fig. 2 we show results for the ratio δπ(q, )/q r as a function of q r for various choices of γ = c/b. One can see that the momentum dependence of δπ(q, ) is essentially linear as in Eq. (4), and the slope (prefactor of the linear dependence) also diverges logarithmically. However, in contrast to Eq. (4), the slope of δπ(q, ) is negative, except for small values of γ. This is simply due to the change of Fermi surface at the hot spots from convex to concave for γ > γ c.5, which is imposed by the momentum dependence of Σ(k, ) as discussed in Sec. III.B. We recall that Eq. (4) was derived with

14 FIG. 2: Numerical results for δπ(q, )/q r computed with a frequency and momentum dependent self-energy for various choices of γ = c/b. The other parameters are v F = m = b =, and the multiplicities M and N have also been set equal to one. The energy scale in the logarithms in Σ(k, ) has been chosen as Λ = 4, and the integration variables k, k r, k t in the polarization function, Eq. (34), were restricted by a UV cutoff equal to one. Note the logarithmic scale on the abscissa. a fermion self-energy depending only on frequency. The proportionality of δπ(q, ) to b at fixed γ, as obtained in Eq. (4), remains true also in the presence of Σ(k, ). The above remarks also hold for q r < (not shown in Fig. 2). The one-loop result for the susceptibility and the fermion self-energy computed with bare fermion propagators is thus far from being self-consistent. Plugging the one-loop self-energy into the polarization function yields a result that deviates drastically from the bare polarization function Π. This is in striking contrast to the case of the Ising nematic QCP, where the one-loop susceptibility and one-loop fermion self-energy computed with bare propagators are also solutions of the self-consistent one-loop equations for the same quantities computed with full propagators (including the fermion self-energy). 6,2

15 5 One might try to find a self-consistent scaling solution of the one-loop equations with a fermion self-energy scaling as ω α with an exponent α 2. However, we have shown that 3 any momentum-independent non-fermi liquid self-energy with α < will spoil the peak of the polarization function at the 2k F wave vector Q, and hence the mechanism for the density wave instability occurring naturally at a wave vector on the 2k F -line. Hence, such a scaling solution doesn t seem to exist. However, in Sec. III.B we have seen that the momentum dependence of the fermion selfenergy computed with bare polarization bubbles leads to a flattening of the Fermi surface near the hot spots. A comparison with one-dimensional Luttinger liquids 2 shows that a flat Fermi surface favors a peak in the susceptibilities at 2k F. Hence, a QCP with a non-fermi liquid self-energy and a 2k F peak in the susceptibility might be stabilized by a Fermi surface flattening around the hot spots. To describe such a QCP requires a solution for Π(q, iq ) and Σ(k, iω), where momentum and frequency dependences of both quantities are taken into account self-consistently. This is beyond the scope of this work. We finally note that vertex corrections do not affect the above results, since there are no singular one-loop vertex corrections for incommensurate density waves. 4 V. CONCLUSION In summary, we have analyzed quantum fluctuation effects at the onset of density wave order with an incommensurate 2k F wave vector Q in two-dimensional metals for the case where Q connects only a single pair of hot spots ±k H on the Fermi surface. We have confirmed the non-fermi liquid scaling proportional to ω 2/3 of the imaginary part of the one-loop fermion self-energy Σ at the hot spots as obtained previously. 7 Computing the momentum and frequency dependence of ImΣ also away from the hot spots we obtained the frequency-momentum region governed by non-fermi liquid behavior. The momentum dependence of the real part ReΣ at zero frequency yields a logarithmic divergence of the renormalized Fermi velocity and a logarithmically diverging renormalization of the tangential momentum dependence of the dispersion at the hot spots. Resumming the perturbative divergences by a simple scaling ansatz promotes the logarithms to power-laws, indicating in particular a flattening of the Fermi surface near the hot spots. Going beyond the leading order perturbation expansion, we have shown that the one-loop

16 6 result computed with bare fermion propagators is far from being self-consistent. Computing the polarization function with propagators dressed by the one-loop self-energy, one obtains a completely different behavior at the 2k F vector Q. In particular, the peak at Q is destroyed, such that the instability wave vector shifts away from the 2k F line. Moreover, a destruction of the 2k F peak found in the bare polarization function is obtained for any momentumindependent self-energy with a non-fermi liquid frequency exponent. Hence, a QCP at the onset of a density wave with an incommensurate 2k F wave vector, which is naturally obtained in mean-field theory, seems to be spoiled by fluctuations. However, a 2k F QCP might still be conceivable in a self-consistent solution with a mutual stabilization of a flattened Fermi surface and a 2k F peak in the susceptibility. Exploring this possibility is not easy and remains for future work. The stability of the 2k F density wave QCP was already discussed by Altshuler et al. 4 They found a strong non-renormalizable divergence in the two-loop contributions to the polarization function (computed with bare fermion propagators), and concluded that fluctuations will destroy the QCP in favor of a first order transition. In view of our results we believe that a continuous quantum phase transition might still occur, either at a wave vector shifted away from the 2k F line, or with 2k F fluctuations being enhanced by a flattened Fermi surface. At the one-loop level, the 2k F density wave QCP is formally equivalent to a QCP separating a normal metal from a Fulde-Ferrell-Larkin-Ovchinikov (FFLO) superconductor, where the one-loop self-energy at the hot spots also scales as ω 2/3 in two dimensions. 22 Our results thus raise the issue of the stability of the FFLO QCP. It would be worthwhile to investigate the fate of the FFLO QCP along the same lines. Acknowledgments We are grateful to Darshan Joshi for a critical reading of the manuscript. T.H. is supported by the Minerva Foundation.

17 7 Appendix A: Effective interaction In this appendix we present two concrete examples for the RPA effective interaction, namely for the charge instability in the spinless fermion model with nearest neighbor interactions, and for the spin instability in the Hubbard model.. Spinless fermion model Spinless lattice fermions with nearest-neighbor interactions exhibit a charge density wave instability, if the interaction V exceeds a certain critical value. The RPA charge susceptibility reads χ(q, iq ) = χ (q, iq ) + V (q)χ (q, iq ), (A) where χ (q, iq ) = Π (q, iq ) is the bare charge susceptibility, and V (q) = 2V (cos q + +cos q d ) is the Fourier transform of the nearest-neighbor interaction (in d dimensions). A charge density wave instability occurs at a critical interaction strength V (Q) = /χ (Q, ), where Q is the ordering wave vector. The effective interaction is obtained from the sum of bubble chains as D(q, iq ) = V (q) V (q)π (q, iq ). (A2) At the QCP, this can be approximated by the form Eq. (8) with g = V (Q). 2. Hubbard model In the Hubbard model, spin- 2 lattice fermions interact via a local interaction U between particles with opposite spin orientation. The RPA spin susceptibility is given by χ(q, iq ) = χ (q, iq ) 2 Uχ (q, iq ), (A3) where χ (q, iq ) = 2Π (q, iq ). Due to the spin-rotation invariance of the model, it is the same for all spin components. A spin density wave instability occurs at a critical interaction strength U = 2/χ (Q, ), where Q is the ordering wave vector. Due to the spin structure and the spin rotation invariance, one can construct several effective interactions by summing bubble chains or particle-hole ladders (see Fig. 3), which

18 8 D, = + + D, = + + D lad = + FIG. 3: Feynman diagrams contributing to the effective interactions in the Hubbard model. all diverge at the QCP. Summing bubble chain diagrams yields an effective interaction between particles with equal spin orientation D, (q, iq ) = U 2 Π (q, iq ) U 2 [Π (q, iq )] 2, (A4) and an effective interaction between particles with opposite spin orientation D, (q, iq ) = U U 2 [Π (q, iq )] 2. (A5) They both diverge at the QCP. Summing particle-hole ladder diagrams, one obtains the effective interaction D lad (q, iq ) = U + UΠ (q, iq ), (A6) which also diverges at the QCP. At the QCP one has UΠ (Q, iq ) =, such that one can approximate the effective interactions as D, (q, iq ) = D(q, iq ), D, (q, iq ) = D(q, iq ), D lad (q, iq ) = 2D(q, iq ), (A7) (A8) (A9) where D(q, iq ) has the form Eq. (8) with g = U/2. The effective interaction can also be written in the manifestly spin-rotation invariant form D 2 2 (q, iq ) as in Eq. (), where the τ z component is obtained from the bubble diagrams, while the τ x and τ y components come from the ladder diagrams. To first order in D, there are two distinct contributions to the fermion self-energy where the divergence of D leads to singular contributions (see Fig. 4). They involve D and D lad

19 9 Σ = + FIG. 4: Feynman diagrams representing RPA contributions to the fermion self-energy in the Hubbard model. in the form Σ(k, iω) = d 2 q dq (2π) 2 d 2 q (2π) 2 2π D lad(q, iq ) G (k q, iω iq ) dq 2π D,(q, iq ) G (k q, iω iq ). (A) (A) Note that here the contribution d 2 q (2π) 2 dq 2π U 2 Π (q, iq ) G (k q, iω iq ) appears twice, and should thus be subtracted once. However, keeping only the singular contributions at the QCP, this regular overcounting term can be ignored, and the singular terms can be combined to Σ(k, iω) = 3 This has the form Eq. (2), with M = 3. d 2 q dq (2π) 2 2π D(q, iq ) G (k q, iω iq ). (A2) Appendix B: Computation of one-loop self-energy In this appendix we derive the asymptotic results for the fermion self-energy, starting from the one-loop integral Eq. (2). We set the global multiplicity factor M equal to one for simplicity. The dispersion relations are expanded as e q = v F q r + q2 t 4m and ξ k q = v F (k r q r ) + 2m (k t q t ) 2. Combining these expansions, we can write the fermion energy as ξ k q = ξ k + e q + 4m (q t 2k t ) 2, and use q t and e q as momentum integration variables. Inserting the propagators with the above expression for ξ k q, we obtain the three-fold integral Σ(k, iω) = deq dqt dq av F 2π 2π 2π ξ k + e q iω + iq + (q 4m t 2k t ) 2 h(e q, q ) be, (B) q + c q2 t 4m

20 2 where h(e q, q ) = 2 e 2 q + q 2 + e q. The real part of Σ(k, iω) requires an ultraviolet regularization (see below). The q t integration can be performed by using the residue theorem. The integrand has four simple poles in the complex plane, q t = 2k t ± 2 m ξ k e q + iω iq, m q t = ±2i h(e q, q ) c be q. Closing the integration contour in the upper or in the lower half-plane yields the same result, but in a different form. To obtain an expression that is manifestly symmetric under k t k t, we average the two expressions to obtain Σ(k, iω) = [ i sgn(ω q ) de q dq 2πN ξk e s=± q + iω iq h(e q, q ) be [ sk q + c m t + ] 2 ξ k e q + iω iq ] + [ i h(eq, q c ) be ] 2 q skt m ξk + e q iω + iq c h(eq, q ) be. (B2) q. Imaginary part: The dominant contributions to ImΣ for small frequencies come from negative e q and q e q. 7 We define dimensionless integration variables q and ẽ q via the relations q = ω q and e q = ( ω / b) 2/3 ẽ q, and dimensionless parameters ξ k and k t via ξ k = ( ω / b) 2/3 ξk and k t = m( ω / b) /3 kt, respectively. In the low frequency limit we can approximate h(q, q ) by q / e q = b /3 ω 2/3 q / ẽ q, and ξ k e q +iω iq ξ k e q = ( ω / b) 2/3 ( ξ k ẽ q ). The neglected terms are suppressed by a factor of order ω /3. Inserting these approximations with dimensionless variables into Eq. (B2), we obtain ImΣ(k, iω) = ω b 2/3 min(, ξ k ) dẽ q 2πN s=± q ẽq ẽ q + γ [ s k t + d q sgn[sgn(ω) q ] ξk ẽ q ξ k ẽ q ] 2, (B3) with γ = c/ b = c/b. The integration over q is now elementary. Using d q sgn(± q ) A q + B = ±2 d q A q + B = ± 2 ln( + A/B), A (B4)

21 2 we obtain ImΣ(k, iω) = sgn(ω) πn with the dimensionless scaling function I( ξ k, k min(, ξk ) dẽ q t ) = ξ ln + k /ẽ q s=± ω b 2/3 I( ξ k, k t ), ( + γ s k t / ẽ q + (B5) ẽ q 3/2 ξ ) 2. (B6) k /ẽ q We now evaluate the scaling function analytically in various limiting cases. For k = k H one has ξ k = and k t =. From I(, ) = 2 dẽ q ln ( + ẽ ) q 3/2 + γ = 4π 3 ( + γ) 2/3 (B7) one obtains Eq. (6) for ImΣ(k H, iω). The same result is valid for b ξ k 3/2 ω and b k t 3 m 3/2 ω, because under these conditions one can approximate ξ k and k t by zero. For k t = and ξ k < one has k t =, ξ k <, and thus [ ] ξk dẽ q ẽ q 3/2 I( ξ k, ) = 2 ξ ln + k /ẽ + γ( ξ. (B8) k /ẽ q ) q Substituting x = ẽ q / ξ k, and expanding for ξ k, one obtains [ ] dx I( ξ k, ) = 2 ξ k ln + ξ k 3/2 x /2 x ( + γ)x γ 2 ξ k dx x ( + γ)x γ = 2π. (B9) + γ ξ k Inserting this into Eq. (B5) yields Eq. (7). For k t = and ξ k > one has k t =, ξ k >, and thus [ ] dẽ q ẽ q 3/2 I( ξ k, ) = 2 ξ ln + k /ẽ + γ( ξ. (B) k /ẽ q ) q Substituting x = ẽ q / ξ k and expanding for ξ k and ξ 3/2 k γ, one obtains [ ] dx I( ξ k, ) = 2ξ k ln + ξ 3/2 x /2 + x k ( + γ)x + γ 2 dx ξ x + ( + γ)x + γ = 4Arcsinh(γ /2 ). (B) + γ k ξ k

22 22 Inserting this into Eq. (B5) yields Eq. (8). Note that the condition ξ k alone is not sufficient to justify the above expansion. Eqs. (7) and (8) hold also for k t with b k t 3 m 3/2 ω, because under this condition one can approximate k t by zero. For ξ k = and k t we have ξ k = and k t, and thus I(, k ẽ q 3/2 t ) = dẽ q ln + ( s=± + γ + s k t / ẽ q ) 2. (B2) Substituting x = ẽ q / k t 2, and expanding for k t and k t 3 γ, one obtains I(, k t ) = k t 2 dx [ ln + ] x /2 s=± k t 3 ( + γ)x + γ + 2sγ x dx k t x ( + γ)x + γ + 2sγ x = 2π γ k t. (B3) s=± Inserting this into Eq. (B5) yields Eq. (9). The result remains valid for ξ k as long as b ξk 3/2 ω since then one can approximate ξ k by zero. 2. Real part We now compute the momentum dependence perpendicular to the Fermi surface of the self-energy at zero frequency. Setting ω = and k t = in Eq. (B2), one obtains Σ(k, ) = de q dq πn h(e q, q ) be q + c(ξ k e q iq ) [ ] isgn(q ) c + ξk e q iq h(eq, q ) be. (B4) q This is a sum of two (real) terms, which we denote by Σ (k, ) and Σ 2 (k, ). For small ξ k, the leading contributions of both terms to δσ(k, ) = Σ(k, ) Σ(k H, ) are proportional to ξ k ln ξ k. Imposing an ultraviolet cutoff e 2 q + q 2 Λ and introducing rescaled integration variables defined by e q = ξ k ẽ q and q = ξ k q, the first term can be written as Σ (k, ) = ξ ( ) k dẽ q d q Θ Λ/ ξ k ẽ πn 2 q + q 2 isgn( q ) sgn(ξk ) ẽ q i q 2 ẽ 2 q + q 2 + ẽ q + ξ k [ c(sgn(ξ k ) ẽ q i q ) bẽ ]. q (B5)

23 23 Dominant contributions for ξ k come from large (UV) and small (IR) integration variables q and ẽ q. To compute the UV contribution, it is convenient to use polar coordinates r and θ such that ẽ q + i q = re iθ. Subtracting Σ (k H, ) and letting ξ k in the integrand, one obtains δσ UV (k, ) = ξ k πn Λ/ ξk ( ) π dr dθ i r sgn(θ) π 2 cos(θ/2) sgn(ξk ) re iθ re iθ. (B6) Note that δσ UV (k, ) does not depend on b and c. For large r, the integrand decays as r, which yields a logarithmic UV divergence. The prefactor is determined by an elementary integration over θ, yielding δσ UV (k, ) = 2N ξ k ln Λ ξ k + O(ξ k). (B7) For ξ k < there is also a contribution of the same order from ẽ q and q ẽ q. Expanding 2 ẽ 2 q + q 2 + ẽ q q / ẽ q, this contribution can be written as Σ IR (k, ) = ξ λ ẽq k dẽ q d q πn λ ẽ q + ẽq q / ẽ q ξ k [ c( + ẽ q ) + bẽ q ] (B8) for small ξ k, where λ is a small positive number. The frequency integration is now elementary, leading to Using dẽ ẽq q +ẽq yielding Σ IR (k, ) = ξ k 2 ẽ q dẽ q πn ln + ẽq λ ẽq ξk [ c( + ẽ q ) + bẽ q ]. (B9) = π/2, the leading contribution for ξ k can be explicitly evaluated, Σ IR (k, ) = ξ k 2N ln Λ ξ k + O(ξ k). (B2) The choice of the arbitrary constants Λ and λ affects only the subleading term of order ξ k. We recall that Σ IR (k, ) contributes only for ξ k <. Hence, δσ (k, ) = 2N [ + Θ( ξ k)] ξ k ln Λ ξ k + O(ξ k). (B2) Making the substitution e q + iq = x + iy, the second contribution to Σ(k, ) can be written as Σ 2 (k, ) = πn Λ dx dy 4(x 2 + y 2 ) 2x b(x2 y 2 ) 2x b(x 2 y 2 ) + c[ξ k (x + iy) 2 ], (B22)

24 24 where Λ is an ultraviolet cutoff with dimension of energy. The imaginary part in the denominator does not contribute to leading order and can thus be ignored. Subtracting Σ 2 (k H, ) and performing the y-integration yields 8 Λ [ ( 2x 2( b + c)x 2 + cξ ) k δσ 2 (k, ) = dx g πn( b + c) 2x bξ g ( ( b + c)x )], k where g(u) = u arctan c/( bu). (B23) Introducing the dimensionless variable x defined by x = b ξ k x, this can be written as δσ 2 (k, ) = 8 ξ Λ b ξk k d x πn( + γ) [ ( 2 x 2( + g γ) b2 ξ k x 2 + γ sgnξ ) k g ( ( + γ) b 2 ξ k x ) ], 2 x sgnξ k (B24) with γ = c/ b. The leading contributions to this integral come from large x. Restricting the integration to x, and expanding for large x and small ξ k, one obtains δσ 2 (k, ) = 8 ξ Λ b ξk k d x g () ( + γ) sgnξ k πn( + γ) 2 x = 2 ( γ πn + γ arctan ) Λ γ ξ k ln b ξk. Summing δσ (k, ) and δσ 2 (k, ) yields Eq. (2). (B25) Appendix C: Computation of polarization bubble In the following we evaluate the polarization bubble Π(q, iq ) defined in Eq. (34) with a self-energy of the form (35). expanding the fermion dispersion yields d 2 k dk Π(q, iq ) = (2π) 2 2π i{k } v F k r k2 t 2m Using radial and tangential momentum coordinates and i{k q } + v F (k r q r ) (kt qt)2 2m. (C) Here we have introduced the short-hand notation {ω} = Csgn(ω) ω α. The k r and k t integrations can be performed by using the residue theorem, Π(q, iq ) = i dk dk t 2π 2π v F = i v F dk 2π m = 4πv F dkt 2π Θ( k ) Θ(k q ) i({k } + {k q }) v F q r k2 t (kt qt)2 2m 2m Θ(k q ) Θ( k ) (k t q t/2) 2 m + e q i({k } + {k q }) dk Θ(k q ) Θ( k ) eq + i({k } + {k q }). (C2)

25 25 Subtracting Π(Q, ) one thus obtains [ m Θ(k q ) Θ( k ) δπ(q, iq ) = dk 4πv F eq + i({k } + {k q }) dk 2i{k } ]. (C3) For the bare bubble one obtains the same expression with {k } = k. The k -integration is then convergent (without UV cutoff) and elementary, yielding the familiar singular contribution δπ (q, iq ) = m 4πv F ( eq + iq + e q iq ). (C4) We now evaluate Eq. (C3) for α < in the static limit q =. Introducing an ultraviolet cutoff Λ, we can write δπ(q, ) as m δπ(q, ) = Re 2πv F Λ Using dimensionless variables k and Λ defined via ( ) dk eq + 2i{k }. (C5) 2i{k } k = C e q k /α, Λ = C e q Λ /α, (C6) we can scale out dimensionful quantities as where δπ(q, ) = δ Π(q, ) = Λ α Re d k k = Re α 2i Λ m eq α 2 δ Π(q, ), 2πv F C α α d k k α 3 2 2i k sgn(e q ) sgn(eq) 2i k 2i k. (C7) (C8) For α > 2 3 the integral is UV convergent, that is, we can set Λ. From d k k α 3 2 = ( π 2i k ± i 2 ) α 2 ( Γ 2 α) ( Γ, (C9) α) where Γ is the Gamma function, one thus obtains Eq. (37), with the coefficients a ± = m 2πv F 2 α παc α [ π ( cos 3 ± 2 )] ( Γ 4 α 2 α) ( Γ. (C) α)

26 26 Note that a + m/(2πv F C) and a for α. For α = 2 the integral in Eq. (C8) is dominated by large k 3, so that we can expand δ Π(q, 3 Λ ) = Re 2 d k 2i sgn(eq) 2i k 3 Re 2 2i Inserting this into Eq. (C7) with α = 2 3 Λ d k sgn(e q ) 4i k = 3 6 sgn(e q) ln( Λ). (C) yields Eq. (39). H. v. Löhneysen, A. Rosch, M. Vojta, and P. Wölfle, Rev. Mod. Phys. 79, 5 (27). 2 J.A. Hertz, Phys. Rev. B 4, 65 (976). 3 A.J. Millis, Phys. Rev. B 48, 783 (993). 4 V. Oganesyan, S. A. Kivelson, and E. Fradkin, Phys. Rev. B 64, 959 (2). 5 W. Metzner, D. Rohe, and S. Andergassen, Phys. Rev. Lett. 9, 6642 (23). 6 M. A. Metlitski and S. Sachdev, Phys. Rev. B 82, 7527 (2). 7 A. Abanov, A. V. Chubukov, and J. Schmalian, Adv. Phys. 52, 9 (23); A. Abanov and A. V. Chubukov, Phys. Rev. Lett. 93, (24). 8 M. A. Metlitski and S. Sachdev, Phys. Rev. B 82, 7528 (2); New J. Phys. 2, 57 (2). 9 C. Castellani, C. Di Castro, and M. Grilli, Phys. Rev. Lett. 75, 465 (995); F. Becca, M. Tarquini, M. Grilli, and C. Di Castro, Phys. Rev. B 54, 2443 (996); A. Perali, C. Castellani, C. Di Castro, and M. Grilli, Phys. Rev. B 54, 626 (996). H. J. Schulz, Phys. Rev. Lett. 64, 445 (99). P. A. Igoshev, M. A. Timirgazin, A. A. Katanin, A. K. Arzhnikov, and V. Yu. Irkhin, Phys. Rev. B 8, 9447 (2). 2 S. Sachdev and R. La Plata, Phys. Rev. Lett., 2722 (23). 3 T. Holder and W. Metzner, Phys. Rev. B 85, 653 (22). 4 B. L. Altshuler, L. B. Ioffe, and A. J. Millis, Phys. Rev. B 52, 5563 (995). 5 D. Bergeron, D. Chowdhury, M. Punk, S. Sachdev, and A.-M. S. Tremblay, Phys. Rev. B 86, 5523 (22). 6 Y. Wang and A. Chubukov, Phys. Rev. B 88, 2456 (23).

27 27 7 T. Holder and W. Metzner, Phys. Rev. B 9, 66(R) (24). 8 J. W. Negele and H. Orland, Quantum Many-Particle Systems (Addison-Wesley, Reading, MA, 987). 9 F. Stern, Phys. Rev. Lett. 8, 546 (967). 2 This was first shown for the formally equivalent problem of electrons coupled to a U() gauge field by J. Polchinski, Nucl. Phys. B 422, 67 (994). 2 T. Giamarchi, Quantum Physics in One Dimension, (Oxford University Press, Oxford, 24). 22 F. Piazza, W. Zwerger, and P. Strack, Phys. Rev. B 93, 852 (26).