
==== Front
Nat Commun
Nat Commun
Nature Communications
2041-1723
Nature Publishing Group UK London

39231956
52124
10.1038/s41467-024-52124-9
Article
Evolution from a charge-ordered insulator to a high-temperature superconductor in Bi2Sr2(Ca,Dy)Cu2O8+δ
Zou Changwei 12
http://orcid.org/0000-0003-4616-4345
Choi Jaewon 3
Li Qizhi 17
Ye Shusen 2
Yin Chaohui 4
http://orcid.org/0000-0002-6982-9066
Garcia-Fernandez Mirian 3
http://orcid.org/0000-0002-3625-880X
Agrestini Stefano 3
Qiu Qingzheng 1
Cai Xinqiang 1
Xiao Qian 1
http://orcid.org/0000-0002-5261-1386
Zhou Xingjiang 4
http://orcid.org/0000-0001-9293-0595
Zhou Ke-Jin 3
Wang Yayu yayuwang@mail.tsinghua.edu.cn

25
http://orcid.org/0000-0002-2657-3590
Peng Yingying yingying.peng@pku.edu.cn

16
1 https://ror.org/02v51f717 grid.11135.37 0000 0001 2256 9319 International Center for Quantum Materials, School of Physics, Peking University, Beijing, China
2 grid.12527.33 0000 0001 0662 3178 State Key Laboratory of Low Dimensional Quantum Physics, Department of Physics, Tsinghua University, Beijing, China
3 https://ror.org/05etxs293 grid.18785.33 0000 0004 1764 0696 Diamond Light Source, Harwell Campus, Didcot, UK
4 grid.458438.6 0000 0004 0605 6806 Beijing National Laboratory for Condensed Matter Physics, Institute of Physics, Chinese Academy of Sciences, Beijing, China
5 grid.12527.33 0000 0001 0662 3178 Frontier Science Center for Quantum Information, Beijing, China
6 https://ror.org/03jn38r85 grid.495569.2 Collaborative Innovation Center of Quantum Matter, Beijing, China
7 Present Address: Shenzhen Pinghu Laboratory, Building C, Chinese Sciences Vally, Industrial Park (iBT), Shenzhen, China
4 9 2024
4 9 2024
2024
15 77394 1 2023
26 8 2024
© The Author(s) 2024
2024
https://creativecommons.org/licenses/by-nc-nd/4.0/ Open Access This article is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License, which permits any non-commercial use, sharing, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if you modified the licensed material. You do not have permission under this licence to share adapted material derived from this article or parts of it. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by-nc-nd/4.0/.
How Cooper pairs form and condense has been the main challenge in the physics of copper-oxide high-temperature superconductors. Great efforts have been made in the ‘underdoped’ region of the phase diagram, through doping a Mott insulator or cooling a strange metal. However, there is still no consensus on how superconductivity emerges when electron-electron correlations dominate and the Fermi surface is missing. To address this issue, here we carry out high-resolution resonant inelastic X-ray scattering and scanning tunneling microscopy studies on prototype cuprates Bi2Sr2Ca0.6Dy0.4Cu2O8+δ near the onset of superconductivity, combining bulk and surface, momentum- and real-space information. We show that an incipient charge order exists in the antiferromagnetic regime down to 0.04 holes per CuO2 unit, entangled with a particle-hole asymmetric pseudogap. The charge order induces an intensity anomaly in the bond-buckling phonon branch, which exhibits an abrupt increase once the system enters the superconducting dome. Our results suggest that the Cooper pairs grow out of a charge-ordered insulating state, and then condense accompanied by an enhanced interplay between charge excitations and electron-phonon coupling.

A major mystery in cuprate superconductors is how Cooper pairs form and condense. Here, the authors use RIXS and STM techniques to reveal that a charge-ordered insulating phase and enhanced bond-buckling phonons play important roles in this process.

Subject terms

Superconducting properties and materials
Electronic properties and materials
Ministry of Science and Technology of China (Grants No. 2019YFA0308401 and No. 2021YFA1401903) and the National Natural Science Foundation of China (Grant No. 11974029).issue-copyright-statement© Springer Nature Limited 2024
==== Body
pmcIntroduction

High-Tc superconductivity (SC) is closely proximal to a charge order (CO) electronic state in the cuprate phase diagram1, which exhibits an incommensurate wavevector QCO that varies with hole concentration (Fig. 1a and ref. 2). While there is growing evidence suggesting that CO competes with SC since its amplitude is weakened below Tc3–6, it might actually be the first emerged phase during the doping of a Mott insulator1, acting as a precursor to SC. Conforming to this scenario, scanning tunneling microscopy (STM) studies on Bi2Sr2−xLaxCuO6+δ (La-Bi2201)7 and Ca2−xNaxCuO2Cl2 (Na-CCOC)8 reveal that a surface-form CO has already become manifest even in the antiferromagnetic regime, with only a small amount of holes (p = 0.07~0.08) introduced into the parent compounds. However, whether the CO exists in insulating cuprates has not been confirmed by bulk-sensitive probes such as resonant inelastic X-ray scattering (RIXS). Thus, it remains unclear whether the CO has a bulk nature in the non-superconducting regime or exists only in a surface form, and whether it extends to even lower doping levels. Furthermore, it is important to determine whether the CO wavevector QCO is commensurate, as reported in refs. 7,8, to understand its relation with the underlying lattice. Answering these questions is critical in delineating a clear story about the relationship between unconventional SC and CO in cuprates.Fig. 1 Demonstrations of the research with combined RIXS and STM techniques.

a QCO of Bi2212 as a function of doping. The blue shaded line is a guide to the eye. Inset: the schematic crystal structure of Dy-Bi2212. b Scheduled measurements of this work, displaying the three paths (black lines and arrows) in the schematic phase diagram. The orange and purple circles denote the specific temperature and doping levels at which the RIXS and STM data were respectively collected. The blue-filled (Bi2212) and open (Bi2201) stars, together with the quadratic fitting, are suggested by ref. 2. The T* line (light gray) is adapted from ref. 13. AF: antiferromagnetic, SM: strange metal, SC: superconducting. c Representative RIXS spectra of the p = 0.04 sample at 17 K. The data are plotted as functions of energy loss distinguishing different excitations (indicated), as well as of the projected in-plane momentum transfer Q// in the (H, 0) direction. Inset: schematics of the RIXS process, displaying the intermediate (∣m >, resonant absorption from Cu 2p to Cu 3d) and final (∣f >, inelastic emission) states. d Representative tunneling spectra dI/dV (r, V) (roughly proportional to the local density-of-states) of the p = 0.04 sample at 5 K, which has extended the lower doping limit of STM's capability in measuring insulating cuprates at liquid-helium temperature. Inset: schematics of STM and the topographic image that clearly resolves the superstructure and Bi atoms of Bi2212.

One more piece of the puzzle that remains elusive is the origin of CO. This CO is typically observed at relatively high doping levels and is intertwined with pseudogap (PG) and SC through a truncated Fermi surface. For example, angle-resolved photoemission spectroscopy (ARPES) experiments have found discontinuous segments known as Fermi arcs9, and momentum-space nesting mechanism has been suggested as a possible explanation for the origin4,10. On a different front, quantum oscillations11 and magnetoresistance12 measurements suggest that the CO is closely connected with a reconstructed Fermi pocket. However, all these previous pictures fail to explain the above-mentioned charge-ordered insulating phase, where the Fermi surface is absent13,14. To address these issues, we have designed a comprehensive set of temperature paths (Fig. 1b) to explore the property, universality, and extension of the CO, and how the SC emerges in the vicinity of its onset doping level.

