Understanding the Ising zigzag antiferromagnetism of FePS3 and FePSe3 monolayers

Abstract

This study investigates the spin-orbital states of FePS3 and FePSe3 monolayers and the origin of their Ising zigzag AFM, using DFT, crystal field level diagrams, superexchange analyses, and parallel tempering MC simulations. Our calculations show that under the trigonal elongation of the FeS6 (FeSe6) octahedra, the $e_g^{\pi }$ doublet of the Fe 3d crystal field levels lies lower than the $a_{1g}$ singlet by about 108 meV (123 meV), which is much larger than the strength of Fe 3d SOC. Then, the half-filled minority-spin $e_g^{\pi }$ doublet of the high-spin Fe${}^{2+}$ ions ($d^{5\uparrow ,1\downarrow }$) splits by the SOC into the lower $L_{z+}$ and higher $L_{z-}$ states. The spin-orbital ground state $d^{5\uparrow }$$L_{z+}^{1\downarrow }$ formally with $S_z$ = 2 and $L_z$ = 1 gives the large z-axis spin/orbital moments of 3.51/0.76 ${\mu }_B$ (3.41/0.67 ${\mu }_B$) for FePS${}_3$ (FePSe${}_3$) monolayer, and both the moments are reduced by the strong (stronger) Fe 3d hybridizations with S 3p (Se 4p) states. As a result, FePS3 (FePSe3) monolayer has a huge perpendicular single-ion anisotropy energy of 19.4 meV (14.9 meV), giving an Ising-type magnetism. Moreover, via the maximally localized Wannier functions, we find that the first nearest neighboring (1NN) Fe-Fe pair has large hopping parameters in between some specific orbitals, and so does the 3NN Fe-Fe pair. In contrast, the 2NN Fe-Fe pair has much smaller hopping parameters and the 4NN Fe-Fe pair has negligibly small ones. Then, a combination of those hopping parameters and the superexchange picture can readily explain the computed strong 1NN ferromagnetic coupling and the strong 3NN antiferromagnetic one but the relatively much smaller 2NN antiferromagnetic coupling. Furthermore, our PTMC simulations give TN of 119 K for FePS3 monolayer and also predict for FePSe3 monolayer the same magnetic structure with a close or even higher TN.

arXiv comment14 pages, 9 figures, 3 tables

References in corpus (19)

• Electric Field Effect in Atomically Thin Carbon Films
• Layer-dependent Ferromagnetism in a van der Waals Crystal down to the Monolayer Limit
• Discovery of intrinsic ferromagnetism in 2D van der Waals crystals
• Wannier90: A Tool for Obtaining Maximally-Localised Wannier Functions
• Maximally localized Wannier functions: Theory and applications
• Exchange Monte Carlo Method and Application to Spin Glass Simulations
• Magnetic 2D materials and heterostructures
• Giant Tunneling Magnetoresistance in Spin-Filter van der Waals Heterostructures
• Ising-Type Magnetic Ordering in Atomically Thin FePS3
• On the origin of magnetic anisotropy in two dimensional CrI${}_3$
• Switching 2D Magnetic States via Pressure Tuning of Layer Stacking
• Pressure-controlled interlayer magnetism in atomically thin CrI3
• Electronic and magnetic properties of single-layer MPX${}_3$ metal phosphorous trichalcogenides
• Interplay between Kitaev interaction and single ion anisotropy in ferromagnetic CrI${}_3$ and CrGeTe${}_3$ monolayers
• Hubbard U and Hund's Exchange J in Transition Metal Oxides: Screening vs. Localization Trends from Constrained Random Phase Approximation
• VI3: a 2D Ising ferromagnet
• Magnetic anisotropy and exchange interactions of two-dimensional FePS${}_3$, NiPS${}_3$ and MnPS${}_3$ from first principles calculations
• Orbital and magnetic ordering in single-layer FePS3: A DFT+U study
• Large orbital magnetic moment in VI3

Cited by in corpus (4)

• Magnetic Anisotropy in Two-dimensional van der Waals Magnetic Materials and Their Heterostructures: Importance, Mechanisms, and Opportunities
• Identifying band structure changes of FePS3 across the antiferromagnetic phase transition
• Giant Anisotropic Magnetoresistance in Magnetic Monolayers CrPX3 (X = S, Se, Te) due to symmetry breaking between the in-plane and out-of-plane crystallographic axes
• Contrasting magnetism in VPS3 and CrI3 monolayers with the common honeycomb S = 3/2 spin lattice