Mechanism for fluctuating pair density wave

<jats:title>Abstract</jats:title><jats:p>In weakly coupled BCS superconductors, only electrons within a tiny energy window around the Fermi energy,<jats:italic>E</jats:italic><jats:sub><jats:italic>F</jats:italic></jats:sub>, form Cooper pairs. This may not be the case in strong coupling superconductors such as cuprates, FeSe, SrTiO<jats:sub>3</jats:sub>or cold atom condensates where the pairing scale,<jats:italic>E</jats:italic><jats:sub><jats:italic>B</jats:italic></jats:sub>, becomes comparable or even larger than<jats:italic>E</jats:italic><jats:sub><jats:italic>F</jats:italic></jats:sub>. In cuprates, for example, a plausible candidate for the pseudogap state at low doping is a fluctuating pair density wave, but no microscopic model has yet been found which supports such a state. In this work, we write an analytically solvable model to examine pairing phases in the strongly coupled regime and in the presence of anisotropic interactions. Already for moderate coupling we find an unusual finite temperature phase, below an instability temperature<jats:italic>T</jats:italic><jats:sub><jats:italic>i</jats:italic></jats:sub>, where local pair correlations have non-zero center-of-mass momentum but lack long-range order. At low temperature, this fluctuating pair density wave can condense either to a uniform<jats:italic>d</jats:italic>-wave superconductor or the widely postulated pair-density wave phase depending on the interaction strength. Our minimal model offers a unified framework to understand the emergence of both fluctuating and long range pair density waves in realistic systems.</jats:p>

section_content: In weakly coupled BCS superconductors, only electrons within a tiny energy window around the Fermi energy, E F, form Cooper pairs. This may not be the case in strong coupling superconductors such as cuprates, FeSe, SrTiO 3 or cold atom condensates where the pairing scale, E B, becomes comparable or even larger than E F. In cuprates, for example, a plausible candidate for the pseudogap state at low doping is a fluctuating pair density wave, but no microscopic model has yet been found which supports such a state. In this work, we write an analytically solvable model to examine pairing phases in the strongly coupled regime and in the presence of anisotropic interactions. Already for moderate coupling we find an unusual finite temperature phase, below an instability temperature T i, where local pair correlations have nonzero center-of-mass momentum but lack long-range order. At low temperature, this fluctuating pair density wave can condense either to a uniform dwave superconductor or the widely postulated pair-density wave phase depending on the interaction strength. Our minimal model offers a unified framework to understand the emergence of both fluctuating and long range pair density waves in realistic systems. 

