Entropic Origin of Pseudogap Physics and a Mott-Slater Transition in Cuprates

We propose a new approach to understand the origin of the pseudogap in the cuprates, in terms of bosonic entropy. The near-simultaneous softening of a large number of different q-bosons yields an extended range of short-range order, wherein the growth of magnetic correlations with decreasing temperature T is anomalously slow. These entropic effects cause the spectral weight associated with the Van Hove singularity (VHS) to shift rapidly and nearly linearly toward half filling at higher T, consistent with a picture of the VHS driving the pseudogap transition at a temperature ~T*. As a byproduct, we develop an order-parameter classification scheme that predicts supertransitions between families of order parameters. As one example, we find that by tuning the hopping parameters, it is possible to drive the cuprates across a transition between Mott and Slater physics, where a spin-frustrated state emerges at the crossover.

Evidence is growing that the 'pseudogap phase' found in cuprates is home to one or more competing phases, including a variety of stripe, spin-, or charge-density wave (S/CDW) phases [1][2][3][4][5][6][7][8][9][10][11][12] . The CDW phase, in particular, has stimulated considerable interest [13][14][15][16][17][18][19][20][21][22] . Many of these phases, including superconductivity, seem to appear at temperatures well below the pseudogap temperature T * , so the exact relation between the pseudogap and these other phases remains elusive. Indeed, the real puzzle is understanding why the pseudogap phase bears so little resemblance to a conventional phase transition. Here we demonstrate that the pseudogap phenomenon arises from strong mode coupling, where a large number of density waves with similar ordering wave-vectors q attempt to soften and condense at the same time. The associated entropy leads to anomalously low transition temperatures and extended ranges of short-range order, characteristic of the pseudogap phase. In an extreme case, this can lead to a condensation bottleneck where long-range order is absent, leaving behind a novel entangled spin-frustrated phase.
To understand how entropy drives short range order, consider the Hubbard model with on-site repulsion U, nearest neighbor hopping energy t, and exchange J = 4t 2 /U. A simple mean field solution correctly captures an antiferromagnetic ground state with ordering vector q = (π, π) and gap 2Δ ~ U. However, the same calculation yields a Neel temperature T N ~ U, instead of the strong coupling result, ~1/U. The resolution of this problem is the electronic entropy: repeating the calculation for magnetic order with any other q-vector, including ferromagnetism (q = 0) also produces a gap of order U. While the (π, π) state is the ground state, the other states are higher in energy by no more than a factor ~J. Hence, at low T the energy wins out, and the system condenses into its ground state, while at higher temperature entropy wins and the system lowers its free energy by maximizing its entropy, forming a mixture of all these competing phases. Thus, longe-range order is lost at a T N ~ J 23 . The problem is how to incorporate this entropic effect into a calculation of the phase diagram. Similar problems have arisen in the past. In Overhauser's theory of CDWs in alkali metals 24 , the q-vectors for all points on a spherical Fermi surface become unstable simultaneously and cannot be handled one at a time. These effects are often referred to as phonon entropy 25 , or more generally as boson entropy, and can be analyzed via vertex corrections which take proper account of mode coupling 26,27 . Here we develop a similar theory for the cuprates via a self-consistent renormalization calculation of the vertex corrections. Similar effects are likely to be present in many other families of correlated materials.
The ultimate origin of the bosonic entropy is localization. We provide two examples of this. First, in the above discussion of AFM, large U means single spins tend to localize on individual copper atoms and the effective hopping is renormalized from t-to J-scales, leading to much flatter dispersions, and to small differences in energy Scientific RepoRts | 7:44008 | DOI: 10.1038/srep44008 phonon from condensing until very low T and giving rise to an extended range of short-range order. We find that similar strong entropic effects arise from electronic bosons (electron-hole pairs), predominantly associated with the VHS. In turn, the entropic effects modify VHS physics, causing the main spectral weight associated with the VHS to shift rapidly towards half filling as T increases, not following the evolution of the density of states (DOS). This anomalous shift suggests that the VHS may play a large role in pseudogap physics, Fig. 1 9 . The key anomaly in pseudogap physics lies not below the pseudogap temperature but above T * , Fig. 1(b), where the resistivity in the strange metal phase varies linearly with T, ρ ~ T. The pseudogap temperature T * indicates the onset of deviations from linear resistivity, while below a lower temperature T ** the resistivity becomes purely quadratic, suggestive of a Fermi liquid regime. While the strange metal phase is suggestive of the quantum critical phase in the Heisenberg quantum phase transition 33 , there is no obvious reason why the resistivity of the Heisenberg model should be linear in T, nor any natural candidate for either the pseudogap (since the Heisenberg model has a well-defined ordered phase), or for the transition at T ** . In contrast, these features arise naturally in our VHS model, Fig. 1(a). For T > T VHS Pauli blocking is completely ineffective, leaving the expected VHS scattering rate τ −1 ~ T. Below T VHS this blocking of the VHS starts to turn on, while coherent Fermi surface features appear, until below a coherent-incoherent crossover, T < T coh , the VHS scattering is suppressed, leaving behind a Fermi liquid τ −1 ~ T 2 . These results are consistent with experimental findings that the pseudogaps in Bi2201 and Bi2212 terminate near the VHS doping x VHS [34][35][36] . Mode coupling and order parameter competition. The bosonic entropy can be incorporated into the DFT-Lindhard function by modifying Moriya's mode-coupling approach 37 to include realistic band-structures via a susceptibility DOS [Results Section (Beyond RPA)]. We use this technique to calculate the growth of the correlation length in LSCO, finding that the Mermin-Wagner theorem is satisfied leading to a broad regime of short-range order. We also find a temperature range where the correlation length does not grow at all. By tuning the hopping parameters, we find that this is associated with a 'supertransition' between two families of order parameters. This is particularly interesting as it involves competition between (π, π) Mott antiferromagnetic order driven by VHS fluctuations and a (π, π − δ) Slater SDW driven by conventional Fermi surface nesting. We are able to define an order parameter DOS that diverges at the transition: Fig. 2 shows the susceptibility near the transition for two different values of t′ /t. At the transition, the susceptibility at (π, π) changes from a maximum to a local minimum, and the susceptibility is quite flat across the entire (π, π)-plateau. This is especially true in frame (b), which is tuned close to the T = 0 transition (note change of scale). It is convenient to introduce a susceptibility density-of-states (DOS) to quantify the number of competing states. This behavior can be approximately captured in a simple 'Mexican hat' model, Fig. 2(c). A similar analysis for Bi2201 finds a transition between (π, π − δ) and antinodal nesting order, consistent with the SDW-to-CDW transition found in most cuprates.

