Molecular intercalation in the van der Waals antiferromagnets FePS${}_{3}$ and NiPS${}_{3}$

Abstract

We have performed electrochemical treatment of the van der Waals antiferromagnetic materials FePS${}_{3}$ and NiPS${}_{3}$ with the ionic liquid EMIM-BF${}_{4}$, achieving significant molecular intercalation. Mass analysis of the intercalated compounds, EMIM${}_{x}$-FePS${}_{3}$ and EMIM${}_{x}$-NiPS${}_{3}$, indicated respective intercalation levels, $x$, of approximately 27% and 37%, and X-ray diffraction measurements demonstrated a massive (over 50%) enhancement of the $c$-axis lattice parameters. To investigate the consequences of these changes for the magnetic properties, we performed magnetic susceptibility and ${}^{31}$P nuclear magnetic resonance (NMR) studies of both systems. For EMIM${}_{x}$-FePS${}_{3}$, intercalation reduces the magnetic ordering temperature from $T_{N}=120$ K to 78 K, and we find a spin gap in the antiferromagnetic phase that drops from 45 K to 30 K. For EMIM${}_{x}$-NiPS${}_{3}$, the ordering temperature is almost unaffected (changing from 148 K to 145 K), but a change towards nearly isotropic spin fluctuations suggests an alteration of the magnetic Hamiltonian. Such relatively modest changes, given that the huge extension of the $c$ axes is expected to cause a very strong suppression any interlayer interactions, point unequivocally to the conclusion that the magnetic properties of both parent compounds are determined solely by two-dimensional (2D), intralayer physics. The changes in transition temperatures and low-temperature spin dynamics in both compounds therefore indicate that intercalation also results in a significant modulation of the intralayer magnetic interactions, which we propose is due to charge doping and localization on the P sites. Our study offers chemical intercalation with ionic liquids as an effective method to control not only the interlayer but also the intralayer interactions in quasi-2D magnetic materials.

I Introduction

Low-dimensional magnetic systems have played a pivotal role not only in advancing our understanding of the quantum properties of materials but also in the realization and exploration of new concepts in many-body physics. Still, with the exception of certain geometrically discretized or strong-frustration scenarios, which lead to quantum disordered magnetic states including quantum spin liquids, most three-dimensional (3D) magnetic systems exhibit long-ranged order. Conversely, as emphasized by Bethe in his seminal work, low dimensions amplify quantum fluctuations and destabilize the conventional order parameters Bethe (1931). In the case of layered compounds, increasing the interlayer spacing by the insertion of large molecules is a very literal means of modulating the dimensionality, and hence the properties, from potentially 3D towards the 2D limit Alias and Sukumaran (1992); Pattayil and Sukumaran (1993). The discovery of 2D magnetic materials suitable for this type of control would hold significant promise for applications in nanoelectronics and spintronics Soumyanarayanan et al. (2016).

For this reason, magnetic van der Waals materials have attracted extensive attention in recent years Burch et al. (2018), although efforts at dimensional manipulation have so far been limited largely to approaching the monolayer limit by exfoliation. In this context, the Ising-type van der Waals magnets Cr${}_{2}$Ge${}_{2}$Te${}_{6}$ and CrI${}_{3}$ exhibit the emergence of intrinsic ferromagnetism with a high transition temperature in few- or monolayer films Gong et al. (2017); Huang et al. (2017); Tian et al. (2016); in particular, CrI${}_{3}$ shows a systematic evolution of the ordering temperature and even the type of magnetic order with the number of layers Huang et al. (2017). MnSe${}_{2}$ remains ferromagnetic (FM) at room temperature in its monolayer form O’Hara et al. (2018), and most surprisingly room-temperature FM order emerges in monolayer VSe${}_{2}$ despite the bulk material being paramagnetic Bonilla et al. (2018).

The transition-metal trisulfide $M$PS${}_{3}$ ($M$ = Mn, Fe, Co, Ni) is a class of antiferromagnetic (AFM) van der Waals materials that has also been studied extensively Le Flem et al. (1982); Ouvrard et al. (1985); Joy and Vasudevan (1992); Chittari et al. (2016); Du et al. (2016); Wang et al. (2018). As Fig. 1(a) shows, the $M$PS${}_{3}$ compounds are isostructural, with the magnetic ions ($M^{2+}$) forming a honeycomb lattice. All exhibit the characteristics of a Mott insulator, displaying high resistivity at room temperature and a band gap well in excess of 1 eV (1.5 eV for FePS${}_{3}$ and 1.6 eV for NiPS${}_{3}$). This gap varies widely on exfoliation, offering an optoelectronic response over a broad frequency range for device applications Du et al. (2016). FePS${}_{3}$ exhibits a phase transition under pressure from insulating to metallic Haines et al. (2018) and NiPS${}_{3}$ in its AFM phase exhibits coherent excitonic states Kang et al. (2020).

Differences in the trigonal distortion of the $M$S${}_{6}$ octahedra and in spin-orbit coupling mean that the anisotropy of intralayer magnetic interactions differs in the $M$PS${}_{3}$ materials. MnPS${}_{3}$ appears to have Heisenberg spin interactions, NiPS${}_{3}$ shows a weak and CoPS${}_{3}$ a stronger easy-plane (XY) anisotropy, and FePS${}_{3}$ has a strong Ising anisotropy Kurosawa et al. (1983); Joy and Vasudevan (1992); Chandrasekharan and Vasudevan (1994); Chatterjee (1995); Rule et al. (2007); Lee et al. (2016); Wildes et al. (2017). Upon cooling, they all order in an AFM pattern, with respective transition temperatures, $T_{N}$, of 78 K, 118 K, 122 K, and 155 K for MnPS${}_{3}$, FePS${}_{3}$, CoPS${}_{3}$ and NiPS${}_{3}$ Le Flem et al. (1982); Brec (1986); Wang et al. (2018).

