求助最多40个原子有机反应的势垒的计算级别选择

Roa

各位老师好，我希望使用高斯计算一个仅含主族元素的有机反应的势垒，体系不算很大（最多40个原子），结构里有氢键，目前打算在B3LYP-D3(BJ)/6-311G**水平下做结构优化与振动分析，再使用双杂化泛函结合def2-TZVP基组计算高精度单点能。但是在查阅文献的时候我看到有文章计算类似反应使用了PBE0/def2-TZVP-D3BJ。请问我这个计算级别是否合适？需要采用文献中的计算方法吗？

cokie

都有参考文献了，直接照抄用 PBE0-D3BJ 多好，PBE0-D3BJ 也不比 B3LYP-D3BJ 慢。写方法选择的时候还可以顺手一引，直接杜绝审稿人刁难你的计算级别选择。高精度单点的选择倒是没有问题，不过个人愚见：在 ORCA 里用 wB97M-V 我感觉精度很好而且比双杂化快得多。

UW_0728.

什么叫def2-TZVP-D3BJ，写对了，PBE0-D3(BJ)/def2-TZVP

优化级别我觉得没啥问题；算单点的话，我觉得有条件用双杂化的话当然双杂化更好（不过只在Gaussian上算的话性价比不高，建议用ORCA开RI算双杂化），我感觉对当前类型问题revDSD-PBEP86-D3(BJ)（ORCA里用D4）就很不错了；你提到的文献中PBE0-D3(BJ)/def2-TZVP算单点也是合理的级别（适合一类体系计算的级别没必要是唯一的，只要你恰当考虑了方法级别本身的特点，以及你自己的软件和硬件条件）。一般来说，如果只用Gaussian算仅由常见主族元素构成的有机体系的话，用M06-2X/def2-TZVP算单点进而得到热力学数据是算势垒方面性价比最高的，这个组合也是最常用的。

Roa

UW_0728. 发表于 2026-2-2 17:12
什么叫def2-TZVP-D3BJ，写对了，PBE0-D3(BJ)/def2-TZVP

优化级别我觉得没啥问题；算单点的话，我觉得有 ...

老师好，这个def2-TZVP-D3BJ文献里确实就是这么写的，我也觉得很奇怪。

感谢您的指导，我先在高斯里试试M06-2X/def2-TZVP

Roa

cokie 发表于 2026-2-2 17:12
都有参考文献了，直接照抄用 PBE0-D3BJ 多好，PBE0-D3BJ 也不比 B3LYP-D3BJ 慢。写方法选择的时候还可以顺 ...

好的，谢谢您的指导

sobereva

认真把此文看了，该说的都说了
简谈量子化学计算中DFT泛函的选择
http://sobereva.com/272（http://bbs.keinsci.com/thread-536-1-1.html）

PBE0算有机体系热力学量是明显不理想的选择，绝对不要随便效仿文献（2/3以上做应用性计算的研究者是计算级别选择方面的外行或准外行）。有的是耗时相仿佛而精度好得多的普通泛函可用，上文说了。参考北京科音中级量子化学培训班（http://www.keinsci.com/KBQC）讲泛函选择的这页幻灯片：

sobereva

Roa 发表于 2026-2-2 17:19
老师好，这个def2-TZVP-D3BJ文献里确实就是这么写的，我也觉得很奇怪。

感谢您的指导，我先在高斯里试 ...

明显作者属于业余水平。这种文章的计算级别的选择毫无参考性可言

Roa

sobereva 发表于 2026-2-2 23:53
认真把此文看了，该说的都说了
简谈量子化学计算中DFT泛函的选择
http://sobereva.com/272（http://bbs.k ...

感谢卢老师的讲解！