过渡态优化SCF不收敛

xujian · 发表于 Post on 2017-3-10 11:05:44

本帖最后由 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优化好的反应物重新画过过渡态，但都不行，请问各位老师，有什么好的建议？

liyuanhe211 · 发表于 Post on 2017-3-10 11:34:30

不保密的话上传输出文件看看

sobereva · 发表于 Post on 2017-3-10 11:49:39

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

xujian · 发表于 Post on 2017-3-10 14:19:42

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

老师，附件中包含了B3LYP下正常结束的输出文件，以及CBS-QB3下用guess=read，以及上述关键词的输出文件，麻烦老师查看，谢谢

PhSOconformer3to5.rar (167.56 KB, 下载次数 Times of downloads: 69)

xujian · 发表于 Post on 2017-3-10 14:20:54

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

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

liyuanhe211 · 发表于 Post on 2017-3-10 15:07:17

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

复制代码

liyuanhe211 · 发表于 Post on 2017-3-10 15:23:01

本帖最后由 liyuanhe211 于 2017-3-10 15:44 编辑

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

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

xujian · 发表于 Post on 2017-3-10 16:30:13

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

老师，我用#p计算过了，用scf=maxcyc=500,感觉SCF后能量第七位一直在变化，看不出是否具有收敛趋势

PhSOconformer3to5.log (865.71 KB, 下载次数 Times of downloads: 6)

xujian · 发表于 Post on 2017-3-10 16:31:04

本帖最后由 xujian 于 2017-3-10 16:34 编辑

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

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

哈哈哈 · 发表于 Post on 2018-5-11 15:48:34

本帖最后由哈哈哈于 2018-5-18 17:41 编辑

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

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

hmlyz · 发表于 Post on 2020-10-26 09:20:49

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

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

snljty · 发表于 Post on 2020-10-26 10:47:50

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

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

moiraaa · 发表于 Post on 2025-11-22 13:18:30

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

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

wal · 发表于 Post on 2025-11-22 13:28:54

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

其实你可以现场试一下的

moiraaa · 发表于 Post on 2025-11-22 14:04:29

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

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

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[新手求助] 过渡态优化SCF不收敛

评分 Rate

评分 Rate

浏览过的版块