SCIENTIFIC REPORTS 非 van der Waals可调二维磁性材料设计平台

SCIENTIFIC REPORTS 非 van der Waals可调二维磁性材料设计平台

SCIENTIFIC REPORTS 16, 19079 (2026)

Cr₂Se₃：非 van der Waals可调二维磁性材料设计平台

A Non-van der Waals Platform for Designing Tunable 2D Magnetic Materials

导读 导读：通过第一性原理 DFT 计算，建立了菱方 Cr₂Se₃ 作为连接性编码的非 vdW 磁性材料设计平台。在位点选择性替代的 Cr₃TMSe₆（TM = V, Mn, W）中，TM 替代作为可控微扰，通过轨道和 SOC 机制系统重塑磁基态、交换网络和磁各向异性。该工作将框架驱动的磁性设计原则从 Heusler 合金扩展到低维非 vdW 体系。

一、前言背景

框架驱动的磁性材料设计

Heusler 合金长期以来是磁性材料设计的模型平台------在刚性结构框架内，不同原子占据特定晶格位点，通过位点选择性化学调控实现丰富的磁性机制。然而，将这种框架驱动的设计范式扩展到低维系统仍面临挑战，因为许多 vdW 磁体缺乏结构复杂性。

菱方 Cr₂X₃（X = S, Se, Te）天然桥接了这两种范式：它们结合了 Heusler 型位点特异性设计的刚性无机晶格和层状非 vdW 几何结构，支持稳健的低维磁性。Cr₂X₃ 由面共享和边共享 CrX₆ 八面体的混合网络构建，具有晶体学不等价金属位点和强共价键。

Cr₂Se₃ 作为磁性设计平台的优势

Cr₂X₃ 可视为 CrX₂（TMD）的自插层版本------vdW 间隙被 Cr 原子填充，提供了 Cr 和 X 原子的结构和化学计量多样性。通过位点选择性化学修饰，磁交换、电子结构和自旋-轨道驱动的各向异性可在强结构约束下系统调控。

核心问题：(1) 过渡金属（V, Mn, W）替代 Cr(III) 位点如何影响基态磁性？(2) 不同 TM 掺杂如何通过轨道和 SOC 机制重塑交换网络和磁各向异性？(3) Cr₂Se₃ 框架能否作为非 vdW 二维磁性材料的通用设计平台？

二、研究方法

第一性原理计算方法

使用 VASP 进行 DFT 计算，PAW 方法，PBE-GGA+U 交换关联泛函。通过线性响应方法确定 Ueff 值：V = 4.43 eV，Mn = 5.87 eV，W = 0.81 eV。平面波截断能 500 eV，力收敛标准 0.002 eV/Å，总能量收敛至 10⁻⁷ eV。

自旋极化计算贯穿始终，SOC 自洽处理用于磁各向异性和轨道分辨性质。声子谱通过 PHONOPY 代码的有限位移法计算（2×2×1 超胞）。AIMD 模拟在 NVT 系综下 300 K 运行 10 ps。

磁交换参数通过四态法总能量映射到经典 Heisenberg 模型提取。CI-NEB 方法计算替代原子的迁移势垒。

经典 Heisenberg 模型：磁交换相互作用的描述

四态法提取交换耦合常数 Jij

稳定性评估的多层次框架

热力学稳定性：形成能 Ef 和凸包上方能量 Ehull------评估合成可行性。

动力学稳定性：声子谱无虚频------确认晶格动力学稳定。

热稳定性：AIMD 模拟在 300 K 下验证结构不退化。

动力学鲁棒性：CI-NEB 迁移势垒 > 1 eV------局部结构重排被动力学抑制，支持非平衡合成（MBE/CVD）。

凸包上方能量：热力学亚稳态的量化指标

【DFT Tip 1】DFT+U的线性响应法：U值不是"调参"而是"计算"

本文使用线性响应方法确定Ueff值（V=4.43 eV, Mn=5.87 eV, W=0.81 eV）------这是最严谨的做法。线性响应法步骤：① 构建超胞，② 对目标原子施加局域势微扰α，③ 自洽计算d轨道占据数n(α)，④ 拟合dn/dα斜率得到U=∂α/∂n。

关键：U值依赖于氧化态、配位环境和自旋态------Cr³⁺在八面体场中Ueff≈3-4 eV，但Mn³⁺因Jahn-Teller畸变Ueff可>5 eV。不能直接使用文献中其他体系的U值。

