内置电场及d带中心理论

说明：本文华算科技主要介绍内置电场与 d 带中心的物理定义、耦合途径及第一性原理计算中的区分方法。

内置电场写成 E(r)=−∇φ(r)，描述体系内部电荷分布产生的静电势梯度。题目沿用“内置电场”，论文中也常写作“内建电场”。它可以来自异质界面两侧的功函数差、极性晶格、Janus 结构上下表面的化学差异、空间电荷区或局域配位破缺。场强具有空间位置和方向，常用 V Å⁻¹ 表示；平面平均静电势 φ(z)、平面平均差分电荷密度 Δρ(z) 与局域介电屏蔽共同限定实际势降。

平衡态中的“内置”指无外接电源仍存在的势场分量。金属—半导体或两种导体接触后，体系的电化学势达到平衡，接触区的 φ(z) 仍可变化。势能平台之间的阶跃对应界面偶极，有限厚度区域内的连续斜率才能换算局域场强。总势差 ΔV 要与势降长度一同报告，相同 ΔV 分布在不同距离内会产生不同 E。

d 带中心 ε_d 是指定原子 d 投影态密度的能量一阶矩，常写成 ε_d=∫E·n_d(E)dE/∫n_d(E)dE。它把一段有宽度、有峰形、可分自旋和轨道分量的 d-PDOS 压缩为一个能量坐标。数值随投影原子、积分能窗、自旋处理及能量零点变化；跨模型比较时，这四项应保持同一口径。

一个量记录实空间中的势梯度，另一个量记录能量空间中的局域 d 态重心。两者经界面电荷迁移、轨道占据变化、晶格弛豫和配位重构发生关联。非对称 InO₅S 中，功函数、电荷重排、电子局域函数与 d 带中心随场强的变化被置于同一原子模型，场由何处产生、局域态怎样响应因而具有共同的结构参照。

场强与 ε_d 的相关曲线属于特定结构族。金属屏蔽、局域轨道对称性和 d 电子填充都会改变曲线斜率；某些位点出现明显电势变化，d-PDOS 仍可能只发生小幅重排。论文中的相关性要绑定原子位点、场方向和能量参照；省略任一项，“上移”或“下移”便缺少可比较的物理坐标。

两种材料接触前具有各自的功函数和电子化学势。接触后，电子迁移持续到整个体系共享费米能级，界面两侧留下空间分离的正负电荷，φ(z) 随 z 产生势阶或斜率。功函数差给出电子迁移倾向，Δρ(z) 定位电荷重排，φ(z) 定位势降；三种量承担不同物理信息。

局域势改变时，所有电子态可能先发生近似同幅的能量平移。把两个模型都以各自 E_F=0 eV 作图，局域势差可能被吸收到零点中；改用真空能级、深层芯态或局域平均静电势对齐后，谱线若重新重合，主要变化来自能量参照移动。谱峰相对局域参照仍有位移，才指向轨道占据、杂化或晶场变化。

界面电荷进入金属位点后，d 轨道的占据数、自旋分裂和峰权重会改变；配位原子随场弛豫，键长、晶场分裂和带宽又产生结构反馈。ε_d 的位移可能来自刚性谱移、谱权重重排、带宽变化或自旋分支重配，相同的质心差值可以对应不同电子图像。

自由载流子常把金属中的静电场压缩到表面数个原子间距，半导体耗尽区和极性绝缘体可保留较长的势降区。场的穿透长度决定哪些金属原子出现 d 态变化，也决定表面位点能否感受到埋藏界面的极化。逐原子或逐原子面计算的 d-PDOS能够检验这种空间衰减，总 DOS 会把界面响应与体相信号混合。

MnO₂ 模型中的定向场使 Mn 3d 与 O 2p 态靠近 E_F，d_x²−y² 分量响应尤为明显，Mn–O 的 −pCOHP 与 ICOHP 也随之改变。这里的 d 带中心变化带有明确的轨道选择和成键后果，信息量高于一条总 d-PDOS 质心数值。

吸附物前线轨道靠近过渡金属位点后，会与对称性匹配的 d 轨道耦合，生成较低能的成键态和较高能的反键态。对一组结构相近、d 电子填充相近的金属表面，ε_d 向 E_F 靠近时，反键态常向较高能量移动，占据减少，吸附键随之增强。真正控制键强的是耦合矩阵元、能级间距和反键态占据，ε_d 负责概括其中一部分趋势。

