评估不同模型化学方法下谐振振动频率的幂律缩放因子

《Journal of Chemical Theory and Computation》：Evaluating Power-Law Scaling Factors for Harmonic Vibrational Frequencies across Model Chemistries

　　

利用最近发展的波数幂律缩放（WPS）方法，推导了基频谐振振动频率和分子热力学校正的缩放因子。WPS方法实现了对振动频率的简单双参数缩放，对于所测试的最佳方法，其误差相比均匀缩放最多改善约64%。研究人员在超过500种模型化学方法上计算了幂律缩放因子和均匀缩放因子，包括波函数理论（HF、MP2、CCSD）以及采用多种原子中心基组的密度泛函理论。实验对比使用包含894个振动模式的数据集，其中已移除C≡C和C≡N三键弯曲模式。B97、B97–1和τ-HCTHh密度泛函方法的计算结果与实验基频非常接近，建议普遍使用。WPS频率随基组尺寸增大的收敛速度快于均匀缩放频率，因此兼具更高的准确性和计算效率。幂律缩放使用统一的缩放因子描述基频振动、焓和熵，并可改善对均匀缩放不敏感的广义梯度近似频率。

该研究发表于《Journal of Chemical Theory and Computation》。量子化学计算的谐振振动频率会因非谐效应、基组不完备和电子相关描述不完整而偏离实验值。Pople等人提出均匀缩放因子进行经验校正，但谐波偏差在全谱范围内并不恒定，已有修正方法可能引入不连续、多项式拐点或过拟合问题。最近发展的波数幂律缩放（WPS）方法采用简单双参数幂律曲线描述缩放因子随波数的单调递减行为，可在整个光谱范围内一致处理非谐性。本研究的目的是为基频振动频率、振动焓和熵提供一套统一的WPS缩放因子，并在大量模型化学方法上评估其适用性。

研究人员使用Q-Chem 6开发版计算优化几何和振动频率；电子结构方法包括Hartree-Fock（HF）、二阶M?ller-Plesset微扰理论（MP2）、耦合簇单双激发（CCSD），以及广义梯度近似（GGA）、meta-GGA、全局杂化GGA、meta-杂化GGA和范围分离杂化GGA等密度泛函理论（DFT）泛函，结合Pople和Dunning系列等原子中心基组。实验对比采用Hanson-Heine修改后的F1频率集，含894个振动，并去除C≡C和C≡N三键弯曲模式。缩放因子通过最小化对数-对数残差拟合得到，并使用Microsoft Excel的LINEST函数和Python库完成交叉验证。

**基频振动**：B97和B97–1杂化泛函在6–311+G(2df,p)基组下WPS RMSD分别为16.1和16.2 cm^–1，优于均匀缩放的25.3和25.0 cm^–1；τ-HCTHh在6–311++G(2df,p)下进一步降至15.8 cm^–1。HF误差近似均匀，WPS无法改进；加入MP2和CCSD相关校正逐步改善。GGA泛函对均匀缩放不敏感，但WPS显著改进，KT3是表现最好的GGA。

**基组测试**：对最佳方法扩展基组测试，B97、B97–1和τ-HCTHh在Def2-TZVPD基组下WPS RMSD分别为15.3、15.3和15.4 cm^–1。WPS频率随基组增大收敛更快，例如B97/6–31+G(d,p)与最大Pople基组仅差约0.7 cm^–1，而均匀缩放差约3.6 cm^–1。最小基组6–31G(d)的均匀缩放误差抵消优势不再适用于WPS，但WPS仍优于均匀缩放。

**热力学性质**：298.15 K下，面向热力学焓、熵专门优化的均匀缩放因子使焓RMSD降低约0.010 kJ/mol、熵降低约0.07–0.08 J/mol/K；而WPS无需专门优化即可给出更优结果，B97/Def2-TZVPD焓RMSD为0.009 kJ/mol、熵RMSD为0.05 J/mol/K。进一步优化WPS参数仅带来极小改进，因此可用同一套缩放因子统一描述基频、焓和熵。

**大分子测试**：对参数集外的萘和苯甲腈测试，苯甲腈在τ-HCTHh/Def2-TZVPD下WPS RMSD比均匀缩放低约34%，去除含C≡N三键弯曲贡献的振动后改进约40%。萘中均匀缩放在低于1000 cm^–1区域略优，WPS在高波数伸缩模式更优；去除第9个弯曲异常模式后，WPS在全部准确度指标上一致改进。

**交叉验证**：将修改后的F1分子集按9:1、8:2、7:3比例随机拆分，重复100次。训练集平均RMSD在15.2–15.3 cm^–1，测试集在14.8–15.5 cm^–1；缩放因子仅在λ₁第三位、λ₂第四位小数变化，表明WPS对数据集缩减稳健，并非过拟合。

**总结与结论**：讨论指出，WPS能在几乎所有情况下显著优于均匀缩放，且对均匀缩放不敏感的GGA频率同样有效。功率缩放因子随基组收敛更快，并能用统一因子同时处理基频、振动焓和振动熵。研究结论为：B97/Def2-TZVPD表现最佳，RMSD 15.3 cm^–1，中位绝对偏差（MED）6.9 cm^–1，较均匀缩放分别改善约40%和64%；推荐B97、B97–1和τ-HCTHh配合Def2-TZVPD基组用于振动性质计算。