# Dipolar spin-waves and tunable band gap at the Dirac points in the 2D magnet ErBr3

Christian Wessler, Bertrand Roessli, Karl W. Krämer, Uwe Stuhr, Andrew Wildes, Hans B. Braun, Michel Kenzelmann

► To cite this version:

Christian Wessler, Bertrand Roessli, Karl W. Krämer, Uwe Stuhr, Andrew Wildes, et al., Dipolar spin-waves and tunable band gap at the Dirac points in the 2D magnet ErBr3. Communications Physics, 2022, 5 (1), pp.185. 10.1038/s42005-022-00965-5. [hal-05442045]

HAL Id: hal-05442045

https://hal.science/hal-05442045v1

Submitted on 17 Feb 2026

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons CC BY 4.0 - Attribution - International License
<figure><img src="image_2.png" /></figure>

# Dipolar spin-waves and tunable band gap at the Dirac points in the 2D magnet ErBr<sub>3</sub>

Christian Wessler<sup>1,2✉</sup>, Bertrand Roessli<sup>1</sup>, Karl W. Krämer<sup>1,3</sup>, Uwe Stuhr<sup>1,3</sup>, Andrew Wildes<sup>4</sup>,
Hans B. Braun<sup>3</sup> & Michel Kenzelmann<sup>1,2</sup>

Topological magnon insulators constitute a growing field of research for their potential use as information carriers without heat dissipation. We report an experimental and theoretical study of the magnetic ground-state and excitations in the van der Waals two-dimensional honeycomb magnet ErBr<sub>3</sub>. We show that the magnetic properties of this compound are entirely governed by the dipolar interactions which generate a continuously degenerate non-collinear ground-state on the honeycomb lattice with spins confined in the plane. We find that the magnon dispersion exhibits Dirac-like cones when the magnetic moments in the ground-state are related by time-reversal and inversion symmetries associated with a Berry phase π as in single-layer graphene. A magnon band gap opens when the dipoles are rotated away from this state, entailing a finite Berry curvature in the vicinity of the K and K' Dirac points. Our results illustrate that the spin-wave dispersion of dipoles on the honeycomb lattice can be reversibly controlled from a magnetic phase with Dirac cones to a topological antiferromagnetic insulator with non-trivial valley Chern number.

<footnote n="1"><sup>1</sup>Laboratory for Neutron Scattering and Imaging, Paul Scherrer Institut, Villigen, Switzerland. <sup>2</sup>Department of Physics, Universität Basel, Basel, Switzerland.</footnote>

<footnote n="3"><sup>3</sup>Department of Chemistry, Biochemistry, and Pharmacy, University of Bern, Bern, Switzerland. <sup>4</sup>Institut Laue-Langevin, Grenoble, France. ✉email: christian.wessler@psi.ch</footnote>
Modern electronics makes use of the electronic charge to transport information. Because of electron scattering in metals and semiconductors, electronic devices experience resistance which results in significant energy loss. Due to Joule heating and dissipation miniaturization of transistors and integrated circuits is about to reach its technological limit. One way of overcoming these limitations is the use of the electron spin to store and process information. Research in the fast-developing field of spintronics mainly focuses on the development of technologies based on spin-polarized currents like spin valves that in the past have revolutionized the capacity of hard disk drives¹,². These efforts led to the emergence of the field of “magnonics” that promises low energy consumption. This technology uses spin waves (“magnons”) to process and store information and exploits the advantage that magnon-based devices can be easily tuned by an external magnetic field³,⁴. Among candidate materials for magnonic applications, magnetic van der Waals crystals are of special interest since a single layer can be exfoliated that allows us to build complex heterostructures⁵,⁶. One example is MnPS₃ that is considered for magnon transport⁷ and detection of anti-ferromagnetic states⁸. A key issue in developing magnonic materials is to control the direction of propagation of the spin waves at the surface of artificial nanostructures. The discovery of magnetic crystals where spin waves are robust at the surface and can propagate only along a single direction, because they are topologically protected, has prompted new ideas for future magnonic applications⁹. Non-trivial topological states can result from spin-orbit interactions that open a gap in the magnon Dirac cone¹⁰,¹¹. The spin-orbit coupling is at the origin of the anti-symmetric Dzyaloshinskii-Moriya (DM) interaction in magnetic crystals¹² that induces the magnon-Hall effect observed in ferromagnetic insulators¹³,¹⁴. Recently, it was proposed that dipolar interactions have a similar effect as the DM interaction¹⁵,¹⁶ that paves the way for the realization of topological artificial nanostructures¹⁷,¹⁸.

