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

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » Gaussian 热化学计算 问题请教
91/1返回列表
查看: 2043  |  回复: 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返回列表
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