Spontaneous decays of magneto-elastic excitations in non-collinear antiferromagnet (Y,Lu)MnO3

Abstract

Magnons and phonons are fundamental quasiparticles in a solid and can be coupled together to form a hybrid quasi-particle. However, detailed experimental studies on the underlying Hamiltonian of this particle are rare for actual materials. Moreover, the anharmonicity of such magnetoelastic excitations remains largely unexplored, although it is essential for a proper understanding of their diverse thermodynamic behaviour and intrinsic zero-temperature decay. Here we show that in non-collinear antiferromagnets, a strong magnon-phonon coupling can significantly enhance the anharmonicity, resulting in the creation of magnetoelastic excitations and their spontaneous decay. By measuring the spin waves over the full Brillouin zone and carrying out anharmonic spin wave calculations using a Hamiltonian with an explicit magnon-phonon coupling, we have identified a hybrid magnetoelastic mode in (Y,Lu)MnO₃ and quantified its decay rate and the exchange-striction coupling term required to produce it.

ARTICLE

Received 3 Jun 2016 | Accepted 8 Sep 2016 | Published 19 Oct 2016

Joosung Oh1,2, Manh Duc Le1,2, Ho-Hyun Nahm1,2, Hasung Sim1,2, Jaehong Jeong1,2, T.G. Perring3, Hyungje Woo3,4, Kenji Nakajima5, Seiko Ohira-Kawamura5, Zahra Yamani6, Y. Yoshida7, H. Eisaki7,S.-W. Cheong8, A.L. Chernyshev9 & Je-Geun Park1,2

Magnons and phonons are fundamental quasiparticles in a solid and can be coupled together to form a hybrid quasi-particle. However, detailed experimental studies on the underlying Hamiltonian of this particle are rare for actual materials. Moreover, the anharmonicity of such magnetoelastic excitations remains largely unexplored, although it is essential for a proper understanding of their diverse thermodynamic behaviour and intrinsic zero-temperature decay. Here we show that in non-collinear antiferromagnets, a strong magnonphonon coupling can signicantly enhance the anharmonicity, resulting in the creation of magnetoelastic excitations and their spontaneous decay. By measuring the spin waves over the full Brillouin zone and carrying out anharmonic spin wave calculations using a Hamiltonian with an explicit magnonphonon coupling, we have identied a hybrid magnetoelastic mode in (Y,Lu)MnO3 and quantied its decay rate and the exchange-striction coupling term required to produce it.

1 Center for Correlated Electron Systems, Institute for Basic Science, Seoul 08826, Korea. 2 Department of Physics and Astronomy, Seoul National University, Kwanak-Gu, Seoul 08826, Korea. 3 ISIS Facility, STFC Rutherford Appleton Laboratory, Didcot OX11 0QX, UK. 4 Department of Physics, Brookhaven National Laboratory, Upton, New York 11973, USA. 5 Materials and Life Science Division, J-PARC Center, Japan Atomic Energy Agency, Tokai, Ibaraki 319-1195, Japan.

6 Chalk River Laboratories, National Research Council, Chalk River, Ontario, Canada K0J 1J0. 7 National Institute of Advanced Industrial Science and Technology, Tsukuba, Ibaraki 305-8565, Japan. 8 Department of Physics and Astronomy, and Rutgers Center for Emergent Materials, Rutgers University, Piscataway, New Jersey 08854, USA. 9 Department of Physics and Astronomy, University of California, Irvine, California 92697, USA. Correspondence and requests for materials should be addressed to J.-G.P. (email: mailto:[email protected]

Web End [email protected] ).

NATURE COMMUNICATIONS | 7:13146 | DOI: 10.1038/ncomms13146 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 1

DOI: 10.1038/ncomms13146 OPEN

Spontaneous decays of magneto-elastic excitations in non-collinear antiferromagnet (Y,Lu)MnO3

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms13146

Spin and lattice constitute two of the four fundamental degrees of freedom in the solid: the other two are charge and orbital. In linearized models that account for many

current understandings of the solid, excitations of spin and lattice, magnon of spin waves and phonon of lattice vibrations, are two principal examples of such quasiparticles1. Although there have been some experimental observations of cross-coupling between magnon and phonon, it is rather rare to actually observe and, more importantly, quantify the magnonphonon coupling for real materials2,3. Nevertheless, it is essential for a proper understanding of their diverse thermodynamic behaviour and intrinsic zero-temperature decay4. Furthermore, it is generally believed that the magnon phonon coupling is important for materials such as multiferroic compounds, geometrically frustrated systems, spin-Peierls systems and Invar materials, to name only a few58.