The weak van der Waals interlayer coupling leads to a cleavage energy close to that of graphite, and thus the $M$PS${}_{3}$ materials are easy to exfoliate. Multiple efforts to thin $M$PS${}_{3}$ samples have shown them to remain structurally stable down to a single layer Du et al. (2016); Kuo et al. (2016). Recent Raman scattering studies of monolayer FePS${}_{3}$ have reported that $T_{N}$ either drops from 117 K to 104 K Wang et al. (2016) or remains unchanged Lee et al. (2016), which raises the possibility of substrate effects on this 2D Ising magnet McCreary et al. (2020). Raman investigation of monolayer NiPS${}_{3}$ indicates a suppression of long-range order relative to the bulk material, which was interpreted in the framework of the Berezinskii-Kosterlitz-Thouless (BKT) phase transition Hu et al. (2023).

An alternative approach to dimensionality reduction is the incorporation of organic cations as spacers in the bulk materials, such that the altered layer separation enables control over the interlayer interactions Zhang et al. (2020, 2022). This is complementary to electrochemical treatments with ionic liquids whose primary effect is to induce protonation, or other electron doping effects, when the sample is placed on the cathode side Cui et al. (2018, 2019). It was reported in the early literature that pyridine (C${}_{5}$H${}_{5}$N) can be intercalated into MnPS${}_{3}$, leading to a transition from an AFM to a weakly FM ground state Alias and Sukumaran (1992), and that the intercalation of alkylamines (C${}_{n}$H${}_{2n+1}$NH${}_{2}$) into FePS${}_{3}$ causes a strong reduction of $T_{N}$ Pattayil and Sukumaran (1993). Much more recently, the intercalation of tetraheptylammonium-bromide (C${}_{28}$H${}_{60}$BrN) into NiPS${}_{3}$ was found to cause a transition from AFM to ferrimagnetic (FIM) order, followed by another transition to AFM order, in effects ascribed to the electron doping of the layers Mi et al. (2022).

In this work, we report the successful interlayer intercalation of EMIM${}^{+}$ into FePS${}_{3}$ and NiPS${}_{3}$, using the ionic liquid C${}_{6}$H${}_{11}$N${}_{2}$BF${}_{4}$ (EMIM-BF${}_{4}$). By measuring changes in the mass and interlayer spacing, we estimate the intercalation levels to be approximately 0.27 EMIM${}^{+}$/f.u. in FePS${}_{3}$ and 0.37 EMIM${}^{+}$/f.u. in NiPS${}_{3}$, and that both procedures dilate the $c$-axis lattice parameter by over 50%. We investigated the magnetic properties of EMIM${}_{x}$-FePS${}_{3}$ and EMIM${}_{x}$-NiPS${}_{3}$, and compared them with FePS${}_{3}$ and NiPS${}_{3}$, by magnetic susceptibility and ${}^{31}$P nuclear magnetic resonance (NMR) measurements.We found that intercalation in FePS${}_{3}$ causes changes in the Curie-Weiss temperatures determined above $T_{N}$, a suppression of $T_{N}$ itself, and a reduced spin gap at the lowest temperatures. For intercalated NiPS${}_{3}$, the magnetic transition temperature barely changes and the spin fluctuations become very isotropic, suggesting an evolution from weakly XY toward Heisenberg-type magnetism. Taken together with a dramatic suppression of any interlayer magnetic interactions, these data provide strong evidence for the systematic modulation of intralayer magnetic interactions by intercalation, for which we deduce that charge doping on the P site should be taken into account.

The structure of this article is as follows. In Sec. II we describe our $M$PS${}_{3}$ samples and intercalation procedure, and in Sec. III report our structural characterization. Section IV reports our measurements of the magnetic properties of pure and intercalated FePS${}_{3}$ and Sec. V our data for pure and intercalated NiPS${}_{3}$. In Sec. VI we discuss the consequences of our results for the understanding of magnetism in $M$PS${}_{3}$ materials and provide a brief conclusion.

II Materials and Methods

Figure 1: Lattice and magnetic structure of $M$PS${}_{3}$. (a) Crystal structure. The figure illustrates two unit cells containing eight $M$PS${}_{3}$ units. (b) Magnetic structure below $T_{N}$. $J_{1}$, $J_{2}$, and $J_{3}$ label magnetic interactions between nearest-, next-nearest-, and next-next-nearest-neighbor sites, respectively; $J^{\prime}$ labels the interlayer interaction, which we note corresponds to ABC stacking of the magnetic ions in the honeycomb layers. Blue and red arrows illustrate the relative orientations of local moments in the ordered phases.

Single crystals of FePS${}_{3}$ and NiPS${}_{3}$ were grown by the method of chemical vapor transport (CVT) Brec (1986). As Fig. 1(a) shows, the $M$PS${}_{3}$ layer is composed of covalently bonded $M$S${}_{6}$ octahedra and double-cone [P${}_{2}$S${}_{6}$]${}^{4-}$ units in a 2:1 ratio Brec et al. (1985). The $M$ ions are thought to be in a robustly divalent $M^{2+}$ state in all compounds, as shown for FePS${}_{3}$ by X-ray photoelectron spectroscopy (XPS) Yu et al. (2019). These ions form the honeycomb lattice, with the P-P pairs located at the centers of the $M$ hexagons Brec (1986); Joy and Vasudevan (1992); Le Flem et al. (1982). The S${}^{2-}$ layers stack in an ABC configuration along the $c$ direction to form a monoclinic structure with space group C2/m Ouvrard et al. (1985).

