镧系和过渡金属的hybrids SCF不收敛

lanthanides · 发表于 Post on 2020-7-19 21:18:59

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

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

感激～～

sobereva · 发表于 Post on 2020-7-19 21:40:29

尝试其它做法
解决SCF不收敛问题的方法
http://sobereva.com/61

并且尝试大核赝势

lanthanides · 发表于 Post on 2020-7-24 10:30:44

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

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

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

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

sobereva · 发表于 Post on 2020-7-25 00:06:24

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

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

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

lanthanides · 发表于 Post on 2020-7-25 15:59:58

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

嗯嗯，好的卢老师。

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

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

lanthanides · 发表于 Post on 2020-7-26 00:21:52

我好像解决这个问题了，我把我找到贴在这防止其他人也有这个问题哈

参考这篇文献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+

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[新手求助] 镧系和过渡金属的hybrids SCF不收敛