结构优化和单点能计算能量差别悬殊

ZPLLL · 发表于 Post on 2023-8-7 21:02:20

在优化EDTA-Fe(II)2-结构的过程中使用的是pbe1pbe泛函和混合基组（genecp(sdd；6-311g（d,p) ），能量为-1223.0980949 hartee；然而使用 b2plypd3/def2TZVPP 泛函和基组对其进行单点能计算时，能量为-2363.2416748 hartee。看过以前的帖子，说不同级别下绝对能量没有可比性《相同DFT不同基组计算的单点能差值较大 - 量子化学 (Quantum Chemistry) - 计算化学公社 (keinsci.com)》。但这个能量悬殊过大，想请教一下大佬这个能量可信吗？；此外，只有在基组为混合基组时（genecp(sdd；6-311g（d,p) ），能量为 -1223.0980949 hartee，在同一个泛函下更换了基组（pbe1pbe/def2tzvp em=gd3bj）之后能量为（-2362.8756447）与b2plypd3/def2TZVPP 的相差不大，能量的变化如此之大是因为赝势基组（sdd）的原因吗？

wzkchem5 · 发表于 Post on 2023-8-7 21:45:03

那叫SDD不叫SSD
参见https://mattermodeling.stackexch ... culations/9962#9962

zjxitcc · 发表于 Post on 2023-8-7 21:45:34

是因为赝势的原因。

至于是否可信，还有其他因素要考虑。比如
（1）你B2PLYPD3能量有没读错？看《谈谈该从Gaussian输出文件中的什么地方读电子能量》http://sobereva.com/488
（2）这个体系是开壳层的嘛？如果是，那波函数稳定性检查了没、自旋布居看了没？

ZPLLL · 发表于 Post on 2023-8-7 21:52:34

wzkchem5 发表于 2023-8-7 21:45
那叫SDD不叫SSD
参见https://mattermodeling.stackexchange.com/questions/9961/effective-core-potential ...

感谢指出书写问题，在计算中使用的是SDD，这里由于编辑的时候打顺手，下次还是要更细心。

ZPLLL · 发表于 Post on 2023-8-7 22:02:54

zjxitcc 发表于 2023-8-7 21:45
是因为赝势的原因。

至于是否可信，还有其他因素要考虑。比如

这个体系是开壳层，感谢解答并指出检查的方向；我后面去查看一下它的波函数和自旋布居

ZPLLL · 发表于 Post on 2023-8-7 22:50:58

zjxitcc 发表于 2023-8-7 21:45
是因为赝势的原因。

至于是否可信，还有其他因素要考虑。比如

您好，再请教一下。赝势基组与其它基组的单点能差别这么大，那么在计算单点能的过程中，是否能选用赝势基组来进行计算呢？

sobereva · 发表于 Post on 2023-8-8 01:53:25

不要自己在标题里手写[gaussian]这种标签，这次给你改了，以后注意

sobereva · 发表于 Post on 2023-8-8 02:12:27

ZPLLL 发表于 2023-8-7 22:50
您好，再请教一下。赝势基组与其它基组的单点能差别这么大，那么在计算单点能的过程中，是否能选用赝势基 ...

赝势基组不明确描述内核电子，而内核电子的能量又远比价层低得多，所以用不用赝势时能量的数量级都有很大不同

ZPLLL · 发表于 Post on 2023-8-8 10:58:49

sobereva 发表于 2023-8-8 02:12
赝势基组不明确描述内核电子，而内核电子的能量又远比价层低得多，所以用不用赝势时能量的数量级都有很大 ...

理解了，感谢老师的解答。

ZPLLL · 发表于 Post on 2023-8-8 10:59:56

sobereva 发表于 2023-8-8 01:53
不要自己在标题里手写[gaussian]这种标签，这次给你改了，以后注意

好的老师，下次我注意

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

[Gaussian/gview] 结构优化和单点能计算能量差别悬殊