3.2 用频率分析验证水分子极小点

这些输入和起始模型为未执行的教学示例。原创配图为示意图，不是计算结果。使用前请验证版本相关语法、授权数据、数值收敛及科学模型。

3.2.1 模型、单位与来源记录

结构长度使用 Å，电子能使用 hartree，振动波数使用 cm⁻¹。核对各输出字段，区分热修正与标准态项。

3.2.2 未执行输入与明确修改片段

按照下方说明确认每个片段的母计算与替换位置；片段不自动构成独立输入。保留空行与文件来源要求。

3.2.2.1 输入块 1

%chk=case10_water_min.chk
%mem=4GB
%nprocshared=2
#p B3LYP/6-31G(d) Opt=(Tight,CalcFC) Freq
   SCF=Tight Integral=UltraFine

Water optimization and minimum verification

0 1
O   0.000000   0.000000   0.000000
H   0.758000   0.000000   0.586000
H  -0.758000   0.000000   0.586000

3.2.3 实作研究

3.2.3.1 问题与直观理解

“优化完成”是否就证明结构稳定？它说明优化算法的收敛标准已经满足，但尚未判断所有振动方向上的局部曲率是否都向上。频率分析检查的正是局部曲率：碗底在各方向向上弯曲，而山口至少有一个下坡方向。对于孤立、非线性的三原子水分子，振动子空间包含 3N−6 = 3 个模式；整分子的平移和转动不是额外的内部振动。

本例只要求能够阅读 Gaussian 输出和笛卡尔坐标。几何优化与 Hessian 计算明确采用相同的方法、基组、DFT 积分网格及 SCF 阈值。用另一种模型在原几何上计算 Hessian，通常不能证明该几何是新模型的驻点。本例对象是中性闭壳层单个水分子，不是自由基、离子、水团簇或液态水样品。

3.2.3.2 完整输入

参见上方输入块 1。

Opt Freq 在几何优化之后计算 Hessian。最终检查点为 case10_water_min.chk，是案例 11、12 的明确前提文件；确认两阶段均完成后应保留。文件名含有 min 或 freq 并不能证明其内容正确，还应同时保存对应日志，因为模型、收敛状态、频率和终止状态等验证信息都在日志里。

3.2.3.3 操作步骤

检查初始分子并运行完整输入，在作业记录中保存 Gaussian 版本、route 行、总电荷、多重度和资源设置。
找到最终优化收敛表，检查最大力、均方根力、最大位移及均方根位移，并结合优化完成信息判断。SCF 能量低不能替代几何收敛。
阅读实际振动频率区块，不要只看最前面的低频诊断行。数出三个内部模式，记录 cm−1 单位的频率和 km mol−1 单位的红外强度，并检查是否有负值振动频率。
在查看器中播放每个模式，根据位移矢量辨认弯曲、对称伸缩和反对称伸缩；实际运动模式比背诵大致频段更可靠。
检查最终终止行，确认热化学数据确实已输出。所有检查通过后再归档检查点；若仍存在有意义的虚频，可沿该模式正反方向做小位移、重新优化并再次计算频率。

3.2.3.4 输出与验收标准

Gaussian 通常以频率列表中的负数表示虚频。它意味着该方向的曲率为负，而不是分子在以“负振动速率”做物理振动。本例中合格的极小点应有三个实的内部振动模式、可忽略的残余梯度和正常完成状态。这证明的是所选模型下的局部极小点，不是所有成键方式、电荷态或电子态中的全局最低点。

在较大且柔软的分子中，很小的负频率需要结合数值精度、位移模式及真实收敛情况分析。不能设定一个通用阈值，把所有较小虚频直接当成正频率。收紧优化和积分网格并重算可以提供证据，手动去掉负号则不能。水分子较刚性，本例出现内部负频率时应立即排查，而非自动套用低频例外。

谐振频率是在一个平衡结构附近得到的模型预测。与实验比较时，可能需要有出处的频率缩放因子、非谐处理，并匹配相态及同位素组成。不能为了让某条峰吻合而随意选择缩放因子；缩放后的光谱也不能修复不稳定的几何。应保留未缩放频率以供核查。