计算化学公社

 找回密码 Forget password
 注册 Register
Views: 17453|回复 Reply: 18
打印 Print 上一主题 Last thread 下一主题 Next thread

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

[复制链接 Copy URL]

588

帖子

0

威望

3310

eV
积分
3898

Level 5 (御坂)

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

88

帖子

0

威望

233

eV
积分
321

Level 3 能力者

19#
发表于 Post on 2025-11-23 17:59:39 | 只看该作者 Only view this author
wal 发表于 2025-11-23 15:35
你可能需要更新multiwfn版本了,我的载入一切正常

好的好的,谢谢老师!

683

帖子

12

威望

2820

eV
积分
3743

Level 5 (御坂)

鸩羽

18#
发表于 Post on 2025-11-23 15:35:24 | 只看该作者 Only view this author
moiraaa 发表于 2025-11-23 13:29
好的,老师,这是我退出时的界面,那个是我的out文件
老师,可不可以帮我再看看502报错,一直是502报错o ...


你可能需要更新multiwfn版本了,我的载入一切正常
某不知名实验组从苞米地里长出来的计算选手

88

帖子

0

威望

233

eV
积分
321

Level 3 能力者

17#
发表于 Post on 2025-11-23 13:29:47 | 只看该作者 Only view this author
本帖最后由 moiraaa 于 2025-11-23 13:32 编辑
wal 发表于 2025-11-22 14:12
那你上传上来看看吧 并没有不支持out文件

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

202511231325383468..png (117.41 KB, 下载次数 Times of downloads: 0)

202511231325383468..png

ZN shuiningjiaoopt.zip

126.65 KB, 下载次数 Times of downloads: 2

ZN shuiningjiaoopt.gjf

7.37 KB, 下载次数 Times of downloads: 2

683

帖子

12

威望

2820

eV
积分
3743

Level 5 (御坂)

鸩羽

16#
发表于 Post on 2025-11-22 14:12:14 | 只看该作者 Only view this author
moiraaa 发表于 2025-11-22 14:04
我试了,但是程序直接退出了,我还去官网看了,好像不支持out文件诶

那你上传上来看看吧 并没有不支持out文件

某不知名实验组从苞米地里长出来的计算选手

88

帖子

0

威望

233

eV
积分
321

Level 3 能力者

15#
发表于 Post on 2025-11-22 14:04:29 | 只看该作者 Only view this author
wal 发表于 2025-11-22 13:28
其实你可以现场试一下的

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

683

帖子

12

威望

2820

eV
积分
3743

Level 5 (御坂)

鸩羽

14#
发表于 Post on 2025-11-22 13:28:54 | 只看该作者 Only view this author
moiraaa 发表于 2025-11-22 13:18
老师,您好,但是Multiwfn不是不能输入out文件吗

其实你可以现场试一下的
某不知名实验组从苞米地里长出来的计算选手

88

帖子

0

威望

233

eV
积分
321

Level 3 能力者

13#
发表于 Post on 2025-11-22 13:18:30 | 只看该作者 Only view this author
snljty 发表于 2020-10-26 10:47
如果你只是想看SCF收敛曲线,用Multiwfn载入单点任务的out文件(输入文件必须有#P),进入主功能100的子 ...

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

1187

帖子

5

威望

2869

eV
积分
4156

Level 6 (一方通行)

12#
发表于 Post on 2020-10-26 10:47:50 | 只看该作者 Only view this author
hmlyz 发表于 2020-10-26 09:20
老师您好,请问您使用的这个脚本有在论坛里发过吗?我只找到了tsmonitor。

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

31

帖子

0

威望

1123

eV
积分
1154

Level 4 (黑子)

11#
发表于 Post on 2020-10-26 09:20:49 | 只看该作者 Only view this author
liyuanhe211 发表于 2017-3-10 15:23
我自己跑了一下你的输入(→_→其实只是为了收集奇怪的错误),从右上角可以看出你的SCF却有收敛趋势,可 ...

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

43

帖子

0

威望

1168

eV
积分
1211

Level 4 (黑子)

10#
发表于 Post on 2018-5-11 15:48:34 | 只看该作者 Only view this author
本帖最后由 哈哈哈 于 2018-5-18 17:41 编辑
liyuanhe211 发表于 2017-3-10 15:07
输入文件中应使用#p代替#来运行,这样可以显示SCF过程,否则像现在输出文件中具体的SCF信息根本没有。
...

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

588

帖子

0

威望

3310

eV
积分
3898

Level 5 (御坂)

9#
 楼主 Author| 发表于 Post on 2017-3-10 16:31:04 | 只看该作者 Only view this author
本帖最后由 xujian 于 2017-3-10 16:34 编辑
liyuanhe211 发表于 2017-3-10 15:23
我自己跑了一下你的输入(→_→其实只是为了收集奇怪的错误),从右上角可以看出你的SCF却有收敛趋势,可 ...

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

588

帖子

0

威望

3310

eV
积分
3898

Level 5 (御坂)

8#
 楼主 Author| 发表于 Post on 2017-3-10 16:30:13 | 只看该作者 Only view this author
liyuanhe211 发表于 2017-3-10 15:07
输入文件中应使用#p代替#来运行,这样可以显示SCF过程,否则像现在输出文件中具体的SCF信息根本没有。
...

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

3097

帖子

29

威望

1万

eV
积分
17221

Level 6 (一方通行)

7#
发表于 Post on 2017-3-10 15:23:01 | 只看该作者 Only view this author
本帖最后由 liyuanhe211 于 2017-3-10 15:44 编辑
xujian 发表于 2017-3-10 14:20
谢谢老师,为了节约公社空间,我只上传在给sob老师的回复中了,麻烦查看,谢谢!

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



评分 Rate

参与人数
Participants 1
eV +1 收起 理由
Reason
sobereva + 1

查看全部评分 View all ratings

3097

帖子

29

威望

1万

eV
积分
17221

Level 6 (一方通行)

6#
发表于 Post on 2017-3-10 15:07:17 | 只看该作者 Only view this author
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,重新运行。如果没有收敛趋势另说。

  1.         Line 341:  SCF Done:  E(UB3LYP) =  -705.035390591     A.U. after  113 cycles
  2.         Line 1002:  SCF Done:  E(UB3LYP) =  -705.035740147     A.U. after   97 cycles
  3.         Line 1245:  SCF Done:  E(UB3LYP) =  -705.035851068     A.U. after   56 cycles
  4.         Line 1490:  SCF Done:  E(UB3LYP) =  -705.036179972     A.U. after   71 cycles
  5.         Line 1722:  SCF Done:  E(UB3LYP) =  -705.036288565     A.U. after  129 cycles
复制代码

评分 Rate

参与人数
Participants 1
eV +1 收起 理由
Reason
sobereva + 1

查看全部评分 View all ratings

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2026-1-24 06:35 , Processed in 0.183001 second(s), 25 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list