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

# Magnon gap excitations in van der Waals antiferromagnet MnPSe$_3$

Dipankar Jana<fnref n="1" />,², D. Vaclavkova<fnref n="1" />, I. Mohelsky<fnref n="1" />, P. Kapuscinski<fnref n="1" />,³, C. W. Cho<fnref n="1" />,⁷, I. Breslavetz<fnref n="1" />, M. Białek⁴,⁵, J.-Ph. Ansermet⁵, B. A. Piot<fnref n="1" />, M. Orlita<fnref n="1" />,⁶, C. Faugeras<fnref n="1" /> & M. Potemski<fnref n="1" />

Magneto-spectroscopy methods have been employed to study the zero-wavevector magnon excitations in MnPSe$_3$. Experiments carried out as a function of temperature and the applied magnetic field show that two low-energy magnon branches of MnPSe$_3$ in its antiferromagnetic phase are gapped. The observation of two low-energy magnon gaps (at 1.70 ± 0.05 meV and 0.09 ± 0.01 meV) implies that MnPSe$_3$ is a biaxial antiferromagnet. A relatively strong out-of-plane anisotropy imposes the spin alignment to be in-plane whereas the spin directionality within the plane is governed by a factor of 2.5 × 10⁻³ weaker in-plane anisotropy.

Two-dimensional (2D) materials beyond graphene constitute an exceptional platform for the exploration of the unique electronic, optical, and magnetic properties of low-dimensional systems<fnref n="1" />⁻⁵. Among currently studied materials are layered, van der Waals magnets, which include a family of metal phosphorus trichalcogenides (MPX$_3$, with M = Mn, Fe, Co, Ni and X = S, Se)⁶⁻¹¹. MPX$_3$ compounds are semiconducting antiferromagnets with diverse spin ordering configurations of the transition metal ions¹². Since MPX$_3$ layers are held together by weak van der Waals forces, typical distances between the metal ions within the layers are shorter than those between ions in different layers, and the magnetic order in MPX$_3$ compounds is governed mainly by spin-spin coupling within the layers. Such quasi-2D or even strictly 2D (in monolayers) magnetic ordering persists in MPX$_3$ materials since they are not necessarily perfect Heisenberg antiferromagnets, but, in addition to the exchange terms, are also characterized by non-zero magnetic anisotropies. A non-negligible spin-orbit interaction within the M-ions is a primary cause of the so-called single-ion magnetic anisotropy. At first sight, the spin-orbit interaction is expected to be rather weak for Mn²⁺ ions in MnPX$_3$, though the nature seems to be against the Mermin–Vagner theorem¹³ (no magnetic order in the 2D limit for Heisenberg antiferromagnets) and MnPX$_3$ compounds do display antiferromagnetic ordering¹⁴,¹⁵. Interestingly enough, the change of the chalcogen atom from sulfur (S) to selenium (Se) imposes the change in the type of magnetic ordering¹⁶,¹⁷. A single-ion anisotropy may still be relevant for ordering spins in MnPSe$_3$¹⁵,¹⁶, whereas the competing anisotropy due to magnetic dipolar interaction is likely imposing the character of spin-ordering in MnPS$_3$¹⁴. From another perspective, the single ion anisotropy in MnPSe$_3$ is considered to be small while interlayer exchange integrals are significantly large¹⁸. Despite these ambiguities, it is known that spins are aligned in the direction nearly perpendicular to the layers' plane in the case of broadly studied MnPS$_3$¹⁹. Instead, the less studied MnPSe$_3$ is an easy-plane antiferromagnet, as deduced from neutron scattering¹⁸,²⁰ and magnetization¹⁵,¹⁶ measurements as well as theoretical calculations⁷. The open questions addressed in the present work are in regards to the character of low energy magnon excitations in MnPSe$_3$. One of the calculated magnon branches⁷ has been concluded to be gapless (being dispersed linearly with the *k*-vector) and the *k* = 0 gap of the second magnon branch is estimated to be ~14 cm⁻<fnref n="1" />. Neutron scattering data¹⁸ confirm the upper 14 cm⁻<fnref n="1" /> gap but lack the resolution to trace the low energy magnon branch in the *k* = 0 limit. On the other hand, the magnetization measurements¹⁶ indicate a spin-flop field at ~0.7 T, pointing out a non-zero onset (at *k* = 0) of the lower magnon branch.

