Stabilization of three-dimensional charge order in YBa{}_{2}Cu{}_{3}O{}_{6+x} via epitaxial growth

Stabilization of three-dimensional charge order in YBaCuO
via epitaxial growth

M. Bluschke Max Planck Institute for Solid State Research, Heisenbergstr. 1, 70569 Stuttgart, Germany Helmholtz-Zentrum Berlin für Materialien und Energie, Wilhelm-Conrad-Röntgen-Campus BESSY II, Albert-Einstein-Str. 15, 12489 Berlin, Germany    A. Frano Department of Physics, University of California, Berkeley, California 94720, USA Advanced Light Source, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA Department of Physics, University of California San Diego, La Jolla, California 92093, USA    E. Schierle Helmholtz-Zentrum Berlin für Materialien und Energie, Wilhelm-Conrad-Röntgen-Campus BESSY II, Albert-Einstein-Str. 15, 12489 Berlin, Germany    D. Putzky Max Planck Institute for Solid State Research, Heisenbergstr. 1, 70569 Stuttgart, Germany    F. Ghorbani Max Planck Institute for Solid State Research, Heisenbergstr. 1, 70569 Stuttgart, Germany    R. Ortiz Max Planck Institute for Solid State Research, Heisenbergstr. 1, 70569 Stuttgart, Germany    H. Suzuki Max Planck Institute for Solid State Research, Heisenbergstr. 1, 70569 Stuttgart, Germany    G. Christiani Max Planck Institute for Solid State Research, Heisenbergstr. 1, 70569 Stuttgart, Germany    G. Logvenov Max Planck Institute for Solid State Research, Heisenbergstr. 1, 70569 Stuttgart, Germany    E. Weschke Helmholtz-Zentrum Berlin für Materialien und Energie, Wilhelm-Conrad-Röntgen-Campus BESSY II, Albert-Einstein-Str. 15, 12489 Berlin, Germany    R. J. Birgeneau Department of Physics, University of California, Berkeley, California 94720, USA    E. H. da Silva Neto Department of Physics, University of California, Davis, California 95616, USA    M. Minola Max Planck Institute for Solid State Research, Heisenbergstr. 1, 70569 Stuttgart, Germany    S. Blanco-Canosa CIC nanoGUNE, 20018 Donostia - San Sebastian, Basque Country, Spain IKERBASQUE, Basque Foundation for Science, 48011 Bilbao, Basque Country, Spain Donostia International Physcis Center, 20018 Donostia - San Sebastian, Basque Country, Spain    B. Keimer Max Planck Institute for Solid State Research, Heisenbergstr. 1, 70569 Stuttgart, Germany

Incommensurate charge order (CO) has been identified as the leading competitor of high-temperature superconductivity in all major families of layered copper oxides, but the perplexing variety of CO states in different cuprates has confounded investigations of its impact on the transport and thermodynamic properties. The three-dimensional (3D) CO observed in YBaCuO in high magnetic fields is of particular interest, because quantum transport measurements have revealed detailed information about the corresponding Fermi surface. Here we use resonant X-ray scattering to demonstrate 3D-CO in underdoped YBaCuO films grown epitaxially on SrTiO in the absence of magnetic fields. The resonance profiles indicate that Cu sites in the charge-reservoir layers participate in the CO state, and thus efficiently transmit CO correlations between adjacent CuO bilayer units. The results offer fresh perspectives for experiments elucidating the influence of 3D-CO on the electronic properties of cuprates without the need to apply high magnetic fields.

Recent nuclear magnetic resonance and X-ray scattering experiments have identified incommensurate charge order (CO) as a generic instability of copper oxide superconductors1, 2, 3, 4, 5. The competition of CO and high-temperature superconductivity was identified as a key factor shaping the phase diagram in the underdoped regime. These results have stimulated a profusion of research on the influence of collective charge fluctuations on the mechanism of superconductivity and on the anomalous transport and thermodynamic properties in the normal state of the underdoped cuprates6, 7, 8, 9. However, many aspects of charge ordering differ considerably between various families of cuprates. Specific CO states reported so far include two-dimensionally (2D) short-range ordered states in hole-doped Bi- and Hg-based materials 3, 4 and electron-doped NdCeCuO 5, 10, longer-range quasi-2D CO with different doping-dependent incommensurabilities in LaSrCuO  11, 12, 13 and YBaCuO 1, 2, 14 (YBCO), and three-dimensional (3D) CO in YBCO in high magnetic fields 15, 16, 17. These differences reflect the distinct lattice architectures and doping mechanisms of individual compound families that are difficult to capture in quantum many-body theories. The generic description of the interplay between CO and superconductivity thus requires a model system where CO correlations can develop with minimal disruption by doping-induced disorder.