Quasiparticles like magnon and phonon, the cornerstone of modern condensed matter physics, are fundamentally the byproducts of linearized theories that ignore all the higher-order terms than quadratic terms and neglect any conceivable interaction among the quasiparticles themselves. As such, they are considered to be stable, except for very few exceptions. For example, for classical spin systems without strong quantum uctuations, magnon breakdown is thought to be unlikely for most of purposes. Therefore, observing and understanding how the otherwise stable quasiparticles break down in these unusual cases are the central theme of condensed matter physics.

One route leading to the breakdown of magnon and phonon is the cubic anharmonicities. Despite the general belief that this nonlinear magnon(-magnon) interaction is rather weak in real materials, recent insights gained mainly through theoretical studies suggest that things should be drastically different for certain cases, namely for non-collinear antiferromagnetic structures4. Unlike collinear magnetic structures that forbid the cubic anharmonicities, it was shown that such interaction is in principle allowed for non-collinear magnetic structures, such as the canonical 120 spin pattern in a two-dimensional (2D) triangular lattice. There have since been several experimental reports9,10 supporting these theoretical postulates. Nevertheless, the full details of the nonlinear interaction still need to be worked out, especially from experiments. We should also point out that this non-collinear magnetic order, in principle, allows a hitherto forbidden magnonphonon coupling that has been less recognized among the community: the rst-order variation of the exchange energy with respect to transverse spin uctuations is non-zero for a non-collinear magnetic order11. As the O(3) symmetry is completely broken in the non-collinear ordered ground state of spins, coupling to phonons necessarily generates a coupling in which a magnon can convert directly into a phonon and vice versa. This is in contrast to the spinlattice coupling in more conventional, collinear magnets, where the coupling usually respects the parity, which necessarily conserves the number of magnons or allows creation (annihilation) of them only in pairs12.

Here we report the direct observation of the magnonphonon coupling and the spontaneous decay of magneto-elastic excitations in the triangular antiferromagnets (Y,Lu)MnO3.

The full magnetic excitation spectra measured by inelastic neutron-scattering experiments show clear deviations from the linear spin wave theory without the magnonphonon coupling: an additional mode at high energies and the downward shift of the bottom mode at the Brillouin zone boundary. This is the most direct and stark evidence of the linear coupling of magnons and phonons, which, in turn, leads to enhancement of the anharmonic effects. We demonstrate that these experimental

anomalies can only be fully resolved by incorporating the magnonphonon coupling and carrying out the nonlinear spin wave analysis. We further reveal that the magneto-elastic excitation leads to signicant broadening of the magnon spectra at the zone boundaries, originating from the decay of the magneto-elastic excitations into the two-magnon continuum.

ResultsFailure of standard spin waves calculation. Hexagonal rare-earth manganite RMnO3 represents a good model system for geometrical frustration on a 2D triangular lattice: the nearest-neighbour antiferromagnetic interaction between S 2 (Mn3 )

spins dominate, whereas the interlayer interaction is relatively weak13,14. We should note that it also exhibits a very large spinlattice coupling when it becomes magnetically ordered5. In a more recent work, we reported experimental evidence for a spontaneous magnon decay and a remarkably large spatial anisotropy for Mn3 ions in the exchange interactions for

LuMnO3. This was attributed to a structural distortion in which groups of three Mn atoms become more closely bound9, such that the intratrimer J1 exchange constants may differ from the intertrimer J2 (see Fig. 1). Similar interpretations were shared by other groups too15,16. However, we should note that the large J1/J2 ratio of 6.4 obtained from tting the data is inconsistent with the value of 1.15, obtained from rst-principles calculations17 that used the experimental atomic positions as reported by neutron diffraction5. This realization prompted us to handle both magnon and phonon, as well as their cross-coupling on an equal footing and to go beyond the linear spin wave analysis.

Figure 1 shows the spin waves measured at the MAPS beamline of the ISIS facility together with the theoretical dispersion relation calculated from the spin Hamiltonian by a linear spin wave theory using the following parameters (see Supplementary Note 2): for YMnO3: J1 4, J2 1.8, D1 0.28, D2 0.02 meV;

