计算化学公社

标题: Gaussian、NWChem和ORCA用MP2在计算负离子VDE时差别 [打印本页]

作者
Author: 李潞琪 时间: 2023-12-4 00:50
标题: Gaussian、NWChem和ORCA用MP2在计算负离子VDE时差别
老师您好，我在学习阶段遇到一些问题，在计算硫酸氢根离子HSO4-的垂直剥离能VDE的时候，用三个不同软件Gaussian、NWChem和ORCA分别在MP2/6-31++G**//MP2/aug-cc-pVTZ计算了HSO4-的VDE，但是得到的结果ORCA和NWChem的结果几乎是一样的，但是Gaussian的结果能相差多达0.7eV。

(, 下载次数 Times of downloads: 2)
所以想请教各位老师，他们的差别来源是什么。 (, 下载次数 Times of downloads: 6)

作者
Author: sobereva 时间: 2023-12-4 02:18
谈谈不同量子化学程序计算结果的差异问题
http://sobereva.com/573（http://bbs.keinsci.com/thread-20122-1-1.html）

作者
Author: zjxitcc 时间: 2023-12-4 08:49
本帖最后由 zjxitcc 于 2023-12-4 08:50 编辑

经典问题。球谐型/Cartesian型基函数不同造成的，原理看这两篇
《谈谈5d、6d型d壳层基函数与它们在Gaussian中的标识》http://sobereva.com/51
《球谐型与笛卡尔型Gaussian函数的转换关系》http://sobereva.com/97

ORCA默认（只支持）5D 7F，NWCHEM默认6D 10F（在我看来这是很不好的默认），Gaussian对6-31++G(d,p)默认6D 7F，而对aug-cc-pVTZ默认5D 7F，所以你用的两个不公平的Gaussian结果在进行对比。

顺带一提，将SDD写在关键词那一行默认6D 10F，而用genecp后将SDD写在坐标底下默认5D 7F，这个估计不少做赝势计算的人会被/曾经被坑。

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