关于过渡金属体系所用泛函的求助

Papillon

请教大家，含有开壳层过渡金属离子的200-300个原子左右的非有机体系（含有几十个H2O分子)，由于计算昂贵，在pbe和blyp两种纯泛函之间抉择。
我两种方法计算发现pbe比blyp反应能垒只高几kcal/mol，但是浏览论坛发现大多推荐pbe方法，是pbe对过渡金属体系的电子和结构描述更加准确吗，还是什么其他原因呢？
请问大家对该体系所用方法有什么建议吗？谢谢大家

sobereva

当前问题没体现出和CP2K有直接联系，别选CP2K分类。这次给你改了，以后注意

如果用CP2K，两三百个原子的杂化计算是小儿科的事，毫无理由为节约时间用纯泛函。前提是没有犯低级错误，相关信息：
CP2K做杂化泛函计算的关键要点和简单例子
http://sobereva.com/690（http://bbs.keinsci.com/thread-41087-1-1.html）

便宜而且相对普适的PBE算是第一性原理领域默认用的泛函，没有理由改成BLYP，本来BLYP算金属有机反应也很烂。PBE算金属有机反应也算不上好，非要纯泛函还不如r2SCAN。然而你也没说清楚具体是什么反应，到底是不是金属有机反应（体系含过渡金属显然不代表参与反应的就是过渡金属），没法给出准确推荐。提问这种问题必须把反应交代得具体。

Papillon

sobereva 发表于 2025-1-11 00:04
当前问题没体现出和CP2K有直接联系，别选CP2K分类。这次给你改了，以后注意

如果用CP2K，两三百个原子的 ...

谢谢sob老师，是我表述不清。
确实是用CP2K，但不是做单点计算，而是用CP2K跑限制性AIMD来看反应势能面，反应体系是过渡金属和S，周围是有几十个水分子的液滴环境

sobereva

Papillon 发表于 2025-1-13 19:50
谢谢sob老师，是我表述不清。
确实是用CP2K，但不是做单点计算，而是用CP2K跑限制性AIMD来看反应势能面 ...

跑AIMD就更不适合在CP2K里用BLYP