Supermetal

We study the effect of electron interaction in an electronic system with a high-order Van Hove singularity, where the density of states shows a power-law divergence. Owing to scale invariance, we perform a renormalization group (RG) analysis to find a nontrivial metallic behavior where various divergent susceptibilities coexist but no long-range order appears. We term such a metallic state as a supermetal. Our RG analysis reveals noninteracting and interacting fixed points, which draws an analogy to the $\phi^4$ theory. We further present a finite anomalous dimension at the interacting fixed point by a controlled RG analysis, thus establishing an interacting supermetal as a non-Fermi liquid.


I. INTRODUCTION
A Bloch electron in a crystal is described by the energy dispersion E k that relates the energy with its wave vector k. For metals, the energy dispersion determines the density of states (DOS) at the Fermi level, which to a large extent governs various thermodynamic properties such as charge compressibility, spin susceptibility, and specific heat. Van Hove's seminal work [1] revealed that the DOS exhibits non-analyticity at an extremum or a saddle point of the energy dispersion, where ∇ k E k = 0. Importantly, Van Hove singularities (VHS) are guaranteed to exist in every energy band by the continuity and the periodicity of E k over the Brillouin zone. The behavior of the DOS at a VHS depends on whether it is at an energy extremum or a saddle point, and also on the dimensionality of the system. For example, at a saddle point in two dimensions with E k = k 2 x − k 2 y , the DOS diverges logarithmically. As the chemical potential crosses the VHS, the topology of Fermi surface changes from electron to hole type, known as an electronic topological transition.
Recently, we have extended the notion of VHS to highorder saddle points, where, besides ∇ k E k = 0, the Hessian matrix D ij = ∂ ki ∂ kj E k satisfies det D(k) = 0 [2]. These high-order saddle points occur where two Fermi surfaces touch tangentially, or at the common intersection of three or more Fermi surfaces [3,4]. An example of the former is E k = k 2 x − k 4 y , and of the latter is E k = k 3 x − 3k x k 2 y . Generally speaking, high-order saddle points can be realized by tuning the energy dispersion with one or more control parameters. At high-order saddle points in two dimensions, the DOS shows a power-law divergence D(E) ∝ |E| − , much stronger than a logarithmic one at ordinary VHS [2,4].
The existence of high-order VHS has recently been identified in a variety of materials including twisted bilayer graphene near a magic angle, trilayer graphene-hexagonal boron nitride heterostructure [2], and Sr 3 Ru 2 O 7 [5]. In particular, a power-law divergent DOS of high-order VHS with exponent −1/4 was found in scanning tunneling spectroscopy measurements [6] on magic-angle twisted bilayer graphene [2].
In the presence of electron-electron interaction, a large DOS near the Fermi level may have important consequences. On the one hand, it may trigger Stoner instability to ferromagnetism. On the other hand, a large DOS may result in strong screening of repulsive interaction, so that a Fermi liquid description remains valid at low energy. For the case of a single conventional VHS with a logarithmically divergent DOS at the Fermi energy, previous works [7][8][9][10][11][12][13][14][15][16][17][18][19][20][21] have shown that repulsive interaction decreases at low energies, likely leading to a marginal Fermi liquid [22][23][24][25][26].
In this work, we study interacting electron systems with a high-order saddle point near the Fermi level. Assuming that electron interaction is weak, dominant contributions to low-energy thermodynamic properties of the system come from those states in the vicinity of the saddle point, from which the DOS divergence originates. This allows us to formulate a continuum field theory of interacting fermions by taking the leading-order energy dispersion relation E k near the saddle point and extending the range of momentum to infinity.
In this field theory, when the high-order VHS is right at the Fermi level, the Fermi surface in k-space becomes scale-invariant. As the VHS approaches the Fermi level, charge and spin susceptibilities exhibit power-law divergence, reminiscent of critical phenomena. Motivated by these observations, we develop a renormalization group (RG) theory for interacting fermions near high-order VHS, which parallels Wilson-Fisher RG approach to the φ 4 theory [27,28]. By introducing a small parameter associated with the DOS divergence, we present a controlled RG analysis and find that short-range repulsive interaction is relevant at the noninteracting fixed point and drives the system into a nontrivial T = 0 interacting fixed point. The former is the analog of the Gaussian fixed point in Fermi systems, and the latter the analog of the Wilson-Fisher fixed point.
The metallic state at the interacting fixed point exhibits scale-invariance in space/time and power-law divergent charge and spin susceptibility, but finite pairing susceptibility. In other words, this is a metal on the verge of charge separation and ferromagnetism. We call such a critical state of metal with various divergent susceptibilities but without any long-range order, a supermetal. This terminology is motivated by a comparison with a arXiv:1905.05188v2 [cond-mat.str-el] 12 Aug 2019 metal and a semimetal. All three are conductors without a band gap at the Fermi level, but differ in the DOS. A semimetal has a vanishing DOS, a metal has a finite DOS, and a supermetal has a divergent DOS.
We further show by a two-loop RG calculation for a high-order saddle point that the fermion field acquires a finite anomalous dimension. Hence the interacting supermetal we found is a non-Fermi liquid, as opposed to a marginal Fermi liquid for the case of a conventional VHS. The singular DOS of supermetal D(E) ∝ |E| − plays a pivotal role by making a non-Fermi liquid possible under weak repulsive interaction. The DOS exponent naturally serves as a small parameter that allows a controlled analysis via perturbative RG calculation.
The outline of the paper is as follows: In Sec. II, we introduce a tight-binding model with a high-order VHS and calculate the power-law divergent DOS, whose exponent is determined from the scaling property of energy dispersion near the high-order saddle point. We show that a high-order VHS appears generically when the energy dispersion around a saddle point is modified by changing just a single hopping parameter.
In Sec. III, we present a mean-field analysis of interacting electrons with a high-order saddle point near Fermi level. We find that in the presence of repulsive contact interaction, as the chemical potential approaches the Van Hove energy, a first-order transition to a ferromagnetic metal occurs, displaying a discontinuous change in spin polarization and charge density.
In Sec. IV, we perform the energy-shell RG analysis step by step. We first define the energy shell as a region of momentum space. Then, the tree-level and oneloop RG equations for the chemical potential and interaction strength are derived in sequence, which resembles the case of the φ 4 theory. We identify the noninteracting fixed point and the nontrivial interacting fixed point, which is the analog of the Wilson-Fisher fixed point in Fermi system. We next consider other relevant perturbations to the system, including Zeeman and pairing fields as well as additional symmetry-allowed terms in the energy dispersion. A discussion about a higher-loop RG analysis follows, while an actual two-loop calculation appears in a later section.
In Sec. V, we combine the results from the mean-field and RG analyses to propose a phase diagram of interacting electrons near a high-order VHS in the parameter space of chemical potential, interaction strength, and detuning of single-particle energy dispersion from the highorder VHS. We show that a supermetal appears on a line in the phase diagram, which can be reached by tuning two parameters. We then perform the scaling analysis for thermodynamic quantities and correlation functions. The generic formalism is first presented, followed by the one-loop result for various exponents of divergent susceptibilities. In addition, we discuss the Ward identity, which results from charge conservation and gives relations among the field renormalization and scaling exponents.
In Sec. VI, we introduce another RG scheme, the field theory approach with a soft UV energy cutoff, which is confirmed to satisfy the Ward identity. Compared to the energy-shell RG analysis, it has the advantage in calculating higher-order perturbative corrections. The oneloop calculation reproduces the energy-shell RG analysis in Sec. IV. Furthermore, the two-loop calculation shows the finite anomalous dimension of the fermion field at a high-order saddle point. This result directly establishes the non-Fermi liquid nature of an interacting supermetal.
In Sec. VII, we evaluate the quasiparticle lifetime at finite temperature due to electron interaction. From a perturbative calculation, we find an unusual temperature dependence in the quasiparticle lifetime, which also implies the non-Fermi liquid behavior.

II. MODEL
A. An example of high-order VHS in two dimensions We consider a tight-binding model on an anisotropic square latticê H = − j t x c † j+x c j + t y c † j+ŷ c j + t y c † j+2ŷ c j + H.c.
(1) t x and t y are the nearest-neighbor hopping amplitudes along the x and y directions, respectively, and t y is the second-nearest neighbor hopping along the y direction. The energy dispersion is obtained as E k = −2t x cos(k x a) − 2t y cos(k y a) − 2t y cos(2k y a), (2) with the lattice constant a.
For |t y | ≥ 4|t y |, there are four VHS points in the Brillouin zone at the high symmetry points: Γ = (0, 0), X = (π/a, 0), Y = (0, π/a), and M = (π/a, π/a). With t x , t y , t y > 0, the energy minimum and maximum are located at Γ and M points, respectively, and X and Y points are the saddle points [ Fig. 1(a)]. Near Y point, the energy dispersions takes the form where k x and k y are rescaled to eliminate the coefficients of k 2 x and k 4 y . The evolution of the Fermi surface by changing λ is shown in Fig. 2. For λ > 0 (t y > 4t y ), there is an ordinary FIG. 1: Lattice model for a high-order VHS. (a) Energy contour plot with ty/tx = 0.8 and t y /tx = 0.2. We find the energy minimum at Γ, the maximum at M , and the two saddle points at X and Y . Y and M are high-order VHS points. At Y , we see that the two Fermi surfaces touch tangentially while they cross linearly at X. saddle point with a logarithmically divergent DOS, where the Fermi surfaces cross at a point at µ = 0. At λ = 0 (t y = 4t y ), the two Fermi surfaces touch tangentially to realize a high-order saddle point. For λ < 0 (t y < 4t y ), the singular point splits into two saddle points and one minimum. Those saddle points are located at (k x , k y ) = (0, ± |λ|/2) with the energy λ 2 /4. We can see that the high-order VHS is realized around the conventional VHS point(s) by controlling the single parameter λ in the energy dispersion [2]. We add that, at t y = 4t y , the energy dispersion near M point becomes k 2 x + k 4 y , describing a high-order energy extremum.
The specific tight-binding model (1) illustrates a general feature of Bloch electron's energy dispersion: the existence of saddle points is mathematically guaranteed [1], and tuning a single parameter can turn an ordinary saddle point into a high-order one [2].
A VHS manifests itself as an analytic singularity in the DOS The DOS for the present model (d = 2) is depicted in Fig. 1(b). We find four singularities in the DOS and each of them is tied to the individual VHS of the model. The band bottom at Γ gives rise to a discontinuity in the DOS and the saddle point at X shows a logarithmic divergence in the DOS. Those two are conventional VHS, known since the original work of Van Hove [1]. Here we focus on the high-order VHS at Y and M . They exhibit distinct behavior: the DOS has a power-law divergence as |E| −1/4 instead of a logarithm. In addition, the divergence at Y is stronger on the electron side by the factor √ 2 than on the hole side. Such an asymmetry is not seen for a conventional VHS with . The parameter λ describes a perturbation around a high-order saddle point. The sign of λ controls the topology of the Fermi surface; there is one VHS point at (kx, ky) = (0, 0) for λ ≥ 0, and it splits into two located at (0, ± |λ|/2) for λ < 0.
a logarithmic divergence at X. The two Fermi surfaces touch tangentially at Y at the Van Hove energy. When the chemical potential µ crosses the Van Hove energy, the Fermi surface topology changes from being closed to open in the k y direction.
In Fig. 1(b), the DOS peaks at the two high-order VHS in our tight-binding model are fitted by the analytical expressions of the DOS calculated from the continuum models in their vicinities. The calculation will be shown in the next subsection. We can see a close fit within a finite energy range. Since the divergent DOS and hence susceptibilities originate from the vicinity of the highorder VHS, the continuum model is expected to capture universal features at low energy. Using the continuum model has the advantage of removing non-universal aspects associated with high-energy regions away from the high-order VHS in the tight-binding model. We will show that infrared (IR) scaling properties are not indeed affected by the UV cutoff in the continuum model.
Before proceeding, we briefly mention the Fermi surfaces in strained Sr 2 RuO 4 [64][65][66]. It has a quasitwo-dimensional electronic structure with a layered perovskite structure. Under uniaxial pressure, a Lifshitz transition occurs on the Brillouin zone boundary [65]. At the transition point, there is one VHS in the Brillouin zone at the Fermi energy. The Fermi surface of the band of interest resembles the one obtained from Eq. (2).

