## Quasi-Two-Dimensional Magnon Identification in Antiferromagnetic FePS3 via Magneto-Raman Spectroscopy

Amber McCreary, 1  Jeffrey R. Simpson, 1,2 Thuc T. Mai, 1,3 Robert D. McMichael, 1 Jason E. Douglas, 4 Nicholas Butch,  Cindi Dennis,  Rolando Valdes Aguilar,  and Angela R. Hight 5 4 3 Walker 1,*

1 Nanoscale Device Characterization Division, Physical Measurement Laboratory, National Institute of Standards and Technology, Gaithersburg, Maryland 20899, USA

2 Department of Physics, Astronomy, and Geosciences, Towson University, Towson, MD 21252, USA

3 Department of Physics, The Ohio State University, Columbus, OH 43210, USA

4 Materials Science and Engineering Division, Material Measurement Laboratory, National Institute of Standards and Technology, Gaithersburg, Maryland 20899, USA

5 NIST Center for Neutron Research, National Institute of Standards and Technology,

Gaithersburg, Maryland 20899, USA

*Email:  angela.hightwalker@nist.gov

KEYWORDS:    Magnon,  spin-wave,  Raman  spectroscopy,  magneto-Raman,  phonon,  2D materials, FePS3, Mott insulator, Ising antiferromagnet

## ABSTRACT:

Recently it was discovered that van der Waals-bonded magnetic materials retain long range magnetic ordering down to a single layer, opening many avenues in fundamental physics and potential applications of these fascinating materials. One such material is FePS3, a large spin (S=2) Mott insulator where the Fe atoms form a honeycomb lattice. In the bulk, FePS3 has been shown to be a quasi-two-dimensional-Ising antiferromagnet, with additional features in the Raman spectra emerging below the Néel temperature ( T N) of approximately 120 K. Using magneto-Raman spectroscopy as an optical probe of magnetic structure, we show that one of these Raman-active modes in the magnetically ordered state is actually a magnon with a frequency of ≈3.7 THz (122 cm -1 ). Contrary to previous work, which interpreted this feature as a phonon, our Raman data shows the expected frequency shifting and splitting of the magnon as a function of temperature and magnetic field, respectively, where we determine the g-factor to be ≈2. In addition, the symmetry behavior of the magnon is studied by polarization-dependent Raman spectroscopy and explained using the magnetic point group of FePS3.

## I. INTRODUCTION

Since the isolation of monolayer graphene in 2004,  there has been a surge of research 1 into van der Waals layered materials, where the strong intralayer and weak interlayer coupling allows for isolation of layers that are only a few atoms thick. These materials exhibit a wide range of electronic properties, including semiconducting, metallic, insulating, superconducting, and charge density waves, 2-4 allowing for device architectures composed of solely twodimensional (2D) materials. While significant research has been dedicated thus far to studying the optical, mechanical, and electrical properties of 2D materials, 3,5-9 exploring magnetism is still

in its infancy, even though 2D magnetic materials provide a solid-state platform to experimentally access fundamental, low-dimensional physics. 10,11 Additionally, any 2D magnetic material would likely still possess the captivating properties of 2D materials, including extremely large mechanical flexibility, 12,13 efficient tuning of transport properties with an electric field, 14-17 relative ease of chemical modification, 18,19 as well as the ability to create van der Waals stacked heterostructures. 20 These myriad of tuning parameters could unlock opportunities for customengineered magnetoelectric and magneto-optical devices, where 2D magnets coupled with other technologically relevant materials could realize unprecedented capabilities in fields such as spintronics. 11,21,22

Until recently, long range magnetic order was believed to be impossible in 2D systems due to thermal fluctuations. 23 In early 2017, however, intrinsic ferromagnetism was observed down to the few-layer and monolayer limit in two different, layered materials with magnetic anisotropy, including Cr2Ge2Te6 and CrI3. 24,25 For CrI3, it was also shown that the interlayer magnetic ordering ( i.e . ferromagnetically or antiferromagnetically coupled) was dependent on the number of layers, 25 and could be controlled by an external electric field. 14,26 Transition metal phosphorus trisulfides X PS3 ( X = Fe, Mn, Ni, etc .) are another class of van der Waals antiferromagnets that are being studied in the 2D limit. Interestingly, although FePS3 ( T N ≈ 120 K), 27 MnPS3 ( T N ≈ 78 K), 27 and NiPS3 ( T N ≈ 155 K) 27 are isostructural, they have different spin structures below the Néel temperature. The choice of transition metal results in varied magnetic phenomena, since the spins align antiferromagnetically in different fashions, including Néel, zigzag, or stripe ordering. 28,29 FePS3, which is a Mott insulator, 30,31 is especially intriguing as a 2D Ising antiferromagnet on a honeycomb lattice. 27,32 In addition, long-distance magnon

transport (several micrometers) has been experimentally observed in MnPS3, demonstrating that these materials are viable candidates for future magnonic devices. 22

Raman spectroscopy, being non-destructive and highly sensitive to minute perturbations, is a powerful technique to study various properties of quantum materials, including the effects of layer number, 33,34 strain, 35 defects/doping, 36 electron-phonon coupling, 37,38 phase transitions, 39 spin-phonon coupling, 40 and magnetic excitations. 41,42 In addition, unlike other measurements that require bulk, large-area crystals, such as neutron diffraction, X-ray diffraction, or magnetic susceptibility, Raman spectroscopy can probe atomically thin flakes with diffraction-limited spatial resolution. The Raman spectra of bulk X PS3 materials has been studied since the 1980's, 28,43,44 and only recently extended to samples in the monolayer limit for NiPS3 45 and FePS3. 46,47 In particular, FePS3 is an interesting candidate to study using Raman spectroscopy because, due to the antiferromagnetic alignment of the spins and the resulting increase in the unit cell, Brillouin zone folding leads to new modes appearing below T N. 46,47

In this work, we show that one of the modes that appears below T N, which has been recently attributed to a phonon, 46 is actually a quantized spin wave, i.e . a magnon. This antiferromagnetic magnon, with a frequency of 3.7 THz (≈122 cm -1 ), softens with temperature more rapidly than typical phonon modes and splits upon application of a magnetic field, as expected for antiferromagnetic magnons. The frequency of the non-degenerate magnons depends linearly on magnetic field, with an effective g-factor g ≈ 2. However, its symmetry behavior with respect to the polarization vectors of the incoming and scattered light contradicts the long accepted work of Fleury and Loudon which states that Raman scattering of first-order magnon excitations involves antisymmetric Raman tensors and can thus only be observed in crosspolarized light configurations. 41 Instead, we observe the magnon in FePS3 in both parallel and

