Probing atomic-scale processes at the ferrihydrite-water interface with reactive molecular dynamics

Interfacial processes involving metal (oxyhydr)oxide phases are important for the mobility and bioavailability of nutrients and contaminants in soils, sediments, and water. Consequently, these processes influence ecosystem health and functioning, and have shaped the biological and environmental co-evolution of Earth over geologic time. Here we employ reactive molecular dynamics simulations, supported by synchrotron X-ray spectroscopy to study the molecular-scale interfacial processes that influence surface complexation in ferrihydrite-water systems containing aqueous molybdate. We validate the utility of this approach by calculating surface complexation models directly from simulations. The reactive force-field captures the realistic dynamics of surface restructuring, surface charge equilibration, and the evolution of the interfacial water hydrogen bond network in response to adsorption and proton transfer. We find that upon hydration and adsorption, ferrihydrite restructures into a more disordered phase through surface charge equilibration, as revealed by simulations and high-resolution X-ray diffraction. We observed how this restructuring leads to a different interfacial hydrogen bond network compared to bulk water by monitoring water dynamics. Using umbrella sampling, we constructed the free energy landscape of aqueous molybdate adsorption at various concentrations and the deprotonation of the ferrihydrite surface. The results demonstrate excellent agreement with the values reported by experimental surface complexation models. These findings are important as reactive molecular dynamics opens new avenues to study mineral-water interfaces, enriching and refining surface complexation models beyond their foundational assumptions.