计算化学公社

标题: 一键拟合任意级别ZPE校正因子的小脚本fitScFact [打印本页]

作者
Author: Voidmio 时间: 2025-4-18 18:59
标题: 一键拟合任意级别ZPE校正因子的小脚本fitScFact
本帖最后由 Voidmio 于 2025-5-3 20:05 编辑

笔者在对有机体系进行能量计算时，经常在完成后才发现自己使用的级别并没有前人进行各项校正因子的拟合，不堪其忧，自己逐个计算拟合又太过麻烦。
发现 sob 社长在自行拟合基频校正因子与ZPE校正因子的简单方法（http://sobereva.com/391）一文中提到了：

顺带一提，以后笔者有时间时候可能还会写个小程序，使得频率校正因子的拟合可以一键完成，届时用户只需填写计算级别就够了。

遂决定完成之。成品可见附件
以Truhlar JCTC, 6, 2872 (2010) 中 ZPVE15/10 数据集与为例。只需已将g16正常写入环境变量，如脚本目录下运行：

./fitScFact.sh "b3lyp 6-31g*"

复制代码

稍等，即可得到输出内容：
(, 下载次数 Times of downloads: 8)
可见结果与帖内人工拟合结果一致。
同时生成 ZPE.txt 以供查看。

亦可在自己新增分子/新建数据集，只需要在数据集目录下添加与我相同格式的 gjf 文件，并将 fitScFact.sh 中 trainingSet 部分改为相应数据集文件夹名称即可。如我使用部分 Z1 数据集：
(, 下载次数 Times of downloads: 5)
可见结果与常用的 0.9806 肥肠接近。

(, 下载次数 Times of downloads: 30)

目前的问题：由于毕业论文压力暂时搁置了，翻出来为我所用的数据集暂时很少，最常使用的 JPC, 100,16502 (1996) 的完整 Z1 数据集害没找到，也没想好频率校正因子相关怎么写比较合适。欢迎各位多多提提意见与想法

作者
Author: happyspider 时间: 2025-4-27 21:49
确实挺方便的，算了一下B3LYP/6-311G**，Z1数据集是0.987712，和http://bbs.keinsci.com/forum.php?mod=viewthread&tid=3805里0.9888差别不大

作者
Author: Voidmio 时间: 2025-5-3 20:00

happyspider 发表于 2025-4-27 21:49
确实挺方便的，算了一下B3LYP/6-311G**，Z1数据集是0.987712，和http://bbs.keinsci.com/forum.php?mod=vie ...

我那个 Z1 不是完整的 Z1，其实是 GJS，所以会有点差距

之前随便翻了翻，没有找到整合好的 Z1 集数据

作者
Author: Voidmio 时间: 2025-5-3 20:02

happyspider 发表于 2025-4-27 21:49
确实挺方便的，算了一下B3LYP/6-311G**，Z1数据集是0.987712，和http://bbs.keinsci.com/forum.php?mod=vie ...

各个数据集间的差别还挺大，所以个人感觉应该针对具体体系再添例子拟合，但业界没有见到这么干的先例，应该也没什么人care

欢迎光临计算化学公社 (http://bbs.keinsci.com/)