Spatially uniform superconducting (SC) order formed from Cooper pairs with zero center-of-mass momentum is the energetically favored ground state in the conventional theory of Bardeen, Cooper and Schrieffer (BCS) 1. Equivalently, the SC instability is signaled by a divergence in the static pair-fluctuation propagator, L(q, Ω = 0), at q = 0 once the pair instability temperature, T i, is achieved 2. On the other hand, a non-uniform order with non-zero center-of-mass momentum Cooper pair can occur when the divergence of the pairfluctuation propagator is shifted to non-zero q. First proposed by Fulde and Farrell (FF) 3 and independently by Larkin and Ovchinnikov (LO) 4, these solutions are stabilized in the presence of explicit timereversal symmetry breaking from an external magnetic field. A modulated order parameter can also be realized in the presence of time-reversal symmetry where the spatial average of the gap vanishes. Termed pair-density waves (PDWs), these states are posited to exist in a variety of systems, including high-temperature cuprate superconductors (for a review, see ref. 5

section_content: While PDWs have been subject to much theoretical [6] [7] [8] [9] [10] [11] [12] [13] [14] [15] [16] [17] [18] [19] and numerical [20] [21] [22] [23] [24] [25] [26] interest, a clear-cut analytically solvable model describing their origin from microscopic ingredients is lacking. From the experimental point of view, the interest for modulated pairing phases has been triggered by increasing experimental evidence for short-ranged PDW order in the underdoped region of the phase diagram of cuprates [26] [27] [28] [29] [30] [31] [32] [33] [34] [35] [36] [37] [38]. In particular, 32 reported the first clear observation via scanning tunneling spectroscopy of a vortex-induced PDW in Bi 2 Sr 2 CaCu 2 O 8 at low temperature. More recent STM experiments provide further evidence in favor of a short-range PDW coexisting with the d-wave superconductivity in the SC phase and evolving into a PDW state in the pseudogap region 26, 38. This phase is characterized by a gap at finite temperatures but lacks long-range order, and can be characterized as a "fluctuating pair density wave", locally pinned by disorder. Such a state also provides an explanation for many other experimental signatures of the cuprates, including the existence of vestigial charge density wave order arising from partial melting of a PDW 5, 15, 39. However, there is currently no microscopic model supporting this picture. Hence it is urgent to seek a unified framework that subsumes both fluctuating and long-range ordered PDW phases under a single paradigm by providing a concrete description of their origin. 

In this work, we show that a Fermi liquid subjected to a finite anisotropic interaction is unstable toward a modulated SC phase in the strong coupling limit. Whether this phase is a 'fluctuating' PDW (FPDW) or long-range order PDW is determined by temperature as well as the coupling strength defined by the ratio α = E B /E F, with E F the Fermi energy and E B the bound state energy for pair formation. 

Our strategy is to solve the self-consistent gap equation for the homogeneous d-wave superconductor and analyze the momentum dependence of the SC fluctuations. The expansion of the static pair propagator L q in powers of momentum transfer q, can reveal, in fact, critical fluctuations of Cooper pairs with finite center-of-mass momentum, that makes the homogeneous solution unstable toward a modulated SC phase. Since we study the action only to second order in the pair field we cannot distinguish between an instability to a FF or LO state. We observe the emergence of a modulated SC state already at intermediate coupling α ~0. 7. The appearance of such a state is linked to the existence of fluctuating terms that lower the momentum rigidity of the Cooper pairs. These terms directly follow from the anisotropy of the pairing interaction that affects the momentum dependence of the pairing susceptibility already in the normal phase. 

Our results are summarized in the phase diagram, Fig. 1. T i is the instability temperature of the homogeneous d-wave state obtained within the mean-field approximation. The analysis of fluctuations allows us to define two different regimes. At weak coupling, α ≪ 1, the uniform d-wave paired state is the ground state; at larger α (strong coupling), SC fluctuations at finite momentum lead to two modulated pairing phases-the T = 0 PDW ground state and a higher temperature FPDW phase that condenses into a PDW ordered phase below a coherence temperature (T c ). The strong and weak coupling regions are separated by the line T = T *. This is the temperature at which the finite momentum fluctuations around the homogeneous solution become critical. As expected in the BCS limit, the instability temperature T i and the coherence temperature T c coincide at weak coupling. In the strong coupling regime, we anticipate that T i and T c decouple since Cooper pairs are formed but with no long range coherence. In this work, we do not perform here any calculation of the coherence temperature T c inside the modulated phase. However, in analogy with results obtained for homogeneous s-wave superconductors in the strong coupling limit 40, we expect that T c < T i for α > 1 as well. The FPDW is found for temperature T c < T < T i and it is characterized by pairs with finite momentum with no coherence. At T = 0, the ground state can be either the uniform d-wave solution or the long-range PDW depending on the value of α. Hence our model captures two key experimentally postulated modulated Cooper phases -a FPDW and a long-range PDW-in a single unified scheme. 

The mechanism we present in this paper predicts spatially modulated pairing phases for α = E B /E F > 1, i. e., in strongly coupled electronic systems, with anisotropic interactions. Examples of lowdensity electronic materials include the Fe-based superconductor FeSe where quantum oscillations 41 as well as transport and scanning tunneling spectroscopy 42 show that both the electron and hole pockets are tiny with Fermi energies comparable or even smaller than the SC gap and for which we find several proposals of BCS-BEC cross-over physics in the literature [43] [44] [45]. Other "mixed-band" superconductors such as O vacancy-or Nb-doped SrTiO 3 have one partially filled band with a large Fermi surface while the Fermi level intersects the other at or close to the band bottom 46. Even if these materials typically have more than one band close to or crossing the Fermi level, the results from our minimal model may eventually provide a suitable startingpoint for the analysis of possible instabilities toward modulated pairing states in dilute multiorbital superconductors. Our results may also be relevant to the recent observation of superconductivity in twistedbilayer graphene 47 where interactions can be large compared to the bandwidth leading to large inter-particle distances 48 and hence possible strongly-coupled Cooper pairing. 

The modulated phases we propose in this work, that include both the long-range ordered PDW as well as the FPDW at finite temperature, are distinct from earlier proposals in literature. Loder et al. 10, considered similar models characterized by nearest neighbor attractive interaction with d-wave symmetry and found Cooper pairing with finite center-of-mass momentum above a critical interaction strength. In refs. 18, 19, a modulated superconducting state is found in models which have correlated pair-hopping interactions. Other models that admit modulated SC ground states were proposed in the context of cold atoms 9 where local interactions were considered in systems with multiple bands. Those references focused on the analysis of the longrange ordered state (mainly at zero temperature) without exploring the FPDW phase. The key contribution of our work is it provides an analytically tractable model where both fluctuating and long-range ordered PDWs can be explained under a single unified framework.
section_name: and references therein).