Dipolar interactions are always present in magnetic solids. They are direction dependent and long-range in contrast to exchange interactions. This property induces frustration in the lattice as in general it is not possible to satisfy all interactions simultaneously. The problem of determining whether a lattice of ordered magnetic dipoles can exist was addressed by Luttinger and Tisza¹⁹. They demonstrated that the face-centered and body-centered cubic lattices display ferromagnetism at low temperature while for the simple cubic lattice the classical ground state is a collinear antiferromagnet. The Mermin-Wagner theorem²⁰ does not apply in the case of long-range potentials and the absence of order in two-dimensional (2D) magnets with dipolar interactions does not generally apply. For a triangular lattice with solely dipolar interactions²¹, indeed long-range order is predicted. For both the square and honeycomb dipolar lattices, the classical ground state exhibits a continuous degeneracy²²,²³. Finally, in the geometrically frustrated Kagome lattice, no ground state satisfies the condition that the spin length is constant on all lattice points and long-range ferrimagnetic order is predicted with numerical simulations²⁴.

In this work, we report an analysis of the magnetic ground state and excitations in the van der Waals magnet ErBr₃ which was shown to undergo a magnetic phase transition to a continuously degenerate non-collinear 2D order below T = 280 mK²⁵. We show that the magnetic ground state in this crystal is explained by dipolar interactions between the Er³⁺ moments. Calculation of the dispersion of the magnetic excitations using the mean-field random-phase approximation (MF-RPA), which includes dipolar interactions and crystalline electric-field (CEF) anisotropy, describes the inelastic neutron-scattering (INS) data. Using our simulations of the magnetic excitations for the dipolar

honeycomb lattice, we show the presence of a Dirac cone at the K-points of the Brillouin zone when the configuration of dipoles possesses inversion (I) and time-reversal (T) symmetry. A gap opens at the Dirac points between the acoustic and lowest optical branch as the degeneracy angle between spins in the two sublattices of the honeycomb structure is varied and the IT symmetry of the magnetic structure is lifted.

## Results

Crystal and magnetic structure of ErBr₃. ErBr₃ crystallizes in the BiI₃ structure²⁶ with the rhombohedral space group R-3 (no. 148) and lattice parameters $a = b = 7.005$ Å and $c = 18.89$ Å at T = 1.5 K²⁵, as shown in Fig. 1a. In this layer-type structure the Br⁻ ions adopt a hexagonal densest packing (AB) and Er³⁺ ions occupy in an ordered way 2/3 of the octahedral voids between every other layer (2/3c). The layer stacking along the c-axis results in a sequence A2/3cBA2/3cBA2/3cB with three Er³⁺ honeycomb layers per unit cell. The Er³⁺ ions are located on site (6c) at (0, 0, ±z), (0, 0, ±z) + (2/3, 1/3, 1/3), and (0, 0, ±z) + (1/3, 2/3, 2/3). In this crystal z = 1/6 and the rare-earth ions form a perfect honeycomb lattice. The magnetic structure of ErBr₃ was determined by powder neutron diffraction²⁵. Magnetic order sets in below T = 280 mK and down to T = 50 mK only 2D order is observed. The magnetic ordering vector $k_0 = (1/3, 1/3, 0)$ triples the volume of the magnetic unit cell to $\sqrt{3}a \times \sqrt{3}a \times c$ with axes rotated by 30° from the crystallographic cell, see Fig. 1b. The ordered magnetic moments are located in the honeycomb plane and amount to $\mu = 4.7$ µB at T = 50 mK. The honeycomb lattice is composed of two trigonal sublattices. The Er³⁺ ions on each sublattice exhibit a 120° anti-ferromagnetic order with opposite chirality, see Fig. 1b. The magnetic ground state has an infinite degeneracy with respect to in-plane rotation of the ordered moments by an arbitrary angle Ψ with opposite signs for the two sublattices²⁵, see Fig. 1c.

In the mean-field approximation, the magnetic ground state is determined by the eigenvector of the largest eigenvalue of the Fourier transform of the exchange interaction tensor²⁷. The Hamiltonian with purely dipolar interactions is written

$$H = - \frac{(g \mu_B)^2 \mu_0}{8\pi} \sum_{i,j} \mathbf{J}_i \cdot \overline{\mathbf{D}}(ij) \cdot \mathbf{J}_j \quad (1)$$

and

$$D_{\alpha,\beta}(ij) = \frac{3(\mathbf{R}_i - \mathbf{R}_j)_\alpha (\mathbf{R}_i - \mathbf{R}_j)_\beta}{|\mathbf{R}_i - \mathbf{R}_j|^5} - \frac{\delta_{\alpha,\beta}}{|\mathbf{R}_i - \mathbf{R}_j|^3} \quad (2)$$

