VASP | 理论计算常见问题解答(102)掺杂建模该先体相还是先表面?吸附能为正、中间体断键、态密度费米能级漂移怎么处理?

Q1:朱老师在晶胞里掺杂后,优化晶胞后再切面,再优化晶面;还是优化原胞后直接切面,面上掺杂后优化晶面就行了?

A:看研究的是体相掺杂还是表面掺杂,不能混用。

①若掺杂元素在材料内部均匀存在:先做体相掺杂优化,再切面建slab;

②若研究表面活性位点掺杂:可先切纯净slab,再替换表面原子优化;

③真实催化表面常需比较:表层、次表层、体相位点的掺杂能;

④最终吸附计算应选最稳定或实验最可能的掺杂表面。

Q2:朱老师,我这个Ag/al2o3模型吸附pms能量算的是正的,这个银是根据xrd找的立方构建的异质结。我是不是可以换一个银团簇来做,不用完全对应xrd卡片?

VASP | 理论计算常见问题解答(102)掺杂建模该先体相还是先表面?吸附能为正、中间体断键、态密度费米能级漂移怎么处理?

A:可以换Ag团簇,但不能为了让能量变负而随意换:

①XRD只能说明Ag主要是立方相,不能说明表面一定是无限平板;

②若实验是负载Ag纳米颗粒,用Ag团簇或Agn更合理;

③需比较Ag(111)、Ag(100)、Ag团簇、Ag/Al₂O₃界面位点;

④正吸附能也可能来自PMS电荷、溶剂、吸附位点或能量公式设置问题。

Q3:朱老师,这个niooh的h总是往下跑。这样可以么?
VASP | 理论计算常见问题解答(102)掺杂建模该先体相还是先表面?吸附能为正、中间体断键、态密度费米能级漂移怎么处理?
A:可以,但要先判断是合理氢键,还是周期性假象:

若 H 只是朝下方 O 偏移,形成 O–H···O 氢键,通常正常;

若 H 直接转移到底层/镜像层,说明层距或真空可能不够;

查 O–H 键长是否约 1.0 Å,H···O 是否合理。

Q4:老师我计算CuO(111)晶面的态密度,自洽计算(K点密度为2 2 1)后DOSCAR中费米能级为0.00149628,然后拷贝一份进行态密度计算(使用自洽计算的电荷密度和波函数,修改ISMEAR = -5,K点密度改为 3 3 1),计算后费米能级变为0.66669597。态密度进行费米能级矫正后,做图时,能量0点落在了导带底,如下图,这是因为我的参数改变导致吗?那这种0点在导带底的态密度是正确的吗?一般费米能级不是都在价带顶吗?是我哪一步有疏漏需要重新计算吗?

VASP | 理论计算常见问题解答(102)掺杂建模该先体相还是先表面?吸附能为正、中间体断键、态密度费米能级漂移怎么处理?
A:主要是参数和费米能级定义导致的,不一定错,半导体DOS更建议以VBM对齐。

ISMEAR=-5 和更密 KPOINTS 会重新确定占据,E-fermi 变化正常;

半导体/绝缘体带隙内没有态,费米能级不唯一,不一定非在价带顶;

若0点落在导带底但费米处DOS≈0,可能只是对齐方式问题,若DOS不为0,则有表面/缺陷态或金属性;

建议DOS步设 ICHARG=11,检查参数,作图时可改用 VBM=0。

Q5老师,加了铁这个ooh总是断开 有啥办法么。试了好多种吸附位置了?

VASP | 理论计算常见问题解答(102)掺杂建模该先体相还是先表面?吸附能为正、中间体断键、态密度费米能级漂移怎么处理?

A:反复断开,说明 Fe 位上 *OOH 很可能本身不稳定,别硬压成中间体。

先确认初始构型:应是 Fe–O–O–H,O–O≈1.45 Å,O–H≈1.0 Å;

若优化后变成 O₂ 脱附或 H 转移,说明该位点更倾向放氧/断键;

可先小步优化:POTIM=0.05–0.1,短暂固定 O–O/O–H 再放开;

若放开仍断,说明 *OOH 不是稳定态,应换活性位或改反应路径分析。

 

声明:如需转载请注明出处(华算科技旗下资讯学习网站-学术资讯),并附有原文链接,谢谢!
(0)
上一篇 1天前
下一篇 2026年3月19日 上午11:58

相关推荐