铁的电子结构计算

我做过一个纯铁的电子结构计算项目，客户要解释为什么 bcc-Fe 是铁磁性的，以及自旋向上和自旋向下的电子态密度为什么不对称。这个案例很经典，但要把磁矩、能带和 DOS 一起讲清楚，还是需要老老实实做一遍自旋极化计算。

我算的体系是体心立方铁，晶格常数取实验值 2.87 Å。铁磁态的关键是 ISPIN = 2，也就是把电子按自旋向上（majority）和自旋向下（minority）分开处理。左侧图是自旋分辨能带。蓝色实线是 majority spin，红色虚线是 minority spin。可以看到在 Γ-H 和 H-N 方向，两条能带明显不重合，尤其是在费米面 E_F 附近，minority 的 d 带更窄、电子数更少。正是这个不对称导致净自旋，从而出现宏观磁矩。

右侧图是自旋分辨态密度（DOS），标注的平均磁矩 m = 2.2 μB/atom。蓝色填充区域是 spin-up DOS，红色是 spin-down。DOS 曲线有几个特征：在 -2 eV 附近 spin-up 有一个强峰，对应成键性 d 电子；-4 eV 以下有 s/p 电子贡献；费米面处 spin-up DOS 明显高于 spin-down，说明自旋向上的导电电子更多。磁矩可以通过积分费米面以下的 DOS 差得到，m = ∫[n↑(E) − n↓(E)]dE，我算出来 2.18 μB，和实验 2.2 μB 很接近。

泛函选择上，PBE 对铁的磁矩和晶格常数已经不错，但要更准的态密度细节可以试 PBEsol。自旋极化计算对初始磁矩很敏感，我在 POSCAR 之后用 MAGMOM 给每个 Fe 原子初始磁矩 2.5 μB，反铁磁初态我也试了一下，结果能量比铁磁态高 0.3 eV/atom，说明基态确实是铁磁。

我这个项目有几个容易踩的坑。第一个是 K 点：铁是金属，能带穿过费米面，K 点密度不能低。我用的 16×16×16 Monkhorst-Pack，否则 DOS 的尖峰会很糙。第二个是 smearing：金属体系我用 Methfessel-Paxton 二阶，SIGMA=0.1 eV，配合 ISMEAR=-5 做精确 DOS 计算。第三个是非自洽 DOS：先自洽优化，再用大 K 点算非自洽 DOS，才能保证 DOS 图平滑。第四个是磁各向异性：我这个项目没算磁晶各向异性能，如果要考虑，需要加 LSORBIT 做自旋轨道耦合，机时会翻几倍。

经验总结：铁的电子结构计算最重要的是自旋极化和 K 点收敛。磁矩 2.2 μB 是一个很强的验证指标，如果算出来偏离超过 0.3 μB，通常是泛函或初始磁矩有问题。这个案例也适合拿给学生教学：能带和 DOS 为什么不对称，本质上是因为交换作用让同自旋电子感受到不同的有效势。

更多 VASP 自旋极化与磁性材料的计算细节，我放在 [铁的电子结构计算](https://www.keyanxueshu.com/category/dft/) 栏目里。[铁的电子结构计算](https://www.keyanxueshu.com/)

从教学角度看，左侧能带图有一个很典型的特征：majority 和 minority 能带在费米面附近的交错。spin-up 的 d 带穿过 E_F，而 spin-down 的 d 带在 E_F 附近有一个赝能隙。这种不对称正是铁磁性的根源——自旋向上电子数多于自旋向下。我们算出来的 m=2.2 μB 和实验值吻合，说明 ISPIN=2 加合理初始磁矩就够用了。但有一个细节：PBE 会略微高估铁的晶格常数和磁矩，如果要做精确对比，建议用 PBEsol 或 SCAN。另外，态密度图里 −2 eV 处的 majority 峰对应 t2g 电子，−1 eV 以下对应 eg 电子，做轨道投影 DOS 时能把这些归属讲得更清楚。