求助gamess HF能量与g16能量不一样

lqhenwunai

我用g16生成gamess的输入文件：
%mem=5000MB
%nprocshared=10
#p HF/cc-pvtz guess=only punch=gamess nosymm

1

0 1
6       0.00321900      -0.00065800     -0.76547100
6       -0.00321900     0.00065800      0.76547100
1       0.87998400      -0.51517000     1.17161100
1       -0.00033800     1.02543800      1.16659100
1       -0.89437500     -0.50725900     1.16414200
1       -0.87998400     0.51517000      -1.17161100
1       0.89437500      0.50725900      -1.16414200
1       0.00033800      -1.02543800     -1.16659100


得到以下的gamess输入文件

$CONTRL SCFTYP=RHF RUNTYP=ENERGY ICHARG= 0 MULT=1 UNITS=BOHR $END
$DATA
1
C1       1
基组和几何坐标...


但是得到的HF能量却差很多

-79.2591675536 （gamess)
-79.2589563034  (G16C)

请问这是哪里出问题了？

zjxitcc

这是能量不一致系列里的常见问题，原因很简单，你可以看输出文件，基函数数目都不一样。GAMESS用的是6D 10F，高斯默认用的是5D 7F。如果你对5D 6D不懂，可以看社长的博文http://sobereva.com/51《谈谈5d、6d型d壳层基函数与它们在Gaussian中的标识》

显然6D 10F下的基函数稍微多一些，能量也更低。如果你想保持严格一致，高斯关键词应当加上6D 10F，一般就能一致了。不过这样GAMESS里你得从头开始迭代波函数，费时。而且有UHF多个解的时候，容易出现高斯与GAMESS收敛到不同解上的问题。

如果想严格保持一致，并且在GAMESS中一两圈就收敛，可以用MOKIT（https://gitlab.com/jxzou/mokit）中的fch2inp小程序从高斯fchk文件（加关键词nosymm int=nobasistransform 6D 10F）直接产生GAMESS .inp文件，然后你用GAMESS算，就能立即得到一模一样的能量了。

biogon

zjxitcc 发表于 2020-6-8 10:10
这是能量不一致系列里的常见问题，原因很简单，你可以看输出文件，基函数数目都不一样。GAMESS用的是6D 10F ...

mokit在win下能编译吗

zjxitcc

biogon 发表于 2020-6-8 11:14
mokit在win下能编译吗

编译的话，你在Win下用Visual Studio自然是能编译。我也在考虑提供windows预编译版，MOKIT中的大部分小程序是可以提供windows版本的，即将更新的一些新功能不行。在近期的更新里我会提供windows预编译版。

biogon

zjxitcc 发表于 2020-6-8 11:26
编译的话，你在Win下用Visual Studio自然是能编译。我也在考虑提供windows预编译版，MOKIT中的大部分小程 ...

win下运行需要的库和linux的一样的吧

sobereva

Multiwfn载入fch文件，用主功能100的子功能2产生GAMESS-US输入文件，直接就可以带着Gaussian产生的波函数当初猜。Multiwfn在windows下双击即可直接使用，无需安装

GAMESS-US的初猜做得很烂，三重态水都能收敛到非稳定波函数上去

zjxitcc

biogon 发表于 2020-6-8 11:14
mokit在win下能编译吗

最新版MOKIT（https://gitlab.com/jxzou/mokit）提供 轨道转换程序 的Windows可执行文件，在\bin目录里解压缩即可使用。这样在Windows下也可以很方便直接从fchk文件生成GAMESS, ORCA, OpenMolcas等软件的输入文件和轨道。