|
|
本帖最后由 Voidmio 于 2025-9-26 16:17 编辑
笔者在对有机体系进行能量计算时，经常在完成后才发现自己使用的级别并没有前人进行各项校正因子的拟合，不堪其忧，自己逐个计算拟合又太过麻烦。
发现 sob 社长在自行拟合基频校正因子与ZPE校正因子的简单方法（http://sobereva.com/391）一文中提到了：
顺带一提，以后笔者有时间时候可能还会写个小程序，使得频率校正因子的拟合可以一键完成，届时用户只需填写计算级别就够了。
遂决定完成之。成品可见附件，亦发布于Github页 (https://github.com/Qemetiel/fitScFact-A-Script-that-can-Auto-fit-Scaling-Factor)
以Truhlar JCTC, 6, 2872 (2010) 中 ZPVE15/10 数据集与为例。只需已将g16正常写入环境变量，如脚本目录下运行：
- ./fitScFact.sh "b3lyp 6-31g*"
稍等，即可得到输出内容：
可见结果与帖内人工拟合结果一致。
同时生成 ZPE.txt 以供查看。
亦可在自己新增分子/新建数据集，只需要在数据集目录下添加与我相同格式的 gjf 文件，并将 fitScFact.sh 中 trainingSet 部分改为相应数据集文件夹名称即可。如我使用部分 Z1 数据集：
可见结果与常用的 0.9806 肥肠接近。
fittingScaleFactor.zip
(1.15 MB, 下载次数 Times of downloads: 47)
|
评分 Rate
-
查看全部评分 View all ratings
|
发表于 Post on 2025-4-18 18:59:31
楼主 Author