cross configurations, thus showing that the established magnon symmetry rule lacks the generality that has been suggested by previous literature. We explain the symmetry of the observed magnon using the magnetic point group with the inclusion of complex tensor elements. Because of the small interlayer exchange coupling, and thus quasi-2D magnetic nature of bulk FePS3, the magnon observed herein is also expected to be quasi-2D. To the best of our knowledge, this work represents the first verification of a magnon in a quasi-2D magnet using magneto-Raman spectroscopy. Furthermore, this work highlights temperature-dependent, magneto-Raman spectroscopy as an important technique to explore properties of magnetic excitations, such as frequencies and lifetimes, and to aid in investigating next-generation, magnonic devices. FePS3 may prove a better candidate for magnon transport than the previously studied MnPS3 as the magnon frequency is an order of magnitude higher, promising faster switching capabilities. 22,48

## II. RESULTS AND DISCUSSIONS

A. Initial Characterization: Sample Details and Magnetization Large crystals of FePS3 were mechanically exfoliated using adhesive tape 49 onto 300 nm SiO2 thermally grown on Si(100). Thicker flakes were selected for Raman spectroscopy measurements due to the more intense scattering signal. Bulk FePS3 crystallizes in the monoclinic structure with the space group C2/m. 50 In the a b -plane, as is shown in Figure 1a, the Fe atoms form a honeycomb lattice. In the center of the Fe honeycomb lattice is a [P2S6] 4unit, where the P2 dimers orient normal to the surface and coordinate with six S atoms. The magnetic structure has been a source of controversy, 32,50-53 where in both proposed structures, the moments of the large spin Fe 2+ ions (S=2) are aligned normal to the a b -plane and ordered in ferromagnetic chains that couple

antiferromagnetically. The proposed structures differed from each other, however, by a 120° rotation of the ferromagnetic chains, as well as whether the layers themselves are ferromagnetically or antiferromagnetically coupled. Recent neutron scattering by D. Lançon et al. 53 of large, single crystals of FePS3 confirmed the magnetic structure consists of ferromagnetic chains oriented along the a -axis, with antiferromagnetically coupled layers along the c -axis, as shown in Figures 1b and 1c.

Using a SQUID magnetometer, the magnetization of the FePS3 crystals ( i.e. before exfoliation) was measured both parallel and perpendicular to the a b -plane with an applied magnetic field of 0.1 T. The magnetization versus temperature from 5 K to 300 K is shown in Figure 1d. Similar to earlier reports, 27,46,47 the strong anisotropy in FePS3 results in higher magnetization perpendicular to the a b -plane, which is along the spin direction in the antiferromagnetic state. In a 2D antiferromagnet, T N is defined to be the temperature where there is a maximum in the slope of the magnetization vs. temperature curve. 27 Thus, by taking the first derivative of the magnetization with respect to temperature (Figure 1e), T N is found to be ≈118 K. The broad maximum above T N in the magnetization vs. temperature for both parallel and perpendicular orientations is typical for low-dimensional magnetic systems and has been attributed to short-range spin-spin correlation. 27

B. Raman Spectroscopy and Emergence of New Modes The majority of the peaks in the Raman spectrum (Figure 1f) can be ascribed to phonons. From the symmetry of bulk FePS3, there are 30 zone-center phonons with irreducible representations: Γ = 8Ag + 6Au + 7Bg + 9Bu, of which only the Ag and Bg modes are Raman active. Figure 1f shows the Raman spectra of bulk FePS3 (laser excitation wavelength λ = 514.5 nm) in the 180  back-scattering geometry at T = ° 135 K and T = 15 K, which corresponds to above and well below the Néel temperature,

respectively. Phonons from the Γ point (ω &gt; 150 cm -1 ) are labeled by whether they are an Ag mode or a combination of Ag/Bg modes based on previous first principle calculations as well as their measured orientation-dependent frequency since a peak that is a combination of Ag/Bg modes will appear to shift in frequency as the light polarization is rotated with respect to the crystal orientation. 47,54 The phonon modes Ag 2 /Bg 2 , Ag 3 , Ag 4 /Bg 4 , and Ag 5 are mostly ascribed to vibrations of the (P2S6) 4unit, and have similar frequencies for the different X PS3 materials due to the lack of contribution from the metal atom. 54 While the Ag and Ag/Bg modes in the spectra are similar in terms of frequency and intensity when comparing above and below the Néel temperature, significant changes are seen in the modes below 150 cm -1 . Above the Néel temperature, there is a broad, asymmetric peak ' ψ ' that exists even up to room temperature. At T = 15 K, however, there are four new modes that appear in this frequency range, labeled:  ψ ,  ψ , 1 2 ψ3, and  ψ . 4

The evolution of the Raman spectra of mode  ψ  into four modes below T N has been observed since the 1980's, 28,43,44 yet there is still uncertainty regarding the nature of the transition and the origin of each component. It is well known that antiferromagnetic ordering in FePS3 leads to a doubling of the unit cell and subsequent zone folding, such that the M point in the paramagnetic Brillouin zone folds into the Γ point of the antiferromagnetically ordered Brillouin zone. From the phonon dispersion calculated by X. Wang et al. , 47 ψ 1 and  ψ 2 appear to be M point phonons present below the Néel temperature due to magnetic-ordering induced zone folding. X. Wang et al. also state that  ψ 3 is a K point phonon, however they do not indicate how it becomes Raman-active since zone folding does not cause the K point to fold into the Γ point. On the contrary, J.-U. Lee et al. 46 attribute these modes to very different origins. From density functional theory calculations, they suggest  ψ  is a Γ point phonon without magnetic ordering, 2

while  ψ  and  ψ 3 4 have both Γ and M point components, and would thus appear above T N. Other older publications also observed  ψ - ψ 1 4, but only comment on possible origins of  ψ 1 and  ψ . 2 44,55 Thus, although the emergence of these modes can provide insight on the transition to a magnetically ordered state in FePS3, there is still substantial debate as to the source of each mode. In this work, we shed light on the origin of  ψ , where we unequivocally show that it is a 4 magnon.

C. Temperature Dependence of Raman Modes: Figure 2a shows the temperature dependence of the Raman spectra from 110 K down to 40 K in steps of 10 K.  ψ  is not observed until 4 temperatures below 100 K, well below T N of ≈118 K. As the temperature is further decreased,  ψ 4 displays a dramatic blueshift (higher in frequency) with temperature, especially compared with the other modes, including  ψ ,  ψ 1 2,  ψ 3 and the Γ-point phonons. This effect is illustrated in Figure 2b, which shows the frequencies of the different modes relative to their frequency at 15 K ( 𝜔𝜔15K ), including standard errors from fitting each peak with a Voigt function. While the frequency changes of the other modes only range from 0.01% to 0.57% (1 cm -1 to 2 cm -1 shifts), as is typical for anharmonic lattice effects, the frequency of  ψ 4 changes by as much as 6.2% (8 cm -1 shift) as the temperature is increased to 100 K. The strong shift in frequency with respect to temperature offers the first indication that this mode is attributed to a magnon, as suggested by T. Sekine et al. in 1990 for FePS3. 56 Similar magnon temperature dependence has been reported for other antiferromagnets, such as FeF2, 57,58 NiF2, 59 and MnF2. 60 The assignment of this mode to a magnon is further supported by recent inelastic neutron scattering measurements that observed a Γ point magnon at ≈15.1 meV (122 cm -1 ), 53,61 which coincides in frequency to our observed mode  ψ . 4