常见错误：① 直接套用Materials Project的U值（它们基于不同赝势和泛函），② 对5d元素（W）使用大U------5d轨道更离域，Ueff通常<1.5 eV，③ 在SOC计算中同时使用U------SOC+Hubbard U的相互作用需谨慎处理（建议对比U+SOC vs 仅SOC的结果）。

【DFT Tip 2】四态法提取磁交换参数J_ij：自旋配置的陷阱

本文使用四态法（E_↑↑, E_↓↓, E_↑↓, E_↓↑）提取J_ij------这是最常用的方法，但有几个陷阱：① 超胞必须足够大以消除自镜像相互作用（至少2×2×1），② 非共线磁序（如螺旋磁序）的J_ij需用广义四态法（含自旋螺旋计算），③ 自旋态必须与基态一致（如Cr³⁺为S=3/2）。

精度验证：提取的J_ij应能重现所有四种自旋配置的DFT能量（误差<1 meV/atom）。如果不满足，可能原因：① 磁矩在计算中偏离了预设值（检查OSZICAR中磁矩），② 存在非Heisenberg项（如双二次交换或DM相互作用），③ 超胞不够大。

Monte Carlo建议：利用提取的J_ij做经典Monte Carlo模拟预测Tc------但Heisenberg模型中的Tc通常比实验高估20-30%，因为忽略了量子涨落和磁各向异性。

【DFT Tip 3】MAE计算中的SOC收敛：2D磁性材料的精度瓶颈

MAE=E_∥-E_⊥通常在μeV~meV量级------这比DFT总能量（~10³ eV）小数百万倍，对数值精度要求极高。本文使用SOC自洽计算MAE，建议：① 能量收敛标准EDIFF≤1×10⁻⁷ eV（比常规-1×10⁻⁵严格100倍），② K点密度≥15×15×5（2D体系需更密），③ 对比非自洽SOC（在标量相对论波函数基础上加SOC）与自洽SOC结果------差异>10%时优先用自洽SOC。

2D体系的特殊考虑：MAE随层数振荡------单层CrI₃的MAE≈0.7 meV/Cr，而块体≈0.1 meV/Cr。本文单层Cr₃TMSe₆的MAE可能显著不同于块体。另外，真空层厚度不足会导致虚假的层间偶极耦合------建议c≥20 Å。

常见错误：① 用非SOC波函数算SOC能量（VASP中LSORBIT=.TRUE.但未从SOC WAVECAR重启），② 比较不同赝势的MAE（MAE对赝势敏感，特别是包含半芯态的体系），③ 忽略MAE的角度依赖性------本文Fig 6展示了MAE(θ)的全角度扫描，揭示了易轴方向。Cr₂Se₃：非vdW可调二维磁性材料设计平台结构框架Cr₂Se₃菱方晶格TM替代V/Mn/W→Cr(III)稳定性评估声子+AIMD+CI-NEB磁基态AFM→FM转变电子结构Cr-Se骨架主导Jahn-Teller畸变Mn³⁺: t₂g³eg¹内在SOCW-5d: 强SOC交换网络J₁-J₁₂重建MAE增强Mn: Se配体SOCMAE增强W: 内在d轨道SOC单层极限维度调控磁性实验展望MBE/CVD合成▸ 框架拓扑: 面共享+边共享八面体→方向性d-p杂化各向异性→位点选择性调控▸ 三种TM路径: V(带边调控) | Mn(JT畸变→Se-SOC) | W(内在5d-SOC)→同一框架多功能▸ 非vdW优势: 强共价键+动力学鲁棒性(CI-NEB>1eV)+三维磁网络→超越vdW磁体

【知识扩展 1】超交换 vs 双交换：磁性半导体的交换机制之争

【理论解释】超交换（Superexchange）：两个磁性离子通过中间非磁性阴离子（如O²⁻、Se²⁻）的p轨道间接耦合------Anderson（1950）和Goodenough-Kanamori规则（1955-1959）定量描述。J_SE∝t²_pd/U，其中t_pd为p-d跳跃积分，U为Hubbard关联能。对于Cr³⁺-Se²⁻-Cr³⁺路径，180°键角→AFM耦合（t₂g-p-t₂g），90°键角→FM耦合（t₂g-p-e_g）。双交换（Double Exchange）：Zener（1951）提出，不同价态离子间电子跃迁（如Mn³⁺-O²⁻-Mn⁴⁺）→FM耦合，J_DE∝t_eff（有效跳跃积分）。

【方法比较】在DFT中区分SE和DE：① 四态法提取的J_ij的符号和大小------AFM J_ij<0指示SE主导，FM J_ij>0指示DE或直接交换，② 能带结构分析------SE体系中带隙明显，DE体系中费米能级附近态密度高，③ Bader电荷分析------混合价态（Mn³⁺/Mn⁴⁺）是DE的标志。

