计算化学公社

标题: 用B2PLYP方法算3d后过渡金属的键能，总出现L9999或者L508错误，请教前辈们是方法的... [打印本页]

作者
Author: Scx 时间: 2018-3-19 20:07
标题: 用B2PLYP方法算3d后过渡金属的键能，总出现L9999或者L508错误，请教前辈们是方法的...
没有加SCF=QC的时候，是出现L9999错误，试了几种改正方法没有效果，加了SCF=QC，又变成了L508错误。计算3d后过渡金属键能，用B3LYP都试了，都没有问题。想着再用B2LYP算一遍，对比一下，但总是出现不收敛错误。

看到有前辈的博客写了这样一句话：（B2PLYP：MP2与B2LYP(B88+HF+LYP)能量相混合，比DFT慢一个数量级，速度在MP2级别。弱相互作用虽然比一般DFT都强，但也不算太好，不如DFT-D，有的稀有气体二聚体预测成非束缚态。算普通体系也经常比某些其它泛函差。）所以想问问是不是B2LYP方法本身不适合的原因还是我的输入有错误？谢谢前辈！

输入文件：
%chk=CoN-M5-B2PLYP.chk
%mem=2GB
%nproc=8
# B2PLYP/gen freq opt(maxcycle=80) IOp(3/125=0360003600,3/78=0640006400,3/76=0350006500,3/77=1000010000,5/33=1)
# empiricaldispersion=gd3bj int(grid=ultrafine) SCF=QC

Title Card Required

0 5
Co             -0.01370958 0.00000000 0.00000000
N                -1.24370958 0.00000000 0.00000000

N  0
6-31G*
****
Co 0
Def2QZVP
****

作者
Author: sobereva 时间: 2018-3-19 21:07
算个这个有什么必要用def2-QZVP，真要是要高精度，N就不可能就用个6-31G*
仔细看
谈谈量子化学中基组的选择
http://sobereva.com/336

508、9999都是因为什么，仔细看，以后问问题别问Lxxx出错怎么办
Gaussian FAQ
http://bbs.keinsci.com/forum.php?mod=viewthread&tid=4829

初学者别写一堆乱七八糟IOp

作者
Author: Scx 时间: 2018-3-19 21:35

sobereva 发表于 2018-3-19 21:07
算个这个有什么必要用def2-QZVP，真要是要高精度，N就不可能就用个6-31G*
仔细看
谈谈量子化学中基组的选 ...

好的，知道了

作者
Author: Scx 时间: 2018-3-19 21:41

Scx 发表于 2018-3-19 21:35
好的，知道了

谢谢

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