for Y0.5Lu0.5MnO3: J1 12.5, J2 0.97, D1 0.18, D2 0.018

meV; for LuMnO3: J1 9, J2 1.4, D1 0.28, D2 0.02 meV.

Despite the apparent success of the linear spin wave calculations, there lies a critical failure: rst, the downward curvature along the AB direction and, second, the additional peaks at B19 meV indicated by a red box in Fig. 1ce. However, most importantly, here we have to use an unphysically large J1/J2 ratio, to explain the additional high-energy peaks. Apart from the large J1/J2 ratio, this analysis of the linear spin waves has another drawback: which is that the calculated dynamical structure factor using the linear spin wave theory as shown in Fig. 1fh always produces stronger intensity at the top mode of the spin waves than at the middle one, in marked contrast with the experimental data.

This discrepancy requires us to adopt a radically different approach to explain the measured full spin waves and go beyond the standard linear spin wave theory. One clue for how to address this problem can be taken from the physical properties: for example, our earlier neutron diffraction data revealed a giant spinlattice coupling for (Y,Lu)MnO3 (ref. 5). This observation was subsequently corroborated by independent measurements on rare-earth hexagonal RMnO3 (refs 1825). More importantly, ultrasound measurements on YMnO3 found marked softening in C11 and C66, supporting our view that there is a strong in-plane deformation below TN26. This observation naturally indicates the importance of magnon and in-plane phonon coupling. This conclusion is also backed up by theoretical calculations17,27.

Calculation of magneto-elastic excitations. Following this idea of a large spinlattice coupling in RMnO3, we took a rst-principles approach to the magnonphonon coupling.

2 NATURE COMMUNICATIONS | 7:13146 | DOI: 10.1038/ncomms13146 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms13146 ARTICLE

a

O Mn

b

C A

J1

Intensity (a.u.)

1

0.5

0

J2

O

B

a

D

b

YMnO3

Y0.5Lu0.5MnO3

LuMnO3

c d e

Intensity (a.u.)

20

Energy (meV) Energy (meV)

15

10

5

25

20

15

10

5

0A B

MK MK MK MK MK

MK

10

89

7

6

5

4

3

2

f

g

h

1

0.25

O M M M M M M

C B D C A B O C B D C A B O C B D C

Figure 1 | Magnetic excitation spectra in (Y,Lu)MnO3. (a) A Mn-O layer in RMnO3 forming a distorted 2D triangular antiferromagnet. (b) Inelastic neutron scattering data on LuMnO3 summed over an energy window of [19.5, 20.5] meV. The arrows in b indicate the reciprocal points where the data shown in ce are cut. (ce) The inelastic neutron-scattering data along the high symmetric directions (red circle and contour map) and tted dispersion (black solid curve) for (c) YMnO3, (d) Y0.5Lu0.5MnO3 and (e) LuMnO3 calculated by linear spin wave theory. (fh) Calculated dynamical spin structure factors using the minimal spin Hamiltonian, equation (4) in the Supplementary Note 2 for (f) YMnO3, (g) Y0.5Lu0.5MnO3 and (h) LuMnO3. For our simulations, we used the momentum and energy resolution of the instrument at the elastic line.

First, we construct the following full model Hamiltonian, which couples the in-plane manganese vibrations directly to the spin system:

HHspin X90i1 oi

by

eOiji ui eOijj uj

Si Sj

1

where eOiji denotes the unit vector connecting the i-th manganese atom and the neighbouring oxygen atoms between the i-th and j-th manganese atoms as shown in Supplementary Fig. 2, is the exchange striction,@J=@r, which is naturally made into a

dimensionless exchangestriction constant a 2d=J28, and d is

the MnO bond length at the equilibrium. Therefore, our Hamiltonian takes into account the modulations of the MnO bond length as a function of Mn displacements.

Before going into detailed discussion, we would like to make a general remark on the related issue. In cases when the spin rotational symmetry is broken completely in the ground state, that is, when the spin structure is non-collinear, the Heisenberg term of the Hamiltonian provides a coupling of the transverse and longitudinal modes on neighbouring sites. That is, in terms of the local site-dependent preferred spin direction of the ordered state, the coupling terms take the form of the type SziSxj and so on.

In the magnon language, they are quantized into the odd terms, producing linear a aw and cubic (awaa and so on) contribu

