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

52220
10.1038/s41467-024-52220-w
Article
Manipulation of anisotropic Zhang-Rice exciton in NiPS3 by magnetic field
Song Feilong 1
Lv Yanpei 12
Sun Yu-Jia 12
Pang Simin 12
Chang Haonan 12
Guan Shan 1
Lai Jia-Min 12
http://orcid.org/0000-0001-9154-4115
Wang Xu-Jie 3
Wu Bang 3
Hu Chengyong 3
http://orcid.org/0000-0001-5276-9151
Yuan Zhiliang 3
http://orcid.org/0000-0002-9831-6796
Zhang Jun zhangjwill@semi.ac.cn

12
1 grid.9227.e 0000000119573309 State Key Laboratory of Superlattices and Microstructures, Institute of Semiconductors, Chinese Academy of Sciences, Beijing, China
2 https://ror.org/05qbk4x57 grid.410726.6 0000 0004 1797 8419 Center of Materials Science and Optoelectronics Engineering, University of Chinese Academy of Sciences, Beijing, China
3 https://ror.org/04nqf9k60 grid.510904.9 0000 0004 9362 2406 Beijing Academy of Quantum Information Science, Beijing, China
7 9 2024
7 9 2024
2024
15 784123 6 2023
29 8 2024
© The Author(s) 2024
2024
https://creativecommons.org/licenses/by/4.0/ 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 effect of external magnetic fields on the behavior of the Zhang-Rice exciton in NiPS3, which captures the physics of spin-orbital entanglement in 2D XY-type antiferromagnets, remains unclear. This study presents systematic study of angle-resolved and polarization-resolved magneto-optical photoluminescence spectra of NiPS3 in the Voigt geometry. We observed highly anisotropic, non-linear Zeeman splitting and polarization rotation of the Zhang-Rice exciton, which depends on the direction and intensity of the magnetic field and can be attributed to the spin-orbital coupling and field-induced spin reorientation. Furthermore, above the critical magnetic field, we detected additional splitting of the exciton peaks, indicating the coexistence of various orientations of Néel vector. This study characterizes orbital change of Zhang-Rice exciton and field-induced spin-reorientation phase transitions in a 2D hexagonal XY-type antiferromagnet, and it further demonstrates the continuous manipulation of the spin and polarization of the Zhang-Rice exciton.

NiPS3 is a van der Waals antiferromagnet with a rich optical response including an exciton of very narrow linewidth, the origin of which is still a topic of active discussion. Herein, Song, Lv and coauthors study the response of the Zhang-Rice exciton to applied magnetic fields.

Subject terms

Two-dimensional materials
Magneto-optics
issue-copyright-statement© Springer Nature Limited 2024
==== Body
pmcIntroduction

Two-dimensional (2D) magnetism, which is associated with strong intrinsic spin fluctuations, has long been the heart of fundamental questions in condensed matter physics1,2. The recent emergence of van der Waals (vdW) magnet provides a new platform to explore spin-correlated physics, such as strong correlation, quantum magnetism, topological spintronics as well as their opto-spintronic and quantum transduction applications3–5. Yet, the study of exciton coupled with the intrinsic magnetic order in 2D magnet which plays an important role in bridging the photon degree with other quasiparticles in condensed matter such as magnon, phonon, etc., remains a challenge due to the lack of suitable materials. Recently, spontaneous circularly polarized photoluminescence (PL) from the Frenkel excitons in ferromagnetic monolayer CrI3 and CrBr36,7, and linear polarization exciton strongly coupled with coherent magnons in CrSBr have been observed5. However, the polarization of the exciton is insensitive to the magnetic field direction, which means that little anisotropic information of magnetic order is reflected from the exciton8.

NiPS3, as a 2D vdW antiferromagnetic insulator, belongs to a class of transition metal phosphorus trichalcogenides (APX3, A: Fe, Mn, Ni and X:S, Se). Recent reports have shown that the X1 exciton in NiPS3 exhibits strong excitonic anisotropy in the optical linear dichroism (LD) spectrum, pump-probe Kerr rotation experiment, PL spectrum, etc9–13. Different from the CrBr3 and CrI36,7, the X1 exciton emission in NiPS3 exhibits an ultra-narrow linewidth (about 350 μeV) and highly linear polarization10,12,14. In addition, the exciton tends to couple with other quasiparticles. For example, the LD spectra show that the exciton is coupled with the A1g phonon and forms the exciton-phonon bound state10. Moreover, the NiPS3 thin flake in a photonic cavity shows that the exciton is strongly coupled with the cavity photon15. However, further applications of the exciton, such as the quantum manipulation and opto-spintronic devices, are hindered by several problems. Firstly, the origin of the X1 exciton remains controversial, like the Zhang-Rice exciton (ZRE)14, the Wannier-like exciton from the recombination of electron and hole10,16, the spin-flip induced d-d emission in the Frenkel exciton12,17, the defect-related exciton18, etc. Then, the relationships of the Néel vector and polarization of the exciton were identified variously10,12,14. Lastly, there is a lack of understanding and observation about energy shift or splitting of the X1 exciton under the magnetic field, which is powerful in determining the characteristics of the exciton such as Zeeman splitting, diamagnetic shift, g-factor, etc19,20. It is essential for further understanding of spin-orbital coupling and magnetic transition such as spin-flop transition and application in opto-spintronics3,21.

