| 查看: 4571 | 回复: 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，请各位指点！非常感谢！ |
回复此楼» 猜你喜欢
中南大学化学化工学院招收电池储能、理论计算方向博士生
已经有16人回复
-
请问哪里可以代烧带压力的炉子？
已经有6人回复
-
物理化学论文润色/翻译怎么收费?
已经有258人回复
-
关于合成MIL-125的求助
已经有0人回复
-
双非本硕光催化方向想读博，请各位老师看看我！
已经有2人回复
-
瑞典林雪平大学范立洲课题组现招聘博士后、岗位制博士生
已经有0人回复
-
中国科大仿生界面材料科学全国重点实验室戴懿涛研究员/江雷院士团队招聘科研助理
已经有0人回复
-
中国科大仿生界面材料科学全国重点实验室戴懿涛研究员/江雷院士团队招聘科研助理
已经有0人回复
-
寻求物理化学学报最新的投稿模板！
已经有4人回复
-
投稿咨询
已经有7人回复
-
自科面上 散金币 祈福
已经有61人回复
» 本主题相关商家推荐: (我也要在这里推广)
11楼2015-12-30 10:48:04
xjyuefan
铁杆木虫 (正式写手)
- 应助: 12 (小学生)
- 金币: 5130.7
- 散金: 100
- 红花: 1
- 帖子: 462
- 在线: 238.7小时
- 虫号: 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
5