自旋轨道耦合指南：来源、晶体场效应及激发态矩阵元解析

说明：本文华算科技主要介绍自旋轨道耦合在哈密顿量中的来源、晶体场对有效耦合的调制、激发态跃迁与自旋弛豫中的矩阵元，以及第一性原理结果所支持的判断范围。

自旋轨道耦合的概念与哈密顿量

非相对论电子哈密顿量常把空间轨道与自旋写成两个自由度。对 Dirac 方程作低能展开后，电子在核电势梯度中的运动会产生相对论修正，自旋磁矩与轨道运动由此发生耦合。在近似球对称势场中，这一项通常写作 H_SOC = ξ(r)L·S。ξ(r) 在原子核附近较大，重元素的相对论效应通常更强，这也是 5d、6p 元素经常出现显著 SOC 响应的原子尺度来源。

L·S = [J² − L² − S²]/2。自由离子中，给定轨道量子数 l 的能级可按总角动量 j 分裂。晶体里的配位电场会先分裂原子轨道，晶体场能标、SOC 常数和交换作用共同决定实际本征态。此时 L 与 S 往往不再各自守恒，波函数要用含两个自旋分量的旋量表示；对选定量子化轴，旋量中的非对角矩阵元负责不同自旋分量之间的混合。

全相对论 Kohn–Sham 方程中的 SOC 项会同时改变本征能量以及波函数的轨道成分、自旋期望值、磁矩方向和跃迁选择规则。标量相对论处理已经包含质量速度项和 Darwin 项，却没有 L·S 引起的方向耦合；全相对论赝势或第二变分 SOC 会把这一矩阵项作用到 Kohn–Sham 态上。两种计算的总能或能带差记录能量修正，旋量投影记录波函数如何混合，二者回答的是不同物理量。

Cu₂N 表面 Fe 单原子的多重态给出了一个可分辨的实例。轴向与横向晶体场先产生能级结构，加入 λ_SO 后，各态获得不同的 ⟨S_z⟩ 和 ⟨L_z⟩；外磁场再解除部分简并。隧穿谱中的两个激发能标分别对应自旋投影改变与轨道投影反转。本征值差确定激发能量，隧穿矩阵元与选择规则控制相应峰强。

原子SOC常数为何会受晶体环境调制

同一种元素处在不同配位中，即使原子 ξ 近似相同，有效 SOC 响应仍可能差很多。晶体场把 d 或 p 轨道分裂成不同表示，成键又会改变轨道占据与空间延展。若占据态主要由实轨道组合构成，⟨L⟩ 可能接近零，轨道磁矩表现为淬灭；若两组轨道在能量上接近且对称性允许耦合，未占据态会通过 SOC 混入占据态，轨道磁矩和磁各向异性能随之增大。

Fe 位于 Cu₂N 的 N 顶位时，局域配位近似保留轴向特征，DFT 自旋密度集中在 Fe 及其近邻区域。这个构型具有接近自由原子值的未淬灭轨道矩，原因包含局域晶体场对轨道简并的保留程度，而非仅由 Fe 的原子序数决定。换成桥位、空位邻位或强共价配位后，晶体场分裂与轨道占据会改变，同一元素可出现另一组磁矩与各向异性能。

在弱 SOC 扰动图像中，某个占据态 |o⟩ 的能量修正常含有 Σ_u|⟨u|H_SOC|o⟩|²/(ε_o−ε_u)。能量分母较小只会提高允许矩阵元的权重；若两个轨道的空间对称性使矩阵元为零，近简并本身不会产生强混合。原子 ξ、晶体场能差、交换劈裂、轨道杂化和占据数共同决定固体中的有效作用，原子常数与能带劈裂量之间不会保持固定比例。

SOC 前后的本征值比较给出能量位移；轨道投影和 ⟨L⟩描述哪些态参与混合；磁矩随量子化轴的变化对应磁各向异性；双群不可约表示可判断交叉点是否允许开隙。若赝势是否含全相对论效应、磁化方向和非共线约束没有写明，几个 meV 的劈裂或各向异性能缺少可复算条件。

SOC如何进入系间窜越与自旋弛豫速率