D. Magnons and Magneto-Raman Spectroscopy of FePS3 Due to the contrasting behavior of phonons and magnons in an applied magnetic field, we utilize our unique magneto-Raman capabilities to investigate the magnetic field dependence of the modes in FePS3 to identify if they have magnetic character. We detail our findings from the magnetic-field dependent Raman spectra of FePS3. Then, by discussing how magnons behave under applied magnetic fields, we assign  ψ 4 as a Raman-active magnon.

An exfoliated flake of bulk FePS3 was cooled to T = 5 K, after which a magnetic field was applied normal to the a b -plane, i.e . parallel/antiparallel to the directions of the spins in the antiferromagnetically ordered state. As seen in Figure 3a, when the magnetic field is increased from 0 T to 9 T, ψ 4 splits into two modes, labeled as ψ 4 (1) and ψ 4 (2) , and the splitting increases with the applied magnetic field. Figure 3b shows the frequencies of ψ 4 (1) and ψ 4 (2) as a function of magnetic field, where the splitting is symmetric for positive and negative magnetic fields. From a linear fit, the extracted slopes for the ψ 4 (1) and ψ 4 (2) branches are 0.93 ± 0.02 cm -1 /T and 0.94 ± 0.01 cm -1 /T, respectively. This splitting only occurs for  ψ 4, where the other  ψ  modes and phonon modes remain as single peaks. As discussed below, this behavior matches well with the expected magnetic field dependence of magnons.

Spin waves, which as quasiparticles are called magnons, are collective excitations of the spins in magnetic materials. The classical dynamics of the magnetization at the Γ-point in an antiferromagnet can be modeled following F. Keffer and C. Kittel 62 as two interacting macrospins (Figure 3c) representing the spin-up sublattice magnetization (M1, pink) and spindown sublattice magnetization (M2, green) that orient in the +  and -  directions in equilibrium, z z respectively. The excitation of a magnon causes the magnetization on both sublattices to precess about their equilibrium directions. Figures 3(d) and 3(e) show side and top view illustrations of

the two normal modes for M1 and M2 precession. The frequencies of the normal modes are given by:

$$\omega _ { k = 0 } = \gamma \{ ( 2 H _ { E } + H _ { A } ) H _ { A } \} ^ { 1 / 2 } \pm \gamma H _ { 0 } \text{ \quad \ \ } ( 1 )$$

where 𝐻𝐸 is the exchange field, 𝐻𝐴 is the anisotropy field, 𝐻0 is an applied field, and 𝛾 = 𝑔𝜇𝐵𝜇0 2𝜋ħ is the gyromagnetic ratio, which for free electrons ( g ≈ 2.0023) 63 is approximately 28.025 GHz/T, or 0.9348 cm -1 /T. Without an applied magnetic field ( 𝐻0 = 0 ), the frequencies of the magnon normal modes are degenerate ( 𝐸 = 𝐸0 ). However, an applied magnetic field 𝐻0 along z will cause a Zeeman splitting, leading to two separate modes with 𝐸 &lt; 𝐸0 (Figure 3d) and 𝐸 &gt; 𝐸0 (Figure 3e), both changing in magnetic field with a slope of |𝛾| . From Equation 1, a splitting of the mode with applied magnetic field offers a conclusive method to determine if a Raman mode is a magnon, as a magnetic field dependence is not expected for a typical phonon.

The behavior of the normal modes of magnons in a magnetic field described above matches with the observed magnetic-field dependence of ψ 4 , where the application of a magnetic field parallel to the spin direction leads to a lifting of the degeneracy of the magnon frequency into two components, ψ 4 (1) and ψ 4 (2) . The frequencies of the two non-degenerate magnon modes obey a linear response as a function of magnetic field, in agreement with Equation 1, with the two slopes of equal (within error) to the predicted slope 𝛾 = 0.93481 cm -1 /T for the free electron limit of g ≈ 2.0023. Using our experimental values of 𝛾 from the slopes of the linear fits in Figure 3b, we estimate that for bulk FePS3, the effective magnon g-factor is g ≈ 1.99 ± 0.05.

Thus, through magneto-Raman spectroscopy, we can unequivocally assign the mode  ψ  at ≈122 4 cm -1 as a magnon. From the full-width at half maximum (FWHM) value of approximately 3 cm -1 , which includes instrument broadening, we estimate the lower bound of the magnon lifetime to be on the order of 10 picoseconds.

E. Magnon Symmetry Behavior and Magnetic Point Group Raman scattering from magnons in both ferromagnets and antiferromagnets has been theoretically discussed by numerous authors, with one of the most cited works by P.A. Fleury and R. Loudon from 1968. 41 In their work, the resulting polarization selection rules of Raman scattering from first order magnon excitations is purely antisymmetric with respect to the polarization of the incoming ( 𝜀𝑖 ) and scattered ( 𝜀𝑠 ) light. Thus, from their theory, a magnon mode can be observed only in the perpendicular (or crossed) polarization configurations ( 𝜀 𝑖 = 𝜀𝑠 ±90° ), regardless of the orientation with respect to the crystallographic axes. This theory has worked well to describe a variety of ferromagnetic and antiferromagnetic materials with Raman-observed magnons, including FeF2, 41 MnF2, 41 YFeO3, 64 and Cd1-xMnxTe. 65

To probe the symmetry behavior of the magnon in FePS3, 𝜀𝑖 and 𝜀𝑠 were selected as shown in the inset of Figure 4a, where incoming light 𝜀𝑖 is at an arbitrary angle φ with respect to the crystallographic b -axis and the scattered light 𝜀𝑠 is at a controlled angle θ with respect to 𝜀𝑖 . Figure 4a shows that the magnon in FePS3, can be observed in parallel polarization ( θ =0 ) in ° addition to perpendicular polarization ( θ =90 ). When ° θ is varied between 0  and 360 , a two-fold ° ° symmetric pattern emerges where the magnon never disappears, but instead has a minimum at θ =60  and maximum at ° θ =150 , although these angle values are dependent on ° φ .

The theory presented in Fleury and Loudon applies rigorously to magnetic ions with orbital angular momentum ground states of L=0, or those with crystal-field quenching of the ground state to L=0, 41 as is the case for the well-studied antiferromagnetic magnon example of FeF2. 66,67 For materials that do not satisfy the L=0 ground state condition, such as FePS3, 68 the magnetic excitations are more difficult to treat; the general 'antisymmetric' nature of the magnon Raman tensor no longer applies. Instead, we use the symmetry properties of the

