24小时热门版块排行榜    

查看: 3848  |  回复: 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,请各位指点!非常感谢!
回复此楼

» 收录本帖的淘帖专辑推荐

高斯

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xjyuefan

铁杆木虫 (正式写手)

两个体系,那分子坐标是分开写吗?
2楼2013-08-11 11:40:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

swfctheochem

铜虫 (初入文坛)

引用回帖:
2楼: Originally posted by xjyuefan at 2013-08-11 11:40:21
两个体系,那分子坐标是分开写吗?

不用分开写啊,每个原子后面有数字描述该原子是属于哪个片段的
3楼2013-08-12 13:33:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangmt

至尊木虫 (著名写手)

我叫MT

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
swfctheochem: 金币+10, ★★★很有帮助, 很有帮助 2013-08-20 21:20:25
搜索"EUMP2",我不信就给你一个能量。
一群自以为正义凛然的年轻人将一切不能以科学解释的事情定性为封建迷信并大刀阔斧地进行消灭,其实这是修养不足学识浅薄的一种体现,也是可恶的偏执和愚蠢的自以
4楼2013-08-13 00:11:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

nunup5

木虫 (小有名气)

搜“SCF Done”一共五个结果,就是你说的五个体系!
5楼2013-08-20 09:04:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

春雨飞飞

银虫 (小有名气)

引用回帖:
4楼: Originally posted by zhangmt at 2013-08-13 00:11:05
搜索"EUMP2",我不信就给你一个能量。

如果搜 EUMP2 的话,也要搜出五个吗?
知识的世界,没有国界,更没有限制与范围,懂得分享,才是智慧。
6楼2013-08-20 10:17:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

swfctheochem

铜虫 (初入文坛)

引用回帖:
5楼: Originally posted by nunup5 at 2013-08-20 09:04:52
搜“SCF Done”一共五个结果,就是你说的五个体系!

SCF Done是在一起的还是分别计算的,要挨个找吗?我找过,在同一个地方没有5个能量,谢谢指点!
7楼2013-08-20 17:10:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

swfctheochem

铜虫 (初入文坛)

引用回帖:
6楼: Originally posted by 春雨飞飞 at 2013-08-20 10:17:48
如果搜 EUMP2 的话,也要搜出五个吗?...

我找到了,EUMP2每轮计算都有有五个,用搜索功能找最后五个就是,原来MP2方法与Hf和DFT的输出结果不同,SCFdone不在一起,而且EUMP2也不在一起,MP2能量应是EUMP2后的,我我用上面说的两种计算相互作用能的方法验证了一下,是一样的。
8楼2013-08-20 21:17:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

swfctheochem

铜虫 (初入文坛)

引用回帖:
4楼: Originally posted by zhangmt at 2013-08-13 00:11:05
搜索"EUMP2",我不信就给你一个能量。

楼主说的对,但没说清楚啊!5个EUMP2不再一起,要挨个找。不过还是谢谢!
9楼2013-08-20 21:21:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liyan620

禁虫 (初入文坛)

本帖内容被屏蔽

10楼2015-04-21 01:16:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 swfctheochem 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 深圳大学硕士招生(2026秋,传感器方向,仅录取第一志愿) +4 xujiaoszu 2026-03-11 8/400 2026-03-16 09:45 by xujiaoszu
[考研] 中科院材料273求调剂 +3 yzydy 2026-03-15 3/150 2026-03-15 21:15 by ms629
[考研] 求老师收留调剂 +4 jiang姜66 2026-03-14 5/250 2026-03-15 20:11 by Winj1e
[考研] 274求调剂 +4 时间点 2026-03-13 4/200 2026-03-15 15:29 by Rambo13
[基金申请] 面上和青基一样限30页不合理 +5 wowsunflower 2026-03-10 7/350 2026-03-14 17:21 by kingkocxr
[考研] 290求调剂 +4 @将就将就看 2026-03-10 8/400 2026-03-14 14:23 by 千千运气
[考研] 材料080500调剂求收留 +3 一颗meteor 2026-03-13 3/150 2026-03-14 10:54 by peike
[考研] 288求调剂 +14 王晓阳- 2026-03-09 19/950 2026-03-14 02:05 by JourneyLucky
[考研] 304求调剂 +6 Mochaaaa 2026-03-12 7/350 2026-03-13 22:18 by 星空星月
[考研] 求材料调剂 085600英一数二总分302 前三科235 精通机器学习 一志愿哈工大 +4 林yaxin 2026-03-12 4/200 2026-03-13 22:04 by 星空星月
[考研] [0860]321分求调剂,ab区皆可 +4 宝贵热 2026-03-13 4/200 2026-03-13 22:01 by 星空星月
[考研] (081700)化学工程与技术-298分求调剂 +12 11啦啦啦 2026-03-11 35/1750 2026-03-13 21:25 by JourneyLucky
[考研] 0703化学求调剂 +7 绿豆芹菜汤 2026-03-12 7/350 2026-03-13 17:25 by njzyff
[考研] 302求调剂 +6 负心者当诛 2026-03-11 6/300 2026-03-13 16:11 by JourneyLucky
[考研] 328化工专硕求调剂 +4 。,。,。,。i 2026-03-12 4/200 2026-03-13 14:44 by JourneyLucky
[考研] 296求调剂 +3 大口吃饭 身体健 2026-03-13 3/150 2026-03-13 10:31 by 学员8dgXkO
[考研] 290求调剂 +3 ADT 2026-03-13 3/150 2026-03-13 10:19 by peike
[考研] 283求调剂,材料、化工皆可 +8 苏打水7777 2026-03-11 10/500 2026-03-13 09:06 by Linda Hu
[考研] 327分求调剂086 +4 西红柿?小帅 2026-03-09 7/350 2026-03-10 14:47 by ruiyingmiao
[硕博家园] 木虫好像不热闹了,是不是? +4 偏振片 2026-03-10 4/200 2026-03-10 09:51 by longwave
信息提示
请填处理意见