Methods
Our calculation is a form of many-body perturbation theory (MBPT) based on Hedin's scheme 31 . The scheme involves four elements: electrons are described by Green's functions G with DFT-based dispersions renormalized by a self-energy Σ ; electronic bosons [electron-hole pairs] are described by a spectral weight [susceptibility] χ , with A 4 > 0 and varying A 2 : A 2 > 0 (red line), χ = 0 (blue line), and < 0 (green line). renormalized by vertex corrections Γ . Neglecting vertex corrections, the self-energy can be calculated as a convolution of G and W = U 2 χ (GW approximation). This approach has been used to solve the energy gap problem in semiconductors, where the Γ correction leads to excitons via the solution of a Bethe-Salpeter equation 38,39 , and in extending DMFT calculations to incorporate more correlations [e.g., DMFT + GW, etc.]. Our approach here is to extend our previous GW calculations [quasiparticle-GW or QPGW 40 ] to approximately include vertex corrections.
In QPGW we introduce an auxiliary function is the bare, or DFT dispersion. G Z behaves like the Green's function of a Landau-type QP -a free electron with renormalized parameters that describes the low-energy dressed electronic excitations. However, this is a non-Fermi liquid type QP, since the frequency-integral of Im(G Z ) is Z and not unity. That is, the Z-QP describes only the coherent part of the electronic dispersion, and is not in a 1:1 correspondence with the original electrons. The importance of such a correction can be readily demonstrated. Since a Z-QP has only the weight Z of a regular electron, the susceptibility [a convolution of two Gs] is weaker by a factor of Z 2 than an ordinary bare susceptibility. To match this effect in the Stoner criterion requires introducing an effective U eff = ZU. In contrast, MBPT calculations in semiconductors typically set the GW-corrected Green's function to 41 , thereby missing the reduced spectral weight of the low-energy, coherent part of the band.
In this paper the DFT-Lindhard susceptibility χ 0 is calculated at intermediate coupling level, using Z-corrected first-principles dispersions to describe the coherent part of the quasiparticle-GW dressed carriers 40 . We calculate χ 0 based on the dressed dispersion  Z = Z 0 , assuming a doping-independent Z = 0.5. We use a single-band model of the cuprates, working in a purely magnetic sector, with the Hubbard U controlling all fluctuations; there is a competition in this model between near-nodal and antinodal nesting (ANN) which mimics the SDW-CDW competition in cuprates, sharing the same nesting vectors 22,32 .

