计算化学公社

标题: gamess输入文件求助 [打印本页]

作者
Author: lqhenwunai 时间: 2020-5-11 11:06
标题: gamess输入文件求助
不好意思请问一下，我的gamess输入文件哪里出问题了？

$contrl maxit=200 scftyp=rhf runtyp=ENERGY ICHARG=0 $end
$DFT DFTTYP=BP86 $END
$statpt nstep=100 $end
$SCF DIRSCF=.FALSE.  FDIFF=.FALSE. $END
$SYSTEM TIMLIM=600000 MEMORY=100000000 $END
$GUESS  GUESS=HUCKEL $END
$BASIS  GBASIS=STO NGAUSS=2 $END
$DATA
ethane canonical
C1  1
  C -0.08025675768537    -0.51288956895687    0.88172622841539
  C 1.14534084103036    0.28405341010292    0.42458076749812
  H -0.97026243070168    0.12358651342552    0.93432047335516
  H -0.29934308838676    -1.33430570967557    0.19085829957482
  H 0.07691556271032    -0.94708713817721    1.87512204623966
  H 0.98810800695018    0.71812894581634    -0.56884020825319
  H 2.03527735757673    -0.35257603926912    0.37212093305213
  H 1.36458677374802    1.10549058277534    1.11537250779105
$END

输出文件虽然显示正常结束，但是
TOTAL ELECTRON NUMBER          =       0.0000000000
能量因此也完全不对。

谢谢。

作者
Author: zjxitcc 时间: 2020-5-11 11:30
本帖最后由 zjxitcc 于 2020-5-11 12:43 编辑

元素名称与Cartesian坐标之间还有一列原子序号，注意看gamess/tests/standard/下的标准示例。
我觉得如今丝毫没有手写GAMESS输入文件的必要（除非特殊目的），直接高斯算完了，用MOKIT（gitlab.com/jxzou/mokit）里的fch2inp小程序从fchk生成.inp文件，根本不会碰到这些格式问题。

作者
Author: lqhenwunai 时间: 2020-5-11 11:55

zjxitcc 发表于 2020-5-11 11:30
元素名称与Cartesian坐标之间还有一列原子序号，注意看gamess/tests/standard/下的标准示例。
我觉得如今 ...

这下对了。谢谢。

作者
Author: lqhenwunai 时间: 2020-5-11 12:53

zjxitcc 发表于 2020-5-11 11:30
元素名称与Cartesian坐标之间还有一列原子序号，注意看gamess/tests/standard/下的标准示例。
我觉得如今 ...

不好意思再问一下，我想对Cl使用def2tzvp基组，结果gamess报错说
*** ERROR, NO MORE THAN 30 GAUSSIANS PER SHELL

请问这该怎么办？

作者
Author: sobereva 时间: 2020-5-11 13:41
初学者建议用Multiwfn直接产生GAMESS-US输入文件，即用Multiwfn载入xyz/mol/pdb/gjf/mol2/fch等任意Multiwfn支持的格式，然后进入主功能100的子功能2，用里面的相应选项产生基本的GAMESS-US输入文件，之后再根据实际需要稍微修改一下关键词就能用

应当上传完整的输入文件，便于别人检查

如果不熟悉写GAMESS-US自定义基组的输入文件，直接用Gaussian生成最省事，还能避免弄错
(, 下载次数 Times of downloads: 56)

作者
Author: snljty 时间: 2020-5-11 13:49

sobereva 发表于 2020-5-11 13:41
初学者建议用Multiwfn直接产生GAMESS-US输入文件，即用Multiwfn载入xyz/mol/pdb/gjf/mol2/fch等任意Multiwf ...

赞！感觉卢老师写的Multiwfn生成各种量化软件输入文件的功能真的帮初学省了巨量时间。拿到一个新软件先生成一个最简单的输入文件自己算一算，然后再自己查手册往后学，很快就入门了。如果拿到一个软件一无所知，面对上千页的手册（有些在线的网络还有问题），得一点儿一点儿找怎么表示哪行是关键词，泛函怎么写基组怎么写，坐标怎么写，资源控制又怎么写等等，头都大了...感觉这个功能让好多软件学习周期最低减少了半个月（ORCA减少的应该远不止）。而且似乎卢老师还在抽时间不断完善这个功能，支持的基本算法越来越多了。这个功能真是帮助太大了。

