24小时热门版块排行榜    

Znn3bq.jpeg
查看: 618  |  回复: 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(金币+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的回帖

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的回帖
相关版块跳转 我要订阅楼主 joy412 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 297,工科调剂? +7 河南农业大学-能 2026-04-14 7/350 2026-04-15 13:26 by 西北望—风沙
[考研] 一志愿华中农业071010,320求调剂 +10 困困困困坤坤 2026-04-14 11/550 2026-04-15 12:45 by 西北望—风沙
[考研] 化学070300 求调剂 +23 哈哈哈^_^ 2026-04-12 23/1150 2026-04-14 16:30 by zhouxiaoyu
[考研] 一志愿085502,267分求调剂 +19 再忙也要吃饭啊 2026-04-08 20/1000 2026-04-14 16:03 by zs92450
[考研] 366求调剂 +11 不知名的小卅 2026-04-11 11/550 2026-04-14 15:50 by zs92450
[基金申请] RY:中国产出的科学垃圾论文,绝对数量和比例都世界第一 +6 zju2000 2026-04-14 17/850 2026-04-14 14:34 by jurkat.1640
[考研] 求调剂 +16 雪逢冬 2026-04-10 16/800 2026-04-14 14:27 by 逆水乘风
[考研] 085600材料与化工349分求调剂 +16 李木子啊哈哈 2026-04-12 17/850 2026-04-14 09:11 by fenglj492
[考研] 本科西工大 324求调剂 +5 wysyjs25 2026-04-10 5/250 2026-04-13 23:08 by pies112
[考研] 2本,初试303,0860求调剂 +6 floriea 2026-04-12 8/400 2026-04-12 18:13 by zhouxiaoyu
[考研] 化工调剂求导师收留!一志愿失利,踏实肯干,有植物提取科研经历 +20 yzyzx 2026-04-09 21/1050 2026-04-12 00:12 by 小小小小啦啦啦
[考研] 343求调剂 +9 王国帅 2026-04-10 9/450 2026-04-11 20:31 by dongdian1
[考研] 调剂求助 +6 果然有我 2026-04-11 7/350 2026-04-11 16:22 by 明月此时有
[考研] 085410-273求调剂 +6 X1999 2026-04-10 6/300 2026-04-11 10:32 by Delta2012
[考研] 中药学调剂 初试324 +4 洋甘菊、 2026-04-10 6/300 2026-04-11 09:41 by gong120082
[考研] 本科211 工科085400 280分求调剂 可跨专业 +11 LZH(等待调剂中 2026-04-10 11/550 2026-04-11 08:39 by zhq0425
[考研] 考研调剂 +26 硕星赴 2026-04-09 27/1350 2026-04-10 22:24 by 猪会飞
[考研] 一志愿京区985,085401电子信息,本科电子信息 +3 阳光开朗的男孩 2026-04-10 3/150 2026-04-10 16:29 by sophia_93
[考研] 青岛科技大学材料学院,环境学院调剂补录4月10日以前都可以 +3 1青科大。 2026-04-09 5/250 2026-04-10 09:58 by 翩翩一书生
[考研] 286求调剂 +19 Faune 2026-04-08 20/1000 2026-04-09 08:36 by 哦哦123
信息提示
请填处理意见