
==== Front
Nat Commun
Nat Commun
Nature Communications
2041-1723
Nature Publishing Group UK London

51643
10.1038/s41467-024-51643-9
Article
Unveiling the spin evolution in van der Waals antiferromagnets via magneto-exciton effects
http://orcid.org/0000-0002-0929-854X
Wang Xingzhi xzwang@xmu.edu.cn

12
http://orcid.org/0000-0002-5829-9119
Tan Qishuo 2
Li Tie 1
http://orcid.org/0000-0001-9451-9919
Lu Zhengguang 34
http://orcid.org/0000-0002-8439-7708
Cao Jun 2
Ge Yanan 1
Zhao Lili 5
http://orcid.org/0000-0003-2365-2952
Tang Jing 2
http://orcid.org/0000-0001-7622-0577
Kitadai Hikari 2
Guo Mingda 6
Li Yun-Mei 1
http://orcid.org/0000-0002-3014-756X
Xu Weigao 5
http://orcid.org/0000-0003-0166-2172
Cheng Ran 67
http://orcid.org/0000-0001-6358-3221
Smirnov Dmitry 3
http://orcid.org/0000-0003-3462-9088
Ling Xi xiling@bu.edu

289
1 https://ror.org/00mcjh785 grid.12955.3a 0000 0001 2264 7233 Department of Physics, Xiamen University, Xiamen, PR China
2 https://ror.org/05qwgg493 grid.189504.1 0000 0004 1936 7558 Department of Chemistry, Boston University, Boston, MA USA
3 https://ror.org/03s53g630 grid.481548.4 0000 0001 2292 2549 National High Magnetic Field Laboratory, Tallahassee, FL USA
4 https://ror.org/05g3dte14 grid.255986.5 0000 0004 0472 0419 Department of Physics, Florida State University, Tallahassee, FL USA
5 https://ror.org/01rxvg760 grid.41156.37 0000 0001 2314 964X Key Laboratory of Mesoscopic Chemistry, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing, China
6 grid.266097.c 0000 0001 2222 1582 Department of Physics and Astronomy, University of California, Riverside, CA USA
7 grid.266097.c 0000 0001 2222 1582 Department of Electrical and Computer Engineering, University of California, Riverside, CA USA
8 grid.189504.1 0000 0004 1936 7558 Division of Materials Science and Engineering, Boston University, Boston, MA USA
9 https://ror.org/05qwgg493 grid.189504.1 0000 0004 1936 7558 The Photonics Center, Boston University, Boston, MA USA
13 9 2024
13 9 2024
2024
15 801124 7 2023
14 8 2024
© The Author(s) 2024
2024
https://creativecommons.org/licenses/by-nc-nd/4.0/ Open Access This article is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License, which permits any non-commercial use, sharing, 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 you modified the licensed material. You do not have permission under this licence to share adapted material derived from this article or parts of it. 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-nc-nd/4.0/.
Among the fascinating phenomena observed in two-dimensional (2D) magnets, the magneto-exciton effect stands out as a pivotal link between optics and magnetism. Although the excitonic effect has been revealed and exhibits a considerable correlation with the spin structures in certain 2D magnets, the underlying mechanism of the magneto-exciton effect remains underexplored, especially under high magnetic fields. Here we perform a systematic investigation of the spin-exciton coupling in 2D antiferromagnetic NiPS3 under high magnetic fields. When an in-plane magnetic field is applied, the exceptional sharp excitonic emission at ~1.4756 eV exhibits a Zeeman-like splitting with g ≈ 2.0, experimentally identifying the exciton as an excitation of dominant triplet-singlet character. By examining the polarization of excitonic emission and simulating the spin evolution, we further verify the correlation between excitonic emission and Néel vector in NiPS3. Our work elucidates the mechanism behind the spin-exciton coupling in NiPS3 and establishes a strategy for optically probing the spin evolutions in 2D magnets.

NiPS3, a van der Waals antiferromagnet exhibits exciton emission with a very sharp linewidth. The exact origin of this is has been a subject of active debate. Here, Wang et al study the behavior of this sharp exciton peak under applied magnetic fields, and find a Zeeman-like splitting, indicating the exciton has triplet-singlet character.

Subject terms

Condensed-matter physics
Materials for optics
https://doi.org/10.13039/100000001 National Science Foundation (NSF) 1945364 Ling Xi https://doi.org/10.13039/100000015 U.S. Department of Energy (DOE) DE-SC0021064 Ling Xi issue-copyright-statement© Springer Nature Limited 2024
==== Body
pmcIntroduction

Since the discovery of monolayer magnets by magneto-optical methods, there has been a tremendous surge of interest in 2D magnetic materials in recent years1–7. These studies have led to the discovery of novel physics phenomena, such as correlated-electron physics8,9, ferromagnetic-superconducting transitions10, and magnon topology11. These discoveries provide an intriguing platform for the development of innovative spintronic devices and quantum information technologies. Among the various 2D magnets, NiPS3 has emerged as a strongly correlated antiferromagnet with complex interactions among spins, electrons, phonons, and magnons8,9,12–15. One of its notable features is the observation of an exceptionally sharp emission with an intrinsic linewidth on the order of hundreds of micro-electron-volts in the antiferromagnetic phase16–18. Furthermore, the excitonic emission shows a strong correlation with the spin order, suggesting the potential for innovative 2D antiferromagnet-based opto-spintronic devices19. The widely accepted explanation for the origin of this emission is Zhang-Rice excitons, a type of many-body exciton16, which is also theoretically understood using density functional theory recently20,21. However, other systematic studies, both experimental and simulation-based, provide substantial evidence to suggest that the emission may originate from either the intrinsic band structures or defect-bound states in NiPS318,22,23. Additionally, there have been conflicting reports regarding the origin of spin-induced linear polarization of exciton emission16–18, which is thought to arise from the spin orientation17 and the ferromagnetic zig-zag chain16,18,22. The controversy surrounding the mechanism of excitonic emission and its coupling to magnetism hinders advanced investigations in 2D magnetism and the development of future opto-magnetic devices24–27. Furthermore, the optical absorption studies on this excitonic transition reveal its resilience to spectral changes in both in-plane and out-of-plane magnetic fields up to 7 T16. Higher magnetic fields are necessary to amplify magneto-exciton effect, which can even result in the rotation or flipping of spins in 2D antiferromagnetic systems.

