24小时热门版块排行榜    

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

joy412

银虫 (小有名气)

[交流] 【求助】求助:利用高斯计算热力学参数

想要计算一下400K时的水杨醛的吉布斯自由能等热力学参数,今天算了一个下午也没搞定,将我的输入文件贴出来,哪位大虾帮忙解释一下是什么地方出错,说明书也看了,就是算不出来,可能是本人比较笨吧
下面是输入的文件
%chk=B11.chk
%mem=6MW
%nproc=1
# freq  b3lyp/6-31g Freq=ReadIsotopes Test

[No Title]

0 1
C
C                  1              B1
C                  2              B2    1              A1
C                  3              B3    2              A2    1              D1
C                  4              B4    3              A3    2              D2
C                  1              B5    2              A4    3              D3
O                  4              B6    3              A5    2              D4
C                  5              B7    4              A6    3              D5
O                  8              B8    5              A7    4              D6
H                  1              B9    6              A8    5              D7
H                  2             B10    1              A9    6              D8
H                  3             B11    2             A10    1              D9
H                  6             B12    1             A11    2             D10
H                  7             B13    4             A12    3             D11
H                  8             B14    5             A13    4             D12

   B1             1.40759949
   B2             1.39208206
   B3             1.40109382
   B4             1.41509317
   B5             1.38808326
   B6             1.38704085
   B7             1.46383284
   B8             1.24411085
   B9             1.08414041
   B10            1.08515672
   B11            1.08318523
   B12            1.08462443
   B13            0.97530533
   B14            1.10863067
   A1           120.65339523
   A2           119.96571234
   A3           120.30730542
   A4           119.33066540
   A5           115.12176998
   A6           122.65596734
   A7           124.70466196
   A8           120.40657685
   A9           119.96134516
   A10          121.84706946
   A11          121.47734423
   A12          113.60332309
   A13          117.10097231
   D1             0.00000000
   D2             0.00000000
   D3             0.00000000
   D4           179.99619537
   D5          -179.99238557
   D6           179.97812872
   D7           179.99317989
   D8          -179.99572647
   D9          -179.99538278
   D10          179.99389045
   D11          179.97449330
   D12           -0.01651788

1 2 1.5 6 1.5 10 1.0
2 3 1.5 11 1.0
3 4 1.5 12 1.0
4 5 1.5 7 1.0
5 6 1.5 8 1.0
6 13 1.0
7 14 1.0
8 9 2.0 15 1.0
9
10
11
12
13
14
15

400.0 1.0
16
16
13
13
13
13
13
13
13
1
1
1
1
1
1
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ZDBWHZ

金虫 (正式写手)

Zero-point correction=                           0.115744 (Hartree/Particle)
Thermal correction to Energy=                    0.127196
Thermal correction to Enthalpy=                  0.128463
Thermal correction to Gibbs Free Energy=         0.070139
Sum of electronic and zero-point Energies=           -420.553934
Sum of electronic and thermal Energies=              -420.542482
Sum of electronic and thermal Enthalpies=            -420.541215
Sum of electronic and thermal Free Energies=         -420.599539

                     E (Thermal)             CV                S
                      KCal/Mol        Cal/Mol-Kelvin    Cal/Mol-Kelvin
Total                   79.817             33.233             91.497
Electronic               0.000              0.000              0.000
Translational            1.192              2.981             41.772
Rotational               1.192              2.981             29.608
Vibrational             77.432             27.271             20.118
Vibration  1             0.806              1.959              3.769
Vibration  2             0.821              1.922              2.923
Vibration  3             0.838              1.882              2.457
Vibration  4             0.897              1.745              1.657
Vibration  5             0.945              1.639              1.319
Vibration  6             1.010              1.504              1.023
Vibration  7             1.016              1.493              1.003
Vibration  8             1.052              1.424              0.887
Vibration  9             1.132              1.277              0.689
Vibration 10             1.190              1.180              0.583
Vibration 11             1.223              1.127              0.532
这只是初步的结果,后续的工作相信阁下应该自己会处理。如还有问题请联系QQ372029567
4楼2009-04-02 17:38:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 5 个回答

ZDBWHZ

金虫 (正式写手)

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
joy412(金币+20,VIP+0):谢谢,我试试 4-2 19:36
yuhuobuku(金币+2):欢迎参加讨论
先不看这么长的输入文件,敢问再算之前有没有进行结构的优化。
2楼2009-04-02 17:13:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ZDBWHZ

金虫 (正式写手)

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
joy412(金币+29,VIP+0):谢谢帮忙算出来,十分感谢 4-5 10:19
在计算频率前,先进性结构优化,然后在计算频率,在关键字栏里直接输入temperature=400 pressure=1
3楼2009-04-02 17:27:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yxb727

铜虫 (小有名气)

★ ★ ★
joy412(金币+1,VIP+0):谢谢,有其他数值的 4-2 19:36
yuhuobuku(金币+2):欢迎参加讨论
其实只算400K是的吉布斯自由能,完全可以不用高斯,有公式哦,楼主肯定还要其他数据吧
5楼2009-04-02 18:04:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见