
==== Front
iScience
iScience
iScience
2589-0042
Elsevier

S2589-0042(24)02006-6
10.1016/j.isci.2024.110781
110781
Article
Magnetic anisotropy and phononic properties of two-dimensional ferromagnetic Fe3GeS2 monolayer
Wei Yu 1
Liu Hui 1
Wang Ke kewang@xupt.edu.cn
12∗
1 Xi’an University of Posts & Telecommunications, Shaanxi 710121, China
∗ Corresponding author kewang@xupt.edu.cn
2 Lead contact

22 8 2024
20 9 2024
22 8 2024
27 9 11078112 6 2024
14 8 2024
16 8 2024
© 2024 The Author(s)
2024
https://creativecommons.org/licenses/by/4.0/ This is an open access article under the CC BY license (http://creativecommons.org/licenses/by/4.0/).
Summary

In 2023, Fe3GeS2 monolayer with Curie temperature of 630 K is predicted, which is promising to be used in next-generation spintronic devices. However, its magnetic anisotropy and phononic properties are still unclear. In this paper, we implemented the first-principles calculations on Fe3GeS2 monolayer, and found its ferromagnetic ground state with robustness to the −1.5%–1.3% biaxial strain. Meanwhile, the out-of-plane magnetic anisotropy dominated by dipolar interaction is found in Fe3GeS2 monolayer. Finally, we studied the phononic properties to identify the dynamical stability of Fe3GeS2 monolayer and highlight the contribution from the anharmonic interaction of optical phonons to the thermal expansion coefficient. We also find two single-phonon modes can be used to design quantum mechanical resonators with a wide cool-temperature range. These results can provide a comprehensive understanding of the magnetism and phonon properties of two-dimensional (2D) Fe3GeS2, beneficial for the application of 2D Fe3GeS2 in spintronics.

Graphical abstract

Highlights

• The ground state of Fe3GeS2 monolayer is robust to the −1.5% to 1.3% biaxial strain

• There are two single-phonon modes (9.08 THz and 11.12 THz) in the phonon spectrum

• The out-of-plane magnetic anisotropy is dominated by dipolar interaction

• The contribution from optical phonons to thermal expansion coefficient is significant

Physics; Magnetism

Subject areas

Physics
Magnetism
Published: August 22, 2024
==== Body
pmcIntroduction

The successful exfoliation of two-dimensional (2D) graphene has opened up substantial research topics on the fascinating properties and potential applications of 2D materials. In the initial stage of studying 2D materials, researchers focused on black phosphorene, transition metal dichalcogenides (TMDs), hexagonal boron nitride, etc., all of which are nonmagnetic (NM). In 2017, the 2D ferromagnetic (FM) CrGeTe3 and CrI3 were mechanically exfoliated from the bulk crystals by Gong et al. and Huang et al.,1,2 which shed light on the applications of 2D materials in spintronics. Subsequently, a large of 2D magnets are predicted and prepared, including CrCl3,3 Fe3GeTe2,4 CrTe2,5 MnBi2Te46 with nontrivial topology, Weyl half-semimetal PtCl3,7 and antiferromagnetic (AFM) MnPS38 and MXenes.9,10 In recent years, these 2D magnets attracted extensive research attention in the field of materials science and condensed matter physics, due to their interesting properties, such as large spin Seebeck coefficient,11 controllable magnetoresistance,12 and quantized anomalous Hall effect.13 These intriguing properties render 2D magnets promising candidates in fabricating spintronic devices. Currently, spintronics devices have been widely applied in the fields of signal transfer,14 data storage,15 biomedicine,16 and energy conversion.17

However, a fatal flaw of 2D magnets is their weak magnetic stability and low magnetic phase-transition temperature (Curie temperature for FM; Néel temperature for AFM), which hinders strongly their application in spintronics. For instance, the Curie/Néel temperatures are ∼60 K for CrGeTe31, ∼45 K for CrI32, and 89 K for FePS3,18 which is much lower than the room temperature. Therefore, in recent years, researchers have devoted themselves to exploring methods to improve the magnetic phase-transition temperature of 2D magnetic materials or design novel 2D magnetic materials with high Curie/Néel temperature. At present, it is recognized that the out-of-plane magnetic anisotropy is key for 2D magnets to break the Mermin-Wagner theorem and to resist thermal fluctuation. The relationship between Curie/Néel temperature and the out-of-plane magnetic anisotropy can be written as:(Equation 1) Tc=−(ZJlf+ZKlf+Azz)S0·13KB,

