说明:本文华算科技主要介绍能带理论中的绝热近似、单电子近似和周期场近似,以及它们在第一性原理电子结构计算中的含义与适用条件。

固体物理教材在建立晶体能带理论时,常用的三大近似是绝热近似、单电子近似、周期场近似。它们分别处理原子核与电子的耦合、多电子之间的作用,以及晶格势的空间重复。晶格振动中的谐近似处理平衡位置附近的势能展开,讨论对象与这组三大近似不同。
由原子核和电子组成的固体,其非相对论哈密顿量包含 H=Tn+Te+Vnn+Vee+Ven。前两项是核与电子的动能,后三项分别是核间排斥、电子间排斥和核—电子吸引。整组核坐标记为 R,整组电子坐标记为 r;电子位置变化会影响其他电子,核位置变化又会改变电子所处的外势。
供体与受体之间的质子—电子模型保留了两类坐标:质子位移改变电子所处的势,电子分布随之调整。论文中的供体和受体固定,质子与电子沿连线运动;这是一维分子模型,用于区分电子自由度与核自由度,晶体则包含大量核和电子。
绝热近似允许在给定核构型下求电子态,再用电子能量描述核的运动。单电子近似把多电子问题组织成有效单粒子方程,周期场近似进一步要求有效势具有晶格平移对称性。三者针对不同物理条件:核的运动快慢、电子关联强弱、缺陷或无序破坏的平移周期,可能分别影响计算结果。

固定核坐标怎样生成势能面?
原子核的质量远大于电子,许多基态过程的核运动时间尺度较长,电子态可以随核构型近似绝热地调整。电子结构计算把 R 作为参数,求解 He(R)Φj(r;R)=Ej(R)Φj(r;R)。这里 He 包含电子动能和相互作用;若将核间排斥计入 Ej,各核构型的本征能量就组成相应的势能面。
沿一个核坐标扫描时,势能面显示为一条曲线。不同颜色的电子态能量随核位置变化,每个离散点对应一次固定核构型的电子本征问题。基态曲线的极小值给出平衡构型候选,曲线的斜率决定核受力,局部曲率决定小振幅振动的恢复力。
Born–Oppenheimer 近似常把总波函数写作 Ψ(r,R)≈Φj(r;R)χ(R),在选定电子态的势能面上处理核运动。结构优化沿能量梯度调整核坐标,Born–Oppenheimer 分子动力学在每个瞬时构型重新求电子态并推进核的轨迹。核波函数也可以在该势能面上作量子力学求解;核的经典处理是另加的选择。
什么会使单一电子态的描述失效?
两个电子态能量接近、电子波函数随核坐标快速变化,或核的运动较快时,不同电子态之间的非绝热耦合可能变得显著。电子与核运动的分离依赖能隙、耦合和核速度,质量差只给出常见的物理依据。光激发后的载流子弛豫、电子转移等过程,常涉及多个电子态上的核运动。
质子—电子模型中,提高初始质子动量后,核波函数在其他电子态上出现非零振幅。振幅在电子态之间转移超出了单一势能面的动力学描述。声子计算和热运动仍可在绝热框架内开展;此处“绝热”描述电子态与核运动的分离,和热力学中系统是否交换热量的用法不同。

有效单粒子方程里保留了什么?
N 个电子的波函数依赖全部电子坐标,多电子相互作用使求解代价迅速增长。单电子近似用轨道和有效单粒子算符描述电子态,把其他电子的作用通过平均场、交换项或交换相关泛函处理。方程中的“单电子”指求解对象,固体中仍有 N 个电子,轨道占据仍须遵守泡利原理。
在常规 Kohn–Sham DFT 中,辅助电子在有效势中满足 [−ℏ2∇2/(2me)+vs(r)]φi=εiφi,并由 n(r)=∑ifi|φi(r)|2 重建密度。辅助体系按基态密度与真实相互作用体系匹配,fi 是轨道占据数。
Hartree 处理保留电子密度产生的经典库仑场,Hartree–Fock 通过反对称波函数引入交换,其交换算符通常非局域。Kohn–Sham 方法引入非相互作用辅助体系,将剩余的多体贡献写入交换相关能;精确泛函原则上可恢复基态密度与能量,实际计算中的 LDA、GGA 等才对未知泛函作近似。辅助轨道互不直接作用,真实体系的电子排斥仍参与能量与密度的确定。
有效势随密度改变,轨道又决定新密度,自洽求解要求二者相容。自洽收敛检验的是所选方程的解;交换相关近似遗漏的物理贡献仍可能影响结果,单纯增加迭代次数无法恢复这些贡献。
单电子能级与实验激发谱怎样区分?
Kohn–Sham 本征值属于辅助体系,而光电子谱探测电子移除后的多体激发。局域电子的库仑排斥与能带展宽相竞争时,Mott 绝缘行为需要电子关联的描述。把普通能带中的空隙直接解释为某种绝缘机制,会遗漏谱权重与激发态能量的变化。
层间二聚化的 Nb3Br8 中,Mott 模型和能带绝缘模型都可产生能隙,但面外动量方向的谱强度分布不同。原研究比较 ARPES 与两类模拟,最高占据带顶的强度重复间隔符合 Mott 相的预期。这里区分两类绝缘态所用的是动量分辨谱,而非一个带隙数值。
能带计算中的平滑色散与实验谱中的展宽、卫星峰涉及不同电子激发信息。随频率变化的多体自能可改变峰位、寿命和谱权重;静态有效单粒子能级没有完整包含这些频率依赖信息。

理想晶体中,原子排列按晶格矢量 T 重复。若电子态没有破坏该平移对称性,有效势满足 veff(r+T)=veff(r)。对理想周期模型,这是明确的对称性条件;用它描述有限、含缺陷或存在热位移的样品时,才涉及忽略实际结构偏离的近似。
石墨和石墨烯上的 MBB-1、MBB-2 分子自组装形成条纹晶格,上排与下排分别对应两种分子。重复排列的结构单元给出空间周期,分子官能团还会改变局部静电环境。STM 图记录高度或隧穿电流对比,计算电势则来自另一个电子结构计算量。
在周期单粒子哈密顿量下,布洛赫定理允许选择 ψnk(r)=eik·runk(r),其中 unk(r+T)=unk(r)。平移改变波函数的相位,概率密度保持周期。按能带序号 n 与晶体波矢 k 求得 En(k),就得到能带色散;周期条件没有预先指定具体的带宽和带隙。
两种超分子晶格可以保持相同的排列周期,却产生不同幅度的电势调制。原图色标单位为 V,表示差分电势 φ,电子势能中的静电贡献按 −eφ 换算。空间重复确定允许的平移,势的幅度、形状及轨道耦合继续影响电子能量。
近自由电子模型把周期势视为对自由电子的微扰,紧束缚模型从局域轨道及其耦合构造能带;两者是周期场问题的求解模型,没有替代绝热、单电子与周期场这三项基本处理。分子和团簇也能用 DFT 求电子结构,但通常不具备晶体的平移周期。
含空位或杂质的超胞计算将局部缺陷连同周围晶格周期排列,恢复的是超胞周期,原胞平移通常无法保持缺陷位置不变。增大超胞会改变镜像间距及缺陷能级展宽;孤立表面的理想 slab 模型则在面内保留晶格平移,常用三维周期超胞另设真空层来分隔相邻 slab。
