说明：本文华算科技主要介绍自旋极化在第一性原理计算中的物理对象、计算来源、图谱读法和材料判断范围。

自旋极化的计算对象是什么？

在普通的非自旋计算里，体系只需要一套电子密度来描述占据；自旋极化计算把这套密度分解为两个自旋分量：上自旋电子密度 ρ↑(r) 和 下自旋电子密度 ρ↓(r)。两者相加给出总电子密度 ρ(r)，两者相减给出 自旋密度 m(r)。局域磁矩来自某个原子区域、投影球或轨道投影中的自旋电子数差值，常写成 μB 为单位的磁矩。

这种拆分直接处理电子占据自由度：同一个 3d 轨道、缺陷态或吸附分子轨道，可以在两个自旋分量中有不同填充。过渡金属中心、氧空位附近的金属原子、吸附的 O2、NO、自由基中间体，本来就可能带未成对电子；若强行让两种自旋占据相等，优化得到的键长、电荷分布和总能量都可能偏离真实低能态。

图1. 双层 silicene 的吸附位、层间自旋方向和吸附体系视图，展示模型中原子位置与自旋排布的关系。 DOI：10.1039/D2RA06955E

计算程序需要一个初始自旋排布来开始迭代，例如给不同原子设定初始磁矩，或把相邻层、相邻金属位点设置成平行和反平行。这个初值只负责把体系带到某个自旋态附近，不能直接当作结果。真正要读的是收敛后的总磁矩、各原子局域磁矩、自旋分辨态密度、自旋分辨能带、自旋密度分布，以及不同磁序之间的总能量差。

判断自旋极化时，总磁矩为零仍可能对应局域磁性。反铁磁材料里，相邻原子的局域磁矩方向相反，积分后的总磁矩可以抵消；若只看一个总数，就会把局域自旋排布漏掉。对缺陷、单原子位点和低维材料来说，局域磁矩的空间位置比总磁矩更能说明电子到底停在哪些原子或键附近。

自旋极化如何改变电子图谱？

自旋极化进入能带和态密度图后，最明显的变化是 交换分裂。原来重合的一组能级，在上、下自旋分量中会发生相对移动；靠近费米能级 EF 的峰、带边或平带，可能只出现在其中一个自旋分量里。若 一个自旋分量穿过 EF，另一个自旋分量仍保持带隙，材料就呈现半金属特征，这类图谱在自旋电子学里很常见。

态密度图常把上自旋画在横轴上方、下自旋画在横轴下方，这只是绘图约定。需要逐项比较 峰位、峰宽、投影轨道来源 和 EF 附近的占据。某个峰只在一个自旋分量增强，通常意味着这个原子或轨道参与了未成对电子占据；若两个自旋分量形状几乎镜像相同，自旋极化对该区域电子结构的影响就弱得多。

图2. 不同非磁吸附原子修饰双层 silicene 后的自旋分辨能带，红色和蓝色曲线对应相反自旋分量。 DOI：10.1039/D2RA06955E

能带给出动量空间中的色散关系，态密度把所有 k 点和能级贡献压缩成能量分布。两张图合读时，可以从 EF 附近的 单个自旋分量穿越 入手，随后用 投影态密度 确认这些态来自吸附原子、近邻 Si 原子，还是远离吸附位的骨架原子。若活性位的局域态靠近 EF，它会影响电子给受能力，也会改变吸附中间体的轨道杂化方式。

图3. Mg 吸附体系和本征双层 silicene 的分层 DOS 对比，显示不同层对自旋分辨态密度的贡献。 DOI：10.1039/D2RA06955E

在金属表面或含 d 电子催化位点，EF 附近的单一自旋态还会改变吸附中间体的成键和反键占据。*O、*OH 与金属 d 轨道杂化时，上、下自旋分量的反键态填充不同，吸附能、键长和振动频率都会随磁序发生偏移，这类变化在磁性氧化物表面尤其明显。

在吸附体系里，局域态常集中在 吸附原子 和 最近邻原子 上，离吸附位越远，态密度形状越接近基底本征特征。这个距离依赖能帮助区分两类情况：一种是吸附原子把未成对电子局限在很小区域，另一种是吸附诱导的自旋极化沿二维骨架扩展。前者更像局域磁中心，后者会改写材料的输运和磁耦合。

