计算化学公社

标题: 用def2TZVP这样来定义金属cu和非金属碘合理吗？ [打印本页]

作者
Author: shengwenyang 时间: 2017-7-11 08:33
标题: 用def2TZVP这样来定义金属cu和非金属碘合理吗？
最近用高斯算一个体系，在碘未出现前换算过后能量是合理的，但是加入碘之后的能量就变得很诡异，怀疑是碘的基组或者极化系数的问题，求大神解疑
  Cu I 0
  def2tzvp
  ****
  Cu 0
  f 1 1.0
  3.525 1.0
  ****
  I 0
  d 1 1.0
  0.289 1.0
  ****

作者
Author: 我本是个娃娃 时间: 2017-7-11 08:45
def2-TZVP的极化函数已经很合理了，没必要加入更多的极化函数。如果十分想要，莫不如用cc-pVTZ

另外，对于铜和碘原子来说，用赝势基组更加的方便，计算结果也不错

作者
Author: shengwenyang 时间: 2017-7-11 08:54

我本是个娃娃发表于 2017-7-11 08:45
def2-TZVP的极化函数已经很合理了，没必要加入更多的极化函数。如果十分想要，莫不如用cc-pVTZ

另外，对 ...

def2-TZVP对前四周期元素是全电子基组，那意思就是cu和碘I在这里是全电子基组？我试过去掉极化函数来算，结果也是不如人意，可能真如娃娃所说，要换成cc-pVTZ会好些

作者
Author: sobereva 时间: 2017-7-11 08:59
def2tzvp对I和Cu本来就有极化函数，你这样多此一举，还可能造成结果更烂

而且不知道你定义赝势没有，这么写对I只定义了赝势基组而没定义赝势，结果会完全错误

这和改用cc-pVTZ没有丝毫关系
cc-pVTZ对I根本没有定义，除非你用-PP版。

作者
Author: 我本是个娃娃 时间: 2017-7-11 08:59

shengwenyang 发表于 2017-7-11 08:54
def2-TZVP对前四周期元素是全电子基组，那意思就是cu和碘I在这里是全电子基组？我试过去掉极化函数来算， ...

试试赝势基组能更好一点，SDD就不错。

如果非要用全电子基组，大可以用赝势基组下优化好的几何后在全电子基组下做一个单点。

作者
Author: shengwenyang 时间: 2017-7-11 09:09

sobereva 发表于 2017-7-11 08:59
def2tzvp对I和Cu本来就有极化函数，你这样多此一举，还可能造成结果更烂

而且不知道你定义赝势没有，这 ...

谢谢Sob老师，我之前查文献知道已经包含极化函数，一开始没给Cu和I加极化函数，能量出现问题后就加回来试试，这是我优化的#p opt(maxstep=1) B3LYP gen pseudo=read scrf=(smd,solvent=dichloromethane)，单点能#p M06 gen pseudo=read EmpiricalDispersion=GD3 scrf=(smd,read,solvent=dichloromethane)，都已经定义赝势，不知道哪里有误，麻烦Sob老师指点

作者
Author: sobereva 时间: 2017-7-11 09:15

shengwenyang 发表于 2017-7-11 09:09
谢谢Sob老师，我之前查文献知道已经包含极化函数，一开始没给Cu和I加极化函数，能量出现问题后就加回来试 ...

参考我4L说的，把你最终写好的完整输入文件贴上来
最省事就是直接写比如B3LYP/def2TZVP，这时候所有原子都用def2TZVP，该用赝势的也直接用赝势了，看
详解Gaussian中混合基组、自定义基组和赝势基组的输入
http://sobereva.com/60

作者
Author: shengwenyang 时间: 2017-7-11 09:30

sobereva 发表于 2017-7-11 09:15
参考我4L说的，把你最终写好的完整输入文件贴上来
最省事就是直接写比如B3LYP/def2TZVP，这时候所有原 ...

Sob老师，这是我的的单点能输入文件，我现在先试一下用def2TZVP来定义全部原子

作者
Author: sobereva 时间: 2017-7-11 10:17

shengwenyang 发表于 2017-7-11 09:30
Sob老师，这是我的的单点能输入文件，我现在先试一下用def2TZVP来定义全部原子

写ExternalIteration、TABS=253.15毫无意义
当前没给I定义赝势，结果完全错误

作者
Author: shengwenyang 时间: 2017-7-11 11:31

sobereva 发表于 2017-7-11 10:17
写ExternalIteration、TABS=253.15毫无意义
当前没给I定义赝势，结果完全错误

谢谢Sob老师，没考虑到还需要给I定义赝势，找到问题根源很感谢Sob的指导，谢谢。

作者
Author: shengwenyang 时间: 2017-7-11 11:31

我本是个娃娃发表于 2017-7-11 08:45
def2-TZVP的极化函数已经很合理了，没必要加入更多的极化函数。如果十分想要，莫不如用cc-pVTZ

另外，对 ...

谢谢娃娃的耐心解答，谢谢

作者
Author: liyuanhe211 时间: 2017-7-11 12:46

我本是个娃娃发表于 2017-7-11 08:45
def2-TZVP的极化函数已经很合理了，没必要加入更多的极化函数。如果十分想要，莫不如用cc-pVTZ

另外，对 ...

def2-TZVP对I本来就是赝势基组吧。
另外对价层描述接近的赝势基组和全电子基组，在后者不单独考虑相对论的情况下，对周期较高的元素结果是不是还没有前者好？

作者
Author: 冰释之川 时间: 2017-7-11 12:52

liyuanhe211 发表于 2017-7-11 12:46
def2-TZVP对I本来就是赝势基组吧。
另外对价层描述接近的赝势基组和全电子基组，在后者不单独考虑相对论 ...

对，第五周期开始的元素要考虑相对论效应，如果全电子基组的话，需做标量相对论

作者
Author: 我本是个娃娃 时间: 2017-7-11 13:02

liyuanhe211 发表于 2017-7-11 12:46
def2-TZVP对I本来就是赝势基组吧。
另外对价层描述接近的赝势基组和全电子基组，在后者不单独考虑相对论 ...

def2-TZVP对碘是赝势基组

PS：李大博士现在越来越谦虚了

作者
Author: beefly 时间: 2017-7-11 16:28
我记得gaussian09.a等几个早期版本的def2系列基组有问题，重元素遗漏了ecp，结果错误。后来的09.d、09.e版本是正常的。

作者
Author: 虎王 时间: 2020-8-2 14:53

sobereva 发表于 2017-7-11 09:15
参考我4L说的，把你最终写好的完整输入文件贴上来
最省事就是直接写比如B3LYP/def2TZVP，这时候所有原 ...

如果是含金银的体系，也可以直接写B3LYP/DEF2TZVP吗？这样对金银是赝式吗？

作者
Author: thanhtam 时间: 2020-8-2 15:10

虎王发表于 2020-8-2 14:53
如果是含金银的体系，也可以直接写B3LYP/DEF2TZVP吗？这样对金银是赝式吗？

可以，是

作者
Author: 虎王 时间: 2020-8-2 15:22

thanhtam 发表于 2020-8-2 15:10
可以，是

能提交，只是出现这个问题 Out-of-memory error in routine...

作者
Author: thanhtam 时间: 2020-8-2 15:32

虎王发表于 2020-8-2 15:22
能提交，只是出现这个问题 Out-of-memory error in routine...

这是另外的问题。内存不够。

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