第一性原理计算中自旋轨道耦合和自旋极化分别改变哪些物理量

说明：本文华算科技主要介绍自旋轨道耦合与自旋极化在第一性原理计算中分别改变的物理量、四组对照计算的归因方法，以及能带、态密度、自旋纹理和磁各向异性中可区分的信号。

自旋极化描述电子密度对自旋自由度的响应。共线计算把总电子密度写成 n(r)=n_↑(r)+n_↓(r)，并用 自旋密度 m(r)=n↑(r)−n↓(r) 表示两类分量的空间差异。若 m(r) 在某些原子附近保持非零，体系含有局域磁矩；对整个晶胞积分为零时，局域磁矩仍可按反铁磁次序互相抵消。

允许 n_↑ 与 n_↓ 分别自洽之后，交换关联势出现自旋依赖项 V_xc^↑ 与 V_xc^↓。自洽交换场 会改变轨道占据、局域磁矩、总能和能带位置。所谓“开启自旋极化”，对应的是放开这组密度自由度，并未预设计算一定收敛到铁磁基态。

反铁磁晶胞常满足总磁矩 M=∫m(r)dr=0，但不同磁性子晶格上的 m(r) 符号相反。零净磁矩与零局域自旋密度属于两个不同陈述，判断前者看晶胞积分，判断后者看原子投影磁矩或实空间自旋密度。

自旋轨道耦合来自电子自旋与轨道运动之间的相对论相互作用。Pauli 形式可写为 H_SOC∝σ·[∇V(r)×p]，原子附近常写成 ξ(r)L·S。SOC 把自旋、轨道与晶格方向关联起来，重元素较大的核势梯度通常带来更强的 ξ，但实际效应还受轨道成分、能级间隔和晶体对称性控制。

引入 H_SOC 后，Kohn–Sham 波函数通常写成双分量旋量。此时单条能带可同时含有 ↑ 与 ↓ 分量，沿选定轴的自旋期望值 ⟨σ_x⟩、⟨σ_y⟩、⟨σ_z⟩ 比纯粹的“上自旋带”“下自旋带”更贴合计算对象。总磁矩仍可为零，非中心对称非磁材料却能在 k 与 −k 处出现相反的自旋纹理。

自旋极化控制磁性密度是否参与自洽，SOC 控制自旋是否与轨道和空间方向耦合。两项设置对应两条可独立改变的计算轴；磁性来源、相对论分裂和方向相关能量只有经过独立对照，才能分别归入交换作用、SOC 或二者的交叉项。

同一结构可形成四种计算状态：A 为非自旋极化且无 SOC，B 为自旋极化且无 SOC，C 为保持非磁约束并加入 SOC，D 为自旋极化并加入 SOC。A→B 检验磁性密度和交换分裂，A→C 检验非磁态中的相对论修正，B→D 则观察既定磁序加入 SOC 后的变化。

四组结果应使用相同晶格、原子位置、电子数、泛函、U 值、k 点和占据展宽。若四种状态各自充分弛豫，能量差会包含结构响应；固定同一几何构型时，能带差异更接近哈密顿量项的受控比较。结构弛豫效应可另列一组结果，避免把键长变化写成纯 SOC 信号。

非磁 WTe₂ 的计算将含 SOC 能带与无 SOC 能带叠加在同一条 Γ—X 路径上。红、蓝计算曲线之间的位移发生在没有磁有序的前提下，属于 SOC 对非磁电子结构的修正；时间反演仍约束 k 与 −k 的成对关系，晶体缺少空间反演时，远离时间反演不变点的能带可以分裂。

B 状态可在 SOC=0 时产生交换分裂。共线铁磁体中，V_xc^↑−V_xc^↓ 使多数与少数自旋带占据不同；反铁磁体还可能保留整体二重简并，却在局域扇区形成相反自旋极化。交换分裂的存在不要求 L·S 项非零。

CuMnAs 的无 SOC 计算给出整体自旋简并能带，α 与 β 两个 MnAs 扇区的投影自旋纹理却符号相反。整体能带简并、晶胞净磁矩为零、局域自旋极化非零 可以同时成立；只读取总磁矩会漏掉反铁磁子晶格和局域投影。

总能归因也可写成有限差分。ΔE_SP=E_B−E_A 包含磁性密度自洽带来的能量变化，ΔE_SOC|mag=E_D−E_B 表示既定磁态的 SOC 修正。δE=E_D−E_B−E_C+E_A 则量化两种作用在所选参考态下的非加和部分。

若只比较 A 与 D，新增局域磁矩、交换分裂、SOC 混合及结构响应会共同改变总能和能带。这样的端点差适合报告完整模型结果，却无法区分主导来源。局域磁矩、能带自旋期望值与三组能量差 应分别对应磁序、相对论项和交叉贡献。

