GKS-EDA计算报错

498746012

各位老师好！想请教一下(FeH2)2的GKS-EDA如何加入赝势基组，supermolecule, monomer1,monomer2都要写赝势基组吗？以下输入文件计算报错 ERROR: $CONTRL KEYWORD PP/ECP WAS GIVEN AN ILLEGAL VALUE RE

$contrl scftyp=rhf runtyp=eda dfttyp=b3lyp
icharg=0 mult=1 maxit=200 pp=read $end
$lmoeda matom(1)=3 3 mcharg(1)=0 0 mmult(1)=1 1
edatyp=gks rdvecm=.t. itrtn=.f. $end
$scf diis=.t. soscf=.f. dirscf=.t. fdiff=.f. $end
$system memddi=1000 mwords=300 $end
$data
GAMESS inp format file produced by MOKIT
C1   1
Fe   26.         1.66648840         0.91958739         0.18516925
   S  3
   1   2.05130280E+01   2.34714123E-01   0.00000000E+00
   2   9.77679200E+00  -8.77527461E-01   0.00000000E+00
   3   4.57359900E+00  -3.41636180E-01   0.00000000E+00
.....

H     1.         2.51062821         0.45945068         0.22467864
   S  3
   1   1.87311370E+01   3.34946043E-02   0.00000000E+00
   2   2.82539436E+00   2.34726953E-01   0.00000000E+00
   3   6.40121692E-01   8.13757326E-01   0.00000000E+00
........

H     1.         1.71503388         1.44488947        -0.61955493
   S  3
   1   1.87311370E+01   3.34946043E-02   0.00000000E+00
   2   2.82539436E+00   2.34726953E-01   0.00000000E+00
   3   6.40121692E-01   8.13757326E-01   0.00000000E+00
.........

Fe   26.        -0.50007313        -1.00682957        -0.02512561
   S  3
   1   2.05130280E+01   2.34714123E-01   0.00000000E+00
   2   9.77679200E+00  -8.77527461E-01   0.00000000E+00
   3   4.57359900E+00  -3.41636180E-01   0.00000000E+00
   S  1
   1   1.94991000E+00   1.00000000E+00   0.00000000E+00
   S  1
.........

H     1.         0.18059841        -0.33345616         0.12122210
   S  3
   1   1.87311370E+01   3.34946043E-02   0.00000000E+00
   2   2.82539436E+00   2.34726953E-01   0.00000000E+00
   3   6.40121692E-01   8.13757326E-01   0.00000000E+00
..........
H     1.        -1.11879779        -0.88421939         0.69906109
   S  3
   1   1.87311370E+01   3.34946043E-02   0.00000000E+00
   2   2.82539436E+00   2.34726953E-01   0.00000000E+00
   3   6.40121692E-01   8.13757326E-01   0.00000000E+00
........

$END
$ECP
1-ECP GEN   10   3
1     ----- f-ul potential -----
    0.00000000E+00    2   1.00000000E+00
2     ----- s-f potential -----
    2.53749588E+02    2   2.09300000E+01
    3.79228450E+01    2   9.44500000E+00
2     ----- p-f potential -----
    1.61036812E+02    2   2.17600000E+01
    2.76512980E+01    2   9.17800000E+00
2     ----- d-f potential -----
   -2.44312760E+01    2   2.59000000E+01
   -1.43425100E+00    2   8.83500000E+00
2 ECP-NONE
3 ECP-NONE
4-ECP GEN   10   3
1     ----- f-ul potential -----
    0.00000000E+00    2   1.00000000E+00
2     ----- s-f potential -----
    2.53749588E+02    2   2.09300000E+01
    3.79228450E+01    2   9.44500000E+00
2     ----- p-f potential -----
    1.61036812E+02    2   2.17600000E+01
    2.76512980E+01    2   9.17800000E+00
2     ----- d-f potential -----
   -2.44312760E+01    2   2.59000000E+01
   -1.43425100E+00    2   8.83500000E+00
5 ECP-NONE
6 ECP-NONE
$END

$vec0
...... !the wave function of supermolecule.
$end
$vec1
...... !the wave function of monomer 1.
$end
$vec2
...... !the wave function of monomer 2.
$end

zjxitcc

阿甘 发表于 2022-7-26 23:44
啊好的，谢谢老师的提醒，我使用的是1.2.3版本，我现在就下载新版本试试看

运行automr -v会显示具体日期，差1天都不是最新版

阿甘

zjxitcc 发表于 2022-7-26 23:38
是否尝试过使用最新版的MOKIT？不要用一个月之前的版本

啊好的，谢谢老师的提醒，我使用的是1.2.3版本，我现在就下载新版本试试看

zjxitcc

阿甘 发表于 2022-7-26 23:32
老师您好，学生最近在使用您写的程序MOKIT中的frag_guess_wfn小程序产生GKS-EDA输入文件时出现数组越界的 ...

是否尝试过使用最新版的MOKIT？不要用一个月之前的版本

阿甘

zjxitcc 发表于 2021-6-11 23:47
利用MOKIT（https://gitlab.com/jxzou/mokit）中的frag_guess_wfn小程序可以可以自动产生GKS-EDA输入文件（ ...

