企业官网建设流程全解析

Li 插层 Fe3GaTe2中的显著磁性增强

J. APPL. PHYS. 2026

Li 插层 Fe₃GaTe₂中的显著磁性增强

Significantly Improved Curie Temperature in van der Waals Fe₃GaTe₂ by Li-Ion Intercalation

导读 导读：本文通过 DFT 计算和蒙特卡洛模拟，预测 Li 离子插层可将范德华磁体 Fe₃GaTe₂ 的居里温度从 365 K 大幅提升至 770 K——超过两倍的增强。电子掺杂同时增强了层间铁磁耦合（自旋极化界面电子中介）和层内自旋交换（Fe 3d 轨道占据优化），为高居里温度 vdW 磁体的设计提供了理论指导。

【一、前言背景】

1.1 范德华磁性材料的发展历程

vdW 层状磁性材料是凝聚态物理近十年的研究热点——从 CrI₃ 单层铁磁性的发现（Nature 2017）开始，该领域迅速发展。

T_C 演进路线：CrI₃（45-61 K, 2015-2017）→ Fe₃GeTe₂（~100-300 K, 2018-2020）→ Fe₄GeTe₂/Fe₅GeTe₂（~250-310 K, 2019-2020）→ Fe₃GaTe₂（~350-380 K, 2022-2024）。

短短数年内，本征 T_C 提升了近一个数量级——这得益于新材料的不断发现和合成技术的进步。

目前已形成庞大的 vdW 磁性材料家族——涵盖 Cr 基、Fe 基、Mn 基等多种体系，为自旋电子学器件提供了丰富的候选材料。

1.2 Fe₃GaTe₂：当前 T_C 最高的 vdW 磁体之一

Fe₃GaTe₂（FGaT）具有六方晶体结构，空间群 P6₃/mmc——每单胞包含两个由 vdW 间隙分隔的单层。

每单层由五个原子亚层组成：Te-FeI-(FeIIGa)-FeI-Te，其中 FeI 和 FeII 为两个不等价 Fe 位点。

实验测量：FGaT 体相 T_C = 350-380 K，单层 T_C ≈ 240-260 K——已接近或超过室温，是 vdW 磁体中的最高水平之一。

晶格常数：a = 4.06-4.09 Å, c = 16.07-16.15 Å——DFT 计算与实验吻合良好，为后续插层研究提供了可靠的结构基础。

1.3 T_C 对器件应用的决定性作用

对于实际器件，T_C 和操作温度之间需要足够的余量以确保磁体维持稳定功能——这是决定器件工作温度上限的关键性能指标。

作为参考：商用 NdFeB 磁体的 T_C ≈ 300°C，但稳定操作温度通常仅限 ~100°C——安全余量约 200°C。

vdW 磁体若要应用于自旋阀、磁隧道结（MTJ）和自旋轨道转矩（SOT）器件，T_C 必须远超室温。

因此，探索大幅提升 T_C 的策略，对于实现基于 vdW 磁体的热稳定室温纳米器件至关重要。

1.4 Li 离子插层策略的历史与成效

Li 离子插层（或注入）是一种广泛应用于 vdW 层状材料的组分调控策略——在磁性增强方面已展现出显著效果。

里程碑一：Deng 等（Nature 2018）实验将 Li 离子嵌入 Fe₃GeTe₂，T_C 从 ~100 K 提升至 >300 K——提升约 3 倍。

里程碑二：Wang 等（ACS Nano 2024）实验将 Li 离子嵌入 CrI₃，T_C 从 ~50 K 提升至 >400 K——提升约 8 倍。

DFT 计算揭示（Huang 等, Appl. Phys. Lett. 2021）：T_C 增强源于电子掺杂——同时强化层间和层内铁磁（FM）耦合。

1.5 Li 插层增强磁性的物理机制

电子掺杂效应：Li 的 2s 电子注入 FGaT 导带——提供界面自旋极化中介载流子，增强跨 vdW 间隙的层间 FM 耦合（类似 RKKY 机制）。

层内电荷重新分布：电子掺杂改变 Fe 3d 轨道占据——部分电子从自旋向下通道翻转到自旋向上通道，增大了局域磁矩和交换劈裂。

协同增强：层间增强的 FM 耦合通过分子场效应进一步稳定层内 FM 序——两个机制相互促进，非简单相加。

Te 介导的间接交换：Li 向 Te 捐赠电子后，Te 对周围 Fe 的电子需求减少——Fe-Fe 直接交换和 Fe-Te-Fe 超交换均被增强。

1.6 本文研究动机与核心目标

既然 Li 插层对低 T_C 的 Fe₃GeTe₂ 和 CrI₃ 有效，那么对已有更高本征 T_C 的 FGaT 进行 Li 插层——有望实现远超室温的稳定铁磁性。

核心目标：通过 DFT 计算构建 Li 插层 FGaT（Li-FGaT）晶格模型——预测 T_C 可达 770 K，并揭示电子掺杂增强磁性的微观机制。

本文还考察了双层模型（不同 Li 浓度和空间分布）的磁性——验证 Li 插层策略在厚度、浓度和空间分布上的普适性。

最终目标：为实验提供可测试的明确预测（Li-FGaT, T_C ~ 770 K），并为高 T_C vdW 磁体的电子掺杂设计提供理论框架。

【二、研究方法】

2.1 DFT 计算框架

采用 VASP（Vienna Ab initio Simulation Package）程序包进行第一性原理计算——基于平面波基组和赝势方法。

交换关联泛函：PBE（Perdew-Burke-Ernzerhof）形式的 GGA 泛函——广泛用于过渡金属化合物的电子结构计算。

范德华相互作用：DFT-D3 方法（Grimme 等, J. Chem. Phys. 2010）——用于准确描述 vdW 层状材料中的层间弱相互作用。

平面波截断能：400 eV——足以收敛 Fe、Ga、Te 等元素的电子波函数。力收敛标准：0.005 eV/Å——确保结构优化充分。

2.2 Hubbard U 校正与参数选择

Hubbard U = 1.3 eV 用于 Fe-3d 电子的关联效应修正（Dudarev 等, Phys. Rev. B 1998）——采用 DFT+U 旋转不变方案。

U 值的选择依据：该 U 值使计算得到的 FGaT 体相 T_C（365 K）与实验测量值（350-380 K）吻合良好——验证了参数设置的可靠性。