In the present work, we demonstrate that magnon gap excitations in MnPSe$_3$ antiferromagnet can be effec-
tively probed using magneto-spectroscopy methods (Raman scattering, far-infrared, and microwave absorption
techniques). While neutron diffraction is a common experimental tool for investigating magnon dispersion
in antiferromagnetic materials¹⁰,¹⁸, the magneto-optical techniques employed in this work represent useful

<footnote n="1">¹Laboratoire National des Champs Magnétiques Intenses, LNCMI-EMFL, CNRS UPR3228, INSA-T, Univ. Grenoble Alpes, Univ. Toulouse, Univ. Toulouse 3, Grenoble and Toulouse, France. ²Institute for Functional Intelligent Materials, National University of Singapore, Singapore 117544, Singapore. ³Institute of Experimental Physics, Faculty of Physics, University of Warsaw, ul. Pasteura 5, 02-093 Warsaw, Poland. ⁴CENTERA Labs, Institute of High Pressure Physics, PAS; CEZAMAT, Warsaw University of Technology, Warsaw, Poland. ⁵Institute of Physics, Ecole Polytechnique Fédérale de Lausanne (EPFL), 1015 Lausanne, Switzerland. ⁶Institute of Physics, Charles University, 121 16 Prague, Czech Republic. ⁷Department of Physics, Chungnam National University, Daejeon 341134, Republic of Korea. ⁸email: marek.potemski@lncmi.cnrs.fr</footnote>
alternatives that offer additional insights<sup>8,9,21,22</sup>. Mostly, these methods address the *k* = 0 excitations only, but they allow one to explore the magnon modes in a broader energy range and with a greater energy resolution. With these experiments, we identify MnPSe<sub>3</sub> as a *de facto* biaxial antiferromagnet, by revealing the energies of two fundamental magnon gaps, both of which are non-zero (14 and 0.7 cm<sup>-1</sup>). This implies that the antiferromagnetic order in MnPSe<sub>3</sub> is governed by two anisotropy terms: one previously identified and directed across the layers' plane (*D*<sub>⊥</sub> ≈ 0.3 cm<sup>-1</sup>)<sup>16</sup>, which imposes the in-plane spin alignment, and another one, weak but non-negligible *D*<sub>‖</sub> term, estimated here to be ≈ 7.5 × 10<sup>-4</sup> cm<sup>-1</sup>, which imposes the spins' alignment along a certain, yet unidentified, direction within the plane.

## Results

The magnetic structure of MnPSe<sub>3</sub> in the antiferromagnetic phase is shown in Fig. 1a. It is assumed<sup>15</sup> that the in-plane alignment of the colinear Mn<sup>2+</sup> magnetic moments is imposed by the (single ion) anisotropy term (*D*<sub>⊥</sub> ≈ 0.3 cm<sup>-1</sup>)<sup>16</sup> along the c-axis. The axis of spin alignment within the a-b plane is chosen arbitrarily. So far, no in-plane anisotropy term was theoretically identified in MnPSe<sub>3</sub>. The net antiferromagnetic order in the planes is governed by the nearest neighbor exchange interactions. The spin alignment across the layers is reported to be antiferromagnetic<sup>18</sup>. Raman scattering is one of the experimental techniques used to probe the low energy excitations in solids and in particular the *k* = 0 spin excitations (magnon gap excitations) in layered antiferromagnets<sup>24</sup>. Such a method revealed the coupled phonon-two magnon excitations in MnPSe<sub>3</sub> antiferromagnet<sup>7</sup> which appear in the spectral range relatively far from the laser line. Instead, in our experiments, we focus on the Raman scattering response of MnPSe<sub>3</sub> in an energy range close to the laser line. The temperature evolution of Raman scattering spectra of MnPSe<sub>3</sub> measured in the range of ± 20 cm<sup>-1</sup> is shown in Fig. 1b. The Raman scattering feature at an energy of *w*<sub>M+</sub> ≈ 14 cm<sup>-1</sup> in the limit of low temperatures, denoted as M+, characteristically softens with rising temperature, is identified as an upper magnon gap excitation of MnPSe<sub>3</sub>. Notably, the *w*<sub>M+</sub> energy coincides well with the previously reported value<sup>7</sup>. The optional approach to trace the temperature evolution of the M+ mode, with the technique of THz monochromatic transmission spectroscopy in the spectral range 270–390 GHz (9–13 cm<sup>-1</sup>), is presented in section S1 of the SM<sup>25</sup>.