tions. In equilibrium, the linear magnon term must vanish, leaving

the anharmonic cubic magnon coupling to be the sole outcome, which is important for magnon decays. However, in the presence of coupling to phonons, the linear a aw terms is activated, as the

local atomic displacements (ui) violate the equilibrium conditions locally, hence the direct coupling of magnons and phonons.

To calculate the full dispersions of all 90 phonon modes for the unit cell with six formula units, we used a rst-principles density functional theory (DFT). We show the full phonon dispersion curves for the three compounds as dashed lines in Fig. 2. We note that the calculated phonon density of states (DOS) is in good agreement with the phonon spectra we measured using powder YMnO3 and LuMnO3 at the AMATERAS beamline of J-PARC (see Supplementary Fig. 1). We then calculated the dynamical spin structure factor within the linear approximation by using the full Hamiltonian above with the explicit magnonphonon coupling: we used the dimensionless exchangestriction coefcients of aB1620.

The exchangestriction constant can also be estimated by using the pressure dependence of the crystal structure and the antiferromagnetic transition reported for YMnO3 (refs 29,30).

Using the experimental data reported in refs 9,10, we came to an estimate of the dimensionless exchange striction a of 14 with the following formula: a d P0@TN P=@PTNP0@d P=@P. Here, d is average MnO bond length, which is approximately one-third of lattice constant a. The experimental parameters used in our estimate are @TN P

@P

k bk X

ij

3 K GPa 1 (ref. 29) and

@d P

@P

0:0057 GPa 1

NATURE COMMUNICATIONS | 7:13146 | DOI: 10.1038/ncomms13146 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 3

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms13146

a

25

20

Energy (meV) Energy (meV) Energy (meV)

15

10

5

0 25

20

15

10

5

0 25

S(q,)

(a.u.) 10

b

c

20

15

10

5

0A B

M K

O

Si Sj:

2

First, we calculate the dynamical spin structure factor by using a standard method. As shown in Fig. 3a, despite the simplication, the calculation results reproduce well the experimental intensity along the CBD direction. For the calculations, we used the following set of the parameters: for YMnO3: J 2.7 meV, g 0.93, :o0 17.5 meV, a 7.2; for

LuMnO3: J 3.2 meV, g 0.93, :o0 19.5 meV, a 8. Next, the

decay rate of the high-energy mode was calculated using an anharmonic spin wave theory within the 1/S approximation. The calculated results summarized in Fig. 3 also show the signicant linewidth broadening for the top mode near the B and D points only for LuMnO3, consistent with the experimental results.

The reason for this is that in LuMnO3, the combination of the higher energies of the magnon and magnetoelastic modes means that more decay channels, including two quasiparticle emission, are kinematically allowed. The different decay channels have different boundaries in the reciprocal space, which also corresponds to logarithmic singularities in the decay rate33, and the largest broadening is observed at momentum transfers where the single-quasiparticle dispersion crosses these boundaries, such as at the B point. In the case of YMnO3, no such crossing occurs;

thus, there are fewer decay channels available explaining why the observed linewidths remain narrow. Here we should stress that the single-magnon branches do not cross the line of singularities, whereas the magnetoelastic mode does; thus, a pure magnon decay is forbidden in this case. Similarly, the intrinsic decay rate of phonons is usually small due to a weak cubic anharmonicity. Thus, the strong hybridization of magnons and phonons provides a new mechanism to enhancing the magnon decays.

DiscussionOur studies using an exchange-striction model indicates that the deviations from the linear spin wave theory must be a common feature for other triangular antiferromagnets with a non-collinear magnetic order. Looking beyond the 2D triangular antiferromagnets, we believe that the idea of magnonphonon coupling can also be important in the studies of a wide variety of

o0 X

k

by

k bk X

ij

eO i ui eO

j

uj

C B D C

M M K

98 7 6 5

4

3

2

1

0.25

M

Figure 2 | Calculated dynamical structure factor of magneto-elastic excitation. The dynamical spin structure factor calculated from the full Hamiltonian, equation (1) (contour map) by taking into account the magnonphonon coupling: the phonon dispersion curves (dashed lines) and the magnon dispersion without the coupling (solid lines) for (a) YMnO3,

(b) Y0.5Lu0.5MnO3 and (c) LuMnO3.

