镧系元素自洽场收敛报错

zhangby

各位老师好，最近在做镧系元素配合物的几何结构优化，发现用PBE0/6-31G(d)优化时，可以收敛，但出现了一个虚频，然后按照虚频振动方向调整初猜结构后，又出现了自洽场不收敛问题，而且把关键词调整为“opt freq genecp em=gd3 int=fine pbe1pbe scf=maxcyc=200 SCF=conver=6”后，报错“No NMR shielding tensors so no spin-rotation constants.”。下面是输入文件和输出文件，请问各位老师有什么方法可以解决这个报错和不收敛的问题呢？

zjxitcc

您搞错了，No NMR shielding tensors so no spin-rotation constants是正常输出内容之一，不是报错；报错信息是SCF没有收敛到正常的收敛限，out文件末尾写得很清楚

1. SCF Error SCF Error SCF Error SCF Error SCF Error SCF Error SCF Error SCF Error
   
2. 
   
3.                                     ERROR!!!!
   
4. 
   
5.     SCF has not converged.  Gradients and post-SCF results would be GARBAGE!!
   
6. 
   
7. SCF Error SCF Error SCF Error SCF Error SCF Error SCF Error SCF Error SCF Error
   

复制代码

这个报错对应的是你写了scf(conver=6)，它只适合初步计算，后续仍应当收敛到正常的收敛限。

您在附件Eu.gjf里没给出基组和赝势信息，估计是在Windows下编辑了一下就直接点鼠标提交Gaussian计算，然后没有及时保存或备份，就再次修改gjf文件，过了十几分钟就不知道自己实际上用的哪个文件做计算、哪个gjf对应哪个out，这是非常不好的习惯。从out文件看，基组和赝势定义有误，例如N原子没有定义相应的基组，首要之急是检查基组和赝势。

在正确定义基组和赝势后，如果还有SCF不收敛问题，可以阅读http://bbs.keinsci.com/thread-45423-1-1.html

zhangby

zjxitcc 发表于 2025-11-6 10:20
您搞错了，No NMR shielding tensors so no spin-rotation constants是正常输出内容之一，不是报错；报错信 ...

好的，非常感谢，我按照您的建议试一下