,

Corresponding Author

Bernhard Keimer

Max Planck Institute for Solid State Research

References

1.
Harris, A.B.; Kumar, D.; Halperin, B.I.; Hohenberg, P.C.
Dynamics of an antiferromagnet at low temperatures – spin-wave damping and hydrodynamics
2.
Tyč, S.; Halperin, B.I.
Damping of spin waves in a two-dimensional Heisenberg antiferromagnet at low temperatures
3.
Kopietz, P.
Magnon damping in the two-dimensional quantum Heisenberg antiferromagnet at short wavelengths

In collaboration with:

D.A. Tennant, M.C.R. Gibson (Helmholtz-Zentrum Berlin)

S.D. Wilson, R.J. Birgeneau (Lawrence Berkeley National Laboratory, USA)

Department "Solid State Spectroscopy"

Lifetimes of antiferromagnetic magnons in two and three dimensions

Authors

Ph. Leininger, T. Keller, and B. Keimer

Departments

Solid State Spectroscopy (Bernhard Keimer)

A high-resolution neutron technique is used to measure magnon lifetimes in the prototypical 2D and 3D antiferromagnets Rb2MnF4 and MnF2, over the full Brillouin zone and a wide range of temperatures. We rederived theories of the lifetime from magnon-magnon scattering, for application beyond asymptotically small regions of q and T. Our computations, plus a small contribution from collisions with domain boundaries, yield excellent quantitative agreement with the data. Comprehensive understanding of magnon lifetimes in simple antiferromagnets provides a foundation for research on complex magnets.

The lifetimes of spin excitations in complex magnets, including quantum-disordered spin chains and ladders and noncollinear antiferromagnets, have recently been the subject of considerable attention. The relative contributions of magnon-magnon scattering and various forms of disorder to the magnon lifetime have been of particular interest, as the same mechanisms are also responsible for the decay of thermal transport currents in magnets. Remarkably, however, our understanding of magnon lifetimes in simple, collinear, large-spin antiferromagnets has remained incomplete for decades. To provide a firm foundation for current research on more complex magnetic systems, we have reexamined this issue both experimentally and theoretically.

Dynamics in Heisenberg antiferromagnets have been studied very extensively by means of theoretical approaches such as hydrodynamics, dynamical scaling, and interacting spin-wave methods. In the dilute magnon gas in two and three dimensions, the dominant damping processes are identified by diagrammatic perturbation theory to be elastic pair-wise collisions of magnons; the transition rate is given by Fermi’s Golden Rule. Analytic formulae for the relaxation as a function of wavevector q and temperature T were derived for this mechanism. Following the discovery of high-temperature superconductivity in doped two-dimensional (2D) antiferromagnets, recent theoretical work has focused specifically on the differences between the behavior of magnons in two and three dimensions. Surprisingly, in recent experimental studies of both 2D and 3D antiferromagnets, the measured linewidths disagreed considerably with predicted values, and anomalous additional line broadening at low temperatures that eluded explanation was observed. These findings provide additional motivation for a detailed investigation of the lifetimes of magnons in 2D and 3D antiferromagnets.

<strong>Fig. 1:</strong> Linewidth data over the full Brillouin zone in (a) Rb<sub>2</sub>MnF<sub>4</sub> and (b) MnF<sub>2</sub>. The curves show self-consistent calculations, with no free parameters, of 4-magnon scattering, added to a small, nearly <em>T</em>-independent contribution, peaked at low <em>q</em>, representing 2-magnon scattering from crystal blocks (with one free parameter, fitted to the 3K data). Zoom Image
Fig. 1: Linewidth data over the full Brillouin zone in (a) Rb2MnF4 and (b) MnF2. The curves show self-consistent calculations, with no free parameters, of 4-magnon scattering, added to a small, nearly T-independent contribution, peaked at low q, representing 2-magnon scattering from crystal blocks (with one free parameter, fitted to the 3K data). [less]

