Publication date: 22nd July 2026
Bismuth vanadate (BiVO4) represents a promising candidate for photoelectrochemical (PEC) water splitting. It has been shown that polarons play an important role in the water-splitting mechanism by introducing charge transition levels (CTLs) in the band gap and modifying the band alignment. Additionally, the
introduction of oxygen vacancies, which can interact with electron polarons, has been reported to improve the PEC efficiency. However, explicit simulations of water interfaces require large numbers of atoms, and calculations of CTLs involve integration of free energies over long timescales, making ab initio methods
prohibitive. In this work, we train a machine-learning potential that captures electron polarons and oxygen vacancies in multiple charge states to investigate their stability at an explicit BiVO4-water interface. By computing finite-temperature free energies, CTLs, and oxygen-vacancy formation energies as a function of distance from the interface, we determine how the aqueous environment modifies charge trapping and defect energetics relative to bulk BiVO4. While the free-energy differences are temperature independent in bulk BiVO4 due to the harmonic response of the VO4 units, we find that the aqueous interface introduces a clear temperature dependence. By further analyzing the local structural distortions associated with electron localization, we provide a microscopic link between interfacial structure and the stability of localized charges. Finally, we observe vacancy migration near the interface, highlighting the strong impact of the aqueous interface on the stability of oxygen vacancies.