magnetic point group of bulk FePS3 to understand the selection rules and polarization dependence of the  ψ  magnon mode shown in Figure 4. 4

When a material orders antiferromagnetically, the symmetry generally lowers because the spin up and spin down sites are no longer equivalent. In this regard, for the magnetic point group, unitary operations are symmetry operations that preserve the spin direction. Other operations that do not preserve the spin direction may become allowed with the application of time reversal symmetry and are referred to as antiunitary operations. The magnetic point group for bulk FePS3, which considers both unitary and antiunitary operations is 2´/m, where the prime means time reversal. The Raman tensors for this magnetic group are given by the co-representations: 69

$$D A ^ { \prime } = \begin{pmatrix} A & B & 0 \\ D & E & 0 \\ 0 & 0 & I \end{pmatrix}, \quad D A ^ { \prime \prime } = \begin{pmatrix} 0 & 0 & C \\ 0 & 0 & F \\ G & H & 0 \end{pmatrix}$$

The experimental 180  backscattering configuration incident on the honeycomb plane can access ° only the upper left 2 × 2 block of these tensors, thus eliminating consideration of 𝐷𝐴′′ . The Raman intensity is proportional to |𝜀 𝑖 ∙ 𝐷𝐴′ ∙ 𝜀 𝑠 | 2 , where 𝜀𝑖 = (sin 𝜑 , cos 𝜑 , 0) and 𝜀𝑠 = (sin 𝜑 + 𝜃 , cos 𝜑 + 𝜃 , 0) are the incident and scatter light polarization vectors, respectively.

It should be noted that if 𝐴 𝐵 𝐶 , , , and 𝐸 are purely real numbers, the intensity profile as a function of θ would necessarily have two nodes, i.e . two angles where the intensity goes to zero, which is contrary to the polar plot we observed in Figure 4b. These nodes, however, disappear from the 𝐷𝐴′ symmetry response when the tensor elements are complex, which is the general case for absorbing materials: 70

$$D A ^ { \prime } = \begin{pmatrix} A & b e ^ { i \delta _ { B } } \\ d e ^ { i \delta _ { D } } & f e ^ { i \delta _ { F } } \end{pmatrix}$$

where 𝛿𝐵 , 𝛿𝐷 , and 𝛿𝐹 are the complex phase factors of the real tensor elements 𝐵 𝐷 , , and 𝐹 with respect to 𝐴 . Considering the complex tensor elements, there are multiple combinations of the

amplitude and phase factors that reproduce the polar plot in Figure 4b. However, there are two points that remain clear. First, at least one of the phase factors must be non-zero in FePS3 to agree with the non-nodal nature of the experimental polar plot. Second, the symmetry response of the magnon M in FePS3 deviates strongly from the antisymmetric nature predicted by Fleury and Loudon 41 for magnetic ions with L=0 ground state. Thus, there is no reason, a-priori , to expect that a magnon will appear only in cross-polarization configurations for these honeycomblike structures, such as FePS3, α -RuCl3, or CrI3, 71 without considering the magnetic point group symmetry and associated Raman tensors.

In principle, one should be able to use the generators of rotations ( i.e. angular momenta 𝐽 𝑥 , 𝐽 𝑦 , and 𝐽 𝑧 ) generally given in the character table as guides for the magnon symmetry behavior. For example, in FePS3 where the spins are oriented along the  -direction (normal to z a b -plane), the excitation of a magnon reduces the magnitude of the magnetization 𝑀𝑧 , which is then converted to the precession of the spins in the x y -( a b -) plane, 𝑀𝑥 and 𝑀𝑦 , as depicted in Figures 3d and 3e. Thus, a magnon in FePS3, which is a magnetization rotation in the x y -plane, can only have the same symmetry as 𝐽 𝑥 and 𝐽 𝑦 , which are generators of rotations. While the character tables that state these transformations are calculated for the non-magnetic space groups, 72-74 we were unable to find a reference that tabulated these for the magnetic space groups. Deriving such transformations is out beyond the scope of this work.

## III. CONCLUSIONS

In conclusion, we examined the temperature- and magnetic field-dependent Raman spectra of bulk FePS3, where new Raman-active modes appear in the antiferromagnetically ordered state. One of these modes ( ψ ), with a frequency of ≈122 cm 4 -1 (15.1 meV), was previously assigned as a phonon mode appearing due to zone folding. However, the strong temperature shift (compared

to other phonon modes) and splitting (and linear shifting) of the mode with applied magnetic field observed herein allows this mode to be unequivocally assigned as a magnon. We also investigated the magnon symmetry behavior, which is not purely antisymmetric, and explained its behavior using the magnetic point group of FePS3 (2´/m). This work will aid in future studies of magnetic excitations in similar magnetic layered materials, such as α -RuCl3 and CrI3, and demonstrates temperature- and magnetic field-dependent Raman spectroscopy as a technique to probe magnon phenomena in 2D materials, including for possible magnon transport applications. The quasi-2D magnetic nature of bulk FePS3, with weak interlayer exchange coupling, indicates that the magnon in bulk FePS3 is also quasi-2D. To the best of our knowledge, this is the first verification of a quasi-2D magnon in a layered material using magneto-Raman spectroscopy.

## IV. EXPERIMENTAL DETAILS

Bulk FePS3 crystals were purchased from 2D Semiconductors  and then exfoliated via the † adhesive tape onto 300 nm SiO2 thermally grown on Si(100). 49 Raman spectra were measured with the  514.5308  nm  excitation  wavelength  of  an  Argon  ion  laser  in  the  180   backscattering ° configuration using a triple grating Raman spectrometer (Horiba JY T64000 , 1800 mm † -1 grating) coupled to a liquid nitrogen cooled CCD detector. Polarization was selected and controlled using ultra  broadband  polarizers  and  achromatic  half  wave  plates.  To  perform  temperature-  and magnetic-field dependent Raman, the sample was placed into an attoDRY1000 cryostat (Attocube Inc. ), where the sample holder is pumped to ≈1x10 † -3 Pa (≈7x10 -6 Torr), backfilled with helium gas, and cooled. Micrometer-sized flakes were studied by focusing the laser with a white light camera onto the sample with a low-temperature, magnetic field compatible objective (50×, N.A. 0.82) and xyz nano-positioners. Integration times were approximately 12 minutes and the laser