section_name: Results

section_content: Let's consider a single band SC system. The kinetic part of the Hamiltonian reads H 0 = ∑ kσ ξ k c kσ c kσ, where ξ k = ϵ k -μ, μ is the chemical potential, ϵ k = k 2 /2m the parabolic dispersion and we further assume 2m = 1. The pairing interaction is given by

g is the constant SC coupling and θ q is defined as

where f k,q = (h k-q/2 + h k+q/2 )/2 is a form factor. In this work h k can be any anisotropic form factor; we consider, e. g.,

y Þ=Λ with dwave form. The pairing energy scale is Λ i. e., the high energy cut-off set by the inverse lattice spacing and much larger than E F. Our results do not depend qualitatively on the exact form of the anisotropy, provided it is strong enough, but they are distinct from the conventional s-wave case f k,q = 1. Note that the interaction Hamiltonian we choose above already assumes an attractive pairing interaction and does not begin from any repulsive Hubbard-type model. Nevertheless, the mechanism we present is 'microscopic' in the sense that the order parameter can be evaluated over the entire phase diagram in terms of the microscopic parameters such as Cooper binding energy and Fermi energy. Such a treatment is distinct from Ginzburg-Landau theory and consistent with terminology used in previous literature 49, 50. 

We use the standard Hubbard-Stratonovich transformation to decouple the interaction term, Eq. ( 1 ) and to derive the effective action in term of the bosonic pairing field Δ (for a detailed derivation see Supplementary Note 1).
section_name: Model

section_content: In standard BCS superconductors, the mean-field value of the pairing field is defined by minimizing the action with respect to the homogeneous q = 0 value of Δ and then solving this equation together with the one for the chemical potential. To study fluctuations of the pairing field around the mean-field value, we analyze the gaussian action obtained by retaining up to the second order in the fluctuating field with arbitrary momentum q given by S

where L À1 q is the static pairing susceptibility. As discussed in 51 and in the Supplementary Notes 1-3 of this work, we need to consider in principle both the real and imaginary part of the fluctuating field. However, here we are interested in the analysis of the static limit of the fluctuation, for which real and imaginary part are decoupled. Hereafter we will analyze the gaussian fluctuation of the amplitude mode only to investigate the possible emergence of a spatially modulated pairing fluctuation out of a homogeneous d-wave SC state. The explicit form of the static pairing susceptibility for the amplitude mode is L À1 q = g À1 + Π q, where the particle-particle propagator reads

