求助：VASP计算Co3O4slab模型的磁矩设置问题

zhj1023

各位老师们好，在用vasp计算四氧化三钴时候，bulk计算可以根据配位环境选择对称原子来设置不同方向的磁矩；
然而在切完晶面之后变成slab模型底部的原子配位的氧原子数改变了（图2），这时候有点不明白应该如何设置磁矩去计算了。
目前尝试过直接前一半设置正的磁矩，后一半原子负的磁矩，结构优化能够收敛，但是后续自洽计算难以收敛（图1，需算了两次250scf还未收敛，EDIFF = 1E-06）。请问大家怎么设置类似体系的磁矩，谢谢！

Boble

这个挺麻烦，建议用脚本来做。Co3O4体相中，O原子不加磁性，八面体配位的Co也不加磁性，四面体配位的Co磁矩在±2.6uB左右，这个很多文献都能查得到。
由于Co3O4是反铁磁性晶体（AFM），所以四面体配位的Co一定要一半数量自旋向上，另一半自旋向下，这样能量最稳定，物理性质如带隙才最接近于实验值
你切了表面之后，slab的中间层中的原子类比于体相，磁矩和体相保持一致，比较麻烦的是表面暴露的配位不饱和原子。
表面保留的不饱和八面体配位的Co一般在±2-3uB之间，不饱和四面体配位的Co±3uB左右，这个得做SCF验证
总之找磁基态是一个很麻烦，以上都是个人经验，不保证对，请仔细斟酌，合理采纳