B. Generalization
From now on, we study a continuum model of fermions with a high-order energy dispersion. For the purpose of a controlled RG analysis later, here we consider the generalized energy dispersion in the d-dimensional k-space The momentum is denoted by where k ± are d ± -dimensional vectors with d + + d − = d, and k ± = |k ± |. Analyticity of the energy dispersion requires n ± to be positive integers. We consider the case of even n ± , so that E k = E −k satisfies time-reversal symmetry. When at least one of n ± is greater than two, this energy dispersion has a high-order VHS at k = 0, which is defined as a point where the Hessian matrix The energy dispersion Eq. (5) follows the scaling relation It then follows from Eqs. (4) and (7) that the DOS satisfies where the DOS singularity exponent is Throughout this work, we consider the case > 0. For example, the high-order VHS introduced in the preceding section corresponds to the case of d + = d − = 1, n + = 2, n − = 4, so that = 1/4. We calculate the prefactors D ± for the dispersion Eq. (5) explicitly and find with the common factor (10b) We note that in calculating the DOS, the d-dimensional momentum integral over k ∈ (−∞, ∞) d is convergent for all E = 0. Also, note that D + = D − for d + /n + = d − /n − . It describes the asymmetry in the DOS above and below E = 0. This is a feature of the high-order saddle points defined by Eq. (5), distinct from conventional saddle points in two dimensions, where the logarithmically divergent DOS peak is symmetric. We also find it useful to consider another generalization with d + = 1 and d − = 2−4 . The original problem in two dimensions corresponds to = 1/4, while the generalized problem is defined in 3−4 dimensions, in a similar spirit as Wilson-Fisher theory in 4 − dimension. Now, the DOS D(E) has a power-law divergence at E = 0 for > 0 with the same form as Eq. (8), but the coefficients are replaced with The nontrivial interacting fixed point to be shown later is controlled by the smallness of . For the model defined by Eq. (5), the exponent can be any rational number between 0 < < 1. By choosing positive integers n ± and d ± judiciously, we can make arbitrarily small in high-dimensional crystals, while keeping the energymomentum dispersion an analytic function.
We now introduce our model of interacting electrons near a high-order VHS: with the density operatorn rσ = c † rσ c rσ . g denotes the coupling constant for the contact interaction between electrons with opposite spins and the summation over the spin index σ =↑ (+), ↓ (−) is implicit. The corresponding action is given by with the fermionic field ψ σ . We set k B = = 1 throughout the paper. Here we formulate the model at temperature T . Temperature T is regarded as the system size L β ≡ 1/T in the imaginary time direction. Later, we shall consider the effects of other interactions and external fields.
From the action, we define the noninteracting Green's function with the fermionic Matsubara frequency ω n = (2n+1)πT (n: integer). The partition function Z is expressed as We are interested in thermodynamic quantities such as specific heat. These are obtained from the free energy density where V is the volume of the system.

III. MEAN-FIELD ANALYSIS
We first consider the effect of interaction in Eq. (13) at T = 0 with a mean-field approximation. We assume repulsive interaction (g > 0) and minimize the energy expectation value Ψ 0 |Ĥ|Ψ 0 with the variational wave function given by This wave function has two independent variational parameters E ↑ and E ↓ , corresponding to Fermi energies for spin-up and spin-down electrons, respectively. W (E σ ) denotes the region in the momentum space where the energy E k is below the variational parameter E σ : W (E σ ) = {k|E k ≤ E σ }. We note that for g < 0 the system becomes unstable against pairing and hence the variational wave function Eq. (18) is inapplicable.
The variational wave function |Ψ 0 gives the exact ground state at g = 0 by choosing the two variational parameters E ↑ = E ↓ = µ. For g = 0, the energy expectation value becomes where the electron density for spin σ at an energy E σ is given by We introduce a lower bound in the energy integral, i.e., the UV cutoff Λ(> 0), which corresponds to the inverse of the microscopic lattice scale. Since D(E) > 0, the electron density n(E σ ) is a monotonic function of E σ . The one-to-one correspondence allows an inverse function of n(E σ ); we define n σ ≡ n(E σ ) − n(0) to write E σ as a function of n σ : Now we can write the energy expectation value Ψ 0 |Ĥ|Ψ 0 as a function of n σ : where we introduce It is convenient to express the energy expectation value with the dimensionless quantities defined bȳ Then, we obtain where the dimensionless function is to be minimized by varyingn ↑ andn ↓ . The function Φ(n σ ) is given bȳ The electron densitiesn ↑ andn ↓ are order parameters in the mean-field analysis. Instead ofn ↑ andn ↓ , we usē where they corresponds to the order parameters for charge and magnetization, respectively. The values of n C andn M are obtained by minimizing the function E(n ↑ ,n ↓ ) with the chemical potentialμ and the coupling constantḡ given. The numerical result for = 1/4 and D r = 1/ √ 2 is shown in Fig. 3. We find discontinuities inn C andn M at the sameμ andḡ, which characterizes a first-order phase transition and defines a critical value of the coupling constantḡ c (μ). Finite magnetizationn M aboveḡ c characterizes a ferromagnetic state with the spin-rotational symmetry broken. The phase boundary in the numerical result well obeysḡ c (μ) ∝ D −1 (μ), as expected from Stoner criterion for ferromagnetism.

IV. ENERGY-SHELL RG ANALYSIS
The mean-field analysis in the previous section leads to a first-order transition to ferromagnetism with discontinuous changes of the charge density and magnetization. The ferromagnetic region shrinks and the discontinuity at the transition weakens as the interaction decreases. Nonetheless, in the mean-field theory, this transition occurs with infinitesimal repulsive interaction atμ = 0 because of the divergent DOS at VHS. However, this is an artifact of the mean-field analysis that neglects longwavelength fluctuations, which becomes increasingly important as the first-order transition becomes weaker. In this section, we perform an energy-shell RG analysis to study the role of these fluctuations near high-order VHS (|μ| 1) with weak repulsion (ḡ 1).