【经典参考】Anderson, Phys. Rev. 79, 350 (1950)------超交换理论；Goodenough, J. Phys. Chem. Solids 6, 287 (1958)------Goodenough-Kanamori规则；Zener, Phys. Rev. 82, 403 (1951)------双交换；Xiang et al., Dalton Trans. 42, 823 (2013)------四态法DFT提取J_ij。

【迁移能力】本文Cr₃TMSe₆的交换机制可通过Bader电荷+轨道分辨PDOS分析：Cr³⁺(3d³)之间为SE→AFM，TM替代引入额外d电子→部分DE→FM。此分析框架适用于所有过渡金属氧化物和硫族化物磁性体系。

三、实验结果

图 1：菱方 Cr₃TMSe₆ 晶体结构。(a) 侧视图和 (b) 俯视图------交替的完全占据 CrSe 层和 TM/Cr 替代层形成非 vdW 堆叠。(c) 面共享和边共享八面体连接性------Cr(I) 和 Cr(II) 位点。(d) CrSe₆ 八面体排列------Se 原子按沿 c 轴方向编号（1-3 在上，4-6 在下），标注键长 la, lb 和键角 θa, θb。

结构与热力学性质

TM 替代后所有结构保持 Cr₂Se₃ 的菱方晶格对称性，晶格参数 a 和 c 仅有适度变化。Mn 替代沿 c 轴膨胀最大，W 替代沿 c 轴收缩。

所有化合物的形成能 Ef 为负（~-0.6 eV），表明替代相的热力学可及性。Cr₃MnSe₆ 接近热力学稳定极限，Cr₃VSe₆ 和 Cr₃WSe₆ 分别为亚稳态和高亚稳态，但在非平衡生长条件下仍实验可行。

TM 替代将所有 Cr₃TMSe₆ 体系的基态稳定为铁磁（FM），不同于原始 Cr₂Se₃ 的补偿反铁磁（AFM）序------替代 Cr(III) 位点有效改变了竞争交换相互作用的平衡。

图 2：声子色散、ELF 和 AIMD 稳定性。(a-d) 声子色散关系和原子投影声子态密度------无虚频，动力学稳定。(e-h) (110) 面 ELF 图------TM 替代系统修改键合特征：V 增加局域化，Mn 显示更显著的离子性键合，W 显示更离域的金属性键合。(i) AIMD 模拟在 300 K 下 10 ps------能量和温度波动可忽略，热稳定性良好。

电子结构与轨道杂化

整个 Cr₃TMSe₆ 系列的前沿电子态仍由 Cr-Se 网络主导，掺杂剂/插层剂作为轨道和能量选择性微扰，而非带边形成物种。

V 对导带边（CBM）的扰动最强------V 3d 轨道与上层 Cr-d 态有效杂化。Mn 的 3d 能级远离带边，在 ±2 eV 内几乎无贡献，但 Mn-Se 键长显著增加（+8%）导致 Cr-Se 键长伸长，产生半金属性。

W 5d 轨道更离域，与面外 Se-pz 态杂化更强------面共享 Cr(I)Se₆ 和 WS₆ 八面体沿 c 轴的 dz²-pz 耦合通道。W 略微增强价带权重但不主导。

图 3：CI-NEB 迁移势垒------Cr₃TMSe₆（TM=V, Mn, W）的动力学鲁棒性。反应坐标连接初始和最终构型。绿色虚线为原始 Cr₂Se₃ 中 Cr(III) 迁移路径。所有势垒 > 1 eV（远超室温热能 kBT ~ 0.026 eV），局部结构重排被动力学抑制。

图 4：Cr₃TMSe₆ 中 Cr(I) 位点周围 CrSe₆ 八面体的结构畸变。(a) 键长不对称性 Δl，(b) 键角不对称性 Δθ，(c) 面外角畸变 Δγ。Mn 替代产生最大局部畸变------高自旋 Mn³⁺（t₂g³eg¹）Jahn-Teller 效应。W 替代仅显示最小偏离------离域 5d 电子密度减少金属-配体各向异性。

磁各向异性能（MAE）的微观机制

所有化合物的易轴在面内。总 MAE 从 V 到 Mn 至 W 单调递增。Mn 替代产生最大的 Se 衍生 MAE------三个 Se 位点贡献相当且显著。W 替代的大 MAE 几乎完全来自 W 位点本身，而非 Se。

Mn 增强 MAE 主要通过畸变激活的 Se 介导 SOC 通道------Jahn-Teller 畸变增强位点依赖杂化。W 通过其强内在 5d SOC 实现大各向异性------dz²-dyz 和 dxy-dx²-y² 通道的强非对角耦合。