作者
Author: lqhenwunai 时间: 2020-5-11 14:26
本帖最后由 lqhenwunai 于 2020-5-11 14:34 编辑

sobereva 发表于 2020-5-11 13:41
初学者建议用Multiwfn直接产生GAMESS-US输入文件，即用Multiwfn载入xyz/mol/pdb/gjf/mol2/fch等任意Multiwf ...

谢谢。我是用avogadro生成的输入文件，以后有机会再学习使用multiwfn。我现在想对Cl原子用def2tzvp 基组。结果报错说*** ERROR, NO MORE THAN 30 GAUSSIANS PER SHELL。我看说明文档里也确实说不能超过30个。
“This program is limited to a total of 2,000 atoms.  The
total number of symmetry unique basis set shells cannot
exceed 5,000, containing no more than 20,000 Gaussian
primitives.  Each contraction must contain no more than 30
Gaussians.”

如果是的话岂不是但凡大一点的基组都不能用了？不会吧。

以下是我的输入文件：
$contrl maxit=200 scftyp=rhf runtyp=ENERGY ICHARG=0 $end
$DFT DFTTYP=BP86 $END
$statpt nstep=100 $end
$SCF DIRSCF=.FALSE.  FDIFF=.FALSE. $END
$SYSTEM TIMLIM=600000 MEMORY=100000000 $END
! $BASIS  EXTFIL=.TRUE. GBASIS=DEF2TZVP $END
$DATA
ethane canonical
C1  1
C 6    0.000000000    0.000000000    -0.796370000
6
1    13575.3496820             0.22245814352E-03
2    2035.2333680             0.17232738252E-02
3    463.22562359          0.89255715314E-02
4    131.20019598          0.35727984502E-01
5       42.853015891          0.11076259931
6       15.584185766          0.24295627626
S 2
1       6.2067138508          0.41440263448
2       2.5764896527          0.23744968655
S 1
1       0.57696339419       1.0000000
S 1
1       0.22972831358       1.0000000
S 1
1       0.95164440028E-01    1.0000000
P 4
1       34.697232244          0.53333657805E-02
2       7.9582622826          0.35864109092E-01
3       2.3780826883          0.14215873329
4       0.81433208183       0.34270471845
P 1
1       0.28887547253          .46445822433
P 1
1       0.10056823671          .24955789874
D 1
1       1.09700000          1.0000000
D 1
1       0.31800000          1.0000000
F 1
1       0.76100000          1.0000000

C 6    0.000000000    0.000000000    0.796370000
S 6
1    13575.3496820             0.22245814352E-03
2    2035.2333680             0.17232738252E-02
3    463.22562359          0.89255715314E-02
4    131.20019598          0.35727984502E-01
5       42.853015891          0.11076259931
6       15.584185766          0.24295627626
S 2
1       6.2067138508          0.41440263448
2       2.5764896527          0.23744968655
S 1
1       0.57696339419       1.0000000
S 1
1       0.22972831358       1.0000000
S 1
1       0.95164440028E-01    1.0000000
P 4
1       34.697232244          0.53333657805E-02
2       7.9582622826          0.35864109092E-01
3       2.3780826883          0.14215873329
4       0.81433208183       0.34270471845
P 1
1       0.28887547253          .46445822433
P 1
1       0.10056823671          .24955789874
D 1
1       1.09700000          1.0000000
D 1
1       0.31800000          1.0000000
F 1
1       0.76100000          1.0000000

Cl  17    1.458260000    -0.841930000    1.400770000
S 7
1    69507.9909450             0.54314897497E-03
2    10426.1568800             0.41990463961E-02
3    2373.2334061             0.21592141679E-01
4    671.56420071          0.84598850094E-01
5    218.41999790          0.24757249724
6       77.572249714          0.47016930228
7       28.888815277          0.37436370716
S 3
1    127.10527185          0.25182166603E-01
2       39.339582961          0.10786112456
3       7.6740679989       -0.27408821574
S 2
1       3.8745627630          1.3213875014
2       1.8385832573          0.68636955368
S 1
1       0.50229057542       1.0000000
S 1
1       0.17962723420       1.0000000
P 5
1    666.50423284          0.23632663836E-02
2    157.64241690          0.18879300374E-01
3       50.262520978          0.87206341273E-01
4       18.536078105          0.25285612970
5       7.2940532777          0.43507154820
P 1
1       2.9433248995          .35026513165
P 1
1       1.0404970818          1.0000000
P 1
1       0.38456415080       1.0000000
P 1
1       0.13069642732       1.0000000
D 2
1       4.61000000          0.20000000
2       1.01100000          1.00000000
D 1
1       0.339000000          1.0000000
F 1