Next, we focus on experiments carried out at low temperature, but as a function of the magnetic field applied either in the direction perpendicular or parallel to the layers' planes (out-of-plane and in-plane configurations, correspondingly). The evolution of the M+ mode with the magnetic field applied in the out-of-plane configuration could be conveniently traced with far-infrared transmission measurements. The results of these experiments are shown in Fig. 2a. We observe that the wM+ increases upon applying the magnetic field. One may anticipate that wM+ initially displays a quadratic B dependence, which turns into a linear one at high magnetic fields. Our approach to trace the M+ by Raman scattering in the out-of-plane configuration was not successful due to the specific arrangement of our experimental set-up, preventing an efficient rejection of the stray laser light. Nevertheless, we performed such experiments by setting an angle of 20° between the direction of the magnetic field and the normal to the surface. The results of these experiments, as shown in section S2 of the SM<sup>25</sup>, are consistent with those presented in Fig. 2a. The performance of our magneto-Raman scattering set-up is superior in the configuration of the in-plane magnetic field. A series of low-temperature (5 K) Raman scattering spectra measured as a function of the in-plane magnetic field (arbitrary direction within the layers' plane) is shown in Fig. 2b. The energy of the M+ mode is found to be unaffected by the in-plane magnetic field. On the other hand, another low energy M− mode emerges in the spectra from the region covered by the stray laser light, when the

<figure><img src="image_2.png" /><figcaption>b</figcaption></figure>

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

**Figure 1.** (a) The magnetic structure of MnPSe<sub>3</sub> in the antiferromagnetic phase. The grey sphere and red arrow represent Mn<sup>2+</sup> ion and spin direction respectively. The figure is created using the VESTA software package<sup>23</sup>. (b) Raman spectra of MnPSe<sub>3</sub> measured at some selected temperatures. M+ resonance corresponds to the upper magnon gap excitation. Spectra are shifted vertically for clarity. Black triangles mark the peak energy wM+ of M+ resonance. The shaded region indicates the spectral range blocked by filtering of the laser stray light. As illustrated in the inset, the wM+ dependence on temperature follows a conventional formula for the antiferromagnets (see the text) as shown with the solid line.
<figure><img src="image_4.png" /><figcaption>**Figure 2.** (a) False color map of normalized far-infrared transmission of MnPSe₃, measured at 4.2 K, as a function of the magnetic field applied along the crystal c-axis. A few representative spectra, measured at a magnetic field strength of 3 T, 6 T, 9 T, 12 T, and 15 T are also plotted. The dashed line corresponds to a fit of w<sub>M+</sub>(B) dependence according to Eq. 1. (b) Raman scattering spectra of MnPSe₃ measured at selected in-plane magnetic fields. M<sub>+</sub> and M<sub>-</sub> resonances are identified as, correspondingly the upper and lower magnon gap excitations. Spectra are shifted vertically for clarity. Black triangles mark the peak energies which are plotted as a function of the magnetic field in the inset. The solid line corresponds to a linear fit to the field dependence of M<sub>-</sub> and M<sub>+</sub> modes while the dotted line corresponds to extrapolation of the fit to 0 T. The shaded region indicates the spectral range blocked by filtering of the laser stray light.</figcaption></figure>

