计算化学公社

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

[新手求助] 想请教关于解释计算结果的问题

[复制链接 Copy URL]

201

帖子

2

威望

1354

eV
积分
1595

Level 5 (御坂)

跳转到指定楼层 Go to specific reply
楼主
本帖最后由 steven 于 2018-10-31 17:38 编辑

各位老师同学好。实验上常温做出来一个产物,是1位产物。我做计算研究1位和三位的竞争问题,但是本人菜鸟一个,这是第一次要发文章。
现在我有两个问题想听取一下大家的建议:


1,我找到的最低的1位路径势垒是 25 kcal/mol,3位是 25.6 kcal/mol.  都超过了21kcal/mol  这俩决速步涉及的TS都是一个脱H的H转移过程。 计算级别M06/6-31g*  溶解自由能直接是用M06/6-31g* SCRF(SMD,solvent=AceticAcid) 得到的单点能加到之前的气相下自由能校正量得到的,相当于用再M06算的溶解自由能。 投稿的话有没有什么比较好的方式能说一下由于计算误差什么的超过了21kcal/mol? 还是直接不提这个事也行?

2,如果用这个势垒算反应速率的话一定很小,而且反应步骤很多,我觉得误差会很大所以没想算速率。 按照决速步来看1位只比3位低了0.6kcal/mol  这样用于说明1位很有优势会有啥问题吗?(或者也像1一样解释一下,比如说可能是计算误差导致差这么小,实际差的可能更大?)

谢谢大家。


561

帖子

0

威望

3410

eV
积分
3971

Level 5 (御坂)

2#
发表于 Post on 2018-10-31 16:29:59 | 只看该作者 Only view this author
首先我就很好奇为啥你用M06而不是M062X…

201

帖子

2

威望

1354

eV
积分
1595

Level 5 (御坂)

3#
 楼主 Author| 发表于 Post on 2018-10-31 16:35:05 | 只看该作者 Only view this author
k64_cc 发表于 2018-10-31 16:29
首先我就很好奇为啥你用M06而不是M062X…

含过渡金属体系

6万

帖子

99

威望

6万

eV
积分
125145

管理员

公社社长

4#
发表于 Post on 2018-10-31 20:15:13 | 只看该作者 Only view this author
基组太烂,要定量靠谱计算势垒,起码得TZVP。用这样的级别算单点8成要挨批,起码涉及反应的那些原子基组得大点。另外计算量允许的情况下,尝试双杂化泛函计算。当前问题的研究对于相对自由能很敏感,相差本来就不大,就更需要单点级别的精度应尽量好
北京科音自然科学研究中心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

201

帖子

2

威望

1354

eV
积分
1595

Level 5 (御坂)

5#
 楼主 Author| 发表于 Post on 2018-10-31 20:27:31 | 只看该作者 Only view this author
本帖最后由 steven 于 2018-10-31 20:45 编辑
sobereva 发表于 2018-10-31 20:15
基组太烂,要定量靠谱计算势垒,起码得TZVP。用这样的级别算单点8成要挨批,起码涉及反应的那些原子基组得 ...

谢谢老师的意见。那我基组用TZVP的话,双杂化泛函您推荐用什么合适一些?(最好是Gaussion09 A.01支持的..)  而且用M06/6-31g*算溶解自由能,应该还合适把?体系是含Pd的有机金属配合物。
比较尴尬的是,因为反应路径算了很多想要重新算一遍计算量极大,而且期限将至不知道能不能算完,而且假如级别改了以后,相对能量仍然是这样变化不大的话,针对我的1和2两点能否给审稿人一个什么解释?  谢谢老师。

6万

帖子

99

威望

6万

eV
积分
125145

管理员

公社社长

