Resolving the odd–even oscillation of water dissociation at rutile TiO2(110)–water interface by machine learning accelerated molecular dynamics

Aqueous rutile TiO2(110) is the most widely studied water–oxide interface, and yet questions about water dissociation are still controversial. Theoretical studies have systematically investigated the influence of the slab thickness on water dissociation energy (Ediss) at 1 monolayer coverage using static density functional theory calculation and found that Ediss exhibits odd–even oscillation with respect to the TiO2 slab thickness. However, less studies have accounted for the full solvation of an aqueous phase using ab initio molecular dynamics due to high computational costs in which only three, four, and five trilayer models of rutile(110)–water interfaces have been simulated. Here, we report Machine Learning accelerated Molecular Dynamics (MLMD) simulations of defect-free rutile(110)–water interfaces, which allows for a systematic study of the slab thickness ranging from 3 to 17 trilayers with much lower costs while keeping ab initio accuracy. Our MLMD simulations show that the dissociation degree of surface water (α) oscillates with the slab thickness and converges to ∼2% as the TiO2 slab becomes thicker. Converting α into dissociation free energy (ΔAdiss) and comparing with dissociation total energy Ediss calculated with a single monolayer of water, we find that the full solvation of the interfaces suppresses surface water from dissociating. It is interesting to note that the machine learning potential trained from the dataset containing exclusively the five trilayer TiO2 model exhibits excellent transferability to other slab thicknesses and further captures the oscillating behavior of surface water dissociation. Detailed analyses indicate that the central plane in odd trilayer slabs modulates the interaction between double trilayers and, thus, the bonding strength between terminal Ti and water, which affects pKa of surface water and water dissociation degree.