KKR方法计算磁性体系交换常数

 

KKR方法计算磁性体系交换常数

一、引言

KKR（Korringa-Kohn-Rostoker）方法计算磁性体系交换常数，通常采用Liechtenstein–Katsnelson–Antropov–Gubanov（LKAG）公式。该方法通过自旋向上/向下电子的Green函数，计算体系总能量对两个磁性原子磁矩相对旋转的二阶响应，从而得到交换常数Jij。

与基于总能量的方法不同，LKAG不需要构建多个磁构型进行能量差计算，而是直接从电子Green函数导出Jij。

二、Heisenberg模型与交换常数定义

磁性体系通常映射为Heisenberg模型：

Heisenberg Hamiltonian

其中i, j表示两个磁性原子，ei、ej为局域磁矩方向的单位矢量，Jij为交换常数。

符号约定：

J positiveJ negative

注：不同程序的Hamiltonian符号约定可能存在差异，实际分析时须参照具体软件定义。

三、LKAG方法核心思想

LKAG方法的关键在于考虑两个局域磁矩发生无穷小旋转：

ei → ei + δei

ej → ej + δej

体系总能量展开为：

Energy expansion

由此得到：

J derivative

交换常数本质上即为总能量对两个局域磁矩相对旋转的二阶导数。KKR方法的优势在于，该二阶导数可用Green函数直接表达。

四、LKAG公式

经典的LKAG公式表达形式为：

LKAG formula

4.1 交换分裂 Δi

Δi表示第i个原子的交换分裂：

Delta i

其中ti↑(E)和ti↓(E)分别为自旋向上和向下的单原子t矩阵。直观上：

Exchange splitting

4.2 Green函数 Gij

Gij(E)表示从原子i传播到原子j的Green函数：

G_ij

LKAG公式计算的物理过程为：磁性原子i → 电子传播 → 原子j → 再传播回i，这一过程对磁矩耦合的贡献。

五、Green函数与交换作用的物理图像

设体系中有两个磁性原子i和j。以3d电子为例，其具有交换分裂，即自旋向上与向下的电子散射势不同。电子从i传播到j（Gij），受到j的磁性散射（Δj），再传播回来（Gji），产生耦合项：

J approximation

该量包含了两个磁矩之间的电子介导耦合。因此，KKR-LKAG特别适合研究以下机制：

  • RKKY相互作用
  • 超交换作用
  • 双交换作用
  • 巡游磁性
  • 长程交换作用

上述机制本质上均与电子传播相关。

六、KKR方法的优势

KKR基于多重散射理论（Multiple Scattering Theory），将晶体视为电子在各原子间不断散射的过程。KKR直接构造实空间Green函数G(r, r′, E)，天然拥有电子从某一原子传播到另一原子的信息。

由于Jij恰好是这种信息的二阶响应，LKAG与KKR的结合具有自然的理论自洽性。

七、实际计算流程

以NiO体系为例，完整计算流程如下：

  1. 建立晶体结构
  2. KKR自洽计算（spin-polarized）
  3. 得到电荷密度与有效势
  4. 得到自旋极化电子结构
  5. 计算单原子t矩阵
  6. 得到散射路径算符τij
  7. 构造Green函数Gij
  8. LKAG积分计算Jij
  9. 按距离/原子对分类整理
  10. 得到J1, J2, J3, ...

7.1 参考磁性状态

计算前需指定参考磁态，例如：

  • 铁磁（FM）：所有磁矩同向排列
  • 反铁磁（AFM）：磁矩交替反向排列

7.2 自洽计算

进行自旋极化KKR计算，分别得到有效势Veff↑(r)和Veff↓(r)，进而得到各原子的局域磁矩：

Magnetic moment

7.3 散射矩阵与Green函数

每个原子对应单原子散射矩阵ti(E)。对于磁性原子，ti↑(E) ≠ ti↓(E)。通过散射路径算符τij(E)：

tau_ij

τij包含电子从i到j的所有多重散射路径，因此可描述复杂的电子传播过程，这也是KKR能够得到长程Jij的原因。

7.4 LKAG积分

J integral

实际程序采用复能量积分或围道积分等数值方法，而非直接从−∞积到EF。该步骤为计算量较大的环节。

7.5 结果分析

最终得到不同近邻距离的交换常数，例如：

J1 = 8.2 meV（最近邻）
J2 = −1.4 meV（次近邻）
J3 = 0.37 meV（第三近邻）

可进一步绘制Jij随距离rij的变化关系：

J vs r

该图对于判断磁性机制具有重要价值。

八、KKR-LKAG与VASP四态法的对比

8.1 VASP四态法

人为建立不同磁构型（如↑↑和↑↓），通过总能量差计算交换常数：

VASP J

该方法需要计算多个有限磁构型的总能量，通过能量差拟合得到J。

8.2 KKR-LKAG

仅需一个参考磁态，通过Green函数直接计算无穷小旋转下的响应：

J proportional

无需构建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（无序体系）分别计算，可直接比较有序与无序状态下的交换相互作用差异，从而判断无序化究竟仅破坏磁长程有序，还是从电子结构层面改变了交换作用本身。

十、常用软件

程序KKRJijCPASOC 
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。

核心公式：

J proportional

对于NiO等磁性体系，通过SPR-KKR或JuKKR等程序，结合LKAG方法及CPA处理无序效应，可系统研究近邻及长程交换相互作用，为理解磁性基态提供定量依据。