因为没看到文献本身，您这个结构的POSCAR也没贴上来，所以我只能猜一下，有以下几种可能的原因。
1，晶格常数本身优化过吗？包括完美的以及引入了缺陷的都要优化晶格常数。
2，如果晶格常数已经完全弛豫了，但和文献里的形成能不一样，那有可能文献是簇模型，您的是周期性模型，形成能的确会不一样。当然如果文献用的软件不是vasp的话形成能也会不一样。
3，文献中同样是用vasp的话，可能您的这个slab的大小和文献里的不一样，那么等同于缺陷的密度不一样，那么形成能也会有差异。
4，计算参数和文献不一样，或者文献Ef是用的考虑了溶剂化的能量算的。
空位缺陷本身结构的差异从图中我用肉眼有些不太看的出来，所谓的“键没有连起来”其实只是可视化软件的视觉效果，不能算作结构不一样的依据。

876449830 发表于 2024-4-18 10:30
二维材料要固定晶轴优化

是把Z轴固定嘛，我固定了算出来三个悬空键中的两个还是没连起来，形成能也不对

二维材料要固定晶轴优化