A. Formalism at zero temperature
Here, we adopt the Wilsonian approach to the RG equations for the action Eq. (14). For clarity, we consider first the action at T = 0, where we will find fixed points. Then, the Matsubara frequency becomes continuous ω n → ω, and the action is written with frequency ω and momentum k as We introduce the shorthand notation k = (k, ω). We impose a UV energy cutoff Λ on this action to remove unphysical UV divergences that appear in electron density of the ground state, etc. We note that the UV cutoff here is imposed on energy, but not on momentum directly. The region in k-space with |E k | ≤ Λ still extends to infinity. Importantly, this UV cutoff does not affect universal scaling properties of IR fixed points in the analysis of high-order VHS, as we shall show. The The thick line is the Fermi surface at the Van Hove energy, which is scale-invariant. The colored region has the energy inside the cutoff Λ, where the red (blue) area corresponds to E > 0 (E < 0). At every RG step of the energy-shell RG scheme, high-energy modes within the energy shell shown in darker colors are integrated out. In the field theory approach, all states below the cutoff Λ are integrated over at once.
UV cutoff merely appears in the prefactors of IR scaling functions.
We use two different energy cutoff schemes in this paper: an energy shell with a hard cutoff and a soft energy cutoff. The former scheme allows the Wilsonian RG approach, which offers a rather simple analysis and understanding. The latter requires a field theoretical analysis, which is apparently complicated, but high-order perturbative corrections become more tractable.
This section focuses on the energy-shell RG scheme, which imposes a constraint on momentum integrals. By converting the momentum integral to an energy integral with the help of the DOS, we write the momentum integral with the cutoff Λ as for an arbitrary function F. We denote the action with the energy cutoff Λ as S Λ , obtained by replacing the momentum integral k by Λ k . The UV energy cutoff designates an unbounded region in k-space, reflecting the extended Fermi surface with scale invariance (Fig. 4). Note that frequency integrals still range from −∞ to +∞. In a high-order VHS, divergences of momentum integrals arise from a singularity at k = 0 but not k → ∞. We will show that this simplifies the energy-shell RG analysis, which includes only the UV energy cutoff Λ. This is in contrast to a conventional VHS with a logarithmic divergence of the DOS [19,20]; it additionally requires a UV momentum cutoff. A further discussion can be found in Sec. VIII.
We now sketch how an RG transformation works with the energy-shell RG scheme. To access the IR behavior, we progressively eliminate UV modes and focus more on remaining modes. In the energy-shell RG scheme, we first split the energy range into two parts; one corresponds to lower energies E k ∈ [−Λ/b, Λ/b] and the other to higher energies E k ∈ [−Λ, −Λ/b) ∨ (Λ/b, Λ] (b > 1). Accordingly, the fermion field ψ is decomposed as where ψ < σ represents the low-energy modes and ψ > σ the high-energy modes. We write a momentum integral in the same way: Due to this division, the action is decomposed into the three parts as The first term S < [ψ < ] consists only of the low-energy modes ψ > and the second term S > [ψ > ] of the highenergy modes ψ > . The last term S <> [ψ < , ψ > ] describes the coupling of the low-and high-energy modes, which arises when the interaction is finite (g = 0). To obtain the effective action without the high-energy modes, we need to integrate out the high-energy modes: Now the high-energy modes are eliminated and the new action has the smaller cutoff Λ/b. One may be tempted to compare S Λ [ψ] and S Λ/b [ψ < ] to look into low-energy properties. However, it is like "comparing apples to oranges" [67] as the two actions are defined in different domains. For a fair comparison, we should make a change of variables (k, ω, and ψ) to restore the cutoff Λ. This procedure, called rescaling, completes the RG step. It results in the change of parameters in the model, which is described by RG equations.
The RG equations describe the flow of the parameters under a scale transformation. When the parameters do not change under a scale transformation, the system reaches an RG fixed point and exhibits scale-invariant properties. Away from a fixed point, the parameters flow. If the flow converges to a fixed point in its vicinity, then the fixed point is called a stable fixed point. If the parameters flow away from a fixed point, then it is an unstable fixed point. The RG equations also tell us how various susceptibilities and correlation lengths diverge as the critical point is approached, and the scaling properties of correlation functions at the critical point.

B. Tree-level analysis
The mixing term S <> can be calculated by expanding the logarithm in powers of the coupling constant g. We first consider the zeroth-order contribution in g. Since the remaining terms are described by tree diagrams without loops, the approximation is referred to as the treelevel analysis.
At tree-level, the effective action with the cutoff Λ/b becomes To compare with S Λ [ψ], we need to change the variables to put the cutoff Λ/b back to Λ. Now we change the variables so that the energy satisfies the relation For the energy dispersion given by Eq. (5), this immediately leads to rescaling of the momentum while the coefficients do not change: To retain the form of the action, we also need to rescale the field ψ, frequency ω, chemical potential µ, and coupling constant g to be When we look at the parameters of the model, the chemical potential µ and the coupling constant g change after an RG step, whereas the coefficients of the energy dispersion A ± do not. The flow of an parameter under an infinitesimal scale transformation (b → 1) is described by a differential equation, namely the RG equation. For µ and g, the RG equations are obtained from Eqs. (37c) and (37d): with l = ln b.
In the present case, we find the noninteracting fixed point at µ = g = 0 in Eq. (38), where the partition function takes a functional form of the Gaussian integral. If the parameters are away from the fixed point, they grow as l increases i.e., in low energies, and flow away from the fixed point. Therefore, the fixed point at µ = g = 0 is unstable and both µ and g are relevant perturbations to the unstable fixed point.
So far we have only considered the contact interaction. However, electron-electron interactions can take a more complicated form. Other types of interactions will be generated under RG even if not present initially, and thus their effects should be considered as well. In general, a finite-range interaction can be expanded in powers of spatial derivatives, with contact interaction being the lowest order term. The next leading term g − (ψ ↑ ∂ r−ψ↓ )(ψ ↓ ∂ r− ψ ↑ ) contains two spatial derivatives, and has a different scaling relation: g − = b −2/n− g − , which has a much smaller exponent than for the contact interaction. As an example, for the energy dispersion (3) in two dimensions, we have = 1/4 and n − = 4, so that g − is irrelevant. It is therefore legitimate to retain only the contact interaction in RG analysis.

C. One-loop analysis
In the presence of interaction, elimination of the highenergy modes gives rise to corrections in the effective action through the mixing of low-and high-energy modes in S <> [ψ < , ψ > ]. When depicted diagrammatically, S <> [ψ < , ψ > ] involves diagrams with loops, corresponding to integrations of the high-energy modes. We here consider perturbative corrections to one-loop order.
The effective action Eq. (35) can be calculated perturbatively with respect to the coupling constant g when it is small. We also treat the chemical potential µ as a perturbation as we are interested in critical phenomena where there is no characteristic scale in the system. Including perturbative corrections, we write down the action in the form where δg is a correction to the coupling constant and (· · · ) consists of interactions with derivatives that may be generated after integrating out the high-energy modes.
As we have discussed above, finite-range interactions are irrelevant, so that we can safely neglect them.
Perturbative corrections to the lowest order, namely to one-loop order, are diagrammatically depicted in Fig. 5(a) and (b), corresponding to Σ and δg, respectively. We find that the one-loop corrections to the selfenergy Σ and the coupling constant δg can be written as We emphasize that the all loop corrections should be evaluated at zero external frequency and momentum. The one-loop corrections are obtained to O(l) as where D(Λ) is the DOS at the cutoff energy and the dimensionless constants c H and c pp are We can see that the particle-hole contribution vanishes identically after the frequency integration, i.e., at T = 0 there is no particle-hole screening coming from states near the cutoff energy Λ. On the other hand, the particleparticle loop has a finite contribution. The Hartree contribution Σ H can be finite only when the DOS is asymmetric on the electron and hole side (D + = D − ), leading to a finite c H at most of order . There is no frequency or momentum dependence in the self-energy to one-loop order, so that the self-energy only renormalizes the chemical potential µ. The field renormalization or renormalization of the energy dispersion does not appear at one-loop order. They appear at two-loop order from the diagram shown in Fig. 5(d), which will be examined with the field theory approach in Sec. VI. There are two fixed points: the fixed point with g * 1 = 0 corresponds to the noninteracting fixed point and the other withḡ * 2 = /cpp to the nontrivial interacting fixed point. The interacting fixed point is stable along a line that connects the two fixed points, whereas the noninteracting fixed point is unstable on the plane.ḡ * 2 is a positive number of order , i.e., the interacting fixed point has weak repulsive interaction with its strength controlled by the DOS singularity exponent . We assume D+ > D− for the RG flow, which makes the DOS larger on the electron side. In such a case, the interacting fixed point is shifted from zero to beμ = O( 2 ) < 0. Note that the coupling constantḡ monotonically grows when the interacting is attractive (ḡ < 0).
With the one-loop corrections obtained, the new parameters µ and g after rescaling are which lead to the RG equations for the chemical potential µ and the coupling constant g. It is convenient to define the dimensionless chemical potentialμ and coupling constantḡ as Then, we obtain RG equations forμ andḡ as Since we are interested in the low-energy behavior, we consider the RG flow by increasing l. The RG flow is shown in Fig. 6. From the RG equations (45) for the coupling constantḡ and the chemical potentialμ, we find two fixed pointsḡ * g * 1 =μ * 1 = 0 corresponds to the noninteracting fixed point. The new fixed point atḡ * 2 > 0,μ * 2 < 0 is the nontrivial interacting fixed point with finite repulsive interaction, whose strength is of order . The smallness of the coupling constant allows a controlled analysis by the DOS singularity exponent about the interacting fixed point.
We can find the similarity to the φ 4 theory in the structure of the RG equation (45): the coefficient of the quadratic term rφ 2 corresponds to the chemical potentialμ and the quartic interaction term φ 4 to the coupling constantḡ. From this viewpoint, our theory can be regarded as the fermionic analog of the φ 4 theory. Like the Wilson-Fisher fixed point, our perturbative RG analysis is analytically controlled thanks to the smallness of the coupling constant on the order of at the interacting fixed point. While the φ 4 theory in three dimensions corresponds to = 1 in Wilson-Fisher RG, in our theory for high-order VHS in two dimensions takes the value of 1/4, given by the DOS exponent.
In the φ 4 theory, the RG flow of r describes the phase transition between ordered and disordered states: the RG flow to r 0 corresponds to the disordered state and r 0 to the ordered state, where the field φ has a finite expectation value associated with spontaneous symmetry breaking. The parameter r is analogous to the chemical potential µ in the present fermionic model, where µ 0 yields the electron Fermi surface and µ 0 the hole Fermi surface. The sign change of µ thus describes a topological transition between electron and hole Fermi liquids, which involves a change of Fermi surface topology without symmetry breaking.
Note that at the interacting fixed pointμ * 2 is nonzero when there is a finite contribution from the Hartree term c H = 0 due to the asymmetry of DOS at E > 0 and E < 0: D(±|E|) = D ± |E| − with D − = D + . This means that in the presence of repulsive interaction, the chemical potential at which scale-invariant Fermi surface appears is shifted from the noninteracting case, similar to the deviation of r at Wilson-Fisher fixed point from the mean-field value. For small , it follows from the expressions for D ± thatμ * 2 is at most of order 2 .