where Jlf, Klf, and Azz are the isotropic magnetic exchange coupling coefficient, out-of-plane magnetic anisotropic parameter, and out-of-plane single-point anisotropic parameter, respectively. KB is the Boltzmann constant, while S0 is the altitude of spin vector on each magnetic lattice. The details of Equation 1 can be referred to a study by Wang et al.19 Hence, many strategies have been employed to enhance the out-of-plane magnetic anisotropy and elevate the Curie/Néel temperature of 2D magnets, such as strain engineering,20,21 charge doping,22,23 surface functionalization,24,25 atomic doping,26,27 intercalation,28,29,30 and external manipulation.31,32 Meantime, researchers also predicted many new 2D magnets with high Curie/Néel temperatures.33,34,35 For example, Wen et al.36 prepared 2D FM CuCr2Te4 flakes with thickness-dependent Curie temperature (260 K–320 K) by the heteroepitaxial growth. Based on Fe3GeTe2 monolayer with the Curie temperature of 150 K,37 Yang et al.38 predicted 60 types of easy exfoliable and highly stable magnetic A3BX2 monolayers by machine learning and high-throughput computation. Among these 60 types of A3BX2 monolayer, Fe3GeS2 monolayer with excellent dynamical and thermal stability owns the highest Curie temperature of 630 K, much higher than room temperature 300 K. Therefore, if it can be successfully prepared in time, it will bring revolutionary progress to spintronic devices based on 2D magnets, making 2D spintronic devices have great application potential in quantum communication, data storage, biomedicine, and other fields. However, in their study, the Curie temperature was estimated using Monte Carlo simulation based on the magnetic exchange coupling, without exploring the magnetic anisotropy in Fe3GeS2 monolayer.

In this paper, we would explore the magnetic anisotropy in Fe3GeS2 monolayer by the first-principles calculations. Besides, in the spintronic devices fabricated by the 2D magnets, the performance and reliability of the devices heavily depend on the thermal conduct and thermal expansion of 2D magnets. Phononic properties of 2D magnets are the basis for analyzing and exploring the thermal conduct and thermal expansion. Therefore, we would also delve into the phononic properties and thermal expansion of Fe3GeS2 monolayer, laying a theoretical foundation for its application in 2D spintronic devices. We found the out-of-plane magnetic anisotropy dominated by dipolar interaction in Fe3GeS2 monolayer, robust to the −1.5%–1.3% biaxial strain. Meantime, we also find two single-phonon modes with frequencies of 9.08 THz and 11.12 THz, which can be used to design quantum mechanical resonators with wide cool-temperature range. Besides, we also calculate the Grüneisen constants of all phonons and thermal expansion coefficient (TEC) of FM Fe3GeS2 monolayer, highlighting the contribution from the anharmonic interaction of optical phonons to the TEC.

Computational details

In our simulations, the calculations of structural relaxation, electronic density of states, and phonon dispersion were obtained by Device Studio program,39 which provides several functions for crystal visualization, modeling and simulation. In this paper, all these calculations were implemented by DS-PAW software integrated in Device Studio program.40 In DS-PAW software, the projected augmented wave (PAW) method was employed to describe the coupling between atomic nuclei and extra-nuclear electrons, while we selected the Perdew-Burke-Ernzerhof (PBE) method of general gradient approximation (GGA) as exchange-correlation functional.41,42 To suppress the non-physical interaction between adjacent layers, a 15 Å vacuum space was imposed along the out-of-plane direction. The convergence limits of energy and force were set as 10−7 eV and 0.001 eV/Å in the relaxation of geometrical structure with cutoff energy of 500 eV and 5×5×1 Monkhorst-Pack (MP) grid. The phonon dispersions were calculated using an MP grid of 9×9×1 for a 2×2×1 supercell, based on the density function perturbation theory (DFPT).43 Besides, it has been reported in previous studies44,45,46 that including the Hubbard “U” into the PBE functional will overestimate the lattice parameters and magnetic moment of based-Fe3GeTe2 materials seriously. Meantime, we calculated the band structures of Fe3GeS2 monolayer with different “U” values, as shown in Figure 1. It can be observed that the metallic state of Fe3GeS2 monolayer is robust to the value of Hubbard “U”. Therefore, the “U” parameter was not included in our calculations.Figure 1 DS-PAW calculated band structures of Fe3GeS2 monolayer with different “U” values

The used "Ueff" values are 0 eV (A), 1 eV (B), 2 eV (C), and 3 eV (D), respectively.

Results and discussion

