Hour-glass magnetic excitations induced by nanoscopic phase separation in cobalt oxides

Abstract

The magnetic excitations in the cuprate superconductors might be essential for an understanding of high-temperature superconductivity. In these cuprate superconductors the magnetic excitation spectrum resembles an hour-glass and certain resonant magnetic excitations within are believed to be connected to the pairing mechanism, which is corroborated by the observation of a universal linear scaling of superconducting gap and magnetic resonance energy. So far, charge stripes are widely believed to be involved in the physics of hour-glass spectra. Here we study an isostructural cobaltate that also exhibits an hour-glass magnetic spectrum. Instead of the expected charge stripe order we observe nano phase separation and unravel a microscopically split origin of hour-glass spectra on the nano scale pointing to a connection between the magnetic resonance peak and the spin gap originating in islands of the antiferromagnetic parent insulator. Our findings open new ways to theories of magnetic excitations and superconductivity in cuprate superconductors.

ARTICLE

Received 7 May 2014 | Accepted 2 Nov 2014 | Published 23 Dec 2014

Y. Drees1,*, Z.W. Li1,*, A. Ricci2,*, M. Rotter1, W. Schmidt3, D. Lamago4,5, O. Sobolev6,7, U. Rtt2,O. Gutowski2, M. Sprung2, A. Piovano8, J.P. Castellan4,5 & A.C. Komarek1

The magnetic excitations in the cuprate superconductors might be essential for an understanding of high-temperature superconductivity. In these cuprate superconductors the magnetic excitation spectrum resembles an hour-glass and certain resonant magnetic excitations within are believed to be connected to the pairing mechanism, which is corroborated by the observation of a universal linear scaling of superconducting gap and magnetic resonance energy. So far, charge stripes are widely believed to be involved in the physics of hour-glass spectra. Here we study an isostructural cobaltate that also exhibits an hour-glass magnetic spectrum. Instead of the expected charge stripe order we observe nano phase separation and unravel a microscopically split origin of hour-glass spectra on the nano scale pointing to a connection between the magnetic resonance peak and the spin gap originating in islands of the antiferromagnetic parent insulator. Our ndings open new ways to theories of magnetic excitations and superconductivity in cuprate superconductors.

DOI: 10.1038/ncomms6731 OPEN

Hour-glass magnetic excitations induced by nanoscopic phase separation in cobalt oxides

1 Max-Planck-Institute for Chemical Physics of Solids, Nthnitzer Str. 40, D-01187 Dresden, Germany. 2 Deutsches Elektronen-Synchrotron DESY, Notkestr. 85, 22603 Hamburg, Germany. 3 Jlich Centre for Neutron Science JCNS, Forschungszentrum Jlich GmbH, Outstation at ILL, BP 156, 6 Rue Jules Horowitz, 38042 Grenoble, France. 4 Laboratoire Lon Brillouin, CEA/CNRS,F-91191 Gif-sur Yvette Cedex, UMR12 CEA-CNRS, Bt 563 CEA Saclay, France. 5 Institute of Solid State Physics, Karlsruhe Institute of Technology, D-76021 Karlsruhe, Germany. 6 Forschungsneutronenquelle Heinz Maier-Leibnitz (FRM-II), TU Mnchen, Lichtenbergstr. 1, D-85747 Garching, Germany. 7 Georg-August-Universitat Gttingen, Institut fr Physikalische Chemie, Tammannstrasse 6, D-37077 Gttingen, Germany. 8 Institut Laue-Langevin (ILL), 6 Rue Jules Horowitz, F-38043 Grenoble, France. * These authors contributed equally to this work. Correspondence and requests for materials should be addressed to A.C.K. (email: mailto:[email protected]

Web End [email protected] ).

NATURE COMMUNICATIONS | 5:5731 | DOI: 10.1038/ncomms6731 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 1

& 2014 Macmillan Publishers Limited. All rights reserved.

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms6731

The spin excitations in high-temperature superconducting (HTSC) cuprates with a quasi two-dimensional layered crystal structure resemble an hour-glass if plotted as a

function of energy and momentum transfer113. Understanding the mechanism of these magnetic excitations and their relation to high-temperature superconductivity remains one of the great challenges in condensed matter physics. Therefore, the recent discovery of such excitation spectra in an isostructural but insulating cobaltate reference system La5/3Sr1/3CoO4 has attracted

enormous attention12. The assumed presence of charge stripes in this cobalt oxide pointed to the presence of (uctuating) stripes in the isostructural cuprates12 and, hence, suggested the importance of charge stripes for the physics in the high-temperature superconductors. In charge stripe phases the doped charges (holes) segregate into stripes, whereas antiferromagnetism is able to recover in the charge-depleted regions between these stripes. However, the very recent nding of an hour-glass dispersion in a checkerboard charge ordered cobaltate La1.6Sr0.4CoO4 has shown

that hour-glass spectra are not necessarily connected to charge stripe phases13. La5/3Sr1/3CoO4 is a compound with a hole-

concentration even further from half-doping that exhibits diagonal magnetic satellite reections around third-integer positions in reciprocal space. For this hole-doping the isostructural nickelates that also exhibit diagonal magnetic satellites at the same positions in reciprocal space exhibit the most stable diagonal charge stripe order. Hence, also for La5/3Sr1/3CoO4 the incommensurate magnetic peaks have been interpreted within a charge stripe scenario12,14.

Here we demonstrate that there are no signatures of charge stripes in La5/3Sr1/3CoO4. Instead, we observe a fractal microstructure with a nanoscopic distribution of undoped and hole-doped regions indicating that hour-glass spectra consist of excitations with distinct origin on the nanometer scale. Whereas the low-energy regime is still dominated by excitations within the nano-sized hole-doped regions, at highest energies only the nanoscopic undoped islands can be excited. In our nanoscale phase separation scenario the energy of the neck of the hour-glass could be interpreted as the spin-gap within undoped islands of type of the antiferromagnetic parent insulator La2CoO4.

