It had been found that spinorbit coupling components were small in FeO systems, [44,45] therefore collinear spin model was used in the magnetene layers. The Brillouin zone integrations for the bulk and slab were carried out on a converged Monkhorst pack k-mesh of sizes 3 × 3 × 3 and 4 × 6 × 1, respectively.