二维半导体加外电场后，带边位置和吸附自由能怎样同步变化？

引言

二维半导体上加外电场，最容易出现两个结果各说各话：能带图显示 VBM 或 CBM 移动，吸附自由能又给出另一个排序。真正要比较的是同一组场强、同一方向、同一真空参照下，带边、表面静电势、吸附构型和 ΔG 是否一起改变。只看能带会忽略吸附态；只看吸附能又会丢掉电荷转移动因。同步扫描的目标，是把电子结构位移和反应端自由能放回同一个模型里。

外电场方向先怎样定义？

二维 slab 计算里，外电场方向和吸附侧必须成对定义。方向不能只写成正负号，而要说明电场从哪一侧指向哪一侧、吸附物位于哪一面、偶极修正是否打开。同一材料翻面后，正场和负场的物理含义会互换，所以目录名、结构文件和图注都应带有 top 或 bottom 这样的方位标记。

第一步不急着放吸附物，而是对干净二维半导体做几组场强的静态或轻度弛豫。每组使用相同真空层、相同面内晶格和相同 k 点。若真空层太薄，外场区与周期镜像会互相干扰；若未做偶极修正，真空平台可能倾斜到难以取参照。此时能带移动既可能来自真实极化，也可能来自数值参照漂移。

实际建模还要先决定原子是否在外场下重新弛豫。只做静态场强扫描，能更清楚地看电子响应；允许离子弛豫，则会把压电、层间距变化和吸附位点重排一并纳入。两者都可以做，但建议把电子项和结构项分成两张表，否则很难判断 ΔG 改变来自带边移动还是来自键长重新优化。

图1显示 g-GaN/Sc₂CO₂ 范德华异质结构的模型与电子结构背景。它适合作为外场扫描前的母体提醒：层间方向、界面电荷转移和真空方向一旦定义不清，后续带边与吸附态排序都会失去可比性。

图1　g-GaN/Sc2CO2 异质结构模型与电子结构背景。 DOI: 10.1039/d2ra00419d。

带边位置用哪个零点对齐？

外电场下的能带本征值不能直接跨场强相减。应先从 LOCPOT 或静电势后处理中提取远离 slab 的真空平台，再把 VBM、CBM 或投影带边换算成相对于同一真空零点的能量。这里真正比较的是 E_VBM − E_vac 和 E_CBM − E_vac，不是 OUTCAR 里孤立的本征值序号。

对二维半导体，带边还要确认来自哪一层、哪一侧或哪类原子轨道。若外场把 VBM 推向上表面，把 CBM 推向下表面，表面吸附物只与其中一侧强耦合。此时全局带隙变化并不能直接代表吸附反应活性；需要用投影能带、PDOS 或部分电荷密度说明吸附侧真正看到的带边。

如果体系是异质结，还要把层分辨带边单独列出。外场可能让总带隙变小，却把反应侧带边推离吸附物；也可能让另一层承担主要光生载流子分离。带边同步表至少应含全局带边、吸附侧投影带边和真空平台，这样才能区分能级移动与空间分离两种效应，并判断吸附物实际接触的是电子富集侧还是空穴富集侧。

图2给出 g-GaN/Sc₂CO₂ 在外电场调控下的带边响应。做自己的 VASP 后处理时，可以把这类图拆成三列：场强、真空对齐后的 VBM/CBM、吸附侧局域态贡献，避免把整体能带移动误判为某个位点活性变化。

图2　外电场调控下的带边与光催化相关能级变化。 DOI: 10.1039/d2ra00419d。

吸附自由能怎样跟随场强算？

吸附自由能必须和干净表面的外场设置严格配对。对每个场强 E，需要分别优化 clean slab、gas 或参照物种、slab+X；若气相盒子不施加同样外场，报告时要明确参照处理方式。对于带偶极的吸附物，要提前枚举不同几何取向，因为同一个初始构型在正负场下可能收敛到不同吸附模式，不能只复制零场 CONTCAR。

常用自由能可以写成电子能差加校正项。若以 H* 为例，至少要记录 ΔE_H*、ZPE、熵项和电位项是否随场强重算。外场本身不是电极电位的完整替代，它提供的是局域势和界面偶极扰动；若要换算到 RHE 或 SHE，还要另建电位参照。

