计算化学公社

标题: 过渡金属Ni4团簇上CH4吸附的问题 [打印本页]

作者
Author: 喂大师兄 时间: 2021-7-4 10:31
标题: 过渡金属Ni4团簇上CH4吸附的问题

我用以下泛函和混合基组对CH4在Ni4团簇上的吸附结构进行了优化，单点能我用的b3lyp/def2tzvp计算的，但是吸附能和文献差的有点多，我算下来吸附能是0.55eV,但文献里用同样Ni4团簇算的吸附能为0.42eV,看文献别人通过模拟CH4在Ni(111)面上的吸附情况，吸附能才0.03eV。因为我是初学者，还有很多不明白的地方，希望大家能帮帮我，在这里先谢过大家了

#p b3lyp/genecp opt=calcall int=ultrafine
Title Card Required

0 1
Ni             -0.82138200 1.15673600 -0.73125000
Ni             -0.40944900 -1.30668500 -0.69902800
Ni             -1.05210100 -0.11311400 1.15960300
Ni                1.07621000 0.25771700 0.25289800
C                3.37492200 0.01870300 0.05212800
H                2.83451800 -0.38341000 0.92750000
H                4.01131600 0.83853800 0.37666600
H                3.96943800 -0.79733500 -0.35284700
H                2.72340300 0.37968300 -0.76633300

Ni 0
lanl2dz
****
C H 0
6-31G**
****

Ni 0
lanl2dz

作者
Author: zjxitcc 时间: 2021-7-4 11:18
文献里用同样Ni4团簇算的吸附能为0.42eV -> 你没讲这文献是什么方法、基组算的，计算细节是什么、与你的是否一样。这些因素造成的差别可能很大。

讲道理这种问题DFT算出来一般不靠谱（SCF多重解、多参考特征显著）。不过靠不靠谱与能否重复是两码事，如果你想严格重复文献的值，你就用严格一样的方法、基组、软件。

作者
Author: 喂大师兄 时间: 2021-7-4 12:23

zjxitcc 发表于 2021-7-4 11:18
文献里用同样Ni4团簇算的吸附能为0.42eV -> 你没讲这文献是什么方法、基组算的，计算细节是什么、与你的是 ...

好的，谢谢你的解答，主要文献里面没有给他计算单点能的泛函和基组，所以我就感到很迷茫。听了你的解答，我下去再用其他基组和泛函再算算

作者
Author: Accelerator 时间: 2021-7-4 12:39
本帖最后由 Accelerator 于 2021-7-4 12:41 编辑

这绿幕是怎么回事看得眼睛疼
全是过渡金属别用B3LYP 建议MN15L
就算用B3LYP也得加上D3BJ
Ni4的自旋多重度是多少就算是单重态很可能也不是闭壳层
opt=calcall没什么必要，有这个机时还不如用tzvp优化文献和你所说的“吸附能”分别是电子能量还是自由能
文献如果真的没给方法部分至少表明作者极不专业，没必要非得重复

作者
Author: 喂大师兄 时间: 2021-7-4 12:41

Accelerator 发表于 2021-7-4 12:39
这绿幕是怎么回事看得眼睛疼
全是过渡金属别用B3LYP 建议MN15L
就算用B3LYP也得加上D3BJ

好的，我再用您说的泛函去算算看，谢谢啊

作者
Author: hebrewsnabla 时间: 2021-7-4 12:41
Ni团簇和Ni(111)面的吸附能没有多少可比性。

文献通常不可能不给方法基组，如果真没给那这篇文献没有参考价值。

作者
Author: granvia 时间: 2021-7-4 13:14

喂大师兄发表于 2021-7-4 12:41
好的，我再用您说的泛函去算算看，谢谢啊

如四楼所说，必须尝试多种可能的自旋多重度。过渡金属体系小团簇一般都要考虑自旋，而且吸附前后基态自旋多重度可能发生改变

作者
Author: wzkchem5 时间: 2021-7-4 16:32

喂大师兄发表于 2021-7-4 05:23
好的，谢谢你的解答，主要文献里面没有给他计算单点能的泛函和基组，所以我就感到很迷茫。听了你的解答， ...

文献的SI里也没有给吗？

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