SOC 对动力学的作用由跃迁矩阵元进入费米黄金规则。对初态 |i⟩ 和末态 |f⟩，速率可写成 k ∝ |⟨f|H_SOC|i⟩|²ρ，其中 ρ 表示满足能量条件的末态与振动相空间。SOC 矩阵元给出自旋变化的耦合幅度，能级差、振动重组和态密度决定跃迁是否拥有足够相空间。相同的 ξ 数值放在不同能级排列中，速率可相差多个数量级。

有机分子的系间窜越发生在单重态与三重态之间。纯自旋波函数在非相对论近似下正交，SOC 通过轨道角动量算符给予有限矩阵元。当跃迁伴随轨道类型改变，例如 n→π^* 与 π→π^* 态之间耦合，矩阵元常比轨道类型相同的组合更大。硝基芳香分子中的 S₁ 与 T₄ 自然跃迁轨道分别具有 ππ^* 和 nπ^* 成分，较小的单重态—三重态能差与有限 SOC 矩阵元共同提高系间窜越速率。

TPA-1N 在晶体中的硝基苯环扭转受分子间作用限制，S₁ 与低能三重态的轨道成分及能量间隔随之区别于气相优化构型。晶体构象的 ξ_S1T3 为 4.85 cm⁻¹，ΔE_S1T3 为 0.17 eV；气相低能三重态变为 T₂，ξ_S1T2 降至 0.91 cm⁻¹，ΔE_S1T2 增至 0.35 eV。几何构型同时改写矩阵元与能量匹配，固定原子组成并不保证相同的系间窜越速率。

固体中的自旋寿命服从相似结构，但末态由布里渊区中的电子态组成，散射源可来自声子、杂质或界面无序。SOC 先让 Bloch 旋量含有少量反向自旋成分，散射势再连接不同 k 点；第一性原理速率包含自旋翻转电子—声子矩阵元与散射态密度。锗烯在 GeH、SiH、GaTe 和 InSe 衬底上的杂化不同，其自旋寿命随温度及杂质面密度出现数量级差异，电场和衬底层数还会改变层间波函数分布。

能量劈裂VS波函数混合VS输运响应

在某个 k 点比较 SOC 前后的本征能量，常会得到两条能带的劈裂 ΔE，该量却没有包含声子占据、杂质浓度或末态数量。自旋弛豫率 τ_s⁻¹ 直接受 |g^↑↓|²D^S 控制，其中 g^↑↓ 是自旋翻转矩阵元，D^S 是满足能量守恒的散射态密度。能带劈裂与自旋寿命属于不同计算对象，前者可参与限制允许跃迁，后者还包含波函数重叠和散射相空间。

不同衬底锗烯的数据中，τ_s⁻¹ 与平均 |g^↑↓|²D^S、sin²(θ^↑↓/2)D^S 接近线性关系；以空间轨道重叠构造横轴时，数据离散更大。这里的 θ^↑↓ 衡量初末态自旋方向差。旋量的自旋混合方式比一个静态劈裂数值更接近实际翻转矩阵元；温度变化还会通过声子占据改变可用散射事件。

非中心对称晶体中的输运响应还可能依赖 SOC 波函数的几何量，分裂能带可产生强烈且动量依赖的 Berry 曲率 Ω(k)；费米面附近速度与 Ω(k) 的加权积分形成 Berry 曲率偶极。时间反演对称会约束相反动量处的 Berry 曲率符号，晶体对称性与费米面形状决定积分后哪些分量保留。非线性霍尔响应由费米面附近的曲率热点和占据共同给出，原子 ξ 没有单调对应的霍尔系数。

NbIrTe₄ 的计算在 Γ 点附近给出 SOC 分裂的部分占据能带，Berry 曲率集中于窄动量区域。化学势移动几十 meV 时，Berry 曲率偶极 D_a 的大小和符号明显变化；在 −20 meV 与费米能级两种占据下，动量空间正负贡献的抵消程度不同。致密 k 点插值和费米能级位置会直接影响这类积分量，即使原子 SOC 常数保持不变，载流子浓度改变仍可反转计算得到的非线性霍尔极性。