Results: DFT-Lindhard susceptibility
Sorting the susceptibility peaks: basins of attraction for cuprates. Fluctuation maps. We define fluctuation maps in x-T space as maps of domains ('basins of attraction') of classes of q-vectors at which χ 0 (q, 0) is maximum. Figure 3 displays fluctation maps for La 2−x Sr x CuO 4 (LSCO) (a) and Bi 2 Sr 2 CuO 6+x (Bi2201) (b), using the full DFT dispersion. We fix an energy scale by taking t = 0.4195 eV, appropriate for LSCO 42 . Most of the colored regions are associated with FS nesting, while the white region is mainly associated with (π, π) VHS nesting, as will be discussed below. This figure confirms our earlier finding that LSCO, Fig. 3(a), has a very different phase diagram from most other cuprates, here represented by Bi2201, Fig. 3(b) 32 . In the phase labeled 'ring' , the susceptibility is approximately constant along a ring of q-vectors surrounding (π, π). We note in particular that antinodal nesting (red-shaded region) emerges naturally in most cuprates at higher doping, but is absent in LSCO. Shown also in Fig. 3 are T VHS (dark red-dashed line) and the shift of the DOS-peak (thin green-dashed line), and a series of mean-field phase boundaries, based on a Stoner criterion, There is a critical value of U = U 0 below which there is no phase transition at x = 0. For U < U 0 , the phase boundary is dome-shaped, but with a peak shifted away from the T = 0 value of x VHS to a value close to the line x VHS (T) (red-dashed line). For U > U 0 the phase boundary rapidly evolves to a strong coupling form where the highest T transition occurs near half-filling. Consistent with this, we will show that VHS phenomena are dominated by entropy effects, and thus scale with the (π, π) susceptibility and not with the peak DOS. In all cases, at sufficiently high T, FS nesting is washed out and the dominant instability is at q = (π, π). In comparing these phase boundaries to experiments in cuprates, it should be kept in mind that the dotted curves refer to constant U, whereas U is likely to be strongly screened near a VHS 43 .
Reference families. The above fluctuation maps allow us to introduce equivalence classes of models for cuprates, and in particular to determine simple reference families that can be used to search parameter space beyond the physical cuprates. To understand differences between various cuprates, and to gauge their proximity to the pure Hubbard model, it is necessary to add a third axis to the usual T-and doping-phase diagrams. Here we find reference families for which differences between physical cuprates' band structures can be parametrized by a single hopping parameter. By a Wannier-type downfolding of DFT results, each cuprate can be modelled by a single CuO 2 band crossing the Fermi level, so their parameter space consists of hopping parameters, traditionally t, t′ , t′ ′ , … . In fact, since t sets an energy scale, the band dispersion is determined by the ratios t′ /t, t′ ′ /t, … To reduce the number of hopping parameters, we introduce a notion of equivalence, wherein a set of materials are equivalent if they have the same fluctuation map. From this equivalence, the DFT parameters can be mapped onto an equivalent reference family depending only on two parameters, namely t′ /t and t′ ′ /t. We find that when the physical cuprates are approximated by their reference families they lie close to a particular cut in this 2D space. . We find that reference family (ii) semi-quantitatively reproduces the DFT derived fluctuation maps for both cuprates. On the other hand, setting t′ ′ = 0 works well for LSCO, Fig. 4(a), but fails for Bi2201 Fig. 4(b). We further find that a reasonable choice of t′ /t is the one which matches the Van Hove singularity (VHS) doping of the full DFT dispersion. We note that both reference families have the additional desirable feature that they evolve from the state with t′ = 0− i.e., the original Hubbard model.
An advantage of using the reference families is that this lets us systematically vary parameter space to study the evolution from one cuprate to another, or to extrapolate the parameter space to better understand the origin of various phases. Figure 5(a) shows the cuts in t′ − t′ ′ -space associated with reference families (i) (t′ ′ = 0, red line) and (ii) (t′ ′ = − t′ /2, blue line). The dotted lines, Fig. 5(a), indicate ground-state transitions at x = 0, from (π, π) to (π, π − δ) (black dots) and from (π, π − δ) to (δ, δ) order (red dots). The black dotted line is the Mott-Slater line discussed in the Results Section (Beyond RPA). The densities-of-states (DOS) shown for the two families in Fig. 5(b,c) provide insight into the origin of the crossovers. Note that the use of these reference cuts demonstrates the clear evolution of the VHS between 2D and 1D behavior.
For reference family (i), Fig. 5(c), the DOS reveals a crossover associated with a change from 2D to 1D nesting. This is seen by the DOS at t′ = − t/2 which shows the characteristic E 1/ divergence expected for a 1D system. To explore the origin of the 1D-behavior we consider the electronic dispersion for reference family (i) of the form E = − 2t(cos(k x a) + cos(k y a)) − 4t′ cos(k x a) cos(k y a). If either k x or k y = 0 and t′ = − t/2 then E = − 2t. Therefore the FS at the VHS energy reduces to just the x and y axes. This shows that the strong 1D nature of the corresponding FS leads to a strong divergence in the DOS at the VHS.
Near t′ = − t/2, the DOS has two peaks, one at the saddle-point VHS found in 2D materials, the other at the leading edge, as found in 1D systems, Fig. 5(c) (violet line). For t′ < − 0.45t (green-black), the leading edge is larger, suggestive of weakly-coupled chains, while for t′ > − 0.28t (green dotted arrow) the leading edge peak turns into a local minimum, Fig. 5(c) inset, at the green dotted arrow. This quasi-1D behavior is suggestive of the nematic phase found in cuprates 45 . Indeed, an earlier VHS renormalization group calculation found that charge order transitions splitting the degeneracy of the VHSs at (0, π) and (π, 0) [which would now be called a nematic transition] are only possible if t′ /t < − 0.276 43 . Remarkably, t′ = − 0.28t approximately coincides with the commensurate-incommensurate transition found in Fig. 2(b) to lead to a severe ξ depression.
For family (ii) a non-zero t′ ′ modifies the dispersion such that it breaks up the 1 dimensionality of the FS and weakens the DOS divergence, as seen in Fig. 5(b). However, a residual 1 dimensionality remains in that there is conventional FS nesting in the ANN region, which is responsible for the CDW. For the red, orange, yellow, and green curves we find a resemblance to that of the DOS for t′ ′ = 0. At the light blue curve there is a topological transition where a pocket centered at (π, 0) splits off from the rest of the Fermi surface. For larger t′ , one now sees two features in the DOS, a peak followed by a step at lower energies. The peak is the topological transition where the pockets split off, and the step is where the (π, 0) pocket disappears. The inset illustrates this progression showing the FS at the DOS peak energy for various t′ .
Origins of structure in the susceptibility: competition between FS and VHS nesting. Separating FS maps and broad background. In order to understand the structure in the fluctuation maps we need to analyze the susceptibilities more deeply. The susceptibility contains sharp features which are a map of the Fermi surface and which quantify the strength of the Fermi surface nesting. In addition, the susceptibility contains a smooth, intense background feature which peaks at (π, π). This peak plays an important role in the cuprates, producing large entropic effects. This peak can shift the balance of the FS nesting to q-vectors closer to the peaks, and in special cases, it can lead to commensurate nesting away from the FS nesting vector at exactly (π, π). Moreover, as T increases, coherent FS features are washed out, leaving behind only the commensurate bulk contribution.
Here we determine that this background can be considered to be a bosonic VHS (b-VHS), in that the electronic excitation comes predominantly from the region of the VHS. It is analogous to the Van Hove excitons found in optical spectra 46 , but present in the intraband susceptibility. The evolution of this feature contains hints as to the origin of the pseudogap.
Consider the expression for the ω = 0 bare susceptibility that enters into the Stoner criterion, where Δ f k,q is the Pauli blocking factor  Fig. 3, that x 0 is the doping at T = 0, but we assume  x x 0 .
Scientific RepoRts | 7:44008 | DOI: 10.1038/srep44008 The denominator is smallest if  k =  k+q , but in this case the numerator vanishes at T = 0 unless  k is at the Fermi level. Even in this case, the term contributes little weight to the susceptibility integral unless the q-shifted FSs are tangential, which is generally satisfied if q = 2k F , giving rise to a folded map of the FS -the ridge of FS nesting 32 .
The Fermi surface nesting features in the susceptibility map for Bi 2 Sr 2 CuO 6 (Bi2201), at x = 0.13 and ω = 5 meV, are illustrated in Fig. 6, with χ ′ 0 given in frame (a), and the corresponding map of χ ″ 0 in frame (b). Due to the energy δ-function, χ ″ 0 contains only the FS contribution 47 , except an extra feature at Γ associated with the DOS. In Fig. 6(d), we plot χ ′ 0 (red line) and χ ″ 0 (blue line) along a cut Γ → (π, 0) → (π, π) → Γ in momentum space. It can be seen that, away from Γ , the peaks in both components fall at the same q-values, and the relative weights of the peaks are comparable. Moreover, the anti-nodal nesting peak near Γ in χ ″ 0 is stronger than the near-(π, π) peak. However, χ ′ 0 is largest on the near-(π, π) peaks, since the FS contribution is riding on the large background contribution to the susceptibility. Here we demonstrate that VHS nesting plays a dominant role in creating this background, thereby providing a deeper understanding of the VHS vs FS nesting competition.
Two approaches are taken to eliminate the FS nesting features to see the background susceptibility by itself. First, raising T (violet long-dashed line in Fig. 6(d)) washes out the sharp FS features first, while leaving the background features behind. Secondly, in frame (c), we plot the difference between χ ′ 0 and a suitably scaled quantity proportional to χ ″ 0 ; this is also plotted as the green curve in frame (d). Since χ ′ 0 and χ ″ 0 are Kramers-Krönig transforms, they have different lineshapes, but since the FS map is nearly one-dimensional, a large part of it is still cancelled in taking the difference. From frame (d) it can be seen that both approaches lead to similar results for the bulk susceptibility, and that this contribution provides the dominant contribution to virtually the whole (π, π) plateau structure. Here the background contribution arises as follows. Even when  k ≠  k+q , there will be a finite contribution to the susceptibility, Eq. 2, reduced by the factor 1/( k −  k+q ), as long as the two energies lie on opposite sides of the Fermi level (for T = 0). This contribution will be especially significant near a VHS, which has a large spectral weight.
Background peak as VHS nesting. Here we demonstrate the connection between the bulk contribution to χ 0 and VHS nesting by deconvolving the bare susceptibility, Eq. 2, into its various  k components. To understand this procedure, recall that the VHS gives rise to susceptibility features near Γ (DOS) and (π, π). The latter are typically Scientific RepoRts | 7:44008 | DOI: 10.1038/srep44008 more intense. Thus, to confirm that the background susceptibility is associated with the VHS, we analyze Eq. 2 for χ 0 , assuming that the electron is associated with  k and the hole with  k+q . We must demonstrate that, when q ~ (π, π), the dominant contribution to χ 0 has k near the VHS at (π, 0) [or (0, π)] and  k near E VHS .
To demonstrate this we divide the Brillouin zone into four patches, around the Γ , (π, 0), (0, π), and (π, π) points, and calculate separately the susceptibility due to each patch, rewriting Eq. 2 as where i runs over the four quadrants Γ , (π, 0), (0, π), and (π, π). In LSCO, the plateau is predominantly composed of equal contributions from the (0, π) and (π, 0) patches, as expected for a VHS. We then explore the distribution of values of  k which contribute to the susceptibility in the (π, 0) patch at x = 0, Fig. 7. There is a clear crossover from a dominant near-FS contribution ( k = 0) at low T to a dominant contribution at high-T which lies along a flat (q-independent) line centered at energy k B T VHS = E F − E VHS . As we vary doping, this flat line shifts towards the Fermi level following T VHS , red-dashed line in Figs 4(a) and 8. This crossover is associated with Pauli blocking, i.e. the factor Δ f in Eq. 5. Since at low T the Fermi function is a step, Δ f blocks the scattering between states far below E F . As T increases f gets smeared out, so that features away from the Fermi surface can contribute to the susceptibility. Since the VHS has a substantial peak, the high-T changes are due to the unveiling of this feature. Note that at low T, Fig. 7(a-c), there is still scattering from one VHS, but to states near E F , and hence contributing intensity at q ≠ (π, π). Correspondingly, the weight near (π, π) comes from states in the tail of the VHS extending to E F . (c) Difference plot of (a) and the scaled version of (b). (d) Cuts of susceptibility from earlier frames plotted along high symmetry lines, with red curve representing data in frame (a), blue curve for frame (b), and green curve for frame (c). Shown also is the curve corresponding to frame (a), but at T = 1000 K (violet long-dashed line) scaled by 1.4 and shifted by − 0.045 eV −1 to match other curves near (π, π).
Scientific RepoRts | 7:44008 | DOI: 10.1038/srep44008 While Δ f enhances features associated with the VHS as T increases, it has the opposite effect on features near the Fermi level. For energy  exactly at E F , f() is exactly 1/2 at any T, so Δ f = 0 if both states are at E F . As T increases, f smooths out, so that if | k |, | k+q | are both < k B T, then f k ≈ f k+q ≈ 1/2, and therefore, Δ f ≈ 0. This explains why the crossover from FS nesting to VHS nesting in Fig. 7 occurs in such a narrow T range (pink shaded region in Fig. 8). Passing between Fig. 7(c,d) the Fermi surface contribution goes through a coherent -incoherent crossover -note that the FS nesting contribution is fully incoherent by 500 K.
In Fig. 8 we superpose this crossover, plotted as a pink shaded region, on the fluctuation map for LSCO, with all calculations for the minimal reference family (t′ ′ = 0). This incoherent-coherent crossover scales with T VHS , with boundary lines given by T VHS /2 and T VHS /3. The same crossover can be seen directly in the susceptibility,   Fig. 3, that x 0 is the doping at T = 0, but we assume  x x 0 .
Scientific RepoRts | 7:44008 | DOI: 10.1038/srep44008 Fig. 9, as the coherent, FS-related features present at low T are washed out, leaving only a featureless (incoherent) (π, π) plateau at high T. These effects of Pauli blocking can explain several additional features of Fig. 8. For instance, hole doping causes the (π, π)-plateau to expand in momentum space, leading to a weakening of VHS nesting. Thus, in Bi2201, the transition between (π, π)-nesting and some form of FS nesting shifts to higher T with increasing x, as seen in Fig. 4(b,d,f). A similar effect is found in LSCO, Fig. 4(a,c,e), for the high-T termination of the incommensurate (π, π − δ) order (green shaded region). However, in LSCO there is also an anomalous reentrant transition, where (π, π) order reappears below a lower domain boundary (blue solid line). To understand this unusual reentrance, consider first what happens to χ 0 when E F = E VHS , Fig. 9(a). While the VHS peak is clearly present at T = 0 (blue line), finite temperature does not lead to simple thermal broadening but actually splits the peak. For small finite T the very sharp VHS at E F is blocked. As T increases, the blocking spreads over a wider range of energy, leading to the splitting of the VHS feature. Thus, raising T causes a blocking first of the sharp VHS peak to produce the lower edge of the domain boundary. Shifting to x < x VHS means that a higher T is necessary to block the sharp VHS peak. Note that the lower branch of the boundary closely follows T VHS , especially for the t′ ′ = 0 reference family. This is also why the conventional VHS theory that relies on the sharp VHS peak only holds at relatively low T as discussed in conjunction with Fig. 3.
Pauli-blocking and the bosonic-VHS. If we think of χ 0 as the propagator of different electron-hole pairs, then the pair associated with q = Q ≡ (π, π) has many similarities with the excitonic-VHS seen in optical studies 46 , and we will refer to it as the bosonic-VHS, or b-VHS for short. We clarify the role of Pauli blocking by describing its effect on the b-VHS in detail, using the t′ ′ = 0 reference family for simplicity. To understand the relationship between the b-VHS and entropy effects, we need to know how Pauli blocking "hides" the states associated with the VHS at low T, and how they recover as T increases. If the renormalized dispersion of a single electron is  k with wave vector k, then an e-h pair at wave vector q has a dispersion ω q (k) . For LSCO, the dominant pairs are those at q = Q, the AF nesting vector. The associated dispersion ω Q (π, π), plotted in Fig. 10(a), resembles the electronic dispersion  k , but with an important distinction: it depends only on  Q− (k), whereas all of the hopping terms that shift the electronic VHS away from half-filling (t′ , t′ ′ ) are contained in  Q+ (k), i.e. the b-VHS is pinned at ω = 0. Since χ ″ 0 is an odd function of ω, χ ω ″ = = Q ( , 0) 0 0 . However, while the b-VHS is pinned at ω = 0, Fig. 10(a), its weight vanishes at T = 0, due to the Pauli-blocking factor, Δ f k,Q = 0 near k = (π, 0) at T = 0, Fig. 10(b). Finite T restores weight, optimally near 1000 K, although Δ f always vanishes exactly at (π, 0). Hence the bosonic entropy only turns on at finite T. Figure 7(g) illustrates how the Pauli window opens with increasing T. As T increases, the Fermi functions spread out, so the FS contribution decreases whereas the bulk contribution grows as Δ f becomes unblocked. For the VHS, we expect full unblocking at a temperature near T ~ T VHS (red short-dashed line in Fig. 8). Since f varies exponentially with T, the unblocking will appear at a temperature T < T VHS , which may vary according to the property studied. This is similar to a Schottky anomaly, where a feature at energy Δ above the ground state produces a peak in the heat capacity at a temperature T m ~ Δ /2.4. The low-T fadeout of the b-VHS is controlled by the Pauli blocking factor, Entropy and a revised VHS scenario. Here we will demonstrate that this b-VHS peak produces large entropic effects. To show this, we calculate several characteristic VHS-dominated properties -the Sommerfeld coefficient of the heat capacity and the susceptibilities at Γ and (π, π). The entropy is calculated from the standard expression Fig. 11. The sum is over the electronic dispersion of LSCO, assuming a paramagnetic phase. At T = 0, γ is proportional to the DOS, and hence diverges at x VHS , Fig. 11. However, as T increases, the peak does not remain near x VHS , but rapidly shifts toward x = 0. We plot the peak position T γ both in the inset of Fig. 11 and as a violet dot-dashed line in Fig. 8. Clearly, since γ is a measure of entropy, it is sensitive to sources of entropy away from the Fermi level, when T is high enough to reduce the Pauli blocking. It has recently been noted that there should be an excess of low-T entropy in the vicinity of a quantum phase transition 48 . Similar behavior is found for χ 0 . The VHS has logarithmically diverging features in the susceptibility at two different q-vectors, an inter-VHS contribution at (π, π), and an intra-VHS contribution at Γ 49,50 . While both peaks cross the Fermi level at x VHS at T = 0, they evolve differently with temperature. This is because the feature at Γ represents the peak in the DOS, and hence it is dominated by near-FS physics. As a result it shifts only weakly with T, green long-dashed line in Fig. 8. In contrast, the inter-VHS (π, π) peak is a bulk contribution and hence is sensitive to Pauli blocking, as discussed above. This leads to the rapid thermal evolution of the (π, π)-susceptibility peak all the way to x = 0 (light blue dot-dashed line), scaling with both T γ and T VHS , Fig. 8. However, this behavior seems to be dominated by a separate aspect of the VHS physics. As discussed in ref. 32, the susceptibility has an approximate electron-hole symmetry, with similar spectra for a nearly empty (x ~ 1) or a nearly full (x ~ − 1) band. The t′ and t′ ′ hopping parameters break this symmetry, and as a result, the low-T susceptibility is roughly symmetrical about x VHS , where it has its largest peak. Here it seems that as T increases, the electron-hole symmetry is restored when T > t′ . Since T VHS ~ t′ , one expects a similar scaling of T vs x for x ≥ 0. But when x < 0 T VHS continues to increase, whereas the peak in χ 0 (q = (π, π)) stays at x = 0, indicative of a restored electron-hole symmetry (see Supplementary Material Fig. 2).
These entropic effects have profound consequences for VHS physics that were overlooked in earlier mean-field studies 49,50 . Most earlier studies of cuprate antiferromagnetism used a model valid only near the FS, or at low-T, overlooking Pauli unblocking of the VHS feature at large T. As a consequence, they predicted that χ 0 would show a maximum for x near x VHS , inconsistent with experiment. Figure 3(a) shows revised RPA phase diagrams for several values of U, based on the DFT-Lindhard susceptibility. For small enough U (U < ~ 1.3 eV), the earlier results are recovered, but when U is larger a new behavior is found, with a maximum transition temperature at half-filling, decreasing with doping as if approaching a quantum critical point. This now starts to resemble the pseudogap phase diagram of the cuprates. The comparison of experimental pseudogap lines 9,51 in Fig. 8 will be discussed in the Discussion Section. Finally, when McMillan 25 introduced the idea of bosonic entropy, he meant that competition between different phonon instabilities (phonons with different q-vectors) slowed down the tendency to CDW ordering, leading to broad T-ranges with only short-range order. In the following section we will show that a similar competition exists for electronic bosons, and that mode coupling allows us to quantify these effects.

