arXiv:cond-mat/0003349v1 [cond-mat.str-el] 21 Mar 2000 w aallfltFrisraepee 6.Znh and functional Zanchi modern with how system shown [6]. model recently pieces a have surface for [7] Fermi par- Schulz study flat in RG parallel demonstrated 1-loop the two been in a has important in This and ticular strong limit. is low-energy surface Fermi the along the since momen- neglected, gential irrelevant are small dependences if energy a inter- even and tum couplings, by effective running accurately the of parameterized that number be is cannot dimension actions one to compared correlation corresponding the - of identi- functions. level. divergencies where from 1-loop instabilities on fied superconducting couplings computed and running and density of defined in number be points small Hove could a which van dom- for between on k-space, focussed processes studies scattering These inant [3–5]. perconductivity n sepce obcm a become an to expected has half-filling is it at and state since ground the ordered planes superconductors, antiferromagnetically for copper-oxide the high-temperature model in prototype of freedom of promising degrees a electronic as interest much ecnutr srte tog h ednytwrsan- towards tendency and the tiferromagnetism strong, rather is perconductors [2]. concentrations smaller slightly oe tre ihsml u neiu cln ap- scaling high- of ingenious discovery the but after simple shortly very with proaches, started model channels footing. equal particle-hole an and on particle-particle the account adrsmain fFymndarm,btrequires but stan- by diagrams, Feynman achieved a of be cannot resummations divergencies dard van these treat- and unbiased of nesting and ment controlled surface A Fermi [3–5]. of singularities consequence Hove ap- a divergences as infrared pear competing where close half-filling, densities for to at down breaks model theory perturbation Hubbard tional 2D the by eomlzto group ao opiaini w-iesoa systems two-dimensional in complication major A h w-iesoa ubr oe 1 a attracted has [1] model Hubbard two-dimensional The lhuhteCuobitrcini h urt su- cuprate the in interaction Coulomb the Although al Gsuiso h w-iesoa Hubbard two-dimensional the of studies RG Early oetmdpnec feetv interactions effective of dependence momentum AS 11.d 11.w 74.20.Mn 71.10.-w, 71.10.Fd, PACS: surface. instab Fermi Pomeranchuk the a to of lead unexpectedly may which More develop interactions fluctuations. spin antiferromagnetic by ae nanwrnraiaingopmto o interacting ex for expected the method confirms group susceptibilities and renormalization interactions new a on based epeetasseai tblt nlssfrtetwo-dime the for analysis stability systematic a present We d d wv uecnutvt n oeacu instability Pomeranchuk and superconductivity -wave wv arn scpue already captured is pairing -wave R)aayi hc ae into takes which analysis (RG) ntttfu hoeicePyi ,Tcnsh Hochschul Technische C, Physik f¨ur Theoretische Institut d wv uecnutrfor superconductor -wave weak ntetodmninlHbadmodel Hubbard two-dimensional the in epegae 5 -25 ahn Germany Aachen, D-52056 55, Templergraben hitp .HlohadWle Metzner Walter and Halboth J. Christoph opig Conven- coupling. Nvme ,2018) 2, (November T c tan- su- 1 nasur atc.Teepce xsec fa of existence expected The hopping lattice. neighbor square next-nearest a Hubbard and on two-dimensional nearest the with to model renormal- method Salmhofer’s group applying by ization obtained as tibilities susceptibil- have for [9]. we equations ities flow instabilities, numerical corresponding physical concrete the derived compute equa- a To flow for resulting suitable evaluation. The particularly are re- of approach. tions version Most theoretic improved an field derived surfaces. this has Fermi [8] curved Salmhofer cently, inter- arbitrary effective for the of treat actions dependence to momentum used tangential be full can methods group renormalization cin eeo hc a edt oeacu [10] Pomeranchuk Fermi the a surface. of symmetry to tetragonal the lead breaking may inter- instability un- scattering which forward More develop strong actions that doubt. find beyond we confirmed expectedly, thereby fluc- spin is antiferromagnetic tuations by driven instability pairing hc ed oa niermgei ntblt o arbi- for instability small antiferromagnetic trarily an to leads which we etnaetnihoso qaeltielasto leads relation lattice dispersion square the a on neighbors next-nearest tween est fsae tteenergy the at states non-interacting of the density in singularities Hove van logarithmic ecie ih-idn lcrn ihalclrepulsion local a > with U tight-binding describes atc site lattice a ǫ prtr o emoswt pnprojection spin with fermions for operators o igepril tts hsdseso eainhas relation dispersion at This points saddle states. single-particle for − k t nti etrw rsn eut o h o fsuscep- of flow the for results present we letter this In h n-adHbadmdl[1] model Hubbard one-band The a h etn property nesting the has ewe ers egbr n namplitude an and neighbors nearest between sec fa of istence .Here 0. ǫ efidta togfradscattering forward strong that find we , k lt raigtettaoa symmetry tetragonal the breaking ility = H em ytm.Teflwo effective of flow The systems. Fermi − soa ubr oe,wihis which model, Hubbard nsional = > U 2 i and , d t c X (cos wv arn ntblt driven instability pairing -wave i i † , σ j k thl-lig[1]. half-filling at 0 and X σ k (0 = n x j Aachen e t σ cos + ij c i = σ c π , i † σ r rainadannihilation and creation are c n ( and ) c j † k ǫ j σ σ k y c + ) + j σ Q − ǫ opn amplitude hopping A . U vH = 4 π, X t ′ j − 4 = ) hc generate which 0), cos ǫ n k j ↑ k t for ′ n x For . j σ ↓ cos Q {↑ ∈ , k ( = y − t d , ′ -wave ↓} t ,π π, 0, = ′ be- (1) (2) on ), The RG equations are obtained as follows (for details, 100 see Salmhofer [8] and Ref. [9]). The infrared singularities sdw (π,π) are regularized by introducing an infrared cutoff Λ > 0 sdw (π,π−δ) into the bare such that contributions from sdw (π−δ,π−δ) momenta with |ǫk − µ| < Λ are suppressed. All Green cdw (π,π) functions of the interacting system will then flow as a 10 sc dx2-y2 function of Λ, and the true theory is recovered in the sc dxy 0 sc s limit Λ → 0. Salmhofer [8] has recently pointed out that χ /