strength of the applied field exceeds *B* = 6 *T*. In this range of the magnetic field, the energy of the M− mode is observed to increase linearly with *B*. Interestingly, at fields around *B* = 14 *T*, the M− mode crosses the M+ mode and *w*<sub>M−</sub> > *w*<sub>M+</sub> in the range of high magnetic fields (see section S3 of the SM²⁵). Of particular interest for defining the character of spin ordering in MnPSe<sub>3</sub> is the inspection of the *M*− versus *B* dependence in the limit of low magnetic fields. To access this limit, the phase-sensitive, resistivity-detected microwave magneto-absorption measurements²² have been carried out. In this experiment, the resonance absorption signals have been detected when exposing a MnPSe<sub>3</sub> specimen to GHz radiation (selected frequencies in the range of 10 GHz to 64 GHz) upon sweeping the magnetic field applied in the arbitrarily chosen direction but remaining along the layers’ plane. As can be seen in Fig. 3a,b, the relatively sharp resonances are observed in the high field range where they disperse almost linearly with *B*. They are recognized as M− magnon modes, indeed merging those observed in Raman scattering experiments at still higher magnetic fields (see Fig. 2b). Importantly, however, the measured GHz resonances significantly broaden, indicating the flattening of the *E*<sub>M−</sub> versus *B* dependence when lowering the magnetic field, and no resonance is observed below 20 GHz. This points towards a small but clearly non-zero (≈ 0.7 cm<sup>−1</sup>) onset for the M− magnon gap at *B* = 0 *T*.

<figure><img src="image_5.png" /><figcaption>**Figure 3.** (a) False color map of GHz absorption phase of MnPSe<sub>3</sub>, measured at 5 K, as a function of the in-plane magnetic field. A few spectra at selected frequencies are also displayed in the same plot showing a change in the absorption phase at the resonance energy. (b) False color map of the simulated (normalized) intensity of M− mode resonance absorption considering the three symmetry orientations (18°, 78°, and 42°) of the spin to the in-plane magnetic field. Peak energies obtained from phase spectra (Fig. 3a) are shown by cyan squares whereas the half-widths of the observed resonances are indicated with horizontal bars.</figcaption></figure>
Discussion

The results presented above are fully consistent with the identification of the observed M+ and M− modes as two low energy magnon gap excitations in MnPSe3 antiferromagnet. The temperature dependence of the energy of the M+ mode presented in Fig. 1b is typical for a large class of antiferromagnets²⁶. The solid line in the inset of Fig. 1b represents the expected temperature softening of the magnon gap given by a conventional formula: $w_{M+}(T) = w_{M+}(1 - T/T_N)^{1/3}$, in which we have set $w_{M+} = 14 \text{ cm}^{-1}$ for the magnon gap in the low temperature limit and $T_N = 74$ K for the Neél temperaturure in MnPSe3²⁰.

Further evidence for the identification of M+ and M− modes comes out from their characteristic dependencies upon the application of the magnetic field. These dependencies are here discussed in reference to a classical theory of the antiferromagnetic resonance in a biaxial antiferromagnet²⁷. First, we focus on the out-of-plane configuration that is the magnetic field applied perpendicularly to the layers’ planes ($B \parallel c$). In this case, all Mn²⁺ spins, regardless of their orientation in the plane, are aligned perpendicularly to B. The field-dependent energies of M+ and M− modes are given by the following expressions²⁷:

$$w_{M-}(B) = w_{M-} \qquad w_{M+}(B) = \sqrt{(w_{M+})^2 + (g\mu_B B)^2} \qquad (1)$$