$\alpha, \beta = a, b, c$ are the crystal directions in Cartesian coordinates and $\mathbf{R}_i$ is the position vector of the i magnetic ion. The dispersion of the eigenvalues of the Fourier transform of the dipolar tensor along $\mathbf{k} = (q, q, 0)$ for the ErBr₃ crystal with six Er-ions in the chemical cell is shown in Supplementary Fig. 1 of Supplementary Note 1. The largest eigenvalue has a maximum at $k_0 = (1/3, 1/3, 0)$ in agreement with the propagation vector of the magnetic structure measured in ErBr₃. This particular non-collinear magnetic ordering cannot be reproduced by frustrated Heisenberg interactions²⁸. The temperature of the magnetic phase transition calculated with mean-field theory²⁹ $T_{MF} = \max\{\lambda_{k_0}/k_B\}/3 \sim 0.3$ K is close to the ordering temperature of the Er moments measured by neutron diffraction. Mean-field theory predicts the same magnetic ground state that is observed in ErBr₃ and agrees with the spin order obtained by Rozenbaum²³ for magnetic dipoles on the bipartite honeycomb lattice. The magnetic ground state is continuously degenerate and the equivalent spin structures are defined in terms of an arbitrary angle $\Psi$, as shown in Fig. 1c. Supplementary Fig. 2 shows the dependence of the ground state energy upon the lattice parameter $c$ of the crystal structure of ErBr₃. The dipolar energy is constant for a lattice parameter $c > 10$ Å
<figure><img src="image_3.png" /><figcaption>**Fig. 1** Crystal and magnetic structure of ErBr₃. a Crystal structure of ErBr₃. Br⁻ ions octahedrally coordinate Er³⁺ ions that form honeycomb layers stacked along the c-axis. b Honeycomb layer of Er³⁺ ions at z = 1/6 with ordered magnetic moments of 4.7 µB (black arrows) for a spin configuration given by angle Ψ = 0. The two trigonal sublattices of the honeycomb are shown as blue and red spheres. The crystallographic and magnetic unit cells are shown as black and green lines, respectively. c The classical magnetic ground state is continuously degenerate and the magnetic moments can be rotated in opposite direction by varying the angle Ψ. In panel c, Ψ = π/6.</figcaption></figure>

that indicates that interlayer interactions become negligibly small above that value. With a lattice parameter c = 18.89 Å and an interlayer distance of c/3 = 6.3 Å the 2D limit is reached in ErBr₃.

**Spin waves.** To calculate the dispersion of the magnetic excitations we used the MF-RPA method that treats the single-ion interactions (crystal-field and mean-field Hamiltonians) exactly while the interactions between the magnetic ions are considered in the random-phase approximation³⁰. The INS measurements of the crystal-field splitting of ErBr₃ are presented in Supplementary Note 2 and the analysis of the data explains the temperature dependence of the single-ion susceptibility shown in Supplementary Fig. 5. In ErBr₃, the first crystal-field level is located at 1.5 meV (see Supplementary Fig. 5), much higher than the energy scale of the dipolar interactions; the spin system can be approximated by a S = 1/2 doublet with effective anisotropic g-factors as detailed in Supplementary Note 3. For ErBr₃, we obtained a large easy-plane anisotropy with gₓ = gᵧ = 9.4 and gz = 0.1.

In MF-RPA the spin-wave excitations appear as poles in the dynamical susceptibility tensor $\bar{\bar{\chi}}(\mathbf{q}, E)$,

$$ \bar{\bar{\chi}}(\mathbf{q}, E) = [\bar{\bar{1}} - \bar{\bar{\chi}}_0(E)\bar{\bar{D}}(\mathbf{q})]^{-1} \bar{\bar{\chi}}_0(E), \quad (3) $$

where $\bar{\bar{\chi}}_0(E)$ is the single-ion susceptibility and $\bar{\bar{D}}(\mathbf{q})$ the Fourier transform of the dipolar interactions. Because the magnetic moments of ErBr₃ in the 2D phase form a non-collinear structure, it is convenient to introduce local coordinates for each spin in the magnetic cell with the quantization axis pointing along the spin direction, so that after transformation the single-ion susceptibility is the same for every magnetic site in the unit cell³⁰. The inelastic neutron cross-section $S(\mathbf{\kappa}, E)$ is then proportional to the imaginary part of the dynamical susceptibility

