求助Gaussian16A03跑过渡态用calcall收敛后fre计算出现NO且没虚频

djjj148

各位老师好，请教一下，我用Gaussian16A03跑过渡态。普通有机分子吸附氧气分子的模型，关键词如下（完整输入输出文件已上传）：%chk=TS-adsO2.chk
# b3lyp/6-31g(d) opt=(ts,calcall,noeigen) fre

Title Card Required

0 2

N                    -1.89954  -1.21148  -1.04412
N                    -0.88402  -4.38165   1.39554

.....（略）
得到的log文件显示过渡态优化收敛了
               Item               Value     Threshold  Converged?
Maximum Force           0.000399     0.000450     YES
RMS        Force            0.000037     0.000300     YES
Maximum Displacement   0.000656     0.001800  YES
RMS     Displacement     0.000146     0.001200    YES
Predicted change in Energy= 4.059102D-07
Optimization completed.
    -- Stationary point found.

但是接下来的fre计算却出现了2个NO！
         Item                         Value     Threshold  Converged?
Maximum Force                 0.000399     0.000450     YES
RMS     Force                    0.000037     0.000300     YES
Maximum Displacement      0.743684     0.001800     NO
RMS     Displacement         0.183600     0.001200     NO
Predicted change in Energy= 1.118061D-03

本人新手，不太理解为何已经用了calcall关键词收敛了，频率计算还会得到这样的结果。恳请各位老师指点一下，谢谢！

djjj148

这是完整的输入输出文件

zjxitcc

都是老生常谈了···花1 min搜一下就自己获得答案了···看sob老师博文最后一段《Gaussian中几何优化收敛后Freq时出现NO或虚频的原因和解决方法》http://sobereva.com/278

PS1：建的分子看上去有点怪怪的，建议确认其合理性····
PS2：是freq不是fre。另外，写了Callcall没必要再写freq，Callcall最后会自动算freq，写了相当于多浪费了些计算时间。

djjj148

zjxitcc 发表于 2020-1-21 10:46
都是老生常谈了···花1 min搜一下就自己获得答案了···看sob老师博文最后一段《Gaussian中几何优化收敛 ...

感谢您的回复。我的确把freq写错了~ 不过我应该是没写清楚，导致您有点误解，再麻烦您一下，如果您有空
PS1：freq的确是写错了，但是我之前是只写了# opt=(ts,calcall,noeigen) b3lyp/6-31g(d)关键词且g16提示收敛了的，而且把输出文件放到gview里看不到一个虚频。
PS2：我之前sob老师的那篇博文我很早就看了，但是里面最后提到的“可能原因有三 (1)优化时用的级别、数值设定（如DFT积分格点）等因素和几何优化时不完全一致 (2)用gview读取优化的输出文件，然后又保存成freq任务的输入文件的时候，由于小数位数有限，因此造成了一点数值误差 (3)几何优化的时候用了GDIIS”，我觉得我的情况都不是那三种，1）我写calcall算完过渡态应该是直接算freq的吧，应该不存在优化的级别和freq的级别不一样的情况。2）我并没有用gview转输出文件。3）我也并没有用gdiis方法算。

djjj148

附上只写了calcall且没写freq的输入输出文件，结果显示过渡态收敛但是没有一个虚频，还清各位老师有空帮看下，谢谢

zjxitcc

djjj148 发表于 2020-1-21 11:46
附上只写了calcall且没写freq的输入输出文件，结果显示过渡态收敛但是没有一个虚频，还清各位老师有空帮看 ...

这次是想问为何没虚频？优化算法并不能 百分百保证 得到想要的结果/结构，想找过渡态却没虚频，说明初始结构不够合理，不够接近真实过渡态。换个更好的结构初猜即可。

PS：看不出你想找啥过渡态，O2与大分子成键的过渡态？O2与大分子不成键、只靠分子间作用力维持一定距离的过渡态？（后者不一定存在）

zjxitcc

djjj148 发表于 2020-1-21 11:42
感谢您的回复。我的确把freq写错了~ 不过我应该是没写清楚，导致您有点误解，再麻烦您一下，如果您有空
...

博文中第（3）点括号内容

djjj148

zjxitcc 发表于 2020-1-21 12:27
这次是想问为何没虚频？优化算法并不能 百分百保证 得到想要的结果/结构，想找过渡态却没虚频，说明初始 ...

我明白了！谢谢您的回复！

Machaophy

可以做势能面扫描，然后在扫描中找到一个合适的过渡态初猜
会靠谱很多

djjj148

Machaophy 发表于 2020-1-25 08:38
可以做势能面扫描，然后在扫描中找到一个合适的过渡态初猜
会靠谱很多

谢谢！