原子结构计算电子数：DFT电子占据与PAW势的实战解析

一个看起来非常简单的问题：单个Fe原子的基态电子组态是什么？化学教科书上写的是[Ar]3d⁶4s²。但用VASP做一个孤立Fe原子的DFT计算，结果可能是3d⁷4s¹，也可能是3d⁶.³4s¹.⁷——取决于你用的赝势、展宽方式和自旋处理。这个差异让项目组意识到，原子结构计算电子数这件事，从理论到计算实践之间隔着好几层需要仔细处理的细节。

电子占据：smearing不是万能的

DFT计算中处理金属性体系的电子占据，标准做法是用Methfessel-Paxton或Gaussian smearing来平滑Fermi面附近的占据数跃变。对于块体Fe这种金属体系，ISMEAR=1（Methfessel-Paxton一阶）配合SIGMA=0.20eV是常规选择，能稳定收敛。

但对孤立原子来说，问题完全不一样。孤立原子的能级是分立的，不存在连续的能带结构。用smearing来处理孤立原子，本质上是在人为地展宽分立的能级，导致分数占据——原本应该整数占据的轨道被”涂抹”成了非整数占据。项目组对孤立Fe原子用ISMEAR=0（Gaussian）配合SIGMA=0.01eV做了测试，得到了3d⁶.⁴4s¹.⁶的电子占据，这显然不是物理上正确的整数占据结果。

解决方案是用ISMEAR=-2（Fermi-Dirac smearing）配合一个极小的SIGMA值（0.001eV），同时关掉对称性（ISYM=0）以避免VASP自动对称化导致能级简并。在这种设置下，孤立Fe原子的电子占据收敛到了3d⁶4s²，与教科书一致。但这个设置在块体计算中是不推荐的——SIGMA太小会导致收敛极其困难。

原子结构计算电子数的一个核心教训是：块体计算和孤立原子/分子计算的参数策略完全不同。把块体计算的那套smearing参数直接搬到原子计算上，电子占据结果大概率是错的。

PAW势的”隐藏”电子

PAW势（Projector Augmented Wave）是VASP中处理芯电子的标准方法。但PAW势有一个容易被忽视的特性：势文件中定义的价电子数不一定等于你在OUTCAR中看到的电子数。

以Fe的PAW_PBE赝势为例，势文件中定义了8个价电子（3d⁶4s²），但VASP实际计算的是Kohn-Sham轨道中的电子密度，而PAW重构会在芯区引入额外的电荷补偿。这意味着你用Bader电荷分析得到的Fe原子上的电子数，不是”价电子数”，而是包含了PAW重构贡献的有效电荷。

项目组做了个对比：对α-Fe₂O₃做Bader电荷分析，Fe原子上的Bader电荷约为+1.7e，O原子约为-1.1e。如果简单地用8（价电子数）减去1.7（Bader电荷）得到6.3作为”Fe的电子数”，这个数字在物理上毫无意义——Bader电荷反映的是电荷转移趋势，不是绝对电子数。

正确的解读方式是：Bader电荷的绝对值意义不大，但相对值有明确的物理含义。Fe₂O₃中两个Fe原子的Bader电荷如果有差异（比如+1.65和+1.75），说明两个Fe位的化学环境不同，可能存在不同的氧化态或配位环境。用Bader电荷做绝对电子数计算是方向性的错误。

态密度积分：看起来简单，做起来容易出错

另一个常被忽视的问题是：从态密度积分得到电子数时，积分上下限和占据判据的选择。项目组在分析Fe掺杂ZnO的电子结构时，需要计算Fe 3d轨道的电子占据数。

操作步骤本身很简单：从VASP的DOSCAR文件中提取Fe的3d投影态密度（通过LORBIT=11计算获得），然后从能量下界积分到费米能级。但这里有两个容易出错的地方：

第一，能量积分下限的选择。如果积分下限设得太高（比如-5eV），会漏掉深层的成键态，导致电子数偏低。项目组建议积分下限至少取到-10eV，确保覆盖所有3d轨道的成键态和反键态中低于费米能级的部分。

第二，投影态密度的归一化问题。LORBIT=11给出的投影态密度是对球内电荷的投影，而球半径（RWIGS）的选择直接影响投影结果。如果RWIGS设得太小，一部分3d电荷会落在投影球外，导致投影态密度偏低，积分得到的电子数自然也偏低。项目组对Fe掺杂ZnO体系的测试表明，当RWIGS从1.3Å增加到1.5Å时，Fe的3d投影电子数从5.8增加到了6.3——差了将近半个电子。

对于需要精确电子占据数的场景（比如确定过渡金属离子的氧化态），不建议依赖投影态密度积分，而应该结合Bader电荷分析和磁矩信息做交叉验证。磁矩对3d占据数非常敏感——Fe³⁺（3d⁵）的高自旋磁矩约为5μB，Fe²⁺（3d⁶）约为4μB。如果Bader电荷显示+1.7但磁矩只有3.2μB，那就需要回头检查计算参数是否合理。

自旋极化：别忘了这个开关

原子结构计算电子数的最后一个常见错误是忘记开自旋极化。对于闭壳层体系（比如NaCl），不开自旋极化影响不大。但对于含3d过渡金属的体系，不开自旋极化（ISPIN=1）意味着强制所有轨道双占据，这会完全改变电子组态。

以孤立Mn原子为例：不开自旋极化时，VASP会把3d轨道按能量从低到高双占据填充，得到一个完全不反映Hund规则的结果。开启自旋极化（ISPIN=2）后，Mn原子的5个3d电子按Hund规则分别占据5个自旋向上的3d轨道，总磁矩5μB，这才是物理正确的基态。

项目组的经验是：只要体系中含d区或f区元素，默认就应该开启自旋极化。即使最终结果显示磁矩很小甚至为零，也应该先开ISPIN=2跑一遍，确认体系确实没有磁矩后再考虑关掉。很多过渡金属氧化物在非磁计算中看起来是导体，开了自旋极化后出现带隙变成半导体——这个差异足以改变对材料电子结构的全部判断。

DFT原子结构计算电子数这件事，拆开来看每个环节都不复杂：选赝势、设smearing、开自旋、做Bader分析。但把这些环节串起来时，任何一个步骤的参数选择偏差都会在最终结果中被放大。电子占据数的精确性取决于你是否理解了每一步近似背后的物理假设——smearing是对温度效应的模拟而非对分立能级的展宽，Bader电荷是相对量而非绝对量，投影态密度依赖于球半径这个人为参数。

回过头看，孤立Fe原子3d⁶4s²这个结果之所以可贵，不是因为它”正确”，而是因为它经过了从赝势选择到占据判据的每一步推敲。

更多DFT计算技术讨论，欢迎访问科研学术网第一性原理栏目。计算资源与技术支持请见科研学术网首页。电子结构的精确计算从来不是一键求解，而是参数选择背后的物理直觉与系统验证。