电子自洽收敛标准：10⁻⁶ eV——确保总能计算精度足以分辨不同磁序之间的微小能量差异（~meV 量级）。

DFT+U 方法在过渡金属氧化物和硫属化合物中广泛应用——能够有效修正 3d 电子的自相互作用误差。

2.3 k 点网格与计算精度

Monkhorst-Pack（Γ 中心）k 点网格：体相单胞 24×24×4，2×2×1 超胞 12×12×4——密集网格确保总能收敛。

单层/双层板模型：24×24×1（单胞）或 12×12×1（2×2×1 超胞）——板模型仅在面向有周期性，z 方向无需 k 点采样。

声子谱计算：使用 PHONOPY 代码基于有限位移法——验证晶格动力学稳定性（无虚频）。

AIMD 模拟：4×4×1 超胞（224 原子），正则系综，Nosé-Hoover 恒温器 500 K——时间步长 1 fs，模拟 5 ps。

2.4 交换参数提取：能量映射方法

采用经典 Heisenberg 自旋哈密顿量描述 Fe 自旋晶格——包含交换项 J_ij 和单离子各向异性（SIA）项 A。

能量映射方法：构建 2×2×1 超胞中六种不同磁序（MO-1 至 MO-6），计算每种磁序的 DFT 总能——建立总能与交换参数的线性方程组（Eqs. 2-7）。

五个交换参数：J₁（FeI-FeI 最近邻）、J₂（FeII-FeII 最近邻）、J₃（FeI-FeII 最近邻）、J₄（FeI-FeI 垂直）、J₅（层间 FeI-FeI）。

SIA 参数 A 通过 MCA 能量提取：E_MCA = E[001] − E[100] = N_Fe·A——其中 N_Fe 为单胞中 Fe 位点数。

2.5 蒙特卡洛模拟设置

使用 MCsolver 软件（Liu 等, J. Phys. Chem. Lett. 2020）——基于 Metropolis 算法的经典自旋蒙特卡洛模拟。

体相系统：60×60×20 超胞——从 10 K 起以 10 K 步长升温，每温度 8×10⁴ 步热平衡 + 1.6×10⁵ 步统计平均。

单层/双层板系统：60×60×1 超胞——自旋涨落显著，需更长的平衡和采样时间（8×10⁵ + 1.6×10⁶ 步/温度）。

初始磁构型：面外 FM 态——通过 Curie-Bloch 方程 M(T) = (1 − T/T_C)^β 拟合 M-T 曲线提取 T_C。

2.6 磁晶各向异性与稳定性验证

磁晶各向异性（MCA）能量：E_MCA = E[001] − E[100]——计算包含自旋轨道耦合（SOC）的 DFT 总能差。

FGaT 体相：E_MCA = −2.82 meV/单胞，A = −0.47 meV/Fe——负值表示垂直磁各向异性（PMA），与实验观察一致。

声子谱验证：无虚频模式 → 晶格动力学稳定——AIMD 验证（500 K, 5 ps）→ 晶格框架保持完好，Li 离子留在最优位点。

Li 插层能：E_int = E(Li-FGaT) − E(FGaT) − E(Li) = −1.69 eV——负值表明 Li 离子在八面体位点能量稳定。

【三、实验结果】

图 1：FGaT 体相。(a) 侧视图和 (b) 俯视图——红色框标示单胞；(c) 五种交换耦合 J₁-J₅ 的示意图；(d) 用于计算交换参数的六种磁序（MO-1 至 MO-6）——深蓝/浅蓝球为自旋向上 FeI/FeII，红/浅粉球为自旋向下 FeI/FeII；(e) 计算得到的五种交换参数；(f) MC 模拟的 M-T 结果——蓝方块/绿圆点为 FGaT 体相/单层数据，红/橙曲线为 Curie-Bloch 拟合。

3.1 FGaT 原始晶体结构与磁性质

FGaT 具有六方结构（P6₃/mmc）——每单胞包含两个单层（2 f.u.），由 vdW 间隙分隔。每单层五个原子亚层：Te-FeI-(FeIIGa)-FeI-Te。

两个不等价 Fe 位点：FeI 位于外层亚层（与 Te 直接键合），FeII 位于中心亚层（与 Ga 共占据）——每个亚层原子占据三个三角形亚晶格位点（site-A/B/C）之一。

计算晶格常数：a = 4.06 Å, c = 16.15 Å——与实验值（a = 4.09 Å, c = 16.07 Å, Lee 等 Nano Lett. 2023）吻合良好。

铁磁性起源于自旋极化 Fe-3d 电子——FGaT 呈 FM 金属性电子结构，每单胞磁矩 13.44 μB。

3.2 FGaT 的交换参数 J₁-J₅

通过能量映射方法，在 2×2×1 超胞中构建六种磁序（MO-1 至 MO-6），利用 DFT 总能建立线性方程组（Eqs. 2-7）求解交换参数。

J₁（FeI-FeI NN）和 J₂（FeII-FeII NN）较小——同一亚层内的最近邻交换对 FM 序的贡献有限。

J₃（FeI-FeII NN）和 J₄（FeI-FeI 垂直）建立层内 FM 序——是主导的层内交换参数。

J₅（层间 FeI-FeI 垂直）对应层间 FM 耦合——证实 FGaT 体相具有层间 FM 耦合，与前期研究一致（Li 等, Appl. Phys. Lett. 2023）。

3.3 FGaT 的 MC 模拟验证

利用 DFT 提取的交换参数和 SIA 参数，进行 MC 模拟获得 M-T 曲线——FGaT 体相 T_C = 365 K，单层 T_C = 260 K。

Curie-Bloch 拟合：M(T) = (1 − T/T_C)^β，β = 0.32——对体相和单层均提供良好拟合。

体相 T_C = 365 K 与实验值（350-380 K, Zhang 等 Nat. Commun. 2022; Wang 等 npj 2D Mater. Appl. 2024）吻合——验证了计算方法和参数设置的可靠性。

单层 T_C = 260 K 与实验值（~240 K, Wang 等）基本一致——进一步确认了 DFT+U（U = 1.3 eV）方案的准确性。

Heisenberg 自旋哈密顿量：J_ij 为交换耦合常数，S_i 为经典自旋单位矢量，A 为单离子各向异性参数——求和遍历所有 Fe 位点对。