$$ S(\mathbf{\kappa}, E) = \frac{1}{\pi 1 - \exp(-E/k_B T)} \frac{1}{N} \frac{1}{\sum_{\alpha, \beta} \sum_{u,v}} \left( \delta_{\alpha, \beta} - \frac{\kappa_{\alpha} \kappa_{\beta}}{|\mathbf{\kappa}|^2} \right) \bar{\bar{\chi}}_{\alpha, \beta}(\mathbf{\kappa}, E) \quad (4) $$

with $\mathbf{\kappa} = \mathbf{\tau} + \mathbf{q}$ the scattering vector and $\mathbf{\tau}$ a reciprocal lattice vector; u, v number the N Er-ions in the magnetic cell and $\alpha, \beta = a, b, c$. Equation (4) is given in the crystal frame and hence, the $\bar{\bar{\chi}}_{\alpha, \beta}(\mathbf{\kappa}, E)$ matrices are evaluated in the rotated local coordinate system.

For the spin-wave calculations, only interactions within a single honeycomb layer were considered as the dipolar interactions between the planes are negligibly small. The spin-wave dispersion is shown in Supplementary Fig. 7. For scattering vectors along (0, q, 0) an acoustic mode and three optic branches are well observable, while the other two spin-wave branches have a considerably lower intensity. Along this direction in reciprocal space $S(\mathbf{\kappa}, E)$ does not change appreciably when Ψ is varied for the optical branches above $E \approx 0.1$ meV and no crossing of the acoustic mode with an optical branch occurs.

The INS measurements in ErBr₃ were performed with an energy resolution of about 70 µeV and the spin waves that have an energy $E < 0.1$ meV are hidden by the incoherent scattering at the elastic position. Fig. 2a shows representative INS measurements along (0, q, 0) at T = 80 mK together with the result of MF-RPA. Additional INS measurements are shown in Supplementary Fig. 6 together with a plot of the magnetic Brillouin zone where the scan directions are indicated. Only an overall scale factor is used to adjust the INS intensity with the MF-RPA calculations. The peaks in the INS spectrum at ~0.2 meV correspond to the dispersion of the highest visible magnon branch. Both the dispersion of the spin waves and the intensity of the INS are well reproduced. In Fig. 2b, we compare the dispersion of the spin waves in ErBr₃ with the MF-RPA calculations. It can be observed that the calculated spectrum of $S(\mathbf{\kappa}, E)$ is in agreement with the INS data.

## Discussion

The electronic band structure of materials with a honeycomb structure, like graphene, exhibits Dirac points at the corners of the Brillouin zone (the K-points in Supplementary Fig. 6c) where the electronic bands cross. Close to the K-point the electronic
<figure><img src="image_4.png" /><figcaption>**Fig. 2** Spin-wave measurements in ErBr<sub>3</sub>. a Representative inelastic neutron scattering measurements of the spin-wave excitations in ErBr<sub>3</sub> compared with the result of mean-field random-phase approximation (blue lines). The red line is a fit of the incoherent scattering measured at momentum transfer **k** = (2/3, -4/3, 0) and at *T* = 6 K, which is well above the temperature of the phase transition *T* = 280 mK and in the paramagnetic phase where there are no propagating coherent spin waves. The sum of both contributions is shown by the black line. Error bars are standard deviations given by the square root of the neutron counts. b False color plot of the inelastic neutron cross-section *S*(*k*, *E*) along **k** = (0, q, 0). In the calculations, we used a linewidth Γ = 0.07 meV in the single-ion susceptibility to approximate the energy resolution of the spectrometer (see Supplementary Note 3). The points represent the spin-wave dispersion in ErBr<sub>3</sub> obtained from a fit with a single Gaussian. Magnetic excitations below the energy *E* = 0.1 meV were not accessible in the experiment. The error bars are a result of the Gaussian fit.</figcaption></figure>

bands have a linear dispersion that corresponds to the dispersion of a relativistic and massless particle³¹. It was shown that the spin waves of a Heisenberg ferromagnet with a honeycomb structure have Dirac cones³²,³³ around the K-point opening the possibility of producing magnetic Dirac materials. In agreement with these predictions topological magnons have been observed as Dirac massless bosons in the ferromagnet CrBr₃³. The acoustic and optical spin-wave branches cross at the K-point of the Brillouin zone boundary and have a linear dispersion analogous to the electronic band structure of graphene. Magnon Dirac cones are not only found in ferromagnetic insulators and topological magnons protected by symmetry have been observed in the 3D antiferromagnets Cu₃TeO₆³⁵⁻³⁷ and CoTiO₃³⁸. In contrast, a magnon band gap is observed at the zone boundary in CrI₃. The origin of the gap is still debated and has been associated either with the DM interaction between next-nearest neighbors on the honeycomb lattice³⁹ or magnon-phonon coupling⁴⁰. In magnetic insulators, the DM interaction plays the role of the spin-orbit interaction present in graphene and leads to avoided band crossing of magnons. According to the bulk-edge correspondence⁴¹, magnons, in the gapped phase, can propagate uni-directionally along the surface or edges which opens the possibility of producing magnetic counterparts of topological insulators.

