关于在高斯中计算金属团簇对气体分子的吸附所选用的基组赝势问题

Zt2020666

最近在算一个金属Al与气体分子的吸附，气体分子是碳氧化物，氮氧化物，查看了一些论坛里的资料，虽然大致有方向，但是具体问题具体分析吧。就这个体系优化应该在什么水平优化呢，尤其是对于金属铝。这里金属铝是建立了簇模型，在簇模型上放置的气体分子，体系是中性。是否可以在B3LYP-D3(BJ)/ma-TZVP水平进行优化？或者用PBE0替换B3LYP？

sobereva

问这种问题始终都应当把计算资源、体系特征（最好给截图，至少说化学组成）给出来。适合的计算级别要在精度、耗时方面做权衡最终综合决定。
诸如氢键复合物，对较大体系6-311G**优化足矣，而对于很小的氢键复合物如水二聚体，同样的计算条件下明显可以轻易达到更高精度，故6-311G**优化就显得略糙，而适合用def2-TZVPP。

snljty2

Zt2020666 发表于 2022-5-20 14:28
能算，现在在算结果，这样得到的能量是否可以直接算吸附能，并分析分子与Al团簇间的相互作用，包括成键相 ...

后者得做EDA研究。前者应该优化后在更高级别下，并考虑BSSE计算吸附能，诸如wB97M-V/def2-TZVPP+CounterPoise校正，或者def2-QZVPP如果算的动，或者ma-def2-TZVPP但是自洽场收敛和辅助基组的线性依赖是问题。

Zt2020666

snljty2 发表于 2022-5-20 13:59
B3LYP-D3(BJ)/6-311G*甚至小道6-31G*优化足以，最大def2-TZVP足以，不需要弥散，团簇不要建太小。考虑DFT-D ...

能算，现在在算结果，这样得到的能量是否可以直接算吸附能，并分析分子与Al团簇间的相互作用，包括成键相互作用，弱相互作用

snljty2

B3LYP-D3(BJ)/6-311G*甚至小道6-31G*优化足以，最大def2-TZVP足以，不需要弥散，团簇不要建太小。考虑DFT-D3的三体校正项可能会有改善，但是程序不一定支持，或者不一定能算解析频率。