应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » Gaussian 热化学计算 问题请教
北京石油化工学院2026年研究生招生接收调剂公告
91/1返回列表
查看: 2205  |  回复: 8
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

cquccc

铜虫 (正式写手)

[求助] Gaussian 热化学计算 问题请教

分别在300K和1000K计算CO2，得到的输出文件按的 - Thermochemistry -
部分的自由能、熵、热容等都完全一致！？为什么不同温度下的这些参数一样呢？不同温度下，至少熵一定是不一样的啊。请教各位了。

输入文件如下：
第一个计算 300K
# opt  b3lyp/6-31g* freq=readisotopes

CO2-300K

0 1
C                  
O 1 1.3
O 1 1.3 2 180.0

300 1.0  0.99
12
16
16

第二个计算 1000K
# opt  b3lyp/6-31g* freq=readisotopes

CO2-1000K

0 1
C                  
O 1 1.3
O 1 1.3 2 180.0

1000 1.0  0.99
12
16
16
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2013-07-11 11:13:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yanrding

金虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
cquccc: 金币+5 2013-07-11 16:00:28
gmy1990: 金币+2 2013-07-12 01:53:46
你先看看输出的log文件里面，到底用了哪个温度，看是不是读入温度的部分错了。

另外，你这里的温度给的都是整数格式，最好写为 300.0 1000.0这样的小数格式。
赞一下 回复此楼
华北电力大学（北京）数理学院，计算化学，凝聚态物理。欢迎合作！招硕士/博士研究生，招老师。
2楼2013-07-11 12:55:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cquccc

铜虫 (正式写手)
把温度改为300.0 1000.0之后仍然一样，不是这个原因。
赞一下 回复此楼
3楼2013-07-11 16:09:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

agent99

禁虫 (正式写手)

物化程序猿

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与，应助指数 +1
cquccc: 金币+5 2013-07-12 08:19:09
我好奇你的优化是怎么进行的？有一个180度的角
赞一下 回复此楼
理论与计算化学，化学物理，生物物理
4楼2013-07-12 00:46:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

agent99

禁虫 (正式写手)

物化程序猿
另外，可否贴出输出文件？
赞一下 回复此楼
理论与计算化学，化学物理，生物物理
5楼2013-07-12 00:47:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangmt

至尊木虫 (著名写手)

我叫MT

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
gmy1990: 金币+3 2013-07-12 01:53:39
cquccc: 金币+20, 很全面很深入，我去试试。谢谢了先！ 2013-07-12 08:28:04
1/ readisotopes 设定从G09开始倾向于废止，如果仅仅是要指定温度/压力/频率校正因子，可以使用temperature pressure scale等关键字直接指定
2/ 温度/压力/频率校正因子的设定，对于优化没有任何影响，因此在opt过程中不需要指定，没有意义
3/ 象# opt freq这类写法，是一种简化，程序在执行时事实上是分为opt 和 freq两个步骤来完成的，按照你贴出来的输入文件，那个readisotopes选项在opt时对程序有效（但对于优化没有意义），而在第二步进行freq时对于程序无效，因此你的程序相当于
# opt
.......
--link1--
# freq

也就是频率分析时根本没有起到任何影响，还是采用默认的298.15K 。。。。。。
4/ 采用temperature/pressure/scale等关键字可以避免readisotopes造成上述的错误
5/ 更合适的做法是先做好优化，取优化后的结构加temperature/pressure/scale等关键字来直接做freq计算，或者显式地使用--link1--来指定多步计算，而不要使用opt freq这种简化的方式
或者取优化后的chk文件，使用freqchk程序，直接输入
freqchk youropt.chk n 300 1.0 0.99
freqchk youropt.chk n 1000 1.0 0.99
等等来分别获取不同温度/压力/校正因子下的频率/能量校正等等热力学数据。
赞一下(1人) 回复此楼
一群自以为正义凛然的年轻人将一切不能以科学解释的事情定性为封建迷信并大刀阔斧地进行消灭，其实这是修养不足学识浅薄的一种体现，也是可恶的偏执和愚蠢的自以
6楼2013-07-12 01:52:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cquccc

铜虫 (正式写手)
引用回帖:
4楼: Originally posted by agent99 at 2013-07-12 00:46:16
我好奇你的优化是怎么进行的？有一个180度的角

分子本身就是线性的，180度没有问题。
赞一下 回复此楼
7楼2013-07-12 08:19:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cquccc

铜虫 (正式写手)
引用回帖:
5楼: Originally posted by agent99 at 2013-07-12 00:47:09
另外，可否贴出输出文件？

由于输出文件太占版面，也不知道哪里上传，将相关部分放置在下面：

300K计算的相关输出文件：
- Thermochemistry -
-------------------
Temperature   298.150 Kelvin.  Pressure   1.00000 Atm.
Atom     1 has atomic number  6 and mass  12.00000
.......
This molecule is a prolate symmetric top.
Rotational symmetry number  2.
Rotational temperature (Kelvin)      0.55466
Rotational constant (GHZ):          11.557186
Zero-point vibrational energy      30436.2 (Joules/Mol)
                                    7.27443 (Kcal/Mol)
