24小时热门版块排行榜    

CyRhmU.jpeg
查看: 1014  |  回复: 8
当前主题已经存档。
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

flyskysky

银虫 (正式写手)

[交流] 【求助】请问如何从频率计算out中得到生成热?

有一个化合物,我要算其生成热,输入文件如下:
%mem=512MB
# opt freq rhf/6-31g

Title Card Required

0 1
C                  0.15906300    0.38436100    0.00004800
N                 -1.08192800   -0.15556300   -0.00005800
O                  1.19623400   -0.24317500   -0.00002500
H                 -1.90151300    0.40881800   -0.00003500
H                  0.13560900    1.46762400    0.00007300
H                 -1.18485100   -1.14826600    0.00027500

优化后得到的out文件中,有这么一段:

Warning -- explicit consideration of   2 degrees of freedom as
           vibrations may cause significant error
Vibrational temperatures:    835.61   892.61  1065.29  1698.60  1699.42
          (Kelvin)           2015.63  2233.12  2616.84  2718.17  4696.26
                             5526.98  5732.59

Zero-point correction=                           0.050243 (Hartree/Particle)
Thermal correction to Energy=                    0.053543
Thermal correction to Enthalpy=                  0.054487
Thermal correction to Gibbs Free Energy=         0.026075
Sum of electronic and zero-point Energies=           -168.804829
Sum of electronic and thermal Energies=              -168.801529
Sum of electronic and thermal Enthalpies=            -168.800585
Sum of electronic and thermal Free Energies=         -168.828997

                     E (Thermal)             CV                S
                      KCal/Mol        Cal/Mol-Kelvin    Cal/Mol-Kelvin
Total                   33.599              9.422             59.798
Electronic               0.000              0.000              0.000
Translational            0.889              2.981             37.339
Rotational               0.889              2.981             21.175
Vibrational             31.821              3.461              1.283
Vibration  1             0.937              1.073              0.484
Vibration  2             0.980              0.989              0.416
                       Q            Log10(Q)             Ln(Q)
Total Bot       0.101457D-11        -11.993716        -27.616553
Total V=0       0.130782D+12         11.116548         25.596797
Vib (Bot)       0.902340D-23        -23.044630        -53.062221
Vib (Bot)  1    0.262171D+00         -0.581415         -1.338757
Vib (Bot)  2    0.235625D+00         -0.627779         -1.445515
Vib (V=0)       0.116315D+01          0.065635          0.151129
Vib (V=0)  1    0.106457D+01          0.027172          0.062566
Vib (V=0)  2    0.105274D+01          0.022320          0.051394
Electronic      0.100000D+01          0.000000          0.000000
Translational   0.118737D+08          7.074585         16.289833
Rotational      0.946953D+04          3.976329          9.155835
***** Axes restored to original set *****
-------------------------------------------------------------------
Center     Atomic                   Forces (Hartrees/Bohr)
Number     Number              X              Y              Z
-------------------------------------------------------------------
    1          6           0.000213175   -0.000103736   -0.000222558
    2          7          -0.000415508   -0.000162012    0.000055919
    3          8           0.000067928    0.000129979    0.000081303
    4          1           0.000079427   -0.000100804    0.000007142
    5          1           0.000025047   -0.000046140    0.000077376
    6          1           0.000029931    0.000282713    0.000000818
-------------------------------------------------------------------
Cartesian Forces:  Max     0.000415508 RMS     0.000156696

这些数据中如何才能求得化合物的生成热呢?
回复此楼
有做顺磁自由基的朋友吗,大家交流一下吧!flyskysky@126.com
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

airytang

银虫 (初入文坛)

★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
lei0736(金币+2,VIP+0):谢谢 7-17 11:35
引用回帖:
Originally posted by flyskysky at 2009-7-17 09:48:
是不是  -168.800585 - 0.054487 - 0.050243 就可以得到我需要的结果了?

如果我比较一个构象 比另一个构象 稳定 多少,我应该用哪一组数据来进行比较啊?

我认为比较两个构象的稳定,用相同方法和基组计算二者的single point energy就可以。单点能越负,构象越稳定。
虫虫
8楼2009-07-17 10:02:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 9 个回答

flyskysky

银虫 (正式写手)

是不是把 Sum of electronic and zero-point Energies ( -168.804829)乘以627.5095 kcal/mol 就可以啦!
有做顺磁自由基的朋友吗,大家交流一下吧!flyskysky@126.com
2楼2009-07-16 23:53:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yjcmwgk

禁虫 (文坛精英)

密度泛函·小卒

优秀版主

★ ★ ★ ★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
lei0736(金币+5,VIP+0):谢谢 7-17 09:29
不是这样的。Sum of electronic and zero-point Energies 是该化合物的energy,而不是生成热。HF和DFT计算中,化合物的energy并无明确意义,只有能量的差值才有意义。很多时候生成热是针对某个反应来说的(这也是生成热的定义,见PS1)。你还要计算你所选定的反应物的energy。

生成热=E(生成物)-E(反应物1)-E(反应物2)-E(反应物3)-……

PS1:生成热的定义:
在一定温度和压力下,由最稳定的单质生成1摩尔纯物质的热效应。现在多称生成焓,因为此生成反应的热效应等于该过程体系焓的增量。标准状态下(各单质与生成的纯物质皆处于标准状态)的生成焓称为标准生成热(standard heat of formation),现行手册上称为摩尔标准生成焓,简称标准生成焓,用符号△fHθm表示
3楼2009-07-17 06:57:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

flyskysky

银虫 (正式写手)

谢谢,非常感谢!
有做顺磁自由基的朋友吗,大家交流一下吧!flyskysky@126.com
4楼2009-07-17 09:20:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复(可上传附件)
信息提示
请填处理意见