镧系Gd3+的结合能高达-4Hartree，是基组错了还是理论方法错了

lph

我想计算一下石墨烯相的氮化碳（g-C3N4）结合Gd3+的能力。
普通原子采用的6-31G*基组，Gd3+采用的Gaussian自带的SDD基组，优化结构主要采用的HF（但得到结构似乎不太合理，有些配位键长度只有1.3埃，文献中相似的健长为2.3埃。换用m62x也没啥改善），
溶剂化的单点能采用的M052x，这样算的的结合能将近-4.3Hartree。

有的文献说用单纯的DFT来计算Gd3+的能量会非常不准，需要用DFT+U方法，但优化出来的健长太短了，会不会基组也选得有问题？

备注：Gd3+电子结构4f7。计算中自旋多重的M=2和M=8多尝试过，似乎对几何构型影响不大。我也试图采用 NH3+Gd3+ 这个4原子小体系来测试一下EMSL中Gd相关的各个基组，但常常出现“ Annihilation of the first spin contaminant:
S**2 before annihilation    15.7637,   after    15.7500
Convergence failure -- run terminated.”这种错误，比较头大。

kekexili_08

还不是基组输入部分给某些原子的基组多输入了一次？ 过短的键长， 过大的结合能， 看起来有这个可能。

sobereva

8成是输入文件有低级错误，诸如没定义赝势、赝势基组定义错误之类，否则不可能配位键键长优化出来成那样。不要赖基组或者理论方法
绝对不要用HF。DFT结果好赖看泛函，无论如何也轮不到HF。量化计算领域没有+U的事。
SCF不收敛看
解决SCF不收敛问题的方法
http://sobereva.com/61

lph

sobereva 发表于 2018-1-30 22:58
8成是输入文件有低级错误，诸如没定义赝势、赝势基组定义错误之类，否则不可能配位键键长优化出来成那样。 ...

多谢！！找到错误了。果然是只写了SDD基组，没写SDD赝势。
原来写的是：
C  N  H  0
6-31G*
****
Gd    0
SDD
****
我改成：
C  N  H  0
6-31G*
****
Gd    0
SDD
****

Gd   0
SDD
来试试