计算化学公社

标题: B3LYP-D3(BJ) 6-31++g(d,p)的sclZPE应该如何选择？ [打印本页]

作者
Author: 花露水 时间: 2023-3-21 21:11
标题: B3LYP-D3(BJ) 6-31++g(d,p)的sclZPE应该如何选择？
通过文献及各位老师的文章中选择了B3LYP-D3(BJ) 6-31++g(d,p)来进行优化和频率计算，使用wB97XD来计算HF，后续打算通过shermo程序进行热力学统计，但是在老师推荐的查找sclZEP的网址中没有找到对应的数值，应该如何取值呢？感谢各位老师指导！

作者
Author: wzkchem5 时间: 2023-3-21 21:22
“使用wB97XD来计算HF”不是规范说法。那个叫做“使用wB97XD来计算单点能”，单点能固然是在"HF="后面读出来的，但它不能叫做HF。
取一个最接近的理论级别下的校正因子就行了。如果被审稿人质疑，就补一个这样的计算：统计一下B3LYP搭配不同double zeta基组，校正因子的范围是多少，取一个最大的，取一个最小的，分别重新计算一遍，如果结果基本不变、不影响文章结论，就证明你这里取校正因子没有取到最合适的值所引入的误差可以忽略不计。
另一个办法是自行拟合校正因子（http://sobereva.com/391），会让结果靠谱一点点，但是改进比较有限（随便取一个类似级别的一般也不会影响文章结论），另一个好处是以后可能有人专门为了用这个校正因子而引你的文章，可以给你的文章增加几个引用数

作者
Author: 花露水 时间: 2023-3-21 21:23
我做的是(H2SO4)1(H2O)n的团簇，也想请问一下各位老师，我使用B3LYP-D3(BJ) 6-31++g(d,p)和wB97XD是否恰当？

作者
Author: wzkchem5 时间: 2023-3-21 21:26

花露水发表于 2023-3-21 14:23
我做的是(H2SO4)1(H2O)n的团簇，也想请问一下各位老师，我使用B3LYP-D3(BJ) 6-31++g(d,p)和wB97XD是否恰当 ...

没必要加弥散函数，尤其是氢原子的弥散函数更没必要；如果整个研究不涉及HSO4-、(SO4)2-等物种，那么非氢原子也可以不加弥散函数。
如果n不是很大，建议用至少6-311G(d,p)优化结构；单点能建议用DLPNO-CCSD(T)或双杂化泛函

作者
Author: 花露水 时间: 2023-3-21 21:27

wzkchem5 发表于 2023-3-21 21:22
“使用wB97XD来计算HF”不是规范说法。那个叫做“使用wB97XD来计算单点能”，单点能固然是在"HF="后面读出 ...

好的，我下次注意规范用词！非常感谢老师！

作者
Author: 花露水 时间: 2023-3-21 21:35

wzkchem5 发表于 2023-3-21 21:26
没必要加弥散函数，尤其是氢原子的弥散函数更没必要；如果整个研究不涉及HSO4-、(SO4)2-等物种，那么非氢 ...

好的，谢谢老师！！！

作者
Author: sobereva 时间: 2023-3-21 23:11
你这种团簇的opt freq用B3LYP-D3(BJ)/6-311G**是首选，性价比最高，准确度也足够一般目的。频率校正因子有现成的

作者
Author: 花露水 时间: 2023-3-22 10:17

sobereva 发表于 2023-3-21 23:11
你这种团簇的opt freq用B3LYP-D3(BJ)/6-311G**是首选，性价比最高，准确度也足够一般目的。频率校正因子有 ...

好的老师！谢谢您！

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