---
title: 'Dirac Magnons: Topological Bosonic Excitations'
url: https://www.emergentmind.com/topics/dirac-magnons
type: topic
---

# Dirac Magnons: Topological Bosonic Excitations

Dirac magnons are bosonic spin-wave excitations whose band structures exhibit linear crossings analogous to Dirac fermions in electronic systems. In the simplest realization, a ferromagnet on a honeycomb lattice yields two magnon branches that touch linearly at the Brillouin-zone corners \(K\) and \(K'\); in three-dimensional antiferromagnets with combined inversion and time-reversal symmetry, \(PT\), a Dirac magnon node can appear as a fourfold crossing of two doubly degenerate magnon branches [1512.04902] [1711.02960]. Across the literature, the term therefore denotes a symmetry-protected bosonic nodal excitation rather than a single lattice type or dimensionality, and it extends naturally to gapped topological Dirac magnons, nodal loops, nodal lines, nodal planes, and related multicomponent crossings [1703.08545] [1803.06349].

## 1. Definition and elementary band-theoretic structure

In honeycomb ferromagnets, the defining structural ingredient is the two-site basis. Linear spin-wave theory then produces a \(2\times 2\) magnon Hamiltonian in sublattice pseudospin space, directly paralleling graphene. For the nearest-neighbor Heisenberg model, one obtains two magnon branches
\[
\varepsilon_{\pm}(\mathbf k)=JS\left(3\pm |\gamma_{\mathbf k}|\right),
\]
with \(\gamma_{\mathbf k}\) the honeycomb structure factor. Because \(\gamma_{\mathbf K}=\gamma_{\mathbf K'}=0\), the two bands touch linearly at \(K\) and \(K'\), yielding Dirac magnons at finite positive energy rather than at a Fermi level [1706.03384].

The same logic can be written in a form convenient for topological perturbations. In the honeycomb ferromagnet with next-nearest-neighbor Dzyaloshinskii–Moriya interaction (DMI), the single-magnon Hamiltonian becomes
\[
M_{\mathbf k}= \begin{pmatrix} 3JS-2DS\,\beta_{\mathbf k} & -JS\,\gamma_{\mathbf k}\ -JS\,\gamma_{\mathbf k}^* & 3JS+2DS\,\beta_{\mathbf k} \end{pmatrix},
\]
with magnon energies
\[
\omega_{\mathbf k}=3JS\pm\sqrt{J^2S^2|\gamma_{\mathbf k}|^2+4D^2S^2\beta_{\mathbf k}^2}.
\]
For \(D=0\), the bands touch at \(K\) and \(K'\); for \(D\neq 0\), a gap
\[
\Delta_g = 6\sqrt{3}\,DS
\]
opens at the Dirac points [2207.03964]. This distinction underlies the common classification into massless and massive Dirac magnons.

The concept is not restricted to honeycomb ferromagnets. The earliest general honeycomb analysis showed that localized spins on a two-dimensional honeycomb lattice naturally generate magnon bands with Dirac-like crossings, with the Dirac points at \(K\) and \(K'\) in the ferromagnetic case and at \(\Gamma\) in the antiferromagnetic case [1512.04902]. A common misconception is therefore that Dirac magnons are intrinsically tied to ferromagnetic honeycomb layers; the literature instead treats them as a broader class of bosonic linear band touchings whose detailed form is fixed by lattice geometry, magnetic order, and symmetry.

## 2. Symmetry protection and nodal taxonomy

The simplest honeycomb Dirac magnons are protected by the two-sublattice structure and the vanishing of the inter-sublattice structure factor at \(K\) and \(K'\). In three-dimensional antiferromagnets, however, the central protecting mechanism can instead be magnetic \(PT\) symmetry. For a collinear \(PT\)-symmetric antiferromagnet with spin \(U(1)\) symmetry, the quadratic magnon Hamiltonian decomposes as
\[
H = H_+ \oplus H_-,
\qquad
H_+ = H_-^\ast,
\]
so a Weyl point in one spin sector is accompanied by an opposite-charge partner in the other sector at the same momentum. The full crossing is then a fourfold Dirac magnon node [1703.08545]. In this setting, Dirac magnons are not merely accidental point degeneracies but composites of symmetry-related Weyl magnons.

Once \(U(1)\) is broken while \(PT\) remains intact, the same three-dimensional crossings need not gap. Instead, the generic descendants are nodal lines with nontrivial \(Z_2\) charge, and the Dirac point becomes the zero-radius limit of a nodal line [1703.08545]. This mechanism is conceptually distinct from the honeycomb DMI gap problem: the same label, “Dirac magnon,” covers both a two-band graphene-like cone and a fourfold \(PT\)-protected three-dimensional crossing.

The taxonomy becomes broader in lower dimensions as well. On the two-dimensional CaVO lattice, a collinear antiferromagnet with \(\mathcal{PT}\) symmetry yields eight doubly degenerate magnon branches whose linear crossings form one-dimensional closed loops in the two-dimensional Brillouin zone. These are Dirac nodal loop magnons rather than isolated Dirac points, and their topology is diagnosed by a parity-based \(\mathbb Z_2\) invariant with \(\nu=1\) [1803.06349]. On the anisotropic square-hexagon-octagon lattice, type-II Dirac magnon points arise without requiring DMI, and their overtilted cones coexist with Dirac nodal loops whose protection is characterized by a \(\mathbb Z_2\) invariant and Wilson-loop analysis [2305.16419].

A further extension concerns the internal spin texture of the magnon eigenmodes. In a kagome \(XY\) antiferromagnet with \(120^\circ\) order, finite-energy Dirac points at \(K\) and \(K'\) are accompanied by a momentum-space magnon spin vortex with winding number
\[
Q=+1,
\]
whereas a related kagome model exhibits a \(Q=-2\) vortex around \(\Gamma\) [1702.02575]. This establishes that Dirac magnons may be characterized not only by linear dispersion but also by a Dirac-like winding of the one-magnon spin expectation value.

## 3. Microscopic Hamiltonians and calculational frameworks

The microscopic starting point is usually an exchange Hamiltonian of Heisenberg type,
\[
H=\sum_{\langle ij\rangle} J_{ij}\,\mathbf S_i\cdot \mathbf S_j,
\]
sometimes augmented by anisotropy, DMI, or Zeeman terms. For Cu\(_3\)TeO\(_6\), the experimentally fitted model is an isotropic Heisenberg Hamiltonian extending to sixth neighbors,
\[
H=\sum_{n=1}^6 \sum_{\langle ij\rangle\in\{n\text{-NN}\}} J_n\, {\bf S}_i\cdot{\bf S}_j ,
\]
with best-fit parameters
\[
J_1=9.07,\quad J_2=0.89,\quad J_3=-1.81,\quad J_4=1.91,\quad J_5=0.09,\quad J_6=1.83~\text{meV},
\]
and \(J_1\) the dominant scale [1711.02960]. For elemental Gd, the basic topological magnon structure is already present in a two-band hcp ferromagnetic model with exchange constants \(J_1\), \(J_2\), and \(J_3\), where the band splitting is controlled by an inter-sublattice structure factor \(g(\mathbf k)\) and nodal manifolds occur where \(g(\mathbf k)=0\) [2107.11372].

The standard calculational route is linear spin-wave theory. After a Holstein–Primakoff expansion about the ordered state, one diagonalizes either an ordinary quadratic bosonic Hamiltonian or, in antiferromagnets, a bosonic Bogoliubov problem. The CaVO antiferromagnet is a representative example: its quadratic spin-wave Hamiltonian contains the usual pairing terms \(a^\dagger b^\dagger+ab\), so diagonalization requires a Bogoliubov transformation rather than a simple unitary one [1803.06349]. For noncollinear orders, one first rotates to local spin frames and then performs the same expansion [1703.08545].

When comparison to inelastic neutron scattering is required, the theory must also produce the neutron cross section. In Cu\(_3\)TeO\(_6\), the measured intensity is compared with
\[
\frac{d^2\sigma}{d\Omega dE}\propto \sum_{\alpha\beta}\left(\delta_{\alpha\beta}-\frac{Q_\alpha Q_\beta}{\bf Q^2}\right) S^{\alpha\beta}({\bf Q},E),
\]
where
\[
S^{\alpha\beta}({\bf Q},E)=\frac{1}{N}\sum_{ij}e^{i{\bf Q}\cdot({\bf r}_i-{\bf r}_j)} \int_{-\infty}^{\infty}\langle S_i^\alpha S_j^\beta(t)\rangle e^{-iEt}\,dt.
\]
These expressions are not topological invariants, but they are the quantitative link between symmetry-based magnon band theory and spectroscopic evidence [1711.02960].

Beyond linear theory, several many-body frameworks have been developed. For honeycomb Dirac magnons with DMI, a spinor Green’s-function formalism yields a Dyson equation up to second order,
\[
G_R=G_R^{(0)}+G_R^{(0)}\Sigma^{(1)}G_R^{(0)}+G_R^{(0)}\Sigma^{(2)}G_R^{(0)},
\]
capturing static self-energy renormalization, momentum-dependent band reshaping, and damping [2207.03964]. A separate pumping-based treatment of a honeycomb ferromagnet constructs a bosonic analog of the Cooper ladder for opposite-momentum Dirac magnons, with a pump-induced pairing scale
\[
\Delta^2 = \frac{J}{4}\left(\frac{\Gamma\sqrt{S}}{3SJ}\right)^2
\]
at the Dirac resonance \(\Omega=3SJ\) [2012.12051].

## 4. Material realizations and spectroscopic evidence

Experimental realizations now span two-dimensional, quasi-two-dimensional, and three-dimensional magnets. The table summarizes representative systems discussed in the literature.

| System | Magnetic setting | Salient Dirac-magnon phenomenology |
|---|---|---|
| Cu\(_3\)TeO\(_6\) [1711.02960] | 3D collinear antiferromagnet with \(PT\) symmetry | Dirac nodes at \(P\) and coexistence with triply degenerate magnons |
| CoTiO\(_3\) [1907.02061] | Stacked honeycomb quantum \(XY\) magnet | Gapless Dirac cone at \(K\) near \(\hbar\omega^\star\sim 8.5\) meV |
| CrCl\(_3\) [2110.10894] | Quasi-2D honeycomb ferromagnet | Massless Dirac magnons; phase interpreted as not topological |
| CrBr\(_3\) [2204.11355] | Honeycomb ferromagnet | Ideal Dirac crossing at \(K\), isospin winding, \(T^2\) thermal renormalization |
| CrI\(_3\) [2605.06291] | Quasi-2D honeycomb ferromagnet | Gapped Dirac magnons near \(K\) with winding feature around the gap |
| Gd [2107.11372] | hcp ferromagnet | Dirac magnons at \(K\), nodal lines along \(K\)-\(H\)-\(K\), nodal planes at half-integer \(\ell\) |

Cu\(_3\)TeO\(_6\) established the three-dimensional antiferromagnetic case experimentally. Inelastic neutron scattering at \(T=5\) K identified a Dirac node at \(P\) near \(E\approx 15\) meV and triply degenerate nodes at \(H\) near \(E\approx 16\) meV and \(18.5\) meV, with the spectrum well described by linear spin-wave theory based on the fitted \(J_1\)-\(J_6\) model [1711.02960]. The paper is explicit that what is directly measured is the magnon dispersion and the nodal crossings in \(({\bf Q},E)\) space; Berry curvature, Chern numbers, and surface arcs are inferred theoretically rather than observed directly.

CoTiO\(_3\) provided a distinct route: a stacked honeycomb lattice magnet with strong easy-plane exchange anisotropy and effective pseudospin-\(1/2\) physics. Its low-energy magnons exhibit a gapless Dirac cone at the honeycomb \(K\) points, within an overall magnon bandwidth of approximately \(12\) meV, and the fitted nearest-neighbor intralayer exchange is almost purely \(XY\),
\[
J_{\parallel,1}=-4.4(9)~{\rm meV}, \qquad J_{\perp,1}=0~{\rm meV},
\]
with smaller interlayer couplings of about \(0.6\) meV [1907.02061].

The chromium trihalides separate the massless and gapped regimes. CrCl\(_3\) is described by a nearly two-dimensional Heisenberg model with fitted in-plane exchanges
\[
J_1=-0.95\ \mathrm{meV},\qquad J_2=-0.024\ \mathrm{meV},\qquad J_3=0.051\ \mathrm{meV},
\]
while the estimated interlayer coupling and anisotropy,
\[
J_c=0.0077\ \mathrm{meV},\qquad D=0.014\ \mathrm{meV},
\]
are too small to matter at the Dirac scale. The Dirac crossing is inferred at \(|\mathbf Q_{1/3,1/3,0}|=0.70\ \text{\AA}^{-1}\) and \(\hbar\omega_D\approx 4.4\) meV, with no observable gap and no anomalous broadening; the phase is therefore interpreted as massless and not topological [2110.10894]. CrBr\(_3\), by contrast, shows an ideal Dirac crossing at \(K\) at
\[
E_D=7.5\pm 0.1~\text{meV},
\]
together with the neutron-scattering signature of isospin winding,
\[
I = I_0 \big[1 \pm \cos(\alpha-\alpha_0)\big],
\]
around the Dirac point [2204.11355].

CrI\(_3\) realizes the gapped case. A recent single-crystal study resolved a Dirac magnon gap near \(K\) centered around \(12\) meV with size
\[
\Delta_K \sim 2~\text{meV},
\]
and, crucially, observed a momentum-space winding of neutron intensity around \(K\): spectral weight lies outside the first Brillouin zone above the gap and inside it below the gap, with angular cuts fit by phase-shifted cosine functions [2605.06291]. In structurally disordered CrI\(_3\), neutron scattering on a mixed rhombohedral/monoclinic sample showed that antiferromagnetic and ferromagnetic stackings can coexist, with the monoclinic antiferromagnetic component shifting magnon energies downward by just under \(1\) meV while leaving a gap at the Dirac-point energy; the paper therefore argued cautiously that the nontrivial magnon topology of rhombohedral CrI\(_3\) may persist in the monoclinic phase as well [2305.11750].

A recurring experimental issue is momentum resolution. In CrBr\(_3\), broad in-plane integration windows in time-of-flight analysis can generate an apparent splitting at \(K\), whereas narrow integration windows recover the continuous Dirac crossing [2204.11355]. In CrCl\(_3\), the evidence is powder averaged and therefore identifies the Dirac crossing through density-of-states fingerprints rather than a directly imaged cone [2110.10894]. These cases illustrate that the phrase “Dirac magnon” can refer either to a fully momentum-resolved cone or to a more indirect but model-consistent spectroscopic signature.

## 5. Interactions, pumping, and finite-temperature evolution

A central issue is the fate of the Dirac cone under magnon–magnon interactions. In ferromagnetic honeycomb layers, extending the theory beyond standard Dyson theory produces strongly momentum-dependent renormalization of the bare magnon bands and strongly momentum-dependent lifetimes, with qualitative agreement to unexplained anomalies in old neutron data for CrBr\(_3\) [1706.03384]. A later experimental reinvestigation of CrBr\(_3\) found that both the thermal band renormalization and the magnon linewidth follow the expected universal \(T^2\) evolution, with fitted exponents
\[
\alpha_{\mathrm{Energy}}=2.25(15),\qquad \alpha_{\mathrm{Width}}=1.95(14),
\]
although the measured momentum dependence lacks the sharp van Hove features predicted by the available low-order interacting spin-wave theory [2204.11355]. In CrI\(_3\), the temperature evolution of several magnon modes likewise fits
\[
E(T)=a+bT^\alpha
\]
with exponents \(2.1(4)\), \(2.4(3)\), and \(2.0(5)\), consistent with \(T^2\)-type renormalization from magnon–magnon interactions [2605.06291].

More elaborate many-body effects arise when DMI and nonequilibrium populations are present. For topological honeycomb Dirac magnons, the first-order self-energy can renormalize the Dirac velocity-like term, shift the Dirac points in momentum space through an interaction-generated field \(g(\mathbf k)\), and modify the effective Haldane mass through an interaction-induced term \(m_H\). The renormalized spectrum takes the form
\[
\omega_{\mathbf k} = 3JS+P \pm \sqrt{\left|(JS-Q)\gamma_{\mathbf k}+g(\mathbf k)\right|^2+\left(2DS\beta_{\mathbf k}-m_H\right)^2},
\]
so interactions reshape not only energies but also Berry curvature and, under strong pumping, can even generate topological flat bands and Chern-number inversion [2207.03964].

Pumping opens a separate nonequilibrium route. A theory of a uniformly pumped honeycomb ferromagnet showed that a direct one-photon process cannot create a single Dirac magnon because the pump carries zero crystal momentum, but a second-order Suhl process can create pairs of opposite-momentum Dirac magnons. The resonant frequency is the Dirac-point energy,
\[
\Omega = 3SJ,
\]
and the resulting instability was interpreted as formation of a Dirac magnon paired state with zero or reduced magnetization, with CrBr\(_3\) and CrCl\(_3\) proposed as candidate materials near the Curie temperature [2012.12051]. This does not describe an equilibrium condensate; it is a pump-induced instability of a driven bosonic Dirac system.

A particularly important current controversy concerns thermally driven topological phase transitions in gapped honeycomb Dirac magnons. A Hartree-based treatment predicted that, in a 2D honeycomb ferromagnet with DMI and Zeeman field, increasing temperature or field can close and reopen the Dirac gap at a critical temperature \(T_c\), reverse the Chern numbers, and induce a sign change of the thermal Hall conductivity \(\kappa_{xy}\) [2104.00266]. A later microscopic theory based on nonlinear spin-wave theory with T-matrix resummation reached the opposite conclusion for CrI\(_3\), CrSiTe\(_3\), and CrGeTe\(_3\): it found a thermally induced reduction of the topological magnon gap but no evidence of thermally driven topological transitions, and it identified two-magnon bound-state resonances as an important source of low-temperature damping [2509.13900]. This suggests that the thermal fate of topological Dirac magnons is approximation-sensitive and that linewidth broadening should not be conflated with gap closing.

## 6. Topological responses, emergent gauge fields, and open directions

Once a gap opens at a Dirac point, the magnon bands acquire Berry curvature and may support Chern-band responses. In Cu\(_3\)TeO\(_6\), a magnetic field is expected to split a Dirac point into two Weyl points with monopole charges \(+1\) and \(-1\), leading to thermal Hall conductivity; similarly, gapping the triply degenerate nodes is expected to produce thermal Hall contributions from bands carrying Chern numbers \(+2\) and \(-2\) [1711.02960]. In honeycomb ferromagnets with DMI, the thermal Hall conductivity is commonly written as
\[
\kappa^{xy} = -\frac{k_B^2T}{(2\pi)^2} \sum_n\int d^2\mathbf k\, c_2\!\left(f(\omega_n(\mathbf k))\right)\Omega_n^{xy}(\mathbf k),
\]
and under interaction-driven band reshaping its sign can be tuned by the pumped magnon population [2207.03964].

Several works extend the response theory beyond thermal Hall transport. In the honeycomb Dirac-magnon continuum coupled to emergent gauge fields, the effective Euclidean Lagrangian
\[
\mathcal{L} = \bar{\psi} \left(i\hbar\tilde{\gamma}^\mu\partial_\mu - m\right)\psi - g \bar{\psi} \tilde{\gamma}^\mu  \psi \mathcal{A}_\mu
\]
yields an induced current
\[
j^\mu(p) = \Pi^{\mu\nu}(p)\,\mathcal{A}_\nu(p),
\]
and in the DC long-wavelength limit the transverse spin conductivity becomes
\[
\sigma^{xy} = \alpha^2 \operatorname{sgn}(m)\frac{\hbar}{4\pi},
\]
a magnonic analog of the quantum Hall effect [2512.13362]. In the optical limit, the same theory predicts a sharp resonance when
\[
\hbar\omega = \Delta = 2m,
\]
interpreted as an interband transition of gapped Dirac magnons [2512.13362].

Elastic degrees of freedom provide another route. In a honeycomb ferromagnet, strain modulates the exchange and generates elastic gauge fields
\[
A_x^{el}=\frac{\beta JS}{v_J}(u_{xx}-u_{yy}), \qquad
A_y^{el}=-\frac{2\beta JS}{v_J}u_{xy},
\]
which couple to the Dirac magnons like valley-odd pseudogauge fields. For deformations producing constant pseudomagnetic fields, the Dirac-magnon spectrum splits into pseudo-Landau levels
\[
E_n=\pm \sqrt{2\hbar |eB|\,v_J^2\,n}.
\]
When DMI opens a topological gap, integrating out the thermally occupied Dirac magnons generates a Chern–Simons term for the elastic gauge field and a phonon Hall viscosity
\[
\eta^H = -\frac{4\hbar\beta^2\tilde C_D}{9\pi a^2},
\]
which vanishes at \(T=0\) and grows with temperature [1711.08653].

Surface and boundary manifestations remain an active theme. In Gd, linear spin-wave theory on finite slabs yields a surface magnon mode associated with the nodal lines, while the neutron intensity around the Dirac cone exhibits a sinusoidal angular modulation that reverses above and below the crossing, reflecting the momentum-dependent internal phase of the magnon eigenvectors [2107.11372]. In CrBr\(_3\) and CrI\(_3\), the experimentally resolved intensity winding around \(K\) similarly probes the pseudospin texture of the Dirac magnons rather than only the eigenvalues [2204.11355] [2605.06291]. A plausible implication is that future progress will depend increasingly on probes sensitive to eigenvectors, not only dispersions.

Open directions stated across the literature are relatively consistent. They include direct detection of surface arc states in three-dimensional topological magnon materials, quantitative tests of field-induced Weyl-magnon scenarios, controlled separation of thermal broadening from genuine topological gap evolution, and exploitation of stacking, strain, pumping, or gauge fields as external knobs for band engineering [1711.02960] [2305.11750] [2509.13900] [2512.13362]. The broader record already makes one point clear: Dirac magnons are not a single phenomenon but a family of symmetry-protected bosonic nodal excitations whose manifestations range from graphene-like honeycomb cones to \(PT\)-protected fourfold nodes, nodal loops, and gauge-responsive topological magnon bands.

Source: https://www.emergentmind.com/topics/dirac-magnons