Here we conduct a systematic magneto-optical study on 2D antiferromagnetic NiPS3 under high magnetic fields to explore the spin-exciton coupling in the materials. The spin-correlated excitonic emission exhibits clear splitting features with applied magnetic fields in the Voigt geometry, accompanied by the spin rotation within the layer plane. The magnetic-field-dependent spectroscopy and polarization features of excitonic emission are further corroborated by precise simulation of spin evolution. Additionally, we also investigate the magneto-optical response of excitonic emission in the Faraday geometry, reaching up to 31 T. The demonstrated magneto-exciton effect in our study provides a deeper understanding of this new type of excitons and their coupling to spins in 2D NiPS3, and will advance fundamental exploration in the field of 2D magnets.

The atomic layers of NiPS3 stack in the C2/m space group, which belongs to the monoclinic crystal system28. When the temperature decreases across its Néel temperature (TN = ~152 K), NiPS3 undergoes a magnetic phase transition and forms a 2D antiferromagnetic XXZ model on a honeycomb lattice (Fig. 1a)12. Each spin is mainly aligned either in parallel or antiparallel along the in-plane a axis, with a slight out-of-plane component resulting in a tilt of the spin from a-b plane at an angle of approximately 7°29. Moreover, the intra-layer spin moments are arranged ferromagnetically in chains along the a axis, but are coupled antiferromagnetically with their neighboring chains. The spin structure of NiPS3 indicates that the a axis of the lattice is the easy axis of the spin structure29,30. Considering the XXZ-type spin structure of NiPS3, the anisotropic energy is small when the spins align along directions within the a-b plane30. Such spin alignment can be achieved by applying in-plane magnetic fields.Fig. 1 Magneto-exciton effect of NiPS3 with external field direction along a axis.

a Top view of lattice and spin structure of an NiPS3 single layer. The crystal orientations (a, b and c axes) are denoted by brown arrows, with the spin orientation primarily along the a axis of the crystal. The black arrow indicates the direction of the applied magnetic field, characterized by the angle θ relative to the a axis of the crystal in the Voigt geometry. b PL spectrum from exfoliated NiPS3 flakes excited by a 532-nm continuous-wave laser at 1.6 K. The sharp excitonic peak is labeled as X, accompanied by two small side peaks labeled as Sα and Sβ. The inset displays an optical microscopy image of exfoliated NiPS3 flakes on a SiO2/Si substrate, with the red point indicating the area where the signal was collected. c Magnetic-field-dependent PL color map of NiPS3 with the external field applied along a axis in the range from −31 to 31 T. d Magnetic-field-dependent PL spectra of NiPS3 with the external field applied along a axis in the range from 0 to 31 T. e Energy shift of the splitting excitonic peaks, X+ (ΔEX+) and X- (ΔEX-), as a function of the magnitude of the applied field. f Full width at half maximum (FWHM) of the excitonic emission as a function of the magnitude of applied field. Error bars represent standard errors from fitting function.

NiPS3 flakes were mechanically exfoliated onto SiO2/Si substrates from as-grown bulk single-crystals28. The photoluminescence spectrum of NiPS3 shows a sharp emission (labeled as X) located at ~1.4756 eV with a linewidth of ~270 μeV at 1.6 K (Fig. 1b). The measured sample has a thickness of ~20 nm, as indicated by the atomic force microscopy measurement (Fig. S1). The PL mapping measurement indicates that its line-width is homogeneous on the flakes (Fig. S2), which contradicts the localized emission from an isolated impurity or a single emitter31,32. In addition to the main excitonic peak, two smaller side peaks labeled as Sα and Sβ were observed at 1.4768 and 1.4783 eV, respectively (Fig. S3). The energies of two side peaks are both ~0.5 meV lower than the reported values16, which could due to the sample variance. In contrast, the energy of the excitonic peak X is consistent with the reported values in previous studies16–18. An important characteristic of the excitonic emission is its near-unity linear polarization, which has been attributed to the anisotropy of the spin structure16–18,33,34. The excitonic emission intensity reaches a maximum when collecting polarization is perpendicular to the a axis17,22. However, within each NiPS3 layer, both the spin orientation and the zigzag ferromagnetic chains align along the a axis and break the three-fold rotation symmetry of the lattice. This results in ambiguity regarding the origin of the linear polarization of X16–18,22. It is worth noting that the orientation of spins in 2D antiferromagnets is of significant interest in the community of 2D magnetism, as it plays a crucial role in presenting information in the antiferromagnet-based devices and technologies24–27. Therefore, these two hypotheses represent distinct underlying physical mechanisms of spin-exciton coupling, which will shape the future development of opto-magnetic and opto-spintronic devices24. One strategy to differentiate the spin orientation from the zig-zag chain direction is to induce spin rotation through external modulations, such as the application of an external magnetic field16,17.

For this purpose, we conducted systematic magneto-optical measurements on NiPS3 bulk crystals in the Voigt geometry. We define θ as the angle between the external field direction and the a axis (Fig. 1a). When the magnetic field is applied along the a axis (θ = 0°), the excitonic peak splits into two branches, labeled as X+ and X-, in the range from −31 to 31 T (Fig. 1c, d). The splitting energy, ΔE = EX+ – EX-, linearly increases with the increase of applied field from 0 to 11 T, reaching a maximum of ~ 2.15 meV. However, when the magnitude of applied field increases above 11 T, the splitting energy drops dramatically and stabilizes at ~0.49 meV (Fig. 1e). The spectral behavior of X implies a wealth of information about the magneto-exciton effect in 2D antiferromagnets. Considering the spin-corrected nature of excitonic emission, we hypothesize that the abrupt change in the splitting energy originates from a field-induced magnetic phase transition from an antiferromagnetic phase to a spin-flop phase35,36. The magnetic phase transition hypothesis is further supported by the field-dependent full width at half maximum (FWHM) of excitonic emission (Fig. 1f). The measured FWHM of X in this set of magneto-optical measurements, even at B = 0 T, is obviously larger than the real value (Fig. 1b), mainly due to the chosen low-resolution grating (see Method section). However, the linewidth still exhibits a clear broadening in the range from 12 to 20 T, suggesting the presence of spin disorder during the magnetic phase transition, which has the same origin as the abrupt change in the splitting energy. In contrast, the original FWHM of excitonic peak at B < 11 T and B > 20 T may suggest a pure magnetic phase. Similar behavior can be also found in the extracted intensity of X (Fig. S4).

