计算化学公社

标题: 在更高级别基组下的单点能比低级别下优化时的单点能高的多，请高手帮忙看一下 [打印本页]

作者
Author: fuzanhui 时间: 2023-10-22 19:22
标题: 在更高级别基组下的单点能比低级别下优化时的单点能高的多，请高手帮忙看一下
(, 下载次数 Times of downloads: 2)
各位大佬，我在使用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

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

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

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

作者
Author: fuzanhui 时间: 2023-10-23 09:00

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

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

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

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