Stable non-Fermi liquid fixed point at the onset of incommensurate charge density wave order
Abstract
We consider the emergence of a non-Fermi liquid fixed point in a two-dimensional metal, at the onset of a quantum phase transition from a Fermi liquid state to an incommensurate charge density wave (CDW) ordered phase. The momentum of the CDW boson is centred at the wavevector , which connects a single pair of antipodal points on the Fermi surface with antiparallel tangent vectors. We employ the dimensional regularization technique in which the co-dimension of the Fermi surface is extended to a generic value, while keeping the dimension of the Fermi surface itself fixed at one. Although the system is strongly coupled at dimension , the interactions become marginal at the upper critical dimension , whose value is found to be . Using a controlled perturbative expansion in the parameter , we compute the critical exponents of the stable infrared fixed point characterizing the quantum critical point. The scalings of the original theory are determined by setting , where the fermion self-energy is seen to scale with the frequency with a fractional power law of , which is the telltale signature of a typical non-Fermi liquid phase.
Contents
I Introduction
While Landau’s Fermi liquid theory has been incredibly successful in describing normal metals, there exists an extensive number of metallic states where the framework fails. Although of widely different origins, these systems are widely known as non-Fermi liquids. A finite density of nonrelativistic fermions interacting with transverse gauge field bosons was the first model, considered by Holstein et al. [1] to study the effects of the electromagnetic fields on a metal, to exhibit a non-Fermi liquid character. Subsequently, it was realized that similar behaviour is found for finite-density fermions coupled with massless order parameter bosons at quantum critical points [2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31] or artificial gauge field(s) emerging in various kinds of scenarios [32, 33, 34, 35, 36, 37, 38, 39]. While the Landau quasiparticles provide a natural single-particle basis for writing down the corresponding low-energy quantum field theories (QFTs) on the merit of being long-lived excitations in a Landau Fermi liquid, they are destroyed by the strong interactions between the soft fluctuations of the Fermi surface and the gapless bosonic quantum fields in the scenarios described above. Since these are fundamentally strongly-interacting theories, it is a challenging task to build a controlled approximation amenable to theoretical analysis. As a result, there have been intensive efforts to devise QFT frameworks to explain the emergent physical characteristics of these non-Fermi liquid systems [1, 40, 35, 41, 42, 43, 39, 17, 44, 38, 7, 45, 25, 2, 26, 27, 46, 47, 48, 49, 50, 19, 28, 20, 21, 22, 23, 24, 29, 30, 31, 51, 52, 53, 54]. In two spatial dimensions, the corresponding theories are genuinely strongly interacting, while in three dimensions, they emerge as marginal Fermi liquids [20, 53]. Analogous situations arise in two-dimensional (2d) and three-dimensional (3d) semimetals, where a non-Fermi liquid state emerges when the chemical potential cuts a band-crossing point giving rise to a Fermi point (rather than a Fermi surface) and a long-ranged (i.e., unscreened) Coulomb potential is switched on [55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65].
For non-Fermi liquids arising at quantum critical points, there are two types of order parameter bosons: (1) cases where the critical bosonic field is centred about zero momentum, inducing the quasiparticles to lose coherence across the entire Fermi surface [2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24]; (2) the momentum of the quantum field describing the bosonic degrees freedom is centred around a finite wavevector equal to , which connects points on the Fermi surface (commonly known as hot-spots), such that the non-Fermi liquid behaviour emerges locally in the vicinity of the hot-spots [2, 26, 27, 28, 29, 30, 31, 52]. Examples in the first category include the Ising-nematic critical point, while the second one include ordering transitions to phases like spin density wave (SDW), charge density wave (CDW) [25, 26, 27, 28, 29, 30, 31], and the Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) states [52].
The CDW (SDW) bosons with give rise to instabilitites involving charge (magnetic) order, where the charge (spin) density spontaneously breaks translational symmetries and develops a density modulation equal to . There are two distinct categorizations depending on the nature of the wavevector: (1) whether it is commensurate or incommensurate; (2) whether is a nesting vector of the Fermi surface or not. While the so-called commensurate wavevectors can be written as a linear combination of the reciprocal lattice vectors with rational coefficients, we cannot do the same for the incommensurate wavevectors. Clearly, there are infinitely many possible rational coefficients, but the quantitative effects of commensurability decrease with the size of their denominators. For the special situation of equalling a nesting vector connecting two points on the Fermi surface with antiparallel Fermi velocities (or, equivalently, tangent vectors), the spin and charge orderings feature a well-known singularity caused by an enhanced phase space for low-energy particle-hole excitations. In an inversion-symmetric crystal with the valence band dispersion , the nesting vectors are given by the condition , where is the Fermi energy. The nesting-vector nature of , coupled with inversion symmetry, implies that , where is the local radius of curvature/reciprocal of curvature of the Fermi surface (i.e., the magnitude of the local Fermi momentum vector). This results from the fact that the two hot-spots are related by inversion symmetry and, hence, both have the same value of . The -wavevector instabilities are ubiquitous in 2d systems exhibiting high-temperature superconductivity — for example, (1) the ground state of the 2d Hubbard model at half-filling of the conduction band exhibits an SDW instability at a -wavevector [66, 67]; (2) d-wave bond charge order, triggered by antiferromagnetic fluctuations, in models for cuprate superconductors occurs naturally at -wavevectors [25, 68]. We would like to emphasize that such a nesting vector causes a partial nesting of the Fermi surface, which is distinct from the perfect-nesting cases when slices (and not just discrete points) of a Fermi surface are connected by the same nesting vector. The evidence of incommensurate CDW orderings has been reported in materials like NbSe2 and TaS2 [69, 70], VSe2 [71], SmTe3 [72], and TbTe3 [73]. In some of these compounds, the CDW transition temperature can be tuned close to zero Kelvin by applying high pressure, unravelling a putative quantum critical point at the onset of the CDW ordering [74]. Such observations indicate the importance of a theoretical understanding of quantum critical points involving incommensurate -wavevector instabilities.
In this paper, we consider a pair of antipodal points on a 1d Fermi surface (of a 2d metal) with parallel tangent vectors interacting with an order parameter boson, whose condensation gives rise to an incommensurate CDW ordered phase [75, 76, 77]. The CDW boson becomes massless right at the quantum critical point, giving rise to strong quantum fluctuations. The physical picture here is that the CDW boson drives the system across a quantum phase transition to an ordered state, where the electron density spontaneously breaks translational symmetries and develops a density modulation with a wave vector that is incommensurate with the underlying reciprocal lattice vectors. For the sake of definiteness, we choose , without any loss of generality.
To deal with the QFT describing the above system, we implement the analytic approach of dimensional regularization [78, 19, 20, 21, 23, 22, 52, 53], in which the co-dimension of the Fermi surface is increased (as a mathematical tool) with the aim to obtain the value of the upper critical dimension . Since is the dimension at which the interactions become marginal, we succeed in formulating a controlled perturbative approximation, although the quasiparticle-description has broken down. The critical exponents and various physical properties can now be calculated in a systematic expansion involving the perturbative parameter . Halbinger et al. [77] have considered this incommensurate CDW problem by employing the same methodology, and have found that the fermion-boson interactions lead to a stable non-Fermi liquid fixed point. Furthermore, their results show that the resulting critical Fermi surface is flattened at the hot-spots. However, they arrived at their conclusions from some inaccurate computations, which, in addition to other results, could not find the correct critical dimension . They carried out their computations assuming that is equal to , predicting that it would be obtained at the two-loop order. Consequently, their answers do not include the all-important frequency depedence of for the fermion self-energy, which is predicted in the earlier (uncontrolled) random phase approximation (RPA) calculations [75, 76]. In order to address these inadequacies, the analysis via the dimensional regularization scheme must be reexamined, if only to put it on a firmer basis.
The paper is organized as follows. In Sec. II, we introduce the effective low-energy Euclidean action in the Matsubara frequency space, describing the fermionic excitations at the two hot-spots of the Fermi surface interacting with the CDW boson fluctuations. The original theory is embedded in spatial dimensions. We also explain how to generalize the theory to a generic value of by increasing the number of dimensions perpendicular to the 1d Fermi surface, as this will allow us to identify the upper critical dimension and, subsequently, to regularize the theory via the renormalization group (RG) procedure. Sec. III is devoted to the derivations of the bosonic and fermionic self-energies, and showing that comes out to be . Using the one-loop results, the RG flow equations are determined in Sec. IV by absorbing the ultraviolet (UV) divergences into singular counterterms. We discuss the nature of the fixed points in the infrared (IR) limit and show that the system flows to a stable non-Fermi liquid point. Finally, we conclude with the relevant discussions in Sec. V, comparing our computations and results with earlier works. The appendix shows a part of the computations of the fermion self-energy.
II Model
We consider the low-energy QFT action describing finite-density fermions confined to two spatial dimensions, and interacting with an incommensurate CDW order parameter with momentum centred at . Thus, the nesting vector connects two hot-spots on the Fermi surface located along the -axis, as shown in Fig. 1. In -dimensions, the effective action describing the electrons near the hot-spots and the CDW order parameter mode is given by [77, 75, 76]
(1) |
where denotes the three-vector comprising the Matsubara space frequency and the spatial momentum vector , , and is the number of spatial dimensions . The fermionic degrees of freedom about the right and left hot-spots are denoted by and , respectively. The fields and refer to the bosonic fluctuations carrying frequency and momenta and , respectively. The bosons are massless as we are considering the quantum critical point. To simplify notations, we have rescaled the fermionic momenta in a way such that the absolute value of the Fermi velocity is unity and the curvature of the Fermi surface is equal to at the hot-spots. Although the bosonic velocity is in general distinct from that of the fermions, we have set the bare velocity of the bosons equal to unity as well. This is because the dynamics of the bosons at the critical point is dominated by the particle-hole excitations of the Fermi surface at low energies, and the actual value of the bosonic velocity does not matter in the low-energy effective theory.
In our action, since the Fermi surface is locally parabolic, we can set the scaling dimensions of and are equal to and , respectively. In order to extract the critical scalings in a controlled approximation, we increase the co-dimensions of the Fermi surface [78, 19, 28] to eventually determine the upper critical dimension , where the fermion self-energy shows a logarithmic singularity. To preserve the analyticity of the theory in momentum space (alternatively, locality in real space) with generic co-dimensions, we introduce the two-component “spinors” [19, 20, 21, 52, 53, 54]
(2) |
and write an action that describes the 1d Fermi surface embedded in a -dimensional momentum space:
(3) |
The -component vector includes the frequency and the -components of the momentum vector due to the added co-dimensions. The original momentum components along the - and -directions have been relabelled as and , respectively. Hence, in the resulting -dimensional momentum space, the set of components represents the the direction perpendicular to the Fermi surface, while is along the parallel direction. Similarly, the vector of matrices has components representing the gamma matrices associated with and the extra co-dimensions. Ultimately, we are interested in continuing to , which implies that, in practice, it is sufficient to consider only the gamma matrices and in our computations.
In the purely bosonic part of the action, only the part of the kinetic term is retained, because is irrelevant under the scaling of the patch-theory formalism [2, 19, 20, 21, 77, 52, 53], where each of has dimension unity and has dimension . Dependence on in the bosonic propagator will be generated via the susceptibility, which is obtained from the dynamics of the strong particle-hole fluctuations. Anticipating this, we have already added the extra term , which will be generated in the RG process, and its has been dictated by the divergent term in the one-loop susceptibility calculated in the following subsection [see Eq. (13)]. In other words, we have simply included a term which will be generated via quantum corrections. This term has a mass dimension equal to unity, similar to the term, with a vanishing engineering dimension for . If we do not include this term, the loop integrations involving the bosonic propagator will turn out to be infrared-divergent, and these divergences will be the mere artifacts of our dropping the irrelevant terms in the minimal local effective action (if it contains only the term). Finally, the engineering dimension of the fermion-boson coupling is equal to — this observation has dictated us to introduce an explicit factor of a mass scale raised to the power , chosen so as to ensure that is dimensionless, which is the usual procedure followed in QFT calculations.
The emergent sliding symmetry in the Ising-nematic case [2, 19, 20, 21] forces the terms proportional to and in Eq. (II) to renormalize in the same way. In other words, the fermion propagator depends on and only through even after loop corrections. However, that is not the case here, and the renormalization process is not guaranteed to retain the sole dependence on . In other words, it is possible that the nature of the renormalized terms may turn out to be such that it leads to a flattening of the Fermi surface at the hot-spots, as found in the RPA calculations of Ref. [76]. Nevertheless, as we will see from our explicit calculations, the flattening does not show up in our one-loop level computations.
III One-loop self-energies and dimensional regularization
The value of tells us that the coupling constant becomes marginal at the upper critical dimension . In other words, is relevant for and irrelevant for . Our aim is to access the interacting phase perturbatively in , using as the perturbative parameter. In particular, this implies that at the end of our systematic -expansion, we have to set for our original two-dimensional theory. Before embarking on deriving the RG flows, in this section, we compute the one-loop self-energies for the bosonic and fermionic degrees of freedom, which will feed into the equations necessary for determining the beta functions of the coupling constants and .
The bare fermion and boson propagators for the action defined in Eq. (II) are given by
(4) |
III.1 One-loop boson self-energy
The one-loop boson self-energy [cf. Fig. 2] is defined by
(5) |
Using the commutation relations between the gamma-matrices, and the identities and , we get
(6) |
Noting that , we first shift , and then use the Feynman parametrization to obtain
(7) |
Changing variables to with the Jacobian factor , we get
(8) |
Since the bare susceptibility at zero temperature diverges at the nesting vector , the well-defined self-energy is given by subtracting off the singular contribution for this momentum, which translates to
(9) |
where (after changing variables as )
(10) |
Note that the zeroes coming from and at are cancelled by the factor present in Eq. (9).
The leading order terms obtained in the limit are found to be:
(11) |
Finally, expanding in for , we get
(12) |
where is the Euler–Mascheroni constant. It is important to note here that the loop integral has an upper critical dimension at and, as a result, it shows a pole in also at . As such, these terms in the one-loop corrected action will be irrelevant, which is reflected by the that they have a scaling dimension equal to at .
From our analysis above, we conclude that the self-energy correction in Eq. (9) is given by
(13) |
where
(14) |
We note that the bare boson propagators are still independent of , the loop integrations involving it are ill-defined unless one resums a series of diagrams that provides a nontrivial dispersion along these frequency and momentum components. Hence, in all loop calculations involving the boson propagators, we include the lowest order fine correction from the one-loop self-energy, which is proportional to . Therefore, both for the and bosonic fields, we use the dressed propagator
(15) |
which is equivalent to rearranging the perturbative loop-expansions such that the one-loop finite part of the boson self-energy, dependent on , is included at the zeroth order. We would like to emphasize that is the so-called Landau-damped term which leads to the signature -dependence of the fermion self-energy, characterizing the non-Fermi liquid behaviour in various quantum critical systems [2, 25, 19, 20, 21, 77, 52, 53]. The Landau-damped part also plays the most significant role in inducing unconventional superconductivity in this kind of non-Fermi liquid systems [48, 49, 79, 23, 24].
III.2 One-loop fermion self-energy
The fermion self-energy [cf. Fig. 2] is given by the integral
(16) |
where
(17) |
and
(18) |
The steps to compute these two parts have been explained in the next two subsections, which can be skipped if the reader is not interested in the tedious intermediate steps. For their benefit, we state here the final result. Setting , we get the singular part to be
(19) |
where the divergence is parametrized by a pole at .
III.2.1 Computation of -dependent part
The leading order dependence of on can be extracted by setting the external momentum components and to zero. Hence, we will evaluate
(20) |
Changing the description to and as integration variables, and dividing into the parts and as , we have
(21) |
and
(22) |
The integral cannot be evaluated exactly and we need to make some reasonable approximations to extract the leading order corrections. We note that the first factor of the integrand tells us that the dominant contribution is concentrated around the region and . As for the second factor, the dominant contribution comes from . Since for small close to zero and , we can substitute in the factor in the overall denominator and the term in the denominator of the second factor, and extend the lower limit of the integral over to , leading to
(23) |
The integral blows up at , which thus gives us the value of the upper critical dimension . The fermion-boson coupling is irrelevant for , relevant for , and marginal for . This allows us to access the strongly interacting non-Fermi liquid state perturbatively, in a controlled approximation, using , where serves as the perturbative/small parameter. In our dimensional regularization scheme, the divergence appears as , with the factor having a pole at . We also note that this term produces the behaviour of the fermion self-energy as at , which matches with the uncontrolled RPA result [75, 76]. We would like to point out that the correct -dependence of could be captured only because we have included the crucial Landau damping term in the dressed bosonic propagator . The -scaling was missed in Ref. [77] due to the non-inclusion of .
III.2.2 Computation of -dependent part
The leading order dependence of on and can be extracted by setting . Hence, we will evaluate
(24) |
where . We note that the first factor of the integrand sets the condition that the dominant contributions must come when . This immediately tells us that, in the denominator of the second factor, the term proportional to is subleading. Hence, to extract the leading order contribution we can set it to zero. In the appendix, we show that the subleading order term turns out to be proportional to and is not singular at . Therefore, we calculate the approximate integral
(25) |
Using the identity
(26) |
for the Cauchy principal value of the integral, the -integral is performed first to obtain the form
(27) |
The factor has a pole at , which shows that this term has (1) a logarithmic divergence at , and (2) a linear divergence at , when we translate the divergences in the language of the Wilsonian cutoff . We need to treat this result carefully by remembering that, in the dimensional regularization procedure, UV divergences of all degrees show up as the poles of -functions. The degree of divergence can be understood by using an explicit UV cut-off in the language of the Wilsonian RG, which is denoted here by . Although this term will play an important role for analyzing UV-stable fixed points, it should be discarded here because we are considering the RG flows in the IR, and this term represents an IR-irrelevant operator for the theory in .
Some more comments are in order. The situation above is similar to having a -term in a scalar field theory in -dimensions. A simple power counting of momenta arising in loop diagrams shows that adding a -interaction vertex to a free scalar field theory gives an upper critical dimension of , while addition of a -vertex has the upper critical dimension value of . Therefore, adding a -vertex to the theory in four spacetime dimensions makes the theory nonrenormalizable.
The above statements can also be quantified in an alternate way. Let us expand the expression in , where . Then we have
(28) |
which indicates that the term has a linear UV divergence at . Another crucial observation is that vanishes identically at . All these observations together dictate that there will be no contribution to the RG flows in the IR from this term and, hence, should not be included in the counterterms. Consequently, the flattening of the Fermi surface at the hot-spots, found in the RPA calculations [76], does not show up in our one-loop results.
III.3 One-loop vertex correction
It is not possible to get a one-loop vertex diagram and, hence, the corresponding correction is zero.
IV Renormalization Group flows under minimal subtraction scheme
Eq. (II) is supposed to be the physical action, defined at an energy scale , consisting of the fundamental Lagrangian with non-divergent quantities. However, we have seen that the loop integrals lead to divergent terms and, in order to cure it, we employ the renormalization procedure, using dimensional regularization as the regularization method. In our dimensional regularization formalism, the UV-divergent terms are the ones arising in the limit. We use the minimal subtraction () renormalization scheme to control the UV divergences [80, 81], which involves cancelling the divergent parts of the loop-contributions via adding appropriate counterterms. More precisely, we adopt the modified minimal subtraction () scheme where, in addition to the divergent term, we absorb the universal term proportional to (that always accompanies the term with the pole) into the corresponding counterterm.
The action, consisting of the counterterms to absorb the singular terms, takes the form:
(29) |
The counterterm-factors are given by the power series
(30) |
such that they cancel the divergent contributions from the Feynman diagrams. Due to the -dimensional rotational invariance in the space perpendicular to the Fermi surface, each term in is renormalized in the same way.
Subtracting from the so-called bare action , we obtain the renormalized action, which is the physical effective action of the theory, re-written in terms of non-divergent quantum parameters. While the bare parameters can be divergent, the physical observables are the renormalized coupling constants, which are determined by the RG equations. The RG flows describe the evolution of the bare couplings as functions of the floating energy scale (i.e., with respect to an increasing logarithmic length scale ). To achieve this objective, we first define the bare (or fundamental) action
(31) |
consisting of the bare quantities, where the superscript “” has been used to denote the bare fields, couplings, frequency, and momenta. We now relate the bare quantities to the so-called renormalized quantities (without the superscript “”) via the multiplicative -factors such that
(32) |
(33) |
and
(34) |
Observing that there exists a freedom to change the renormalization of the fields and the renormalization of momenta without affecting the action, we have exploited it by requiring , which is equivalent to measuring the scaling dimensions of all the other quantities relative to the scaling dimension of . now represents the renormalized action (also known as the Wilsonian effective action) because it consists of the renormalized quantities. Basically, we have written the fundamental action of our theory in two different ways [82], which allows us to subtract off the divergent parts (represented by ).
IV.1 RG flow equations from one-loop results
At one-loop order, the divergent contributions are obtained from Eqs.(13) and (19). These lead to
(35) |
To this leading order correction, we find that , and they do not get any correction from the loop integrals.
Because , we define a single dynamical critical exponent for the fermions as
(36) |
This applies to our one-loop level calculations where the -part as a whole is not renormalized. Furthermore, the anomalous dimensions for the fermions and the bosons are given by
(37) |
respectively. We also define the beta functions for the two coupling constants as
(38) |
The sole purpose of the introduction of the ad hoc mass scale is to regularize the theory, thus eliminating the infinities emerging from the loop integrals of Feynman diagrams. However, since physical quantities must be independent of , as is not really a parameter of the fundamental theory, the bare parameters must be independent of it as well. Imposing this condition, as well as the requirement that the regular (i.e., non-singular) parts of the final solutions are of the forms
(39) |
in the limit , we get the following differential equations:
(40) |
The above set of equations have been obtained by (1) demanding that ; (2) plugging in the values from Eqs. (IV.1) and (39); (3) expanding in powers of ; and (4) matching the coefficients of the regular powers of in the resulting equations.
Solving the equations, we get
(41) | ||||
(42) | ||||
(43) |
where
(44) |
Since we are interested in the behaviour at the IR energy scales, we determine the RG flows with respect to the logarithmic length scale , which are given by the derivatives
(45) |
for the two coupling constants and , respectively.
IV.2 Stable fixed points
The fixed points of the theory are obtained as the points where the two beta functions, viz. and , go to zero simultaneously. The nomenclature originates from the fact that these are the equilibrium points of the differential equations describing the RG flows in the space of the two coupling constants, and . Due to the complicated form of the beta functions, it is not feasible to obtain closed-form expressions for these fixed points. Hence, we find their values numerically for a given value of . In order to determine the stability of a fixed point, one needs to figure out whether the flow lines (in the IR), given by the vector field with the components in the - plane, are towards or away from it. Accordingly, they are classified as stable or unstable. The values of , , and at the stable fixed points, as functions of , are shown in Fig. 3. In the range , we find precisely one stable fixed point for each value of . The stable nature of these fixed points has been illustrated by the RG flows for and in Fig. 4, which is obvious from the fact that the RG flow trajectories in each subfigure converge towards the fixed point while flowing from a high value of the floating mass scale to lower and lower values.
Since at the one-loop order, we do not find any flattening of the Fermi surface at the hot-spots, unlike the interpretation of Ref. [77], where the authors performed the RG using the pole of at . Furthermore, we would like to point out that both the fermionic and the bosonic anomalous dimensions remain equal, each taking a negative value (see the last panel of Fig. 3) in the range , unlike the results found in Ref. [77].
V Discussions and outlook
In this paper, we have revisited the QFT of the quantum critical point emerging at the continuous phase transition from a normal metal phase to an ordered phase involving an incommensurate CDW modulation. In our one-loop computations, we have used a dressed boson propagator by including the Landau-damping correction , which is instrumental in inducing the non-Fermi liquid behaviour with a characteristic scaling of for the fermion self-energy [19, 20, 21, 52, 53, 54, 83]. In Ref. [77], the authors computed the fermion self-energy by assuming the ad hoc form of the boson-self energy to be , which they argued is generated in RG due to symmetry arguments. However, a careful analysis shows that the RG procedure generates a term proportional to and, additionally, terms with Landau damping [cf. Eq. (13)]. As a consequence of not including the all-important Landau-damped term in the dressed bosonic propagator, they did not obtain any contribution in proportional to and, hence, missed the crucial -dependence at . They concluded that the correct frequency-dependence would show up at two-loop order. Furthermore, Halbinger et al. included the part proportional to in the one-loop fermion self-energy in the counterterms, which leads to their conclusion of the flattening of the Fermi surface at the hot-spots (as found in the RPA calculations of Ref. [76]). However, that term has a factor of , which first blows up at — hence, it should not contribute to the beta-functions for the flows towards IR, when is equal to . Their integrals are actually somewhat similar to the scenario for our calculation, and we have argued in detail why that term should not be included while computing the RG flows.
While we have addressed here the behaviour of a system with two hot-spots in the incommensurate CDW setting in two spatial dimensions, it will be worthwhile to extend it to the case of three dimensions [84], where we expect a marginal Fermi liquid behaviour, analogous to the results found in Refs. [20, 21, 53]. Another interesting direction to investigate is the scenario when the Fermi surface harbours two pairs of hot-spots [85]. Last, but not the least, we would like to compute the nature of superconducting instabilities in the presence of these critical CDW bosons, utilizing the RG set-ups constructed in Refs. [23, 59].
Acknowledgments
We are grateful for useful discussions with Sung-Sik Lee, Ashoke Sen, and Dimitri Pimenov. We thank Peng Rao for participating in the initial stages of the calculations. This research, leading to the results reported, has received funding from the European Union’s Horizon 2020 research and innovation programme under the Marie Skłodowska-Curie grant agreement number 754340.
Appendix: Details of calculating the -dependent part of
In this appendix, we compute the -dependent part of the fermion self-energy contribution from [cf. Eqs. (16) and (18)]. Starting with Eq. (24), we define
(46) |
where
(47) |
and
(48) |
We note that . Since the part is included in , whose computation has been detailed in the main text itself, we focus only on here.
To simply the calculations a bit, we set to zero to obtain
(49) |
We observe that the first factor of the integrand forces the dominant contribution to be concentrated around the region and . As for the second factor, the dominant contribution comes from . Hence, we can substitute in the factor (in the overall denominator) and the term (in the denominator of the second factor), and extend the lower limit of the integral over to . This leads to
(50) |
Hence, to leading order, the integral vanishes, and there remains no -dependent term.
References
- Holstein et al. [1973] T. Holstein, R. E. Norton, and P. Pincus, de Haas-van Alphen effect and the specific heat of an electron gas, Phys. Rev. B 8, 2649 (1973).
- Metlitski and Sachdev [2010a] M. A. Metlitski and S. Sachdev, Quantum phase transitions of metals in two spatial dimensions. I. Ising-nematic order, Phys. Rev. B 82, 075127 (2010a).
- Oganesyan et al. [2001] V. Oganesyan, S. A. Kivelson, and E. Fradkin, Quantum theory of a nematic Fermi fluid, Phys. Rev. B 64, 195109 (2001).
- Metzner et al. [2003] W. Metzner, D. Rohe, and S. Andergassen, Soft Fermi surfaces and breakdown of Fermi-liquid behavior, Phys. Rev. Lett. 91, 066402 (2003).
- Dell’Anna and Metzner [2006] L. Dell’Anna and W. Metzner, Fermi surface fluctuations and single electron excitations near Pomeranchuk instability in two dimensions, Phys. Rev. B 73, 045127 (2006).
- Kee et al. [2003] H.-Y. Kee, E. H. Kim, and C.-H. Chung, Signatures of an electronic nematic phase at the isotropic-nematic phase transition, Phys. Rev. B 68, 245109 (2003).
- Lawler et al. [2006] M. J. Lawler, D. G. Barci, V. Fernández, E. Fradkin, and L. Oxman, Nonperturbative behavior of the quantum phase transition to a nematic Fermi fluid, Phys. Rev. B 73, 085101 (2006).
- Rech et al. [2006] J. Rech, C. Pépin, and A. V. Chubukov, Quantum critical behavior in itinerant electron systems: Eliashberg theory and instability of a ferromagnetic quantum critical point, Phys. Rev. B 74, 195126 (2006).
- Wölfle and Rosch [2007] P. Wölfle and A. Rosch, Fermi liquid near a quantum critical point, Journal of Low Temperature Physics 147, 165 (2007).
- Maslov and Chubukov [2010] D. L. Maslov and A. V. Chubukov, Fermi liquid near Pomeranchuk quantum criticality, Phys. Rev. B 81, 045110 (2010).
- Quintanilla and Schofield [2006] J. Quintanilla and A. J. Schofield, Pomeranchuk and topological Fermi surface instabilities from central interactions, Phys. Rev. B 74, 115126 (2006).
- Yamase and Kohno [2000] H. Yamase and H. Kohno, Instability toward formation of quasi-one-dimensional Fermi surface in two-dimensional t-J model, Journal of the Physical Society of Japan 69, 2151 (2000).
- Yamase et al. [2005] H. Yamase, V. Oganesyan, and W. Metzner, Mean-field theory for symmetry-breaking Fermi surface deformations on a square lattice, Phys. Rev. B 72, 035114 (2005).
- Halboth and Metzner [2000] C. J. Halboth and W. Metzner, d-wave superconductivity and Pomeranchuk instability in the two-dimensional Hubbard model, Phys. Rev. Lett. 85, 5162 (2000).
- Jakubczyk et al. [2008] P. Jakubczyk, P. Strack, A. A. Katanin, and W. Metzner, Renormalization group for phases with broken discrete symmetry near quantum critical points, Phys. Rev. B 77, 195120 (2008).
- Zacharias et al. [2009] M. Zacharias, P. Wölfle, and M. Garst, Multiscale quantum criticality: Pomeranchuk instability in isotropic metals, Phys. Rev. B 80, 165116 (2009).
- Kim et al. [2008] E.-A. Kim, M. J. Lawler, P. Oreto, S. Sachdev, E. Fradkin, and S. A. Kivelson, Theory of the nodal nematic quantum phase transition in superconductors, Phys. Rev. B 77, 184514 (2008).
- Huh and Sachdev [2008] Y. Huh and S. Sachdev, Renormalization group theory of nematic ordering in d-wave superconductors, Phys. Rev. B 78, 064512 (2008).
- Dalidovich and Lee [2013] D. Dalidovich and S.-S. Lee, Perturbative non-Fermi liquids from dimensional regularization, Phys. Rev. B 88, 245106 (2013).
- Mandal and Lee [2015] I. Mandal and S.-S. Lee, Ultraviolet/infrared mixing in non-Fermi liquids, Phys. Rev. B 92, 035141 (2015).
- Mandal [2016a] I. Mandal, UV/IR mixing in non-Fermi liquids: higher-loop corrections in different energy ranges, The European Physical Journal B 89, 278 (2016a).
- Eberlein et al. [2016] A. Eberlein, I. Mandal, and S. Sachdev, Hyperscaling violation at the Ising-nematic quantum critical point in two-dimensional metals, Phys. Rev. B 94, 045133 (2016).
- Mandal [2016b] I. Mandal, Superconducting instability in non-Fermi liquids, Phys. Rev. B 94, 115138 (2016b).
- Mandal [2024] I. Mandal, Erratum: Superconducting instability in non-Fermi liquids [Phys. Rev. B 94, 115138 (2016)], Phys. Rev. B 109, 079902 (2024).
- Metlitski and Sachdev [2010b] M. A. Metlitski and S. Sachdev, Quantum phase transitions of metals in two spatial dimensions. II. Spin density wave order, Phys. Rev. B 82, 075128 (2010b).
- Abanov and Chubukov [2004] A. Abanov and A. Chubukov, Anomalous scaling at the quantum critical point in itinerant antiferromagnets, Physical Review Letters 93, 255702 (2004).
- Abanov and Chubukov [2000] A. Abanov and A. V. Chubukov, Spin-Fermion model near the quantum critical point: One-loop renormalization group results, Physical Review Letters 84, 5608 (2000).
- Sur and Lee [2015] S. Sur and S.-S. Lee, Quasilocal strange metal, Phys. Rev. B 91, 125136 (2015).
- Mandal [2017] I. Mandal, Scaling behaviour and superconducting instability in anisotropic non-Fermi liquids, Annals of Physics 376, 89 (2017).
- Schlief et al. [2017] A. Schlief, P. Lunts, and S.-S. Lee, Exact critical exponents for the antiferromagnetic quantum critical metal in two dimensions, Phys. Rev. X 7, 021010 (2017).
- Lunts et al. [2017] P. Lunts, A. Schlief, and S.-S. Lee, Emergence of a control parameter for the antiferromagnetic quantum critical metal, Phys. Rev. B 95, 245109 (2017).
- Baskaran and Anderson [1988] G. Baskaran and P. W. Anderson, Gauge theory of high-temperature superconductors and strongly correlated Fermi systems, Phys. Rev. B 37, 580 (1988).
- Ioffe and Larkin [1989] L. B. Ioffe and A. I. Larkin, Gapless fermions and gauge fields in dielectrics, Phys. Rev. B 39, 8988 (1989).
- Lee [1989] P. A. Lee, Gauge field, Aharonov-Bohm flux, and high-Tc superconductivity, Phys. Rev. Lett. 63, 680 (1989).
- Lee and Nagaosa [1992] P. A. Lee and N. Nagaosa, Gauge theory of the normal state of high-Tc superconductors, Phys. Rev. B 46, 5621 (1992).
- Blok and Monien [1993] B. Blok and H. Monien, Gauge theories of high-Tc superconductors, Phys. Rev. B 47, 3454 (1993).
- Ubbens and Lee [1994] M. U. Ubbens and P. A. Lee, Superconductivity phase diagram in the gauge-field description of the t-J model, Phys. Rev. B 49, 6853 (1994).
- Nayak and Wilczek [1994a] C. Nayak and F. Wilczek, Non-Fermi liquid fixed point in 2 + 1 dimensions, Nuclear Physics B 417, 359 (1994a).
- Chakravarty et al. [1995] S. Chakravarty, R. E. Norton, and O. F. Syljuåsen, Transverse gauge interactions and the vanquished Fermi liquid, Physical Review Letters 74, 1423 (1995).
- Reizer [1989] M. Y. Reizer, Relativistic effects in the electron density of states, specific heat, and the electron spectrum of normal metals, Phys. Rev. B 40, 11571 (1989).
- Halperin et al. [1993] B. I. Halperin, P. A. Lee, and N. Read, Theory of the half-filled Landau level, Phys. Rev. B 47, 7312 (1993).
- Polchinski [1994] J. Polchinski, Low-energy dynamics of the spinon-gauge system, Nuclear Physics B 422, 617 (1994).
- Altshuler et al. [1994] B. L. Altshuler, L. B. Ioffe, and A. J. Millis, Low-energy properties of fermions with singular interactions, Phys. Rev. B 50, 14048 (1994).
- Nayak and Wilczek [1994b] C. Nayak and F. Wilczek, Renormalization group approach to low temperature properties of a non-Fermi liquid metal, Nuclear Physics B 430, 534 (1994b).
- Lee [2009] S.-S. Lee, Low-energy effective theory of Fermi surface coupled with U(1) gauge field in dimensions, Phys. Rev. B 80, 165102 (2009).
- Mross et al. [2010] D. F. Mross, J. McGreevy, H. Liu, and T. Senthil, Controlled expansion for certain non-Fermi-liquid metals, Phys. Rev. B 82, 045121 (2010).
- Jiang et al. [2013] H.-C. Jiang, M. S. Block, R. V. Mishmash, J. R. Garrison, D. N. Sheng, O. I. Motrunich, and M. P. A. Fisher, Non-Fermi-liquid d-wave metal phase of strongly interacting electrons, Nature (London) 493, 39 (2013).
- Chung et al. [2013] S. B. Chung, I. Mandal, S. Raghu, and S. Chakravarty, Higher angular momentum pairing from transverse gauge interactions, Phys. Rev. B 88, 045127 (2013).
- Wang et al. [2014] Z. Wang, I. Mandal, S. B. Chung, and S. Chakravarty, Pairing in half-filled Landau level, Annals of Physics 351, 727 (2014).
- Sur and Lee [2014] S. Sur and S.-S. Lee, Chiral non-Fermi liquids, Phys. Rev. B 90, 045121 (2014).
- Lee [2018] S.-S. Lee, Recent developments in non-Fermi liquid theory, Annual Review of Condensed Matter Physics 9, 227–244 (2018).
- Pimenov et al. [2018] D. Pimenov, I. Mandal, F. Piazza, and M. Punk, Non-Fermi liquid at the FFLO quantum critical point, Phys. Rev. B 98, 024510 (2018).
- Mandal [2020] I. Mandal, Critical Fermi surfaces in generic dimensions arising from transverse gauge field interactions, Phys. Rev. Research 2, 043277 (2020).
- Mandal and Fernandes [2023] I. Mandal and R. M. Fernandes, Valley-polarized nematic order in twisted moiré systems: In-plane orbital magnetism and crossover from non-Fermi liquid to Fermi liquid, Phys. Rev. B 107, 125142 (2023).
- Abrikosov [1974] A. A. Abrikosov, Calculation of critical indices for zero-gap semiconductors, Journal of Experimental and Theoretical Physics 39, 709 (1974).
- Moon et al. [2013] E.-G. Moon, C. Xu, Y. B. Kim, and L. Balents, Non-Fermi-liquid and topological states with strong spin-orbit coupling, Phys. Rev. Lett. 111, 206401 (2013).
- Nandkishore and Parameswaran [2017] R. M. Nandkishore and S. A. Parameswaran, Disorder-driven destruction of a non-Fermi liquid semimetal studied by renormalization group analysis, Phys. Rev. B 95, 205106 (2017).
- Mandal and Nandkishore [2018] I. Mandal and R. M. Nandkishore, Interplay of Coulomb interactions and disorder in three-dimensional quadratic band crossings without time-reversal symmetry and with unequal masses for conduction and valence bands, Phys. Rev. B 97, 125121 (2018).
- Mandal [2018] I. Mandal, Fate of superconductivity in three-dimensional disordered Luttinger semimetals, Annals of Physics 392, 179 (2018).
- Mandal and Freire [2021] I. Mandal and H. Freire, Transport in the non-Fermi liquid phase of isotropic Luttinger semimetals, Phys. Rev. B 103, 195116 (2021).
- Freire and Mandal [2021] H. Freire and I. Mandal, Thermoelectric and thermal properties of the weakly disordered non-Fermi liquid phase of Luttinger semimetals, Physics Letters A 407, 127470 (2021).
- Mandal and Freire [2022] I. Mandal and H. Freire, Raman response and shear viscosity in the non-Fermi liquid phase of Luttinger semimetals, Journal of Physics: Condensed Matter 34, 275604 (2022).
- Roy et al. [2018] B. Roy, M. P. Kennett, K. Yang, and V. Juričić, From birefringent electrons to a marginal or non-Fermi liquid of relativistic spin- fermions: An emergent superuniversality, Phys. Rev. Lett. 121, 157602 (2018).
- Mandal [2021] I. Mandal, Robust marginal Fermi liquid in birefringent semimetals, Physics Letters A 418, 127707 (2021).
- Mandal and Freire [2024] I. Mandal and H. Freire, Transport properties in non-Fermi liquid phases of nodal-point semimetals, arXiv e-prints (2024), arXiv:2404.08635 [cond-mat.str-el] .
- Schulz [1990] H. J. Schulz, Incommensurate antiferromagnetism in the two-dimensional Hubbard model, Phys. Rev. Lett. 64, 1445 (1990).
- Igoshev et al. [2010] P. A. Igoshev, M. A. Timirgazin, A. A. Katanin, A. K. Arzhnikov, and V. Y. Irkhin, Incommensurate magnetic order and phase separation in the two-dimensional Hubbard model with nearest- and next-nearest-neighbor hopping, Phys. Rev. B 81, 094407 (2010).
- Sachdev and La Placa [2013] S. Sachdev and R. La Placa, Bond order in two-dimensional metals with antiferromagnetic exchange interactions, Phys. Rev. Lett. 111, 027202 (2013).
- J.A. Wilson and Mahajan [1975] F. D. S. J.A. Wilson and S. Mahajan, Charge-density waves and superlattices in the metallic layered transition metal dichalcogenides, Advances in Physics 24, 117 (1975).
- Scholz et al. [1982] G. Scholz, O. Singh, R. Frindt, and A. Curzon, Charge density wave commensurability in 2H-TaS2 and AgxTaS2, Solid State Communications 44, 1455 (1982).
- Chen et al. [2018] P. Chen, W. W. Pai, Y.-H. Chan, V. Madhavan, M. Y. Chou, S.-K. Mo, A.-V. Fedorov, and T.-C. Chiang, Unique gap structure and symmetry of the charge density wave in single-layer VSe2, Phys. Rev. Lett. 121, 196402 (2018).
- Gweon et al. [1998] G.-H. Gweon, J. D. Denlinger, J. A. Clack, J. W. Allen, C. G. Olson, E. DiMasi, M. C. Aronson, B. Foran, and S. Lee, Direct observation of complete Fermi surface, imperfect nesting, and gap anisotropy in the high-temperature incommensurate charge-density-wave compound SmTe3, Phys. Rev. Lett. 81, 886 (1998).
- Fang et al. [2007] A. Fang, N. Ru, I. R. Fisher, and A. Kapitulnik, STM studies of TbTe3: Evidence for a fully incommensurate charge density wave, Phys. Rev. Lett. 99, 046401 (2007).
- Feng et al. [2015] Y. Feng, J. van Wezel, J. Wang, F. Flicker, D. M. Silevitch, P. B. Littlewood, and T. F. Rosenbaum, Itinerant density wave instabilities at classical and quantum critical points, Nature Physics 11, 865 (2015).
- Holder and Metzner [2014] T. Holder and W. Metzner, Non-Fermi-liquid behavior at the onset of incommensurate charge- or spin-density wave order in two dimensions, Phys. Rev. B 90, 161106 (2014).
- Sýkora et al. [2018] J. Sýkora, T. Holder, and W. Metzner, Fluctuation effects at the onset of the density wave order with one pair of hot spots in two-dimensional metals, Phys. Rev. B 97, 155159 (2018).
- Halbinger et al. [2019] J. Halbinger, D. Pimenov, and M. Punk, Incommensurate density wave quantum criticality in two-dimensional metals, Phys. Rev. B 99, 195102 (2019).
- Senthil and Shankar [2009] T. Senthil and R. Shankar, Fermi surfaces in general codimension and a new controlled nontrivial fixed point, Phys. Rev. Lett. 102, 046406 (2009).
- Metlitski et al. [2015] M. A. Metlitski, D. F. Mross, S. Sachdev, and T. Senthil, Cooper pairing in non-Fermi liquids, Phys. Rev. B 91, 115111 (2015).
- ’t Hooft [1973] G. ’t Hooft, Dimensional regularization and the renormalization group, Nuclear Physics B 61, 455 (1973).
- Weinberg [1973] S. Weinberg, New approach to the renormalization group, Phys. Rev. D 8, 3497 (1973).
- Srednicki [2007] M. Srednicki, Quantum field theory (Cambridge University Press, 2007).
- Rao and Piazza [2023] P. Rao and F. Piazza, Non-Fermi-liquid behavior from cavity electromagnetic vacuum fluctuations at the superradiant transition, Phys. Rev. Lett. 130, 083603 (2023).
- Schäfer et al. [2017] T. Schäfer, A. A. Katanin, K. Held, and A. Toschi, Interplay of correlations and kohn anomalies in three dimensions: Quantum criticality with a twist, Phys. Rev. Lett. 119, 046402 (2017).
- Sýkora and Metzner [2021] J. Sýkora and W. Metzner, Fluctuation effects at the onset of density wave order with two pairs of hot spots in two-dimensional metals, Phys. Rev. B 104, 125123 (2021).