计算化学公社

标题: 过渡态优化SCF不收敛 [打印本页]

作者
Author: xujian 时间: 2017-3-10 11:05
标题: 过渡态优化SCF不收敛
本帖最后由 xujian 于 2017-3-10 11:06 编辑

在B3LYP/Def2-TZVPP水平下得到的过渡态，用CB3-QB3优化，出现SCF不收敛的问题（L502，Convergence failure）,尝试过SCF(novaracc,noincfock) int=acc2e=12，SCF=vshift=500，以及用B3LYP的波函数做初猜，guess=read。也直接用CBS-QB3优化好的反应物重新画过过渡态，但都不行，请问各位老师，有什么好的建议？

作者
Author: liyuanhe211 时间: 2017-3-10 11:34
不保密的话上传输出文件看看

作者
Author: sobereva 时间: 2017-3-10 11:49
你是已经在B3LYP/Def2-TZVPP下找到了过渡态结构，然后又要用CBS-QB3继续优化过渡态并算出热力学数据？如果是这样，输入输出文件压缩后传上来

作者
Author: xujian 时间: 2017-3-10 14:19

sobereva 发表于 2017-3-10 11:49
你是已经在B3LYP/Def2-TZVPP下找到了过渡态结构，然后又要用CBS-QB3继续优化过渡态并算出热力学数据？如果 ...

老师，附件中包含了B3LYP下正常结束的输出文件，以及CBS-QB3下用guess=read，以及上述关键词的输出文件，麻烦老师查看，谢谢
(, 下载次数 Times of downloads: 69)

作者
Author: xujian 时间: 2017-3-10 14:20

liyuanhe211 发表于 2017-3-10 11:34
不保密的话上传输出文件看看

谢谢老师，为了节约公社空间，我只上传在给sob老师的回复中了，麻烦查看，谢谢！

作者
Author: liyuanhe211 时间: 2017-3-10 15:07

xujian 发表于 2017-3-10 14:20
谢谢老师，为了节约公社空间，我只上传在给sob老师的回复中了，麻烦查看，谢谢！

输入文件中应使用#p代替#来运行，这样可以显示SCF过程，否则像现在输出文件中具体的SCF信息根本没有。
应当重新带着#p运行文件，观察SCF收敛过程

在你的输出中（失败前）一共有4次成功的SCF，所需的循环数均甚多，可能该体系/构型就是不太容易SCF收敛，最后一次失败前也已收敛到1E-7的水平，故有可能遇到了少见的增加scf=maxcyc可以解决的问题，应带着#p观察最后一次scf过程是否却有收敛趋势，如果确实有，应当增加scf=maxcyc至某个大于128的数字，如scf=maxcyc=200，重新运行。如果没有收敛趋势另说。

Line 341: SCF Done: E(UB3LYP) = -705.035390591 A.U. after 113 cycles
Line 1002: SCF Done: E(UB3LYP) = -705.035740147 A.U. after 97 cycles
Line 1245: SCF Done: E(UB3LYP) = -705.035851068 A.U. after 56 cycles
Line 1490: SCF Done: E(UB3LYP) = -705.036179972 A.U. after 71 cycles
Line 1722: SCF Done: E(UB3LYP) = -705.036288565 A.U. after 129 cycles

复制代码

作者
Author: liyuanhe211 时间: 2017-3-10 15:23
本帖最后由 liyuanhe211 于 2017-3-10 15:44 编辑

xujian 发表于 2017-3-10 14:20
谢谢老师，为了节约公社空间，我只上传在给sob老师的回复中了，麻烦查看，谢谢！

我自己跑了一下你的输入（→_→其实只是为了收集奇怪的错误），从右上角可以看出你的SCF却有收敛趋势，可以考虑从当前结构开始计算(如果当前结构合理)，增加 SCF=maxcyc=N的数值（命令如SCF=(novaracc,noincfock,maxcyc=200)）尝试，如果还不收敛，当前收敛结果也比较高了，稍微升高SCF收敛阈值也可以。

(, 下载次数 Times of downloads: 78)

作者
Author: xujian 时间: 2017-3-10 16:30

