Computer simulation of electron transfer at hematite surfaces

Molecular dynamics simulations in combination with ab initio calculations were carried out to determine the rate of electron transfer at room temperature in bulk hematite (α-Fe2O3) and at two low-index surfaces, namely the (012) and (001) surfaces. The electron transfer reactions considered here involve the II/III valence interchange between nearest-neighbor iron atoms. Two electron transfer directions were investigated, namely the basal plane and c direction electron transfers. Electron transfer rates obtained in bulk hematite were in good agreement with ab initio electronic structure calculations thus validating the potential model. The surfaces were considered both in vacuum and in contact with an equilibrated aqueous solution. The reorganization energy is found to increase significantly at the first surface layer and this value is little affected by the presence of water. In addition, in the case of the (012) surface, the electronic coupling matrix element for the topmost basal plane transfer was calculated at the Hartree-Fock level and was found to be weak compared to the corresponding electron transfer in the bulk. Therefore, most surfaces show a decrease in the rate of electron transfer at the surface. However, where iron atoms involved in the electron transfer reaction are directly coordinated to water molecules, water lowers the free energy of activation to a great extent and provides a large driving force for electrons to diffuse toward the bulk thus opposing the intrinsic surface effect. The surfaces considered in this work show different electron transfer properties. Hematite has been shown to exhibit anisotropic conductivity and thus different surfaces will show different intra- and inter-layer rates depending on their orientation. Moreover, the calculations of electron transfers at the hydroxyl- and iron-terminated (001) surfaces revealed that surface termination has a significant effect on the electron transfer parameters in the vicinity of the surface. Finally, our findings indicate that undercoordinated terminal iron atoms could act as electron traps at the surface.