Although the spin orientation has been verified to align along the a axis of the NiPS3 crystal29, the XXZ nature of the spin structure still indicates a weak magnetic anisotropy in the a-b plane12,37. To further explore the field-induced energy splitting of the excitonic peak, we conducted series of magneto-PL measurements by applying the magnetic field with the magnitudes ranging from 0 to 12 T, while along different directions in the a-b plane of NiPS3. The energy splitting effect of X from the exfoliated NiPS3 flake exhibits distinct features at various angles (θ) (Fig. 2a). In contrast to the result when θ = 0°, the abrupt change in the spectra is absent with other six θ values from 15° to 90°, where the splitting energy shows a non-linear dependence on the applied field and a gradual variation in the range from 0 to 12 T. The energy splitting effects become inconspicuous as the field direction deviates from the easy axis, and no obvious peak splitting can be observed when the magnetic field is applied along the b axis (i.e., θ = 90°). The FWHM of X+ and X- is set at the same value during the fitting based on their shared origins. A clear broadening is observed from the bulk NiPS3 crystal with an external field up to 31 T along the a axis (Fig. 1f). Another notable characteristic of peak X in the spin-flop phase is the asymmetric behaviors in the peak intensities of the two splitting branches. In most of cases, the intensity of the low-energy branch is weaker than that of the high-energy branch, with an intensity ratio IX-/IX+ below 1. However, the intensity ratio dramatically increases at B > 10 T only when θ = 0° (Fig. S5). Similar results are also observed in another set of magneto-PL study conducted on bulk NiPS3 with magnetic fields ranging from −14 to 14 T in the Voigt geometry (Fig. S6).Fig. 2 Energy splitting of excitonic emission in the Voigt geometry.

a Magneto-optical measurements of exfoliated NiPS3 flakes in the Voigt geometry with the external field applied along different angles within the a-b plane. b Schematic of excitonic energy states in the antiferromagnetic (AFM) phase, featuring the quantum number difference Δms=±1 in the transition. c Schematic of excitonic energy states in the canted antiferromagnetic (CAFM) phase, featuring the quantum number difference Δms=±δ in the transition, where δ < 1.

The excitonic peak X can be well understood in the framework of Zhang–Rice exciton, which originates from the transition from Zhang–Rice singlet (ZRS) to Zhang–Rice triplet (ZRT) states16. ZRT state represents a state where spin directions at the Ni and S sites are the same, while ZRS state represents a state where spin directions at the Ni and S sites are opposite. The magnetic behaviors of X remain a variety of controversy within the framework of Zhang-Rice exciton. For instance, the spectral features of sharp emission have been reported to be insensitive to the magnetic field, which were considered to conform to the Zhang-Rice exciton nature16. However, a recent study presented a distinct phenomenon: a clear splitting effect of the X transition under in-plane fields38. Contrary to previous perspectives on the robustness of the spectral features of Zhang-Rice excitons, this remarkable response to magnetic fields was treated as clear evidence to rule out the Zhang-Rice picture as the origin38. Unlike the previous viewpoints on the magnetic response of Zhang-Rice excitons, in this work, we argue that the spectral evolution of X peak reveals its nature of singlet-triplet transition, which is consistent with the framework of Zhang-Rice exciton.

When the exciton transitions from singlet to triplet states, the spin orientations experience a spin flip, resulting in a change in quantum number from ms = 0 to ms = ±1 (Fig. S7a, b). Consequently, the triplet state with ms = ±1 splits under the external magnetic field, leading to the splitting of X (Fig. 2b). Meanwhile, with the application of an in-plane magnetic field, the spins rotate in-plane, and the spin structure of NiPS3 changes from the antiferromagnetic to canted antiferromagnetic phases when θ ≠ 0° (Fig. S7c, d). In the canted antiferromagnetic phase, the spin orientations at the Ni and S sites vary in two directions with an angle smaller than 180°. Therefore, the excitonic transition in the canted antiferromagnetic phase results in a quantum number difference Δms = ± δ = ± |M1 – M2 | /2 = ± |L | /2, where M1 and M2 represent two different spin momenta in NiPS3; and L = M1 – M2 is the Néel vector. Considering the antiferromagnetic nature of the spin structure, we assume the magnitude of each spin momentum as unity, i.e., M1 = M2 = 1, and we have δ < 1 in the canted antiferromagnetic phase (Fig. 2c). Compared to the split triplet state with ms = +δ, the ms = –δ state has lower energy and consequent higher carrier density with an application of positive magnetic field, resulting in a preference for transition from singlet and a higher peak intensity of X+ in the region of low magnetic field (Fig. S5). In contrast, the reversed intensity ratio between two splitting branches at B > 10 T and θ = 0° suggests that the energy diagrams described in Fig. 2b, c is no longer applicable under in the case of high applied fields. A stronger peak of low-energy transition potentially indicates a splitting of excited state, and further implies an entirely different energy diagram in the spin-flop phase.

The energy splitting of X branches (X+ and X−) is proposed to be described by the following equation,1 EX+/X−=EX0±gsμBBcosφΔms=EX0±12gsμBB⋅L,

where gs ≈ 2.0 is the Landé g-factor for electrons, and φ is defined as the angle between B and L. When θ ≠ 0, with the increase of B, the spin momenta M1 and M2 experience an in-plane rotation, which is determined by the competition between super-exchange coupling and magnetic anisotropy constant36 (Fig. 3a). Therefore, φ is also expected to change with the increase of B in the Voigt geometry, resulting in a variation in the splitting energy. Here we define k = gsΔms cos φ to describe the slope of magnetic-field-dependent energy splitting of X. With B = 0, we have φ = θ and Δms = ±1 for the transition from singlet to triplet states. As a result, on the premise that the spin rotation could be neglected at a low field (e.g. B ≤ 4 T), the energy shift of X branches is nearly linearly dependent on the field magnitude for various applied field directions, where k = gs cos θ. Thus, we extract the k values from experimental data using k =ΔE/2μBB under B = 2, 3 and 4 T (Fig. 3b). The values of k match well with the calculated value, where we obtain gs ≈ 1.90 ± 0.02 by fitting. The result indicates that the splitting energy of peak X is determined by the effective field applied along the Néel vector direction.Fig. 3 Simulated spin evolution in 2D antiferromagnetic NiPS3.