liyuanhe211 发表于 2017-3-10 15:07
输入文件中应使用#p代替#来运行，这样可以显示SCF过程，否则像现在输出文件中具体的SCF信息根本没有。
...

老师，我用#p计算过了，用scf=maxcyc=500,感觉SCF后能量第七位一直在变化，看不出是否具有收敛趋势
(, 下载次数 Times of downloads: 6)

作者
Author: xujian 时间: 2017-3-10 16:31
本帖最后由 xujian 于 2017-3-10 16:34 编辑

liyuanhe211 发表于 2017-3-10 15:23
我自己跑了一下你的输入（→_→其实只是为了收集奇怪的错误），从右上角可以看出你的SCF却有收敛趋势，可 ...

好的，我再试试，刚看到您的这条回复，可以忽略我上一条，也可以帮我看看

作者
Author: 哈哈哈 时间: 2018-5-11 15:48
本帖最后由哈哈哈于 2018-5-18 17:41 编辑

liyuanhe211 发表于 2017-3-10 15:07
输入文件中应使用#p代替#来运行，这样可以显示SCF过程，否则像现在输出文件中具体的SCF信息根本没有。
...

老师您好，我有些没看懂。麻烦您解答一下。就是您说所需的循环数均甚多是怎么来判定的多？ #p 看收敛趋势是怎么看，看四个收敛标准代表的值在下降？麻烦您了。

作者
Author: hmlyz 时间: 2020-10-26 09:20

liyuanhe211 发表于 2017-3-10 15:23
我自己跑了一下你的输入（→_→其实只是为了收集奇怪的错误），从右上角可以看出你的SCF却有收敛趋势，可 ...

老师您好，请问您使用的这个脚本有在论坛里发过吗？我只找到了tsmonitor。

作者
Author: snljty 时间: 2020-10-26 10:47

hmlyz 发表于 2020-10-26 09:20
老师您好，请问您使用的这个脚本有在论坛里发过吗？我只找到了tsmonitor。

如果你只是想看SCF收敛曲线，用Multiwfn载入单点任务的out文件（输入文件必须有#P），进入主功能100的子功能6即可。

作者
Author: moiraaa 时间: 2025-11-22 13:18

snljty 发表于 2020-10-26 10:47
如果你只是想看SCF收敛曲线，用Multiwfn载入单点任务的out文件（输入文件必须有#P），进入主功能100的子 ...

老师，您好，但是Multiwfn不是不能输入out文件吗

作者
Author: wal 时间: 2025-11-22 13:28

moiraaa 发表于 2025-11-22 13:18
老师，您好，但是Multiwfn不是不能输入out文件吗

其实你可以现场试一下的

作者
Author: moiraaa 时间: 2025-11-22 14:04

wal 发表于 2025-11-22 13:28
其实你可以现场试一下的

我试了，但是程序直接退出了，我还去官网看了，好像不支持out文件诶

作者
Author: wal 时间: 2025-11-22 14:12

moiraaa 发表于 2025-11-22 14:04
我试了，但是程序直接退出了，我还去官网看了，好像不支持out文件诶

那你上传上来看看吧并没有不支持out文件
(, 下载次数 Times of downloads: 0)

作者
Author: moiraaa 时间: 2025-11-23 13:29
本帖最后由 moiraaa 于 2025-11-23 13:32 编辑

wal 发表于 2025-11-22 14:12
那你上传上来看看吧并没有不支持out文件

好的，老师，这是我退出时的界面，那个是我的out文件
老师，可不可以帮我再看看502报错，一直是502报错o(╥﹏╥)o

作者
Author: wal 时间: 2025-11-23 15:35

moiraaa 发表于 2025-11-23 13:29
好的，老师，这是我退出时的界面，那个是我的out文件
老师，可不可以帮我再看看502报错，一直是502报错o ...

你可能需要更新multiwfn版本了，我的载入一切正常
(, 下载次数 Times of downloads: 0)

作者
Author: moiraaa 时间: 2025-11-23 17:59

wal 发表于 2025-11-23 15:35
你可能需要更新multiwfn版本了，我的载入一切正常

好的好的，谢谢老师！

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