Mottness induced superfluid phase fluctuation with increased density

Recent observation of diminishing superfluid phase stiffness upon increasing carrier density in cuprate high-temperature superconductors is unexpected from the quantum density-phase conjugation of superfluidity. Here, through analytic estimation and verified via variational Monte Carlo calculation of an emergent Bose liquid, we point out that Mottness of the underlying carriers can cause a stronger phase fluctuation of the superfluid with increasing carrier density. This effect turns the expected density-increased phase stiffness into a dome shape, in good agreement with the recent observation. Specifically, the effective mass divergence due to"jamming"of the low-energy bosons reproduces the observed nonlinear relation between phase stiffness and transition temperature. Our results suggest a new paradigm, in which unconventional superconductivity in some strongly correlated materials is described by physics of bosonic superfluidity, as opposed to pairing-strength limited Cooper pairing.


I. INTRODUCTION
Superfluid phase stiffness is known to suffer from low density, due to the canonical conjugation between carrier density and phase φ ∆n∆φ ∼h [1]. Since the density fluctuation is very limited at low density, the phase must then fluctuate significantly, leaving only a small portion in the well-defined phase of the superfluid component. This quantum mechanical limitation on superfluid density provides a natural explanation [1,2] for the observed reduction of the superconducting transition temperature T c in the underdoped cuprates [cf. δ < 15% of Fig. 1(a)]. As shown schematically in Fig. 1(a), in the underdoped region (δ < 15%), T c decreases as the phase stiffness is weakened at lower carrier density (gray dotted line). In that case the superconducting transition temperature T c would be controlled by the phase stiffness, as opposed to the strength of the pairing.
Unexpectedly, recent measurement [3] of penetration depth λ in the overdoped cuprates (δ > 15%) found a surprising reduction of superfluid stiffness upon increasing carrier density. Correspondingly, the low-temperature superfluid phase stiffness (∝ λ −2 ) forms a dome shape against doping, qualitatively different from the simple proportionality to carrier density [4] expected in the pairing strength-limited BCS theory. Furthermore, the low-temperature phase stiffness even scales with T c in the overdoped regime, obeying the same universal Uemura relation [5][6][7] well-known in the underdoped regime. Since this relation implies a phase coherence-limited superfluidity, this new data apparently reveals that even in the overdoped side, the dominant physics is still the phase fluctuation. Nonetheless, it is not obvious at all why superfluid stiffness should reduce at high carrier density.
In fact, several previous studies already raised similar doubts against the common belief of BCS-like amplitude fluctuation in the overdoped cuprates, and faced the same puzzle * corresponding email: weiku@sjtu.edu.cn how the superfluidity can suffer at higher density. As shown in Fig. 1(b), ARPES observed [8,9] a nearly doping independent superconducting gap ∆ near momentum (π/2, π/2) over a wide doping range(∼7-20%), in great contrast to the expected proportionality to T c in the BCS theory. Similarly, Fig. 1(d) shows that even earlier, the observed momentumdependence of superconducting gap-induced weight transfer deviates qualitatively from the expected d x 2 −y 2 form of the superconducting order parameter, also from the underdoped regime all the way to the overdoped regime. Both observations suggest more commonality between the underdoped regime and the overdoped one than previously expected. Consistently, variational Monte-Carlo (VMC) calculations [10][11][12], which ignore phase fluctuation, all found that the pairing still exists beyond 25% doping [black dashed line in Fig. 1(a)]. In this case, the entire superconducting dome would be deep inside the pairing region of the phase diagram, and thus the demise of superconductivity in the overdoped region must be controlled by additional phase fluctuation of unknown origin.
So, the key scientific question is how it is possible to host such a strong superfluid phase fluctuation at such high carrier density? Particularly in the overdoped cuprates, how can the phase quantum-fluctuates even stronger with increased carrier density, as to wipe out superconductivity at roughly δ > 25%? Why is this critical doping level seemingly universal across different families of cuprates and the new nickelate superconductors [13,14]? What determines this special doping level? What is the origin of the striking similarity between the underdoped and overdoped regime, especially near the quantum critical points (the end of the dome) δ ∼ 5% and 25%? Finally, what essential nature do these puzzles reveal about the unconventional high-temperature superconductivity in the cuprates?
Here we show that all these important questions can be answered naturally in the scenario of the emergent Bose liquid [16][17][18][19][20][21]. Essentially, the high-energy local Coulomb repulsion that prevents double occupation of carriers in each lattice site at low energy (so-called "Mottness") leads to a rather large hard core of the superfluid carriers at even lower temper- Beyond the antiferromagnetic (AFM) phase at very low doping, the superconductivity exists in a "dome" region of the phase diagram, much smaller than the one predicted by the VMC calculation. This suggests that below the classical phase coherent temperature, T θ , additional phase fluctuation must exist in the region in gray.
(b) Observed superconducting gap size (triangles) near momentum (π/2,π/2) [8,9] displays a doping dependence very different from that of superconducting temperature (purple line). Normalized spectral weight transfer ∆W from condensed EBL (c) and ARPES [15] (d)both show a d-wave form distinct from the cos(k x ) − cos k y form of the superconducting order parameter, against expectation of BCSlike theories.
ature. Such a large hard core would unavoidably cause a serious jamming in carriers' kinetic motion, and in turn damage the superfluid phase coherence at higher density, a trend opposite to the common lore. Below we demonstrate this general mechanism of stronger superfluid fluctuation at higher density via simple close-pack estimation and Gutzwiller approximation through statistical counting and variational Monte Carlo, all showing a consistent destruction of superfluidity at around 25-30% carrier density. Specifically for the cuprates and nickelates, this classically intuitive mechanism explains the robust upper limit of doping level in the entire families of cuprate and nickelate superconductors. Combined with the persistent consistency of optical [21], transport [19,21], quasiparticle [16,18,20], and superconducting [17] properties of the exact same model with corresponding experimental observations, our results reinforce the notion of a new paradigm that the high-temperature superconductivity in the cuprates (and likely nickelates) is dominated by physics of bosonic superfluidity, as opposed to pairing-strength limited Cooper pairing.
The contents of this paper are organized as follows. Section II introduces the emergent Bose liquid model and the basic assumptions. In Section III, we discuss the main results of this study. Section V gives the connection of our model to the case of the cuprate superconductors. 0.0 0.2 (π,π) (0,0) (π,0) (π,π) (c) Model of an emergent Bose liquid. (a) Hopping of a bondcentered boson (solid ellipse) to nearest τ and next nearest τ neighboring bonds (empty ellipse), giving rise to a two-band dispersion (right panel), colored in red/blue to reflect the dand s-wave nature. (b) d-wave boson corresponding to the low-energy effective one-orbital description with one band dispersion (right panel), requiring 16 forbidden bonds. (c) Extended hard-core constraint among site-centered d-wave bosons, up to third neighboring sites. (d) Simple close-pack estimation of the highest boson density with minimal non-zero kinetic energy allowing hopping within two-sites on average.

