沿着负频的振动方向调整一下原子位置不就行了？这个应该不难吧。

避免freq最后出现NO乃至虚频，或者说让优化更精确的做法有二：
(1)几何优化时用更严的收敛限，即opt里面写上tight。对于DFT来讲，建议同时搭配int=ultrafine使用更高精度的积分格点，效果会明显更好（注：从G16开始int=ultrafine已成为默认）。这样来解决出现NO和虚频对计算耗时增加不太多，不过也有可能会反复震荡很难达到这个收敛限。
注意：如果写opt=tight freq，则freq最后也会用tight标准来判断是否收敛。实际上只要能满足默认判断标准即可，不需要非得满足tight的判断标准。

(2)使用精确的Hessian矩阵。在opt里写上calcall，那么几何优化过程中每一步都会精确计算Hessian，此时几何优化和freq所用的Hessian都是一致的，因此对收敛的判断结果也是完全相同的。也就是说，只要opt最后是YES则freq最后也必然是YES，并且多数情况也可以确保无虚频。（很值得一提的是，用calcall的时候，优化任务最后会自动做振动分析，也就是说，此时实际上完全没有必要再单独用freq关键词做振动分析了。）
不过，每一步优化都精确计算Hessian太耗时。如果你之前可以已经在默认情况下优化过，可以取最后的结构重新优化，并在opt里写上calcfc，这样只在优化第一步的时候精确计算Hessian，而之后还是近似计算Hessian。这样收敛后也往往可以避免freq之后出现NO。但是如果初始结构离实际极小点太远则这么做无效，因为等到收敛时Hessian矩阵可能又偏离精确Hessian比较远了。

参考http://sobereva.com/278，