ResultsHard X-ray diffraction measurements. Using 100 keV hard X-rays at the synchrotron we rst studied the charge correlations in La5/3Sr1/3CoO4. In a disordered charge stripe scenario the corresponding charge stripe ordering superstructure reections should exhibit similar peak widths as the magnetic reections15. However, in our highly precise synchrotron measurements we were not able to detect any indications for charge stripe ordering reections (for L 0 and L 3440) in La5/3Sr1/3CoO4 with an unprecedented

accuracy of about 10 9 of a strong Bragg reection, see Fig. 1a,b. A charge stripe-ordered La5/3Sr1/3NiO4 reference sample with same hole-concentration exhibits strong charge stripe ordering peaks in both equivalent scans. Note that the huge Co2 ions and the small

Co3 ions exhibit distinctly bigger differences in their ionic radii than the corresponding Ni ions16 such that any occurring stripe phases should be much easier to detect in cobaltates than in nickelates. Furthermore, magnetic stripes and charge stripes can not be decoupled due to the nonmagnetic properties of the Co3 ions17,18 and due to the highly insulating localized nature18 of this cobaltate system. Instead of charge stripe ordering reections we, surprisingly, observe weak checkerboard charge ordering intensities at half-integer H and K positions in reciprocal space that survive up to room temperature, see Fig. 1c,d. This important nding indicates that the robust checkerboard charge ordering correlations still persist also far away from half-doping down to

one-third doping in the La2 xSrxCoO4 system and that there are

no disordered charge stripe phases in this entire incommensurate magnetic regime of the La2 xSrxCoO4 system (1/3rxo1/2).

Instead, checkerboard charge ordered regions with roughly 0.7(2) nm correlation length appear. The remaining Co2 ions have to be distributed somewhere into the cobalt oxygen planes as the hole-doping in La5/3Sr1/3CoO4 is far below half-doping. Hence,

checkerboard charge ordered regions have to be interspersed with La2CoO4-like undoped islands within a nano-phase separation scenario.

Resonant micro X-ray diffraction measurements. We analysed the real-space distribution of the charge correlations in La5/3Sr1/3CoO4 by means of resonant scanning micro X-ray

diffraction (mXRD) measurements. The mXRD technique has been already successfully applied for the study of phase separation in the isostructural cuprates19,20. Here, we additionally applied an incident X-ray energy of 7,719 eV to measure any Co2/Co3 charge ordering correlations in resonance for a high contrast of charge ordering reections. Moving the sample under a 2 mm focussed beam with an x-y translator, we scanned a sample area of 202 210 mm2 collecting 10,605 different diffraction

patterns at B300 K. For each scanned point of the sample the(0.5 0.5 2) peak prole has been extracted and the full-width half-maximum (FWHM) along ab-direction has been evaluated to obtain the domain size of the checkerboard charge ordered regions. The spatial map of this domain size along ab-direction is shown in Fig. 1e. To quantify21 the distribution of undoped islands and checkerboard charge ordered regions, we calculated the probability density function (PDF) of the measured domain sizes, see Fig. 1f. Similar to that observed in the isostructural cuprates19,20 for oxygen ordering superstructure reections, the PDF tail presents a perfect power-law behaviour and can be tted by the function PDF(x) x a, where x is the ab-domain size of

checkerboard charge ordering and a 1.99(8) is the critical

exponent of the power-law. This particular behaviour highlights a scale-free distribution of the checkerboard charge ordered domain size in La5/3Sr1/3CoO4. From the non power-law region

of the PDF we can estimate the expected values of the checkerboard charge ordering ab-domain size that amounts to B0.8(2) nm. Thus, our mXRD measurements indicate a fractal microstructure up to a maximum observed domain size of

B39 nm and corroborate the presence of a disordered checkerboard charge ordered phase interspersed with undoped islands and, hence, our nano phase separation scenario.

Neutron scattering measurements. Studying the magnetic excitations in La5/3Sr1/3CoO4 in detail, we surprisingly observe a hitherto unknown high-energy magnon band above the entire classical hour-glass dispersion, see Fig. 2a. The polarization analysis conrms that this additional high-energy mode is magnetic in origin. Note that this is the rst observation of an additional separated high-energy magnon band appearing completely above an entire hour-glass dispersion. Usually, one expects optical modes in a magnon dispersion without hour-glass shape or the optical magnon mode merges with the low-energy dispersion such that both together are forming an hour-glass-shaped dispersion2224. Also a 40% Sr-doped cobaltate within the checkerboard charge ordered regime exhibits this same additional high-energy mode, see Fig. 2c. Our polarized neutron measurements uncover the interesting nature of these magnetic excitations, whereas all low-energy spin excitations within the entire classical hour-glass dispersion are mainly in-plane excitations (see Fig. 2b), the additional yet unknown high-energy mode appears to be an out-of-plane polarized

2 NATURE COMMUNICATIONS | 5:5731 | DOI: 10.1038/ncomms6731 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

& 2014 Macmillan Publishers Limited. All rights reserved.

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms6731 ARTICLE

La5/3Sr1/3CoO4 at 12 K

106

La5/3Sr1/3CoO4 at 12 K

Intensity (a.u.)Intensity (photons s1 )

Intensity (photons s1 )

106

105

104

103

102

LSNO at 24K

(H 3H 0) scan

LSNO at 20K

(H H 3) scan

105

104

103

102

Hard X-ray

3.0 2.0

Hard X-ray

2.0 1.0

2.5 1.5

H (r.l.u.)

H (r.l.u.)

4

3

2

1

0

La Sr CoO

Intensity (a.u.)

La5/3Sr1/3CoO4

Counts (a.u.)

Hard X-ray

LSNO

400

300

200

100

0

(H H 2) scan 285 K Fit12 K

Total fit LTO CBCO

0.0 0.5 1.0

H 1H 0 (r.l.u.)

LTO

+

CBCO CBCO

TStruct

1.0 1.2 1.4 1.6 1.8 2.0 0 100 200 300