Both FePS${}_{3}$ and NiPS${}_{3}$ order magnetically with a “zig-zag” pattern, as shown in Fig. 1(b). The Fe${}^{2+}$ ions in FePS${}_{3}$ take their high-spin, $S=2$ state and order with the moments normal to the $ab$ plane Joy and Vasudevan (1992). These moments have FM alignment parallel to the $a$ axis and AFM alignment along the $b$ and $c$ axes, with propagation wave vector $\bm{k}=[0~{}1~{}\frac{1}{2}]$ Kurosawa et al. (1983); Rule et al. (2007); Wildes et al. (2012); Lançon et al. (2016); Coak et al. (2021); Paul et al. (2023). For NiPS${}_{3}$, the Ni${}^{2+}$ ($S=1$) moments are oriented primarily along the $a$ axis with a small component perpendicular to the $ab$ plane; the zig-zag chains also lie along the $a$ axis, with AFM alignment along $b$ but an FM repeat along the $c$ axis, resulting in $\bm{k}=[0~{}1~{}0]$ Kurosawa et al. (1983); Wildes et al. (2015); Kim and Park (2021); Wildes et al. (2022).

Figure 2(a) shows the configuration that we adopt for electrochemical treatment in order to achieve interlayer intercalation. The ionic liquid EMIM-BF${}_{4}$ was packed in a container with two platinum electrodes placed approximately 20 mm apart and subjected to a 3 V potential difference Cui et al. (2019, 2018). A single crystal of FePS${}_{3}$ or NiPS${}_{3}$ was attached to the cathode and covered with silver paint. The ionic liquid was heated to a temperature around 60${}^{\circ}$C. After 24 hours the mass of the samples had changed significantly, but longer treatment times led to no further change, and so one could declare the intercalation process to be complete.

X-ray diffraction (XRD) measurements were performed using Cu${}_{\alpha}$ and Cu${}_{\beta}$ radiation to determine the $c$-axis lattice parameters. The d.c. magnetic susceptibility was measured in a Magnetic Property Measurement System (MPMS) with a field of 100 Oe and at temperatures down to 1.8 K in field-cooled (FC) and zero-field-cooled (ZFC) conditions. ${}^{31}$P has nuclear spin $I=1/2$ and a Zeeman factor of ${}^{31}\gamma=17.235$ MHz/T. The ${}^{31}$P NMR measurements were performed on single crystals using a top tuning circuit and the spectra collected by spin-echo pulse sequences. For broad spectra, frequency sweeps were used to acquire the full spectrum. NMR Knight shifts were calculated from $K_{n}=(f/^{31}\gamma H-1){\times}100\%$, where $f$ is the resonance frequency of the spectral peaks in the paramagnetic phase and the average frequency in the ordered phase. Spin-lattice relaxation rates, $1/^{31}T_{1}$, were measured by the standard magnetization inversion-recovery method. These constitute a sensitive probe of low-energy spin fluctuations, because $1/T_{1}T={\rm lim}_{\omega\to 0}\Sigma_{q}A^{2}_{\rm hf}(q)\,{\rm Im}\,\chi(q,\omega)/\omega$, where $\chi(q,\omega)$ is the dynamical susceptibility, $A_{\rm hf}(q)$ is the hyperfine coupling, and $\omega$ is the NMR measurement frequency, which for electronic spins lies in the zero-energy limit.

III Lattice structure and doping of EMIM${}_{x}$-FePS${}_{3}$ and EMIM${}_{x}$-NiPS${}_{3}$

Figure 2: Sample intercalation and structural characterization of FePS${}_{3}$ and NiPS${}_{3}$. (a) Configuration for electrochemical treatment with EMIM-BF${}_{4}$. Two platinum electrodes are placed in the ionic liquid, set to a potential difference of 3 V, and the sample is attached to the negative electrode. (b) Single-crystal XRD data for FePS${}_{3}$ and EMIM${}_{x}$-FePS${}_{3}$. (c) Single-crystal XRD data for NiPS${}_{3}$ and EMIM${}_{x}$-NiPS${}_{3}$.

Figure 2(b) shows single-crystal XRD data for our FePS${}_{3}$ samples before and after 24-hour electrochemical treatment with the ionic liquid EMIM-BF${}_{4}$, and Fig. 2(c) the analogous data for NiPS${}_{3}$. We focus on the (0 0 $L$) Bragg peaks in order to extract the interlayer ($c$-axis) dimension of the unit cell. Studying the intralayer structure, meaning the $a$ and $b$ parameters, requires powder XRD measurements, but it has been found that the grinding process introduces disorder within the layers that precludes a meaningful analysis. The positions of the (0 0 $L$) XRD peaks change appreciably after the intercalation treatment, and the $c$-axis lattice parameters are shown and compared in Table I. The $c$-axis dimension increases from 6.72 Å to 10.40 Å for FePS${}_{3}$ and from 6.64 Å to 10.34 Å for NiPS${}_{3}$. Such an extremly large (over 50%) increase in interlayer spacing can be expected to have a very strong effect on reducing the dimensionality to the 2D limit.

Because the sample is attached to the cathode, and because of the large change in $c$-axis parameter, we deduce that EMIM${}^{+}$, rather than H${}^{+}$, is intercalated into the materials. A similar observation of intercalation by a large organic molecule has also been reported in NiPS${}_{3}$ using a different type of ionic liquid Mi et al. (2022). We note that our intercalation in EMIM${}_{x}$-NiPS${}_{3}$ was not complete by volume, in that a small portion of the XRD pattern of EMIM${}_{x}$-NiPS${}_{3}$ coincided with that of pristine NiPS${}_{3}$ [visible in Fig. 2(c)]. However, the volume ratio of residual NiPS${}_{3}$ is rather small, as we establish later from the absence of detectable signals in the magnetic susceptibility and NMR spectra.