In this paper, the model of crystal Fe3GeS2 is established by using the bulk Fe3GeTe2 as a template. The bulk Fe3GeS2 is hexagonal, with the space group of P63/mmc (No. 194). The Fe3GeS2 monolayer is exfoliated from the bulk Fe3GeS2, and its space group is P-6m2 (No. 187). The top and side views of Fe3GeS2 monolayer are shown in Figures 2A and 2B, where the Fe, Ge, and S atoms are represented by the brown, grayish purple and gold balls, respectively. In Fe3GeS2 monolayer, there are five sub-layers, where Fe3Ge substructure is sandwiched by two S layers. The red dashed line marks the unit cell of Fe3GeS2 monolayer in Figure 2A. In this unit cell, there are six atoms, including three Fe atoms, one Ge atom, and two S atoms. In the unit cell of Fe3GeS2 monolayer, these three Fe atoms can be divided in two inequivalent types Fe1 (including the top and bottom Fe atoms) and Fe2, based on their positions and the crystal symmetry. Hence, there are four possible magnetic phases in the unit cell of Fe3GeS2 monolayer, including FM, AFM-1, AFM-2, and NM magnetic configurations, as shown in Figure 3. When both the spins of two Fe1 atoms are parallel to that of Fe2 atom, it results in an FM phase, as shown in Figure 3A. In AFM-1 phase, the spin of top Fe1 is antiparallel to the bottom Fe1, but is parallel to that of Fe2. If the coupling between Fe2 and two Fe1 atoms is AFM, while the coupling between two Fe1 atoms is FM, it leads to AFM-2 phase. To determine the magnetic ground state, we calculated the energies of Fe3GeS2 monolayer at these four possible magnetic phases, as shown in Figure 4. We can find the Fe3GeS2 monolayer with FM phase owns the lowest energy at lattice constant of 3.95 Å. In Figures 3C and 3D, it can be found that the energies of AFM-1, AFM-2 and NM phases are ∼1.322 eV, ∼1.323 eV, and ∼1.346 eV higher than the FM phase as lattice constant is of 3.95 Å. Meantime, with the lattice constant increasing, these energy differences also elevate significantly. When the lattice constant is smaller than 3.95 Å, all of the energy differences between AFM-1, AFM-2, NM, and FM phases are still positive, indicating the energy of FM phase is lowest among these four possible magnetic phases. According to the least-energy principle, it can be drawn that the magnetic ground state of Fe3GeS2 monolayer is FM, and the optimized lattice constant is a = b = 3.95 Å. Furthermore, this FM magnetic ground state of Fe3GeS2 monolayer is robust to the biaxial strain within the range of −1.5%–1.3%.Figure 2 Geometrical structure and high-symmetry path of Fe3GeS2 monolayer

(A and B) are the top and side views of geometrical structure, and (C) is the high-symmetry path in the irreducible Brillouin zone. In (A) and (B), the Fe, Ge, and S atoms are represented by the brown, grayish purple and gold balls, respectively.

Figure 3 Four possible magnetic phases in the Fe3GeS2 unit cell

Four possible magnetic phases include ferromagnetic (FM) (A), antiferromagnetic-1 (AFM-1) (B), antiferromagnetic-2 (AFM-2) (C), and nonmagnetic (NM) (D) phases.

Figure 4 The curve of energy at possible magnetic phases versus strain

(A) is the energy of Fe3GeS2 monolayer at FM phase, while (B–D) are the energy differences between FM and, AFM-1, AMF-2 and NM phases obtained by DS-PAW, respectively.

For 2D magnets, magnetic anisotropy is essential to resist thermal fluctuation and keep the long-range magnetic order. Magnetic anisotropy, introduced by the spin-orbit coupling, is composed of magnetocrystalline anisotropic energy (C-MAE) and dipolar magnetic anisotropic energy (D-MAE). In this paper, C-MAE is calculated by the XXZ model47:(Equation 2) EC−MAE=EZZ−EXX,

where the EXX (EZZ) is the energy of 2D magnet with the spin along x- (z-) direction. In the XXZ model, the energy with spin along the x-direction is considered equal to the spin along the y-direction, resulting in-plane magnetic isotropy. Thus, the C-MAE only considers the energy difference of 2D magnets as the spin is along x- and z-directions. The dipolar interaction energy in magnet can be calculated by48:(Equation 3) Ed=μ0(gμB)22∑i≠j|Rij|2Si·Sj−3(Rij·Si)(Rij·Sj)|Rij|5,

where g of 2 is the gyromagnetic ratio, μB and μ0 are the Bohr magnon and the vacuum permeability, respectively. Rij is the position vector between the i- and j-magnetic lattices, while the spin angular momentum of the i- (j-) magnetic lattice is Si (Sj). According to Equation 3, it can be found that the dipolar interaction energy relies on the relative position of spin pair and their spin orientations. Early in 2002, Politi et al.49 reported the dipolar coupling between single-domain FM particles could induce long-range FM order. Based on the strong and tunable dipole-dipole interaction, Young et al.50 realized fast two-qubit entangling gates, providing significant speedups for quantum algorithms. Meantime, Utesov51 proposed that dipolar forces originating from the dipole-dipole interaction can lead to biaxial anisotropy in the reciprocal space of antiferromagnets with skyrmion, and the dipolar forces is considered as the critical ingredient to stabilize the nanometer-sized skyrmions in antiferromagnets.52 Hence, it is recognized that dipolar forces can result in the complicated sequences of magnetic phase transition in magnets under an external magnetic field.53,54 To obtain the D-MAE, we make a difference in the dipole interaction energy of spin angular momentum along the z- and x-directions. The obtained magnetic anisotropic energies are plotted in Figure 5. Interestingly, the C-MAE keeps positive revealing an in-plane magnetocrystalline anisotropy. At the optimized lattice constant (3.95 Å), the C-MAE is 0.28 meV. When the lattice constant is stretched to 3.97 Å, the C-MAE of Fe3GeS2 monolayer decreases to 0.01 meV indicating a nearly isotropic magnetism. As the lattice constant is compressed to 3.90 Å, the C-MAE increases to 0.56 meV, suggesting an enhancement of in-plane magnetocrystalline anisotropy. For D-MAE, it keeps negative, revealing an out-of-plane magnetic dipolar anisotropy. At the optimized lattice constant (3.95 Å), the D-MAE is −236.19 meV, while the sum of C-MAE and D-MAE is −235.91 meV, revealing an out-of-plane magnetic anisotropy in the optimized Fe3GeS2 monolayer. Furthermore, the D-MAE is within the range of −293.50∼-97.68 meV, with the lattice constant changing from 3.89 to 4.00 Å. These results suggest the Fe3GeS2 monolayer owns out-of-plane magnetic anisotropy dominated by the dipolar interaction, and this out-of-plane magnetic anisotropy is robust to the biaxial strain within the range of −1.5%–1.3%. It is worth emphasizing that this out-of-plane magnetic anisotropy dominated by the dipolar interaction is rare in previous synthesized 2D ferromagnets, such CrGeTe3, CrI3, Fe3GeTe2, which is promising to realize dipolar-induced magnon chirality.Figure 5 Magnetic anisotropy energy of Fe3GeS2 monolayer with strain