Temperature (K)

Domain size x (nm)

H (r.l.u.)

100.0 nm

La5/3Sr1/3CoO4+ at 300 K

1

10.00 nm

0.1

0.01

PDF

1.000 nm

1E3

Resonant XRD

0.1000 nm

50 m 0.01

0.1 1 10 100

1E4

Figure 1 | Charge correlations Synchrotron radiation single-crystal X-ray diffraction measurements using 100 keV hard X-rays at beamline P07 at PETRA III. (a) Measured intensities of scans along (H 3 H 0) and b scans along (H H 3) for La5/3Sr1/3CoO4 and a roughly one-third hole-doped LSNO reference sample. Strong fundamental Bragg reections are of the order of 1010 photons/s. Intensities in a,b have been normalized to the monitor (incident beam) for a comparison of cobaltate and nickelate reference sample. The dashed lines are Lorentzian peaks with same integrated intensity as for the charge stripe ordering peaks of the LSNO reference sample but with peak widths that have been multiplied by the width ratio of the incommensurate magnetic peaks of both compounds measured with neutrons (inset of d); dotted lines are additionally normalized to the (200) peak intensities. No diagonal charge stripe correlations like in the nickelates could be observed for La5/3Sr1/3CoO4. (c) Weak checkerboard charge ordering reections appear at half-integer positions for even values of L already at room-temperature. The correlation length that could be obtained from Lorentzian ts indicates an average nanometer size of checkerboard charge ordered regions (B0.7(2) nm). Below the structural transition tetragonal-orthorhombic additional sharp octahedral tilting reections appear on top of the broad checkerboard charge ordering reections. (d) The temperature dependence of the total integrated intensities at half-integer peak positions. The inset shows neutron scans across the magnetic satellites for La5/3Sr1/3CoO4 and B1/3 hole-doped LSNO.

(e) Spatial map of the checkerboard charge ordered domain size distribution on the crystal surface along the ab-direction obtained from mXRD measurements at beamline P10 at PETRA-III. The distribution appears to be inhomogeneous and checkerboard charge-ordered regions with varying domain size can be observed. (f) The PDF36 of the checkerboard charge ordered domain size along the ab-direction. The behaviour of the PDF tail obeys a perfect power-law (PDF(x) x a) indicating a scale-free distribution in the checkerboard charge ordered domain sizes, thus corroborating our nano-phase-

separation scenario. Only B1% (B11%) of the measured sample exhibits domain sizes bigger than 50 (10) unit cells. Error bars are s.d.

magnon band, see Fig. 2a. As the low-energy excitations in La5/3Sr1/3CoO4 basically resemble on the excitations in the half-doped cobaltate La1.5Sr0.5CoO4 and as the high-energy excitations

resemble on the excitations of the undoped parent compound La2CoO4 including the out-of-plane polarized additional high-energy mode, our observations suggest that the hour-glass dispersion is not one single dispersion, but, that it consists of two dispersions with distinct origin, see Fig. 3. The rst dispersion arises from the low-energy magnetic excitations of the hole-doped regions that are governed by frustration and that

exhibit chiral or non-collinear magnetic phases and the second one arises from the high-energy magnetic excitations of the undoped islands strongly resembling on the dispersion of the undoped parent compound (La2CoO4) that exhibits a spin gap and an additional high-energy magnon mode. This model is shown schematically in Fig. 3. The intensity enhancement at the magnon merging point or neck of the hour-glass (we would like to call this by denition also simply resonance peak for the non-superconducting cobaltates throughout this article) that can be traced with energy scans at the planar antiferromagnetic

NATURE COMMUNICATIONS | 5:5731 | DOI: 10.1038/ncomms6731 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 3

& 2014 Macmillan Publishers Limited. All rights reserved.

ARTICLE

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms6731

5

300

250

200

200

x=1/3

NSFxx

SFzz ~2K

x=1/3

SFzz

SFyy

Unpolarized

SFyy ~2K

SFxx ~2K

4

Counts (a.u.)

Counts (a.u.)

E(meV)

Counts (a.u.)

Counts (a.u.)

E(meV)

3

150

100

50

0

0

2

1

x=0 reference material

40

Energy (meV)

10 20 30

0 1 2 3 4 5

0

Energy (meV) Energy (meV)

300

100

0

50

40

30

20

10

0 0.0 0.5

x=0.4 :

SFxx

NSFxx

x=1/3 :

1 x=1/3 SFzz

SFzz

SFyy

1

Scan

29 meV

1 20

H (r.l.u.)

K(r.l.u.)

10 20 30 0.0 0.5 0

50

40

30

20

10

0

0

K (r.l.u.)

x=1/3

(1.5 K 0) scans

x=0.4 (1.5 K 0) scans

Phonon

Phonon

(1K K 0 ) scan

(1K K 0 ) scan

0.0 0.5 1.0

K (r.l.u.)

1.0

MinNormalizec counts (a.u.)Max

K (r.l.u.)

Figure 2 | Spin excitations (a) Energy scans at the planar AFM wavevector studied by polarized neutrons. In-plane polarized excitations that can be detected in the Spin-Flip (SF) channel ||z (z> scattering plane) emerge at the neck of the hour-glass (red data points). Surprisingly, an additional high-energy mode appears at elevated energies. The polarization analysis indicates the out-of-plane character of these high-energy excitations as it can be detected in the SF channel ||y (blue data points; y within scattering plane, y>Q). In addition, also the results of unpolarized neutron scattering experiments of La5/3Sr1/3CoO4 measured at the IN8 spectrometer (dark cyan data points) and La2CoO4 taken from ref. 25 (dark yellow data points) are shown. As can be seen, the small additional 20 meV feature belongs to a phonon, whereas the resonance-like increase of intensity is located B15 meV that is close to the in-plane anisotropy gap of La2CoO4. (b) Study of the nature of the low-energy excitations by polarized neutrons measured by constant-Q scans at (1/2 e

1/2 e 0) at B2 K at the IN12 spectrometer. The in-plane character of these low-energy magnetic excitations is revealed. (c) Also La1.6Sr0.4CoO4 exhibits the same energy scan at the planar AFM wavevector measured with polarized neutrons in the channel ||Q that detects in-plane and out-of-plane excitations simultaneously. (d) For the outwards dispersing branches above the neck of the hour-glass also in-plane excitations can be detected in the SF channel ||z. The inset shows the corresponding scan in reciprocal space within a 29 meV constant-E map. (e,f) Comparison of strength of the high-energy magnetic intensities for different Sr-doping measured at the IN8 spectrometer. Intensities have been normalized to the Phonon peak intensity (black arrrow). Intensity error bars are statistical error bars calculated by the square root of intensity.