Table 1: Lattice parameters of FePS${}_{3}$, EMIM${}_{x}$-FePS${}_{3}$, NiPS${}_{3}$, and EMIM${}_{x}$-NiPS${}_{3}$, and the estimated doping, $x$.

	FePS${}_{3}$	EMIM${}_{x}$-FePS${}_{3}$	NiPS${}_{3}$	EMIM${}_{x}$-NiPS${}_{3}$
$c$ (Å)	6.723(1)	10.401(1)	6.642(1)	10.343(1)
$x$	0	$27\pm 3$%	0	$37\pm 0.4$%

The mass of the FePS${}_{3}$ crystal increased from $0.18\pm 0.01$ mg before to $0.21\pm 0.01$ mg after intercalation. From this we estimate the intercalant concentration and the charge doping level to be $x=27\pm 3\%$, given the molecular mass of 183 g/mol for FePS${}_{3}$ and 111 g/mol for EMIM${}^{+}$. For the NiPS${}_{3}$ crystal, the mass changed from $3.86\pm 0.01$ mg to $4.71\pm 0.01$ mg, and thus we extracted the doping level as $x=37\pm 0.4\%$, using the molecular mass of 186 g/mol for NiPS${}_{3}$. Our NMR studies revealed in addition the presence of F${}^{-}$ ions in the intercalated samples, but at an unknown small concentration (data not shown), which suggests that the actual doping $x$ may vary with position across the samples.

IV Experimental data on FePS${}_{3}$ and EMIM${}_{x}$-FePS${}_{3}$

IV.1 Magnetic susceptibility

Our susceptibility measurements were conducted on single crystals of FePS${}_{3}$ and EMIM${}_{x}$-FePS${}_{3}$ with a small field applied parallel and perpendicular to the crystalline $ab$ plane, and data under FC and ZFC conditions are shown in Fig. 3. In both pure and intercalated FePS${}_{3}$, $\chi(T)$ is much bigger at temperatures above $T_{N}$ when the field is perpendicular to the $ab$ plane than when it lies in the $ab$ plane, but much smaller below $T_{N}$. Such an anisotropy should indicate the Ising nature of the system, which is attributed to combination of trigonal distortion and spin-orbit coupling Chatterjee (1995); Chandrasekharan and Vasudevan (1994); Joy and Vasudevan (1992); Amirabbasi and Kratzer (2023). With $H\perp ab$, $\chi$ does not reach zero in the zero-temperature limit, which may be the consequence of a small field misalignment.

Figure 3: Magnetic susceptibility of FePS${}_{3}$ and EMIM${}_{x}$-FePS${}_{3}$. $\chi(T)$ measured in a field of 100 Oe applied parallel and perpendicular to the $ab$ plane, under ZFC and FC conditions. (a) Susceptibility of FePS${}_{3}$ shown as a function of temperature. (b) $1/\chi(T)$ of FePS${}_{3}$ in the high-temperature PM phase. (c) Susceptibility of EMIM${}_{x}$-FePS${}_{3}$ shown as a function of temperature. (d) $1/\chi(T)$ of EMIM${}_{x}$-FePS${}_{3}$ in the high-temperature PM phase. $T_{N}$ in panels (a) and (c) labels the AFM transition temperature, which is located at the peak of $\chi(T)$. Solid lines in panels (b) and (d) are linear fits to the data made over the temperature range shown, with $\theta$ obtained from the intercept on the $x$ axis (i.e. at $y=0$).

On cooling in FePS${}_{3}$ [Fig. 3(a)], $\chi(T)$ exhibits a peak followed by a sharp drop at 120 K in both field orientations, indicating the Néel transition at $T_{N}\simeq 120$ K as found in previous reports Joy and Vasudevan (1992); Le Flem et al. (1982); Okuda et al. (1983). There is no appreciable hysteresis with FC and ZFC conditions, which supports a second-order AFM transition, in contrast to earlier reports that this transition may be of first order Jernberg et al. (1984); Ferloni and Scagliotti (1989). Turning to EMIM${}_{x}$-FePS${}_{3}$, $T_{N}$ taken from the peak position in $\chi(T)$ is reduced to 80 K by intercalation [Fig. 3(c)], which is a value still 2/3 of that in the pristine sample. These values of $T_{N}$ will be verified by the NMR measurements to follow, and already contain one of our primary results, that there is little to no connection between $T_{N}$ and a hypothetical interlayer coupling. The FC and ZFC data for EMIM${}_{x}$-FePS${}_{3}$ already differ at temperatures above $T_{N}$, but this difference is affected only minimally by the onset of magnetic order at $T_{N}$; thus we believe that the difference between the FC and ZFC datasets is not an intrinsic effect, but an extrinsic one caused by impurities that enter the sample during chemical intercalation.

