计算化学公社

标题: 请问我这个gamess自定义基组哪里出问题了？ [打印本页]

作者
Author: 413 时间: 2019-11-29 09:25
标题: 请问我这个gamess自定义基组哪里出问题了？
本帖最后由 413 于 2019-11-29 09:39 编辑

大家好，我按照之前一个大神的帖子整理了一个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

作者
Author: zjxitcc 时间: 2019-11-29 13:12
本帖最后由 zjxitcc 于 2019-11-29 13:54 编辑

“从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天有效）

作者
Author: 413 时间: 2019-11-29 15:16

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

fort.7 ...

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

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

作者
Author: 413 时间: 2019-12-1 19:29

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%

作者
Author: 413 时间: 2019-12-1 19:49
本帖最后由 413 于 2019-12-1 20:05 编辑

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.

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

作者
Author: 413 时间: 2019-12-1 20:33

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

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

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