计算化学公社

标题: 寻找过渡态/柔性扫描中SCF不收敛求助 [打印本页]

作者
Author: wrtgcn 时间: 2020-1-25 09:41
标题: 寻找过渡态/柔性扫描中SCF不收敛求助
计算一个气体小分子在催化剂表面吸附后断键过程的过渡态。用QST2法生成的初猜结果不是很好，第一个结构的单点能都收敛失败；而用TS法根据末态手动摆的结构可能也不是太好，无法计算到过渡态。无奈下想用柔性扫描找个初猜。用PM6法计算第一个点就SCF不收敛；提高到HF/3-21G后，跑了10步左右再次出现SCF不收敛（共18步）。之前一般在结构很差时才会碰到SCF不收敛的情况。请问大家，在过渡态/柔性扫描过程中，碰到这样的情况，有什么好的办法吗？谢谢大家。

作者
Author: lijiayisjtu 时间: 2020-1-25 12:27
解决SCF不收敛问题的方法 http://sobereva.com/61 看完基本搞定

作者
Author: Machaophy 时间: 2020-1-25 15:09
HF/3-21G 计算的方法不可信吧，没有用过HF。。。。。。。。。。。。
套用Sob老师的话，最穷的级别B3LYP/6-31g*,这个都用不起就别做量化计算了
推荐M06-2X，这个找过渡态性能好些
你不贴文件出来，不知道你怎么做的呀

作者
Author: sobereva 时间: 2020-1-25 23:16
所有可能的原因和解决办法在http://sobereva.com/61那篇文章里都提了。
没有具体文件没法给出确切的建议。
HF/3-21G算了也白算
PM6都SCF不收敛，有较大可能输入文件就存在显著硬伤。如果催化剂里有过渡金属，改用xtb靠谱得多，Gaussian结合xtb找过渡态又快又容易，见《将Gaussian与Grimme的xtb程序联用搜索过渡态、产生IRC、做振动分析》（http://sobereva.com/421）

作者
Author: wrtgcn 时间: 2020-1-26 06:21

Machaophy 发表于 2020-1-25 15:09
HF/3-21G 计算的方法不可信吧，没有用过HF。。。。。。。。。。。。
套用Sob老师的话，最穷的级别B3LYP/6- ...

想用HF/3-21G扫描个靠谱的初始结构，再提高级别做过渡态计算。之前用B3LYP/6-31g*做柔性扫描，28个小时才跑了10个点不到，每个点的SCF收敛的很慢。

作者
Author: wrtgcn 时间: 2020-1-26 06:39

sobereva 发表于 2020-1-25 23:16
所有可能的原因和解决办法在http://sobereva.com/61那篇文章里都提了。
没有具体文件没法给出确切的建议。 ...

谢谢社长。体系中没有过渡金属，但是气体分子是臭氧，有没有可能是其断键过程中的电子分布比较特殊，导致SCF难以收敛？

作者
Author: sobereva 时间: 2020-1-26 18:54

wrtgcn 发表于 2020-1-26 06:39
谢谢社长。体系中没有过渡金属，但是气体分子是臭氧，有没有可能是其断键过程中的电子分布比较特殊，导致 ...

O3属于很难算的体系，静态相关强，HF的结果是离谱的，要是再牵扯其断键，那事情就更复杂了。
SCF难收敛跟这个不无关系

作者
Author: wrtgcn 时间: 2020-1-27 01:49

sobereva 发表于 2020-1-26 18:54
O3属于很难算的体系，静态相关强，HF的结果是离谱的，要是再牵扯其断键，那事情就更复杂了。
SCF难收敛 ...

感谢您的回复。如果确实想要计算O3在材料表面吸附后断键过程的过渡态，您在方法或相应关键词上有什么建议吗？我在这一步上卡了好长时间了，一直没能找到。

作者
Author: sobereva 时间: 2020-1-27 02:21

wrtgcn 发表于 2020-1-27 01:49
感谢您的回复。如果确实想要计算O3在材料表面吸附后断键过程的过渡态，您在方法或相应关键词上有什么建议 ...

用稍微像样的级别，诸如至少B3LYP/6-31G*。HF完全不能接受
解决SCF不收敛的做法我博文里都写了，没有可补充的

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