求助：几何优化时加入溶剂模型报错

beidou12138

老师们好！我在使用m062x-def2svp对一个分子进行几何优化，输入指令如下：#p opt freq m062x/def2svp scrf(smd,solvent=benzene)
现在发现几何优化会震荡，于是我尝试加入了ultrafine，随后会出现如下报错：
>>>>>>>>>> Convergence criterion not met. Error on total polarization charges = ********
SCF Done:  E(RM062X) =  -10997.6975194     A.U. after  129 cycles            
NFock=128  Conv=0.91D-01     -V/T= 2.3970 KE= 7.872480186971D+03 PE=-5.349930099180D+04 EE= 1.816038597129D+04
SMD-CDS (non-electrostatic) energy       (kcal/mol) =     -25.85 (included in total energy above)
Convergence failure -- run terminated. Error termination via Lnk1e in /G16/g16/l502.exe at Sat Aug  1 05:20:22 2026.  
我后面尝试加入XQC和caclfc均会出现此类报错，中途使用gview查看，发现会出现图片中的情况，最后还是会出现上述的报错。

我的结构是126个原子，之前使用过不加入溶剂模型的时候是可以收敛的，这次使用的初猜结构使用的是没有加溶剂模型的结果。
求助各位老师这个应该如何处理呢？

fatpig

smd用于优化，报错几率会增大很多。
建议换成pcm...

陈AG

http://sobereva.com/164

wal

结合你的描述，看这个哗啦一下降30万Ha的能量曲线，像是溶剂引起的Inv3 failed in PCMMkU的一种隐蔽变体，可以参考这个帖子尝试更改溶剂表面类型试试

这个问题也有可能不报错，但造成 SCF 能量明显异常。这个错误比较隐蔽，当出现能量明显不对的时候应该检查。

例如我在一次优化中发现有两次 SCF 迭代后的能量（用>>>标出）低了好几十个 Hartree，明显是不对的，SCF 迭代所需的次数也明显较多。这也是上面所述的溶剂空腔问题。

        Line 68580:  SCF Done:  E(RB-P86) =  -3825.13599635     A.U. after   13 cycles
        Line 69843:  SCF Done:  E(RB-P86) =  -3825.13662431     A.U. after   12 cycles
        Line 71105:  SCF Done:  E(RB-P86) =  -3825.13636208     A.U. after   12 cycles
>>> Line 73188:  SCF Done:  E(RB-P86) =  -3845.72449495     A.U. after   54 cycles
        Line 74770:  SCF Done:  E(RB-P86) =  -3825.13665330     A.U. after   20 cycles
        Line 76261:  SCF Done:  E(RB-P86) =  -3825.13647203     A.U. after   23 cycles
        Line 77804:  SCF Done:  E(RB-P86) =  -3825.13594059     A.U. after   26 cycles
        Line 79343:  SCF Done:  E(RB-P86) =  -3825.13482223     A.U. after   26 cycles
        Line 81098:  SCF Done:  E(RB-P86) =  -3825.13408631     A.U. after   35 cycles
>>> Line 83672:  SCF Done:  E(RB-P86) =  -3888.91182114     A.U. after   73 cycles
        Line 85289:  SCF Done:  E(RB-P86) =  -3825.13665328     A.U. after   22 cycles
        Line 86819:  SCF Done:  E(RB-P86) =  -3825.13647358     A.U. after   25 cycles
        Line 88367:  SCF Done:  E(RB-P86) =  -3825.13595614     A.U. after   26 cycles
        Line 89942:  SCF Done:  E(RB-P86) =  -3825.13490799     A.U. after   27 cycles

http://bbs.keinsci.com/forum.php ... =6682&fromuid=63020

sobereva

没事甭用XQC，我反复强调过无数次了
认真看下文
解决SCF不收敛问题的方法
http://sobereva.com/61

不建议M06-2X优化，下文说了
谈谈量子化学研究中什么时候用B3LYP泛函优化几何结构是适当的
http://sobereva.com/557（http://bbs.keinsci.com/thread-17899-1-1.html）

若无特殊情况，opt freq目的甭用SMD，始终用scrf默认的IEFPCM，下文明确说了
谈谈隐式溶剂模型下溶解自由能和体系自由能的计算
http://sobereva.com/327（http://bbs.keinsci.com/thread-3345-1-1.html）
另外，苯的极性很弱，多数情况下，用苯作为溶剂和在真空下优化出来的结构没有明显差异，可以用真空下优化的结果作为合理近似