同一 ε_d 可对应宽度和填充数不同的谱形：窄带在 ε_d 附近形成高峰，宽带把态权重分散到较大能区；自旋极化位点的多数与少数自旋反键态占据也不相同。O^*、CO^*、H^* 的轨道对称性差异明显，d_z²、d_xz/d_yz 或 d_xy/d_x²−y² 的分轨道中心可直接对应实际成键原子对。

FeN/Fe₃N 界面把这一路径具体化：Fe₃N 向 FeN 一侧的电荷极化形成场，界面 Fe d 态重心和 E_F 附近态密度随之变化，COOH^* 与 CO^* 获得不同的电荷稳定。自由能图把电子态变化对应到 CO₂→COOH^* 的能量代价，避免用 ε_d 单独替代整条反应路径。

电场还可作用于吸附态自身的偶极矩 μ 和极化率 α。小场近似下，场致能量项可写成 ΔE_field≈−μ·F−½αF²。当反应物、过渡态和产物具有不同 μ 或 α，同一场强对各状态的稳定幅度不同；即使金属 d-PDOS 变化很小，吸附能或能垒仍可发生可测差异。直接场项与 d 态介导项在同一反应中可以共存。

InO₅S、InO₆ 与 InS₆ 的 d 带中心差异进入 NO 还原后，作用结果体现在 In–N 成键、N–O 断键与多步质子电子转移中。−ICOHP 和能量分辨 COHP 指向具体原子对，自由能曲线给出每个中间体的热力学位置。d 带中心适合提出吸附趋势，COHP 和自由能负责核对趋势对应哪根键、哪一步反应。

界面模型与分离组分模型应保持相同原子数、晶胞和冻结几何，所得 Δρ(r)、Δρ(z) 与差分静电势负责识别接触诱导项。显式外场模型可在同一构型上使用 +F、0、−F 三个状态；冻结离子得到电子响应，弛豫离子后的差额对应结构响应。场反向后能量的一次项随符号反转，二次极化项保持同号，可用来辨认偶极与极化率贡献。

d 态响应应标明具体原子和轨道。位点分辨、自旋分辨 d-PDOS 经统一能量参照后计算 ε_d，再与带宽、Bader 电荷、轨道占据及 COHP 对照。金属载体会屏蔽远离表面的场，纳米尖端与界面又产生强烈空间起伏；用整颗粒 d 态求一个平均 ε_d，容易掩盖尖端富电子位点和界面缺电子位点的相反响应。

+F 与 −F 若主要改变吸附态能量，而清洁表面的 d-PDOS 经局域势对齐后近乎重合，直接偶极—场项占比突出；位点 Bader 电荷、分轨道 PDOS 与 COHP 若随 F 成组变化，电子结构介导项已有数值支撑；差异只在离子弛豫后出现时，键长和配位变化承担主要贡献。三套信号对应三种可检验的物理来源。

电化学工作态叠加电极电位、溶剂取向和双电层离子产生的运行场。真空 slab 得到的内置分量适合描述固体接触后的自发极化，溶液中的活性位能量则随电位和覆盖度变化。无偏压界面与恒电位界面的 ε_d、吸附能应分别标注模型条件，无偏压势场不宜外推到任意工作电位。

Ir/NiPS₃ 中，Ir 与载体接触区向 NiPS₃ 失电子，Ir 尖端区域却因几何场增强和簇内电荷迁移呈富电子状态。平面平均势给出总体势降，原子分区和局域场矢量揭示颗粒内部的非均匀分布。同一界面含有多种局域电子环境，吸附物所在原子决定哪一组 d 态进入成键。

局域场矢量、KPFM 表面势和原位振动信号为势场分布给出实验约束；位点分辨能垒把尖端、界面和载体的反应差异写入同一能量坐标。计算报告若同时列出场方向与强度、ε_d 的投影对象和能量参照、吸附态偶极、COHP 及自由能路径，电场直接稳定吸附态、场诱导 d 态重排和几何弛豫三类贡献便可获得各自的数值证据。