M062X与PWPB95计算单点能结果差异大

sai77 · 发表于 Post on 2024-3-15 13:15:54

本帖最后由 sai77 于 2024-3-15 13:18 编辑

150个原子的有机体系，包含C H O N P S元素，没有大共轭基团，结构优化级别为B3LYP(D3BJ)/def2SVP。
对比的是手性产物R/S构型的过渡态能量差值，但是计算单点用 PWPB95(D3BJ)/def2TZVP 和 M062X(D3)/def2TZVP 给出了完全相反的结果：M062X的结果与晶体结构吻合，而双杂化泛函PWPB95直接定性错误。
只看单点能的话都和实验结果吻合，但是用Shermo处理加上自由能校正后就出现差异了。
请问为什么会出现这种结果呢？之前用PWPB95泛函表现都不错，这次的结果真的出乎意外。

wzkchem5 · 发表于 Post on 2024-3-15 16:46:37

之前PWPB95泛函表现不错的体系，和当前的体系，最大的分子结构差别是什么？

sai77 · 发表于 Post on 2024-3-15 21:14:35

wzkchem5 发表于 2024-3-15 16:46
之前PWPB95泛函表现不错的体系，和当前的体系，最大的分子结构差别是什么？

差别有点大，之前的体系是过渡金属催化的有机反应

wzkchem5 · 发表于 Post on 2024-3-15 23:15:00

sai77 发表于 2024-3-15 14:14
差别有点大，之前的体系是过渡金属催化的有机反应

构象搜索做得充分吗？
可能是有的关键构象没搜到，M062X的泛函误差和构象搜索不充分导致的误差碰巧抵消了

sai77 · 发表于 Post on 2024-3-16 00:02:29

wzkchem5 发表于 2024-3-15 23:15
构象搜索做得充分吗？
可能是有的关键构象没搜到，M062X的泛函误差和构象搜索不充分导致的误差碰巧抵消 ...

体系不算很柔性，构象应该不多，劣势结构只找了一个，优势结构的构象找了五个才找到一个比劣势能量低的。
准备换wB97M-V再算单点看看

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[Gaussian/gview] M062X与PWPB95计算单点能结果差异大