计算化学公社

标题: GAMESS软件求助 [打印本页]

作者
Author: Tinker 时间: 2014-12-10 23:49
标题: GAMESS软件求助
本帖最后由 Tinker 于 2014-12-10 23:51 编辑

求助各位大大，有两个问题，

1. 第一个问题是关于GAMESS里边怎么使用混合基组，比如高斯里边我可以写，
B3LYP/GENECP 。。。。

下边写
Ir 0
Lanl2dz
****
C H 0
6-31g*
****

Ir 0
Lanl2dz
，

而在GAMESS，他们给我的输入文件是这样的在分子信息后边写价电子基组，之后再写个$ECP去写赝势
IRIDIUM 77.0    9.277565000       1.202463000       6.804552000
S 3
  1    2.3500000          -1.6784642
  2    1.5820000             2.0952553
  3    0.5018000             0.4162934
S 4
  1    2.3500000             1.6464467
  2    1.5820000          -2.2748150
  3    0.5018000          -1.0494357
  4    0.2500000             1.2167791
S 1
  1    0.0598000             1.0000000
P 3
  1    2.7920000          -0.3889212
  2    1.5410000             0.9077516
  3    0.5285000             0.4691443
P 2
  1    0.5100000          -0.1170669
  2    0.0980000             1.0489002
P 1
  1    0.0290000             1.0000000
D 2
  1    1.2400000             0.5087022
  2    0.4647000             0.5862102
D 1
  1    0.1529000             1.0000000

$ECP
  IR-ECP GEN    60    4
  5    ----- g-ul potential    -----
   -0.1578014       1 823.5880147
  -1517.5270446       2 364.6613336
-316.5306529       2    55.7082801
-91.8880941       2    12.0464544
   -9.2241773       2    3.5120610
  6    ----- s-ul potential    -----
   3.1578014       0 188.0490770
   26.8322577       1 340.4194712
800.4250007       2 128.2373673
369.4050683       2    33.8644961
242.4171899       2    4.7560005
-118.2173282       2    3.9649974
  5    ----- p-ul potential    -----
   2.1578014       0 289.7291139
   61.9678610       1    87.4633789
269.0581986       2    30.4363766
231.1654793       2    4.0553412
-133.6952667       2    3.5525341
  5    ----- d-ul potential    -----
   3.1578014       0 136.4017106
   45.9349803       1    95.0776925
359.0344668       2    49.2258410
176.4740119       2    15.0874145
   54.5155286       2    4.0405764
  5    ----- f-ul potential    -----
   3.9546197       0 127.3507908
   52.9773655       1    66.2364374
274.8643383       2    34.4299229
137.2047338       2    10.1995721
   14.8633305       2    2.5409702
O-ECP NONE
C-ECP  NONE
C-ECP  NONE
C-ECP  NONE
H NONE
H NONE
H NONE
H NONE
H NONE
H NONE
H NONE
H NONE
H NONE
O-ECP NONE
$END

一定要通过这种方法写混合基组吗？

还有一个问题，最近在练习使用GAMESS做MCSCF的计算，我师兄的做法是，是先做HF计算，得到的轨道信息作为MCSCF的初始轨道，贴在MCSCF输入文件的下边，是否一定需要这样？

作者
Author: jjjspring 时间: 2014-12-11 00:15
GAMESS里已经定义好的ECP很少，大部分要像你所说的方法输入，具体有哪些你可以看手册。

是的。需要guess=moread 然后用$VEC提供初始轨道，就像Gaussian里面guess=read，不过轨道信息在check文件中，而不是写在输入文件里。做casscf也常常需要guess=（read，alter），gamess貌似没有alter这个命令，可以手动交换。

作者
Author: Tinker 时间: 2014-12-11 01:42

jjjspring 发表于 2014-12-11 00:15
GAMESS里已经定义好的ECP很少，大部分要像你所说的方法输入，具体有哪些你可以看手册。

是的。需要guess ...

谢谢！我再试试。

找了两个说明先贴在这里
http://www.spec.org/hpc96/docs/R ... s/gamess/input.html
http://myweb.liu.edu/~nmatsuna/gamess/input/ORMAS.html

作者
Author: beefly 时间: 2015-1-2 09:35
基组部分：可以设置rungms脚本中的外部基组文件变量，然后就可以在$BASIS中直接引用自定义基组名称。不过近两年的rungms脚本把这个环境变量丢掉了，需要用两年前的rungms脚本。

ECP部分：没办法，目前还得对每个原子（不是元素！）挨个定义

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