能量映射方程组（示例）：E_bulk(1) = E₀ + 48J₁ + 24J₂ + 48J₃ + 8J₄ + 8J₅——六种磁序对应六个方程，联立求解 J₁-J₅。

图 2：Li-FGaT 体相。(a) 侧视图和俯视图——粉色球（site-A）为插层 Li 离子；(b) 声子能带结构——无虚频，确认动力学稳定；(c) AIMD 模拟中温度和能量随时间的演化——时间步长 1 fs；(d) 5 ps AIMD 模拟后的晶格结构快照——Li 周围 Te 配位结构显示为白色八面体，FeII（或 Ga）周围 FeI 配位结构显示为浅蓝（或浅绿）三棱柱。

3.4 Li-FGaT 晶格模型构建

每单胞两个 vdW 间隙中各插入一个 Li 离子——Li 位于 site-A（六个界面 Te 原子包围的八面体中心），此为最有利的插层位点。

验证：site-A 比 site-B（等效 site-C）能量上更优 0.49 eV/Li——符合已有研究对金属离子在 FGaT 中插层位点的判断（Huang 等, Appl. Phys. Lett. 2022）。

Li 插层引起 ~3% 的 c 轴晶格膨胀：a = 4.13 Å, c = 16.64 Å——vdW 间隙从 2.94 Å 增至 3.14 Å（+0.20 Å），单层厚度从 5.13 Å 增至 5.18 Å。

vdW 间隙膨胀是 c 轴参数增加的主要原因——与预期一致，Li 离子嵌入 vdW 间隙自然推开相邻 Te 层。

3.5 Li-FGaT 的结构稳定性验证

声子谱计算：无虚频模式——Li-FGaT 晶格在动力学上稳定，无结构畸变或局部离子位移的趋势。

Li 插层能：E_int = −1.69 eV——负值足够大，表明 Li 离子在八面体位点能量稳定，无聚集驱动力。

AIMD 模拟（500 K, 5 ps）：体系能量和温度在稳定范围内波动——晶格框架保持完好，Li 离子留在最优插层位点。

Li 周围保持 Te 八面体配位，FeII（或 Ga）周围保持 FeI 三棱柱配位——说明 Li-FGaT 在远超室温的温度下（500 K）维持结构完整性。

Li 插层能定义：E_int = E(Li-FGaT) − E(FGaT) − E(Li)——负值越大，Li 越稳定。

电子密度差定义：Δρ_1D = ρ_1D(Li-FGaT) − ρ_1D(Li⁰-FGaT)——Li⁰-FGaT 指移除 Li 但保持原子位置不变的 FGaT 晶格。

图 3：Li-FGaT 体相。(a) 元素分辨投影态密度（DOS）——Fermi 能级设为 0 eV；(b) 计算得到的五种交换参数（J₁-J₅）——水平轴在 −140 至 −90 处断开；(c) MC 模拟的 M-T 结果——粉色方块为 MC 结果，蓝色曲线为 Curie-Bloch 拟合（β=0.35）。

3.6 Li-FGaT 的电子结构与磁矩

Li-FGaT 呈 FM 金属性电子结构——Fermi 能级附近态密度主要由 Fe-3d 轨道贡献，Te-5p 轨道参与杂化。

每单胞磁矩：14.91 μB——比原始 FGaT（13.44 μB）增加 1.47 μB，说明 Li 插层增强了体系的净磁化强度。

磁矩增加主要来自 FeII 位点的自旋翻转——部分电子从自旋向下通道转移到自旋向上通道（Stoner 型机制）。

Li 本身非磁性——磁矩增加完全源于电子掺杂对 FGaT 电子结构的调控效应。

3.7 Li-FGaT 的交换参数增强

与原始 FGaT 类似，Li-FGaT 的层内 FM 序由 J₃ 和 J₄ 决定，J₁ 和 J₂ 很小，层间为 FM 耦合（J₅ > 0）。

关键发现：主导交换参数 J₃、J₄、J₅ 均显著增强——Li 插层使层内和层间 FM 耦合同时大幅增强。

MCA 计算：Li-FGaT 保持 PMA——但 SIA 参数 A 减小至 −0.167 meV/Fe（原始 FGaT 为 −0.47 meV/Fe）。

A 减小归因于插层诱导的巡游电子对面向电子轨道 PMA 贡献的抑制效应（Mo 等, J. Appl. Phys. 2024）——这是电子掺杂的副作用。

3.8 Li-FGaT 的 T_C = 770 K

基于 DFT 计算参数进行 MC 模拟——Li-FGaT 体相的 T_C = 770 K，是原始 FGaT（365 K）的 2.1 倍。

Curie-Bloch 拟合：β = 0.35——与原始 FGaT（β = 0.32）略有不同，反映 Li 插层后自旋系统的维度特征变化。

770 K 远超已有 vdW 磁体的 T_C 记录——Li-Fe₃GeTe₂（>300 K）和 Li-CrI₃（>400 K）均无法达到此水平。

这一预测为实验提供了明确的可测试目标——如果实验验证，770 K 将创造 vdW 磁体 T_C 的新纪录。

Curie-Bloch 拟合：M(T) = (1 − T/T_C)^β——β=0.32（FGaT）或 β=0.35（Li-FGaT）。

MCA 能量定义：E_MCA = E[001] − E[100] = N_Fe·A——负值表示 PMA，N_Fe 为单胞中 Fe 位点数。

图 4：Li-FGaT 的电荷分布。(a) Li 插层诱导电子沿 c 轴的分布 Δρ_1D——电子主要局域在界面 Te 亚层之间的 vdW 间隙；(b) Bader 分析——每个 Li 原子捐赠 ~0.84 e 给 FGaT，其中 ~93% 被界面 Te 原子获取（~0.39 e/Te）；(c) 自旋极化的界面电子——每个 Te 贡献 0.13 μB；(d) (110) 面内插层诱导的电子密度变化——黄色/蓝色区域为电子增益/损失。

3.9 电子掺杂增强磁性的微观机制

电子掺杂分布：每个 Li 捐赠 ~0.84 e 给 FGaT——~93% 被界面 Te 原子获取（~0.39 e/Te），Fe 的电子占据数变化很小。