power  was  kept  below  300  µW  to  reduce  local  heating  of  the  sample.  Spectra  with  applied magnetic field were corrected for Faraday rotation in the objective using half wave plates external to the cryostat. Magnetization vs. temperature measurements were performed on a Superconducting Quantum Interference Device-based Vibrating Sample Magnetometer (SQUIDVSM; Quantum Design, Inc. † ). The piece of the unexfoliated FePS3 crystal was mounted on to a quartz holder using GE varnish (LakeShore Cryotronics, Inc. ), and then measured under vacuum † (&lt; 1 kPa). The magnetization was measured under an applied magnetic field µ0H = 0.1 T from 5 K to 300 K and back down to 5 K. The temperature was varied in steps of 5 K between 5 K - 300 K, 1 K between 30 K - 150 K, and then 5 K between 150 K - 300 K.

† Certain commercial equipment, instruments, or materials are identified in this manuscript in order to  specify  the  experimental procedure  adequately. Such identification is  not intended to imply recommendation or endorsement by the National Institute of Standards and Technology, nor is it intended to imply that the materials or equipment are necessarily the best available for the purpose.

## V. ACKNOWLEDGEMENTS

A.M., T.T.M., and A.R.H.W. would like to acknowledge the National Institute of Standards and Technology (NIST)/National Research Council Postdoctoral Research Associateship Program and NIST-STRS (Scientific and Technical Research and Services) for funding. Work at The Ohio State University was supported by the Center for Emergent Materials, an NSF MRSEC under grant DMR-1420451. We would also like to acknowledge Prof. Natalia Drichko for insightful discussions and Hans Nembach for a careful readthrough of the paper.

Figure 1 . (a) Top view of the crystal structure of FePS3, where the Fe atoms form a honeycomb lattice. The a- and b-axes are labeled. (b) Magnetic structure of the Fe atoms in the a b -plane below T N. The spins are collinear and normal to the surface, where pink indicates spin 'up' (pointing out of the page) and green indicates spin 'down' (pointing into the page). The ferromagnetic chains are aligned along the a -axis. 50 (c) Below T N, the layers are stacked antiferromagnetically, resulting in an increase in the unit cell along the c-axis. (d) Magnetization ( M ) vs. temperature ( T ) for a bulk crystal of FePS3, where the applied magnetic field ( B = 0.1 T) is parallel (orange) or perpendicular (cyan) to the a b -plane. (e) dM/dT of (d) with an inflection point at T N ≈ 118 K. (f) Raman spectra above ( T = 135 K) and below ( T = 15 K) T N.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00608v1-figure-e560622d-4edf-4610-b3ff-85525f34dd17.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00608v1-figure-3d2597f2-4954-412e-8425-4b84a5f82a42.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00608v1-figure-649ffc5c-b4aa-4c85-9a7a-6f7ff4e95264.jpg)

Raman Shift (cm')

Figure 2. (a) Temperature-dependent Raman spectra of bulk FePS3 below T N. The  ψ  mode appears 4 below 100 K and its frequency dramatically increases with decreasing temperature. (b) Frequencies of the various modes (relative to frequency at 15 K). While the other  ψ  modes and phonon modes only shift by 1 cm -1 to 2 cm -1 , as is expected by lattice anharmonic effects, the frequency of  ψ  shifts by as much as 8 cm 4 -1 before it is no longer measurable above 100 K.

Figure 3. (a) Magnetic-field dependent Raman of FePS3 at T = 5K, showing the splitting of ψ 4 into two components, ψ 4 (1) and ψ 4 (2) , where the frequency of ψ 4 (1) ( ψ 4 (2) ) decreases (increases) with increasing magnetic field. (b) Frequency vs. magnetic field of ψ 4 (1) and ψ 4 (2) with the slopes of the linear fits for the two branches. (c) If we consider two magnetization sublattices M1 (pink) and M2 (green) in FePS3, in the ground state M 1 and M 2 point in the + z and -  directions. (d, e) The magnon spectrum includes two normal z modes with net moments Mnet in the x,y plane that precess in opposite directions. The degeneracy of these modes is broken in an applied field 𝐻0𝑧̂ with the energies of the modes in (d) and (e) respectively decreasing and increasing with 𝐻0 . In terms of this simple model, the modes in (d) and (e) would represent ψ 4 (1) and ψ 4 (2) in (b), respectively.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00608v1-figure-adf3e6bc-769f-4d8b-b08a-43bbb472d632.jpg)

Figure 4 . (a) Raman spectra ( T = 10K) showing intensity of magnon ψ 4 in parallel (θ=0 ) and cross ° (θ=90 ) polarization configurations at some arbitrary angle φ. Inset shows a schematic defining angles θ ° and φ with respect to the incoming (εi) and scattered (εs) light polarization and the a - and b - crystal axes. (b) Polar plot of intensity of the magnon as a function of θ at a constant angle φ. The intensity scale is represented by the radial lines that span from 0.1 to 1.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00608v1-figure-febce55c-faf5-4c9f-8c2f-84584591aa58.jpg)

## REFERENCES:

- [1] K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, Y. Zhang, S. V. Dubonos, I. V. Grigorieva, and A. A. Firsov . Electric Field Effect in Atomically Thin Carbon Films. Science 306 , 666-669 (2004).
- [2] Q. H. Wang, K. Kalantar-Zadeh, A. Kis, J. N. Coleman, and M. S. Strano. Electronics and Optoelectronics of Two-Dimensional Transition Metal Dichalcogenides . Nature Nanotechnol . 7 , 699-712 (2012).
- [3] K. S. Novoselov, A. Mishchenko, A. Carvalho, and A. H. C. Neto. 2D Materials and Van der Waals Heterostructures . Science 353 , aac9439 (2016).
- [4] D. S. Schulman, A. J. Arnold, A. Razavieh, J. Nasr, and S. Das. The Prospect of TwoDimensional Heterostructures: A Review of Recent Breakthroughs . IEEE Nanotechnol. Mag . 11 , 6-17 (2017).
- [5] Z. Lin, A. McCreary, N. Briggs, S. Subramanian, K. Zhang, Y. Sun, X. Li, N. J. Borys, H. Yuan, S. K. Fullteron-Shirey, et al. 2D Materials Advances: From Large Scale Synthesis and Controlled Heterostructures to Improved Characterization Techniques, Defects and Applications . 2D Mater. 3 , 042001 (2016).
- [6] K. F. Mak and J. Shan. Photonics and Optoelectronics of 2D Semiconductor Transition Metal Dichalcogenides . Nat. Photonics 10 , 216-226 (2016).
- [7] J. R. Schaibley, H. Yu, G. Clark, P. Rivera, J. S. Ross, K. L. Seyler, W. Yao, and X. Xu . Valleytronics in 2D Materials . Nat. Rev. Mater . 1 , 16055 (2016).
- [8] S. Manzeli, D. Ovchinnikov, D. Pasquier, O. V. Yazyev, and A. Kis. 2D Transition Metal Dichalcogenides . Nat. Rev. Mater. 2 , 17033 (2017).
- [9] D. Akinwande, C. J. Brennan, J. S. Bunch, P. Egberts, J. R. Felts, H. Gao, R. Huang, J. S. Kim, T. Li, Y. Li, et al. A Review on Mechanics and Mechanical Properties of 2D Materials-Graphene and Beyond . Extreme Mech. Lett . 13 , 42-77 (2017).
- [10] K. S. Burch, D. Mandrus, and J. G. Park. Magnetism in Two-Dimensional van der Waals Materials . Nature 563 , 47-52 (2018).
- [11] C. Gong, and X. Zhang. Two-Dimensional Magnetic Crystals and Emergent Heterostructure Devices . Science 363 , 706-717 (2019).
- [12] C. Lee, X. D. Wei, J. W. Kysar, and J. Hone. Measurement of the Elastic Properties and Antrinsic Strength of Monolayer Graphene . Science 321 , 385-388 (2008).
- [13] A. Castellanos-Gomez, M. Poot, G. A. Steele, H. S. J. van der Zant, N. Agrait, and G. Rubio-Bollinger . Elastic Properties of Freely Suspended MoS2 Nanosheets . Adv. Mater . 24 , 772-775 (2012).
- [14] B. Huang, G. Clark, D. R. Klein, D. MacNeill, E. Navarro-Moratalla, K. L. Seyler, N. Wilson, M. A. McGuire, D. H. Cobden, D. Xiao, et al. Electrical Control of 2D Magnetism in Bilayer CrI3 . Nat. Nanotechnol . 13 , 544-548 (2018).
- [15] S. W. Jiang, J. Shan, J. and K. F. Mak. Electric-Field Switching of Two-Dimensional van der Waals Magnets . Nat. Mater . 17 , 406-410 (2018).
- [16] J. S. Ross, P. Klement, A. M. Jones, N. J. Ghimire, J. Yan, D. G. Mandrus, T. Taniguchi, K. Watanabe, K. Kitamura, W. Yao, et al. Electrically Tunable Excitonic Light-Emitting Diodes Based on Monolayer WSe2 p-n Junctions . Nat. Nanotechnol . 9 , 268-272 (2014).
- [17] D. Akinwande, N. Petrone, and J. Hone. Two-Dimensional Flexible Nanoelectronics . Nat. Commun . 5 , 5678 (2014).

