计算化学公社

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

[新手求助] 在SMD模型下使用M06-2X/6-31++G**计算相互作用能是否需要考虑BSSE校正能

[复制链接 Copy URL]

22

帖子

0

威望

262

eV
积分
284

Level 3 能力者

跳转到指定楼层 Go to specific reply
楼主
我想计算SMD模型下溶质和溶剂分子之间的相互作用能,首先用B3LYP-D3/6-311+G**进行结构优化,获得气相下的最优结构。再在气相最优构型下,使用SMD模型和M06-2X/6-31++G**计算溶质与溶剂之间的相互作用能。在计算相互作用能的过程中我需要加入BSSE校正能?还有在优化和计算能量过程中选用的计算方法合理吗?望各位赐教,谢谢!

28

帖子

0

威望

220

eV
积分
248

Level 3 能力者

14#
发表于 Post on 2024-10-31 17:29:09 | 只看该作者 Only view this author
sobereva 发表于 2024-10-31 13:14
下文明确说了,即便是阴离子体系用SMD计算溶解自由能也不该带弥散函数
谈谈隐式溶剂模型下溶解自由能和 ...

感谢sob老师,之前是我读的不仔细

5万

帖子

99

威望

5万

eV
积分
112351

管理员

公社社长

13#
发表于 Post on 2024-10-31 13:14:38 | 只看该作者 Only view this author
qaqfdmmj 发表于 2024-10-31 10:43
sob老师,这里如果计算的是阴离子体系的话,这里6-31G*是要改成6-31+G*吗

下文明确说了,即便是阴离子体系用SMD计算溶解自由能也不该带弥散函数
谈谈隐式溶剂模型下溶解自由能和体系自由能的计算
http://sobereva.com/327http://bbs.keinsci.com/thread-3345-1-1.html
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

28

帖子

0

威望

220

eV
积分
248

Level 3 能力者

12#
发表于 Post on 2024-10-31 10:43:11 | 只看该作者 Only view this author
sobereva 发表于 2016-11-30 17:59
你要计算溶剂化能,就用B3LYP-D3/6-311G**优化得到的结构,在M052X/6-31G*下用SMD溶剂模型来算就完了,这个 ...

sob老师,这里如果计算的是阴离子体系的话,这里6-31G*是要改成6-31+G*吗

22

帖子

0

威望

262

eV
积分
284

Level 3 能力者

11#
 楼主 Author| 发表于 Post on 2017-3-11 18:41:01 | 只看该作者 Only view this author

好的,谢谢您!

5万

帖子

99

威望

5万

eV
积分
112351

管理员

公社社长

10#
发表于 Post on 2017-3-11 18:35:27 | 只看该作者 Only view this author
C41014007 发表于 2017-3-11 16:56
请问在用M052X/6-31G*方法计算完溶剂化自由能之后,能够直接用该方法下计算溶剂化自由能时得到的波函数文 ...

可以。
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

22

帖子

0

威望

262

eV
积分
284

Level 3 能力者

9#
 楼主 Author| 发表于 Post on 2017-3-11 16:56:16 | 只看该作者 Only view this author
sobereva 发表于 2016-11-30 17:59
你要计算溶剂化能,就用B3LYP-D3/6-311G**优化得到的结构,在M052X/6-31G*下用SMD溶剂模型来算就完了,这个 ...

请问在用M052X/6-31G*方法计算完溶剂化自由能之后,能够直接用该方法下计算溶剂化自由能时得到的波函数文件进行波函数分析(进行静电势、AIM、RDG分析等)吗?

5万

帖子

99

威望

5万

eV
积分
112351

管理员

公社社长

8#
发表于 Post on 2016-12-1 14:26:22 | 只看该作者 Only view this author
C41014007 发表于 2016-11-30 20:49
因为之前看有文献中使用SMD模型计算溶剂化能的时候使用的是M06-2X/6-31++G**,所以就改了


看了SMD原文里的对比你就知道了
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

22

帖子

0

威望

262

eV
积分
284

Level 3 能力者

7#
 楼主 Author| 发表于 Post on 2016-11-30 20:49:35 | 只看该作者 Only view this author
happyrabbit 发表于 2016-11-30 15:40
为何要改变计算方法?B3LYP-D3计算若相互作用很好啊,单点计算的时候增加弥散基组就行

因为之前看有文献中使用SMD模型计算溶剂化能的时候使用的是M06-2X/6-31++G**,所以就改了

5万

帖子

99

威望

5万

eV
积分
112351

管理员

公社社长

6#
发表于 Post on 2016-11-30 17:59:53 | 只看该作者 Only view this author
你要计算溶剂化能,就用B3LYP-D3/6-311G**优化得到的结构,在M052X/6-31G*下用SMD溶剂模型来算就完了,这个级别是最适合结合SMD算溶解自由能的级别
用隐式溶剂模型算溶解自由能和BSSE完全没有丝毫联系
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

41

帖子

0

威望

829

eV
积分
870

Level 4 (黑子)

5#
发表于 Post on 2016-11-30 15:40:48 | 只看该作者 Only view this author
为何要改变计算方法?B3LYP-D3计算若相互作用很好啊,单点计算的时候增加弥散基组就行

22

帖子

0

威望

262

eV
积分
284

Level 3 能力者

4#
 楼主 Author| 发表于 Post on 2016-11-30 13:41:12 | 只看该作者 Only view this author
liyuanhe211 发表于 2016-11-30 11:52
仿佛没一句对的→_→

SMD计算的是(平均)溶剂化自由能,不是“溶质与溶剂之间的相互作用能”。后者应通 ...

不好意思,表达有误。是想计算溶剂化自由能。因为之前已用B3LYP-D3/6-311G**方法对气相结构进行优化,在cosmo模型下使用PWPB95-D3/def2-QZVPP进行了能量计算。但后来了解到用SMD模型更为合适,在SMD模型下使用M06-2X/6-31G*计算能量也可以,所以就想用原来B3LYP-D3/6-311G**优化出来的气相结构,在SMD模型下用M06-2X/6-31++G**进行能量计算。现在看来这样做是错的了。

3098

帖子

28

威望

1万

eV
积分
16892

Level 6 (一方通行)

3#
发表于 Post on 2016-11-30 11:52:30 | 只看该作者 Only view this author
仿佛没一句对的→_→

SMD计算的是(平均)溶剂化自由能,不是“溶质与溶剂之间的相互作用能”。后者应通过把溶剂分子显式的放在体系里,计算弱相互作用实现。先用3-zeta基组优化再用2-zeta计算能量毫无道理,包括方法、基组、BSSE的问题看 sobereva.com/214。

如果只是计算溶剂化自由能,看 sobereva.com/327

1102

帖子

18

威望

6643

eV
积分
8105

Level 6 (一方通行)

計算化学の社畜

2#
发表于 Post on 2016-11-30 11:45:40 | 只看该作者 Only view this author
本帖最后由 冰释之川 于 2016-11-30 11:47 编辑

优化用了3-zeta基组;算相互作用能为何用2-zeta基组?6-31++G**基组还不如用6-311+G**,给氢加弥散绝大多数时候没啥意义。
关于BSSE,参考《谈谈BSSE校正与Gaussian对它的处理》http://www.sobereva.com/46  BSSE要不要做,要看你的研究体系和研究内容,还有就是基组本身质量(完备性程度)
Stand on the shoulders of giants

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

GMT+8, 2024-11-23 10:01 , Processed in 0.187414 second(s), 22 queries , Gzip On.

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