计算化学公社

标题: gaussian09 用Def2TZVPP计算过渡金属 [打印本页]

作者
Author: 高雪皎 时间: 2014-11-12 09:42
标题: gaussian09 用Def2TZVPP计算过渡金属
亲爱的社员们：
我看到一篇文献，用GAMESS程序优化了含有镧系金属的体系，用的是杂化泛函B3lyp，对于轻原子用6-311g*,镧系金属La用了赝势，def2yzvpp基组。我想用相同的方法和基组用gaussian09重复一下计算，但是总是报错，说金属超出了def2yzvpp的计算范围，但是我在输出文件中并没有找到金属原子赝势的部分，所以我觉得可能是我的输入文件写的有问题，昨天我查了好久，高斯手册上关于def2tzvpp只有很少的一段话，并没有给出例子。社员们有用过这个基组的吗，能不能告诉我哪里出错了，非常谢谢大家。我把文献（参看计算方法部分，已高亮显示）和输入输出文件附上。

作者
Author: sobereva 时间: 2014-11-12 10:13
def和def2系列基组都是没有镧系的，那篇文章写得有问题。推荐用SDD

作者
Author: jiangning198511 时间: 2014-11-12 10:17

sobereva 发表于 2014-11-12 10:13
def和def2系列基组都是没有镧系的，那篇文章写得有问题。推荐用SDD

def2-TZVPP Basis set for La in Gaussian-format

****
La 0
S  2  1.00
  5.0873990000    -0.44174952534
  4.2709780000    0.85812466843
S  1  1.00
  2.1823656045    1.0000000000
S  1  1.00
0.48966533785    1.0000000000
S  1  1.00
0.23301675348    1.0000000000
S  1  1.00
0.55719523051E-01 1.0000000000
S  1  1.00
0.22854708808E-01 1.0000000000
P  4  1.00
  6.0000000000    -0.11397960924E-01
  3.6808191615    0.14675038550
  2.3265462081    -0.35581819167
0.64342629633    0.45834955229
P  1  1.00
0.33584281962    1.0000000000
P  1  1.00
0.16519051916    1.0000000000
P  1  1.00
0.354E-01          1.0000000000
D  3  1.00
  1.2675288018    -0.17569274035
0.89395340284    0.25139922933
0.33095767339    0.44603267053
D  1  1.00
0.13461377501    1.0000000000
D  1  1.00
0.51441629839E-01 1.0000000000
F  1  1.00
0.38941 1.0
F  1  1.00
0.12043 1.0
G  1  1.00
0.17907 1.0
! F. Weigend, R. Ahlrichs; Phys. Chem. Chem. Phys. 7, 3297 (2005).
****

LA    0
LA-ECP    3    46
f POTENTIAL
  1
2    4.02860000       -36.01001600
s-f POTENTIAL
  3
2    3.30990000          91.93217700
2    1.65500000          -3.78876400
2    4.02860000          36.01001600
p-f POTENTIAL
  3
2    2.83680000          63.75948600
2    1.41840000          -0.64795800
2    4.02860000          36.01001600
d-f POTENTIAL
  3
2    2.02130000          36.11617300
2    1.01070000          0.21911400
2    4.02860000          36.01001600
! M. Dolg, H. Stoll, A. Savin, H. Preuss; Theor. Chim. Acta 75, 173 (1989).
! M. Dolg, H. Stoll, H. Preuss; Theor. Chim. Acta 85, 441 (1993).

作者
Author: jiangning198511 时间: 2014-11-12 10:23
这个基组是在TURBOMOLE上找到的，没有收录在 EMSL-BASIS
http://www.cosmologic.de/basis-sets/basissets.php

作者
Author: jiangning198511 时间: 2014-11-12 12:27
这系列的基组只有La有，其他镧系元素没有def2,只有def

作者
Author: truemantwo 时间: 2014-11-12 14:27

