KKR方法计算磁性体系交换常数KKR方法计算磁性体系交换常数一、引言KKRKorringa-Kohn-Rostoker方法计算磁性体系交换常数通常采用Liechtenstein–Katsnelson–Antropov–GubanovLKAG公式。该方法通过自旋向上/向下电子的Green函数计算体系总能量对两个磁性原子磁矩相对旋转的二阶响应从而得到交换常数Jij。与基于总能量的方法不同LKAG不需要构建多个磁构型进行能量差计算而是直接从电子Green函数导出Jij。二、Heisenberg模型与交换常数定义磁性体系通常映射为Heisenberg模型其中i, j表示两个磁性原子ei、ej为局域磁矩方向的单位矢量Jij为交换常数。符号约定注不同程序的Hamiltonian符号约定可能存在差异实际分析时须参照具体软件定义。三、LKAG方法核心思想LKAG方法的关键在于考虑两个局域磁矩发生无穷小旋转ei →ei δeiej →ej δej体系总能量展开为由此得到交换常数本质上即为总能量对两个局域磁矩相对旋转的二阶导数。KKR方法的优势在于该二阶导数可用Green函数直接表达。四、LKAG公式经典的LKAG公式表达形式为4.1 交换分裂 ΔiΔi表示第i个原子的交换分裂其中ti↑(E)和ti↓(E)分别为自旋向上和向下的单原子t矩阵。直观上4.2 Green函数 GijGij(E)表示从原子i传播到原子j的Green函数LKAG公式计算的物理过程为磁性原子i → 电子传播 → 原子j → 再传播回i这一过程对磁矩耦合的贡献。五、Green函数与交换作用的物理图像设体系中有两个磁性原子i和j。以3d电子为例其具有交换分裂即自旋向上与向下的电子散射势不同。电子从i传播到jGij受到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函数GijLKAG积分计算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无序体系分别计算可直接比较有序与无序状态下的交换相互作用差异从而判断无序化究竟仅破坏磁长程有序还是从电子结构层面改变了交换作用本身。十、常用软件程序KKRJijCPASOCSPR-KKR✓✓✓✓JuKKR/KKRhost✓✓✓✓KKRnano✓✓✓✓AkaiKKR✓△✓✓QuestaalGF/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处理无序效应可系统研究近邻及长程交换相互作用为理解磁性基态提供定量依据。 KKR方法计算磁性体系交换常数
KKR方法计算磁性体系交换常数KKR方法计算磁性体系交换常数一、引言KKRKorringa-Kohn-Rostoker方法计算磁性体系交换常数通常采用Liechtenstein–Katsnelson–Antropov–GubanovLKAG公式。该方法通过自旋向上/向下电子的Green函数计算体系总能量对两个磁性原子磁矩相对旋转的二阶响应从而得到交换常数Jij。与基于总能量的方法不同LKAG不需要构建多个磁构型进行能量差计算而是直接从电子Green函数导出Jij。二、Heisenberg模型与交换常数定义磁性体系通常映射为Heisenberg模型其中i, j表示两个磁性原子ei、ej为局域磁矩方向的单位矢量Jij为交换常数。符号约定注不同程序的Hamiltonian符号约定可能存在差异实际分析时须参照具体软件定义。三、LKAG方法核心思想LKAG方法的关键在于考虑两个局域磁矩发生无穷小旋转ei →ei δeiej →ej δej体系总能量展开为由此得到交换常数本质上即为总能量对两个局域磁矩相对旋转的二阶导数。KKR方法的优势在于该二阶导数可用Green函数直接表达。四、LKAG公式经典的LKAG公式表达形式为4.1 交换分裂 ΔiΔi表示第i个原子的交换分裂其中ti↑(E)和ti↓(E)分别为自旋向上和向下的单原子t矩阵。直观上4.2 Green函数 GijGij(E)表示从原子i传播到原子j的Green函数LKAG公式计算的物理过程为磁性原子i → 电子传播 → 原子j → 再传播回i这一过程对磁矩耦合的贡献。五、Green函数与交换作用的物理图像设体系中有两个磁性原子i和j。以3d电子为例其具有交换分裂即自旋向上与向下的电子散射势不同。电子从i传播到jGij受到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函数GijLKAG积分计算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无序体系分别计算可直接比较有序与无序状态下的交换相互作用差异从而判断无序化究竟仅破坏磁长程有序还是从电子结构层面改变了交换作用本身。十、常用软件程序KKRJijCPASOCSPR-KKR✓✓✓✓JuKKR/KKRhost✓✓✓✓KKRnano✓✓✓✓AkaiKKR✓△✓✓QuestaalGF/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处理无序效应可系统研究近邻及长程交换相互作用为理解磁性基态提供定量依据。