The magnetocrystalline anisotropic energy (C-MAE) (A) and magnetic dipole anisotropic energy (D-MAE) (B) of Fe3GeS2 monolayer obtained by DS-PAW.

To identify the dynamic stability of optimized Fe3GeS2 monolayer, we calculated its phonon dispersion based on the DFPT, shown in Figure 6. There is little imaginary frequency in the phonon dispersion and phonon density of states (PDOS), suggesting the dynamic stability of our optimized Fe3GeS2 monolayer. In phonon dispersion, there are 18 phonon branches including 3 acoustic branches and 15 optical branches, because there are six atoms in the unit cell of Fe3GeS2 monolayer. These three acoustic branches are named as longitudinal acoustic (LA), transverse acoustic (TA), out-of-plane acoustic (ZA) phonons, and their eigenvectors describe the translation of the Fe3GeS2 monolayer along x-, y-, and z-directions, respectively. It has been reported that the thermal conductivity usually is dominated by the anharmonic interaction of acoustic phonons,55,56 which is fundament for the performance and reliability of the 2D devices. We mark the LA, TA, and ZA branches by the red, green, and purple solid lines in phonon dispersion, as shown in Figure 6A. To estimate the anharmonicity of acoustic phonons in Fe3GeS2 monolayer, we calculated the Grüneisen constant (γq,s) by57:(Equation 4) γq,s=−a2ωq,s·dωq,sda,

where ωq,s is the phonon frequency of the s-branch at the wave vector of q, and a is the lattice constant. In this paper, the Grüneisen constant was calculated by applying a 2% biaxial strain to lattice at 0 K, and the obtained Grüneisen constants of acoustic phonons are presented in Figure 7A. As known, a larger absolute value of Grüneisen constant suggests a stronger anharmonicity. In Figure 6A, the absolute Grüneisen constant of ZA phonon is larger than LA and TA phonons, indicating the strong anharmonic interaction with other phonon modes. Meantime, the Grüneisen constant of ZA remains positive with the wave vector changing, suggesting a softening phonon frequency of ZA mode with lattice expansion. At the long-wavelength limit (near Γ point), the Grüneisen constants of LA and TA phonons are 0.18 and 4.42, respectively, which is much smaller than ZA phonon (37.24). These results are different from other 2D materials where ZA phonon owns remarkable negative Grüneisen constant, such as graphene,58 biphenylene,59 and black phosphorene.60 As the wave vector increases, the Grüneisen constant of LA phonon also remains positive, but there is a negative Grüneisen constant (−0.51) of TA phonon.Figure 6 The phonon dispersion and phonon density of states of Fe3GeS2 monolayer obtained by DS-PAW

(A) is phonon dispersion, while (B) shows the phonon density of states.

Figure 7 Grüneisen constants of phonons

(A) is Grüneisen constants of acoustic phonons along the high-symmetry path, while (B and C) are Grüneisen constants for fifteen optical phonons at the center of Brillouin zone and all phonons, respectively.

Compared with acoustic phonons, optical phonon modes can be used to measure thermal conductivity,8 evaluate the thickness of 2D materials,61 and design quantum mechanical resonator.62 Here, we also investigate and analyze the optical phonons in Fe3GeS2 monolayer. In the phonon dispersion of Fe3GeS2 monolayer, there are fifteen optical phonons named from O01 to O15, as shown in Figure 6A. We can find two sing-phonon optical modes O12 and O15 with frequencies of 9.08 THz and 11.12 THz, respectively, which is rare in other 2D materials, because PDOS always is composed of the mixed contribution from coupled phonon modes. In sing-phonon optical mode, there is no degeneracy, rendering it promising to design quantum mechanical resonators. In designing quantum mechanical resonators, how to cool the resonator to its ground state is a critical problem. Generally, the cooling temperature Tcool should be smaller than hf/kB,63 and h, f, kB are Planck’s constant, phonon frequency, and Boltzmann’s constant, respectively. For O12 and O15 modes, the corresponding maximum cooling temperatures Tcool are 432.63 K and 529.83 K, higher than room temperature and can be realized by standard cryogenic methods. Meantime, the Grüneisen constants of fifteen optical phonons at the center of Brillouin zone (Γpoint) were also calculated by Equation 4, as shown in Figure 7B. Obviously, the O15 mode owns a negative Grüneisen constant of −0.4, suggesting a frequency softening with lattice expansion. The largest Grüneisen constant of 5.85 occurs at O04 mode, but is still smaller than the ZA phonon at the long-wavelength limit.

