Nonequilibrium Ab Initio Molecular Dynamics Simulation of Water Splitting at Fe₂O₃-Hematite/Water Interfaces in an...

In the exploration of the optimal material for achieving the photoelectrochemical dissociation of water into hydrogen, hematite (α-Fe₂O₃) emerges as a highly promising candidate for proof-of-concept demonstrations. Recent studies suggest that the concurrent application of external electric fields could enhance the photoelectrochemical (PEC) process. To delve into this, we conducted nonequilibrium ab initio molecular dynamics (NE-AIMD) simulations in this study, focusing on hematite-water interfaces at room temperature under progressively stronger electric fields. Our findings reveal intriguing evidence of water molecule adsorption and dissociation, as evidenced by an analysis of the structural properties of the hydrated layered surface of the hematite-water interface. Additionally, we scrutinized intermolecular structures using radial distribution functions (RDFs) to explore the interaction between the hematite slab and water. Notably, the presence of a Grotthuss hopping mechanism became apparent as the electric field strength increased. A comprehensive discussion based on intramolecular geometry highlighted aspects such as hydrogen-bond lengths, H-bond angles, average H-bond numbers, and the observed correlation existing among the hydrogen-bond strength, bond-dissociation energy, and H-bond lifetime. Furthermore, we assessed the impact of electric fields on the librational, bending, and stretching modes of hydrogen atoms in water by calculating the vibrational density of states (VDOS). This analysis revealed distinct field effects for the three characteristic band modes, both in the bulk region and at the hematite-water interface. We also evaluated the charge density of active elements at the aqueous hematite surface, delving into field-induced electronic charge-density variations through the Hirshfeld charge density analysis of atomic elements. Throughout this work, we drew clear distinctions between parallel and antiparallel field alignments at the hematite-water interface, aiming to elucidate crucial differences in local behavior for each surface direction of the hematite-water interface.