计算化学公社

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

[新手求助] 请问不同基组优化出来的几何结构,键长键角相差多少算基组影响很小

[复制链接 Copy URL]

97

帖子

0

威望

569

eV
积分
666

Level 4 (黑子)

跳转到指定楼层 Go to specific reply
楼主
各位老师好。假如我用cc-pvdz和cc-pvtz优化几何结构,来判断选用哪个基组。请问二者优化出来的键长键角等参数,相差值在多大范围以内,就可以算基组的影响较小呢?
以及还想问各位老师一个问题,计算绝热能的时候,计算值与实验值相差在多少以内是比较好的呢,假如以kcal/mol为单位?
谢谢各位老师。

1万

帖子

0

威望

9902

eV
积分
22161

Level 6 (一方通行)

2#
发表于 Post on 2023-2-15 05:24:12 | 只看该作者 Only view this author
视你需要的精度而定,没有固定的答案。
和实验值比的时候,如果差别小于实验误差范围,或远小于你的理论计算结果和实验结果之差,那么这个差别一定是可以算较小的;如果差别大到会改变你文章的结论,那一定是算较大的;如果差别大于实验误差范围,却又不足以导致文章结论改变,那么没有一定的规则,如果你机时紧张就用那个便宜的理论级别,否则就用精确的那个理论级别。
举例:你用CASPT2/cc-pVDZ和CASPT2/cc-pVTZ优化一个键长,结果分别为R1和R2(埃);实验的键长为R3。
如果R1=1.200, R2=1.205, R3=1.205±0.010,则基组影响较小,因为两者的差别小于实验误差;
如果R1=1.200, R2=1.205, R3=1.255±0.001,则基组影响较小,因为两者的差别远小于理论方法的误差,此时应该考虑换用比CASPT2更精确的方法,或增大活性空间等。
如果R1=1.200, R2=1.250, 且计算的目的在于判断两个互相矛盾的实验值R3a=1.205±0.001和R3b=1.245±0.001哪个更合理,那么基组影响较大,因为改变理论方法会改变文章结论。
如果R1=1.200, R2=1.210, 仍要判断以上两个实验值哪个更合理,则应视机时情况等综合判断,不同人答案不一样。因为虽然两者的差别既不远小于理论方法误差,又不小于实验误差,但也不影响文章结论。
以上所有讨论不仅适用于键长键角,也适用于能量以及其他几乎任何物理量
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
本团队长期招收研究生,有意者可私信联系

6万

帖子

99

威望

6万

eV
积分
125234

管理员

公社社长

3#
发表于 Post on 2023-2-15 06:31:35 | 只看该作者 Only view this author
做这种测试毫无意义,看下文
谈谈量子化学研究中什么样的benchmark才有意义
http://sobereva.com/554http://bbs.keinsci.com/thread-17826-1-1.html

从基组造成的误差角度来讲,有最基本基组常识的人都知道cc-pVTZ的误差明显小于cc-pVDZ,不需要做任何测试。即便你对于某个体系,由于理论方法、基组、溶剂效应、计算条件与实验条件可比性等等各种因素的误差抵消原因,十分巧合地cc-pVDZ算出来的几何结构表面上比cc-pVTZ和实验更相符,也并不能说明任何问题,因为没有任何普遍意义。而且如果以此为依据选了cc-pVDZ做其它计算,被审稿人知道了大概率要挨批。
北京科音自然科学研究中心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

97

帖子

0

威望

569

eV
积分
666

Level 4 (黑子)

4#
 楼主 Author| 发表于 Post on 2023-2-15 15:49:06 | 只看该作者 Only view this author
wzkchem5 发表于 2023-2-15 05:24
视你需要的精度而定,没有固定的答案。
和实验值比的时候,如果差别小于实验误差范围,或远小于你的理论计 ...

好的,谢谢老师,明白了

97

帖子

0

威望

569

eV
积分
666

Level 4 (黑子)

5#
 楼主 Author| 发表于 Post on 2023-2-15 15:49:25 | 只看该作者 Only view this author
sobereva 发表于 2023-2-15 06:31
做这种测试毫无意义,看下文
谈谈量子化学研究中什么样的benchmark才有意义
http://sobereva.com/554(ht ...

好的,谢谢老师

本版积分规则 Credits rule

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

GMT+8, 2026-2-28 14:04 , Processed in 1.054354 second(s), 20 queries , Gzip On.

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