第一性原理计算中的自旋极化：自旋分解、自旋密度与磁电子结构解析

说明：本文华算科技主要介绍自旋极化在第一性原理计算中如何把电子密度拆成上、下自旋分量，以及自旋密度、磁矩、能带和模型条件之间的相互指向。

自旋极化在DFT中改变的计算对象是啥

在自旋极化 DFT 里，Kohn-Sham 电子密度不再写成单一的 ρ(r)，而写成 ρ_↑(r) 与 ρ_↓(r) 两套空间分布。总电子密度来自 ρ_↑+ρ_↓，自旋密度来自 ρ_↑-ρ_↓。当两套分量处处相等，体系没有净自旋密度；当某些原子、缺陷或吸附物附近出现差值，计算会给出局域磁矩、交换分裂和自旋分辨能带。

图1. FeTAP、ZnTAP、RhTAP 和 AgTAP 的自旋密度以及 O₂ 吸附后的电荷密度差。DOI：10.1039/d6ra03135h

自旋极化把电子占据从一套能级扩展为两套可变占据，交换-相关能也按自旋分量处理。闭壳层分子、非磁金属或普通半导体常常收敛到 ρ_↑=ρ_↓；含未成对 d/f 电子、自由基、氧空位、吸附 O₂ 或过渡金属单原子位点时，上、下自旋占据不对称会进入总能计算，几何、电荷和吸附能随之改变。初始磁矩只是在自洽迭代开始时给电子一个分配方向，收敛后的 ρ_↑ 与 ρ_↓ 才是物理读数。

图1里的红、绿等值面对应两种自旋分量的差值区域。FeTAP 和 RhTAP 的自旋密度主要围绕金属中心和邻近配位原子分布，O₂ 吸附后等值面位置和形状发生变化，说明吸附物与金属 d 轨道、配体 p 轨道之间存在自旋相关杂化。自旋密度图读的是未成对电子的空间位置，不是普通总电荷密度的颜色深浅。

DFT 输出的总磁矩通常是整个计算胞内 ρ_↑-ρ_↓ 的积分，局域磁矩则依赖原子球、Bader 分区或投影方式。一个体系总磁矩接近整数 μ_B，并不代表每个原子都有同样大小的自旋；总磁矩接近零，也可能来自反平行排列的局域磁矩相互抵消。自旋极化的第一层读法应回到密度分量和空间分布。

在计算输入中，常见做法是给 Fe、Co、Ni、Mn、Cr 等原子设置初始磁矩，或给 O₂、NO、自由基吸附物指定自旋多重度。这个设置不等同于结论，它只是帮助自洽过程进入某个电子占据盆地。若不同初始磁矩收敛到不同总能和不同局域磁矩，体系存在多个可竞争自旋态，后续讨论应列出能量差而不是只保留最低能文件。

上下自旋通道为什么会形成不同能带？

自旋分量进入 Kohn-Sham 方程后，上、下自旋电子看到的有效势不完全相同。过渡金属 d 轨道附近，交换作用会把同一类 d 态分裂成两个自旋通道；晶体场再把 d 轨道分成 e_g、t_2g 或更低对称性的组合。DOS 中靠近 E_F 的上、下自旋峰位差，常用来识别未成对电子是否参与导电、吸附和磁交换。

图2. Mn、Fe、W 和 Co 掺杂 SnO 的结构以及自旋分辨 DOS 和 PDOS。DOI：10.1039/d6na00140h

Mn、Fe、W、Co 掺杂 SnO 的 DOS 面板把总态密度和 d 轨道投影分开画出。黑线给出总 DOS，彩色曲线给出 d_z2、d_xz、d_xy、d_x2-y2 等轨道贡献；上下半轴分别代表两套自旋通道。峰位不对称、峰强不对称和 E_F 附近的单通道态，都属于自旋极化在 DOS 中留下的直接信号。

图3. FeTAP、ZnTAP、RhTAP 和 AgTAP 在 O₂ 吸附前后的自旋分辨能带和 PDOS。DOI：10.1039/d6ra03135h

TMTAP 单原子位点吸附 O₂ 后，Fe 3d、Rh 4d、Ag 4d 与 O 2p 态在 E_F 附近重新分布。读此类图时，哪条曲线更高只是表面现象，吸附前后自旋通道的占据、带隙和轨道杂化是否同步变化更值得关注。若 O₂ 的 π* 轨道保留未成对电子，单线态、三线态或电荷转移态的排序会直接影响吸附能和反应路径。

半金属材料常出现一个自旋通道金属性、另一个自旋通道半导体性的图像；普通磁性半导体则可能上下通道都有带隙，只是带边位置不同。自旋极化升高或降低某一组 d 态时，应同时看该态的占据、相邻 p 态混合、E_F 位置和总磁矩，单独读一条 d 峰的位移很容易把轨道能级、杂化强度和费米能级重排合并成同一个判断。