6#
发表于 Post on 2018-10-31 23:13:36 | 只看该作者 Only view this author
steven 发表于 2018-10-31 20:27
谢谢老师的意见。那我基组用TZVP的话,双杂化泛函您推荐用什么合适一些?(最好是Gaussion09 A.01支持的.. ...

溶解自由能那个级别没问题
双杂化用B2PLYP-D3(BJ)。此时尽量用不低于def2-TZVP的基组

不需要重新优化结构,只需要算个高精度单点,不会花多少时间。而且用ORCA的话跑双杂化速度极快
先算完了看情况再说
北京科音自然科学研究中心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

201

帖子

2

威望

1354

eV
积分
1595

Level 5 (御坂)

7#
 楼主 Author| 发表于 Post on 2018-11-1 16:58:40 | 只看该作者 Only view this author
sobereva 发表于 2018-10-31 23:13
溶解自由能那个级别没问题
双杂化用B2PLYP-D3(BJ)。此时尽量用不低于def2-TZVP的基组

谢谢老师。之前的都是用Gaussian09算的(而且我没用过ORCA怕折腾不明白),如果双杂化用不起的话,单点计算从M06/6-31g*换成 M06L/def2TZVP 或 BP86/def2TZVP 您看行吗?

6万

帖子

99

威望

6万

eV
积分
125145

管理员

公社社长

8#
发表于 Post on 2018-11-1 16:59:34 | 只看该作者 Only view this author
steven 发表于 2018-11-1 16:58
谢谢老师。之前的都是用Gaussian09算的(而且我没用过ORCA怕折腾不明白),如果双杂化用不起的话,单点计 ...

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

201

帖子

2

威望

1354

eV
积分
1595

Level 5 (御坂)

9#
 楼主 Author| 发表于 Post on 2018-11-1 20:36:37 | 只看该作者 Only view this author
本帖最后由 steven 于 2018-11-1 22:15 编辑

非常非常感谢老师。
1.我试用了几个泛函,最后用了  M06L/def2-TZVP  EmpiricalDispersion=GD3  计算单点。总体上自由能垒确实下降了。
1位决速步能垒17.9 kcal/mol      3位决速步能垒 19.1 kcal/mol
都小于21了,但是还是也就差了1kcal/mol,  我把这个作为论据老师您认为可行吗?(实验为室温,产物都是1位)

2.另外若被审稿人问到为何优化振动分析和溶解自由能都是用的M06泛函(ultrafine,无D3矫正), 而单点除却基组的不同外,把泛函改成了 M06L (没加ultrafine,加了D3矫正) ,请教老师这该怎样回答合适?  (加没加ultrafine应该没事,但是泛函本身变了,而且多了D3矫正)


3.优化和振动分析没加色散矫正,我觉得M06不加不会成为问题吧?

谢谢老师。


6万

帖子

99

威望

6万

eV
积分
125145

管理员

公社社长

10#
发表于 Post on 2018-11-2 00:31:37 | 只看该作者 Only view this author
steven 发表于 2018-11-1 20:36
非常非常感谢老师。
1.我试用了几个泛函,最后用了  M06L/def2-TZVP  EmpiricalDispersion=GD3  计算单 ...

虽然1kcal/mol看起来不大,但是对k的影响很大,可以作为证据

优化和振动分析最好也用M06L,或者算势垒的时候用M06,应当统一,免得引起麻烦。或者势垒用双杂化来算
参看
量子化学研究中切换泛函应当注意的问题
http://sobereva.com/415

加不加色散校正应统一,见上文
北京科音自然科学研究中心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

201

帖子

2

威望

1354

eV
积分
1595

Level 5 (御坂)

11#
 楼主 Author| 发表于 Post on 2018-11-2 01:43:36 | 只看该作者 Only view this author
本帖最后由 steven 于 2018-11-2 01:51 编辑
sobereva 发表于 2018-11-2 00:31
虽然1kcal/mol看起来不大,但是对k的影响很大,可以作为证据

优化和振动分析最好也用M06L,或者算势垒 ...

完了,都用M06L优化振动分析一遍应该来不及了。只能再用M06重新算势垒,就用fine是不是也行呢老师?优化振动分析是ultrafine,我试了如果单点也用ultrafine,1和3位势垒差的更小…

或者我现在可以把M06L的D3矫正去掉。这样被审稿人问起的概率大吗…?

感谢老师!

6万

帖子

99

威望

6万

eV
积分
125145

管理员

公社社长

12#
发表于 Post on 2018-11-2 05:26:12 | 只看该作者 Only view this author
steven 发表于 2018-11-2 01:43
完了,都用M06L优化振动分析一遍应该来不及了。只能再用M06重新算势垒,就用fine是不是也行呢老师?优化 ...

原理上单点应当也用ultrafine
D3要用就所有过程都用,要么就都别用
北京科音自然科学研究中心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

201

帖子

2

威望

1354

eV
积分
1595

Level 5 (御坂)

13#
 楼主 Author| 发表于 Post on 2018-11-2 09:17:59 | 只看该作者 Only view this author
sobereva 发表于 2018-11-2 05:26
原理上单点应当也用ultrafine
D3要用就所有过程都用,要么就都别用

好的,谢谢老师的耐心指导。

201

帖子

2

威望

1354

eV
积分
1595

Level 5 (御坂)

14#
 楼主 Author| 发表于 Post on 2018-11-2 17:08:33 | 只看该作者 Only view this author
steven 发表于 2018-11-2 09:17
好的,谢谢老师的耐心指导。

再次打扰老师。。
1.  我几何优化和振动分析是 M06/6-31G* int=ultrafine    溶解自由能是  M06/6-31G* int=ultrafine SCRF(SMD,solvent=AceticAcid)。 我要做反应中一个中间体的波函数分析,是不是应该用 M06/6-31G* int=ultrafine SCRF(SMD,solvent=AceticAcid) 的波函数? 谢谢老师!

2. 我把泛函和D3矫正都统一了,这下 1位 和 3位 就差了 0.6 kcal/mol   应该作为依据也可以吧?

感谢老师。

6万

帖子

99

威望

6万

eV
积分
125145

管理员

公社社长

15#
发表于 Post on 2018-11-2 17:49:29 | 只看该作者 Only view this author
steven 发表于 2018-11-2 17:08
再次打扰老师。。
1.  我几何优化和振动分析是 M06/6-31G* int=ultrafine    溶解自由能是  M06/6-31G*  ...

1 是
2 可以作为依据,但是对于普通泛函来说,能正确分辨0.6 kcal/mol数量级的差异并不容易,看审稿人信不信了。如果能做个双杂化计算进一步验证的话更有说服力。
北京科音自然科学研究中心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

本版积分规则 Credits rule

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

GMT+8, 2026-2-22 09:15 , Processed in 0.165055 second(s), 20 queries , Gzip On.

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