计算化学公社

标题: 求助在Gaussian中输入6-31g(d)和gen下面再写6-31g(d)的区别 [打印本页]

作者
Author:
630380802    时间: 2020-8-26 19:47
标题: 求助在Gaussian中输入6-31g(d)和gen下面再写6-31g(d)的区别
请教各位老师一个低级的问题,在高斯输入文件的#行直接输入6-31g(d) (举例)和在#行写gen后再在下面写6-31g(d)(对所有原子),这两种输入方式有什么区别么?为什么输出能量结果不一样?谢谢各位老师。

作者
Author:
zjxitcc    时间: 2020-8-26 19:54
本帖最后由 zjxitcc 于 2020-8-26 19:59 编辑

经典问题,我今天还在github PySCF下回答了一个这样的问题https://github.com/pyscf/pyscf/issues/688

直接写6-31g(d),高斯默认用的是6D 7F,是Cartesian/spherical harmonic混合型基函数。而在#行写genecp,底下写6-31g(d),高斯默认用的是5D 7F,全是球谐型基函数。前者会比后者能量稍低一些,因为基函数稍多一些。而cc-pVnZ和def2系列默认都是用球谐型的。如果不懂5D, 6D的含义,看sob老师博文《谈谈5d、6d型d壳层基函数与它们在Gaussian中的标识》http://sobereva.com/51

实际做计算的时候,一般无需注意这个问题。但是若有的文件直接用了6-31g(d),有的文件写的是gen或genecp,要注意这个问题,需要保持一致。要么全用球谐型(关键词5D 7F),要么全用Cartesian型(关键词6D 10F)。

作者
Author:
630380802    时间: 2020-8-26 20:49
zjxitcc 发表于 2020-8-26 19:54
经典问题,我今天还在github PySCF下回答了一个这样的问题https://github.com/pyscf/pyscf/issues/688。

...

感谢感谢!
作者
Author:
sobereva    时间: 2020-8-27 01:16
如置顶的新社员必读贴和论坛首页的公告栏所示,求助帖必须在帖子标题明确体现出此帖内容是求助或提问,并反映出帖子具体内容,避免有任何歧义,仔细看http://bbs.keinsci.com/thread-9348-1-1.html。我已把你的不恰当标题“高斯 基组”改了,以后务必注意,下次将扣分处理
作者
Author:
dingniu2    时间: 2020-8-28 21:48
zjxitcc 发表于 2020-8-26 19:54
经典问题,我今天还在github PySCF下回答了一个这样的问题https://github.com/pyscf/pyscf/issues/688。

...

你好,请问关键词5D 7F写在哪里,是写在命令行还是坐标后面,是直接写5D 7F中间带一个空格吗?谢谢!
作者
Author:
hebrewsnabla    时间: 2020-8-28 22:38
dingniu2 发表于 2020-8-28 21:48
你好,请问关键词5D 7F写在哪里,是写在命令行还是坐标后面,是直接写5D 7F中间带一个空格吗?谢谢!

关键词那一行
作者
Author:
SS123    时间: 2022-3-4 16:58
dingniu2 发表于 2020-8-28 21:48
你好,请问关键词5D 7F写在哪里,是写在命令行还是坐标后面,是直接写5D 7F中间带一个空格吗?谢谢!

请问是直接在route行写5D 7F吗?
作者
Author:
sobereva    时间: 2022-3-12 06:48
SS123 发表于 2022-3-4 16:58
请问是直接在route行写5D 7F吗?


作者
Author:
翟yingxin123    时间: 2024-3-4 13:45
SS123 发表于 2022-3-4 16:58
请问是直接在route行写5D 7F吗?

同学你好,#b3lyp/6-31d 5D 7F
这5D怎么运行不了,我想用5D,现在试了几种直接写怎么都不对




欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3