where $w_{M-}$ and $w_{M+}$ denote the M− and M+ magnon energies at $B = 0$ T and g is the effective g-factor. We thus expect that the $w_{M-}$ remains unaffected by the out-of-plane magnetic field, while $w_{M+}$ increases with B, initially quadratically but linearly at high fields. The $w_{M+}(B)$ dependence given in the above formulas is in accordance with the results presented in Fig. 2a: The dashed line in Fig. 2a has been drawn when setting $w_{M+} = 14 \text{ cm}^{-1}$ and g = 2 for the effective g-factor. We also note that when replacing B by $B \cos(20°)$, the above formula describes well the Raman scattering data presented in the Fig. S2 of the SM²⁵.

The analysis of the data obtained in the in-plane magnetic field configuration is more complex. As illustrated in Fig. S4 of the SM, the apparent shapes of $w_{M-}(B)$ and $w_{M+}(B)$ dependencies are different for different orientations of the magnetic field with respect to the crystal axis along which the spins are initially (at $B = 0$ T) aligned²⁷. Unfortunately, this axis is not identified in our work, and moreover, the samples used in our experiments are thick and likely composed of grains with different orientations of the in-plane crystal axis. The Raman scattering data shown in Fig. 2b correspond to the high field range, above the spin-flop field $B_{sf}$. At sufficiently high fields (see section S4 of the SM²⁵), independent of the actual orientation of the in-plane magnetic field, the $w_{M+}(B)$ and $w_{M-}(B)$ energies, satisfy the following conditions:

$$\sqrt{(g\mu_B B)^2 - (w_{M-})^2} \le w_{M-}(B) \le \sqrt{(g\mu_B B)^2 + (w_{M-})^2}$$

$$\sqrt{(w_{M+})^2 - (w_{M-})^2} \le w_{M+}(B) \le w_{M+}$$

Expecting a rather small value for the spin-flop field ($B_{sf} = 0.7$ T)¹⁶ and thus a small value for the $w_{M-} = 0.7 \text{ cm}^{-1}$, one finds that in the high field range $w_{M+}$ is not affected by the applied field ($w_{M+}(B) = w_{M+}$) whereas, $w_{M-}(B)$ depends linearly on the applied field ($w_{M-}(B) = g\mu_B B$). This accounts for the observation, as shown in the inset to Fig. 2b, of linear with B-field dependence of $w_{M-}(B)$.

Less evident is the interpretation of the data presented in Fig. 3a, which results from our approach to trace the $w_{M-}(B)$ dependence in the range of low in-plane magnetic fields, on either side of the spin-flop field. As previously stated, access to this range can be achieved through an alternative experimental approach, resembling the multi-frequency magnetic resonance technique. This method necessitates the utilization of relatively thick specimens, unavoidably composed of domains with six-fold in-plane rotational symmetry of the crystal axis. The crystal symmetry confines the spin orientation in three directions, each rotated by 60° in a multi-domain sample. Additionally, the magnetic field is applied in an arbitrary in-plane direction relative to the magnetic easy axis. The $w_{M-}(B)$ dependence is different for different orientations of the crystal axis with respect to the magnetic field (see Fig. S4b of the SM). Consequently, the data presented in Fig. 3a can be effectively reproduced when “averaging” the $w_{M-}(B)$ dependencies expected for the three differently oriented domains. False color maps shown in Fig. S5 of the SM have been constructed by summing up the $w_{M-}(B)$ dependencies distinctive of each of the three symmetry domains for a specific magnetic field orientation relative to the spin orientation of one domain. Each of those three dependencies has been calculated assuming $w_{M-} = 0.7 \text{ cm}^{-1}$ for the energy of the M− mode at zero magnetic field and consistently with the previous analysis, g = 2 for the effective g-factor. Most of the traces shown in Fig. S5 concentrate in the upper energy limit given by $w_{M+}(B) = \sqrt{(w_{M+})^2 + (g\mu_B B)^2}$. The trace map which reproduces the experimental data well is presented in Fig. 3b. The resonance is broad in the limit of low magnetic fields but relatively narrow at higher fields, where independently of a specific configuration $w_{M-}(B)$ disperses linearly with the magnetic field and matches the data traced with magneto-Raman scattering (see Fig. 2b).