吸附端的目录可以按 Efield_m0p2、Efield_0p0、Efield_p0p2 分组，每组内固定 clean、ads、freq 和 post 四类目录。这样得到的 ΔG(E) 才能与带边 E_VBM/CBM(E) 一一对应。

对带电或强极性中间体，频率和熵校正也不宜完全照搬零场结果。若资源有限，可先对零场、最强正场和最强负场各做一次频率，检查 ZPE 改变量是否远小于电子能差。当校正项变化接近目标 ΔG 差值时，只比较 TOTEN 会把自由能趋势放大或反向，也会让不同吸附物之间的排序失去热力学含义。

ΔG_X(E) = E_slab+X(E) − E_slab(E) − μ_X + ΔZPE(E) − TΔS(E)　　式(1)

带边移动和 ΔG 不一致时看什么？

如果 CBM 上移而 H* 自由能没有明显改善，先检查 H* 态是否真的与导带边耦合。吸附反应位可能由边缘金属 d 态、缺陷态或表面悬挂键控制，而不是由二维半导体的体相带边控制。带边给出电子供受能力边界，吸附自由能给出局域成键结果，两者相符才支持同一机制解释。

还要把电荷差分、Bader 电荷和局域势放在一起看。外场可增强界面极化，却不一定增加净电荷转移；也可能改变吸附键长，却把电子重新分布在远离反应位的一侧。若 ΔG 变化主要来自分子偶极取向，应把不同取向作为独立构型比较，而不是把它们混入同一条场强曲线。

另一个常见分歧来自半导体带隙误差。低成本扫描看趋势，高精度计算定参照，这是一种更清楚的分层安排。PBE 下的带边响应趋势可以用于内部比较，但绝对带边若要对接水氧化还原电位或实验电位，通常需要 HSE06、GW 近似或经验校正辅助；两层任务分开后，外场响应才不会被泛函误差掩盖。

图3展示 ReS₂ 在应变与电场工程下的光催化响应。对 VASP 数据而言，更有用的不是照搬材料结论，而是学习它把能带、带边位置和反应可行性放在同一组变量里扫描。

图3　ReS2 在应变与电场工程下的光催化响应。 DOI: 10.1039/d1ra03821d。

怎样做成可复核的 VASP 流程？

推荐用零场结构作为母体，先完成干净 slab 的场强扫描，再挑选两到三个代表性场强做吸附构型枚举。每个场强都要在同一套场强表里对应 POSCAR、INCAR、KPOINTS、OUTCAR、LOCPOT、vasprun.xml 和吸附态振动输入，否则后期很难发现方向或参照写反的问题。

图4　有向电场增强 MoS2 边缘的 DFT 电荷响应和 ΔGH*。 DOI: 10.1021/acsnano.4c06341。

INCAR 中常用 LDIPOL、IDIPOL 和 EFIELD 控制外场与偶极修正，具体写法要与 VASP 版本说明核对。二维 slab 通常把 IDIPOL 设在真空方向；如果坐标轴不是 z，就不能沿用别人的模板。结构弛豫时还要观察吸附物是否被场强拉出表面势阱，出现脱附趋势时应把该点标为不稳定构型，而不是强行读能量。

最后把结果画成同一横轴的四联图：场强、E_VBM/E_CBM、吸附侧 Bader 电荷、ΔG_X。材料筛选阶段应优先看四联同步图；若四条曲线在某个场强同时出现拐点，就检查结构重排或吸附态身份变化，若只有 ΔG 改变，说明外场主要影响局域成键。

图4和图5来自 oriented electric field 调控 MoS₂ 边缘 HER 的研究。图4把异质结边缘的电荷响应与 ΔG_H* 放在同一框架，图5再连接电化学表现；这正是同步判断的核心：外场改变电子结构，电子结构必须能落到吸附自由能与实验可观察量上。

图5　有向电场增强边缘位点后的 HER 表现。 DOI: 10.1021/acsnano.4c06341。

【本章要点总结】

外电场扫描先固定方向、吸附侧和偶极修正，不能只写正负场强。

带边要用真空平台或局域势对齐，VBM/CBM 本征值不能裸比较。

ΔG(E) 需要 clean、ads 和校正项同场强配对，构型取向也要重新枚举。

带边位移与吸附自由能一致时，才适合把外场效应解释为反应活性变化。