求助，理论上无单电子的slab出现mag不为0的问题

XIEH

最近对一个正交晶系的{013}晶面进行了计算，之前都是对{100}、{010}、{001}这种低指数且较规整晶面进行计算，符合mag=0
但是{013}切面切出来由于角度的原因，导致有些最小重复单元缺损，比如最底下的层会出现O原子没有金属原子配位或者底层金属原子配位数不足，
分析OUTCAR的magnetization发现体系mag不为0的原因就来自于底层这些配位数较低的原子

这样的表面能否进行计算？是否科学合理？如何解决？
目前想法是在底部也添加同样厚度的真空层，但不确定是否可行，且考虑到计算资源过大，烦请各位大佬支支招

XIEH

@卡开发发 请教一下卡开发发老师

卡开发发

XIEH 发表于 2025-6-18 19:15
@卡开发发 请教一下卡开发发老师

看看具体模型？

XIEH

卡开发发 发表于 2025-6-19 09:46
看看具体模型？

抱歉，晚上才看到。
像这样的，V2O5{013}

卡开发发

XIEH 发表于 2025-6-19 20:45
抱歉，晚上才看到。
像这样的，V2O5{013}

你要不然还是把完整结构传一下。

XIEH

卡开发发 发表于 2025-6-19 22:38
你要不然还是把完整结构传一下。

我前天又重新切了一个面013re，一天多了还没优化到力小于0.05 eV/Ang

卡开发发

XIEH 发表于 2025-6-21 09:50
我前天又重新切了一个面013re，一天多了还没优化到力小于0.05 eV/Ang

切面切的有问题：
1、底部有O原子，几乎不与V相连；
2、不符合化学计量比。

XIEH

卡开发发 发表于 2025-6-21 10:03
切面切的有问题：
1、底部有O原子，几乎不与V相连；
2、不符合化学计量比。

感谢老师指点