II. MODEL
In the emergent Bose liquid (EBL) model [16][17][18][19][20][21], the charge carriers are assumed to be mostly in a tightly bound two-body state involving the nearest neighboring atoms, below some energy scale, say 150meV, due to various high-energy physics, for example, bi-polaronic correlation [22], two-dimensional short-range anti-ferromagnetic correlation [23,24] and/or cancellation of topological spin current [25]. (For example, in cuprates this bosonic two-body state would consist of doped holes of opposite spin residing in neighboring Cu sites, using the half-filled Mott insulator with strong short-range antiferromagnetic correlation as reference.) Correspondingly, the probability for the fermionic carriers in an unbound one-body state is relatively low below this energy scale. Therefore, most physical properties of the system are dominated by the emerged bond-centered bosonic carriers [c.f. solid ellipse in Fig. 2(a)] that reside in the corresponding bond lattice of the original fermionic lattice [black cross in Fig. 2(a)]. Note that this assumption of EBL does not require a complete depletion of the probability for low-energy fermion (or fermionic spectral weight), but only its overall smallness. (More detailed considerations on low-energy fermions are available in Appendix A.) In addition, in many strongly correlated materials, intraatomic electronic repulsion is very strong, well above the energy scale of this effective low-energy model, resulting in an exclusion of double occupation [26][27][28] (so-called "Mottness") in the low-energy sector of the fermionic carriers. This gives rise to an 'extended hard-core constraint' of EBL, namely forbidden occupation of the vicinity of a boson by another, as illustrated by the shadowed area in Fig. 2(a).
The kinetic process of the EBL is dominated by the pivoting motion [16] in the emerged bond-centered lattice. For the specific case of a square lattice of fermions, for example the Zhang-Rice singlets [29] in the cuprates, the bond lattice is a checkerboard lattice in Fig. 2(a) with two orbitals in the unit cell corresponding to the vertical and horizontal bonds: where b l denotes the annihilation of a boson located at bond l, and τ ll = τ or τ is the fully dressed kinetic process involving the nearest and second-nearest neighboring bonds. In this study, we employ the exact same set of parameters τ and τ previously extracted from ARPES dispersion data [16,30,31] (cf. Appendix C.) It is important to emphasize that the low-energy effective Hamiltonian of EBL is under the above-mentioned extended hard-core constraint, which prevents phase separation related to clustering of the emergent bosons. At low temperature, it also helps to establish the superfluid phase stiffness in the superfluid phase. As demonstrated below, this constraint is the key factor that generates the surprising physical behavior of interest in this study. Now, since the focus of our study is on the phase stiffness of superfluidity, which involves only the lowest energy excitation, it is advantageous to further integrate out irrelevant (higher-energy) kinetic degree of freedom of Eq. 1. Figure 2(a) shows that within the bosonic band structure, the lower-energy band is dominated by d-wave symmetry due to the positive τ . We proceed to decouple the d-wave boson from the s-wave one via a canonical transformation. Practically, this is performed approximately via construction of the d-wave bosonic Wannier orbital corresponding to the lowerenergy band in Fig. 2(a), and then representing the effective Hamiltonian in this decoupled subspace. (See Appendix B for more details.) Essentially, this is equivalent to assuming the formation of a well-defined lowest energy quasi-particle with a d-wave structure responsible for the superfluidity.
Finally, we arrive at the resulting lowest-energy one-orbital Hamiltonian that controls directly the physics of phase fluctuation,H corresponds to a larger d-wave Wannier orbital, denoted by d † i creation operator, now centered at sites i in the original square lattice [cf.: Fig. 2(b)]. Note that due to the shift of the Wannier center, the d-wave boson is not strictly local, but instead with a power-law decay similar to the Zhang-Rice singlet [29]. Correspondingly, the d-wave boson can now hop beyond the first and second neighbors with further renormalizedτ ii given in Tab II.
Most importantly, these low-energy superfluid related dwave bosons now acquire an even more extended hard-core constraint [shadowed area in Fig. 2 [16,17] containing additional mass divergence around δ ∼ 5% doping.

III. MOTTNESS INDUCED SUPERFLUID PHASE FLUCTUATION WITH INCREASED DENSITY
Concerning the main puzzles of the strong phase fluctuation and diminishing superfluidity density at high carrier density, we will now show that the EBL also provides a simple resolution through its "extended hard-core constraint". Such a large extended hard-core implies that the essential kinetic processes will be easily suppressed at moderate density due to blocking and jamming between the low-energy d-wave bosons. Figure 2(d) shows a simple estimation of the special case, in which each of the bosons reach minimum non-zero mobility on average, being able to hop back and forth between two sites. The corresponding density, 2 holes (1 boson) / 8 atomic sites = 25%, marks the approximate maximum doping level of superfluidity, above which the d-wave bosons become nearly impossible to move and thus unable to maintain phase coherence. Interestingly, this 25% coincides very well with the experimentally observed end of the superconducting dome in overdoped cuprates in general [32]. Therefore, this mechanism offers a simple and natural explanation of the key puzzle of strong phase fluctuation in the overdoped cuprates.
We employ the well-known Gutzwiller g-factor approximation [33,34]τ ii → g ii τ ii as the simplest way to capture this renormalization of the kinetic process. The g-factors are calculated via two numerical approaches (See Appendix D). First, we count statistically the probability of each hopping under the constraint at various doping level. Figure 3(a) shows that the resulting g-factors are approximately g ∼ 1 − δ /0.3. Second, we calculate the g-factor using VMC based on a noninteracting BEC wavefunction under the constraint. The resulting g-factors in Correlation between low-temperature superfluid stiffness λ −2 and superconducting transition temperature T c near 5% quantum phase transition point [17,37], normalized by λ −2 UD at T c = 17K. (b)Superfluid phase stiffness λ −2 vs doping level δ showing a dome shape [3], normalized by λ −2 OP at δ = 0.16. Theoretical overestimation in the overdoped region is expected from the Gutzwiller treatment. (c) same as (a) for the overdoped region [3], normalized by λ −2 OD at T c = 11K. (d)-(f) Renormalized band structure of EBL near ∼ 5%, ∼ 15%, and ∼ 25% doping, showing similar effective mass divergence of low-energy band near transition points, but with distinct higher energy features.
As expected, both results decrease rapidly and become rather small beyond 25% doping and diminish around 30%.
We now demonstrate the "dome" shape of T c with our simple picture. Estimated from the standard BEC [35] of a uniform boson gas, the condensation temperature T c ∝ n 2/3 /m * is a simple function of density n and effective mass m * . Obviously, followingτ, 1/m * and T c are also renormalized by g. Figure 3(c) plots T c ∝ δ 2/3 (in black) and the renormalized T c → gT c (in red). Indeed gT c is strongly suppressed at higher doping and eventually vanishes around 30%, where the average m * diverges. This δ ∼ 30% upper bound of superfluidity is in excellent agreement with Tl 2 Ba 2 CuO 6+δ [36] and well consistent with the common ∼ 25% limit of typical cuprates superconductors. If one further incorporates the previously proposed doping dependent bosonic mass that diverges at around 5% due to a level crossing [17], the gT c in Fig. 3(d) reproduces the experimentally observed dome shape quite well. In short, the phase fluctuation can indeed grow in the overdoped regime when the suppression of kinetic processes overcomes the increasing density.
Similarly, the low-temperature limit of the phase stiffness λ −2 would also suffer from the jamming reduced kinetic process. Phase stiffness λ −2 ∝ n s /m * is proportional to the superfluid density n s . Near the above-mentioned mass divergence dictated quantum phase transition points, n s ∝ (1/m * ) β δ must vanish with 1/m * with a exponent β > 0, leading to λ −2 ∝ (1/m * ) 1+β δ . Figure 4(b) shows that this phase stiffness indeed forms a dome shape, vanishing at around δ ∼ 5% and ∼ 30%. Particularly, due to the extended hard-core constraint of the EBL, the phase stiffness reduces in the overdoped regime against the growing density, as observed in recent experiments [3,4]. (Our results' deviation from exper-iment merely reflects the above-mentioned overestimation of phase coherence at high doping in our simple Gutzwiller treatment.) Furthermore, the effective mass divergence of EBL has an important physical consequence in the relation between λ −2 and T c [17]. Estimated from T c ∝ δ 2/3 /m * for a free boson gas, one finds λ −2 ∝ T 1+β c . Therefore, distinct from the simple linear relation expected from standard phase fluctuation scenario, λ −2 vs T c must show a zero slope as T c approaches zero (since β > 0). This provides a natural explanation of the super-linear relation in the experimental observations shown in Fig. 4(a) and (c). In fact, with β = 1 this formula describes very well the observed relation in the entire low-T c region in both underdoped and overdoped sides. Not surprisingly, the resulting scaling λ −2 ∝ T 2 c agrees with that of the 3+1D XYmode that describes quantum phase fluctuation of superfluidity upon suppressed kinetic energy [38,39].
Such an apparent "symmetry" near the quantum phase transition points between the underdoped and overdoped sides is puzzling within the current lore, given that one is supposedly governed by fluctuation in the phase and the other in the amplitude. Especially, the observed physical properties of the normal-state above T c in these two regions are quite distinct [32,40,41]. In our EBL picture, the contrast in nonsuperfluid properties is easily understood from the rising importance of the jamming effect discussed above. Concerning the higher-energy features, Fig. 4(d)-(f) show that in the overdoped region the renormalized s-wave (in blue) and d-wave bosons (in red) suffer a significant renormalization, quite distinct from the underdoped and optimally doped region. Yet, the renormalized low-energy coherent d-wave bosonic band near the underdoped (d) and the overdoped transition points becomes heavier (flatter) in a similar fashion (through very different microscopic mechanism though.) In other words, the low-energy phenomenon of superfluidity experiences a similar loss of phase stiffness due to the effective mass divergence despite the distinctly different underlying higher-energy physics.
From this perspective, the non-superconducting "normal state" above 25% doping should have interesting properties. While energetically most favorable locally, the pure d-wave boson would suffer from severe jamming and loss of kinetic energy. It would thus tend to morph into a p-wave boson whose cigar shape allows them to remain mobile at a much larger doping (See Appendix E). Therefore, the EBL is expected to behave as an unconventional Bose metal consisting of fluctuating dand p-wave bosons. Alternatively, since at very high doping, the electronic correlation will eventually become weaker, the emergent bosons might lose their binding strength and start to decompose into weakly bound fermions. Concerning this scenario, we note two very relevant experimental observations: 1) robust paramagnon dispersion at 40% doping [23] indicating that the local magnetic correlation remains strong and 2) linear resistivity at 30% doping [3] indicating that the normal states is still not a fermionic quasiparticle based "normal metal". Both observations suggest that decomposition of the emergent boson is not immediate after the disappearance of superconductivity.
It is also interesting to realize that in every energy scale our EBL experiences important influence of the "Mottness" of the underlying doped holes through the suppression of their double occupation at each Cu site due to large local interaction. First, with the help of kinetic energy, it lays the foundation of near neighboring antiferromagnetic [42] and bipolaronic [22] correlations that provide a strong tendency to bind doped holes into pairs [26,43]. Then, it forces the bound pairs to occupy two atomic sites and thus indirectly establishes d-wave form [17]. Finally, it produces the extended hard-core constraint of the EBL that ultimately induces jamming of lowenergy d-wave bosons at high enough density. This jamming in turn suppresses the superfluid phase coherence especially with higher density.
Since the dominant physics of jamming inherits from the Mottness at very high energy scale, details of the intermediate-and low-energy scale is therefore not as essential, as long as the generic assumptions of our EBL model remains valid. Interestingly, the recently discovered Nd 1−x Sr x NiO 2 also demonstrates a superconducting critical doping at around 25% [13,14], similar to our analysis above. This strongly suggests that the nickelate superconductors can also be described by an EBL, with strong jamming-induced phase fluctuation in the overdoped regime. Future experimental verification of this prediction would be highly valuable.
In short, we show that emergent Bose liquid built from carriers under strong intra-atomic repulsion (Mottness) provides a generic mechanism that reduces superfluid stiffness at higher density due to jamming of its kinetic process at around 25-30% carrier density. This mechanism naturally explains the recent observation of strong phase fluctuation in the overdoped cuprates [3,44] and likely to be found also in the nickelates. Assuming that this mechanism is indeed responsible for the puzzling observation in the cuprates, it further suggests that in the entire doping range of the superconducting dome, the unconventional superconductivity can be described by bosonic superfluidity for example of an EBL, without resorting to a crossover into BCS-like amplitude-fluctuating descriptions.