(AFM) wavevector exhibits a surprisingly similar energy dependence like in equivalent scans for La2CoO4, see Fig. 2a (dark yellow dots). This clearly points to a related origin of this resonance peak energy in La5/3Sr1/3CoO4 and the spin gap in

La2CoO4. This is also consistent with the temperature dependence of this energy scan and a beginning of a closing of this spin gap with increasing temperature, see Fig. 4a and also compare ref. 25 for the temperature dependence of the long-range ordered undoped compound La2CoO4 (with higher magnetic ordering temperature). In contrast the high-energy part of the hour-glass-shaped dispersion in La5/3Sr1/3CoO4 does not exceed

to exactly that high energies as the dispersion for pure La2CoO4. This might be owed to the structural changes and the amount of disorder in the charge sector with extremely short-ranged La2CoO4-like correlations. However, this difference in the maximum energy clearly proves that we did not simply measure the magnetic signal of a macroscopically phase-separated sample with conventional phase separation. Our polarized neutron study underlines our ndings as the hitherto unknown high-energy magnon band clearly belongs to an out-of-plane excitation, whereas the outwards dispersing contributions from the neck of the hour-glass are mostly in-plane excitations,

4 NATURE COMMUNICATIONS | 5:5731 | DOI: 10.1038/ncomms6731 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

& 2014 Macmillan Publishers Limited. All rights reserved.

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms6731 ARTICLE

Min Counts (a.u.) Max

45 40 35 30 25 20

Energy

E(meV)

15 10 5

1.0

J

0.2

1+K K 0 (r.l.u.)

0.4

0.6

0.8

H

J

K

Figure 3 | The two hour-glass dispersions (a) Our nano-phase separation model within La5/3Sr1/3CoO4. The doping of additional electrons or Co2 ions (green spheres) at Co3 sites (black spheres) of an ideal checkerboard charge ordered structure of the half-doped AFM parent compound induces strong frustration as we already pointed out in ref. 13. However, due to the large ratio of the exchange interactions J/J0 the non-collinearity of the spin structure can be expected basically in the hole-doped regions (blue) and not within the undoped islands (red) where the spins should be aligned in a much more collinear manner. Neutrons (yellow waves) scattered from these different regions will behave differently. (b) Schematic presentation of the magnetic spectra in La2 xSrxCoO4 together with measured neutron scattering intensity behind. The blue contributions indicate the low-energy excitations

originating from the hole-doped regions (blue areas in a) with small J,0 whereas the high-energy signal arises from excitations within the undoped regions (red areas in a) with large J. Here, the upper dispersion (red) is shifted somewhat to higher energies such that the lower dispersion (blue) becomes clearly visible, although both dispersions should inter-penetrate each other. The yellow waves indicate neutrons that were either scattered from the undoped islands or from the hole-doped regions. (c) One-dimensional mechanical analogy of decoupled oscillations of balls within a one-dimensional chain that are connected with weak (blue) or strong (red) force constants.

see Fig. 2d, and, thus, fully compatible to the La2CoO4 excitation spectrum25. Also the 40% Sr-doped cobaltate exhibits qualitatively the same energy scan at the planar AFM wavevector, see Fig. 2c. This is fully in agreement with our scenario that the spin gap within the undoped islands is responsible for the resonance-like increase of intensity at the neck of the hour-glass, as the spin anisotropy gap energy should be much less doping dependent than a dispersion that starts from different incommensurate magnetic peak positions in samples with different hole-doping. Also the comparison of magnetic intensities of high and low energy magnetic signal for two different hole-concentrations corroborates our whole scenario, as the additional high energy mode is the strongest for low Sr-doping where larger fractions of undoped regimes can be expected, compare Fig. 2e,f. Furthermore, the temperature dependence of the hour-glass dispersion is also consistent with our nano phase separation model, see Fig. 4. The temperature dependence of the high energy excitations is different from the strongly temperature affected low energy part. Fitting the magnetic peaks of the constant-E scans at 5 meV and at 29 meV (shown in Fig. 4b,c) with two symmetric gaussians indicates a different temperature dependence in both energy regions. In our model, this could be explained by a different decrease of the spin wave stiffness of upper and lower dispersions with rising temperature, see Fig. 4d. Furthermore, the peak widths are unchanged with temperature, see inset in Fig. 4d. This unaffected peak width reects an almost temperature independent dynamic correlation length. In our nano phase separation model this would be naturally implied by the limiting factor of the limited static domain size of hole-doped checkerboard charge ordered and undoped islands. For the

low-energy excitations this dynamic correlation length would be of the order of B0.6 nm that would nicely correspond to the measured average correlation length of the static checkerboard charge order which amounts to B0.7 nm.

DiscussionThe physical reason for the occurrence of these two distinct sets of dispersions that we schematically present in Fig. 3b can be understood as well. The magnetic exchange interaction J0 between two Co2 ions in the hole-doped (checkerboard charge ordered) regions is much smaller than the comparably large magnetic interactions J within the undoped islands of La2CoO4 type (J44J0), see Fig. 3a. Thus, the excitations within the hole-doped regions stay at low energies, whereas the excitations of the undoped islands (red area in Fig. 3a) exceed to high energies. The reason that the excitations within both types of regions are that strongly de-coupled is the huge difference in the exchange couplings J and J0 that is of the order of B9.7 and B1.4 meV, respectively17,25, and also the existence of a sizeable spin gap within the phase with larger J. For a simplied mechanical analogy see Fig. 3c. The neck of the hour-glass dispersion results basically from the existence of a spin gap within the undoped islands. Note that such a spin gap of these dimensions has been observed in La2CoO4 (ref. 25), see Fig. 2a (dark yellow dots).