Over the past decade, this role has been served by the YBaCuO system where resistivity, magnetization, and specific heat experiments in magnetic fields high enough to substantially degrade superconductivity have revealed pronounced quantum oscillations even in the underdoped regime 18, 19, 20, 21, 9. The oscillation period inferred from quantum transport data implies a small electron pocket, and thus indicates a reconstruction of the hole-like Fermi surface due to translational symmetry breaking. The size of the experimentally observed electron pocket is consistent with the one expected from CO-induced folding of the Fermi surface 22. The experimental discovery of 3D-CO with in-plane correlation lengths comparable to the cyclotron radius of the electron orbit in high magnetic fields has provided further support for models relating the fermiology of YBCO to CO. However, several important issues remain unresolved – including especially the lack of evidence for hole pockets (the counterpart of electron pockets in the reconstructed Fermi surface) in quantum oscillation data. At the same time, the need to impose high magnetic fields precludes experiments (such as angle-resolved photoemission) that could reveal the full Fermi surface geometry in the 3D-CO state.

Here we report a resonant X-ray scattering (RXS) study of moderately doped epitaxial thin films of YBCO, which has revealed Bragg reflections indicative of 3D-CO in zero magnetic field. The high-field 3D-CO state in bulk YBCO has thus far only been investigated by non-resonant hard X-ray diffraction, which probes the spatial distribution of the entire charge density. RXS with soft X-rays tuned to the Cu -absorption edge, on the other hand, is a spectroscopic measurement which selectively probes spatial variations of the electronic states near the Fermi level 23, 6. Based on spectroscopic information, RXS can also distinguish contributions of different crystallographic sites to charge-ordered superstructures 23. However, the technical challenges associated with producing high magnetic fields have thus far prohibited RXS studies of the 3D-CO state in bulk YBCO. Our thin-film data now demonstrate that Cu atoms in the charge-reservoir layers between the CuO sheets participate in the CO state. The RXS results thus yield fresh insight into the mechanisms stabilizing the unique 3D-CO state in YBCO and point out new opportunities to explore the influence of CO on the transport and thermodynamic properties of the cuprates in the absence of high magnetic fields.

We investigated underdoped YBCO films of different thicknesses that were grown epitaxially on SrTiO (STO) substrates using pulsed laser deposition (see Methods and Table 1). The in-plane lattice parameters of the nearly cubic perovskite structure of STO are  Å, larger than those in the orthorhombic structure of bulk YBCO ( Å  and  Å). STO thus applies tensile strain at least to the first layers deposited in the growth process. Epitaxial strain generally relaxes along the growth direction, and in our 50 nm thick films the volume-averaged in-plane lattice constants observed were  Å. A more profound difference between the crystal structure of our YBCO films and that of bulk YBCO is indicated by the absence of any of the conventional CuO chain ordering patterns24 associated with the charge-reservoir layer in our scattering experiments. The local orthorhombic distortion of the unit cell caused by the dopant oxygens in the charge-reservoir layer is therefore uncorrelated, and the crystal structure is essentially tetragonal to any non-local probe. This change in lattice symmetry is presumably a consequence of the epitaxial relationship with the STO substrate which enforces a tetragonal structure in the deposited layers; prior work has shown that such effects can persist over tens of nm in the growth direction 25, 26. We will discuss the influence of this distinct crystal symmetry on the CO pattern below.

Figure 1: Reciprocal space mapping. Three YBCO films grown on STO (001) are studied; (a) 50 nm, = 53 K; (b) 50 nm, = 55 K; (c) 20 nm, = 53 K. Scattering data were collected using polarized photons tuned to the Cu absorption edge. The upper panel shows the high and low temperature scattering intensity along the reciprocal space direction ( 0 1). The lower panel exhibits reciprocal space maps in the plane, these are produced by subtracting the scattered intensity detected at 300 K from the equivalent map measured at low temperature. The linear intensity scale is the same for all three plots in the upper panel, and as indicated by the color bars in the lower panel. The color maps were smoothed using a Gaussian filter.

.1 Results

2D and 3D charge-order correlations. Reciprocal space mapping of epitaxial YBCO films performed with resonant X-rays tuned to the Cu absorption edge reveals a rod of diffracted intensity in the plane peaked at incommensurate in-plane momentum transfer (Fig. 1), corresponding to the well studied 2D-CO in the CuO planes of bulk underdoped YBCO 2, 27. (The reciprocal-space coordinates , , are quoted in reciprocal lattice units, r.l.u., based on the tetragonal structure of the YBCO films.) The diffraction rod extending in the direction indicates 2D short-range order without correlation in the crystallographic -direction, qualitatively analogous to bulk YBCO in the absence of magnetic fields. The in-plane correlation lengths are comparable to those previously observed in bulk YBCO 27, indicating that the absence of oxygen order in the films does not have a major impact on the 2D-CO, in agreement with prior work on oxygen-disordered bulk YBCO 28.

