过渡金属插层范德华材料展现出丰富的磁学性质强烈依赖于插层离子的超晶格有序性。本文通过 DFT 构建成对晶格模型结合蒙特卡洛模拟确定热力学相图。使用抗磁离子 Zn2 和 Sc3 作为磁性离子的代理揭示了 √3×√3 有序在远宽于传统预期的范围内占主导——即使在 x1/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×2x0.25和 √3×√3x0.33决定磁交换机制磁交换RKKY 型振荡交换——微妙的结构变化导致巨大的磁性差异1.2 为什么需要理论模型尽管理论上广泛预期特定化学计量比如 x1/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 σiJ(r)距离依赖的成对耦合常数由 DFT 拟合得到μ参考化学势相对于空晶格这对应于一个二维三角晶格气模型相互作用超出最近邻2.2 DFT 计算参数计算设置GGAU 水平Dudarev 型 U线性响应获得刚性晶格近似弹性贡献 ~0.4 meV远小于电子效应 ~200-300 meVKawasaki 动力学蒙特卡洛模拟固定插层密度序参量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) 分离Zn2d10 闭壳层 → 抗磁Sc3d0 闭壳层 → 抗磁图 1 | DFT 结构参数(a) 2H-TaS2、2H-Zn0.5TaS2 和 2H-Sc0.5TaS2 双层结构参数(b) 三角畸变对晶体场分裂的影响(c) Zn2 和 Sc3 的预期 d 轨道电子排布四、插层驱动的电子效应4.1 电子态密度DOSZn 插层刚性带图像在低插层密度下基本成立——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接近 1Sc 的 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 插层x1/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 耦合常数拟合结果拟合结果表 1Zn 插层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 有序极为脆弱仅局限于相图的狭窄区域在 x1/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 有序在高温下占主导。这一实验证据支持了相图的关键结论在 x0.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 有序仅存在于狭窄相区。关键科学意义① 挑战传统假设x1/4 不必然意味着 2×2 有序——√3×√3 共存或主导是热力学结果② 屏蔽机制金属主体的高效介电屏蔽使静电耦合几乎完全消失——耦合来自电荷密度重排而非长程静电③ 构型熵驱动短程排斥 构型熵 → √3×√3 有序的广泛存在④ 实验启示为靶向特定超晶格相提供材料选择和热处理的指导局限性与未来方向① 成对模型对高价态离子Sc3精度有限 → 需团簇展开② 刚性晶格近似 → 需考虑面外弹性效应③ 仅建模了抗磁代理离子 → 需纳入磁性自由度④ 当前模型可推广至其他插层 TMD 体系如 NbS2、TiSe2九、支持信息 (Supporting Information)S1. PAW 球投影 DOSZn 插层PAW 球面谐波投影的态密度展示 Zn 插层 Ta18S36Zn9x 在不同插层密度和不同离子排序下的 DOS。图 S1 | PAW 球投影 DOSTa18S36Zn9x 在不同插层密度 x 和不同离子排序下的电子态密度S2. PAW 球投影 DOSSc 插层PAW 球面谐波投影的态密度展示 Sc 插层 Ta18S36Sc9x 在不同插层密度和不同离子排序下的 DOS。图 S2 | PAW 球投影 DOSTa18S36Sc9x 在不同插层密度 x 和不同离子排序下的电子态密度S3. Bader 电荷球投影 DOS使用 Bader 电荷分析半径的球面谐波投影 DOS与 PAW 球投影结果定性一致。图 S3 | Bader 电荷球投影 DOSTa18S36Zn9x 的态密度Fermi 能级以灰线标注S4. DFT 能量 vs a0ZnTa2S4 的 DFT 能量作为 a0 的函数及二次拟合用于确定弹性弹簧常数 k。图 S4 | DFT 能量 vs a0ZnTa2S4 的能量随晶格常数的变化及二次拟合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 角分辨光电子能谱]