In this work, we conduct angle-resolved and polarization-resolved magneto-optical experiments of NiPS3 in the Voigt geometry. With applying an in-plane magnetic field (μ0H), we observe the Zeeman splitting of the linear polarized PL of the X1 exciton. With rotating the sample in the ab-plane and changing the μ0H, we observe that the Zeeman splitting is highly anisotropic and non-linearly dependent on the μ0H. With the ZRE model, we attribute the anisotropy and non-linearity to magnetic field-induced spin reorientation transition (SRT) and read out the Néel vector in NiPS3. With measuring the polarization of splitting PL peaks, we observe the identical anisotropic behavior of the polarization and confirm the spin-polarization coupling during the SRT process, which reflects the localism and orbital change of ZRE and provide an explanation for the short lifetime of ZRE. We accomplish the accurate determination of the spin orientation based on the polarization and splitting energy of the ZRE. Based on this spin read-out technique, we observe the coexistence of various orientations of Néel vector near the critical magnetic field approximately parallel with the a-axis. Our work demonstrates that the ZRE in NiPS3 is an ideal platform for the study of the spin-orbital entanglement and the magnetic phase transition in the XY-type antiferromagnets. Besides, we also provide a method to continuously manipulate and detect the spin and the polarization of the ZRE, which enables the potential applications of antiferromagnets in quantum information and opto-spintronic devices.

Results and discussion

Bulk NiPS3 has a monoclinic structure with point group C2h, whereas monolayer NiPS3 has a trigonal structure with point group D3d22,23. Figure 1a indicates the crystalline axis of the honeycomb lattice structure with three-fold rotational symmetry in monolayer NiPS3. Below the Néel temperature (TN = 155 K), NiPS3 becomes an antiferromagnetic insulator10,16, in which the Néel vector is confirmed almost aligned in-plane and along the a-axis24. The interlayer coupling in NiPS3 is ferromagnetic, as shown in Fig. 1b. X-ray absorption studies suggest that NiPS3 is a charge-transfer antiferromagnetic insulator that exhibits strong spin-charge correlation effects below the Néel temperature25.Fig. 1 Crystal structure and antiferromagnetic exciton state of NiPS3.

a Schematic of the spin orientation in a single layer of NiPS3. b Monoclinic stacking with ferromagnetic interlayer coupling in a bulk crystal. c Temperature-dependent reflection contrast spectra. The red curve denotes the reflection contrast spectra below the Néel temperature, while the blue curve is above Néel temperature. These peaks in the reflection contrast spectra are assigned to the X1 and X3 exciton. d Polarized PL spectra collected from the parallel and perpendicular directions to PL polarization, the linewidth of the PL spectra is about 260 μeV. The red line is the PL spectrum collected parallel to the PL polarization, while the blue line is perpendicular to the PL polarization. In the spectra, there is a main peak position at 1.476 eV (labeled as X1) and a shoulder peak position at 1.479 eV (labeled as X2).

With the experimental setup as shown in Method and SI Fig. S1, we measured the temperature-dependent and magnetic field-dependent spectra. Figure 1c is the temperature-dependent reflection contrast spectra of the NiPS3 thin flakes, in which there are two sharp peaks at 1.476 eV and 1.498 eV below the Néel temperature. The low-energy one (labeled as X1) corresponds to the reported spin-correlated exciton. The high-energy one (labeled as X3) was earlier attributed to the two-magnon sideband or Rydberg exciton of the X114,16, while subsequent temperature-dependent experiments identified the X3 as the phonon sideband of the X112,23,26. The temperature dependence of the reflection contrast spectra indicates that these two peaks are coupled to the long-range magnetic order as the two peaks get considerably weakened and disappear above 150 K when the temperature increases. And the temperature-dependent polarization-resolved PL spectra (as shown in SI Fig. S2) also indicate similar magnetic behavior with 2D XY-type systems12,27,28. Figure 1d shows the polarization-resolved PL spectra of NiPS3 thin flakes at 4 K, from which we can easily find that the PL is almost linearly polarized. The sharp PL peak (labeled as X1) is observed at 1.4758 eV with a side peak (labeled as X2) at 1.4788 eV. As shown in SI Fig. S3, by varying the intensity of the laser and the sample position, the PL intensity of X1 is linearly dependent on the laser power and the profile of the X1’s PL is insensitive to the sample position. Hence, the defect-related origin of X1 could be excluded. The X1 exciton shows an ultra-narrow linewidth in our experiments, about 260 μeV (Fig. 1d), similar to those of other groups10,12,14, which indicates the good crystallization of our sample and the suitability for studying the exciton properties with the magneto-optical experiment29,30.

