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

Li 插层 Fe3GaTe2中的显著磁性增强J. APPL. PHYS. 2026Li 插层 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. 2021T_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 插层 FGaTLi-FGaT晶格模型——预测 T_C 可达 770 K并揭示电子掺杂增强磁性的微观机制。本文还考察了双层模型不同 Li 浓度和空间分布的磁性——验证 Li 插层策略在厚度、浓度和空间分布上的普适性。最终目标为实验提供可测试的明确预测Li-FGaT, T_C ~ 770 K并为高 T_C vdW 磁体的电子掺杂设计提供理论框架。【二、研究方法】2.1 DFT 计算框架采用 VASPVienna Ab initio Simulation Package程序包进行第一性原理计算——基于平面波基组和赝势方法。交换关联泛函PBEPerdew-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——采用 DFTU 旋转不变方案。U 值的选择依据该 U 值使计算得到的 FGaT 体相 T_C365 K与实验测量值350-380 K吻合良好——验证了参数设置的可靠性。电子自洽收敛标准10⁻⁶ eV——确保总能计算精度足以分辨不同磁序之间的微小能量差异~meV 量级。DFTU 方法在过渡金属氧化物和硫属化合物中广泛应用——能够有效修正 3d 电子的自相互作用误差。2.3 k 点网格与计算精度Monkhorst-PackΓ 中心k 点网格体相单胞 24×24×42×2×1 超胞 12×12×4——密集网格确保总能收敛。单层/双层板模型24×24×1单胞或 12×12×12×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 离子在八面体位点能量稳定。【三、实验结果】图 1FGaT 体相。(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 等基本一致——进一步确认了 DFTUU 1.3 eV方案的准确性。Heisenberg 自旋哈密顿量J_ij 为交换耦合常数S_i 为经典自旋单位矢量A 为单离子各向异性参数——求和遍历所有 Fe 位点对。能量映射方程组示例E_bulk(1) E₀ 48J₁ 24J₂ 48J₃ 8J₄ 8J₅——六种磁序对应六个方程联立求解 J₁-J₅。图 2Li-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 晶格。图 3Li-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——比原始 FGaT13.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是原始 FGaT365 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.32FGaT或 β0.35Li-FGaT。MCA 能量定义E_MCA E[001] − E[100] N_Fe·A——负值表示 PMAN_Fe 为单胞中 Fe 位点数。图 4Li-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/TeFe 的电子占据数变化很小。自旋极化界面电子Te 获取的掺杂电子受金属性电子结构自旋极化影响——每个 Te 贡献 0.13 μB为层间自旋交换提供中介载流子。FeII 磁矩增加尽管 Fe 的电子占据数变化很小FeII 的磁矩显著增加——部分电子从自旋向下通道翻转到自旋向上通道。层内电子重新分布Te 获得额外电子后减少对周围 Fe 的电子抽取——Fe-Fe 直接交换和 Fe-Te-Fe 超交换通道均被增强。图 5MC 模拟的 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)₂每间隙一个 LiT_C 490 K、Li₁/₄-(FGaT)₂每四个间隙一个 LiT_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_imcag-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 DFTU 方法的可靠性与局限U 1.3 eV 的选择使原始 FGaT 体相 T_C365 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 核心结论DFTMC 计算预测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】DFTU 中 U 值的系统验证Fe₃GaTe₂ 的 DFTU 计算中 U1.3 eV 是关键参数——该值使原始 FGaT 的 T_C365 K与实验350-380 K吻合。常见错误① 直接使用文献中其他 Fe 基材料的 U 值如 Fe₃GeTe₂ 的 U2.0 eV→ T_C 严重偏离实验② 未验证 U 值对磁矩和交换参数的敏感性。建议U 值在 1.0-2.0 eV 范围内测试比较磁矩、交换参数 J_ij 和最终 T_C 的变化——U1.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 偏低。建议IVDW11DFT-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/2NGXF×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_C365 K (实验: 350-380 K) → Li-FGaT: T_C770 K (2.1× 增强)↓层间增强自旋极化界面电子Te 获取 ~0.39 e → 0.13 μB/Te层内增强Fe 3d 轨道占据优化自旋翻转 → FeII 磁矩增加↓Li-FGaTT_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 磁体中的 RKKYWang 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、PBEU、SCAN的 J_ij 差异③ 逐步增大 MC 超胞20×20×10 → 60×60×20确认 T_C 收敛。【建议】DFTMC 的 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 排斥可能导致相分离。【解决方案】① 使用 NEBNudged 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 浓度x0, 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)₂每间隙一个 Lix 1.0和原始 FGaT 双层x 0对比——正文 Fig. 5 展示了三种模型的 M-T 曲线。图 S3Li₁/₄-(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 耦合。图 S4ag-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 磁性