3. 磁性材料电子结构基础
在DFT自旋计算中，磁性的根源并不在“自旋”这个标签本身，而在电子结构。要理解为什么Fe、Co、Ni是铁磁金属，为什么MnO是反铁磁绝缘体，为什么Fe₃O₄同时存在两种磁性子晶格，必须回到能带、态密度和交换作用这三个基本概念上。本节从能带论出发，逐步建立从电子结构到宏观磁有序的完整逻辑链。
一、能带论：自旋极化的起点
能带论把晶体中电子的能量描述为波矢k的函数E(k)。在非磁性计算中，自旋向上和自旋向下的能带完全简并，即E↑(k)=E↓(k)。一旦体系存在净磁矩，自旋向上和自旋向下的电子感受到不同的交换势，两条能带发生相对位移，这就是交换劈裂（exchange splitting）。
交换劈裂的物理来源是库仑交换作用：自旋平行的电子因泡利不相容原理而彼此回避，库仑排斥能降低。在Stoner模型框架下，交换劈裂能可写为ΔE = I·m，其中I为Stoner交换积分，m为磁矩。这意味着磁矩越大，劈裂越显著。对Fe、Co、Ni，ΔE约为1～2 eV量级，正是这个劈裂使自旋向上和向下的d带发生错位，从而产生净自旋磁矩。
二、态密度：磁性的“指纹”
态密度DOS(E)描述单位能量区间内的电子态数目。在自旋极化计算中，我们分别得到DOS↑(E)和DOS↓(E)。两者的差异直接反映磁性的强弱：
- 铁磁体：DOS↑与DOS↓在费米面附近明显不对称，自旋向上和向下的d带占据数不同，净磁矩m = ∫(DOS↑−DOS↓)dE。
- 反铁磁体：两个子晶格的DOS↑和DOS↓互为镜像，总态密度在费米面附近对称，净磁矩为零，但局域磁矩不为零。
- 亚铁磁体：两个子晶格的磁矩大小不等或方向相反，总DOS不对称但净磁矩小于各子晶格磁矩之和。
在DFT后处理中，通常用dos.x或projwfc.x（Quantum ESPRESSO）或VASP的DOSCAR输出，观察费米能级附近自旋向上/向下d态的占据差异。例如bcc Fe的典型结果：自旋向上d带几乎全满，自旋向下d带部分占据，费米面处DOS↓远大于DOS↑，净磁矩约2.2 μB/atom。
三、Stoner模型与巡游磁性判据
Stoner模型给出了巡游电子系统出现铁磁性的判据：I·N(E_F) > 1，其中N(E_F)为费米面处非磁性态密度。当交换积分I与费米面态密度的乘积超过1时，顺磁态失稳，系统自发极化。这一判据解释了为什么Fe、Co、Ni是铁磁金属：它们的d带在费米面处DOS很高，且I足够大。
Stoner模型还预言了磁矩与劈裂的关系：m = ΔE/I。在计算中，可以通过固定磁矩做约束DFT，扫描不同m下的总能量，找到能量最低点对应的平衡磁矩。这一方法在Fe、Co、Ni及其合金的磁矩预测中非常实用。
四、巡游磁性与局域磁性
巡游磁性（itinerant magnetism）指磁矩由扩展的能带电子共同贡献，典型代表是3d过渡金属Fe、Co、Ni。其特点是磁矩非整数、对结构敏感、存在自旋涨落。局域磁性（localized magnetism）指磁矩定域在特定原子或离子周围，典型代表是4f稀土和3d氧化物中的金属离子，其磁矩接近整数、由局域d/f轨道占据决定。
在DFT计算中，区分两者的实用判据是：观察磁矩是否主要来自某个原子球内的d/f轨道投影，以及磁矩对晶格常数、体积的敏感程度。巡游磁体的磁矩随体积变化明显，局域磁体的磁矩则相对稳定。
五、铁磁、反铁磁与亚铁磁的电子结构特征
三种磁有序在电子结构上的差异，可以通过自旋分辨的能带和DOS来识别：
- 铁磁（FM）：自旋向上和向下能带整体劈裂，费米面处自旋极化率P = (N↑−N↓)/(N↑+N↓)不为零。典型如bcc Fe、fcc Co、fcc Ni。
- 反铁磁（AFM）：两个子晶格的自旋方向相反，能带在布里渊区中发生折叠，总DOS在费米面附近对称。典型如MnO、NiO，其带隙约2～4 eV，属于电荷转移型绝缘体。
- 亚铁磁（FiM）：两个子晶格磁矩不等，总DOS不对称但净磁矩小于各子晶格之和。典型如Fe₃O₄，其中四面体A位Fe³⁺与八面体B位Fe³⁺/Fe²⁺的磁矩反平行，净磁矩约4 μB/f.u.。
六、从电子结构到磁有序的判断流程
在实际计算中，判断一个材料属于哪类磁性，可以按以下步骤进行：先做非磁性计算得到N(E_F)，用Stoner判据I·N(E_F)>1初步判断是否可能铁磁；再分别设置FM、AFM、FiM初始磁矩做自旋极化计算，比较总能量；最后分析自旋分辨DOS和原子投影磁矩，确认磁有序类型和磁矩大小。总能量最低的磁态即为基态磁结构。
上图给出了从非磁性计算到磁有序类型判断的完整逻辑链。需要强调的是，Stoner判据只是初步筛选，最终判断必须依赖自旋极化总能量比较和DOS分析。对于强关联体系如MnO、NiO，还需引入DFT+U或杂化泛函才能正确描述带隙和磁矩。
七、计算实践中的关键参数与判断依据
在VASP或Quantum ESPRESSO中进行自旋极化计算时，以下参数直接影响磁性结果：
| 参数/设置 | 作用 | 典型值/建议 |
|---|---|---|
| ISPIN=2（VASP） | 开启自旋极化 | 必须设置，否则无法得到磁矩 |
| MAGMOM | 初始磁矩 | Fe: 4～5 μB；AFM用正负号区分 |
| k点密度 | 影响费米面附近DOS精度 | 金属体系≥9×9×9 |
| 展宽方法 | 影响DOS和能量收敛 | 金属用Methfessel-Paxton，绝缘体用Gaussian |
| DFT+U | 修正强关联d/f电子 | MnO: U_eff≈4～5 eV |
判断磁态是否合理，可以检查三个量：原子投影磁矩是否接近实验值或化学直觉；自旋分辨DOS在费米面处的极化方向是否与预期一致；总能量是否低于非磁性态。若磁矩远小于预期或DOS完全对称，通常说明初始磁矩设置不当或计算未收敛到磁性基态。
总结来说，能带论提供了自旋极化的框架，态密度给出了磁性的定量指纹，Stoner模型给出了巡游铁磁的判据，而铁磁、反铁磁、亚铁磁的区分最终落在自旋分辨DOS的对称性和原子投影磁矩的分布上。掌握这条从电子结构到磁有序的分析链，是后续自旋计算和磁各向异性研究的基础。