Up to now, number of phenomena of this sharp exciton emission have been fully researched with magneto-optical experiments12,31,32, except for the Zeeman splitting. The magnetic field will induce the Zeeman splitting of the excitons in traditional semiconductors when its component along the spin direction of the excitons is nonzero30,33,34. Hence, we apply a magnetic field in the ab-plane of NiPS3 to detect the magneto-optical properties of X1. We observe that the X1 exciton splits into two peaks when μ0H is parallel to the a-axis, as shown in Fig. 2a. The splitting energy is linearly dependent on the μ0H with a slope of 0.22 meV/T, as shown in SI Figure S4. When the magnetic field is increased from 0 T to 9 T, the PL peaks of ZRE keep linearly polarized and the polarization keeps parallel with the external field, as shown in the Fig. 2c. This polarization behavior indicates that the PL of ZRE is a π light which corresponds to the unchanged total angular momentum from the excited state to the ground state35,36. The Zeeman splitting of the X1 can be explained by a model as shown in Fig. 2b. The ground state (Zhang-Rice triplet, ZRT) has a high-spin state of S = 1 while the spin of the Zhang-Rice singlet (ZRS) is zero14. The π light of ZRE means the contribution of spin-orbital coupling and the change in orbitals from ZRS to ZRE. As a result, the difference of gyromagnetic ratio (Δg) between the ZRT and ZRE leads the energy of ZRE in different spin chains to shift in an opposite direction under the magnetic field along the a-axis, as discussed in Methods, SI Note 4, and SI Fig. S6. The orbital change means that the spin-correlated exciton does not originate from d-d transition and the exciton could obey optical selection rules, which could provide potential explanation for the short lifetime of ZRE10,12, as discussed in the SI Note 5. The a-axis could be determined with an angle-resolved measurement as discussed in Methods and SI Note 1.Fig. 2 Zeeman splitting of the excitons in NiPS3 thin flakes.

a Magnetic field-dependent PL spectra and Zeeman splitting of ZRS exciton in NiPS3 thin flakes. b A schematic view of the Zeeman splitting of the ZRE under the magnetic field parallel to the a-axis. c Contour plot spectra of the polarized PL spectra with a magnetic field at 0 T and 9 T (left panel and right panel respectively). The direction of the magnetic field is parallel to a-axis.

Since the ZRE is coupled with magnetic order and its PL spectrum exhibits an ultra-narrow linewidth, the Zeeman splitting could accurately indicate the information of the Néel vector. Similar protocols to detect and manipulate Néel vector have been accomplished with the Zeeman effect of magnons in bulk crystals such as MnF2 and Cr2O321,37. As shown in Fig. 3a, we systematically measured the magnetic field-dependent PL spectra of NiPS3 thin flakes. The θ denotes the angle between μ0H and the a-axis. In the first panel (θ = 4°) of the Fig. 3a, a typical Zeeman splitting phenomenon with the magnetic field below 10 T can be observed, where the peaks of the exciton linearly split with μ0H. When the θ increases from 4° to 86°, the Zeeman splitting shows highly anisotropic and nonlinear dependences on the μ0H, which could be illustrated by the ZRE model as discussed above. The anti-paralleled magnetic moments of ZRE have different Zeeman energies±ΔgμBμ0Hcosφ, where φ represents the angle between the μ0H and the Néel vector. The nonlinearity and anisotropy of Zeeman splitting indicate that φ depends on μ0H and θ. In other word, a direct conclusion from the spectra based on the ZRE model is that the Néel vector will rotate to a stable orientation under the in-plane magnetic field. The similar spin-rotating phenomenon, which is called as the magnetic field-induced spin reorientation transition (SRT), have been reported by Néel and Yosida in the research of the tetragonal MnF238,39.Fig. 3 Angle-dependent Zeeman splitting of ZRE.

a Contour plot of the magnetic field-dependent PL spectra of ZRE under different θ. θ is the angle between the a-axis and the magnetic field. b The Zeeman splitting energy as a function of the angle θ and magnetic field. The dots are the experimental data, while the lines are the corresponding fitted results of the same-colored data. c Magnetic field-dependent Néel vector rotation. The dots are calculated from the experimental data, while the lines are the corresponding fitted results. The inset cartoon schematically illustrates the angle φ between magnetic field and the Néel vector L.

To further research the dependence of Zeeman splitting on the μ0H, we extracted the Zeeman splitting from the spectra in Fig. 3a and fitted them with a formula 2ΔgμBμ0Hcosθ−12arctansin2θcos2θ−Hc2H2, as shown in Fig. 3b, where μ0Hc is selected as 10.75 T, consistent with the reported experimental result12. We should note here that the fit formula, which could be directly obtained from our data as shown in Methods, is just the same as Néel’s formula in the discussion of MnF239. In Fig. 3c, the φ−θ represents the rotation angle of the Néel vector, which is extracted from the Fig. 3b and fitted with φ=θ−12arctansin2θcos2θ−Hc2H2. With Monte Carlo simulations for a 2D hexagonal XY-type antiferromagnet (as shown in Methods and SI Fig. S5), the simulated Zeeman splitting and φ are the same as the curves calculated with Néel’s formula in Fig. 3 b and c, reflecting the anisotropy akin to a tetragonal lattice under a magnetic field for the zig-zag antiferromagnets. Due to the narrow linewidth of NiPS3, we obtain small errors in the extracted energies and the angle φ, which indicate the potential application of ZRE for accurate spin readout. With the μ0H ranging from 0 T to 12 T, we realize the manipulation of the Néel vector to almost any direction in the ab-plane of NiPS3, as shown in Fig. 3c.

