Nature Communications

Nat Commun. 8: 1152

Breakdown of magnons in a strongly spin-orbital coupled magnet

http://orcid.org/0000-0002-0025-1446Stephen M. Winter

winter@physik.uni-frankfurt.de

¹, Kira Riedl¹, Pavel A. Maksimov², Alexander L. Chernyshev², http://orcid.org/0000-0001-6383-3200Andreas Honecker³, Roser Valentí¹

1. 0000 0004 1936 9721grid.7839.5Institut für Theoretische Physik, Goethe-Universität Frankfurt, Max-von-Laue-Strasse 1, 60438 Frankfurt am Main, Germany

2. 0000 0001 0668 7243grid.266093.8Department of Physics and Astronomy, University of California, Irvine, CA 92697 USA

3. 0000 0001 2290 0120grid.7901.fLaboratoire de Physique Théorique et Modélisation, CNRS UMR 8089, Université de Cergy-Pontoise, 95302 Cergy-Pontoise Cedex, France

© The Author(s) 2017

DOI: 10.1038/s41467-017-01177-0

Published online: 27 October 2017

Abstract

The description of quantized collective excitations stands as a landmark in the quantum theory of condensed matter. A prominent example occurs in conventional magnets, which support bosonic magnons—quantized harmonic fluctuations of the ordered spins. In striking contrast is the recent discovery that strongly spin-orbital-coupled magnets, such as α-RuCl₃, may display a broad excitation continuum inconsistent with conventional magnons. Due to incomplete knowledge of the underlying interactions unraveling the nature of this continuum remains challenging. The most discussed explanation refers to a coherent continuum of fractional excitations analogous to the celebrated Kitaev spin liquid. Here, we present a more general scenario. We propose that the observed continuum represents incoherent excitations originating from strong magnetic anharmonicity that naturally occurs in such materials. This scenario fully explains the observed inelastic magnetic response of α-RuCl₃ and reveals the presence of nontrivial excitations in such materials extending well beyond the Kitaev state.

Abstract

The observation of magnetic excitation continua in putative Kitaev spin liquid compounds is consistent with predictions of exotic quantum behavior, attracting significant interest. Here, the authors show the observed continua could instead be a consequence of realistic non-Kitaev interactions.

Introduction

From magnons in ordered magnets to phonons in periodic crystals, the appearance of bosonic collective excitations is ubiquitous in condensed phases of matter¹. For this reason, special attention is given to those states that support more exotic collective modes, for which the conventional paradigm breaks down. In the context of magnetic phases, the breakdown of magnons is commonly thought to require closeness to an unconventional state such as a quantum spin liquid²–⁴. A notable example occurs in Kitaev’s exactly solvable honeycomb model⁵, for which strongly anisotropic and bond-dependent interactions fractionalize conventional spin excitations into Majorana spinons and fluxes. This Kitaev state has recently risen to prominence due to the suggestion that it may be realized in heavy metal 4d ⁵ and 5d ⁵ insulators via a specific interplay between the crystal field and strong spin-orbit coupling⁶, and, consequently, a variety of candidate materials based on Ir⁴⁺ and Ru³⁺ have been intensively explored⁷. Encouragingly, evidence of a continuum of magnetic excitations that is inconsistent with conventional magnons was found in the majority of such materials, including the two-dimensional (2D) honeycomb Na₂IrO₃ ⁸, ⁹ and α-RuCl₃ ¹⁰–¹⁴, as well as the three-dimensional (3D) analogs β-,γ-Li₂IrO₃ ¹⁵, despite all of them having magnetically ordered ground states.

While the observed excitation continua in these systems have been interpreted in terms of signatures of the Kitaev state, the low-symmetry crystalline environment of the real materials also allows various additional interactions beyond Kitaev’s model¹⁶–¹⁸, which are thought to be large based on both experimental¹⁹, ²⁰ and theoretical¹⁸, ²¹, ²² considerations. In this sense, understanding the mechanism for the breakdown of magnons and the appearance of a broad continuum of magnetic excitations remain a key challenge.

In this work, we study a representative case α-RuCl₃, which forms a layered 2D honeycomb lattice and displays zigzag magnetic order below T N ~7 K¹², ¹³, ²³. We specifically address the recent inelastic neutron scattering (INS) measurements, which have revealed low-energy magnons²⁴ coexisting with an intense excitation continuum¹². The latter continuum possesses a distinctive six-fold star shape in momentum space, and large intensity at the 2D Γ-point over a wide energy range E = 2–15 meV¹². To resolve the nature of this continuum, we take two complementary approaches. We first theoretically investigate the neutron spectra over a range of relevant magnetic interactions in order to determine the correct spin Hamiltonian for α-RuCl₃, which has been a subject of intense recent discussion¹⁸, ²⁵–²⁸. Second, we identify the conditions that lead to the breakdown of conventional magnons in the presence of strongly anisotropic and frustrated interactions, revealing that nontrivial excitations naturally persist well beyond the Kitaev spin liquid.

Results

The model

Based on previous ab initio studies¹⁸, ²⁵–²⁸, the largest terms in the spin Hamiltonian of α-RuCl₃ are generally expected to include nearest neighbor Heisenberg J ₁, Kitaev K ₁, and off-diagonal Γ ₁ couplings, supplemented by a third neighbor Heisenberg J ₃ term:

1. H=∑i,jJ1Si⋅Sj+K1SiγSjγ+Γ1SiαSjβ+SiβSjα+∑i,jJ3Si⋅Sj,[Image: 41467_2017_1177_Article_Equ1.gif]

where ⟨i, j⟩ and ⟨⟨⟨i, j⟩⟩⟩ refer to summation over first and third neighbor bonds, respectively (see Fig. 1). The bond-dependent variables {α, β, γ} distinguish the three types of first neighbor bonds, with {α, β, γ} = {y, z, x}, {z, x, y}, and {x, y, z} for the X-bonds, Y-bonds, and Z-bonds, respectively. The third neighbor interactions are bond-independent. The phase diagram of this model has been discussed previously¹⁷, ¹⁸, ²⁶, ²⁹, and is further detailed in Supplementary Note 1; here we review the key aspects.

Fig. 1View larger versionFig. 1. From material to model. Within the honeycomb ab-layer of α-RuCl₃ are illustrated the RuCl₆ octahedra, magnetic zigzag ordering pattern, and definition of the underlying magnetic interactions. Crystal axes are labeled with respect to the C2/m structure