两种 MAE 增强机制提供不同功能路径：畸变敏感各向异性可用于应变/电场控制磁结构，重元素内在 SOC 提供高稳定性磁存储。

磁各向异性能：面内与面外能量差

图 5：Cr₃TMSe₆ 的电子能带结构。(a) Cr₃VSe₆，(c) Cr₃MnSe₆，(e) Cr₃WSe₆------自旋向上（红色虚线）和自旋向下（蓝色实线）。气泡大小定性表示 3d-3d/3d-5d 轨道杂化强度。(b,d,f,h) 总态密度和元素投影 PDOS。(i-l) 掺杂位点的轨道分辨 PDOS------dz²(pz) 红色、dxy+dx²-y²(py) 灰色、dxz+yz(px) 蓝色。

图 6：磁各向异性能的角度依赖性和轨道分辨 SOC 矩阵。(a) MAE 随磁化方向角度的变化。(b) Cr、Se 和 TM 对总 MAE 的相对原子贡献。(c-f) p 轨道和 d 轨道的 SOC 能量矩阵------揭示 Mn 通过 Se 配体 SOC 和 W 通过内在 d 轨道 SOC 增强 MAE 的不同机制。

磁交换网络的调控

原始 Cr₂Se₃ 具有复杂的交换网络------多个显著的 AFM 相互作用跨越层内和层间路径，导致 frustrated AFM 基态。TM 替代从根本上重构交换拓扑：改变多个交换通道的符号和大小。

Cr₃VSe₆ 和 Cr₃MnSe₆ 中多个竞争 AFM 相互作用被削弱或反转，降低磁阻挫程度，稳定 FM 序。Cr₃WSe₆ 显示独特的交换层次------J₅ 和 J₈ 主导，扩展 5d 轨道和增强金属-配体杂化沿 c 轴抑制竞争路径。

TM 替代诱导的杂化增强 FM 交换通道（双交换机制）并抑制 AFM 超交换------这是 AFM → FM 转变的微观起源。

超交换（J ∝ t²pd）与双交换（J ∝ teff）竞争

图 7：Cr₃TMSe₆ 中 Cr 和 TM 位点间磁交换路径（J₁-J₁₂）的示意图------左边三维单胞表示，右边二维表示。Cr(I) 浅蓝、Cr(II) 青色、Cr(III)/TM 红色。层内和层间交换耦合共同构成三维连接性磁网络。

图 8：Cr₂X₃ 框架：从基础设计到实验展望。(左) Cr₂Se₃ 作为连接性编码的结构骨架------Cr(III) 层作为化学活性替代位点，MBE 层-by-层生长。(右) 未来展望------非 vdW 系统的结构完整性和动力学鲁棒性支持 MBE/CVD 大尺寸合成，利用可调各向异性机制和重构交换拓扑工程化 Curie 温度和自旋动力学。

【DFT Tip 4】2D磁性材料的声子谱计算：虚频排查与修复

本文所有Cr₃TMSe₆声子谱无虚频------但2D磁性材料中虚频极其常见。常见原因：① 结构未充分弛豫（力标准>0.001 eV/Å），② 超胞太小（声子力常数在2×2×1超胞中可能未收敛→建议3×3×1），③ 磁性基态设置错误（如AFM用FM计算声子），④ 真空层太薄（<15 Å）。

修复虚频策略：① 如果虚频仅在Γ点附近→用Phonopy的NAC修正（非极性材料不需要），② 如果虚频在整个BZ→结构弛豫未收敛，重新弛豫（EDIFFG≤-0.001），③ 如果虚频在特定q点→检查该q点对应的实空间原子位移模式，手动调整原子位置后重新弛豫。

2D体系的特殊处理：对于2D材料，Phonopy中需设置DIM=2×2×1（而不包括真空方向），且phonopy.conf中需添加DIM=...正确指定超胞。常见的错误：在3D Phonopy中计算2D材料→真空方向产生虚假的声学支折返。

【DFT Tip 5】AIMD模拟的实用参数：10 ps够不够？

本文使用300 K下10 ps的AIMD验证热稳定性，能量波动<22 meV/atom。但10 ps对以下情况可能不够：① 相变动力学（需>100 ps），② 离子扩散（需>50 ps），③ 磁矩翻转（需>20 ps+自旋翻转采样）。

AIMD参数建议：① 时间步长1-2 fs（轻元素H需0.5 fs），② NVT系综+Nosé-Hoover恒温器，③ 前2-5 ps用于平衡（丢弃），后5-10 ps用于统计，④ 只对Γ点采样（KPOINTS=1×1×1）以降低计算成本。