Furthermore, we varied the μ0H and θ to measure the polarization of ZRE to study the anisotropic behavior of it. To research the rotation of polarization, we defined the polarization-rotation angle β, as shown in Fig. 4a. Figure 4b shows that the polarization of ZRE is unchanged below the critical magnetic field and rotates about 90° above the critical magnetic field, when θ is 0°. Figure 4c shows that the polarization continuously rotates when θ is 43° and Fig. 4d shows that nearly no polarization rotates when θ is 86°, when the magnetic field increases from 0 T to 12 T. In Fig. 4e, we extract the polarization angle for different θ under the μ0H, increasing from 0 T to 12 T. We find that the polarization angle can be described by the same formula that characterizes the reorientation of the Néel vector under the same μ0H, as shown in Fig. 4e, indicating the spin-polarization coupling of ZRE during the SRT process. Figure 4f is the comparison of the extracted β and the φ−θ at 9 T and 12 T, which confirms that the rotation of the Néel vector and polarization is identical under the magnetic field below and above the μ0Hc. Interestingly, the distorted sine function in Fig. 4f at 9 T is similar to that in magnetic transport experiments40.Fig. 4 Anisotropic rotation of the polarization of ZRE under the magnetic field from 0 T to 12 T.

a The β is defined as rotation angle of polarization in ZRE with and without magnetic field. b–d The rotation of the polarization of ZRE at a magnetic field of 0 T, 9 T and 12 T under θ is about 0°, 45°, and 90°. e The scattering points are measured polarizations and the curves are the same as curves in Fig. 3c. f Anisotropic rotation of the polarization of ZRE and the Néel vector when the magnetic field is 9 T and 12 T. The red and blue small balls represent the data of φ−θ and β respectively. The dash and solid line are calculated from Néel’s formula.

To research the behaviors of the SRT at the critical magnetic field, we conducted magnetic field-dependent polarization measurements for ZRE when θ = 1°, as shown in Fig. 5a. When the μ0H increases from 0 T to 9.5 T, the exciton peaks split into two branches without polarization rotation. However, when the magnetic field increases to 10 T, the energy and polarization of the exciton peaks that existed below 9.5 T remain almost unchanged, while two new peaks appear with smaller splitting energy and polarizations nearly perpendicular to the original peaks. For clarity, the new emerging peaks are indicated by green arrows (dots and curves) and pre-existing peaks are indicated by red arrows (dots and curves), as shown in Fig. 5a (Fig. 5 b–e). As shown in the Fig. 5a, d, and e, when increasing the magnetic field from 10 T to 12 T, the intensity of the pre-existing peaks become weaker, while the intensity of the new emerging peaks becomes stronger, indicating the transformation of different spin configurations. Besides, the polarization and energy of the pre-existing peaks gradually get closer towards the new emerging peaks, while the polarization and energy of the new emerging peaks remain essentially unchanged, indicating that the spin-polarization coupling is still maintained. The multiple peaks reflect the change of magnetic order in the transition region near the critical field. This phenomenon could originate from inhomogeneity-induced various spin-flop critical fields41 or the transition to an intermediate phase42. If the phenomenon originates from inhomogeneity, our findings can help to further investigate the origin of inhomogeneity and thus domain engineering43 in antiferromagnetic materials. If the phenomenon originates from the intermediate phase, the neutron scattering could help to further investigation42 and the magnetic field dependence of the peak positions and polarization in Fig. 5 could contribute to the study of magnetic phase transitions in XXZ-type or XY-type antiferromagnets (AFM)23. Further determination of the origin needs more experimental evidence and is beyond the scope of our discussion. Our results show that the PL polarization and energy of ZRE is a non-destructive and convenient method to detect the critical behavior of NiPS3.Fig. 5 Field-induced evolution of multiple spin configurations of ZRE.

a Polarization-resolved PL spectra of ZRE when θ is approximately 1°. The four arrows represent the different spin orientations extracted from the four exciton peaks. Red arrows represent the pre-existing peaks, while green arrows represent the new emerging peaks. b–e Polarizations extracted from the four peaks. Dots are extracted polarization from experimental data, while the lines are fitting curves. Red dots represent the pre-existing peaks, while green dots represent the new emerging peaks.