a Schematic of spin evolution under an in-plane magnetic field. The blue, red, and orange arrows represent the spin momenta M1, M2, and Néel vector L, with angles φ1, φ2, and φ with respect to the positive direction of the applied field (depicted as a black arrow). b Magnetic-field-direction-dependent k values, which are derived from the experimental data using the formula k=ΔE/2μBB. The dashed line represents the fitting line based on k=gscosθ. Error bars represent standard errors from fitting function. c Simulated k value as a function of the magnitude B and direction θ of the magnetic field. k is defined as gsΔmscosφ. d Simulated energy splitting ΔE as a function of the magnitude B and direction θ of the magnetic field. e Magnetic-field-dependent energy splitting extracted from experimental data by ΔE=EX+−EX−, compared with calculated results. The experimental data of θ = 0o is compared with calculated data of θ = 2o.

In order to quantitatively understand the excitonic transition under an external field, the spin states with a given applied in-plane field, described by B and θ, have been simulated. We consider a two-sublattice antiferromagnet characterized by two macrospin, S1 and S2, which correspond to the magnetic momentums M1 and M2 with unity magnitudes in NiPS3. M1 and M2 are antiparallel when the magnetic field is zero, representing the initial antiferromagnetic phase. The general Hamiltonian of the antiferromagnetic system can be written as36:2 H=JM1⋅M2−A2M1⋅x^2+M2⋅x^2−gμBB⋅M1+M2,

where J is the antiferromagnetic exchange coupling, and A is the easy axis anisotropy with x^ refers to the easy axis of spin structure. The ground state can be obtained by minimizing the Hamiltonian. In an approximation, we can use this model to simulate the spin structure evolution of NiPS3 under an in-plane field. To simplify the calculation, we convert this equation into a function of three angles θ, φ1, and φ2 (see Methods). φ1 and φ2 are defined as the angles of M1 and M2 concerning the positive direction of the external field (Fig. 3a). The constant k = gs Δms cos φ and splitting energy ΔE = 2kμBB are simulated as a function of B and θ (Fig. 3c–d and Fig. S8). Also, ΔE are extracted from the experimental data in Fig. 2a through ΔE = EX+ – EX–. We compare the experimental results with the calculated values (Fig. 3e). Of particular note, the experimental data with θ = 0° matches better with the simulated curve of θ = 2° (Fig. S9). This minor angle difference may result from experimental errors in aligning the sample’s a axis with the magnetic field direction. For all the directions of the applied field, the simulated results strongly agree with the experimental data at low field region of B < 10 T, while they start to deviate at higher field region. The deviation may be attributed to emergence of complex local spin structures in NiPS3 induced by the high field, causing considerable variations on excitonic energy levels20,39. The small total energy differences among different magnetic configurations at the level of several to tens of meV might be reached by the application of an external field20.

In addition to the spectral features under in-plane fields, the polarization of the excitonic emission also carries rich spin information in NiPS3. Previous studies have reported that the intensity of X is nearly independent of the incident laser polarization but strongly dependent on the collection light polarization16–18. The strong coupling between exciton dipole and local spin originates from the spin-flip during the singlet-triplet transition, resulting in the exciton dipole being perpendicular to the spin orientation17. In light of this, we measure the polarization-dependent photoluminescence spectra by rotating the analyzer angle with varying direction (θ) and magnitude (B) of the applied field (Fig. 4a and Fig. S10). Across all configurations, two splitting branches exhibit clear linear polarizations, with the maximal intensity consistently occurring at the same polarization angle (Fig. S11). The polarization angles of both branches further rotate with the increase of the in-plane field (Fig. 4b). The measured rotation of the polarization angle of X can be attributed to two roots: the rotation of the Néel vector induced by the in-plane magnetic field and the Faraday rotation of collection light passing through the objective lens along the magnetic field direction. To exclude the contribution from Faraday rotation, we conduct a controlled experiment described in Methods and Fig. S12. By removing the influence of Faraday rotation, the measured polarization of the X can be used as a probe to determine the Néel vector orientation in NiPS317. The Néel vector orientation can be calculated as φ = φ1 + φ2 + 90°, where φ1 and φ2 are also functions of B and θ (Fig. S13). Notably, the field-dependent polarization angle of the X with different θ exhibits excellent agreement with the simulation results (Fig. 4c), confirming the effectivity of our simulation in capturing both the energy splitting ΔE and the Néel vector orientation φ. This enables us to obtain a detailed understanding of the spin evolution in 2D antiferromagnetic NiPS3, as illustrated by typical spin configurations at various applied magnetic fields (Fig. 4d). Our magneto-optical investigation of the excitonic transition offers a multi-angle approach to study the unique spin-exciton coupling and magnetic behaviors in NiPS3.Fig. 4 Exciton polarization evolution and spin evolution probe in NiPS3.

a Collection-polarization-dependent PL color map with applied fields along the a axis under different magnetic fields at 5 K. b Polar plots of the intensity of X- branch as a function of the collection polarization angle at selected magnitudes of the applied field. In each panel, the black dashed line represents the polar fitting without an external field. The polarization angle of the X- rotates with the increase of magnitude of the applied field. c PL polarization angle of X as a function of the in-plane magnetic field with different θ, compared with calculated PL polarization (labeled with solid curves) by φCal-X = φL + 90°, where φL is the simulated Néel vector direction. The experimental polarization angles are the value after excluding the Faraday rotation effect. The experimental data of θ = 0o is compared with calculated data of θ = 2o. In the four configurations, the applied field direction is along 90o, which is denoted by the black dashed line. Error bars represent standard errors from fitting function. d Simulated spin evolution under magnetic field with different magnitudes B and directions θ. The red and blue arrows represent the spin momenta M1 and M2, respectively.

In contrast to the sensitive response of the exciton to the magnetic field in the Voigt geometry, the excitonic emission becomes robust to the field in the Faraday geometry, where the applied field is nearly perpendicular to the spin direction in NiPS3 (Fig. 5a). We notice that X does not exhibit a clear variation at B < 10 T (Fig. 5b), which is consistent with previous results in the literature16. The anisotropy between the in-plane and out-of-plane directions should lead to considerable difficulty in rotating the spins by an out-of-plane field in the present experimental conditions. However, considering the small out-of-plane component of the spin momentum of NiPS3, the magneto-optical response of the excitonic transition is still expected to be observed, even though a strong external field is required. As a result, X exhibits an energy shift of ~ 97 μeV with a quenched intensity at B = 31 T, compared with the case of zero magnetic field (Fig. 5c). Moreover, the field-dependent peak position and FWHM exhibit nonlinear relations (Fig. S14), which are similar to the magneto-exciton behavior of bright excitons in transition metal dichalcogenides in the Voigt geometry40. This similarity suggests that the slight peak shift of X possibly originates from the diamagnetic shift and implies the existence of a dark state at higher energy that coupling with the X state. The X emission energy (EX) can be fitted well using the following formula40,3 EX=EX0+ΔE02−ΔE022+gFμBB22,

