The fate of chemical and radioactive wastes in the environment is related to the ability of natural phases to attenuate and immobilize contaminants through chemical sorption and precipitation processes. Our understanding of these complex processes at the atomic level is provided by a few experimental and analytical methods such as X-ray absorption and NMR spectroscopies. However, due to complexities in the structure and composition of clay and other hydrated minerals, and the inherent uncertainties of the experimental methods, it is important to apply theoretical molecular models for a fundamental atomic-level understanding, interpretation, and prediction of these phenomena. In this effort, we have developed a general force field, CLAYFF, suitable for the simulation of hydrated and multicomponent mineral systems and their interfaces with aqueous solutions. Interatomic potentials were derived from parametrizations incorporating structural and spectroscopic data for a variety of simple hydrated compounds. A flexible SPC-based water model is used to describe the water and hydroxyl behavior. Metal-oxygen interactions are described by a LennardJones function and a Coulombic term with partial charges derived by Mulliken and ESP analysis of DFT results. Bulk structures, relaxed surface structures, and intercalation processes are evaluated and compared to experimental and spectroscopic findings for validation. Our approach differs from most others in that we treat most interatomic interactions as nonbonded. This allows us to effectively use the force field for a wide variety of phases and to properly account for energy and momentum transfer between the fluid phase and the solid, while keeping the number of parameters small enough to allow modeling of relatively large and highly disordered systems. Simulations of clay, hydroxide, and oxyhydroxide phases and their interfaces with aqueous solutions combine energy minimization and molecular dynamics methods to describe the structure and behavior of water, hydroxyl, surface species, and intercalates in these systems. The results obtained to date demonstrate that CLAYFF shows good promise to evolve into a widely adaptable and broadly effective force field for molecular simulations of fluid interfaces with clays and other clay-related phases, as well as other inorganic materials characterized by complex, disordered, and often ill-determined structure and composition.

Molecular dynamics (MD) computer simulations of liquid water adsorbed on the muscovite (001) surface provide a greatly increased, atomistically detailed understanding of surface-related effects on the spatial variation in the structural and orientational ordering, hydrogen bond (H-bond) organization, and local density of H2O molecules at this important model phyllosilicate surface. MD simulations at constant temperature and volume (statistical NVT ensemble) were performed for a series of model systems consisting of a two-layer muscovite slab (representing 8 crystallographic surface unit cells of the substrate) and 0 to 319 adsorbed H2O molecules, probing the atomistic structure and dynamics of surface aqueous films up to 3 nm in thickness. The results do not demonstrate a completely liquid-like behavior, as otherwise suggested from the interpretation of X-ray reflectivity measurements and earlier Monte Carlo simulations. Instead, a more structurally and orientationally restricted behavior of surface H2O molecules is observed, and this structural ordering extends to larger distances from the surface than previously expected. Even at the largest surface water coverage studied, over 20% of H2O molecules are associated with specific adsorption sites, and another 50% maintain strongly preferred orientations relative to the surface. This partially ordered structure is also different from the well-ordered 2-dimensional ice-like structure predicted by ab initio MD simulations for a system with a complete monolayer water coverage. However, consistent with these ab initio results, our simulations do predict that a full molecular monolayer surface water coverage represents a relatively stable surface structure in terms of the lowest diffusional mobility of H2O molecules along the surface. Calculated energies of water adsorption are in good agreement with available experimental data.

Naturally occurring clay minerals provide a distinctive material for carbon capture and carbon dioxide sequestration. Swelling clay minerals, such as the smectite variety, possess an aluminosilicate structure that is controlled by low-charge layers that readily expand to accommodate water molecules and, potentially, CO2. Recent experimental studies have demonstrated the efficacy of intercalating CO2 in the interlayer of layered clays, but little is known about the molecular mechanisms of the process and the extent of carbon capture as a function of clay charge and structure. A series of molecular dynamics simulations and vibrational analyses have been completed to assess the molecular interactions associated with incorporation of CO2 and H2O in the interlayer of montmorillonite clay and to help validate the models with experimental observation. An accurate and fully flexible set of interatomic potentials for CO2 is developed and combined with Clayff potentials to help evaluate the intercalation mechanism and examine the effect of molecular flexibility on the diffusion rate of CO2 in water.