D. Relevant perturbations
We have identified the two fixed points: the noninteracting and interacting fixed points. With the chemical potential tuned at the fixed points, the noninteracting fixed point atḡ * 1 is an unstable fixed point and the interacting fixed point atḡ * 2 is a stable fixed point. The chemical potential is a relevant perturbation around both fixed points. We have included the chemical potential even in the analysis of the simplest case above as it can be generated by interaction due to the absence of particle-hole symmetry in the single-particle DOS.
In addition to the chemical potential, we consider other relevant perturbations to the fixed points, including the magnetic field h and the s-wave pairing field ∆. Those relevant perturbations add the following terms to the ac-tion at criticality: Finite temperature is also a relevant perturbation. Its effect is taken account of via Matsubara frequencies; see Appendix A. We further consider other relevant perturbations. For an energy dispersion x − k 4 y , the fermion bilinear terms with derivatives ∂ x , ∂ y , ∂ 2 y , ∂ 3 y , ∂ x ∂ y are also relevant perturbations. Perturbations to the system are subject to symmetry constraints: Particle conservation forbids the pairing term, spin-rotational symmetry nonzero h, and reflection symmetry odd-derivative terms in x or y. With all three symmetries present, only two terms µψψ andψ∂ 2 y ψ are allowed as perturbations to the system with E k = k 2 x −k 4 y ; see Fig. 2. This means that we need to tune two parameters to reach the critical metallic state governed by the interacting RG fixed point shown earlier.
In our RG analysis so far, the starting point is the single-particle dispersion at the high-order VHS, where the termψ∂ 2 y ψ is absent. To one-loop order, this term is not generated from the interaction since the self-energy Σ is independent of momentum. However, it may be generated at higher-loop order. As we shall show later, this means that in the presence of interacting, the critical state is reached when theψ∂ 2 y ψ term is present in the single-particle dispersion and its coefficient is tuned to a particular value.
The pairing field can be introduced by proximity to an external superconductor, or it can be regarded as a test field for studying s-wave pairing susceptibility. Likewise, the magnetic field h can be externally introduced or regarded as a test field for the spin susceptibility. In this viewpoint, the chemical potential is conjugate to the particle number, and hence it is related to the charge compressibility.
Corrections to the perturbations h and ∆ are calculated similarly as those for µ and g at T = 0. To consider a correction to the pairing field ∆, we include the particle-particle loop diagram, where the one-loop diagram is shown in Fig. 5(c). We include the corrections to write the magnetic field h + δh and the pairing field ∆ + δ∆.
Integrating out the high-energy modes is followed by rescaling. The parameters of the model should be rescaled at tree level as h = bh and ∆ = b∆. Those parameters are relevant and thus their values increase as we proceed with RG steps. When the perturbative corrections are included, the new parameters after an RG step are To one-loop order, the correction terms are expressed as The one-loop correction Π pp is obtained in Eq. (41b). Then, the parameters change as With the dimensionless quantities we reach the RG equations We confirm that the perturbations h and ∆ are relevant around the two fixed points, given in Eqs. (46) and (47). Finite temperature is also a relevant perturbation, which scales in the same manner as energy and frequency. All low-energy fixed points are found at T = 0, and thus we focus on zero temperature in the main part. The oneloop RG equations at finite temperature are presented in Appendix A. The physical consequences, i.e., scaling properties of thermodynamic quantities, are discussed in the next section.

E. Structure of higher-order RG
So far, we have made the energy-shell RG analysis to one-loop order. We now illustrate how it works in the case with higher-order corrections. Again, for clarity we consider here the minimal case at T = 0 without symmetry-breaking fields. Inclusion of other relevant contributions such as T , h, and ∆ is straightforward.
Higher-order perturbative corrections give rise to the frequency and momentum dependence in the self-energy Σ in Eq. (39), while the one-loop corrections are independent of frequency or momentum and depend only on the DOS as we have seen. We expand the self-energy with respect to the frequency and momentum to find corrections to the field, energy dispersion, and chemical potential.
As we shall show later in Sec. VI, the momentum dependence may give corrections to the energy dispersion. In that case, one has to be wary of the generation of relevant corrections in the single-particle energy dispersion even when they are initially absent. We represent such a term as λk n , where the coefficient λ transforms under Eq. (7) as λ = b a λ (a > 0). For instance, for the case of E k = k 2 x − k 4 y , this term corresponds to λk 2 y (a = 1/2), which is the only relevant perturbation to the energy dispersion.
In order to keep track of such relevant term(s), we include λk n in the energy dispersion: At least one such relevant perturbation term exists for a high-order VHS, and if present, turns a high-order saddle point in the noninteracting single-particle dispersion into ordinary one. For the generalized dispersion Eq. (11), there is only one relevant perturbation to find For other types of dispersion, it is possible to have multiple relevant terms. An extension to a case with multiple relevant terms is straightforward.
Then, the expansion of the self-energy is generally given by where irrelevant high-order terms are safely neglected. After integrating out the high-energy modes within the energy shell, we obtain the effective action The next step in the energy-shell RG analysis is to rescale the momentum and restore the energy cutoff Λ/b to Λ; see Eqs. (36a) and (36b). However, the effective action S Λ/b still evidently has a different form from S Λ . To recover the form of the action, we rescale the other quantities as follows: Here we introduce the scaling exponents γ A± , γ µ , and γ ψ . Note that there is an ambiguity in defining ω and ψ as the factor (1 − Σ ω ) can be imposed on either ω or ψ . We choose to scale ω linearly in b and hence the factor (1 − Σ ω ) contributes to the field renormalization. For γ A± = 0, if we continue to rescale momentum according to Eq. (36b) and the coefficients A ± according to Eq. (57b), the cutoff energy Λ/b is not mapped to Λ. To remedy this issue, we rescale momentum as so that E k = bE k is satisfied. In this way, the coefficients A ± do not change under rescaling. At a fixed point, the parameters in the action are determined to satisfy scale invariance; i.e., they do not vary under rescaling [Eqs. (57c)-(57e)]. To reach a fixed point, initial values of the relevant perturbations µ and λ should be tuned so that they cease to flow when the coupling constant g reaches the fixed point value.
Rescaling of the magnetic field h and the pairing field ∆ can be considered similarly. Including the field renormalization, we obtain where we define the exponents γ h and γ ∆ . We shall show later that the Ward identity requires γ µ = γ h = 1.

V. ANALYSIS
In this section, we combine the results obtained from the mean-field and RG analyses to present a phase diagram of interacting electrons near a high-order VHS. We then describe scaling properties for thermodynamic quantities and correlation functions near the supermetal critical point. In addition, we discuss the Ward identity, which results from charge conservation and gives relations among scaling exponents of electronic specific heat, magnetic susceptibility, and charge compressibility.

A. Phase diagram
As we have discussed in Sec. II, realization of a highorder VHS requires tuning of the energy dispersion in addition to the chemical potential. There is at least one parameter for a relevant perturbation in the energy dispersion to be tuned; see Sec. IV D. Therefore, to present a global phase diagram near high-order VHS, we need three axes for the coupling constant g, the chemical potential µ, and a tuning parameter λ of the energy dispersion. All those three are relevant perturbations at noninteracting fixed point; we use dimensionless parameters defined byḡ = g/Λ ,μ = µ/Λ, andλ = λ/Λ a .  We now incorporate the results from the mean-field analysis (Fig. 3) and the RG analysis (Fig. 6). The meanfield analysis is expected to be qualitatively correct for relatively largeḡ, while the RG analysis is valid for small µ andḡ. Based on these considerations, we propose a global phase diagram of interacting electrons near highorder VHS shown in Fig. 7.
For largeḡ, mean-field calculation reveals that an itinerant ferromagnetic metal exists over a wide range of chemical potential, and the ferromagnetic transition is first order. Due to the finite correlation length, we expect these results to be qualitatively correct and continue to hold in the presence of a small λ.
On the other hand, for smallḡ, there is no spontaneous symmetry breaking or long-range order. When the DOS is not divergent, the system is either a electron or a hole Fermi liquid depending on the sign ofμ. These two Fermi liquid states are indistinguishable by symmetry but differ in the Fermi surface topology. A transition between electron and hole Fermi liquids, i.e., a Lifshitz transition, occurs as the chemical potential crosses the VHS. This transition occurs on a surface in the three-dimensional phase diagram.
Our RG analysis reveals that by tuning both µ and λ, a multicritical line on the Lifshitz transition surface can be reached, where the system displays various di-vergent susceptibilities and scale-invariant Fermi surface. We coin a term, supermetal, to describe such an unusual metallic state. At the end of this multicritical line g = µ = λ = 0, the noninteracting supermetal exhibits divergent charge, spin and pairing susceptibilities determined by the power-law divergent DOS. In contrast, for g = 0, the interacting supermetal displays universal critical properties governed by the nontrivial interacting fixed point, located atḡ * = /c pp ,μ * = 0,λ * = 0 to first order in . As we shall show in next subsection, at this fixed point, while the charge compressibility and spin susceptibility diverge, the s-wave pairing susceptibility remains finite. We shall also show later by a two-loop RG calculation that the electron Green's function has the scaling form G(ω) ∝ 1/|ω| 1−η with η > 0, thus establishing the non-Fermi liquid nature of an interacting supermetal.
The interacting fixed point is stable along the multicritical line and unstable in two other directions. One of the unstable direction (roughly speakingλ) lies within the Lifshitz transition surface, while the other direction (roughly speakingμ) drives the system into electron or hole Fermi liquid. Note that a finiteλ converts a highorder VHS to a conventional VHS (λ > 0) or splits it to two conventional VHS points (λ < 0). For the latter case for Eq. (3), over a finite range of the chemical potential (Fig. 2) there is an extra small pocket around k = 0 in addition to large Fermi surfaces.
Since the relevant perturbations µ and λ introduce an intrinsic momentum scale to the system, the resulting Fermi liquids may be unstable to superconductivity at very low temperature via the Kohn-Luttinger mechanism associated with non-analyticity of susceptibility at momentum 2k F [68]. This scenario is neglected in the phase diagram (Fig. 7). In contrast, being a quantum critical state of metal at T = 0, the interacting supermetal is immune from the Kohn-Luttinger mechanism since its Fermi surface is scale-invariant without any intrinsic scale.
Finally, we conjecture how the ferromagnetic transition at largeḡ andμ and the Lifshitz transition at small g andμ meet together. A plausible scenario is that the multicritical line of supermetal meets the first-order ferromagnetic transition line at a tricritical point between electron Fermi liquid, hole Fermi liquid and ferromagnetic metal. The nature of this tricritical point is an interesting open question.
B. Scaling analysis

