We have shown that two models, a fully planar and a convex one, can explain the common Mössbauer doublets D1, D1b, D2, and the operando doublet D4, however with different influences of the axial ligands. While two different model sizes in terms of the equatorial carbon-based ligand environment were studied, the differences in the resulting Mössbauer parameters were too subtle to distinguish between these model sizes. This is consistent with Mössbauer spectroscopy being a local probe, and previous work that studied model size influences .