Pyramidal charged domain walls in ferroelectric BiFeO₃

Domain structures play a crucial role in the electric, mechanical and other properties of ferroelectric materials. In this study, we uncover the physical origins of the enigmatic zigzag domain structure in the prototypical multiferroic material BiFeO$_3$. Using phase-field simulations within the Landau-Ginzburg-Devonshire framework, we demonstrate that spatially-homogeneous defect charges result in domain structures that closely resemble those observed experimentally. The acquired understanding of the underlying physics of pyramidal-domain formation may enable the engineering of new materials with self-assembled domain structures exhibiting defined domain periodicity at the nanometre scale, opening avenues for advanced applications.