| 查看: 3368 | 回复: 10 | |||
swfctheochem铜虫 (初入文坛)
|
[求助]
关于gaussian09中BSSE校正
|
|
在论坛上看到以下内容： 要计算A、B两个分子的相互作用能，在Gaussian03中使用counterpoise=2关键字（可简写为counter=2），会计算5个体系，输出的能量按照如下顺序： E_AB：A、B基组下AB复合物的能量 E_A,bAB：A、B基组下A的能量 E_B,bAB：A、B基组下B的能量 E_A：A基组下A的能量 E_B：B基组下B的能量 最后会输出corrected energy (E_corrected)和BSSE energy (E_BSSE)。E_BSSE是校正能，E_BSSE = (E_A - E_A,bAB) + (E_B - E_B,bAB) ；E_corrected就是消除了因单体基组重叠造成的能量降低后的AB复合物能量，E_corrected = E_AB + E_BSSE。 BSSE校正后的真实的相互作用能可以这样计算：E_interaction = E_corrected - (E_A + E_B)。也可以这样计算，是等价的：E_interaction = E_AB - E_A,bAB - E_B,bAB。 我用gaussian09计算两个分子的相互作用能，输入格式如下 # MP2/6-31+G** Counterpoise=2 Opt Counterpoise with Cartesian 0,1 分子描述 发现输出结果中只有如下两项： Counterpoise: corrected energy = -734.521269484568 Counterpoise: BSSE energy = 0.003653491328 这两项应该就是上面说的corrected energy (E_corrected)和BSSE energy (E_BSSE)。那么前面五项在哪里啊？没有E_A和E_B也算不了最终的相互作用能啊！难道我还要单独计算E_A和E_B吗？按上面的说法应该一次计算就可得到结果啊！哪里出错了？是03和09输出格式的区别吗？另外out文件的最后给出MP2=-385.2025592，这似乎是其中一个分子的总能量，那么另一个分子的能量在哪里？直接用Counterpoise=2 Opt，并没有做频率分析，如果单独算A和B的能量，需要做频率分析考虑零点能吗？第一次计算BSSE，请各位指点！非常感谢！ |
回复此楼» 收录本帖的淘帖专辑推荐
 高斯 |
» 猜你喜欢
纽约州立大学水牛城分校张培鸿教授招聘博士后研究员
已经有0人回复
-
美国西北大学 生物医学实验室 招聘优秀博士后！
已经有3人回复
-
物理化学论文润色/翻译怎么收费?
已经有243人回复
-
北京化工大学普通博士招生---生物方向
已经有2人回复
-
细胞培养遇到黑点怎么办？这里有办法！
已经有0人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 电荷与自旋多重度的问题
已经有5人回复
- Gaussian计算的苯紫外光谱和教材上的差别太大，求助
已经有10人回复
- 关于ZPE和BSSE校正，转载一个个人感觉很好的帖子。
已经有29人回复
- 在Gaussian09计算中考虑标量相对论和旋轨耦合
已经有9人回复
- 求助！gauss03中BSSE校正的问题
已经有19人回复
- gaussian09 linux安装问题
已经有5人回复
- 有关opt+freq+bsse校正输出文件的意思
已经有3人回复
- BSSE校正中 片段2的变化如何改写
已经有8人回复
- 请问在Gaussian计算中怎样考虑弱相互作用？
已经有23人回复
- #p ccsd(t)/cc-pwCVQZ iop(5/13=1) opt freq 提示输入文件有错 用的是GS09
已经有8人回复
- 请教关于零点能校正的问题
已经有5人回复
- 这句话怎么个理解？这个能量怎么计算？
已经有6人回复
- Gaussian 09 中溶剂化的荧光发射
已经有11人回复
- gaussian关键词scan使用
已经有11人回复
- Linux Gaussian09 无法安装。。求高人指点
已经有10人回复
- 【gchcaj个人文集】浅谈Orca计算BSSE校正
已经有5人回复
- 【lihb个人文集】浅谈BSSE矫正之相互作用能的计算
已经有187人回复
- 【sobereva个人文集】乱弹DFT-D（从sobereva的博客中转载过来的）
已经有27人回复
- 【讨论】请问现在有自己写量化代码的虫子吗？
已经有4人回复
- 【求助】siesta计算中，如何进行BSSE修正？
已经有9人回复
- 【求助】BSSE计算出错
已经有19人回复
- 【停止更新】GaussianFAQ实践版（V27.01）【我将在近期对这个帖子进行大修】
已经有568人回复
xjyuefan
铁杆木虫 (正式写手)
- 应助: 12 (小学生)
- 金币: 5125.7
- 散金: 100
- 红花: 1
- 帖子: 461
- 在线: 238.4小时
- 虫号: 338536
- 注册: 2007-04-04
- 性别: GG
- 专业: 配位化学
2楼2013-08-11 11:40:21
swfctheochem
铜虫 (初入文坛)
- 应助: 2 (幼儿园)
- 金币: 13.8
- 散金: 10
- 帖子: 42
- 在线: 43.3小时
- 虫号: 1859455
- 注册: 2012-06-14
- 专业: 理论和计算化学
3楼2013-08-12 13:33:57
zhangmt
至尊木虫 (著名写手)
我叫MT
- QC强帖: 5
- 应助: 99 (初中生)
- 金币: 6961.8
- 散金: 10406
- 红花: 49
- 帖子: 1761
- 在线: 763.2小时
- 虫号: 880392
- 注册: 2009-10-22
- 性别: GG
- 专业: 理论和计算化学
4楼2013-08-13 00:11:05
nunup5
木虫 (小有名气)
- 应助: 4 (幼儿园)
- 金币: 5617.1
- 散金: 178
- 红花: 2
- 帖子: 239
- 在线: 1341小时
- 虫号: 1091965
- 注册: 2010-09-07
- 性别: MM
- 专业: 理论和计算化学
5楼2013-08-20 09:04:52
6楼2013-08-20 10:17:48
swfctheochem
铜虫 (初入文坛)
- 应助: 2 (幼儿园)
- 金币: 13.8
- 散金: 10
- 帖子: 42
- 在线: 43.3小时
- 虫号: 1859455
- 注册: 2012-06-14
- 专业: 理论和计算化学
7楼2013-08-20 17:10:05
swfctheochem
铜虫 (初入文坛)
- 应助: 2 (幼儿园)
- 金币: 13.8
- 散金: 10
- 帖子: 42
- 在线: 43.3小时
- 虫号: 1859455
- 注册: 2012-06-14
- 专业: 理论和计算化学
8楼2013-08-20 21:17:23
swfctheochem
铜虫 (初入文坛)
- 应助: 2 (幼儿园)
- 金币: 13.8
- 散金: 10
- 帖子: 42
- 在线: 43.3小时
- 虫号: 1859455
- 注册: 2012-06-14
- 专业: 理论和计算化学
9楼2013-08-20 21:21:38
|
本帖内容被屏蔽 |
10楼2015-04-21 01:16:58
50