Having established two low energy magnon gaps in MnPSe3, we identify this antiferromagnet as a biaxial system with two anisotropy fields ($B_{D\perp}$ and $B_{D\parallel}$) or energies ($D_{\perp}$ and $D_{\parallel}$): $D_{\perp} = g\mu_B B_{D\perp}$ and $D_{\parallel} = g\mu_B B_{D\parallel}$. In agreement with previous reports¹⁶, the out-of-plane anisotropy $D_{\perp} = 0.3 \text{ cm}^{-1}$ restricts the in-plane spin alignment and accounts for the upper magnon gap: $w_{M+} = \sqrt{2E_{exch}D_{\perp}} = 14 \pm 0.4 \text{ cm}^{-1}$, where $E_{exch}$ is the mean exchange energy. On the other hand, the reported here, non-zero gap $w_{M-} = 0.7 \pm 0.1 \text{ cm}^{-1}$ for the lower magnon branch indicates the presence of the additional in-plane anisotropy $D_{\parallel}(w_{M-} = \sqrt{2E_{exch}D_{\parallel}})$. Although rather weak: $D_{\parallel} = D_{\perp} \times (w_{M-}/w_{M+})^2 = 7.5 (\pm 2.5) \times 10^{-4} \text{ cm}^{-1}$, $D_{\parallel}$ imposes spins to be aligned along certain in-plane crystal axis. This axis remains unidentified in the present work which calls for further, experimental and theoretical studies of MnPSe3 antiferromagnet.
## Concluding remarks

We have employed Raman scattering, far-infrared transmission, and microwave absorption measurements, carried as a function of temperature and of the magnetic field, to identify magnon gap excitations in MnPSe₃ van der Waals antiferromagnet. The observation of two distinct magnon gap excitations supports the identification of MnPSe₃ as a biaxial antiferromagnet with non-negligible in-plane anisotropy. Relevant parameters (fundamental magnon gap energies, g-factor, spin-flop field, anisotropy energies) are also estimated.

## Methods

We utilized commercially available thick and large area (≈ 25 mm² × 500 µm) MnPSe₃ crystals for FIR-transmission and Raman scattering measurements and relatively thin and small size crystals for GHz absorption measurements. Magneto-optical experiments were performed in both the Faraday and Voigt configurations with the magnetic field (from a superconducting solenoid magnet) applied perpendicular and parallel to the layers' plane respectively. In the case of Raman scattering measurements, the sample was excited by a semiconductor-based laser (λ = 515 nm) at a tilt of 20° with respect to the surface normal?. The laser beam was focused on the sample using a microscope objective with a numerical aperture of 0.83. The scattered signal was collected using the same objective and was analyzed by a 0.7m grating spectrometer equipped with a liquid nitrogen-cooled charge-coupled device (CCD). A set of reflection-based Bragg filters was employed in both the excitation and collection paths to effectively eliminate any unwanted laser line and to reject backscattered laser. For far infrared magneto-transmission measurements²¹, the radiation emitted from a mercury lamp was modulated using a Bruker Vertex 80v Fourier-transform spectrometer. The modulated radiation was delivered to the sample through light-pipe optics and subsequently detected by a composite bolometer placed directly behind the sample. Magneto-absorption experiments in the GHz frequency range were performed utilizing the modulation technique²². The radiation emitted by an Agilent gigahertz source was transmitted to the sample through a coaxial cable. The sample was positioned in the proximity of a resistive thermometer, which served as a bolometer. When the sample effectively absorbs the incident radiation, it results in a reduction of the radiation reaching the resistor. This causes a change in the resistor's temperature, which can be detected by monitoring the variation in its resistance. The alteration in resistivity is detected through changes in modulation response amplitude, while the delay in the resistivity response, owing to the finite thermalization time, is observed through the modulation response phase. The change in the modulation response phase with respect to the off-resonance phase is plotted in Fig. 3a.

