计算化学公社

标题: 镧系Gd3+的结合能高达-4Hartree，是基组错了还是理论方法错了 [打印本页]

作者
Author: lph 时间: 2018-1-30 21:57
标题: 镧系Gd3+的结合能高达-4Hartree，是基组错了还是理论方法错了
我想计算一下石墨烯相的氮化碳（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.”这种错误，比较头大。

作者
Author: kekexili_08 时间: 2018-1-30 22:09
还不是基组输入部分给某些原子的基组多输入了一次？过短的键长，过大的结合能，看起来有这个可能。

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

作者
Author: lph 时间: 2018-1-31 13:00

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
来试试

欢迎光临计算化学公社 (http://bbs.keinsci.com/)