Generic case
Scale invariance at the fixed points enables us to extract various scaling relations. Since the partition function Z is invariant under the scale transformation, the free energy density F , defined in Eq. (17), reflects the scaling of the factor T /V : where the volume V scales according to Eq. (58) and temperature scales the same manner as energy and frequency. For convenience, we rewrite the exponent as By explicitly showing the parameters of F , we obtain the scaling relation of the free energy density Here, the scaling exponents γ A± , γ µ , γ h , and γ ∆ correspond to the values at a fixed pointḡ * . We later see γ µ = γ h = 1, but we keep them in the following scaling analysis. The coupling constant g itself does not appear in the scaling relation of the free energy density F , but the effect is imprinted on γ ∆ and˜ as the fixed point properties. We shall see that γ A± are at most of order 2 at the interacting fixed point and thus˜ is also a small positive quantity. We then consider the critical exponents of the charge compressibility κ, magnetic susceptibility χ, heat capacity per unit volume C V , and s-wave pairing susceptibility χ BCS . From Eq. (62), we find We also examine the pair correlation function with O = DψDψOe −S /Z. From the comparison between χ BCS and C(r, τ ), we obtain the scaling form where ν is an arbitrary energy scale andĉ is a scaling function.

One-loop results
To one-loop order, we find from the RG equations (45a) and (53) the exponents at the fixed points with˜ = . Most exponents in Eqs. (63)- (65) are the same at the noninteracting and interacting fixed points, which is identical to that of the DOS in the noninteracting state. The difference is found when the pairing field ∆ is involved. The exponent for the pairing field γ ∆ renders different exponents for the pairing susceptibility χ BCS : The s-wave pairing susceptibility remains finite at the interacting fixed point whereas it diverges at the noninteracting fixed point. We also find a difference in the pair correlation function It shows a faster decay at the interacting fixed point, reflecting the suppressed pairing susceptibility.

C. Ward identity
In Sec. V B 1, we mentioned the relations γ µ = 1 and γ h = 1. They result from the conservation laws for charge and spin. The Ward identity (more generally the Ward-Takahashi identity) describes a conservation law [69,70]. The identity is regarded as the quantum analog to Noether's theorem. We present how the Ward identity works in our present analysis. The identity should hold even after an RG analysis, and thus it can be used to check the validity of an RG scheme, or specifically a choice of a cutoff. It also gives relations between the exponents for thermodynamic quantities Eqs. (63)- (65). Now we investigate the structure of the self-energy and vertex corrections. To be concrete, we look into the expansion of the self-energy Eq. (55) to find a relation between Σ ω and δµ. The Ward identity concludes at T = 0. The identity is based on charge conservation or the U(1) gauge invariance; the action and correlation functions are invariant under the transformations ψ → e iθ(r,τ ) ψ andψ →ψe −iθ(r,τ ) with a smooth scalar function θ(r, τ ). In the present model, charge conservation holds for each spin separately, thus leading to Then, Eqs. (57d) and (59a) yield The result of the energy-shell RG analysis to one-loop order in Sec. IV satisfies the Ward identity, which means that the conservation laws are correctly taken account of. We notice that a frequency shell instead of the energy shell violates the Ward identity.
Furthermore, Eq. (76) makes some ratios among scaling relations Eqs. (63)-(65) constant as functions of temperature T . One is the Wilson ratio R W between the electronic specific heat C V and the magnetic susceptibility χ and the other is the ratio R C between the charge compressibility κ and C V : In the following, we sketch the derivation of the Ward identity from the diagrammatic point of view. A detailed derivation is given in Appendix B. In the present analysis, the Ward identity relates the frequency derivative of the self-energy and the vertex function corresponding to the coupling term α σψσ ϕψ σ in the action. α σ is the spin-dependent coupling constant and the ϕ is a bosonic field. We write the vertex function as Γ (2,ασ) σ (ω + ω , ω), where we focus only on the frequency dependence. The vertex function modifies the coupling term to be α σ Γ (2,ασ) σ (ω + ω , ω)ψ σ (ω + ω )ϕ(ω )ψ σ (ω). The derivation of the Ward identity makes use of the equality G −1 0 (k, ω + ω ) − G −1 0 (k, ω) = iω , or equivalently (78) This equation is diagrammatically shown in Fig. 8(a). It relates the noninteracting vertex function Γ (2,ασ) σ (ω + ω , ω) = 1 and the noninteracting Green's function G 0 . Now we add corrections to the self-energy, depicted in Fig. 8(b) as shaded blobs. The dressed vertex function is obtained from the dressed self-energy by attaching the external scalar field ϕ to every internal fermion line. As a result, we find the Ward-Takahashi identity The full Green's function G(k, ω) is given by with the full self-energy Σ. Taking the zero frequency limit ω → 0, we obtain the Ward identity The vertex function Γ (2,ασ) σ (ω, ω) gives the quantum correction to the coupling α σ to be α σ + δα σ = α σ Γ (2,ασ) σ (ω, ω). α σ represents the chemical potential with α ↑ = α ↓ = −µ and the magnetic field with α ↑ = −α ↓ = h. Since the right-hand side of Eq. (81) is independent of spin σ, we confirm Eqs. (74) and (75).

VI. FIELD THEORY APPROACH
This section focuses on the RG analysis from the field theory approach. To begin with, we briefly argue the two RG schemes: the energy-shell RG analysis and the field theory approach. We then confirm that the two methods give the same result at one-loop order. We also perform a two-loop analysis of the self-energy (two-point function) from the field theory approach to show the anomalous dimension and the correction to the energy dispersion.

A. RG schemes
An objective of RG analyses is to track the flow of parameters in a theory under a scale transformation. Here, we illustrate two different RG schemes: the Wilsonian approach, including the preceding energy-shell RG analysis, and the field theory approach. The common feature is to divide the integration manifold (frequency and momentum in the present case) into two parts and integrate out modes belonging to one of them. The two schemes differ in intervals of integrations. The first scheme involves an integration within a hard shell. In the energyshell RG analysis, fluctuations inside the thin energy shell E ∈ [−Λ, −Λ/b) ∨ (Λ/b, Λ] are eliminated. This mode elimination followed by rescaling enables us to keep track of the change of parameters under a scale transformation. On the other hand, in the field theory approach, we integrate out all low-energy fluctuations below the cutoff Λ. Then, we deduce the RG flow of parameters by comparing results at different cutoffs Λ and Λ . The two schemes have advantages in different aspects. In the Wilsonian approach, the frequency-momentum space is progressively integrated over, so the interpretation of the RG procedure is rather simple. The inclusion of low-energy modes results in a theory at low energies with different parameters. In spite of its simple interpretation, higher-loop calculations are not easy with the Wilsonian approach. In a one-loop calculation, we have only one shell to be concerned about. However, higher-loop diagrams consist of many internal lines (virtual states), so that we have to take care of shells for each of them. On the other hand, the field theory approach does not require such error-prone steps as it deals with all modes below the cutoff at once. This makes higher-loop calculations more tractable. Although not as intuitive as the Wilsonian approach, the field theory approach leads to the same results about critical phenomena. More descriptions about the comparison between the two schemes can be found in e.g. Ref. [67]. A brief review of the field theory approach is given in Appendix C.

B. Soft cutoff
In the field theory approach, we calculate the connected N -point correlation function G (N ) or the oneparticle irreducible N -point function Γ (N ) . If we face a UV divergence in calculating them, we need to cure the divergence to obtain physically meaningful results. There are several ways to do so; we here choose to employ the UV energy cutoff Λ to make a comparison to the preceding energy-shell RG analysis.
The functions G (N ) and Γ (N ) can be obtained perturbatively with the noninteracting Green's function G 0 . We introduce the UV energy cutoff by suppressing the high-energy contributions in G 0 . We define the noninteracting Green's function with the energy cutoff G 0Λ (k, ω n ) as with the UV energy cutoff factor Note that the cutoff factor smoothly varies from 0 to 1 and thus works as a soft energy cutoff. This is in contrast to the energy-shell RG analysis, where the interval of an energy integration is cut off abruptly at Λ/b and Λ. We can interpret the modified Green's function as a Green's function with an energy-dependent quasiparticle weight K Λ (E). The weight fades away in the highenergy limit E → ±∞ to eliminate UV divergences, while K Λ (E) → 1 for energies much lower than the cutoff Λ. One may be tempted to see the modified Green's function in a different way. For example, it can be rewritten as It may be viewed as a variation of the Pauli-Villars regularization, where the additional term cures a UV divergence but vanishes in the limit Λ → ∞. However, we cannot think of it as a propagator with a large mass term since we cannot add a mass term for the electronic energy dispersion which is continuous and unbounded.
It should be noted that the cutoff factor K Λ (E) does not depend on frequency. It potentially causes a violation of the Ward identity, which would result in wrong conclusions. For example, if one chooses a cutoff factor of the form Λ 2 /(Λ 2 + E 2 + ω 2 n ), it invalidates the Ward identity. The absence of the frequency in the cutoff factor ensures the Ward identity.

Structure of the RG analysis
To derive RG equations and see scaling properties, we calculate the one-particle irreducible N -point function Γ R . The energy scale at the reference point is referred to as the renormalization scale. The procedure of fixing the model to the reference is equivalent to setting the initial parameters in the Wilsonian approach.
We first analyze the case with T = h = ∆ = 0. We impose the renormalization conditions where the condition for Γ (4) should be considered at k 1 + k 2 = k 1 + k 3 = k 1 + k 4 = 0. The subscript 0 denotes quantities at the renormalization scale. The interaction dresses the two-point and four-point functions and they acquire cutoff-dependent corrections. We here use the energy dispersion Eq. (54), which includes a relevant perturbation to a high-order VHS, since such a term could be generated under the RG analysis at two-loop order or higher; see the discussion in Sec. IV E. Then, the two-point and four-point functions at the cutoff Λ can be expressed as where the corrections Z ψ , Z A± , Z λ , Z µ , and Z g are calculated perturbatively. The N -point functions at the renormalization scale and the cutoff Λ are related by We note the structure of the RG analysis is general, so that an analysis of other energy dispersions such as Eq. (11) is straightforward.
The last equation leads to the RG equations. Since the left-hand side does not depend on the cutoff Λ, we obtain the differential equation leading to the Callan-Symanzik equation [71][72][73]. We obtain the Callan-Symanzik for the one-particle irreducible N -point function The beta functions and γ ψ are defined by Since the renormalized values are given bȳ we can rewrite the beta functions as Those equations show that the field renormalization gives additional effects to the beta functions and hence the scaling properties.