Results: Beyond RPA susceptibility
Mode coupling formalism. Above, we have established the outlines of a classification scheme for phase transitions, analogous to the spectrum-generating algebras (SGAs) of nuclear physics. We summarize this development in Supplementary Material Section I.A, and compare it to SGAs in Supplementary Material Section I.C. In the present section we show how competition between condensing phases affects the transition temperature, and how to quantify the number of competing modes. A few striking examples of this competition will be discussed.
From χ 0 we calculate an interacting χ of modified random-phase approximation (RPA) form, where λ is a mode-coupling parameter which must be calculated self-consistently. The imaginary part of χ can be thought of as the density-of-states (DOS) of electronic bosons, electron-hole (e-h) pairs, which may become excitons or excitonic resonances when a long-range Coulomb interaction is turned on. However, static (ω = 0) instabilities depend on a Stoner criterion, and hence on The self-consistent parameter λ is found from a Matsubara sum of the susceptibility i n q 0 , n where N is the number of q-points, and the summation in Eq. 9 can be transformed: The term λ 1 introduces a small, nonsingular correction 23 to Eq. 9, which we neglect. In order to explore the role of suseptibility plateaus for realistic band dispersions, the self-consistency equation is evaluated numerically. Note that since χ − 0 1 has dimensions of energy, a plot of χ − q ( , 0) 0 1 resembles a dispersion map. Hence we can define a SDOS: , is the variable of integration, and we introduce a corresponding susceptibility density of states (SDOS) N − . Eq. 9 then becomes U c 0 c with δ = U c − U sp , U c = 1/χ 0 (Q 0 , 0). In Section II of Supplementary Materials, we show that Eq. 14 is closely related to excitonic Bose condensation, with λ proportional to the effective number of bosons. Once N − has been calculated, Eq. 14 can be evaluated numerically. The singular part of the integral is treated analytically, and the remainder numerically, with u and U approximated as constants, u = 0.8eV −1 , U = 2 eV. It is convenient to fix the SDOS at some temperature T′ , and then solve Eq. 14 for T(ξ th , T′ ), with self-consistency requiring T(ξ th , T′ ) = T′ . Figure 12 illustrates how integrating over the inverse susceptibility, Fig. 12(a), leads to the SDOS, Fig. 12(b), at x = 0 where the near-(π, π) plateau is dominant. Here we use hopping parameters appropriate to the DFT dispersion of Bi2201, t = 419.5, t′ = − 108.2 and t′ ′ = 54.1 meV (Subsection Mode coupling in Bi2201). This yields a doping phase diagram that is qualitatively similar to that of most cuprates, except LSCO. By comparing Fig. 12(a,b), one can see how features in χ − 0 1 translate into features in N − . Thus, the intense, flat-topped peak in N − at small values of χ − 0 1 represents the near-(π, π) plateau. Its broad leading edge (smaller χ − 0 1 ) is controlled by anisotropy of the plateau edge between (π, π − δ) and (π − δ, π − δ), while the sharp trailing edge corresponds to the local maximum of χ − 0 1 at (π, π). For the ANN peak, its leading edge scarcely leaves any feature in N − , but a local maximum translates into a large peak in N − . Finite temperature, T = 100K (blue line), rounds off the cusp in Fig. 12(c), leading to a step in N − , Fig. 12(d), but otherwise has little effect. Instead of such a self-consistent calculation, many treatments of phase transitions, including path-integral approaches, spin-fermion models, and Hertz-Millis 52,53 quantum critical theories, approximate the susceptibility denominator by a T-independent Ornstein-Zernicke(OZ) 54 form, in terms of various deviations from the critical point (in q, ω, and a 'tuning parameter' which is proportional to ξ −2 , where ξ is the correlation length). Here z is a dynamic exponent and = − q q Q 0 is the difference between the actual wave number q and Q 0 , the wave number where the susceptibility has a peak, If the OZ parameters are assumed to be T-and doping-independent 23,37,55,56 , the resulting physics becomes quite simple. For 2D materials, the Mermin-Wagner (MW) theorem 57 is satisfied, and the RPA transition at T mf turns into a pseudogap onset at T* ~ T mf , with a crossover to long-range order when interlayer coupling is strong enough -in short, not much changes from the mean-field results. When the OZ form of χ is assumed, 2 , then N − becomes a constant, which we denote by N a , and and N a = a 2 /4πB. Eq. 14 leads to long range order at T = 0 only, with correlation length ξ π δ = a N / 1/ 4 th a given by 23 where T 2 = πA 2 λ/6Γ ua 2 , A i is the coefficient of q i , and q c is a wave number cutoff. The present G Z − W Z − Γ Z model is the simplest model which captures the essential physics of the pseudogap. For a more quantitative comparison with experiment, two additional problems must be solved. First, an extension to fully self-consistent G − W − Γ may be needed to see how the susceptibility responds to the gap-opening. Secondly, the term in Eq. 15 proportional to ω z must be included to describe quantum fluctuations. Since we are primarily interested in the opening of the pseudogap at higher temperatures, we ignored it in our analysis.
Scientific RepoRts | 7:44008 | DOI: 10.1038/srep44008 Strong mode coupling leads to an extended range of short-range order. Thus, mode coupling modifies the Stoner criterion to where λ satisfies Eq. 14. In the OZ approximation, this becomes This OZ approximation misses the essential and strongly T-dependent mode-coupling physics, and we develop an alternative formulation of the theory. In particular, the inverse curvature of the susceptibility 1/A 2 can diverge due to a competition between conventional Fermi-surface nesting and Van Hove singularity (VHS) nesting, and this parameter plays a special role in the theory. In particular, we can get a semiquantitative picture of the correlation length changes by simply accounting for the T dependence of A 2 (T). We define the order-parameter DOS (OPDOS) as the threshold SDOS, OP c c 2 2 2 and use this in the OZ formalism, Eqs 16-18. While this neglects the fact that q c can also have a strong T-dependence, and the contribution away from threshold is not negligible, the simplification should work when ξq c ≫ 1. We believe that this OPDOS quantifies the number of competing modes and provides an essential link to the SGA theory. The connection between Eq. 14 and entropy is further discussed in Supplementary Material Section III. Figure 13 illustrates the profound effects that strong mode coupling has in LSCO, as well as the complete inability of the OZ approximation to capture this physics. The SDOS, Fig. 13(a), contains VHS-like features characteristic of conventional DOSs. However, the singular behavior of Eq. 20 involves only features near the threshold, −X U c , which evolve strongly with T, see inset to Fig. 13(a). For T > 0, the threshold behavior is always a step at X − = U c , indicative of a parabolic peak in χ 0 with curvature inversely proportional to the step height, Eq. 21. For analyzing critical phenomena, we define a threshold correlation length ξ th (red solid lines in Fig. 13(b)) in terms of the initial curvature of χ 0 at Q 0 , Eq. 17. We then evaluate ξ th using the full Eq. 14, but the results can be interpreted in terms of an OZ form with a strongly T-dependent step height N OP (T). However, we caution that the experimentally measured ξ 1/2 is typically taken from the susceptibility half-width. This ξ 1/2 is generally larger than ξ th due to the fast fall-off of χ near the edge of the (π, π)-plateau. Thus, for the example of Fig. 13(c), we calculated a ξ 1/2 (filled blue circles) from the full half-width of the renormalized χ, and found that typically ξ ξ 2 th 1/2 , where the solid red line in Fig. 13(c) plots ξ th vs 1/T; the values of ξ 1/2 are in good agreement with experiment (green dot-dot-dashed line) 58 . We note that while the OZ form is used in quantum critical theory 35,36 , it is explicitly stated that it is to be used only in a limited T-range, and only in the absence of FS nesting 36 , both of which are violated here.  Figure 13(b) illustrates how strong mode coupling slows down the correlation length divergence. Undoped LSCO (red solid line) shows two regions I and III of exponential growth of ξ th with decreasing T, separated by an anomalous region II where ξ th actually decreases with decreasing T. While the model captures the MW-like divergence at low T (region III), our main interest is in the higher-T behavior. In the high-T limit (region I), the leading-edge parabolic curvature is quite small, and if it were T-independent, as in the OZ approximation, the growth in ξ th would follow the blue dashed line (Eq. 18), but thermal broadening causes the curvature to decrease with increasing T, leading to the faster growth of the red solid line. The origin of the anomalous region II for undoped LSCO, Fig. 13(b), can be readily understood from Fig. 13(a), where for T < 1500 K there is excess SDOS weight near U c , leading to a strong peak (inset) as T → 0. This feature represents the development of the (π, π)-plateau. In Supplementary Material Section III we present a simple calculation of ξ th , showing how a range of power-law behavior, ξ~T 1/ th 1/2 (dotted line in Fig. 13(b)), could arise. While this power-law behavior captures the average T dependence of the anomaly, it cannot reproduce the nonmonotonic variation of the correlation length.
Exploring parameter space. The reference families can also be used to interpolate between physical cuprates. Here we use interpolation to understand the origin of the condensation bottleneck found in region II. We find that it is caused by proximity to a much stronger transition, a crossover between Mott and Slater physics, and that LSCO is more correlated than other cuprates, lying close to this crossover. Remarkably, this conclusion echoes an early neutron scattering study, which claimed that LSCO is close to an instability that is off of the physical parameter plane 59 .
We study two important cuts in t′ /t − t′ ′ /t space, a minimal cut (t′ ′ = 0) and the Pavarini-Andersen [PA] cut (t′ ′ = − t′ /2) -the latter seems to best capture the physics of the cuprates. By tuning t′ we unveil the origin of the condensation bottleneck as a localization-delocalization crossover tied to the crossover from (π, π)to FS-nesting. As a byproduct, we gain insight into why LSCO is so different from other cuprates, and how cuprates evolve from the pure Hubbard limit (t′ = 0). Figure 14(a) shows the T-evolution of ξ th for several values of t′ along the minimal cut, including the data of Fig. 13(b). For t′ /t > − 0.17, the system is characterized by commensurate (π, π) order with an effective Neel temperature T N ~ 1000 K [the correlation length grows so rapidly that interlayer correlations will drive a transition to full 3D order]. Similarly, for t′ /t < − 0.345 there is incommensurate (π, π − δ) order with T N about a factor of 10 smaller. But for intermediate t′ /t the commensurate-incommensurate transition is highly anomalous, with correlation length orders of magnitude smaller than expected. Figure 14(b) shows that a similar evolution follows along the PA cut in parameter space. The reason for this anomalous behavior can be seen by looking at the leading edge SDOS in the crossover regime, shown for three values of t′ along the minimal cut in Fig. 14(c-e). For ease in viewing, these curves have been shifted to line up the SDOS at (π, π) at all T. It is seen that the anomalous collapse of ξ th is associated with a rapid growth of the step height N OP ~ 1/A 2 , culminating in a near-divergence at ′ = − . t t 0 27 c , where the leading edge curvature A 2 goes to zero, Fig. 2(b). Note that ξ th drops by 9 orders of magnitude at T = 200 K when t′ /t changes from − 0.17 to − 0.27, then grows by a similar amount at 100 K when t′ /t changes from − 0.27 to − 0.345. This transition may be hard to detect experimentally, since over this same range, the ξ 1/2 will be frozen at the value corresponding to the half-width of the (π, π) plateau. The green shaded region in Fig. 14(c) shows that the range of the anomalous growth (II) of ξ 1/2 in Fig. 13(b) coincides with the range of rapid growth of the SDOS leading edge; such behavior is absent if a T-independent OZ form (black dashed line) is assumed.
From the relationship between ξ th and the step height A 2 , Eq. 18, we see that a diverging A 2 will cause ξ th → 0. The divergence arises when the susceptibility at (π, π) crosses over from a maximum to a local minimum. In the latter case, the maximum intensity of χ(q, 0) is spread along a 'ring' in q-space surrounding (π, π), so that the inverse susceptibility resembles a 'Mexican hat' , Fig. 2 Fig. 15(a). Note that the SDOS-divergence resembles a 1D VHS in the conventional DOS, even though here the 1-D direction is the radial direction away from (π, π). The divergence of χ along a ring can be thought of as a 2D analog of Overhauser's effect 24 . Similar effects are found in Bi2201, next Subsection.
The strength of the anomaly in ξ th at ′ = − . t t 0 27 c suggests that it is not an ordinary quantum critical point. Indeed, from Fig. 3 it can be seen to be a 'supertransition' between two domains of attraction, one commensurate at (π, π), the other incommensurate at (π.π − δ) (compare Fig. 3(a,b) at x = 0). Such a supertransition is known from spectrum-generating algebras (SGAs) to be highly anomalous. The diverging OPDOS at the supertransition, also characteristic of SGAs, indicates that many different q-vectors compete simultaneously, frustrating the divergence of any particular mode. This can be clearly seen in Fig. 2(b), where the susceptibility is nearly flat. This is the electronic analog of McMillan's phonon entropy: if many phonons are simultaneously excited, the transition is suppressed to anomalously low temperatures. At this transition the ground state is infinitely degenerate, representing an anomalous spin frustrated state arising in the absence of disorder, with the frustration arising from strong mode coupling, although whether it represents a spin liquid or spin glass requires further analysis. Note also that when the commensurate-incommensurate transition is at T > 300 K, ξ th (T) has a sharp downward cusp at the transition, while above the transition ξ th is strongly suppressed, ξ th < a, arrows in Fig. 14(a,b).
For LSCO, the anomaly in Fig. 13(b) (large red dot in Fig. 8(c)) seems to fall near the coherent-incoherent crossover (pink shaded region), consistent with an early neutron scattering result which indicated that doped LSCO is close to a magnetic QCP 59 . This proximity to a novel disorder-free spin-glass or spin-liquid QCP may play a role in stripe physics. Notably, in LSCO, commensurate (π, π) order disappears rapidly by ~2% doping, being replaced by incommensurate magnetic fluctuations and low-T spin-glass effects. Moreover, the commensurate-incommensurate transition involves a highly disordered regime separating two well-ordered phases ( Fig. 14(a,b)). Since ordering tends to lower the free energy of the electronic system, the disordered Scientific RepoRts | 7:44008 | DOI: 10.1038/srep44008 regime represents a state of high free energy, and doping across this region can lead to a regime of [nanoscale] phase separation (NPS), which in LSCO is manifest as the stripe phase 60 . It should be noted that also in Cr the commensurate-incommensurate AF transition is first order 61 . The very different situation in most other cuprates is discussed in the next Subsection, where it is shown that high-T magnetic fluctuations give way at low T to a fragile CDW order, consistent with experiment.
Finally, we note that the t′ of the Mott-Slater transition is also the point at which the DOS at the bottom of the band changes from a local minimum to a local maximum (green dotted arrow in the inset to Fig. 5(c) -see also Fig. 2), and close to the point at which a nematic transition [splitting the degeneracy of the VHSs at (0, π) and (π, 0)] becomes possible 43 . Mode coupling in Bi2201. We noted earlier that the OZ form should not be used in quantum critical theory when FS nesting is present 53 . The problem can however be addressed by using SDOS. The issue is that at T = 0 χ 0 is nonanalytic, with cusp-like peaks, Fig. 12, so the leading q-dependence is not parabolic. Indeed, since the dominant instabilities tend to be associated with double-nesting, the peak shape can be quite anisotropic, ~q in some directions and ~q 1/2 in others 32 . There is a simple relation between the q-dependence of the peak susceptibility and the form of N − (X − ) at threshold. If the SDOS varies as a power law, . For a parabolic onset , or a step in N − at threshold (p = 0 or N − = const > 0 for X − > U c ). Instead, for Bi2201 there are cusps in χ 0 near (π, π − δ) and an antinodal cusp at (δ, δ). For both cusps, a smooth power law is found, with p ~ 2, or 2 /3 at threshold. However, this leads to new problems. When a term q p for p ≤ 1 is substituted into Eq. 14, the integral converges at threshold, so the Mermin-Wagner effect appears to be absent. That this is not the case can be seen from Fig. 12: at any finite T, the cusp is rounded off, and χ 0 varies parabolically, which means A 2 is finite. The cusps still produce a strong effect in that A 2 has an exceptionally large T-dependence, A 2 → ∞ as T → 0. This leads to an extra strong divergence of ξ th , Eq. 18, at T → 0. But, as we increase T from zero, this also means that the corresponding ξ th decreases rapidly, leading to a very fragile phase. We have separately confirmed that the order is also sensitive to ω and disorder.
In Bi2201 and most other cuprates, a larger ratio |t′ /t| expands the (π, π)-plateau, leading to a rich phase diagram, Fig. 3(b). In this case the b-VHS is dominant only at high T, while FS-related susceptibility features cause a crossover of the (π, π) susceptibility from a maximum at high T to a local minimum at lower T. This leads to a different kind of strong mode coupling at intermediate T, best exemplified by the ring phase, Fig. 4(d,f). In this phase, the leading-edge SDOS has a 'Mexican hat'-like divergence (previous Subsection), Fig. 15(a) near 2000 K, leading to a slow growth in ξ − T th 1 , Fig. 15(c), quite similar to that observed in doped YBCO 62,63 . Near the onset of ring order, T = 1625 K, the correlation length has a nonmonotonic T-dependence, similar to that seen in LSCO, Fig. 14(a,b). For each curve, a partial SDOS near the ANN peak is superposed to show the SDOS associated with the ANN peak whose onset is marked by arrows. (c) Correlation length vs T. Blue dotted line represents ξ 1/2 ~ T −1 . Arrow shows dip in ξ 1/2 associated with commensurate-incommensurate transition. (d) Susceptibility of χ(q) peak vs T for q near the (π, π)-plateau (red line) or for ANN nesting (blue line).
Scientific RepoRts | 7:44008 | DOI: 10.1038/srep44008 While many pseudogap-like features are similar to LSCO, Bi2201 has some novel properties associated with the ANN instability, which can become a CDW when charge order is included in the model. Since the ANN susceptibility peak has a much larger intrinsic curvature, mode-coupling effects are weak. Hence the ANN susceptibility peak grows rapidly at low temperature, and as doping increases, it can surpass the (π, π)-peak, thereby sharply cutting off the growth of near-(π, π) fluctuations. This is illustrated in Fig. 15 for x = 0.20. Notably, since this peak is not associated with strong mode coupling it has a weak SDOS with a very strong T-dependence, Fig. 15(b,d), crossing the AFM peak near T = 100 K. The large subsequent growth of the ANN correlation length would ensure the rapid establishment of a true long-range 3D ANN order. However, corresponding to the large T-dependence of χ ANN , there is a similar strong disorder dependence -particularly in the presence of short-range near-(π, π) order, which should cause the ANN correlation length to saturate.