Our main finding is the discovery of an additional temperature dependent diffraction peak centered at momentum transfer q = ( 0 1) in multiple YBCO films (Fig. 1 and Fig. 2c). The reflection intensity is sharply peaked in all three reciprocal-space directions, indicating 3D phase coherence, and is also present for positive in-plane momentum transfer q = (0.329 0 1). While it is clearly observable for incident photons with polarization (electric field ), the peak at q is completely suppressed with polarized X-rays (Fig. 3f). Since magnetic scattering is typically stronger for incident polarized X-rays, vanishing intensity of the 3D diffraction feature for incident polarization (/) leads us to conclude a non-magnetic origin and attribute the resonant diffraction peak at q to 3D-CO. Moreover the dramatic polarization dependence is not unexpected. In bulk YBaCuO at zero field, the ratio of scattering intensities / recorded at various values along the 2D-CO rod is observed to vanish as the momentum transfer along approaches   29. Considering all of these factors, we identify the 3D-CO reflection observed in thin films of YBCO as analogous to the behavior of bulk YBCO in high magnetic fields.

XAS and RXS energy profiles. RXS further provides spectroscopic information about the 3D-CO via analysis of the photon energy dependence. The spectral profile of the 3D-CO reflection near the Cu edge is plotted in Fig. 3a and compared with X-ray absorption spectroscopy (XAS) of the same film, measured with X-ray polarization parallel to the planar Cu-O bonds. Both the RXS and XAS processes involve the promotion of a core 2 electron into an unoccupied electronic state near the Fermi level. Since the XAS technique does not provide spatially resolved information, the measured lineshape represents the sum over all absorption processes associated with all ions in the material. In the YBCO structure there are two distinct Cu sites: Cu(I) located in the charge-reservoir layer containing the CuO chains, and Cu(II) in the CuO planes. As a function of oxygen doping the Cu(I) sites evolve from a ground-state configuration, associated with empty CuO chains, to a configuration associated with filled chains 30. Here indicates a hole on the ligand oxygens. Since the dopant oxygen ions are intercalated into the charge reservoir layer and not into the CuO planes, the Cu(II) sites receive a comparatively light hole doping and evolve between the and states as a function of increasing oxygen content. While the main maximum of the Cu XAS can be associated with transitions to the state of the Cu(II) sites, transitions to the and states of the Cu(I) sites occur for slightly higher photon energies, 932.7 and 934.3 eV respectively 30. The latter of these is clearly visible in the XAS profile, whereas the transition associated with the state is hardly recognizable. Inspection of the 3D-CO resonance profile reveals a main maximum coinciding with the XAS maximum, a strong shoulder around 933 eV, and a second significantly weaker shoulder at even higher energy, corresponding to the Cu(I) states. The enhancement of the resonant scattering intensity at energies associated with transitions to states of both the Cu(I) and Cu(II) ions indicates the participation of both Cu sites in the 3D-CO.

For comparison we have measured the resonance profiles of a number of other Bragg peaks associated with various ordering tendencies intrinsic to underdoped YBCO. Fig. 3b shows the resonance profile of the 2D-CO measured in one of our YBCO films. In agreement with previous studies 31, the 2D-CO reflection resonates solely at the energy corresponding to transitions into the states of the Cu(II) ions. In Fig. 3c we display the resonant enhancement of the (001) structural Bragg peak measured in the same film. The same XAS spectrum is plotted in Fig. 3a–c to facilitate comparison with the resonant scattering profiles.

Finally, Fig. 3d shows the resonant enhancement of an ortho-II oxygen superstructure peak measured in bulk single crystal YBCO, along with XAS measured in the same crystal with . Although the temperature dependence of the 3D-CO peak is inconsistent with that of a structural Bragg reflection (discussed below), one might be tempted to suggest that this peak originates from an ordered oxygen superstructure in the charge reservoir layer. In order to exclude this possibility we consider the polarization dependent resonant enhancement of the ortho-II oxygen superstructure peak displayed in Fig. 3d, which is consistent with previous resonant scattering studies of ortho-II, -III and -VIII ordered YBCO 30, 31. Since the dopant oxygens order into Cu-O chains in the charge reservoir layer, it is expected that this superstructure modulates the atomic scattering form factors of the Cu(I) ions. Accordingly the resonance profile of the ortho-II superstructure reflection is peaked at photon energies associated with the and states of the Cu(I) ions. The fact that the ortho-II peak does not resonate at the energy of the XAS maximum (associated with the Cu(II) sites), together with the polarization dependence, which actually demonstrates stronger scattering from the Cu-O chain with , leads us to conclude that the 3D-CO observed in our YBCO films does not correspond to a Cu-O chain ordering superstructure.