(amputated) Green functions obtained by expanding the χ sc xs effective action of the theory in powers of normal ordered monomials of fermion fields obey differential flow equa- 1 tions with a structure that is particularly convenient for a power counting analysis to arbitrary loop order. With the bare interaction as initial condition at the highest scale Λ0 = max |ǫk − µ|, these flow equations determine 0.1 the exact flow of the effective interactions as Λ sweeps 0.01 0.1 1 over the entire Brillouin zone down to the . Λ/t The effective low-energy theory can thus be computed FIG. 1. The flow of the ratio of interacting and directly from the microscopic model without introducing non-interacting susceptibilities for t′ = −0.01t, U = t and any ad hoc parameters. µ = −0.055t. For a weak coupling stability analysis it is sufficient to truncate the exact hierarchy of flow equations at 1- pair susceptibilities with form factors [2] loop level. The effective 2-particle interaction then re- duces to the one-particle irreducible 2-particle vertex ΓΛ, Λ 1 (s-wave) and its flow is determined exclusively by Γ itself (no  1 √2 (cos kx + cos ky) (extended s-wave) higher many-particle interactions enter). Flow equations d(k)=  1 (3)  (cos kx − cos ky) (d-wave dx2 y2 ) for susceptibilities are obtained by considering the exact √2 − sin k sin k (d-wave d ). RG equations in the presence of suitable external fields,  x y xy which leads to an additional 1-particle term in the bare  interaction, and expanding everything in powers of the Some of these susceptibilities diverge together with the external fields to sufficiently high order [9]. vertex function at the scale Λc. Depending on the choice One cannot solve the flow equations with the full en- of U, t′ and µ, the strongest divergence is found for the ergy and momentum dependence of the vertex function, commensurate or incommensurate spin susceptibility or Λ since Γ has three independent energy and momentum for the pair susceptibility with dx2 y2 symmetry. In Fig. − variables. The problem can however be much simpli- 1 we show a typical result for the flow of susceptibili- fied by ignoring dependences which are irrelevant in the ties as a function of Λ. Note the threshold at Λ ≈ 0.03t low energy limit, namely the energy dependence and the below which the amplitudes for various scattering pro- momentum dependence normal to the Fermi surface (for cesses, especially umklapp scattering, renormalize only details, see Ref. [9]). This approximation is exact for the very slowly. The flow of the antiferromagnetic spin sus- bare Hubbard vertex, and asymptotically exact in the ceptibility is cut off at the same scale. The pairing sus- low-energy regime. The remaining tangential momen- ceptibility with dx2 y2 -symmetry is obviously dominant − tum dependence is discretized for a numerical evaluation. here (note the logarithmic scale). Following the flow of Most of our results where obtained for a discretization the susceptibilities one can see that the dx2 y2 -pairing − with 16 points on the Fermi surface (yielding 880 ”run- correlations develop in the presence of pronounced but ning couplings”), and we have checked that increasing the short-range antiferromagnetic spin-correlations, in agree- number of points does not change our results too much. ment with earlier ideas on d-wave superconductivity in We have computed the flow of the vertex function for the Hubbard model [2]. many different model parameters t′ and U (t just fixes the In Fig. 2 we show the (µ,U) phase diagram for t′ = absolute energy scale) and densities close to half-filling. −0.01t obtained by identifying the dominant instability In all cases the vertex function develops a strong momen- from the flow for many different values of µ and U. For tum dependence for small Λ with divergencies for several µ =4t′ the Fermi surface touches the saddle points (0, π) momenta at some critical scale Λc > 0, which vanishes and (π, 0), while µ =4t′+0.01t corresponds to half-filling. exponentially for U → 0. To see which physical insta- Note that for U → 0 the pairing instability always dom- bility is associated with the diverging vertex function we inates, because the BCS channel dominates the flow in have computed commensurate and incommensurate spin the limit Λ → 0. A spin density wave is the leading in- susceptibilities χS(q) with q = (π, π), q = (π − δ, π) and stability for U → 0 only in the special case with perfect q = (1−δ)(π, π), where δ is a function of density [11], the nesting, t′ = 0 and µ = 0 (see the (µ,U) phase diagram commensurate charge susceptibility χC(π, π), and singlet computed from the 1-loop flow for t′ = 0 in Ref. [9]).