From a systematically angle-resolved and polarization-resolved magneto-optical PL spectra measurement, we observed non-linear and anisotropic Zeeman splitting originated from magnetic field-induced SRT. Mapping the energy splitting to the direction of the Néel vector, we confirmed that the anisotropic linearly polarized PL emission of ZRE in NiPS3 is coupled with Néel vector when varying the in-plane magnetic field, which indicates the spin-orbital coupling and orbital change from ZRS to ZRT. With the advantage of the ultra-narrow linewidth, we found that the orientation of the Néel vector could be continuously read out and manipulated by an external in-plane magnetic field, which opens a path for detecting and manipulating the spin-orbit-entangled exciton state in a 2D magnetic material and promotes the development of opto-spintronic devices and magnetic quantum information technology. With this spin-readout technique, we probe the coexistence of various orientations of the Néel vector, which could be attributed to the inhomogeneity or the transition to an intermediate phase. Our findings provide new insights for the localism, orbital change, and short lifetime of ZRE, as well as AFM phase transition in NiPS3. We thus foresee that NiPS3 is a promising 2D platform for exploring magnetic phase transition and strongly correlated systems with the optical methods. After this article was submitted, we became aware of another paper44 reporting related to the Zeeman splitting of NiPS3.

Methods

Sample fabrication and PL measurement

The sample in our experiment was mechanically exfoliated from the high-quality NiPS3 bulk single crystals grown by chemical vapor transport method. It was then transferred to a marked SiO2/Si substrate. A He-Ne laser (633 nm) was used to excite the NiPS3 sample, which was placed in a closed cycle cryostat (attodry 2100). The PL spectra of the exciton in NiPS3 were polarization-resolved and collected by a spectrometer with a liquid nitrogen-cooled CCD. A magnetic field up to 12 T was applied in Voigt geometry. A low-temperature objective was used to focus the excitation laser onto the sample and collect the PL signal. An electronically controlled rotary stage enables us to rotate the sample in a low-temperature and high-magnetic-field environment. The magneto-optical experiments were performed in Voigt geometry. The sample was positioned vertically inside the magnetic cell with the surface parallel to the applied magnetic field for the Voigt geometry. A mirror was set between the objective and sample with an angle of 45° to change the optical path by 90°. For both the polarization-dependent and magnetic-field-dependent experiments, the PL signal was filtered by an analyzer and then collected and measured by a spectrometer with a CCD. Two samples were measured in our experiments, as shown in SI Note 1.

Experimental details

The images of the sample inside the magnetic cell were captured with a camera, as shown in the SI Fig. S7. The polarization-resolved spectra, when θ ranges from -30° to 90° (SI Fig. S8), show that the Zeeman splitting when θ=0° is the largest. Combining the images, the Zeeman splitting, and the polarization, we could determine a-axis as discussed in the SI Note 1. By flipping the magnetic field direction, the Faraday rotation effect could be excluded in the research of polarization rotation, as shown in the SI Fig. S9. Before and after flipping the μ0H, the behaviors of SRT were consistent when the μ0H is higher than the μ0Hc, as shown in Fig. S10.

ZRE model and data analysis

Although the total angular momentum of the ground state and excited state is the same (JZRT=JZRS≡J0), there is a distinction in the quantum numbers of spin and orbital angular momentum between the ground and excited states, as discussed in the SI Note 4. Consequently, the effective gyromagnetic ratios of the ground state and excited state are typically different, that is, gZRT≠gZRS. In antiferromagnetic materials, the localized magnetic moment is represented by μBgJ0 or −μBgJ0, where μB denotes Bohr magneton. Under the external magnetic field B, the exciton ground state energy is given by EZRTB=EZRT0±μBgZRTJZRT⋅B, while the corresponding exciton excited state energy is EZRSB=EZRS0±μBgZRSJZRS⋅B. Therefore, the peak positions of the ZRE is EZREB=EZRE0±μBgZRT−gZRSJ0⋅B=EZRE0±μBΔgJ0⋅B. In the external magnetic field, the “spin reorientation” means that the localized magnetic moment (±μBgJ0) rotates in the plane, leading to anisotropic and nonlinear Zeeman effects, as stated in the original version. The peak positions in Fig. 3b are extracted from PL spectra in Fig. 3a with peaks fitted by the Lorentzian function. The angles (φ−θ) are extracted from the splitting energy in Fig. 3b with formula φ−θ=arccosΔEθ,μ0H2gμB, whose errors include one part propagated from errors of peak positions and another part resulted from angle measurement in experiments. The fit curves in Fig. 3b and Fig. 3c are both plotted based on the following steps. Firstly, as shown in SI Fig. S11a or c, we plot the dependence of φ−θ on θ under different μ0H. The dependence is well fitted with formula φ−θ = −12arctan(sin2θcos2θ−A) and the fitting parameter A are various in different μ0H. Then, as shown in SI Fig. S11b or d, we plot the dependence of A on μ0H, which is fitted with A=Hc2H2 and Hc=10.75T. As a result, we could plot the dependence of φ−θ on μ0H for any θ. Finally, combined with the Zeeman splitting energy ±gμBμ0Hcosφ, we could plot the dependence of splitting energy on μ0H for any θ. The formula φ−θ =12arctansin2θcos2θ−H2Hc2−θ, obtained from our data here, is different from Néel’s formula, while it is consistent in physics for the C2 symmetry of NiPS3. More discussions could refer to SI Note 2.

