Broken translational and rotational symmetry via charge stripe order in underdoped YBa2Cu3O6+y

See allHide authors and affiliations

Science  20 Mar 2015:
Vol. 347, Issue 6228, pp. 1335-1339
DOI: 10.1126/science.1258399

Picking out the elusive stripes

Copper-oxide superconductors have periodic modulations of charge density. Typically, the modulation is not the same for the whole crystal, but breaks up into small nanosized domains. Bulk experiments show that the density is modulated along both axes in the copper-oxide plane, but it is not clear whether this is true only on the scale of the whole crystal or also locally, for each domain. Comin et al. analyzed the charge order in the compound YBa2Cu3O6+y, using resonant x-ray scattering, and found that it was consistent with a local unidirectional, so-called stripy, ordering.

Science, this issue p. 1335


After the discovery of stripelike order in lanthanum-based copper oxide superconductors, charge-ordering instabilities were observed in all cuprate families. However, it has proven difficult to distinguish between unidirectional (stripes) and bidirectional (checkerboard) charge order in yttrium- and bismuth-based materials. We used resonant x-ray scattering to measure the two-dimensional structure factor in the superconductor YBa2Cu3O6+y in reciprocal space. Our data reveal the presence of charge stripe order (i.e., locally unidirectional density waves), which may represent the true microscopic nature of charge modulation in cuprates. At the same time, we find that the well-established competition between charge order and superconductivity is stronger for charge correlations across the stripes than along them, which provides additional evidence for the intrinsic unidirectional nature of the charge order.

Recent studies of Y-based copper oxides have highlighted the importance of a charge-ordered electronic ground state, also termed a charge density wave (CDW), for the phenomenology of high-temperature superconductors (114). Experiments on the family of YBa2Cu3O6+y (YBCO) compounds have yielded a wealth of experimental results that have enabled advancements in our understanding of CDW instabilities and their interplay with superconductivity (911, 1522).

YBCO is a layered copper oxide–based material in which hole doping is controlled by the oxygen stoichiometry in the chain layer, characterized by uniaxial CuO chains running along the crystallographic b axis. In addition to ordering within the chain layer—attained via the periodic alternation of fully oxygenated and fully depleted CuO chains—recent experiments have extensively shown the presence of charge ordering in the CuO2 planes, with an incommensurate wave vector Q ≈ 0.31 reciprocal lattice units (23), corresponding to a period of approximately three unit cells in real space (911). Although the stripy nature of La-based cuprates is long established (13), the local symmetry of the CDW in YBCO has not yet been resolved. Both charge stripes (in the presence of 90° rotated domains) and a checkerboard pattern are consistent with the globally bidirectional structure of the CDW, which is characterized by wave vectors along both the a and b axes, at Qa ≈ (0.31, 0) and Qb ≈ (0, 0.31), respectively (10, 11, 2426). This leaves open the fundamental question of whether stripes are the underlying charge instability in the whole class of hole-doped cuprates.

We used resonant x-ray scattering (RXS) to study the local density correlations of the charge-ordered state and the interaction of this state with superconductivity (SC) in underdoped YBCO. The RXS technique, which is now at full maturity, represents a unique combination of diffraction (to probe reciprocal space) and resonant absorption (allowing element specificity and therefore site selectivity). RXS directly measures the structure factor S(Qx, Qy), where Qx and Qy represent the momenta along the reciprocal axes H and K, respectively. The structure factor is linked to the density-density correlation function and therefore to the CDW order parameter in momentum space (27). To reconstruct the two-dimensional (2D) structure factor with high resolution, we used a specifically devised RXS probing scheme whereby a charge-ordering peak is sliced along different directions in the (Qx, Qy) plane, parameterized by the azimuthal angle α (Fig. 1). The resulting 2D shape of the CDW peaks rules out checkerboard order and is instead consistent with a stripy nature of charge modulations in YBCO (28). We carried out RXS measurements in the vicinity of the CDW wave vectors Qa ≈ (0.31, 0) and Qb ≈ (0, 0.31) for three detwinned, oxygen-ordered YBCO samples: YBa2Cu3O6.51 (Y651, with hole doping p ≈ 0.10), YBa2Cu3O6.67 (Y667, p ≈ 0.12), and YBa2Cu3O6.75 (Y675, p ≈ 0.13).

Fig. 1 Charge order topology in momentum space.