Rb2MnF4 and MnF2 are near-ideal realizations of 2D and 3D antiferromagnets, respectively. Due to their structural and chemical similarity, they are a well-suited pair of compounds for comparison of linewidths in 2D and 3D. The magnetic interactions are dominated by the antiferromagnetic exchange coupling J between S=5/2 Mn2+ spins; the dispersion relations in the two systems are very similar. We used the high-resolution NRSE-TAS (neutron resonance spin-echo triple-axis spectroscopy) technique to measure magnon linewidths in the ab-plane of these antiferromagnets (see Figs. 1– 3). Temperatures ranged from 3K up to 0.8TN (Rb 2MnF4) and 0.6TN (MnF 2); the wavevectors spanned the full magnetic Brillouin zone.

We initially compared the results with analytical approximations for the magnon linewidth originating from 4-magnon (two-in/two-out) collisions. The corresponding expressions derived by Harris et al. [1] and Tyc and Halperin [2] apply to restricted regions of small wavevectors at low temperatures. Their general low-temperature expressions, which contain no adjustable parameters and involve summations over the Brillouin zone, are intended for use at long wavelengths. Umklapp processes, which become increasingly important as q and T increase, are not treated properly in the scattering matrix elements.

We rederived the scattering matrix elements for general q, treating Umklapp processes correctly , and then computed the corresponding linewidths Γmag(q,T) numerically. It became clear that "on-shell" computations do not give correct line broadenings, particularly near the zone boundaries, because of the curvature of the dispersion relation and the large density of states present there. For this reason, we evaluated Γmag(q,T) self-consistently. We incorporated experimentally-measured values of the renormalized magnon energies into the calculation, and hence our results can be applied at higher temperatures than those considered by Harris et al. and Tyc and Halperin. Through such calculations, we determined that incorrect treatment of Umklapp processes yields inaccurate linewidth results, even at small wavevectors and low temperatures. Independent of how these are treated, our numerical results differ considerably in magnitude from the analytical approximations, and thus numerical evaluations are indispensable even for the limiting cases of small q and low T.

<strong>Fig. 2:</strong> Linewidths at 3K for Rb<sub>2</sub>MnF<sub>4</sub> (red) and MnF<sub>2</sub> (blue), with the calculated 4-magnon scattering subtracted off (this is less than 0.4&micro;eV for Rb<sub>2</sub>MnF<sub>4</sub> and 0.05&micro;eV for MnF<sub>2</sub>, for all <em>q</em>). The dashed curves show fits to 2-magnon scattering (see text); the fits shown by the solid curves also allow for a constant offset. Zoom Image
Fig. 2: Linewidths at 3K for Rb2MnF4 (red) and MnF2 (blue), with the calculated 4-magnon scattering subtracted off (this is less than 0.4µeV for Rb2MnF4 and 0.05µeV for MnF2, for all q). The dashed curves show fits to 2-magnon scattering (see text); the fits shown by the solid curves also allow for a constant offset. [less]

An additional potential scattering process consists of a non-momentum-conserving collision of a magnon with a boundary (e.g. a grain boundary, antiferromagnetic domain wall, or microcrack in the crystal). A rough estimate of the associated magnon lifetime for such a process is τb=L/(2v(q)), where L is the size of the grain, domain, or microcrystallite and v(q) is the magnon velocity. The corresponding additive contribution Γb=hτb-1 to the overall linewidth exhibits a peak at low q that is similar to that in the 3K data (see Fig. 2). We fitted these data to the corresponding expression after first subtracting off the calculated Γmag. Using the measured magnon dispersion data, we obtained an average L of 0.58±0.04µm for Rb2MnF4 and 0.50±0.04µm for MnF2. The latter is consistent with the antiferromagnetic domain size of 0.63±0.05µm obtained from independent high-resolution neutron Larmor diffraction measurements performed on the same MnF2 crystal at 5K.