Cl  17    -1.458260000    -0.841930000    1.400770000
S 7
1    69507.9909450             0.54314897497E-03
2    10426.1568800             0.41990463961E-02
3    2373.2334061             0.21592141679E-01
4    671.56420071          0.84598850094E-01
5    218.41999790          0.24757249724
6       77.572249714          0.47016930228
7       28.888815277          0.37436370716
S 3
1    127.10527185          0.25182166603E-01
2       39.339582961          0.10786112456
3       7.6740679989       -0.27408821574
S 2
1       3.8745627630          1.3213875014
2       1.8385832573          0.68636955368
S 1
1       0.50229057542       1.0000000
S 1
1       0.17962723420       1.0000000
P 5
1    666.50423284          0.23632663836E-02
2    157.64241690          0.18879300374E-01
3       50.262520978          0.87206341273E-01
4       18.536078105          0.25285612970
5       7.2940532777          0.43507154820
P 1
1       2.9433248995          .35026513165
P 1
1       1.0404970818          1.0000000
P 1
1       0.38456415080       1.0000000
P 1
1       0.13069642732       1.0000000
D 2
1       4.61000000          0.20000000
2       1.01100000          1.00000000
D 1
1       0.339000000          1.0000000
F 1

Cl  17    0.000000000    1.683860000    1.400770000
S 7
1    69507.9909450             0.54314897497E-03
2    10426.1568800             0.41990463961E-02
3    2373.2334061             0.21592141679E-01
4    671.56420071          0.84598850094E-01
5    218.41999790          0.24757249724
6       77.572249714          0.47016930228
7       28.888815277          0.37436370716
S 3
1    127.10527185          0.25182166603E-01
2       39.339582961          0.10786112456
3       7.6740679989       -0.27408821574
S 2
1       3.8745627630          1.3213875014
2       1.8385832573          0.68636955368
S 1
1       0.50229057542       1.0000000
S 1
1       0.17962723420       1.0000000
P 5
1    666.50423284          0.23632663836E-02
2    157.64241690          0.18879300374E-01
3       50.262520978          0.87206341273E-01
4       18.536078105          0.25285612970
5       7.2940532777          0.43507154820
P 1
1       2.9433248995          .35026513165
P 1
1       1.0404970818          1.0000000
P 1
1       0.38456415080       1.0000000
P 1
1       0.13069642732       1.0000000
D 2
1       4.61000000          0.20000000
2       1.01100000          1.00000000
D 1
1       0.339000000          1.0000000
F 1

Cl  17    -1.458260000    0.841930000    -1.400770000
S 7
1    69507.9909450             0.54314897497E-03
2    10426.1568800             0.41990463961E-02
3    2373.2334061             0.21592141679E-01
4    671.56420071          0.84598850094E-01
5    218.41999790          0.24757249724
6       77.572249714          0.47016930228
7       28.888815277          0.37436370716
S 3
1    127.10527185          0.25182166603E-01
2       39.339582961          0.10786112456
3       7.6740679989       -0.27408821574
S 2
1       3.8745627630          1.3213875014
2       1.8385832573          0.68636955368
S 1
1       0.50229057542       1.0000000
S 1
1       0.17962723420       1.0000000
P 5
1    666.50423284          0.23632663836E-02
2    157.64241690          0.18879300374E-01
3       50.262520978          0.87206341273E-01
4       18.536078105          0.25285612970
5       7.2940532777          0.43507154820
P 1
1       2.9433248995          .35026513165
P 1
1       1.0404970818          1.0000000
P 1
1       0.38456415080       1.0000000
P 1
1       0.13069642732       1.0000000
D 2
1       4.61000000          0.20000000
2       1.01100000          1.00000000
D 1
1       0.339000000          1.0000000
F 1

