求助稀土金属配合物自旋多重度的判断

zzzzhang199311

计算对象：中性，总电子数309个，稀土金属Gd配合物。
使用Gaussian对该配合物进行计算，关键词：opt freq b3pw91/genecp nosymm；尝试了2,4,6,8,10,12等自旋多重度。其中只有自旋多重度为8时顺利完成计算。使用其它自旋多重度计算，第一步时就报错：

Annihilation of the first spin contaminant:
S**2 before annihilation     3.6631,   after     8.6136
Convergence failure -- run terminated.
Error termination via Lnk1e in /public/home/gws_dmcui_02/g09/l502.exe


问题：在论坛里看到sob老师说，基态比激发态更容易收敛，我能根据这个来判断，该稀土金属配合物基态时的自旋多重度就是8吗？（想偷个懒）


   还请大佬们指点迷津！

zjxitcc

不能判断。至少应该算2~3种自旋多重度。先用化学直觉判断最可能的那个，然后再算2种（比最可能的高一点的，低一点的）。SCF不收敛就去解决SCF不收敛。波函数稳定性也没判断，这离下结论还有好几天的工作量。。

bomsaude

对于Gd3+中性配合物，就是Gd的6s2 4f8失去3个电子，变成6s0 4f7，结构化学中的半充满稳定结构。 你现在自旋多重度尝试2,4,6,8等，说明你采用了小核如SDD之类的基组把7个4f电子考虑进去了，这时候7个f电子自旋一样，自旋多重度考虑首先考虑8，是常识，但自旋多重度考虑10,12，应该没有什么意义吧。另外，选2,4,6出错，感觉未必是自旋多重度的问题。
如果用大核计算的话，那7个f电子就不考虑了，自旋多重度就是1了。

sobereva

多看看下文的讨论
使用Gaussian做镧系金属配合物的量子化学计算
http://sobereva.com/581（http://bbs.keinsci.com/thread-21274-1-1.html）

基态电子态往往更容易收敛，决不代表容易收敛的就一定是基态电子态

zzzzhang199311

zjxitcc 发表于 2021-4-16 16:33
不能判断。至少应该算2~3种自旋多重度。先用化学直觉判断最可能的那个，然后再算2种（比最可能的高一点的， ...

谢谢老师，还以为能偷个懒呢，现在老老实实全算了吧

zzzzhang199311

bomsaude 发表于 2021-4-19 00:30
对于Gd3+中性配合物，就是Gd的6s2 4f8失去3个电子，变成6s0 4f7，结构化学中的半充满稳定结构。 你现在自旋 ...

感谢回复，已经用8的自旋多重度做了后续计算，发现能量已经不对劲了。确实是选错了，哎.......

zzzzhang199311

sobereva 发表于 2021-4-21 17:41
多看看下文的讨论
使用Gaussian做镧系金属配合物的量子化学计算
http://sobereva.com/581（http://bbs.ke ...

前两天因为错误发帖被关小黑屋了，现在才来赶来回复，谢谢sob老师！