实用检查：如果AIMD中总能量漂移>0.1 eV/atom/ps→结构不稳定或时间步长太大。如果温度波动>±50 K→恒温器参数不当。对于磁性体系，额外检查：AIMD中磁矩是否保持------如果磁矩在300 K下消失→居里温度低于室温，材料不适合室温应用。

【DFT Tip 6】CI-NEB迁移势垒计算：收敛标准与插点策略

本文CI-NEB给出>1 eV的迁移势垒，确认动力学稳定性。CI-NEB的关键参数：① 插点数量（images）------本文8个，对于简单路径足够，复杂路径需10-15个，② 弹簧常数（SPRING=-5）------默认值通常OK，③ 收敛标准（EDIFFG≤-0.02 eV/Å）------比几何弛豫宽松。

常见CI-NEB失败原因：① 初始路径不合理------插点直接线性插值导致原子重叠，② 插点陷入局部极小值------用improved tangent方法或增加插点数，③ 磁性体系的自旋态在路径中改变------需固定MAGMOM。

对于2D体系：CI-NEB中超胞的真空层必须足够（>15 Å），否则周期性镜像原子会引入虚假相互作用，导致势垒偏低或偏高。另外，CI-NEB仅考虑单原子迁移------协同迁移（多个原子同时移动）需要更复杂的过渡态搜索方法。

【DFT Tip 7】应变工程在2D磁性中的DFT实现：弹性极限与MAE响应

本文在单层极限下进一步研究了应变效应------2D磁性对应变极其敏感。VASP中施加应变：① 手动修改POSCAR的晶格矢量（固定应变），② 或使用ISIF=3弛豫面内晶格+固定c轴。建议：面内应变ε=±2%, ±4%, ±6%扫描，每个应变下重新弛豫原子位置。

应变效应分析：① 带隙vs应变→半导体-金属转变，② MAE vs应变→易轴翻转（面内↔面外），③ J_ij vs应变→FM/AFM转变。关键参数：弹性极限（应力-应变曲线的线性范围）------对于2D硫族化物通常ε<10%。

常见错误：① 施加应变后不做原子弛豫（仅缩放晶格）→高估应力，② 用ISIF=3做2D应变计算（会同时改变真空层），③ 忽略泊松效应------单轴应变会引起正交方向的被动应变。

【科研经验 1】为什么DFT预测的磁性基态与实验不符？

问题：DFT+U预测Cr₂Se₃为FiM基态，但实验可能报告不同磁序。原因：① Ueff值的选择直接影响磁基态------U偏大→AFM稳定，U偏小→FM稳定，② 四态法提取的J_ij对U值极其敏感（J_ij随U变化可达±50%），③ DFT忽略了自旋-轨道耦合的量子涨落效应（在2D体系中尤其重要）。

解决方案：① 对Ueff做收敛性测试（Ueff=2,3,4,5 eV→观察磁基态是否变化），② 用HSE06对比验证（混合泛函通常给出更准确的磁基态，但成本高），③ 如果有实验数据（如中子衍射磁结构），用实验磁序反推Ueff。

建议：DFT预测的磁基态应视为"候选"而非"结论"------尤其对于竞争磁序能量差<10 meV/f.u.的体系，不同泛函结果可能不同。

【科研经验 2】为什么声子谱无虚频但AIMD显示结构不稳定？

问题：声子谱（T=0 K的简谐近似）无虚频，但AIMD（T=300 K）显示结构退化。原因：① 非简谐效应------高温下声子-声子散射导致有效势软化，② 声子谱在软模附近能量极小值很浅→热涨落即可跨越势垒，③ 磁性体系的磁-声子耦合------磁矩波动导致动态晶格不稳定。

解决方案：① 在多个温度下做AIMD（300 K, 500 K, 800 K）→确定结构失稳温度，② 计算声子谱的Grüneisen参数------大的负γ指示强非简谐性，③ 用SSCHA（随机自洽简谐近似）方法计算有限温度下的有效声子谱。

建议：对于2D磁性材料，声子谱（T=0 K）和AIMD（T=300 K）的一致性验证是必需的------单靠声子谱不足以确认实际工作温度下的稳定性。

【知识扩展 2】Jahn-Teller效应：从配位化学到磁性材料设计

