| 查看: 3338 | 回复: 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 |
回复此楼» 猜你喜欢
085404（计算机技术）293分 希望调剂到材料相关专业
已经有0人回复
-
ELISA试验中不容忽视的细节盘点(二)
已经有0人回复
-
物理化学论文润色/翻译怎么收费?
已经有138人回复
-
华南师范大学广州市光谱分析与功能探针重点实验室诚招青年英才和博士后岗位研究人员
已经有11人回复
-
华南师范大学广州市光谱分析与功能探针重点实验室诚招青年英才和博士后岗位研究人员
已经有10人回复
-
细胞培养，这22个细节一定要注意！（一）
已经有0人回复
-
CSC访学博后项目获批，外方学校暂停合作该怎么办？
已经有51人回复
-
求标准卡片Cu0.4In0.4Zn1.2S2 PDF#32-0340
已经有2人回复
-
量子点电致发光器件（QLED）：蓝光量子点（QDs）和纳米氧化锌电子传输材料
已经有0人回复
-
需要一个圆偏振光诱导化学反应，请问使用圆偏振光作为光源，应该需要哪些设备搭建
已经有0人回复
paramecium86
版主 (著名写手)
- QC强帖: 5
- 应助: 586 (博士)
- 贵宾: 0.01
- 金币: 26316.1
- 红花: 348
- 帖子: 2021
- 在线: 1440.5小时
- 虫号: 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