有机小分子不对称催化体系，如何选择计算精度

typ_gauss

本人对于一个有机小分子不对称催化体系进行了研究，采用如下方法计算，SMD(solvent)/B3LYP-D3(BJ)/6-311+G**(SDD for Br)//B3LYP-D3(BJ)/6-31G*(SDD for Br) level of theory at 218.15 K. 实验能量差为1.4 kcal/mol，计算为1.7kcal/mol，审稿人说计算精度不够，这个真的精度不够吗？如果只把单点能换成SMD(solvent)M062X-D3/def2TZVPP(SDD for Br)计算可以解决吗？原子近200个，重新换方法优化气相结构感觉有点浪费时间。
顺便想一问下，当拿到一个需要计算的反应体系，请问一般如何判断需要什么精度呢？谢谢

KazusaT

几kcal/mol已经属于B3LYP-D3(BJ)的误差范围了吧，能跟实验对的好不一定说明能量精度高
坚持使用B3LYP算有机体系能量的人的下场

像这种体系（几十个原子），用wB97XD或B3LYP-D3(BJ)结合6-311G**或def-TZVP优化，然后M06-2X或wB97M-V或双杂化泛函结合def2-TZVP算单点，妥妥的，耗时不很高精度又满足一般研究需要，哪个审稿人也挑不出毛病（真碰见有找茬的可以轻易怼回去）。

sobereva

审稿人缺乏基本的计算化学常识。如此小的误差，哪来的精度不够，直接怼回去
当前的结果和实验差异那么小，已经很可能有明显的误差幸运抵消的运气成分了，换成原理上精度更高的级别反倒误差还可能更大

看下文以及提到的benchmark文章了解计算精度的常识
谈谈量子化学研究中什么样的benchmark才有意义
http://sobereva.com/554（http://bbs.keinsci.com/thread-17826-1-1.html）