Vibrational temperatures:    920.93   920.93  1974.07  3505.34
          (Kelvin)

Zero-point correction=                           0.011593 (Hartree/Particle)
Thermal correction to Energy=                    0.014240
Thermal correction to Enthalpy=                  0.015184
Thermal correction to Gibbs Free Energy=        -0.009126
Sum of electronic and zero-point Energies=           -188.569348
Sum of electronic and thermal Energies=              -188.566700
Sum of electronic and thermal Enthalpies=            -188.565756
Sum of electronic and thermal Free Energies=         -188.590066

                     E (Thermal)             CV                S
                      KCal/Mol        Cal/Mol-Kelvin    Cal/Mol-Kelvin Total                    8.936              6.983             51.165
......................................

1000K计算的输出文件
- Thermochemistry -
-------------------
Temperature   298.150 Kelvin.  Pressure   1.00000 Atm.
Atom     1 has atomic number  6 and mass  12.00000
........
Zero-point correction=                           0.011593 (Hartree/Particle)
Thermal correction to Energy=                    0.014240
Thermal correction to Enthalpy=                  0.015184
Thermal correction to Gibbs Free Energy=        -0.009126
Sum of electronic and zero-point Energies=           -188.569348
Sum of electronic and thermal Energies=              -188.566700
Sum of electronic and thermal Enthalpies=            -188.565756
Sum of electronic and thermal Free Energies=         -188.590066

                     E (Thermal)             CV                S
                      KCal/Mol        Cal/Mol-Kelvin    Cal/Mol-Kelvin
Total                    8.936              6.983             51.165
...........................
赞一下 回复此楼
8楼2013-07-12 08:26:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cquccc

铜虫 (正式写手)
实践证明，zhangmt
的说法完全正确。即先优化，之后在freq计算时候设定温度等。
谢谢各位！谢谢zhangmt！
赞一下 回复此楼
9楼2013-07-12 08:49:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 cquccc 的主题更新
91/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 085600材料与化工301分求调剂院校 +16 刺痛jk 2026-04-06 17/850 2026-04-06 19:47 by 电池阿白
[考研] 308求调剂 +13 倘若起风了呢 2026-04-05 13/650 2026-04-06 14:20 by 蒋皓禹
[考研] 070300化学279求调剂 +16 哈哈哈^_^ 2026-03-31 19/950 2026-04-06 14:14 by 无际的草原
[考研] 295求调剂 +14 xndjjj 2026-04-04 15/750 2026-04-06 13:24 by 272061663
[考研] 301求调剂 +7 细胞相关蛋白 2026-04-03 7/350 2026-04-06 11:47 by lijunpoly
[考研] 0855求调剂材料 +11 红桃灼灼 2026-04-04 12/600 2026-04-06 10:26 by 蓝云思雨
[考研] 材料专硕283求调剂 +17 试试看呗 2026-04-04 18/900 2026-04-06 09:24 by 286640313
[考研] 材料专硕322分 +11 哈哈哈吼吼吼哈 2026-04-02 11/550 2026-04-04 23:37 by 永字号
[考研] 材料调剂 +10 懒羊羊轻置玉臀 2026-04-02 11/550 2026-04-04 21:56 by laoshidan
[考研] 材料科学与工程考研 +10 拯救皮特托先生 2026-04-02 10/500 2026-04-03 23:57 by userper
[考研] 总分328生物与医药考数学求调剂 +7 aaadim 2026-04-02 9/450 2026-04-03 22:53 by syh9288
[考研] 366求调剂 +7 sbdnd 2026-04-03 7/350 2026-04-03 12:40 by cymywx
[考研] 326求调剂 +3 9ahye 2026-04-02 4/200 2026-04-03 08:43 by Jaylen.
[考研] 求调剂 +7 Aniyaio 2026-04-02 7/350 2026-04-02 16:42 by zzsw+
[考研] 318求调剂，计算材料方向 +10 吸喵有害笙命 2026-04-01 11/550 2026-04-02 16:29 by oooqiao
[考研] 282求调剂 +13 呼吸都是减肥 2026-04-01 13/650 2026-04-02 14:10 by baoball
[考研] 求调剂，一志愿南京师范大学计算机专硕，初试373，六级通过， +3 计算机追梦人 2026-04-01 3/150 2026-04-02 07:57 by fxue1114
[考研] 生物学327，求调剂 +5 书上的梅子 2026-04-01 6/300 2026-04-02 06:47 by ilovexiaobin
[考研] 求调剂：一志愿：南京大学 专业：0705 总分320 ，本科985，四六级已过 +3 lfy760306 2026-03-31 3/150 2026-04-01 01:57 by Creta
[考研] 070300化学专业279调剂 +10 哈哈哈^_^ 2026-03-31 10/500 2026-03-31 23:13 by liu823948201
信息提示
关闭
请填处理意见
关闭