如何把Gaussian输入文件转成gamess-us的输入文件

498746012 · 发表于 Post on 2021-3-25 09:25:14

各位老师好，想请教一下如何把Gaussian计算开壳层单重态收敛的chk文件转成gamess-us的输入文件，我也想用gamess-us计算开壳层单重态并且得到同样的能量，如下是高斯的输入文件
%nproc=36
%mem=30GB
%chk=1.chk
# 5d 7f
# guess=mix ub3lyp/genecp em=gd3 int=(ultrafine,acc2e=11) scf=xqc

Title Card Required

0 1
坐标

C H 0
6-31G**
****
Pd 0
SDD
****

Pd 0
sdd

zjxitcc · 发表于 Post on 2021-3-25 09:32:43

本帖最后由 zjxitcc 于 2021-5-7 09:31 编辑

高斯输入文件关键词加上int=nobasistransform nosymm，算完后获得fchk文件

formchk xxx.chk xxx.fchk

复制代码

从fchk文件生成GAMESS .inp文件，SCF可迅速收敛（注意你用的是UB3LYP，GAMESS与高斯DFT积分格点设置不同，不能1圈收敛，但能迅速收敛；若用的是UHF，可1圈收敛）

fch2inp xxx.inp -uhf

复制代码

fch2inp小程序免费开源，可在https://gitlab.com/jxzou/mokit下载。

2021-05-07更新：前段时间MOKIT代码已更新，用户现不需要指定-uhf，会自动识别fch文件里有无Beta轨道。

biogon · 发表于 Post on 2021-3-25 09:33:15

用mokit 还需要nosymm int=nobasistransform关键词

498746012 · 发表于 Post on 2021-3-25 09:38:24

zjxitcc 发表于 2021-3-25 09:32
高斯输入文件关键词加上int=nobasistransform nosymm，算完后获得fchk文件

从fchk文件生成GAMESS .inp文 ...

谢谢老师的耐心帮助

498746012 · 发表于 Post on 2021-3-25 09:38:48

biogon 发表于 2021-3-25 09:33
用mokit 还需要nosymm int=nobasistransform关键词

谢谢老师

498746012 · 发表于 Post on 2021-3-25 09:56:24

本帖最后由 498746012 于 2021-3-25 09:57 编辑

zjxitcc 发表于 2021-3-25 09:32
高斯输入文件关键词加上int=nobasistransform nosymm，算完后获得fchk文件

从fchk文件生成GAMESS .inp文 ...

老师打扰了，我的mkl在/share/intel/2019u5下边，而which gfortran命令后是/usr/bin/gfortran，gfortran的版本是4.8.5，我在装mokit的时候make all报错,这种情况该如何处理呀？ifort -c file_op.f90 -o file_op.o -O2 -fpp -traceback
make: ifort: Command not found
make: *** [file_op.o] Error 127

snljty · 发表于 Post on 2021-3-25 09:58:42

498746012 发表于 2021-3-25 09:56
老师打扰了，我的mkl在/share/intel/2019u5下边，而which gfortran命令后是/usr/bin/gfortran，gfortran ...

看Makefile，有ifort就加载ifort的环境，用ifort编译。没有就改Makefile里面指定的编译器。

498746012 · 发表于 Post on 2021-3-25 09:59:46

本帖最后由 498746012 于 2021-3-25 10:02 编辑

snljty 发表于 2021-3-25 09:58
看Makefile，有ifort就加载ifort的环境，用ifort编译。没有就改Makefile里面指定的编译器。

谢谢老师啦，老师我的服务器工程师没有装ifort，只装了gfortan。老师我把ifor这几行注释掉，把gfortan这几行的注释去掉就可以吗# --- gfortran ---
#F90 = gfortran
#FFLAGS = -O2 -cpp -fbacktrace
#MKL_FLAGS = -lmkl_rt
#F2PY = f2py
#F2_FLAGS =
#F2_MKL_FLAGS = --link-lapack_opt