(ref. 30). Below we point out that our own data for the magnon excitation spectrum imply the value of the magneto-striction of the same order. This should be contrasted with the cuprates family, where estimates for an equivalent quantity are substantially smaller, aB27(refs 31,32).

By comparing with the experimental data, we obtained the best tting results with the following sets of the parameters: J1 J2 2.5 meV, D1 0.28 meV, D2 0.02 meV and a 16

for YMnO3; J1 J2 2.7 meV, D1 0.28 meV, D2 0.02 meV

and a 20 for Y0.5Lu0.5MnO3; and J1 J2 3 meV,

D1 0.28 meV, D2 0.02 meV and a 16 for LuMnO3

(see Supplementary Table 1). We note that according to the DFT calculations17 the relative difference between J1 and J2 is theoretically about 1020% at maximum. Therefore, we judge that this choice of J1 J2 in our analysis is good enough to

capture the essential underlying physics of (Y,Lu)MnO3, which is the magnonphonon coupling.

The results shown as colour contour plot in Fig. 2 reproduce the overall features of the observed spectra of the experimental data in Fig. 1. It clearly shows that the high-energy signals located at 18B20 meV come from a direct coupling between the magnon and the optical phonons, that is, a magneto-elastic mode. The intensity of this high energy magneto-phonon modes becomes stronger for Lu-enriched compounds due to the larger Mn phonon DOS at this energy for the Lu-enriched compounds, consistent with the experimental results shown in Fig. 1. This conclusion on the relevance of the magnonphonon coupling for RMnO3 is also supported by the fact that our calculated spin waves successfully explains the downward curvature of the bottom magnon branch along the AB direction.

In fact, an estimate of the exchange striction from the splitting of the high-energy hybrid modes in Fig. 2 yields the values in the same range, aB1020 (see Supplementary Note 3). We also note that our polarized neutron-scattering data are in good agreement with our calculations (see Supplementary Fig. 4 and Supplementary Note 4).

Spontaneous decay of the hybrid mode. In addition to the magnonphonon hybridization, non-collinear spin structures allow three magnon interactions as discussed above, which can lead to spontaneous magnon decay into two magnon states when the kinematic conditions are satised. The magneto-elastic excitations have, by denition, both magnon and phonon characters. Therefore, the above mechanism can also lead to the decay of magneto-elastic excitations inside the two quasi-particles continuum of magnon. Indeed, we observe signicant broadening of the top mode in LuMnO3 near the B and D points as shown in

Fig. 3 and Supplementary Fig. 3: less strong broadening has been seen for other two compounds.

To calculate the decay rate directly and compare it with the experimental data, we simplify the problem by assuming a dispersionless optical phonon mode at B20 meV, where the strongest coupling has been observed in our data. Next, our model Hamiltonian reads as follows:

HJ X

ij

SxiSxj SyiSyj gSziSzj

4 NATURE COMMUNICATIONS | 7:13146 | DOI: 10.1038/ncomms13146 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms13146 ARTICLE

a

b

C

FWHM

(meV)

FWHM

(meV)

FWHM (meV)

Energy (meV) Energy (meV)

2

0

25

20

15

10

5

0 25

20

15

10

5

0C B DC B D 0

4

5

4

3

2

1

0

2

1

1.5

0.5

Intensity (a.u.)

4

3.5

3

O B

2.5

c

C

2

1.5

1

0.5

O B

Figure 3 | Linewidth broadening of magneto-elastic excitation. (a) The neutron-scattering data along the CBD direction (middle) and the calculated dynamical structure factor within the linear spin wave theory (bottom). Observed linewidth broadening of the top mode (square) together with the calculated result from the 1/S approximation (line) and the experimental resolution (dashed line) (top) for YMnO3 (left) and LuMnO3 (right). (b) Observed linewidth of the top mode for LuMnO3 (contour map) and (c) calculated intrinsic broadening of the top mode using the model Hamiltonian, equation (2).

The experimental linewidth of the top mode was estimated by using multi-Gaussian functions and the typical results are shown in Supplementary Fig. 3.

the 3d transition-metal magnetic compounds. For example, similar analysis might shed a light on the investigation of spin phonon coupling mechanism and anharmonic effects in many other important non-collinear magnets that exhibit a rather large spinlattice coupling, such as spinel34 and invar materials35, which should be dominated by the exchangestriction as discussed in ref. 36. This is in contrast to the previous studies on magnonphonon coupling focusing on materials with strong spinorbit coupling such as rare-earth elements37.