2 2.5 π π

2 free free 0 deformed 0 deformed 1.5

U/t −π −π 1 (a)−π0 π (b) −π 0 π FIG. 4. Schematic plot of Fermi surface deformations 0.5 sdw breaking the square symmetry; the deformed surface may be sc dx2-y2 closed (a) or open (b). 0 -0.1 -0.05 0 0.05 0.1 Previous investigations of Fermi surface deformations (4t'-µ)/t within standard perturbation theory have yielded only ′ FIG. 2. The (µ, U) phase diagram for t = −0.01t near very small shifts even for sizable interaction strengths half-filling (marked by the dashed vertical line); the symbols [12], but in these studies the possibility of a spontaneous represent the parameter values for which the flow has been breaking of the point group symmetry of the square lat- computed and whether the dominant instability is magnetic tice has not been taken into account. (squares) or superconducting (circles); the solid line sepa- To analyze systematically the stability of the Fermi rates the spin-density wave regime from the superconducting surface shape, we define a susceptibility regime. δskF ′ κkF k = (4) F ′ δµkF

0.08 which measures the size of Fermi surface shifts δskF for sdw small momentum dependent shifts of the chemical po- 0.07 tential δµk′ at points k′ on the Fermi surface. The sc dx2-y2 F F ′ 0.06 matrix κkF kF defines a linear integral operator acting on functions of k . A simple consideration in the spirit 0.05 F of phenomenological shows that the t / c 0.04 corresponding inverse operator is given by Λ 1 0.03 ′ c (κ ) = v F δ(k − k )+2f ′ (5) − kF kF k F F′ kF kF 0.02 c where vkF is the Fermi velocity and f ′ is the Landau 0.01 kF kF function in the charge (spin-symmetric) channel. It is ′ 0 now obvious that the matrix κkF kF is symmetric. The -0.08 -0.04 0 0.04 0.08 1 Fermi surface is stable, if all eigenvalues of κ (or κ− ) (4t'-µ)/t are positive. Note that Landau’s energy functional can 1 FIG. 3. The critical energy scale Λc as a function of the be written as a quadratic form in δskF , with κ− as ker- chemical potential µ for U = 1.5t and t′ = −0.01t. The dif- nel [13], and negative eigenvalues would imply that this ferent symbols indicate whether the leading instability is a energy can be lowered by a suitable deformation of the spin-density wave or d-wave pairing instability. Fermi surface. In isotropic Fermi liquids such instabil- ities occur for strongly negative Landau parameters, as How the critical energy scale Λc varies as a function first pointed out by Pomeranchuk [10]. of the chemical potential (i.e. as a function of density) We have computed the renormalization group flow of 1 is shown in Fig. 3 for an interaction strength U = 1.5t. the eigenvalues and eigenvectors of the operator κ− from cΛ the flow of the Landau function f ′ , which is given di- Obviously Λc is maximal for a chemical potential at the kF kF van Hove energy. Note that Λc must not be interpreted rectly by the vertex function in the forward scattering as a transition temperature for spin density wave forma- channel [14]. For various choices of the model parame- tion or superconductivity, but rather as an energy scale ters we have always found a Fermi surface instability at a P where bound particle-hole or particle-particle pairs are scale Λc above the scale Λc where the vertex function di- formed. verges. In all cases the instability favors a deformation of Since some of the forward scattering interactions grow the Fermi surface which breaks the point group symme- strong for small Λ, while the Fermi velocity is very small try of the square lattice, as shown schematically in Fig. near the saddle points, the Fermi surface may be signifi- 4. The instability is mainly driven by a strong attrac- cantly deformed by interactions, especially for µ ≈ ǫvH . tive interaction between particles (or holes) on opposite

3 corners of the Fermi surface near the saddle points and In summary, we have shown that modern renormal- a repulsive interaction between particles on neighboring ization group methods can be used to establish the ex- corners. pected d-wave pairing instability in the two-dimensional The above diagnosis of Pomeranchuk instabilities Hubbard beyond doubt. Note that for small bare inter- would be rigorous for a normal Fermi liquid with finite actions and in a parameter regime where only particle- renormalized interactions in the infrared limit. In the particle pairing fluctuations grow strong, the strong cou- present system, however, the vertex function diverges at pling problem associated with the formation of a super- a finite scale and possible Pomeranchuk type instabilities conducting state can be treated rigorously [16]. Further- compete with magnetic and superconducting instabili- more, we have pointed out that a Pomeranchuk instabil- ties. Since we have no quantitative theory of the strong ity breaking the tetragonal symmetry of the Fermi sur- coupling physics near and below the scale Λc, we can only face is likely to occur for a suitable choice of the model list and discuss two possible scenarios: parameters, especially for a Fermi surface close to the i) Energy gaps due to particle-particle or particle-hole saddle points in the absence of perfect nesting. binding may stop the flow of forward scattering interac- tions before a Pomeranchuk instability sets in. Acknowledgments: ii) The Pomeranchuk instability is not blocked by bind- We are very grateful to Maurice Rice, Manfred Salmhofer ing phenomena. In that case one would have a finite and Eugene Trubowitz for valuable discussions. This temperature phase transition with a spontaneous break- work has been supported by the Deutsche Forschungs- ing of the (discrete) tetragonal symmetry of the square gemeinschaft under Contract No. Me 1255/4 and within lattice, and subsequent continuous symmetry breaking the Sonderforschungsbereich 341. associated with magnetic order or superconductivity in the ground state. Which of the two scenarios is realized depends on the choice of the model parameters. The Pomeranchuk insta- bility occurs more easily if the Fermi surface is close to the saddle points of ǫk. On the other hand, nesting raises the scale for particle-hole binding (leading ultimately to [1] See, for example, The Hubbard Model, ed. A. Montorsi magnetic order). The best candidate is therefore the (World Scientific 1992). Hubbard model with a sizable t (reducing nesting) and [2] For a review on d-wave superconductivity, see D. J. ′ Scalapino, Phys. Rep. 250, 329 (1995). µ = ǫ . vH [3] H. J. Schulz, Europhys. Lett. 4, 609 (1987). We emphasize that the Pomeranchuk instability does [4] I. Dzyaloshinskii, Sov. Phys. JETP 66, 848 (1987). not cut off the singularity in the Cooper channel since [5] P. Lederer, G. Montambaux, and D. Poilblanc, J. Phys. it does not break the reflection invariance. Hence, at (Paris) 48, 1613 (1987). sufficiently large doping away from half-filling, d-wave [6] A. T. Zheleznyak, V. M. Yakovenko, and I. E. superconductivity will set in in any case, with an or- Dzyaloshinskii, Phys. Rev. B 55, 3200 (1997). der parameter that may be slightly distorted away from [7] D. Zanchi and H. J. Schulz, Z. Phys. B 103, 339 (1997); perfect d-wave symmetry. The Pomeranchuk instability Europhys. Lett. 44, 235 (1998); cond-mat 9812303. would also not destroy the umklapp scattering route to [8] M. Salmhofer, Commun. Math. Phys. 194, 249 (1998); an insulating spin liquid discussed recently by Furukawa Renormalization (Springer 1998). et al. [15]. [9] C.J. Halboth and W. Metzner, Phys. Rev. B 61, 4364 To our knowledge a Pomeranchuk instability has not (2000). yet been observed in numerical solutions of the two- [10] I.J. Pomeranchuk, Soviet. Phys. JETP 8, 361 (1958). 64 dimensional Hubbard model. Of course this may be [11] See, for example, H. J. Schulz, Phys. Rev. Lett. , 1445 due to finite size limitations or too high temperatures (1990), where also the choice of δ as a function of doping in Monte Carlo simulations. It would thus be interest- is discussed. [12] C.J. Halboth and W. Metzner, Z. Phys. B 102, 501 ing to compute the Fermi surface susceptibility κ ′ kF kF (1997), and references therein. numerically. [13] P. Nozi`eres, Theory of interacting Fermi systems (Ben- In real systems a Pomeranchuk instability as in Fig. jamin 1964). 4 may lead to an orthorhombic lattice distortion, as a [14] For the relation between RG and Fermi liquid theory see, consequence of the coupling of electronic and lattice de- for example, W. Metzner, C. Castellani and C. Di Castro, grees of freedom. High temperature superconductors in- Adv. Phys. 47, No. 3 (1998). deed exhibit structural phase transitions between tetrag- [15] N. Furukawa, T. M. Rice, and M. Salmhofer, Phys. Rev. onal and orthorhombic phases. It would be interesting Lett. 81, 3195 (1998); C. Honerkamp, M. Salmhofer, N. to clarify whether a Pomeranchuk instability might drive Furukawa, and T.M. Rice, cond-mat/9912358. (at least to a significant extent) the transition into the [16] J. Feldman, J. Magnen, V. Rivasseau, and E. Trubowitz, orthorhombic phase in these materials. Europhys. Lett. 24, 437 (1993).

4