In 2D devices, the difference in TEC between the substrates and 2D material can induce strain inevitably, and this strain increases with the temperature. At high temperature, this induced strain could destroy the geometrical structure of 2D devices, leading to the performance degradation of 2D devices. Based on these Grüneisen constants, the TEC α can be calculated by64:(Equation 5) α=kBV0·B2D∑q,sγq,s(hωq,s/kBT)2exp(hωq,s/kBT)(exp(hωq,s/kBT)−1)2,

where V0 and B2D are the volume of Fe3GeS2 primitive cell and the bulk modulus of Fe3GeS2 monolayer. The bulk modulus of Fe3GeS2 monolayer is calculated as 2,165.67 N/m by using the Yang’s modulus and Poisson ratio in a study by Long and Yang.38 The calculated TEC of Fe3GeS2 monolayer is plotted in Figure 8. At 20 K, the TEC contributed by acoustic phonons is 1.39 × 10−5 1/K, which is 97.89% of that considering the contribution from both acoustic and optical phonons. As temperature increases to 300 K, the TEC by only considering the contribution from acoustic phonons is 2.80 × 10−5 1/K. Meantime, the TEC by considering the contribution from both acoustic and optical phonons is up to 11.84 × 10−5 1/K whose 76.35% comes from the optical phonons. Obviously, the contribution from optical phonons grows remarkably with temperature, because the Grüneisen constants of optical phonons are comparable to acoustic phonons, as shown in Figure 7C. With temperature increasing, more and more optical phonons are excited, and then the anharmonic interaction between optical phonons becomes strengthening. Therefore, the contribution from optical phonons to TEC of Fe3GeS2 monolayer is enhanced significantly, which has also been observed in SiP2 monolayer.65Figure 8 The calculated thermal expansion coefficient (TEC) of Fe3GeS2 monolayer

Conclusions

According to first-principles calculations implemented in Device Studio program, we find the magnetic ground state of Fe3GeS2 monolayer is the FM, and this FM state is robust to the −1.5%–1.3% biaxial strain. Meantime, we also find the out-of-plane magnetic anisotropy is dominated by dipolar interaction in Fe3GeS2 monolayer. This out-of-plane magnetic anisotropy is beneficial for Fe3GeS2 monolayer to resist thermal fluctuation and sustain long-range FM order. In phonon dispersion, there are two single-phonon modes with frequencies of 9.08 THz and 11.12 THz, which can be used to design quantum mechanical resonator with a wide cool temperature range. Besides, the Grüneisen constants of all phonons and TEC of FM Fe3GeS2 monolayer were calculated, which highlights the contribution from the anharmonic interaction of optical phonons to the TEC. Our study provides theoretical support for the application of FM Fe3GeS2 monolayer in future spintronic devices.

Limitations of the study

This study implemented first-principles calculations by Device Studio program to investigate the magnetic anisotropy and phononic properties of ideal 2D FM Fe3GeS2 monolayer with high Curie temperature of 630 K, and found the out-of-plane magnetic anisotropy dominated by dipolar interaction, two single-phonon modes with frequencies of 9.08 THz and 11.12 THz, and the large contribution from optical phonons to thermal expansion.

Nonetheless, this study is not without its limitations that warrant further attention. Generally, the first-principles calculation and Monte Carlo simulation often overestimate the magnetic phase-transition temperature of 2D magnets, underscoring the necessity for additional experimental corroboration. The largest limitation of the study is that 2D FM Fe3GeS2 has not been successfully prepared, and experimental corroboration of its interesting properties cannot be carried out. This limitation presents a pivotal direction for future research endeavors on 2D FM Fe3GeS2.

Resource availability

Lead contact

Further information and requests for resources and reagents should be directed to and will be fulfilled by the lead contact, Ke Wang (kewang@xupt.edu.cn).

Materials availability

This study did not generate new unique materials.

Data and code availability

• Data reported in this paper will be shared by the lead contact upon request.

• This paper does not report original codes.

• Any additional information required to reanalyze the data reported in this paper is available from the lead contact upon request.

Acknowledgments

Y.W. acknowledges the 10.13039/501100007128 Natural Science Foundation of Shaanxi Province (2024JC-YBMS-325 ). K.W. acknowledges the support from 10.13039/501100001809 National Natural Science Foundation of China (NNSFC) (12204373 ). We gratefully acknowledge HZWTECH for providing computation facilities.

Author contributions

Conceptualization, Y.W. and H.L.; methodology, Y.W. and K.W.; investigation, Y.W. and K.W.; software, K.W.; writing – original draft, Y.W.; writing – review and editing, K.W. and H.L.; funding acquisition and project administration, Y.W. and K.W.

Declaration of interests

The authors declare no competing interests.

STAR★Methods

Key resources table

REAGENT or RESOURCE	SOURCE	IDENTIFIER	
Software and algorithms	
	
