高斯Anharmonic计算：为什么谐振近似会误导你

做振动光谱计算的人，大多在某个阶段被谐振近似的结果搞失望过。算出来的红外频率，和实验值稳定地差5-10%，氢能键体系甚至能差到20%。不是你的计算水平不够，是谐振近似这个前提，在有些体系里根本站不住脚。

谐振近似到底忽略了什么

谐振近似把每个振动模式近似成简谐振子，势能面用二次函数描述。实际分子的势能面，在平衡位置附近确实近似抛物形，但一旦偏离平衡位置，非谐性就出来了——键的弯曲模式尤其明显。

非谐性的直接后果是：

• 基频和泛频的能量不是整数倍关系
• 费米共振（Fermi resonance）会导致谱峰分裂
• 热力学函数（熵、热容）在高温下偏差明显

Gaussian里做anharmonic计算，本质上是在谐振频率的基础上，用微扰论（VPT2，Vibrational Perturbation Theory 2nd order）修正非谐贡献。

Gaussian里的实际操作

关键关键词是Anharmonic。基本输入长这样：

# B3LYP/6-311+G(d,p) Opt Freq Anharmonic

有几个注意点：

第一，计算量会爆炸。谐振频率计算只需求 Hessian 矩阵（二阶导）。Anharmonic 需要三阶和四阶导信息，计算量是同基组下谐振计算的数十倍。对超过20个重原子的分子，基本不现实。

第二，需要足够大的基组。6-31G(d) 这个级别做 anharmonic 可靠性很差，至少要用 6-311+G(d,p) 或更大的基组，最好加上弥散函数。

第三，VPT2 的收敛性问题。有些体系的微扰级数不收敛，Gaussian 会报 VPT2 failed to converge 的警告。遇到这种情况，可以尝试 Anharmonic=VCI（振动组态相互作用），精度更高但更贵。

什么时候必须用Anharmonic

常见坑点

坑一：把Anharmonic频率直接当实验值比。Anharmonic修正虽然大幅改善精度，但仍有系统偏差——主要来自交换关联泛函误差和基组不完全。拿Anharmonic算出来的O-H拉伸频率跟实验红外光谱比，偏差一般在20-50 cm⁻¹，这属于正常水平。

坑二：忽略了溶剂效应。气相Anharmonic计算对气相实验光谱有价值，但液相光谱需要考虑溶剂。Gaussian里可以用SMD隐式溶剂加上Anharmonic，但计算量进一步放大，需要权衡。

坑三：热力学函数的温度外推。Anharmonic对热力学函数的修正主要体现在高温区。如果只关心298K的热力学数据，谐振近似+频率校正（quasi-harmonic approximation，用QuasiHarmonic关键词）通常够用，不必硬上完整的Anharmonic计算。

实用建议

对于日常计算，我的策略是分层处理：

1. 先跑谐振频率，看振动模式，确认没有虚频，基频位置和实验值的偏差在可接受范围（<5%）的话，直接出结果。
2. 对关键模式（比如涉及氢键的拉伸模式）做Anharmonic修正，不必全分子做，只针对你关心的振动自由度。
3. 如果要做完整的非谐光谱，考虑换用专门的程序（比如Molpro、CFOUR），它们对高阶导的处理比Gaussian更稳健。

热力学函数方面，QuasiHarmonic 近似（在Gaussian里用 Freq=QuasiHarmonic 调用）用经验公式修正低频模式的熵贡献，对大多数有机分子已经够用，计算量只比普通频率计算多一步后处理。

Anharmonic计算不是万能药。它的价值在于，让你知道谐振近似的误差来自哪里，而不是盲目相信一个精确到小数点后一位的频率数值。搞清楚误差来源，比追求绝对精度更有意义。