计算化学公社

标题: 柔性长链脂肪酸优化不收敛求助 [打印本页]

作者
Author: chanyeolhgh 时间: 2020-10-19 09:01
标题: 柔性长链脂肪酸优化不收敛求助
Sob老师好，在进行几何优化的过程中，总是出现508错误，按照您的解决SCF不收敛问题的方法，尝试过SCF(novaracc,noincfock) int=acc2e=12、SCF=XQC、下面是SCF=XQC的输出文件，下一步具体该用什么关键词比较好呢？恳请老师赐教。感谢老师。输出文件：
(, 下载次数 Times of downloads: 3)
输入文件：
(, 下载次数 Times of downloads: 23)

作者
Author: 喵星大佬 时间: 2020-10-19 10:04
本帖最后由喵星大佬于 2020-10-19 10:06 编辑

我寻思你这个酸，质子还在上面，加个负电荷是要做甚。。。。。
变成开壳层了还带弥散函数本来也难收敛。。。。
如果是要算绝热电子亲和能的话，手动微调一下结构，或者设置片段组合波函数把初猜自旋给羧基部分，可能会好收敛点
或者先把他变回闭壳层做个点单，把收敛的波函数作为初猜扔进去

作者
Author: chanyeolhgh 时间: 2020-10-19 12:15

喵星大佬发表于 2020-10-19 10:04
我寻思你这个酸，质子还在上面，加个负电荷是要做甚。。。。。
变成开壳层了还带弥散函数本来也难收敛。。 ...

是的，老师。想要做绝热电子亲合能老师您所说的设置片段组合波函数把初猜给羧基部分具体要怎么做呢？小的水平不够，还望老师指点，感谢老师。

作者
Author: wzkchem5 时间: 2020-10-19 13:20

chanyeolhgh 发表于 2020-10-19 12:15
是的，老师。想要做绝热电子亲合能老师您所说的设置片段组合波函数把初猜给羧基部分具体要怎么做呢？ ...

这个稍微有点麻烦，个人建议一开始先不要尝试做这个。
先试试读取闭壳层波函数，然后加一个比较大的vshift

作者
Author: chanyeolhgh 时间: 2020-10-19 14:05

wzkchem5 发表于 2020-10-19 13:20
这个稍微有点麻烦，个人建议一开始先不要尝试做这个。
先试试读取闭壳层波函数，然后加一个比较大的vshi ...

好的老师感谢指导

作者
Author: GoldenBaby 时间: 2020-10-19 15:15

chanyeolhgh 发表于 2020-10-19 12:15
是的，老师。想要做绝热电子亲合能老师您所说的设置片段组合波函数把初猜给羧基部分具体要怎么做呢？ ...

片段组合波函数可以参考这篇博文http://sobereva.com/82

作者
Author: sobereva 时间: 2020-10-19 23:54
尽量少用xqc，博文已经说了。往往白浪费很多时间，最后依然不收敛。

你可以考虑泛函改用wB97XD试试，并去掉xqc。
并且opt加上maxstep=10,notrust，让每一步结构变化别那么厉害
另外，鉴于你最后一步在启动qc之前SCF已经收敛得还算比较充分了，可以带上scf=conver=6，等任务收敛之后，再用guess=read结合默认的SCF收敛限进一步做opt freq

作者
Author: chanyeolhgh 时间: 2020-10-20 13:33

sobereva 发表于 2020-10-19 23:54
尽量少用xqc，博文已经说了。往往白浪费很多时间，最后依然不收敛。

你可以考虑泛函改用wB97XD试试，并 ...

好的感谢老师

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