Figure 2: Temperature and doping dependence. (a,b) Temperature dependence of the scattering intensity along the ( 0 1) reciprocal space direction for 2 films with  = 53 K and 55 K, respectively. The background (dashed line) was achieved by slightly detuning the detector angle away from the Bragg condition, and then subtracted before integrating to produce the normalized data points in the upper panel (blue circles). For comparison the upper panel of (a) also shows the temperature dependent intensity of the 2D-CO rod collected (red triangles). The intensity scales of the two lower panels are identical and may be compared. (c) Scattering intensity at 12 K along ( 0 1) for 5 YBCO films of equal thickness but varying doping, as indicated by variation in the superconducting (see Methods). The commensurate momentum transfer is indicated by the dashed vertical line. All data were collected using polarized X-rays tuned to the Cu resonance.

In contrast to both the (001) Bragg peak and the ortho-II superstructure peak, neither the 2D-CO nor the 3D-CO were detected at the Cu resonance (951.4 eV). The 3D-CO peak intensity recorded at the Cu resonance is at least 8 times greater than the noise level, from which we derive a lower limit on the / branching ratio of at least 8. A similar lower limit on the / branching ratio can be extracted from RXS measurements of the 2D charge correlations in bulk YBCO27. This lower limit may help to distinguish the nature of the ordering type when modeling the resonant scattering profile across the Cu edges. Additionally this lower limit further discriminates the 2D- and 3D-CO reflections from the (001) Bragg peak and ortho-II superstructure peak, which both exhibit / branching ratios of less than 8.

Temperature dependence. The intensity of the 2D-CO rod evolves upon cooling as expected for bulk YBCO, onsetting below 120 K and achieving a maximum intensity at the superconducting (red triangles in the upper panel of Fig. 2a for the 50 nm film with  = 53 K). In stark contrast, the 3D-CO peak observed at q exhibits a monotonic intensity increase upon cooling, and no saturation (Fig. 2a upper panel, blue circles). Although a significant intensity remains at room temperature, the intensity variation between 12 K and 300 K is too strong to be accounted for by the Debye-Waller factor, and therefore makes a purely structural origin of this reflection unlikely. In fact, the greater thermal stability of the 3D-CO compared to the 2D-CO may be a simple consequence of the increased dimensionality of the CO correlations. The temperature dependence of the 3D-CO peak in a second film with 55 K, shown in Fig. 2b, is much more pronounced with almost no intensity remaining at room temperature. Additionally, its temperature dependence demonstrates an apparent maximum at , however this maximum likely arises from contamination of the -scan with intensity from the CO rod which in this film occurs at approximately the same in-plane momentum transfer as the 3D-CO and with comparable intensity. Unlike the integrated intensity, the correlation length is already saturated at 300 K and found to be between 8 and 10 nm, depending on the sample, showing no significant in-plane out-of-plane anisotropy. Among the series of samples studied, only the 20 nm film shows significant deviation in the correlation length of the 3D-CO (see Table 1).

Film thickness dependence. In order to test the extent to which CO correlations are affected by tensile strain imposed by the STO substrate, which expands the in-plane lattice parameters of YBCO, we have prepared a thinner 20 nm YBCO film with  K. Since epitaxial strain is expected to relax over the first 10-20 nm in the growth direction, a purely strain-induced mechanism for 3D phase coherence would be expected to exist only in the near-substrate region. Therefore reducing the film thickness from 50 nm to 20 nm (Fig. 1c), might be expected to leave the absolute intensity of the 3D-CO unchanged. In fact, however, both the intensity and correlation length of the 3D-CO reflection were found to decrease in the 20 nm thin film. Comparing the 20 and 50 nm films with  K, we find that the 3D-CO intensity is equal in both films after normalizing to the film thickness, but that the coherence length is significantly reduced in the 20 nm film (see Table 1). We conclude that the 3D-CO exists throughout the entire thickness of the 50 nm films, and that confining the 3D-CO by reducing the film thickness causes a decrease in its correlation length.