To summarize, we mapped out the spin waves and phonon excitations of (Y,Lu)MnO3 over the Brillouin zone. By carrying out the spin wave calculations using the full Hamiltonian with both magnons and phonons on an equal footing and their coupling, we have not only demonstrated in our inelastic neutron-scattering data a clear sign of magnonphonon coupling, but also have quantied the coupling strength directly. Our work provides the rare experimental test and quantication of magnonphonon coupling in real materials and opens a new window of opportunities in other materials such as 2D triangular lattice and other frustrated systems, where such couplings, hitherto hidden, have been long suspected.

Methods

Sample preparation. We synthesized powder samples using a solid-state reaction method following the recipe as described in the literature38. We then grew single crystals of Y1 xLuxMnO3 (with typical size of 5 5 40 mm3) by using a

commercial optical oating zone furnace (Crystal Systems, Japan). Our subsequent powder and single crystal XRD conrmed that all our samples are prepared in high quality. We also measured the bulk properties (susceptibility and heat capacity) of all the samples to further conrm the quality by using a commercial set-up (MPMS5XL and PPMS9, Quantum Design USA).

Inelastic neutron scattering. Inelastic neutron-scattering experiments have been performed on single crystal samples using the MAPS time-of-ight spectrometer at

ISIS, UK, and a triple axes spectrometer with polarization analysis at Chalk River, Canada. In the time-of-ight experiments, incident energies were chosen at 40 meV for LuMnO3, 35 meV for Y0.5Lu0.5MnO3 and 30 meV for YMnO3 to adjust to the

slight variations in the energy scales for each samples. The chopper frequency has been set to 250 Hz, which gave us a full width at half maximum energy resolution of 0.43B1.36 meV depending on the energy transfer. The measurements have been performed at 4 K for YMnO3 and Y0.5Lu0.5MnO3, and 13 K for LuMnO3. We used the Horace programme for our data analysis39. In the triple axes spectrometer experiments, spin-polarized neutrons have been produced by using vertical focusing Heusler monochromator and analyser with the nal energy xed at13.7 meV. To measure phonon DOS, inelastic neutron-scattering experiments have also been performed on the powder samples with the incident energy of 42 meV using the AMATERAS beamline at J-PARC, Japan.

Theoretical calculations. We carried out rst-principles calculations of phonon using a DFT U method with U 4 eV. We used the PHONOPY code based on

the force constant method40. In addition, the force constants were constructed by means of a supercell approach based on the density functional perturbation theory41, implemented in the VASP code42. Detailed discussion is given in the Supplementary Note 1.

For the spin waves calculations, we used a rotating framework with the direction of easy axis anisotropy being rotated from parallel to perpendicular to the crystallographic axes. To make our calculations simpler and transparent, we ignored the interlayer exchange coupling as it is known to be more than 100 times smaller than the in-plane coupling9,14. Using this approximation, we can have the following minimal Hamiltonian:

J1

intra

X Si Sj J2 X

inter

Si Sj D1 X

i

Szi

2; 3

where J1 and J2 represent intra- and inter-trimer exchange constants and D1 and D2 are two magnetic anisotropies. We then calculated the spin wave dispersion using the standard linear spin wave theory43. We give detailed description of our spin wave calculations for the full Hamiltonian with the magnon and phonon coupling and the magnonmagnon nonlinear interaction in the Supplementary Notes 2 and 3.

Data availability. All relevant data that support the ndings of this study are available from the corresponding author on request.

2 D2 X

i

Sni

NATURE COMMUNICATIONS | 7:13146 | DOI: 10.1038/ncomms13146 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 5

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms13146

References

1. Landau, L. D. The theory of a Fermi liquid. Sov. Phys. JETP 3, 920925 (1957).2. Pimenov, A. et al. Possible evidence for electromagnons in multiferroic manganites. Nat. Phys. 2, 97100 (2006).

3. Ogawa, N. et al. Photodrive of magnetic bubbles via magnetoelastic waves. Proc. Natl Acad. Sci. USA 112, 89778981 (2015).

4. Zhitomirsky, M. E. & Chernyshev, A. L. Colloquium: spontaneous magnon decays. Rev. Mod. Phys. 85, 219242 (2013).

5. Lee, S. et al. Giant magneto-elastic coupling in multiferroic hexagonal manganites. Nature 451, 805809 (2008).