(A) Schematic representation of the momentum structure of charge modulations in YBCO. Left inset: Selected momentum scans of the CDW peak along the b axis at Qb = (0, 0.31), for different azimuthal angles (α = 0°, 45°, 90°). Continuous lines represent Lorentzian fits; horizontal bars denote the linewidth ΔQ (HWHM). Right inset: Same as for the left inset, but for the CDW peak along the a axis at Qa ≈ (0.31, 0). (B) Color map of a series of Q-scans (normalized to the peak height) slicing the Qb peak between α = –90° and α = 90°; black bars represent the linewidth ΔQ, which is largest at α = 0°. (C to E) Polar plots of ΔQ as a function of the azimuthal angle α for Qa (red) and Qb (blue) in YBa2Cu3O6.51 (Y651), YBa2Cu3O6.67 (Y667), and YBa2Cu3O6.75 (Y675), respectively. Concentric gray circles are spaced by 0.01 Å–1; continuous lines are fits to an elliptic profile. Bottom right insets: CDW peaks represented as solid ellipses and compared with their rotated versions (hollow ellipses) for each doping.

In our experimental scheme, the CDW peaks are scanned in a radial geometry via control of the azimuthal angle α (29) (Fig. 1A). At the Cu-L3 edge, the measured signal is mainly sensitive to periodic variations in the Cu-2p → 3d transition energy (30, 31), which is a scalar quantity, even though the detailed contribution of a pure charge modulation versus ionic displacements to the RXS signal cannot be decoupled (32). In addition, the poor coherence of the CDW across the CuO2 planes (11, 24, 26) qualifies this electronic ordering as a 2D phenomenon, thus motivating our focus on the structure factor in the (Qx, Qy) plane. Representative scans of the CDW peak for different α values and at the superconducting critical temperature TTc are shown in the insets to Fig. 1A, for the Qb and Qa CDW peaks of a Y667 sample. A change in the peak half width at half maximum (HWHM) ΔQ between α = 0° and α = 90° is already apparent but is even better visualized in the color map of Fig. 1B, which shows the sequence of Q-scans versus azimuthal angle and the corresponding variation of ΔQ for Qb in the range α = –90° to 90°. This same procedure is applied to all three YBCO doping levels, for both the Qa and Qb CDW peaks; polar plots of ΔQ versus α are shown in Fig. 1, C to E, for Y651, Y667, and Y675, respectively. With the aid of the ellipse fits to the CDW profiles (continuous lines), four key aspects of these data stand out: (i) All peaks show a clear anisotropy between the two perpendicular directions α = 0° and α = 90°; (ii) for each doping, the Qa and Qb peaks have different shapes, and in the case of Y651 and Y667 this is more evident as the peaks are elongated along two different directions; (iii) the departure from an isotropic case, quantified by the elongation of the CDW ellipses, increases toward the underdoped regime and is opposite to the evolution of orthorhombicity, which is instead maximized at optimal doping [see (27) and fig. S4 for a more detailed discussion]; and (iv) the peak elongation at Qa and Qb, evolving from biaxial (Y651 and Y667) to uniaxial (Y675), is inconsistent with the doping independence of the uniaxial symmetry of the CuO chain layer, which rules out the possibility that the observed CDW peak structure is exclusively controlled by the crystal’s orthorhombic structure [however, the uniaxial anisotropy observed for Y675 might reflect a more pronounced interaction between the Cu-O planes and chains in this compound, possibly caused by the increase in orthorhombicity upon hole doping (27)].

The observed 2D peak shape indicates that the four-fold (C4) symmetry is broken at both the macro- and nanoscale, which is consistent with the emergence of a stripe-ordered state. In fact, under C4 symmetry the electronic density would be invariant under a 90° rotation in real space (xy, y → –x), which is equivalent to a 90° rotation in momentum space (QxQy, Qy → –Qx). Instead, the CDW structure factor S(Qx, Qy) is clearly not invariant under such operation, as shown in the insets of Fig. 1, C to E, which compare the original S(Qx, Qy) factors to their 90° rotated versions S(Qy, –Qx). This finding demonstrates an unambiguous breaking of global C4 symmetry in all investigated samples and might elucidate the origin of the anisotropy observed in the Nernst effect (20) and in optical birefringence measurements (33).

