VASP理论计算常见问题（85）：CIF 查找、MAGMOM 设置、过渡态 ZPE 修正等计算关键问题答疑

A：MP给的两个NiMoO₄，严格说都不是多数电化学文献里常说的高温β-NiMoO₄结构；

SpringerMaterials 有专门条目 “β-NiMoO₄(Ni[MoO₄] ht)”，页面写明可下载 CIF；

很多电化学文章把 β 指向 JCPDS 45-0142。

Q2：朱老师，在算Ni（111）表面的吸附能时，MAGMOM参数该如何填写，假如是4*4*4的晶胞，全部默认是1嘛？还是怎么写？

A：建议进行参数测试，选择最合理的那个；

算 Ni(111) 表面吸附能时，MAGMOM 可以把它当成“初始磁矩猜测”：写得不需要非常精确，但要能把体系推到正确的铁磁（FM）解上；

如果不写 MAGMOM，它会用内部默认，有时会让 Ni 收敛到低磁矩甚至非磁解。所以做 Ni 表面吸附能，建议还是显式写上。

Q3：朱老师，请问过渡态计算，能量如果进行zpe修正，是不是还需要给过渡态单独算一下zpe？

A：是的：

如果想把过渡态能垒也做成 ZPE 修正后的能垒，那需要给过渡态（TS）单独算 ZPE；

TS 的振动频率里会有且只应有 1 个虚频，只把所有实频振动模加入 ZPE，那个虚频不参与 ZPE 求和，它对应反应坐标，不是稳定振动。

Q4：朱老师，文献中的d带中心一般是指整体团簇的d带中心平均值呢还是某个特定吸附原子的d带中心呢？

A：一般是平均;

指的是某一类“活性金属位点”的 d 带中心；

建议计算吸附位点原子，或直接成键的那几个金属原子的 d 带中心。

Q5：朱老师，做吸附能计算 那个切面应该用完美晶胞的来切面再在110面上创造氧空位吧这样对吗 ？

A：对的

先用无缺陷bulk去切 (110) 表面，再在表面上造氧空位；

氧空位是表面缺陷，它的结构，电荷重排强烈依赖表面环境；

如果在 bulk 里先造空位再切 slab，经常会变成缺陷位置不受控、甚至缺陷跑到不想要的层、并且更难保证终止面一致。