计算化学公社

标题: Gaussian结构优化某一步后突然发散 [打印本页]

作者
Author: fresco604 时间: 2022-3-30 10:28
标题: Gaussian结构优化某一步后突然发散
本帖最后由 fresco604 于 2022-3-30 10:31 编辑

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

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

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

(, 下载次数 Times of downloads: 16)
(, 下载次数 Times of downloads: 10)
(, 下载次数 Times of downloads: 14)
(, 下载次数 Times of downloads: 15)

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

作者
Author: fresco604 时间: 2022-3-30 14:29

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

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

作者
Author: RAL 时间: 2022-3-30 15:07

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

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

作者
Author: sobereva 时间: 2022-3-31 15:08

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）

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

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

作者
Author: fresco604 时间: 2022-3-31 20:23

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

谢谢老师指点！

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