Here T is the temperature and V the volume. Note that since 2m = 1, energies have dimensions of 2-D V -1, and L À1 q is therefore dimensionless. 

The static susceptibility can be expanded in the hydrodynamic limit as

The instability temperature is defined as the highest temperature at which the susceptibility diverges, i. e., c 0 = g À1 + Π 0 | T = T i = 0, as we assume that the minimum of the action, Eq. ( 3 ), is associated with the homogeneous order parameter. The coefficient c 2 = ð∂ 2 L À1 q =∂q 2 | q = 0 Þ=2 provides instead information about the momentum rigidity of the fluctuating Cooper pairs i. e., the energy needed to move the center-ofmass momentum of the Cooper pairs from zero to a finite value. A negative momentum rigidity, c 2 < 0, implies that finite momentum fluctuations can lower the energy of the system making the homogeneous SC solution unstable. This means that the highest temperature at which the pairing susceptibility, Eq. ( 5 ), diverges is actually associated to a critical mode with finite momentum. 

In what follows we analyze the momentum-dependence of the static susceptibility, Eq. ( 5 ), looking for a sign change of the momentum rigidity parameter c 2 and using it as a proxy to identify possible spatially modulated SC regions in the phase diagram. It is worth noticing that c 2 is directly affected by the momentum properties of the pairing susceptibility i. e., the pairing symmetry. From Eq. ( 4 ), it is easy to verify that the anisotropy of the interactions affects the momentum dependence of the propagator not only in the SC phase via the symmetry of the SC order parameter, but also above the instability temperature T i where Δ = 0 due to the overall form factor f 2 k,q at the numerator. This reflects in a strong momentum dependence of the contributions to the rigidity parameter depending on the symmetry of the pairing interaction. We discuss below how this affects the development of critical finite-momentum fluctuations. 

The mean-field analysis for the homogeneous d-wave superconductor is shown in Fig. 2. In panels (a)-(b) we report the selfconsistent numerical mean-field solutions for the pairing function Δ and the chemical potential μ as a function of temperature T for three representative cases of the pairing strength α = E B /E F = 0. 5, 1. 0, 2. 0, where for simplicity the weak-coupling expression E B = Λe -2/g is used at all α. In panels (c)-(d) we show the same mean-field results at T = T i and T = 0 as a function of α. The change of sign of the chemical potential with increasing coupling strength is well-known from the BCS-BEC crossover problem [52] [53] [54] [55] [56]. In the weak-coupling regime, the pairs are loosely bound and we recover the BCS expression μ ~EF. As the interaction increases, all fermions strongly bind in pairs and μ becomes negative and proportional to -E B. In both the weak and strong coupling limits, the curves are similar to those derived for s-wave superconductors in 56, showing that the d-wave symmetry of the pairing interaction does not affect the mean-field results qualitatively. 

We first study the SC fluctuations above the instability temperature by analyze the static pairing susceptibility in the hydrodynamic limit, Eq. ( 5 ). The mass term c 0 is positive and vanishes as the temperature approaches the instability temperature as expected from a Ginzburg-Landau description of the transition. 

The analysis of the momentum rigidity of the fluctuating pairs above T i is shown in Fig. 3. The weak coupling region is characterized by a standard regime of fluctuations with c 2 > 0. Here Cooper pairs with zero center-of-mass momentum are stable. Increasing α, the momentum rigidity for the d-wave pairing interaction (continuous line) monotonically decreases and becomes negative at intermediate coupling, α > 0. 7, as shown in Fig. 3(a). This means that finite momentum critical fluctuations grow, increasing the coupling strength up to a critical value of the interaction for which the homogeneous SC solution can become unstable toward a modulated phase. Notice that c 2 becomes very small and eventually changes sign in the crossover between weak and strong coupling where also the chemical potential changes sign from positive to negative, see inset Fig. 3(a). The result changes qualitatively for the isotropic s-wave interaction (dashed line) where the rigidity parameter decreases but remains positive even at strong coupling for the set of model parameter of our study (This result differs from the analysis of 55, in which a sign change of the rigidity parameter is found for a SC system with s-wave pairing symmetry at strong coupling. In this case, however, the authors first perform a strong coupling expansion of the pairing susceptibility and only subsequently expand the approximated result in powers of q. ). 