IV. KEY CHARACTERISTICS OF THE BOSON IN THE EBL MODEL
The EBL model assumes a high-energy binding between the nearest neighboring carriers such that the unbinding processes can be integrated out from the low-energy Hilbert space. Therefore, the bosons in EBL model are in a tightly bound two-body state of fermionic carriers of the nearest neighboring position (instead of momentum) with energy spanning far away from the chemical potential (since the binding is assumed strong). Most importantly, with strong shortrange magnetic correlation in the background, the two-body state is not necessarily a spin singlet. Due to their emergent nature, they contain the following key characteristics that distinguish them from other boson-like carriers proposed in the literature, such as bosonized Cooper pairs [45][46][47], resonating valence bond (RVB) [48][49][50]: 1. They are in a tightly bound two-body state of fermionic carriers (not necessarily a spin singlet, given a strong shortrange magnetic correlation.) 2. They consist of fermionic state with energy spanning far away from the chemical potential, since the binding is assumed strong.
3. They are small in size and have a well-defined location, with underlying fermions residing in the nearest neighboring position (instead of a fixed total momentum as in a Cooper pair).
4. They have a rod-like shape with two ends at nearest neighboring atoms.
5. Correspondingly, they reside in the bonds between atomic sites.
6. The bosonic lattice is therefore the bond lattice, different from the underlying fermionic lattice.
7. For systems of 2D or higher dimension, the multiple directions of bonds imply multi-orbital nature of EBL.
8. The bosons in EBL carry a limited pivoting motion, moving one end at a time. This paves ways to frustration in the kinetic energy and can therefore heavily renormalize the kinetic energy of the first bosonic band (to ∼ 30 meV for the cuprates.) 9. They inherit an extended hard-core constraint from the large on-site repulsion between the underlying fermionic carriers.

