求助SDD赝势基组作者官网上基组的含义和选择

mt13

请教各位老师，在SDD赝势基组作者官网http://www.tc.uni-koeln.de/PP/clickpse.en.html上，碘的ECP28MDF有多种基组，比如ECP28MDF_VDZ (VDZ-PP cc-pVDZ-PP)、ECP28MDF_AVDZ (AVDZ-PP aug-cc-pVDZ-PP)等等。
1、如何知道Gaussian中SDD赝势搭配的是什么赝势基组呢，在Gaussian用户指南上好像只能看到某种元素是MDF或者MWB，看不到默认的赝势基组；
2、ECP28MDF_VDZ基组大小是不是和6-31G**差不多，ECP28MDF_VTZ是不是和6-311G**差不多；
3、ECP28MDF_AVDZ (AVDZ-PP aug-cc-pVDZ-PP)这个基组的意思，是带着弥散的基组吧；
4、ECP28MDF_AVDZ (AVDZ-PP aug-cc-pVDZ-PP)是不是和https://www.basissetexchange.org/上的aug-cc-pVDZ-PP完全一样？
    ECP28MDF_AVDZ：
I 0
S 8 1.00
0.24497900E+04 0.41900000E-03
0.35980800E+03 0.22400000E-02
0.14405800E+02 0.39722300E+00
0.90763200E+01 -0.93224900E+00
0.20881000E+01 0.93713800E+00
0.10349800E+01 0.39208600E+00
0.31628400E+00 0.12485000E-01
0.12171900E+00 -0.13290000E-02
S 8 1.00
0.24497900E+04 0.17500000E-03
0.35980800E+03 0.10570000E-02
0.14405800E+02 0.16900000E+00
0.90763200E+01 -0.42179300E+00
0.20881000E+01 0.63886400E+00
0.10349800E+01 0.32011500E+00
0.31628400E+00 -0.81442800E+00
0.12171900E+00 -0.48979800E+00
S 1 1.00
0.31628400E+00 1.0
S 1 1.00
0.12171900E+00 1.0
。。。。。。


    aug-cc-pVDZ-PP：
I     0
S    8   1.00
      2.449790D+03           4.190000D-04
      3.598080D+02           2.240000D-03
      1.440580D+01           3.972230D-01
      9.076320D+00          -9.322490D-01
      2.088100D+00           9.371380D-01
      1.034980D+00           3.920860D-01
      3.162840D-01           1.248500D-02
      1.217190D-01          -1.329000D-03
S    1   1.00
      3.162840D-01           1.000000D+00
S    8   1.00
      2.449790D+03           1.750000D-04
      3.598080D+02           1.057000D-03
      1.440580D+01           1.690000D-01
      9.076320D+00          -4.217930D-01
      2.088100D+00           6.388640D-01
      1.034980D+00           3.201150D-01
      3.162840D-01          -8.144280D-01
      1.217190D-01          -4.897980D-01
S    1   1.00
      1.217190D-01           1.000000D+00
S    1   1.00
      4.200000D-02           1.0000000
。。。。。。

beefly

gaussian内置的sdd基组比较老，大概是Basis Set Exchange里的Stuttgart RLC或者Stuttgart RSC 1997，不包括cc-pvnz-pp型的基组

mt13

beefly 发表于 2021-12-7 18:03
gaussian内置的sdd基组比较老，大概是Basis Set Exchange里的Stuttgart RLC或者Stuttgart RSC 1997，不包括 ...

谢谢老师，虽然比较老，用起来也没问题吧？

sobereva

1 搭配的是赝势提出时标配的Dolg搞的赝势基组。aug-cc-pVnZ-PP系列是其他人搞的基于Stuttgart赝势的赝势基组。

2 把基组定义打印出来对比。

3 是

4 打印出来对比便知

coolrainbow

使用SDD基组时，建议直接拷贝用GENECP读取。Gaussian自带的无法确定是什么年代的东西了，除非打印出来看

mt13

sobereva 发表于 2021-12-8 08:07
1 搭配的是赝势提出时标配的Dolg搞的赝势基组。aug-cc-pVnZ-PP系列是其他人搞的基于Stuttgart赝势的赝势基 ...

谢谢老师

mt13

coolrainbow 发表于 2021-12-8 09:14
使用SDD基组时，建议直接拷贝用GENECP读取。Gaussian自带的无法确定是什么年代的东西了，除非打印出来看