Monte Carlo simulations

We calculated the zig-zag antiferromagnetic spin evolution at 1.6 K by Monte Carlo simulations. We treated the spins classically and included the short-range exchange interactions up to third neighbors with ignoring the anisotropic exchange interaction. At each magnetic field, total steps are set as 30000. For all the numerical data from Monte Carlo simulations, the numerical uncertainties are smaller than or comparable to the size of symbols. More details are indicated in SI Note 3.

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-52220-w.

Acknowledgements

We are grateful for high-field measurements supported by Xiamen University (X. Wang) and helpful discussion about the theoretical calculations with University of Science and Technology of China (J. Zhao). Authors acknowledge the funding support from Chinese Academy of Sciences - the Scientific and Technological Research council of TÜRKİYE Joint Research Projects (172111KYSB20210004), the CAS Interdisciplinary Innovation Team, National Natural Science Foundation of China (12074371), Research Equipment Development Project of Chinese Academy of Sciences (YJKYYQ20210001).

Author contributions

J.Z. supervised the project; J.Z. and F.S. conceived the ideas; F.S. and Y.L. prepared the samples and carried out the experiments, assisted by H.C., Y.S., X.W. and B.W. F.S., Y.L. and J.Z. analyzed the data. Y.L. fitted the anisotropic Zeeman splitting data. S.G. and Y.L. conducted the Monte Carlo simulations. F.S., Y.L. and J.Z. wrote the paper with input from all authors. All authors discussed the results and commented on the manuscript.

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 authors declare that data generated in this study are provided in the paper and the Supplementary Information file. Further datasets are available from the corresponding author upon request. Source data are provided in this paper.

Code availability

The authors declare that code supporting the findings of this study is available within the paper and the Supplementary Information file. Further code is available from the corresponding author upon request.

Competing interests

The authors declare that they have no competing interests.

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

These authors contributed equally: Feilong Song, Yanpei Lv.
==== Refs
References