The total calculated linewidth Γ=Γmagb agrees well with the data (Figs. 1–3). In Rb2MnF4, the agreement worsens with increasing T at intermediate q (0.025–0.15r.l.u.); in MnF2, for all q, and the discrepancy is largest at intermediate q (0.2–0.4r.l.u.). Additional scattering channels, such as 6-magnon (three-in/three-out) scattering, may become important at higher temperatures.

<strong>Fig. 3:</strong> (a) Linewidth data vs. <em>T</em> in Rb<sub>2</sub>MnF<sub>4</sub> (symbols), measured by NRSE-TAS at wavevectors along the (<em>qq</em>0) direction. The curves are calculated as described in Fig. 1. (b) Temperature-dependent component of the magnon linewidths in Rb<sub>2</sub>MnF<sub>4</sub> (open symbols) and MnF<sub>2</sub> (filled symbols). The fitted contribution from 2-magnon scattering from crystal blocks at 3K (dashed curves from Fig. 2) has been subtracted off. The lines (solid for Rb<sub>2</sub>MnF<sub>4</sub> and dashed for MnF<sub>2</sub>) are power-law fits to the data. The data and curves are color-coded by <em>q/q<sub>ZB</sub></em>. Inset: the exponent corresponding to such power-law fits in Rb<sub>2</sub>MnF<sub>4</sub> (red symbols) and MnF<sub>2</sub> (blue symbols). Zoom Image
Fig. 3: (a) Linewidth data vs. T in Rb2MnF4 (symbols), measured by NRSE-TAS at wavevectors along the (qq0) direction. The curves are calculated as described in Fig. 1. (b) Temperature-dependent component of the magnon linewidths in Rb2MnF4 (open symbols) and MnF2 (filled symbols). The fitted contribution from 2-magnon scattering from crystal blocks at 3K (dashed curves from Fig. 2) has been subtracted off. The lines (solid for Rb2MnF4 and dashed for MnF2) are power-law fits to the data. The data and curves are color-coded by q/qZB. Inset: the exponent corresponding to such power-law fits in Rb2MnF4 (red symbols) and MnF2 (blue symbols). [less]

Selected linewidth data from which the fitted curves to the 3K data (Fig. 2) have been subtracted are compared in the main panel of Fig. 3b. The exponents from (phenomenological) power-law fits for these and additional values of q/qZB , where qZB is the zone-boundary wavevector, are shown in the inset. The q-dependence of the exponents is non-monotonic; the values for the two compounds agree within error at low q, and diverge near the respective zone boundaries. The exponent for Rb2MnF4 there is 2.6±0.1, which is not too different from Kopietz’ analytical result of 3 for the isotropic case (derived from approximations made to simplify the general expression of Tyc and Halperin at the zone boundary) [3], though the region of validity of his approximations lies lower in temperature than our data.

In summary, comprehensive comparison of experimental data and numerical results suggests that the processes that determine the magnon lifetime in 2D and 3D antiferromagnets are two-in/two-out magnon collisions and boundary scattering. The difference between 2D and 3D systems is thus quantitative, not qualitative: it originates from the difference between near-cylindrical and ellipsoidal (near-spherical) surfaces of constant magnon energy with respect to the conservation of energy and momentum in 4-magnon collisions. The seminal theories of 4-magnon scattering considered here could not previously be compared realistically and quantitatively with experimental data. This limitation is relieved here through proper inclusion of the Umklapp scattering and the use of numerical evaluations, permitting the first complete experimental and theoretical description of magnon linewidths over broad ranges of wavevector and temperature in two and three dimensions. The quantitative understanding of prototypical antiferromagnets that we have obtained provides a solid basis for investigations of the influence of noncollinearity, disorder, and doping in more complex magnetic systems. Comprehensive, momentum-resolved linewidth data that are described accurately by theory also enable the first calculation of the magnon-mediated thermal conductivity without adjustable parameters.


 
loading content