V. APPLICABILITY OF EBL TO THE CUPRATES
Discovered almost thirty years ago [51], the exotic phenomenon of high-temperature superconductivity (HT-SC) in the cuprates still remains puzzling to researchers. Conventional superconductivity, a state of the matter that shows no resistance in conducting current, is well described by the standard "BCS" theory [52] via weakly bound "Cooper pairs" that fluctuate in amplitude. On the other hand, the HT-SC in cuprates shows qualitatively different behavior. For example, in the weakly hole doped ("underdoped") region, the isotope effect increases dramatically and yet the corresponding superconducting transition temperature T c decreases [53] instead. The observed superconducting gap ∆ 0 in the cuprates is typically significantly larger than the canonical value of twice the transition temperature ∼ 2k B T c in BCS theory [54,55]. During the phase transition, the measured specific heat shows two clear kinks 10K apart [56], qualitatively different from a standard second-order phase transition from the BCS theory. Furthermore, in the underdoped region, the low-temperature specific heat shows no T 2 contribution expected from the observed d-wave quasiparticles, but only a dominant T 3 instead [56]. In addition, the observed upper critical field H c2 does not saturate at low temperature and sometimes even exceeds the Pauli limit [57]. These qualitative distinct features indicate clearly that the HT-SC in the cuprates is of a different nature.
The assumptions of the EBL model are actually quite applicable given the strong high-energy spin and lattice correlations observed in the cuprates [22,24,42]. Particularly, the EBL model is in good consistency with the recent observations of bosonic behaviors in the cuprates [58][59][60], and is obviously the most intuitive one to account for the observed dominant role of phase fluctuation in the superconducting phase transition of the cuprates [1,3]. Below we revisit several recent studies that demonstrated persistent agreement of optical [21], transport [19,21], quasiparticle [16,18,20], and superconducting [17] properties of EBL with corresponding experimental observations. The fact that these successful applications of EBL are all via the same Hamiltonian with the same set of parameters strongly supports the scenario in which charge-related physics in the cuprates can be understood intuitively by a simple EBL.

