24小时热门版块排行榜    

查看: 1603  |  回复: 5
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

suosuosky

铁杆木虫 (著名写手)

[求助] 求助:请问MP2计算团簇到底哪个才是体系的能量?

我计算的是一个二元团簇,输入文件是:
...
#p opt freq rmp2/6-311++g(d,p) counterpoise=2
...

输出文件中:
1. 中间有一个能量数值是:
...
Zero-point correction=                           0.137249 (Hartree/Particle)
Thermal correction to Energy=                    0.148166
Thermal correction to Enthalpy=                  0.149111
Thermal correction to Gibbs Free Energy=         0.095875
Sum of electronic and zero-point Energies=           -951.897947
Sum of electronic and thermal Energies=              -951.887029
Sum of electronic and thermal Enthalpies=            -951.886085
Sum of electronic and thermal Free Energies=         -951.939321
...

2. 最后有个能量数值是:
....
\HF=-526.8069254\MP2=-526.9546706\RMSD=3.217e-15\RMSF=1.433e-05\ZeroPo
......


请问上面的1和2给出的能量为什么相差这么多呢,哪个才是团簇体系真正的能量?
其中的2是否经过零点能和BSSE校正了呢?


谢谢大家!


[ Last edited by suosuosky on 2013-10-30 at 18:44 ]
回复此楼
follow_my_dream
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

suosuosky

铁杆木虫 (著名写手)

引用回帖:
5楼: Originally posted by wolfzhong at 2013-10-31 08:55:58
Sum of electronic and zero-point Energies是MP2+zpe能量
未BSSE校正的能量

bsse校正在后面吧,没有连用过opt 和counterpoise。

谢谢。
我刚才算了一下,"Sum of electronic and zero-point Energies“好像是已经包含BSSE校正的。
因为我用"Counterpoise: corrected energy"的能量数值加上"Zero-point correction"的能量数值就等于"Sum of electronic and zero-point Energies“的数值。
follow_my_dream
6楼2013-10-31 14:29:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 6 个回答

fubest

铁杆木虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与,应助指数 +1
gmy1990: 金币+2 2013-10-30 23:42:09
suosuosky: 金币+3, ★★★很有帮助, 谢谢!那你计算过团簇结合能了吗?是用的是这个能量减去两个单体的优化后的能量吗? 2013-10-31 08:41:40
我也是算的团簇用的是Sum of electronic and zero-point Energies=           -951.897947这个值
2楼2013-10-30 21:53:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wolfzhong

木虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
gmy1990: 金币+3 2013-10-30 23:42:03
suosuosky: 金币+5, ★★★很有帮助, 谢谢! 2013-10-31 08:42:18
Sum of electronic and zero-point Energies=           -951.897947
这是实际体系的能量。

因为你用了关键词couterpoise=2
是做BSSE校正的。对应给出片段1和片段2的能量(未校正的),和校正后片段1和2的能量【里面顺序不太记得了】,以及最后计算的BSSE能量。

MP2=-526.9546706这个能量是其中某个片段的MP2能量,应该是不包括zpe corrected energy(片段没算频率,没zpe)。
3楼2013-10-30 23:23:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

suosuosky

铁杆木虫 (著名写手)

引用回帖:
3楼: Originally posted by wolfzhong at 2013-10-30 23:23:23
Sum of electronic and zero-point Energies=           -951.897947
这是实际体系的能量。

因为你用了关键词couterpoise=2
是做BSSE校正的。对应给出片段1和片段2的能量(未校正的),和校正后片段1和2的能量 ...

谢谢!那请问:Sum of electronic and zero-point Energies这个能量是MP2能量吗,是已经经过ZPL和BSSE校正后的能量吗?
如果计算团簇结合能,是用这个能量减去优化后的两个单体的能量吗?
follow_my_dream
4楼2013-10-31 08:44:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 hKe7U1kad5Sg 2026-10-31 7/350 2026-11-01 15:04 by BCdvBAEpt3zm
[基金申请] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 xvJS1BYDxvBS 2026-10-31 9/450 2026-11-01 14:55 by BCdvBAEpt3zm
[考研] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 xvJS1BYDxvBS 2026-10-31 10/500 2026-11-01 14:50 by BCdvBAEpt3zm
[论文投稿] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 xvJS1BYDxvBS 2026-10-31 7/350 2026-11-01 14:50 by BCdvBAEpt3zm
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +5 1uaZt2kAZj2S 2026-10-31 6/300 2026-11-01 14:42 by BCdvBAEpt3zm
[硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 1uaZt2kAZj2S 2026-10-30 13/650 2026-11-01 14:29 by BCdvBAEpt3zm
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-30 14/700 2026-11-01 14:24 by BCdvBAEpt3zm
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 xvJS1BYDxvBS 2026-10-31 3/150 2026-11-01 12:53 by BCdvBAEpt3zm
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 xvJS1BYDxvBS 2026-10-31 5/250 2026-11-01 12:51 by BCdvBAEpt3zm
[考博] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 1uaZt2kAZj2S 2026-10-31 6/300 2026-11-01 12:40 by BCdvBAEpt3zm
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 1uaZt2kAZj2S 2026-10-31 7/350 2026-11-01 12:40 by BCdvBAEpt3zm
[考研] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 1uaZt2kAZj2S 2026-10-30 16/800 2026-11-01 12:14 by BCdvBAEpt3zm
[硕博家园] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 1uaZt2kAZj2S 2026-10-30 13/650 2026-11-01 12:14 by BCdvBAEpt3zm
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 X5bHxOULxFqb 2026-10-30 18/900 2026-11-01 12:11 by BCdvBAEpt3zm
[教师之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 X5bHxOULxFqb 2026-10-30 15/750 2026-11-01 12:11 by BCdvBAEpt3zm
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 hKe7U1kad5Sg 2026-10-31 5/250 2026-11-01 11:14 by BCdvBAEpt3zm
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 xvJS1BYDxvBS 2026-10-31 9/450 2026-11-01 10:53 by BCdvBAEpt3zm
[论文投稿] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +5 X5bHxOULxFqb 2026-10-30 26/1300 2026-11-01 10:22 by BCdvBAEpt3zm
[博后之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-30 15/750 2026-11-01 05:22 by YXM1FDCSE8Bn
[文学芳草园] 闲聊 +3 myrtle 2026-10-25 6/300 2026-10-29 15:28 by myrtle
信息提示
请填处理意见