-axis correlation and - anisotropy. The CO rods observed in our YBCO films exhibit no enhancement of the scattering intensity towards (Fig. 1) indicating that the antiphase correlation between the CO in neighboring CuO planes, previously established for bulk underdoped YBCO 32, 16, is not present in our films. While one might be tempted to suggest that the CO rods observed in Fig. 1 are simply tails of the 3D-CO peak at , the in-plane separation of the rod and the 3D charge peak in Fig. 1a disfavors this interpretation. A more feasible scenario to explain the observed diffraction pattern involves a microscopic phase separation of regions hosting a 2D-CO (rod) and regions hosting a 3D-CO (peak at ). Alternatively, the 2D and 3D-CO may coexist spatially, but propagate along distinct in-plane crystallographic directions. In this scenario both reciprocal-space features appear along in Fig. 1 due to the presence of spatially separated domains in which the and directions are exchanged. Due to the effective tetragonality of the crystal structure, a scenario involving distinct charge ordering instabilities along distinct in-plane crystallographic directions would imply a nematic spontaneous symmetry breaking. Hard X-ray diffraction studies 16 of bulk single crystal YBCO in static magnetic fields reveal the abrupt emergence of a sharp reflection at q = (0 1) above 15 T, indicating a field-induced 3D phase coherence of the CO propagating along the direction, whereas along the CO remains two-dimensional. However, in ortho-ordered single crystal YBCO, the equivalence of the and directions is already broken by the CuO chain order in the charge-reservoir layers.

Figure 3: Photon energy and polarization dependence. (ad) Comparison of various energy dependent scattering intensities across the Cu edge. The data in (ac) were measured on a 50 nm YBCO thin film grown on SrTiO (no chain order), whereas the data in (d) was collected from a single crystal of ortho-II ordered YBCuO. (e) Scattering geometry used to study the 3D-CO reflection at q. (f) Polarization dependence of the 3D-CO reflection as measured in (e), and with photon energy tuned to the Cu edge.
(K)   thickness (nm) (r.l.u.) peak  area (arb. units) FWHM along   (r.l.u.) in-plane correlation length (nm)
  41   50   - - - -
  50   50   -0.329(1) 41(6) 0.016(2) 8(1)
  53   50   -0.329(1) 100(6) 0.012(1) 10(1)
  55   50   -0.329(1) 56(6) 0.015(2) 8(1)
  62   50   - - - -
  53   20   -0.324(1) 40(10) 0.040(6) 3(1)
Table 1: Fitting of the 3D-CO peaks. Fit parameters for Lorentzian fits to the 3D-CO peak measured at 12 K in a series of YBCO thin films grown on STO (100). The raw data for the 50 nm films is shown in Fig. 2c, and for the 20 nm film in the top panel of Fig. 1c. Since the 3D-CO is still observable in the 50 nm films at 300 K, fluorescence backgrounds of the type shown in Fig. 2a–b (dashed lines) were subtracted from the 12 K data before fitting. In contrast the 3D-CO reflection in the 20 nm film is broader, making necessary the subtraction of the 300 K data (top panel of Fig. 1c, red points) before fitting.

.2 Discussion

Our RXS study has demonstrated 3D-CO in YBCO thin films in zero magnetic field and uncovered spectroscopic evidence for the participation of Cu(I) sites in the charge reservoir layer in the CO state. The electronically active charge reservoir distinguishes YBCO from most other cuprates, where the dopant ions have closed electronic shells. Its prior manifestations include a uniaxial in-plane anisotropy of the transport properties due to the contributions of partially occupied electronic states in the CuO chains to the Fermi surface. The RXS energy profile of Fig. 3a now reveals a spatial modulation of the valence electron density in the charge-reservoir layer with the same incommensurate wave vector that characterizes the 2D-CO in the CuO layers. The charge-modulated CuO chain layer can thus effectively transmit CO correlation between adjacent CuO bilayers, bridging the large distance between CuO layers that characterizes most other cuprate families. The phenomenology of the 3D-CO state in the YBCO films shares key features with the one previously observed by non-resonant X-ray scattering in bulk YBCO in high magnetic fields 16, including particularly its coexistence with the 2D-CO state with the same incommensurability. It is therefore highly likely that the charge modulation of Cu chain states directly observed in our thin-film data also pertains to 3D-CO in bulk YBCO where RXS experiments in high magnetic fields have thus far not proven feasible. Our results therefore go a long way towards explaining why 3D-CO has so far only been observed in YBCO, and not in the other cuprate families.