where EX0 and ΔE0 are the energy of X transition and the energy difference between X and higher dark state without the applied field; gF is the Landé g-factor for the transition in the Faraday geometry; μB is Bohr magneton; B is the magnitude of the applied field. We obtain ΔE0 = ~ 1.5 meV by fitting, which is matches the energy of Sα peak. However, two side peaks from the measured sample cannot be well resolved on the spectra, so we did not observe their field-dependent behaviors to further confirm whether Sα is the dark state for the bright excitonic state. Further investigation is necessary to verify this hypothesis.Fig. 5 Magneto-exciton effect of NiPS3 in the Faraday geometry.

a Schematic of the Faraday geometry with the applied field direction perpendicular to the a-b plane of NiPS3 crystal. The red and blue arrows represent the spin momenta, suggesting the XXZ-type antiferromagnetic spin structure of NiPS3. The black dashed lines represent the honeycomb spin structure, and the gray solid lines show the edge of the unit cell. b Magnetic-field-dependent PL spectra of NiPS3 in the Faraday geometry under 532-nm laser excitation at 5 K. The black dashed line indicates the peak position at B = 0 T for comparison. c Magnetic-field-dependent PL color map of NiPS3 in the Faraday geometry under 532-nm laser excitation at 5 K. The black dashed line represents the fitting curve of the magnetic-field-dependent peak position of X emission.

Uncovering the underlying physics of antiferromagnetic spin states, especially in the 2D cases, poses a significant challenge due to the absence of net magnetism41–43. In this work, we observed a pronounced magneto-exciton effect in NiPS3, leading to substantial variations in the optical spectra. The exceptional sharp emission exhibits a Zeeman-like splitting under an in-plane magnetic field with a Landé g-factor around 2, providing strong experimental evidence to identify the exciton mode as an excitation of dominant triplet-singlet character. Moreover, the field-induced polarization rotation of exitonic emission is unambiguously proven to be correlated with the Néel vector rotation, which has been further supported by a simulation of spin evolution. Our investigation into the magneto-exciton coupling in NiPS3 sheds light on new facets of magnetism and opto-magnetism at the extreme nanoscale limit. The unique characteristics of antiferromagnets, such as terahertz resonance, multilevel states, and the absence of stray fields, make them highly attractive for spintronic devices. Leveraging the strong correlation between spins and excitons, optical methods could be employed to probe and manipulate the spin orders, paving the way for the further advanced studies to achieve a deeper understanding of fundamental physics in the field of 2D magnetism.

Note: Recently, a related work by Jana et al. appeared, reporting similar magneto-PL results38.

Methods

Sample preparation and characterization

NiPS3 single crystals were grown using a chemical vapor transport (CVT) method28. NiPS3 flakes were mechanically exfoliated onto Si substrates with a 285-nm SiO2 layer from a bulk single-crystal. The morphology and thickness of exfoliated flakes were characterized by optical microscopy (Olympus BX53) and atomic force microscopy (Bruker NW4) in a tapping mode. The crystal orientations (a and b axes) of NiPS3 sample were determined by polarized PL measurements, where the excitonic emission intensity reaches a maximum with collecting polarization along the b axis17,22.

Magneto-optical measurements

The magneto-PL experiments on exfoliated NiPS3 flakes were performed in a closed cycle cryostat (Attocube attoDRY2100). The magnetic field was applied in the Voigt geometry from 0 to 12 T. A 532-nm laser was used to excite the sample. PL signals were collected by a low-temperature objective, propagate in a free beam throughout collecting path, and measured by a spectrometer (Horiba iHR550) equipped with a liquid-nitrogen-cooled charge-coupled device (CCD). A high-resolution 1800g/mm grating was chosen to disperse the PL signal with a spectral resolution of ~70 μeV. The magneto-PL experiments of bulk NiPS3 crystals were performed in a 14 T DC magnet and a 31 T DC magnet in the National High Magnetic Field Laboratory. A 532-nm laser was used to excite the sample. In the Faraday geometry, the sample was horizontal placed in the center of the magnetic cell with the surface perpendicular to the applied magnetic field. In the Voigt geometry, a mirror was set between the objective and the sample with 45° to guide the optical path by 90°, and the sample was vertically placed in the center of the magnetic cell with the surface parallel to the applied magnetic field. A converted 10× objective in the magnetic cell was used to focus the laser beam on the sample and collect the PL signal. Both the incident and collection lights propagate in a free-beam in the magnetic cell. For the polarization-dependent PL measurement, the collected signal was filtered by an analyzer, subsequently collected by a multi-mode optical fiber, and measured by a spectrometer (Princeton Instruments, IsoPlane 320) equipped an electrically cooled CCD. A 1800g/mm grating was chosen to disperse the PL signal in the measurement in the Faraday geometry and a 1200 g/mm grating was selected in the measurement in the Voigt geometry.

Determination of the Faraday rotation of the PL signal

An ideal method to exclude the contribution of Faraday rotation is to conduct a control experiment with θ = 0°, where the Néel vector is considered unchanged before the field-induced phase transition to the spin-flop phase happens17. The Faraday rotation at different applied fields can be obtained by measuring the rotation of the polarization angle of excitonic peak X with θ = 0°. However, accurately aligning the NiPS3 sample with the a axis of the crystal along the direction of the applied field is challenging, the Néel vector will exhibit a considerable rotation (φ) even with a small value of θ, as shown in our simulation (Fig. S12a). The effect caused by the experimental error in θ can be estimated by the simulated value of dφ/dθ. In the range from 0 to 10 T, the value of dφ/dθ is relatively small when θ is around 30°. Therefore, in this work, we mount two samples with θ = 30° and −30°, and measured their polarization rotations of X emission, α1 and α2. In this case, these two samples will have the opposite Néel vector rotation and the same Faraday rotation. Thus, the contribution from the Néel vector rotation due to the small error in sample alignment will be canceled out. We obtain a Faraday rotation φF = (α1 + α2)/2 = ~2.305°/T in this work, which applies to all θ values for the same magnetic cell and optical path (Fig. S12b).

Simulation of spin rotation with the external magnetic field

To calculate the splitting energy with given B and θ, we consider a two-sublattice antiferromagnet featured by two macrospin, S1 and S2, manipulated by in-plane magnetic fields along different angles. S1 and S2 are antiparallel when the magnetic field is zero. The free energy density of two dimensionless vectors m1 and m2 can be written as36:4 E=Jm1⋅m2−A2m1⋅x^2+m2⋅x^2−MSH0⋅(m1+m2),

