请问我这个gamess自定义基组哪里出问题了？

413

大家好，我按照之前一个大神的帖子整理了一个gamess使用自定义基组的测试输入，想请大家帮我看看哪里出问题了。
我也按照另一个大神教我的方法，从gaussian里使用pouch命令输出fort.5，但好像不能输出ecp，于是还是自己写了个输入文件
基组是从basissetexchange直接copy过来的lanl2dz的gamess-US的格式错误提示如下：
......
INPUT CARD>2         0.6274579110E+03       0.1399462700E-01                              
INPUT CARD>3         0.1429020930E+03       0.6858655181E-01                              
INPUT CARD>4         0.4023432930E+02       0.2322408730E+00                              
INPUT CARD>5         0.1282021290E+02       0.4690699481E+00                              
INPUT CARD>6         0.4390437010E+01       0.3604551991E+00                              
INPUT CARD>L   3                                                                           
INPUT CARD>1         0.1162636186E+02      -0.1149611817E+00       0.6757974388E-01        
INPUT CARD>2         0.2716279807E+01      -0.1691174786E+00       0.3239072959E+00        
    1000000 WORDS OF MEMORY AVAILABLE


     RUN TITLE
     ---------
   structure                                                                     

THE POINT GROUP OF THE MOLECULE IS C1      
THE ORDER OF THE PRINCIPAL AXIS IS     1
*** ERROR, ILLEGAL BASIS FUNCTION TYPE=N        NGAUSS=   7
ERROR OCCURED READING ATOM NO.   1, INPUT ATOM NO.    1, NAME=AG        
CHARGE= 47.0 X,Y,Z=    -5.30792526    -3.81937926     3.31861567
    LUCY: DO NOT WORRY, CHARLIE BROWN, WE LEARN MORE
          FROM LOSING THAN WE DO FROM WINNING.
CHARLIE: THEN THAT MAKES ME THE SMARTEST PERSON IN
          THE WORLD!

          ILLEGAL          BASIS FUNCTION REQUESTED.
EXECUTION OF GAMESS TERMINATED -ABNORMALLY- AT Fri Nov 29 09:36:54 2019
                580000  WORDS OF DYNAMIC MEMORY USED
CPU     0: STEP CPU TIME=     2.24 TOTAL CPU TIME=        2.2 (    0.0 MIN)
TOTAL WALL CLOCK TIME=        4.6 SECONDS, CPU UTILIZATION IS  49.34%
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
输入如下:
$CONTRL scftyp=rhf runtyp=energy dfttyp=revTPSS ICHARG=0 MULT=1 $END
$basis basnam=metal,ligN,ligH,ligH  $end
$DATA
  structure
  C1 1
Ag  47.00                 -5.30792526   -3.81937926    3.31861567
N    7.00                 -0.77529728    6.45282669    3.34426673
H    1.00                 -0.71804967    0.99897898    7.22207219
H    1.00                  1.15071639    6.75932473    2.56687569
$END
$metal
S   4
1         2.9500000              1.0141125
2         2.1490000             -1.2413971
3         0.6684000             -0.4901427
4         0.0997000              1.1128375
S   3
1         2.9500000             -1.7910564
2         2.1490000              2.0244570
3         0.6684000              0.6072839
S   1
1         0.0347000              1.0000000
P   3
1         6.5530000             -0.1079117
2         1.5650000              0.7403645
3         0.5748000              0.3721008
P   2
1         0.9085000             -0.0418371
2         0.0833000              1.0087586
P   1
1         0.0252000              1.0000000
D   3
1         3.3910000              0.1396938
2         1.5990000              0.4744421
3         0.6282000              0.5156311
D   1
1         0.2108000              1.0000000

$end
$ligN
S   6
1         0.4173511460E+04       0.1834772160E-02
2         0.6274579110E+03       0.1399462700E-01
3         0.1429020930E+03       0.6858655181E-01
4         0.4023432930E+02       0.2322408730E+00
5         0.1282021290E+02       0.4690699481E+00
6         0.4390437010E+01       0.3604551991E+00
L   3
1         0.1162636186E+02      -0.1149611817E+00       0.6757974388E-01
2         0.2716279807E+01      -0.1691174786E+00       0.3239072959E+00
3         0.7722183966E+00       0.1145851947E+01       0.7408951398E+00
L   1
1         0.2120314975E+00       0.1000000000E+01       0.1000000000E+01
D   1
1         0.8000000000E+00       1.0000000

$end
$ligH
S   3
1         0.1873113696E+02       0.3349460434E-01
2         0.2825394365E+01       0.2347269535E+00
3         0.6401216923E+00       0.8137573261E+00
S   1
1         0.1612777588E+00       1.0000000

$end
$ligH
S   3
1         0.1873113696E+02       0.3349460434E-01
2         0.2825394365E+01       0.2347269535E+00
3         0.6401216923E+00       0.8137573261E+00
S   1
1         0.1612777588E+00       1.0000000

$end
$ecp
Ag-ECP GEN    28    3
5     ----- f-ul potential -----
     -0.0587930       0     568.7006237
    -20.1145146       1     162.3579066
   -104.2733114       2      51.1025755
    -40.4539787       2      16.9205822
     -3.4420009       2       6.1669596
5     ----- s-f potential -----
      2.9861527       0      76.0974658
     35.1576460       1      15.3327359
    450.1809906       2      18.7715345
   -866.0248308       2      13.3663294
    523.1110176       2       9.8236948
5     ----- p-f potential -----
      4.9640671       0      56.3318043
     21.5028219       1      69.0609098
    546.0275453       2      19.2717998
   -600.3822556       2      12.5770654
    348.2949289       2       8.7956670
