Gaussian、NWChem和ORCA用MP2在计算负离子VDE时差别

李潞琪 · 发表于 Post on 2023-12-4 00:50:49

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

所以想请教各位老师，他们的差别来源是什么。

HSO4.rar (3.17 MB, 下载次数 Times of downloads: 6)

sobereva · 发表于 Post on 2023-12-4 02:18:59

谈谈不同量子化学程序计算结果的差异问题
http://sobereva.com/573（http://bbs.keinsci.com/thread-20122-1-1.html）

zjxitcc · 发表于 Post on 2023-12-4 08:49:23

本帖最后由 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，这个估计不少做赝势计算的人会被/曾经被坑。

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[Gaussian/gview] Gaussian、NWChem和ORCA用MP2在计算负离子VDE时差别