Figure 3a shows the dispersion of the spin waves for ErBr<sub>3</sub> along the (q, q, 0) direction. As the magnetic ground state is degenerate, we have performed the calculations for different values of Ψ for the magnetic moments of the Er-sublattices. The dispersion of the spin waves for the dipolar honeycomb lattice shows a Dirac cone at the K-point when the spins in ErBr<sub>3</sub> form a periodic vortex structure with Ψ = 0 (Fig. 3a). When Ψ is varied an anti-crossing of the lowest optical branch with the acoustic branch appears at the K-point of the magnetic Brillouin zone. The size of the energy gap between the two branches depends on Ψ as shown in Fig. 3b and reaches a maximum for Ψ = π/6. The splitting between the acoustic and optic branch decreases again for values of Ψ larger than π/6; finally, the two branches form a Dirac cone again when the arrangement of the Er moments corresponds to a periodic vortex structure with Ψ = π/3 (Fig. 3a).

Furthermore, the variation of Ψ from 0 to 2π produces six distinct spin structures on the honeycomb net corresponding to three vortex positions with opposite chiralities where the Dirac cones are present in the spin-wave spectrum at the K and K' points.

For Ψ = n × π/3 (n an integer), the acoustic and the lowest optical spin-wave branch cross at the K-point with a linear slope, and a Dirac cone is found at an energy of E ≈ 0.083 meV, as shown in Fig. 3c. We calculated a Berry phase γ = ±π (see Fig. 4a and Supplementary Note 4) that reflects the degeneracy of the two magnon branches at the Dirac points and shows that the topology of the magnon dispersion in ErBr₃ mimics the topology of the electronic band structure in graphene. Our results demonstrate that Dirac cones are not only present in honeycomb magnets with Heisenberg interactions, but also in van der Waals crystals like ErBr₃, where long-range dipolar forces are the primary magnetic interaction.

As in Cu₃TeO₆, Dirac points are present in the spectrum of magnetic excitations due to the invariance of the magnetic structure under combined inversion and time-reversal symmetry operations³⁵⁻³⁷. When Ψ ≠ n × π/3 (n an integer), the symmetry of the spin structure is broken which generates a finite Berry curvature as shown in Fig. 4b. In Fig. 4c and Supplementary Fig. 8, we show the evolution of the Berry curvature calculated for the magnon acoustic branch which is well separated from the other branches. The Berry curvatures at the K and K' points have the same magnitude but opposite signs as is the case in graphene⁴². For degeneracy angles close to Ψ = 0, the Berry curvature Ω(**k**) is peaked around the K and K' points and diverges when Ψ → 0. With increasing values of Ψ we observe that Ω(**k**) continuously broadens and weakens but remains centered around the K and K' points. The Berry curvature finally disappears at Ψ = π/6 when the gap between the acoustic and optical branch has its maximum value. Upon increasing the angle Ψ further, the Berry curvature reappears close to the K and K' points albeit with opposite signs. The Berry curvature of the dipolar honeycomb magnet hence depends on the size of the spin gap and can be manipulated by varying the angle Ψ. The Chern number *C* is equal to 0 when the Berry curvature is integrated over the entire Brillouin zone. Although *C* = 0 indicates a topologically trivial
<figure><img src="image_5.png" /><figcaption>**Fig. 3** Magnon gap in the dipolar honeycomb lattice. **a** Simulation of the inelastic neutron cross-section S(κ, E) along the momentum transfer κ = (q, q, 0) for selected Ψ. Note the opening of a gap at κ = (2/3, 2/3, 0) between the acoustic and the lowest optical mode at the energy E ≈ 0.083 meV when the spin configuration given by angle Ψ is varied. The color bar represents the calculated intensity of the inelastic neutron cross-section. **b** Dependence of the magnon band gap as a function Ψ at κ = (2/3, 2/3, 0). **c** Calculation of the spin-wave dispersion close to the Dirac point for Ψ = -5.6 × 10⁻⁴ π, 0, and +5.6 × 10⁻⁴ π, respectively.</figcaption></figure>
<figure><img src="image_6.png" /><figcaption>**Fig. 4** Berry phase and Berry curvature in the dipolar honeycomb lattice.
a The Berry phase associated with the Dirac cone ($\Psi = 0^\circ$) is either $+\pi$ (represented in yellow) or $-\pi$ (blue). The Berry curvature $\Omega(k)$ in ErBr$_3$ calculated for the acoustic magnon branch: **b** $\Psi = -5^\circ$, **c** $\Psi = 5^\circ$, $\Psi = 30^\circ$ and $\Psi = 55^\circ$. $\Psi$ is the angle representing the spin configuration. The Brillouin zone boundary is outlined with dashed lines. The color bar shows the magnitude of the Berry curvature.</figcaption></figure>

