DFT如何计算交换耦合系数？

说明：本文华算科技主要介绍多层磁性材料中交换耦合系数 J 的物理含义、DFT 总能差怎样映射到 Heisenberg 模型、层间堆垛和应变怎样改变 J，以及报告计算结果时需要限定的模型条件。

交换耦合系数J是什么？

多层磁性材料中，每个磁性原子或磁性离子附近可以形成局域磁矩。DFT 给出的直接结果是不同自旋排布下的总能、磁矩、电子态和力；交换耦合系数 J 是把这些量投影到自旋模型后的参数，常用于描述两个局域磁矩倾向于平行还是反平行排列。

常见写法是 Heisenberg 型哈密顿量，例如 H = -Σ_ij J_ijS_i·S_j。在这一约定下，J_ij > 0 对应 FM 耦合，J_ij ij J_ijS_i·S_j，符号会反过来。读 J 值时，哈密顿量符号、是否包含 1/2、S 是单位矢量还是自旋量子数 S，都要与数值绑定。

双层或多层材料比单层多出一类变量：层间磁矩之间的配对关系。同一层内可以有 J₁、J₂、J₃，相邻层之间可能写成 J_z、J_c、J₁¹² 或 J₂¹²。若只把所有层间作用压成一个有效 J_z，文章应说明这是平均后的模型参数，不能替代每一条 Cr-I-I-Cr、M-X-X-M 或 M-X-M 路径的真实差异。

多层材料的 J 还与磁性构型有关。双层 CrI₃、CrCl₃ 这类范德华磁体里，层内 FM 和层间 AFM 可以同时出现；三层材料还会产生 ↑↓↑、↑↑↓、↑↑↑ 等排布。DFT 计算 J 时比较的是同一个结构模型中的一组自旋态能量，不是单独看某一次自旋极化计算的总磁矩。

怎样把FM/AFM总能差换算成 J？

能量映射法从几个固定自旋构型出发。先选定晶胞、堆垛、层数、磁性原子数和需要求解的配对项，再计算 FM、A-type AFM、G-type AFM 或更多非等价自旋态的总能。每个总能都可以写成 E₀ 加上一串 J_ijS_i·S_j 项，不同构型给出一组线性方程。未知 J 的数量越多，需要的独立磁构型越多。

以最简双层模型为例，若只保留层间最近邻 J_z，并采用 H = -ΣJ_ije_i·e_j 的单位矢量写法，那么一对层间磁矩从 FM 变成 AFM 时，能量项由 -J_z 变成 +J_z，每一对贡献变化 2J_z。若超胞中有 N_pair 个等价层间配对，ΔE = E_AFM – E_FM = 2N_pairJ_z。换成量子自旋 S 时，需要乘 S² 或 S(S+1) 形式，取决于映射约定。

一个双层原胞里有几个 Cr、每个 Cr 跨层连到几个 Cr、是否双计数、能量是每个原胞还是每个磁性原子，都会直接改变 J 的分母。报告 J 时写成 meV/Cr、meV/pair 或 meV/unit cell，含义并不相同。只给“AFM 比 FM 低 8 meV”还不能复现 J，必须同时给出磁性原子数和配对计数。

含多个未知量时，FM/AFM 两个能量差不够。层内 J₁、J₂、J₃ 与层间 J_z 共同存在时，可以设计多种超胞自旋排布，让每个构型中的平行/反平行配对数不同，再解线性方程；构型数多于未知数时，用最小二乘拟合。拟合残差若达到 J 本身同一量级，说明简单 Heisenberg 模型缺项，可能需要加入双二次交换、DMI、Kitaev 项、各向异性交换或更远邻配对。

层数、堆垛和应变为什么会改写J？

层间交换通常比层内交换小得多，因为相邻层之间没有强共价骨架，磁性原子间隔更长，电子轨道重叠依赖卤素、硫族元素或有机配体连接。堆垛改变的是磁性原子之间的相对横向位置，同一层距下也会改变 M-X-X-M 或 M-X-M 的角度、距离和中间配体轨道参与程度。

压缩或拉伸应变改写键长、键角、层间距和晶场分裂后，CrCl₃ 这类体系中的 FM 与 AFM 能量差会移动，GGA 与 GGA+U 对临界应变的估计也不同。这里能读出的不是“某个应变一定触发相变”，而是在指定泛函、U、vdW 修正和堆垛下，E_AFM-E_FM 的符号和斜率已经对结构变量敏感。

交换参数还能分解为轨道贡献。Cr 的 t_2g、e_g 轨道与 Cl 或 I 的 p 轨道形成不同的超交换路径，某些轨道组合给 FM 贡献，另一些轨道组合给 AFM 贡献。轨道分解图的价值在于区分 J 的来源：同一个总 J 可能来自几个相反符号通道的抵消，轻微应变或 U 值变化就可能改变抵消后的净值。

近邻配对求和也要谨慎。图中的实线和虚线分别对应不同的配对汇总方式，层内最近邻、次近邻、第三近邻和层间多条路径对总交换的权重不同。若文章只取一个 J_z，读者需要知道它是最近邻层间配对、所有层间配对求和，还是为了后续 Monte Carlo 模拟而重标定的有效参数。

CrI₃ 中的 Cr-I-I-Cr 层间超超交换说明了另一点：层间 J 并非只由层距决定。I 原子子晶格破坏某些镜面对称后，两个 Cr-Cr 距离相同的层间路径也可能给出不同的 J。多层材料的交换耦合需要检查堆垛、配体连接路径和电子关联处理，单独改变真空层或层距后比较总能，容易把几何效应和电子结构效应误归入同一个 ΔE。

DFT计算J时哪些条件必须同步报告？

J 的数值和方法设置紧密相连。含 Cr、Mn、Fe、Co、Ni 等元素的层状磁体通常要检查自旋极化、U 值、vdW 修正、SOC 和非共线设置。U 改变 d 轨道占据和 p-d 能级差，vdW 修正改变层间距，SOC 决定磁各向异性和 DMI；交换耦合常数需要绑定结构、电子关联和自旋方向。

DFT+U、QSGW 和 QSGW+U 对自旋激发谱的差异提醒我们：能量映射得到的 J 只是某种电子结构方法下的模型参数。某些体系里，DFT 可能高估交换耦合；U 的加入有时降低线性响应中的最近邻 J，却可能在总能映射里给出相反趋势。比较 T_C、磁各向异性能和自旋波时，J 定义要保持一致。

结构弛豫口径也要写明。若每个磁构型都独立弛豫，ΔE 会混入磁致伸缩和构型弛豫能；若固定同一个晶格和原子坐标，J 更接近给定结构下的自旋模型参数。研究磁性基态时可以比较完全弛豫后的能量，研究交换常数时更常把结构固定在同一参考构型，再翻转自旋排布。

双层的 J_z 不能直接替代三层、四层或块体中的 J_c，因为层数和周期条件会改变层间配对；周期性边界会让最上层与最下层在 z 方向重复相邻，真空层不足时还会引入非物理层间作用。J 的单位、符号约定、配对数、每层磁矩数、U 值、vdW 方法、SOC 状态和结构是否固定应与表格或图注同页出现。

把 J 用于 Monte Carlo、Vampire、UppASD 或自旋波计算时，需要确认程序读取的自旋是单位矢量、磁矩 μ_B 还是量子自旋 S。DFT 表格中的 meV/pair 若直接填入自旋动力学软件，可能需要乘以或除以磁矩大小、配位数或双计数因子。最终比较对象若是 T_C、T_N、磁滞回线或 magnon dispersion，交换参数表应保留到能复算这些量的程度。