The real-space representation of charge order branches off into two possible scenarios: (i) a biaxial anisotropy, where both x- and y-elongated domains (34) are present (Fig. 2, A and C); (ii) a uniaxial anisotropy, where only y-elongated (or, equivalently, x-elongated) domains are found (Fig. 2, B and D). Note that these domains need not necessarily lie in the same layer, but they need to be present at the same time within the bulk of the material (e.g., they can be present in alternating layers while still leading to the same momentum structure). The momentum-space representation of the order parameter—and therefore of the electronic density fluctuations—is shown in the corresponding panels in Fig. 2, E to H, where Sa(Q) (red ellipses) and Sb(Q) (blue ellipses) represent the structure factor associated to the charge modulations along a and b, respectively. The profile of a single structure-factor peak is the result of two contributions: the underlying CDW symmetry as well as its 2D correlation length, which can also be anisotropic.

Fig. 2 Domain mesostructure in real and reciprocal space.

(A and B) Stripy domains along a (red stripes) and b (blue stripes) in the presence of biaxial (A) and uniaxial (B) anisotropic correlations. The domains are visualized in the same layer, although a situation in which they are present in alternating layers in a 90°-rotated arrangement is equally possible. (C and D) Checkerboard domains in the presence of biaxial (C) and uniaxial (D) anisotropic correlations. (E to H) Corresponding structure factors in reciprocal (Q) space. In case of the simultaneous presence of both CDW components (checkerboard), the imprinted correlations must be equal for the density wave along a and b, thus imposing an equivalent peak structure at Qa and Qb as seen in (G) and (H).

As a net result, the anisotropy of a single peak in Q-space cannot be used to discriminate between different CDW symmetries. Instead, the latter can be resolved by probing the 2D CDW structure factor—that is, by comparing the CDW peak shapes for Qa and Qb. Inspection of the diagrams in Fig. 2, E to H, reveals a common trait of checkerboard structures in momentum space, in that the following conditions (Fig. 2, G and H) must always hold by symmetry: Embedded Image = Embedded Image and Embedded Image = Embedded Image. That is, the peak broadening along a given direction must coincide for Δa and Δb (see bottom of Fig. 2, E to H, for case-specific conditions on the peak linewidths). Intuitively, this follows from the fact that for the checkerboard case, the charge modulations along a and b axes must be subject to the same correlation lengths within each domain—irrespective of its orientation—and must therefore lead to an equivalent peak broadening along the same direction in reciprocal space, in contrast to our findings for the CDW linewidths (27). From this symmetry analysis, we can conclude that for both uniaxial and biaxial anisotropy it is in principle possible (35) to discriminate between a pure checkerboard and a pure stripe charge order, even in the presence of a distribution of canted domains (see tables S1 to S3 for a complete classification). Ultimately, the inequivalence of the peak broadening ΔQ along different directions for all studied YBCO samples, combined with the macroscopic C4 symmetry breaking, provides clear evidence for the unidirectional (stripe) intrinsic nature of the charge order (28).

Having established the underlying stripelike character of charge modulations in the CuO2 planes, we turn to the temperature dependence of the longitudinal and transverse correlation lengths, respectively parallel and perpendicular to the specific ordering wave vector. These can be extracted by inverting the momentum HWHM ΔQ, as illustrated in Fig. 3A. Longitudinal correlations are then given by ξ|| = ΔQ||–1, where ΔQ|| represents the momentum linewidth in the direction parallel to the ordering wave vector; transverse correlations are given by ξ = ΔQ–1, where ΔQ represents the momentum linewidth in the direction perpendicular to the ordering wave vector.

We subsequently studied the temperature dependence of ξ|| and ξ for both the Qa and Qb ordering wave vectors (Fig. 3, B to F). We observed a rise of correlation lengths below the CDW onset near 150 K, followed by their suppression below the SC transition temperature Tc; this confirms the competition between these two orders, in agreement with recent energy-integrated as well as energy-resolved RXS studies (10, 11, 24, 25, 30). However, the drop in the correlation lengths below Tc (Δξ) was in all instances larger for the longitudinal correlations, or Δξ|| > Δξ. In particular, the discrepancy between Δξ|| and Δξ, although small for Y675, was quite substantial for the more underdoped Y667 and Y651. This anisotropy provides additional evidence for the unidirectional nature of the charge ordering and thus the breaking of C4 symmetry, because a bidirectional order would be expected to exhibit an isotropic drop in correlation length across Tc. The anisotropy has an opposite doping trend from the crystal orthorhombicity, whose associated anomalies across Tc increase with hole doping [as revealed, for example, by lattice expansivity measurements (27, 36)].