Discussion
Cuprate pseudogap. In Figs 1 and 8, and Results Section (DFT-Lindhard), we discussed a revised VHS scenario in which a number of VHS related features evolve nearly linearly in doping and scale with T VHS . We noted that this scaling resembles the pseudogap crossover T * (x). Here we provide additional evidence for this identification. Note that T VHS is typically close to the measured T * in most cuprates, including YBCO and Hg-cuprates 9 (yellow-green line in Fig. 8) and LSCO (magenta dotted line) 51 . Furthermore, the pseudogaps in Bi2201 and Bi2212 terminate near the point T VHS → 0, i.e., at the conventional VHS [34][35][36]64,65 .
The phase diagram of Fig. 8 has a number of further consequences for cuprate physics. First, the VHS onset near T VHS combined with the coherent -incoherent crossover (pink shaded region in Fig. 8) can explain the anomalous transport properties found near the pseudogap, Fig. 1. This parallels transport, where for T > T * , the resistivity ρ varies linearly with T 9 , a behavior expected near a VHS 66 . For lower T, the resistivity is mixed, but ρ ~ T 2 , as expected for a coherent Fermi liquid below a T coh < T VHS . The brown long-dashed line in Fig. 8 represents the corresponding T coh found in YBCO and Hg-cuprates 9 , in excellent agreement with the coherent-incoherent crossover.
We note further that the compressibility should have a peak similar to γ that extrapolates to the VHS doping as T → 0, consistent with our Fig. 8. In a related cluster-DFT calculation 67 , the pseudogap is associated with a Widom line, related to a peak in the compressibility. This is further discussed in Supplementary Material Section IV. Notably, the Widom line terminates in a first-order phase transition near the VHS. The similarity of this phase diagram to the excitonic liquid-gas transition found in semiconductors 68 suggests that the b-VHS may be best described as a VHS-exciton 69,70 .
This picture bears a resemblance to the Barzykin-Pines model of the cuprate pseudogap 71 , identifying T VHS and the coherent -incoherent crossover with T * and T coh < T*/3 in their model. Since their model is related to Kondo lattice physics 72 , the question arises as to whether a similar mode-coupling calculation in heavy-fermion compounds could lead to a similar anomalous entanglement at the f-electron incoherent-to-coherent transition.
Classifying order parameters. Competition between order parameters also arises in a totally different field, nuclear physics, where order parameters are formally classified in terms of a spectrum generating algebra (SGA) (Supplementary Material Section I.C). Order parameters are grouped into families (Lie subalgebras), and an order parameter density-of-states (OPDOS) is defined. In addition to conventional phase transitions, there are 'supertransitions' (called excited-state quantum critical points) where the system is driven across a boundary between two families and the OPDOS diverges. As a byproduct of this work, we develop a similar classification of order parameters for condensed matter systems, where similar phenomena can be observed (Supplementary Material Section I.A).
Strong coupling physics. The condensation bottleneck that we identified in Fig. 14 appears to be a previously unexpected transition between Mott and Slater physics, shedding light on the controversies between Anderson and Scalapino 73,74 . On the Slater-side of the line (larger |t′ |), the cuprate physics is dominated by FS nesting with Neel temperatures (strictly, ξã / 1000 th ) near 100 K, while on the Mott side, the Fermi surface plays no role, and the magnetic transition is at T ~ 1000 K, and always lies at (π, π) due to VHS nesting. Thus, the phenomenological strong-coupling Yang-Rice-Zhang model 75 can only work on the Mott side of the transition, since it requires (π, π) nesting. However, at the crossover, the correlation length remains nanoscale, leading to strong mixing of the two competing phases. Among the physical cuprates, all fall into the Slater domain, except for LSCO, which lies close to the Mott phase. Notably, LSCO shows the strongest signs of nanoscale phase separation in the form of stripe physics [Results Section (Beyond RVB)]. This result is consistent with a recent first-principles calculation of the Hubbard U, which finds LSCO more strongly correlated than other cuprates 76 . We recall that our self-energy formalism 40 is able to reproduce most spectral features of the cuprates in terms of competition between d-wave superconductivity, a (π, π) ordered AFM phase, and an ANN CDW. We reproduce not only the photoemission dispersions, limited to the lower Hubbard bands, but optical and x-ray spectra, which depend sensitively on the Mott gap. In a related 3-band model, we reproduced the Zhang-Rice result that the first doped holes are predominantly of oxygen character, with overall dispersions at half-filling agreeing with subsequent DMFT results 40 . While our earlier calculations predicted long-range (π, π) AFM order to persist to too high T, this is corrected in the present calculations which have only short range order and reproduce the upper and lower Hubbard band dispersions with broadening ~1/ξ th 23 . Moreover, the calculated ξ 1/2 values (filled blue dots in Fig. 13(c)) agree with experiment in LSCO (green dot-dot-dashed line) 58 . The experimental data are shown only above 300 K, since at lower T interlayer coupling drives a transition to long-range order. In the Heisenberg model, ξ 1/2 is a function only of T/J, and these data have been used to adduce the value of the exchange J. Hence, our mode coupling calculation successfully reproduces the crossover from the U-scale to the J-scale physics, as was found earlier for electron-doped cuprates 23 .
For more insight into the strong coupling limit, we note that in the two-particle self-consistent approach 77 , λ is determined by a sum rule involving double occupancy, which is fixed by assuming = ↑ ↓ U U n n n / / ( /2) , sp 2 which leads to a saturation of U sp as U → ∞ (or ↑ ↓~ñ n U J 1/ ). In our calculation, this saturation arises naturally, since U sp can never exceed U c .