# --- ifort ----
F90 = ifort
FFLAGS = -O2 -fpp -traceback
MKL_FLAGS = -mkl
F2PY = f2py
F2_FLAGS = --fcompiler=intelem --compiler=intelem
F2_MKL_FLAGS = --link-lapack_opt

zjxitcc · 发表于 Post on 2021-3-25 10:00:43

498746012 发表于 2021-3-25 09:56
老师打扰了，我的mkl在/share/intel/2019u5下边，而which gfortran命令后是/usr/bin/gfortran，gfortran ...

你只需要fch2inp小程序，无需执行make all，执行make fch2inp即可。如果你没有ifort编译器，可以打开Makefile将前几行ifort的注释掉，且去掉gfortran的注释。就可以用gfortran了

498746012 · 发表于 Post on 2021-3-25 10:07:11

zjxitcc 发表于 2021-3-25 10:00
你只需要fch2inp小程序，无需执行make all，执行make fch2inp即可。如果你没有ifort编译器，可以打开Make ...

谢谢老师，我按照您说的，执行完make fch2inp，显示如下是不是就装成功了呀？
gfortran -c file_op.f90 -o file_op.o -O2 -cpp -fbacktrace
gfortran -c string_manipulate.f90 -o string_manipulate.o -O2 -cpp -fbacktrace
gfortran -c read_fch.f90 -o read_fch.o -O2 -cpp -fbacktrace
gfortran -c fch2inp.f90 -o fch2inp.o -O2 -cpp -fbacktrace
gfortran file_op.o string_manipulate.o read_fch.o fch2inp.o -o ../bin/fch2inp -O2 -cpp -fbacktrace

zjxitcc · 发表于 Post on 2021-3-25 10:18:35

498746012 发表于 2021-3-25 10:07
谢谢老师，我按照您说的，执行完make fch2inp，显示如下是不是就装成功了呀？
gfortran -c file_op.f90 ...

是的。

498746012 · 发表于 Post on 2021-3-25 10:58:49

zjxitcc 发表于 2021-3-25 10:18
是的。

谢谢老师成功了，老师打扰了还想问一下如何在生成的输入文件中加入B3LYP-D3, 不是D3BJ，还有坐标呀，输入文件没有B3LYP-D3和坐标
$CONTRL SCFTYP=UHF RUNTYP=ENERGY ICHARG=0 MULT=1 NOSYM=1 ICUT=11 PP=READ
MAXIT=200 ISPHER=1 $END
$SYSTEM MWORDS=500 $END
$SCF DIRSCF=.TRUE. $END
$GUESS GUESS=MOREAD NORB=767 $END
$DATA
GAMESS inp format file produced by MOKIT, nbf=767

zjxitcc · 发表于 Post on 2021-3-25 11:15:08

498746012 发表于 2021-3-25 10:58
谢谢老师成功了，老师打扰了还想问一下如何在生成的输入文件中加入B3LYP-D3, 不是D3BJ，还有坐标呀，输入 ...

你可以自行修改$CONTRL里的内容，来达到你的计算需求（泛函、D3等等），怎么写方法名称自己看GAMESS手册。坐标和基组都有了（认真看文件），不用你再写，也不能去修改坐标和基组。

498746012 · 发表于 Post on 2021-3-25 11:16:23

zjxitcc 发表于 2021-3-25 11:15
你可以自行修改$CONTRL里的内容，来达到你的计算需求（泛函、D3等等），怎么写方法名称自己看GAMESS手册 ...

谢谢老师您辛苦啦

498746012 · 发表于 Post on 2021-3-25 12:01:10

zjxitcc 发表于 2021-3-25 11:15
你可以自行修改$CONTRL里的内容，来达到你的计算需求（泛函、D3等等），怎么写方法名称自己看GAMESS手册 ...

老师gamess输入文件如何写5d,7f呀

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[GAMESS-US] 如何把Gaussian输入文件转成gamess-us的输入文件

浏览过的版块