计算化学公社

标题: 用Def2-TZVP优化含有Al、O原子的团簇时,O原子也用了Al的基组输入,是否需要重算? [打印本页]

作者
Author:
jiameiye    时间: 2016-6-13 13:37
标题: 用Def2-TZVP优化含有Al、O原子的团簇时,O原子也用了Al的基组输入,是否需要重算?
请问,用Def2-TZVP优化含有Al、O原子的团簇时,O原子也用了Al的基组输入,见下文,这样的结果是否需要重算?O和Al 分开在文件末尾各用各的输入?还是在路径中直接写B3LYP/Def2TZVP?
谢谢!

Al  O   0
S   7   1.00
  37792.5507720              0.57047888709E-03
   5668.0682165              0.44093016538E-02
   1289.8582841              0.22630967411E-01
    364.86596028             0.88025644295E-01
    118.57631515             0.25223701612
     42.024867605            0.45960547169
     15.499501629            0.33277886014
S   3   1.00
     75.208026598            0.19250560190E-01
     23.031408972            0.87906743952E-01
      3.6348797649          -0.34246704535
S   2   1.00
      1.6065049957           1.5106266058
      0.76103394581          0.58071016470
S   1   1.00
      0.16556708849          1.0000000
S   1   1.00
      0.60041577113E-01            1.0000000
P   5   1.00
    452.52303192             0.23110812466E-02
    107.08195049             0.18568641823E-01
     34.131021255            0.87216237035E-01
     12.587037428            0.26902101523
      4.9811919704           0.52128324272
P   1   1.00
      2.0070350900           1.0000000
P   1   1.00
      0.80083714374          1.0000000
P   1   1.00
      0.20178927472          1.0000000
P   1   1.00
      0.57895550392E-01            1.0000000
D   2   1.00
      1.57000000             0.20000000
      0.333000000            1.00000000
D   1   1.00
      0.111000000            1.0000000
F   1   1.00
      0.244000000            1.0000000
****









作者
Author:
liyuanhe211    时间: 2016-6-13 13:38
Gaussian或ORCA都自带def2tzvp自定义它干啥
O用Al的基组的结果显然不能用,不过我倒是好奇结果长什么样。
作者
Author:
jiameiye    时间: 2016-6-13 13:49
liyuanhe211 发表于 2016-6-13 13:38
Gaussian或ORCA都自带def2tzvp自定义它干啥
O用Al的基组的结果显然不能用,不过我倒是好奇结果长什么样。

因为还有别的重金属原子用赝势,写输入文件时候一下就给把O和Al合并写了!
上述这样情况需要自定义吗?然后O和Al各写自己的全电子Def2TZVP,重金属分别写基组部分和赝势部分?谢谢
作者
Author:
jiameiye    时间: 2016-6-13 13:54
liyuanhe211 发表于 2016-6-13 13:38
Gaussian或ORCA都自带def2tzvp自定义它干啥
O用Al的基组的结果显然不能用,不过我倒是好奇结果长什么样。

结果?正常结束,几何结构看上去也是对的。没分析结果文本信息呢。
作者
Author:
sobereva    时间: 2016-6-13 13:57
jiameiye 发表于 2016-6-13 13:49
因为还有别的重金属原子用赝势,写输入文件时候一下就给把O和Al合并写了!
上述这样情况需要自定义吗? ...


你直接写genecp,后面写比如
Al O
def2TZVP
****
Fe
SDD
****

Fe
SDD

就完了。仔细看此文
详解Gaussian中混合基组、自定义基组和赝势基组的输入
http://sobereva.com/60

图省事的话索性全都用def2,即B3LYP/def2TZVP。不过对于过渡金属def2TZVP比SDD耗时高不少。

程序能算不报错、表面上看起来正确不代表真的定量正确。给O用Al的基组显然结果不准。

作者
Author:
jiameiye    时间: 2016-6-13 16:43
sobereva 发表于 2016-6-13 13:57
你直接写genecp,后面写比如
Al O
def2TZVP

好的,谢谢社长耐心解答回复,我去看帖子。
作者
Author:
jiameiye    时间: 2016-6-13 16:54
sobereva 发表于 2016-6-13 13:57
你直接写genecp,后面写比如
Al O
def2TZVP

社长好,写成下列这样可以吗?

Al O Fe 0
def2TZVP
****

Fe 0
def2TZVP


作者
Author:
sobereva    时间: 2016-6-13 17:24
jiameiye 发表于 2016-6-13 16:54
社长好,写成下列这样可以吗?

Al O Fe 0

def2你别这么写,否则程序可能没法正确判断
作者
Author:
jiameiye    时间: 2016-6-13 19:15
sobereva 发表于 2016-6-13 17:24
def2你别这么写,否则程序可能没法正确判断

社长,那这个def2TZVP该怎么样写?谢谢,又来麻烦您了,我还是有些困惑,

我现在是这么写的:
Al O 0
Def2TZVP
****
Fe 0
EMSL上Def2TZVP的基组自定义片段
****

Fe 0
EMSL上Def2TZVP的赝势自定义片段

可不可以简化成下列格式?
Al O 0
Def2TZVP
****
Fe 0
Def2TZVP
****

Fe 0
Def2TZVP
作者
Author:
sobereva    时间: 2016-6-13 20:03
jiameiye 发表于 2016-6-13 19:15
社长,那这个def2TZVP该怎么样写?谢谢,又来麻烦您了,我还是有些困惑,

我现在是这么写的:


前面已经回复过了:“图省事的话索性全都用def2,即B3LYP/def2TZVP。不过对于过渡金属def2TZVP比SDD耗时高不少。”
只写B3LYP/def2TZVP就完事了

作者
Author:
jiameiye    时间: 2016-6-13 20:16
sobereva 发表于 2016-6-13 17:24
def2你别这么写,否则程序可能没法正确判断

社长好,我是想所有原子(Ag,Au,Pd,Pt,O,Al)都用def2TZVP,不使用混合基组,
这样写#B3LYP/def2TZVP,是否可以实现Al 、O用全电子,其它重原子用赝势?
作者
Author:
sobereva    时间: 2016-6-13 23:11
jiameiye 发表于 2016-6-13 20:16
社长好,我是想所有原子(Ag,Au,Pd,Pt,O,Al)都用def2TZVP,不使用混合基组,
这样写#B3LYP/def2TZVP,是 ...

def2对于前四周期都是全电子基组,对于后头的都是赝势基组
作者
Author:
jiameiye    时间: 2016-6-14 10:17
sobereva 发表于 2016-6-13 23:11
def2对于前四周期都是全电子基组,对于后头的都是赝势基组

好的,明白了,谢谢社长。




欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3