A. Doping dependent Superconducting gap
Previously this EBL model was shown to offer a simple explanation for the disappearance of superfluidity in cuprates at δ ∼ 5%, as the effective mass of the emergent boson diverges [17]. (The same study also suggested a scenario for the lack of superfluidity below 5% doping.) In a related study, a second kind of "superconducting gap" was found in the quasi-particle spectrum, resulting from coherent kinetic scattering against the Bose-Einstein condensation (BEC) of the EBL [16]: where f d (k) = 1 2 [cos(k x ) − cos(k y )] gives the momentum dependent d-wave form factor of the order parameter and n 0 (T ) the temperature-dependent condensation density. ε k − ε 0 denotes the quasi-particle energy measured from the low-energy band center, and signifies the kinetic origin of the gap (as opposed to the pairing potential). Beyond the "Fermi arc", where ε k moves from the chemical potential toward the band center, it provided additional momentum dependence [16] that reproduced nicely the ARPES measurements [15] at various doping levels [cf. Fig. 1(c). ] We first observe that this model actual "predicted" naturally the above-mentioned nearly doping independent super-conducting gap in ARPES measurements [8,9]. Since the quasi-particles on the Fermi arc are at the chemical potential µ, ε k − ε 0 is given by µ − ε 0 ∝ 1 − 2δ upon hole doping into the quasi-2D singlet band of the cuprates. Away from the ends of the dome, in the absence of strong quantum fluctuation, the condensation density n 0 (T = 0) is roughly proportional to the superfluid density, n s (T = 0) ∝ T c . Together, this gives a weakly doping dependent near-node gap scale √ T c , very similar to the experiments shown in Fig. 5, but in great contrast to the stronger dome shape of T c . Particularly, notice that the 1 − 2δ factor shifts the maximum of V ∆ to a significantly lower doping δ ∼ 0.12, compared to T c . Similar to the above-mentioned deviation from d-wave form of the gap, this again reflects the kinetic nature of the anomalous scattering gap absent in typical weak-coupling pairing scenario.

