Abstract
Hematite (alpha-Fe2O3) is widely used as a catalytic electrode material in photoelectrochemical water oxidation, where its surface compositions and stabilities can strongly impact the redox reaction process. Here, its surface configurations in environmental or electrochemical conditions are assessed via density functional theory (DFT) calculations conducted at the Perdew, Burke, and Ernzerhof (PBE)+U level. The most energetically favorable surface domains of alpha-Fe2O3 (0001) and (1 (1) over bar 02) are predicted by constructing the surface phase diagrams in the framework of first-principle thermodynamics. The relative surface stabilities are investigated as a function of partial pressures of oxygen and water, temperature, solution pH, and electrode potential not only for perfect bulk terminations but also for defect-containing surfaces having various degrees of hydroxylation and hydration. In order to assess the impact on the redox reactions of the surface planes as well as of the extent of surface hydration/hydroxylation, the thermodynamics of the four-step oxygen evolution reaction (OER) mechanism are examined in detail for different models of the alpha-Fe2O3 (0001) and (1 (1) over bar 02) surfaces. Importantly, the results underline that the nature of the surface termination and the degree of near-surface hydroxylation give rise to significant variations in the OER overpotentials.