图4. Mg 吸附体系中不同距离 Si 原子的平均 DOS，对比最近邻、次近邻和远离吸附位原子的贡献。 DOI：10.1039/D2RA06955E

磁矩与自旋密度有何差异？

磁矩是一个积分后的数，自旋密度 是空间中的分布。把 ρ↑(r)-ρ↓(r) 画成等值面，可以看到哪一侧、哪一层、哪个配位环境积累了净自旋。对于过渡金属氧化物、硫化物和二维材料，d 轨道与邻近 p 轨道的杂化 会让自旋密度从金属中心延伸到配位原子上；这时只报金属原子的磁矩，信息是不够的。

自旋密度图和 差分电荷密度图 看的是两个变量。差分电荷密度强调成键或吸附前后的电荷再分布，自旋密度强调两个自旋分量的占据差。一个区域可以有明显电荷转移，却没有显著净自旋；也可以电荷转移很小，但未成对电子集中在某个 p 轨道或 d 轨道上。

图5. Mg 吸附双层 silicene 的差分电荷密度与自旋密度，分别展示电荷重排和净自旋分布。 DOI：10.1039/D2RA06955E

同一个磁矩数值换成 Bader 分区、Mulliken 布居、PAW 投影球 或轨道投影后，数值口径会发生变化，绝对值不宜横向混用；在同一套方法、同一类模型中比较趋势更有意义。例如同一个缺陷位点在不同电荷态下磁矩从 1 μB 变到 0 μB，这说明未成对电子被填充或抽空；若只是不同软件给出 0.82 和 0.95 μB，不能据此判断磁性强弱发生本质变化。

铁磁、反铁磁和亚铁磁构型 不能只靠一个初始态来代表。相同化学组成下，平行排列、反平行排列和层间反平行排列可能都能收敛，各自对应不同的局域磁矩、带隙和总能量。若能量差只有几 meV，k 点、U 值、SOC、结构畸变和温度扰动都会影响排序；若能量差达到几十到几百 meV，磁序判断才更容易在计算设置变化下保持稳定。

图6. 双层 silicene 在 100 K 和 300 K 自旋极化 AIMD 中的温度与磁矩波动。 DOI：10.1039/D2RA06955E

自旋极化如何影响能量计算？

自旋极化会改变 总能量排序，所以它会影响形成能、吸附能、缺陷形成能、反应能垒和自由能修正。以 O2、NO、*OOH、过渡金属位点和氧空位为例，参考态、初态、中间体和终态 若采用不兼容的自旋处理，能量差会混入人为误差。催化计算里常见的

把 PBE、DFT+U、HSE06 和 SOC 放到同一个体系里比较，常会看到同一组 d 态在 EF 附近移动。PBE 可能低估局域电子相关，DFT+U 会把部分 d 态推离 EF，HSE06 常拉大带隙，SOC 会把自旋方向和轨道角动量耦合起来。低维磁性材料还受真空层、层间距、覆盖度和电荷态影响。读自旋极化结果时，需要把自旋分量和这些模型条件放在同一张表里，而不是只截取某一张漂亮的 DOS 图。

图7. 反铁磁双层 silicene 在 PBE+D2 与 HSE06 条件下的能带对比，显示方法选择对自旋分辨带结构的影响。 DOI：10.1039/D2RA06955E

AIMD 轨迹 里，键长和层间距随时间抖动，局域磁矩也跟着涨落。磁矩随时间变化，反映的是原子振动、键长波动和电子占据共同作用后的响应；它和 0 K 优化结构上的单点磁矩不是同一个问题。若材料应用场景在室温、外场或界面环境中，静态磁矩只能给出参考态，自旋有序能否保留，需要检查结构扰动后的图谱。

图8. 自旋极化 AIMD 后双层 silicene 在 100 K 和 300 K 下的侧视与正视结构。 DOI：10.1039/D2RA06955E

一份有用的自旋极化计算结果，至少能回答四个具体问题：哪几个原子带 局域磁矩，哪个自旋分量在 EF 附近贡献主要态，不同磁序的总能量差有多大，结构变量改变时自旋密度如何移动。把这些结果连起来，才能判断材料的本征活性、缺陷态、界面电场调控或磁输运是否真的和自旋占据有关。

读图顺序可以从模型里的未成对电子来源开始，随后检查自旋分辨能带和 DOS，再用自旋密度定位空间分布，最后比较磁序和能量差。