方案总结

下面盘点从不同网站中摘抄出的解决方法，并对相关情况进行总结。

方案1：Gaussian计算消除虚频实用技巧

地址：http://www.goscience.cn/article/articleInfo?type=2&id=331

1. 降低对称性： 例如用D3h对称性优化NH3的时候会得到一个平面分子，会有一个很大的虚频，垂直于分子平面，使用C3v对称性即可消除虚频，此降低对称性操作在Gaussview中即可完成。或者直接在路径行加入命令nosymm，就可以搞定很多出现小虚频的情况。当然，用Gaussian降低对称性的内部高级指令IOp(2/15＝3)也行；
2. 降低结构的收敛精度，opt＝loose；
3. 使用opt=tight或者verytight，提高收敛精度；
4. 将虚频对应的振动坐标（严格地应称之为笛卡尔简正位移坐标）的1/2（或者其他因子，具体需要经验）叠加到原始坐标（优化好结构的笛卡尔坐标）。很多时候此招甚灵，我将其称之为“头疼医头，脚疼医脚”；
5. 在GaussView里沿虚频振动方向（这在Gaussview是可视的）有针对性地调整初始结构，重新优化；
6. 计算力常数(opt=calcfc)，重新开始优化。

评论：目前尝试了3，添加opt=tight，没啥效果，可见整体性虚频不适合这么整；同样试了方法6，也没啥效果；看很多评论都说nosymm没啥用，但是很多又写了，我到底要不要试一下呢？

方案2：小木虫转帖

地址：http://muchong.com/t-4150220-1

查不到原地址的垃圾中文互联网……

在进行Gaussian计算时，经常会出现多个虚频个，每个都只有－30~－20左右，但是如果结构改动过大，则一不小心，虚频就会更大。消除虚频的主要方法是改变构型。其次在计算上还可以尝试：nosymm；加大循环次数；提高收敛度；iop(1/8=1)等。通常对于势能面很平，较小的虚频很难消除。
为了解决以上的这种现象,按如下的步骤一步步调整：

1. 在优化时采用Scf(tight)的选项，增加收敛的标准。再去计算频率。如果还有虚频，参见下一步。
2. 对称性的影响，很多情况下的虚频是由于分子本身的对称性造成的。这样，在优化时，如果必要，要将对称性降低，还有，输入文件有时是用内坐标。建议如果有虚频的话，将内坐标改成直角坐标优化。但是有时这一方法也并没有多大作用。
3. 如果上述方法还有虚频，看一下虚频，找到强度较大的，将在频率中产生的原子的振动坐标加到相应的输入文件中。这样，重新计算。直到虚频没有。比较好的做法是分子结构输入用直角坐标，将这些产生虚频的振动坐标直接与原始坐标相加。加的时候不一定要将坐标直接相加，可以选择权重如1/3相加等。
4. 实际上，如果分子柔性较大，很难找到最低点，这是电子结构计算的问题，这种情况下，需要动力学的东西，用构象搜寻的办法解决。如：模拟退火，最陡下降法，淬火法等。将得到的能量最低的构象做一般的电子结构计算，这样，应当没有问题。

如果上面几步仍然无法解决问题，则有可能是因为计算还没有得到最稳定的结构，那么或者是分子结构有问题，或者是计算错了，再者就是游离出现在计算的范畴。
最后如果是寻找过渡态，如果只有一个虚频，而且该虚频对应某单键旋转震动，那么产生虚频的原因就是初始结构比较接近重叠构象。这时把优化过的结构绕振动方向旋转几十度再重新优化就有可能得到正确的结果。优化中等大小的有机分子（50-100原子）的时候，应该经常碰到这种情况。
当出现两个以上虚频时，可以先看一下振动模式。如果多余的虚频，和想要的振动模式发生在不同原子上，可以先固定过渡态模式所对应的原子，然后做opt到minimal,一般可以把多余的虚频去掉，再放开做全优化。或许会有所收获。