To characterize the modulated SC state and check its stability, we expand the static susceptibility to higher order in momentum

We report the coefficients of the momentum expansion at T i in Fig. 3(b). Results are shown as a function of α for the coupling regime in which c 2 ≲ 0. We need to expand the susceptibility up n = 6 to find c 6 > 0, since for our set of model parameters c 4 < 0 as in the conventional BCS case. 

We analyze the momentum dependence of the static susceptibility at T i in Fig. 4, where we show the expansion of Eq. ( 6 ) up to sixth order for different values of α. At the instability temperature, c 0 = 0 by definition and the minimum of the function is determined by the higher order coefficients. At weak coupling, where c 2 is large and positive, the minimum of L À1

q is located at zero momentum. As the pairing interaction increases c 2 becomes small and eventually changes sign at α ~0. 7. Here, since c 4 < 0, the minimum shifts discontinuously to a finite momentum Q, i. e., by increasing the interactions the modulated phase emerges at T i via a first order transition from the homogeneous d-wave SC solution, in analogy with the results found at T = 0 in 10, 18. The non-zero value of Q at α ~0. 7 signals the formation of the FPDW state with finite momentum pairing but no long range coherent order. Note that the finite order parameter jump Q is a nonuniversal quantity and depends on microscopic details of the chosen model, as is a feature of any generic first order transition. The momentum characterizing the modulated phase shifts toward larger values increasing the coupling parameters. In the strong coupling regime, α ≫ 1, the minimum occurs at q= ffiffiffi ffi Λ p ≫ 1, (not shown), for this range of the interaction the analysis of the momentum characterizing the modulated phase requires the implementation of a non perturbative approach. 

To further prove the stability and quantitatively characterize the modulated phases we expand the polynomial form, Eq. ( 6 ), around its minima Q so that we can define the susceptibility of the modulated SC phase

Here c 0 and c 2 are the mass and the momentum rigidity parameter associated with the modulated SC phase. In the Supplementary Notes 2, 3, 4, we show that c 0 = aðT À T i Þ with T i = T i + δT and, δT > 0, and that c 2 > 0, thus demonstrating stability of the modulated phases. 

The sign change of the momentum rigidity parameter discussed at T = T i can be traced down in temperature (dashed line in Fig. 1 ). At T = 0 the homogeneous d-wave state becomes unstable, now toward a PDW, for a slightly higher value of the coupling where the chemical potential μ also changes sign (see Fig. 2d ). The stability of the PDW phase requires expanding up to the sixth-order, c 4 < 0, c 6 > 0 as we show in the Supplementary Notes 2-4. 

The results of our numerical study are summarized in the phase diagram of Fig. 1. We characterized the SC region below T i by the sign of the momentum rigidity parameter (dashed line). The sign change of the c 2 coefficient at strong coupling signals the presence of critical SC fluctuations at finite momentum that make the d-wave homogeneous state unstable toward either an FPDW or PDW. The pink and purple regions indicate the FPDW and the long-range ordered PDW state at high and low temperatures respectively. We leave for future work the explicit calculation of the coherence temperature below which the FPDW condenses. The color gradient indicates approximately the expected T c (α) behavior based on previous analysis of the coherence energy scale for the homogeneous s-wave SC state 40. 

Analytical calculations of the momentum rigidity can be easily performed within a simplified model in which the chemical potential is used as parameter. Both at T i and T = 0, we find qualitatively the same results discussed within the numerical study. In particular, within the analytical calculations sketched in the Supplementary Note 5, the momentum rigidity parameter follows the chemical potential behavior, i. e., c 2 (μ) < 0 for μ < 0. This relation is qualitatively in agreement with the numerical study performed computing self-consistently μ(α), as one can see from the inset of Fig. 3a. 

