计算化学公社

标题: 求助gamess HF能量与g16能量不一样 [打印本页]

作者
Author: lqhenwunai 时间: 2020-6-8 10:02
标题: 求助gamess HF能量与g16能量不一样
我用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)

请问这是哪里出问题了？

作者
Author: zjxitcc 时间: 2020-6-8 10:10
本帖最后由 zjxitcc 于 2020-6-8 10:17 编辑

这是能量不一致系列里的常见问题，原因很简单，你可以看输出文件，基函数数目都不一样。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算，就能立即得到一模一样的能量了。

作者
Author: biogon 时间: 2020-6-8 11:14

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

mokit在win下能编译吗

作者
Author: zjxitcc 时间: 2020-6-8 11:26
本帖最后由 zjxitcc 于 2020-6-8 11:27 编辑

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

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

作者
Author: biogon 时间: 2020-6-8 11:32

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

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

作者
Author: sobereva 时间: 2020-6-8 12:54
Multiwfn载入fch文件，用主功能100的子功能2产生GAMESS-US输入文件，直接就可以带着Gaussian产生的波函数当初猜。Multiwfn在windows下双击即可直接使用，无需安装

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

作者
Author: zjxitcc 时间: 2020-8-8 18:16
本帖最后由 zjxitcc 于 2020-8-8 18:32 编辑

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

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

(, 下载次数 Times of downloads: 37)

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