Propensity of hydroxide and hydronium ions for the air–water and graphene–water interfaces from <i>ab initio</i> and force field simulations

Understanding acids and bases at interfaces is relevant for a range of applications from environmental chemistry to energy storage. We present combined ab initio and force-field molecular dynamics simulations of hydrochloric acid and sodium hydroxide highly concentrated electrolytes at the interface with air and graphene. In agreement with surface tension measurements at the air-water interface, we find that HCl presents an ionic surface excess, while NaOH displays an ionic surface depletion, for both interfaces. We further show that graphene becomes less hydrophilic as the water ions concentration increases, with a transition to being hydrophobic for highly basic solutions. For HCl, we observe that hydronium adsorbs to both interfaces and orients strongly toward the water phase, due to the hydrogen bonding behavior of hydronium ions, which donate three hydrogen bonds to bulk water molecules when adsorbed at the interface. For NaOH, we observe density peaks of strongly oriented hydroxide ions at the interface with air and graphene. To extrapolate our results from concentrated electrolytes to dilute solutions, we perform single ion-pair ab initio simulations, as well as develop force-field parameters for ions and graphene that reproduce the density profiles at high concentrations. We find the behavior of hydronium ions to be rather independent of concentration. For NaOH electrolytes, the force-field simulations of dilute NaOH solutions suggest no hydroxide adsorption but some adsorption at high concentrations. For both interfaces, we predict that the surface potential is positive for HCl and close to neutral for NaOH.