计算化学公社

标题: 【求助】l9999 不收敛之后我继续这样做有意义吗？ [打印本页]

作者
Author: Adolphin 时间: 2019-4-26 13:19
标题: 【求助】l9999 不收敛之后我继续这样做有意义吗？
# opt=(calcfc,ts,noeigen,z-matrix) freq b3lyp/6-31++g(d,p) 基组
再次优化时我用上一次最后的坐标开始优化，已经四次了，请问我这样优化下去有意义吗？本人萌新请求老师解答。

作者
Author: DoubeeTwT 时间: 2019-4-26 13:39
没意义（或者说没必要）

Problem1:
你的结果似乎已经在TS附近了但是时常会跳开，建议加上maxstep=5限制最大步长（如有必要，可以加上notrust固定最大步长为5而不是由Gaussian自己给出）。

Problem2:
你的方法与基组不匹配，b3lyp是一个差劲的方法你却给他一个6-31++g(d,p)那么高精度的基组。。。。建议改用m062x。

Problem3:
建议检查你的体系是否需要那么高精度的基组（体系大小、是否有必要加弥散函数）。如果没那么高的要求建议先在m062x/6-31g(d)下做优化，做完后再用你的高精度算一个单点能，而不是再如此高精度下优化。

作者
Author: Adolphin 时间: 2019-4-26 17:24

DoubeeTwT 发表于 2019-4-26 13:39
没意义（或者说没必要）

Problem1:

谢谢您的回复，我去换一下计算试试。

作者
Author: sobereva 时间: 2019-4-27 00:53
仔细看下文，你重复拿最后一步结构还用之前的关键词优化根本毫无意义！
量子化学计算中帮助几何优化收敛的常用方法
http://sobereva.com/164

基组和泛函选用问题好好看此文及文中引用的博文
坚持使用B3LYP算有机体系的人的下场
http://bbs.keinsci.com/thread-12773-1-1.html

顺带一提，6-31++(d,p)绝对称不上高精度基组，不需要弥散函数情况时，精度和中等偏低的6-31G**基本没太大差别
谈谈量子化学中基组的选择
http://sobereva.com/336（http://bbs.keinsci.com/thread-3545-1-1.html）

一般情况建议B3LYP-D3(BJ)/6-31G*优化，或者更好点结合6-311G*优化就够了，跟氢关系密切时再给氢加极化。

对于优化过渡态，notrust是默认的，不需要自己写。

作者
Author: Adolphin 时间: 2019-4-27 11:26

sobereva 发表于 2019-4-27 00:53
仔细看下文，你重复拿最后一步结构还用之前的关键词优化根本毫无意义！
量子化学计算中帮助几何优化收敛的 ...

谢谢sob老师，我用m062x已经解决了。

作者
Author: Adolphin 时间: 2019-4-27 11:27

DoubeeTwT 发表于 2019-4-26 13:39
没意义（或者说没必要）

Problem1:

谢谢，已经解决了

作者
Author: DoubeeTwT 时间: 2019-4-27 17:08

sobereva 发表于 2019-4-27 00:53
仔细看下文，你重复拿最后一步结构还用之前的关键词优化根本毫无意义！
量子化学计算中帮助几何优化收敛的 ...

谢谢，sob大大指出不足，还需要多学习~

作者
Author: DoubeeTwT 时间: 2019-4-27 17:09

Adolphin 发表于 2019-4-27 11:27
谢谢，已经解决了

解决了就好，sob老师还指出了我说的有瑕疵的地方，抱歉。。。（雾~

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