自旋极化界面电子：Te 获取的掺杂电子受金属性电子结构自旋极化影响——每个 Te 贡献 0.13 μB，为层间自旋交换提供中介载流子。

FeII 磁矩增加：尽管 Fe 的电子占据数变化很小，FeII 的磁矩显著增加——部分电子从自旋向下通道翻转到自旋向上通道。

层内电子重新分布：Te 获得额外电子后，减少对周围 Fe 的电子抽取——Fe-Fe 直接交换和 Fe-Te-Fe 超交换通道均被增强。

图 5：MC 模拟的 M-T 结果——FGaT 双层、Li₁-(FGaT)₂ 和 Li₁/₄-(FGaT)₂ 模型的对比。空心符号为 MC 结果，实线为 Curie-Bloch 拟合（β=0.32）。

【四、对比分析】

4.1 Li-FGaT vs 其他 Li 插层 vdW 磁体

Li-Fe₃GeTe₂：T_C 从 ~100 K → >300 K（提升 ~3 倍）——起始 T_C 低，提升倍数大但绝对值有限（~300 K）。

Li-CrI₃：T_C 从 ~50 K → >400 K（提升 ~8 倍）——起始 T_C 极低，Li 插层效果最显著（倍数最大）。

Li-Fe₃GaTe₂（本工作）：T_C 从 365 K → 770 K（提升 ~2.1 倍）——起始 T_C 最高，提升后绝对值最高（770 K），远超已有报道。

趋势：起始 T_C 越高，Li 插层提升的倍数越小但绝对值越大——电子掺杂在已有强 FM 体系中"叠加"效果更显著，暗示交换耦合强度存在上限。

4.2 双层模型：Li 浓度效应

构建三个双层模型：原始 FGaT 双层（T_C = 300 K，与实验一致）、Li₁-(FGaT)₂（每间隙一个 Li，T_C = 490 K）、Li₁/₄-(FGaT)₂（每四个间隙一个 Li，T_C = 395 K）。

高 Li 浓度（Li₁）对应更强的 T_C 增强——490 K vs 395 K——表明磁性增强程度与 Li 掺杂浓度正相关。

即使低 Li 浓度（Li₁/₄），T_C 仍从 300 K 提升至 395 K——说明 Li 插层策略在低掺杂水平下依然有效。

双层模型更接近实际器件构型——vdW 异质结器件通常基于少层而非体相材料，验证了该策略在低维极限下的有效性。

4.3 Li 离子聚集 vs 分散排列

构造对比模型 ag-Li₁/₄-(FGaT)₂——Li 离子以聚集方式排列（与 Li₁/₄-(FGaT)₂ 的分散排列对比），浓度相同。

能量比较：ag-Li₁/₄ 能量比分散排列高 3.8 meV/Li——Li 离子倾向于分散排列，聚集在能量上不利，但可能以一定概率出现。

层间磁耦合能 E_imc：ag-Li₁/₄ 为 9.0 meV，分散 Li₁/₄ 为 11.2 meV，原始 FGaT 双层仅为 3.3 meV。

关键结论：即使在能量不利的聚集构型下，Li 插层仍能有效增强层间 FM 耦合（E_imc 是原始的 2.7 倍）——磁性增强效果具有空间分布鲁棒性。

4.4 实验验证前景

Li 离子插层实验已在 Fe₃GeTe₂ 和 CrI₃ 中成功实现——电化学插层法可在室温下将 Li⁺ 嵌入 vdW 间隙，浓度精确可控。

磁化测量（SQUID-VSM）可直接验证 T_C 提升——770 K 远超通常 SQUID 测量范围（~400 K），但可通过高温磁化率或高温霍尔效应间接验证。

XRD 可确认 Li 插层后的 c 轴晶格膨胀（~3%）——层间距增大是 Li 成功嵌入的明确标志，是实验验证的第一步。

声子谱测量（拉曼/非弹性中子散射）可验证 Li-FGaT 的动力学稳定性——与 DFT 计算的声子谱对比，确认结构完整性。

4.5 器件应用潜力

770 K 的 T_C 意味着 Li-FGaT 可在极端高温环境下保持铁磁功能——适用于航空航天（>200°C）、汽车电子（~150°C）、深井勘探等高温场景。

vdW 材料天然适合异质结集成——Li-FGaT 可与石墨烯、hBN 等 2D 材料堆叠构建自旋阀和磁隧道结（MTJ）。

Li 插层可在器件制备后通过电化学方法实现——为后处理调控器件性能提供灵活手段，类似于离子液体栅控技术。

高 T_C 确保器件在宽温度范围内的稳定操作——这是 vdW 磁体从实验室走向实际应用的关键一步。

【五、讨论】

5.1 层间 FM 耦合增强：自旋极化界面电子

Li 插层引入的掺杂电子主要局域在 vdW 间隙中的界面 Te 亚层之间——形成二维自旋极化电子气。

这些自旋极化的界面电子充当层间中介载流子——类似于 RKKY 机制中传导电子介导的间接交换，增强了跨 vdW 间隙的 Fe-Fe 交换（J₅）。

界面 Te 获取 ~0.39 e/Te 并贡献 0.13 μB/Te——自旋极化程度约为 33%，足以有效介导层间 FM 耦合。

这一机制解释了为何 Li 插层对层间 FM 耦合的增强效果显著——自旋极化界面电子直接桥接了上下两个 FGaT 单层的 Fe 自旋。

5.2 层内自旋交换增强：Fe 3d 轨道占据优化

FeII 磁矩的显著增加源于自旋翻转——部分电子从自旋向下通道转移到自旋向上通道（Stoner 型机制）。

自旋翻转的驱动力：Te 从 Li 获得额外电子后，减少了对周围 Fe 的电子抽取——Fe-Te 键合区域的电子密度降低（region-B2）。

为维持电荷平衡，FeI-FeII 连接区域（region-Y1）获得电子——增加的电子占据发生在自旋向上通道，增强了自旋极化。

FeI-FeII 连接区域的额外电子 + FeII 增强的自旋极化——共同促进了 J₃ 的增强，优化了层内自旋交换通道。

5.3 电子重新分布的空间特征

Te 从 Li 获得额外电子后，由于 Te 价轨道的有限容量——某些外围区域（region-B1）的电子密度相应降低。