Device Studio (DS-PAW) 2022	HONGZHIWEI TECHNOLOGY	https://cloud.hzwtech.com/web/product-service?id=6	

Experimental model and study participant details

Our study does not use experimental models typical in the life sciences.

Method details

The detailed numerical setting in first-principles calculations have been presented in computational details.

Quantification and statistical analysis

Our study does not include statistical analysis or quantification.
==== Refs
References

1 Gong C. Li L. Li Z. Ji H. Stern A. Xia Y. Cao T. Bao W. Wang C. Wang Y. Discovery of intrinsic ferromagnetism in two-dimensional van der Waals crystals Nature 546 2017 265 269 28445468
2 Huang B. Clark G. Navarro-Moratalla E. Klein D.R. Cheng R. Seyler K.L. Zhong D. Schmidgall E. McGuire M.A. Cobden D.H. Layer-dependent ferromagnetism in a van der Waals crystal down to the monolayer limit Nature 546 2017 270 273 28593970
3 Cai X. Song T. Wilson N.P. Clark G. He M. Zhang X. Taniguchi T. Watanabe K. Yao W. Xiao D. Atomically Thin CrCl3: An In-Plane Layered Antiferromagnetic Insulator Nano Lett. 19 2019 3993 3998 31083954
4 Fei Z. Huang B. Malinowski P. Wang W. Song T. Sanchez J. Yao W. Xiao D. Zhu X. May A.F. Two-dimensional itinerant ferromagnetism in atomically thin Fe3GeTe2 Nat. Mater. 17 2018 778 782 30104669
5 Sun X. Li W. Wang X. Sui Q. Zhang T. Wang Z. Liu L. Li D. Feng S. Zhong S. Room temperature ferromagnetism in ultra-thin van der Waals crystals of 1T-CrTe2 Nano Res. 13 2020 3358 3363
6 Deng Y. Yu Y. Shi M.Z. Guo Z. Xu Z. Wang J. Chen X.H. Zhang Y. Quantum anomalous Hall effect in intrinsic magnetic topological insulator MnBi2Te4 Science 367 2020 895 900 31974160
7 You J.-Y. Chen C. Zhang Z. Sheng X.-L. Yang S.A. Su G. Two-dimensional Weyl half-semimetal and tunable quantum anomalous Hall effect Phys. Rev. B 100 2019 064408
8 Kargar F. Coleman E.A. Ghosh S. Lee J. Gomez M.J. Liu Y. Magana A.S. Barani Z. Mohammadzadeh A. Debnath B. Phonon and Thermal Properties of Quasi-Two-Dimensional FePS3 and MnPS3 Antiferromagnetic Semiconductors ACS Nano 14 2020 2424 2435 31951116
9 He J. Ding G. Zhong C. Li S. Li D. Zhang G. Cr2TiC2-based double MXenes: novel 2D bipolar antiferromagnetic semiconductor with gate-controllable spin orientation toward antiferromagnetic spintronics Nanoscale 11 2019 356 364
10 Kumar H. Frey N.C. Dong L. Anasori B. Gogotsi Y. Shenoy V.B. Tunable Magnetism and Transport Properties in Nitride MXenes ACS Nano 11 2017 7648 7655 28558192
11 Marfoua B. Khan I. Hong J. Giant spin Seebeck effect in two-dimensional ferromagnetic CrI3 monolayer Nanotechnology 31 2020 455404
12 Liu W. Wang Y. Han Y. Tong W. Fan J. Pi L. Hao N. Zhang L. Zhang Y. Anisotropic magnetoresistance behaviors in the layered ferromagnetic Cr2Ge2Te6 J. Phys. D Appl. Phys. 53 2020 025101
13 Li J. Li Y. Du S. Wang Z. Gu B.-L. Zhang S.-C. He K. Duan W. Xu Y. Intrinsic magnetic topological insulators in van der Waals layered MnBi2Te4-family materials Sci. Adv. 5 2019 eaaw5685
14 Heinz B. Brächer T. Schneider M. Wang Q. Lägel B. Friedel A.M. Breitbach D. Steinert S. Meyer T. Kewenig M. Propagation of Spin-Wave Packets in Individual Nanosized Yttrium Iron Garnet Magnonic Conduits Nano Lett. 20 2020 4220 4227 32329620
15 Kent A.D. Worledge D.C. A new spin on magnetic memories Nat. Nanotechnol. 10 2015 187 191 25740126
16 Cortie M.B. Cortie D.L. Timchenko V. Heat transfer from nanoparticles for targeted destruction of infectious organisms Int. J. Hyperther. 34 2018 157 167
17 Arnold D.P. Review of Microscale Magnetic Power Generation IEEE Trans. Magn. 43 2007 3940 3951
18 Olsen T. Magnetic anisotropy and exchange interactions of two-dimensional FePS3, NiPS3 and MnPS3 from first principles calculations J. Phys. D Appl. Phys. 54 2021 314001
19 Wang K. Ren K. Hou Y. Cheng Y. Zhang G. Physical insights into enhancing magnetic stability of 2D magnets J. Appl. Phys. 133 2023 110902
20 Wu L. Zhou L. Zhou X. Wang C. Ji W. In-plane epitaxy-strain-tuning intralayer and interlayer magnetic coupling in CrSe2 and CrTe2 monolayers and bilayers Phys. Rev. B 106 2022 L081401
21 Gu Y. Zhang S. Zou X. Tunable magnetism in layered CoPS3 by pressure and carrier doping Sci. China Mater. 64 2020 673 682
22 Liu M. Zhang L. Liu J. Wan T.L. Du A. Gu Y. Kou L. Density Functional Theory Studies on Magnetic Manipulation in NiI2 Layers ACS Appl. Electron. Mater. 5 2023 920 927
23 Ji S. Wang Y. Yao R. Quan C. Yang J. Li X. Dimension-dependent magnetic anisotropy for tunable anomalous Hall effect in transition metal dichalcogenides Phys. Rev. B 107 2023 035438
24 Niu K. Qiu G. Wang C. Li D. Niu Y. Li S. Kang L. Cai Y. Han M. Lin J. Self-Intercalated Magnetic Heterostructures in 2D Chromium Telluride Adv. Funct. Mater. 33 2023 2208528
25 Dey D. Ray A. Yu L. Intrinsic ferromagnetism and restrictive thermodynamic stability in MA2N4 and Janus VSiGeN4 monolayers Phys. Rev. Mater. 6 2022 L061002
26 Tang C. Ostrikov K.K. Sanvito S. Du A. Prediction of room-temperature ferromagnetism and large perpendicular magnetic anisotropy in a planar hypercoordinate FeB3 monolayer Nanoscale Horiz. 6 2021 43 48 33237052
27 Zhang S. Liang X. Zhao H. Chen Y. He Q. Liu J. Lv L. Yang J. Wu H. Chen L. Tuning the magnetic properties of Fe3GeTe2 by doping with 3d transition-metals Phys. Lett. 396 2021 127219
28 Zhang L. Tang C. Sanvito S. Gu Y. Du A. Hydrogen-Intercalated 2D Magnetic Bilayer: Controlled Magnetic Phase Transition and Half-Metallicity via Ferroelectric Switching ACS Appl. Mater. Interfaces 14 2022 1800 1806 34962753
29 Guo Y. Liu N. Zhao Y. Jiang X. Zhou S. Zhao J. Enhanced Ferromagnetism of CrI3 Bilayer by Self-Intercalation Chin. Phys. Lett. 37 2020 107506
30 Popov Z.I. Sukhanova E.V. Kvashnin D.G. Metallocene inspired 2D metal intercalated carbon allotropes: Stability and properties via DFT calculations Carbon 184 2021 714 720
31 Cheng H.-X. Zhou J. Wang C. Ji W. Zhang Y.-N. Nonvolatile electric field control of magnetism in bilayer CrI3 on monolayer In2Se3 Phys. Rev. B 104 2021 064443
32 Zhuo W. Lei B. Wu S. Yu F. Zhu C. Cui J. Sun Z. Ma D. Shi M. Wang H. Manipulating Ferromagnetism in Few-Layered Cr2Ge2Te6 Adv. Mater. 33 2021 2008586
33 Rhone T.D. Bhattarai R. Gavras H. Lusch B. Salim M. Mattheakis M. Larson D.T. Krockenberger Y. Kaxiras E. Artificial Intelligence Guided Studies of van der Waals Magnets Adv. Theory Simul. 6 2023 2300019
34 Lu H. Yang T. Sun Z. Robertson J. Zhao W. Automatic Calculation of the Transition Temperatures for two-dimensional Heisenberg type Magnets Preprint at arXiv 2023 10.48550/arXiv:2312.04497
35 Kothakonda M. Zhu Y. Guan Y. He J. Kidd J. Zhang R. Ning J. Gopalan V. Xie W. Mao Z. Sun J. High-Throughput Screening Assisted Discovery of a Stable Layered Anti-Ferromagnetic Semiconductor: CdFeP2Se6 Adv. Funct. Mater. 33 2023 2210965
36 Wang H. Wen Y. Zhao X. Cheng R. Yin L. Zhai B. Jiang J. Li Z. Liu C. Wu F. He J. Heteroepitaxy of 2D CuCr2Te4 with Robust Room-temperature Ferromagnetism Adv. Mater. 35 2023 2211388
37 Zhang L. Huang X. Dai H. Wang M. Cheng H. Tong L. Li Z. Han X. Wang X. Ye L. Han J. Proximity-Coupling-Induced Significant Enhancement of Coercive Field and Curie Temperature in 2D van der Waals Heterostructures Adv. Mater. 32 2020 2002032
38 Long D.-B. Yang L.-M. Transmutation engineering makes a large class of stable and exfoliable A3BX2 compounds with exceptional high magnetic critical temperatures and exotic electronic properties ACS Appl. Mater. Interfaces 15 2023 24549 24569 37171271
39 Hongzhiwei Technology Device Studio, Version 2022A, China https://cloud.hzwtech.com/web/product-service?id=6 2022
40 Blöchl P.E. Projector augmented-wave method Phys. Rev. B 50 1994 17953 17979
41 Perdew J.P. Burke K. Ernzerhof M. Generalized gradient approximation made simple Phys. Rev. Lett. 77 1996 3865 3868 10062328
42 Kresse G. Hafner J. Ab initio simulation of the metal/nonmetal transition in expanded fluid mercury Phys. Rev. B 55 1997 7539 7548
43 Gonze X. Adiabatic density-functional perturbation theory Phys. Rev. 52 1995 1096 1114
44 Li X. Lü J.-T. Zhang J. You L. Su Y. Tsymbal E.Y. Spin-Dependent Transport in van der Waals Magnetic Tunnel Junctions with Fe3GeTe2 Electrodes Nano Lett. 19 2019 5133 5139 31276417
45 Jiang M.-C. Guo G.-Y. Large magneto-optical effect and magnetic anisotropy energy in two-dimensional metallic ferromagnet Fe3GeTe2 Phys. Rev. B 105 2022 014437
46 Yang X. Zhou X. Feng W. Yao Y. Strong magneto-optical effect and anomalous transport in the two-dimensional van der Waals magnets FenGeTe2 (n=3, 4, 5) Phys. Rev. B 104 2021 104427
47 Lado J.L. Fernández-Rossier J. On the origin of magnetic anisotropy in two dimensional CrI3 2D Mater. 4 2017 035002
48 Shen K. Magnon Spin Relaxation and Spin Hall Effect Due to the Dipolar Interaction in Antiferromagnetic Insulators Phys. Rev. Lett. 124 2020 077201
49 Politi P. Pini M.G. Dipolar interaction between two-dimensional magnetic particles Phys. Rev. B 66 2002 214414
50 Young J.T. Bienias P. Belyansky R. Kaufman A.M. Gorshkov A.V. Asymmetric Blockade and Multiqubit Gates via Dipole-Dipole Interactions Phys. Rev. Lett. 127 2021 120501
51 Utesov O.I. Thermodynamically stable skyrmion lattice in a tetragonal frustrated antiferromagnet with dipolar interaction Phys. Rev. B 103 2021 064414
52 Bogdanov A.N. Panagopoulos C. Physical foundations and basic properties of magnetic skyrmions Nat. Rev. Phys. 2 2020 492 498
53 Utesov O.I. Syromyatnikov A.V. Cascades of phase transitions in spiral magnets caused by dipolar forces Phys. Rev. B 95 2017 214420
54 Utesov O.I. Syromyatnikov A.V. Phase competition in frustrated anisotropic antiferromagnet in strong magnetic field J. Magn. Magn Mater. 527 2021 167732
55 Freedman J.P. Leach J.H. Preble E.A. Sitar Z. Davis R.F. Malen J.A. Universal phonon mean free path spectra in crystalline semiconductors at high temperature Sci. Rep. 3 2013 2963 24129328
56 Regner K.T. Sellan D.P. Su Z. Amon C.H. McGaughey A.J.H. Malen J.A. Broadband phonon mean free path contributions to thermal conductivity measured using frequency domain thermoreflectance Nat. Commun. 4 2013 1640 23535661
57 Cai Y. Lan J. Zhang G. Zhang Y.-W. Lattice vibrational modes and phonon thermal conductivity of monolayer MoS2 Phys. Rev. B 89 2014 035438
58 Kong B.D. Paul S. Nardelli M.B. Kim K.W. First-principles analysis of lattice thermal conductivity in monolayer and bilayer graphene Phys. Rev. B 80 2009 033406
59 Wang K. Ren K. Zhang D. Cheng Y. Zhang G. Phonon properties of biphenylene monolayer by first-principles calculations Appl. Phys. Lett. 121 2022 042203
60 Cai Y. Ke Q. Zhang G. Feng Y.P. Shenoy V.B. Zhang Y.-W. Giant Phononic Anisotropy and Unusual Anharmonicity of Phosphorene: Interlayer Coupling and Strain Engineering Adv. Funct. Mater. 25 2015 2230 2236
61 Ghosh S. Bao W. Nika D.L. Subrina S. Pokatilov E.P. Lau C.N. Balandin A.A. Dimensional crossover of thermal transport in few-layer graphene Nat. Mater. 9 2010 555 558 20453845
62 Zhou H. Cai Y. Zhang G. Zhang Y.-W. Unusual phonon behavior and ultra-low thermal conductance of monolayer InSe Nanoscale 10 2017 480 487 29231225
63 O’Connell A.D. Hofheinz M. Ansmann M. Bialczak R.C. Lenander M. Lucero E. Neeley M. Sank D. Wang H. Weides M. Quantum ground state and single-phonon control of a mechanical resonator Nature 464 2010 697 703 20237473
64 Huang L.F. Gong P.L. Zeng Z. Correlation between structure, phonon spectra, thermal expansion, and thermomechanics of single-layer MoS2 Phys. Rev. B 90 2014 045409
65 Hou Y. Ren K. Wei Y. Cui Z. Wang K. Anisotropic phonon properties in SiP2 monolayer: A first-principles study Vacuum 222 2024 113094
