计算化学公社

标题: 求助Gaussian16A03跑过渡态用calcall收敛后fre计算出现NO且没虚频 [打印本页]

作者
Author: djjj148 时间: 2020-1-21 09:06
标题: 求助Gaussian16A03跑过渡态用calcall收敛后fre计算出现NO且没虚频
各位老师好，请教一下，我用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关键词收敛了，频率计算还会得到这样的结果。恳请各位老师指点一下，谢谢！

作者
Author: djjj148 时间: 2020-1-21 09:06
这是完整的输入输出文件

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

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

作者
Author: djjj148 时间: 2020-1-21 11:42

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方法算。

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

作者
Author: zjxitcc 时间: 2020-1-21 12:27
本帖最后由 zjxitcc 于 2020-1-21 12:30 编辑

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

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

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

作者
Author: zjxitcc 时间: 2020-1-21 12:31

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

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

作者
Author: djjj148 时间: 2020-1-21 13:17

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

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

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

作者
Author: djjj148 时间: 2020-2-2 21:39

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

谢谢！

欢迎光临计算化学公社 (http://bbs.keinsci.com/)