Cl  17    0.000000000    -1.683860000    -1.400770000
S 7
1    69507.9909450             0.54314897497E-03
2    10426.1568800             0.41990463961E-02
3    2373.2334061             0.21592141679E-01
4    671.56420071          0.84598850094E-01
5    218.41999790          0.24757249724
6       77.572249714          0.47016930228
7       28.888815277          0.37436370716
S 3
1    127.10527185          0.25182166603E-01
2       39.339582961          0.10786112456
3       7.6740679989       -0.27408821574
S 2
1       3.8745627630          1.3213875014
2       1.8385832573          0.68636955368
S 1
1       0.50229057542       1.0000000
S 1
1       0.17962723420       1.0000000
P 5
1    666.50423284          0.23632663836E-02
2    157.64241690          0.18879300374E-01
3       50.262520978          0.87206341273E-01
4       18.536078105          0.25285612970
5       7.2940532777          0.43507154820
P 1
1       2.9433248995          .35026513165
P 1
1       1.0404970818          1.0000000
P 1
1       0.38456415080       1.0000000
P 1
1       0.13069642732       1.0000000
D 2
1       4.61000000          0.20000000
2       1.01100000          1.00000000
D 1
1       0.339000000          1.0000000
F 1

Cl  17    1.458260000    0.841930000    -1.400770000
S 7
1    69507.9909450             0.54314897497E-03
2    10426.1568800             0.41990463961E-02
3    2373.2334061             0.21592141679E-01
4    671.56420071          0.84598850094E-01
5    218.41999790          0.24757249724
6       77.572249714          0.47016930228
7       28.888815277          0.37436370716
S 3
1    127.10527185          0.25182166603E-01
2       39.339582961          0.10786112456
3       7.6740679989       -0.27408821574
S 2
1       3.8745627630          1.3213875014
2       1.8385832573          0.68636955368
S 1
1       0.50229057542       1.0000000
S 1
1       0.17962723420       1.0000000
P 5
1    666.50423284          0.23632663836E-02
2    157.64241690          0.18879300374E-01
3       50.262520978          0.87206341273E-01
4       18.536078105          0.25285612970
5       7.2940532777          0.43507154820
P 1
1       2.9433248995          .35026513165
P 1
1       1.0404970818          1.0000000
P 1
1       0.38456415080       1.0000000
P 1
1       0.13069642732       1.0000000
D 2
1       4.61000000          0.20000000
2       1.01100000          1.00000000
D 1
1       0.339000000          1.0000000
F 1

$END

作者
Author: sobereva 时间: 2020-5-11 14:49

lqhenwunai 发表于 2020-5-11 14:26
谢谢。我是用avogadro生成的输入文件，以后有机会再学习使用multiwfn。我现在想对Cl原子用def2tzvp 基组 ...

内容太多的时候不要直接贴，应当上传输入文件

明显是你的基组定义写得有问题。直接用我幻灯片里说的通过Gaussian的punch=gamess产生的def2-TZVP输入文件，完全没问题：
(, 下载次数 Times of downloads: 34)

作者
Author: lqhenwunai 时间: 2020-5-11 15:15

sobereva 发表于 2020-5-11 14:49
内容太多的时候不要直接贴，应当上传输入文件

明显是你的基组定义写得有问题。直接用我幻灯片里说的通 ...

谢谢！

作者
Author: 火星哥 时间: 2022-8-4 16:57

sobereva 发表于 2020-5-11 14:49
内容太多的时候不要直接贴，应当上传输入文件

明显是你的基组定义写得有问题。直接用我幻灯片里说的通 ...

老师我用Multiwfn将高斯输出的转成GAMESS的输入文件 $VEC后的数据很多很多，这正常吗

作者
Author: sobereva 时间: 2022-8-5 08:18

火星哥发表于 2022-8-4 16:57
老师我用Multiwfn将高斯输出的转成GAMESS的输入文件 $VEC后的数据很多很多，这正常吗

那是初猜轨道信息，当然正常
如果不想产生，别让Multiwfn写入初猜就完了

作者
Author: 火星哥 时间: 2022-8-26 11:28

sobereva 发表于 2022-8-5 08:18
那是初猜轨道信息，当然正常
如果不想产生，别让Multiwfn写入初猜就完了

老师，咱们培训课会讲这一部分吗

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