jiangning198511 发表于 2014-11-12 10:23
这个基组是在TURBOMOLE上找到的，没有收录在 EMSL-BASIS
http://www.cosmologic.de/basis-sets/basissets. ...

呃，这个页面报错：

Page not Found (404)

作者
Author: chrinide 时间: 2014-11-12 15:49
本帖最后由 chrinide 于 2014-11-12 15:54 编辑

检查了下TB的basen文件下 la 的基组确实有def2系列的基组 EMSL上的基组没有更新。
TB的基组下载地址改了： http://www.cosmologic-services.de/basis-sets/basissets.php

#  available basis sets and corresponding atomic La(2D) energies:
#
#    Basis                UHF(noneq) / a.u.    ROHF(equiv) / a.u.
# ------------------------------------------------------------------
# ecp-46-mwb-SVP          -31.08076818692    -31.07033330186
# def-SVP                         -31.08076818692    -31.07033330186
# def-SV(P)                      -31.08076818692    -31.07033330186
# def2-SV(P)                   -31.08076818692    -31.07033330186
# def2-SVP                      -31.08342531180    -31.07033330186
# ecp-46-mwb-TZVP       -31.08122307866    -31.07033330186
# def-TZVP                      -31.08122307866    -31.07033330186
# ecp-46-mwb-TZVext    -31.08139200128    -31.07040696627
# ecp-46-mwb-QZV          -31.08164671966    -31.07058686661
# ecp-46-mwb-TZVPP       -31.08298007337    -31.07033330186
# def-TZVPP                   -31.08298007337    -31.07033330186
# ecp-46-mwb-TZVPext    -31.08405204339    -31.07040696627
# def2-TZVP                   -31.08405204339    -31.07040696627
# ecp-46-mwb-TZVPPext -31.08496756403    -31.07040696627
# def2-TZVPP                   -31.08496756403    -31.07040696627
# ecp-46-mwb-QZVP       -31.08539566442    -31.07058686661
# def-QZVP                      -31.08539566442    -31.07058686661
# def2-QZVP                      -31.08539566442    -31.07058686661
# ecp-46-mwb-QZVPP       -31.08553371306    -31.07058686661
# def-QZVPP                      -31.08553371306    -31.07058686661
# def2-QZVPP                   -31.08553371306    -31.07058686661
# ecp-46-mwb-QZVPPh    -31.08557507678    -31.07058686661
# ecp-46-mwb-QZV5f4g3h -31.08576635399    -31.07058686661
# Hay-Wadt MB (n+1) ECP -30.76610762512    -30.75558910789
# Hay-Wadt VDZ (n+1) ECP -30.76610762512    -30.75558910789
# LANL2DZ ECP                -30.76610762512    -30.75558910789
# ------------------------------------------------------------------

作者
Author: 高雪皎 时间: 2014-11-12 17:39

jiangning198511 发表于 2014-11-12 12:27
这系列的基组只有La有，其他镧系元素没有def2,只有def

哦，原来是这样，只有La的有，没有其他镧系金属的，那就不能用这个方法算了，谢写您啦！

作者
Author: 高雪皎 时间: 2014-11-12 17:50

chrinide 发表于 2014-11-12 15:49
检查了下TB的basen文件下 la 的基组确实有def2系列的基组 EMSL上的基组没有更新。
TB的基组下载地址改了 ...

我在网上下载了Gd的def2TZVPP基组，是不是要用extrabasis关键字来把这个基组写到输入文件中？我用的是混合基组，只有金属用这个基组，输入文件该怎么写呢，您可以告诉一下我吗，输入文件在附件里？谢谢啦

作者
Author: sobereva 时间: 2014-11-13 00:45

高雪皎发表于 2014-11-12 17:50
我在网上下载了Gd的def2TZVPP基组，是不是要用extrabasis关键字来把这个基组写到输入文件中？我用的是混 ...

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

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