镧系金属配合物单点能报错，对于赝势基组格式不确定

声聚变

参考sob老师对于自定义基组输入的贴文，以下是我的输入代码，不确定赝势基组定义部分是否正确，目前运行报错

%nprocshared=4
%mem=8GB
%chk=D:\G16W\data\121\OR2-stableopt.chk
#p pbe1pbe/def2TZVP em=gd3bj int=fine stable=opt

opted

0 5
C                  3.24302200   -0.32509500    1.38762300
H                 -3.92731800    7.41612900   -1.12252000

Eu    0
S  5  1.00
  70059.420000      0.97000000000E-04
  10776.235000      0.70700000000E-03
  2482.4901000      0.33620000000E-02
  702.15260000      0.10206000000E-01
  216.79260000      0.16597000000E-01
S  1  1.00
  41.495000000       1.0000000000
S  1  1.00
  29.496900000       1.0000000000
S  1  1.00
  15.030900000       1.0000000000
S  1  1.00
  3.7772000000       1.0000000000
S  1  1.00
  1.9219000000       1.0000000000
S  1  1.00
0.71310000000       1.0000000000
S  1  1.00
0.30970000000       1.0000000000
S  1  1.00
0.54900000000E-01   1.0000000000
S  1  1.00
0.22700000000E-01   1.0000000000
P  6  1.00
  1545.4788969     -0.21477687577E-03
  348.79762213     -0.16803101931E-02
  105.52943529     -0.57348847854E-02
  25.617997454     -0.75899902959E-01
  15.784430330      0.24468326999
  3.3852871239     -0.26307923005
P  3  1.00
  7.4940301111     -0.10983416002
  4.5438093966     -0.24311706000
  2.0144742222     -0.18613405999
P  1  1.00
  1.1083886524       1.0000000000
P  1  1.00
0.51142648878       1.0000000000
P  1  1.00
0.22561415288       1.0000000000
P  1  1.00
0.76000181463E-01   1.0000000000
P  1  1.00
0.27117030486E-01   1.0000000000
D  6  1.00
  389.08410000      0.53400000000E-03
  117.76020000      0.42580000000E-02
  44.986900000      0.17160000000E-01
  22.859900000      0.12427000000E-01
  7.8331000000      0.27041600000
  4.1134000000      0.44527200000
D  1  1.00
  2.0999000000       1.0000000000
D  1  1.00
  1.0043000000       1.0000000000
D  1  1.00
0.34220000000       1.0000000000
D  1  1.00
0.10510000000       1.0000000000
F  5  1.00
  128.23410000      0.34240000000E-02
  46.268500000      0.34313000000E-01
  20.736600000      0.12319400000
  9.4479000000      0.25182400000
  4.3149000000      0.35129200000
F  1  1.00
  1.9121000000       1.0000000000
F  1  1.00
0.78970000000       1.0000000000
F  1  1.00
0.28590000000       1.0000000000
G  1  1.00
0.47520000000       1.0000000000

EU     0
EU-ECP     5     28
h POTENTIAL
  1
2      1.00000000            0.00000000
s-h POTENTIAL
  1
2     23.47138400          607.65933100
p-h POTENTIAL
  1
2     16.77247900          264.38547600
d-h POTENTIAL
  1
2     13.98134300          115.38137500
f-h POTENTIAL
  1
2     23.96288800          -49.40079400
g-h POTENTIAL
  1
2     21.23245800          -26.74827300

54687

1. 这叫输入文件，不叫输入代码；
2. 你的坐标信息缺失，只有C和H的。

声聚变

54687 发表于 2025-8-12 22:35
1. 这叫输入文件，不叫输入代码；
2. 你的坐标信息缺失，只有C和H的。

感谢指正，中间的原子坐标我省略了，只保留了第一行和最后一行

Ferrocene

def2TZVP对Eu有定义，不需要自定义基组和赝势；
如果想用自定义基组，基组那里写genecp，文件末尾要对所有元素都写基组定义

声聚变

Ferrocene 发表于 2025-8-12 22:58
def2TZVP对Eu有定义，不需要自定义基组和赝势；
如果想用自定义基组，基组那里写genecp，文件末尾要对所有 ...

sob老师帖文意思不是要单独定义吗，我在坐标后面写的是从turbmole网站下载的赝势定义

sobereva

仔细阅读下文
谈谈赝势基组的选用
http://sobereva.com/373（http://bbs.keinsci.com/thread-5625-1-1.html）

你都没写genecp，当然不会读你自己在坐标后面定义的基组和赝势

声聚变

sobereva 发表于 2025-8-13 00:04
仔细阅读下文
谈谈赝势基组的选用
http://sobereva.com/373（http://bbs.keinsci.com/thread-5625-1-1.ht ...

感谢sob老师解惑，已正常运行