| [18]   | S. Zhang, H. M. Hill, K. Moudgil, C. A. Richter, A. R. Hight Walker, S. Barlow, S. R. Marder, C. A. Hacker, and S. J. Pookpanratana . Controllable, Wide-Ranging n-Doping and p-Doping of Monolayer Group 6 Transition-Metal Disulfides and Diselenides . Adv. Mater . 30 , 1802991 (2018).   |
|--------|-----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------|
| [19]   | Y. Q. Guo, K. Xu, C. Z. Wu, J. Y. Zhao, and Y. Xie. Surface Chemical-Modification for Engineering the Intrinsic Physical Properties of Inorganic Two-Dimensional Nanomaterials . Chem. Soc. Rev . 44 , 637-646 (2015).                                                                        |
| [20]   | M. Y. Li, C. H. Chen, Y. M. Shi, and L. J. Li. Heterostructures Based on Two- Dimensional Layered Materials and Their Potential Applications . Mater. Today 19 , 322- 335 (2016).                                                                                                             |
| [21]   | A. V. Chumak, V. I. Vasyuchka, A. A. Serga, and B. Hillebrands, B. Magnon Spintronics . Nat. Phys . 11 , 453-461 (2015).                                                                                                                                                                      |
| [22]   | W. Y. Xing, L. Qiu, X. Wang, Y. Yao, Y. Ma, R. Cai, S. Jia, X. C. Xie, and W. Han . Magnon Transport in Quasi-Two-Dimensional van der Waals Antiferromagnets . Phys. Rev. X 9 , 011026 (2019).                                                                                                |
| [23]   | N. D. Mermin and H. Wagner. Absence of Ferromagnetism or Antiferromagnetism in One- or Two-Dimensional Isotropic Heisenberg Models . Phys. Rev. Lett . 17 , 1133-1136 (1966).                                                                                                                 |
| [24]   | C. Gong, L. Li, Z. Li, H. Ji, A. Stern, Y. Xia, T. Cao, W. Bao, C. Wang, Y. Wang, et al. Discovery of Intrinsic Ferromagnetism in Two-Dimensional van der Waals Crystals . Nature 546 , 265-269 (2017).                                                                                       |
| [25]   | B. Huang, G. Clark, E. Navarro-Moratalla, D. Klein, R. Cheng, K. Seyler, D. Zhong, E. Schmidgall, M. A. McGuire, D. H. Cobden, et al. Layer-Dependent Ferromagnetism in a van der Waals Crystal Down to the Monolayer Limit . Nature 546 , 270-273 (2017).                                    |
| [26]   | S. W. Jiang, L. Z. Li, Z. F. Wang, K. F. Mak, and J. Shan. Controlling Magnetism in 2D CrI 3 by Electrostatic Doping . Nat. Nanotechnol . 13 , 549-553 (2018).                                                                                                                                |
| [27]   | P. A. Joy and S. Vasudevan. Magnetism in the Layered Transition-Metal Thiophosphates MPS 3 (M = Mn, Fe, and Ni) . Phys. Rev. B 46 , 5425-5433 (1992).                                                                                                                                         |
| [28]   | M. Balkanski, M. Jouanne, G. Ouvrard, and M. Scagliotti. Effects due to Spin Ordering in Layered MPX 3 Compounds Revealed by Inelastic Light-Scattering . J. Phys. C-Solid State 20 , 4397-4413 (1987).                                                                                       |
| [29]   | B. L. Chittari, Y. Park, D. Lee, M. Han, A. H. MacDonald, E. Hwang, and J. Jung . Electronic and Magnetic Properties of Single-layer MPX 3 Metal Phosphorous Trichalcogenides . Phys. Rev. B 94 , 184428 (2016).                                                                              |
| [30]   | C. R. S. Haines, M. J. Coak, A. R. Wildes, G. I. Lampronti, C. Liu, P. Nahai-Williamson, H. Hamidov, D. Daisenberger, and S. S. Saxena . Pressure-Induced Electronic and Structural Phase Evolution in the van der Waals Compound FePS 3 . Phys. Rev. Lett. 121 , 266801 (2018).              |
| [31]   | M. Tsurubayashi, K. Kodama, M. Kano, K. Ishigaki, Y. Uwatoko, T. Watanabe, K. Takase, and T. Takano . Metal-Insulator Transition in Mott-Insulator FePS 3 . AIP Adv. 8 , 101307 (2018).                                                                                                       |
| [32]   | K. C. Rule, G. J. McIntyre, S. J. Kennedy, and T. J. Hicks. Single-Crystal and Powder Neutron Diffraction Experiments on FePS 3 : Search for the Magnetic Structure . Phys. Rev. B 76 , 134402 (2007).                                                                                        |