老师您好，学生最近在使用您写的程序MOKIT中的frag_guess_wfn小程序产生GKS-EDA输入文件时出现数组越界的报错，请问您可以告知我如何修改源代码然后重新编译吗？以下是输出信息：
%mem=10000MB
%nprocshared=24
natom=80
nfrag0=2
nfrag=5
method=m052x, basis=6-31+g(d), wfn_type=3
i=  1, frags(i)%e =    -1331.144661210, frags(i)%ssquare=  0.78
i=  2, frags(i)%e =    -1331.144662060, frags(i)%ssquare=  0.78
i=  3, frags(i)%e =    -1331.148736200, frags(i)%ssquare=  0.78
i=  4, frags(i)%e =    -1331.148737050, frags(i)%ssquare=  0.78
ERROR in subroutine direct_sum_frag_mo2super_mo: the 2nd dimension of array mo_a out of range!
k3=644,nif0=543,na0=101
k4=1086,nif=1082,i=2,n=2

498746012

hdhxx123 发表于 2021-6-14 21:03
frag_guess_wfn比fch2inp对解决GAMESS的收敛问题要好得多

谢谢老师

hdhxx123

498746012 发表于 2021-6-14 15:06
谢谢老师，我的版本是gcc version 4.8.5. 老师能不能只用fch2inp，来生成GKS-EDA的读取轨道的输入文件呀[ ...

frag_guess_wfn比fch2inp对解决GAMESS的收敛问题要好得多

498746012

zjxitcc 发表于 2021-6-14 12:08
要把信息一次性给足。例如你没说gfortran版本。手册里写了需用>=4.8.5版本的gfortran，低的话就不必尝试 ...

谢谢老师，我的版本是gcc version 4.8.5. 老师能不能只用fch2inp，来生成GKS-EDA的读取轨道的输入文件呀

zjxitcc

498746012 发表于 2021-6-14 09:45
老师想问一下，您开发的软件，我make all,报如下错误该怎么办呢？我服务器上没有ifort，只有gfortran,
g ...

要把信息一次性给足。例如你没说gfortran版本。手册里写了需用>=4.8.5版本的gfortran，低的话就不必尝试了，可以选择安装一下intel编译器，十分简单《Linux下安装Intel oneAPI》。

PS1: MOKIT产生GKS-EDA输入文件功能确实只需要frag_guess_wfn和fch2inp两个小程序。
PS2: 感觉你这报错信息不是最初的，上面找找可能还有报错。

hdhxx123

498746012 发表于 2021-6-14 09:45
老师想问一下，您开发的软件，我make all,报如下错误该怎么办呢？我服务器上没有ifort，只有gfortran,
g ...

1.装个ifort，方法网上都能查到，这里帖一个

1. https://mp.weixin.qq.com/s/7pQETkrDO1C83vQjKQqI4w

复制代码

2.用gfortran编译，看makefile里开头的是可以选择的

498746012

zjxitcc 发表于 2021-6-12 18:01
不是你这么问题的，你看这问答 一来一回，一天时间就过去了，本来20 min就能算出结果了。

（1）当你说 ...

老师想问一下，您开发的软件，我make all,报如下错误该怎么办呢？我服务器上没有ifort，只有gfortran,
gfortran file_op.o string_manipulate.o read_gms_inp.o add_bgcharge_to_inp.o -o ../bin/add_bgcharge_to_inp -O2 -cpp -fbacktrace
gfortran -c do_hf.f90 -o do_hf.o -O2 -cpp -fbacktrace
do_hf.f90:9.22:

use mr_keyword, only: hf_prog, readuhf, readrhf, skiphf, mem, nproc, basis, &
                      1
Error: Symbol 'hf_prog' referenced at (1) not found in module 'mr_keyword'
do_hf.f90:68.25:

select case(TRIM(hf_prog))
                         1
Error: Symbol 'hf_prog' at (1) has no IMPLICIT type
do_hf.f90:68.13:

select case(TRIM(hf_prog))
             1
Error: Argument of SELECT statement at (1) cannot be UNKNOWN
make: *** [do_hf.o] Error 1

498746012

zjxitcc 发表于 2021-6-12 18:01
不是你这么问题的，你看这问答 一来一回，一天时间就过去了，本来20 min就能算出结果了。

（1）当你说 ...

谢谢老师耐心的回复。

zjxitcc

498746012 发表于 2021-6-12 15:46
老师我试过了还是不行，还是没有基组信息，估计frag_guess_wfn程序运行的时候会调用mokit-master其他的小 ...

不是你这么问题的，你看这问答 一来一回，一天时间就过去了，本来20 min就能算出结果了。

（1）当你说“不行”时，应该给出报错信息；如果自己找不到出错，就应该上传输入输出文件（可压缩后上传）。平时就算是高斯报错，你肯定也是遵循这个步骤。
（2）“只是我没有装”->如果你怀疑这个问题，直接执行make all全部编译就行了，也就1 min的事。
（3）文中Ti的例子里，cc-pVTZ对Ti而言是全电子基组，因此你看不到赝势数据是正常的。LANL2DZ对Ti才有赝势。
（4）由于你提供的信息过少，我只能提供一些模糊的回答：
(i) 确保{gks}顶格写了
(ii) 确保使用的是https://gitlab.com/jxzou/mokit页面上最新的程序，不要从Releases那里去下载，那里会滞后一些。

498746012

zjxitcc 发表于 2021-6-12 13:42
Title Card那一行你得写{gks}，这是frag_guess_wfn小程序识别关键词的地方。虽然这行对Gaussian是没用的 ...

老师我试过了还是不行，还是没有基组信息，估计frag_guess_wfn程序运行的时候会调用mokit-master其他的小程序，只是我没有装

498746012

zjxitcc 发表于 2021-6-12 13:42
Title Card那一行你得写{gks}，这是frag_guess_wfn小程序识别关键词的地方。虽然这行对Gaussian是没用的 ...

谢谢老师，感恩