24小时热门版块排行榜    

查看: 3369  |  回复: 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的回帖

zhangmt

至尊木虫 (著名写手)

我叫MT

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
swfctheochem: 金币+10, ★★★很有帮助, 很有帮助 2013-08-20 21:20:25
搜索"EUMP2",我不信就给你一个能量。
一群自以为正义凛然的年轻人将一切不能以科学解释的事情定性为封建迷信并大刀阔斧地进行消灭,其实这是修养不足学识浅薄的一种体现,也是可恶的偏执和愚蠢的自以
4楼2013-08-13 00:11:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 swfctheochem 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 挑数据是造假吗 +9 电话的建设 2024-12-21 9/450 2024-12-23 08:55 by Xmc_金虫
[考博] 申博的科研计划书怎么写? +6 爱喝风的龙卷茶 2024-12-19 11/550 2024-12-23 07:55 by vincent_hpax
[硕博家园] 纠结退学ing +17 幽皮皮皮 2024-12-22 26/1300 2024-12-23 07:23 by 张开翅膀的鱼
[论文投稿] MDPI投稿连续拒稿是什么原因?文章质量不好吗?还是其他原因? 2+9 多听多看多学 2024-12-17 18/900 2024-12-22 19:20 by jurkat.1640
[考博] 2025申博求助 +7 125814 2024-12-22 7/350 2024-12-22 17:07 by 毕生所学
[找工作] 化学博士 +5 实名心动国民萌 2024-12-17 7/350 2024-12-22 15:33 by 大发财树
[有机交流] 装柱干法上样 +4 刘松垚 2024-12-20 4/200 2024-12-22 14:24 by 88817753
[硕博家园] 脂质氧化 +3 song794 2024-12-22 4/200 2024-12-22 14:12 by song794
[教师之家] 咨询一下,是不是教授之间比较忌讳谈论文造假,因为可能很多人都是造假 +7 akslis2024 2024-12-21 8/400 2024-12-22 12:30 by Quakerbird
[考博] 想找一位年轻的985博导 +6 艾弗森迪尔 2024-12-21 7/350 2024-12-22 09:33 by 飞翔和星辰
[硕博家园] 博一 +10 幽皮皮皮 2024-12-17 16/800 2024-12-22 07:37 by Ermito
[教师之家] 今年的科研绩效只奖励Q1了 +12 holypower 2024-12-20 17/850 2024-12-21 16:13 by willbuilder
[教师之家] 你们都降薪了吗? +10 红枣葡萄干 2024-12-18 10/500 2024-12-20 17:32 by 开心就好56
[考博] 有偿聘请兼职硕士 +10 大发财树 2024-12-19 10/500 2024-12-20 17:30 by mopsite
[论文投稿] Transportation Research Part C 投稿模板 5+3 Yetzirah 2024-12-20 4/200 2024-12-20 16:32 by 北京莱茵润色
[教师之家] 某老师在组会上总是针对我带的学生怎么办? +15 lice_1987 2024-12-17 23/1150 2024-12-20 04:50 by 胖胖的大海
[论文投稿] Molecular Crystals and Liquid Crystals投稿8个月没有信息 +4 书剑如风 2024-12-17 12/600 2024-12-20 04:40 by 胖胖的大海
[论文投稿] 求推荐5分二区以上肿瘤接受生信结合实验的期刊。。。 +3 lice_1987 2024-12-17 6/300 2024-12-19 17:01 by Andy_124
[考博] 有未就业的博士吗 +4 大发财树 2024-12-17 4/200 2024-12-18 14:13 by 她在学习
[论文投稿] 小论文投稿求助 100+3 容城侯 2024-12-16 4/200 2024-12-17 17:07 by 大发财树
信息提示
请填处理意见