IOp(8/117=-99)无法跳过L801的SCFChk检查

Kmetsch

在优化一个初始结构偏离较大且开壳层的体系时，我尝试使用社长的解决SCF不收敛问题的方法 - 思想家公社的门口：量子化学·分子模拟·二次元中第六点的放宽SCF收敛限的方式来进行预优化，并按照令Gaussian 16中SCF未收敛到默认收敛限也能继续做后续计算的方法 - 思想家公社的门口：量子化学·分子模拟·二次元的方法，添加关键词IOp(7/127=-99,8/117=-99)来防止G16自动检查SCF收敛限的报错。然而G16任务依然在L801处报错，out文件中显示信息如下：
Leave Link  502 at Thu Sep  5 09:47:16 2024, MaxMem=  2147483648 cpu:             912.2 elap:              29.2
(Enter /usr/local/theochem/g16c01-i2018amd/g16/l801.exe)
DoSCS=F DFT=T ScalE2(SS,OS)=  1.000000  1.000000
SCFChk:  SCF convergence  4.14D-07 required  1.00D-08




SCF Error SCF Error SCF Error SCF Error SCF Error SCF Error SCF Error SCF Error


                                    ERROR!!!!


    SCF has not converged.  Gradients and post-SCF results would be GARBAGE!!


SCF Error SCF Error SCF Error SCF Error SCF Error SCF Error SCF Error SCF Error




Error termination via Lnk1e in /usr/local/theochem/g16c01-i2018amd/g16/l801.exe at Thu Sep  5 09:47:16 2024.

显然在计算过程中，G16进行了SCFChk的收敛限检查，并转为tight的收敛判断从而终止了任务。但我查找这个任务的前一段输出，发现在L701中并未启用SCFChk，任务正常推进了

(Enter /usr/local/theochem/g16c01-i2018amd/g16/l701.exe)
... and contract with generalized density number  0.
Compute integral first derivatives.

zjxitcc

完全不需要这几个绕过SCF检查的IOP，可以以更为强大的关键词和收敛算法强硬收敛，除非是本身结构不合理，例如MoS2材料非要截成团簇来计算。

wal

看看是不是过完link1丢了

Kmetsch

zjxitcc 发表于 2026-8-8 15:08
完全不需要这几个绕过SCF检查的IOP，可以以更为强大的关键词和收敛算法强硬收敛，除非是本身结构不合理，例 ...

请问这里的

更为强大的关键词和收敛算法

是指哪个呢，因为我的体系是含过渡金属的开壳层反应中间体，SCF收敛非常难搞TAT求教

Kmetsch

wal 发表于 2026-8-8 15:12
看看是不是过完link1丢了

是的，8/117=-99在Link 1 中没有被载入。请问这会是什么原因导致的呢？
#p opt freq b3lyp/genecp scrf=(solvent=water) em=gd3bj int=ultrafine S
CF(conver=6,novaracc,vshift=500) IOp(7/127=-99,8/117=-99)
----------------------------------------------------------------------
1/18=20,19=15,26=6,38=1/1,3;
2/9=110,12=2,17=6,18=5,40=1/2;
3/5=7,11=2,16=1,17=8,25=1,30=1,70=2201,71=1,72=1,74=-5,75=-5,124=41/1,2,3;
4//1;
5/5=2,6=6,10=500,38=1,53=1/2;
6/7=2,8=2,9=2,10=2,28=1/1;
7/127=-99/1,2,3,16;
1/18=20,19=15,26=6/3(2);
2/9=110/2;
99//99;
2/9=110/2;
3/5=7,6=1,11=2,16=1,17=8,25=1,30=1,70=2205,71=1,72=1,74=-5,75=-5,82=7,124=41/1,2,3;
4/5=5,16=3,69=1/1;
5/5=2,6=6,10=500,38=1,53=1/2;
7/127=-99/1,2,3,16;
1/18=20,19=15,26=6/3(-5);
2/9=110/2;
6/7=2,8=2,9=2,10=2,19=2,28=1/1;
99/9=1/99;
Leave Link    1 at Fri Aug  7 15:09:55 2026, MaxMem=  2147483648 cpu:               3.9 elap:               0.2

zjxitcc

Kmetsch 发表于 2026-8-8 15:56
请问这里的
是指哪个呢，因为我的体系是含过渡金属的开壳层反应中间体，SCF收敛非常难搞TAT求教

只需仔细阅读这两个链接中的计算技巧，完全不需要那些毒害用户的、绕过SCF检查的IOP。
(1) http://bbs.keinsci.com/forum.php ... 20&fromuid=2632
(2) 谈谈片段组合波函数与自旋极化单重态
http://sobereva.com/82