Our spin wave calculations based on our nano-phase separation scenario, which we modelled with Monte Carlo simulations using the Metropolis algorithm, are able to reproduce all basic features of the hour-glass dispersion including our additional high-energy mode and, thus, corroborate our experimentally found nano phase

NATURE COMMUNICATIONS | 5:5731 | DOI: 10.1038/ncomms6731 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 5

& 2014 Macmillan Publishers Limited. All rights reserved.

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms6731

E-Scan

5 meV diagonal

10 K 100 K

29 meV vertical

3

20 K

10 K 60 K70 K80 K

90 K 100 K

Counts (a.u.)

I (a.u.)

2

1

0

0

0.0

30 K40 K50 K

10 K100 K

I (a.u.)

10 20 30 0.0 0.5 1.0

40

E (meV)

K (r.l.u.)

0.05

Width (r.l.u.)

0.2

0.0

0 100

Temp. (K)

10 K100 K

Peak shift (r.l.u.)

5 meV 29 meV

0.00

0.5

H (r.l.u.) Temperature (K)

1.0

0 50 100

Max Min

E(meV)

40

35

30

25

20

15

10

5

40

35

30

25

20

15

10

5

E(meV)

E(meV)

I (a.u.)

E(meV)

0.0 0.5 1.0

K (r.l.u.)

0.0 0.5 1.0

K (r.l.u.)

Figure 4 | Temperature dependence of the hour-glass dispersion (a) The temperature dependence of energy scans at the planar AFM wavevector. (b) Temperature dependence of constant-E scans at 29 meV scanning vertical through the planar AFM wavevector. (c) Temperature dependence of constant-E scans at 5 meV scanning diagonal through the magnetic satellite positions. (d) Temperature dependence of the magnetic peak positions relative to the peak positions at lowest temperature (10 K) if two Gaussians are tted to the data in b,c. The inset shows the corresponding peak widths. (eh) Colour-contour plots of neutron scattering intensity measured as a function of energy and momentum transfer for low (10 K: e,g) and higher temperatures (100 K: f,h). This temperature dependence is consistent with our nano phase separation model (see text). Intensities in bh are multiplied by

k k

1 exp ho = kBT = f2Co Q . Error bars are s.d. (obtained from gaussian ts). Intensity error bars are statistical error bars calculated

by the square root of intensity.

separation model, see Methods section and Fig. 5. First of all, chiral or non-collinear magnetic structures arise due to frustration in analogy to a proposed scenario for LSCO26,27. More important, our calculations are able to reproduce the surprising additional high-energy mode that arises from out-of-plane excitations exactly as observed in our polarized neutron scattering experiments (see also the movie in the Supplementary Material). Furthermore, our simulations indicatein close relation to our mechanical analogy

shown in Fig. 3cthat the hole-doped (checkerboard charge ordered) regions can be only excited at lower energies, whereas at high energies only the undoped islands of type of the antiferro-magnetic parent insulator La2CoO4 can be excited, see Fig. 5.

Hence, our spin wave simulations corroborate our experimental nding of a microscopically split origin of upper and lower part of magnetic excitations within the so-called hour-glass dispersion.

6 NATURE COMMUNICATIONS | 5:5731 | DOI: 10.1038/ncomms6731 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

& 2014 Macmillan Publishers Limited. All rights reserved.

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms6731 ARTICLE

Elastic map

(spin correlations)

1.0

K(r.l.u.)

0.5

Intensity

MinMax

0.0 1.0 1.5 2.0

H (r.l.u.)

(1.5 0.5 0)

~42 meV

Excited spins

coupled with J

50

High energy

magnetic

excitation

H 0.5 0

45

Max

40

35

30

Energy (meV)

MinIntensity

(1.5 0.5 0)

~15 meV

Excited spins

coupled with J

25

Resonance peak

(1.29 0.29 0)

~4 meV

20

15

H H 1 0

10

Excited spins

coupled with J

5

1 2

H (r.l.u.)

Low energy magnetic excitation

Figure 5 | Spin wave simulations (a) Spin correlations within a 30 30

lattice calculated for our nano phase separation model. A non-collinear magnetic structure arises from frustration. Black spheres: nonmagnetic Co3 ions, coloured spheres: Co2 ions, red shaded areas: undoped islands (Co2 ions coupled with at least one J), blue shaded areas: checkerboard charge ordered islands (Co2 ions coupled only with J0). The direction of spins is additionally indicated by a colour coding. (b) Fourier transformation of the non-excited spin congurations in Q-space. The vertical crosses indicate wavevectors expected for an ideal checkerboard charge ordered system and the diagonal crosses indicate the expected incommensurate magnetic peak positions for our nominal hole-concentration according to an experimentally observed linear hole-doping dependence13. (c) Simulated spin excitation spectrum. The basic features of the hour-glass dispersion can be reproduced. (df) Calculated magnetic excitations. The movement of the spins is indicated by plotting strongest displacements of the spins for this spin excitation on top of the static spin conguration such that direction and colour of that spin are differing when spins are involved in the excitation. In addition, the dashed cyan lines mark regions of excited spins. Our theoretical simulations yield a purely out-of-plane character of the novel high-energy spin excitations, see also our Supplementary Video File; (here in d, the out-of-plane oscillation has been articially transformed into an in-plane excitation). As can be seen, the checkerboard charge ordered regions (coupled with J0ooJ) can be excited only at low energies, whereas at high energies only the undoped islands (coupled with J44J0) can be excited. At the crossing of the hour-glass the undoped islands become in resonancehence the occurrence of a resonant peakand exhibit the strongest movements.

