Publications

Ilgen, Anastasia G.; Stewart, Thomas A.; Cygan, Randall T.; Dewers, Thomas D.; Aman, Michael A.; Espinoza, Nicolas E.

Abstract not provided.

Journal of Physical Chemistry C

Greathouse, Jeffery A.; Cygan, Randall T.; Fredrich, J.T.; Jerauld, G.R.

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.

International Journal of Greenhouse Gas Control

Ilgen, A.G.; Cygan, Randall T.

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%.

Journal of Physical Chemistry C

Cygan, Randall T.; Daemen, Luke L.; Ilgen, Anastasia G.; Krumhansl, James L.; Nenoff, T.M.

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.

Teich-McGoldrick, Stephanie T.; Cygan, Randall T.; Cygan, Randall T.; Hart, David B.

Abstract not provided.

Sandia journal manuscript; Not yet accepted for publication

Gordon, Margaret E.; Cygan, Randall T.; Teich-McGoldrick, Stephanie T.; Rodriguez, Mark A.; Meserole, Stephen M.

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.

Cygan, Randall T.; Tenney, Craig M.

Abstract not provided.

Clays and Clay Minerals

Greathouse, Jeffery A.; Geatches, Dawn L.; Pike, Darin Q.; Greenwell, H.C.; Johnston, Cliff T.; Wilcox, Jennifer; Cygan, Randall T.

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.

Greathouse, Jeffery A.; Hart, David B.; Bowers, Geoffrey M.; Kirkpatrick, R.J.; Cygan, Randall T.

Abstract not provided.

Ilgen, Anastasia G.; Cygan, Randall T.; Stewart, Thomas A.; Dewers, Thomas D.

Abstract not provided.

Alam, Todd M.; Wheeler, Jill S.; Spoerke, Erik D.; Cygan, Randall T.; Bell, Nelson S.; Ingersoll, David I.

Abstract not provided.

Materials Research Society Symposium Proceedings

Alam, Todd M.; Bell, Nelson S.; Wheeler, Jill; Spoerke, Erik D.; Cygan, Randall T.; Ingersoll, David I.

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.

Journal of Physical Chemistry A

Zaunbrecher, Laura K.; Teich-McGoldrick, Stephanie T.; Cygan, Randall T.; Elliott, W.C.

Abstract not provided.

Environmental Science and Technology

Ilgen, Anastasia G.; Cygan, Randall T.

Abstract not provided.

Spoerke, Erik D.; Bell, Nelson S.; Edney, Cynthia E.; Wheeler, Jill S.; Cygan, Randall T.; Alam, Todd M.; Ihlefeld, Jon I.; Clem, Paul G.; Ingersoll, David I.

Abstract not provided.

Journal of Physical Chemistry C

Zeitler, Todd R.; Greathouse, Jeffery A.; Gale, Julian D.; Cygan, Randall T.

We introduce a nonbonded three-body harmonic potential energy term for Mg-O-H interactions for improved edge surface stability in molecular simulations. The new potential term is compatible with the Clayff force field and is applied here to brucite, a layered magnesium hydroxide mineral. Comparisons of normal mode frequencies from classical and density functional theory calculations are used to verify a suitable spring constant (k parameter) for the Mg-O-H bending motion. Vibrational analysis of hydroxyl librations at two brucite surfaces indicates that surface Mg-O-H modes are shifted to frequencies lower than the corresponding bulk modes. A comparison of DFT and classical normal modes validates this new potential term. The methodology for parameter development can be applied to other clay mineral components (e.g., Al, Si) to improve the modeling of edge surface stability, resulting in expanded applicability to clay mineral applications. © 2014 American Chemical Society.

Cygan, Randall T.; Bishop, Joseph E.; Martinez, Mario J.; Rinehart, Alex R.; Yoon, Hongkyu Y.; Dewers, Thomas D.; Tenney, Craig M.; Matteo, Edward N.; Ilgen, Anastasia G.; Arnold, Bill W.; Newell, Pania N.

