在ORCA里面有必要增加SCF迭代次数吗

魔卡少女晗

这里已经循环到125cycle，然后ORCA提示波函数不收敛。
116  -2540.4066240555   0.000008569255 0.00016353  0.00000066  0.0106526 0.9200
117  -2540.4066194433   0.000004612129 0.00015358  0.00000062  0.0095098 0.9200
118  -2540.4066165093   0.000002934076 0.00014220  0.00000058  0.0095216 0.9200
119  -2540.4066161669   0.000000342337 0.00012845  0.00000052  0.0036936 0.9200
120  -2540.4066197852  -0.000003618256 0.00011680  0.00000048  0.0036915 0.9200
121  -2540.4066258284  -0.000006043168 0.00133950  0.00000546  0.0029269 0.0000
122  -2540.4066514811  -0.000025652726 0.00000899  0.00000003  0.0476809 0.9200
123  -2540.4066525558  -0.000001074733 0.00000890  0.00000002  0.0434154 0.9200
124  -2540.4066535046  -0.000000948799 0.00001000  0.00000002  0.0386176 0.9200

               *****************************************************
               *                      ERROR                        *
               *        SCF NOT CONVERGED AFTER 125 CYCLES         *
               *****************************************************

但是从Max-DP     RMS-DP的数值来看，已经离收敛标准很近了，各位老师能不能帮我看下，如果增加循环次数，是不是有希望达到收敛呢？

hdhxx123

你用的是DIIS的话倒数第二列看力[F,P]离收敛限还远着，收敛方法可以参考帖子http://bbs.keinsci.com/thread-20748-1-1.html和李老师的帖子http://bbs.keinsci.com/thread-4570-1-1.html以及sob老师博文里说的把高斯收敛的波函数传到ORCA帮助收敛

魔卡少女晗

hdhxx123 发表于 2021-5-13 11:52
你用的是DIIS的话倒数第二列看力[F,P]离收敛限还远着，收敛方法可以参考帖子http://bbs.keinsci.com/thread ...

好的好的，非常感谢老师！！我试试把高斯的波函数导进去。另外还有一个问题想请教您，计算BSSE的时候，复合物AB的SCF收敛了，但是片段A的波函数不收敛，能否将片段A的坐标复制出来，用高斯得到波函数，然后也加进orca输入文件中JOB2的命令行呢？

hdhxx123

魔卡少女晗 发表于 2021-5-13 14:35
好的好的，非常感谢老师！！我试试把高斯的波函数导进去。另外还有一个问题想请教您，计算BSSE的时候，复 ...

理论上说，泛函不变，基组zeta数不变时， 高斯得到稳定的波函数传到orca里收敛也可以视作稳定波函数。不清楚BSSE计算中对波函数有没有额外处理，如果只是单点能相减的话原则上是没有问题的，当然如果你都单独放到高斯里收敛完再传到orca的行为肯定不会被质疑，最多耗时长了一点

魔卡少女晗

hdhxx123 发表于 2021-5-13 14:42
理论上说，泛函不变，基组zeta数不变时， 高斯得到稳定的波函数传到orca里收敛也可以视作稳定波函数。不 ...

嗯嗯我明白啦，谢谢老师！

LJF

从能量上来看，应该是小于两百个原子的体系？RMS-DP最容易达到收敛标准，但Max-DP的数值其实还是不小，而且从这几个数值来看，应该还在震荡。