1. 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
2. Mermin ND Wagner H Absence of ferromagnetism or antiferromagnetism in one- or two-dimensional isotropic Heisenberg models Phys. Rev. Lett. 1966 17 1133 1136 10.1103/PhysRevLett.17.1133
Mermin, N. D. & Wagner, H. Absence of ferromagnetism or antiferromagnetism in one- or two-dimensional isotropic Heisenberg models. Phys. Rev. Lett. 17, 1133–1136 (1966).10.1103/PhysRevLett.17.1133
3. Diederich GM Tunable interaction between excitons and hybridized magnons in a layered semiconductor Nat. Nanotechnol. 2023 18 23 28 10.1038/s41565-022-01259-1 36577852
Diederich, G. M. et al. Tunable interaction between excitons and hybridized magnons in a layered semiconductor. Nat. Nanotechnol. 18, 23–28 (2023).36577852 10.1038/s41565-022-01259-1
4. 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
5. 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
6. Seyler KL Ligand-field helical luminescence in a 2D ferromagnetic insulator Nat. Phys. 2018 14 277 281 10.1038/s41567-017-0006-7
Seyler, K. L. et al. Ligand-field helical luminescence in a 2D ferromagnetic insulator. Nat. Phys. 14, 277–281 (2018).10.1038/s41567-017-0006-7
7. Zhang Z Direct photoluminescence probing of ferromagnetism in monolayer two-dimensional CrBr3 Nano. Lett. 2019 19 3138 3142 10.1021/acs.nanolett.9b00553 30945874
Zhang, Z. et al. Direct photoluminescence probing of ferromagnetism in monolayer two-dimensional CrBr3. Nano. Lett. 19, 3138–3142 (2019).30945874 10.1021/acs.nanolett.9b00553
8. Wilson NP Interlayer electronic coupling on demand in a 2D magnetic semiconductor Nat. Mater. 2021 20 1657 1662 10.1038/s41563-021-01070-8 34312534
Wilson, N. P. et al. Interlayer electronic coupling on demand in a 2D magnetic semiconductor. Nat. Mater. 20, 1657–1662 (2021).34312534 10.1038/s41563-021-01070-8
9. Zhang Q Observation of giant optical linear dichroism in a zigzag antiferromagnet FePS3 Nano. Lett. 2021 21 6938 6945 10.1021/acs.nanolett.1c02188 34428905
Zhang, Q. et al. Observation of giant optical linear dichroism in a zigzag antiferromagnet FePS3. Nano. Lett. 21, 6938–6945 (2021).34428905 10.1021/acs.nanolett.1c02188
10. 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
11. Afanasiev D Controlling the anisotropy of a van der Waals antiferromagnet with light Sci. Adv. 2021 7 130 10.1126/sciadv.abf3096
Afanasiev, D. et al. Controlling the anisotropy of a van der Waals antiferromagnet with light. Sci. Adv. 7, 130 (2021).10.1126/sciadv.abf3096
12. 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
13. Tan Q Observation of three-state nematicity and domain evolution in atomically thin antiferromagnetic NiPS3 Nano Lett. 2024 24 7166 7172 10.1021/acs.nanolett.4c00772
Tan, Q. et al. Observation of three-state nematicity and domain evolution in atomically thin antiferromagnetic NiPS3. Nano Lett. 24, 7166–7172 (2024).10.1021/acs.nanolett.4c00772
14. 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
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. Ho C-H Hsu T-Y Muhimmah LC The band-edge excitons observed in few-layer NiPS3 npj 2D Mater. Appl. 2021 5 8 10.1038/s41699-020-00188-8
Ho, C.-H., Hsu, T.-Y. & Muhimmah, L. C. The band-edge excitons observed in few-layer NiPS3. npj 2D Mater. Appl. 5, 8 (2021).10.1038/s41699-020-00188-8
17. Gnatchenko SL Kachur IS Piryatinskaya VG Vysochanskii YM Gurzan MI Exciton-magnon structure of the optical absorption spectrum of antiferromagnetic MnPS3 Low. Temp. Phys. 2011 37 144 148 10.1063/1.3556660
Gnatchenko, S. L., Kachur, I. S., Piryatinskaya, V. G., Vysochanskii, Y. M. & Gurzan, M. I. Exciton-magnon structure of the optical absorption spectrum of antiferromagnetic MnPS3. Low. Temp. Phys. 37, 144–148 (2011).10.1063/1.3556660
18. 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
19. Wu S Anisotropies of the g-factor tensor and diamagnetic coefficient in crystal-phase quantum dots in InP nanowires Nano Res. 2019 12 2842 2848 10.1007/s12274-019-2522-5
Wu, S. et al. Anisotropies of the g-factor tensor and diamagnetic coefficient in crystal-phase quantum dots in InP nanowires. Nano Res. 12, 2842–2848 (2019).10.1007/s12274-019-2522-5
20. van Bree J Anisotropy of electron and hole g tensors of quantum dots: an intuitive picture based on spin-correlated orbital currents Phys. Rev. B 2016 93 035311 10.1103/PhysRevB.93.035311
van Bree, J. et al. Anisotropy of electron and hole g tensors of quantum dots: an intuitive picture based on spin-correlated orbital currents. Phys. Rev. B 93, 035311 (2016).10.1103/PhysRevB.93.035311
21. Vaidya P Subterahertz spin pumping from an insulating antiferromagnet Science 2020 368 160 165 10.1126/science.aaz4247 32273462
Vaidya, P. et al. Subterahertz spin pumping from an insulating antiferromagnet. Science 368, 160–165 (2020).32273462 10.1126/science.aaz4247
22. Susner MA Chyasnavichyus M McGuire MA Ganesh P Maksymovych P Metal thio- and selenophosphates as multifunctional van der Waals layered materials Adv. Mater. 2017 29 1602852 10.1002/adma.201602852
Susner, M. A., Chyasnavichyus, M., McGuire, M. A., Ganesh, P. & Maksymovych, P. Metal thio- and selenophosphates as multifunctional van der Waals layered materials. Adv. Mater. 29, 1602852 (2017).10.1002/adma.201602852
23. 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
24. 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
25. 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
26. 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
27. Bramwell ST Holdsworth PCW Magnetization and universal sub-critical behaviour in two-dimensional XY magnets J. Phys. Condens. Matter 1993 5 L53 10.1088/0953-8984/5/4/004
Bramwell, S. T. & Holdsworth, P. C. W. Magnetization and universal sub-critical behaviour in two-dimensional XY magnets. J. Phys. Condens. Matter 5, L53 (1993).10.1088/0953-8984/5/4/004
28. Hu L Observation of a magnetic phase transition in monolayer NiPS3 Phys. Rev. B 2023 107 L220407 10.1103/PhysRevB.107.L220407
Hu, L. et al. Observation of a magnetic phase transition in monolayer NiPS3. Phys. Rev. B 107, L220407 (2023).10.1103/PhysRevB.107.L220407
29. Kopteva NE Zeeman and Davydov splitting of Frenkel excitons in the antiferromagnet CuB2O4 Phys. Rev. B 2022 105 024421 10.1103/PhysRevB.105.024421
Kopteva, N. E. et al. Zeeman and Davydov splitting of Frenkel excitons in the antiferromagnet CuB2O4. Phys. Rev. B 105, 024421 (2022).10.1103/PhysRevB.105.024421
30. Gantz L Controlling the dark exciton spin eigenstates by external magnetic field Phys. Rev. B 2016 94 045426 10.1103/PhysRevB.94.045426
Gantz, L. et al. Controlling the dark exciton spin eigenstates by external magnetic field. Phys. Rev. B 94, 045426 (2016).10.1103/PhysRevB.94.045426
31. Basnet R Wegner A Pandey K Storment S Hu J Highly sensitive spin-flop transition in antiferromagnetic van der Waals material MPS3(M=Ni and Mn) Phys. Rev. Mater. 2021 5 064413 10.1103/PhysRevMaterials.5.064413
Basnet, R., Wegner, A., Pandey, K., Storment, S. & Hu, J. Highly sensitive spin-flop transition in antiferromagnetic van der Waals material MPS3(M=Ni and Mn). Phys. Rev. Mater. 5, 064413 (2021).10.1103/PhysRevMaterials.5.064413
32. Bougamha F 31P NMR investigation of quasi-two-dimensional magnetic correlations in T2P2S6 (T=Mn, Ni) Phys. Rev. B 2022 105 024410 10.1103/PhysRevB.105.024410
Bougamha, F. et al. 31P NMR investigation of quasi-two-dimensional magnetic correlations in T2P2S6 (T=Mn, Ni). Phys. Rev. B 105, 024410 (2022).10.1103/PhysRevB.105.024410
33. Lyons TP Giant effective Zeeman splitting in a monolayer semiconductor realized by spin-selective strong light–matter coupling Nat. Photon. 2022 16 632 636 10.1038/s41566-022-01025-8
Lyons, T. P. et al. Giant effective Zeeman splitting in a monolayer semiconductor realized by spin-selective strong light–matter coupling. Nat. Photon. 16, 632–636 (2022).10.1038/s41566-022-01025-8
34. Yu JH Giant Zeeman splitting in nucleation-controlled doped CdSe:Mn2+ quantum nanoribbons Nat. Mater. 2010 9 47 53 10.1038/nmat2572 19915554
Yu, J. H. et al. Giant Zeeman splitting in nucleation-controlled doped CdSe:Mn2+ quantum nanoribbons. Nat. Mater. 9, 47–53 (2010).19915554 10.1038/nmat2572
35. Wang G In-plane propagation of light in transition metal dichalcogenide monolayers: optical selection Rules Phys. Rev. Lett. 2017 119 047401 10.1103/PhysRevLett.119.047401 29341750
Wang, G. et al. In-plane propagation of light in transition metal dichalcogenide monolayers: optical selection Rules. Phys. Rev. Lett. 119, 047401 (2017).29341750 10.1103/PhysRevLett.119.047401
36. Andelkovic Z Laser cooling of externally produced Mg ions in a penning trap for sympathetic cooling of highly charged ions Phys. Rev. A 2013 87 033423 10.1103/PhysRevA.87.033423
Andelkovic, Z. et al. Laser cooling of externally produced Mg ions in a penning trap for sympathetic cooling of highly charged ions. Phys. Rev. A 87, 033423 (2013).10.1103/PhysRevA.87.033423
37. Li J Spin current from sub-terahertz-generated antiferromagnetic magnons Nature 2020 578 70 74 10.1038/s41586-020-1950-4 31988510
Li, J. et al. Spin current from sub-terahertz-generated antiferromagnetic magnons. Nature 578, 70–74 (2020).31988510 10.1038/s41586-020-1950-4
38. Yosida K On the antiferromagnetism of single crystals Prog. Theor. Phys. 1951 6 691 701 10.1143/ptp/6.5.691
Yosida, K. On the antiferromagnetism of single crystals. Prog. Theor. Phys. 6, 691–701 (1951).10.1143/ptp/6.5.691
39. Néel L Propriétés magnétiques de l'état métallique et énergie d’interaction entre atomes magnétiques Ann. Phys. 1936 11 232 279 10.1051/anphys/193611050232
Néel, L. Propriétés magnétiques de l'état métallique et énergie d’interaction entre atomes magnétiques. Ann. Phys. 11, 232–279 (1936).10.1051/anphys/193611050232
40. Fina I Anisotropic magnetoresistance in an antiferromagnetic semiconductor Nat. Commun. 2014 5 4671 10.1038/ncomms5671 25204755
Fina, I. et al. Anisotropic magnetoresistance in an antiferromagnetic semiconductor. Nat. Commun. 5, 4671 (2014).25204755 10.1038/ncomms5671
41. Stöhr, J. & Siegmann, H. C. Magnetism. Solid-State Sciences. Springer, Berlin, Heidelberg 11, 479–520 (2006).
42. Smeets JPM Frikkee E de Jonge WJM Intermediate phase in a spin-flop system with coupled order parameters Phys. Rev. Lett. 1982 49 1515 1518 10.1103/PhysRevLett.49.1515
Smeets, J. P. M., Frikkee, E. & de Jonge, W. J. M. Intermediate phase in a spin-flop system with coupled order parameters. Phys. Rev. Lett. 49, 1515–1518 (1982).10.1103/PhysRevLett.49.1515
43. Thong HC Domain engineering in bulk ferroelectric ceramics via mesoscopic chemical inhomogeneity Adv. Sci. 2022 9 2200998 10.1002/advs.202200998
Thong, H. C. et al. Domain engineering in bulk ferroelectric ceramics via mesoscopic chemical inhomogeneity. Adv. Sci. 9, 2200998 (2022).10.1002/advs.202200998
44. 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