Nonetheless, there are significant differences between the 3D-CO state in bulk YBCO and the one we have reported in YBCO films. Whereas the former state is only realized in high magnetic fields and temperatures below K, the latter state is stable in zero field and temperatures up to at least 300 K. Based on the observed dependence of CO correlations on film thickness, it is unlikely that tensile strain induced by the STO substrate is primarily responsible for the unexpected robustness of the 3D-CO state in the films. Likewise, pinning to lattice defects such as stacking faults and dislocations generated by steps on the substrate surface is unlikely to play a central role in generating the 3D-ordered state, because the corresponding Bragg reflections are consistently observed in multiple films, with in-plane correlation lengths of order 8 nm (Table 1). However, we cannot rule out that extended defects may promote phase separation into regions with and without 3D order.

The most salient difference between the lattice structures of bulk and thin-film YBCO is the absence of ortho ordering of the dopant oxygen ions in the films, such that their lattice symmetry is tetragonal. The difference in lattice symmetry, which is presumably enforced by epitaxy with the STO substrate, can have a substantial influence on the electronic structure. The CuO chains in bulk YBCO give rise to a separate Fermi surface sheet and cause a complex folding of the CuO plane-derived sheets, which degrade the nesting properties of the Fermi surface. In contrast, the Fermi surface of the effectively tetragonal structure of the films is expected to be more strongly nested, albeit broadened by disorder due to the randomly distributed oxygen dopants. The enhanced nesting, with nearly parallel Fermi surface segments in the Cu-O bond direction, is expected to favor a charge density wave instability. Detailed electronic-structure calculations are required to assess whether this effect is sufficient to explain the strongly enhanced 3D-CO instability we observed.

Alternatively, one might consider scenarios in which charge transfer 33, orbital hybridization 34, or phonon hybridization 35 at the YBCO-STO interface nucleate the 3D-CO instability which then propagates throughout the entire film. A related long-range electronic proximity effect has recently been reported for 2D-CO at (La,Ca)MnO-YBCO interfaces 36. We also note that incommensurate 3D-CO with comparable ordering temperatures is observed in insulating layered nickelates of composition SrNiO (where = rare earth) 37, where NiO octahedral tilts and rotations are believed to facilitate the coupling between electronic charge order and the lattice structure. However, analogous soft structural degrees of freedom are not apparent in the YBCO structure.

We end our discussion by pointing out parallels to other recent observations on cuprates. In particular, recent RXS experiments on overdoped bulk (Bi,Pb)SrCuO revealed 2D-CO diffraction features with large in-plane correlation lengths that persist up to at least room temperature 38. These findings were discussed in terms of enhanced nesting properties of the Fermi surface in proximity to the Lifshitz point, where a 2D van Hove singularity crosses the Fermi level. Together with our results, these observations illustrate the delicate stability conditions for CO in the cuprates, and offer new targets for theoretical research aimed at a realistic description of the phase behavior of these materials.

Independent of its microscopic origin, the 3D-CO state we have discovered in epitaxial YBCO films is much more accessible to experimental probes than the high-field state in bulk YBCO, and thus provides fresh perspectives for experimental exploration. To assess the influence of 3D-CO on the Fermi surface and macroscopic properties, it will be particularly interesting to perform angle-resolved photoemission and transport experiments on YBCO films and compare the results with those for bulk YBCO.

.3 Methods

Thin film growth and characterization. The YBaCuO thin films were grown by pulsed laser deposition on (100) oriented SrTiO (STO) substrates at the Max-Planck-Institute for Solid State Research in Stuttgart. A KrF excimer laser was used to ablate material from a commercially purchased stoichiometric YBaCuO target, with a 2 Hz pulse rate and an energy density of 1.6 J/cm. The STO substrate was kept under an oxygen partial pressure of 0.5 mbar and held at a temperature of 730 C during the deposition process. This procedure yields (001) oriented films of YBCO with in-plane twinning. Following the deposition process a series of 6 films were annealed for 15 to 60 minutes at 500 C in oxygen partial pressures ranging between 0.08 and 3 mbar, resulting in a series of superconducting critical temperatures (see Table 1). Whereas in bulk YBCO the hole doping level can be determined by measuring the out-of-plane lattice parameter 39, the tensile-strain-induced reduction of the lattice spacing precludes a quantitative evaluation of the doping using this method. Similarly, a universal relationship between hole doping and is difficult to establish for thin-film YBCO, as the critical temperature depends sensitively on the substrate, strain and growth conditions. Accordingly we have chosen to study a series of YBCO films all grown on STO substrates with constant growth parameters. Although we are unable to accurately determine the hole doping in the CuO planes, we attribute relative variations in among our systematically grown films to changes in the hole doping achieved during the annealing procedure. Furthermore we note that the lattice parameters of our films were all found to be within the range 11.69 - 11.73 Å, consistent with bulk underdoped YBCO.