## Data availability

The datasets generated during and/or analysed during the current study are available from the corresponding author on reasonable request.

Received: 28 July 2023; Accepted: 10 July 2024

Published online: 30 July 2024

## References

1. Bhimanapati, G. R. et al. Recent advances in two-dimensional materials beyond graphene. ACS Nano 9, 11509–11539. https://doi.org/10.1021/acsnano.5b05556 (2015).

2. Manzeli, S., Ovchinnikov, D., Pasquier, D., Yazyev, O. V. & Kis, A. 2D transition metal dichalcogenides. Nat. Rev. Mater. 2, 1–15. https://doi.org/10.1038/natrevmats.2017.33 (2017).

3. Mak, K. F. & Shan, J. Photonics and optoelectronics of 2D semiconductor transition metal dichalcogenides. Nat. Photon. 10, 216–226. https://doi.org/10.1038/nphoton.2015.282 (2016).

4. Geng, D. & Yang, H. Y. Recent advances in growth of novel 2D materials: Beyond graphene and transition metal dichalcogenides. Adv. Mater. 30, 1800865. https://doi.org/10.1002/adma.201800865 (2018).

5. Molle, A. et al. Buckled two-dimensional Xene sheets. Nat. Mater. 16, 163–169. https://doi.org/10.1038/nmat4802 (2017).

6. Xing, W. et al. Magnon transport in quasi-two-dimensional van der Waals antiferromagnets. Phys. Rev. X 9, 011026. https://doi.org/10.1103/PhysRevX.9.011026 (2019).

7. Mai, T. T. et al. Magnon-phonon hybridization in 2D antiferromagnet MnPSe₃. Sci. Adv. 7, eabj3106. https://doi.org/10.1126/sciadv.abj3106 (2021).

8. Liu, S. et al. Direct observation of magnon-phonon strong coupling in two-dimensional antiferromagnet at high magnetic fields. Phys. Rev. Lett. 127, 097401 (2021).

9. Cui, J. et al. Chirality selective magnon-phonon hybridization and magnon-induced chiral phonons in an atomically thin antiferromagnet. Nat. Comm. 14, 3396. https://doi.org/10.1038/s41467-023-39123-y (2023).

10. Wildes, A. R. et al. Spin wave spectra of single crystal CoPS₃. Phys. Rev. B 107, 054438. https://doi.org/10.1103/PhysRevB.107.054438 (2023).

11. Kang, S. et al. Coherent many-body exciton in van der Waals antiferromagnet NiPS₃. Nature 583, 785–789. https://doi.org/10.1038/s41586-020-2520-5 (2020).

12. Brec, R. Review on structural and chemical properties of transition metal phosphorous trisulfides MPX₃. Solid State Ion. 22, 3–30. https://doi.org/10.1016/0167-2738(86)90055-X (1986).

13. Mermin, N. D. & Wagner, H. Absence of ferromagnetism or antiferromagnetism in one-or two-dimensional isotropic Heisenberg models. Phys. Rev. Lett. 17, 1133. https://doi.org/10.1103/PhysRevLett.17.1133 (1966).

14. Joy, P. & Vasudevan, S. Magnetism in the layered transition-metal thiophosphates MPS₃ (M = Mn, Fe and Ni). Phys. Rev. B 46, 5425. https://doi.org/10.1103/PhysRevB.46.5425 (1992).

15. Jeevanandam, P. & Vasudevan, S. Magnetism in MnPSe₃; a layered 3d⁵ antiferromagnet with unusually large XY anisotropy. J. Phys.: Condens. Matter 11, 3563. https://doi.org/10.1088/0953-8984/11/17/314 (1999).

16. Basnet, R. et al. Controlling magnetic exchange and anisotropy by nonmagnetic ligand substitution in layered MPX₃ (M = Ni, Mn; X = S, Se). Phys. Rev. Res. 4, 023256. https://doi.org/10.1103/PhysRevResearch.4.023256 (2022).

