计算化学公社

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

[新手求助] 用不同基组算单点能时个别数据能量差距较大

[复制链接 Copy URL]

13

帖子

0

威望

217

eV
积分
230

Level 3 能力者

用ub3lyp/6-31g(d,p)优化,用ub3lyp/6-311++g(d,p)计算单点能,后又提高精度换成了ub3lyp def2tzvp计算单点能。发现换了基组以后计算的能量大部分没有差别,但个别数据出现了14kcal/mol的能量差,后又用ub3lyp def2qzvpp对其进行计算,发现结果与ub3lyp/6-311++g(d,p)计算的能量一致,请问这种情况该怎么办呢?




6万

帖子

99

威望

6万

eV
积分
125156

管理员

公社社长

2#
发表于 Post on 2024-7-1 17:58:06 | 只看该作者 Only view this author
不说具体算的东西问这个毫无意义
倘若是必须带弥散的情况,更贵的def2-TZVP结果很可能还不如6-311++g(d,p)
仔细看
谈谈弥散函数和“月份”基组
http://sobereva.com/119

北京科音自然科学研究中心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

1万

帖子

0

威望

9888

eV
积分
22133

Level 6 (一方通行)

3#
发表于 Post on 2024-7-1 19:50:52 | 只看该作者 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?user=XW6C6eQAAAAJ
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1133/1776.htm
GitHub:https://github.com/wzkchem5
本团队长期招收研究生,有意者可私信联系

13

帖子

0

威望

217

eV
积分
230

Level 3 能力者

4#
 楼主 Author| 发表于 Post on 2024-7-1 20:31:33 | 只看该作者 Only view this author
本帖最后由 Jiangxue 于 2024-7-1 20:36 编辑
sobereva 发表于 2024-7-1 17:58
不说具体算的东西问这个毫无意义
倘若是必须带弥散的情况,更贵的def2-TZVP结果很可能还不如6-311++g(d,p) ...

计算的是一个铁卟啉环和药物分子的加成,输入和输出文件如下,麻烦老师帮我看看
C:\Users\Jiangxue\Desktop

13

帖子

0

威望

217

eV
积分
230

Level 3 能力者

5#
 楼主 Author| 发表于 Post on 2024-7-1 20:36:04 | 只看该作者 Only view this author
本帖最后由 Jiangxue 于 2024-7-1 20:49 编辑
Jiangxue 发表于 2024-7-1 20:31
计算的是一个铁卟啉环和药物分子的加成,输入和输出文件如下,麻烦老师帮我看看

输入文件输出文件[/ur ...



4s2FBJTC4epo-4ts1-RCspdef2.out (123.73 KB, 下载次数 Times of downloads: 0) 2FBJTC4NIF-4ts2-RCspdef2.gjf (6.13 KB, 下载次数 Times of downloads: 10)
4s2FBJTC4epo-ts1irc-RCsp1.log (102.13 KB, 下载次数 Times of downloads: 2) 4s2FBJTC4epo-ts1irc-PC1sp.log (102.22 KB, 下载次数 Times of downloads: 0)
4s2FBJTC4epo-4ts1-RCspqzvp.log (231.83 KB, 下载次数 Times of downloads: 0) 4s2FBJTC4epo-4ts1-PCspqzvp.log (231.96 KB, 下载次数 Times of downloads: 0)

6万

帖子

99

威望

6万

eV
积分
125156

管理员

公社社长

6#
发表于 Post on 2024-7-2 19:33:56 | 只看该作者 Only view this author

光从gjf文件上没看出有问题
普通泛函结合def2-TZVP算这类体系的电子能量完全够用,没必要考虑其它的
北京科音自然科学研究中心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

13

帖子

0

威望

217

eV
积分
230

Level 3 能力者

7#
 楼主 Author| 发表于 Post on 2024-7-3 16:54:21 | 只看该作者 Only view this author
wzkchem5 发表于 2024-7-1 19:50
检查是不是波函数收敛到了错误的地方

请问怎么检查是不是波函数收敛到了错误的地方呢?是检查计算的opt输出结果是否收敛吗?我的计算结果opt输出能量曲线图降低到了最低点,收敛结果有3个yes,用6-311++g(d,p)计算的结果和def2qzvpp保持一致,但用def2tzvp计算出的结果差了很多,反复检查了输入文件和初始结构,都没发现问题,不知道该如何解决这一问题

1万

帖子

0

威望

9888

eV
积分
22133

Level 6 (一方通行)

8#
发表于 Post on 2024-7-3 20:48:24 | 只看该作者 Only view this author
Jiangxue 发表于 2024-7-3 09:54
请问怎么检查是不是波函数收敛到了错误的地方呢?是检查计算的opt输出结果是否收敛吗?我的计算结果opt输 ...

不是,是检查最后一帧的布居分析(电荷和自旋)以及前线轨道,是否符合化学常识
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?user=XW6C6eQAAAAJ
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1133/1776.htm
GitHub:https://github.com/wzkchem5
本团队长期招收研究生,有意者可私信联系

本版积分规则 Credits rule

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

GMT+8, 2026-2-24 23:28 , Processed in 0.176687 second(s), 23 queries , Gzip On.

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