| 查看: 2041 | 回复: 8 | ||
[求助]
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 |
回复此楼» 猜你喜欢
招贤纳博（已结束）
已经有29人回复
-
表征测试机构后台私聊我
已经有0人回复
-
物理化学论文润色/翻译怎么收费?
已经有299人回复
-
需要一个圆偏振光诱导化学反应，请问使用圆偏振光作为光源，应该需要哪些设备搭建
已经有0人回复
-
Coordination Chemistry Reviews
已经有1人回复
-
请问现在还有电池材料方向的博导招人吗
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
利用声子谱准确地计算气固反应的焓变及熵变
已经有4人回复- 【请教】 2D Gaussian Fitting based on Matlab
已经有4人回复
- 小白问题：吉布斯自由能和活化能
已经有6人回复
- 求gaussian计算反应速率常数，指前因子方法。
已经有25人回复
- Gaussian计算和其他方法得到的热力学数据哪个更准确呢
已经有11人回复
- 化学键的力常数如何获得
已经有17人回复
- 【求助】在gaussian过渡态优化计算出现的问题~请教高手~
已经有3人回复
- 【求助】绝对熵值的计算
已经有8人回复
- 【求助】Gaussian计算完后能量看哪个值啊？
已经有23人回复
- 【求助】Linux下gaussian09安装问题请教
已经有6人回复
- 【求助】Gaussview可以直接导入Chem3D保存的Gauss/Gif文件吗？
已经有7人回复
- 【活动】各种量子力学计算软件使用经验谈（回复就有奖励，讲述详细的奖励5-10金币）
已经有30人回复
- 【求助】前线轨道研究反应选择性的问题
已经有24人回复
2楼2013-07-11 12:55:53
3楼2013-07-11 16:09:09
agent99
禁虫 (正式写手)
物化程序猿
- 应助: 75 (初中生)
- 贵宾: 0.001
- 金币: 3630.7
- 散金: 70
- 红花: 8
- 沙发: 1
- 帖子: 392
- 在线: 177.5小时
- 虫号: 1247746
- 注册: 2011-03-28
- 性别: GG
- 专业: 理论和计算化学
4楼2013-07-12 00:46:16
agent99
禁虫 (正式写手)
物化程序猿
- 应助: 75 (初中生)
- 贵宾: 0.001
- 金币: 3630.7
- 散金: 70
- 红花: 8
- 沙发: 1
- 帖子: 392
- 在线: 177.5小时
- 虫号: 1247746
- 注册: 2011-03-28
- 性别: GG
- 专业: 理论和计算化学
5楼2013-07-12 00:47:09
zhangmt
至尊木虫 (著名写手)
我叫MT
- QC强帖: 5
- 应助: 99 (初中生)
- 金币: 6961.8
- 散金: 10406
- 红花: 49
- 帖子: 1761
- 在线: 763.2小时
- 虫号: 880392
- 注册: 2009-10-22
- 性别: GG
- 专业: 理论和计算化学
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
gmy1990: 金币+3 2013-07-12 01:53:39
cquccc: 金币+20, 很全面很深入，我去试试。谢谢了先！ 2013-07-12 08:28:04
感谢参与，应助指数 +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 等等来分别获取不同温度/压力/校正因子下的频率/能量校正等等热力学数据。 |
6楼2013-07-12 01:52:18
7楼2013-07-12 08:19:40
|
由于输出文件太占版面，也不知道哪里上传，将相关部分放置在下面： 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
9楼2013-07-12 08:49:47