The superconducting were determined using a mutual inductance device, and defined as the onset temperature of the Meissner effect. Electrical transport measurements were performed using a Quantum Design Physical Properties Measurement System (PPMS), indicating consistent with those determined from mutual inductance measurements, and revealing residual resistivities as low as 50 cm. Transmission electron microscopy of one of the same YBCO films studied here can be found in the supplementary information of Ref. 36.

Resonant soft X-ray absorption and scattering. The resonant X-ray scattering experiments were performed using linearly polarized photons produced by the elliptical undulator of the UE46-PGM1 beamline at the BESSY II synchrotron of the Helmholtz-Zentrum-Berlin. The sample environment consists of an ultra-high-vacuum two-circle diffractometer in which the sample is mounted directly to a liquid Helium flow cryostat. Careful adjustment of the sample position with respect to the X-ray beam was accomplished to a precision of 5 m enabling reproducible background intensities when cycling the temperature. This precision is necessary in order to accurately detect charge reflections with signals as weak as 1 of the background intensity. Scattered photons were detected using a standard photodiode, and the reciprocal space maps were produced by scanning the angle of the sample with respect to the incident beam for a series of detector angles. The photodiode does not discriminate photon energies, implying that the measured intensities represent an integration over all elastic and inelastic scattering processes.