17. Han, H. et al. Emergent mixed antiferromagnetic state in MnPS₃(1−x)Se₃x. Appl. Phys. Lett. 122, 033101. https://doi.org/10.1063/5.0135557 (2023).

18. Calder, S., Haglund, A., Kolesnikov, A. I. & Mandrus, D. Magnetic exchange interactions in the van der Waals layered antiferromagnet MnPSe₃. Phys. Rev. B 103, 024414. https://doi.org/10.1103/PhysRevB.103.024414 (2021).
19. Goossens, D. J., Wildes, A., Ritter, C. & Hicks, T. J. Ordering and the nature of the spin flop phase transition in MnPS₃. *J. Phys.: Condens. Matter* **12**, 1845. https://doi.org/10.1088/0953-8984/12/8/327 (2000).

20. Wiedenmann, A., Rossat-Mignod, J., Louisy, A., Brec, R. & Rouxel, J. Neutron diffraction study of the layered compounds MnPSe₃ andFePSe₃. *Solid State Commun.* **40**, 1067–1072. https://doi.org/10.1016/0038-1098(81)90253-2 (1981).

21. Jana, D. et al. Magnon gap excitations and spin-entangled optical transition in van der Waals antiferromagnetNiPS₃. *Phys. Rev. B* **108**, 115149. https://doi.org/10.1103/PhysRevB.108.115149 (2023).

22. Cho, C. et al. Microscopic parameters of the van der Waals CrSBr antiferromagnet from microwave absorption experiments. *Phys. Rev. B* **107**, 094403. https://doi.org/10.1103/PhysRevB.107.094403 (2023).

23. Momma, K. & Izumi, F. VESTA 3 for three-dimensional visualization of crystal, volumetric and morphology data. *J. Appl. Crystallogr.* **44**, 1272–1276. https://doi.org/10.1107/S0021889811038970 (2011).

24. Balkanski, M., Jouanne, M., Ouvrard, G. & Scagliotti, M. Effects due to spin ordering in layered MPX₃ compounds revealed by inelastic light scattering. *J. Phys. C: Solid State Phys.* **20**, 4397. https://doi.org/10.1088/0022-3719/20/27/017 (1987).

25. See Supplemental Material at [URL will be inserted by publisher] for details on THz monochromatic transmission measurements, Raman scattering measurements in out-of-plane configuration, Raman scattering measurements in in-plane configuration till 30 T and simulation for in-plane magnetic field dependence of magnon gap of MnPSe₃.

26. Eibschütz, M., Shtrikman, S. & Treves, D. Internal field in orthoferrites and the one third power law. *Solid State Commun.* **4**, 141–145. https://doi.org/10.1016/0038-1098(66)90213-4 (1966).

27. Nagamiya, T., Yosida, K. & Kubo, R. Antiferromagnetism. *Adv. Phys.* **4**, 1–112. https://doi.org/10.1080/00018735500101154 (1955).

## Acknowledgements

The work has been supported by the EC Graphene Flagship Core 3 project. M.P and M.B acknowledge the support from the EU ERC TERAPLASM Grant (No. 101053716) and CENTERA projects (MAB/2018/9 and FENG.02.01-IP.05-T004/23) founded within the IRA program of the FNP-Poland, co-financed by the EU, SG OP and FENG Programmes.

## Author contributions

D.J., J.-Ph., B.A.P., M.O., C.F., and M.P. conceived the project and experiments. D.J., D.V., I.M., P.K., C.W.C., I.B., and M.B. conducted the experiments and analyzed the results. The manuscript, initially written by D.J. and M.P., has been reviewed by all authors.

## 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/s41598-024-67356-4.

**Correspondence** and requests for materials should be addressed to M.P.

**Reprints and permissions information** is available at 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_6.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 licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence 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 licence, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2024