In the limit J ₁ = Γ ₁ = J ₃ = 0, the ground state is a gapless Z ₂ spin liquid for either positive or negative K ₁, as demonstrated in Kitaev’s seminal work⁵. Small perturbations from the pure K ₁ limit may induce various magnetically ordered states, such as the zigzag antiferromagnetic (AFM) state observed in α-RuCl₃ and shown in Fig. 1. The simplest perturbation is the introduction of a finite J ₁, which yields the well-studied (J ₁, K ₁) nearest neighbor Heisenberg-Kitaev (nnHK) model. This model hosts zigzag order in the region K ₁ > 0, J ₁ < 0, as discussed in Supplementary Note 1. Accordingly, previous analysis of the powder INS experiments within the context of the nnHK model¹³ suggested that K ₁ ~+7 meV, and ∣J₁∕K₁∣~0.3 - 0.7 for α-RuCl₃. On this basis, the excitation continua observed experimentally were initially interpreted in terms of proximity to the AFM K ₁ > 0 spin liquid¹², ¹³. However, the further consideration of finite Γ ₁ and J ₃ interactions in Eq. (1) significantly expands the experimentally relevant region, as both interactions generally stabilize zigzag order. Indeed, recent ab initio studies¹⁸, ²⁵–²⁸ have suggested that the zigzag order in α-RuCl₃ emerges from J₁~0,  K₁ < 0,  Γ₁ > 0, and J ₃ > 0, with ∣Γ₁∕K₁∣~0.5 - 1.0 and ∣J₃∕K₁∣~0.1 - 0.5, as reviewed in Supplementary Note 2. That is, K ₁ is ferromagnetic, and supplemented by significant Γ ₁ and J ₃ interactions. Such interactions would represent large deviations from both Kitaev’s original model and the region identified by initial experimental analysis. Before discussing the origin of the excitation continua, it is therefore crucial to first pinpoint the relevant interaction parameters.

In order to address this issue directly, we have computed the neutron scattering intensity ℐ(k, ω) for a variety of interactions within the zigzag ordered phase via both linear spin-wave theory (LSWT) and exact diagonalization (ED). For the latter case, we combine results from various periodic 20-site and 24-site clusters compatible with the zigzag state in order to probe a wider range of k-points (see “Methods” section). Full results for the complete range of models are presented in Supplementary Note 5. Here, we highlight the key results for two representative sets of interactions. Within the (J ₁, K ₁) nnHK model, we focus on Model 1 (J ₁ = −2.2, K ₁ = +7.4 meV; ∣J₁∕K₁∣ = 0.3), which lies on the border of the region initially identified in ref. ¹², close to the spin liquid. Beyond the nnHK model, we consider Model 2 (J ₁ = −0.5, K ₁ = −5.0, Γ ₁ = +2.5, J ₃ = +0.5 meV) for which parameters have been guided by recent ab initio studies¹⁸, ²⁵–²⁸, and further optimized to improve agreement with the experimental spectra. Results for Models 1 and 2 are first presented in Figs. 2 and 3, which show detailed ω-dependence and k-dependence of ℐ(k, ω), along with the evolution of the spectra upon changing parameters toward the K ₁ > 0 or K ₁ < 0 spin liquid regions.

Fig. 2View larger versionFig. 2. Neutron scattering intensity ℐ(k, ω) within the nnHK model. a–c Detailed results for Model 1 (J ₁ = −2.2, K ₁ = +7.4 meV): a ℐ(k, ω) computed via LSWT; results are averaged over the three zigzag ordering wavectors, parallel to the X-bonds, Y-bonds, and Z-bonds. Inset: Definition of Brillouin zone and high-symmetry k-points. b ED results, combining data from several 20-site and 24-site periodic clusters (see “Methods”). c ED k-dependence of ℐ(k, ω) integrated over the indicated energies, as obtained from a single 24-site cluster respecting all symmetries of Eq. (1) (see “Methods”). d Comparison of Γ-point intensities for the K ₁ = +7.7 meV AFM spin liquid (exact results⁵³, ⁵⁴), Model 1 (ED), and the experimental data for α-RuCl₃ ¹². e Evolution of the ED Γ-point intensity with decreasing ∣J₁∕K₁∣, showing absence of low-energy intensity close to the K ₁ > 0 spin liquid. The top three interaction sets correspond to zigzag order, while the bottom is the K ₁ > 0 spin liquid. For all spectra, a Gaussian broadening of 0.5 meV has been applied

Fig. 3View larger versionFig. 3. Neutron scattering intensity ℐ(k, ω) within the extended model. a–c Detailed results for Model 2 (J ₁ = −0.5, K ₁ = −5.0, Γ ₁ = +2.5, J ₃ = +0.5 meV): a ℐ(k, ω) computed via LSWT; results are averaged over the three zigzag domains with ordering wavectors parallel to the X-bonds, Y-bonds, and Z-bonds. b ED results, combining data from several 20-site and 24-site periodic clusters (see “Methods”). c ED k-dependence of ℐ(k, ω) integrated over the indicated energies, as obtained from a single 24-site cluster respecting all symmetries of Eq. (1) (see “Methods”). d Comparison of Γ-point intensities for the K ₁ = −7.7 meV FM spin liquid (exact results⁵³, ⁵⁴), Model 2 (ED), and the experimental data for α-RuCl₃ ¹². e Evolution of the ED Γ-point intensity with decreasing ∣Γ₁∕K₁∣, showing significant broadening at finite Γ ₁. The top three interaction sets correspond to zigzag order, while the bottom is the K ₁ < 0 spin liquid. For all spectra, a Gaussian broadening of 0.5 meV has been applied

Nearest neighbor Heisenberg-Kitaev model