In all cases shown in Fig. 3, $\chi(T)$ above $T_{N}$ can be well fitted by the Curie-Weiss (CW) form, $\chi=\dfrac{N\mu_{0}\mu_{\rm eff}^{2}}{3k_{B}(T-\theta)}$, where $\mu_{\rm eff}$ is the effective paramagnetic (PM) moment and $\theta$ is the CW temperature. In FePS${}_{3}$ [Fig. 3(b)], the CW temperatures are found to be $\theta_{\perp}=22\pm 1$ K and $\theta_{\parallel}=-104\pm 2$ K for the two field orientations. In our convention, the positive $\theta_{\perp}$ indicates the presence of FM interactions, as one may expect from the zig-zag ground-state order; nevertheless, this ordered state is dominated by AFM interactions and by the Ising anisotropy, so our result reflects the magnetic frustration in the system. The CW fit to the perpendicular-field data returns a PM moment $\mu_{\rm eff}=6.4\pm 0.2\mu_{B}$/Fe, which is close to the literature value Joy and Vasudevan (1992); as noted in early studies, such a large $\mu_{\rm eff}$ for divalent Fe${}^{2+}$ ions can be explained only by taking the spin-orbit coupling into account Joy and Vasudevan (1992).

For EMIM${}_{x}$-FePS${}_{3}$, the susceptibility shows the same type of anisotropy as the pristine sample [Fig. 3(c)]. The PM moments we extract, $6.0\pm 0.2\mu_{B}$/Fe, are largely unaltered, indicating that intercalation causes no significant changes to the ionic properties. The CW temperatures we obtain, $\theta_{\perp}=26\pm 1$ K and $\theta_{\parallel}=-18\pm 3$ K [Fig. 3(d)], suggest that no dramatic changes take place in the signs or energy scales of the couplings, but offer no straightforward interpretation as alterations to the dominant AFM or frustrating FM interactions. Because the susceptibility measurements contain significant impurity contributions whose subtraction is difficult to benchmark, we turn next from a bulk to a microscopic probe in order to gain deeper insight into the magnetic properties after intercalation.

IV.2 NMR spectra

Figure 4: ${}^{31}$P NMR spectra of FePS${}_{3}$ and EMIM${}_{x}$-FePS${}_{3}$. (a) Spectra of FePS${}_{3}$ over a range of temperatures above and below $T_{N}$, taken at a fixed in-plane field of 4 T. Successive datasets are offset vertically for clarity. (b) Spectra of EMIM${}_{x}$-FePS${}_{3}$ taken under the same conditions. The spectral intensities at 30 and 60 K are multiplied by a factor of two; dashed lines represent Gaussian fits to the data.

We measured ${}^{31}$P NMR spectra for FePS${}_{3}$ and EMIM${}_{x}$-FePS${}_{3}$ over a wide range of temperatures on both sides of $T_{N}$, and show our results in Fig. 4. In FePS${}_{3}$, the single NMR peak observed at temperatures above 120 K splits into two nearly symmetrical peaks, presenting a clear signature of AFM ordering at $T_{N}$ [Fig. 4(a)]. We note here that the field was applied in the $ab$ plane, where from Fig. 1(b) one expects the ${}^{31}$P nucleus to experience hyperfine fields from the two different sites that are in opposite directions relative to the external field. We then fitted the peak or peaks in each NMR spectrum with Gaussian profiles to extract the full width at half-maximum height (FWHM) and the peak splitting below $T_{N}$. As Fig. 5(a) shows, the single NMR peak above 120 K broadens rapidly at $T_{N}$, from approximately 18 to over 50 kHz where it splits. The line splitting, ${\Delta}f$, approaches 2.8 MHz at low temperatures and, as shown in Fig. 5(b), from its shape can serve as an order parameter for the AFM phase.

Figure 5: NMR line widths and splittings in FePS${}_{3}$ and EMIM${}_{x}$-FePS${}_{3}$. (a) FWHM of the primary spectral peak, or peaks, in FePS${}_{3}$ [Fig. 4(a)] shown as a function of temperature. (b) Splitting, $\Delta f$, of NMR peaks in FePS${}_{3}$ as a function of temperature. (c) FWHM of the broad NMR spectrum of EMIM${}_{x}$-FePS${}_{3}$ [Fig. 4(b)] shown as a function of temperature.

For EMIM${}_{x}$-FePS${}_{3}$, the single NMR peak remains rather narrow when cooled down to 80 K [Fig. 4(b)]. Here it shows no evidence of a line splitting, but instead undergoes a large increase in line width [Fig. 5(c)]. This broadening can be taken to indicate the onset of magnetic order at $T_{N}\simeq 80$ K, consistent with the susceptibility data of Fig. 3, and Gaussian peak fits show the FWHM broadening from 300 kHz above $T_{N}$ to 3 MHz below $T_{N}$. We draw attention to the fact that the low-temperature FWHM in EMIM${}_{x}$-FePS${}_{3}$ has the same order as ${\Delta}f$ in FePS${}_{3}$, which suggests a distribution of ordered moments in the system. Because we observe only one peak, we propose that these results reflect an incommensurate magnetic order in the 2D layers, with a systematic distribution of local-moment magnitudes and/or directions, as opposed to a multi-domain structure. Such an incommensurate order could arise from rather small alterations of the intralayer interactions in FePS${}_{3}$, whose pre-intercalation compromise configuration is the Ising zig-zag order shown in Fig. 1(b).

IV.3 Spin-lattice relaxation rates and spin gaps

Figure 6: Spin-lattice relaxation rates of FePS${}_{3}$ and EMIM${}_{x}$-FePS${}_{3}$. (a) $1/^{31}T_{1}$ shown as a function of temperature, measured with an in-plane field of 4 T. (b) $1/^{31}T_{1}$ shown in semilog form as a function of $1/T$. The straight solid lines represent linear fits to the data in the low-temperature regime. $\Delta$ is a gap value extracted from the gradient of each line.

