Roles of Hydration and Magnetism on the Structure of Ferrihydrite from First Principles

We report density functional theory calculations aimed at predicting thermodynamically stable structures for ferrihydrite across a range of possible compositions determined by the amount of structural water. Based on an assumed formula unit of Fe₅O₈H + nH₂O, we performed ab initio calculations with evolutionary searching to find the lowest enthalpy structures as a function of the water content up to n = 2. This is the most exhaustive search for the ferrihydrite structure conducted so far; more than 5000 unique configurations were generated and evaluated over five states of hydration. Among them, the Michel akdaliate model was generated, along with several energetically comparable new structures at higher states of hydration. However, a direct comparison between calculated and experimental pair distribution function and X-ray diffraction patterns for the 50 lowest energy structures shows that none beyond the Michel model could be associated with ferrihydrite. Nevertheless, this energetically comparable structure set provides a novel basis for analyzing and understanding the effects of hydration and magnetism on the topology of ferrihydrite, from which we conclude that any tetrahedral Fe should be viewed as a metastable structural defect, created either as a result of the rapid kinetics of crystal growth or to accommodate a local magnetic stress between neighboring Fe atoms.