We begin by analyzing the spectra ℐ(k, ω) within the zigzag phase of the (J ₁, K ₁) nnHK model, starting with Model 1 (Fig. 2). Despite proximity to the spin liquid, the ED calculations on Model 1 (Fig. 2b) show sharp dispersive modes appearing over the majority of the Brillouin zone that are consistent with the conventional magnons of LSWT (Fig. 2a). Indeed, for energies below the main spin-wave branch (ω = 1.3–2.3 meV), intensity is localized around the M-points and Y-points, corresponding to the pseudo-Goldstone modes of the zigzag order (Fig. 2c). ED calculations show clear spin-wave cones emerging from such points and extending to higher energies. Large deviations from LSWT are observed only for the highest energy excitations, which appear near the 2D Γ-point for energies ω > 5 meV. Here, the ED calculations display a broad continuum-like feature centered at ω~K₁ that resembles the response expected for the K ₁ > 0 Kitaev spin liquid, as shown in Fig. 2d. However, comparison with the experimental ℐ(Γ, ω) shows poor agreement, while the experimental intensity extends from 2 to 15 meV, the ED results for Model 1 predict intensity only at high energies >5 meV. Indeed, the evolution of the Γ-point intensity with ∣J₁∕K₁∣ is shown in Fig. 2e. On approaching the K ₁ > 0 spin liquid by decreasing ∣J₁∕K₁∣, excitations at the Γ-point shift to higher energy, such that none of the parameters in the vicinity of the spin liquid reproduce the experimental ω-dependence of ℐ(Γ, ω). Similar conclusions can also be drawn from recent Density Matrix Renormalization Group (DMRG) studies of the nnHK model³⁰. We therefore conclude that the broad features observed experimentally in ℐ(Γ, ω) at relatively low energies¹² are incompatible with the nnHK model with J ₁ < 0 and K ₁ > 0.

Extended ab initio guided model

In order to treat the effect of interactions beyond the nnHK model, we consider now the ab initio guided Model 2. In contrast to Model 1, ED calculations on Model 2 (Fig. 3b) show large deviations from standard LSWT (Fig. 3a) over a wide range of k and ω. This model reproduces many of the experimental spectral features¹², ²⁴. In particular, sharp single-magnon-like peaks appear only near the pseudo-Goldstone modes at the M-points and Y-points. Elsewhere in the Brillouin zone, broad continuum-like features are observed within the ED resolution. As demonstrated in Fig. 3c, we find significant intensity at low energies (ω < 2.3 meV), at both the Γ-points and (M,Y)-points. For the intermediate energy region (ω = 5.5–8.5 meV), ℐ(k) resembles the six-fold star shape observed in ref. ¹². At higher energies (ω > 10.5 meV) scattering intensity is mainly located at the Γ-point, also in accord with experiment. Furthermore, the ED results for the Γ-point intensity ℐ(Γ, ω) show a broad range of excitations peaked around 4 and 6 meV, and extending up to ~15 meV (Fig. 3d). Therefore, ED calculations on Model 2 reproduce all of the main experimental spectral features, validating the range of interactions indicated by ab initio calculations. The only aspect that is not quantitatively reproduced within the Model 2 is the magnitude of the gap at the M-point (~0.8 meV at the level of LSWT vs. ~2 meV experimentally¹³, ²⁴). This discrepancy may result from deviations from C ₃ symmetry, which are allowed within the C2/m space group¹⁸, ³¹, but not considered here for simplicity (see Supplementary Fig. 11). Interestingly, the spectral features at the Γ-point become dramatically sharper on approaching the K ₁ < 0 spin liquid, as shown in the evolution of ℐ(Γ, ω) with the ratio ∣Γ₁∕K₁∣ (Fig. 3e). This result reveals that the broad continuum may not be directly associated with a proximity to the Kitaev state.

Magnon stability beyond LSWT

To gain further insight into the reason for such a drastic contrast between the stability of magnons in Models 1 and 2, it is useful to consider possible magnon decay channels in the zigzag ordered phase. At the level of LSWT, the spin-wave Hamiltonian is truncated at quadratic order, and can be written H2= ∑k,mϵk,mak,m†ak,m in terms of magnon creation (annihilation) operators a † (a), where ϵk,m denotes the dispersion of the mth magnon band. In this harmonic approximation, magnons represent sharp, well-defined excitations. However, when higher-order anharmonic terms are included, the total magnon number Ntot= ∑k,mak,m†ak,m is typically not a conserved quantity, such that the stability of magnons is not guaranteed beyond quadratic order. Quantum fluctuations associated with the higher-order anharmonic decay terms may mix sharp single-magnon modes with the multi-magnon continuum³²–³⁴. Similar considerations also apply to the breakdown of other collective modes, such as phonons in anharmonic crystals³⁵, ³⁶. From this perspective, a large decay rate is expected for any single-magnon mode that is energetically degenerate with the multi-particle continuum, unless there are specific symmetries guaranteeing that the two do not couple. It is therefore useful to consider the prerequisites for magnon breakdown in the presence of the strongly anisotropic interactions of Eq. (1).

Magnon decay channels for the nnHK model

We first examine the stability of magnons in the nnHK model. For pure J ₁ and K ₁ interactions, the total spin projections Stotγ= ∑iSiγ are conserved along the cubic axes γ = {x, y, z} modulo two. Since the ordered moment also lies along one of the cubic axes in the zigzag phase²⁰, ³⁷ (see Fig. 4c), the possible magnon decay channels are restricted. In the local picture, the relevant quantum fluctuations are local singlet SixSjx↑↓=↓↑ and triplet SixSjx↑↑=↓↓ fluctuations shown in Fig. 4a, with ΔStotz=0 and 2, respectively. In the magnon picture, the Hamiltonian can only contain even-order terms (i.e., ℋ = ℋ₂ + ℋ₄ + ...), analogous to conventional Heisenberg antiferromagnets with collinear ordered spins³², ³⁴. For example, the fourth-order decay process due to ℋ₄ mixes the one-magnon states with the three-magnon continuum (ΔN tot = ±2), where

2. H4= ∑1-4V1234a1†a2†a3†a4†δk1+k2+k3-k4+H.c.[Image: 41467_2017_1177_Article_Equ2.gif]

Here, the bold index (n ≡ k n, m n) labels both momentum and magnon band. This process is pictured in Fig. 4b. As noted above, the effect of such terms depends crucially on the availability of low-energy three-magnon states in which to decay.

Fig. 4View larger versionFig. 4. Magnon decay channels for the nnHK model. a Local picture of quantum fluctuations away from zigzag order. The energy cost for the left process vanishes on approaching the spin liquid ∣J₁∕K₁∣→0. b Momentum space picture for the corresponding fourth-order decay process due to ℋ₄. c Ordered moment direction for Model 1 (J ₁ = −2.2, K ₁ = +7.4 meV), corresponding to the zigzag domain with ordering wavevector Q = Y. d LSWT dispersions ϵk,m, and three-magnon density of states (DOS) for Model 1 for the same zigzag domain as c. The dashed line indicates the bottom of the three-magnon continuum E3min, which is coincident with the lowest magnon band