The spin-lattice relaxation rates, $1/^{31}T_{1}$, of pure and intercalated FePS${}_{3}$ are shown in Fig. 6. FePS${}_{3}$ exhibits the characteristic behavior on decreasing temperature that $1/^{31}T_{1}$ first increases to a maximum at $T=120$ K, which reflects the AFM transition at $T_{N}$, followed by a rapid drop [Fig. 6(a)]. Far below $T_{N}$, the exponential fall suggests the form $1/^{31}T_{1}\propto e^{-{\Delta}/k_{B}T}$, where $\Delta$ is the spin gap expected in an ordered Ising magnet. Plotting $\ln(1/^{31}T_{1})$ as a function of $1/T$ reveals that the data below 50 K follow a straight line [Fig. 6(b)], whose gradient gives a gap $\Delta=45\pm 3$ K. Here we comment that inelastic neutron scattering (INS) measurements of the anisotropic spin-wave spectra have reported a gap of 15 meV Lançon et al. (2016), which is significantly larger than our NMR value. Because the NMR $1/T_{1}$ is a very sensitive, low-energy local probe that sums contributions from all of momentum space, it is possible that our results reveal a process missed by both triple-axis and time-of-flight INS for reasons of energy resolution or reciprocal-space coverage; however, the higher temperatures at which we obtained an NMR signal also leave open the possibility of probing a finite-energy transition, and we hope that future studies may resolve this apparent contradiction.

Similar magnetic ordering and spin-gap behavior are observed in EMIM${}_{x}$-FePS${}_{3}$, as Fig. 6(a) shows. Here $T_{N}$ is determined quite precisely at 78 K from the sharp peak in $1/^{31}T_{1}$. Analysis of the low-temperature response, shown in Fig. 6(b), gives the equivalent spin gap as $\Delta=30\pm 3$ K. This suppression of the gap upon intercalation scales linearly with the suppression of $T_{N}$, which would be consistent with a straightforward reduction of the relevant intralayer energies.

V Experimental data on NiPS${}_{3}$ and EMIM${}_{x}$-NiPS${}_{3}$

V.1 Magnetic susceptibility

Figure 7: Magnetic susceptibility of NiPS${}_{3}$ and EMIM${}_{x}$-NiPS${}_{3}$. $\chi(T)$ measured in an out-of-plane field of 100 Oe under FC and ZFC conditions.

Figure 8: ${}^{31}$P NMR spectra of NiPS${}_{3}$ and EMIM${}_{x}$-NiPS${}_{3}$. (a) Spectra of NiPS${}_{3}$ measured in an out-of-plane field of 3 T over a wide range of temperatures above and below $T_{N}$. Successive datasets are offset vertically for clarity. (b) As in panel (a) with an in-plane field. The low-temperature spectra are fitted to double-Lorentzian functions. (c) and (d) Spectra of EMIM${}_{x}$-NiPS${}_{3}$ measured under the same conditions of field and temperature.

The magnetic susceptibilities of NiPS${}_{3}$ and EMIM${}_{x}$-NiPS${}_{3}$ are shown in Fig. 7. For both compounds, the peak in $\chi(T)$ occurs around 250 K, which makes it impossible to reach a high-temperature regime in which to extract the CW behavior. Such high intralayer energy scales are found by high-precision INS measurements of the spin dynamics of NiPS${}_{3}$ Wildes et al. (2022), and the broad peak in $\chi$ indicates the development of 2D AFM correlations. We note that the AFM ordering transition is essential invisible in $\chi(T)$, and we will show by NMR that $T_{N}=148$ K (below). Detailed measurements on NiPS${}_{3}$ showed that the susceptibility is isotropic above $T_{N}$ and anisotropic below it Wildes et al. (2015). The most dramatic feature in our measurements (Fig. 7) is the onset of a large difference between the FC and ZFC data below 25 K, which may be caused by quenched disorder or a change of magnetic structure.

For EMIM${}_{x}$-NiPS${}_{3}$, $\chi(T)$ is remarkably similar to the pristine compound, with only a small deviation setting in at temperatures below 200 K. This minimal alteration despite the massive increase in interlayer spacing is a strong indication that all the magnetic properties of NiPS${}_{3}$ are intrinsically 2D. For the intercalated system, a difference between the ZFC and FC data develops gradually with decreasing temperature, unaffected by the onset of magnetic order (shown below by NMR), and we believe this to be a consequence of impurities introduced by the intercalation, similar to the situation in EMIM${}_{x}$-FePS${}_{3}$ (Fig. 3).

V.2 NMR spectra

Figure 9: NMR line widths and splittings in NiPS${}_{3}$ and EMIM${}_{x}$-NiPS${}_{3}$. (a) FWHM of the primary spectral peak, or peaks, in NiPS${}_{3}$ [Figs. 8(b)] shown as a function of temperature. (b) Splitting, $\Delta f$, of NMR peaks in NiPS${}_{3}$ as a function of temperature. (c) FWHM of the broad NMR spectrum of EMIM${}_{x}$-NiPS${}_{3}$ [Figs. 8(d)] shown as a function of temperature. The arrows mark the transition temperatures, $T_{N}$ and $T^{*}$, extracted for the two compounds.

Figure 8 shows the ${}^{31}$P NMR spectra of NiPS${}_{3}$ and EMIM${}_{x}$-NiPS${}_{3}$ over a wide range of temperature for both field orientations. For NiPS${}_{3}$ in an out-of-plane field, the spectra above 148 K have a single peak characteristic of the PM phase, which broadens significantly on further cooling from 148 K to 140 K [Fig. 8(a)]. By contrast, the double-peak feature developing at all temperatures below 148 K in an in-plane field demonstrates clearly the onset of AFM order and fixes $T_{N}=148\pm 2$ K. In this case, Lorentzian functions gave a slightly better fit to the spectra, although this choice made little difference for the purpose of extracting the FWHM of each peak and the line splitting, $\Delta f$, below $T_{N}$ [shown in Figs. 9(a) and 9(b)]. The FWHM is about 40 kHz at temperatures above $T_{N}$ and saturates around 0.2 MHz far below $T_{N}$, while ${\Delta}f$ also tends to level off around 0.2 MHz at low temperatures in the ordered state.