吸附反应中还会遇到另一个细节：表面自旋态和分子自旋态可能同时改变。O₂、NO、过氧中间体、金属-氧自由基都可能把未成对电子带入反应坐标。若只比较非自旋极化吸附能，三线态 O₂ 到表面吸附态的自旋转变能会被漏掉。这类体系应把气相分子自旋、表面磁矩和吸附后自旋密度放在同一张能量表里。

自旋密度和局域磁矩记录哪类信号？

自旋密度图把 ρ_↑-ρ_↓ 放回实空间，适合观察未成对电子集中在金属中心、空位邻近原子、边缘态还是吸附分子上。局域磁矩则把这部分差值积分到某个原子或分区。二者分别对应“空间分布”和“积分数值”，一个给位置，一个给大小，读法不能互相替代。

图4. Fe 和 Mn 掺杂 Cs₂SnI₆ 的位点分辨磁矩。DOI：10.1039/d6ra02174c

Fe、Mn 掺杂 Cs₂SnI₆ 的位点磁矩图中，主要磁矩集中在 TM9 附近，I、Sn、Cs 位点贡献很小，胞间区域仍有可见贡献。这样的分布告诉我们，掺杂原子保留了强 d 电子局域特征，配体 p 态则通过杂化承担少量自旋极化。总磁矩相近的两个模型，局域磁矩分布仍可能差别很大。

缺陷体系尤其要把总磁矩、局域磁矩和自旋密度拆开读。氧空位可能把电子留在邻近金属 d 轨道，也可能形成更分散的缺陷态；边缘纳米带可能两条边反平行排列，总磁矩为零但每条边都有局域自旋；吸附自由基可能把未成对电子转移到表面。自旋极化的数值结果必须对应到原子位点或空间区域，否则很难判断活性中心来自金属、缺陷还是吸附物。

自旋密度等值面还受绘图阈值影响，同一份 CHGCAR 或 cube 文件用 0.002 e Å^-3 和 0.01 e Å^-3 画出来，等值面体积可能差很多。定量比较时，应固定等值面阈值、积分分区和超胞大小，再讨论不同缺陷、掺杂位点或吸附构型之间的自旋差异。

图5. Co 掺杂 SnO 在 U=0 eV 和 U=3 eV 下的自旋分辨能带及光学电导。DOI：10.1039/d6na00140h

Co 掺杂 SnO 的能带在 U=0 eV 和 U=3 eV 下出现变化，原因来自局域 d 电子的库仑排斥增强，而不是参数数值本身带来“更真实”的结果。U 会改变 d 态位置、带隙和自旋分裂。当磁矩主要来自局域 d/f 轨道时，泛函、U、SOC 和初始磁矩设置都会进入自旋态排序。

模型条件怎样改变自旋态排序？

自旋极化计算给出的能量差通常很小，几 meV 到几十 meV 就足以改变铁磁、反铁磁和非磁解的排序。含过渡金属卤化物、二维磁体或稀磁半导体时，磁序、晶格常数、U、SOC、k 点和超胞尺寸都会牵动交换参数。把一次收敛结果当成唯一自旋态，往往会忽略亚稳态。

图6. 单层 CrI₃ 中带隙和 Cr 局域磁矩随 U 参数变化。DOI：10.1021/acsomega.6c02634

单层 CrI₃ 的 DFT+U+SOC 结果中，U 从 0 eV 增加到 5 eV 时，带隙下降而 Cr 局域磁矩缓慢上升。这个趋势提醒我们，磁矩大小和电子结构间隙并不一定同向变化；局域电子更“像原子态”时，能带宽度、交换分裂和超交换路径会一起调整。

图7. 单层 CrI₃ 中 J₁、J₂、J₃ 交换参数随 U 参数变化。DOI：10.1021/acsomega.6c02634

同一篇 CrI₃ 计算中，J₁、J₂、J₃ 随 U 的曲线并不保持同一种变化方式。J₁ 在中等 U 区间增强后下降，J₃ 还会跨过零点。局域磁矩存在，只代表体系有自旋自由度；铁磁还是反铁磁、交换强弱和居里温度，还要看自旋之间的耦合路径。

做材料判断时，自旋极化至少要回答三件事：未成对电子分布在哪些原子或键区；上、下自旋 DOS 和能带在 E_F 附近如何分裂；不同磁序或初始磁矩收敛到哪些总能。若体系含重元素、非共线磁结构或强 SOC，标量自旋极化还不足以描述磁各向异性和自旋纹理。自旋极化的含义来自电子密度分解，判断则回到能量、图谱和模型条件的共同约束。

非共线磁计算把自旋从单一 z 方向扩展为三维矢量，自旋轨道耦合再把电子自旋和晶体方向联系起来。对 CrI₃、Fe₃GeTe₂、重元素掺杂二维材料或拓扑磁体，只看上、下自旋 DOS 不能覆盖磁各向异性、Dzyaloshinskii-Moriya 相互作用和自旋纹理。常规自旋极化结果适合回答未成对电子和磁矩问题，涉及磁易轴、反常霍尔或自旋输运时，SOC 与非共线设置要单独列入模型说明。