Molecular Dynamics Simulations of the Interfacial Region between Boehmite and Gibbsite Basal Surfaces and High Ionic Strength Aqueous Solutions

Abstract

Classical molecular dynamics (MD) simulations were used to study the interactions of up to 2 M NaCl and NaNO3 aqueous solutions with the presumed inert boehmite (010) and gibbsite (001) surfaces. The force field parameters used in these simulations were validated against density functional theory calculations of Na+ and Cl– hydrated complexes adsorbed at the boehmite (010) surface. In all the classical MD simulations and regardless of the ionic strength or the nature of the anion, Na+ ions were found to preferably form inner-sphere complexes over outer-sphere complexes at the aluminum (oxy)hydroxide surfaces, adsorbing closer to the surface than both water molecules and anions. In contrast, Cl– ions were predicted to distribute preferably in outer-sphere positions. The resulting asymmetry in adsorption strengths offers molecular-scale evidence for the observed isoelectric point (IEP) shift to higher pH at high ionic strength for aluminum (oxy)hydroxides. As such, the MD simulations also provided clear evidence against the assumption that the basal surfaces of boehmite and gibbsite are inert to background electrolytes. Finally, the MD simulations indicated that the different affinities of NO3– and Cl– for the surfaces might have macroscopic consequences, such as difference in the sensitivity of the IEP to the electrolyte concentration.