Finally, the insulating cobaltates appear to be a copper-free reference system that enables us to understand some most basic underlying mechanisms of hour-glass spectra and that overcomes the obstacle of uctuating charges, thus unraveling a direct connection between the occurrence of hour-glass magnetic spectra and nano phase separated microstructures. Our ndings

might help to understand the observations of isotropic high-energy magnetic excitations in the isostructural La2 xSrxCuO4

(LSCO) materials7,11,28 or the recent observations of similar high-energy magnetic excitations in superconducting LSCO samples and the undoped parent materials probed by resonant inelastic X-ray scattering (RIXS)29,30. Whereas it seems to be reasonable to attribute the spin gap within the undoped islands of the cobaltates to an anisotropy gap like in La2CoO4 (ref. 25), in the cuprates the higher energy scales within the excitation spectra would suggest a different physics. Spin gaps with higher energies are known for singlet-triplet excitations of cuprate materials due to nite size effects of the undoped regions3134.

Methods

Crystal synthesis. A.C.K. has grown the La2 xSrxCoO4 and La2 xSrxNiO4

(LSNO) single crystals by the optical oating zone technique using a CSC four-mirror image furnace following the route described in ref. 13. The crystals are single phase as ascertained by high-resolution powder X-ray diffraction measurements. The single crystals are single domain as ascertained by Laue X-ray and single crystal neutron diffraction techniques. Owing to the measured incommensurability in these neutron measurements the oxygen content in La5/3Sr1/3CoO4d is close to stoichiometric dB0.01. The crystal quality was further

ascertained by WDX measurements indicating close to nominal stoichiometry within our single crystals La1.675(6)Sr0.325(2)CoO4d (averaged over 10 measurements). Finally, we performed a highly accurate single crystal X-ray diffraction measurement at a Bruker D8 VENTURE single-crystal X-ray diffractometer equiped with a Photon large area CMOS detector. The sample has a spherical shape with 95(5) mm diameter. In total, 93,028 reections have been collected using Mo Ka radiation and a bent graphite monochromator for B3 times intensity enhancement. The internal R-value amounts to 2.20% with a redundancy of 50.4. The R and weighted R-value in the renement on F2 amount to 1.66 and4.19% and the goodness of t amounts to 1.27. The resulting structural parameters within space group I4/mmm are z(La/Sr) 0.362234(1), Uiso(La/Sr) 0.005885(2),

Uiso(Co) 0.004560(4), Uiso(O1) 0.009831(25), z(O2) 0.173122(24),

Uiso(O2) 0.020418(36), a 3.8552(3) , c 12.6267(9) , and the renement

of the Sr-occupancy indicates that x amounts to 33.2(2)% that is close to stoichiometric. The nominally 40% Sr doped reference sample was already described in ref. 13 and has a composition of La1.591(3)Sr0.404(1)CoO4.01(1) according to ICP-OES and titration measurements.

Z.W.L. has also grown another La5/3Sr1/3CoO4d single crystal following the route described in ref. 12, which reveals larger additional oxygen contents than the sample grown by the author A.C.K. and was only used for the measurement of the temperature dependence of the magnetic excitations and the mXRD experiments. Comparing the magnetic incommensurabilities the excess of oxygen amounts to d 0.035(1) in this sample. However, the dispersion matches

reasonably well with that of ref. 12, thus indicating a similar incommensurability. Complementary inductively coupled plasma optical emission spectroscopy (ICP-OES) and WDX measurements indicate close to nominal stoichiometry of that sample La1.675(11)Sr0.326(3)CoO4d and La1.663(7)Sr0.330(6)CoO4d, respectively.

Synchrotron measurements. The synchrotron radiation single crystal X-ray diffraction measurements have been performed at beamline P07 at PETRA III at DESY using 100 keV hard X-rays monochromatized by a bent Si(111) double monochromator providing an energy resolution of 0.2%. The signal was analyzed with a Si-Ge(111) gradient crystal analyser and detected with a point detector. The measured samples were very B(0.51.5)3 mm3 sized single crystals. The monitor-normalized (200) reection intensities of La5/3Sr1/3CoO4 and the LSNO reference sample amount to 2.7 1010 photons s 1 and 9.7 109 photons s 1.

The mXRD experiments were carried out at the Coherence Beamline P10 at PETRA III. The X-ray beam is monochromatized using a cryogenically cooled Si(111) double crystal monochromator. An X-ray energy of 7,719 eV with a bandwidth of dE/EB1.4 10 4 was selected. About 1 m upstream of the focusing

optics situated at 85.6 m from the undulator source a slit of a size of 150 75 mm2

selected a highly coherent fraction of the X-ray beam. This collimated coherent X-ray beam was focused using the in-vacuum Beryllium compound refractive lens (CRL) transfocator of the P10 beamline to a size of about 2 2 mm2 at the sample

position 2.2 m downstream of the transfocator centre. A pair of JJ X-ray in-vacuum slits located at 220 and 800 mm upstream of the sample was used to clean up the small angle scattering signal of the CRLs from focused X-ray beam. The incident ux on the sample was about 12 1011 photons s 1. The scattered signal from

the sample was recorded at a sample to detector distance of 25 cm using a Pilatus 300 K detector.

Neutron scattering measurements. For the elastic and inelastic neutron scattering experiments one large La5/3Sr1/3CoO4, La1.6Sr0.4CoO4 or LSNO single crystal has been mounted with an [100]/[010] orientation in the scattering plane. The

NATURE COMMUNICATIONS | 5:5731 | DOI: 10.1038/ncomms6731 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 7

& 2014 Macmillan Publishers Limited. All rights reserved.

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms6731

