Polaronic structure of excess electrons and holes for a series of bulk iron oxides

Iron oxides such as hematite (α-Fe2O3) play an important role in diverse fields ranging from biogeochemistry to photocatalysis. Here we perform calculations of both the electron and electron hole polaron structures and associated reorganisation energies for a series of bulk iron oxides: hematite (α-Fe2O3), lepidocrocite (γ-FeOOH), goethite (α-FeOOH) and white rust (Fe(OH)2). Through the use of gap-optimized hybrid functionals and large supercells under periodic boundary conditions, we remove some of the complications and uncertainties present in earlier cluster model calculations. It is found that while the electron hole polaron in these materials generally localises onto a single iron site, the electron polaron localises across two iron sites of the same spin layer as a consequence of the lower reorganisation energy for electrons compared to holes. An exception to these trends is the hole of goethite, which according to our calculations does not form a localised polaron