关于比较结合能时切换泛函的问题

johngu111

各位老师好，想请教一下如果我想比较一个主族元素（例如Al离子）和镧系元素（例如La离子）与同一阴离子间的结合能差异。1.优化结构的时候对Al用了B3LYP泛函，对La用了PBE0泛函，那之后算单点都用PBE0泛函得出的结论合理吗？2.如果算单点用对Al用M06泛函，对La用PBE0泛函，这样是不是不行？3.优化结构的时候用不同泛函合理吗，还是也必须统一泛函？

wzkchem5

必须用同一个泛函。如果两个配合物的最优泛函不同，必须进行折衷，挑一个对两个配合物的误差都不是特别大的泛函。
例如这里统一用PBE0就行了，Al用PBE0虽没有特别的优势，但也并没有劣势

johngu111

wzkchem5 发表于 2024-4-9 02:26
必须用同一个泛函。如果两个配合物的最优泛函不同，必须进行折衷，挑一个对两个配合物的误差都不是特别大的 ...

谢谢老师，还有如果我想比较不同金属离子的软硬酸碱程度，例如定量比较Al离子和La离子哪个硬酸属性更强，用什么参数比较合适呢，静电势？全局电子硬度？LUMO轨道能量？

sobereva

johngu111 发表于 2024-4-9 07:46
谢谢老师，还有如果我想比较不同金属离子的软硬酸碱程度，例如定量比较Al离子和La离子哪个硬酸属性更强， ...

电子硬度或者极化率

johngu111

sobereva 发表于 2024-4-9 11:40
电子硬度或者极化率

谢谢sob老师，电子硬度是可以取Multiwfn生成的CDFT参数里面的Hardness吗？

sobereva

johngu111 发表于 2024-4-9 12:22
谢谢sob老师，电子硬度是可以取Multiwfn生成的CDFT参数里面的Hardness吗？

是

johngu111

sobereva 发表于 2024-4-10 03:10
是

老师我用multiwfn尝试计算电子硬度的时候出现一个问题，是计算四价Sn离子的时候，无法输出正常的N-1波函数文件，Gaussian报错信息如下：
20 basis functions,    21 primitive gaussians,    21 cartesian basis functions
     0 alpha electrons       -1 beta electrons
       nuclear repulsion energy         0.0000000000 Hartrees.
Bad numbers of electrons.
然后我用Multiwfn单独跑了一下Sn4+的HOMO轨道能量发现是0，为什么会这样呢？那我如何能计算这个离子的电子硬度呢？附件里是Gaussian计算的N-1输入和输出文件，还有Sn4+的fch文件，输入文件中用的是拷贝的LANL08(d)基组。请问老师这个计算哪里有问题？

paramecium86

johngu111 发表于 2024-4-10 10:08
老师我用multiwfn尝试计算电子硬度的时候出现一个问题，是计算四价Sn离子的时候，无法输出正常的N-1波函 ...

Sn 一共50个电子， 输入文件用了ECP46 赝势 也就是说Sn只考虑了4个电子。而N-1是+5的charge 自然就会出现-1个 beta电子的情况。这种情况换ECP28赝势就可以了。

johngu111

paramecium86 发表于 2024-4-10 10:35
Sn 一共50个电子， 输入文件用了ECP46 赝势 也就是说Sn只考虑了4个电子。而N-1是+5的charge 自然就会出现 ...

谢谢，明白了，我再去试试

johngu111

sobereva 发表于 2024-4-10 03:10
是

老师我还想到一个问题，就比如我想对比Al、Fe、Sn、La、Ce离子等多种不同族、周期也差的比较远的元素的极化率或者电子硬度，靠前的Al这种一般都是全电子基组，而Sn、La、Ce等靠后的元素是赝势基组，这样是否会对结果产生较大的影响，如果想获得比较好一点的结果该用什么精度的基组比较好呢？def2-SVP这种是否可行？

sobereva

johngu111 发表于 2024-4-10 15:11
老师我还想到一个问题，就比如我想对比Al、Fe、Sn、La、Ce离子等多种不同族、周期也差的比较远的元素的极 ...

没任何不良影响
def2-SVP算极化率是开玩笑。至少def2-SVPD
仔细看
谈谈量子化学中基组的选择
http://sobereva.com/336（http://bbs.keinsci.com/thread-3545-1-1.html）