计算化学公社

标题: 镧系金属配合物单点能报错，对于赝势基组格式不确定 [打印本页]

作者
Author: 声聚变 时间: yesterday 21:54
标题: 镧系金属配合物单点能报错，对于赝势基组格式不确定
参考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

作者
Author: 54687 时间: yesterday 22:35
本帖最后由 54687 于 2025-8-12 22:39 编辑

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

作者
Author: 声聚变 时间: yesterday 22:47

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

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

作者
Author: Ferrocene 时间: yesterday 22:58
def2TZVP对Eu有定义，不需要自定义基组和赝势；
如果想用自定义基组，基组那里写genecp，文件末尾要对所有元素都写基组定义

作者
Author: 声聚变 时间: yesterday 23:20
本帖最后由声聚变于 2025-8-12 23:23 编辑

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

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

作者
Author: sobereva 时间: 12 hour ago
仔细阅读下文
谈谈赝势基组的选用
http://sobereva.com/373（http://bbs.keinsci.com/thread-5625-1-1.html）

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

作者
Author: 声聚变 时间: 2 hour ago

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

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

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