24小时热门版块排行榜    

查看: 2674  |  回复: 5

hyshun7017

新虫 (小有名气)

[求助] Gaussian计算结果能量的处理

举例说明:

用CCSD(T)//6-311++G**//MP2/6-31+G**计算得到结果如下:

MP2/6-31+G**结构优化得到的单点能为:
MP2 = -188.1179592                                                    Single point energy

MP2/6-31+G** FREQ得到的结果是:
Zero-point correction=                           0.011399 (Hartree/Particle)        ZPE
Thermal correction to Energy=                    0.014052
Thermal correction to Enthalpy=                  0.014997
Thermal correction to Gibbs Free Energy=        -0.009341
Sum of electronic and zero-point Energies=           -188.106560                        E0
Sum of electronic and thermal Energies=              -188.103907                          E
Sum of electronic and thermal Enthalpies=            -188.102963                          H
Sum of electronic and thermal Free Energies=         -188.127300                        G

CCSD(T)//6-311++G**//MP2/6-31+G**计算得到的单点能如下:
CCSD(T) = -188.2207503                                     Single point energy

问题是:
1.    请问在进行论文发表的时候应该用哪个数值呢?
2.     如果有相应的实验结果的话是不是就需要用如下的Gibbs Free Energy进行能量处理呀?
CCSD(T) [-188.2207503]  + Thermal correction to Gibbs Free Energy [-0.009341] =  -188.2300913

3.   如果没有实验结果对照的话是需要用哪个能量进行能垒图作图?

谢谢指教!!不甚感激!!
回复此楼

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

过渡态计算

» 猜你喜欢

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

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

hyshun7017

新虫 (小有名气)

找了一下,好像不可以修改。

我问的问题是针对化学反应过渡态能垒图的计算的。

谢谢!!!
hyshun
2楼2012-03-29 21:20:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

nedved33

金虫 (小有名气)

同求啊,个人觉得应该是内能E
3楼2012-03-29 23:52:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xchem6628

银虫 (初入文坛)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
小红豆: 金币+5, 3x 2012-03-30 12:24:47
hyshun7017: 金币+2, 谢谢 2012-03-30 13:38:15
你好,我们一般画化学反应势能曲线时用的是这个值:
E=CCSD(T) [-188.2207503]  +  Zero-point correction [0.011399],
而且一般能量结果:以Hartree为单位的,保留5位小数,计算相对能量Er,以kJ mol-1为单位的,保留1位小数。
4楼2012-03-30 11:41:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hyshun7017

新虫 (小有名气)

引用回帖:
4楼: Originally posted by xchem6628 at 2012-03-30 11:41:35:
你好,我们一般画化学反应势能曲线时用的是这个值:
E=CCSD(T)   +  Zero-point correction ,
而且一般能量结果:以Hartree为单位的,保留5位小数,计算相对能量Er,以kJ mol-1为单位的,保留1位小数。

E=CCSD(T) [-188.2207503]  +  Zero-point correction [0.011399]
就是通常所说的E0,0 K时分子的内能。

比较困惑的是:
举例来说:以下所列文献全部出之于 K. N. Houk:

1. Org. Lett. 2003, 5, 1249  和  Org. Lett. 2006, 6, 4273  这两篇文献,作者进行过渡态能量的比较时,使用的是 E0 的能量的对比。

2. 在 Org. Lett. 2006, 8, 1513 这篇文献中,对于过渡态能量的对比,作者却是用Enthalpy 进行比较的。为什么此篇文章不用 E0 进行能量的对比?

希望能够得到解答,谢谢!!!
hyshun
5楼2012-03-30 13:35:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xchem6628

银虫 (初入文坛)

引用回帖:
5楼: Originally posted by hyshun7017 at 2012-03-30 13:35:36:
E=CCSD(T)   +  Zero-point correction
就是通常所说的E0,0 K时分子的内能。

比较困惑的是:
举例来说:以下所列文献全部出之于 K. N. Houk:

1. Org. Lett. 2003, 5, 1249  和  Org. Lett. 2006, 6,  ...

你好,这个没有规定必须用E0进行对比,还有用G0,S等来做对比的,具体情况具体分析
6楼2012-04-01 09:41:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 hyshun7017 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 Vi50GxzrFcSG 2026-08-07 6/300 2026-08-08 15:07 by oEVWOejN9taj
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 WQyTGMbfH7Cx 2026-08-07 4/200 2026-08-08 15:02 by oEVWOejN9taj
[基金申请] 这样的filecode谁见过 +5 布布和一二 2026-08-08 8/400 2026-08-08 14:54 by archvillain
[硕博家园] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 HEQlVqMTIA7d 2026-08-07 6/300 2026-08-08 14:47 by oEVWOejN9taj
[博后之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 HEQlVqMTIA7d 2026-08-07 6/300 2026-08-08 14:42 by oEVWOejN9taj
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 HEQlVqMTIA7d 2026-08-07 6/300 2026-08-08 14:39 by oEVWOejN9taj
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 HEQlVqMTIA7d 2026-08-07 5/250 2026-08-08 14:27 by oEVWOejN9taj
[基金申请] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 HEQlVqMTIA7d 2026-08-07 4/200 2026-08-08 14:27 by oEVWOejN9taj
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 HEQlVqMTIA7d 2026-08-07 5/250 2026-08-08 14:22 by oEVWOejN9taj
[论文投稿] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 KXLV3nuBVBY7 2026-08-07 6/300 2026-08-08 14:01 by oEVWOejN9taj
[基金申请] 关于代码变化问题,想知道的进来 +11 且听虎啸 2026-08-07 14/700 2026-08-08 12:50 by 布布和一二
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 Vi50GxzrFcSG 2026-08-07 6/300 2026-08-08 10:19 by 3OOjAIS77qg2
[基金申请] 基金中了 +14 laoda193707 2026-08-06 14/700 2026-08-08 00:23 by 实验小白ha
[基金申请] 化学口download_prp&fileCode的固定段好像这几天一直没变,有变的大神么? +3 Tide man 2026-08-07 4/200 2026-08-07 22:39 by Tide man
[基金申请] filecode与中标关系的预测 +3 布布和一二 2026-08-07 3/150 2026-08-07 15:09 by gltch
[基金申请] filecode变化情况 +6 布布和一二 2026-08-07 22/1100 2026-08-07 14:45 by 且听虎啸
[基金申请] 关于豆爷回答的JTJC与%2F数量 +4 yang182083 2026-08-06 6/300 2026-08-07 12:43 by 虫友是什么虫
[基金申请] 8月时间戳变的,举个手。玩一下,释放压力 +9 archvillain 2026-08-04 11/550 2026-08-05 20:06 by wlwhappy
[考博] 【2027博士申请】纳米药物递送方向 20+3 13586093586 2026-08-03 4/200 2026-08-05 09:59 by lfy8008
[基金申请] 有没有H口的?有收到消息的吗? +3 超级海虾 2026-08-04 3/150 2026-08-04 17:26 by 学教育滴
信息提示
请填处理意见