The density of three-magnon states for Model 1 is shown in Fig. 4d, based on the one-magnon dispersions obtained in LSWT. At each k-point, the lowest energy three-magnon state aq1†aq2†aq3†0 (with q ₁ + q ₂ + q ₃ = k) is obtained by placing two particles in the pseudo-Goldstone modes at opposite M-points (q ₁ + q ₂ = 0), and the third particle at q ₃ = k, with total energy E3min(k)=ϵk,1+2ϵM,1. This implies E3min(k)≥ϵk,1. That is, the three-magnon states lie above the lowest one-magnon band at every k-point. As a result, every magnon in the lowest band remains kinetically stable, due to the absence of low-energy three-particle states in which to decay. Precisely this condition ensures the stability of low-energy magnons in conventional isotropic antiferromagnets, and explains the sharp magnon-like peaks observed in Fig. 2b for Model 1. Strong spectral broadening in the nnHK model can occur only for high-lying excitations with ϵk,m>ϵk,1=E3min, where the density of three-magnon states is finite, such as at the 2D Γ-point. On approaching the spin liquid (at J ₁/K ₁ = 0), this condition is relaxed due to the vanishing dispersion of the lowest magnon band (i.e., ϵk,1→0), which corresponds to a vanishing energy cost the singlet fluctuations shown on the left of Fig. 4a. The relevant fluctuations in the limit J ₁/K ₁→0 therefore correspond to ΔN tot = ±2. For other values of J ₁/K ₁, the majority of magnons are expected to remain stable due to the absence of low-energy three-magnon states.

Magnon decay channels for the extended model

In Model 2, the character of the quantum fluctuations away from zigzag order is notably different (Fig. 5). The finite Γ ₁ interaction reduces the local symmetry and leads to rotation of the ordered moments away from the cubic axes²⁰, ³⁷ (Fig. 5c). In the local picture, this allows additional single-spin fluctuations SixSjz↑↑=↓↑ (Fig. 5a), which correspond to odd-order anharmonic terms ℋ₃, ℋ₅, ... in the magnon Hamiltonian, where³³, ³⁴:

3. H3= ∑1-3Λ123a1†a2†a3δk1+k2-k3+H.c.[Image: 41467_2017_1177_Article_Equ3.gif]

At lowest order, such terms mix the single-magnon states with the two-magnon continuum (ΔN tot = ±1), via the scattering process depicted in Fig. 5b. The density of two-magnon states is shown in Fig. 5d, for the zigzag domain with Q = Y. In this case, at each k-point the lowest energy two-magnon state aq1†aq2†0 is obtained by placing one particle in the pseudo-Goldstone mode at an M-point, and the second particle at q ₂ = k − M, with total energy E2min(k)=ϵk-M+ϵM≠E3min. It should be emphasized that this condition differs from that of a conventional Heisenberg antiferromagnet (for which E2min=E3min)³⁴. In the case of Model 2, the difference is directly related to the strong anisotropic K ₁ and Γ ₁ interactions, which shift the pseudo-Goldstone modes to the M-points, such that only high-energy magnons remain at the Γ-point or ordering wavevector Q ³⁸. This shift therefore leads to an offset of the low-energy even and odd magnon states in k-space such that E2min(k)<ϵk,1 over a wide region of the Brillouin zone; there are many two-magnon states with equal or lower energy than the one-magnon states. Provided there is a finite Γ ₁, the spontaneous decay of single magnons into the two-particle continuum is therefore allowed even for the lowest magnon band. The decay rate is expected to be particularly large near the zone center, which represents a minimum in E2min. Similar kinematic conditions may also occur in other systems³⁴, ³⁹. For Model 2, the pseudo-Goldstone magnons near the M-points remain coherent due to the absence of low-energy two particle states in which to decay (Fig. 5d). This explains the experimental observation of sharp magnon-like modes near the M-points²⁴. In contrast, the magnon bands in the remainder of the Brillouin zone directly overlap with the two-particle continuum. It is therefore natural to anticipate a large decay rate even for the lowest magnon bands.

Fig. 5View larger versionFig. 5. Magnon decay channels for the extended model. a Local picture of additional quantum fluctuations away from zigzag order induced by Γ ₁ interactions. b Momentum space picture of the third-order decay process ℋ₃. c Ordered moment direction for Model 2 (J ₁ = −0.5, K ₁ = −5.0, Γ ₁ = +2.5, J ₃ = +0.5 meV) with zigzag ordering wavevector Q = Y, parallel to the Z-bond. d LSWT dispersions ϵk,m, and two-magnon DOS for Model 2 with Q = Y. Dashed lines indicate the bottom of the two-magnon and three-magnon continuum (E2min(k) and E3min(k), respectively)

To confirm this idea, we have computed the three-magnon interactions and decay rates for all magnon bands for Model 2 using the self-consistent imaginary Dyson equation (iDE) approach⁴⁰. Within this approach, it is assumed that the real part of the magnon self-energy is already captured by the LSWT parameters, while the imaginary part is obtained self-consistently (see “Methods” and Supplementary Note 3). The iDE approach therefore represents an extension of LSWT, in which the one-magnon excitations are broadened according to the momentum and band-dependent decay rate γ k,n, while other contributions to the neutron intensity from multi-magnon excitations are also absent⁴¹. As a result, comparison of LSWT, ED, and iDE results (Fig. 6) allows for the identification of the origin of different contributions to the spectra.

Fig. 6View larger versionFig. 6. Effects of two-magnon decays in ℐ(k, ω) for extended model. Results are shown for Model 2 computed via a LSWT, b self-consistent iDE approach, and c ED. Results in a and b are averaged over the different zigzag domains. The white and pink dashed lines indicate the bottom of the two-magnon continuum, E2min(k) for the different zigzag domains. In the iDE results, the effects of two-magnon decays strongly broadens any magnon bands overlapping with the two-magnon continuum