【理论解释】Jahn-Teller（JT）定理（1937）：非线性分子在简并电子基态下会发生几何畸变以降低能量。在八面体场中：d⁴（HS Mn³⁺: t₂g³eg¹）和d⁹（Cu²⁺: t₂g⁶eg³）产生强JT畸变------eg轨道简并导致轴向伸长/压缩。本文Mn³⁺（3d⁴）替代Cr³⁺（3d³）→引入JT活性→八面体畸变（Δl, Δθ显著增大）→通过Se配体SOC增强MAE。

【方法比较】在DFT中检测JT畸变：① 八面体键长不对称性Δl=la-lb，② 键角不对称性Δθ=θa-θb，③ 轨道分辨PDOS中eg轨道的分裂。JT畸变程度：Mn³⁺（强）> Cr²⁺（中）> Cu²⁺（强）> Ni³⁺（弱）。

【经典参考】Jahn & Teller, Proc. R. Soc. A 161, 220 (1937)------原始定理；Bersuker, Chem. Rev. 101, 1067 (2001)------JT效应综述；Kugel & Khomskii, Sov. Phys. Usp. 25, 231 (1982)------轨道有序与JT效应。

【迁移能力】JT效应不仅是Mn³⁺的专利------任何d⁴/d⁷（HS）或d⁹体系都有JT活性。在磁性材料设计中，可主动利用JT畸变来增强MAE：选择JT活性离子（Mn³⁺, Cu²⁺, Cr²⁺）替代非JT活性位点→畸变驱动SOC增强→MAE可调。此策略适用于钙钛矿、尖晶石和层状硫族化物。

从非vdW框架到低维磁性材料基因组

本文已完成了Cr₂Se₃框架的DFT基础表征------结构、稳定性、电子结构、磁交换和MAE。但作为非vdW磁性设计平台，以下计算可以进一步挖掘其潜力，构建完整的低维磁性材料设计链条。

Further 1 | Monte Carlo模拟精确预测Curie温度

为什么值得算：本文仅定性讨论了TM替代对磁性的影响，但缺乏Tc的定量预测。利用已提取的J_ij参数，做经典Monte Carlo（Metropolis算法+Heisenberg模型）模拟→预测Tc。

能回答的问题：Cr₃VSe₆、Cr₃MnSe₆、Cr₃WSe₆的Tc分别是多少？是否高于室温？维度效应（块体→单层）对Tc的影响？

适合体系：所有Cr₃TMSe₆（块体+单层），以及扩展到Cr₂S₃和Cr₂Te₃。

输入：四态法提取的J_ij（已完成），自旋S值（Cr³⁺: S=3/2, Mn³⁺: S=2, V³⁺: S=1, W⁴⁺: S=0），磁性各向异性常数。

Further 2 | 自旋波色散与磁振子谱

为什么值得算：本文的交换网络包含12个J_ij参数------复杂的磁阻挫结构可能产生非常规磁振子（如Dirac磁振子、拓扑磁振子）。计算自旋波色散ω(q)可以揭示磁激发的拓扑性质。

能回答的问题：Cr₃TMSe₆中是否存在Dirac磁振子或Weyl磁振子？TM替代如何改变磁振子能隙？磁振子-声子耦合（磁振子极化子）？

适合体系：所有具有复杂交换网络的磁性体系，特别是存在竞争AFM/FM交换的体系。

输入：J_ij矩阵（已完成），SpinW或TBiLab代码，线性自旋波理论（LSWT）。

Further 3 | 光学性质与磁光Kerr效应（MOKE）

为什么值得算：本文的电子结构显示半金属到半导体的转变------这意味着丰富的光学响应。计算介电函数ε(ω)、吸收系数α(ω)和MOKE旋转角θ_K。

能回答的问题：Cr₃TMSe₆的MOKE信号是否足够强以用于磁光存储？半金属性的光学指纹是什么？

适合体系：所有Cr₃TMSe₆（块体+单层），特别是Mn替代体系（大MAE→大MOKE）。

输入：SOC自洽的WAVECAR（已完成），VASP介电函数计算（LOPTICS=.TRUE.），K点密度≥15×15×5。

Further 4 | 磁振子热导率与自旋Seebeck效应

为什么值得算：非vdW体系的三维磁网络连接性赋予其独特的热输运性质。计算磁振子热导率κ_mag(T)和自旋Seebeck系数S_s。

能回答的问题：Cr₃TMSe₆的磁振子热导率是否可与YIG竞争？Mn替代（大畸变→强磁振子-声子散射）vs W替代（强SOC→大自旋Seebeck）哪种更适合热电器件？

适合体系：所有磁性绝缘体/半导体，特别是具有大MAE和强SOC的体系。

输入：磁振子色散ω(q)（从Further 2），弛豫时间近似，Boltzmann输运方程。