B. Transport, quasiparticle and superconducting properties
Here, we recite our recent studies of various non-trivial physical properties of EBL in its non-superfluid and superfluid phase, which all demonstrate excellent agreement with experimental observations of the cuprates, within a single set of parameters and the exact same Hamiltonian of EBL.
First, the transport properties of EBL in the non-superfluid phase are generically non-Fermi liquid like, not surprisingly. Its optical conductivity [21] contains a featureless continuum extending a broad (> 300 meV) energy range, with a midinfrared peak around 100 meV, nicely reproducing the puzzling experimental spectra. The corresponding DC resistivity [21] demonstrates no saturation at high temperature (the "bad metal" behavior), a linear temperature dependence at a large temperature range all the way to very low temperature (the "strange metal" behavior), and even sometimes an insulator-like upturn at lower temperature (the "weak insulator" behavior). These exotic non-Fermi liquid transport behaviors all exist at different temperature range of the same system, just like what is observed in the cuprate samples. When τ > τ (corresponding to extremely underdoped δ < 5% cuprates), EBL has been numerically [17] and analytically [19] proven to be a homogeneous Bose metal that cannot superflow and acquires finite resistivity with temperature and/or disorder without the protection of an energy gap, resembling the low-temperature pseudogap phase of the hole doped cuprates. Under moderate magnetic field, not surprisingly EBL displays no quantum oscillation (since it has no Fermi surface). The corresponding Hall coefficient is highly temperature dependent just like the cuprates observation, and does not give directly the carrier density [61].
Second, upon scattering against EBL, the residual fermionic quasi-particle (QP) demonstrates strong non-Fermi liquid behaviors as well. The EBL density of state contains a clear ∼ 50 meV structure and a large continuum extended to 300 meV [18], resembling perfectly the Eliashberg function α 2 F extracted from ARPES. The corresponding QP selfenergy demonstrates two well-defined features at ∼ 25 meV (that render low-energy QP heavier) and ∼ 50 − 70 meV (that generates a "kink" in the QP dispersion), in excellent agreement with ARPES observation. Due to the non-vanishing density of the EBL at low temperature (contrary to phonon and magnon), the normal-state QP demonstrates non-Fermi liquid scattering rate [18,20] at the chemical potential, just like the ARPES observations. Furthermore, since it is the same set of bosons the QP scatters against in the non-superfluid and the superfluid phase, the incoherent normal-state scattering rate is trivially proportional to the low-temperature coherent superconducting QP gap [18], explaining the ARPES observed puzzling connection between the incoherent and the coherent features. Specifically in the Bose metal phase that corresponding to the pseudogap phase, scattering against EBL generates a Fermi arc with an obvious pseudogap [20] in the anti-nodal region that fills up (instead of closes) at higher temperature, same as the observations by ARPES and scanning tunneling spectroscopy (STS). For fixed momentum, the resulting pseudogap displays a strong asymmetric gap edge in energy, while for a fixed site it has a perfect symmetric gap edge [20], reconciling the puzzling qualitative inconsistency between ARPES and STS.
Third, in the superfluid phase, coherent scattering against EBL generates a second kind of QP superconducting gap [16] whose energy scale corresponds to the effective kinetic energy, instead of the pairing energy. This generates additional momentum dependence beyond that of the d-wave superconducting order parameter, similar to the weight transfer observation by ARPES. Furthermore, this analytically and quantitatively reproduces the correspondence between the QP dispersion and the superconducting gap observed by ARPES.
Finally, at low temperature, EBL was shown to host a dwave superfluid [17], scattering against which the QP display a nodal direction between the second neighboring Cu-Cu directions [16], consistent with the observation. The EBL superfluidity suffers from strong phase fluctuation associated with the mass divergence [17] due to a local level crossing at τ → τ , corresponding to the observed quantum fluctuation near 5% hole doping. And, in this article, we show that a similar mass divergence has to take place at around 25-30% carrier density, but this time due to the robust jamming effect induced by the Mottness of the electronic structure. This effect leads to an unusual physical trend of weaker superconducting phase stiffness at higher density, as recently observed experimentally.
We stress that these persistent successful applications of the EBL model are all via the exact same Hamiltonian with the same set of parameters as the ones used in this article. Combined with recent observations of bosonic behaviors in the cuprates [58][59][60], these results lay strong support for a consistent EBL-based scenario for the cuprates. Further experimental confirmations, for example direct observation of the bosonic quasi-particles and their dispersion [c.f. Fig. 2(a)] via one-photon-in-two-electron-out photo-emission process, or observation of the EBL-predicted second superconducting dome [19] below 5% doping under a strong (110) uniaxial pressure, would offer stronger verification (or falsification) of the suitability of the EBL in describing the low-energy physics of charges in the cuprates.