magnon insulator phase, when the integration of $\Omega(k)$ is restricted around a single Dirac point we numerically obtain for $\Psi = 1.5^\circ C_V \approx \pm 1/2$ at the K and K' points that is a prerequisite for the quantum valley Hall effect$^{43,44}$. A non-trivial Berry curvature is also at the origin of the magnon valley transports that were recently reported in refs.$^{45,46}$. For the dipolar honeycomb lattice, a topological phase transition occurs at $\Psi = 0$ where the spin gap closes at the Dirac points and $C_V$ changes sign. As the Berry curvature is antisymmetric with respect to the transition point, a topological magnon can propagate along the interface of two spin-domains located on opposite sides of the phase boundary$^{47}$. It is worth restating that these considerations hold for purely dipolar interactions. If exchange interactions are present, long-range ordering will occur with the moments normal to the plane$^{48}$.

In conclusion, we have shown that ErBr$_3$ serves as a model crystal for the dipolar antiferromagnet on the honeycomb lattice. The classical magnetic ground state forms a degenerate non-collinear vortex structure and we have found that the lowest magnon bands evolve from Dirac magnon excitations to a topological phase with non-trivial Berry curvature by varying the phase $\Psi$ that defines the spin configuration. The application of an external magnetic field lifts the degeneracy of the ground state and can be used to select a particular spin configuration by varying the direction of the magnetic field$^{22}$. As van der Waals

magnets can be cleaved this could be relevant for future investigations of magnon currents in honeycomb lattices. Our results suggest that if artificial dipolar honeycomb lattices can be built, it would have the potential to be extremely useful for magnonic devices since the energy at which the magnon band crossing occurs can be adapted by varying the size of the magnetic moment or the lattice spacing that allows us to engineer the magnon band structure$^{49,50}$.

**Methods**

**Single crystal of ErBr₃ and characterization.** The single crystal of ErBr₃ was synthesized according to the method described previously in ref. 25. The magnetic susceptibility was measured with an MPMS-5XL SQUID system (Quantum Design).

**Neutron experimental setup.** The CEF splitting was determined with the thermal-neutron spectrometer Eiger at the SINQ spallation source. The spectrometer was operated in the constant final energy mode with $k_f = 2.662$ Å⁻¹. To maximize the intensity the monochromator was double focused and the analyzer was horizontally focused. With this configuration the energy resolution is about $\Delta E = 0.8$ meV at the elastic position. Contamination by higher-order neutron wavelengths was eliminated by a PG002-filter installed in the scattered beam. The dispersion of the magnetic excitations was measured with the cold-neutron three-axis spectrometer IN14 at the Institut Laue-Langevin. The spectrometer was operated in a similar configuration as for the previous measurements albeit with $k_f = 1.15$ Å⁻¹ that resulted in an improved energy resolution of 70 µeV. A cold Be-filter, which scatters off neutrons with wavevectors longer than $k_f = 1.55$ Å⁻¹, was used to reduce the background and avoid spurious scattering. For these measurements, a single crystal of approximately 0.5 cm³ was mounted inside a dilution refrigerator and cooled down to the base temperature of $T = 80$ mK that is well below the ordering temperature of the Er³⁺ moments at $T = 280$ mK.

**Data availability**

All data needed to evaluate the conclusions in the paper are present in the paper and/or the Supplementary information. Additional data related to this paper may be requested from the authors.

Received: 18 February 2022; Accepted: 1 July 2022;

Published online: 16 July 2022

**References**

1. Parkin, S., Jiang, X., Kaiser, C., Panchula, A. & Samant, M. Magnetically engineered spintronic sensors and memory. *Proc. IEEE* **91**, 661 (2003).

2. Hirohata, A. et al. Review on spintronics: principles and device applications. *J. Magn. Magn. Mater.* **509**, 166711 (2020).

3. Kruglyak, V. V., Demokritov, S. O. & Grundler, D. Magnonics. *J. Phys. D: Appl. Phys.* **43**, 264001 (2010).

4. Chumak, A. V., Vasyuchka, V. I., Serga, A. A. & Hillebrands, B. Magnon spintronics. *Nat. Phys.* **11**, 453 (2015).

