| 查看: 910 | 回复: 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有说明，也没有看太明白，恳请高手答疑！ |
回复此楼» 猜你喜欢
求ENDNOTE破解软件安装包
已经有4人回复
-
研究生招生
已经有0人回复
-
物理化学论文润色/翻译怎么收费?
已经有286人回复
-
河海大学 （211.双一流高校）钙钛矿光电实验室丁勇团队欢迎有志青年！
已经有3人回复
-
黑龙江省原子能研究院化学课题组招收2026级化学专业或化学相关专业硕士生
已经有0人回复
-
山西大学资源与环境工程研究所化工分离/环境功能材料小组招收2026年博士生
已经有3人回复
-
26博士申请
已经有0人回复
-
一志愿南师大0703化学 275求调剂
已经有6人回复
-
A 区过线就可申！重庆双非公办高校 0856调剂，接受跨专业！纳米材料 / 润滑热门方向
已经有0人回复
-
欢迎调剂长江师范学院汤老师课题组的硕士研究生
已经有0人回复
-
中国电机工程学报的投稿系统
已经有0人回复
wxjbuilder
木虫 (知名作家)
- 应助: 1 (幼儿园)
- 贵宾: 0.05
- 金币: 2986.4
- 散金: 1761
- 红花: 2
- 沙发: 1
- 帖子: 5392
- 在线: 273.2小时
- 虫号: 442356
- 注册: 2007-10-27
- 性别: GG
- 专业: 理论和计算化学
3楼2008-10-19 14:12:52
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
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