VI. SUMMARY
In summary, within the scenario emergent Bose liquid with extended hard-core constraint, we propose a generic mechanism to explain the puzzling reduction of superfluid phase stiffness and transition temperature at high carrier density. Essentially, the strong intra-atomic repulsion (the Mottness) leads to a large extended hard-core constraint for the emergent bosons that can be easily jammed in their motion at high density and thereby host weakened superfluid stiffness. Our results are well consistent with the experimental observations in some superconductors, such as the cuprates and the nickelates. Furthermore, together the persistent consistency of the EBL with many key spectral and transport observations of the cuprates, this new finding adds credibility to the notion that EBL model can provide a good description for the charge channel of the low-energy electronic structure of the cuprates.
Note added: After initial submission of our manuscript, the shot noise experiment proposed in our original manuscript was performed [59], whose data indicated indeed a 2e-charge quanta, subject to further independent experimental confirmation. In the final version of this manuscript, we therefore remove this suggestion and instead incorporate this new development in the reference. It's a common question whether one is allowed to employ a bosonic picture like the EBL model to systems with experimentally observed low-energy fermionic quasi-particles, such as the cuprates [40,62]. As shown in Figure 6(c), if one would consider only the BEC limit corresponding to an overwhelming binding energy E B , the low-energy fermions indeed should be completely depleted with a clean gap of the scale of E B − W around the chemical potential in the fermionic onebody spectral function. (Here W denotes the bandwidth of unbound fermions.) However, in general, with intermediate binding strength (E B slightly larger than W /2) illustrated in Figure 6(b) [63] one often finds a small probability of the unbound fermions (in blue) being occupied and forming the Fermi surface. Even though, fermions in the system are actually predominately part of the well-defined bound bosonic state, as reflected by the larger spectral weight in red below the chemical potential. In this scenario, the rare gapless fermionic quasi-particles' contributions to low-energy physical properties are generally overwhelmed by those of the more likely bosonic carriers. Thus, a bosonic description of the lowenergy physics, such as the EBL model used in this study, would be the most convenient even in the presence of rare gapless unbound fermions. Importantly, in contrast to the weak influence of the rare unbound fermions on the bound fermion, the latter is however expected to have significant impact on the former. Particularly at low temperature, contrary to the depletion of typical bosonic excitations such as phonon and magnon, the number (or more precisely the probability) of these bound bosons is basically fixed. Therefore, scattering against the bosonic carrier can lead to many unusual non-Fermi-liquid behaviors of the rare unbound fermions that cannot result from phonon or magnon, such as kink structure [18], Fermi arc and pseudogap [20] (c.f. Section V of the main text.) In essence, the unbound fermions can sensitively reflect the influence of bosonic carriers and serve as good "probes" to the more essential bosonic carriers.

VII. ACKNOWLEDGMENTS
Appendix B: Construction of the lowest-energy orbital and the corresponding one-orbital Hamiltonian Since the focus of this study is on the superfluid stiffness, our main interest is the lowest-energy fluctuation around the superfluid. We therefore aim to construct a low-energy quasiparticle whose local structure absorbs the influence of the high-energy orbitals. Conceptually, this is similar to integrating out τ to construct a local Wannier orbital with a procedure similar to that of the Zhang-Rice singlet, where the nearest neighbor hopping processes between two oxygen are integrated [29]. As long as the energy of resulting local object is lower enough than other state, this approximation is controlled.
In practice, it is easily achieved numerically by simply constructing the d-wave Wannier function from the lowerenergy band in Fig. 2(a). The resulting energy of d-wave boson is ε d = τ − 2τ = −41.2meV which is lower than the other states, such as p-wave boson with energy ε p = −τ = −25.8meV (with the parameter τ = 33.5meV, τ = 25.8meV,  [16,30,31] first and second nearest neighboring hopping parameters τ and τ of bosonic carriers in Eq. 1 at different doping level δ of La 2−δ Sr δ CuO 4 . Note that these parameters are kept fixed in all existing studies [16][17][18][19][20][21]61]  see Table I). With the Wannier orbital, the lowest-energy effective one-band Hamiltonian Eq. 2 can be obtained by representing Eq. 1 within the subspace of d-wave boson, similar in spirit to the construction of the t-J model from the Zhang-Rice singlet [29].

Appendix C: Parameters of EBL Hamiltonian
The key physics of interest in this study, namely enhanced superfluid phase fluctuation at higher density, is dominated by the extended hard-core constraint inheriting from the highenergy physics of Mottness. Therefore, it is completely independent of the parameters of Eq. 1. Nonetheless, this section provides the parameters for interested readers. Table I lists the parameters of Eq. 1 used in this study. These parameters were initially extracted [16] from ARPES experiments [30,31] on La 2−δ Sr δ CuO 4 , and kept fixed in all the studies [16][17][18][19][20][21]61] using the same EBL model. Other families of cuprates might have slightly different parameters, but we don't expect any significant deviation. Table II gives examples ofτ ii of the d-wave boson at δ = 15%, corresponding to the band in Fig. 2(b). They are derived from integrating out the s-wave boson corresponding to the higher-energy band in Fig. 2(a).
states. Specifically, the definition of the Gutzwiller g-factor involves the ratio of thermal average of the hopping process g ii = d † i d i c / d † i d i 0 evaluated using low-energy states with and without the constraint (denoted by subscripts c and 0.) In our case, d † i is the creation operator of local d-wave boson. Therefore, this g-factors should in general be energy and temperature dependent, especially around the scale of the interaction that induces this constraint. On the other hand, below this energy scale where the constraint is enforced, the gfactor should be quite energy-and temperature-independent, and thus not very sensitive to the detailed structure of the lowlying states. This makes such statistical counting a rather good approximation of the actual g-factor. Indeed, as shown in the main text, the resulting g-factor resembles the VMC calculation very well. This energy independence (within the scale of our low-energy Hamiltonian) also allows us to apply on an approximate ground state in our VMC calculation.