Solutions
The Callan-Symanzik equation can be solved by the method of characteristics; see Appendix C. The beta functions describe the RG flows of the parameters: dλ dl = β λ (ḡ,λ). (95d) l = ln Λ 0 /Λ denotes the RG scale, measured relative to the renormalization scale Λ 0 . Those RG equations are to be compared with those obtained by the energy-shell RG analysis in Sec. IV. In general, they are coupled differential equations and zeros of the beta functions determine fixed points. We can write the beta functions β µ , β A± , and β λ around a fixed point withḡ * as γ µ (ḡ * ), γ A± (ḡ * ), and γ λ (ḡ * ) give the exponents in the scaling region. Recall that γ µ = 1 is required by the Ward identity, regardless ofḡ. From the beta functions around the fixed point, we observe the scaling properties Since the energy dispersion does not receive correction at one-loop order, we have γ A± (ḡ * ) = O( 2 ) and γ λ (ḡ * ) = a + O( 2 ). The shift of the chemical potential and the generation of the relevant perturbation to the energy dispersion are also seen from the beta functions. When β µ (ḡ, 0) = 0, the chemical potential is displaced from zero under the RG analysis, while it does not alter the scaling behavior ofμ. Similarly, a finite relevant perturbationλ is generated if β λ (ḡ, 0) = 0, even when it is initially absent.
The function γ ψ is ascribed to the anomalous dimension η when it is computed at a fixed point. To see this, we solve the Callan-Symanzik equation (91); see Appendix C for details. The solution of the two-point function is given by We now examine the behavior in the critical region as a function of ω, k + , and k − . We assume the two-point function is a function of A + k n+ + , A − k n− − , ω in the scaling region. Since those three quantities, Λ, and Γ (2) Λ have the dimension of energy, the two-point function can be written as whereΓ (2) Λ is a dimensionless scaling function. Here, we do not need to assume homogeneity forΓ  98), l is an arbitrary quantity; to inspect the scaling behavior in terms of ω, we set l = ln(ω 0 /ω) and k + = k − = 0. The momentum dependence is considered in the same manner with l = ln(k ±,0 /k ± ) n± and Eq. (97). We then find where γ ψ (ḡ(l)) = η is used. It confirms the scaling relation of the two-point correlation function Eq. (69) along with the relation G = [Γ (2) ] −1 .

D. One-loop calculations (h = ∆ = 0)
We calculate the two-point and four-point functions to obtain the beta functions and γ ψ (ḡ). This is accomplished by evaluating the perturbative corrections to the two-point and four-point functions (Figs. 9 and 10). As the corrections to the coupling constant g, there are three possible one-loop diagrams shown in Fig. 10. To determine the perturbative correction δg, all diagrams should be evaluated with zero momentum transfer q = 0, which is required by the renormalization condition Eq. (86). The three one-loop diagrams in Fig. 10(a) correspond to the BCS, density-density, and exchange channels (from left to right). Out of the three, the density-density channel does not contribute because of the Pauli exclusion principle for the contact interaction. This contribution is allowed when we assume the density-density interaction in finite rangeψ σψσ ψ σ ψ σ with arbitrary spins σ, σ or when there is an additional valley of orbital degree of freedom. (For reference, we note that the three channels are referred to as the BCS, ZS (zero sound), and ZS in Ref. [67]; or s-, t-, and u-channels with the Mandelstam variables.) To one-loop order, the two-point and four-point functions give corrections to the chemical potential and the coupling constant, but not to field or the energy dispersion as we have seen in the energy-shell RG analysis. One-loop diagrams can be represented by Σ H , Π pp , and Π ph like Eq. (40). Then, the two-point and four-point functions become which lead to Here we calculate the perturbative corrections with the soft cutoff K Λ . The actual calculations for the beta functions require the Λ-derivatives instead of the corrections themselves. We thus obtain the one-loop corrections as follows: As a result, we obtain the beta functions Eqs. (92a) and (92b) Note that the tree-level scaling terms appear from the definitions of the dimensionless parametersḡ 0 = g 0 D(Λ) andμ 0 = µ/Λ. The beta functions are to be compared with the result from the energy-shell RG analysis Eq. (45). To confirm, we first evaluate the coefficients c H ,c pp ,c ph at T = 0: The zeros of the beta function β(ḡ) give the two fixed pointsḡ * We now find the noninteracting and interacting fixed points from the field theory approach. Although the value ofḡ * 2 differs in the two schemes, resulting exponents for the thermodynamic quantities are not suffered from the difference as the exponents are not directly dependent on the coupling constantḡ at fixed point. We explicitly confirm this in the next subsection by calculating the beta functions for the magnetic field h and pairing field ∆.

E. RG equations for h and ∆
The beta functions for the magnetic field h and pairing field ∆ can be obtained from the corresponding vertex functions Γ (2,h) and Γ (2,∆) , respectively. Perturbative corrections to them are depicted in Figs. 9(b) and 9(c). We impose the renormalization conditions where the vertex functions with the cutoff Λ are expressed as To obtain the beta functions to one-loop order, it is sufficient to consider the Callan-Symanzik equations without corrections to the energy dispersion and the chemical potential: where the beta functions for the magnetic field and pairing field are defined by Using the relations the beta functions become They are related to the exponents γ h and γ ∆ when evaluated at a fixed point: We calculate the vertex functions for h and ∆ to oneloop order and find The vertex functions lead to the beta functions Now we confirm by taking T → 0 that the exponent for the pairing field ∆ is the same independent of the RG schemes. Particularly at the interacting fixed point, we obtain β ∆ (ḡ * 2 ) = (1 − )∆. This is consistent with the result from the energy-shell RG analysis. The coefficient c pp , which determines the value of the coupling constant at the interacting fixed point, does not appear to the exponent of the pairing field.

F. Two-loop calculations
So far we have calculated the perturbative corrections from the field theory approach to confirm that the two distinct RG schemes conclude the same physical results. An advantage of the field theory approach is considerable when we deal with higher-order corrections. In the following, we consider the two-loop corrections to the twopoint correlation function at T = 0 for the anomalous dimension and the correction to the energy dispersion.
The field renormalization is seen from the frequency dependence of the self-energy. The linear term Σ ω in Eq. (55) is given by We expand Σ with respect to the coupling constant g.
On the other hand, the zero-frequency part is related to corrections to the chemical potential and the energy dispersion: The corrections δµ, δA ± , and δλ are obtained as We have used the fact that the one-loop correction, i.e., the Hartree contribution, does not yield the frequency or momentum dependence. The renormalization condition Eq. (85) reads The field renormalization γ ψ is expressed from Eq. (92e) as and the beta functions for the chemical potential and the coefficients of the energy dispersions are obtained from Eqs. (94b)-(94d) as We now calculate the two-loop correction to the self-energy Σ (2) . For the case of the contact interaction, there is only one connected two-loop diagram, i.e., the sunrise diagram shown in Fig. 5(d) and 9(a) as the rightmost term.
The frequency and momentum dependent contribution appears from this diagram, calculated from We use the shorthand notations p = (p, ω p ) and p = dωp (2π) and the momentum-dependent part (125b) Here we denote the dimensionless quantities by adding bars; we defineω = ω/Λ,p + = p + /Λ 1/n+ , andp − = p − /Λ 1/n− . The momentum is scaled by Λ so that the energy becomes dimensionless: Ek = E k /Λ. ± p = p Θ(±E p ) stands for the momentum integral within the positive (negative) energy domain. The constraints on the momentum integrals emerge after the frequency integrals. They can be evaluated by identifying the position of poles on the complex plane, leading to the restricted regions of the momentum integrals.
We expect finite results for the two-loop results Eqs. (125a) and (125b) at a saddle point of an energy dispersion because of the constraints on the momentum integrals ± pq ∓ l . The two-loop contributions vanish at a band edge since there is no sign change in the energy dispersion. Now we scrutinize the frequency-dependent part Σ (2) ω , which is responsible to the field renormalization and hence the anomalous dimension. As we have discussed, the contribution vanishes at a band edge and thus an anomalous dimension does not arise. It can be finite only at an energy saddle point. In addition, it is worth pointing out that the integrand of Eq. (125a) is guaranteed to be positive. Therefore, if there exists a finite volume that satisfies the constraint of the momentum integrals, we find a finite result: C (2) > 0. The constraints on the momentum integrals can be rephrased as follows: There exists a momentum l = p + q such that sgn(E l ) = − sgn(E p ) = − sgn(E q ). Such a momentum l in general exists near a saddle point because the energy dispersion near a saddle point comprises two or more filled Fermi seas and the area is not convex. We do not further evaluate the expression of the two-loop correction as its value depends on the explicit form of the energy dispersion.
Equation (125a) defines a numerical factor C (2) , which is independent of the cutoff Λ. From Eq. (122), we find the field renormalization This quantity gives the anomalous dimension when evaluated at a fixed point. It can be finite at the interacting fixed point to become A finite anomalous dimension concludes a non-Fermi liquid behavior at the interacting fixed point. This happens at a saddle point of an energy dispersion with a powerlaw DOS singularity. The uniform component of Eq. (125b) adds a correction to the beta function for the chemical potential Eq. (105b), but it does not change the structure of the RG flow for smallḡ. Here, we focus on the momentum dependence, which is absent to one-loop order. Similarly to C (2) , it becomes finite only at a saddle point of an energy dispersion but not at a band edge. The momentum dependence of C (2) k leads to the beta functions From Eqs. (96b) and (96c), we can identify the scaling exponents γ A± and γ λ . The former affect the exponents of susceptibilities via Eqs. (58) and (61). We can see that a finite C (2) (> 0) and hence an anomalous dimension has a negative contribution to γ A± . When γ A+ and γ A− are negative, we have˜ > , leading to stronger divergences with respect to T , µ, and h; see Eqs. (63)- (65). When β λ (ḡ,λ = 0) is finite, the relevant perturbation to the energy dispersion λk n is generated under the RG analysis. It is observed if ∂C (2) k /∂k n | k=0 does not vanish when it is evaluated with λ = 0. Then, the beta function for λ has the form β λ = (a + c 1ḡ 2 )λ + c 2ḡ 2 + O(ḡ 3 ), where c 1 and c 2 are determined by Eq. (128b). This is analogous to the shift of the chemical potential when the Hartree term Σ H is finite, but it occurs at different order inḡ. Generation ofλ curves the scale-invariant line in the phase diagram [ Fig. 7(b)] in theλ direction at order g 2 , while a change in theμ direction can happen at order g. Lastly, we note that the discussion from Eq. (125) is general for any energy dispersion with a power-law divergent DOS, including Eq. (11) with a relevant perturbation λk 2 − .

