说明:本文华算科技主要介绍电子结构调控的物理含义、结构与轨道之间的传递关系,以及电子态变化怎样进入带隙、吸附和反应能量判断。




在第一性原理计算中,电子并不面对一个固定不变的能级表。原子核的位置和电荷、其他电子产生的库仑作用以及交换关联项,共同构成 Kohn–Sham 方程中的有效势 Veff(r)。电子态生成电子密度 n(r),电子密度又进入有效势,两者需要迭代到自洽。收敛后得到本征值 εi、波函数 ψi 和占据数。调控的直接对象是电子态的能量、空间形状与占据关系,结构或成分变化通过有效势作用于这些量。
元素替换会改写核电荷和价电子数,原子位移会改变键长、键角与轨道重叠,晶格位点交换会改变局域静电环境,外加电场则在空间中叠加势能梯度。自洽计算完成后,价带顶、导带底、局域态、交换分裂或成键—反键间隔可能重新排列。电子结构调控描述扰动变量—有效势—电子态之间的可重复比较。比较对象通常包含未扰动参考态和调控态,两者应采用相同电子数约束、能量基准与计算精度。
同一化学式在阳离子排序发生改变时也会产生不同电子结构。ABX2 型硫属化物可看作 NaCl 母结构的有序超结构,A、B 阳离子的排列方式不同,B 位周围会形成三角锥、跷跷板或近八面体配位。图1中的 NaAsS2、NaSbS2 和 RbBiS2 具有不同局域几何;电子所经历的周期势、键轴连接和轨道对称性随之改变。结构变量必须对应明确的电子态响应,此处的响应可沿键轴、带边和孤对电子分布继续追踪。





电子结构变化常从局域配位进入。配位数相同,电子环境仍可能不同。把图1中的三种结构暂时约束到高对称 NaCl 位点,可以沿 B–X 方向辨认两类原子序列。图2右侧把三条互相垂直的键轴单独画出:一类沿轴重复 B–X–B–X,另一类在 A、B 阳离子之间交替。键轴序列决定哪些 p 轨道能够正面对接,也决定局域畸变释放能量的方式。

B–X–B 型键轴倾向把相邻原子拉成一长一短的 B–X 对,形成常见的两中心—两电子键。A–X–B 型键轴中,一个 B 位两侧对应两个 X 阴离子,X–B–X 片段可形成三中心—四电子组合:电子占据较低能的成键态和中间非键态,较高能的反键态空置。图3画的是轨道计数关系。改变阳离子排序会改变两类键轴的数量,随后改变 B 位的配位畸变、孤对电子方向和带边轨道。

图4把轨道示意放回相同 NaAsS2 化学式的自洽电子密度中核验,三种竞争结构的电子局域函数采用相同的 0.82 等值面。三角锥配位中,As 位孤对电子朝远离三条短键的方向伸出;跷跷板配位同时保留定向孤对与多中心成键特征;近八面体环境中的局域区更靠近中心。ELF 在这里没有承担“键越强颜色越亮”的判断,它检验孤对电子与局域成键方向是否和键轴模型一致。同一等值面参数下的位置与形状对比,才能对应配位畸变造成的电子局域变化。





局域轨道进入周期晶体后展宽成能带。RbBiS2 的投影能带与态密度显示,最低三条导带含有明显 Bi 6p 与 S 3p 成分,最高三条价带主要来自 S 3p,紧邻其下的三条价带同时含有 Bi p 和 S p 成分。三组能带依次带有反键、非键和成键特征。图5把能带、PDOS、分子轨道计数和带分辨电子密度并列呈现,能量顺序、元素贡献和实空间形状由同一个成键模型描述。