We choose lightly doped Bi2Sr2Ca1−xDyxCu2O8+δ (Dy-Bi2212, x = 0.4) to study utilizing the state-of-art RIXS and STM techniques. Dy-Bi2212 has proved to be a feasible system for our purposes that can achieve hole concentration as low as the Mott insulating limit15. Importantly, its crystal structure (Fig. 1a inset) allows an atomically flat cleavage that is ideal for both STM and RIXS. A recent RIXS study on underdoped Dy-Bi2212 (p ≈ 0.09) has identified dispersive excitations associated with the CO, which is proposed to induce a peak of electron-phonon coupling (EPC) strength at the in-plane bond-stretching (BS) phonon branch16. This EPC peak grows once the system enters the SC dome6, providing intriguing connections between SC, CO, and EPC. However, the implications of these findings remain contentious among different cuprate compounds when the inelastic spectral weight is solely attributed to the BS phonon mode16–18. Moreover, the extension of the EPC “anomaly" to lower doping levels, especially to the insulating phase diagram, remains largely unexplored. In this work, we investigate this lightly doped region with p = 0.04 and p = 0.06 (Methods and Supplementary Fig. S1) across the edge of the superconducting dome (p = 0.05), where the Fermi surface is found missing in both samples from ARPES studies13.

Results

RIXS and STM measurements

Figure 1 c shows a series of high-resolution RIXS spectra (ΔE = 36 meV, ΔQ = 0.006 r.l.u., see Methods)19 taken in the insulating sample (p = 0.04) at 17 K. The incident X-ray was tuned to the Cu-L3 edge, corresponding to resonant absorption and emission between Cu 2p and Cu 3d orbitals (Fig. 1c inset and Supplementary Fig. 1), which provides energy- and momentum-resolved information of elementary excitations, such as the paramagnons and phonons. The same sample was cleaved at an ultra-high vacuum and immediately loaded into the STM stage. As displayed in Fig. 1d, we succeeded in performing single-particle spectroscopic and atomically-resolved topographic measurements in Dy-Bi2212 with p = 0.04, by exploiting the reliable methodology we have previously developed (Methods). As the sample is spatially inhomogeneous, the spectra evolve from a charge-transfer gap (black curve) to a PG lineshape (red curve), similar to that exhibited in superconducting samples20.

Observation of charge order

