请教给位大大，为什么我的计算中Gaussian与DMol3结果差别这么大？

tjugonghao

对C-Cu键进行了扫描，发现Gaussian与DMol3计算的得出的势能曲线差异很大：
1.都选用B3LYP，Gaussian与DMol3结果差异大，这是为什么呢？
2.用DMol3计算，rPBE和B3LYP差异很小，说明差异主要来自于软件的计算方法不同？那我该选择哪个软件计算呢？

卡开发发

这种小分子用Gaussian就行，绝对没争议，DMol3的杂化泛函速度和精度都是很坑的。

ggdh

基组不一样。
另外b3lyp 有vwn3 和vwn5两个版本。
不知道dmol3的是哪个版本。

tjugonghao

卡开发发 发表于 2017-2-15 10:22
这种小分子用Gaussian就行，绝对没争议，DMol3的杂化泛函速度和精度都是很坑的。

谢谢！还有就是有点诧异dmol3里面rpbe和b3lyp计算结果如此接近

tjugonghao

ggdh 发表于 2017-2-15 10:24
基组不一样。
另外b3lyp 有vwn3 和vwn5两个版本。
不知道dmol3的是哪个版本。

确实不一样，dmol里面dnp4.4，高斯用的是6311g++，就是感概差异很大。vmn3和5区别大吗？gaussian貌似用的是3

卡开发发

tjugonghao 发表于 2017-2-15 10:32
谢谢！还有就是有点差异dmol3里面rpbe和b3lyp计算结果如此接近

要弄清差异的话，或者你可以同样比较一下PBE看看？当然导致结果差异的因素还有可能有其他的，轨道截断、积分精度等不合理的设置都会引入一些误差。

但如果同样使用B3LYP，抛开精度的问题，DMol3的效率实在是太低了。小卒前辈@yjcmwgk 以前好像有过相关的测试。

winterzen

tjugonghao 发表于 2017-2-15 10:37
确实不一样，dmol里面dnp4.4，高斯用的是6311g++，就是感概差异很大。vmn3和5区别大吗？gaussian貌似用的 ...

Gaussian中使用的基组没有加极化函数吗？
DNP是个2-zeta大小的基组，你在Gaussian中用的是个3-zeta大小的，先都用2-zeta的计算看看差别有多大，另外别忘了加极化函数。

tjugonghao

卡开发发 发表于 2017-2-15 10:46
要弄清差异的话，或者你可以同样比较一下PBE看看？当然导致结果差异的因素还有可能有其他的，轨道截断、 ...

谢谢！我用pbe试试

tjugonghao

winterzen 发表于 2017-2-15 10:54
Gaussian中使用的基组没有加极化函数吗？
DNP是个2-zeta大小的基组，你在Gaussian中用的是个3-zeta大小 ...

好建议~目前都加了极化函数（d,p），忘记写了，不过没加弥散。

winterzen

tjugonghao 发表于 2017-2-15 10:57
好建议~目前都没加极化函数，请问为什么要加呢？我看文献中计算CH3Ni小分子也没加

见社长的博文：谈谈量子化学中基组的选择
http://sobereva.com/336

tjugonghao

winterzen 发表于 2017-2-15 11:01
见社长的博文：谈谈量子化学中基组的选择
http://sobereva.com/336

不好意思啊，刚刚看错了，看成你问我加没加弥散了。。。我加了极化