The predicted neutron scattering intensity within the iDE approach (Fig. 6b) captures many of the most notable features that are observed in the ED and experimental data, showing a significant improvement over the LSWT results (Fig. 6a). First, there is an almost complete washout of the two high-energy one-magnon modes due to strong decays. This implies that the higher-energy features >4 meV appearing in ED are primarily multi-magnon in character (including the 6 meV peak at the Γ-point). The appearance of these higher-energy features in the inelastic neutron response may arise partly from direct contributions from the broadened two-magnon continuum via the longitudinal component of the structure factor, which is not included in the iDE approach (see Supplementary Note 3). Second, the broadening of the two lower magnon bands in the iDE results and the resultant variation of their intensities are in a close agreement with the ED—particularly in a wide region near the Γ-point (see also Supplementary Fig. 5). These are precisely the features with which the LSWT results were most incompatible. Over much of the Brillouin zone—and especially for the higher magnon bands—the computed γ k,n is on the same scale as the one-magnon bandwidth, confirming the absence of coherent magnons.

Discussion

The general requirements for strong two-magnon decays are less restrictive than a proximity to a spin liquid state. Indeed, a large decay rate is ensured by the following three conditions: large anisotropic interactions, deviation of the ordered moments away from the high-symmetry axes, and strong overlap of the one-magnon states with the multi-magnon continuum (see Supplementary Note 3). Of these, the first two conditions ensure that the scattering vertex Λ123 is large—of the order of the underlying interactions, i.e., Λ123~OK1,Γ1. For α-RuCl₃, the strong overlap with the multi-magnon continuum is ensured by shifting of the low-energy magnons away from the Γ-point. Since the bottom of the two-magnon continuum must always have an energetic minimum at the Γ-point, the shifting of the pseudo-Goldstone modes to a finite momentum ensures the remaining higher-energy magnons are degenerate with the continuum near the zone center. Experimentally, these conditions are also likely to be satisfied by the zigzag ordered Na₂IrO₃ ⁹, and spiral magnets α-Li₂IrO₃, β-Li₂IrO₃, and γ-Li₂IrO₃ ⁴²–⁴⁴. This picture is also consistent with recent indications that the magnetically disordered phase observed at high pressure in β-Li₂IrO₃ ⁴⁵ is driven primarily by large Γ ₁ interactions⁴⁶.

With this in mind, there are two general scenarios that can explain the observed continuum excitations in α-RuCl₃ and the iridates A ₂IrO₃. In the first scenario, which has been advanced by several studies, the excitations can be treated as free particles with a small number of flavors. Such excitations are weakly interacting and have well-defined dispersions, but possess quantum numbers (e.g., ΔS tot = ±1/2) or topological properties inconsistent with the experimental neutron scattering selection rules (i.e., ΔS tot = 0, ±1). The appearance of the broad continuum in energy therefore results only from the fact that these fractional excitations must be created in multiples. If they could have been created individually, they would have represented long-lived and coherent quasiparticles with sharply peaked energies. This scenario indeed describes the Kitaev spin liquid, where the special symmetries of the Hamiltonian allow an exact description in terms of two flavors of particles: non-interacting Majorana spinons and localized fluxes⁵. Such excitations are long-lived, but belong to nontrivial topological sectors, and therefore cannot be created individually by any local operations. For the Kitaev spin liquid, the predicted continuum therefore represents coherent multiparticle excitations.

In contrast, upon moving away from the pure Kitaev point, the relevant symmetries that protect the spinons and fluxes are lifted both by additional magnetic interactions and by spontaneous symmetry breaking of the magnetic order. This tends to confine spinons into gauge neutral objects such as magnons⁴⁷, ⁴⁸. Despite this latter tendency, we have argued that coherent magnons are unlikely to appear at large Γ ₁ due to the strong anharmonicity in the magnon Hamiltonian. While this leaves open the possibility that nearly free Majorana spinons persist into the zigzag ordered phase in some regions of the Brillouin zone, a more general scenario is that the observed continua represent fully incoherent excitations. In this second scenario, the excitations are not describable in terms of any type of free particles with small decay rates and well-defined dispersions. The broad continua therefore reflect the absence of coherent quasiparticles altogether, rather than particular experimental selection rules related to fractionalization. At present, it is not clear which of these scenarios applies to the iridates and α-RuCl₃, although a key role must be played by both the Kitaev K ₁ and off-diagonal couplings such as Γ ₁. In any case, the study of these materials calls into question the stability of magnetic quasiparticles in the presence of strongly anisotropic interactions.

In summary, we have shown that all main features of the magnetic excitations in α-RuCl₃ ¹², ¹³, ²⁴ are consistent with strongly anisotropic interactions having signs and relative magnitudes in agreement with ab initio predictions. The ferromagnetic Kitaev coupling (K ₁ < 0) is supplemented by a significant off-diagonal term (Γ ₁ > 0) that plays a crucial role in establishing both the zigzag order and the observed continua. In the presence of such interactions, the conventional magnon description breaks down even deep in the ordered phase, due to strong coupling of the one-magnon and two-magnon states. This effect is expected to persist over a large range of the phase diagram suggesting that the observed continua in α-RuCl₃ and the iridates A ₂IrO₃ represent a rich and general phenomenon extending beyond the Kitaev spin liquid. For this class of strongly spin-orbital-coupled magnets, the presence of complex and frustrated anisotropic interactions leads naturally to dominant anharmonic effects in the inelastic magnetic response. Fully describing the dynamics of these and similar materials therefore represents a formidable challenge that is likely to reveal aspects not found in conventional isotropic magnets.

Methods

Exact diagonalization

The neutron scattering intensity was computed via:

4. I(k,ω)∝f2(k)∫dt∑μ,νδμ,ν-kμkν∕k2×∑i,jSiμ(t)Sjν(0)e-ik⋅ri-rj-iωt,[Image: 41467_2017_1177_Article_Equ4.gif]

where f(k) is the atomic form factor of Ru³⁺ from ref. ⁴⁹. ED calculations were performed using the Lanczos algorithm⁵⁰, on several 20-site and 24-site clusters with periodic boundary conditions. Such periodic clusters are detailed in Supplementary Note 4. Excitations were computed using the continued fraction method⁵¹. Further details and additional results are presented in the Supplementary Notes 4 and 5; these extensive calculations go beyond previous ED studies¹⁶, ¹⁷, ²⁰, ²⁶, ²⁹, which focused mainly on the static properties, or a limited portion of the phase diagram. ED results shown for the high-symmetry Γ, M, Y, X, and Γ′ points were averaged over all clusters. The ED k-dependence of ℐ(k, ω), integrated over the energy windows E = 1.3–2.3, 5.5–8.5, and 10.5+ meV (Figs. 2c, 3c), was obtained from a single 24-site cluster respecting all symmetries of the model. The discrete ED spectra were Gaussian broadened by 0.5 meV, consistent with the width of experimental features¹². The intensities were also averaged over the same range of out-of-plane momentum as in the experiment¹².

