24小时热门版块排行榜    

查看: 2079  |  回复: 8
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

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
回复此楼

» 猜你喜欢

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

已阅   回复此楼   关注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的回帖
查看全部 9 个回答

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的回帖
信息提示
请填处理意见