The XAS spectra were measured at the same endstation, and collected in the total electron yield mode.


  • 1 Wu, T. et al. Magnetic-field-induced charge-stripe order in the high-temperature superconductor YBaCuO. Nature 477, 191–194 (2011).
  • 2 Ghiringhelli, G. et al. Long-range incommensurate charge fluctuations in (Y,Nd)BaCuO. Science 337, 821–825 (2012).
  • 3 da Silva Neto, E. H. et al. Ubiquitous interplay between charge ordering and high-temperature superconductivity in cuprates. Science 343, 393–396 (2014).
  • 4 Tabis, W. et al. Charge order and its connection with Fermi-liquid charge transport in a pristine high- cuprate. Nat. Commun. 5, 5875 (2014).
  • 5 da Silva Neto, E. H. et al. Charge ordering in the electron-doped superconductor NdCeCuO. Science 347, 282–285 (2015).
  • 6 Comin, R. & Damascelli, A. Resonant X-ray scattering studies of charge order in cuprates. Annu. Rev. Conden. Ma. P. 7, 369–405 (2016).
  • 7 Kloss, T., Montiel, X., de Carvalho, V. S., Freire, H. & Pépin, C. Charge orders, magnetism and pairings in the cuprate superconductors. Rep. Prog. Phys. 79, 084507 (2016).
  • 8 Sebastian, S. E. & Proust, C. Quantum oscillations in hole-doped cuprates. Annu. Rev. Conden. Ma. P. 6, 411–430 (2015).
  • 9 Yu, F. et al. Magnetic phase diagram of underdoped YBaCuO inferred from torque magnetization and thermal conductivity. Proc. Natl. Acad. Sci. USA 113, 12667–12672 (2016).
  • 10 da Silva Neto, E. H. et al. Doping-dependent charge order correlations in electron-doped cuprates. Sci. Adv. 2, e1600782 (2016).
  • 11 Cheong, S.-W. et al. Incommensurate magnetic fluctuations in . Phys. Rev. Lett. 67, 1791–1794 (1991).
  • 12 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–563 (1995).
  • 13 Yamada, K. et al. Doping dependence of the spatially modulated dynamical spin correlations and the superconducting-transition temperature in . Phys. Rev. B 57, 6165–6172 (1998).
  • 14 Blanco-Canosa, S. et al. Momentum-dependent charge correlations in YBaCuO superconductors probed by resonant X-ray scattering: evidence for three competing phases. Phys. Rev. Lett. 110, 187001 (2013).
  • 15 Gerber, S. et al. Three-dimensional charge density wave order in YBaCuO at high magnetic fields. Science 350, 949–952 (2015).
  • 16 Chang, J. et al. Magnetic field controlled charge density wave coupling in underdoped YBaCuO. Nat. Commun. 7, 11494 (2016).
  • 17 Jang, H. et al. Ideal charge-density-wave order in the high-field state of superconducting YBCO. Proc. Natl. Acad. Sci. USA 113, 14645–14650 (2016).
  • 18 Doiron-Leyraud, N. et al. Quantum oscillations and the Fermi surface in an underdoped high- superconductor. Nature 447, 565–568 (2007).
  • 19 Jaudet, C. et al. de Haas - van Alphen oscillations in the underdoped high-temperature superconductor . Phys. Rev. Lett. 100, 187005 (2008).
  • 20 Sebastian, S. E. et al. A multi-component Fermi surface in the vortex state of an underdoped high- superconductor. Nature 454, 200–203 (2008).
  • 21 Riggs, S. C. et al. Heat capacity through the magnetic-field-induced resistive transition in an underdoped high-temperature superconductor. Nat. Phys. 7, 332–335 (2011).
  • 22 Sebastian, S. E. et al. Normal-state nodal electronic structure in underdoped high- copper oxides. Nature 511, 61–64 (2014).
  • 23 Fink, J., Schierle, E., Weschke, E. & Geck, J. Resonant elastic soft X-ray scattering. Rep. Prog. Phys. 76, 056502 (2013).
  • 24 Zimmermann, M. v. et al. Oxygen-ordering superstructures in underdoped studied by hard X-ray diffraction. Phys. Rev. B 68, 104515 (2003).
  • 25 Butko, V. Y., Logvenov, G., Božović, N., Radović, Z. & Božović, I. Madelung strain in cuprate superconductors - a route to enhancement of the critical temperature. Adv. Mater. 21, 3644–3648 (2009).
  • 26 Frano, A. et al. Layer selective control of the lattice structure in oxide superlattices. Adv. Mater. 26, 258–262 (2014).
  • 27 Blanco-Canosa, S. et al. Resonant X-ray scattering study of charge-density wave correlations in YBaCuO. Phys. Rev. B 90, 054513 (2014).
  • 28 Achkar, A. J. et al. Impact of quenched oxygen disorder on charge density wave order in . Phys. Rev. Lett. 113, 107002 (2014).
  • 29 Achkar, A. J. et al. Orbital symmetry of charge-density-wave order in LaBaCuO and YBaCuO. Nat. Mater. 15, 616–620 (2016).
  • 30 Hawthorn, D. G. et al. Resonant elastic soft X-ray scattering in oxygen-ordered YBaCuO. Phys. Rev. B 84, 075125 (2011).
  • 31 Achkar, A. J. et al. Distinct charge orders in the planes and chains of ortho-III-ordered YBaCuO superconductors identified by resonant elastic X-ray scattering. Phys. Rev. Lett. 109, 167001 (2012).
  • 32 Forgan, E. M. et al. The microscopic structure of charge density waves in underdoped YBaCuO revealed by X-ray diffraction. Nat. Commun. 6, 10064 (2015).
  • 33 Chen, H. & Millis, A. Charge transfer driven emergent phenomena in oxide heterostructures. J. Phys.: Condens. Matter 29, 243001 (2017).
  • 34 Chakhalian, J. et al. Orbital reconstruction and covalent bonding at an oxide interface. Science 318, 1114–1117 (2007).
  • 35 Driza, N. et al. Long-range transfer of electron-phonon coupling in oxide superlattices. Nat. Mater. 11, 675–681 (2012).
  • 36 Frano, A. et al. Long-range charge-density-wave proximity effect at cuprate/manganate interfaces. Nat. Mater. 15, 831–834 (2016).
  • 37 Hücker, M. et al. Unidirectional diagonal order and three-dimensional stacking of charge stripes in orthorhombic and . Phys. Rev. B 74, 085112 (2006).
  • 38 Peng, Y. Y. et al. Re-entrant charge order in overdoped (Bi,Pb)SrCuO outside the pseudogap regime. arXiv:1705.06165v2 (2017).
  • 39 Liang, R., Bonn, D. A. & Hardy, W. N. Evaluation of CuO plane hole doping in YBaCuO single crystals. Phys. Rev. B 73, 180505 (2006).

.4 Acknowledgements

We thank the German Science Foundation (DFG) for financial support under grant No. SFB/TRR 80. S. B.-C. acknowledges IKERBASQUE and the Spanish Ministry of Economy, Industry and Competitiveness under the María de Maeztu Units of Excellence Programme - MDM-2016-0618. The work at the University of California, Berkeley and Lawrence Berkeley National Laboratory was supported by the Office of Science, Office of Basic Energy Sciences, Materials Sciences and Engineering Division, of the U.S. Department of Energy (DOE) under Contract No. DE-AC02-05-CH11231 within the Quantum Materials Program (KC2202).

.5 Author contributions

G.C. and G.L. grew the samples and characterized them with help from M.B., D.P. and R.O.; M.B., A.F., E.S., M.M., E.d.S.N. and S.B.-C. conceived the experiments and performed them with the help of E.W., F.G. and H.S.; M.B. performed the data analysis with the help of A.F. and E.d.S.N.; E.d.S.N., M.M., S.B.-C., R.B. and B.K. supervised the project. M.B. and B.K. wrote the manuscript together with contributions from all authors.

.6 Additional information

Competing financial interests: The authors declare no competing financial interests.

Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description