where m1 = S1/S1; m2 = S2/S2; H0 is the external magnetic field; x^ refer to the easy-axis; J and A are the antiferromagnetic exchange coupling and the easy axis anisotropy, respectively.

We can then simplify the free energy density as a function of three angles:5 E=ℏ2ωJcosθ1−θ2−ωA2cos2θ1+cos2θ2−ℏωH2cosϕ−θ1+cosϕ−θ2

where ℏωJ = J, ℏωA = A, ℏωH = H0Ms; θ1, θ2, and ϕ are the angles of m1, m2, and in-plane magnetic field B to the easy axis. Here we scale all parameters into angular frequencies for calculation. We compensate for the calculated results by multiplying constants to match the experimental data.

Supplementary information

Supplementary Information

Peer Review File

Source data

Source Data

Supplementary information

The online version contains supplementary material available at 10.1038/s41467-024-51643-9.

Acknowledgements

This material is based upon work supported by the National Science Foundation (NSF) under Grant No. (1945364). Work by X.L. and J.C. was supported by the U.S. Department of Energy (DOE), Office of Science, Basic Energy Sciences (BES) under Award DE-SC0021064. X.L. acknowledges the membership of the Photonics Center at Boston University. X.W. acknowledges the support of the National Key R&D Program of China (2022YFA1602704), the National Natural Science Foundation of China (62275225), the Fundamental Research Funds for the Central Universities under grant (20720220034), and the 111 Project (B16029). Z.L. and D.S. acknowledge the support from the U.S. Department of Energy (No. DE-FG02-07ER46451) for magneto-spectroscopy measurements performed at the National High Magnetic Field Laboratory, which is supported by National Science Foundation through NSF/DMR-1644779 and the State of Florida. M.G. and R.C. are supported by Air Force Office of Scientific Research under Grant No. (FA9550-19-1-0307).

Author contributions

X.W. and X.L. conceived the experiment. X.W., T.L., Z.L., and D.S. conducted the magneto-PL measurements. Q.T., J.C., Y.G., and H.K. prepared samples and carried out the sample characterization. Q.T. and J.T. performed theoretical calculations under the supervision of M.G. and R.C.; L.Z. and W.X. performed the PL mapping measurement. X.W. and X.L. performed the analysis and interpretation of the data with significant input from Y.L. X.W., and X.L. wrote the manuscript with contributions from all authors.

Peer review

Peer review information

Nature Communications thanks the anonymous reviewers for their contribution to the peer review of this work. A peer review file is available.

Data availability

The data represented in Figs. 1–5 are provided with the paper as source data. All other data that support results in this article are available from the corresponding authors on reasonable requests. Source data are provided with this paper.

Competing interests

The authors declare no competing interests.

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
==== Refs
References