Further 5 | 非共线磁结构与DM相互作用

为什么值得算：本文假设共线FM/AFM磁序，但非vdW体系中DM相互作用（Dzyaloshinskii-Moriya）可能导致螺旋磁序或斯格明子。计算DM矢量D_ij和螺旋磁序的能量。

能回答的问题：Cr₃TMSe₆中是否存在非共线磁基态？DM相互作用是否足以稳定斯格明子？W替代（强SOC）是否增强DM？

适合体系：所有具有SOC和反演对称性破缺的磁性体系（Cr₃TMSe₆满足面内反演对称性破缺）。

输入：四态法提取DM矢量（需非共线自旋配置），VASP非共线磁计算（LNONCOLLINEAR=.TRUE.），广义Bloch定理。

Further 6 | 高通量筛选Cr₂X₃框架：从Cr₂Se₃到Cr₂S₃和Cr₂Te₃

为什么值得算：本文仅覆盖Cr₂Se₃+3种TM（V, Mn, W），但扩展相空间巨大：X=S, Se, Te（3种）×TM=3d,4d,5d系列（~30种）×替代浓度（25%, 50%, 75%）=~270种候选。高通量DFT筛选可发现最优Tc和MAE组合。

能回答的问题：Cr₂Te₃框架的Tc是否更高（Te的SOC更强）？哪种TM替代产生最高的MAE？是否存在多功能材料（高Tc+大MAE+半金属性）？

适合体系：所有Cr₂X₃衍生化合物。

输入：自动弛豫+四态法J_ij+MAE+MC模拟，Atomate或AFlow自动化。

单层极限与维度调控

所有 Cr₃TMSe₆ 体系均可从块体剥离至单层。单层中所有体系保持 FM（Cr₂Se₃ 为 FiM），大多数单层变为金属性，Cr₃WSe₆ 显示半金属性。

维度缩减/扩展为调控电子和磁性提供了额外自由度------Cr-Se 框架保持完整，但带特性和磁各向异性发生显著变化，强化了 Cr₂Se₃ 作为多功能低维材料工程平台的多样性。

四、对比分析

Cr₂Se₃ vs Heusler 合金 vs vdW 磁体

Heusler 合金：刚性立方框架，位点选择性化学调控，但限于三维。Cr₂Se₃：类 Heusler 的位点特异性设计 + 层状非 vdW 几何 = 低维磁性工程。

vdW 磁体（CrI₃, Fe₃GeTe₂）：层间弱耦合，环境稳定性差，易降解。Cr₂Se₃：强共价键，非 vdW 堆叠，结构鲁棒性更高。

Cr₂Se₃ 的核心优势：面共享和边共享八面体连接性定义了方向性 d-p 杂化各向异性，为位点选择性调控电子和磁性行为提供了结构编码平台。

V、Mn、W 替代的差异化机制

V：中等 MAE，弱 SOC 和适度杂化------CBM 扰动最强。

Mn：最大八面体畸变（Jahn-Teller），通过畸变激活 Se 配体 SOC 通道增强 MAE------离子性键合，层间分离最大。

W：最小畸变，但最强内在 SOC------大 MAE 几乎完全来自 W 位点，W-Se 面外杂化增强层间耦合------金属性键合。

三种 TM 提供不同功能路径：应变敏感（Mn）、SOC 增强（W）、带边调控（V）------同一框架内实现多功能性。

五、讨论

非 vdW 二维磁体的设计范式

Cr₂Se₃ 建立了连接性编码的无机设计平台：框架拓扑（面共享 vs 边共享八面体）决定 d-p 杂化方向性，位点化学（TM 替代）作为可控微扰。

与 vdW 磁体相比，非 vdW 体系具有更强的层间耦合和三维磁网络连接性------这既是优势（更高稳定性）也是挑战（更复杂的交换网络）。

单层极限进一步扩展了设计空间：维度效应可显著改变带特性和磁各向异性，为非 vdW 异质结中的邻近效应和界面工程提供新机会。

自旋电子学应用前景

结构可调的配体介导各向异性（Mn 替代）→ 应变/电场控制的磁结构，电压可调磁存储。

内在重元素 SOC（W 替代）→ 高稳定性磁存储应用。

同一非 vdW 框架内可调性 vs 鲁棒性的权衡 → 多功能自旋电子学器件。

未来方向：实验合成验证（MBE/CVD），Monte Carlo 模拟精确预测 Tc，探索 Cr₂S₃ 和 Cr₂Te₃ 的扩展相空间。

六、总结

核心结论

(1) Cr₂Se₃ 是连接性编码的非 vdW 磁性设计平台------面共享和边共享八面体定义方向性杂化各向异性。

