带电吸附态常出现在电催化模型里，例如 CO₂⁻、NO₃⁻、OH⁻、O₂⁻ 或带额外电子的极性中间体。直接把吸附物放到 slab 上，再用 NELECT 增减电子，VASP 可以给出一个自洽解；但这个解是否能拿来比较吸附能，要看参照态、电荷补偿和电位口径是否一致。带电吸附态可以计算，但不能把它当成普通中性吸附能直接排序，后续比较必须先限定电荷和电位口径。这个边界贯穿后面的模型选择。
图1 电位依赖的界面特异吸附与电荷转移示意。图中给出了吸附物、电极和界面电荷分布随电位变化的关系，可作为带电吸附态讨论的物理背景。 DOI: 10.1038/s41467-026-69559-x。
建模前需要先区分两种对象。第一种是吸附物本身在化学式上带电，例如 NO₃⁻ 或 OH⁻，但放到周期 slab 中时总体系是否带净电荷还取决于是否加入反离子、显式水层或电子数调整。第二种是中性吸附物在电极表面得到电子，形成具有 CO₂δ⁻ 或 O₂δ⁻ 特征的吸附态，此时总体系可能仍是中性，也可能通过 NELECT 固定为带额外电子的计算。
如果只看吸附物名字，很容易把两类情况混在一起。NO₃⁻ 从溶液来到表面，电荷来源包含溶液参照和界面补偿；CO₂ 在表面形成弯曲吸附态，则可能是中性总体系中的局部电荷转移。局部带电不等于周期超胞带净电荷，Bader 电荷、差分电荷密度和局域态分布只能说明电子密度偏移，不能直接等同于外加一个电子。这个区别决定了后面能量参照怎样设置。
VASP 常规静态计算更接近固定电子数。设置 NELECT 后，周期体系由均匀背景电荷保持整体可计算；体系总能包含了这个边界条件下的静电响应。电化学实验常讨论固定电位，电极可以和外部电子库交换电子。固定电荷和固定电位的自然变量不同：前者固定电子数，后者固定电子化学势，因此同一吸附态在两种口径下得到的稳定性排序可能不同。后处理时要先说明采用哪一种口径。
图2 的 CO₂ 电还原界面说明，阳离子、水和局域电场会同步改变 CO₂ 活化。若模型里没有显式离子或电位控制，只通过 NELECT 增加电子，得到的是某个固定电子数下的近似构型；它可以帮助判断电子注入后的结构响应，但不等于完整恒电位吸附自由能。
图2 Cu(100) 表面 CO₂ 电还原中碱金属阳离子与界面水的作用。图中展示了局域电场、阳离子和吸附构型对 CO₂ 活化的共同影响。 DOI: 10.1038/s41467-024-44896-x。
若 slab+X、slab 和气相/溶液吸附物都能定义在同一电荷口径下，可以写成式(1)。上标 q 表示整个周期体系的净电荷状态，Ecorr 表示因带电周期边界、偶极、真空层或其他修正带来的能量项。这个表达只在三个参照都明确时才有意义。
Eadsq=Eslab+Xq−Eslabq′−EXq−q′+Ecorr 式(1)
在实际任务中，q 和 q′ 不能后处理时再临时指定。若 slab+X 带 −1 电荷，而参照 slab 用中性计算，差值里不仅包含吸附，还包含向电极注入电子的能量、背景电荷响应和电势零点移动。固定电荷吸附能比较的是同一电荷边界下的相对稳定性；若参照电荷口径混用，它就不是独立的溶液离子吸附自由能。这个判断要写在结果说明中。
还有一类常见情况是只给吸附态加电子，却没有给参照 slab 建立对应的固定电荷状态。这样得到的能量差会把“吸附导致的稳定化”和“充电导致的能量变化”叠在一起。若要做固定电荷比较，slab、slab+X 和吸附物参照应先写成同一个电荷守恒反应。
对金属表面而言，额外电子往往在整个 slab 中离域分布；对半导体或氧化物表面而言，额外电子可能局域在缺陷态、吸附物 π* 轨道或表面金属 d 态。计算后应先确认电子实际进入哪里决定了固定电荷吸附能的物理含义，再把能量差和 Bader 电荷、PARCHG 或差分电荷密度一起核对。这个检查决定固定电荷结果是否仍服务于目标吸附态。
如果目标是电极电位下的吸附稳定性，应使用 grand potential 或近似的电位修正形式。式(2)把总能 E、电子数 N 和电子化学势 μe 联系起来，Ω 是固定电位下更合适的热力学势。实际计算中可通过恒电位方法、充电序列外推或电荷-电位曲线近似得到 ΔΩ。
Ω(U)=EDFT−Neμe(U) 式(2)
电位改变时不只是吸附态能量变了，电子库贡献也变了。若只在单个 NELECT 下算一个带电吸附态，缺少 Ne 与 U 的对应关系，就不能声称得到某个 RHE 电位下的吸附自由能。进行恒电位比较时，需要把电子化学势写进热力学势，否则不同材料之间的电位零点会混在能量差里。这个处理会改变材料排序的解释范围。
电荷-电位曲线提供了一种折中处理。对同一个 slab，可以计算多个 NELECT 下的功函数或界面电势，再把每个吸附态映射到相近电位区间。这个方法仍然是近似，但比只比较同一个额外电子数更接近电极问题，因为它把电容差异和电势漂移纳入了同一张结果图。
图3 的电场驱动 CO₂ 极化案例表明，局域电场会改变 CO₂ 取向和电子响应，因此电位或电场扫描比单点带电计算更接近机制判断。若目标中间体对电场非常敏感，还需要检查吸附构型是否随场强发生改变，而不是只把单个构型的能量做线性修正。
图3 强酸环境下电场驱动的 CO₂ 极化机制。图中给出 CO₂ 在界面场中的取向和电荷响应，说明固定电荷计算需要和电势响应分开解释。 DOI: 10.1038/s41467-026-68435-y。
第一步建议从中性 slab+X 做起。优化后检查吸附物是否保持目标构型，读取键长、Bader 电荷、差分电荷密度、PDOS 或 PARCHG，判断局部电荷转移是否已经足以形成目标吸附态。许多 CO₂、O₂ 和 NOₓ 中间体在中性总体系中就能出现明显极化，因此中性模型应先作为带电模型的结构参照，再进入固定电荷序列。后续所有带电结果都应和这个参照对齐。
若中性模型不能稳定目标中间体，再进入固定电荷序列。典型做法是以同一几何初态分别设置 NELECT = N、N+0.25、N+0.50、N+1.00，观察电子主要进入吸附物、表面态还是真空区。若额外电子分布在真空平台或非目标表面态，说明该带电吸附态在当前模型中缺少稳定物理承载对象。
固定电荷计算至少要同步检查真空层、LDIPOL/IDIPOL、LOCPOT 平面平均电势和表面偶极方向。带电 slab 的能量会随真空厚度、偶极修正和背景电荷处理变化；如果这些变化和吸附能差值同量级，吸附能排序就不应作为定量结论。对不对称 slab，IDIPOL 方向必须沿真空方向，且上下表面不能在比较中更换。
下面的 INCAR 只作为静态电荷序列参考示例，参数需要按体系收敛测试。NELECT 的数值应由中性体系价电子数加减得到，LDIPOL 与 IDIPOL 只在模型方向明确时使用。若是金属表面，ISMEAR、SIGMA、k 点和真空层应和中性参照保持一致。
ISTART = 1
ICHARG = 1
ENCUT = 520
EDIFF = 1E-6
NELECT = N + 0.50
LDIPOL = .TRUE.
IDIPOL = 3
LVTOT = .TRUE.
图4 的界面吸附状态随电位变化而改变，说明单个带电点只能描述局部条件。对于可发表的电化学判断，更稳妥的结果矩阵应包含中性态、多个固定电荷态、对应功函数或电势指标、吸附构型身份、关键键长和局部电荷分布。这样才能判断趋势来自吸附物电荷变化，还是来自整个 slab 的电势漂移。
图4 界面特异吸附下的电化学性能和吸附状态对比。不同电位区间的界面物种说明吸附态稳定性不应脱离电荷边界条件单独排序。 DOI: 10.1038/s41467-026-69559-x。
固定电荷结果适合回答“额外电子进入后结构怎样响应”和“目标中间体是否存在局部稳定构型”。例如 CO₂ 在额外电子或强局域场下弯曲，C─O 键长增加，表面金属─O 或金属─C 相互作用增强，这些现象可以作为 CO₂ 活化的计算证据。但它们不能自动变成某个实验电位下的绝对吸附自由能。
若比较同一个 slab、同一个吸附物、同一个 NELECT 序列，固定电荷吸附能可用于内部趋势。例如 N、N+0.25、N+0.50 三个电子数下，CO₂ 弯曲角、Bader 电荷和 Eadsq同步变化，可以说明电子注入增强了该构型的局部稳定性。这种趋势只在同一模型族内部成立，不能和另一个不同电荷、电位或真空设置的体系直接混排。跨材料比较需要重新建立电位标尺。
当结论需要比较不同材料在同一电极电位下的吸附自由能时，固定电荷单点就不够。不同材料的电容、功函数和表面态不同，同样增加 0.5 个电子会对应不同电位偏移。此时应把 U 作为可比变量，通过电荷-电位曲线、恒电位 DFT、隐式溶剂/显式离子补偿或外电场扫描组织结果。这样得到的趋势才更接近同一电极条件。
若吸附物是溶液离子，例如 NO₃⁻ 或 OH⁻，还需要同时定义溶液参照。直接把孤立 NO₃⁻ 放在真空盒中算气相能量，通常不能代表水溶液离子化学势。更可取的处理是使用 CHE、实验化学势、溶剂化校正或反应式组合，让带电物种通过中性反应循环进入自由能表达。
结果表达可以分成三层。第一层报告模型边界：NELECT 改变量、是否使用背景电荷、真空层、偶极修正和电势读取方式。第二层报告结构身份：吸附物是否保持目标构型、键长和局部电荷是否支持目标吸附态。第三层报告可比范围：该结果用于同一模型族内的电荷响应趋势，还是已经通过 ΔΩ(U) 转换为同一电位下的自由能。
图5 的 CO₂ 活化案例把几何构型、电子描述符和吸附态响应联系在一起。对带电吸附态而言，也应采用这种多证据判断：弯曲角、键长、Bader 电荷、PDOS/PARCHG、差分电荷和能量项需要互相支持。只有电荷边界、吸附身份和热力学参照同时明确，带电吸附态结果才适合进入材料比较。否则它只能作为局部响应的辅助证据。
图5 半导体氧化物上 CO₂ 活化的结构与描述符。图中把吸附构型、电子描述符和 CO₂ 活化联系起来，可用于说明电荷态结论需要几何与电子证据互校。 DOI: 10.1038/s41467-022-28042-z。
如果目标是探索某个吸附态在额外电子下是否稳定，可以直接做固定电荷序列。此时重点不是给出绝对吸附自由能，而是观察局部结构、电荷分布和电子态怎样随 NELECT 改变。适合的对象包括 CO₂ 活化、O₂ 超氧/过氧特征、NOₓ 中间体电子注入和强极性吸附物的局域响应。
如果目标是比较不同催化剂的电化学活性，直接固定一个 NELECT 值通常不够。不同 slab 的电势响应不同，同一个电子数改变不代表同一个工作电位。此时应把功函数、LOCPOT 真空平台、双电层补偿或 grand potential 纳入计算流程，再讨论电位下 ΔG 的材料排序。
要点：带电吸附态可以用 VASP 计算，但固定电荷结果不能直接等同于恒电位吸附自由能。
要点：局部带电和周期超胞带净电荷不同，Bader 电荷与差分电荷只能说明电子密度偏移。
要点：固定电荷吸附能必须明确 slab、slab+X 和吸附物参照的电荷口径。
要点：材料间电位下比较应优先使用 grand potential、电荷-电位曲线或显式补偿模型。