6. Buschow, K. H. J. (ed.) Handbook of Magnetic Materials Vol.13 (Elsevier, Amsterdam, 2001).

7. Hase, M., Terasaki, I. & Uchinokura, K. Observation of the spin-Peierls transition in linear Cu2 (spin-1/2) chains in an inorganic compound

CuGeO3. Phys. Rev. Lett. 70, 36513654 (1993).8. Kim, D. J. The electron-phonon interaction and itinerant electron magnetism. Phys. Rep. 171, 129229 (1988).

9. Oh, J. et al. Magnon breakdown in a two dimensional triangular lattice Heisenberg antiferromagnet of multiferroic LuMnO3. Phys. Rev. Lett. 111, 257202 (2013).

10. Toth, S. et al. Magnetic soft modes in the distorted triangular antiferromagnet a-CaCr2O4. Phys. Rev. Lett. 109, 127203 (2012).

11. Kim, J. H. & Han, J. H. Coupling of phonons and spin waves in a triangular antiferromagnet. Phys. Rev. B 76, 054431 (2007).

12. Dixon, G. S. Lattice thermal conductivity of antiferromagnetic insulators. Phys. Rev. B 21, 28512864 (1980).

13. Sim, H., Oh, J., Jeong, J., Le, M. D. & Park, J.-G. Hexagonal RMnO3: a model system for 2D triangular lattice antiferromagnets. Acta. Cryst. B 72, 319 (2016).

14. Fabrges, X. et al. Spin-lattice coupling, frustration, and magnetic order in multiferroic RMnO3. Phys. Rev. Lett. 103, 067204 (2009).

15. Sato, T. J. et al. Unconventional spin uctuations in the hexagonal antiferromagnet YMnO3. Phys. Rev. B 68, 014432 (2003).

16. Lewtas, H. J., Boothroyd, A. T., Rotter, M. & Prabhakaran, D. Magnetic excitations in multiferroic LuMnO3 studied by inelastic neutron scattering.

Phys. Rev. B 82, 184420 (2010).17. Solovyev, I. V., Valentyuk, M. V. & Mazurenko, V. V. Magnetic structure of hexagonal YMnO3 and LuMnO3 from a microscopic point of view. Phys. Rev. B 86, 054407 (2012).

18. Litvinchuk, A. P., Iliev, M. N., Popov, V. N. & Gospodinov, M. M. Raman and infrared-active phonons in hexagonal HoMnO3 single crystals: magnetic ordering effects. J. Phys. Condens. Matter 16, 809819 (2004).

19. Souchkov, A. B. et al. Exchange interaction effects on the optical properties of LuMnO3. Phys. Rev. Lett. 91, 027203 (2003).

20. Dela Cruz, C. et al. Strong spin-lattice coupling in multiferroic HoMnO3: thermal expansion anomalies and pressure effect. Phys. Rev. B 71, 060407R (2005).

21. Park, J. et al. Doping dependence of spin-lattice coupling and two-dimensional ordering in multiferroic hexagonal Y1 xLuxMnO3 (0oxo1). Phys. Rev. B 82,

054428 (2010).22. Chatterji, T., Ouladdiaf, B., Henry, P. F. & Bhattacharya, D. Magnetoelastic effects in multiferroic YMnO3. J. Phys. Condens. Matter 24, 336003 (2012).

23. Petit, S. et al. Spin phonon coupling in hexagonal multiferroic YMnO3. Phys. Rev. Lett. 99, 266604 (2007).

24. Sharma, P. A. et al. Thermal conductivity of geometrically frustrated, ferroelectric YMnO3: extraordinary spin-phonon interactions. Phys. Rev. Lett.

93, 177202 (2004).25. Pailhs, S. et al. Hybrid Goldstone modes in multiferroic YMnO3 studied by polarized inelastic neutron scattering. Phys. Rev. B 79, 134409 (2009).

26. Poirier, M. & Lalibert, F. Magnetoelastic coupling in hexagonal multiferroic YMnO3 using ultrasound measurements. Phys. Rev. B 76, 174426 (2007).

27. Varigon, J., Petit, S., Gele, A. & Lepetit, M. B. An ab initio study of magneto-electric coupling of YMnO3. J. Phys. Condens. Matter 25, 496004 (2013).