(2) TM 替代（V, Mn, W）在 Cr(III) 位点系统调控磁基态、交换网络和磁各向异性------所有体系从 AFM 转为 FM。

(3) Mn 通过 Jahn-Teller 畸变激活 Se 配体 SOC 增强 MAE，W 通过强内在 5d SOC 实现大 MAE------两种不同机制。

(4) 所有 Cr₃TMSe₆ 动力学、热力学和动力学稳定------CI-NEB 势垒 > 1 eV，支持非平衡合成。

(5) 单层极限提供额外维度调控自由度------Cr-Se 框架保持完整，带特性和磁各向异性显著变化。

展望

实验合成 Cr₃TMSe₆ 薄膜并验证磁性预测。

Monte Carlo 模拟精确预测各体系的 Curie/Néel 温度。

将设计范式扩展到 Cr₂S₃ 和 Cr₂Te₃，探索非 vdW 异质结中的邻近效应。

七、支撑信息

SI 图概览（S1-S5）与附表（S1-S3）

SI 提供 5 张补充图和 3 张补充表：S1------Hubbard U 参数的线性响应计算，S2------Cr₂Se₃/Cr₃TMSe₆ 菱方单胞中 Cr 和 TM 原子的四种可能自旋排列，S3------位点和轨道分辨投影态密度（PDOS），S4------原子分辨磁各向异性能（MAE）贡献，S5------单层 Cr₃TMSe₆ 的自旋密度分布和自旋极化能带结构。附表 S1------晶格参数与磁基态汇总，S2------Wyckoff 位置，S3------不同相的相对能量。

图 S1：Hubbard U 参数的线性响应计算。(a)-(c) 替代 TM 位点 d 轨道占据数随外加局域势的变化------线性响应斜率确定 Ueff 值（V=4.43 eV, Mn=5.87 eV, W=0.81 eV）。

图 S2：Cr₂Se₃/Cr₃TMSe₆ 菱方单胞中 Cr 和 TM 原子的四种可能自旋排列。(a) 铁磁（↑↑↑），(b) 亚铁磁-II（↑↑↓），(c) 亚铁磁-I（↑↓↑），(d) 反铁磁（↑↓↓）------确定 FM 为基态。

图 S3：Cr₂Se₃ 和 Cr₃TMSe₆（TM=V, Mn, W）的位点和轨道分辨投影态密度（PDOS）。(a-d) Cr₂Se₃，(e-h) Cr₃VSe₆，(i-l) Cr₃MnSe₆，(m-p) Cr₃WSe₆------展示 Cr-d、TM-d 和 Se-p 轨道的能量分布。

图 S4：原子分辨磁各向异性能（MAE）贡献。(a) Cr₃VSe₆，(b) Cr₃MnSe₆，(c) Cr₃WSe₆。蓝色、绿色和红色柱分别代表 Cr、Se 和 TM（V/Mn/W）的 MAE 贡献------揭示 Mn 通过 Se 配体和 W 通过内在 d 轨道 SOC 增强 MAE 的不同机制。

图 S5：单层 Cr₃TMSe₆ 的自旋密度分布和自旋极化能带结构。(a, b) 自旋密度分布（青色和紫色代表相反自旋符号），(c-f) TM=(c)Cr, (d)V, (e)Mn, (f)W 的自旋极化能带结构------展示单层极限下的电子和磁性特征。

主要参考文献

核心引用

1 Gebredingle Y, You S, Kim H, Kim N. Sci. Rep. 16, 19079 (2026) --- 本工作。

2 Liu X, Gebredingle Y, Guo X, Zhang F, Kim N. Adv. Funct. Mater. 34, 2316834 (2024) --- Cr 硫族化物薄膜综述。

3 Zhao X et al. Nature 581, 171 (2020) --- 自插层工程化共价键二维材料。

4 Gebredingle Y, Joe M, Lee C. ACS Appl. Nano Mater. 5, 10383 (2022) --- Cr₂X₃ 自旋电子学。

5 Gebredingle Y, Joe M, Lee C. Phys. Rev. Mater. 6, 054405 (2022) --- Cr₂S₃ 双极半金属 FiM。

6 Xiang HJ et al. Phys. Rev. B 84, 224429 (2011) --- 四态法磁交换参数提取。

7 Henkelman G et al. J. Chem. Phys. 113, 9901 (2000) --- CI-NEB 方法。

Gebredingle Y et al. | Scientific Reports 16, 19079 (2026) | Cr₂Se₃ · 非vdW · 二维磁性 · 位点选择性替代