1. Huang B Layer-dependent ferromagnetism in a van der Waals crystal down to the monolayer limit Nature 2017 546 270 273 10.1038/nature22391 28593970
Huang, B. et al. Layer-dependent ferromagnetism in a van der Waals crystal down to the monolayer limit. Nature 546, 270–273 (2017).28593970 10.1038/nature22391
2. Gong C Discovery of intrinsic ferromagnetism in two-dimensional van der Waals crystals Nature 2017 546 265 269 10.1038/nature22060 28445468
Gong, C. et al. Discovery of intrinsic ferromagnetism in two-dimensional van der Waals crystals. Nature 546, 265–269 (2017).28445468 10.1038/nature22060
3. Wang X Raman spectroscopy of atomically thin two-dimensional magnetic iron phosphorus trisulfide (FePS3) crystals 2D Mater. 2016 3 031009 10.1088/2053-1583/3/3/031009
Wang, X. et al. Raman spectroscopy of atomically thin two-dimensional magnetic iron phosphorus trisulfide (FePS3) crystals. 2D Mater. 3, 031009 (2016).10.1088/2053-1583/3/3/031009
4. Lee J-U Ising-type magnetic ordering in atomically thin FePS3 Nano Lett. 2016 16 7433 7438 10.1021/acs.nanolett.6b03052 27960508
Lee, J.-U. et al. Ising-type magnetic ordering in atomically thin FePS3. Nano Lett. 16, 7433–7438 (2016).27960508 10.1021/acs.nanolett.6b03052
5. Burch KS Mandrus D Park J-G Magnetism in two-dimensional van der Waals materials Nature 2018 563 47 52 10.1038/s41586-018-0631-z 30382199
Burch, K. S., Mandrus, D. & Park, J.-G. Magnetism in two-dimensional van der Waals materials. Nature 563, 47–52 (2018).30382199 10.1038/s41586-018-0631-z
6. Mak, K. F., Shan, J. & Ralph, D. C. Probing and controlling magnetic states in 2D layered magnetic materials. Nat. Rev. Phys. 1, 646–661 (2019).
7. Gong C Zhang X Two-dimensional magnetic crystals and emergent heterostructure devices Science 2019 363 eaav4450 10.1126/science.aav4450 30765537
Gong, C. & Zhang, X. Two-dimensional magnetic crystals and emergent heterostructure devices. Science 363, eaav4450 (2019).30765537 10.1126/science.aav4450
8. Wang X Electronic Raman scattering in the 2D antiferromagnet NiPS3 Sci. Adv. 2022 8 eabl7707 10.1126/sciadv.abl7707 35030029
Wang, X. et al. Electronic Raman scattering in the 2D antiferromagnet NiPS3. Sci. Adv. 8, eabl7707 (2022).35030029 10.1126/sciadv.abl7707
9. Kim SY Charge-spin correlation in van der Waals antiferromagnet NiPS3 Phys. Rev. Lett. 2018 120 136402 10.1103/PhysRevLett.120.136402 29694193
Kim, S. Y. et al. Charge-spin correlation in van der Waals antiferromagnet NiPS3. Phys. Rev. Lett. 120, 136402 (2018).29694193 10.1103/PhysRevLett.120.136402
10. Li Y High-temperature quantum anomalous hall insulators in lithium-decorated iron-based superconductor materials Phys. Rev. Lett. 2020 125 086401 10.1103/PhysRevLett.125.086401 32909795
Li, Y. et al. High-temperature quantum anomalous hall insulators in lithium-decorated iron-based superconductor materials. Phys. Rev. Lett. 125, 086401 (2020).32909795 10.1103/PhysRevLett.125.086401
11. Viñas Boström E Direct optical probe of magnon topology in two-dimensional quantum magnets Phys. Rev. Lett. 2023 130 026701 10.1103/PhysRevLett.130.026701 36706407
Viñas Boström, E. et al. Direct optical probe of magnon topology in two-dimensional quantum magnets. Phys. Rev. Lett. 130, 026701 (2023).36706407 10.1103/PhysRevLett.130.026701
12. Kim K Suppression of magnetic ordering in XXZ-type antiferromagnetic monolayer NiPS3 Nat. Commun. 2019 10 345 10.1038/s41467-018-08284-6 30664705
Kim, K. et al. Suppression of magnetic ordering in XXZ-type antiferromagnetic monolayer NiPS3. Nat. Commun. 10, 345 (2019).30664705 10.1038/s41467-018-08284-6
13. Tan Q Charge-transfer-enhanced d–d emission in antiferromagnetic NiPS3 Appl. Phys. Rev. 2022 9 041406 10.1063/5.0107065
Tan, Q. et al. Charge-transfer-enhanced d–d emission in antiferromagnetic NiPS3. Appl. Phys. Rev. 9, 041406 (2022).10.1063/5.0107065
14. Ergeçen E Magnetically brightened dark electron-phonon bound states in a van der Waals antiferromagnet Nat. Commun. 2022 13 98 10.1038/s41467-021-27741-3 35013277
Ergeçen, E. et al. Magnetically brightened dark electron-phonon bound states in a van der Waals antiferromagnet. Nat. Commun. 13, 98 (2022).35013277 10.1038/s41467-021-27741-3
15. Dirnberger F Spin-correlated exciton–polaritons in a van der Waals magnet Nat. Nanotechnol. 2022 17 1060 1064 10.1038/s41565-022-01204-2 36097046
Dirnberger, F. et al. Spin-correlated exciton–polaritons in a van der Waals magnet. Nat. Nanotechnol. 17, 1060–1064 (2022).36097046 10.1038/s41565-022-01204-2
16. Kang S Coherent many-body exciton in van der Waals antiferromagnet NiPS3 Nature 2020 583 785 789 10.1038/s41586-020-2520-5 32690938
Kang, S. et al. Coherent many-body exciton in van der Waals antiferromagnet NiPS3. Nature 583, 785–789 (2020).32690938 10.1038/s41586-020-2520-5
17. Wang X Spin-induced linear polarization of photoluminescence in antiferromagnetic van der Waals crystals Nat. Mater. 2021 20 964 970 10.1038/s41563-021-00968-7 33903748
Wang, X. et al. Spin-induced linear polarization of photoluminescence in antiferromagnetic van der Waals crystals. Nat. Mater. 20, 964–970 (2021).33903748 10.1038/s41563-021-00968-7
18. Hwangbo K Highly anisotropic excitons and multiple phonon bound states in a van der Waals antiferromagnetic insulator Nat. Nanotechnol. 2021 16 655 660 10.1038/s41565-021-00873-9 33707746
Hwangbo, K. et al. Highly anisotropic excitons and multiple phonon bound states in a van der Waals antiferromagnetic insulator. Nat. Nanotechnol. 16, 655–660 (2021).33707746 10.1038/s41565-021-00873-9
19. Bae YJ Exciton-coupled coherent magnons in a 2D semiconductor Nature 2022 609 282 286 10.1038/s41586-022-05024-1 36071189
Bae, Y. J. et al. Exciton-coupled coherent magnons in a 2D semiconductor. Nature 609, 282–286 (2022).36071189 10.1038/s41586-022-05024-1
20. Klaproth T Origin of the magnetic exciton in the van der Waals antiferromagnet NiPS3 Phys. Rev. Lett. 2023 131 256504 10.1103/PhysRevLett.131.256504 38181357
Klaproth, T. et al. Origin of the magnetic exciton in the van der Waals antiferromagnet NiPS3. Phys. Rev. Lett. 131, 256504 (2023).38181357 10.1103/PhysRevLett.131.256504
21. Yan M Jin Y Voloshina E Dedkov Y Electronic correlations in FexNiyPS3 Van der Waals materials: insights from angle-resolved photoelectron spectroscopy and DFT J. Phys. Chem. Lett. 2023 14 9774 9779 10.1021/acs.jpclett.3c02688 37882477
Yan, M., Jin, Y., Voloshina, E. & Dedkov, Y. Electronic correlations in FexNiyPS3 Van der Waals materials: insights from angle-resolved photoelectron spectroscopy and DFT. J. Phys. Chem. Lett. 14, 9774–9779 (2023).37882477 10.1021/acs.jpclett.3c02688
22. Kim DS Anisotropic excitons reveal local spin chain directions in a van der Waals antiferromagnet Adv. Mater. 2023 35 2206585 10.1002/adma.202206585
Kim, D. S. et al. Anisotropic excitons reveal local spin chain directions in a van der Waals antiferromagnet. Adv. Mater. 35, 2206585 (2023).10.1002/adma.202206585
23. Liu, S., Malik, I. A., Zhang, V. L. & Yu, T. Lightning the spin: harnessing the potential of 2D magnets in opto-spintronics. Adv. Mater. 2306920 (2023) 10.1002/adma.202306920.
24. Němec P Fiebig M Kampfrath T Kimel AV Antiferromagnetic opto-spintronics Nat. Phys. 2018 14 229 241 10.1038/s41567-018-0051-x
Němec, P., Fiebig, M., Kampfrath, T. & Kimel, A. V. Antiferromagnetic opto-spintronics. Nat. Phys. 14, 229–241 (2018).10.1038/s41567-018-0051-x
25. Železný J Wadley P Olejník K Hoffmann A Ohno H Spin transport and spin torque in antiferromagnetic devices Nat. Phys. 2018 14 220 228 10.1038/s41567-018-0062-7
Železný, J., Wadley, P., Olejník, K., Hoffmann, A. & Ohno, H. Spin transport and spin torque in antiferromagnetic devices. Nat. Phys. 14, 220–228 (2018).10.1038/s41567-018-0062-7
26. Wang Z Determining the phase diagram of atomically thin layered antiferromagnet CrCl3 Nat. Nanotechnol. 2019 14 1116 1122 10.1038/s41565-019-0565-0 31712666
Wang, Z. et al. Determining the phase diagram of atomically thin layered antiferromagnet CrCl3. Nat. Nanotechnol. 14, 1116–1122 (2019).31712666 10.1038/s41565-019-0565-0
27. Han J Birefringence-like spin transport via linearly polarized antiferromagnetic magnons Nat. Nanotechnol. 2020 15 563 568 10.1038/s41565-020-0703-8 32483320
Han, J. et al. Birefringence-like spin transport via linearly polarized antiferromagnetic magnons. Nat. Nanotechnol. 15, 563–568 (2020).32483320 10.1038/s41565-020-0703-8
28. Du K Weak Van der Waals stacking, wide-range band gap, and Raman study on ultrathin layers of metal phosphorus trichalcogenides ACS Nano 2016 10 1738 1743 10.1021/acsnano.5b05927 26607168
Du, K. et al. Weak Van der Waals stacking, wide-range band gap, and Raman study on ultrathin layers of metal phosphorus trichalcogenides. ACS Nano 10, 1738–1743 (2016).26607168 10.1021/acsnano.5b05927
29. Wildes AR Magnetic structure of the quasi-two-dimensional antiferromagnet NiPS3 Phys. Rev. B 2015 92 224408 10.1103/PhysRevB.92.224408
Wildes, A. R. et al. Magnetic structure of the quasi-two-dimensional antiferromagnet NiPS3. Phys. Rev. B 92, 224408 (2015).10.1103/PhysRevB.92.224408
30. Lançon D Ewings RA Guidi T Formisano F Wildes AR Magnetic exchange parameters and anisotropy of the quasi-two-dimensional antiferromagnet NiPS3 Phys. Rev. B 2018 98 134414 10.1103/PhysRevB.98.134414
Lançon, D., Ewings, R. A., Guidi, T., Formisano, F. & Wildes, A. R. Magnetic exchange parameters and anisotropy of the quasi-two-dimensional antiferromagnet NiPS3. Phys. Rev. B 98, 134414 (2018).10.1103/PhysRevB.98.134414
31. He Y-M Single quantum emitters in monolayer semiconductors Nat. Nanotechnol. 2015 10 497 502 10.1038/nnano.2015.75 25938571
He, Y.-M. et al. Single quantum emitters in monolayer semiconductors. Nat. Nanotechnol. 10, 497–502 (2015).25938571 10.1038/nnano.2015.75
32. Srivastava A Optically active quantum dots in monolayer WSe2 Nat. Nanotechnol. 2015 10 491 496 10.1038/nnano.2015.60 25938570
Srivastava, A. et al. Optically active quantum dots in monolayer WSe2. Nat. Nanotechnol. 10, 491–496 (2015).25938570 10.1038/nnano.2015.60
33. Belvin CA Exciton-driven antiferromagnetic metal in a correlated van der Waals insulator Nat. Commun. 2021 12 4837 10.1038/s41467-021-25164-8 34376692
Belvin, C. A. et al. Exciton-driven antiferromagnetic metal in a correlated van der Waals insulator. Nat. Commun. 12, 4837 (2021).34376692 10.1038/s41467-021-25164-8
34. Afanasiev D Controlling the anisotropy of a van der Waals antiferromagnet with light Sci. Adv. 2021 7 eabf3096 10.1126/sciadv.abf3096 34078601
Afanasiev, D. et al. Controlling the anisotropy of a van der Waals antiferromagnet with light. Sci. Adv. 7, eabf3096 (2021).34078601 10.1126/sciadv.abf3096
35. Long G Persistence of magnetism in atomically thin MnPS3 crystals Nano Lett. 2020 20 2452 2459 10.1021/acs.nanolett.9b05165 32142288
Long, G. et al. Persistence of magnetism in atomically thin MnPS3 crystals. Nano Lett. 20, 2452–2459 (2020).32142288 10.1021/acs.nanolett.9b05165
36. Guo M Zhang H Cheng R Manipulating ferrimagnets by fields and currents Phys. Rev. B 2022 105 064410 10.1103/PhysRevB.105.064410
Guo, M., Zhang, H. & Cheng, R. Manipulating ferrimagnets by fields and currents. Phys. Rev. B 105, 064410 (2022).10.1103/PhysRevB.105.064410
37. Joy PA Vasudevan S Magnetism in the layered transition-metal thiophosphates MPS3 (M = Mn, Fe, and Ni) Phys. Rev. B 1992 46 5425 5433 10.1103/PhysRevB.46.5425
Joy, P. A. & Vasudevan, S. Magnetism in the layered transition-metal thiophosphates MPS3 (M = Mn, Fe, and Ni). Phys. Rev. B 46, 5425–5433 (1992).10.1103/PhysRevB.46.5425
38. Jana D Magnon gap excitations and spin-entangled optical transition in the van der Waals antiferromagnet NiPS3 Phys. Rev. B 2023 108 115149 10.1103/PhysRevB.108.115149
Jana, D. et al. Magnon gap excitations and spin-entangled optical transition in the van der Waals antiferromagnet NiPS3. Phys. Rev. B 108, 115149 (2023).10.1103/PhysRevB.108.115149
39. Mellado, P. Spin model for the honeycomb NiPS3. Appl. Phys. Lett. 123, 242403 (2023).
40. Lu Z Magnetic field mixing and splitting of bright and dark excitons in monolayer MoSe2 2D Mater. 2020 7 015017 10.1088/2053-1583/ab5614
Lu, Z. et al. Magnetic field mixing and splitting of bright and dark excitons in monolayer MoSe2. 2D Mater. 7, 015017 (2020).10.1088/2053-1583/ab5614
41. Sun Z Giant nonreciprocal second-harmonic generation from antiferromagnetic bilayer CrI3 Nature 2019 572 497 501 10.1038/s41586-019-1445-3 31367036
Sun, Z. et al. Giant nonreciprocal second-harmonic generation from antiferromagnetic bilayer CrI3. Nature 572, 497–501 (2019).31367036 10.1038/s41586-019-1445-3
42. Chen W Direct observation of van der Waals stacking–dependent interlayer magnetism Science 2019 366 983 987 10.1126/science.aav1937 31753996
Chen, W. et al. Direct observation of van der Waals stacking–dependent interlayer magnetism. Science 366, 983–987 (2019).31753996 10.1126/science.aav1937
43. Ni Z Imaging the Néel vector switching in the monolayer antiferromagnet MnPSe3 with strain-controlled Ising order Nat. Nanotechnol. 2021 16 782 787 10.1038/s41565-021-00885-5 33875873
Ni, Z. et al. Imaging the Néel vector switching in the monolayer antiferromagnet MnPSe3 with strain-controlled Ising order. Nat. Nanotechnol. 16, 782–787 (2021).33875873 10.1038/s41565-021-00885-5
