求助高斯过渡金属络合物反应的吉布斯自由能问题

张博

各位老师好，我在计算一个过渡金属络合物（Ru基）和N2分子的反应吉布斯自由能变化值，过渡金属络合物我先用PBE0-D3BJ泛函和SDD（过渡金属）、6-31G*（其他原子）基组做几何优化和振动分析，随后用PBE0-D3BJ泛函和def2-TZVP基组做单点能计算，这样合理么？另外，为计算反应的吉布斯自由能变化值，请问N2分子我应该用过渡金属几何优化的基组还是单点能的基组来计算单点能呢？

18217265596

涉及反应的原子6-31G*实在是太低了。

张博

18217265596 发表于 2025-6-9 10:00
涉及反应的原子6-31G*实在是太低了。

我是想着前面几何优化先用低一点的基组，后面做单点能计算再提高基组，这样有问题么

ljh123

张博 发表于 2025-6-9 10:34
我是想着前面几何优化先用低一点的基组，后面做单点能计算再提高基组，这样有问题么

那也不能太低，改用6-311g**稳健些

张博

ljh123 发表于 2025-6-9 10:36
那也不能太低，改用6-311g**稳健些

好的感谢，那请问反应的吉布斯自由能计算，其中其他的简单小分子自由能计算要用前面的低基组还是用后面的高基组呢，我记得是要统一基组才能做反应自由能变化的计算。

sai77

张博 发表于 2025-6-9 10:41
好的感谢，那请问反应的吉布斯自由能计算，其中其他的简单小分子自由能计算要用前面的低基组还是用后面的 ...

那当然要统一基组
其他不含金属的小分子优化时在相同基组下加上5D关键词就可以了

18217265596

张博 发表于 2025-6-9 10:34
我是想着前面几何优化先用低一点的基组，后面做单点能计算再提高基组，这样有问题么

如果计算资源不够 涉及反应中心的原子也要尽可能给足够高的基组，比如tzvp，原子数也不多，不会对耗时产生太大的影响，远离反应中心的原子可以用低基组，比如6-31g*

Sally-L

sai77 发表于 2025-6-9 11:44
那当然要统一基组
其他不含金属的小分子优化时在相同基组下加上5D关键词就可以了

请问络合物中的金属和配体是分开优化的，那么计算单点能的时候是统一用一套泛函和基组还是金属和配体分开

sai77

Sally-L 发表于 2025-6-9 16:39
请问络合物中的金属和配体是分开优化的，那么计算单点能的时候是统一用一套泛函和基组还是金属和配体分开

不明白你的意思
opt+freq：所有分子泛函基组统一，含金属的PBE0-D3BJ/genecp(6-31G*+SDD)，不含金属的PBE0-D3BJ/6-31G*(5D)
单点能：所有分子泛函基组统一，PBE0-D3BJ/def2TZVP

Sally-L

sai77 发表于 2025-6-10 09:10
不明白你的意思
opt+freq：所有分子泛函基组统一，含金属的PBE0-D3BJ/genecp(6-31G*+SDD)，不含金属的PB ...

谢谢