Gaussian结构优化某一步后突然发散

fresco604 · 发表于 Post on 2022-3-30 10:28:07

本帖最后由 fresco604 于 2022-3-30 10:31 编辑

求助各位老师，我有一个150原子左右的配合物体系结构优化
使用 Gaussian 16

#p opt b3lyp/6-31g(d,p) em=gd3bj

在37步左右突然开始发散，请问这种情况应该如何处理？能不能提取发散之前的结构（第37步的结构），加入什么关键词，就可以抑制后面的发散？谢谢！

妙角不脆 · 发表于 Post on 2022-3-30 14:09:01

这个还没出震荡，你可以让他继续跑下去，说不定跑到后面就收敛了。或者你要觉得这一桢结构就挺接近的，可以提取出来再跑。

fresco604 · 发表于 Post on 2022-3-30 14:29:51

妙角不脆发表于 2022-3-30 14:09
这个还没出震荡，你可以让他继续跑下去，说不定跑到后面就收敛了。或者你要觉得这一桢结构就挺接近的，可以 ...

想问下如果振荡以后的结构不收敛的话，是不是应该提取振荡以前的结构（比如第37步左右）继续跑，是不是需要加入一些关键词避免再次振荡不收敛？比如int=ultrafine, scf=tight之类的？

RAL · 发表于 Post on 2022-3-30 15:07:37

fresco604 发表于 2022-3-30 14:29
想问下如果振荡以后的结构不收敛的话，是不是应该提取振荡以前的结构（比如第37步左右）继续跑，是不是需 ...

能量波动前和波动后的结构很可能是不同的，有很多时候是从势能面平缓但不是极小值点的部分走向真正的极小值点

sobereva · 发表于 Post on 2022-3-31 15:08:25

fresco604 发表于 2022-3-30 14:29
想问下如果振荡以后的结构不收敛的话，是不是应该提取振荡以前的结构（比如第37步左右）继续跑，是不是需 ...

int=ultrafine对G16是默认的，scf=tight从G09开始就是默认的，都是多余关键词
这里还专门说了
常见的多余的和被滥用的Gaussian关键词
http://sobereva.com/331（http://bbs.keinsci.com/thread-3460-1-1.html）

sobereva · 发表于 Post on 2022-3-31 15:09:48

看结构怎么变化的，如果是结构突然改变极大，以突跃的前一步当初始结构继续优化，加上calcfc以及maxstep=10,notrust
参考
量子化学计算中帮助几何优化收敛的常用方法
http://sobereva.com/164

对于过渡金属配合物不建议用B3LYP，此文明确说了
简谈量子化学计算中DFT泛函的选择
http://sobereva.com/272（http://bbs.keinsci.com/thread-536-1-1.html）

fresco604 · 发表于 Post on 2022-3-31 20:23:41

sobereva 发表于 2022-3-31 15:09
看结构怎么变化的，如果是结构突然改变极大，以突跃的前一步当初始结构继续优化，加上calcfc以及maxstep=10 ...

谢谢老师指点！

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[新手求助] Gaussian结构优化某一步后突然发散

评分 Rate

浏览过的版块