| 查看: 862 | 回复: 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有说明，也没有看太明白，恳请高手答疑！ |
回复此楼» 猜你喜欢
夏威夷大学生物胶体与界面实验室诚招全奖博士生和博士后
已经有1人回复
-
200年来人们一直无法破解的难题： 气液临界点状态方程
已经有7人回复
-
物理化学论文润色/翻译怎么收费?
已经有149人回复
-
细胞培养，这22个细节一定要注意！（一）
已经有0人回复
-
关于引力延迟下双星系统角动量不守恒的疑惑---求高手解答！
已经有10人回复
-
求助立方晶体In(OH)3的CIF卡片
已经有2人回复
-
CSC访学博后项目获批，外方学校暂停合作该怎么办？
已经有51人回复
-
[电子教材]天津大学物理化学（第七版）上、下册
已经有1人回复
-
阴离子交换膜电解二氧化碳还原
已经有0人回复
-
分子动力学模拟合作需求
已经有0人回复
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
★ ★ ★
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
6楼2008-10-20 22:19:11