过渡金属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

zjxitcc

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

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

喂大师兄

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

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

Accelerator

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

喂大师兄

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

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

hebrewsnabla

Ni团簇和Ni(111)面的吸附能没有多少可比性。

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

granvia

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

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

wzkchem5

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

文献的SI里也没有给吗？