5     ----- d-f potential -----
      3.0467486       0      53.4641078
     23.3656705       1      40.1975457
    777.2540117       2      11.9086073
  -1238.8602423       2       9.7528183
    608.0677121       2       8.1788997
N-ecp none
H-ecp none
H-ecp none
$end

zjxitcc

“从gaussian里使用pouch命令输出fort.5”有诸多错误：（1）是punch=gamess；（2）文件叫fort.7

fort.7确实无法产生赝势，这是高斯的bug。自己到EMSL网站上拷贝一下GAMESS格式的赝势即可。高斯log文件里可以看到赝势数据，一般与EMSL上一模一样。可以对照之。

这种半自动半手动的写法 出错的概率 比 手动在GAMESS里自定义基组出错的概率 低多了。


我上传了我做的GAMESS输入文件，从高斯收敛的RHF传给GAMESS，一圈即能收敛。注意Ag的原子序数要写47，开头要写PP=READ。
链接: https://pan.baidu.com/s/1HBlogP2gnNi53Kct4Bclmw
提取码: 8snx

（1天有效）

413

zjxitcc 发表于 2019-11-29 13:12
“从gaussian里使用pouch命令输出fort.5”有诸多错误：（1）是punch=gamess；（2）文件叫fort.7

fort.7 ...

sorry, 原贴中误写了诸多错误，谢谢您的更正。

我已完全弃暗投明，谢谢谢谢

413

zjxitcc 发表于 2019-11-29 13:12
“从gaussian里使用pouch命令输出fort.5”有诸多错误：（1）是punch=gamess；（2）文件叫fort.7

fort.7 ...

不好意思哈，再请教一下
我后面自己补上ECP信息后碰到了这个错误提示，请问您碰到过吗？
          ECP POTENTIALS
          --------------
**** ERROR READING VARIABLE IZCORE   CHECK COLUMN 13
5     ----- f-ul potential -----
....V....1....V....2....V....3....V....4....V....5....V....6....V....7....V....8
EXECUTION OF GAMESS TERMINATED -ABNORMALLY- AT Sun Dec  1 19:26:39 2019
                580000  WORDS OF DYNAMIC MEMORY USED
CPU     0: STEP CPU TIME=     0.05 TOTAL CPU TIME=        0.0 (    0.0 MIN)
TOTAL WALL CLOCK TIME=        0.1 SECONDS, CPU UTILIZATION IS  55.56%

413

413 发表于 2019-12-1 19:29
不好意思哈，再请教一下
我后面自己补上ECP信息后碰到了这个错误提示，请问您碰到过吗？
          ECP ...

找到原因了，应该是因为Ag 19没换成Ag 47。
请问我这个是因为内存不足？
我在算的是一个Ag的配合物，有十几个Ag原子, Ag原子用的lanl2dz，还有一些轻原子应该可以忽略。
$CONTRL SCFTYP=RHF RUNTYP=ENERGY ICHARG= 0 MULT=1 UNITS=BOHR NOSYM=1 PP=READ MAXIT=200 dfttyp=revTPSS $END
$DATA
#p RHF genecp nosymm int=nobasistransform punch=gamess                          
C1       1
Ag        47.        0.0000000000        0.0000000000        3.9020349639
    S    3
    1  0.295000000E+01  -.179105558E+01  0.000000000E+00
    2  0.214900000E+01  0.202445653E+01  0.000000000E+00
    3  0.668400000E+00  0.607283758E+00  0.000000000E+00

...........


INPUT CARD>Ag        47.       -3.7112021231       -3.7112021231        3.7112021231      
INPUT CARD>    S    3                                                                     
INPUT CARD>    1  0.295000000E+01  -.179105558E+01  0.000000000E+00                        
INPUT CARD>    2  0.214900000E+01  0.202445653E+01  0.000000000E+00                        
INPUT CARD>    3  0.668400000E+00  0.607283758E+00  0.000000000E+00                        
INPUT CARD>    S    4                                                                     
INPUT CARD>    1  0.295000000E+01  0.101411151E+01  0.000000000E+00                        
INPUT CARD>    2  0.214900000E+01  -.124139640E+01  0.000000000E+00                        
INPUT CARD>    3  0.668400000E+00  -.490142462E+00  0.000000000E+00                        
INPUT CARD>    4  0.997000000E-01  0.111283746E+01  0.000000000E+00                        
INPUT CARD>    S    1                                                                     
INPUT CARD>    1  0.347000000E-01  0.100000000E+01  0.000000000E+00                        
INPUT CARD>    P    3                                                                     
INPUT CARD>    1  0.655300000E+01  -.107911702E+00  0.000000000E+00                        
INPUT CARD>    2  0.156500000E+01  0.740364515E+00  0.000000000E+00                        
INPUT CARD>    3  0.574800000E+00  0.372100807E+00  0.000000000E+00                        
DDI Process 11: semget return an error.
semget errno=ENOSPC -- check system limit for sysv semaphores.
DDI Process 17: semget return an error.
semget errno=ENOSPC -- check system limit for sysv semaphores.
DDI Process 20: semget return an error.
DDI Process 26: semget return an error.
semget errno=ENOSPC -- check system limit for sysv semaphores.

zjxitcc

这样贴也看不出太多信息。只能笼统地说：
（1）在$SYSTEM里指定内存，具体看手册；
（2）$CONTRL一行不要超过80个字符，写不下的可以回车切换到下一行（同样，最好不要顶格写）。

413

zjxitcc 发表于 2019-12-1 20:20
这样贴也看不出太多信息。只能笼统地说：
（1）在$SYSTEM里指定内存，具体看手册；
（2）$CONTRL一行不要 ...

好的，感谢，我是为了帖子里看的方便放到了一行
我看下手册里指定内存的用法，有问题再请教，谢谢谢谢