The strategy implemented here to investigate how finite momentum fluctuations become critical at strong coupling is based on the analysis of the momentum rigidity parameter. This method presents two main advantages with respect to other theoretical approaches. On the one hand, as already discussed, it allows us to explore the finite temperature regime and analyze the FPDW state. On the other, it provides a physical understanding of the importance of the anisotropy of the pairing interactions in the development of the modulated phase. As one can see in Eq. ( 4 ), the symmetry of the pairing interactions dramatically affects the momentum dependence of the propagator not only in the SC phase, but also in the normal one when Δ = 0 due to the overall form factor f 2 k,q. This is reflected in a strong momentum dependence of the contribution to the momentum rigidity parameter. In fact, after performing analytically the Matsubara summation, the computation of the c 2 coefficient reduces to an integral over the Brillouin zone c 2 = 1 V P k I 2 ðkÞ. The expression for I 2 is given in the Supplementary Note 3, but here we show here in Fig. 5 2D maps of I 2 (k) for both s-wave and d-wave at T = 0 and T = T i. In the isotropic s-wave case, the contributions to the momentum rigidity coming from different momenta, I 2 (k), are positive at any (k x, k y ). Whereas, in the d-wave case the contributions to the momentum rigidity coming from the nodal regions are negative and dominate the overall sign of the c 2 coefficient.
section_name: Pair susceptibility