We now focus on the lower-energy range of our high-resolution RIXS and STM data, which is most relevant to CO, PG, EPC, and SC. Intriguingly, in this insulating Dy-Bi2212 we observe an evident bulk CO through careful examination of the energy-momentum RIXS intensity maps (Fig. 2a–c), which has never been reported before. It is directly visualized that at the quasi-elastic region a broad peak along the (H, 0) direction (i.e., Cu-O bond direction) associated with CO has already formed at 273 K and gradually becomes stronger as the temperature decreases (along path #1 of Fig. 1b). The CO intensity profile can be well-fitted by a Gaussian peak (blue shaded) plus a Lorentzian background (Methods), which is centered at an incommensurate QCO  = 0.31 r.l.u. at 17 K. The elastic intensity peak is absent along the diagonal (H, H) direction (gray squares in Fig. 2c), consistent with CO characteristics in other cuprate materials2. Previous studies have only proposed CO down to p = 0.07 using surface-specific STM7,8,21, while a bulk-form CO was only identified in superconducting Bi2212 down to p = 0.09 (Tc = 45 K)16. Our findings substantially extend the lower doping limit of CO to p = 0.04 in the insulating regime, suggesting that CO in insulating cuprates is incommensurate with the underlying lattice, contrary to some previous reports7,8.Fig. 2 Observation of charge order and the corresponding electronic structure in the insulating sample (p = 0.04).

a–c RIXS intensity maps with decreasing temperature, revealing the existence of CO in the quasi-elastic region centered around zero energy loss (indicated by the blue arrow). The gray circles in the lower panel of each figure represent the in-plane H-dependence of RIXS intensity, which is summed over the energy range between the two white dashed lines in the upper panel. The black curve and blue shaded areas are the total fits and the Gaussian component characterizing the CO peak, respectively. The CO is observed only along the (H, 0) or antinodal direction, and it is absent along the (H, H) or nodal direction (gray squares in c). d dI/dV map at V = 30 mV, displaying a very weak CO pattern along the Cu-O or antinodal directions (white arrows) at 77 K. e, Representative dI/dV spectra at points indicated in the inset, which is a zoomed-in view of the white dashed box in d. f, dI/dV map at V = 30 mV, displaying the CO pattern at 5 K. g The PG map at the same area as f, which displays the spatial distribution of the fitted PG magnitude (see Methods) from the tunneling spectra. h The spatial evolution of dI/dV spectrum along a line cutting through CO-poor and CO-rich regions is indicated by the red arrows in f and g, which reveals the particle-hole asymmetric PG by subtracting the CT gap background. The non-dispersive nature of the dI/dV maps in a wide energy range confirms the existence of static CO (see Supplementary Fig. 3 and Fig. 5). Scale bar: 5 nm.

From a real-space perspective, the checkerboard CO feature is particularly weak at a high temperature of 77 K as revealed by the dI/dV map (see Methods) in Fig. 2d. At this temperature, only a small portion of space with large local density-of-states (bright) holds charge ‘puddles’ that are poorly connected along the Cu-O bond directions (white arrows). This observation aligns with the broad full-width at half-maximum (FWHM) of the CO peak revealed by RIXS (e.g., 0.072 r.l.u. at 90 K, corresponding to a very short correlation length ξ ≈ 4.4 a0). A closer examination in Fig. 2e reveals that the dark region, where no CO develops, is characterized by a full charge-transfer gap. In contrast, the CO is accompanied by a density-of-states modulation, especially near the energy scale of the ‘V-shaped’ PG. Figure 2f, g are the low-temperature (T = 5 K) real-space patterns of the CO and PG magnitude (Methods) in the same area. It becomes immediately apparent that the CO is spatially correlated with smaller, or more precisely, sharper PG regions. The spectral linecut in Fig. 2h illustrates a smooth transition from the charge-transfer gap (which has been subtracted and shown as the red curve) in CO-poor locations to PG (CO-rich). Interestingly, the PG size is asymmetric at positive and negative bias, distinct from SC. This asymmetry is not unexpected, as a broken translational-symmetry phase, such as the CO here or the PG itself20,22, generally breaks particle-hole symmetry. Our result has, for the first time, demonstrated such asymmetry in both the occupied and unoccupied states.

To comprehend the evolution of a charge-ordered insulator into a high-temperature superconductor (along path #2 of Fig. 1b), we next move to the p = 0.06 sample with Tc = 20 K. The dI/dV map at 5 K (Fig. 3a) shows that the CO pattern is stronger than that in the insulator (Fig. 2f), which is also visualized by comparing the CO peak intensity and area from the RIXS data (Figs. 3c and 2c). This observation appears inconsistent with the notion that CO competes with SC, but rather plays a reinforced role. Remarkably, both STM (Fig. 3b) and RIXS (Fig. 3c) provide virtually identical QCO = 0.30 r.l.u. (≈3.3 a0 periodicity) along the (H, 0) directions, reconciling the charge modulations at the surface and in the bulk. The QCO of Bi2212 extended to p = 0.04 and p = 0.06 follows the extrapolated trend shown in Fig. 1a, in agreement with the behavior observed in most cuprates2.Fig. 3 Charge order and phonons in the superconducting sample (p = 0.06).

dI/dV map at V = 30 mV (a) and corresponding Fourier transformation (b), revealing a CO with QCO = 0.30 ± 0.01 r.l.u. along the antinodal direction (blue arrow). Scale bar: 5 nm. c The RIXS intensity map (upper) and integrated quasi-elastic CO peak (lower), with the color scale identical to that in Fig. 2. The red, blue, and black arrows indicate the characteristic energy scales of CO, bond buckling (BB) phonon and bond stretching (BS) phonon, respectively. d Joint fit of the RIXS spectrum at H = 0.12. e Subtracted spectra highlighting the BB (red shaded) and BS (black shaded) phonon modes. The black and red triangles indicate the fitted energies of BS and BB phonons, respectively.

Interplay between CO and phonons

Regarding the SC, it has been proposed that EPC plays a significant role in realizing high Tc, irrespective of the exact pairing mechanism23,24. However, the phonon modes are inaccessible by STM in the strongly underdoped regime, since the prominent PG prevents the extraction process, such as a strong-coupling Eliashberg approach25. In contrast, the RIXS cross-section offers the capability to simultaneously map the dispersion and EPC strength for specific phonon branches26, which also have an intricate interplay with CO16. Consequently, we conduct a more in-depth investigation to explore the potential connection between CO, EPC, and SC, utilizing the inelastic spectral weight of RIXS.

A previous study on antiferromagnetic NBCO demonstrated that Cu L3-edge RIXS can simultaneously detect the Cu-O bond stretching (BS) mode at approximately 70 meV and the out-of-plane A1g bond buckling (BB) phonon at about 35 meV27. When there is no charge order, the electron coupling to the BS phonon follows a sin2(πH) trend, while the BB phonon follows a cos2(πH) trend27, consistent with the model calculations of momentum-dependent EPC26. Figure 3d, e shows a raw RIXS spectrum and the inelastic parts after subtracting a resolution-limited elastic peak, respectively. The highly asymmetric lineshapes in the spectra indicate the necessity of considering two phonons in data fitting. With improved energy resolution and superb momentum accuracy in our experiment, we use two Gaussian functions (Methods) to extract the BS mode (shown in black) and the BB mode (in red). Our analyses show that both modes are nearly dispersionless within the energy resolution, which is distinct from previous one-phonon fitting models that presented a BS mode “softening”16–18. However, it is important to note that our present study is insufficient to rule out the softening behavior of lower-energy phonons. As energy resolution is further improved, softened phonons have been revealed by IXS on YBCO28 and higher-resolution RIXS in LSCO29.

Figure 4 a, b compare the decomposed CO (lower panels), BB phonon (middle panels), and BS phonon (upper panels) in the insulating and superconducting Dy-Bi2212 samples (Methods and Supplementary Fig. 8). The corresponding energy-integrated intensity profiles are shown in Fig. 4c. Intriguingly, the BB phonon exhibits a broad peak centered around 0.32 r.l.u., which is likely due to an interaction between the phonon and the dispersive CO excitation16. Both the CO and BB phonon intensities are higher in the superconducting sample (see also Fig. 4d). In contrast, the momentum-dependent EPC of the BS phonon lacks a peak structure and is well captured by a sin2(πH) profile within a non-interacting model27. Therefore, our new analysis suggests that the intensity anomaly observed in previous studies should be attributed to the BB phonon mode, which varies with material, doping, and temperature.Fig. 4 The doping evolution of RIXS spectra near the entrance of high-temperature superconductivity.

The RIXS intensity map highlighting each fitted excitation is shown for the p = 0.04 insulating (a) and p = 0.06 superconducting (b) samples for comparison. Upper panels are obtained by subtracting the fitted elastic and bond buckling (BB) phonon peaks from the raw data. Middle panels: subtracting the fitted elastic and bond stretching (BS) phonon peaks. Lower panels: subtracting the fitted BB and BS phonon peaks. White symbols superimposed on the maps indicate the fitted dispersion of the phonons, and the error bars represent the 95% confidence intervals of the fit. c, Comparison of RIXS intensity profiles integrated within the white dashed boxes in a and b. Note that a peak feature of the BS phonon branch reported at H = 0.36 r.l.u. on underdoped Dy-Bi2212 (ref. 16) is absent in our analysis. d Raw RIXS spectra (light-colored) at selective momenta indicated by the black dashed lines in c, highlighting the doping comparisons of different components (red and blue lines).

Doping- and T-dependence of CO and EPC

To gain more insights into what happens when the system enters the superconducting phase, we consider both the doping and systematic temperature evolutions following the three paths depicted in Fig. 1b. The short-range charge correlation already exists at 273 K (Fig. 5a), which is accompanied by a very broad BB phonon intensity peak (Fig. 5b). Therefore, the weak CO and its interplay with EPC set the stage for the formation and condensation of Cooper pairs. Upon cooling along paths #1 and #3, both CO and BB phonon gradually gain more weight, but they follow very different trends (see Fig. 5c,d): the CO intensity grows faster in the superconducting sample than in the insulating one. In contrast, the BB phonon displays a concurrent and smooth increase, except for an abrupt jump (paths #2 and #3) when the SC appears. These results agree with previous reports on the separation of CO intensity and excitations6,17. The intensity profile of the BS phonon, on the other hand, maintains a sinusoidal trend at different dopings, temperatures (Fig. 5e), and materials (refs. 26,27), indicating a negligible correlation with CO excitations and SC.Fig. 5 The temperature and doping evolutions of CO and EPC.

a CO peak profiles obtained by integrating the quasi-elastic maps (Supplementary Fig. 8). The curves are fit to the peaks (Methods). b Normalized phonon intensity profiles of the p = 0.04 (upper panel) and p = 0.06 (lower panel) samples obtained by integrating the BB phonon maps (Supplementary Fig. 8 and Fig. 10). The colored curve is the fit to each profile using a Gaussian function plus a linear background. c, Elastic peak area in a retrieved by the fits. The error bars represent the 95% confidence intervals. d, BB phonon peak area obtained by integrating the spectral weight within the black dashed boxes in b. e, BS phonon intensity profiles at the two Bi2212 samples and all temperatures investigated, with a sinusoidal line (gray shaded) overlaid. The predicted sin2(πH) trend agrees with and substantially extends a recent RIXS result on La2−xSrxCuO4 cuprates18.

Figure 6 delineates our idea of how the SC emerges below Tc. Our findings reveal that while the SC emerges from the CO phase (along paths #2 and #3 in Fig. 1b), the CO wavevector QCO continuously changes with doping, irrespective of the SC (Fig. 6a). In concert with this observation, previous studies have noted that the CO intensity could respond differently to SC3,30,31. It is the BB phonon intensity anomaly, or the dispersive CO excitations6 triggered by the underlying CO, that induces an abrupt enhancement in the SC state (Fig. 6b). The enhanced EPC strength of the BB phonon may drive or accompany the condensation of Cooper pairs, as evidenced by the coherent superconducting gap everywhere (Fig. 6c). Therefore, the CO and EPC play distinct roles but work together in establishing the macroscopic SC.Fig. 6 Intimate correlation between CO, EPC and SC.

a Summarized plot of QCO as a function of hole concentration, displaying a monotonously decreasing trend with increasing doping. The black symbols are the same as Fig. 1a, adapted from refs. 5,6,16. Note that this decreasing trend may only be at work below p = 0.1653,54. b In the superconducting Dy-Bi2212 (p = 0.06), the dispersive excitations (white dash-dotted) of the short-range CO cross the BB phonon dispersion and induce a strong peak at this branch. The intensity map is obtained by first subtracting the BS phonon mode from Fig. 3c and is then smoothed. c A superconducting gap of  ~ 20 meV is found in different locations on the p = 0.06 sample, indicated by different colors (crosses in Fig. 3a). The shaded blue and red arrows indicate the intimate connections between RIXS and STM experiments. On one hand, the CO branch revealed by the quasi-elastic peak of RIXS gives the same QCO and its doping evolution as STM. On the other hand, the inelastic BB phonon strength abruptly increases once the Cooper pairs condense, corresponding to the SC gap revealed by STM.

Discussion

The combined RIXS and STM results here touch on the important issues regarding the relations between PG, CO, EPC, and SC, all of which are at play in the present study. While a short-range CO has been found pervading the hole-doped phase diagram, where SC acts as the ground state, its extension to the insulating regime remains elusive, with only two STM works (refs. 7 and 8) reporting surface CO. Our work, for the first time, has not only provided concrete evidence of CO in an insulating cuprate through a combined bulk and surface study but also substantially extended the lower doping limit of CO to p = 0.04 with a broad temperature range. Our observations reveal an incommensurate CO wave vector surrounding  ~ 0.31 r.l.u. However, our RIXS results and traditional Fourier-transformed dI/dV (proportional to the local density-of-states) channel cannot differentiate between truly incommensurate charge order and discommensurations within a locally commensurate phase resulting from disorder or topological defects29,32,33. Therefore, the true nature of the CO, whether it is genuinely incommensurate or merely appears to be incommensurate due to certain factors, remains a topic of controversy.

Another hotly debated issue pertains to the origin of CO. In a weak coupling k-space Fermi-surface nesting scenario, the CO wave vector would continuously decrease with increased doping4,10. Alternatively, in a strong coupling r-space scenario, the charge order would be commensurate with the underlying lattice. While we have noted a continuous decrease in QCO with doping, the nesting scenario is not applicable in the heavily underdoped regime where the Fermi surface is entirely gapped13. However, our findings can be reconciled through discommensuration in the r-space strong correlation scenario, where smaller random commensurate domains can cause a larger mismatch of commensurate wave vectors at low doping, as proposed by phase-resolved STM studies32. One may also suspect if this incipient CO is caused by antiferromagnetic excitations in the insulating regime, as suggested in other hole-doped cuprates34,35. However, our RIXS results reveal a smooth paramagnon dispersion and damping across QCO (Fig. 1c and Supplementary Fig. 11), showing an absence of intertwining between spin excitations and CO.

We move forward to the initial question of how the unconventional SC emerges. In the lightly doped cuprates, we are close to a non-Fermi-liquid Bose-Einstein condensate regime. Although our experiments do not give direct evidence on incoherent Cooper pairs, we do suggest that CO is indispensable for their formation. Signatures of “pre-formed” Cooper pairs have been widely observed by STM36,37, ARPES38, and transport measurements39,40. The coexistence of pre-formed Cooper pairs and CO outside the superconducting dome suggests that the charge-ordered electronic state is energetically favorable for the formation of Cooper pairs, akin to the description provided by the stripe model41. This agrees with the observation that two holes are distributed on the peak of the checkerboard plaquette, as visualized by a periodical superconducting coherence peak modulation atop the CO37,42.

In cuprates, another intriguing correspondence is between the short-range nature of CO and nanoscale electronic inhomogeneity43, indicative of the weakly-overlapped Cooper pairs in a real-space picture. The CO, though still short-ranged, gains more amplitude with cooling (along paths #1 and #3) or doping (along path #2). Working in conjunction with the CO, it is the enhanced EPC strength that accompanies the eventual condensation of Cooper pairs (Fig. 6b, c). Remarkably, the BS phonon is invariable against temperature and doping (Fig. 5e), agreeing well with the model calculations26, RIXS measurements in Nd1+xBa2−xCu3O7−δ27, and ARPES measurements in Bi221244. In contrast, the coupling to the BB phonon is a more relevant process, which has been proposed to promote the magnitude of Tc when the long-range Coulomb interaction is plausible24. Similar observations sensitive to the Cooper pairs condensation include the antinodal kink of the buckling B1g phonon observed by ARPES45, and the spin resonance mode observed by inelastic neutron scattering46. Whether these are the driving forces or natural consequences of the itinerant Cooper pairs remains an open question.

Methods

Sample growth and characterizations

High-quality Bi2Sr2Ca0.6Dy0.4Cu2O8+δ single crystals were grown by the traveling solvent floating zone method. The high Dy concentration was intentionally introduced during the growth to achieve the highly underdoped regime. The as-grown sample displays a superconducting transition with Tconset = 40 K as revealed by the magnetic susceptibility curve, corresponding to a hole doping level of p = 0.067 (0.067 ± 0.008). The doping and error bar were estimated using the empirical formula from ref. 47, where the Tcmax was set to 96 K for the Bi2212 system. One as-grown sample was divided into two crystals, which were used in the present study after different post-annealing processes in a vacuum atmosphere. Notably, the STM measurements were carried out in the very same samples after the RIXS experiments. The first annealing condition (Ann1) was T = 350 ∘C for 2 days, which yields a superconducting transition with Tc ~ 20 K and p = 0.062 (0.062 ± 0.006). The second one (Ann2) was T = 450 ∘C for 3 days, which completely suppressed SC showing no evidence of the Meissner effect. We estimate the doping level of Ann2 to be p ~ 0.04 (0.04 ± 0.01) for two reasons: (i) It is below the onset doping of the superconducting dome (p ~ 0.05) and is highly insulating as revealed by the resistivity measurement, and (ii) Previous reports have proposed that it is impossible to carry out liquid-helium-temperature STM measurements on cuprates at or below p = 0.037, even with our delicate methodology (see below). Therefore, the sample prepared by Ann2 is located somewhere below but in the vicinity of the SC onset. The magnetization and resistivity curves are shown in Supplementary Fig. 1.

RIXS measurements

The RIXS experiments were performed at the I21 RIXS beamline of the Diamond Light Source19. The momentum resolution is estimated to be 0.01 Å−1 based on the scattering geometry employed in our study19. For Bi2212, the lattice parameters are a = b = 3.83 Å, corresponding to ΔQ ~ 0.006 r.l.u., which is sufficient for the CO analysis. The Dy-Bi2212 samples were cleaved to expose a fresh surface and then immediately transferred into the ultrahigh vacuum (UHV) chamber. Before the RIXS measurements, we did careful alignments for our samples utilizing the (002) Bragg peak and the (1-10) superstructure peak, which are critical for an accurate determination of the scattering vector q. All RIXS spectra in our study were collected with the incident photon energy tuned to the Cu L3 resonant edge (see Supplementary Fig. 1), where we can achieve an ultra-high energy resolution ΔE ~ 36 meV. To enhance the charge-related excitations, such as CO and phonon branches, all data, except for those shown in Supplementary Fig. 9, were taken with the vertical (σ) polarization of the incident X-ray. To maximize the in-plane momentum-transfer H along the (H, 0) direction, a large scattering angle 2θ = 154° was chosen.

STM measurements

The STM experiments were performed on the same crystals immediately after their return from the RIXS measurements at the Diamond Light Source. We used a commercial (CreaTec) STM system working at low-temperature and UHV environments. Each Dy-Bi2212 single crystal was cleaved in situ at T = 77 K and then immediately transferred into the STM stage stabilized at 77 K or 5 K.

To carry out STM measurements in insulators and to obtain high-quality tunneling spectra over a large bias range (e.g. Fig. 1d), it is crucial to have a stable metallic tip. In our experiment, an electrochemically etched tungsten tip was transferred into the UHV chamber and was then electron-beam bombarded for 15 min to remove the oxides. Before the measurements, the tip was carefully treated and calibrated on an Au(111) single crystal. The basic requirements include: (i) an atomically sharp apex, (ii) metallic DOS which reveals a well-defined surface state and quasi-particle interference pattern of the Au(111), and (iii) stable enough under a point-contact current as high as 10 uA between the tip and Au(111). The tip was then slowly approached onto the cleaved Dy-Bi2212 sample surface with a large magnitude of negative bias voltage, since the occupied state has a much larger DOS due to the tunneling asymmetry seen in CT gap and PG (e.g. Fig. 1d). More details of this methodology can be found in our previous paper48.

Set-point conditions for Fig. 1: tunneling current I = 30 pA, bias voltage V0 = − 1.0 V, ac − modulation Vm = 30 meV. To mitigate the effect of V0 on the low-energy density-of-states, which is known as the set-point effect, a higher V0 is utilized to reveal the CO signal. In our experiment, set-point conditions for Fig. 2 (the p = 0.04 sample): I = 100 pA, V = − 0.5 V, Vm = 6 mV. Set-point conditions for Fig. 3 and Fig. 5 (the p = 0.06 sample): I = 150 pA, V = − 0.3 V, Vm = 4 mV.

RIXS Data analysis

We first determined the zero-energy-loss of each RIXS spectrum by comparing it with that recorded on a nearby carbon tape with the same scattering geometry. After that, the data were normalized to the incident beamline flux and then corrected for the self-absorption effect using the procedures described in ref. 49. The data in Supplementary Fig. 8 were normalized to the areas of the dd-excitations.

The quasi-elastic CO peak profiles shown in Figs. 2a–c, 3c, and 5a were fitted to a Gaussian peak with a Lorentzian background centered at H = 0. From this analysis, the CO wavevector (e.g., Fig. 6a) and amplitude (e.g., Fig. 5c) were respectively characterized by the position and area of the Gaussian component.

The inelastic phonon branches were extracted through a fitting model to each RIXS spectrum, as exemplified in Fig. 3d. It involves resolution-restricted Gaussians for the quasi-elastic, BB phonon, and BS phonon peaks, and a higher-energy background of an anti-symmetric Lorentzian convoluted with the energy resolution. From this analysis, we decomposed the raw data to CO, BB phonon, and BS phonon intensity maps shown in Fig. 4 and Supplementary Fig. 8, which further generate Fig. 5a,b. The adequacy of our model and the comparison with other one-phonon fitting models are demonstrated in Supplementary Fig. 7 and discussed in the accompanying text.

STM data analysis

For characterizing the bi-directional checkerboard CO, the dI/dV map is widely employed for STM analysis4,7,8,10,21,37, which provides a direct comparison with the resonant X-ray scattering signal4,10. From the STM point of view, the defining feature of the CO is an in-phase density-of-states distribution between the particle (-E) and hole (+E) state, distinguishing it from other electronic modulations that are studied by the “Z-map"20,32,36,50,51. In practice, the “Z-map" has been effective in mitigating the CO signal for QPI and uni-directional modulation studies.

The PG map ΔPG(r) shown in Fig. 2g was obtained by fitting the tunneling spectrum at each pixel dI/dV(r, E). The fit was done on the positive-E side (unoccupied state) with an evident peak, where a s-wave Dynes formula works pretty well for the broad feature in this highly underdoped regime50. We note that the density-of-states, or tunneling rate, asymmetry seen in the PG lineshape is generic for cuprates20, which is due to the inequivalence between electron injection and extraction and reflects, thereby, the strong electron-electron correlation effect52.

The STM momentum resolution is estimated to be ΔQ = 2π/23.8 nm = 0.026 Å−1 ~ 0.016 r.l.u. (where 23.8 nm represents the real-space size of Fig. 3a), which is sufficient for the CO analysis. To directly compare the QCO between Fourier-transformed STM and RIXS (Fig. 3), we have corrected the thermal and piezoelectric drifts of the p = 0.06 dataset. The methodology used for this correction, Lawler-Fujita’s algorithm, has been well documented in ref. 51.

Supplementary information

Supplementary Information

Peer Review File

Source data

Source data

Supplementary information

The online version contains supplementary material available at 10.1038/s41467-024-52124-9.

Acknowledgements

We thank J. Wang, Y. Li, Z. Y. Wang, S. Johnston and Z. Y. Weng for their helpful discussions. Y. Peng is grateful for the financial support from the Ministry of Science and Technology of China (Grants No. 2021YFA1401903 and No. 2019YFA0308401), the Beijing Natural Science Foundation (Grant No. JQ24001), and the National Natural Science Foundation of China (Grants No. 12374143 and 11974029). Y. Wang acknowledges the support by the MOST of China grant No. 2017YFA0302900. C. Zou benefits from the Project funded by China Postdoctoral Science Foundation (Grant No. BX2021019). X. J. Z. is supported by the National Natural Science Foundation of China (Grant No. 11888101). We acknowledge Diamond Light Source for providing the beamtime at the I21-RIXS beamline under Proposal MM28855.

Author contributions

C.Z. and Y.P. proposed and designed the research. C.Z., J.C., K.J.Z. and Y.P. carried out the RIXS experiments with the help of Q.L., M.G.F., S.A., Q.Q., X.C., and Q.X., C.Z., S.Y. and Y.W. carried out the STM experiments. C.Z. and Q.L. analyzed the data. C.Y. and X.Z. provided the Dy-Bi2212 single crystals. C.Z. performed the vacuum annealing. C.Z. and S.Y. characterized the samples. C.Z. and Y.P. prepared the manuscript. All authors have read and approved the final version of the manuscript.

Peer review

Peer review information

Nature Communications thanks Marc Janoschek, and the other, anonymous, reviewer(s) for their contribution to the peer review of this work. A peer review file is available.

Data availability

The data used to support the findings of this work are available from the corresponding author upon request. Source data are provided with this paper.

Competing interests

The authors declare no competing interests.

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
==== Refs
References

1. Keimer B Kivelson SA Norman MR Uchida S Zaanen J From quantum matter to high-temperature superconductivity in copper oxides Nature 2015 518 179 186 10.1038/nature14165 25673411
Keimer, B., Kivelson, S. A., Norman, M. R., Uchida, S. & Zaanen, J. From quantum matter to high-temperature superconductivity in copper oxides. Nature 518, 179–186 (2015).25673411 10.1038/nature14165
2. Comin R Damascelli A Resonant X-ray scattering studies of charge order in cuprates Annu. Rev. Condens. Matter Phys. 2016 7 369 405 10.1146/annurev-conmatphys-031115-011401
Comin, R. & Damascelli, A. Resonant X-ray scattering studies of charge order in cuprates. Annu. Rev. Condens. Matter Phys. 7, 369–405 (2016).10.1146/annurev-conmatphys-031115-011401
3. Ghiringhelli G Long-range incommensurate charge fluctuations in (Y,Nd)Ba2Cu3O6+x Science 2012 337 821 825 10.1126/science.1223532 22798406
Ghiringhelli, G. et al. Long-range incommensurate charge fluctuations in (Y,Nd)Ba2Cu3O6+x. Science 337, 821–825 (2012).22798406 10.1126/science.1223532
4. Neto EHDS Ubiquitous interplay between charge ordering and high-temperature superconductivity in cuprates Science 2014 343 393 396 10.1126/science.1243479 24356110
Neto, E. H. D. S. et al. Ubiquitous interplay between charge ordering and high-temperature superconductivity in cuprates. Science 343, 393–396 (2014).24356110 10.1126/science.1243479
5. Boschini F Dynamic electron correlations with charge order wavelength along all directions in the copper oxide plane Nat. Commun. 2021 12 1 8 10.1038/s41467-020-20824-7 33397941
Boschini, F. et al. Dynamic electron correlations with charge order wavelength along all directions in the copper oxide plane. Nat. Commun. 12, 1–8 (2021).33397941 10.1038/s41467-020-20824-7
6. Lee WS Spectroscopic fingerprint of charge order melting driven by quantum fluctuations in a cuprate Nat. Phys. 2021 17 53 57 10.1038/s41567-020-0993-7
Lee, W. S. et al. Spectroscopic fingerprint of charge order melting driven by quantum fluctuations in a cuprate. Nat. Phys. 17, 53–57 (2021).10.1038/s41567-020-0993-7
7. Cai P Visualizing the evolution from the Mott insulator to a charge-ordered insulator in lightly doped cuprates Nat. Phys. 2016 12 1047 1051 10.1038/nphys3840
Cai, P. et al. Visualizing the evolution from the Mott insulator to a charge-ordered insulator in lightly doped cuprates. Nat. Phys. 12, 1047–1051 (2016).10.1038/nphys3840
8. Hanaguri T A ‘checkerboard’ electronic crystal state in lightly hole-doped Ca2−xNaxCuO2Cl2 Nature 2004 430 1001 1005 10.1038/nature02861 15329714
Hanaguri, T. et al. A ‘checkerboard’ electronic crystal state in lightly hole-doped Ca2−xNaxCuO2Cl2. Nature 430, 1001–1005 (2004).15329714 10.1038/nature02861
9. Damascelli A Hussain Z Shen Z-X Angle-resolved photoemission studies of the cuprate superconductors Rev. Mod. Phys. 2003 75 473 541 10.1103/RevModPhys.75.473
Damascelli, A., Hussain, Z. & Shen, Z.-X. Angle-resolved photoemission studies of the cuprate superconductors. Rev. Mod. Phys. 75, 473–541 (2003).10.1103/RevModPhys.75.473
10. Comin R Charge order driven by fermi-arc instability in Bi2Sr2−xLaxCuO6+δ Science 2014 343 390 392 10.1126/science.1242996 24356115
Comin, R. et al. Charge order driven by fermi-arc instability in Bi2Sr2−xLaxCuO6+δ. Science 343, 390–392 (2014).24356115 10.1126/science.1242996
11. Doiron-Leyraud N Quantum oscillations and the fermi surface in an underdoped high-Tc superconductor Nature 2007 447 565 568 10.1038/nature05872 17538614
Doiron-Leyraud, N. et al. Quantum oscillations and the fermi surface in an underdoped high-T c superconductor. Nature 447, 565–568 (2007).17538614 10.1038/nature05872
12. Fang, Y. et al. Fermi surface transformation at the pseudogap critical point of a cuprate superconductor. Nat. Phys. 10.1038/s41567-022-01514-1 (2022).
13. Vishik I Phase competition in trisected superconducting dome Proc. Natl. Acad. Sci. USA 2012 109 18332 18337 10.1073/pnas.1209471109 23093670
Vishik, I. et al. Phase competition in trisected superconducting dome. Proc. Natl. Acad. Sci. USA 109, 18332–18337 (2012).23093670 10.1073/pnas.1209471109
14. Peng Y Disappearance of nodal gap across the insulator-superconductor transition in a copper-oxide superconductor Nat. Commun. 2013 4 2459 10.1038/ncomms3459 24051514
Peng, Y. et al. Disappearance of nodal gap across the insulator-superconductor transition in a copper-oxide superconductor. Nat. Commun. 4, 2459 (2013).24051514 10.1038/ncomms3459
15. Ruan W Relationship between the parent charge transfer gap and maximum transition temperature in cuprates Sci. Bull. 2016 61 1826 1832 10.1007/s11434-016-1204-x
Ruan, W. et al. Relationship between the parent charge transfer gap and maximum transition temperature in cuprates. Sci. Bull. 61, 1826–1832 (2016).10.1007/s11434-016-1204-x
16. Chaix L Dispersive charge density wave excitations in Bi2Sr2CaCu2O8+δ Nat. Phys. 2017 13 952 956 10.1038/nphys4157
Chaix, L. et al. Dispersive charge density wave excitations in Bi2Sr2CaCu2O8+δ. Nat. Phys. 13, 952–956 (2017).10.1038/nphys4157
17. Wang Q Charge order lock-in by electron-phonon coupling in La1.675Eu0.2Sr0.125CuO4 Sci. Adv. 2021 7 27
Wang, Q. et al. Charge order lock-in by electron-phonon coupling in La1.675Eu0.2Sr0.125CuO4. Sci. Adv. 7, 27 (2021).
18. Lin JQ Strongly correlated charge density wave in La2−xSrxCuO4 evidenced by doping-dependent phonon anomaly Phys. Rev. Lett. 2020 124 207005 10.1103/PhysRevLett.124.207005 32501068
Lin, J. Q. et al. Strongly correlated charge density wave in La2−xSrxCuO4 evidenced by doping-dependent phonon anomaly. Phys. Rev. Lett. 124, 207005 (2020).32501068 10.1103/PhysRevLett.124.207005
19. Zhou KJ I21: an advanced high-resolution resonant inelastic X-ray scattering beamline at Diamond Light Source J. Synchrotron Radiat. 2022 29 563 580 10.1107/S1600577522000601 35254322
Zhou, K. J. et al. I21: an advanced high-resolution resonant inelastic X-ray scattering beamline at Diamond Light Source. J. Synchrotron Radiat. 29, 563–580 (2022).35254322 10.1107/S1600577522000601
20. Kohsaka Y An intrinsic bond-centered electronic glass with unidirectional domains in underdoped cuprates Science 2007 315 1380 1385 10.1126/science.1138584 17289939
Kohsaka, Y. et al. An intrinsic bond-centered electronic glass with unidirectional domains in underdoped cuprates. Science 315, 1380–1385 (2007).17289939 10.1126/science.1138584
21. Yu Y Ma L Cai P Zhong R Zhang Y High-temperature superconductivity in monolayer Bi2Sr2CaCu2O8+δ Nature 2019 575 156 163 10.1038/s41586-019-1718-x 31666697
Yu, Y., Ma, L., Cai, P., Zhong, R. & Zhang, Y. High-temperature superconductivity in monolayer Bi2Sr2CaCu2O8+δ. Nature 575, 156–163 (2019).31666697 10.1038/s41586-019-1718-x
22. Hashimoto M Particle-hole symmetry breaking in the pseudogap state of Bi2201 Nat. Phys. 2010 6 414 418 10.1038/nphys1632
Hashimoto, M. et al. Particle-hole symmetry breaking in the pseudogap state of Bi2201. Nat. Phys. 6, 414–418 (2010).10.1038/nphys1632
23. Sakai T Poilblanc D Scalapino DJ Hole pairing and phonon dynamics in generalized two-dimensional t-J Holstein models Phys. Rev. B 1997 55 8445 8451 10.1103/PhysRevB.55.8445
Sakai, T., Poilblanc, D. & Scalapino, D. J. Hole pairing and phonon dynamics in generalized two-dimensional t-J Holstein models. Phys. Rev. B 55, 8445–8451 (1997).10.1103/PhysRevB.55.8445
24. Johnston S Systematic study of electron-phonon coupling to oxygen modes across the cuprates Phys. Rev. B 2010 82 064513 10.1103/PhysRevB.82.064513
Johnston, S. et al. Systematic study of electron-phonon coupling to oxygen modes across the cuprates. Phys. Rev. B 82, 064513 (2010).10.1103/PhysRevB.82.064513
25. Levy de Castro G Berthod C Piriou A Giannini E Fischer O Preeminent role of the van Hove singularity in the strong-coupling analysis of scanning tunneling spectroscopy for two-dimensional cuprate superconductors Phys. Rev. Lett. 2008 101 267004 10.1103/PhysRevLett.101.267004 19437664
Levy de Castro, G., Berthod, C., Piriou, A., Giannini, E. & Fischer, O. Preeminent role of the van Hove singularity in the strong-coupling analysis of scanning tunneling spectroscopy for two-dimensional cuprate superconductors. Phys. Rev. Lett. 101, 267004 (2008).19437664 10.1103/PhysRevLett.101.267004
26. Devereaux TP Directly characterizing the relative strength and momentum dependence of electron-phonon coupling using resonant inelastic X-ray scattering Phys. Rev. X 2016 6 041019
Devereaux, T. P. et al. Directly characterizing the relative strength and momentum dependence of electron-phonon coupling using resonant inelastic X-ray scattering. Phys. Rev. X 6, 041019 (2016).
27. Braicovich L Determining the electron-phonon coupling in superconducting cuprates by resonant inelastic x-ray scattering: Methods and results on Nd1+xBa2−xCu3O7−δ Phys. Rev. Res. 2020 2 023231 10.1103/PhysRevResearch.2.023231
Braicovich, L. et al. Determining the electron-phonon coupling in superconducting cuprates by resonant inelastic x-ray scattering: Methods and results on Nd1+xBa2−xCu3O7−δ. Phys. Rev. Res. 2, 023231 (2020).10.1103/PhysRevResearch.2.023231
28. Tacon ML Bosak A Souliou SM Dellea G Keimer B Inelastic X-ray scattering in YBa2Cu3O6.6 reveals giant phonon anomalies and elastic central peak due to charge-density-wave formation Nat. Phys. 2013 10 52 58 10.1038/nphys2805
Tacon, M. L., Bosak, A., Souliou, S. M., Dellea, G. & Keimer, B. Inelastic X-ray scattering in YBa2Cu3O6.6 reveals giant phonon anomalies and elastic central peak due to charge-density-wave formation. Nat. Phys. 10, 52–58 (2013).10.1038/nphys2805
29. Huang HY Quantum fluctuations of charge order induce phonon softening in a superconducting cuprate Phys. Rev. X 2021 11 041038
Huang, H. Y. et al. Quantum fluctuations of charge order induce phonon softening in a superconducting cuprate. Phys. Rev. X 11, 041038 (2021).
30. Miao H Formation of incommensurate charge density waves in cuprates Phys. Rev. X 2019 9 031042
Miao, H. et al. Formation of incommensurate charge density waves in cuprates. Phys. Rev. X 9, 031042 (2019).
31. Wen JJ Huang H Lee SJ Jang H Lee JS Observation of two types of charge-density-wave orders in superconducting La2−xSrxCuO4 Nat. Commun. 2019 10 3269 10.1038/s41467-019-11167-z 31332190
Wen, J. J., Huang, H., Lee, S. J., Jang, H. & Lee, J. S. Observation of two types of charge-density-wave orders in superconducting La2−xSrxCuO4. Nat. Commun. 10, 3269 (2019).31332190 10.1038/s41467-019-11167-z
32. Mesaros A Commensurate 4a0-period charge density modulations throughout the Bi2Sr2CaCu2O8+x pseudogap regime Proc. Natl. Acad. Sci. USA 2016 113 12661 12666 10.1073/pnas.1614247113 27791157
Mesaros, A. et al. Commensurate 4a0-period charge density modulations throughout the Bi2Sr2CaCu2O8+x pseudogap regime. Proc. Natl. Acad. Sci. USA 113, 12661–12666 (2016).27791157 10.1073/pnas.1614247113
33. Vinograd, I., Zhou, R., Hirata, M., Wu, T. & Julien, M. H. Locally commensurate charge-density wave with three-unit-cell periodicity in YBa2Cu3Oy. Nat. Commun. 12 (2021).
34. Miao H High-temperature charge density wave correlations in La1.875Ba0.125CuO4 without spin-charge locking Proc. Natl. Acad. Sci. USA 2017 114 12430 12435 10.1073/pnas.1708549114 29114049
Miao, H. et al. High-temperature charge density wave correlations in La1.875Ba0.125CuO4 without spin-charge locking. Proc. Natl. Acad. Sci. USA 114, 12430–12435 (2017).29114049 10.1073/pnas.1708549114
35. Peng YY Magnetic excitations and phonons simultaneously studied by resonant inelastic x-ray scattering in optimally doped Bi1.5Pb0.55Sr1.6La0.4CuO6+δ Phys. Rev. B 2015 92 064517 10.1103/PhysRevB.92.064517
Peng, Y. Y. et al. Magnetic excitations and phonons simultaneously studied by resonant inelastic x-ray scattering in optimally doped Bi1.5Pb0.55Sr1.6La0.4CuO6+δ. Phys. Rev. B 92, 064517 (2015).10.1103/PhysRevB.92.064517
36. Lee J Spectroscopic fingerprint of phase-incoherent superconductivity in the underdoped Bi2Sr2CaCu2O8+δ Science 2009 325 1099 1103 10.1126/science.1176369 19713522
Lee, J. et al. Spectroscopic fingerprint of phase-incoherent superconductivity in the underdoped Bi2Sr2CaCu2O8+δ. Science 325, 1099–1103 (2009).19713522 10.1126/science.1176369
37. Ye, S. et al. The emergence of global phase coherence from local pairing in underdoped cuprates. Nat. Phys. 10.1038/s41567-023-02100-9 (2023).
38. He Y Superconducting fluctuations in overdoped Bi2Sr2CaCu2O8+δ Phys. Rev. X 2021 11 031068
He, Y. et al. Superconducting fluctuations in overdoped Bi2Sr2CaCu2O8+δ. Phys. Rev. X 11, 031068 (2021).
39. Wang Y Field-enhanced diamagnetism in the pseudogap state of the cuprate Bi2Sr2CaCu2O8+δ superconductor in an intense magnetic field Phys. Rev. Lett. 2005 95 247002 10.1103/PhysRevLett.95.247002 16384409
Wang, Y. et al. Field-enhanced diamagnetism in the pseudogap state of the cuprate Bi2Sr2CaCu2O8+δ superconductor in an intense magnetic field. Phys. Rev. Lett. 95, 247002 (2005).16384409 10.1103/PhysRevLett.95.247002
40. Tallon JL Storey JG Loram JW Fluctuations and critical temperature reduction in cuprate superconductors Phys. Rev. B 2011 83 092502 10.1103/PhysRevB.83.092502
Tallon, J. L., Storey, J. G. & Loram, J. W. Fluctuations and critical temperature reduction in cuprate superconductors. Phys. Rev. B 83, 092502 (2011).10.1103/PhysRevB.83.092502
41. Tranquada JM Dean MPM Li Q Superconductivity from charge order in cuprates J. Phys. Soc. Jpn. 2021 90 111002 10.7566/JPSJ.90.111002
Tranquada, J. M., Dean, M. P. M. & Li, Q. Superconductivity from charge order in cuprates. J. Phys. Soc. Jpn. 90, 111002 (2021).10.7566/JPSJ.90.111002
42. Ruan W Visualization of the periodic modulation of Cooper pairing in a cuprate superconductor Nat. Phys. 2018 14 1178 1182 10.1038/s41567-018-0276-8
Ruan, W. et al. Visualization of the periodic modulation of Cooper pairing in a cuprate superconductor. Nat. Phys. 14, 1178–1182 (2018).10.1038/s41567-018-0276-8
43. Zeljkovic I Imaging the impact of single oxygen atoms on superconducting Bi2+ySr2−yCaCu2O8+x Science 2012 337 320 323 10.1126/science.1218648 22822144
Zeljkovic, I. et al. Imaging the impact of single oxygen atoms on superconducting Bi2+ySr2−yCaCu2O8+x. Science 337, 320–323 (2012).22822144 10.1126/science.1218648
44. Lanzara A Evidence for ubiquitous strong electron-phonon coupling in high-temperature superconductors Nature 2001 412 510 514 10.1038/35087518 11484045
Lanzara, A. et al. Evidence for ubiquitous strong electron-phonon coupling in high-temperature superconductors. Nature 412, 510–514 (2001).11484045 10.1038/35087518
45. Cuk T Coupling of the B1g Phonon to the antinodal electronic states of Bi2Sr2Ca0.92Y0.08Cu2O8+δ Phys. Rev. Lett. 2004 93 117003 10.1103/PhysRevLett.93.117003 15447370
Cuk, T. et al. Coupling of the B 1g Phonon to the antinodal electronic states of Bi2Sr2Ca0.92Y0.08Cu2O8+δ. Phys. Rev. Lett. 93, 117003 (2004).15447370 10.1103/PhysRevLett.93.117003
46. Liu DZ Zha Y Levin K Theory of neutron scattering in the normal and superconducting states of YBa2Cu3O6+x Phys. Rev. Lett. 1995 75 4130 4133 10.1103/PhysRevLett.75.4130 10059822
Liu, D. Z., Zha, Y. & Levin, K. Theory of neutron scattering in the normal and superconducting states of YBa2Cu3O6+x. Phys. Rev. Lett. 75, 4130–4133 (1995).10059822 10.1103/PhysRevLett.75.4130
47. Presland M Tallon J Buckley R Liu R Flower N General trends in oxygen stoichiometry effects on Tc in Bi and Tl superconductors Phys. C 1991 176 95 105 10.1016/0921-4534(91)90700-9
Presland, M., Tallon, J., Buckley, R., Liu, R. & Flower, N. General trends in oxygen stoichiometry effects on Tc in Bi and Tl superconductors. Phys. C 176, 95–105 (1991).10.1016/0921-4534(91)90700-9
48. Ye C Visualizing the atomic-scale electronic structure of the Ca2CuO2Cl2 Mott insulator Nat. Commun. 2013 4 1365 10.1038/ncomms2369 23340405
Ye, C. et al. Visualizing the atomic-scale electronic structure of the Ca2CuO2Cl2 Mott insulator. Nat. Commun. 4, 1365 (2013).23340405 10.1038/ncomms2369
49. Minola M Collective nature of spin excitations in superconducting cuprates probed by resonant inelastic X-Ray scattering Phys. Rev. Lett. 2015 114 217003 10.1103/PhysRevLett.114.217003 26066453
Minola, M. et al. Collective nature of spin excitations in superconducting cuprates probed by resonant inelastic X-Ray scattering. Phys. Rev. Lett. 114, 217003 (2015).26066453 10.1103/PhysRevLett.114.217003
50. Kohsaka Y Visualizing the emergence of the pseudogap state and the evolution to superconductivity in a lightly hole-doped mott insulator Nat. Phys. 2012 8 604 10 10.1038/nphys2321
Kohsaka, Y. et al. Visualizing the emergence of the pseudogap state and the evolution to superconductivity in a lightly hole-doped mott insulator. Nat. Phys. 8, 604–10 (2012).10.1038/nphys2321
51. Lawler MJ Intra-unit-cell electronic nematicity of the high-Tc copper-oxide pseudogap states Nature 2010 466 347 351 10.1038/nature09169 20631795
Lawler, M. J. et al. Intra-unit-cell electronic nematicity of the high-T c copper-oxide pseudogap states. Nature 466, 347–351 (2010).20631795 10.1038/nature09169
52. Anderson P Ong N Theory of asymmetric tunneling in the cuprate superconductors J. Phys. Chem. Solids 2006 67 1 5 10.1016/j.jpcs.2005.10.132
Anderson, P. & Ong, N. Theory of asymmetric tunneling in the cuprate superconductors. J. Phys. Chem. Solids 67, 1–5 (2006).10.1016/j.jpcs.2005.10.132
53. Lu H Identification of a characteristic doping for charge order phenomena in Bi-2212 cuprates via RIXS Phys. Rev. B 2022 106 155109 10.1103/PhysRevB.106.155109
Lu, H. et al. Identification of a characteristic doping for charge order phenomena in Bi-2212 cuprates via RIXS. Phys. Rev. B 106, 155109 (2022).10.1103/PhysRevB.106.155109
54. Du Z Periodic atomic displacements and visualization of the electron-lattice interaction in the cuprate Phys. Rev. X 2023 13 021025
Du, Z. et al. Periodic atomic displacements and visualization of the electron-lattice interaction in the cuprate. Phys. Rev. X 13, 021025 (2023).