28. Bramwell, S. T. Temperature dependence of the isotropic exchange constant.J. Phys. Condens. Matter 2, 75277536 (1990).29. Lancaster, T. et al. Magnetism in geometrically frustrated YMnO3 under hydrostatic pressure studied with muon spin relaxation. Phys. Rev. Lett. 98, 197203 (2007).

30. Kozlenko, D. P. et al. High-pressure effect on the crystal and magnetic structures of the frustrated antiferromagnet YMnO3. JETP Lett. 82, 212216 (2005).

31. Haiger, P. S. et al. Quantum and thermal ionic motion, oxygen isotope effect, and superexchange distribution in La2CuO4. Phys. Rev. B 89, 085113 (2014).

32. Chernyshev, A. L. & Brenig, W. Thermal conductivity in large-J two-dimensional antiferromagnets: role of phonon scattering. Phys. Rev. B 92, 054409 (2015).

33. Chernyshev, A. L. & Zhitomirsky, M. E. Spin waves in a triangular lattice antiferromagnet: decays, spectrum renormalization, and singularities. Phys. Rev. B 79, 144416 (2009).

34. Ji, S. et al. Spin-lattice order in frustrated ZnCr2O4. Phys. Rev. Lett. 103, 037201 (2009).

35. Schilfgaarde, M. V., Abrikosov, I. A. & Johansson, B. Origin of the Invar effect in iron-nikel alloys. Nature 400, 4649 (1999).

36. Tchernyshyov, O., Moessner, R. & Sondhi, S. L. Order by distortion and string modes in pyrochlore antiferromagnets. Phys. Rev. Lett. 88, 067203 (2002).

37. Jensen, J. & Houmann, J. G. Spin waves in terbium. II. Magnon-phonon interaction. Phys. Rev. B. 12, 320331 (1975).

38. Chandra Sekhar, M., Lee, S., Choi, G., Lee, C. & Park, J.-G. Doping effects of hexagonal manganites Er1 xYxMnO3 with triangular spin structure. Phys. Rev.

B 72, 014402 (2005).39. Ewings, R. A. et al. HORACE: software for the analysis of data from single crystal spectroscopy experiments at time-of-ight neutron instruments. Nucl. Instrum. Meth. A 834, 132142 (2016).

40. Parlinski, K., Li, Z. Q. & Kawazoe, Y. First-principles determination of the soft mode in cubic ZrO2. Phys. Rev. Lett. 78, 40634066 (1997).

41. Baroni, S., De Gironcoli, S., Dal Corso, A. & Giannozzi, P. Phonons and related crystal properties from density-functional perturbation theory. Rev. Mod. Phys. 73, 515562 (2001).

42. Kresse, G. & Furthmu+ller, J. Efcient iterative schemes for ab initio total-energy calculations using a plane-wave basis set. Phys. Rev. B 54, 1116911186 (1996).

43. White, R. M. et al. Diagonalization of the antiferromagnetic magnon-phonon interaction. Phys. Rev. 139, A450A454 (1965).

Acknowledgements

We thank Kisoo Park and W.J.L. Buyers for useful discussions. The work at theIBS CCES was supported by the research programme of Institute for Basic Science (IBS-R009-G1). Experiment at J-PARC was carried out through the proposal number (2011B0064). The work of ALC was supported by the U.S. DOE, Ofce of Science, Basic Energy Sciences under Award Number DE-FG02-04ER46174. The workat Rutgers University was supported by the DOE under Grant Number DOE: DE-FG02-07ER46382.

Authors contributions

J.G.P. conceived and supervised the project. H.S., Y.Y., H.E. and S.W.C. grew the single crystals and H.S. carried out bulk characterization. J.O., M.D.L., J.J., T.G.P., H.W., K.N., S.O.K., Z.Y. and J.G.P. were involved in the neutron scattering experiments. H.H.N. carried out rst-principles phonon calculation and J.O. made all the spin-wave calculations under the guidance of J.G.P. and A.L.C. J.G.P., J.O. and A.L.C. wrote the paper after discussion with all the authors.

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/

How to cite this article: Oh, J. et al. Spontaneous decays of magneto-elastic excitations in non-collinear antiferromagnet (Y,Lu)MnO3. Nat. Commun. 7, 13146 doi: 10.1038/ncomms13146 (2016).

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/

r The Author(s) 2016

6 NATURE COMMUNICATIONS | 7:13146 | DOI: 10.1038/ncomms13146 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications