计算化学公社

标题: orca计算结果和高斯不同问题 [打印本页]

作者
Author: yjf1997 时间: 2021-7-20 16:45
标题: orca计算结果和高斯不同问题
我用orca和高斯计算50原子体系的有机分子，泛函用了b3lyp和pbe，但算下来结果差了有点多，请问是输入文件的问题吗？orca不加RI时，pbe的结果相差0.0225eV，b3lyp相差了0.0533eV。orca加了RI之后差的就更多了，pbe差了0.0452eV，b3lyp差了有0.0200eV。加了RI会对结果产生这么大的误差吗？这个误差在合理范围内吗？

作者
Author: wzkchem5 时间: 2021-7-20 17:17
RI的误差在这个量级是合理的，因为它远小于泛函的误差，没有必要做得比这个更精确。
至于高斯和orca剩余的误差，你可以先试试加大格点，看看是不是格点的影响。参见http://sobereva.com/573

作者
Author: yjf1997 时间: 2021-7-20 17:52

wzkchem5 发表于 2021-7-20 17:17
RI的误差在这个量级是合理的，因为它远小于泛函的误差，没有必要做得比这个更精确。
至于高斯和orca剩余的 ...

加了grid7之后，b3lyp误差只有0.00014eV了，但pbe还是有0.014eV。其他分子我用pbe算出来误差都有点大，是pbe在orca上计算方法和高斯有什么不一样吗？会不会和orca上的b3lpy要加/g一样，高斯上的pbe在orca上是另一个写法？

作者
Author: wzkchem5 时间: 2021-7-20 18:01

yjf1997 发表于 2021-7-20 10:52
加了grid7之后，b3lyp误差只有0.00014eV了，但pbe还是有0.014eV。其他分子我用pbe算出来误差都有点大，是 ...

可能是PBE的LDA部分用得不一样。ORCA的PBE的LDA部分用的是PW91-LDA（你的输出文件里写了），高斯可能是VWN5，不过我不确定

作者
Author: hebrewsnabla 时间: 2021-7-20 18:12
给高斯用int=nobasistransform可以进一步减小差别。

但是PBE的差别还是不会像b3lyp那么小的，也许泛函细节确实有差别

作者
Author: yjf1997 时间: 2021-7-20 18:39
谢谢各位老师的解答

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