FeI-FeI 配对（J₄）：插层导致 FeI 亚层外侧电子密度增加——为减少 Coulomb 排斥，FeI-FeI 间隙区域（region-B3）电子密度降低。

FeI-Ga 连接区域（region-Y2）电子密度增加——提供额外的间接交换通道，增强 J₄。

这一分析视角——将交换参数增强与层内电子空间重新分布关联——为理解 vdW 磁体中磁耦合调控提供了新的视角。

5.4 DFT+U 方法的可靠性与局限

U = 1.3 eV 的选择使原始 FGaT 体相 T_C（365 K）与实验值（350-380 K）吻合——验证了方法的可靠性。

MC 模拟的 T_C 对交换参数 J_ij 敏感——J_ij 来自六种磁序的 DFT 总能差（Eqs. 2-7），精度受 U 值和 k 网格影响。

声子谱无虚频 + AIMD 500 K 稳定——Li-FGaT 是动力学和热力学稳定的，MC 模拟结果具有物理意义。

局限性：MC 模拟使用经典自旋近似，忽略了量子涨落——实际 T_C 可能略低于 770 K，但定性结论（Li 插层大幅提升 T_C）是稳健的。

5.5 研究局限性

理想化模型：本文聚焦于理想化的 Li-FGaT 体相模型——每个 vdW 间隙恰好一个 Li 离子，实际实验中 Li 浓度和分布可能存在不均匀性。

双层模型已部分弥补此局限——验证了不同 Li 浓度和空间分布下的磁性增强效果，但更复杂的浓度梯度效应尚未考虑。

经典自旋近似：MC 模拟基于经典 Heisenberg 模型——忽略了量子涨落、自旋波激发和磁振子-电子耦合等量子效应。

常温稳定性：虽然 AIMD 验证了 500 K 的结构稳定性，但更长时间尺度的 Li 离子扩散和可能的相分离未考虑。

5.6 未来研究方向

实验验证：通过电化学 Li 插层 + 高温磁化测量直接验证 770 K 的 T_C 预测——这是本工作提出的核心可测试目标。

Li 浓度系统优化：研究不同 Li 浓度（x）下 Li_xFe₃GaTe₂ 的 T_C——寻找最优掺杂浓度，平衡 T_C 提升和结构稳定性。

其他碱金属插层：Na、K、Rb 等碱金属离子——离子半径差异影响插层效率和晶格膨胀程度，可能产生不同的磁性增强效果。

多体量子效应：超越经典自旋近似，考虑磁振子-电子耦合和量子涨落——更精确预测 T_C 和自旋动力学行为。

【六、总结】

6.1 核心结论

DFT+MC 计算预测：Li 离子插层可将 Fe₃GaTe₂ 体相 T_C 从 365 K 提升至 770 K——超过两倍的增强，创造了 vdW 磁体 T_C 的预测新高。

T_C 提升的双重机制：(1) 层间——自旋极化界面电子充当 RKKY 型中介载流子，增强跨 vdW 间隙的 FM 耦合（J₅）；(2) 层内——电子掺杂诱导 Fe 3d 轨道占据优化，增强层内自旋交换（J₃、J₄）。

Li-FGaT 结构稳定性：声子谱无虚频（动力学稳定），AIMD 500 K 验证（热力学稳定），Li 插层能 −1.69 eV（能量稳定）。

双层模型验证：Li 插层对 FGaT 磁性的增强效果在厚度（体相→双层）、浓度（Li₁→Li₁/₄）和空间分布（分散→聚集）三个维度上具有普适性。

6.2 展望

实验验证：电化学 Li 插层 + 高温磁化测量——直接验证 T_C 预测，是本工作最直接的后续方向。

材料扩展：将 Li 插层策略应用于其他高 T_C vdW 磁体（如 Fe₅GeTe₂、Fe₄GeTe₂）——探索更高 T_C 的可能性。

器件集成：构建 Li-FGaT 基自旋阀和 MTJ——验证高温下的自旋输运性能，推动 vdW 磁体器件化。

理论深化：超越经典自旋近似，引入量子涨落和自旋波激发效应的 MC 模拟——更精确预测 T_C 和理解自旋动力学。

6.3 本文贡献

首次提出 Li 插层 Fe₃GaTe₂ 的晶格模型——通过系统的 DFT 计算和 MC 模拟，预测 T_C 可达 770 K。

揭示了电子掺杂增强磁性的双重微观机制——层间自旋极化界面电子中介 + 层内 Fe 3d 轨道占据优化，为磁性调控提供了新的理论视角。

通过双层模型验证了 Li 插层策略的普适性——在厚度、浓度和空间分布三个维度上均有效，增强了实验可行性。

为实验提供了明确的可测试目标——Li-FGaT 体系，T_C ~ 770 K，为高 T_C vdW 磁体的电子掺杂设计提供了理论框架。

【【DFT Tips】】

【DFT Tip】DFT+U 中 U 值的系统验证

Fe₃GaTe₂ 的 DFT+U 计算中 U=1.3 eV 是关键参数——该值使原始 FGaT 的 T_C（365 K）与实验（350-380 K）吻合。

常见错误：① 直接使用文献中其他 Fe 基材料的 U 值（如 Fe₃GeTe₂ 的 U=2.0 eV）→ T_C 严重偏离实验；② 未验证 U 值对磁矩和交换参数的敏感性。

建议：U 值在 1.0-2.0 eV 范围内测试，比较磁矩、交换参数 J_ij 和最终 T_C 的变化——U=1.3 eV 是 FGaT 的最优值。

【DFT Tip】能量映射方法的磁序选择

能量映射法通过构建不同磁序的 DFT 总能差提取交换参数 J_ij——磁序的选择必须覆盖所有独立交换通道。

常见错误：① 磁序数量不足 → 方程组欠定，J_ij 不唯一；② 磁序包含冗余组合 → 方程组病态，J_ij 误差大。

建议：n 个交换参数至少需要 n 个独立磁序；本文用 6 种磁序（MO-1 至 MO-6）提取 5 个 J_ij + 1 个参考能量，是最优设计。

【DFT Tip】vdW 磁性材料的 DFT-D3 修正

Fe₃GaTe₂ 的层间耦合由 vdW 力主导——PBE 泛函无法描述 vdW 相互作用，层间距和层间交换 J₅ 都会被严重低估。

