优化铀酰卤化物阴离子配合物结构报错

AndrewX

在计算不同构型铀酰卤化物阴离子配合物的结构时，加int=ultrafine或SCF=QC才能正常结束，否则就报错L401( Inaccurate quadrature in CalDSu)或L502(Convergence failure -- run terminated)，在阅读Sob关于解决SCF不收敛的博文后了解到必须加QC才能收敛时需要慎重考虑该结构是否稳定，但该如何进一步处理此种情况呢？另外，针对Inaccurate quadrature in CalDSu报错，int=ultrafine的确有奇效，但正如Sob提到的，该关键词会改变能量，我认为这几种结构应当能够不加多余关键词正常优化完，麻烦老师们帮忙看看该如何处理。PS: 尝试过将加了qc或ultrafine后正常结束的输出文件另存为新的输入不加qc和ultrafine接着算，但就会出现上述报错，相关输入文件和输出文件见附件。

wzkchem5

scf=qc做完以后应当用stable关键字检查稳定性。不要为了绕过stable就避开scf=qc，因为那些只有scf=qc才能收敛的体系，即使不用scf=qc也有很大概率得到不稳定波函数，所以stable永远是需要做的，除非你算过类似的体系发现波函数是稳定的，而你当前体系的波函数和那个体系定性相符，此时可以不做stable。sob老师说的是，“必须加QC才能收敛时需要慎重考虑该结构是否稳定”，不是说“QC收敛的波函数需要慎重考虑该结构是否稳定”，也就是说你用qc收敛了、又不加qc接着qc收敛的波函数重新收敛一次，结果仍然“需要慎重考虑该结构是否稳定”，不加qc接着算并不能”洗白“qc的不稳定倾向

红米饭1234

先用xtb跑一下在放在高斯里面优化会快一点

AndrewX

wzkchem5 发表于 2021-10-11 19:33
scf=qc做完以后应当用stable关键字检查稳定性。不要为了绕过stable就避开scf=qc，因为那些只有scf=qc才能收 ...

我之前都是加了qc收敛了之后再去掉qc另存新的输入文件重新算，感谢wzkchem5老师的纠正！但如何理解“当前体系的波函数和那个体系定性相符”呢？例如我在附件中的这几个文件，只替换了其中某个原子，或者换掉一个基团，整体的电荷和多重度一样，这样算是定性相符的体系吗？

wzkchem5

AndrewX 发表于 2021-10-11 15:29
我之前都是加了qc收敛了之后再去掉qc另存新的输入文件重新算，感谢wzkchem5老师的纠正！但如何理解“当前 ...

如果替换的基团和原来的基团很像，比如用甲基替换配体上的氢，这种应该问题不大。但是如果比如说把金属上配的Cl换成CN，那就不好说了，因为配体场的强度变了很多。具体不好一概而论，如果不确定就统一都做stable看看，如果scf很难收敛的话，做stable的时间应该比scf本身要快。
所谓波函数定性相符，至少需要看电荷布居、自旋布居等等，保险起见再看个自旋密度分布图。光是整体的电荷和多重度一样，什么也不能说明。

AndrewX

wzkchem5 发表于 2021-10-11 23:27
如果替换的基团和原来的基团很像，比如用甲基替换配体上的氢，这种应该问题不大。但是如果比如说把金属上 ...

明白了，感谢老师解答！