| 查看: 5472 | 回复: 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 |
回复此楼» 猜你喜欢
Tc数据拟合
已经有1人回复
-
寻求物理化学学报最新的投稿模板！
已经有4人回复
-
物理化学论文润色/翻译怎么收费?
已经有243人回复
-
韩国釜山大学储能电池与功能材料实验室诚招博后（电池/电化学/器件制备/MOF合成）
已经有187人回复
-
祈福！
已经有79人回复
-
投稿咨询
已经有6人回复
-
自科面上 散金币 祈福
已经有61人回复
-
散金求好运
已经有103人回复
-
博士生招生, 材料模拟方向，University of South Carolina
已经有0人回复
-
Postdoc position, materials simulations/machine learning, U. of South Carolina
已经有0人回复
-
祈福青基
已经有34人回复
paramecium86
版主 (著名写手)
- QC强帖: 5
- 应助: 588 (博士)
- 贵宾: 0.01
- 金币: 26695.9
- 红花: 349
- 帖子: 2029
- 在线: 1461.6小时
- 虫号: 2466443
- 注册: 2013-05-16
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 量子化学
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
感谢参与，应助指数 +1
maoxiao: 金币+100, ★★★★★最佳答案 2024-09-18 14:37:06
感谢参与，应助指数 +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
40