VII. QUASIPARTICLE DECAY RATE
The preceding RG analyses focused on the real part of the self-energy or equivalently the two-point function. They give rise to the corrections to the action, which are captured through the RG equations. On the other hand, the imaginary part of the self-energy describes the damping of the quasiparticle, which is the focus of this section. It is generated by the interaction in the present model. Unlike the real part of the self-energy, the imaginary part can be calculated without a cutoff; we do not employ an RG method in this section, but integrate over the entire frequency and momentum space at once.
We calculate the quasiparticle decay rate Γ(k, ω), obtained from the retarded self-energy as The retarded self-energy Σ R (k, ω) is calculated from the self-energy Σ(k, ω n ), with the analytic continuation of the Matsubara frequency to the real frequency: iω n = ω + iδ (δ: infinitesimal positive quantity). In the presence of the contact interaction, a finite imaginary part of the self-energy Σ R emerges at two-loop order and higher. The one-loop correction, or the Hartree term Σ H , does not yield a finite imaginary component. Here we consider the two-loop diagram (the sunrise diagram) [Fig. 5(d)] to calculate the quasiparticle decay rate Γ. Like Eq. (124), it is given by but we do not need a cutoff for the imaginary part. The calculation of Σ (2) is standard and can be found in e.g. Ref. [74]; we also show the derivation in Appendix D and just present the result here. The quasiparticle decay rate to two-loop order is given by Σ (2) after the analytic continuation: This relation holds for an arbitrary energy dispersion E k . We extract the temperature dependence by introducing dimensionless quantities in terms of temperature T : we definep ± = p ± /T 1/n± , so that the energy dispersion satisfies Ep = E p /T . Here we are interested in the lowfrequency limit with ω T . By substituting k = 0 and ω = 0, we obtain the temperature dependence [4] The integral gives a finite constant without a cutoff.
In the Fermi liquid theory, when the temperature is much smaller than the Fermi energy F T , the decay rate is proportional to T 2 . This result relies on the existence of the Fermi surface with finite DOS. On the other hand, the decay rate Eq. (132) is distinct from the Fermi liquid results, reflecting the divergent DOS at µ = 0. The behavior is different also from the case for a conventional VHS with a logarithmic DOS, which shows a (marginal) Fermi liquid behavior [23][24][25][26]. We note that the result does not depend on whether the power-law divergent DOS is located at a saddle point or a band edge of the energy dispersion. This is in contrast to the anomalous dimension, which can only be found at a saddle point as we have discussed in Sec. VI F.

VIII. SUMMARY AND DISCUSSIONS
We now summarize our main results, compare supermetal with normal metal and other non-Fermi liquid systems, and discuss experimental signatures of supermetal.

A. Summary
We have analyzed electron interaction effects near a high-order VHS with a scale-invariant Fermi surface with a power-law divergent DOS. Scale invariance of the system allows an RG analysis to search for fixed points and a scaling analysis of thermodynamic quantities and correlation functions around the fixed points.
The one-loop RG analysis finds that electron interaction around high-order VHS offers a fermionic analog of the φ 4 theory. We have identified the two RG fixed points: the noninteracting and interacting fixed points. The latter is an analog of the Wilson-Fisher fixed point in the φ 4 theory. Like the φ 4 theory, the noninteracting fixed point is unstable and the interacting fixed point is stable in terms of the RG flow of the coupling constant.
We performed a controlled RG analysis up to twoloop order about the interacting fixed point owing to the smallness of the DOS singularity exponent . We reveal that the quantum critical metal at the interacting fixed point is a non-Fermi liquid that exhibits a finite anomalous dimension of electrons, and power-law divergent charge and spin susceptibilities. We term such a metallic state with various divergent susceptibilities but yet without a long-range order as a supermetal. In this regard, the noninteracting fixed point can be viewed as a noninteracting supermetal and the interacting fixed point as an interacting supermetal.
A supermetal appears at the topological transition between electron and hole Fermi liquids. An interacting supermetal is a multicritical state reached by tuning two parameters-chemical potential and detuning of energy dispersion from high-order saddle point. Combining the RG and mean-field analyses, we conjecture a global phase diagram where a supermetal is at the border between electron/hole Fermi liquids and on the verge of becoming ferromagnetic.

B. Comparison with normal metal and other non-Fermi liquids
It is worth drawing a comparison between a supermetal and a normal metal. Being at finite density, a normal metal is characterized by a Fermi surface with a characteristic momentum scale. The RG theory of metals with a closed Fermi surface, commonly referred to as Shankar's RG [67], requires a judicious RG procedure that only consider electrons within a small energy shell around the Fermi surface. Then, Fermi liquid appears as the RG fixed point in the limit that the energy range is taken to zero. It is characterized by an infinite number of marginal coupling constants, i.e., Landau forward scattering parameters in all angular momentum channels. Moreover, this Fermi liquid fixed point is only stable when BCS interactions in all angular momentum channels are repulsive [68].
These behaviors of a normal metal should be contrasted with the case of a supermetal. Our theory is formulated with a large UV energy cutoff on the order of bandwidth. The supermetal fixed point is characterized by a single coupling constant-the contact interaction, with all other interactions being irrelevant and without suffering from the Kohn-Luttinger instability to superconductivity.
We have shown that a high-order saddle point with repulsive interaction exhibits the non-Fermi liquid behavior. Non-Fermi liquids are realized also in e.g., onedimensional systems, other kinds of quantum critical metals, and doped Mott insulators. In a one-dimensional electronic system, there is no quasiparticle, but instead, collective charge and spin waves are the elementary excitations [29][30][31]. Electron interaction as a forward scattering renormalizes the velocities of the charge and spin modes separately, thus leading to a non-Fermi liquid.
Great efforts have been devoted to search for generalized Luttinger liquids in dimensions greater than one. While significant progress has been made in quasi-one dimensional systems [75,76], to our knowledge results are limited on Luttinger liquid type behavior in metals with truly two-dimensional Fermi surface.
Another situation for a non-Fermi liquid arises around a quantum critical point where strong electron interaction drives a phase transition from a metallic state to a symmetry-breaking ordered state at T = 0 [32][33][34][35]. Seminal works by Hertz [36], Moriya [37], and Millis [38] deal with the quantum critical phenomenon in itinerant magnets, which describe the coupling between electrons with a finite Fermi surface and bosonic fluctuations of an order parameter near a magnetic transition. In their theories, low-energy modes of electrons are integrated out to yield a nonlocal singular effective action for bosonic modes. This challenging problem has invoked intense work and considerable progress [39][40][41][42][43][44][45][46][47].
Unlike these non-Fermi liquids, the supermetal we found near a high-order VHS is obtained under weak electron interaction. It relies on the singular DOS instead of singular interaction. This feature enabled us to develop an analytically controlled theory of supermetal with local interaction, using a small parameter-the DOS exponent.
Quantum criticality in nodal semimetals also hosts a non-Fermi liquid [48,49]. The vanishing DOS allows unscreened long-range Coulomb interaction, which is dressed by the gapless electron spectrum. Possible platforms include graphene [50][51][52], pyrochlore iridates [53,54], and topological phase transition in two dimensions [55,56]. The RG procedure for semimetals has a similar spirit as ours for supermetal, both of which are Wilson-Fisher type instead Shankar type. Unlike semimetals though, a supermetal has divergent DOS at Fermi level and an extended Fermi surface.
We also note that a saddle point of an energy dispersion in two dimensions E k = k 2 x − k 2 y gives a conventional VHS with a logarithmic divergence mentioned above. This is also scale invariant; however, it requires an additional care in an RG analysis [19,20]. As we mentioned in Sec. IV, momentum integrals in a perturbative RG calculation suffer from a singularity at k → ∞. This cannot be regularized by the UV energy cutoff Λ since every energy contour extends to k → ∞, so that a UV momentum cutoff is needed in addition. This is also related to the non-analyticity of the DOS. As a result, the UV cutoff is not eliminated from RG equations. It occurs as a sequel that the low-energy physics is affected by the UV scale Λ.
We have listed several other non-Fermi liquid systems. Supermetal fixed points are regarded as multicritical points in the phase space spanned by the coupling constant, the chemical potential, and a parameter for the energy dispersion (Fig. 7). Importantly, our present analysis does not suffer from the difficulties related to a closed Fermi surface, the presence of a length scale, a logarithmic DOS, and a singular bosonic fluctuation. A small DOS singularity exponent guarantees that only the short-range interaction is relevant.

C. Discussion
Our analysis is for the case of a single high-order VHS in the Brillouin zone at the energy range in focus. In reality, materials may have multiple high-order VHS points at the same energy in the Brillouin zone, related by symmetry. In that case, additional interactions involving different VHS should be included and their presence may lead to symmetry-breaking instabilities [4], as opposed to quantum criticality. However, for a certain parameter range before an ordering instability takes place, there could exist a scaling region where thermodynamic or transport quantities follow scaling properties. For ex-ample, when temperature T and the carrier density n are control parameters, a physical quantity Q follows the scaling relation around the noninteracting fixed point, where the exponent a is determined by a dimensional analysis of Q and F is a scaling function.