5. Geim, A. K. & Grigorieva, I. V. Van der Waals heterostructures. *Nature* **499**, 419 (2013).

6. Gibertini, M., Koperski, M., Morpurgo, A. F. & Novoselov, K. S. Magnetic 2D materials and heterostructures. *Nat. Nanotechnol.* **14**, 408 (2019).

7. Xing, W. et al. Magnon transport in quasi two-dimensional van der Waals antiferromagnets. *Phys. Rev. X* **9**, 011026 (2019).

8. Zhang, Y. et al. MnPS₃ spin-flop transition-induced anomalous Hall effect in graphite flake via van der Waals proximity coupling. *Nanoscale* **12**, 23266 (2020).

9. Wang, X. S., Zhang, H. W. & Wang, X. R. Topological magnonics: a paradigm for spin-wave manipulation and device design. *Phys. Rev. Appl.* **9**, 024029 (2018).

10. Mook, A., Plekhanov, K., Klinovaja, J. & Loss, D. Interaction-stabilized topological magnon insulator in ferromagnets. *Phys. Rev. X* **11**, 021061 (2021).

11. Mook, A., Henk, J. & Mertig, I. Edge states in topological magnon insulators. *Phys. Rev. B* **90**, 024412 (2014).

12. Moriya, T. Anisotropic superexchange interaction and weak ferromagnetism. *Phys. Rev.* **120**, 91 (1960).

13. Onose, Y. et al. Observation of the magnon Hall effect. *Science* **329**, 297 (2010).

14. Zhang, L., Ren, J., Wang, J. S. & Li, B. Topological magnon insulator in insulating ferromagnet. *Phys. Rev. B* **87**, 144101 (2013).
15. Shindou, R., Ohe, J., Matsumoto, R., Murakami, S. & Saitoh, E. Chiral spin-wave edge modes in dipolar magnetic thin films. *Phys. Rev. B* **87**, 174402 (2013).

16. Shen, K. Magnon spin relaxation and spin Hall effect due to the dipolar interaction in antiferromagnetic insulators. *Phys. Rev. Lett.* **124**, 077201 (2020).

17. Liu, J., Wang, L. & Shen, K. Spin-orbit coupling and linear crossings of dipolar magnons in van der Waals antiferromagnets. *Phys. Rev. B* **102**, 144416 (2020).

18. Liu, J., Wang, L. & Shen, K. Dipolar spin waves in uniaxial easy-axis antiferromagnets: a natural topological nodal-line semimetal. *Phys. Rev. Res.* **2**, 023282 (2020).

19. Luttinger, J. M. & Tisza, L. Theory of dipole interaction in crystals. *Phys. Rev.* **70**, 954 (1946).

20. Mermin, N. D. & Wagner, H. Absence of ferromagnetism or antiferromagnetism in one- or two-dimensional isotropic Heisenberg models. *Phys. Rev. Lett.* **17**, 1133 (1966).

21. Malozovsky, Y. M. & Rozenbaum, V. M. Orientational ordering in two-dimensional systems with long-range interaction. *Phys. A* **175**, 127 (1991).

22. Zimmerman, G. O., Ibrahim, A. K. & Wu, F. Y. Planar classical dipolar system on a honeycomb lattice. *Phys. Rev. B* **37**, 2059 (1988).

23. Rozenbaum, V. M. Ground state and vibrations of dipoles on a honeycomb lattice. *Phys. Rev. B* **51**, 1290 (1995).

24. Maksymenko, M., Chandra, V. R. & Moessner, R. Classical dipoles on the kagome lattice. *Phys. Rev. B* **91**, 184407 (2015).

25. Krämer, K. W. et al. Noncollinear two- and three-dimensional magnetic ordering in the honeycomb lattices of $ErX_3$ ($X=Cl, Br, I$). *Phys. Rev. B* **60**, R3724 (1999).

26. Braekken, H. Die Kristallstruktur der Trijodide von Arsen, Antimon und Wismut. *Z. Krist.* **74**, 67 (1930).

27. Bertaut, E. F. In *Magnetism*, vol. III (eds Rado Suhl, G. T. & Suhl, H.) 150, (Academic Press, 1963).

28. Rastelli, E., Carbognani, A., Regina, S. & Tassi, A. Order by thermal disorder in 2D planar rotator model with dipolar interactions. *Eur. Phys. J. B* **9**, 641 (1999).

29. Enjalran, M. & Gingras, M. J. P. Theory of paramagnetic scattering in highly frustrated magnets with long-range dipole-dipole interactions: the case of the $Tb_2Ti_2O_7$ pyrochlore antiferromagnet. *Phys. Rev. B* **70**, 174426 (2004).