For EMIM${}_{x}$-NiPS${}_{3}$ in both field orientations, the NMR peaks also broaden significantly towards low temperatures [Fig. 8(c,d)]. However, we were not able to resolve any line splitting, even in an in-plane field. Again we used Lorentzian fits, to the peaks obtained in in-plane fields, in order to extract the FWHM shown in Fig. 9(c). The sudden increase in FWHM occurring at 145 K, which we label $T^{*}$, is characteristic of a static magnetic order. However, the absence of any line splitting at any lower temperatures suggests that this order could be of a type different from that in pristine NiPS${}_{3}$.

V.3 Spin-lattice relaxation rates

Figure 10: Spin-lattice relaxation rates of NiPS${}_{3}$ and EMIM${}_{x}$-NiPS${}_{3}$. $1/^{31}T_{1}$ shown as a function of temperature, measured with in- and out-of-plane fields of 3 T.

The spin-lattice relaxation rates of NiPS${}_{3}$ and EMIM${}_{x}$-NiPS${}_{3}$, measured in both field orientations, are presented in Fig. 10. We comment that it was not possible to obtain a reliable $1/^{31}T_{1}$ signal at temperatures significantly below 100 K, in contrast to the situation in FePS${}_{3}$ (Fig. 6), but our data characterize well the magnetic transitions in both materials, which are clearly evident from the kinks in $1/^{31}T_{1}$. For NiPS${}_{3}$ in both orientations, the rapid drop in $1/^{31}T_{1}$ determines precisely the AFM ordering temperature $T_{N}=148$ K, a value slightly lower than that obtained from early susceptibility data Le Flem et al. (1982); Joy and Vasudevan (1992); Ziolo et al. (1988). Below $T_{N}$, $1/^{31}T_{1}$ for an in-plane field drops significantly faster than with an out-of-plane field, which is consistent with XY magnetic order Rosenblum et al. (1994); Kim et al. (2019); Jana et al. (2023) in that the gapless (Goldstone) modes should be preserved only when the field is applied out of the plane.

For EMIM${}_{x}$-NiPS${}_{3}$, our $1/^{31}T_{1}$ data for the two field orientations are barely distinguishable above $T^{*}$, where they decrease with decreasing temperature from 200 K to 160 K and then increase over the range from 160 K to 145 K. Below this they drop more rapidly, allowing us to fix the putative magnetic ordering temperature as $T^{*}=145$ K. Below $T^{*}$, the relaxation rates under in- and out-of-plane fields are truly indistinguishable, presenting the first of two distinctive contrasts between NiPS${}_{3}$ and EMIM${}_{x}$-NiPS${}_{3}$. Such isotropic behavior suggests that intercalation in this system genuinely alters the intralayer magnetic interactions to a Heisenberg form. Second, although $T^{*}$ is reduced by only 3 K on intercalation, $1/^{31}T_{1}$ at $T=100$ K (far below $T^{*}$) is enhanced by a factor of five. This apparent strong enhancement of low-energy spin fluctuations points once again to the possibility of a magnetically ordered state different from that of pure NiPS${}_{3}$ [Fig. 1(b)], candidates for which would include an incommensurate order. Further measurements are therefore required to investigate the magnetic structure of the intercalated system.

Table 2: Magnetic interactions of FePS${}_{3}$ and NiPS${}_{3}$ Lançon et al. (2016); Wildes et al. (2022); $J_{i}$ indicates a Heisenberg interaction between neighboring ions of differing separations [Fig. 1(b)] and $D_{z}$ the leading magnetic anisotropy fitted by INS.

	$J_{1}$ (meV)	$J_{2}$ (meV)	$J_{3}$ (meV)	$J^{\prime}$ (meV)	$D_{z}$ (meV)
FePS${}_{3}$	$-1.46$	0.04	0.96	0.0073	$-2.66$
NiPS${}_{3}$	$-2.6$	0.2	13.5	$-0.3$	0.21

VI Discussion

Our susceptibility data show that the intercalation of FePS${}_{3}$ and NiPS${}_{3}$ does not cause any depletion of the available magnetic moments. Thus the changes we observe, specifically the reduction of $T_{N}$ in EMIM${}_{x}$-FePS${}_{3}$ and the isotropic spin response of EMIM${}_{x}$-NiPS${}_{3}$, should be attributed to changes of the magnetic interactions. To discuss these, in Table II we list the values of the magnetic interactions deduced by INS experiments for FePS${}_{3}$ Lançon et al. (2016) and NiPS${}_{3}$ Wildes et al. (2022). The authors of both studies adopted a minimal model consisting of the three intralayer Heisenberg interactions illustrated in Fig. 1(b), an interlayer Heisenberg interaction ($J^{\prime}$), and single-ion anisotropy terms, of which only one was found to be significant. The $J^{\prime}$ term was found to be very small in FePS${}_{3}$ but potentially significant in NiPS${}_{3}$. The $J_{3}$ term, for third-neighbor $M^{2+}$ ions, was found to be completely dominant in NiPS${}_{3}$ and to play an essential role in establishing the zig-zag order of FePS${}_{3}$.

