VASP | 理论计算常见问题解答(105):催化体系 DFT 计算常见问题解答汇总

Q1:朱老师,请问在计算Ag催化剂吸附pms bader电荷的时候,总体的电荷转移了0.8e。银是用ag6团簇做的,因为一个e正好对应ag→ag+的电荷数,一个pms也就对应一个ag,导体电荷相当于平均分了,这个bader电荷是看整体0.8e还是说六个的平均值0.15e呢?

A:这个结果应优先看整体转移量,平均值只能作为辅助说明:

PMS吸附活化时,0.8e代表Ag6整体向PMS供电子;

六个Ag平均成0.15e,不能等同单个Ag真实氧化;

若讨论Ag价态变化,应看配位Ag的局域电荷变化;

Bader电荷不宜直接对应Ag变Ag正离子的整数电荷。

Q2:朱老师,我想在催化剂pms体系下加一个cl-,看看对整个体系的影响,加氯原子,INCAR里面的necelect需要加一个电荷吗?

A:如果要模拟Cl-,需要在NELECT中多加一个电子,但要注意带电slab误差:

只加Cl原子不改NELECT,实际算的是中性氯物种;

若目标是Cl-影响,NELECT应比默认值多一个电子;  

体系会有均匀背景电荷,能量不能直接随意比较;

建议配合大真空层和偶极矫正,重点比较相同设置结果。

Q3:朱老师,这个报错什么情况啊 ?
VASP | 理论计算常见问题解答(105):催化体系 DFT 计算常见问题解答汇总
A:这个报错一般是电子步对角化不稳定,优先排查结构和重启文件:

先检查原子是否过近,吸附物是否已经异常跑偏;

删除WAVECAR和CHGCAR后,用干净电荷重新启动;

若在结构优化中出现,适当减小POTIM再继续测试;

若反复出现,可改用ALGO等于Normal提高稳定性。

Q4:朱老师,我在计算双原子催化剂吸附S8的结构优化过程中,S8环会断成好几块,这是什么原因呀?应该怎么调整呢?

A:S8断裂可能是真实强吸附解离,也可能是优化设置诱发的假断裂:

先看优化轨迹,若S S键逐步拉长,更像真实活化解离;

若一开始剧烈断裂,常见原因是初始距离过近或步长过大;

可以先把S8放远预优化,再逐步靠近双原子活性位点;

建议减小POTIM,并分阶段固定部分S8或催化剂骨架。

Q5:朱老师,台阶图的高低和bader转移电子数的多少有关系吗?

A:台阶图高低和Bader电荷有关联,但不能简单画等号:

台阶图高低通常表示自由能差,不是电荷转移量本身;

Bader电荷反映电子重新分布,只是影响能量的因素之一;

电荷转移多可能增强吸附,但不一定让能垒或能阶降低;

判断催化性能仍应以自由能变化和关键步骤能垒为主。

声明:如需转载请注明出处(华算科技旗下资讯学习网站-学术资讯),并附有原文链接,谢谢!
赞 (0)
上一篇 2026年6月11日 上午10:02
下一篇 2026年6月16日 上午10:34

相关推荐