| [33]   | X. -L. Li, X. -F. Qiao, W. -P. Han, Y. Lu, Q. -H. Tan, X. -L. Liu, and P. -H. Tan. Layer Number Identification of Intrinsic and Defective Multilayered Graphenes Up to 100 Layers by the Raman Mode Intensity from Substrates . Nanoscale 7 , 8135-8141 (2015).                            |
|--------|--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------|
| [34]   | L. B. Liang, J. Zhang, B. G. Sumpter, Q. -H. Tan, P. -H. Tan, and V. Meunier. Low- Frequency Shear and Layer-Breathing Modes in Raman Scattering of Two-Dimensional Materials . ACS Nano 11 , 11777-11802 (2017).                                                                          |
| [35]   | T. M. G. Mohiuddin, A. Lombardo, R. R. Nair, A. Bonetti, G. Savini, R. Jalil, N. Bonini, D. M. Basko, C. Galiotis, N. Marzari, et al. Uniaxial Strain in Graphene by Raman Spectroscopy: G Peak Splitting, Gruneisen Parameters, and Sample Orientation . Phys. Rev. B 79 , 205433 (2009). |
| [36]   | R. Beams, L. G. Cancado, and L. Novotny. Raman Characterization of Defects and Dopants in Graphene. J. Phys.-Condens. Mat . 27 , 083002 (2015).                                                                                                                                            |
| [37]   | E. del Corro, A. Botello-Mendez, Y. Gillet, A. L. lias, H. Terrones, S. Feng, C. Fantini, D. RHodes, N. Pradhan, L. Balicas et al. Atypical Exciton-Phonon Interactions in WS 2 and WSe 2 Monolayers Revealed by Resonance Raman Spectroscopy . Nano Lett . 16 , 2363-2368, (2016).        |
| [38]   | A. McCreary, J. R. Simpson, Y. Wang, D. Rhodes, K. Fujisawa, L. Balicas, M. Dubey, V. H. Crespi, M. Terrones, and A. R. Hight Walker . Intricate Resonant Raman Response in Anisotropic ReS 2 . Nano Lett. 17 , 5897-5907 (2017).                                                          |
| [39]   | K. N. Zhang, C. Bao, Q. Gu, X. Ren, H. Zhang, K. Deng, Y. Wu, Y. Li, J. Feng, and S. Zhou . Raman Signatures of Inversion Symmetry Breaking and Structural Phase Transition in Type-II Weyl Semimetal MoTe 2 . Nat. Commun. 7 , 13552 (2016).                                              |
| [40]   | Y. Tian, M. J. Gray, H. W. Ji, R. J. Cava, K. S. Burch. Magneto-Elastic Coupling in a Potential Ferromagnetic 2D Atomic Crystal . 2D Mater . 3 , 025035 (2016).                                                                                                                            |
| [41]   | P. A. Fleury and R. Loudon. Scattering of Light By One- and Two-Magnon Excitations . Phys. Rev. 166 , 514-530 (1968).                                                                                                                                                                      |
| [42]   | A. Sethi, T. Byrum, R. D. McAuliffe, S. L. Gleason, J. E. Slimak, D. P. Shoemaker, and S. L. Cooper. Magnons and Magnetodielectric Effects in CoCr 2 O 4 : Raman Scattering Studies . Phy. Rev. B 95 , 174413 (2017).                                                                      |
| [43]   | M. Scagliotti, M. Jouanne, M. Balkanski, G. Ouvrard, and G. Benedek. Raman- Scattering in Antiferromagnetic FePS 3 and FePSe 3 Crystals . Phys. Rev. B 35 , 7097-7104 (1987).                                                                                                              |
| [44]   | M. Bernasconi, G. L. Marra, G. Benedek, L. Miglio, M. Jouanne, C. Julien, M. Scagliotti, and M. Balkanski. Lattice-Dynamics of Layered MPX 3 (M=Mn, Fe, Ni, Zn, X=S, Se) Compounds . Phys. Rev. B 38 , 12089-12099 (1988).                                                                 |
| [45]   | C. -T. Kuo, M. Neumann, K. Balamurugan, H. J. Park, S. Kang, H. W. Shiu, J. H. Kang, B. H. Hong, M. Han, T. W. Noh, and J. -G. Park . Exfoliation and Raman Spectroscopic Fingerprint of Few-Layer NiPS 3 van der Waals Crystals . Sci. Rep . 6 , 20904 (2016).                            |
| [46]   | J. -U. Lee, S. Lee, S. H. Ryoo, S. Kang, T. Y. Kim, P. Kim, C. -H. Park, J. -G. Park, and H. Cheong. Ising-Type Magnetic Ordering in Atomically Thin FePS 3 . Nano Lett . 16 , 7433-7438 (2016).                                                                                           |
| [47]   | X. Z. Wang, K. Du, Y. Y. F. Liu, P. Hu, J. Zhang, Q. Zhang, M. H. S. Owen, X. Lu, C. K. Gan, and P. Sengupta. Raman Spectroscopy of Atomically Thin Two-Dimensional Magnetic Iron Phosphorus Trisulfide (FePS 3 ) Crystals . 2D Mater . 3 , 031009 (2016).                                 |

