计算化学公社

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

[新手求助] 求助关于计算BSSE矫正能!

[复制链接 Copy URL]

105

帖子

0

威望

867

eV
积分
974

Level 4 (黑子)

跳转到指定楼层 Go to specific reply
楼主
大家好,想请问一下,我想计算A、B相互作用能,想考虑BSSE矫正,不知道是不是下面这样做呢?1、首先在某一个基组Q下opt+freq优化A,没有虚频;
2、接着同一个基组Q下opt+freq优化B,同样没有虚频;
3、还是统一基组Q下将优化完的A和B几个构型凑在一块(相互作用的部位)     opt+freq  并添加counterpoise=2
4、然后;作用能= 矫正的二聚体corrected energy - 矫正的片段1(A)能量 - 矫正的片段(B)的能量
是这样的吗?第1和2步需要吗??

评分 Rate

参与人数
Participants 1
eV +1 收起 理由
Reason
teao + 1 赞!

查看全部评分 View all ratings

6万

帖子

99

威望

5万

eV
积分
120165

管理员

公社社长

2#
发表于 Post on 2014-12-8 22:06:45 | 只看该作者 Only view this author
一般计算考虑了BSSE的相互作用能的时候,就是这么几步
(1)对AB复合物进行优化。(优化过程中不建议考虑BSSE,不仅意义不大,而且由于没有解析梯度,耗时增加甚巨。)
(2)从复合物结构中提取A的结构,算A的单点能E_A
(3)从复合物结构中提取B的结构,算B的单点能E_B
(4)计算考虑了BSSE校正后的复合物能量E_AB_corrected
然后E_AB_corrected - E_A - E_B即是相互作用能。

下面这个帖子有详细讨论。实际上Gaussian里用counterpoise关键词对优化后的复合物做单点能计算的时候,已经把(2)、(3)、(4)一次性全都计算出来了。
谈谈BSSE校正与Gaussian对它的处理
http://sobereva.com/46

在计算A和B的能量前是否对A和B单独进行优化看具体考察什么。比如对于A,它在优化后结构下的能量和在复合物结构下的能量之差是所谓的变形能,体现了与B相互作用前结构发生变形造成的能量改变。多数相互作用能的研究都是不考虑变形能的。但如果要计算诸如反应热,那是肯定要考虑的。
北京科音自然科学研究中心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

145

帖子

0

威望

290

eV
积分
435

Level 3 能力者

3#
发表于 Post on 2014-12-9 20:10:07 | 只看该作者 Only view this author
在计算ab复合物结构的时候加上关键词counter再定义分子片段,这个也可以吧?我一直是这样做的

6万

帖子

99

威望

5万

eV
积分
120165

管理员

公社社长

4#
发表于 Post on 2014-12-9 20:51:27 | 只看该作者 Only view this author
小菜鸟求关照 发表于 2014-12-9 20:10
在计算ab复合物结构的时候加上关键词counter再定义分子片段,这个也可以吧?我一直是这样做的

上面回帖中我说的“实际上Gaussian里用counterpoise关键词对优化后的复合物做单点能计算的时候,已经把(2)、(3)、(4)一次性全都计算出来了。”就是你说的这个过程
北京科音自然科学研究中心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

145

帖子

0

威望

290

eV
积分
435

Level 3 能力者

