Abstract
The charge density wave in the high-temperature superconductor YBa2Cu3O7−x (YBCO) has two different ordering tendencies differentiated by their c-axis correlations. These correspond to ferro- (F-CDW) and antiferro- (AF-CDW) couplings between CDWs in neighbouring CuO2 bilayers. This discovery has prompted several fundamental questions: how does superconductivity adjust to two competing orders and are either of these orders responsible for the electronic reconstruction? Here we use x-ray diffraction to study YBa2Cu3O6.67 as a function of magnetic field and temperature. We show that regions with F-CDW correlations suppress superconductivity more strongly than those with AF-CDW correlations. This implies that an inhomogeneous superconducting state exists, in which some regions show a fragile form of superconductivity. By comparison of F-CDW and AF-CDW correlation lengths, it is concluded that F-CDW ordering is sufficiently long-range to modify the electronic structure. Our study thus suggests that F-CDW correlations impact both the superconducting and normal state properties of YBCO.
Similar content being viewed by others
Introduction
Many theories and experiments point to an inhomogeneous nature of the superconductivity of high-temperature cuprate superconductors (HTC). Indeed, inhomogeneity in the sense of phase separation, intertwining of competing order parameters, stripes or pair-density-waves1,2,3,4,5 may be at the heart of cuprate superconductivity. Underdoped high-temperature superconductors exhibit competing tendencies towards charge-density-wave (CDW) and superconducting (SC) orders6,7,8,9,10. This means that the effects of quenched disorder may induce defects that disrupt or suppress the CDW order, and in turn lead to a resurgence of superconductivity. The application of a magnetic field provides another interesting control parameter because the introduction of vortices into the superconducting state tips the energy balance in favour of CDW order.
In this context, YBa2Cu3O6+x (YBCO) is an important model system where large onset temperatures of both superconductivity and CDW order are present8,9,11,12. As a function of both magnetic field and uniaxial stress, two intimately related CDW ordering tendencies have been realised13,14,15,16. These correspond to different ordering patterns along the c-axis13,14, as schematically illustrated in Fig. 1. Magnetic field (or uniaxial stress) induces an in-phase ferro-coupled CDW (F-CDW) along the c-axis14,15 on top of the original out-of-phase bi-axial antiferro-coupled CDW (AF-CDW) order. Quenched disorder, that in YBCO is naturally introduced through imperfections in the chain layer, may have important implications for the CDW order15,17,18. It is believed that disorder locally favours the AF-CDW order and that the evolution of the CDW correlations with magnetic field may be understood as a crossover transition in which the CDW coupling along the c-axis varies14,18. The introduction of normal vortices with field leads to a more ordered CDW. This suggests disorder associated with vortices is of secondary importance.
a Representative CDW stacking sequences along the c-axis at low and high magnetic field for a fixed point in the a-b plane. A and B represent the two possible phases in the CDW modulation in each bilayer14. b Extension to the b-c plane. Spatially separated regions where the F-CDW (ℓ = 1) and AF-CDW (ℓ = 1 ∕ 2) correlations in a are present. c, d Idealised diffraction intensity for the F-CDW and AF-CDW correlations at high field.
More recently, the case where competing tendencies towards CDW and SC order in the presence of topological defects due to quenched disorder and of vortices in the SC state has been considered17. It was demonstrated theoretically that a fragile SC state can occur at low temperatures and at high magnetic field. This state is based on regions in which the CDW order is weakened due to defects where locally superconducting halos can form. These regions can then couple to form a state with global SC phase coherence. Given the richness of theoretical possibilities when superconductivity competes with CDW order in the presence of vortices and weak disorder, it is of great interest to scrutinise the magnetic field-induced phase competition in YBCO.
Here we present a comprehensive study of the F-CDW and AF-CDW correlations in YBa2Cu3O6.67. By varying the magnetic field strength B – applied along the crystallographic c-axis – and temperature T, new insights into in-plane correlation lengths in the phase space (B, T) and competition of the correlations with superconductivity are obtained. F-CDW correlations are traced from low fields. From the diffraction intensities of F- and AF-CDW correlations versus temperature and magnetic field, the competition with superconductivity is studied. In particular, we extract the characteristic temperatures below which the respective CDW correlations are suppressed. They are significantly different at high magnetic fields, where the F-CDW correlation length becomes large, and the F-CDW correlations are suppressed at a significantly lower temperature scale. These different temperatures indicate where superconducting order becomes strong enough to influence the F-CDW and AF-CDW correlations. This observation suggests that superconductivity in the F-CDW regions is more strongly suppressed by the CDW. In this fashion, inhomogeneous superconductivity emerges due to competition with CDW correlations which are locally different. Consequently, regions with weaker or more fragile superconductivity are created. Another implication of our study is that in-plane correlation lengths can be deduced across the (B, T) phase space. We extract the in-plane correlation lengths at the field-dependent temperature scales set by the phase competition between superconductivity and the two CDW ordering tendencies. It is demonstrated that the F-CDW in-plane correlation length exceeds that of the AF-CDW order down to magnetic fields as small as 5 T. This result raises the question as to whether F-CDW correlations may be related to the Fermi surface reconstruction. Together with the recent finding that uniaxial strain induces F-CDW order16, our results suggest that F-CDW order in YBCO plays an important role for both the electronic structure and in rendering superconductivity into an inhomogeneous state that includes a more fragile flavour.
Results
Charge-density-waves in YBCO
Charge density wave (CDW) ordering has been found to be a universal property of hole doped cuprates6,7,8,9,19,20,21. YBCO has a bilayer structure with relatively close (d = 3.3 Å) pairs of CuO2 planes, separated by layers containing CuO chains running along the b-axis. In YBCO, the CDW results in weak scattering centred at QCDW(ℓ) = τ + qCDW + ℓc⋆, where τ is a reciprocal lattice point of the unmodulated structure and qCDW = (δa, 0, 0) or (0, δb, 0)8,9,11. For the YBa2Cu3O6.67 (Tc = 67 K) composition studied here, two charge-density-wave components have been identified13,14,15. The first, found in zero magnetic field, develops below TCDW = 140(10) K with δa = 0.305(2) and δb = 0.314(2). The phase of this CDW in neighbouring bilayers has an anti-parallel correlation with only weak correlation along the c-axis, manifested by a broad intensity peak with an ℓ ≈ 1 ∕ 2 modulation. The intensity of these rods of scattering is enhanced by the application of a magnetic field along the c-axis9,22,23. Figure 2a displays this rod in the (0, k, ℓ) scattering plane for B = 14 T. A second magnetic-field-induced CDW exists along the b-axis direction13,14,15. This CDW has the same phase across neighbouring bilayers (ferro-coupled) and hence is manifested by a ℓ = 1 reflection on top of the ℓ ~ 1 ∕ 2 scattering rod (see Fig. 1c, d and Fig. 2a, b). In this paper we use the term “CDW order” loosely to mean static CDW correlations with a finite correlation length which develop at low temperatures and high magnetic fields.
a, b X-ray diffraction intensity maps around (0, 4 − δb, ℓ) with B = 14 T and 16.5 T, respectively. c, d, and e are ℓ- and k-dependent intensity profiles along the vertical and horizontal white dashed lines in a, b. f Delineation of the scattering plane and probed reciprocal space (grey). The in-plane incommensurability δb is defined by the horizontal arrow. g Diffraction intensity map displayed in false colour, as a function of ℓ and magnetic field. Source data are provided as a Source Data file.
Modelling14 of the c-axis correlations in YBCO suggests that magnetic field changes the inter-bilayer (IB) coupling of the CDW. This is an indirect effect in that the field suppresses the superconductivity which changes the IB coupling. The CDW will also be pinned by defects in the chain layer. Thus, we can view the evolution of the CDW with field as a crossover transition. At low field B ≲ 14 T, the CDW is pinned and the system displays weakly correlated ℓ = 1 ∕ 2 order. At very high field (B ≫ 14 T), in-phase coupling of the bilayers is favoured, which is sufficient to overcome the pinning and forms a state with long range l = 1 order along the b-axis only.
Two-component analysis
The two CDW orderings can be analysed separately as a function of temperature and magnetic field. Here we employ three separately independent but consistent methodologies for data analysis. First, as shown in Figs. 2d, e and 3d, e, the in-plane correlations and diffraction intensities can be inferred from Gaussian fits versus Qb at ℓ = 1 ∕ 2 and ℓ = 1. The correlation lengths are defined as \({\xi }_{b}^{\ell =1/2,1}=1/{\sigma }_{b}^{\ell =1/2,1}\) where \({\sigma }_{b}^{\ell =1/2}\) and \({\sigma }_{b}^{\ell =1}\) are the respective standard deviations. Second, along the (0, 4 − δb, ℓ) direction, the out-of-plane correlation lengths \({\xi }_{c}^{\ell =1/2,1}\) can be deduced by fitting with a double Gaussian function (Fig. 3c). The centre in ℓ of the broad AF-CDW signal rises slightly with magnetic field14 (the values that were fitted and at high field were all ~0.6), however, for simplicity we label the correlation lengths as \({\xi }_{c}^{\ell =1/2,1}\). We find that the temperature dependence of the broad peak is insensitive to the exact ℓ value (see supplementary Fig. 1), and all h- and k-scans were carried out at ℓ = 1 ∕ 2. Third, as the ℓ ≈ 1 ∕ 2 CDW component has essentially no magnetic field dependence in the narrow range 14–16 T (Fig. 2c), it is possible – by subtraction of Gaussian fits of 14 T diffraction intensities – to isolate the field and temperature dependence of the ℓ = 1 diffraction intensity.
a, b Raw diffraction intensities along (0, k, 1 ∕ 2) and (0, k, 1) for temperatures and magnetic fields as indicated. c Raw diffraction intensity profile along the (0, 4 − δb, ℓ) direction, fitted with two Gaussian functions on a linear background. The dark grey line is their sum. d, e Magnetic-field-dependent evolution of the respective correlation lengths for ℓ = 1 ∕ 2 and 1. The dashed lines and shaded area are guides to the eye. Error bars are standard deviations determined by counting statistics and Gaussian fits, respectively. Source data are provided as a Source Data file.
Superconductivity and CDW phase competition
By subtracting the Gaussian fit of 14 T diffraction intensity from those measured at high fields (see Supplementary Figs. 2 and 3), the temperature dependence of the ℓ = 1 CDW component at high field is deduced. For both B = 14.5 and 15 T, the diffraction intensity is maximum at a finite temperature of T ≈ 25 K (see Fig. 4). It thus demonstrates that the CDW intensity at ℓ = 1 has a non-monotonic temperature dependence. This effect can also be analysed from in-plane intensity scans. In Fig. 5a–f, we plot the peak intensity I of Gaussian fits to k-scans through (0, 4 − δb, ℓ) as a function of temperature and at magnetic fields as indicated. At ℓ = 1 ∕ 2 (Fig. 5), the zero-field temperature dependence displays a maximum at \(T={T}_{\max }^{\ell =1/2}={T}_{{\rm{c}}}\) – consistent with previous reports8,9. As the magnetic field is increased, \({T}_{\max }^{\ell =1/2}(B)\) moves to lower temperatures. This effect is found along both a− and b−axis directions (Fig. 5a, b).
a–c ℓ-dependent diffraction intensity after subtracting the Gaussian fitted B = 14 T profiles at temperatures and magnetic fields as indicated. In this fashion, the relative field effects above 14 T are extracted. Raw 14 T ℓ-scans are shown in Supplementary Fig. 2. Data at different fields have been given an arbitrary vertical shift for the sake of presentation. Source data are provided as a Source Data file.
a–f Peak intensities versus temperature measured at ℓ = 1 ∕ 2 and 1 charge-density-wave modulation vectors versus temperature. a, b Intensity for (2 − δa, 0, 1 ∕ 2) and (0, 4 − δa, 1 ∕ 2) with magnetic fields as indicated. Data in a and zero-field data in b are re-plotted from ref. 9. Non-zero magnetic field data in b is ℓ ~ 0.6 amplitudes obtained from two-Gaussian fits of ℓ-scans. c Diffraction intensities at (0, 4 − δa, 1). d Same data as in c but zoomed on the low-field curves. For visibility, these curves have been given an arbitrary vertical shift – as indicated by the dashed base-lines. e, f Comparisons of ℓ = 1 and 1 ∕ 2 temperature dependencies and with that the associated maximum temperature scales \({T}_{\max }\) for the high- and low- magnetic field regimes respectively. Solid lines are fits to a skewed Gaussian (see text) on a linear background. Vertical bars are the \({T}_{\max }\) temperature scales defined by the maximum of the second derivative of the skewed Gaussian. Error bars are standard deviations due to counting statistics. Source data are provided as a Source Data file.
The behaviour of the ℓ = 1 intensity is very different. For low fields, B ≲ 10 T, the field induces a relatively weak response at ℓ = 1. It shows a maximum at similar temperature as the ℓ = 1 ∕ 2 component (see Fig. 5d). For 10 ≲ B ≲ 16 T, there is a rapid increase of I(ℓ = 1) at low temperature (Fig. 5c). The rapid temperature changes below 40 K associated with the crossover field may obscure any maximum due to competition with superconductivity. Then at higher fields B ≥ 14.5 T, a maximum again develops at \({T}_{\max }^{\ell =1}\) which is considerably lower than \({T}_{\max }^{\ell =1/2}\) for the same field. Direct comparisons of the high- and low-field diffraction intensities at ℓ = 1 ∕ 2 and ℓ = 1 versus temperature are shown in Fig. 5e, f. To extract \({T}_{\max }^{\ell =1/2}\) and \({T}_{\max }^{\ell =1}\) in a systematic fashion where the maximum is less well defined, the diffraction intensity versus temperature curves are fitted using a skewed Gaussian function \(\varGamma (T)=G(T)[1+{\mathrm{erf}}(\alpha T/\sqrt{2})]\) on a linear background (solid lines in Fig. 5). Here G(T) is a Gaussian function and \({\rm{erf}}(\alpha T)\) is the error function. The characteristic temperatures, shown in Fig. 6c versus magnetic field, are derived from the criterion \({T}_{\max }=\max (-{d}^{2}\varGamma /d{T}^{2})\) – this selects either a maximum or saturation point in the diffraction intensity versus temperature. We view \({T}_{\max }\) as a characteristic temperature scale for competition between CDW and superconducting order.
a Intensity map of (0, 4 – δb, 1) (logarithmic colour scale) as a function of temperature and magnetic field. The open circles indicate the temperature scale \({T}_{\max }^{\ell =1}\) below which suppression of the F-CDW order is observed (Figs. 4a–c). The curved dashed line is a guide to the eye. b Intensity map (linear colour scale) focusing on the high-magnetic-field region. The thin dashed lines indicate contour lines. c \({T}_{\max }\) scales inferred from the ℓ = 1 ∕ 2 (red) and ℓ = 1 (blue) charge-density-wave, respectively. Solid lines are guides to the eye. For the ℓ = 1 ∕ 2 \({T}_{\max }\) temperature scale, triangular and circular data points stem from data recorded along the a- (ref. 9) and b-axis (this work) respectively. As the \({T}_{\max }\) temperature scales indicate the onset of competition with superconductivity, these lines are labelled \({B}_{{\rm{c}}2}^{\ell =1}\) and \({B}_{c2}^{\ell =1/2}\) in the T → 0 limit. Error bars indicate the uncertainty related to the determination of the \({T}_{\max }\) temperature scales. d Correlation lengths \({\xi }_{b}^{\ell =1/2}\) and \({\xi }_{b}^{\ell =1}\) at the respective \({T}_{\max }\) temperature values. Source data are provided as a Source Data file.
Correlation lengths
Next, we turn to discuss the CDW correlation lengths. Fig. 3d, e show correlation lengths ξb and ξc for the lowest temperatures we measured: T = 22 and 3 K, respectively. The out-of-plane correlation length \({\xi }_{c}^{\ell =1/2} \sim 10\) Å is essentially magnetic-field independent while the in-plane correlation length \({\xi }_{b}^{\ell =1/2} \sim 60\) Å increases weakly. This is in strong contrast to the ℓ = 1 order that gains a dramatic out-of-plane “coherence” above 14 T (Fig. 3e). In fact, \({\xi }_{c}^{\ell =1}\) undergoes a fourfold increase from 14 to 16 T. Prior to this change, \({\xi }_{b}^{\ell =1}\) increases significantly for field strengths above 10 T (Fig. 3d). Thus at high-fields, the correlation length of the ℓ = 1 order exceeds that of the ℓ = 1 ∕ 2 CDW component in both in- and out-of-plane directions. It is also interesting to discuss the in-plane correlations at the \({T}_{\max }\) temperature scales. For ℓ = 1 ∕ 2, the correlation length \({\xi }_{b}^{\ell =1/2}\) at \({T}_{\max }\), is about 20% longer than that at T = 22 K in the studied field range (compare Figs. 3d and 6d). The ℓ = 1 in-plane correlations extend down to zero-magnetic-field (Fig. 6d). At \({T}_{\max },\,{\xi }_{b}^{\ell =1}\) exceeds \({\xi }_{b}^{\ell =1/2}\) for B > 5 T (Fig. 6d), while at 22 K this is only seen at 10 T (Fig. 3d). It is of great interest to note the interplay between the observed peak intensities and the correlation lengths. Our estimates of the k-space integrated intensity for the ℓ = 1 peak show approximately constant values versus field above 14 T (See Supplementary Fig. 4). This implies that the F-CDW develops longer range correlations with increasing field, but does not increase in amplitude or occupy a significantly larger fraction of the sample volume. The integrated intensity of AF-CDW order does not significantly increase over a range 0–16 T, which leads us to the similar conclusion for the development of AF-CDW.
Discussion
The zero-field maximum at Tc of the ℓ = 1 ∕ 2 CDW diffraction intensity has previously been the subject of different interpretations9,24,25,26. Our view is that the intensity maximum is the result of phase competition between ℓ = 1 ∕ 2 CDW order and superconductivity. As diffraction intensities are proportional to the square of the CDW order parameter23,27, temperature (or field) dependence reflects the evolution of the CDW order. Upon cooling below Tc, the superconducting order gradually grows, resulting in a partial suppression of the CDW. The balance of this competition can be tipped in favour of the CDW order by application of an external c-axis magnetic field that acts to suppress the superconducting order parameter. With increasing field-strength, the \({T}_{\max }\) temperature scale shifts to lower values (Fig. 6c) delineating a line in the magnetic field phase diagram that provides a firm lower bound for the border between superconductivity and normal state. An approximate extrapolation to the T → 0 limit yields an upper critical field consistent with that inferred from most bulk measurements28,29,30,31. It should be stressed that there is still an ongoing debate about the exact field scale of Hc232. A similar interpretation can be applied to the temperature dependence of the ℓ = 1 response. As shown in Fig. 6c, the ℓ = 1 maximum temperature scale \({T}_{\max }^{\ell =1}\) is different from \({T}_{\max }\) for ℓ = 1 ∕ 2. This implies the existence of two superconducting onset temperatures at a given applied field or for a fixed temperature two upper critical fields \({B}_{\max }^{\ell =1/2}\) and \({B}_{\max }^{\ell =1}\). The different competing temperature and magnetic field scales observed at ℓ = 1 ∕ 2 and ℓ = 1 lend support to our starting hypothesis, namely that there are two different CDW order constituents. In what follows, we discuss the relationship between the CDW and the Fermi surface reconstruction.
It is interesting to consider the implications of our correlation length measurements for the electronic band structure33,34. Although the unfolded structure has not been resolved, the Fermi surface reconstruction (FSR) revealed by quantum oscillation35,36,37 and transport38,39,40 experiments is undoubtedly linked to the appearance of CDW order. It has, however, been debated whether the reconstruction is triggered by the F-CDW or AF-CDW or whether the CDW order is triggered by FSR. The fact that ARPES does not see any FSR41 has been used as an argument for a field-induced electronic transformation, rather than the reconstruction being already present at zero-field. However, the sign change of the Hall effect persists down to as low as 5 T42, casting doubt on whether the high-field F-CDW order is responsible for the reconstruction. Recent detailed considerations of CDW correlations concluded that the ℓ = 1 CDW order has about the right in-plane correlation length to explain the quantum oscillation experiments whereas the AF-CDW correlations are too short33. We note that at B = 16 T the in-plane correlation lengths of the F-CDW and AF-CDW correlations are 300 and 100 Å, respectively. Another signature of FSR is the Hall coefficient which changes sign from positive to negative on cooling38,39,40,43. It has been demonstrated that for YBCO p ≈ 0.12, the characteristic sign-changing temperature TB = 66 K ≃ Tc is essentially independent of magnetic field38. The sign change in the Hall effect persists down to much lower fields, B = 5 T, than the quantum oscillations mentioned above. Thus it is interesting to consider whether the F-CDW correlations might be involved in the FSR at lower fields. We have previously demonstrated14 that for high magnetic fields, B = 16.5 T, \({\xi }_{b}^{\ell =1}\ > \ {\xi }_{b}^{\ell =1/2}\) for T ≈ 66 K ≈ TB. We can see from Fig. 6d that this condition persists at \({T}_{\max }\) as B → 0 and \({T}_{\max }\to {T}_{B}\). In fact, we find that \({\xi }_{b}^{\ell =1}\) is longer than \({\xi }_{b}^{\ell =1/2}\) down to 5 T where \({T}_{\max }\approx 55\) K. Thus, F-CDW correlations with \({\xi }_{b}^{\ell =1}\ \ge \ 100\) Å \(\approx {\xi }_{b}^{\ell =1/2}\) exist for temperatures T ≲ TB and to lowest measured magnetic fields. Considering the (B, T) phase space, both the F-CDW and AF-CDW correlations could be associated with the Fermi surface reconstruction inferred from Hall effect experiments. However, in both cases the correlation lengths are shorter than the expected threshold for quantum oscillations to be observed. Therefore, either the two-dimensionality of the CuO2 enables the electronic structure to fold also for intermediate correlation lengths or the Fermi surface reconstruction has an entirely different causality and charge order is consequence rather than a trigger.
A central part of this work is the experimental probing of ℓ = 1 in-plane correlations as a function of magnetic field and temperature. With this, the interaction between ℓ = 1 CDW order and superconductivity can be discussed. The ℓ = 1 order displays different competing field and temperature scales, compared to the extensively studied ℓ = 1 ∕ 2 order8,9,22,23,44,45. The \({T}_{\max }^{\ell =1}\) and \({T}_{\max }^{\ell =1/2}\) temperature scales indicate the point where superconductivity is strong enough to partially suppress the respective CDW orders. In presence of large enough magnetic fields B ≥ 14 T, we find that \({T}_{\max }^{\ell =1}\ <\ {T}_{\max }^{\ell =1/2}\) (Fig. 5e). Extrapolating \({T}_{\max }^{\ell =1}\ \to \ 0\), a field scale significantly smaller than Bc2, inferred from transport measurements28,29,40, is reached. In contrast, \({T}_{\max }^{\ell =1/2}\) may track Bc2 at high fields. This suggests that superconductivity occuring in regions with ℓ = 1 CDW order is weaker than in ℓ = 1 ∕ 2 regions. This is especially the case when the ℓ = 1 correlations develop coherence along the c-axis. The emergence of ferro-coupled long-range ℓ = 1 CDW order, therefore, induces a fragile state of superconductivity. In the presence of quenched disorder, such a fragile state of superconductivity has, in fact, been predicted theoretically17. Disorder may therefore play an important role for the organisation of the observed phase diagram17,18 and the interaction between CDW order and superconductivity. The exact role of different types of disorder (oxygen-chain vacancies, impurities and vortices) on the AF-CDW order and the fragile superconducting state should be clarified in future experiments. The regime with competition induced weakening of superconductivity has been the setting for predictions of pair-density-wave formations3,46,47,48,49,50. The existence of such a state should manifest itself through a QPDW = QCDW∕2 modulation. Although STM studies have reported evidence for PDW order in the halos of vortices51, there are no x-ray diffraction results supporting the existence of a PDW order in YBCO.
Finally, we comment on the S-shape of BDOS versus temperature reported by density-of-states sensitive probes (specific heat, thermal conductivity, and NMR)52. For B > BDOS, the density-of-states is magnetic-field independent and the S-shape, refers to an effect where BDOS, in addition to standard Bardeen-Cooper-Schrieffer (BCS) behaviour, displays an sudden increase in the T → 0 limit. Our results do not exclude the possibility that \({B}_{\max }^{\ell =1}\) would increase steeply in the T → 0 limit. In fact, it is possible that the ℓ = 1 ∕ 2 and ℓ = 1 superconducting flavours couple in the T → 0 limit to form a uniform condensate. If so, \({B}_{\max }^{\ell =1}\) and \({B}_{\max }^{\ell =1/2}\) should merge for T = 0. This would lead \({B}_{\max }^{\ell =1}\) to rise in the low-temperature limit and give it an S-shaped dependence reminiscent of that reported by thermodynamic and spin susceptibility probes of Bc252.
Methods
Experimental details
A high quality single crystal of Y Ba2Cu3O6.67 (Tc = 67 K) with ortho-VIII oxygen ordering was grown and detwinned as described in ref. 53. Hard x-ray (100 keV) diffraction experiments were carried out, with a horizontal 17 T magnet54, at PETRA III’s P07 triple-axis diffractometer at DESY (Hamburg, Germany). The YBCO sample was mounted with the magnetic field along the c-axis direction giving access to the Q = (h, 0, ℓ) scattering plane. The setup is identical to that described in ref. 14, with the exception of an improved sample cooling power allowing a base temperature of T ≈ 3 K to be reached. Scattering vectors are specified in orthorhombic notation with reciprocal lattice units (r.l.u) (2 π∕a, 2 π∕b, 2 π∕c) where a ≃ 3.82 Å, b ≃ 3.87 Å, and c ≃ 11.7 Å.
References
Zaanen, J. & Gunnarsson, O. Charged magnetic domain lines and the magnetism of high-T c oxides. Phys. Rev. B 40, 7391–7394 (1989).
Kivelson, S. A. et al. How to detect fluctuating stripes in the high-temperature superconductors. Rev. Mod. Phys. 75, 1201 (2003).
Berg, E. et al. Dynamical layer decoupling in a stripe-ordered high-T c superconductor. Phys. Rev. Lett. 99, 127003 (2007).
Machida, K. Magnetism in La2CuO4 based compounds. Physica C 158, 192 (1989).
Emery, V. J. & Kivelson, S. A. Frustrated electronic phase separation and high-temperature superconductors. Physica C 209, 597 (1993).
Tranquada, J. M., Sternlieb, B. J., Axe, J. D., Nakamura, Y. & Uchida, S. Evidence for stripe correlations of spins and holes in copper oxide superconductors. Nature 375, 561 (1995).
Wu, T. et al. Magnetic-field-induced charge-stripe order in the high temperature superconductor YBa2Cu3Oy. Nature 477, 191 (2011).
Ghiringhelli, G. et al. Long-range incommensurate charge fluctuations in (Y,Nd)Ba2Cu3O6+x. Science 337, 821 (2012).
Chang, J. et al. Direct observation of competition between superconductivity and charge density wave order in YBa2Cu3Oy. Nat. Phys. 8, 871 (2012).
LeBoeuf, D. et al. Thermodynamic phase diagram of static charge order in underdoped YBa2Cu3Oy. Nat. Phys. 9, 79 (2013).
Achkar, A. J. et al. Distinct charge orders in the planes and chains of ortho-III ordered YBa2Cu3O6+d identified by resonant elastic x-ray scattering. Phys. Rev. Lett. 109, 167001 (2012).
Wu, T. et al. Incipient charge order observed by nmr in the normal state of YBa2Cu3Oy. Nat. Commun. 6, 6438 (2015).
Gerber, S. et al. Three-dimensional charge density wave order in YBa2Cu3O6.67 at high magnetic fields. Science 350, 949 (2015).
Chang, J. et al. Magnetic field controlled charge density wave coupling in underdoped YBa2Cu3O6+x. Nat. Commun. 7, 11494 (2016).
Jang, H. et al. Ideal charge-density-wave order in the high-field state of superconducting YBCO. Proc. Natl Acad. Sci. USA 113, 14645 (2016).
Kim, H.-H. et al. Uniaxial pressure control of competing orders in a high-temperature superconductor. Science 362, 1040–1044 (2018).
Yu, Y. & Kivelson, S. A. Fragile superconductivity in the presence of weakly disordered charge density waves. Phys. Rev. B 99, 144513 (2019).
Caplan, Y. & Orgad, D. Dimensional crossover of charge-density wave correlations in the cuprates. Phys. Rev. Lett. 119, 107002 (2017).
Comin, R. et al. Charge order driven by Fermi-arc instability in Bi2Sr2−xLaxCuO6+δ. Science 343, 390 (2014).
daSilvaNeto, E. H. et al. Ubiquitous interplay between charge ordering and high-temperature superconductivity in cuprates. Science 343, 393 (2014).
Tabis, W. et al. Charge order and its connection with Fermi-liquid charge transport in a pristine high-Tc cuprate. Nat. Commun. 5, 5874 (2014).
Blackburn, E. et al. X-ray diffraction observations of a charge-density-wave order in superconducting ortho-II YBa2Cu2O6.54 single crystals in zero magnetic field. Phys. Rev. Lett. 110, 137004 (2013).
Hücker, M. et al. Competing charge, spin and superconducting orders in underdoped YBa2Cu3Oy. Phys. Rev. B 90, 054514 (2014).
Liu, Y.-H., Konik, R. M., Rice, T. M. & Zhang, F.-C. Giant phonon anomaly associated with superconducting fluctuations in the pseudogap phase of cuprates. Nat. Commun. 7, 10378 EP (2016).
Hayward, L. E., Hawthorn, D. G., Melko, R. G. & Sachdev, S. Angular fluctuations of a multicomponent order describe the pseudogap of YBa2Cu3O6+x. Science 343, 1336 (2014).
Nie, L., Sierens, L. E. H., Melko, R. G., Sachdev, S. & Kivelson, S. A. Fluctuating orders and quenched randomness in the cuprates. Phys. Rev. B 92, 174505 (2015).
Jaramillo, R. et al. Breakdown of the Bardeen-Cooper-Schrieffer ground state at a quantum phase transition. Nature 459, 405 (2009).
Ramshaw, B. J. et al. Vortex lattice melting and H c2 in underdoped YBa2Cu3Oy. Phys. Rev. B 86, 174501 (2012).
Grissonnanche, G. et al. Direct measurement of the upper critical field in a cuprate superconductor. Nat. Commun. 5, 3280 (2014).
Marcenat, C. et al. Calorimetric determination of the magnetic phase diagram of underdoped ortho ii in YBa2Cu3O6.54 single crystals. Nat. Commun. 6, 7927 (2015).
Zhou, R. et al. Spin susceptibility of charge-ordered YBa2Cu3Oy across the upper critical field. Proc. Natl Acad. Sci. USA 114, 13148 (2017).
Yu, F. et al. Magnetic phase diagram of underdoped YBa2Cu3Oy inferred from torque magnetization and thermal conductivity. Proc. Natl Acad. Sci. USA 113, 12667 (2016).
Gannot, Y., Ramshaw, B. J. & Kivelson, S. A. Fermi surface reconstruction by a charge-density-wave with finite correlation length. Phys. Rev. B 100, 045128 (2019).
Briffa, A. K. R. et al. Fermi surface reconstruction and quantum oscillations in underdoped YBa2Cu3O7−x modeled in a single bilayer with mirror symmetry broken by charge density waves. Phys. Rev. B 93, 094502 (2016).
Doiron-Leyraud, N. et al. Quantum oscillations and the fermi surface in an underdoped high-T c superconductor. Nature 447, 565 (2007).
Sebastian, S. E. & Proust, C. Quantum oscillations in hole-doped cuprates. Ann. Rev. Condens. Matter Phys. 6, 411 (2015).
Yao, H., Lee, D. H. & Kivelson, S. Fermi-surface reconstruction in a smectic phase of a high-temperature superconductor. Phys. Rev. B 84, 012507 (2011).
LeBoeuf, D. et al. Electron pockets in the Fermi surface of hole-doped high-Tc superconductors. Nature 450, 533 (2007).
LeBoeuf, D. et al. Lifshitz critical point in the cuprate superconductor YBa2Cu3Oy from high-field Hall effect measurements. Phys. Rev. B 83, 054506 (2011).
Chang, J. et al. Nernst and Seebeck coefficients of the cuprate superconductor YBa2Cu3O6.67 : a study of Fermi surface reconstruction. Phys. Rev. Lett. 104, 057005 (2010).
Hossain, M. A. et al. In situ doping control of the surface of high-temperature superconductors. Nat. Phys. 4, 527 (2008).
Cyr-Choinière, O. et al. Anisotropy of the Seebeck coefficient in the cuprate superconductor YBa2Cu3Oy : fermi-surface reconstruction by bidirectional charge order. Phys. Rev. X 7, 031042 (2017).
Laliberté, F. et al. Fermi-surface reconstruction by stripe order in cuprate superconductors. Nat. Commun. 2, 432 (2011).
Blanco-Canosa, S. et al. Momentum-dependent charge correlations in YBa2Cu3O6+δ. Phys. Rev. Lett. 110, 187001 (2013).
Blanco-Canosa, S. et al. Resonant x-ray scattering study of charge-density wave correlations in YBa2Cu3O6+x. Phys. Rev. B 90, 054513 (2014).
Agterberg, D. F. et al. The physics of pair density waves. Preprint at http://arxiv.org/abs/1904.09687 (2019).
Agterberg, D. F. & Tsunetsugu, H. Dislocations and vortices in pair-density-wave superconductors. Nat. Phys. 4, 639 (2008).
Pépin, C., de Carvalho, V. S., Kloss, T. & Montiel, X. Pseudogap, charge order, and pairing density wave at the hot spots in cuprate superconductors. Phys. Rev. B 90, 195207 (2014).
Dai, Z., Zhang, Y., Senthil, T. & Lee, P. A. Pair-density waves, charge-density waves, and vortices in high-T c cuprates. Phys. Rev. B 97, 174511 (2018).
Wang, Y. et al. Pair density waves in superconducting vortex halos. Phys. Rev. B 97, 174510 (2018).
Edkins, S. D. et al. Magnetic field-induced pair density wave state in the cuprate vortex halo. Science 764, 976 (2019).
Kačmarčík, J. et al. Unusual interplay between superconductivity and field-induced charge order in YBa2Cu3Oy. Phys. Rev. Lett. 121, 167002 (2018).
Liang, R., Bonn, D. A. & Hardy, W. N. Growth of YBCO single crystals by the self-flux technique. Philos. Mag. 92, 2563 (2012).
Holmes, A. T., Walsh, G. R., Blackburn, E., Forgan, E. M. & Savey-Bennett, M. A 17 T horizontal field cryomagnet with rapid sample change designed for beamline use. Rev. Sci. Instrum. 83, 023904 (2012).
Acknowledgements
We acknowledge informative discussions with Steve Kivelson. This work was supported by the EPSRC (grant numbers EP/R011141/ & EP/J016977/1), the U.S. Department of Energy (DOE), under Contract No. DE-AC02-98CH10886, the Wolfson Foundation, the Royal Society, the Office of Basic Energy Sciences, Division of Materials Science and Engineering, the Danish Agency for Science, Technology and Innovation under DANSCATT, the Leverhulme Trust, and the Swiss National Science Foundation. Parts of this research were carried out at beamline P07 at DESY, a member of Helmholtz Association HGF.
Author information
Authors and Affiliations
Contributions
R.L., D.A.B. and W.N.H. grew and characterised the YBa2Cu3O6.67 single crystal. E.B., O.I., A.T.H., N.B.C., M.H., E.M.F., O.G., U.R., M.v.Z., S.M.H. and J. Chang executed the x-ray diffraction experiments. J. Choi, O.I., S.G., S.M.H. and J. Chang analysed the data. All authors contributed to the manuscript. J. Choi and O.I. equally contributed to this work.
Corresponding authors
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Peer review information Nature Communications thanks Jun-Sik Lee and the other, anonymous, reviewer(s) for their contribution to the peer review of this work. Peer reviewer reports are available.
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary information
Source data
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Choi, J., Ivashko, O., Blackburn, E. et al. Spatially inhomogeneous competition between superconductivity and the charge density wave in YBa2Cu3O6.67. Nat Commun 11, 990 (2020). https://doi.org/10.1038/s41467-020-14536-1
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41467-020-14536-1
This article is cited by
-
Using strain to uncover the interplay between two- and three-dimensional charge density waves in high-temperature superconducting YBa2Cu3Oy
Nature Communications (2024)
-
Unified energy law for fluctuating density wave orders in cuprate pseudogap phase
Communications Physics (2022)
-
Locally commensurate charge-density wave with three-unit-cell periodicity in YBa2Cu3Oy
Nature Communications (2021)
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.