inelastic neutron scattering experiments at the 1T.1 and 2T spectrometers have been performed in constant-kf operation mode with xed focusing pyrolytic graphite (PG) monochromator and analyser using two PG lters for the suppression of higher order contamination. At the IN8 spectrometer an experimental set-up with double focusing PG monochromator and analyzer as well as two PG lters was used. At the IN22 spectrometer a set-up with xed vertically focusing (125 cm radius) and at horizontal (111) Heusler monochromator and vertically focusing and at horizontal (111) Heusler analyser was used. For polarization analysis a CRYOPAD polarimeter for spherical polarization analysis was installed. The ip-ping ratio amounts to B16. A PG lter was used for suppression of higher order contamination. At the IN12 spectrometer that is situated at the end of a horizontal elliptic focusing guide a double focusing PG monochromator and a horizontally focusing Heusler analyser were used. The incident neutron beam was polarised with a polarising cavity. For polarization analysis a set-up with Helmholtz coils for linear polarization analysis was installed. The ipping ratio amounts to B17.4. A velocity selector was used for elimination of higher order contributions (suppression better than 10 4). The temperature dependence of the magnetic excitations in

La5/3Sr1/3CoO4d was measured at the PUMA spectrometer using a double

focusing PG monochromator and analyser as well as two PG lters for elimination of higher order contamination.

Numerical simulations. The Monte-Carlo simulations for simulating disordered checkerboard charge ordered structures are based on the Metropolis algorithm similar as in ref. 15. Only one potential for Coulomb repulsion and another for steric repulsion of two adjacent Co2 ions with large ionic radii has been considered for our simulation of the charge distribution on a N N lattice (N 30,

50) and statistical/random distribution of the charges at each site with a certain probability (here 40%) in the initial step for about a dozen congurations.

To model the magnetism we introduce AFM exchange interactions between the Co2 ions i,j with spin operators ^Si and ^Sj, respectively (the Co3 ions were treated as nonmagnetic): Hex 12 Pi;j Jij

^Si^Sj. According to our nano-phase separation model, we only assume short-range interactions J(100) 5.8 meV

and J(200) 0.85 meV. To simplify the numerical calculation we consider only

the Kramers ground state doublet (|4) of each Co2 ion. For the easy plane anisotropy of each Co2 ion we assume an in-plane magnetic saturation moment of mx my 2mBh j

^Sx j i 2imBh j

^Sy j i 3mB and an out of plane

saturation moment of mz 2mBh j

^Sz j i 1:3mB. Thus, the exchange

interaction Hex can be projected onto the following spin1/2 Hamiltonian, which is of the same form as that used in ref. 15: H 12 Pi;j Jij Pa x;y;z

mm ^sia^sja. Some

representative N N congurations with N 30 and N 50 were taken as the