VHS and excitons.
While we have demonstrated the important role of the VHS in pseudogap physics, we have not addressed the deeper question of the physical significance of the VHS. We suggest that the significance is primarily topological: the saddle-point VHS is associated with the crossover of a Fermi surface from electron-like to hole-like, and hence represents the point nearest to electron-hole symmetry. This explains its strong role in transport, Fig. 1. While electron-electron scattering is strong, it usually plays only a small role in transport, since if an electron-electron scattering conserves energy and momentum, it does not change the net current. In contrast, any electron-hole collision has a large effect on current 78 . Now strong electron-hole scattering is indicative of a tendency to exciton formation, and we believe our results can most easily be understood in terms of a VHS-driven excitonic instability. Thus, the large entropic effects are related to a broad spread of the electrons in k-space, consistent with a localization of the real-space electron-hole separation in an exciton. A two-band excitonic insulator transition is analogous to a superconducting BCS-BEC transition; we suggest that in our one-band case this becomes the Mott-Slater transition we have found, with the Mott (BEC) phase dominated by preformed excitons. It is often stated that excitonic effects are unimportant in metals due to screening. However, we are studying the Hubbard model, where excitons could be present near half-filling, and their evolution with doping could contribute to a metal-insulator transition.

Conclusions
The present paper has provided significant progress toward DFT corrected MBPT. (1) We have demonstrated the significant role of bosonic entropy in driving short-range order. Our computational scheme incorporates the DFT-Lindhard function, and should capture strong correlation effects, at least for T N ~ 1/U. (2) We have shown that the DFT-Lindhard function can be used to classify phase transitions, and may yield a condensed matter analog of SGAs. (3) We have shown that entropy effects in a one-band model involve a non-FS contribution associated with the VHS. In turn, these entropic effects modify the VHS physics, causing the spectral weight associated with the VHS [in the heat capacity and susceptibility] to shift rapidly towards half-filling as T increases, instead of following the evolution of the (q = 0) DOS. This large shift suggests that the role of the VHS in pseudogap physics needs to be reexamined.
Finally, we comment on the limitations of the present mode-coupling framework. In particular, one would like to go beyond the present results to full inclusion of vertex corrections and excitonic physics. Another issue is the evolution of the Fermi surface in the regime of short-range order. As the FS shrinks from a large to a small size, the susceptibility will also change, and drive changes in the fluctuation maps and the associated properties discussed in this study such as the appearance of a second phase out of a regime of short-range order of a primary phase. With a DFT-corrected Lindhard function, one should be able to obtain insight into these issues.