计算化学公社

标题: 势能面计算的能量与CCSD/TZ计算的能量差别较大怎么解决？ [打印本页]

作者
Author: PESPES 时间: 2021-12-3 15:13
标题: 势能面计算的能量与CCSD/TZ计算的能量差别较大怎么解决？
使用高斯做势能面的计算，venus-nwchem做动力学模拟。体系30个原子左右，研究的是E2反应。首先用B3LYP、m06、m062x、bhandhlyp、bhandh、mp2、b971等方法，aug-cc-pvdz基组做了势能面的计算，得到的结果与CCSD(T)/aug-cc-pvTZ值比较，其中B3LYP、m062x、bhandhlyp这三种方法计算的平均偏差和标准偏差比较小，之后选择了B3LYP这个泛函去跑动力学，因为他用的时间更短一点，之后在轨迹中得到一些新的反应路径，然后又把新的路径补充到势能面上去，然而对于新的路径的过渡态b3lyp计算的值小于CCSD(T)/aug-cc-pvTZ的结果，也就是低估了势垒，最后计算的偏差反而是后两种方法更好一些。另外反应的放热性B3LYP计算的值和实验值比较吻合的。
1.我轨迹已经用B3LYP方法跑完了，速率常数与实验值吻合。

2.存在的问题就是B3LYP对新的路径的势垒描述不好，从而导致整体结果偏差，如何解决？

3.B3LYP对能量描述不好是否会影响到后续的动力学计算？

作者
Author: sobereva 时间: 2021-12-4 01:44
都用M06-2X就完了
也可以考虑wB97XD

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