说明：本文华算科技主要介绍 PDOS 峰位移动所对应的能量参照、局域势、晶场与交换分裂、轨道杂化和谱权重变化，并说明哪些情况下不能把曲线横向位移直接解释为电荷转移。
PDOS峰位移动是什么量的变化？
分波态密度把每个 Kohn–Sham 态投影到指定原子、轨道、层或自旋分量。若第 i 个本征态的能量为 εi，它在投影轨道 φα 上的权重为 Piα，绘图时常把离散能级用展宽函数 g 展开：
曲线峰顶来自一组能量相近、投影权重较大的电子态。峰顶不等同于单个原子轨道能级，峰形还受能带色散、简并度、投影半径和展宽宽度影响。比较前应先固定峰位采用局部最大值、峰包重心还是带边起始位置；三个量可能沿不同方向变化。
图1. RuO2 的能带投影、Ru 中心 Wannier 轨道与 ARPES 能谱。DOI：10.1038/s41467-020-20252-7
RuO2 的例子把这一计算关系画得很清楚：能带上的颜色来自 Kohn–Sham 态对 Ru 中心 Wannier 轨道的投影，d∥ 与 dxz、dyz 在动量空间占据不同区域。把投影权重沿整个布里渊区求和才得到相应 PDOS。曲线上的一个峰往往汇集多条能带，单看峰顶无法追踪某一条能带。
横轴以何处为零会直接影响“移动”的表述：两次计算各自把费米能级 EF 设为 0 eV 时，载流子浓度改变会让化学势移动，所有横坐标都随新的零点重新标定。此时图上峰位相对 0 eV 的变化，包含能带相对化学势的变化。若要比较不同表面、掺杂浓度或带电体系的绝对能量，还需对齐真空能级、深层芯能级或远离缺陷区域的平均静电势。
同一能量参照仍不够。不同赝势、交换关联泛函、U 值、晶胞尺寸、磁构型、投影方式和展宽参数都会改变峰形。尤其是窄峰，展宽从 0.05 eV 改成 0.20 eV 便可能让相邻峰合并，最大值转到新的位置。只有投影定义与数值设置相同，横向差值才适合讨论材料因素。
哪些电子结构变化会推动峰位移动？
掺杂、界面偶极、外电场和缺陷电荷会改变局域静电势，随之调整某个原子附近电子态的能量。电子流入某个片段时，库仑作用和屏蔽会改变局域势，但峰向低能移动不能单独充当得电子的计量结果；电荷差分、Bader 电荷或占据矩阵才负责判断电荷集中于哪些空间区域与轨道。
图2. 外电场下 TCNQ/石墨烯的隧穿谱、分波态密度及 HOMO、LUMO 能量变化。DOI：10.1038/s41699-019-0087-5
TCNQ/石墨烯体系给出了可辨认的非刚性移动。外电场反向后，分子的 HOMO、LUMO 峰与石墨烯 PDOS 朝相反方向变化；计算把它归因于场依赖的分子—石墨烯电荷转移。若只是全体系化学势平移，两组谱线相对彼此的距离应近似不变。不同投影对象出现相反位移，说明局域势和界面偶极发生了重新分配。
键长、键角和配位环境改变会重排晶场能级，例如八面体拉伸可能降低一组 d 轨道、抬高另一组，对称性下降还会解除原有简并。此类变化常表现为特定轨道选择性移动或分裂，而非所有原子和轨道整齐平移。应变计算中需要同时检查优化后的局域结构，才能把峰位差与晶场变化对应起来。
图3. 体相与外延应变 RuO2 的计算能带、态密度及 ARPES 能谱。DOI：10.1038/s41467-020-20252-7
外延应变下的 RuO2 平带向 EF 靠近，7 nm 高应变薄膜相对 48 nm 较松弛薄膜的位移约为 120 ± 20 meV。计算能带、态密度和 ARPES 能量分布给出相同趋势，近费米谱权重随 d∥ 态靠近而增加。这里可讨论的是应变对轨道相对能量和占据的调节，数值仍依赖所用 DFT+U 模型。
磁性体系还要分开观察两个自旋分量。交换作用可把原本简并的峰推向不同能区；改变磁矩、缺陷浓度或 U 值时，峰间距也会变化。若把上下自旋相加，两个峰可能合成宽峰，视觉上像峰顶偏移。自旋分辨 PDOS 与局域磁矩应成组核对，否则容易把交换分裂误写成整体能级下降。
峰位变化为何常伴随分裂、展宽和谱权重转移？
两个轨道发生耦合后，体系的本征态由原来的局域轨道线性组合而成，成键态与反键态分居不同能区。投影到某个原子时，原有峰可能分成两支，也可能在连续能带中展宽。此时峰顶移动只是耦合后谱函数变化的一部分，还应观察峰宽、峰面积、肩峰及不同原子峰形的同步变化。
图4. 不同金属单原子接触 C60 后的 C60 投影态密度。DOI：10.1038/ncomms4659
C60 与 Cu、Au、Pd、Fe、Ti、Al 单原子形成接触后，分子共振相对孤立 C60 的虚线向低能移动，Al 和 Ti 对应的位移较大；Fe、Ti、Pd 附近还出现与金属 d 共振重合的附加结构。同一张 PDOS 同时记录电荷转移和杂化，仅用“峰向低能移动，所以成键增强”会漏掉接触元素、峰宽和新共振携带的信息。
缺陷会打破局域对称性，并改变邻近原子的配位数和电子数。原先等价的 O 原子变成多类局域环境后，某一 p 轨道可能穿过 EF，另一轨道仍留在价带深处。跨越费米能级还涉及占据变化，判断时要记录峰位于占据侧还是未占据侧，而不能只报“左移”或“右移”。
图5. 本征与不同空位 Bi0.5Na0.5TiO3 中邻位 O 原子的轨道分辨 PDOS。DOI：10.1038/s41598-020-62889-w
Bi0.5Na0.5TiO3 中引入 Bi、Na 或 Ti 空位后，邻位 O 的 px、py、pz 峰形和上下自旋不再保持本征体系中的一致性。论文据此把缺陷诱导磁性定位到空位附近 O 2p 态。这里的证据来自轨道选择、自旋不对称与费米能级穿越的组合，并非某个峰的单一位移。
向同一晶胞注入 1 或 2 个电子后，Ti 3d 态向 EF 附近重排，两个自旋分量分开，部分多数自旋态开始占据。每幅子图又各自把 EF 设为零，因此能量参照变化和电子态重构叠加在横轴上。若研究目标是绝对带边移动，应另做静电势对齐；若研究占据和磁性，零点设在 EF 更便于观察跨越行为。
图6. 向 Bi0.5Na0.5TiO3 晶胞注入 0、1、2 个电子时 Ti 3d 轨道的自旋分辨 PDOS。DOI：10.1038/s41598-020-62889-w
有些峰在改性前后已不再代表同一组本征态。峰 A 与峰 B 接近后发生避免交叉，原先属于原子 1 的谱权重可能转移到原子 2，追踪“最高峰”便会换到另一支能带。峰的轨道成分是否连续比横坐标差值更先需要确认，可借助投影能带、波函数局域分布或 Wannier 轨道追踪同一电子态。
怎样确认峰位移动对应哪种物理过程？
比较计算应保留同一晶胞取向、原子轨道投影、能量网格、展宽和磁序，再分别输出以 EF 为零与经静电势对齐的两套横轴。若前者明显移动、后者基本重合，主要变化来自化学势；若对齐后某个原子或轨道仍移动，则需考察局域势、晶场和耦合。排除坐标基准造成的表观位移，可以减少后续机理判断的歧义。
窄峰可记录峰顶，宽峰或多峰包更适合计算能量重心 εc＝∫E D(E)dE/∫D(E)dE，积分窗口应在各体系中保持一致。再把峰宽、积分谱权重和 EF 附近占据列在同一张表中。峰顶、重心和面积给出的变化方向可能不同，三者分开报告能避免用一个数字概括整段谱形。
图7. 裸露与 FAS17 修饰 LNMO 表面的价带 XPS 和原子轨道 PDOS。DOI：10.1038/s41598-018-30135-z
LiNi0.5Mn1.5O4−δ 表面固定 FAS17 分子后，计算得到靠近 EF 的 Ni 3d、Mn 3d 和 O 2p 峰向低能移动，价带 XPS 的边缘变化与之对应。计算与实验共享的是能量变化趋势，峰宽和绝对位置仍会受光电子矩阵元、仪器展宽、表面缺陷及 DFT 能级误差影响，不能逐峰强行配对。
机理核对可按信号类型组合。电荷转移需要电荷差分和原子电荷；晶场变化需要键长、键角、轨道分裂与局域对称性；杂化需要两侧投影峰形、投影能带以及 COHP/COOP 或 Wannier 耦合；磁性需要自旋分辨 PDOS、局域磁矩和自旋密度。PDOS 负责定位能量区间与轨道来源，其他量负责检验空间分布、成键符号和磁矩来源。
对比对象也限定结论范围。同一成分的应变序列可较干净地归因于结构变化；不同化学元素掺杂会一并改变电子数、局域势、尺寸失配和磁交换，峰位差无法自动分解为单一贡献。可设置固定几何但替换元素、固定元素但改变几何、以及中性补偿电荷等对照计算，分别观察各因素造成的谱线变化。对照模型越接近单变量，峰位与物理因素之间的联系越清楚。
实际判读时，把“哪一条投影、相对什么零点、用什么峰位定义、是否发生分裂或谱权重转移、有哪些独立证据”写在结果旁边。表述需要指向可核查的计算对象，例如“真空能级对齐后，界面 O 2p 峰包重心降低 0.35 eV，并伴随金属 3d 权重进入同一能区”。数值、轨道、参照和伴随信号齐全，才足以界定 PDOS 峰位移动所对应的电子结构变化。