常见错误：不加 DFT-D3 → 层间距过大 → 层间 FM 耦合 J₅ 被低估 → MC 模拟的 T_C 偏低。

建议：IVDW=11（DFT-D3 BJ damping），这是 vdW 磁体计算的标配。对于 Li 插层体系，vdW 修正对插层能的计算也至关重要。

【DFT Tip】Li 插层位点的能量筛选

Li 在 vdW 间隙中有多个候选位点（site-A、site-B、site-C）——必须通过 DFT 总能比较确定最优位点。

常见错误：仅凭几何直觉选择插层位点 → 可能选到亚稳态位点，插层能偏正 ~0.5 eV。

建议：在 2×2×1 超胞中测试所有高对称位点，计算插层能 E_int = E(Li-FGaT) − E(FGaT) − E(Li)——选择 E_int 最低的位点。

【DFT Tip】声子谱计算的超胞收敛

声子谱计算验证 Li-FGaT 的动力学稳定性——有限位移法需要足够大的超胞以避免周期性镜像的虚假相互作用。

常见错误：使用 1×1×1 单胞计算声子谱 → 某些声子模式在 Γ 点以外出现虚频，实际可能是超胞不够大导致的伪影。

建议：使用至少 2×2×1 超胞（~120 原子），验证虚频是否随超胞增大而消失——若虚频在 Γ 点，则可能是真实的不稳定性。

【DFT Tip】Monte Carlo 模拟的收敛性判断

MC 模拟的 T_C 对热平衡步数和采样步数敏感——不足的 MC 步数会导致 T_C 偏高或偏低。

常见错误：① 热平衡步数不足 → 初始磁构型的记忆效应，T_C 偏高；② 临近 T_C 时磁化涨落大 → 采样不足导致拟合误差。

建议：T_C 附近的温度点加倍采样步数（2×），验证 M-T 曲线在两次独立 MC 运行中的可重复性。

【DFT Tip】SOC 对磁晶各向异性的影响

Fe₃GaTe₂ 的垂直磁各向异性（PMA）完全由 SOC 决定——MCA 计算必须开启 SOC 并比较不同磁化方向的能量。

常见错误：不加 SOC 计算 MCA → 能量差为零，无法确定易轴方向；SOC 计算后未检查磁矩是否收敛到预期方向。

建议：LSORBIT=.TRUE.，SAXIS 分别沿 [001] 和 [100]，比较总能量——负值 E_MCA 表示 PMA，正值表示面内易磁化。

【DFT Tip】电子掺杂效应的 Bader 分析

Bader 电荷分析是量化 Li→FGaT 电子转移的关键工具——每个 Li 捐赠 ~0.84 e 给 FGaT。

常见错误：① 使用默认 FFT 网格 → Bader 电荷误差 ±0.1 e；② 将 Li 视为中性原子而非 Li⁺ → 电荷分析逻辑错误。

建议：使用 LAECHG=.TRUE. 生成 AECCAR0/1/2，NGXF×2 增加 FFT 精度；Bader 分析在自旋极化框架下进行。

【DFT Tip】AIMD 模拟的恒温器选择

AIMD 验证 Li-FGaT 在 500 K 的热力学稳定性——Nosé-Hoover 恒温器是 NVT 系综的标准选择。

常见错误：① 时间步长过大（>2 fs）→ 轻元素（Li）的动力学失真；② 模拟时间过短（<2 ps）→ 未达到热平衡。

建议：时间步长 1 fs，模拟至少 5 ps——前 2 ps 为平衡阶段，后 3 ps 为统计采样阶段。

【DFT Tip】层间磁耦合能 E_imc 的提取

E_imc = (E_AFM − E_FM)/u 量化层间 FM 耦合强度——Li 插层后 E_imc 从 3.3 meV 增至 11.2 meV。

常见错误：① 超胞中面内单胞数 u 计算错误 → E_imc 偏差因子为 2；② AFM 构型选择不当 → 层间 AFM 耦合被低估。

建议：AFM 构型应为相邻层反平行，FM 构型为所有层平行——两个构型在其他方面（层内磁序、原子位置）完全一致。

原始 FGaT 模型P6₃/mmc, 2 f.u./cellLi 插层位点筛选site-A/B/C 能量比较Li-FGaT 晶格弛豫c 轴膨胀 ~3%↓声子谱验证无虚频 → 动力学稳定AIMD 验证 (500 K)5 ps → 热力学稳定插层能E_int = −1.69 eV↓6 种磁序 DFTMO-1 至 MO-6能量映射J₁-J₅ 线性方程组MCA 计算 (SOC)SIA 参数 A↓Monte Carlo 模拟 → M-T 曲线 → Curie-Bloch 拟合 → T_C = 770 KFGaT: T_C=365 K (实验: 350-380 K) → Li-FGaT: T_C=770 K (2.1× 增强)↓层间增强：自旋极化界面电子Te 获取 ~0.39 e → 0.13 μB/Te层内增强：Fe 3d 轨道占据优化自旋翻转 → FeII 磁矩增加↓Li-FGaT：T_C = 770 K — vdW 磁体 T_C 新纪录预测实验验证：电化学 Li 插层 + 高温磁化测量 | 材料扩展：Na/K/Rb 插层结构构建 → 稳定性验证 → 交换参数提取 → MC 模拟 → 机制分析 → 实验预测

【【知识扩展】】

知识扩展①：能量映射方法与 Heisenberg 模型

【理论解释】能量映射法（Energy Mapping）将 DFT 总能映射到经典 Heisenberg 自旋哈密顿量 H = Σ J_ij S_i·S_j——通过构建不同磁序并计算 DFT 总能差，建立线性方程组求解交换参数 J_ij。这是连接第一性原理和自旋模型的标准方法。

【方法比较】交换参数提取方法：① 能量映射（最常用，适合共线磁序）；② 磁力定理（Lichtenstein 公式，适合非共线/海森堡体系）；③ 自旋波色散拟合（实验或 DFT 计算 magnon 谱）。能量映射法简单、稳健，但要求磁序为共线或接近共线。

【经典参考】经典论文：Xiang et al., Dalton Trans. 42, 823 (2013) — 能量映射法在 2D 磁体中的系统应用；Lichtenstein et al., JMMM 67, 65 (1987) — 磁力定理的原始推导。MC 工具：MCsolver (Liu et al., JPCL 2020), VAMPIRE, UppASD。

