Metal–water interface formation: Thermodynamics from ab initio molecular dynamics simulations

Metal–water interfaces are central to many electrochemical, (electro)catalytic, and materials science processes and systems. However, our current understanding of their thermodynamic properties is limited by the scarcity of accurate experimental and computational data and procedures. In this work, thermodynamic quantities for metal–water interface formation are computed for a range of FCC(111) surfaces (Pd, Pt, Au, Ag, Rh, and PdAu) through extensive density functional theory based molecular dynamics and the two-phase entropy model. We show that metal–water interface formation is thermodynamically favorable and that most metal surfaces studied in this work are completely wettable, i.e., have contact angles of zero. Interfacial water has higher entropy than bulk water due to the increased population of low-frequency translational modes. The entropic contributions also correlate with the orientational water density, and the highest solvation entropies are observed for interfaces with a moderately ordered first water layer; the entropic contributions account for up to ∼25% of the formation free energy. Water adsorption energy correlates with the water orientation and structure and is found to be a good descriptor of the internal energy part of the interface formation free energy, but it alone cannot satisfactorily explain the interfacial thermodynamics; the interface formation is driven by the competition between energetic and entropic contributions. The obtained results and insight can be used to develop, parameterize, and benchmark theoretical and computational methods for studying metal–water interfaces. Overall, our study yields benchmark-quality data and fundamental insight into the thermodynamic forces driving metal–water interface formation.