30. Jensen, J. & Macintosh, A. R. *Rare Earth Magnetism* (Clarendon Press, Oxford, 1991).

31. Wehling, T. O., Black-Schaffer, A. M. & Balatsky, A. V. Dirac materials. *Adv. Phys.* **63**, 1 (2014).

32. Fransson, J., Black-Schaffer, A. M. & Balatsky, A. V. Magnon Dirac materials. *Phys. Rev. B* **94**, 075401 (2016).

33. Boyko, D., Balatsky, A. V. & Haraldsen, J. T. Evolution of magnetic Dirac bosons in a honeycomb lattice. *Phys. Rev. B* **97**, 014433 (2018).

34. Pershoguba, S. S. et al. Dirac magnons in honeycomb ferromagnets. *Phys. Rev. X* **8**, 011010 (2018).

35. Li, K. et al. Dirac and nodal line magnons in three-dimensional antiferromagnets. *Phys. Rev. Lett.* **119**, 247202 (2017).

36. Bao, S. et al. Discovery of coexisting Dirac and triply degenerate magnons in a three-dimensional antiferromagnet. *Nat. Commun.* **9**, 2591 (2018).

37. Yao, W. et al. Topological spin excitations in a three-dimensional antiferromagnet. *Nat. Phys.* **14**, 1011 (2018).

38. Yuan, B. et al. Dirac magnons in a honeycomb lattice quantum XY magnet $CoTiO_3$. *Phys. Rev. X* **10**, 011062 (2020).

39. Chen, L. et al. Topological spin excitations in honeycomb ferromagnet $CrI_3$. *Phys. Rev. X* **8**, 041028 (2018).

40. Delugas, P. et al. Magnon-phonon interactions open a gap at the Dirac point in the spin-wave spectra of $CrI_3$ 2D magnets. Preprint at arXiv:2105.04531 (2021).

41. Hasan, M. Z. & Kane, C. L. Colloquium: topological insulators. *Rev. Mod. Phys.* **82**, 3045 (2010).

42. Xiao, D., Yao, W. & Niu, Q. Valley-contrasting physics in graphene: magnetic moment and topological transport. *Phys. Rev. Lett.* **99**, 236809 (2007).

43. Yao, W., Yang, S. A. & Niu, Q. Edge states in graphene: from gapped flat-band to gapless chiral modes. *Phys. Rev. Lett.* **102**, 096801 (2009).

44. Qian, K., Apigo, D. J., Prodan, C., Barlas, Y. & Prodan, E. Topology of the valley-Chern effect. *Phys. Rev. B* **98**, 155138 (2018).

45. Ghader, D. Valley-polarized domain wall magnons in 2D ferromagnetic bilayers. *Sci. Rep.* **10**, 16733 (2020).

46. Zhai, X. & Blanter, Y. M. Topological valley transport of gapped Dirac magnons in bilayer ferromagnetic insulators. *Phys. Rev. B* **102**, 075407 (2020).

47. Mook, A., Henk, J. & Mertig, I. Topologically nontrivial magnons at an interface of two kagome ferromagnets. *Phys. Rev. B* **91**, 224411 (2015).

48. Pich, C. & Schwabl, F. Spin-wave dynamics of two-dimensional isotropic dipolar honeycomb antiferromagnets. *J. Magn. Magn. Mater.* **148**, 30–31 (1995).

49. Streubel, R. et al. Spatial and temporal correlations of XY macro spins. *Nano. Lett.* **18**, 7428–7434 (2018).

50. Lenk, B., Ulrichs, H., Garbs, F. & Münzenberg, M. The building blocks of magnonics. *Phys. Rep.* **507**, 107 (2011).

## Acknowledgements

We would like to thank C. Mudry for the useful discussions. The financial support by the Swiss National Science Foundation under grant no. SNF 200020_172659 is gratefully acknowledged.

## Author contributions

The single crystal was grown by K.W.K. The inelastic neutron scattering experiments were performed by H.B.B., K.W.K., B.R., U.S., A.W., and C.W. The susceptibility measurements were done by K.W.K. Modeling and simulations were performed by B.R. and C.W. in discussion with H.B.B., K.W.K., and M.K. All authors contributed to the writing of the manuscript.

## Competing interests

The authors declare no competing interests.

## Additional information

**Supplementary information** The online version contains supplementary material available at https://doi.org/10.1038/s42005-022-00965-5.

**Correspondence** and requests for materials should be addressed to Christian Wessler.

**Peer review information** *Communications Physics* thanks the anonymous reviewers for their contribution to the peer review of this work.

**Reprints and permission information** is available at http://www.nature.com/reprints

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

<figure><img src="image_7.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2022