【迁移能力】能量映射法适用于任何可构建多种磁序的体系：vdW 磁体（CrI₃、Fe₃GeTe₂）、氧化物磁体、金属间化合物。关键限制：磁矩必须局域化，巡游磁体（如 Fe、Co）不适合。

知识扩展②：RKKY 机制与层间磁耦合

【理论解释】Ruderman-Kittel-Kasuya-Yosida (RKKY) 机制描述局域磁矩通过传导电子的自旋极化间接交换——交换强度 J(R) ∝ cos(2k_F R)/R³ 振荡衰减。在 vdW 磁体中，层间 FM/AFM 耦合常由跨 vdW 间隙的 RKKY 型间接交换主导。

【方法比较】层间耦合的调控策略：① 电子掺杂（Li/Na 插层）→ 增加传导电子浓度 → 增强 RKKY 耦合；② 层间距调控（压力/应变）→ 改变 RKKY 振荡相位 → 可能翻转 FM/AFM；③ 插层磁性离子 → 直接交换替代间接交换。

【经典参考】经典论文：Ruderman & Kittel, PR 96, 99 (1954) — RKKY 原始理论；Kasuya, PTP 16, 45 (1956)；Yosida, PR 106, 893 (1957)。vdW 磁体中的 RKKY：Wang et al., PRB 104, 064430 (2021) — CrI₃ 层间耦合的 RKKY 解释。

【迁移能力】RKKY 分析适用于所有金属性层状磁体：Fe₃GeTe₂、Fe₃GaTe₂、Fe₅GeTe₂ 等。关键：需有可极化的传导电子（金属性），且费米波矢 k_F 与层间距匹配。半导体/绝缘体（如 CrI₃）的层间耦合更适用超交换机制。

【【科研经验】】

科研经验①：为什么 MC 模拟的 T_C 与实验不符？

【问题】MC 模拟的 T_C 系统性地高于实验值（如 FGaT 预测 365 K vs 实验 350-380 K 吻合，但其他体系可能偏差 50-100 K）。

【原因】① 经典自旋近似忽略了量子涨落——量子涨落在低维体系中降低 T_C 约 10-20%；② 交换参数 J_ij 来自 DFT 总能差，DFT 泛函误差（~5-10 meV）可导致 T_C 偏差 ~50 K；③ 有限尺寸效应——MC 超胞不够大，T_C 偏高。

【解决方案】① 使用量子 MC 或自旋波修正因子（~0.7-0.8）校正经典 T_C；② 比较不同泛函（PBE、PBE+U、SCAN）的 J_ij 差异；③ 逐步增大 MC 超胞（20×20×10 → 60×60×20），确认 T_C 收敛。

【建议】DFT+MC 的 T_C 预测应视为半定量——定性趋势（Li 插层增强 T_C）是稳健的，但绝对 T_C 值（770 K）有 ±50-100 K 的不确定性。

科研经验②：为什么插层能 E_int 为负但实验难以实现？

【问题】DFT 计算显示 Li 插层能 E_int = −1.69 eV（强放热），但实验上电化学 Li 插层可能不完全或需要特定条件。

【原因】① DFT 计算的是 0 K 静态总能，未考虑插层动力学（Li 扩散势垒）；② 电化学插层涉及溶剂化效应和界面反应，DFT 未包含；③ 高 Li 浓度下 Li-Li 排斥可能导致相分离。

【解决方案】① 使用 NEB（Nudged Elastic Band）计算 Li 扩散势垒——评估插层动力学可行性；② 构建不同 Li 浓度的凸包图（Convex Hull）——判断热力学稳定浓度范围；③ 考虑溶剂化效应（隐式溶剂模型 VASPsol）。

【建议】DFT 插层能是必要非充分条件——负值表明热力学允许，但实验实现需结合电化学知识。双层低浓度模型（Li₁/₄）更接近实际实验条件。

进一步计算①：Li 扩散势垒的 NEB 计算

【为什么值得算】插层能 E_int 只说明热力学可行性，NEB 计算 Li 在 vdW 间隙中的扩散势垒——评估插层动力学可行性。

【能回答的问题】Li 在 FGaT 层间的扩散势垒是多少？是否存在快速扩散通道？Li 倾向于聚集还是分散？

【适合体系】所有插层/离子电池材料。

【输入】VASP NEB 计算（初始态和末态结构），3-5 个中间镜像。

进一步计算②：不同碱金属插层的比较（Na, K, Rb）

【为什么值得算】Li 是最小碱金属，Na/K/Rb 离子半径更大——插层效率和磁性增强效果可能不同，甚至有质的差异。

【能回答的问题】Na/K/Rb 插层的 T_C 是否高于 Li？离子半径与 T_C 增强的关系？哪种碱金属是最优选择？

【适合体系】所有 vdW 层状材料。

【输入】Na/K/Rb 插层 FGaT 的 DFT 计算 + MC 模拟（与 Li 类似流程）。

【

进一步计算③：自旋波谱与磁振子能隙

【为什么值得算】自旋波谱（magnon dispersion）揭示低能磁激发——磁振子能隙决定自旋波对 T_C 的量子修正，也是非弹性中子散射的直接可观测量。

【能回答的问题】FGaT 和 Li-FGaT 的磁振子能隙分别是多少？自旋波刚度如何变化？量子涨落对 T_C 的修正比例？

【适合体系】所有有序磁体。

【输入】DFT 交换参数 J_ij + 自旋波计算（SpinW 或 McPhase）。

进一步计算④：Li 浓度依赖的 T_C 相图

【为什么值得算】本文仅计算了 Li₁ 和 Li₁/₄ 两个浓度——系统构建 Li_xFGaT 的 T_C(x) 相图，寻找最优掺杂浓度。

【能回答的问题】T_C 随 Li 浓度如何变化？是否存在最优浓度（T_C 峰值）？过掺杂是否导致 T_C 下降？

【适合体系】所有掺杂/插层调控的磁性体系。

【输入】不同 Li 浓度（x=0, 0.25, 0.5, 0.75, 1.0）的 DFT + MC 模拟。

进一步计算⑤：电子结构中的拓扑特征