5#
发表于 Post on 2014-12-10 09:17:16 | 只看该作者 Only view this author
sobereva 发表于 2014-12-9 20:51
上面回帖中我说的“实际上Gaussian里用counterpoise关键词对优化后的复合物做单点能计算的时候,已经把(2 ...

嗯呢 明白啦

30

帖子

0

威望

153

eV
积分
183

Level 3 能力者

6#
发表于 Post on 2016-9-28 16:42:49 | 只看该作者 Only view this author
sobereva 发表于 2014-12-8 22:06
一般计算考虑了BSSE的相互作用能的时候,就是这么几步
(1)对AB复合物进行优化。(优化过程中不建议考虑BSS ...

老师  你说的E_AB_corrected - E_A - E_B即是相互作用能   这里的  E-A 、E-B是不是就是输出文件里的  E(MCBS)1和E(MCBS)2   可以这样理解吗?   这样减出来的值就是文献里常说的Eint吗?

6万

帖子

99

威望

5万

eV
积分
120165

管理员

公社社长

7#
发表于 Post on 2016-9-28 16:50:18 | 只看该作者 Only view this author
zhangyina 发表于 2016-9-28 16:42
老师  你说的E_AB_corrected - E_A - E_B即是相互作用能   这里的  E-A 、E-B是不是就是输出文件里的  E ...




一般直接取counterpoise任务最后输出的complexation energy =      xxx kcal/mole (corrected)就完了,就是这个值,省得还自己加减了
北京科音自然科学研究中心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

30

帖子

0

威望

153

eV
积分
183

Level 3 能力者

8#
发表于 Post on 2016-9-28 17:10:51 | 只看该作者 Only view this author
sobereva 发表于 2016-9-28 16:50


一般直接取counterpoise任务最后输出的complexation energy =      xxx kcal/mole (corrected)就 ...

我用高斯09 a.02  做的输出结果  只有BSSE  和  校正后的复合物的能量  

30

帖子

0

威望

153

eV
积分
183

Level 3 能力者

9#
发表于 Post on 2016-9-28 17:25:04 | 只看该作者 Only view this author
sobereva 发表于 2016-9-28 16:50


一般直接取counterpoise任务最后输出的complexation energy =      xxx kcal/mole (corrected)就 ...

用mp2计算时   E(MCBS 1)  E(MCBS 2)应该是EUMP2后面的值对吧?那就算Eint时  也减得是这个数值吗?好像不对  那应该 是什么?

6万

帖子

99

威望

5万

eV
积分
120165

管理员

公社社长

10#
发表于 Post on 2016-9-28 19:07:04 | 只看该作者 Only view this author
zhangyina 发表于 2016-9-28 17:25
用mp2计算时   E(MCBS 1)  E(MCBS 2)应该是EUMP2后面的值对吧?那就算Eint时  也减得是这个数值吗? ...

MP2计算当然得取EUMP2的值,scf done后面是HF能量
北京科音自然科学研究中心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

30

帖子

0

威望

153

eV
积分
183

Level 3 能力者

11#
发表于 Post on 2016-9-28 19:41:14 | 只看该作者 Only view this author
本帖最后由 zhangyina 于 2016-9-28 21:54 编辑
sobereva 发表于 2016-9-28 19:07
MP2计算当然得取EUMP2的值,scf done后面是HF能量

E_AB_corrected - E_A - E_B即是相互作用能  用MP2计算时 E-A  E_B   应该是EUMP2 后面的值  还是  scf done 后的HF的能量?如果都不是  那应该是什么?  老师  下面是我做的mp2计算的输出文件  请问 在这种情况下 Eint  是该如何算?

43

帖子

0

威望

1638

eV
积分
1681

Level 5 (御坂)

12#
发表于 Post on 2016-9-28 19:45:03 | 只看该作者 Only view this author
sobereva 发表于 2016-9-28 16:50


一般直接取counterpoise任务最后输出的complexation energy =      xxx kcal/mole (corrected)就 ...

两个片段还好,如果设置八个片段,如何从结果文件中计算两两片段之间的相互作用能呢?是不是很复杂了

6万

帖子

99

威望

5万

eV
积分
120165

管理员

公社社长

13#
发表于 Post on 2016-9-28 22:21:48 | 只看该作者 Only view this author
zycz138 发表于 2016-9-28 19:45
两个片段还好,如果设置八个片段,如何从结果文件中计算两两片段之间的相互作用能呢?是不是很复杂了

不复杂啊。
你要算哪两部分的相互作用能,就把哪两部分在counterpoise时各设成一个片段就完了
北京科音自然科学研究中心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, 2025-8-16 17:40 , Processed in 0.181856 second(s), 22 queries , Gzip On.

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