计算化学公社

标题: Gamess如何读取赝势基组的自定义$ECP部分 [打印本页]

作者
Author: molx79 时间: 2015-4-3 09:05
标题: Gamess如何读取赝势基组的自定义$ECP部分
本帖最后由 molx79 于 2015-4-3 09:10 编辑

本人才学疏浅，求各位老师、大师不吝赐教。我做的体系：SI2....HCHO的能量分解。其中Ｉ原子选取的基组是ａｕｇ－ｃｃ－ｐｖｔｚ－ｐｐ，其余的原子选用的基组是ａｕｇ－ｃｃ－ｐｖｔｚ。如果Ｉ原子选用系统自带的SBKJC就能算成：
$ECP
I-ECP SBKJC
I-ECP SBKJC
O-ecp none
S-ecp none
C-ecp none
H-ecp none
H-ecp none
$END    。

由于我高斯计算的Ｉ采用ａｕｇ－ｃｃ－ｐｖｔｚ－ｐｐ。所以我把网页的https://bse.pnl.gov/bse/portal查到的Ｉ的：
I-ECP GEN    28    4  1    ----- g-ul potential    -----    0.00000000    2    1.00000000
3    ----- s-ul potential    -----
49.98964900    2    40.03337600
281.00655600    2    17.30057600
61.41673900    2    8.85172000
4    ----- p-ul potential    -----
67.41623900    2    15.72014100
134.80769600    2    15.20822200
14.56654800    2    8.29418600
28.96842200    2    7.75394900
4    ----- d-ul potential    -----
35.53875600    2    13.81775100
53.33975900    2    13.58780500
9.71646600    2    6.94763000
14.97750000    2    6.96009900
4    ----- f-ul potential    -----
-20.17661800    2    18.52295000
-26.08807700    2    18.25103500
-0.22043400    2 　　7.55790100
-0.22164600    2    7.59740400
按照ｇａｍｅｓｓ说明书的内容，改为：
I-ECP GEN    28    4
1 ----- IODINE U(g) -----
0.00000000    2    1.00000000
3    ----- IODINE U(s)-u(g) -----
49.98964900    2    40.03337600
281.00655600    2    17.30057600
61.41673900    2    8.85172000
4    ----- IODINE U(p)-u(g)    -----
67.41623900    2    15.72014100
134.80769600    2    15.20822200
14.56654800    2    8.29418600
28.96842200    2    7.75394900
4    ----- IODINE U(d)-u(g)    -----
35.53875600    2    13.81775100
53.33975900    2    13.58780500
9.71646600    2    6.94763000
14.97750000    2    6.96009900
4    ----- IODINE U(f)-u(g)    -----
-20.17661800    2    18.52295000
-26.08807700    2    18.25103500
-0.22043400    2    7.55790100
-0.22164600    2    7.59740400
但是Ｇａｍｅｓ读取不了，给的错误是
**** ERROR READING VARIABLE IZCORE CHECK COLUMN 16 1    ----- IODINE U(g)    -----                                           ....V....1....V....2....V....3....V....4....V....5....V....6....V....7....V....8 EXECUTION OF GAMESS TERMINATED -ABNORMALLY- AT Fri Apr  3 08:40:01 2015             580000  WORDS OF DYNAMIC MEMORY USED请大家不吝赐教，拜托

作者
Author: molx79 时间: 2015-4-3 09:07
有的乱码。我重发一下自定义的赝势部分：
I-ECP GEN    28    4
  1 ----- IODINE U(g) -----
0.00000000    2    1.00000000
  3    ----- IODINE U(s)-u(g) -----
49.98964900    2    40.03337600
281.00655600    2    17.30057600
61.41673900    2    8.85172000
  4    ----- IODINE U(p)-u(g)    -----
67.41623900    2    15.72014100
134.80769600    2    15.20822200
14.56654800    2    8.29418600
28.96842200    2    7.75394900
  4    ----- IODINE U(d)-u(g)    -----
35.53875600    2    13.81775100
53.33975900    2    13.58780500
9.71646600    2    6.94763000
14.97750000    2    6.96009900
  4    ----- IODINE U(f)-u(g)    -----
-20.17661800    2    18.52295000
-26.08807700    2    18.25103500
-0.22043400    2    7.55790100
-0.22164600    2    7.59740400

作者
Author: brothers 时间: 2015-4-3 12:31
顶一下，来学习~

静待高手出现

作者
Author: molx79 时间: 2015-4-3 12:51

brothers 发表于 2015-4-3 12:31
顶一下，来学习~ 静待高手出现

昨天晚上看说明书，就知道了把： g-ul potentia转换成----- IODINE U(g) -----
但是还是有问题

作者
Author: brothers 时间: 2015-4-3 14:42

molx79 发表于 2015-4-3 12:51
昨天晚上看说明书，就知道了把： g-ul potentia转换成----- IODINE U(g) -----
但是还是有问题

这一行只是block名字，可以随意设置的看sob的帖子这么说的

作者
Author: molx79 时间: 2015-4-4 14:39

brothers 发表于 2015-4-3 14:42
这一行只是block名字，可以随意设置的看sob的帖子这么说的

已计算成功，因为是两个I，还需要写第一个I的ECP，第二个I写为：I-ECP，仅此而已。不知道我描述清楚没。需有需要，我把输入文件发给您。

作者
Author: brothers 时间: 2015-4-4 21:41

molx79 发表于 2015-4-4 14:39
已计算成功，因为是两个I，还需要写第一个I的ECP，第二个I写为：I-ECP，仅此而已。不知道我描述清楚没。 ...

懂了~我试试

作者
Author: 小靖靖 时间: 2023-4-24 11:20

molx79 发表于 2015-4-4 14:39
已计算成功，因为是两个I，还需要写第一个I的ECP，第二个I写为：I-ECP，仅此而已。不知道我描述清楚没。 ...

您好，我跟您遇到了同样的问题，我有三个I，可以交流一下，您是怎么解决你嘛？谢谢

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