section_content: A consistent explanation for the occurrence of both static and fluctuating Cooper pairs with finite momentum in the phase diagram of materials such as cuprates has been a long-standing problem. This is primarily because an identification of the microscopic ingredients driving such exotic pairing has been elusive. The results in this paper point toward a simple and unified framework that naturally promotes both fluctuating and static pair-density wave (FPDW and PDW) phases over their zero momentum counterparts. Figure 1 summarizes the main conclusions of our work, supported not only by numerical evaluations but also transparent analytical estimates (see Supplementary Note, where E B is the pair binding energy for two electrons on the Fermi surface in the presence of an attractive interaction, and E F is the Fermi energy. For the specific set of parameters presented here, below a critical value of α ~0. 7, only uniform zero momentum d-wave pairing is favored. In the approximate range of 0. 7 ≲ α ≲ 1. 5, the FPDW phase, characterized by a negative momentum rigidity c 2 and positive c 6 (see Fig. 3 ), is stable over a range of temperatures below the instability temperature T i. However, in this range of α a uniform d-wave pair is still favored at zero temperature. For α ≳ 1. 5, the PDW phase is more stable than a uniform solution at T = 0 and a finite momentum pair exists for all temperatures below T i. The modulation wave vector Q of the paired phases is determined by the ratio α and acquires a jump with increasing α as in a first order transition (see Fig. 4 ). 

It is important to note that while in our paper the critical value of α for which c 2 changes sign appears to be of order unity, which is nominally outside the range of standard BCS weak-coupling theory, we emphasize that this value is a finite non-zero number that does not take a universal value. Depending on parameters, we can easily produce critical values of α of order 0. 3, which might be considered weakcoupling (see Supplementary Fig. 2 ). Hence, there appears to be no physical reason that necessarily constrains FPDW and PDW phases to be in the strong coupling regime of α ~1. 

Furthermore, we note the key role of the existence of a lattice momentum cut-off. If the cut-off were taken to infinity, the modulated Fig. 4 | Momentum dependence of the sixth order expansion of L -1 at T i (dimensionless units). At weak coupling, α = 0. 22, we find the homogeneous dwave SC. The momentum rigidity c 2 is large and positive, and the minimum of the inverse of the susceptibility is at q = 0. At intermediate coupling, α ~0. 7, c 2 vanishes and the minimum of L -1 appears at a finite Q of order 1. Same set of parameters of Fig. 2. solution vanishes as is outlined in the Supplementary Note 5. This is further highlighted in the recent exact two-body and variational wave function solution 57 where it was shown that in the absence of a cut-off in the interaction, the modulated solution on a lattice loses to the homegeneous solution. 

Recently refs. 58, 59 discussed the possible existence of a diagonal pair density wave order in the cuprates. To obtain a diagonal pair density wave order in our analysis, the interaction must contain a dominant B 2g pair-fluctuation interaction. This would imply that the form factor of the kind sin k x sin k y would be dominant in the fluctuations instead of the cos k x À cos k y form factor in Eq. (1). With this form factor, our result continues to hold but with the dominant instability now occurring along the diagonals. Note, however, that our theory does not justify which of the two (B 1g or B 2g ) symmetry channels are the dominant fluctuations. In either of the two cases, one can obtain the finite momentum pairing instability without requiring a strong coupling treatment of the theory. 

The FPDW and PDW phases are stabilized by contributions to the fluctuation free energy arising from momenta close to the nodal regions in the Brillouin zone. These contributions, which also should drive strong anisotropy in the phase stiffness near T i, are suppressed (enhanced) at weak (strong) coupling thus leading to a modulated phase above a critical pairing strength. This simplified picture is confirmed from our numerical calculations (Fig. 5 ). Finally, while our work primarily focuses on the instability temperature T i in the strong coupling limit, the behavior of the condensation temperature T c and the fluctuations around the PDW ground state in this setting are open problems that will require further investigations. Our work does not consider the competing effects of a nematic superconducting phase that has been phenomenologically found to suppress the PDW at T > 0 in 2D 7, 11. In addition, even if allowed by our model, we have not addressed the possible coexistence at low T of a PDW and a homogeneous d-wave superconductor, as suggested by cuprate experiments 5, 26. Our results as such set the stage for future microscopic descriptions of modulated superconductivity in strongly coupled materials. 

Note added: During review of the manuscript, another work based on repulsive interactions yielding PDWs became available 60.
section_name: Discussion

section_content: We used standard many-body field theoretic methods for all computations. Details of these methods are included in the Supplementary Information provided.
section_name: Methods

section_name: Acknowledgements

section_content: We thank P. Abbamonte, B. M. Andersen, S. Caprara, A. Chubukov, E. Fradkin, S. A. Kivelson, M. Granath and A. Toschi for useful discussions and suggestions. C. S., L. F., and P. J. H. are supported by the DOE grant number DE-FG02-05ER46236. L. F. acknowledges support by the European Union 's Horizon 2020 research and innovation programme through the Marie Skłodowska-Curie grant SuperCoop (Grant No 838526 ).
section_name: Acknowledgements

section_name: Data availability

section_content: All data generated or analyzed during this study are included in this published paper (and its Supplementary Information files)
section_name: Data availability

section_content: All codes used to generate or analyze the results of this study are available from the corresponding author (C. S. ) on reasonable request.
section_name: Code availability

section_name: Author contributions

section_content: C. S. and L. F. contributed equally to this work. C. S., L. F., and P. J. H discussed and contributed to the project design. C. S. performed the analytical calculation, L. F. performed the numerical analysis. C. S., L. F., and P. J. H. contributed to the data analysis, to the interpretation of the theoretical results, and to the writing of the text.
section_name: Author contributions

section_content: The authors declare no competing interests.
section_name: Competing interests

section_name: Additional information

section_content: The online version contains supplementary material available at https://doi. org/10. 1038/s41467-023-38956-x. 

Correspondence and requests for materials should be addressed to Chandan Setty, Laura Fanfarillo or P. J. Hirschfeld. 

Peer review information Nature Communications thanks the anonymous reviewers for their contribution to the peer review of this work. A peer review file is available. 

Reprints and permissions information is available at http://www. nature. com/reprints Publisher's note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations. 

Open Access This article is licensed under a Creative Commons Attribution 4. 0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons. org/ licenses/by/4. 0/.
section_name: Supplementary information