Considering the increase of the $c$-axis lattice parameter that we achieve, in excess of 50%, it is natural to expect that the interlayer couplings should be reduced dramatically. Thus our primary conclusion is that the $J^{\prime}$ terms have no role in the magnetism of either material, where 2D ordering remains well established in both in EMIM${}_{x}$-FePS${}_{3}$ and in EMIM${}_{x}$-NiPS${}_{3}$. In the light of theoretical studies demonstrating that quasi-2D magnets with very small $J^{\prime}$ may nevertheless have significant $T_{N}$ values Yasuda et al. (2005), this result was not fully clear for FePS${}_{3}$ (Table II). Our study demonstrates that it is also true for NiPS${}_{3}$, despite the 3 K interlayer coupling scale extracted in Ref. Wildes et al. (2022).

It is then important to understand how the intralayer magnetic interactions are modified by intercalation. In principle, the addition of EMIM${}^{+}$ should lead to electron doping in both materials, with the electrons presumably located on the Fe or P ions, and reaching a doping around 1/3 e${}^{-}$/f.u. (0.27 in EMIM${}_{x}$-FePS${}_{3}$, 0.37 in EMIM${}_{x}$-NiPS${}_{3}$). Considering the unchanged paramagnetic moments we extracted for EMIM${}^{+}$-FePS${}_{3}$ and the unchanged magnetic properties of EMIM${}^{+}$-NiPS${}_{3}$, a change of the $M^{2+}$ valence seems to be excluded. Thus we propose that the electronic doping affects the P sites, by transforming the less favorable P${}^{4+}$ ions Yu et al. (2019) to P${}^{3+}$. With the P ions located at the centers of the $M$ hexagons, we expect that this doping should be very effective in modifying the intralayer superexchange interactions, particularly the strong $J_{3}$ terms. We remark in this context that we found the intercalated samples to remain highly insulating, which is a prerequisite to regard the doped charges as localized, such that they can affect the magnetic interactions rather than drive a metallic transition.

Considering FePS${}_{3}$ in more detail, the $T_{N}$ we measure for bulk EMIM${}_{x}$-FePS${}_{3}$ is lower than the value $T_{N}=104$ K obtained for monolayer FePS${}_{3}$ Wang et al. (2016). Thus our intercalation is more than a simple isolation of the 2D units, and its effects on the competing intralayer interactions require more detailed microscopic modelling. In fitting their INS results, the authors of Ref. Lançon et al. (2016) assumed an extremely minimal model where all of the anisotropy was ascribed to one single-ion term, and further anisotropic interactions such as Dzyaloshinskii-Moriya or Ising-type spin-spin interactions were neglected. The modelling of these terms is a complex task that was approached only recently by correlated density functional theory (DFT) methods Amirabbasi and Kratzer (2023), while further experimental information on the magnetic anisotropy has been obtained by X-ray photoelectron microscopy Lee et al. (2023). With the ongoing interest in FePS${}_{3}$ as a candidate spintronic material Soumyanarayanan et al. (2016), and for band engineering by ultrafast coherent light Mertens et al. (2023), a deeper understanding of the key interaction terms remains a priority.

In EMIM${}_{x}$-NiPS${}_{3}$, the fact that $T^{*}$ is only 3 K lower than $T_{N}$ in NiPS${}_{3}$ demonstrates immediately that $J^{\prime}$ terms are irrelevant. Nevertheless, the absence of line splitting at low temperatures in the ${}^{31}$P spectra of EMIM${}_{x}$-NiPS${}_{3}$ [Fig. 8(d)], in contrast to the behavior in NiPS${}_{3}$ [Fig. 8(b)], and the enhancement of $1/^{31}T_{1}$ (Fig. 10), lead us to suggest that the ordering at $T^{*}$ may be of a different type, such as incomplete (short-ranged) or incommensurate. Recent theoretical work has indeed proposed that NiPS${}_{3}$ is close to incommensurate magnetic order Mellado (2023), and the large overlap of the two NMR peaks we observe in Fig. 8(b) may support this scenario. Another possibility could be weak anisotropy caused by intralayer strain arising from a lattice mismatch with the size of the intercalated molecules. Finally, a BKT transition is also possible if the system retains an easy-plane anisotropy, however weak, and in the intercalated system it is likely that quenched disorder would help to pin true long-range order from the quasi-long-ranged BKT order. We note here that a BKT phase has already been reported in monolayer NiPS${}_{3}$ Hu et al. (2023).

In summary, we performed electrochemical treatment of FePS${}_{3}$ and NiPS${}_{3}$ with the ionic liquid EMIM-BF${}_{4}$, obtaining a significant intercalation of the large EMIM${}^{+}$ ions into both compounds. The most important consequence of intercalation is a large and uniform increase of the $c$-axis lattice parameters, which unavoidably causes a strong reduction of interlayer couplings. The fact that we do not observe order-of-magnitude changes in the magnetic properties of either system demonstrates conclusively that the magnetic properties of the pristine materials are intrinsically 2D, being determined almost exclusively by intralayer physics. The changes we do observe are a consequence of modulated intralayer magnetic interactions. In FePS${}_{3}$, the drop in magnetic ordering temperature suggests a change in the competing FM and AFM interactions while the system remains in the Ising limit; the effects in NiPS${}_{3}$ are far more subtle, with the interactions changing from weakly XY-type to almost isotropic (Heisenberg-type). We propose the charge doping that accompanies intercalation, which appears to be localized around the P ions and reaches a level of 0.3-0.4 e${}^{-}$/f.u., as the mechanism most likely to drive these changes. Our study of chemical intercalation with ionic liquids paves the way not only for approaching the 2D limit in van der Waals magnetic materials but also for tuning their intralayer magnetic interactions.