24小时热门版块排行榜    

Znn3bq.jpeg
查看: 620  |  回复: 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

金虫 (正式写手)

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

金虫 (正式写手)

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的回帖

yxb727

铜虫 (小有名气)

★ ★ ★
joy412(金币+1,VIP+0):谢谢,有其他数值的 4-2 19:36
yuhuobuku(金币+2):欢迎参加讨论
其实只算400K是的吉布斯自由能,完全可以不用高斯,有公式哦,楼主肯定还要其他数据吧
5楼2009-04-02 18:04:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 22408 312求调剂 +18 门路摸摸 2026-04-14 19/950 2026-04-15 14:23 by 巴塞罗那2015
[考研] 求助调剂,跨调 +18 X十甫寸Y 2026-04-11 19/950 2026-04-14 19:26 by Art1977
[考研] 化工学硕294分,求导师收留 +32 yzyzx 2026-04-12 36/1800 2026-04-14 17:45 by lhj2009
[考研] 考研调剂 +13 长弓傲 2026-04-13 14/700 2026-04-14 14:44 by zs92450
[考研] 食品与营养(0955)271求调剂 +15 升格阿达 2026-04-12 16/800 2026-04-14 13:18 by 浮若_安生
[考研] 0856专硕求调剂 希望是a区院校 +24 好好休息好不好 2026-04-09 27/1350 2026-04-13 22:22 by pies112
[考研] 344 材料专业 求调剂211 无地域要求 +8 hualkop 2026-04-11 8/400 2026-04-12 22:24 by fqwang
[硕博家园] 新一代电子信息294求调剂 不挑学校 +7 Ytyt11 2026-04-09 8/400 2026-04-12 16:57 by ajpv风雷
[考研] 0854调剂 +12 长弓傲 2026-04-09 13/650 2026-04-12 09:56 by 逆水乘风
[考研] 一志愿西北工业大学289 085602 +33 yang婷 2026-04-10 34/1700 2026-04-12 08:11 by Art1977
[考研] 一志愿985机械学硕380求调剂 +5 关关雎鸠10 2026-04-11 5/250 2026-04-11 10:10 by 知念。A
[考研] 考研调剂 +26 硕星赴 2026-04-09 27/1350 2026-04-10 22:24 by 猪会飞
[考研] 263能源动力专硕求调剂 +3 加大号饭盒袋 2026-04-10 3/150 2026-04-10 22:23 by 286640313
[考研] 22408 366分,本科211,一志愿西工大 +4 Rubt 2026-04-09 4/200 2026-04-10 19:51 by chemisry
[考研] 265求调剂 +12 风说她早忘了 2026-04-10 13/650 2026-04-10 18:56 by chemisry
[考研] 求调剂 材料与工程 324分 专硕 +19 翩翩一书生 2026-04-10 21/1050 2026-04-10 11:41 by wp06
[考研] 一志愿中科大070300化学,314分求调剂 +12 wakeluofu 2026-04-09 12/600 2026-04-10 09:57 by liuhuiying09
[考研] 085601初试330分找调剂 +10 流心奶黄包l 2026-04-09 10/500 2026-04-10 08:14 by Sammy2
[考研] 材料专硕(0856) 339分求调剂 +9 哈哈哈鹅哈哈哈 2026-04-09 10/500 2026-04-09 20:01 by Orcid
[考研] 086004 求调剂 309 +7 Yin DY 2026-04-08 7/350 2026-04-09 13:59 by Delta2012
信息提示
请填处理意见