计算化学公社

标题: 求助SDD赝势基组作者官网上基组的含义和选择 [打印本页]

作者
Author: mt13 时间: 2021-12-7 10:10
标题: 求助SDD赝势基组作者官网上基组的含义和选择
本帖最后由 mt13 于 2021-12-7 10:25 编辑

请教各位老师，在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
。。。。。。

作者
Author: beefly 时间: 2021-12-7 18:03
gaussian内置的sdd基组比较老，大概是Basis Set Exchange里的Stuttgart RLC或者Stuttgart RSC 1997，不包括cc-pvnz-pp型的基组

作者
Author: mt13 时间: 2021-12-7 22:15

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

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

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

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

3 是

4 打印出来对比便知

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

作者
Author: mt13 时间: 2021-12-8 23:21

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

谢谢老师

作者
Author: mt13 时间: 2021-12-8 23:22

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

谢谢老师

作者
Author: mt13 时间: 2021-12-9 14:09
本帖最后由 mt13 于 2021-12-9 22:25 编辑

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)吗？

作者
Author: mt13 时间: 2021-12-10 20:47

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

麻烦sob老师解答一下

作者
Author: sobereva 时间: 2021-12-11 04:34

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

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

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

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

作者
Author: mt13 时间: 2021-12-11 11:08

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

2 可以。也 ...

谢谢老师

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