吸附态频率计算结束后，OUTCAR 或 vasprun.xml 里偶尔会出现一个很小的虚频。小虚频本身只说明当前 Hessian 在某个方向上有负曲率，是否影响吸附态自由能要看这个方向对应什么运动。很多人会立刻把结果删掉，也有人直接把虚频当作 0 处理并继续做 ΔG，两种处理都可能把问题简化过度。
VASP 吸附态频率中的小虚频，关键在于它对应数值噪声、表面软模、吸附物摆动，还是已经落到真实反应方向。活动原子范围、ZPE/振动熵处理和扰动后重新优化结果，会共同决定这个吸附态还能不能作为局域极小点参与后续热力学校正。
频率计算给出的不是一个抽象数字，而是一组归一化振动模式。判断小虚频要先看振动模式动画，再看频率绝对值。若虚频对应整个 slab 的平移、底层原子的集体漂移或真空方向轻微摆动，它往往和固定层、力收敛、数值精度或低频软模有关；若虚频对应 O-H 断裂、C-S 断裂、H 迁移或吸附物离位，则可能不是稳定吸附态。
Fe(110) 上 *OH 在不同表面位点的优化构型见图1。对于这类吸附物，真实稳定态应在局域位点附近形成正定曲率；若小虚频沿着 Fe-O 键伸缩或 O-H 断裂方向，后续 ΔG 修正不能直接沿用。若虚频只是 *OH 绕吸附键的受限转动，且重新优化后回到同一构型，则可按低频处理方案评估其对熵项的影响。
图1 Fe(110) 表面不同位点的 *OH 吸附优化结构。图中吸附位点差异说明虚频方向需要和局域构型一起判断。 DOI: 10.1021/acsomega.6c05647。
没有一个对所有体系都成立的固定阈值。阈值只能作为筛查入口，不能替代振动模式和力收敛检查。金属 slab、大真空层、低覆盖度吸附态和柔性有机分子常见几十 cm⁻¹ 以内的低频扰动；强成键小分子、过渡态附近结构或未充分优化的吸附态，即使虚频数值不大，也可能指向真实反应坐标。
若要判断该虚频能否保留，力收敛、模式方向和扰动构型要一起看。优化后的最大残余力先达到设定标准；虚频模式再区分吸附物运动、表面原子软模或整个模型漂移；随后沿虚频方向正反向轻微扰动并重新优化。若两侧都回到同一吸附构型，说明它更可能是数值或受限低频问题；若滑向另一个位点、断键或形成新中间体，当前结构不应作为稳定吸附态参与自由能修正。
处理这个频率结果时，不只写“有一个小虚频”。小虚频的处理结论必须来自模式证据，而不是来自一句口头经验。虚频数值、位移方向、受影响原子、重新扰动后的结构去向，以及重新优化前后的最大力，都要一起决定该模式属于可修正的低频扰动，还是说明结构本身需要重算。
吸附态频率通常不让整块 slab 都参与有限差分，因为底层原子的低频模式会显著增加计算量，也会把整体平移或表面声子混入吸附物热校正。底层原子通常保持固定，频率计算只覆盖吸附物和邻近表层原子。这样得到的是局域 Hessian，适合估计吸附态 ZPE 和局域振动熵，但不等同于完整表面声子谱。
层状 FeS 结构和不同表面见图2。对于层状或低维表面，表面原子本身存在柔软的剪切、弯曲和层间位移模式；活动原子范围若过大，低频模式会明显增多，熵项可能被放大。若活动原子范围过小，吸附物和表面协同振动又会被截断。活动原子范围决定 Hessian 维度，也决定小虚频来自哪一类运动。
图2 层状 FeS 的晶体结构。层状表面的柔性位移会影响吸附态频率中低频和小虚频的来源。 DOI: 10.1021/acs.jpcc.7b08711。
频率计算前的优化应使用比普通筛选更严格的力收敛，例如 EDIFFG 按体系测试到 −0.02 eV Å⁻¹ 或更严；静态和频率步骤保持相同 ENCUT、k 点、赝势、ISPIN、U 值和 vdW 设置。IBRION=5 或 6、NFREE=2、POTIM 取 0.01–0.03 Å 作为起点，具体位移幅度要看轻元素和表面柔性。
若小虚频只在一组 POTIM 或一个 k 点下出现，优先检查有限差分步长、电子收敛和力噪声。小虚频是否可接受，取决于模式稳定性，而不是单次输出是否好看。若改变 POTIM、提高 EDIFF、增加 FFT 精度或重新优化后虚频消失，原结果不适合作为最终自由能输入；若多组设置都保留同一方向的虚频，则要把它看作结构或路径问题。
对于含 H、O、N 等轻元素的吸附物，POTIM 过大可能让有限差分跨过局域极小点附近的谐振区，POTIM 过小又可能被力噪声淹没。可用两组位移幅度做对照，并保持同一活动原子列表。若频率符号随位移幅度反复改变，优先回到结构优化和电子收敛，而不是直接把该频率纳入热校正。
ZPE 来自所有实频振动的 1/2hν 求和。若只有一个很小虚频，且模式对应整体漂移或受限转动，许多工作会把该低频模式从谐振子 ZPE/熵中剔除，或改用受限平动/转动近似；但这个处理必须在同一批吸附态中统一执行。若虚频对应键断裂或迁移方向，则 ZPE 没有稳定极小点基础。
FeS 表面噻吩脱硫的反应能垒路径见图3。吸附态小虚频若指向反应坐标，说明当前结构更接近鞍点或未收敛构型，不能直接进入 ΔG 计算。过渡态需要且只应有一个沿反应坐标的虚频；稳定吸附态则应没有指向反应坐标的虚频。
图3 FeS 表面噻吩直接脱硫反应路径。图中能垒路径说明稳定吸附态与过渡态的频率判据不同。 DOI: 10.1021/acs.jpcc.7b08711。
低频模式对振动熵很敏感，因为谐振子公式中低 ν 会放大热熵贡献。一个 20–40 cm⁻¹ 的摆动模式，对 ZPE 的影响很小，对 −TΔS 的影响却可能明显。若不同吸附态低频处理不一致，自由能排序会被熵项牵动，尤其是大分子、显式水层和弱吸附体系。
处理方式应在同一套算例中保持一致：对于所有吸附态，统一设定低频截断或受限平动近似；频率修正的可比性来自统一规则，不是把单个异常模式单独修掉。对于气相分子，保留标准平动、转动和振动热化学；对于 slab，通常只把活动区域的局域频率用于吸附态修正。
Ru(0001) 的表面声子色散和声子态密度见图4。金属表面和低维材料本来就存在低频声子，局域吸附态频率计算若把较多表面原子纳入活动区域，容易出现靠近 0 的模式。这类模式不一定代表吸附物不稳定，但固定层、超胞尺寸和整体漂移约束仍要重新核对。
图4 Ru(0001) 表面声子色散和声子态密度。图中低频表面模式说明吸附态频率需要区分局域振动和表面软模。 DOI: 10.1021/acs.jpclett.7b00784。
若虚频位移主要分布在远离吸附物的表层原子，可先重新优化 slab、增加固定层约束或缩小频率活动区。若位移集中在吸附物与活性位之间，则要沿模式方向扰动并重新优化。保留结果的前提是结构回到同一局域极小点，且低频处理方式对关键 ΔG 排序没有改变。
CH₄/Pt(111) 吸附能随 Pt-C 距离变化的曲线见图5。稳定吸附态应位于势能面极小点附近；若频率中出现沿脱附方向的虚频，说明当前距离附近不是局域极小点，需要重新优化吸附高度或检查 vdW、k 点和固定层设置。势能曲线或离散位移扫描能帮助判断虚频是否来自真实脱附方向。
图5 CH₄/Pt(111) 吸附能随 Pt-C 距离变化的曲线。势能面极小点可用于辅助判断脱附方向虚频是否来自未充分优化结构。 DOI: 10.1021/acs.jctc.3c01308。
数值噪声型小虚频对应精度提高后的频率更新；表面软模型小虚频对应统一的低频处理和活动原子范围说明；反应坐标型虚频则要回到结构优化或 NEB，不能把它当稳定吸附态自由能。只有确认结构是局域极小点，吸附态自由能才有热力学校正基础。
要点：小虚频要结合振动模式、残余力和扰动后优化结果判断，不能只按数值大小处理。
要点：活动原子范围会改变 Hessian 维度，吸附物局域振动、表面软模和整体漂移要分开看。
要点：ZPE 对小低频不太敏感，振动熵更容易被低频模式放大，因此同一批算例要统一低频规则。
要点：若虚频指向断键、迁移或脱附方向，当前结构不能直接作为稳定吸附态参与自由能修正。