Fig. 3 Interplay between stripe order and the superconducting state.

(A) Schematics of the CDW structure factor in Q-space, illustrating how longitudinal (ξ||) and transverse (ξ) correlation lengths are extracted from the data by taking the inverse of the peak HWHM ΔQ. (B to F) Temperature dependence of ξ|| (green markers) and ξ (orange markers) for the various samples (where data are available) and CDW peaks investigated. Continuous lines are guides to the eye; double arrows highlight the drop in correlation lengths Δξ|| and Δξ below Tc. (G and H) Cartoons representing the evolution of the stripy nanodomains from above Tc (G) to below Tc (H), illustrating how longitudinal correlations undergo a larger suppression in the presence of superconducting order.

The inferred real-space representation of the evolution across Tc is schematically illustrated in Fig. 3, G and H, where nanodomains are used to pictorially represent a charge-ordered state with finite correlation lengths. We conclude that the largest change occurs along the direction perpendicular to the stripes. This reflects the tendency of the SC order to gain strength as temperature is lowered, primarily at the expense of longitudinal CDW correlations; the implication is that the mechanism responsible for the density fluctuations across the periodically modulated stripes might be the main one competing with the Cooper pairing process.

Our results may explain many common aspects between CDW physics in YBCO and the stripy cuprates from the La-based family, such as thermoelectric transport (37, 38), strength of the order parameter (39), out-of-equilibrium response (40, 41), and energy-dependent RXS response (30, 31). The nanoscopic nature of the stripe instability and the presence of both a- and b-oriented domains also clarify why this broken symmetry has been difficult to disentangle from a native bidirectional order (10, 11, 42), therefore requiring a tailored experimental scheme to resolve the 2D CDW structure factor S(Qx, Qy). The pronounced directionality in the competition between superconductivity and stripe order reveals the underlying interplay between particle-particle and particle-hole pairing in high-temperature superconductors and provides insights for an ultimate understanding of these materials.

Supplementary Materials

Materials and Methods

Supplementary Text

Figs. S1 to S9

Tables S1 to S3

References (4448)

References and Notes

  1. The small but finite doping dependence of the ordering wave vector is here neglected, because it is not the focus of our investigation. For more details, see (24, 25, 43).
  2. See supplementary materials on Science Online.
  3. Because of our inability to measure magnetic ordering, our concept of stripes only applies to the charge modulations, unlike the canonical definition used in La-based cuprates, which also implies concomitant magnetic order.
  4. However, x-ray diffraction experiments estimate the upper bound for ionic displacements to be around 7 × 10–4 a, or ~2.7 × 10–3 Å (11).
  5. We use the concept of domain to indicate a coherence subregion (i.e., a portion of space where density modulations are largely in phase); in this sense, the size of the domains is intended to reflect the magnitude of correlation lengths.
  6. Unless Embedded Image = Embedded Image, in which case the patterns shown in Fig. 2, F and H, become indistinguishable. This situation, however, is not encountered in our study.
  7. The data-fitting routines used for this analysis implicitly assume that the errors are normally distributed with zero mean and constant variance and that the fit function is a good description of the data. The coefficients and their sigma values (representing uncertainties) are estimates of what one would get if the same fit were performed an infinite number of times on the same underlying data (but with different noise each time) and then the mean and standard deviation were calculated for each coefficient.
  8. Acknowledgments: We thank P. Abbamonte, E. Blackburn, L. Braicovich, J. Geck, G. Ghiringhelli, B. Keimer, S. Kivelson, M. Le Tacon, and S. Sachdev for insightful discussions. Supported by the Max Planck–UBC Centre for Quantum Materials, fellowships from the Killam, Alfred P. Sloan, and Alexander von Humboldt Foundations, and a Natural Sciences and Engineering Research Council (NSERC) Steacie Memorial Fellowship (A.D.); the Canada Research Chairs Program (A.D., G.A.S.); and NSERC, CFI, and CIFAR Quantum Materials. All of the experiments were performed at beamline REIXS of the Canadian Light Source, which is funded by the CFI, NSERC, NRC, CIHR, the Government of Saskatchewan, WD Canada, and the University of Saskatchewan. R.C. acknowledges the receipt of support from the CLS Graduate Student Travel Support Program. E.H.d.S.N. acknowledges support from the CIFAR Global Academy.
View Abstract

Navigate This Article