显式水 AIMD 的主要作用，是记录界面水、吸附物和表面位点在有限温度下的共同演化。将轨迹中的构型用于静态吸附计算前，需要确认水层分布已进入可取样区间，目标中间体在该时间段内保持相同成键关系，抽取后的结构在静态优化中仍对应同一化学对象。
界面水轨迹的第一层判断是能量、温度和关键几何量是否进入平台区。温度围绕目标值波动并不等于结构已经平衡；吸附物高度、金属-吸附物键长、水层质心和氢键数量也要随时间检查。可用于静态吸附的构型应来自平衡后的轨迹片段，而不是升温、重排或质子转移正在发生的片段。
以 Fe–N₄/C 界面水为例，图1显示 *OOH、*O 和 *OH 周围水分子会直接参与质子转移和氢键重排。如果轨迹中 *OOH 已经解离，后续静态计算就不能继续把该构型命名为 *OOH 吸附态；它对应的是新的 *O 与溶液 OH⁻ 环境。AIMD 先回答界面实际停留在哪类状态，静态优化只是在这些状态附近找局部极小点。
图1 Fe–N₄/C 等单原子模型上界面水诱导中间体转化的 AIMD 快照。图中展示了 *OOH、*O 和 *OH 附近的水分子重排。 DOI: 10.1038/s41467-022-29357-7。
轨迹诊断至少保留三个时间序列：Etot 或保守量、目标吸附键长 dM–X、以及水层相对表面的 z 向分布。若 dM–X 长时间漂移，或水层整体向真空侧迁移，说明模型尚未进入可取样区间。若只在一个瞬间出现有利氢键，后续静态优化很可能回到另一类构型，这种构型不适合作为代表结构。
轨迹长度还要和问题尺度匹配。只判断局域水分子是否能稳定围绕 *OH，几皮秒的平衡片段可能提供初步线索；若要比较两种水网络的出现概率，单条短轨迹就不够。可以从不同初始水层排布启动两到三条独立 AIMD，检查它们是否收敛到相近的水密度和氢键统计。多个初始水层给出相近统计时，构型选择才更有代表性。
MSD(t)=O(t)-rO(0)|²> 式(1)
式(1)可用于水氧原子的均方位移。若研究对象是局域氢键网络，MSD 不需要无限增大才算有效；它更适合判断水层是否发生整体扩散、交换或脱离界面。对于 20-30 ps 的短轨迹，MSD、密度和 RDF 应联合使用，不能单凭一条曲线断定界面结构已经充分采样。
构型选择不应只按固定时间间隔进行。可先根据几个可计算坐标分组：吸附物高度 R、吸附物与水的氢键数、第一层水氧密度峰位置、以及目标质子的转移坐标。图2把 OH⁻ 相对 Fe 位点的 z 坐标和氧的氢配位数放在一起，展示了轨迹中可区分的状态区域。构型集合应覆盖概率高且物理含义清楚的状态，而不是只挑自由能较低或图像整齐的瞬间。
每个候选状态至少选取 3-5 个相隔足够远的轨迹构型，再分别做静态优化。如果这些构型优化后回到同一局部极小点，说明该水环境可用一个静态构型近似；如果优化后分裂成多类构型，就要把它们作为不同微观状态处理。对 *OH、*OOH、*COOH 这类极性中间体，氢键方向和供受体身份比水分子数量更重要。
选择轨迹构型时还要避开周期边界带来的假近邻。水分子可能从盒子顶部跨周期靠近表面，视觉上看似形成氢键，实际距离需要按最小镜像重新计算。对于倾斜晶胞或非正交表面，先用未包裹坐标重构连续轨迹，再统计 O···H 距离和角度。构型坐标必须与周期边界处理一致，否则静态结构会继承错误的水分子相对位置。
Pi = Ni / Nsample， = Σ PiAi 式(2)
式(2)给出离散状态的占比和加权平均。Pi 来自平衡轨迹中的样本数，Ai 可以是吸附能、键长或 Bader 电荷等静态优化后得到的量。若轨迹很短或状态间转换不足，Pi 只能作为构型出现频率，不应写成严格热力学权重。
图2 准吸附 OH⁻ 物种在 AIMD 轨迹中的空间演化和概率分布。R 与配位数用于区分可优化的局部状态。 DOI: 10.1038/s41467-022-29357-7。
静态吸附能要比较同一水层背景下的吸附态和参照态。若吸附态含 12 个水分子，干净表面参照却不含水，Eads 会混入建水层的能量代价。更合理的做法是保留相同 slab、相近水层数量和同一晶胞，在抽取的水环境中分别构造吸附态与对应参照。显式水吸附能的可比性来自相同界面背景。
图3给出水密度和 OH⁻ 概率分布，说明界面水存在稳定分层和局部极小点。静态优化时，第一层水和关键氢键应保留，远离活性位、只提供体相水背景的分子可以按研究目标处理。若优化后水层坍缩、目标氢键断裂或质子自动转移，应把该结果归入新的界面状态，而不是强行维持原命名。
图3 界面水密度、OH⁻ 概率分布和近似自由能剖面。局部极小值可作为抽取静态构型的候选。 DOI: 10.1038/s41467-022-29357-7。
吸附能可按下式记录，其中 H₂O 层必须在左右两侧口径一致：
Eadswater = Eslab+water+X – Eslab+water – EX,ref 式(3)
式(3)中的 X 可以是 *OH、*OOH 或 *COOH 的来源参照，EX,ref 按对应反应式选取。若水层在吸附态和参照态发生不同重排，应同时报告结构差异，例如第一层水数目、关键 O···H 距离和氢键网络。水层重排本身是结果的一部分，不能只把它视为静态计算中的背景噪声。
参照态也要从同一水环境中构造。若从吸附态删除 X 后水层立即塌陷，得到的 Eslab+water 就不再对应原吸附环境；可以保留同一时间点的水层构型做单点参照，再另行优化水层评估松弛贡献。这样能区分吸附相互作用、界面水重排和参照态弛豫三部分，而不是把它们压缩进一个总吸附能。
AIMD 后静态优化之前，先检查吸附物的成键拓扑。*OOH 若已断裂为 *O 与溶液 OH⁻，它不再是 OOH 吸附态；*OH 若通过水链获得额外质子，则可能接近 H₂O 脱附或准吸附 OH⁻ 状态。图4中的 RDF 和代表性快照显示，背侧水分子也会改变金属位点的局域配位环境。目标物种归属必须由键长、配位数和质子位置共同确定。
图4 背侧轴向水的径向分布函数和代表性 AIMD 快照。RDF 可用于判断界面水是否形成稳定局域结构。 DOI: 10.1038/s41467-022-29357-7。
静态优化后还要做回读。若初始构型来自轨迹中的局部极小点，优化后目标距离 R、配位数和氢键模式不应偏离该状态太多；若偏离明显，说明该构型只是热涨落中的过渡结构。对于质子转移问题，可以沿着 O–H 距离差 ξ = ddonor–H – dacceptor–H 检查状态身份。ξ 改变符号时，反应对象通常已经转换。
图5中的能量分布说明，不同金属位点形成准吸附 OH⁻ 的倾向不同。把这些状态带回静态吸附计算时，应按金属位点、氧配位、氢键和质子位置分组，而不是只比较一个最低能构型。静态构型应代表轨迹中可辨认的状态族，这样后续吸附自由能才有清楚的结构来源。
图5 不同金属单原子位点形成准吸附 OH⁻ 状态的能量分布。静态吸附判断应与轨迹中真实出现的状态对应。 DOI: 10.1038/s41467-022-29357-7。
如果目标是继续做 NEB，轨迹构型选择标准还要更严格。初态和终态应来自同一水层状态族，旁观水分子的数量、位置和氢键角色尽量一致；否则 NEB 会同时描述反应事件和水层重新排布。对于质子转移路径，可以先固定几个关键水分子位置做路径预筛，再释放水层做端点回算，确认势垒变化是否来自目标反应坐标。
目录可分为 `01_relax_dry`、`02_build_water`、`03_aimd_equil`、`04_aimd_sample`、`05_static_snapshots` 和 `06_compare_adsorption`。`03_aimd_equil` 保存升温或平衡段，`04_aimd_sample` 只放用于构型选择的平衡轨迹，`05_static_snapshots` 按状态编号保存 POSCAR、OUTCAR 和 CONTCAR。这样的归档能把热轨迹、构型选择依据和静态结果分开。
AIMD 参数应写成示例边界：NVT 可采用 300 K 或目标反应温度，POTIM 常取 0.5-1.0 fs 量级，NSW 由所需采样时间决定，氢参与转移时需检查时间步是否过大。电子步可使用 Γ 点或较稀 k 点取样，但被选中的静态构型要用正式 k 点和能量截断重新计算。
最终交付给后续吸附或 NEB 的结构，应附带四项说明：来自哪段平衡轨迹、属于哪个水构型状态、静态优化后是否保留原有成键关系、以及与无水或单水模型相比的能量差。若这些信息缺失，显式水模型只会增加自由度，并不能提高结论可信度。
静态优化后的构型可再进行短时间 AIMD 复核：把结构重新放入目标温度，观察关键氢键和吸附物成键关系是否能保留若干皮秒。若结构一升温就发生质子转移、脱附或水层逃逸，说明该静态极小点对热涨落不够稳定；若多个候选结构都能短时保留，则应结合轨迹出现频率和静态自由能排序。静态吸附构型需要同时经受优化稳定性和短程热稳定性检验。
静态能量读取也要统一：所有候选构型使用相同 ENCUT、k 点、展宽、偶极修正和电子收敛标准，能量取同一类最终静态输出。这样比较得到的是水构型差异，而不是后处理设置差异。
要点：AIMD 构型只有来自平衡轨迹、且吸附物身份清楚时，才适合进入静态吸附计算。
要点：构型选择应依据水密度、RDF、氢键和关键反应坐标，而不是只按固定时间间隔。
要点：显式水吸附能必须保持界面背景一致，并把水层重排作为结果记录。