| [48]   | A. R. Wildes, B. Roessli, B. Lebech, and W. Godfrey. Spin Waves and the Critical Behaviour of the Magnetization in MnPS 3 . J. Phys.: Condens. Matter 10 , 6417-6428 (1998).                                                                                |
|--------|-------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------|
| [49]   | K. S. Novoselov, D. Jiang, F. Schedin, T. J. Booth, V. V. Khotkevich, S. V. Morozov, and A. K. Geim . Two-Dimensional Atomic Crystals . P. Natl. Acad. Sci. USA 102 , 10451- 10453 (2005).                                                                  |
| [50]   | G. Leflem, R. Brec, G. Ouvard, A. Louisy, and P. Segransan. Magnetic-Interactions in the Layer Compounds MPX 3 (M=Mn, Fe, Ni; X=S, Se) . J. Phys. Chem. Solids 43 , 455- 461 (1982).                                                                        |
| [51]   | K. Kurosawa, S. Saito, and Y. Yamaguchi. Neutron-Diffraction Study on MnPS 3 and FePS 3 . J. Phys. Soc. Jpn. 52 , 3919-3926 (1983).                                                                                                                         |
| [52]   | K. C. Rule, S. J. Kennedy, D. J. Goossens, A. M. Mulders, and T. J. Hicks. Contrasting Antiferromagnetic Order Between FePS 3 and MnPS 3 . Appl. Phys. A-Mater . 74 , S811- S813 (2002).                                                                    |
| [53]   | D. Lancon, H. C. Walker, E. Ressouche, B. Ouladdiaf, K. C. Rule, G. J. McIntyre, T. J. Hicks, H. M. Ronnow, and A. R. Wildes. Magnetic Structure and Magnon Dynamics of the Quasi-Two-Dimensional Antiferromagnet FePS 3 . Phys. Rev. B 94 , 214407 (2016). |
| [54]   | A. Hashemi, H. P. Komsa, M. Puska, and A. V. Krasheninnikov. Vibrational Properties of Metal Phosphorus Trichalcogenides from First-Principles Calculations . J. Phys. Chem. C 121 , 27207-27217 (2017).                                                    |
| [55]   | M. Scagliotti, M. Jouanne, M. Balkanski, and G. Ouvrard. Spin Dependent Phonon Raman-Scattering in Antiferromagnetic FePS 3 Layer-Type Compound . Solid State Commun . 54 , 291-294 (1985).                                                                 |
| [56]   | T. Sekine, M. Jouanne, C. Julien, and M. Balkanski. Light-Scattering Study of Dynamic Behavior of Antiferromagnetic Spins in the Layered Magnetic Semiconductor FePS 3 . Phys. Rev. B 42 , 8382-8393 (1990).                                                |
| [57]   | P. A. Fleury, S. P. S. Porto, L. E. Cheesman, and H. J. Guggenheim. Light Scattering by Spin Waves in FeF 2 . Phys. Rev. Lett. 17 , 84-87 (1966).                                                                                                           |
| [58]   | R. C. Ohlmann and M. Tinkham. Antiferromagnetic Resonance in FeF 2 at Far-Infrared Frequencies . Phys. Rev . 123 , 425-434 (1961).                                                                                                                          |
| [59]   | P. A. Fleury. Paramagnetic Spin Waves and Correlation Functions in NiF 2 . Phys. Rev . 180 , 591-593 (1969).                                                                                                                                                |
| [60]   | F. M. Johnson and A. H. Nethercot. Antiferromagnetic Resonance in MnF 2 . Phys. Rev . 114 , 705-716 (1959).                                                                                                                                                 |
| [61]   | A. R. Wildes, K. C. Rule, R. I. Bewley, M. Enderle and T. J. Hicks. The Magnon Dynamics and Spin Exchange Parameters of FePS 3 . J. Phys.-Condens. Mat . 24 , 416004 (2012).                                                                                |
| [62]   | F. Keffer and C. Kittel. Theory of Antiferromagnetic Resonance . Phys. Rev . 85 , 329-337 (1952).                                                                                                                                                           |
| [63]   | P. J. Mohr, D. B. Newell, and B. N. Taylor. CODATA Recommended Values of the Fundamental Physical Constants: 2014 . J. Phys. Chem. Ref. Data 45 , 043102 (2016).                                                                                            |
| [64]   | R. M. White, R. J. Nemanich, and C. Tsang. Raman-Scattering from Magnons in Rare- Earth Ortho-Ferrites . J. Magn. Magn. Mater. 15-18 , 773-774 (1980).                                                                                                      |
| [65]   | S. Venugopalan, A. Petrou, R. R. Galazka, A. K. Ramdas, and S. Rodriguez. Raman- Scattering by Phonons and Magnons in Semimagnetic Semiconductors - Cd 1-x Mn x Te . Phys. Rev. B 25 , 2681-2696 (1982).                                                    |

- [66] D. I. Khomskii. Transition Metal Compounds . (Cambridge University Press, 2014).
- [67] J. W. Stout and S. A. Reed. The Crystal Structure of MnF2, FeF2, CoF2, NiF2 and ZnF2 . J . Am. Chem. Soc. 76 , 5279-5281 (1954).
- [68] P. Jernberg, S. Bjarman, and R. Wappling. FePS3: A First-Order Phase Transition in a '2D' Ising Antiferromagnet . J. Magn. Magn. Mater. 46 , 178-190 (1984).
- [69] A. P. Cracknell. Scattering Matrices for Raman Effect in Magnetic Crystals . J. Physics Part C Solid 2 , 500-511 (1969).
- [70] H. B. Ribeiro, M. A. Pimenta, C. J. S. de Matos, R. L. Moreira, A. S. Rodin, J. D. Zapata, E. A. T. de Souza, and A. H. C. Neto . Unusual Angular Dependence of the Raman Response in Black Phosphorus . ACS Nano 9 , 4270-4276 (2015).
- [71] W. Jin, H. H. kim, Z. Ye, S. Li, P. Rezaie, F. Diaz, S. Siddiq, E. Wauer, B. Yang, C. Li, S. Tian, et al. Raman Fingerprint of Two Terahertz Spin Wave Branches in a TwoDimensional Honeycomb Ising Ferromagnet . Nat. Commun . 9 , 5122 (2018).
- [72] M. I. Aroyo, A. Kirov, C. Capillas, J. M. Perez-Mato, and H. Wondratschek. Bilbao Crystallographic Server: II. Representations of Crystallographic Point Groups and Space Groups. Acta Crystallogr . A 62 , 115-128 (2006).
- [73] M. I. Aroyo, J. M. Perez-Mato, C. Capillas, E. Kroumova, S. Ivantchev, G. Madariaga, A. Kirov, H. Wondratschek. Bilbao Crystallographic Server: I. Databases and Crystallographic Computing Programs . Z. Kristallogr . 221 , 15-27 (2006).
- [74] M. I. Aroyo, A. Kirov, C. Capillas, J. M. Perez-Mato, and H. Wondratschek . Crystallography Online: Bilbao Crystallographic Server . Bulg. Chem. Commun . 43 , 183-197 (2011).
Amber McCreary Jeffrey R. Simpson Thuc T. Mai Robert D. McMichael Jason E. Douglas Nicholas Butch Cindi Dennis Rolando Valdes Aguilar Angela R. Hight Walker
Quasi-Two-Dimensional Magnon Identification in Antiferromagnetic FePS3 via Magneto-Raman Spectroscopy
Recently it was discovered that van der Waals-bonded magnetic materials
retain long range magnetic ordering down to a single layer, opening many
avenues in fundamental physics and potential applications of these fascinating
materials. One such material is FePS3, a large spin (S=2) Mott insulator where
the Fe atoms form a honeycomb lattice. In the bulk, FePS3 has been shown to be
a quasi-two-dimensional-Ising antiferromagnet, with additional features in the
Raman spectra emerging below the Neel temperature of approximately 120 K. Using
magneto-Raman spectroscopy as an optical probe of magnetic structure, we show
that one of these Raman-active modes in the magnetically ordered state is
actually a magnon with a frequency of of approximately 3.7 THz (122 cm-1).
Contrary to previous work, which interpreted this feature as a phonon, our
Raman data shows the expected frequency shifting and splitting of the magnon as
a function of temperature and magnetic field, respectively, where we determine
the g-factor to be approximately 2. In addition, the symmetry behavior of the
magnon is studied by polarization-dependent Raman spectroscopy and explained
using the magnetic point group of FePS3.