Linear spin-wave theory

LSWT results shown in Figs. 1 and 2 were obtained with the aid of SpinW⁵². Following the approach with the ED data, the discrete LSWT spectra were as well Gaussian broadened by 0.5 meV and the intensities were also averaged over the same range of out-of-plane momentum as with ED and in the experiment¹².

Imaginary self-consistent Dyson equation approach

In order to calculate magnon decay rates γ k,n, we have evaluated three-magnon interaction vertices by performing rotation to local reference frames of spins. The obtained value of the real–space interaction is quite large, about ~3 meV. Next, the Born approximation calculation of the decay rates results in unphysical divergencies³⁴, thus the self-energy Σk,n needs to be regularized. We have used the so-called iDE approach: a self-consistent solution on the imaginary part of the Dyson’s equation, Σk,n(ϵk,n + iγk,n) =  - iγk,n, see ref. ⁴⁰. We have obtained the regularized broadening for the magnon spectrum and have calculated the transverse part of the dynamical structure factor, shown in Fig. 6, by adding the calculated decay rates to experimental resolution of 0.25 meV. The spectral function is approximated as a Lorentzian. More technical details can be found in the Supplementary Note 3.

Code availabilty

Custom computer codes used in this study are available from the corresponding author upon reasonable request. Documentation of the codes is not available.

Data availability

Data are available from the corresponding author upon reasonable request.

Supplementary Material

Open In Web Browser

Footnotes

Electronic supplementary material

Supplementary Information accompanies this paper at doi:10.1038/s41467-017-01177-0.

Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Acknowledgements

The authors acknowledge useful discussions with J. Chaloupka, A. Banerjee, S.E. Nagler, A.A. Tsirlin, R. Moessner, F. Pollmann, and M. Zhitomirsky. S.M.W. acknowledges support through an NSERC Canada Postdoctoral Fellowship. R.V. and K.R. acknowledge support by the Deutsche Forschungsgemeinschaft through grant SFB/TR 49. The work of P.A.M. and A.L.C. was supported by the U.S. Department of Energy, Office of Science, Basic Energy Sciences under award no. DE-FG02-04ER46174.

Author contributions

R.V. and S.M.W. conceived the project. K.R., S.M.W. and A.H. performed and analyzed the ED calculations. P.A.M. and A.L.C. performed and analyzed the iDE results. All authors contributed equally to the manuscript.

Notes

Competing interests

The authors declare no competing financial interests.

Articles from Nature Communications are provided here courtesy of Nature Publishing Group

PMC Copyright Notice

The articles available from the PMC site are protected by copyright, even though access is free. Copyright is held by the respective authors or publishers who provide these articles to PMC. Users of PMC are responsible for complying with the terms and conditions defined by the copyright holder.

Users should assume that standard copyright protection applies to articles in PMC, unless an article contains an explicit license statement that gives a user additional reuse or redistribution rights. PMC does not allow automated/bulk downloading of articles that have standard copyright protection.

See the copyright notice on the PMC site, https://www.ncbi.nlm.nih.gov/pmc/about/copyright/, for further details and specific exceptions.

References

Lifshitz, E. M. & Pitaevskii, L. P. Statistical PhysicsPart 2 (Pergamon Press, Oxford, 1980).

Balents L, author. Spin liquids in frustrated magnets. Nature. 2010;464:199–208. [Cross Ref] [PubMed]

Coldea R, Tennant DA, Tsvelik AM, Tylczynski Z, authors. Experimental realization of a 2D fractional quantum spin liquid. Phys. Rev. Lett. 2001;86:1335–1338. [Cross Ref] [PubMed]

Han T-H, et al., authors. Fractionalized excitations in the spin-liquid state of a kagome-lattice antiferromagnet. Nature. 2012;492:406–410. [Cross Ref] [PubMed]

Kitaev A, author. Anyons in an exactly solved model and beyond. Ann. Phys. 2006;321:2–111. [Cross Ref]

Jackeli G, Khaliullin G, authors. Mott insulators in the strong spin-orbit coupling limit: from Heisenberg to a quantum compass and Kitaev models. Phys. Rev. Lett. 2009;102:017205 [Cross Ref] [PubMed]

Rau JG, Lee EK-H, Kee H-Y, authors. Spin-orbit physics giving rise to novel phases in correlated systems: iridates and related materials. Annu. Rev. Condens. Matter Phys. 2016;7:195–221. [Cross Ref]

Gretarsson H, et al., authors. Magnetic excitation spectrum of Na₂IrO₃ probed with resonant inelastic X-ray scattering. Phys. Rev. B. 2013;87:220407 [Cross Ref]

Hwan Chun S, et al., authors. Direct evidence for dominant bond-directional interactions in a honeycomb lattice iridate Na₂IrO₃. Nat. Phys. 2015;11:462–466. [Cross Ref]

Sandilands LJ, Tian Y, Plumb KW, Kim Y-J, Burch KS, authors. Scattering continuum and possible fractionalized excitations in α-RuCl₃. Phys. Rev. Lett. 2015;114:147201 [Cross Ref] [PubMed]

Nasu J, Knolle J, Kovrizhin DL, Motome Y, Moessner R, authors. Fermionic response from fractionalization in an insulating two-dimensional magnet. Nat. Phys. 2016;12:912–915. [Cross Ref]

Banerjee A, et al., authors. Neutron scattering in the proximate quantum spin liquid α-RuCl₃. Science. 2017;356:1055–1059. [Cross Ref] [PubMed]

Banerjee A, et al., authors. Proximate Kitaev quantum spin liquid behaviour in a honeycomb magnet. Nat. Mater. 2016;15:733–740. [Cross Ref] [PubMed]

Do, S.-H. et al. Incarnation of Majorana fermions in Kitaev quantum spin lattice. Preprint at https://arxiv.org/abs/1703.01081 (2017).

Glamazda A, Lemmens P, Do S-H, Choi YS, Choi K-Y, authors. Raman spectroscopic signature of fractionalized excitations in the harmonic-honeycomb iridates β- and γ-Li₂IrO₃. Nat. Commun. 2016;7:12286 [Cross Ref] [PubMed]

