24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

cquccc

铜虫 (正式写手)

应助: 2 (幼儿园)
金币: 17
散金: 385
红花: 1
帖子: 597
在线: 294.8小时
虫号: 1085708
注册: 2010-08-31
性别: GG
专业: 理论和计算化学

[求助] 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

回复此楼

» 猜你喜欢

法国里昂大学招生CSC 资助博士生已经有0人回复
法国里昂大学接收CSC 资助访问学者已经有0人回复
物理化学论文润色/翻译怎么收费? 已经有159人回复
如何掌握一门外语？已经有0人回复
企业为什么一定要做“专精特新”认定？如何成为专精特新企业？已经有5人回复
求助ASIQwin新版本的软件已经有0人回复
北京工业大学材料科学与工程学院《院士团队》招收工程专业博士研究生（计算优先）已经有1人回复
期待科研合作，共同发表论文已经有214人回复
如何掌握一门外语？已经有25人回复
吉林大学化学学院理论化学研究所刘靖尧教授课题组招收2026级申请考核制博士研究生已经有2人回复
国家级人才团体课题组招收2026届博士已经有1人回复

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

利用声子谱准确地计算气固反应的焓变及熵变已经有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人回复

1楼 2013-07-11 11:13:53

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

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/ 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
等等来分别获取不同温度/压力/校正因子下的频率/能量校正等等热力学数据。