Abstract not provided.

Molecular Simulation

Teich-McGoldrick, Stephanie T.; Greathouse, Jeffery A.; Cygan, Randall T.

Anthropogenic activities have led to an increased concentration of uranium on the Earth’s surface and potentially in the subsurface with the development of nuclear waste repositories. Uranium is soluble in groundwater, and its mobility is strongly affected by the presence of clay minerals in soils and in subsurface sediments. We use molecular dynamics simulations to probe the adsorption of aqueous uranyl (UO₂²⁺) ions onto the basal surface of muscovite, a suitable proxy for typically ultrafine-grained clay phases. Model systems include the competitive adsorption between potassium counterions and aqueous ions (0.1 M and 1.0 M UO₂Cl₂ , 0.1 M NaCl). We find that for systems with potassium and uranyl ions present, potassium ions dominate the adsorption phenomenon. Potassium ions adsorb entirely as inner-sphere complexes associated with the ditrigonal cavity of the basal surface. Uranyl ions adsorb in two configurations when it is the only ion species present, and in a single configuration in the presence of potassium. Finally, the majority of adsorbed uranyl ions are tilted less than 45° relative to the muscovite surface, and are associated with the Si₄Al₂ rings near aluminum substitution sites.

Environmental Science and Technology

Tenney, Craig M.; Cygan, Randall T.

Capture and subsequent geologic storage of CO2 in deep brine reservoirs plays a significant role in plans to reduce atmospheric carbon emission and resulting global climate change. The interaction of CO2 and brine species with mineral surfaces controls the ultimate fate of injected CO2 at the nanoscale via geochemistry, at the pore-scale via capillary trapping, and at the field-scale via relative permeability. We used large-scale molecular dynamics simulations to study the behavior of supercritical CO2 and aqueous fluids on both the hydrophilic and hydrophobic basal surfaces of kaolinite, a common clay mineral. In the presence of a bulk aqueous phase, supercritical CO2 forms a nonwetting droplet above the hydrophilic surface of kaolinite. This CO2 droplet is separated from the mineral surface by distinct layers of water, which prevent the CO2 droplet from interacting directly with the mineral surface. Conversely, both CO2 and H2O molecules interact directly with the hydrophobic surface of kaolinite. In the presence of bulk supercritical CO2, nonwetting aqueous droplets interact with the hydrophobic surface of kaolinite via a mixture of adsorbed CO2 and H2O molecules. Because nucleation and precipitation of minerals should depend strongly on the local distribution of CO2, H2O, and ion species, these nanoscale surface interactions are expected to influence long-term mineralization of injected carbon dioxide. © 2014 American Chemical Society.

The Journal of Physical Chemistry

Cygan, Randall T.; Greathouse, Jeffery A.

Abstract not provided.

Journal of Physical Chemistry. C

Tenney, Craig M.; Cygan, Randall T.

Graph theoretic tools were used to identify and classify clusters of ions and solvent molecules in molecular dynamics simulations of lithium-ion battery electrolytes. Electrolytes composed of various concentrations of LiPF6 dissolved in ethylene carbonate (EC), dimethylene carbonate (DMC), or a 1:1 EC/DMC mixture were simulated at multiple temperatures using classical molecular dynamics. Contrary to Nernst–Einstein theory but consistent with experiment, pure DMC systems had the greatest diffusivity but the lowest conductivity. This disagreement with Nernst–Einstein theory is explained by the observed clustering behavior, which found that systems with pure EC as a solvent formed ion clusters with nonzero charge, whereas systems with pure DMC as a solvent formed primarily neutral clusters.

Elements

Cygan, Randall T.

Abstract not provided.

Cygan, Randall T.; Spoerke, Erik D.

Abstract not provided.

Greathouse, Jeffery A.; Cygan, Randall T.

Abstract not provided.

Greathouse, Jeffery A.; Teich-McGoldrick, Stephanie T.; Jove Colon, Carlos F.; Cygan, Randall T.

Abstract not provided.