Chaloupka J, Jackeli G, Khaliullin G, authors. Kitaev-Heisenberg model on a honeycomb lattice: possible exotic phases in iridium oxides A₂IrO₃. Phys. Rev. Lett. 2010;105:027204 [Cross Ref] [PubMed]

Rau JG, Lee EK-H, Kee H-Y, authors. Generic spin model for the honeycomb iridates beyond the Kitaev limit. Phys. Rev. Lett. 2014;112:077204 [Cross Ref] [PubMed]

Winter SM, Li Y, Jeschke HO, Valent R, authors. Challenges in design of Kitaev materials: magnetic interactions from competing energy scales. Phys. Rev. B. 2016;93:214431 [Cross Ref]

Sears JA, et al., authors. Magnetic order in α-RuCl₃: a honeycomb-lattice quantum magnet with strong spin-orbit coupling. Phys. Rev. B. 2015;91:144420 [Cross Ref]

Chaloupka J, Khaliullin G, authors. Magnetic anisotropy in the Kitaev model systems Na₂IrO₃ and RuCl₃. Phys. Rev. B. 2016;94:064435 [Cross Ref]

Kimchi I, Coldea R, Vishwanath A, authors. Unified theory of spiral magnetism in the harmonic-honeycomb iridates α, β, and γLi₂IrO₃. Phys. Rev. B. 2015;91:245134 [Cross Ref]

Kim H-S, Lee EK-H, Kim Y-B, authors. Predominance of the Kitaev interaction in a three-dimensional honeycomb iridate: from ab initio to spin model. EPL. 2015;112:67004 [Cross Ref]

Cao HB, et al., authors. Low-temperature crystal and magnetic structure of α-RuCl₃. Phys. Rev. B. 2016;93:134423 [Cross Ref]

Ran K, et al., authors. Spin-wave excitations evidencing the Kitaev interaction in single crystalline α-RuCl₃. Phys. Rev. Lett. 2017;118:107203 [Cross Ref] [PubMed]

Kim H-S, Kee H-Y, authors. Crystal structure and magnetism in α-RuCl₃: an ab initio study. Phys. Rev. B. 2016;93:155143 [Cross Ref]

Yadav R, et al., authors. Kitaev exchange and field-induced quantum spin-liquid states in honeycomb α-RuCl₃. Sci. Rep. 2016;6:37925 [Cross Ref] [PubMed]

Wang W, Dong Z-Y, Yu S-L, Li J-X, authors. Theoretical investigation of the magnetic dynamics and superconducting pairing symmetry in α-RuCl₃. Phys. Rev. B. 2017;96:115103 [Cross Ref]

Hou YS, Xiang HJ, Gong XG, authors. Unveiling magnetic interactions of ruthenium trichloride via constraining direction of orbital moments: potential routes to realize a quantum spin liquid. Phys. Rev. B. 2017;96:054410 [Cross Ref]

Katukuri VM, et al., authors. Kitaev interactions between j = 1/2 moments in honeycomb Na₂IrO₃ are large and ferromagnetic: insights from ab initio quantum chemistry calculations. New J. Phys. 2014;16:013056 [Cross Ref]

Gohlke, M., Verresen, R., Moessner, R. & Pollmann, F. Dynamics of the Kitaev-Heisenberg model. Preprint at https://arxiv.org/abs/1701.04678 (2017).

Johnson RD, et al., authors. Monoclinic crystal structure of α-RuCl₃ and the zigzag antiferromagnetic ground state. Phys. Rev. B. 2015;92:235119 [Cross Ref]

Harris AB, Kumar D, Halperin BI, Hohenberg PC, authors. Dynamics of an antiferromagnet at low temperatures: spin-wave damping and hydrodynamics. Phys. Rev. B. 1971;3:961–1024. [Cross Ref]

Chernyshev AL, Zhitomirsky ME, authors. Magnon decay in noncollinear quantum antiferromagnets. Phys. Rev. Lett. 2006;97:207202 [Cross Ref] [PubMed]

Zhitomirsky ME, Chernyshev AL, authors. Colloquium: spontaneous magnon decays. Rev. Mod. Phys. 2013;85:219–242. [Cross Ref]

Ruvalds J, Zawadowski A, authors. Two-phonon resonances and hybridization of the resonance with single-phonon states. Phys. Rev. B. 1970;2:1172–1175. [Cross Ref]

Kosevich, A. M. The Crystal Lattice: Phonons, Solitons, Dislocations, Superlattices (Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim, 2005).

Sizyuk Y, Wölfle P, Perkins NB, authors. Selection of direction of the ordered moments in Na₂IrO₃ and α-RuCl₃. Phys. Rev. B. 2016;94:085109 [Cross Ref]

Chaloupka J, Khaliullin G, authors. Hidden symmetries of the extended Kitaev-Heisenberg model: implications for the honeycomb-lattice iridates A₂IrO₃. Phys. Rev. B. 2015;92:024413 [Cross Ref]

Hwang K, Kim YB, authors. Theory of triplon dynamics in the quantum magnet BiCu₂PO₆. Phys. Rev. B. 2016;93:235130 [Cross Ref]

Maksimov PA, Zhitomirsky ME, Chernyshev AL, authors. Field-induced decays in XXZ triangular-lattice antiferromagnets. Phys. Rev. B. 2016;94:140407 [Cross Ref]

Mourigal M, Fuhrman WT, Chernyshev AL, Zhitomirsky ME, authors. Dynamical structure factor of the triangular-lattice antiferromagnet. Phys. Rev. B. 2013;88:094407 [Cross Ref]

Biffin A, et al., authors. Noncoplanar and counterrotating incommensurate magnetic order stabilized by Kitaev interactions in γ-Li₂IrO₃. Phys. Rev. Lett. 2014;113:197201 [Cross Ref] [PubMed]

Biffin A, et al., authors. Unconventional magnetic order on the hyperhoneycomb Kitaev lattice in β-Li₂IrO₃: full solution via magnetic resonant X-ray diffraction. Phys. Rev. B. 2014;90:205116 [Cross Ref]

Williams SC, et al., authors. Incommensurate counterrotating magnetic order stabilized by Kitaev interactions in the layered honeycomb α-Li₂IrO₃. Phys. Rev. B. 2016;93:195158 [Cross Ref]

Takayama T, et al., authors. Hyperhoneycomb iridate β-Li₂IrO₃ as a platform for Kitaev magnetism. Phys. Rev. Lett. 2015;114:077202 [Cross Ref] [PubMed]

