Narasimhan Loganathan

On the Hydrogen Bonding Structure at the Aqueous Interface of Ammonium-Substituted Mica: A Molecular Dynamics Simulation

Molecular dynamics (MD) computer simulations were performed for an aqueous film of 3 nm thickness adsorbed at the (001) surface of ammonium-substituted muscovite mica. The results provide a detailed picture of the near-surface structure and topological characteristics of the interfacial hydrogen bonding network. The effects of deuterium=hydrogen isotopic substitution in N(H=D)4+ on the dynamics and consequently on the convergence of the structural properties have also been explored. Unlike many earlier simulations, a much larger surface area representing 72 crystallographic unit cells was used, which allowed for a more realistic representation of the substrate surface with a more disordered distribution of aluminium=silicon isomorphic substitutions in muscovite. The results clearly demonstrate that under ambient conditions both interfacial ammonium ions and the very first layer of water molecules are H-bonded only to the basal surface of muscovite, but do not form H-bonds with each other. As the distance from the surface increases, the H-bonds donated to the surface by both N(H=D)4+ and H2O are gradually replaced by the H-bonds to the neighbouring water molecules, with the ammonia ions experiencing one reorientational transition region, while the H2O molecules experiencing three such distinct consecutive transitions. The hydrated N(H=D)4+ ions adsorb almost exclusively as inner-sphere surface complexes with the preferential coordination to the basal bridging oxygen atoms surrounding the aluminium=silicon substitutions.

Tipping Point for Expansion of Layered Aluminosilicates in Weakly Polar Solvents: Supercritical CO2

Layered aluminosilicates play a dominant role in the mechanical and gas storage properties of the subsurface, are used in diverse industrial applications, and serve as model materials for understanding solvent-ion-support systems. Although expansion in the presence of H2O is well-known to be systematically correlated with the hydration free energy of the interlayer cation, particularly in environments dominated by nonpolar solvents (i.e., CO2), uptake into the interlayer is not well-understood. Using novel high-pressure capabilities, we investigated the interaction of dry supercritical CO2 with Na-, NH4-, and Cs-saturated montmorillonite, comparing results with predictions from molecular dynamics simulations. Despite the known trend in H2O and that cation solvation energies in CO2 suggest a stronger interaction with Na, both the NH4- and Cs-clays readily absorbed CO2 and expanded, while the Na-clay did not. The apparent inertness of the Na-clay was not due to kinetics, as experiments seeking a stable expanded state showed that none exists. Molecular dynamics simulations revealed a large endothermicity to CO2 intercalation in the Na-clay but little or no energy barrier for the NH4- and Cs-clays. Indeed, the combination of experiment and theory clearly demonstrate that CO2 intercalation of Na-montmorillonite clays is prohibited in the absence of H2O. Consequently, we have shown for the first time that in the presence of a low dielectric constant, gas swelling depends more on the strength of the interaction between the interlayer cation and aluminosilicate sheets and less on that with solvent. The finding suggests a distinct regime in layered aluminosilicate swelling behavior triggered by low solvent polarizability, with important implications in geomechanics, storage, and retention of volatile gases, and across industrial uses in gelling, decoloring, heterogeneous catalysis, and semipermeable reactive barriers.