计算化学公社

标题: 求助，稀土金属使用Stuttgart赝势基组SCF不收敛，有什么解决方法吗 [打印本页]

作者
Author: naive 时间: 2023-7-6 16:58
标题: 求助，稀土金属使用Stuttgart赝势基组SCF不收敛，有什么解决方法吗
(, 下载次数 Times of downloads: 6)

作者
Author: 筱婷397 时间: 2023-7-6 17:02
可以看sob老师的博文
1.SCF不收敛解决办法：http://sobereva.com/61
2.Gaussian对镧系金属配合物的量子化学计算：http://sobereva.com/581

作者
Author: zjxitcc 时间: 2023-7-6 17:29
本帖最后由 zjxitcc 于 2023-7-6 17:30 编辑

（1）你这体系是闭壳层单重态么？怎么用的是RDFT？
（2）B3LYP并不适合用于过渡金属、镧系锕系金属配合物。建议阅读《简谈量子化学计算中DFT泛函的选择》http://sobereva.com/272
（3）必要时应恰当使用片段组合波函数作为SCF初始猜测。
另外，建议晒出除坐标外的所有信息，便于检查关键词、基组和赝势是否合理。

作者
Author: naive 时间: 2023-7-6 17:42

筱婷397 发表于 2023-7-6 17:02
可以看sob老师的博文
1.SCF不收敛解决办法：http://sobereva.com/61
2.Gaussian对镧系金属配合物的量子化 ...

好的，谢谢，我先去了解一下

作者
Author: naive 时间: 2023-7-6 17:51

zjxitcc 发表于 2023-7-6 17:29
（1）你这体系是闭壳层单重态么？怎么用的是RDFT？
（2）B3LYP并不适合用于过渡金属、镧系锕系金属配合物 ...

感谢问答，我是新手，刚使用这个软件，第一个问题不太了解，我是使用师兄做过的来类比使用的。
第二个问题，配合物是使用的B3LYP，不知道会不会有影响，另外再请麻烦能帮我看一下基组合理吗，下面相关是设置
#p opt freq rb3lyp/genecp geom=connectivity

Title Card Required
...

H O Si 0
6-31+G*
****
Pr 0
S 3 1.00
5.7618400 -0.3811710
...
D 1 1.00
0.0461390 1.0
****

Pr 0
ECP48MWB 4 48
G-Komponente
1
2 1.000000 0.000000
S-G

作者
Author: zjxitcc 时间: 2023-7-6 23:22
本帖最后由 zjxitcc 于 2023-7-6 23:24 编辑

naive 发表于 2023-7-6 17:51
感谢问答，我是新手，刚使用这个软件，第一个问题不太了解，我是使用师兄做过的来类比使用的。
第二个问 ...

体系的电荷和自旋多重度与实验真实情况、研究问题有关，与计算程序无关，这两个数如果写错了，计算就完全没意义了，所以一定要结合实验背景、化学常识等确认一下电荷和自旋多重度。

如果体系是中性的，又没有明显带负电荷的物种/基团，可以不用弥散函数（这会给计算带来很多好处）；如果有明显带负电荷的物种/基团，可以只给那几个重要原子加弥散函数，没必要给所有非金属原子都用加弥散函数。

作者
Author: sobereva 时间: 2023-7-7 10:58

naive 发表于 2023-7-6 17:51
感谢问答，我是新手，刚使用这个软件，第一个问题不太了解，我是使用师兄做过的来类比使用的。
第二个问 ...

不适合B3LYP，仔细看前面的人回帖里我的博文

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