Kim H-S, Kim YB, Kee H-Y, authors. Revealing frustrated local moment model for pressurized hyperhoneycomb iridate: paving the way toward a quantum spin liquid. Phys. Rev. B. 2016;94:245127 [Cross Ref]

Mandal S, Bhattacharjee S, Sengupta K, Shankar R, Baskaran G, authors. Confinement-deconfinement transition and spin correlations in a generalized Kitaev model. Phys. Rev. B. 2011;84:155121 [Cross Ref]

Schaffer R, Bhattacharjee S, Kim Y-B, authors. Quantum phase transition in Heisenberg-Kitaev model. Phys. Rev. B. 2012;86:224417 [Cross Ref]

Cromer DT, Waber JT, authors. Scattering factors computed from relativistic Dirac-Slater wave functions. Acta Cryst. 1965;18:104–109. [Cross Ref]

Lanczos C, author. An iteration method for the solution of the eigenvalue problem of linear differential and integral operators. J. Res. Nat. Bur. Stand. 1950;45:255–282. [Cross Ref]

Dagotto E, author. Correlated electrons in high-temperature superconductors. Rev. Mod. Phys. 1994;66:763–840. [Cross Ref]

Toth S, Lake B, authors. Linear spin-wave theory for single-Q incommensurate magnetic structures. J. Phys. Condens. Matter. 2015;27:166002 [Cross Ref] [PubMed]

Knolle J, Kovrizhin DL, Chalker JT, Moessner R, authors. Dynamics of fractionalization in quantum spin liquids. Phys. Rev. B. 2015;92:115127 [Cross Ref]

Knolle J, Kovrizhin DL, Chalker JT, Moessner R, authors. Dynamics of a two-dimensional quantum spin liquid: signatures of emergent Majorana fermions and fluxes. Phys. Rev. Lett. 2014;112:207203 [Cross Ref]

[Back]

Fig. 1.

From material to model. Within the honeycomb ab-layer of α-RuCl₃ are illustrated the RuCl₆ octahedra, magnetic zigzag ordering pattern, and definition of the underlying magnetic interactions. Crystal axes are labeled with respect to the C2/m structure

[Back]

[Back]

Fig. 2.

Neutron scattering intensity ℐ(k, ω) within the nnHK model. a–c Detailed results for Model 1 (J ₁ = −2.2, K ₁ = +7.4 meV): a ℐ(k, ω) computed via LSWT; results are averaged over the three zigzag ordering wavectors, parallel to the X-bonds, Y-bonds, and Z-bonds. Inset: Definition of Brillouin zone and high-symmetry k-points. b ED results, combining data from several 20-site and 24-site periodic clusters (see “Methods”). c ED k-dependence of ℐ(k, ω) integrated over the indicated energies, as obtained from a single 24-site cluster respecting all symmetries of Eq. (1) (see “Methods”). d Comparison of Γ-point intensities for the K ₁ = +7.7 meV AFM spin liquid (exact results⁵³, ⁵⁴), Model 1 (ED), and the experimental data for α-RuCl₃ ¹². e Evolution of the ED Γ-point intensity with decreasing ∣J₁∕K₁∣, showing absence of low-energy intensity close to the K ₁ > 0 spin liquid. The top three interaction sets correspond to zigzag order, while the bottom is the K ₁ > 0 spin liquid. For all spectra, a Gaussian broadening of 0.5 meV has been applied

[Back]

[Back]

Fig. 3.

Neutron scattering intensity ℐ(k, ω) within the extended model. a–c Detailed results for Model 2 (J ₁ = −0.5, K ₁ = −5.0, Γ ₁ = +2.5, J ₃ = +0.5 meV): a ℐ(k, ω) computed via LSWT; results are averaged over the three zigzag domains with ordering wavectors parallel to the X-bonds, Y-bonds, and Z-bonds. b ED results, combining data from several 20-site and 24-site periodic clusters (see “Methods”). c ED k-dependence of ℐ(k, ω) integrated over the indicated energies, as obtained from a single 24-site cluster respecting all symmetries of Eq. (1) (see “Methods”). d Comparison of Γ-point intensities for the K ₁ = −7.7 meV FM spin liquid (exact results⁵³, ⁵⁴), Model 2 (ED), and the experimental data for α-RuCl₃ ¹². e Evolution of the ED Γ-point intensity with decreasing ∣Γ₁∕K₁∣, showing significant broadening at finite Γ ₁. The top three interaction sets correspond to zigzag order, while the bottom is the K ₁ < 0 spin liquid. For all spectra, a Gaussian broadening of 0.5 meV has been applied

[Back]

[Back]

Fig. 4.

Magnon decay channels for the nnHK model. a Local picture of quantum fluctuations away from zigzag order. The energy cost for the left process vanishes on approaching the spin liquid ∣J₁∕K₁∣→0. b Momentum space picture for the corresponding fourth-order decay process due to ℋ₄. c Ordered moment direction for Model 1 (J ₁ = −2.2, K ₁ = +7.4 meV), corresponding to the zigzag domain with ordering wavevector Q = Y. d LSWT dispersions ϵk,m, and three-magnon density of states (DOS) for Model 1 for the same zigzag domain as c. The dashed line indicates the bottom of the three-magnon continuum E3min, which is coincident with the lowest magnon band

[Back]

[Back]

Fig. 5.

Magnon decay channels for the extended model. a Local picture of additional quantum fluctuations away from zigzag order induced by Γ ₁ interactions. b Momentum space picture of the third-order decay process ℋ₃. c Ordered moment direction for Model 2 (J ₁ = −0.5, K ₁ = −5.0, Γ ₁ = +2.5, J ₃ = +0.5 meV) with zigzag ordering wavevector Q = Y, parallel to the Z-bond. d LSWT dispersions ϵk,m, and two-magnon DOS for Model 2 with Q = Y. Dashed lines indicate the bottom of the two-magnon and three-magnon continuum (E2min(k) and E3min(k), respectively)

[Back]

[Back]

Fig. 6.

Effects of two-magnon decays in ℐ(k, ω) for extended model. Results are shown for Model 2 computed via a LSWT, b self-consistent iDE approach, and c ED. Results in a and b are averaged over the different zigzag domains. The white and pink dashed lines indicate the bottom of the two-magnon continuum, E2min(k) for the different zigzag domains. In the iDE results, the effects of two-magnon decays strongly broadens any magnon bands overlapping with the two-magnon continuum

[Back]