计算化学公社

标题: 镧系和过渡金属的hybrids SCF不收敛 [打印本页]

作者
Author: lanthanides 时间: 2020-7-19 21:18
标题: 镧系和过渡金属的hybrids SCF不收敛
求助～
我有个Ru和Ln的hybrids，已经拿到了单晶数据，我在此基础上用PBE0，B3lyp配合 Ru和Ln SDD，其他的就用6-311G(d). 不过都是前几轮SCF后就不收敛了。简化了结构后也不行，加了Vshift=400 也不行。请问有镧系金属大神嘛，能否给些建议。

我单独算过Ln系部分和Ru部分，都能拿到优化结果。只不过当这两个合在一起后问题就变得比较棘手。

感激～～

作者
Author: sobereva 时间: 2020-7-19 21:40
尝试其它做法
解决SCF不收敛问题的方法
http://sobereva.com/61

并且尝试大核赝势

作者
Author: lanthanides 时间: 2020-7-24 10:30

sobereva 发表于 2020-7-19 21:40
尝试其它做法
解决SCF不收敛问题的方法
http://sobereva.com/61

卢老师，我这几天尝试了大核赝势，时间缩短了很多，而且能够收敛。不过分子结构有些变化，比如两个苯环之间的二面角没了，之前是28°左右。

我的目标主要是研究HOMO-LUMO，非线性值。如果我用大核赝势的话，结果可以采用么，还是我应该在大核基础上改成小核继续算呢？

我的产物是一系列的镧系，如果我都用大核赝势的话，结果的趋势有可比性么？

作者
Author: sobereva 时间: 2020-7-25 00:06

lanthanides 发表于 2020-7-24 10:30
卢老师，我这几天尝试了大核赝势，时间缩短了很多，而且能够收敛。不过分子结构有些变化，比如两个苯环之 ...

用大核赝势优化后收敛的结构，改成小核赝势继续做，看能否成功。结构经过初步优化后，SCF通常会更容易收敛。

大核赝势可以用，但注意有些大核赝势针对不同氧化态分别拟合了赝势，得恰当选取。

作者
Author: lanthanides 时间: 2020-7-25 15:59

sobereva 发表于 2020-7-25 00:06
用大核赝势优化后收敛的结构，改成小核赝势继续做，看能否成功。结构经过初步优化后，SCF通常会更容易收 ...

嗯嗯，好的卢老师。

我还有个关于大核赝势拟合赝势的问题，比如Eu的SDD大核赝势有ECP52MWB， ECP53MWB。。52，53代表的是内核电子数，我现在做的是Eu3+用的是ECP52MWB，这个ECP53MWB代表的是Eu2+么？

我在网上搜了搜，没找到相关的帖子，可能是太基础了

作者
Author: lanthanides 时间: 2020-7-26 00:21
我好像解决这个问题了，我把我找到贴在这防止其他人也有这个问题哈

参考这篇文献M. Dolg, H. Stoll, A. Savin, H. Preuss, Theor. Chim. Acta 75, 173 (1989)
’trivalent rare earth: Q = 11 and m = 0 for La through m = 14 for Lu; divalent rare earth: Q = 10 and m = 1 for La through m = 14 for Yb‘

然后根据表格ECP52MWB对应的是Q=11，Eu3+; ECP53MWB对应的是Q=10，Eu2+

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