在更高级别基组下的单点能比低级别下优化时的单点能高的多，请高手帮忙看一下

fuzanhui

LaOF.rar (63.32 KB, 下载次数 Times of downloads: 2)

2023-10-22 19:08 上传 Uploaded

点击下载
Click to download
优化及单点能的输入输出文件

各位大佬，我在使用pbe1pbe方法对LaOF分子进行优化时，O、F原子使用6-31G*基组，La原子使用SDD基，优化收敛，HF=-610.455057.

然后使用更高级别的def2QZVPP基组进行单点能计算，输入文件如下： 计算得到的单点能却为：HF=-206.5767804，这相差也太悬殊了，不知道哪里出了问题，请大家帮忙看一下，输入输出文件在附件中，不胜感激。

%chk=LaOF_energy.chk
#p pbe1pbe/genecp em=gd3bj

Title Card Required

0 3
F                 -1.63926000   -0.94412100    0.00000000
O                  1.84416800   -0.82295300    0.00000000
La                 0.00000000    0.26457400    0.00000000

F O 0
def2qzvpp
****
La     0
S    3   1.00
     26.077272000            0.23198684846D-02
      5.7127838000          -0.29787719611
      4.4057848000           0.70022722515
S    1   1.00
      2.1315325829           1.0000000
S    1   1.00
      0.49385808956          1.0000000
S    1   1.00
      0.23595687888          1.0000000
S    1   1.00
      0.63593489661D-01      1.0000000
S    1   1.00
      0.31545394754D-01      1.0000000
S    1   1.00
      0.15777655296D-01      1.0000000
P    4   1.00
      6.0000000000          -0.11397960924D-01
      3.6808191615           0.14675038550
      2.3265462081          -0.35581819167
      0.64342629633          0.45834955229
P    1   1.00
      0.33584281962          1.0000000
P    1   1.00
      0.16519051916          1.0000000
P    1   1.00
      0.56715971797D-01      1.0000000
P    1   1.00
      0.22100824776D-01      1.0000000
D    3   1.00
      1.2000000000          -0.35803482414
      1.0239524401           0.40483734122
      0.37508490282          0.37649639849
D    1   1.00
      0.17161940245          1.0000000
D    1   1.00
      0.78780423402D-01      1.0000000
D    1   1.00
      0.35000047794D-01      1.0000000
F    1   1.00
      0.6558000              1.0000000
F    1   1.00
      0.2568300              1.0000000
F    1   1.00
      0.1005800              1.0000000
G    1   1.00
      0.3321800              1.0000000
G    1   1.00
      0.1353800              1.0000000
****

LA     0
LA-ECP     3     46
f potential
  1
2      4.0286000            -36.0100160
s-f potential
  3
2      3.3099000             91.9321770
2      1.6550000             -3.7887640
2      4.0286000             36.0100160
p-f potential
  3
2      2.8368000             63.7594860
2      1.4184000             -0.6479580
2      4.0286000             36.0100160
d-f potential
  3
2      2.0213000             36.1161730
2      1.0107000              0.2191140
2      4.0286000             36.0100160

sobereva

不同赝势赝化的内核电子数不一样，存在这种差异时，总能量不在一个数量级很正常

并且注意SDD对同一种元素有不同核尺寸的版本，同为大核赝势时还可能有针对不同价态的版本，不能简单说个SDD，参看
谈谈赝势基组的选用
http://sobereva.com/373（http://bbs.keinsci.com/thread-5625-1-1.html）

详解Gaussian中混合基组、自定义基组和赝势基组的输入
http://sobereva.com/60

fuzanhui

sobereva 发表于 2023-10-23 00:20
不同赝势赝化的内核电子数不一样，存在这种差异时，总能量不在一个数量级很正常

并且注意SDD对同一种元 ...

谢谢Sob老师的回复，我先学习一下看看，谢谢