Acknowledgments
We acknowledge helpful discussions with A. V. Chubukov, E. Fradkin, S.-S. Lee We work on the RG equations at finite temperature T = 0 with the energy-shell RG analysis to one-loop level. Temperature has a dimension of energy, so that it is one of relevant perturbations. The one-loop corrections Σ H , Π pp , and Π ph , shown in Eq. (41) at T = 0, should be calculated at T = 0. To order l, we obtain where we introduce the dimensionless temperatureT = T /Λ. Again, all quantities are evaluated at zero external frequency and momentum, so that the results depend only on the DOS. The temperature-dependent dimensionless coefficients c µ , c pp , and c ph are Unlike the calculation at T = 0, Π ph becomes finite for T = 0, while the correction to the field or the energy dispersion remains absent to one-loop order.
In the analysis at zero temperature, we rescale the frequency as Eq. (37b). At finite temperature, rescaling of the Matsubara frequency leads to rescaling of temperature [38]. Temperature obeys the same scaling relation as that for the frequency: T = bT .
Including the temperature-dependent factors c µ (T ), c pp (T ), and c ph (T ), we obtain the changes of the parameters at an RG step as Then, we reach the RG equations Since temperature T is relevant and its fixed point is located at T = 0, the fixed points of the parameters are found at T = 0, as we discussed in the main part. The RG equations are consistent with the beta functions derived by the field theory approach in Sec. VI D and VI E.
We show the derivation of the Ward-Takahashi identity, based on a diagrammatic discussion by Peskin and Schroeder [86]. It consists of two parts: equalities for a through-going line from an initial state to a final state and an internal loop. Suppose that we calculate a diagram with the same numbers of incoming and outgoing electron lines. Then we can decompose the diagram into lines that connect an incoming and an outgoing line, and internal loops. For the present analysis with the contact interaction, we choose each decomposed diagram to include only one spin species. We set T = 0 in the following discussion.
To derive the Ward-Takahashi identity, we insert an external line for the bosonic field ϕ(k) with k = (k, ω k ). The coupling between an electron and the bosonic field is given by where we introduce the spin-dependent coupling constant α σ .
Step 1. For a line connecting an initial state and a final state, the corresponding equation contains a product of the noninteracting Green's functions. When we consider an electron line with (M + 1) electron line segments, the product is written as with p J = p J−1 + q J , p 0 = p, p M = p , and {q} = (q 1 , q 2 , · · · , q M ) [ Fig. 11(a)]. Remember that all connected line segments has the spin index σ. Now we insert a vertex to the J-th line, connecting a bosonic line carrying momentum and frequency k [ Fig. 11(b)]. We denote it as L σ M,J (p , p; {q}; k), given by For k = (0, ω k ), the following equality holds: or equivalently Using the equality, we obtain If we insert a vertex on the (J − 1)-st line, we then have We find a cancellation of the second term in the brackets in Eq. (B6) by the first term in the brackets in Eq. (B7) when we sum the two. Summing all possible insertions of vertices, we find ↑ . It consists of a line L ↑ and a loop R ↓ . The wavy lines represent a bosonic field; in the present case, internal wavy lines can be regarded as a Hubbard-Stratonovich field for the contact interaction.
Step 2. We apply a similar argument for an internal loop. The difference from the discussion for a line is that internal frequency and momentum that run through the loop have to be integrated out. Thus, the equation corresponding to the loop diagram is given by Step 3. By combining the results of Step 1 and Step 2, we deduce the Ward-Takahashi identity. We introduce the N -point function Γ (N ) ({p }, {p}) and the associated vertex function Γ (N,ασ) ({p }, {p}; k), which is obtained by attaching an external bosonic line with frequency and momentum k to the N -point function Γ (N ) . {p} stands for the set of frequencies and momenta (p 1 , · · · , p N ). Frequency and momentum should be conserved before and after the process: p J = p J for Γ (N ) and p J = p J + k for Γ (N,gσ) . The N -point function can be decomposed into N/2 lines L σ and loops R σ . [For example, Fig. 11(d) shows a decomposition of a sunrise diagram for a two-point function (N = 2) into L ↑ and R ↓ .] Internal frequencies and momenta {q} are to be integrated out. Therefore, from Eqs. (B8) and (B11), we obtain (B12) We use the notation {p} σ J,+k = (p 1 , · · · , p J + k, · · · , p N ), where the addition of k occurs only when p J is associated with an electron with spin σ. {p} σ J,+k reflects the spin-dependent coupling α σ . The sign (−1) on the lefthand side is solely up to the definition of the vertex function. Equation (B12) is the Ward-Takahashi identity for the N -point function. For the present model, internal bosonic lines mediate the contact interaction and external lines represent an external field as a perturbation, such as the chemical potential (α ↑ = α ↓ = −µ) and the magnetic field (α ↑ = −α ↓ = h).
Here, we considered the case where the external bosonic field carries finite frequency only, so that the right-hand side of Eq. (B4) is simply proportional to the frequency. The extension to a case with a finite momentum is straightforward, but the right-hand side becomes not as simple as that for frequency, which depends on the energy dispersion.

Ward identity for two-point functions
The Ward-Takahashi identity originates from conservation laws: the expression Eq. (B12) obeys due to the charge conservation for each spin. The identity yields a relation between a self-energy and a vertex function, thus leading to a relation between the field renormalization and an exponent for a thermodynamic quantity. We can see this from a two-point function (N = 2). The Ward-Takahashi identity for N = 2 is written as where the subscript σ is added to show the external lines correspond to electrons with spin σ. The two-point function is written with the self-energy Σ σ as Therefore, in the limit ω k → 0, we obtain the Ward identity, which relates the vertex function Γ (2,ασ) σ and the self-energy Σ σ : Appendix C: Brief review of the field theory approach to RG analyses Here we describe a field theory approach to RG equations, in light of the Wilsonian approach. We derive the Callan-Symanzik equation and the beta functions, which show how scale-dependent parameters affect physical quantities. We partly owe the following descriptions to the references [85][86][87][88].
For the sake of clarity, we consider a theory with a scalar field φ and a set of dimensionless parameters {ḡ ρ }, where we write the action as S[φ;ḡ]. The partition function is given by When the model suffers from UV divergences, i.e., perturbative loop corrections have UV divergences, we need to cure them to extract meaningful information. Those UV divergences can be regularized by removing UV modes from the model. To this end, we decompose the field φ depending on the energy range to which they contribute: φ Λ Λ accounts for the energy between Λ and Λ. Then, we redefine the partition function as It does not obviously have a UV divergence because no UV modes are included. Here the energy scale Λ 0 works as a UV energy cutoff. (This is equivalent to impose the effective action to be finite at Λ 0 without a cutoff. To this end, divergent counterterms should be introduced to cure UV divergences.) The scale Λ 0 is an arbitrary energy scale to regularize UV divergences. The next thing we should check is how a change of the characteristic energy scale affects the theory. To see it, we define the effective action at an energy scale Λ(< Λ 0 ) as We require that the effective action S eff have the same form as the action S. Then, the partition function can be written as This is simply rewriting of the partition function with the effective action at the scale Λ. We relate the partition functions at different scales to find This equality tells us that we have the same partition function defined at different energy scales Λ and Λ , together with the changes of the weight Z and the parametersḡ ρ . We then aim to calculate the N -point correlation function with the cutoff Λ and the parametersḡ ρ (Λ): When all momenta k a correspond to energies below Λ and Λ , we find φ Λ 0 (k a ) = φ Λ 0 (k a ), which enables us to relate the N -point correlation functions at different scales as Z −N/2 Λ φ Λ 0 (k 1 ) · · · φ Λ 0 (k N ) Λ;ḡ(Λ) = Z −N/2 Λ φ Λ 0 (k 1 ) · · · φ Λ 0 (k N ) Λ ;ḡ(Λ ) . The scale dependence of this equality can be written in the form of a differential equation: Note that the repeated index is summed over. This equation is called the Callan-Symanzik equation [71][72][73] for the connected N -point correlation function G (N ) with the beta functions β ρ and the field renormalization γ defined by The correlation functions are obtained from perturbative calculations. Actually, it is rather straightforward to calculate the one-particle irreducible N -point function Γ (N ) instead of the N -point correlation function G (N ) . When a model involves a quartic interaction φ 4 without a cubic term φ 3 , Γ (2) and Γ (4) are given by Γ (4) (k 1 , k 2 , k 3 , k 4 ) = G (4) (k 1 , k 2 , k 3 , k 4 ) G (2) (k 1 )G (2) (k 2 )G (2) (k 3 )G (2) (k 4 ) .

(C12)
For the definition of Γ (N ) from the effective action, see the references [85][86][87][88]. Roughly speaking, Γ (N ) corresponds to the coefficient of the φ N term in the effective action. Again, using the fact that φ Λ 0 (k a ) = φ Λ 0 (k a ) holds when the energy corresponding to the momentum k a is smaller than Λ and Λ , we find the relation So far, we have compared G (N ) or Γ (N ) at different cutoffs Λ and Λ , so that the dependence on Λ is explicit.
However, this comparison is still theoretical; i.e., this is a comparison of different systems. Our aim is to compare the two theories with the same cutoff Λ. For this sake, we rescale the coordinate to change the cutoff. Suppose we have the scaling relations Λ = bΛ, (C15a) which can be interpreted in the following way. We integrate out fluctuations between the range of (Λ/b, Λ] (corresponding to Z N/2 Λ;ḡ(Λ) /Z N/2 Λ/b;ḡ(Λ/b) ) and rescale the field and parameters (multiplying the factor b d Γ (N ) ) to obtain the new action with a different coupling constant but with the same cutoff Λ.
We can also write down the Callan-Symanzik equation to describe the momentum dependence, instead of the cutoff Λ. We differentiate Eq. (C17) with respect to b and then set b = 1 to obtain d kj k a,j ∂ ∂k a,j + β ρ (ḡ) Λ;ḡ ({k a,j }) = 0.
The deviation η corresponds to the anomalous dimension.
Appendix D: Two-loop self-energy for the quasiparticle lifetime The quasiparticle damping is captured by a finite imaginary part of the self-energy Σ. In a series of perturbative expansions, the lowest order correction appears at second order, which is diagrammatically shown in Fig. 5(d). The result is shown in Eq. (130) and here we calculate it explicitly.
DOS; while χ ph has a logarithmic divergence reflecting the DOS, χ pp exhibits a double logarithmic divergence. The additional logarithm in the BCS channel appears in the presence of time-reversal symmetry, which can be regarded as the Fermi surface nesting with itself.