无 SOC 的共线计算中，哈密顿量可按 ↑、↓ 两个通道分块。交换作用常使相同轨道成分的两组能带发生能量错位，DOS 也可分别画成 D_↑(E) 与 D_↓(E)。两通道占据积分之差给出自旋磁矩；峰位错开对应交换分裂的能量尺度。

DOS 曲线在费米能级 E_F 附近不对称时，可定位贡献磁矩的原子与轨道。反铁磁体的总 DOS 可能恢复上下通道对称，原子投影 DOS 仍在不同子晶格间交换颜色。总 DOS 对称只约束晶胞求和结果，并未消去每个磁性原子的局域交换场。

WTe₂ 没有磁有序，空间反演破缺与较强 SOC 仍产生动量相关自旋极化。自旋分辨 ARPES 在不同 k 点测得 S_y、S_z 分量的符号变化，计算旋量期望值呈相同趋势。SOC 型自旋纹理满足时间反演配对，k 与 −k 的自旋方向相反，晶胞净磁矩保持为零。

含 SOC 能带用红、蓝着色时，颜色通常表示某个量子化轴上的 ⟨σ_i⟩，并非严格的电子数标签。反交叉附近两条带会交换轨道与自旋成分，颜色连续过渡。SOC 旋量中的“上、下自旋”依赖投影轴，改变磁化方向或输出轴即可改变同一波函数的投影数值。

TMCM-MnCl₃ 在固定反铁磁构型下分别计算无 SOC 与含 SOC 能带。无 SOC 时 Γ 点附近保留二重简并，加入 SOC 后部分简并解除并出现小尺度分裂。磁序在两组计算中保持相同，新增分裂才可归入 SOC 与该磁序共同允许的对称性变化。

交换作用更常表现为同类轨道的多数、少数自旋能级整体错开；SOC 对近简并且角动量矩阵元较大的态影响更强，常见信号包括反交叉开隙、简并解除和自旋方向随 k 改变。轨道成分与对称性标签 能区分整体交换位移与局部 SOC 混合。

非共线自旋纹理本身也不是 SOC 的专属信号。Mn₃IrSi 的非共线磁序可由几何阻挫和交换作用形成，无 SOC 的第一性原理结果仍给出动量空间自旋纹理。SOC=0 的对照仍保留纹理时，纹理来源应归入磁序和自旋空间对称性；SOC 只对其作较小修正。

同一能带点的能量、轨道权重和三分量自旋期望值若随 SOC 开启共同变化，可把反交叉附近的混合描述得具体。若仅有两条分裂曲线而缺少 SOC=0 参照与磁序参照，“交换分裂”或“SOC 分裂”的命名仍缺少受控差分，图形相似并不能替代哈密顿量对照。

无 SOC 且忽略磁偶极相互作用时，整体旋转所有自旋不会改变电子总能，因为自旋空间与晶格坐标相互独立。SOC 加入后，L·S 把磁矩方向关联到晶轴，沿 [001] 与 [100] 的总能可出现差值。磁晶各向异性能量来自含 SOC 的方向差分，仅有自旋极化无法给出这项能量。

常用定义为 E_MAE=E_[100]−E_[001]，符号约定必须随文给出。两次计算要保持磁性构型、结构、k 网格和收敛阈值一致，能量精度需低于目标 μeV—meV 尺度。易磁化轴对应所用定义下能量较低的方向，方向差不能从单次总磁矩大小推算。

Ir₂ 的多数、少数自旋分子轨道在无 SOC 时已因自旋极化分开；加入 SOC 后，平行与垂直 Ir—Ir 轴的磁化方向产生不同能级移动。交换分裂提供自旋不等占据，SOC 再把占据轨道与空间方向关联，两项共同形成较大的磁各向异性能量。

轨道磁矩 m_L、磁晶各向异性、Dzyaloshinskii–Moriya 相互作用和各向异性交换都涉及自旋与空间方向的耦合，计算模型通常采用 自旋极化加 SOC 的非共线旋量。其中自旋极化建立磁序和交换场，SOC 决定轨道磁矩、易轴及允许的反对称交换项。

对轻元素闭壳层非磁体系，目标若是平衡键长、体模量或远离简并点的普通能带，A 状态可能已覆盖主要物理；含未成对电子的缺陷、自由基、过渡金属或稀土体系应检验 B 状态。重元素能带分裂、谷极化、轨道磁矩与磁各向异性 对 C 或 D 状态更敏感。

磁性重元素、磁性二维材料与非共线磁体常需要比较 B 和 D。若 SOC 开启后局域自旋磁矩变化很小，而 MAE、轨道磁矩或局部反交叉明显变化，可把前者归入稳定磁序，把后者归入相对论方向耦合。报告磁序、SOC 状态与磁化方向 可使总能差、能带颜色和 DOS 投影对应同一计算定义。