计算化学公社

 找回密码 Forget password
 注册 Register
Views: 3221|回复 Reply: 13
打印 Print 上一主题 Last thread 下一主题 Next thread

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

[复制链接 Copy URL]

26

帖子

0

威望

439

eV
积分
465

Level 3 能力者

在优化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)的原因吗?




1万

帖子

0

威望

8993

eV
积分
20741

Level 6 (一方通行)

2#
发表于 Post on 2023-8-7 21:45:03 | 只看该作者 Only view this author
Zikuan Wang
山东大学光学高等研究中心 研究员
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员
Google Scholar: https://scholar.google.com/citations?hl=zh-CN&user=XW6C6eQAAAAJ&view_op=list_works&sortby=pubdate
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1034/1702.htm
本团队长期招收研究生,有意者可私信联系

4104

帖子

4

威望

8869

eV
积分
13053

Level 6 (一方通行)

MOKIT开发者

3#
发表于 Post on 2023-8-7 21:45:34 | 只看该作者 Only view this author
是因为赝势的原因。

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

26

帖子

0

威望

439

eV
积分
465

Level 3 能力者

4#
 楼主 Author| 发表于 Post on 2023-8-7 21:52:34 | 只看该作者 Only view this author
wzkchem5 发表于 2023-8-7 21:45
那叫SDD不叫SSD
参见https://mattermodeling.stackexchange.com/questions/9961/effective-core-potential ...

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

26

帖子

0

威望

439

eV
积分
465

Level 3 能力者

5#
 楼主 Author| 发表于 Post on 2023-8-7 22:02:54 | 只看该作者 Only view this author
zjxitcc 发表于 2023-8-7 21:45
是因为赝势的原因。

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

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

26

帖子

0

威望

439

eV
积分
465

Level 3 能力者

6#
 楼主 Author| 发表于 Post on 2023-8-7 22:50:58 | 只看该作者 Only view this author
zjxitcc 发表于 2023-8-7 21:45
是因为赝势的原因。

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

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

6万

帖子

99

威望

5万

eV
积分
120110

管理员

公社社长

7#
发表于 Post on 2023-8-8 01:53:25 | 只看该作者 Only view this author
不要自己在标题里手写[gaussian]这种标签,这次给你改了,以后注意
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

6万

帖子

99

威望

5万

eV
积分
120110

管理员

公社社长

8#
发表于 Post on 2023-8-8 02:12:27 | 只看该作者 Only view this author
ZPLLL 发表于 2023-8-7 22:50
您好,再请教一下。赝势基组与其它基组的单点能差别这么大,那么在计算单点能的过程中,是否能选用赝势基 ...

赝势基组不明确描述内核电子,而内核电子的能量又远比价层低得多,所以用不用赝势时能量的数量级都有很大不同
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

26

帖子

0

威望

439

eV
积分
465

Level 3 能力者

9#
 楼主 Author| 发表于 Post on 2023-8-8 10:58:49 | 只看该作者 Only view this author
sobereva 发表于 2023-8-8 02:12
赝势基组不明确描述内核电子,而内核电子的能量又远比价层低得多,所以用不用赝势时能量的数量级都有很大 ...

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

26

帖子

0

威望

439

eV
积分
465

Level 3 能力者

10#
 楼主 Author| 发表于 Post on 2023-8-8 10:59:56 | 只看该作者 Only view this author
sobereva 发表于 2023-8-8 01:53
不要自己在标题里手写[gaussian]这种标签,这次给你改了,以后注意

好的老师,下次我注意

58

帖子

0

威望

369

eV
积分
427

Level 3 能力者

11#
发表于 Post on 2025-5-29 15:18:07 | 只看该作者 Only view this author
sobereva 发表于 2023-8-8 02:12
赝势基组不明确描述内核电子,而内核电子的能量又远比价层低得多,所以用不用赝势时能量的数量级都有很大 ...

老师,那我比较Yb和Li这两个原子的单点能的时候,Li是不是也要采用赝势基组呀

168

帖子

1

威望

2217

eV
积分
2405

Level 5 (御坂)

黄昏粼粼

12#
发表于 Post on 2025-5-29 17:10:36 | 只看该作者 Only view this author
陈AG 发表于 2025-5-29 15:18
老师,那我比较Yb和Li这两个原子的单点能的时候,Li是不是也要采用赝势基组呀

比较Yb和Li的单点能有何意义?
愿乘长风,破万里浪

6万

帖子

99

威望

5万

eV
积分
120110

管理员

公社社长

13#
发表于 Post on 2025-5-29 17:16:07 | 只看该作者 Only view this author
陈AG 发表于 2025-5-29 15:18
老师,那我比较Yb和Li这两个原子的单点能的时候,Li是不是也要采用赝势基组呀


元素都不同根本没有可比性
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

58

帖子

0

威望

369

eV
积分
427

Level 3 能力者

14#
发表于 Post on 2025-5-29 17:34:13 | 只看该作者 Only view this author
pikachuupup 发表于 2025-5-29 17:10
比较Yb和Li的单点能有何意义?

老师,我主要想看这个Yb金属配合物结构优化和单点能基组如何选用
http://bbs.keinsci.com/forum.php ... 3&fromuid=35653
(出处: 计算化学公社)

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2025-8-14 17:15 , Processed in 0.164329 second(s), 23 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list