magnetic unit cell in a periodic lattice, and the spin structure and spin dynamics were calculated by the Dynamical Matrix diagonalization method35 using the McPhase program package (http://creativecommons.org/licenses/by/4.0/

Web End =http://www.mcphase.de). Finally, the spectra of all calculated congurations have been averaged.

References

1. Tranquada et al. Quantum magnetic excitations from stripes in copper oxide superconductors. Nature 429, 534538 (2004).

2. Matsuda, M. et al. Magnetic dispersion of the diagonal incommensurate phase in lightly doped La2 xSrxCuO4. Phys. Rev. Lett. 101, 197001 (2008).

3. Arai, M. et al. Incommensurate spin dynamics of underdoped superconductor YBa2Cu3O6.7. Phys. Rev. Lett. 83, 608611 (1999).

4. Bourges, P. et al. The spin excitation spectrum in superconducting YBa2Cu3O6.85. Science 288, 12341237 (2000).

5. Christensen, N. B. et al. Dispersive excitations in the high-temperature superconductor La2 xSrxCuO4. Phys. Rev. Lett. 93, 147002 (2004).

6. Hayden, S. M., Mook, H. A., Dai, P., Perring, T. G. & Dogan, F. The structure of the highenergy spin excitations in a high-transition-temperature superconductor. Nature 429, 531534 (2004).

7. Vignolle, B. et al. Two energy scales in the spin excitations of the high-temperature superconductor La2 xSrxCuO4. Nat. Phys. 3, 163167

(2007).8. Xu, G. et al. Testing the itinerancy of spin dynamics in superconducting Bi2Sr2CaCu2O8d. Nat. Phys. 5, 642646 (2009).

9. Lipscombe, O. J., Vignolle, B., Perring, T. G., Frost, C. D. & Hayden, S. M. Emergence of coherent magnetic excitations in the high temperature underdoped La2 xSrxCuO4 superconductor at low temperatures. Phys. Rev.

Lett. 102, 167002 (2009).10. Hinkov, V. et al. Spin dynamics in the pseudogap state of a high-temperature superconductor. Nat. Phys. 3, 780785 (2007).

11. Fujita, M. et al. Progress in neutron scattering studies of spin excitations in high- t c cuprates. J. Phys. Soc. Jpn 81, 011007 (2012).

12. Boothroyd, A. T., Babkevich, P., Prabhakaran, D. & Freeman, P. G. An hourglass magnetic spectrum in an insulating, hole-doped antiferromagnet. Nature 471, 341344 (2011).

13. Drees, Y., Lamago, D., Piovano, A. & Komarek, A. C. Hour-glass magnetic spectrum in a stripeless insulating transition metal oxide. Nat. Commun. 4, 2449 (2013).

14. Lancaster, T. et al. Stripe disorder and dynamics in the hole-doped antiferromagnetic insulator La5/3Sr1/3CoO4. Phys. Rev. B 89, 020405 (2014).

15. Andrade, E. C. & Vojta, M. Disorder, cluster spin glass, and hourglass spectra in striped magnetic insulators. Phys. Rev. Lett. 109, 147201 (2012).

16. Shannon, R. D. Revised effective ionic radii and systematic studies of interatomic distances in halides and chalcogenides. Acta Cryst. A32, 751767 (1976).

17. Helme, L. M. et al. Magnetic order and dynamics of the charge-ordered antiferromagnet La1.5Sr0.5CoO4. Phys. Rev. B 80, 134414 (2009).

18. Chang, C. F. et al. Spin blockade, orbital occupation, and charge ordering in La1:5Sr0:5CoO4. Phys. Rev. Lett. 102, 116401 (2009).

19. Fratini, M. et al. Scale-free structural organization of oxygen interstitials in La2CuO4y. Nature 466, 841844 (2010).

20. Poccia, N. et al. Optimum inhomogeneity of local lattice distortions in La2CuO4y. Proc. Natl Acad. Sci. 109, 1568515690 (2012).

21. Ricci, A. et al. Multiscale distribution of oxygen puddles in 1/8 doped YBa2Cu3O6.67. Scientic Rep. 3, 2383 (2013).

22. Krger, F. & Scheidl, S. Spin dynamics of stripes. Phys. Rev. B 67, 134512 (2003).

23. Vojta, M. & Ulbricht, T. Magnetic excitations in a bond-centered stripe phase: Spin waves far from the semiclassical limit. Phys. Rev. Lett. 93, 127002 (2004).

24. Yao, D. X., Carlson, E. W. & Campbell, D. K. Magnetic excitationsof stripes and checkerboards in the cuprates. Phys. Rev. B 73, 224525 (2006).

25. Babkevich, P., Prabhakaran, D., Frost, C. D. & Boothroyd, A. T. Magnetic spectrum of the two-dimensional antiferromagnet La2CoO4 studied by inelastic neutron scattering. Phys. Rev. B 82, 184425 (2010).

26. Sushkov, O. P. & Kotov, V. N. Theory of incommensurate magnetic correlations across the insulator-superconductor transition of underdoped La2 xSrxCuO4. Phys. Rev. Lett. 94, 097005 (2005).

27. Lscher, A., Milstein, A. I. & Sushkov, O. P. Structure of the spin-glassstate of La2 xSrxCuO4: The spiral theory. Phys. Rev. Lett. 98, 037001037004

(2007).28. Matsuda, M., Granroth, G. E., Fujita, M., Yamada, K. & Tranquada, J. M. Energy-dependent crossover from anisotropic to isotropic magnetic dispersion in lightly doped La1.96Sr0.04CuO4. Phys. Rev. B 87, 054508 (2013).

29. Dean, M. P. M. et al. Persistence of magnetic excitations in La2 xSrxCuO4 from

the undoped insulator to the heavily overdoped non-superconducting metal. Nat. Mater. 12, 10191023 (2013).30. Dean, M. P. M. et al. Spin excitations in a single La2CuO4 layer. Nat. Mater. 11, 850854 (2012).

31. Lorenzo, J. E. et al. Macroscopic quantum coherence of the spin tripletin the spin-ladder compound Sr14Cu24O41. Phys. Rev. Lett. 105, 097202 (2010).

32. Deng, G. et al. Spin-gap evolution upon Ca doping in the spin-ladder series Sr14 xCaxCu24O41 studied by inelastic neutron scattering. Phys. Rev. B 88,

014504 (2013).33. Eccleston, R. S., Azuma, M. & Takano, M. Neutron-scattering and susceptibility study of spin chains and spin ladders in (Sr0.8Ca0.2)14Cu24O41. Phys. Rev. B 53,

R14721R14724 (1996).34. Azuma, M., Takano, M. & S. Eccleston, R. Disappearance of the spin gapin a Zn-doped 2-leg ladder compound Sr(Cu1 xZnx)2O3. J. Phys. Soc. Jpn 67,

740743 (1998).35. Rotter, M., Duc Le, M., Bothroyd, A. T. & Blanco, J. A. Dynamical matrix diagonalization for the calculation of dispersive excitations. J. Phys.: Condens. Matter 24, 213201 (2012).

36. Barabsi, A.-L. & Albert, R. Emergence of scaling in random networks. Science 286, 509512 (1999).

Acknowledgements

We thank O. Stockert, C. Schssler-Langeheine, P. Thalmeier and L.H. Tjeng for helpful discussions. We thank A. Todorova for assistance. We thank L.P. Regnault for help at the IN22 spectrometer. For sample-characterization using powder X-ray diffraction and EDX/WDX measurements we thank H. Borrmann, Y. Prots and S. Hckmann as well as the team of U. Burkhardt, respectively.

Author contributions

A.C.K. has conceived all experiments, performed the single-crystal X-ray diffraction experiments, performed the Monte-Carlo simulations, analysed the results and written the article. A.C.K. and Z.W.L. synthesized all studied materials. A.C.K., Y.D., Z.W.L., O.G. and U.R. have performed the hard X-ray experiments. A.R., Z.W.L. and M.S. have performed the mXRD measurements. Y.D., Z.W.L., W.S., D.L., J.P.C. and A.P. have measured the novel high-energy magnetic excitations with unpolarized and polarized

8 NATURE COMMUNICATIONS | 5:5731 | DOI: 10.1038/ncomms6731 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

& 2014 Macmillan Publishers Limited. All rights reserved.

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms6731 ARTICLE

How to cite this article: Drees, Y. et al. Hour-glass magnetic excitations induced by nanoscopic phase separation in cobalt oxides. Nat. Commun. 5:5731doi: 10.1038/ncomms6731 (2014).

This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the articles Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/

Web End =http://creativecommons.org/licenses/by/4.0/

neutrons. A.C.K., Z.W.L. and O.S. have performed the temperature-dependent neutron scattering experiments. M. R. and A.C.K. have performed the spinwave calculations.

Additional information

Supplementary Information accompanies this paper at http://www.nature.com/naturecommunications

Web End =http://www.nature.com/ http://www.nature.com/naturecommunications

Web End =naturecommunications

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

Reprints and permission information is available online at http://npg.nature.com/reprintsandpermissions/

Web End =http://npg.nature.com/ http://npg.nature.com/reprintsandpermissions/

Web End =reprintsandpermissions/

NATURE COMMUNICATIONS | 5:5731 | DOI: 10.1038/ncomms6731 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 9

& 2014 Macmillan Publishers Limited. All rights reserved.