过渡金属插层范德华材料展现出丰富的磁学性质，强烈依赖于插层离子的超晶格有序性。本文通过 DFT 构建成对晶格模型，结合蒙特卡洛模拟确定热力学相图。使用抗磁离子 Zn2+ 和 Sc3+ 作为磁性离子的代理，揭示了 √3×√3 有序在远宽于传统预期的范围内占主导——即使在 x=1/4 处，√3×√3 畴也优于 2×2。Raman 光谱实验证实了 2×2 有序的脆弱性。
来源：J. Am. Chem. Soc. 147, 13629-13641 (2025) · UC Berkeley / LBNL
一、研究背景：过渡金属插层 TMD 的超晶格相图
1.1 过渡金属插层范德华材料
范德华（vdW）材料插层过渡金属（TM）离子后展现出丰富的磁学性质，强烈依赖于插层离子的结构有序性。2H-TaS2 是最具代表性的插层主体材料之一，当 1/4 和 1/3 的八面体间隙位被占据时（对应体相组成 2H-M0.25TaS2 和 2H-M0.33TaS2），插层离子可形成 2×2 和 √3×√3 超晶格，展现出截然不同的磁学性质。
插层 TMD 体系的关键要素
主体晶格：2H-TaS2 提供八面体间隙位，形成三角晶格
插层离子：Fe2+、Co2+、V3+、Cr3+ 等磁性离子
超晶格类型：2×2（x=0.25）和 √3×√3（x=0.33），决定磁交换机制
磁交换：RKKY 型振荡交换——微妙的结构变化导致巨大的磁性差异
1.2 为什么需要理论模型？
尽管理论上广泛预期特定化学计量比（如 x=1/4）对应特定超晶格（如 2×2），但实验上 x≈0.25 的样品常偏离 2×2 对称性，出现有序-无序转变和 √3×√3 共存畴。STEM 显示 2H-Cr0.25TaS2 中存在 2×2 有序畴（数百纳米）与少量 √3×√3 区域共存，而单晶 XRD 仅显示 2×2 有序。这表明需要有限温度下的热力学相图理论来指导实验。
本工作的策略：用抗磁离子代理磁性离子
磁性转变温度（~10-150 K）远低于合成温度（~1000-1400 K）→ 磁性耦合在决定高温位点占据能方面作用不大
因此可以用抗磁离子 Zn2+ 和 Sc3+ 作为 Fe2+/Co2+ 和 Cr3+/V3+ 的代理
Zn/Sc 与 Fe/Cr 具有相似的离子半径、功函数和氧化态 → 结构行为相似，但无磁性复杂性
二、计算方法：DFT + 晶格模型 + 蒙特卡洛
2.1 成对晶格哈密顿量
定义占据变量 σi ∈ {0, 1}（0=空位，1=占据），插层晶格能量为：
晶格哈密顿量
E(σ) = E0 + Σij J(|ri−rj|) σiσj + μ Σi σi
J(r)：距离依赖的成对耦合常数，由 DFT 拟合得到
μ：参考化学势（相对于空晶格）
这对应于一个二维三角晶格气模型，相互作用超出最近邻
2.2 DFT 计算参数
计算设置
GGA+U 水平（Dudarev 型 U，线性响应获得）
刚性晶格近似（弹性贡献 ~0.4 meV，远小于电子效应 ~200-300 meV）
Kawasaki 动力学蒙特卡洛模拟（固定插层密度）
序参量：RMS(γi) 表征 √3×√3 长程序，RMS(φi) 表征 2×2 长程序，⟨α̅⟩ 表征短程序
三、插层驱动的结构效应
3.1 结构弛豫
DFT 结构优化显示：Zn 完全插层导致面内 a0 膨胀 2%、面外 c 膨胀 3%；Sc 插层导致 a0 几乎不变、c 大幅膨胀 23%。Zn 和 Sc 的面内膨胀均很小，意味着面内弹性效应可忽略。弹性哈密顿量能量尺度 ~0.4 meV，而电子效应 ~200-300 meV，相差几个数量级。
三角畸变与晶体场
TaS2 八面体位存在三角畸变（键角 > 90°）
影响 d 轨道晶体场分裂：dz² 与 (dxz, dyz) / (dx²−y², dxy) 分离
Zn2+：d10 闭壳层 → 抗磁
Sc3+：d0 闭壳层 → 抗磁
图 1 | DFT 结构参数：(a) 2H-TaS2、2H-Zn0.5TaS2 和 2H-Sc0.5TaS2 双层结构参数；(b) 三角畸变对晶体场分裂的影响；(c) Zn2+ 和 Sc3+ 的预期 d 轨道电子排布
四、插层驱动的电子效应
4.1 电子态密度（DOS）
Zn 插层：刚性带图像在低插层密度下基本成立——Zn 态与主体态几乎无杂化，与 FexTaS2 中观测到的更强杂化形成对比。Fermi 能级稳步上升反映约 1 个电子/Zn 离子的电荷转移。Sc 插层：刚性带假设明显失效，Sc 与 S 之间更强的杂化导致较大偏离——类似于 Cr 插层 TaS2 中观测到的强 Cr-Ta 杂化。
DOS 关键发现
Zn 插层：Zn 态孤立，几乎无杂化 → 刚性带图像近似成立
Sc 插层：Sc 轨道更大，与 S 能级对齐更好 → 杂化增强 → 刚性带假设失效
轨道特征：Fermi 能级附近主要为 Ta dz²，更高能态为 (dx²−y², dxy) 和 (dxz, dyz)
Bader 电荷：Zn 的 Bader 电荷 ~0.83-0.90（接近 1+），Sc 的 Bader 电荷 ~1.47-1.86
图 2 | 电子态密度：(a) 2H-ZnxTaS2 在不同插层密度下的 DOS；(b) 2H-ScxTaS2 在不同插层密度下的 DOS。Fermi 能级以灰线标注
五、电荷密度重排与晶格哈密顿量拟合
5.1 电荷密度差的实空间分布
电荷转移后，掺杂主体的电子仍然局域在插层离子周围的弥散壳层中，而非分散到整个晶格。Zn 离子周围：vdW 界面电荷大幅增加，TMD 晶格电荷分布变化极小（~0.0015 au/Bohr³）。Sc 离子周围：行为类似，但电荷密度变化略微更分散。
关键含义
J(r > 2a0) 的电荷密度变化可忽略 → 仅需考虑短程耦合
J(r = √3a0) 主要由周围 TaS2 晶格电荷变化驱动
耦合常数 J(r) 不代表静电相互作用，而是局域电荷密度变化的粗粒化描述
图 3 | 电荷密度差：(a) 稀 Sc 插层（x=1/18）和 Zn 插层后的电荷密度差，红色/橙色/黄色/青色等值面分别对应 0.025/0.01/0.0015/−0.0015 au/Bohr³；(b) 代表性非稀占据构型的电荷密度差
5.2 静电屏蔽的定量证据
计算以 Zn 和 Sc 位点为中心、半径递增的球内总电荷：Zn 离子在半径 1.58 Å 处完全电荷补偿，Sc 离子在 1.88 Å 处完全补偿。这意味着 Zn 离子间静电耦合在 J(r ≥ a0) 完全屏蔽，Sc 离子间在 J(r ≥ √3a0) 完全屏蔽。
图 4 | 球内总电荷：(a) Zn 和 (b) Sc 插层位点周围递增半径球内总电荷，标注 PAW 球内电荷和 Bader 电荷；(c-d) 不同插层构型下的对比
5.3 耦合常数拟合结果
拟合结果（表 1）
Zn 插层：RMSE = 9 meV，拟合质量良好 → 成对模型合理
Sc 插层：RMSE = 52 meV，拟合质量较差 → 需更高阶团簇模型
J(a0) 比 J(r > a0) 大至少一个数量级 → 最近邻排斥主导
Sc 的 J(a0) 显著大于 Zn —— 因更强的静电排斥
六、超晶格相图
6.1 蒙特卡洛模拟结果
基于拟合的晶格气参数，蒙特卡洛模拟得到的温度-插层密度相图（图 5）揭示：
相图核心发现
√3×√3 有序远优于 2×2：√3×√3 在更宽的插层密度和温度范围内持续存在
2×2 有序极为脆弱：仅局限于相图的狭窄区域
在 x=1/4 处：传统上预期 2×2 有序，但实际 √3×√3 畴占主导
短程序：即使无序相中，⟨α̅⟩ 仍显著低于随机分布预期值——最近邻排斥持续存在
6.2 物理起源：构型熵
由于 J(a0) 远大于 J(r > a0)，最近邻排斥主导了短程序。2×2 和 √3×√3 有序的偏好主要由构型熵驱动——而非能量差异。√3×√3 有序在更宽范围存在是因为它提供了更多的构型自由度。线扫描（图 6）显示：即使在无序相中，⟨α̅⟩ 仍远低于零（随机分布的预期值），表明最近邻空位偏好在所有合成温度下均存在。
图 5 | 温度-插层密度相图：(a) 整体相图；(b) 2×2 长程序 RMS(φi)；(c) √3×√3 长程序 RMS(γi)；(d) 短程序 ⟨α̅⟩
图 6 | ⟨α̅⟩ vs x 的线扫描：Zn 插层 2H-TaS2 双层在不同温度下的序参量，对比随机分布的 ⟨α̅⟩=0 预期
七、实验验证：Raman 光谱
7.1 V0.25NbS2 的淬火实验
为验证理论预测的 2×2 有序的脆弱性，对 V0.25NbS2 晶体进行 Raman 光谱测量。~200 cm−1 处的尖锐峰是 2×2 有序平面内激发的特征峰。淬火后该峰消失——表明 2×2 有序被破坏，√3×√3 有序在高温下占主导。这一实验证据支持了相图的关键结论：在 x=0.25 处，√3×√3 畴在高温下优于 2×2。
Raman 证据
~200 cm−1 峰：2×2 有序的特征标志
淬火后峰消失 → 2×2 → √3×√3 转变
支持高温下 √3×√3 优于 2×2 的理论预测
图 7 | V0.25NbS2 晶体在 900°C 淬火前后的 Raman 光谱。~200 cm−1 的尖锐峰是 2×2 有序平面内激发的特征峰
八、讨论与展望
本工作建立了 TM 插层 TMD 超晶格有序的热力学相图理论框架，揭示了 √3×√3 有序在远宽于传统预期的范围内占主导，而 2×2 有序仅存在于狭窄相区。
关键科学意义
① 挑战传统假设：x=1/4 不必然意味着 2×2 有序——√3×√3 共存或主导是热力学结果
② 屏蔽机制：金属主体的高效介电屏蔽使静电耦合几乎完全消失——耦合来自电荷密度重排而非长程静电
③ 构型熵驱动：短程排斥 + 构型熵 → √3×√3 有序的广泛存在
④ 实验启示：为靶向特定超晶格相提供材料选择和热处理的指导
局限性与未来方向
① 成对模型对高价态离子（Sc3+）精度有限 → 需团簇展开
② 刚性晶格近似 → 需考虑面外弹性效应
③ 仅建模了抗磁代理离子 → 需纳入磁性自由度
④ 当前模型可推广至其他插层 TMD 体系（如 NbS2、TiSe2）
九、支持信息 (Supporting Information)
S1. PAW 球投影 DOS：Zn 插层
PAW 球面谐波投影的态密度，展示 Zn 插层 Ta18S36Zn9x 在不同插层密度和不同离子排序下的 DOS。
图 S1 | PAW 球投影 DOS：Ta18S36Zn9x 在不同插层密度 x 和不同离子排序下的电子态密度
S2. PAW 球投影 DOS：Sc 插层
PAW 球面谐波投影的态密度，展示 Sc 插层 Ta18S36Sc9x 在不同插层密度和不同离子排序下的 DOS。
图 S2 | PAW 球投影 DOS：Ta18S36Sc9x 在不同插层密度 x 和不同离子排序下的电子态密度
S3. Bader 电荷球投影 DOS
使用 Bader 电荷分析半径的球面谐波投影 DOS，与 PAW 球投影结果定性一致。
图 S3 | Bader 电荷球投影 DOS：Ta18S36Zn9x 的态密度，Fermi 能级以灰线标注
S4. DFT 能量 vs a0
ZnTa2S4 的 DFT 能量作为 a0 的函数及二次拟合，用于确定弹性弹簧常数 k。
图 S4 | DFT 能量 vs a0：ZnTa2S4 的能量随晶格常数的变化及二次拟合
S5. 蒙特卡洛轨迹（一）
选定温度和位点占据率的蒙特卡洛轨迹及运行平均值。
图 S5 | 蒙特卡洛轨迹：选定温度 T 和位点占据率 f 的 MC 运行平均值
S6. 蒙特卡洛轨迹（二）
更多温度-占据率组合的蒙特卡洛轨迹及运行平均值。
图 S6 | 蒙特卡洛轨迹：选定温度 T 和位点占据率 f 的 MC 运行平均值
S7. 蒙特卡洛轨迹（三）
更多温度-占据率组合的蒙特卡洛轨迹及运行平均值。
图 S7 | 蒙特卡洛轨迹：选定温度 T 和位点占据率 f 的 MC 运行平均值
S8. U 值选择对电荷密度的影响
不同 U 值（0, 1, 2, 2.9）对 2H-Sc1/18TaS2 电荷密度的影响——Sc 体系对 U 值的依赖比 Zn 更强。
图 S8 | U 值对电荷密度的影响：2H-Sc1/18TaS2 在不同 U 值下的 Bader 电荷
S9. Vq 和 VR
方程 3 给出的 Vq（左）和方程 2 给出的 VR（右），通过逆 Fourier 变换评估。
图 S9 | Vq（左）和 VR（右）：通过逆 Fourier 变换评估的相互作用势
S10. V0.25NbS2 SEM 图像
V0.25NbS2 晶体的扫描电子显微镜图像，展示晶体形貌。
图 S10 | V0.25NbS2 晶体的 SEM 图像
S11. V0.25NbS2 Raman 光谱（多晶体）
四个不同 V0.25NbS2 晶体在淬火前后的 Raman 光谱，每个晶体取三个采集点，展示了一定的不均匀性。
图 S11 | V0.25NbS2 晶体的 Raman 光谱：淬火前后对比，四个不同晶体，每个三个采集点
参考文献
[1] Craig IM, Kim BJ, Fender SS, et al. J. Am. Chem. Soc.147, 13629-13641 (2025).[2] Morosan E, et al.Phys. Rev. B76, 014408 (2007). [FexTaS2磁性][3] Ko K-T, et al.Nat. Commun.2, 567 (2011). [Fe0.25TaS2磁各向异性][4] Checkelsky JG, et al.Phys. Rev. B86, 180402(R) (2012). [FxNbS2磁性][5] Nair NL, et al.Nat. Mater.19, 153 (2020). [Cr1/3NbS2手性磁体][6] Hardy WJ, et al.Phys. Rev. B91, 054426 (2015). [Fe0.25TaS2磁畴][7] Park S, et al.Nat. Commun.15, 8797 (2024). [Cr1/3TaS2晶体场][8] Fan S, et al.Phys. Rev. Lett.119, 126402 (2017). [插层 TMD 角分辨光电子能谱]
