计算化学公社

标题: 请教给位大大，为什么我的计算中Gaussian与DMol3结果差别这么大？ [打印本页]

作者
Author: tjugonghao 时间: 2017-2-15 09:06
标题: 请教给位大大，为什么我的计算中Gaussian与DMol3结果差别这么大？
对C-Cu键进行了扫描，发现Gaussian与DMol3计算的得出的势能曲线差异很大：
1.都选用B3LYP，Gaussian与DMol3结果差异大，这是为什么呢？
2.用DMol3计算，rPBE和B3LYP差异很小，说明差异主要来自于软件的计算方法不同？那我该选择哪个软件计算呢？

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

作者
Author: ggdh 时间: 2017-2-15 10:24
基组不一样。
另外b3lyp 有vwn3 和vwn5两个版本。
不知道dmol3的是哪个版本。

作者
Author: tjugonghao 时间: 2017-2-15 10:32
本帖最后由 tjugonghao 于 2017-2-15 10:43 编辑

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

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

作者
Author: tjugonghao 时间: 2017-2-15 10:37

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

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

作者
Author: 卡开发发 时间: 2017-2-15 10:46

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

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

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

作者
Author: winterzen 时间: 2017-2-15 10:54

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

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

作者
Author: tjugonghao 时间: 2017-2-15 10:55

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

谢谢！我用pbe试试

作者
Author: tjugonghao 时间: 2017-2-15 10:57
本帖最后由 tjugonghao 于 2017-2-15 11:01 编辑

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

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

作者
Author: winterzen 时间: 2017-2-15 11:01

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

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

作者
Author: tjugonghao 时间: 2017-2-15 11:03

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

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

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