计算化学公社

标题: 有机小分子不对称催化体系,如何选择计算精度 [打印本页]

作者
Author:
typ_gauss    时间: yesterday 23:37
标题: 有机小分子不对称催化体系,如何选择计算精度
本人对于一个有机小分子不对称催化体系进行了研究,采用如下方法计算,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个,重新换方法优化气相结构感觉有点浪费时间。
顺便想一问下,当拿到一个需要计算的反应体系,请问一般如何判断需要什么精度呢?谢谢

作者
Author:
KazusaT    时间: 2 hour ago
几kcal/mol已经属于B3LYP-D3(BJ)的误差范围了吧,能跟实验对的好不一定说明能量精度高
坚持使用B3LYP算有机体系能量的人的下场
像这种体系(几十个原子),用wB97XD或B3LYP-D3(BJ)结合6-311G**或def-TZVP优化,然后M06-2X或wB97M-V或双杂化泛函结合def2-TZVP算单点,妥妥的,耗时不很高精度又满足一般研究需要,哪个审稿人也挑不出毛病(真碰见有找茬的可以轻易怼回去)。


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

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







欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3