In this study, the diffusion of water and ions in the interlayer region of smectite clay minerals represents a direct probe of the type and strength of clay–fluid interactions. Interlayer diffusion also represents an important link between molecular simulation and macroscopic experiments. Here we use molecular dynamics simulation to investigate trends in cation and water diffusion in montmorillonite interlayers, looking specifically at the effects of layer charge, interlayer cation and cation charge (sodium or calcium), water content, and temperature. For Na-montmorillonite, the largest increase in ion and water diffusion coefficients occurs between the one-layer and two-layer hydrates, corresponding to the transition from inner-sphere to outer-sphere surface complexes. Calculated activation energies for ion and water diffusion in Na-montmorillonite are similar to each other and to the water hydrogen bond energy, suggesting the breaking of water–water and water–clay hydrogen bonds as a likely mechanism for interlayer diffusion. A comparison of interlayer diffusion with that of bulk electrolyte solutions reveals a clear trend of decreasing diffusion coefficient with increasing electrolyte concentration, and in most cases the interlayer diffusion results are nearly coincident with the corresponding bulk solutions. Trends in electrical conductivities computed from the ion diffusion coefficients are also compared.
During the Frio-I Brine Pilot CO2 injection experiment in 2004, distinct geochemical changes in response to the injection of 1600tons of CO2 were recorded in brine samples collected from the monitoring well. Previous geochemical modeling studies have considered dissolution of calcite and iron oxyhydroxides, or release of adsorbed iron, as the most likely sources of the increased ion concentrations. In this modeling study we explore possible alternative sources of the increasing calcium and iron, based on the data from the detailed petrographic characterization of the Upper Frio Formation "C". Particularly, we evaluate whether dissolution of pyrite and oligoclase (anorthite component) can account for the observed geochemical changes. Due to kinetic limitations, dissolution of pyrite and anorthite cannot account for the increased iron and calcium concentrations on the time scale of the field test (10 days). However, dissolution of these minerals is contributing to carbonate and clay mineral precipitation on the longer time scales (1000 years). We estimated that during the field test dissolution of calcite and iron oxide resulted in ~0.02wt.% loss of the reservoir rock mass. The reactive transport models were constructed for 25 and 59°C temperature and using Pitzer and B-dot activity correction methods. These models predict carbonate minerals, dolomite and ankerite, as well as clay minerals kaolinite, nontronite and montmorillonite, will precipitate in the Frio Formation "C" sandstone as the system progresses toward chemical equilibrium during a 1000-year period. Cumulative uncertainties associated with using different thermodynamic databases, activity correction models (Pitzer vs. B-dot), and extrapolating to reservoir temperature, are manifested in the difference in the predicted mineral phases. However, these models are consistent with regards to the total volume of mineral precipitation and porosity values which are predicted to within 0.002%.
The study of mineral-water interfaces is of great importance to a variety of applications including oil and gas extraction, gas subsurface storage, environmental contaminant treatment, and nuclear waste repositories. Understanding the fundamentals of that interface is key to the success of those applications. Confinement of water in the interlayer of smectite clay minerals provides a unique environment to examine the interactions among water molecules, interlayer cations, and clay mineral surfaces. Smectite minerals are characterized by a relatively low layer charge that allows the clay to swell with increasing water content. Montmorillonite and beidellite varieties of smectite were investigated to compare the impact of the location of layer charge on the interlayer structure and dynamics. Inelastic neutron scattering of hydrated and dehydrated cation-exchanged smectites was used to probe the dynamics of the interlayer water (200-900 cm-1 spectral region) and identify the shift in the librational edge as a function of the interlayer cation. Molecular dynamics simulations of equivalent phases and power spectra, derived from the resulting molecular trajectories, indicate a general shift in the librational behavior with interlayer cation that is generally consistent with the neutron scattering results for the monolayer hydrates. Both neutron scattering and power spectra exhibit librational structures affected by the location of layer charge and by the charge of the interlayer cation. Divalent cations (Ba2+ and Mg2+) characterized by large hydration enthalpies typically exhibit multiple broad librational peaks compared to monovalent cations (Cs+ and Na+), which have relatively small hydration enthalpies.
The swelling properties of smectite clay minerals are relevant to many engineering applications including environmental remediation, repository design for nuclear waste disposal, borehole stability in drilling operations, and additives for numerous industrial processes and commercial products. We used molecular dynamics and grand canonical Monte Carlo simulations to study the effects of layer charge location, interlayer cation, and temperature on intracrystalline swelling of montmorillonite and beidellite clay minerals. For a beidellite model with layer charge exclusively in the tetrahedral sheet, strong ion-surface interactions shift the onset of the two-layer hydrate to higher water contents. In contrast, for a montmorillonite model with layer charge exclusively in the octahedral sheet, weaker ion-surface interactions result in the formation of fully hydrated ions (two-layer hydrate) at much lower water contents. Clay hydration enthalpies and interlayer atomic density profiles are consistent with the swelling results. Water adsorption isotherms from grand canonical Monte Carlo simulations are used to relate interlayer hydration states to relative humidity, in good agreement with experimental findings.
Heterogeneous nucleation of methane hydrates has been examined using molecular simulation, experimental bulk synthesis, and scanning probe microscopy. Theoretical nucleation rates were determined using molecular dynamics simulations as a function of clay surface represented by hydrophobic and hydrophilic systems. Methane hydrates were synthesized with and without Na-montmorillonite in a bulk reactor pressure assembly. X-ray diffraction and Raman spectroscopy confirm the nucleation and growth of the synthesized hydrates. Various kinetic pathways were explored to produce methane or isobutene clathrates in an ultra-high vacuum apparatus at very low temperatures but scanning probe microscopy only indicates the formation of ice.
In geologic settings relevant to a number of extraction and potential sequestration processes, nanopores bounded by clay mineral surfaces play a critical role in the transport of aqueous species. Solution structure and dynamics at clay-water interfaces are quite different from their bulk values, and the spatial extent of this disruption remains a topic of current interest. We have used molecular dynamics simulations to investigate the structure and diffusion of aqueous solutions in clay nanopores approximately 6 nm thick, comparing the effect of clay composition with model Na-hectorite and Na-montmorillonite surfaces. In addition to structural properties at the interface, water and ion diffusion coefficients were calculated within each aqueous layer at the interface, as well as in the central bulk-like region of the nanopore. The results show similar solution structure and diffusion properties at each surface, with subtle differences in sodium adsorption complexes and water structure in the first adsorbed layer due to different arrangements of layer hydroxyl groups in the two clay models. Interestingly, the extent of surface disruption on bulk-like solution structure and diffusion extends to only a few water layers. A comparison of sodium ion residence times confirms similar behavior of inner-sphere and outer-sphere surface complexes at each clay surface, but 1% of sodium ions adsorb in ditrigonal cavities on the hectorite surface. The presence of these anhydrous ions is consistent with highly immobile anhydrous ions seen in previous nuclear magnetic resonance spectroscopic measurements of hectorite pastes.
Methane hydrates are extremely important naturally-occurring crystalline materials that impact climate change, energy resources, geological hazards, and other major environmental issues. Whereas significant experimental effort has been completed to understanding the bulk thermodynamics of methane hydrate assemblies, little is understood on heterogeneous nucleation and growth of methane hydrates in clay-rich environments. Controlled synthesis experiments were completed at 265-285 K and 6.89 MPa to examine the impact of montmorillonite surfaces in clay-ice mixtures to nucleate and form methane hydrate. The results suggest that the hydrophilic and methane adsorbing properties of Namontmorillonite reduce the nucleation period of methane hydrate formation in pure ice systems.
Organic dyes such as methylene blue (MB) are often used in the characterization of clays and related minerals, but details of the adsorption mechanisms of such dyes are only partially understood from spectroscopic data, which indicate the presence of monomers, dimers, and higher aggregates for varying mineral surfaces. A combination of quantum (density functional theory) and classical molecular simulation methods was used to provide molecular detail of such adsorption processes, specifically the adsorption of MB onto kaolinite basal surfaces. Slab models with vacuum-terminated surfaces were used to obtain detailed structural properties and binding energies at both levels of theory, while classical molecular dynamics simulations of aqueous pores were used to characterize MB adsorption at infinite dilution and at higher concentration in which MB dimers and one-dimensional chains formed. Results for the neutral MB molecules are compared with those for the corresponding cation. Simulations of the aqueous pore indicate preferred adsorption on the hydrophobic siloxane surface, while charge-balancing chloride ions adsorb at the aluminol surface. At infinite dilution and in the gas-phase models, MB adsorbs with its primary molecular plane parallel to the siloxane surface to enhance hydrophobic interactions. Sandwiched dimers and chains are oriented perpendicular to the surface to facilitate the strong hydrophobic intermolecular interactions. Compared with quantum results, the hybrid force field predicts a weaker MB adsorption energy but a stronger dimerization energy. The structure and energetics of adsorbed MB at infinite dilution are consistent with the gas-phase binding results, which indicate that monomer adsorption is driven by strong interfacial forces rather than by the hydration properties of the dye. These results inform spectroscopic studies of MB adsorption on mineral surfaces while also revealing critical areas for development of improved hybrid force fields.
High temperature solid state sodium (23Na) magic angle spinning (MAS) NMR spin lattice relaxation times (T1) were evaluated for a series of NASICON (Na3Zr2PS12O12) materials to directly determine Na jump rates. Simulations of the Ti temperature variations that incorporated distributions in Na jump activation energies, or distribution of jump rates, improved the agreement with experiment. The 23Na NMR T1 relaxation results revealed that distributions in the Na dynamics were present for all of the NASICON materials investigated here. The 23Na relaxation experiments also showed that small differences in material composition and/or changes in the processing conditions impacted the distributions in the Na dynamics. The extent of the distribution was related to the presence of a disordered or glassy phosphate phase present in these different sol-gel processed materials. The 23Na NMR T1 relaxation experiments are a powerful tool to directly probing Na jump dynamics and provide additional molecular level details that could impact transport phenomena.