gaussian和molpro计算结果是否有轻微的不同？

诸葛壹次心

我用gaussian直接算的数据可以和http://cccbdb.nist.gov/匹配到小数点后六位，但是使用molpro的计算结果和gaussian在第四位就有不同，大概相差0.23kcal/mol，请问这个是正常的么？我用的是最简单的hf/6-31g(d)算单点能（用其他方法的时候貌似差异就更大了）。

sobereva

换个基组，或者调节一下收敛限再试。不同程序的HF能量应当符合得非常好。

诸葛壹次心

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. }
    
12. 
    
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
    
22. 
    
23. basis=6-31g(d)                  !double-zeta basis set
    
24. hf
    

复制代码

sobereva

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

有可能molpro和Gaussian对6-31G*基组的定义有微小差别。Gaussian是肯定标准的。
你可以把基组定义打印出来然后对照。

jiangning198511

建议检查一下轨道占据，Molpro非常容易收敛，但可能会占据到错误的态上，我对比过两个软件，HF的结果一致到小数点第六位

beefly

在molpro中，默认都用球函数。

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

把两个程序中的球函数/直角函数设置改为一致，结果就可以重复了。

诸葛壹次心

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. }
    
13. 
    
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
    
23. 
    
24. cartesian
    
25. basis=6-31g(d)                                                                  
    
26. hf

复制代码

lastzealot

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

用笛卡尔坐标就把直角函数转化为球形函数了？

sobereva

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

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