【为什么值得算】Fe₃GaTe₂ 是金属性铁磁体，其能带可能具有非平庸拓扑（如外尔点、节线）——拓扑与磁性的共存是当前热点。

【能回答的问题】FGaT 和 Li-FGaT 中是否存在外尔点？Li 插层是否改变拓扑特征？反常霍尔电导（AHC）的大小？

【适合体系】所有金属性铁磁体。

【输入】DFT 能带结构 + Wannier90 拟合 + Berry 曲率计算。

进一步计算⑥：单层/少层 Li-FGaT 的磁性

【为什么值得算】器件应用通常基于少层而非体相材料——验证 Li 插层策略在单层和双层 FGaT 中的有效性。

【能回答的问题】单层 Li-FGaT 的 T_C 是多少？少层体系的 T_C 是否仍远超室温？维度降低对 Li 插层效应的影响？

【适合体系】所有 vdW 磁体的少层体系。

【输入】单层/双层 Li-FGaT 的 DFT 计算（板模型 + 真空层）+ MC 模拟。

【支撑信息 (Supporting Information)】

【S1 · 额外磁序的示意图】

S1.1 MO-7、MO-8 和 MO-9 磁序

展示三种额外磁序（MO-7、MO-8 和 MO-9）的示意图——用于验证仅用 J₅ 描述层间 FM 耦合的充分性。

正文中仅使用了六种磁序（MO-1 至 MO-6）提取交换参数——引入额外磁序是为了检验额外层间交换参数（J₆、J₇）的必要性。

磁序选择的理由：MO-1 至 MO-6 覆盖了五种交换参数的所有独立线性组合——确保方程组（Eqs. 2-7）可解。

额外磁序 MO-7/8/9 的层间自旋排列与 MO-1 至 MO-6 不同——用于检验 J₅ 是否足以描述所有层间耦合模式。

图 S1：三种额外磁序（MO-7、MO-8 和 MO-9）的示意图——用于验证交换参数提取方案的充分性。

【S2 · 两组交换参数的 M-T 曲线比较】

S2.1 交换参数提取方案的验证

比较使用两组交换参数获得的 MC 模拟 M-T 曲线——第一组仅包含 J₁-J₅，第二组额外包含 J₆ 和 J₇。

两条 M-T 曲线几乎完全重合——T_C 差异在统计误差范围内，说明额外层间交换参数对 T_C 的影响可忽略。

该结果确认：仅用 J₅ 描述层间 FM 耦合是充分的——简化了交换参数提取方案，只需六种磁序即可准确描述。

这一验证为后续 Li-FGaT 的交换参数计算提供了方法学基础——无需引入额外磁序，节省计算资源。

图 S2：使用两组交换参数获得的 MC 模拟 M-T 曲线比较——验证仅用 J₅ 描述层间 FM 耦合的充分性。

【S3 · Li₁/₄-(FGaT)₂ 双层模型结构】

S3.1 低 Li 浓度双层模型

展示 Li₁/₄-(FGaT)₂ 模型的侧视图和俯视图——为清晰显示 Li 离子，上层 FGaT 单层被部分移除。

该模型每四个 vdW 间隙中插入一个 Li 离子——代表较低的 Li 插层浓度（x = 0.25），更接近实验上可能实现的掺杂水平。

FGaT 双层模型包含真空层——模拟 vdW 异质结中少层 FGaT 的实际环境，上下表面无悬挂键。

与 Li₁-(FGaT)₂（每间隙一个 Li，x = 1.0）和原始 FGaT 双层（x = 0）对比——正文 Fig. 5 展示了三种模型的 M-T 曲线。

图 S3：Li₁/₄-(FGaT)₂ 模型的侧视图和俯视图——为清晰显示 Li 离子，上层 FGaT 单层被部分移除。

【S4 · ag-Li₁/₄-(FGaT)₂ 聚集模型结构】

S4.1 Li 离子聚集排列的双层模型

展示 ag-Li₁/₄-(FGaT)₂ 模型的侧视图和俯视图——Li 离子以聚集（aggregated）方式排列，附有 Li 离子周围局部配位环境的放大视图。

与 S3 中 Li 离子均匀分散排列模型对比——评估 Li 离子排列方式对磁性的影响。

能量上分散排列更优（低 3.8 meV/Li）——但聚集排列可能以一定概率出现，尤其在非平衡插层条件下。

层间磁耦合能 E_imc：聚集 9.0 meV > 原始 3.3 meV——即使聚集构型，Li 插层仍显著增强 FM 耦合。

图 S4：ag-Li₁/₄-(FGaT)₂ 模型的侧视图和俯视图——附有 Li 离子周围局部配位环境的放大视图。

层间磁耦合能定义：E_imc = (E_AFM − E_FM)/u——u 为超胞中面向单胞数，正值越大表示层间 FM 耦合越强。

【主要参考文献】

核心引用

[1] Song Z, Mo Y, Zeng R, Liu L, Jiang X, Huang X, Liu J-M. J. Appl. Phys. 140, 014302 (2026) — 本工作。

[2] Deng Y, et al. Nature 563, 94 (2018) — Li 插层 Fe₃GeTe₂ 提升 T_C 至 >300 K。

[3] Wang Z, et al. ACS Nano 18, 23058 (2024) — Li 插层 CrI₃ 提升 T_C 至 >400 K。

[4] Zhang G, et al. Nat. Commun. 13, 5067 (2022) — Fe₃GaTe₂ 的室温铁磁性实验发现。

[5] Wang M, et al. npj 2D Mater. Appl. 8, 22 (2024) — FGaT 单层/双层 T_C 实验测量。

[6] Huang X, et al. Appl. Phys. Lett. 119, 012405 (2021) — DFT 揭示 Li 插层 Fe₃GeTe₂ 的磁性增强机制。

[7] Huang X, et al. Phys. Rev. B 113, 184409 (2026) — 能量映射方法在 vdW 磁体中的应用。

[8] Liu L, et al. J. Phys. Chem. Lett. 11, 7893 (2020) — MCsolver 蒙特卡洛模拟软件。

Z. Song, Y. Mo, R. Zeng, L. Liu, X. Jiang, X. Huang, J.-M. Liu | J. Appl. Phys. 140, 014302 (2026) | Fe₃GaTe₂ · Li 插层 · 居里温度 · vdW 磁性