| 查看: 919 | 回复: 5 | |||
| 当前主题已经存档。 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[交流]
【求助】freqchk工具的用法
|
|||
|
gaussian 附带一些小工具，比如说freqchk工具，可以根据某一个温度下频率计算生成的chk文件，来快速的得到其他温度下的能量。 我现在有疑问就是，在运行freqchk工具的时候有个参数： Scale factor for frequencies during thermochemistry? [0=>1/1.12] 对应着输出中的： Thermochemistry will use frequencies scaled by 0.8929 这个参数可以默认此值，也可以不需要换算，两种情况，对于结果输出的能量值是不同的。 这个换算是什么意思？是不是根据你得到chk文件的那个输入文件中调用的方法和基组相对应？计算能量一般是用到很高级别的方法，我用的是组合算法，比如说G3MP2那么我的这个因子又应该怎么设置呢？ 手册P206有说明，也没有看太明白，恳请高手答疑！ |
回复此楼» 猜你喜欢
杰青/长江团队招收调剂
已经有0人回复
-
中南林业科技大学罗勇锋教授团队2026年招收材料类调剂研究生（第二批）
已经有3人回复
-
物理化学论文润色/翻译怎么收费?
已经有79人回复
-
投稿期刊求助
已经有1人回复
-
电子科技大学材料学院王翀教授拟招收工程博士一名（化学、材料方向）
已经有3人回复
-
湖南农业大学能源材料创新团队招收2026化学专业研究生（学硕）调剂
已经有8人回复
-
江西师大化学工程学院调剂今晚0点开放！ 缺额100个，速冲！
已经有0人回复
-
奥尔堡Xinxin Xiao课题组招收一名岗位制博士（酶基生物电化学）
已经有0人回复
-
招调剂
已经有0人回复
-
26年电池方向博士申请
已经有1人回复
-
运动康复求调剂
已经有0人回复
★ ★ ★
zzgyb(金币+3,VIP+0):谢谢你的参与，欢迎再次光临计算模拟版！
zzgyb(金币+3,VIP+0):谢谢你的参与，欢迎再次光临计算模拟版！
|
对，这个换算是根据你得到chk文件的那个输入文件中调用的方法和基组相对应的。比如，如果你用的是HF/6-31G(d)，那么scal=0.8929。如果你把scal设为了1，那么freqchk输出的就是你原始chk文件中未经scaled的Zero-point energy和Thermal correction to energy。 G2等复合方法包含多步能量计算，但主要是对电子能量的计算，频率计算（也就是能量热修正的计算）只在其中的一步进行，且那一步中输出的Zero-point energy和Thermal correction to energy都是未经scaled，到最后一步计算时才对它们进行scale并加到总能量中去，进行scale的scale factor也是根据频率计算的那一步方法确定的。 所以freqchk读G2等复合方法的chk文件时，如果把scal设为了1，输出的也是复合方法中频率计算那一步中输出的未经scaled的Zero-point energy和Thermal correction to energy，你必须把scal设为和复合方法中频率计算所用方法相应的scale factor才能得到正确结果。 举个例子说，如果你有了G2方法的298K的chk文件,但想知道443K的能量，那你就把298K的chk文件作为freqchk的输入文件，freqchk的输入参数中温度则设为443K，scale factor设为0.8929（因为G2中算频率的是HF/6-31G(d)）. |
5楼2008-10-20 22:06:10
wxjbuilder
木虫 (知名作家)
- 应助: 1 (幼儿园)
- 贵宾: 0.05
- 金币: 2986.4
- 散金: 1761
- 红花: 2
- 沙发: 1
- 帖子: 5392
- 在线: 273.2小时
- 虫号: 442356
- 注册: 2007-10-27
- 性别: GG
- 专业: 理论和计算化学
★ ★
zzgyb(金币+2,VIP+0):谢谢你的参与，欢迎再次光临计算模拟版！
zzgyb(金币+2,VIP+0):谢谢你的参与，欢迎再次光临计算模拟版！
|
scale factor 是为了消除系统误差 一般计算考虑的是简谐正则振动，实际分子振动还要考虑非简谐正则振动，所以要乘上频率校正因子。 一般来说，不同的方法和基组有不同的校正因子，http://www.bsc.ustc.edu.cn/~dxl/gaussian/VFS.mht 在《Exploring Chemistry with Electronic Structures》频率计算那一章里P64，对于G2,zpe校正是0.8929，G3没有列出。这个得要查查文献看别人如何校正的了 |
2楼2008-10-19 14:03:02
wxjbuilder
木虫 (知名作家)
- 应助: 1 (幼儿园)
- 贵宾: 0.05
- 金币: 2986.4
- 散金: 1761
- 红花: 2
- 沙发: 1
- 帖子: 5392
- 在线: 273.2小时
- 虫号: 442356
- 注册: 2007-10-27
- 性别: GG
- 专业: 理论和计算化学
3楼2008-10-19 14:12:52
4楼2008-10-19 14:47:05
