24小时热门版块排行榜    

CyRhmU.jpeg
查看: 3339  |  回复: 1

maoxiao

木虫 (正式写手)

[求助] 求助大神关于SCF不收敛(Convergence failure -- run terminated)的问题已有1人参与

在Gaussion 计算出现下面问题,求助大神怎么解决?
Error on total polarization charges =  0.01305
SCF Done:  E(RB3LYP) =  -1557.98733960     A.U. after  129 cycles
            NFock=128  Conv=0.17D-04     -V/T= 2.0514
Convergence failure -- run terminated.
Error termination via Lnk1e in D:\Gaussion09\G09W\l502.exe at Wed Sep 11 07:12:33 2024.
Job cpu time:       0 days  7 hours 45 minutes  7.0 seconds.
File lengths (MBytes):  RWF=    499 Int=      0 D2E=      0 Chk=     28 Scr=      1
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

paramecium86

版主 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
感谢参与,应助指数 +1
maoxiao: 金币+100, ★★★★★最佳答案 2024-09-18 14:37:06
可以参考我这个帖子里的教程 2.2节   https://muchong.com/t-14242735-1


######################################################################################
A. 首先检查输入结构是不是合理。绝大部分情况下都是因为结构跟 “真实” 结构差的很远,导致的SCF难以收敛。尽量给一个能给出的最好结构 (比如先用UFF粗略优化一下,再用PM6 在优化一下。然后再用DFT计算之类的)
B. 检查体系的自旋多重度是否设置的合理。尤其是过渡金属体系。如果多重度设置的偏离 “真实” 情况很远,则非常难以收敛
C. 检查基组是否合理。尤其是包含重原子,是否合理的设置了混合基组之类的都会影响。具体参考本文1.1节设置混合基组和赝势的方法。
D. 检查当前给的初始结构是否真实存在于 “气相” 中。因为有时候气相和溶液相结构很不一样。这就会导致初始结构以为很准其实是很差的情况。如果是溶液中的结构则需要加上溶剂模型SCRF,见本文1.2 .有时候加了溶剂模型就可以帮助收敛。
E. 以上确定是没有问题的,这时候就考虑加入一些关键词帮助收敛
SCF=(XQC,MAXCONVENTIONALCYCLES=64)
这个关键词可以应用在绝大部分情形下。但是绝不要一上来就用,一定要检查之前的步骤一都没问题了再去使用。G09的话就用SCF=XQC就好。 G16可以写这一串。
SCF=NoVarAcc
高斯在计算SCF过程中,最开始的几步为了加速会降低积分精度。这样有时候就会让波函数不小心进入了一个错误的区间。很难再优化下去。所以可以用这个关键词,取消这种加速降低积分精度的方法。
SCF(VSHIFT=500) 这个一般可以用在过渡金属体系。这种体系,可能HOMO-LUMO间距过小导致波函数不容易收敛。可以计算中人为太高这个值。这个不会影响结果只会影响收敛过程。一般设定在500上下比较合理。

INT(Acc2E=12, UltraFine)
如果用G09进行计算,也可以考虑这个关键词提高了DFT的积分精度和格点精度。这个仅仅适用在某些DFT计算中。比如M062X 这类对于格点要求比较高的泛函中。G16中这个已经是默认值不需要写。

实际计算可以用以上提到的一个或者几个关键词一起帮助收敛。
如果体系实在是太过巨大,这些都做了还是很难收敛。可以使用

SCF(CONVER=6)
这个关键词降低SCF收敛限到10-6, 而默认是10-8。但是慎用。这样的结构可能得到的根本不准。最好是优化完成后,读取波函数和结构恢复到默认的收敛限再进行一次优化。
2楼2024-09-12 10:28:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 maoxiao 的主题更新
信息提示
请填处理意见