调控带边时还要区分能级位置与能带色散。局域非键态接近价带顶,可能抬高价带边;轨道重叠增强会增大带宽,并改变载流子有效质量。同样大小的带隙可以对应陡峭或平坦的带边,也可以由不同元素轨道构成。带隙数值、带边轨道成分和色散应分别比较,光吸收阈值、迁移能力与缺陷俘获倾向才有各自的电子依据。
带隙的变化来自带边相对位置,不取决于某条峰是否单独移动。图6比较 24 组 ABX2 化合物在三类结构中的 HSE06 带隙。硫化物的带隙通常高于同组成骨架的硒化物,平均差约 0.37 eV,与 S 3p 价带位置较低相符。A 位从 Li 增大到 Rb 时,NaAsS2 型和 NaSbS2 型结构中的带隙可增加 1~2 eV;较大的 A 位离子增大共价 B–X 链之间的非键接触距离。RbBiS2 型结构的带隙对 A 位变化弱得多。同一种成分变量会因键连拓扑不同而产生不同电子响应,跨结构类型外推时必须保留这一差别。

催化体系中,目标量换成吸附轨道的占据和反应中间体能量。PtFe@FeSAs–N–C 中的 PtFe 合金与 Fe–N4 位点形成磁性耦合环境。图7显示,电位改变时有效磁矩和 Pt L3 边白线强度同步变化;轨道图把这种响应解释为 Fe 与 Pt 的 dz2 自旋电荷注入。Pt 5d 占据从接近 5d8 向 5d9 调整后,Pt–O* 与 Pt–OH* 的成键、反键占据不再相同,OH* 脱附所对应的键级下降。此处可用一个具体量描述电子结构调控:与反应物成键的轨道有多少电子、占据的是哪一种自旋分量,该占据随后影响吸附键级和脱附能量。





一个完整判断至少包含结构、电子态和性质三个比较对象。结构层要写清改变的是位点排序、键长、配位、应变、界面距离还是外场;电子层要指出响应发生在什么轨道、能量区间、自旋分量或空间区域;性质层再比较带隙、有效质量、吸附自由能、反应势垒、功函数或磁矩。若只看到结构已经变化,没有计算电子态,结论停留在调控手段;若只看到 PDOS 峰位不同,却没有对应的结构控制变量,因果关系无法定位到具体原子环境。
PtFe 与 Pt–FeN4 模型把电子态响应连接到氧还原反应能量。图8左上角的平面平均静电势给出 PtFe(111) 的真空能级与费米能级差,其余面板比较 PtFe(111)、Pt(111)、FeN4 以及 Pt(111)/FeN4 的反应自由能。电位为 0.85 V 时,不同模型的 OH* 质子化能量和决速步发生变化。轨道占据要传到指定反应路径的自由能差,才能用于催化性质判断;“电子更多”或“d 带中心移动”没有给出这一步能量对应。

三道比较还要共享计算口径。带隙对比应使用相同的泛函、晶胞和 k 点路径;PDOS 对比应统一能量零点、投影方式和展宽;自旋占据要检查初始磁矩与收敛终态;吸附体系要保持覆盖度、表面尺寸、溶剂和电位处理一致。图6采用 PBE 优化结构后再用 HSE06 计算带隙,图5用 PBE 分析轨道组成,两者回答的问题不同。方法变化造成的数值差异应从材料调控效应中剔除。
成分变化往往伴随几何弛豫,可以补做两组受控计算。其一把调控态的原子种类放回参考几何,检查核电荷与电子数造成的变化;其二保留调控后的几何,却采用参考态成分,观察键长与键角的贡献。两者与全弛豫结果比较,能够区分化学效应和几何效应。这种拆分尤其适合掺杂、空位、合金化和界面应变并存的模型。
计算模型还需对应可存在的结构。候选阳离子排序要比较总能、形成条件和动力学稳定;缺陷或掺杂要限定化学势与浓度;外场调控要给出场强范围;催化表面要考虑工作电位下的吸附覆盖和可能重构。电子态、结构稳定区间与目标性质应属于同一个模型,此时得到的关系才可用于检验材料机制。
