计算化学公社

标题: gaussian和molpro计算结果是否有轻微的不同？ [打印本页]
作者
Author: 诸葛壹次心    时间: 2014-12-21 23:58
标题: gaussian和molpro计算结果是否有轻微的不同？
我用gaussian直接算的数据可以和http://cccbdb.nist.gov/匹配到小数点后六位，但是使用molpro的计算结果和gaussian在第四位就有不同，大概相差0.23kcal/mol，请问这个是正常的么？我用的是最简单的hf/6-31g(d)算单点能（用其他方法的时候貌似差异就更大了）。


作者
Author: sobereva    时间: 2014-12-22 02:42
换个基组，或者调节一下收敛限再试。不同程序的HF能量应当符合得非常好。
作者
Author: 诸葛壹次心    时间: 2014-12-22 03:19
本帖最后由 诸葛壹次心 于 2015-1-25 11:22 编辑
sobereva 发表于 2014-12-22 02:42
换个基组，或者调节一下收敛限再试。不同程序的HF能量应当符合得非常好。

还真是，我换了几个基组，sto-3g, 3-21g, 6-31g, cc-pvtz等都有很好的一致性，唯独6-31(d)不行，这可真是个奇葩。。。
我有在用6-31g(d)的时候减小了收敛阀值，但是结果没有变化，不知道是不是我的脚本写得不对？

  1. ***,CH4
  2. gthresh,energy=1.d-8
  3. print,basis,orbitals
  4. angstrom
  5. geometry={                !define the nuclear coordinates
  6. C1              
  7. H2                  C1            B1
  8. H3                  C1            B2    H2            A1
  9. H4                  C1            B3    H2            A2    H3            D1    0
  10. H5                  C1            B4    H2            A3    H4            D2    0
  11. }

  13.    B1=             1.07000001
  14.    B2=             1.07000000
  15.    B3=             1.07000000
  16.    B4=             1.07000000
  17.    A1=           109.47120241
  18.    A2=           109.47125087
  19.    A3=           109.47121844
  20.    D1=          -119.99998522
  21.    D2=          -120.00000058

  23. basis=6-31g(d)                  !double-zeta basis set
  24. hf
复制代码

作者
Author: sobereva    时间: 2014-12-22 04:49
诸葛壹次心 发表于 2014-12-22 03:19
还真是，我换了几个基组，sto-3g, 3-21g, 6-31g, cc-pvtz等都有很好的一致性，唯独6-31(d)不行，这可真是 ...

有可能molpro和Gaussian对6-31G*基组的定义有微小差别。Gaussian是肯定标准的。
你可以把基组定义打印出来然后对照。
作者
Author: jiangning198511    时间: 2014-12-22 12:49
本帖最后由 jiangning198511 于 2014-12-22 12:50 编辑

建议检查一下轨道占据，Molpro非常容易收敛，但可能会占据到错误的态上，我对比过两个软件，HF的结果一致到小数点第六位
作者
Author: beefly    时间: 2015-1-2 09:24
在molpro中，默认都用球函数。

在Gaussian中，cc-pvnz或用户自定义基组默认用球函数，而pople型基组默认用直角函数。对于s，p函数两种选择不影响结果，但是从d函数开始结果就不一样了，这就是为啥仅有6-31(d)结果不一致的原因。这个默认设置主要是为了重复古代的计算结果，其实却画蛇添足。

把两个程序中的球函数/直角函数设置改为一致，结果就可以重复了。
作者
Author: 诸葛壹次心    时间: 2015-1-28 08:49
本帖最后由 诸葛壹次心 于 2015-1-28 09:06 编辑
beefly 发表于 2015-1-2 09:24
在molpro中，默认都用球函数。

在Gaussian中，cc-pvnz或用户自定义基组默认用球函数，而pople型基组默认 ...

果然如此！非常感谢！
我在molpro中加上了cartesian就和gaussian一致了~
PS: 抱歉回复晚了，上次提问之后就陷入了无穷无尽的考试周当中。。。

  1. ***,CH4
  2. gthresh,energy=1.d-8
  3. print,basis,orbitals
  4. geomtyp=cartesian
  5. angstrom
  6. geometry={                                                                      !define the nuclear coordinates
  7.   C1
  8.   H2                  C1            B1
  9.   H3                  C1            B2    H2            A1
  10.   H4                  C1            B3    H2            A2    H3            D1    0
  11.   H5                  C1            B4    H2            A3    H4            D2    0
  12. }

  14.     B1=             1.07000001
  15.     B2=             1.07000000
  16.     B3=             1.07000000
  17.     B4=             1.07000000
  18.     A1=           109.47120241
  19.     A2=           109.47125087
  20.     A3=           109.47121844
  21.     D1=          -119.99998522
  22.     D2=          -120.00000058

  24. cartesian
  25. basis=6-31g(d)                                                                  
  26. hf
复制代码

作者
Author: lastzealot    时间: 2016-12-20 11:48
诸葛壹次心 发表于 2015-1-28 08:49
果然如此！非常感谢！
我在molpro中加上了cartesian就和gaussian一致了~
PS: 抱歉回复晚了，上次提问之 ...

用笛卡尔坐标就把直角函数转化为球形函数了？
作者
Author: sobereva    时间: 2016-12-20 12:41
lastzealot 发表于 2016-12-20 11:48
用笛卡尔坐标就把直角函数转化为球形函数了？

用什么高斯函数形式和用什么坐标完全两码事



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