KKR方法计算磁性体系交换常数
KKR方法计算磁性体系交换常数
KKR方法计算磁性体系交换常数
一、引言
KKR（Korringa-Kohn-Rostoker）方法计算磁性体系交换常数，通常采用Liechtenstein–Katsnelson–Antropov–Gubanov（LKAG）公式。该方法通过自旋向上/向下电子的Green函数，计算体系总能量对两个磁性原子磁矩相对旋转的二阶响应，从而得到交换常数Jij。
与基于总能量的方法不同，LKAG不需要构建多个磁构型进行能量差计算，而是直接从电子Green函数导出Jij。
二、Heisenberg模型与交换常数定义
磁性体系通常映射为Heisenberg模型：
其中i, j表示两个磁性原子，ei、ej为局域磁矩方向的单位矢量，Jij为交换常数。
符号约定：
注：不同程序的Hamiltonian符号约定可能存在差异，实际分析时须参照具体软件定义。
三、LKAG方法核心思想
LKAG方法的关键在于考虑两个局域磁矩发生无穷小旋转：
ei →ei + δei
ej →ej + δej
体系总能量展开为：
由此得到：
交换常数本质上即为总能量对两个局域磁矩相对旋转的二阶导数。KKR方法的优势在于，该二阶导数可用Green函数直接表达。
四、LKAG公式
经典的LKAG公式表达形式为：
4.1 交换分裂 Δi
Δi表示第i个原子的交换分裂：
其中ti↑(E)和ti↓(E)分别为自旋向上和向下的单原子t矩阵。直观上：
4.2 Green函数 Gij
Gij(E)表示从原子i传播到原子j的Green函数：
LKAG公式计算的物理过程为：磁性原子i → 电子传播 → 原子j → 再传播回i，这一过程对磁矩耦合的贡献。
五、Green函数与交换作用的物理图像
设体系中有两个磁性原子i和j。以3d电子为例，其具有交换分裂，即自旋向上与向下的电子散射势不同。电子从i传播到j（Gij），受到j的磁性散射（Δj），再传播回来（Gji），产生耦合项：
该量包含了两个磁矩之间的电子介导耦合。因此，KKR-LKAG特别适合研究以下机制：
RKKY相互作用
超交换作用
双交换作用
巡游磁性
长程交换作用
上述机制本质上均与电子传播相关。
六、KKR方法的优势
KKR基于多重散射理论（Multiple Scattering Theory），将晶体视为电子在各原子间不断散射的过程。KKR直接构造实空间Green函数G(r,r′, E)，天然拥有电子从某一原子传播到另一原子的信息。
由于Jij恰好是这种信息的二阶响应，LKAG与KKR的结合具有自然的理论自洽性。
七、实际计算流程
以NiO体系为例，完整计算流程如下：
建立晶体结构
KKR自洽计算（spin-polarized）
得到电荷密度与有效势
得到自旋极化电子结构
计算单原子t矩阵
得到散射路径算符τij
构造Green函数Gij
LKAG积分计算Jij
按距离/原子对分类整理
得到J1, J2, J3, ...
7.1 参考磁性状态
计算前需指定参考磁态，例如：
铁磁（FM）：所有磁矩同向排列
反铁磁（AFM）：磁矩交替反向排列
7.2 自洽计算
进行自旋极化KKR计算，分别得到有效势Veff↑(r)和Veff↓(r)，进而得到各原子的局域磁矩：
7.3 散射矩阵与Green函数
每个原子对应单原子散射矩阵ti(E)。对于磁性原子，ti↑(E) ≠ ti↓(E)。通过散射路径算符τij(E)：
τij包含电子从i到j的所有多重散射路径，因此可描述复杂的电子传播过程，这也是KKR能够得到长程Jij的原因。
7.4 LKAG积分
实际程序采用复能量积分或围道积分等数值方法，而非直接从−∞积到EF。该步骤为计算量较大的环节。
7.5 结果分析
最终得到不同近邻距离的交换常数，例如：
J1 = 8.2 meV（最近邻）
J2 = −1.4 meV（次近邻）
J3 = 0.37 meV（第三近邻）
可进一步绘制Jij随距离rij的变化关系：
该图对于判断磁性机制具有重要价值。
八、KKR-LKAG与VASP四态法的对比
8.1 VASP四态法
人为建立不同磁构型（如↑↑和↑↓），通过总能量差计算交换常数：
该方法需要计算多个有限磁构型的总能量，通过能量差拟合得到J。
8.2 KKR-LKAG
仅需一个参考磁态，通过Green函数直接计算无穷小旋转下的响应：
无需构建FM、AFM1、AFM2等多个构型，一个自洽参考态即可得到大量Jij。
8.3 方法特性差异
LKAG本质上是δθ → 0极限下的计算，得到的是参考磁态附近的交换参数。若体系具有强非Heisenberg特性（磁矩大小随磁构型显著变化），则Jij可能不是固定不变的参数，此时需谨慎使用。
九、无序体系：KKR-CPA
对于无序体系（如非化学计量比的NiO），若采用超胞方法建立随机构型并逐一计算Jij，工作量巨大。KKR结合相干势近似（CPA）可有效处理此类问题。
9.1 CPA基本思想
以Ni0.8O为例，Ni占据概率cNi = 0.8，空位概率c□ = 0.2。CPA不建立具体的随机超胞，而是构建一个有效介质，代表80% Ni + 20% 空位的统计平均，在此有效介质中计算Jij。
9.2 有序与无序的比较
通过KKR-LKAG（有序体系）和KKR-CPA（无序体系）分别计算，可直接比较有序与无序状态下的交换相互作用差异，从而判断无序化究竟仅破坏磁长程有序，还是从电子结构层面改变了交换作用本身。
十、常用软件
程序 | KKR | Jij | CPA | SOC | |
|---|---|---|---|---|---|
SPR-KKR | ✓ | ✓ | ✓ | ✓ | |
JuKKR/KKRhost | ✓ | ✓ | ✓ | ✓ | |
KKRnano | ✓ | ✓ | ✓ | ✓ | |
AkaiKKR | ✓ | △ | ✓ | ✓ | |
Questaal | GF/LMTO | ✓ | ✓ | ✓ |
10.1 SPR-KKR
由Hubert Ebert团队发展的自旋极化相对论KKR程序，支持Green函数、自旋极化DFT、相对论/自旋轨道耦合、CPA无序、交换参数Jij、DMI、磁矩、磁各向异性等计算。官方手册包含Exchange coupling parameters (Jij)的完整章节及Fe示例，通过TASK JXC和CLURAD控制交换参数计算。
10.2 JuKKR
德国于利希研究中心发展的全电子KKR程序，支持ASA、全势、标量相对论、SOC、非共线磁性、CPA、Jij、DMI及并行计算。其KKRhost模块适合体相及合金CPA计算，KKRnano模块针对大规模体系（数千原子）的交换耦合计算。
10.3 AkaiKKR
经典的KKR-CPA程序，KKR与CPA结合自然，支持SCF、DOS、能带、合金、空位、杂质、磁性、SOC及CPA计算，适合入门KKR-CPA方法。
10.4 Questaal
包含LMTO、Green函数、ASA、全势、CPA及磁性计算功能，提供基于Liechtenstein公式的Heisenberg交换参数计算教程。
十一、总结
KKR-LKAG方法的本质是计算两个磁性原子之间电子传播的二阶磁能响应。该方法利用KKR的Green函数Gij，将电子从i到j的传播与两个原子的交换分裂Δi、Δj相结合，直接得到交换常数Jij。
核心公式：
对于NiO等磁性体系，通过SPR-KKR或JuKKR等程序，结合LKAG方法及CPA处理无序效应，可系统研究近邻及长程交换相互作用，为理解磁性基态提供定量依据。