求助gaussian结合shermo计算pyridine的Gibbs自由能结果与文献存在误差

李唯一

各位老师好，我使用与文献相同的计算方法计算pyridine分子的gibbs自由能结果与文献略有误差，参数设置都是完全一致的，因为我算的是pyridine分子，所以并没有设置混合基组SDD，文献中算出的pyridine分子gibbs自由能为-155755.3493kcal/mol，我用shermo算出来的gibbs自由能换算成kcal/mol并且也加上了1.89kcal/mol，结果为-155755.8463kcal/mol，想问问这些误差是怎么导致的，下面我也上传了我的gaussian和shermo输入文件信息以及文献中所用的方法，恳请各位老师指点，文献链接https://doi.org/10.1021/acs.macromol.3c00959

sobereva

谈谈重复不出来计算化学文献里的数据的可能原因
http://sobereva.com/678（http://bbs.keinsci.com/thread-38911-1-1.html）
谈谈不同量子化学程序计算结果的差异问题
http://sobereva.com/573（http://bbs.keinsci.com/thread-20122-1-1.html）

仔细把以上文章看了

李唯一

Uus/pMeC6H4-/キ 发表于 2025-10-28 10:47
6-31G(d)基组有个坑是“按genecp的混合基组写法”来用与“直接写在关键词行”来用有不同表现，基函数形式不 ...

谢谢老师的指点

Uus/pMeC6H4-/キ

6-31G(d)基组有个坑是“按genecp的混合基组写法”来用与“直接写在关键词行”来用有不同表现，基函数形式不同、结果差异明显且不可比较，需对照手册https://gaussian.com/basissets/?tabid=1与讨论帖http://bbs.keinsci.com/thread-19095-1-1.html排查。

编辑：文献除了截图外还应给出doi号等便于检索的信息。泛函选得不合适，看http://sobereva.com/272了解有过渡金属的情况应该怎么用。