Variational Monte Carlo (VMC) calculation of Gutzwiller g-factor
Due to the extremely large many-body Hilbert space of bosonic systems, we can only afford to evaluate the Gutzwiller g-factor using an approximate wave function. Luckily, as argued above, the results should be quite insensitive to the choice of wave function.
We calculate the ratio between the expectation values of the hopping term d † i d i in the "Gutzwiller-projected" Bose-Einstein condensation (BEC) state and that in the mean-field BEC state as the Gutzwiller g-factor, g ii , Here, the P G represents the extended hard-core constraint imposed on the bosonic hole-pairs, and the BEC mean-field state |BEC is expressed by the following normalized formula, where M and N represent the total site number and boson numbers respectively, with N = Mδ /2.
The denominator of Eq.(D1) can be easily obtained as, The numerator of Eq.(D1) can be evaluated as, where P α ≡ | α |P G | BEC | 2 represents the weight of each configuration |α in the Gutzwiller-projected wave function P G |BEC , and B α = ∑ β α d † i d i β · β |P G |BEC α|P G |BEC represents the measurement for the configuration |α . The weight P α and the measurement B α for any configuration |α are easily obtained. From Eq.(D2), one easily finds that the weights P α for all the configurations |α which are permitted by the extended hard-core constraint are equal, while those for the constraint-prohibited configurations are zero. The value B α of any constraint-permitted configuration |α is chosen as follows: defining |β = d † i d i |α , if configuration |β is permitted by the extended hard-core constraint then B α = 1; otherwise B α = 0. These formulae provide an appropriate start-point for the following Monte-Carlo calculations.
In the Monte-Carlo calculation, we start from an arbitrary configuration |α 1 . We then randomly select a particle in configuration that let the particle hop to any other hard-core-constraint-permitted position on the lattice to obtain a second configuration |α 2 . This completes an update. Continuing the updates, we obtain a configuration series |α 1 , |α 2 , |α 3 , · · · · · · . After about 10 4 Monte-Carlo steps thermalization can be realized. In successive Monte-Carlo steps, one begins to extract the measurement B α for the configuration |α , and to average these B α to obtain the numerator defined by Eq.(D4). To avoid auto-correlation, we take one measurement after each N Monte-Carlo steps. About 10 5 measurements are performed to attain convergence. The final result for the Gutzwiller g-factor is given by Eq.(D1), Eq.(D3) and Eq.(D4). a p-wave emergent boson would have a total kinetic energy 2ε p = −2τ , half of which is also gained by forming local p-wave symmetry. Therefore, when τ > τ , 2ε d < 2ε p and d-wave form is dominant.
As discussed in the main text, as the density increases upon higher doping, due to the large hard core the d-wave boson starts to suffer from a significant jamming effect. Represented by the Gutzwiller factor g d , a crude estimation of the jamming-induced renormalized kinetic energy can be made: (1 + g d )ε d . Similarly, for the p-wave boson, the energy reduces to (1 + g p )ε p .
Interestingly, as shown in Figure. 7(a) due to the narrower shape of the p-wave bosons, they naturally suffer a lot less from the jamming effect than the d-wave bosons. Consequently, the p-wave boson can survive a much larger doping, approximately δ = 50% based on the estimation in Figure. 7(b). Indeed, Figure. 7(c) shows a much slower decay in g p against δ from a simple statistical counting. Therefore, one can expect that at high-enough doping, the d-wave symmetry should no longer be dominant, since the bosons would be forced to adapt to p-wave whenever they need to pass each other in short distance. Figure. 7(d) shows that for a fixed τ and τ (from optimal doping), (1 + g p )ε p will eventually become lower than (1 + g d )ε d . Therefore, assuming that the emergent Bose liquid remains intact beyond the superconducting dome, the system should become a Bose metal consisting of a mixture of pand d-wave bosons. Notice that such a Bose metal should experience a different structure of main phase fluctuation compared with the recently proven [19] homogeneous Bose metal in the underdoped region, and is an interesting candidate for another homogeneous Bose metal. Such a strong phase fluctuating EBL should present non-Fermi liquid behaviors as well, as found experimentally in Ref. [3,23]. The emergence of p-wave boson is also consistent with the observation of a similar direction of the charge order found in beyond underdoped and beyond overdoped systems [64].