24小时热门版块排行榜    

查看: 500  |  回复: 1

gongxd325

至尊木虫 (文坛精英)

[交流] 【讨论】振动分析和ReadIso得到的H和G为何不同? 已有1人参与

对CH4采用优化后直接作振动分析的方法(1)以及ReadIsotopes的办法(2),计算1atm、298K的热力学参数,得到的结果不同,为什么?
(1)
% chk=ch4
# B3LYP/3-21g opt freq

CH4

0  1
6    -2.014320    0.304667   -0.022383
1    -1.657666   -0.704143   -0.022383
1    -1.657647    0.809065    0.851269
1    -1.657647    0.809065   -0.896034
1    -3.084320    0.304680   -0.022383

Sum of electronic and thermal Enthalpies=             -40.252193
Sum of electronic and thermal Free Energies=          -40.273319

(2)
% chk=ch4
# B3LYP/3-21g Freq=(ReadIso,ReadFC) Geom=Check Guess=Read

Repeat at 298K

0 1

298.15 1.0
12
2
2
2
2

Sum of electronic and thermal Enthalpies=             -40.264270
Sum of electronic and thermal Free Energies=          -40.286811

振动分析得到的结果(括弧中的)和不同温度下ReadIso所得结果在一起,发现有问题,变化不连续,是ReadIso计算有问题?
           H           G
200  -40.265628  -40.279649
298  -40.264270  -40.286811
(298  -40.252193  -40.273319)
300  -40.264242  -40.286951
400  -40.262609  -40.294759
500  -40.260684  -40.303013
600  -40.258480  -40.311681
700  -40.256022  -40.320739

另外,采用ReadIsotopes的办法(2),G09的结果都是298.15K的,没有其他温度的。是G09改了输入格式?Manual给的例子和G03一样的呀!

[ Last edited by gongxd325 on 2010-6-8 at 15:47 ]
回复此楼

» 猜你喜欢

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

一个人做点好事并不难难的是一辈子做好事!!!
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

雪狼乖乖

木虫 (正式写手)



小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
458337楼: Originally posted by gongxd325 at 2010-06-08 15:36:56:
对CH4采用优化后直接作振动分析的方法(1)以及ReadIsotopes的办法(2),计算1atm、298K的热力学参数,得到的结果不同,为什么?
(1)
% chk=ch4
# B3LYP/3-21g opt freq

CH4

0  1
6    -2.014320    ...

两年前的帖子.......是不是因为第二个计算的是同位素D的原因造成的
你的输入是 12  2  2  2  2
2楼2012-05-14 02:12:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 gongxd325 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[公派出国] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 CyTWnZiM0cS6 2026-08-02 5/250 2026-08-02 20:52 by 3VGOm6paSgqv
[公派出国] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 CyTWnZiM0cS6 2026-08-02 3/150 2026-08-02 20:44 by 3VGOm6paSgqv
[教师之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 CyTWnZiM0cS6 2026-08-02 5/250 2026-08-02 20:43 by 3VGOm6paSgqv
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 AS6fX1tr1on4 2026-08-02 9/450 2026-08-02 20:32 by 3VGOm6paSgqv
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 AS6fX1tr1on4 2026-08-02 5/250 2026-08-02 20:31 by 3VGOm6paSgqv
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 AS6fX1tr1on4 2026-08-02 7/350 2026-08-02 20:24 by 3VGOm6paSgqv
[基金申请] 这种情况还有戏吗 +5 drbart 2026-07-27 12/600 2026-08-02 20:21 by drbart
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 AS6fX1tr1on4 2026-08-02 6/300 2026-08-02 20:14 by 3VGOm6paSgqv
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 AS6fX1tr1on4 2026-08-02 14/700 2026-08-02 20:11 by 3VGOm6paSgqv
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 AS6fX1tr1on4 2026-08-02 7/350 2026-08-02 20:11 by 3VGOm6paSgqv
[基金申请] 面上提前没消息,有中的吗 +11 archvillain 2026-08-02 13/650 2026-08-02 20:09 by symmetry
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 AS6fX1tr1on4 2026-08-02 13/650 2026-08-02 20:04 by 3VGOm6paSgqv
[考博] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 cBGhfJFuNMca 2026-08-02 21/1050 2026-08-02 19:52 by 3VGOm6paSgqv
[博后之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 AS6fX1tr1on4 2026-08-02 8/400 2026-08-02 18:24 by ca3DVAx4OB19
[基金申请] 准备明年的基金了 +7 Tide man 2026-07-27 7/350 2026-08-02 16:43 by gltch
[基金申请] 2026年国自然面上资助率 +15 布布和一二 2026-07-30 21/1050 2026-08-02 10:56 by iaeyuan
[教师之家] 基础研究怎么拉横向,学校到款任务越来越多,难以完成 拉横向,都有哪些途径啊 +8 锦衣卫寒战 2026-07-28 10/500 2026-08-02 10:00 by 锦衣卫寒战
[考博] 2027年申博 50+3 射雕英雄胜 2026-07-30 3/150 2026-08-02 09:36 by lfy8008
[高分子] HXDI做水性聚氨酯乳液,是不是特别容易出渣 15+3 yuyusuv 2026-07-29 3/150 2026-07-31 09:21 by huizingga
[有机交流] 产物和副产物价值比较 30+4 小汤02 2026-07-27 5/250 2026-07-30 12:44 by czyzsu
信息提示
请填处理意见