Monte Carlo and Molecular Dynamics Simulations of Electrical Double-Layer Structure in Potassium−Montmorillonite Hydrates

Abstract

Monte Carlo and molecular dynamics simulations of interlayer molecular structure in the one-, two-, and (hypothetical) three-layer hydrates of K-montmorillonite were performed concurrently in order to elucidate counterion speciation and water structure in the electrical double layer of this clay mineral. Calculated layer spacings, interlayer water potential energies, and counterion mobilities were in agreement with available experimental data. In the one-layer hydrate, both outer-sphere and inner-sphere surface complexes of K+ were observed, the latter always near sites of tetrahedral charge substitution, with the counterion species exchanging readily on the simulation time scale (up to 200 ps). In the two- and three-layer hydrates, the surface complexes persisted, but an incipient diffuse layer of counterions also was observed, with all three types of surface species engaging in sluggish exchange. Water molecules in the one-layer hydrate resided at the interlayer midplane, whereas in the two-layer hydrat...