计算化学公社

标题: 关于Gaussian 16与Gaussian 09在优化基态构型的差异求助 [打印本页]

作者
Author: songxiangan 时间: 2020-9-15 10:33
标题: 关于Gaussian 16与Gaussian 09在优化基态构型的差异求助
各位老师同学好，Gaussian 16与Gaussian 09在优化基态构型时采用同一泛函与基组的情况下差异在什么范围内是正常的，我用16优化的结果算出的HOMO能级与用09算出的HOMO能级有零点几个电子伏特的差别，LOMO能级也有类似的差别，试了好几个分子都是这样，请问这种差别是正常的吗？

作者
Author: naonao5205 时间: 2020-9-15 10:38
g16默认积分格点双电子积分精度不同详见g09default关键词。

作者
Author: zjxitcc 时间: 2020-9-15 10:43
如果你的操作无误，说明你的体系很可能对于积分格点的精细程度(fine, ultrafine) 较为敏感，或者你用的泛函对此较为敏感（例如M062X）。没有“差异在什么范围内是正常的”这种说法，只要你的操作无误，直接用G16的结果就行，不用管G09的。如果对自己的计算过程不自信，可以贴出计算详细的计算步骤或过程。

至于G16与G09的区别，论坛上以往的帖子很多，搜一下就有了。想要让G16严格与G09对比，需要在g16输入文件里加上关键词g09defaults。见官网http://gaussian.com/g09def/

作者
Author: songxiangan 时间: 2020-9-15 10:43

naonao5205 发表于 2020-9-15 10:38
g16默认积分格点双电子积分精度不同详见g09default关键词。

不好意思，请问这个g09default关键词在哪里看？我试过在16里边用int(fine,acc2e=10)，结果还是有较大差异

作者
Author: hebrewsnabla 时间: 2020-9-15 10:46
本帖最后由 hebrewsnabla 于 2020-9-15 10:57 编辑

如果仅仅是格点引起的问题，（相同结构下的）单点能相差0.01~0.1eV还算正常，但是优化的结果就不好说了，也许优化出的结构相差比较大，能量会差的更多。

作者
Author: hebrewsnabla 时间: 2020-9-15 10:47
本帖最后由 hebrewsnabla 于 2020-9-15 10:48 编辑

songxiangan 发表于 2020-9-15 10:43
不好意思，请问这个g09default关键词在哪里看？我试过在16里边用int(fine,acc2e=10)，结果还是有较大差异

直接加上这个关键词就行了。

楼上上不是都把链接给你了么……

作者
Author: songxiangan 时间: 2020-9-15 11:14

zjxitcc 发表于 2020-9-15 10:43
如果你的操作无误，说明你的体系很可能对于积分格点的精细程度(fine, ultrafine) 较为敏感，或者你用的泛 ...

老师您好，想要让G16严格与G09对比，您看下这样写是否正确
%chk=C:\G16W\keti\DP3.chk
%nprocshared=20
%mem=10GB
# opt b3lyp/6-31g(d) g09defaults

......

......

作者
Author: sobereva 时间: 2020-9-15 14:35

songxiangan 发表于 2020-9-15 11:14
老师您好，想要让G16严格与G09对比，您看下这样写是否正确
%chk=C:\G16W\keti\DP3.chk
%nprocshared=20 ...

是

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