谢谢老师

mt13

1、再请教各位老师，碘元素ECP46MWB赝势，BSE网站和SDD基组作者官网上赝势基组和赝势定义完全一样。但是在Gaussian16 A.03中，我用IOp(3/24=1)或者GFInput打印基组，发现赝势定义和BSE、SDD网站上的一样，但是赝势基组部分的一些系数不一样，甚至IOp(3/24=1)和GFInput显示的赝势基组的一些系数也不同。是不是G16的赝势基组有问题？是不是更稳妥的做法，是像楼上老师说的使用赝势基组时，建议直接拷贝用GENECP读取，就不用G16自带的了？

BSE、SDD网站上:
I 0
S 3 1.00
2.1227650 1.1040280
1.7704810 -1.5353260
0.3130840 0.751603
S 1 1.00
0.1240710 1.0
P 3 1.00
2.4328870 0.4422320
2.1372490 -0.5838090
0.3145460 0.6266060
P 1 1.00
0.1049450 1.0
P 1 1.00
0.0326410 1.0

GFInput显示的：
S   3 1.00       0.000000000000
      0.2122765000D+01  0.2063429899D+01
      0.1770481000D+01 -0.2869526473D+01
      0.3130840000D+00  0.1404747074D+01
S   1 1.00       0.000000000000
      0.1240710000D+00  0.1000000000D+01
P   3 1.00       0.000000000000
      0.2432887000D+01  0.7732753148D+00
      0.2137249000D+01 -0.1020833156D+01
      0.3145460000D+00  0.1095666871D+01
P   1 1.00       0.000000000000
      0.1049450000D+00  0.1000000000D+01
P   1 1.00       0.000000000000
      0.3264100000D-01  0.1000000000D+01

IOp(3/24=1)显示的：
*     I      0.00000  0.00000  3.78941                                                                                            *
+                                               44   S           1.00                                                               
+                                                                      0.212276D+010.258628D+010.000000D+000.000000D+000.000000D+00
*                                                                     0.177048D+01-.313898D+010.000000D+000.000000D+000.000000D+00*
*                                                                     0.313084D+000.419038D+000.000000D+000.000000D+000.000000D+00*
*                                              45   S           1.00                                                              *
+                                                                      0.124071D+000.148992D+000.000000D+000.000000D+000.000000D+00
*                                          46- 48   P           1.00                                                              *
+                                                                      0.243289D+010.000000D+000.334909D+010.000000D+000.000000D+00
*                                                                     0.213725D+010.000000D+00-.376022D+010.000000D+000.000000D+00*
*                                                                     0.314546D+000.000000D+000.367895D+000.000000D+000.000000D+00*
*                                          49- 51   P           1.00                                                              *
+                                                                      0.104945D+000.000000D+000.851417D-010.000000D+000.000000D+00
*                                          52- 54   P           1.00                                                              *
+                                                                      0.326410D-010.000000D+000.197763D-010.000000D+000.000000D+00

2、我计算的是Cr配合物，配体有I。对于Cr，我用的是G16自带的SDD；对于碘，我不想用G16自带的SDD，因为是大核MWB赝势，这时可以对I用lanl08(d)吗（同时使用SDD和lanl08(d)是否可以）？
3、对于过渡金属，Lanl08(f)和Lanl08不如lanl2TZ、Lanl2TZ(f)吗？

mt13

sobereva 发表于 2021-12-8 08:07
1 搭配的是赝势提出时标配的Dolg搞的赝势基组。aug-cc-pVnZ-PP系列是其他人搞的基于Stuttgart赝势的赝势基 ...

麻烦sob老师解答一下

sobereva

mt13 发表于 2021-12-9 14:09
1、再请教各位老师，碘元素ECP46MWB赝势，BSE网站和SDD基组作者官网上赝势基组和赝势定义完全一样。但是在G ...

1 收缩系数有的地方给的是归一化后的，有的是没归一化的，这无所谓，实际计算时都会归一化

2 可以。也可以用def2-SVP、def2-TZVP

3 lanl08比lanl2TZ好也更大，lanl08(f)比lanl2TZ(f)好也更大

mt13

sobereva 发表于 2021-12-11 04:34
1 收缩系数有的地方给的是归一化后的，有的是没归一化的，这无所谓，实际计算时都会归一化

2 可以。也 ...

谢谢老师