Publications

Physical Chemistry Chemical Physics

Zeitler, Todd Z.; Greathouse, Jeffery A.; Cygan, Randall T.

Abstract not provided.

Journal of Physical Chemistry C

Cygan, Randall T.

Abstract not provided.

Tenney, Craig M.; Criscenti, Louise C.; Cygan, Randall T.

Abstract not provided.

Journal of Physical Chemistry Letters

Greathouse, Jeffery A.; Cygan, Randall T.

Abstract not provided.

Cygan, Randall T.

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, carbon dioxide. Recent experimental studies have demonstrated the efficacy of intercalating carbon dioxide 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 in the interlayer of montmorillonite clay and to help validate the models with experimental observation.

Muller, Richard P.; Schultz, Peter A.; Cygan, Randall T.; Frischknecht, Amalie F.; Larson, Richard S.; Kanouff, Michael P.; Hewson, John C.; Moffat, Harry K.

Abstract not provided.

Cygan, Randall T.

This work compares the sorption and swelling processes associated with CO2-coal and CO2-clay interactions. We investigated the mechanisms of interaction related to CO2 adsortion in micropores, intercalation into sub-micropores, dissolution in solid matrix, the role of water, and the associated changes in reservoir permeability, for applications in CO2 sequestration and enhanced coal bed methane recovery. The structural changes caused by CO2 have been investigated. A high-pressure micro-dilatometer was equipped to investigate the effect of CO2 pressure on the thermoplastic properties of coal. Using an identical dilatometer, Rashid Khan (1985) performed experiments with CO2 that revealed a dramatic reduction in the softening temperature of coal when exposed to high-pressure CO2. A set of experiments was designed for -20+45-mesh samples of Argonne Premium Pocahontas No.3 coal, which is similar in proximate and ultimate analysis to the Lower Kittanning seam coal that Khan used in his experiments. No dramatic decrease in coal softening temperature has been observed in high-pressure CO2 that would corroborate the prior work of Khan. Thus, conventional polymer (or 'geopolymer') theories may not be directly applicable to CO2 interaction with coals. Clays are similar to coals in that they represent abundant geomaterials with well-developed microporous structure. We evaluated the CO2 sequestration potential of clays relative to coals and investigated the factors that affect the sorption capacity, rates, and permanence of CO2 trapping. For the geomaterials comparison studies, we used source clay samples from The Clay Minerals Society. Preliminary results showed that expandable clays have CO2 sorption capacities comparable to those of coal. We analyzed sorption isotherms, XRD, DRIFTS (infrared reflectance spectra at non-ambient conditions), and TGA-MS (thermal gravimetric analysis) data to compare the effects of various factors on CO2 trapping. In montmorillonite, CO2 molecules may remain trapped for several months following several hours of exposure to high pressure (supercritical conditions), high temperature (above boiling point of water) or both. Such trapping is well preserved in either inert gas or the ambient environment and appears to eventually result in carbonate formation. We performed computer simulations of CO2 interaction with free cations (normal modes of CO2 and Na+CO2 were calculated using B3LYP / aug-cc-pVDZ and MP2 / aug-cc-pVDZ methods) and with clay structures containing interlayer cations (MD simulations with Clayff potentials for clay and a modified CO2 potential). Additionally, interaction of CO2 with hydrated Na-montmorillonite was studied using density functional theory with dispersion corrections. The sorption energies and the swelling behavior were investigated. Preliminary modeling results and experimental observations indicate that the presence of water molecules in the interlayer region is necessary for intercalation of CO2. Our preliminary conclusion is that CO2 molecules may intercalate into interlayer region of swelling clay and stay there via coordination to the interlayer cations.

Greathouse, Jeffery A.; Cygan, Randall T.

Abstract not provided.

Greathouse, Jeffery A.; Cygan, Randall T.

The behavior of water confined in porous materials influences macroscopic phenomena such as solute and water mobility, ion exchange, and adsorption. While the properties of bulk water are generally understood, that of nanoconfined water remains an active area of research. We used molecular simulation and inelastic neutron scattering (INS) to investigate the effect of local structure on the vibrational behavior of nanoconfined water. We focus specifically on the nanosized pores found in the 2:1 phyllosilicate clay minerals palygorskite and sepiolite. These are charge neutral, Mg-rich trioctahedral clays with idealized formulas Mg{sub 5}Si{sub 8}O{sub 20} (OH){sub 2} {center_dot} 8H{sub 2}O and Mg{sub 8}Si{sub 12}O{sub 30} (OH){sub 2} {center_dot} 12H{sub 2}O for palygorskite and sepiolite, respectively. The regular pattern of inverted phyllosilicate layers results in narrow channels with effective van der Waals dimensions of 3.61 {angstrom} x 8.59 {angstrom} (palygorskite) and 4.67 {angstrom} x 12.29 {angstrom} (sepiolite). These clay minerals represent a unique opportunity to study water adsorbed at 'inner edge' sites of uncoordinated Mg{sup 2+}. INS spectra taken at 90 K reveal a large shift in the water librational edge between palygorskite (358 cm{sup -1}) and sepiolite (536 cm{sup -1}), indicating less restricted water motion in the smaller-pore palygorskite. The librational edge of the reference sample (ice I{sub h}) is similar to sepiolite, which confirms the unique water behavior in palygorskite. We used both classical molecular dynamics (CMD) simulations and more rigorous density functional theory (DFT) calculations to investigate the hydrogen bonding environment and vibrational behavior of structural water, defined as those water molecules coordinated to Mg{sup 2+} along the pore walls. These waters remain coordinated throughout the 1-ns timescale of the CMD simulations, and the resulting vibrational spectra indicate a similar shift in the water librational edges seen in the INS spectra. The DFT-optimized structures indicate differences in hydrogen bonding between palygorskite and sepiolite, which could explain the librational shift. Corner-sharing silicate tetrahedra in palygorskite are tilted with respect to the crystallographic a-axis due to the induced strain of layer inversion. As a result, only two short (1.9 {angstrom}) hydrogen bonds form between each water and the framework. In contrast, the relatively unstrained sepiolite structure, each water forms three hydrogen bonds with the framework, and at greater distances (2.0 {angstrom} - 2.5 {angstrom}).

Cygan, Randall T.

Abstract not provided.

Hansen, Francis D.; Gaither, Katherine N.; Sobolik, Steven R.; Cygan, Randall T.; Hardin, Ernest H.; Rechard, Robert P.; Freeze, Geoffrey A.; Sassani, David C.; Brady, Patrick V.; Stone, Charles M.; Martinez, Mario J.; Dewers, Thomas D.

This report evaluates the feasibility of high-level radioactive waste disposal in shale within the United States. The U.S. has many possible clay/shale/argillite basins with positive attributes for permanent disposal. Similar geologic formations have been extensively studied by international programs with largely positive results, over significant ranges of the most important material characteristics including permeability, rheology, and sorptive potential. This report is enabled by the advanced work of the international community to establish functional and operational requirements for disposal of a range of waste forms in shale media. We develop scoping performance analyses, based on the applicable features, events, and processes identified by international investigators, to support a generic conclusion regarding post-closure safety. Requisite assumptions for these analyses include waste characteristics, disposal concepts, and important properties of the geologic formation. We then apply lessons learned from Sandia experience on the Waste Isolation Pilot Project and the Yucca Mountain Project to develop a disposal strategy should a shale repository be considered as an alternative disposal pathway in the U.S. Disposal of high-level radioactive waste in suitable shale formations is attractive because the material is essentially impermeable and self-sealing, conditions are chemically reducing, and sorption tends to prevent radionuclide transport. Vertically and laterally extensive shale and clay formations exist in multiple locations in the contiguous 48 states. Thermal-hydrologic-mechanical calculations indicate that temperatures near emplaced waste packages can be maintained below boiling and will decay to within a few degrees of the ambient temperature within a few decades (or longer depending on the waste form). Construction effects, ventilation, and the thermal pulse will lead to clay dehydration and deformation, confined to an excavation disturbed zone within a few meters of the repository, that can be reasonably characterized. Within a few centuries after waste emplacement, overburden pressures will seal fractures, resaturate the dehydrated zones, and provide a repository setting that strongly limits radionuclide movement to diffusive transport. Coupled hydrogeochemical transport calculations indicate maximum extents of radionuclide transport on the order of tens to hundreds of meters, or less, in a million years. Under the conditions modeled, a shale repository could achieve total containment, with no releases to the environment in undisturbed scenarios. The performance analyses described here are based on the assumption that long-term standards for disposal in clay/shale would be identical in the key aspects, to those prescribed for existing repository programs such as Yucca Mountain. This generic repository evaluation for shale is the first developed in the United States. Previous repository considerations have emphasized salt formations and volcanic rock formations. Much of the experience gained from U.S. repository development, such as seal system design, coupled process simulation, and application of performance assessment methodology, is applied here to scoping analyses for a shale repository. A contemporary understanding of clay mineralogy and attendant chemical environments has allowed identification of the appropriate features, events, and processes to be incorporated into the analysis. Advanced multi-physics modeling provides key support for understanding the effects from coupled processes. The results of the assessment show that shale formations provide a technically advanced, scientifically sound disposal option for the U.S.

Altman, Susan J.; Tidwell, Vincent C.; Boney, Craig M.; Ballantine, Marissa D.; Cygan, Randall T.; Ho, Clifford K.

Abstract not provided.

Greathouse, Jeffery A.; Nenoff, T.M.; Cygan, Randall T.

Inelastic neutron scattering, density functional theory, ab initio molecular dynamics, and classical molecular dynamics were used to examine the behavior of nanoconfined water in palygorskite and sepiolite. These complementary methods provide a strong basis to illustrate and correlate the significant differences observed in the spectroscopic signatures of water in two unique clay minerals. Distortions of silicate tetrahedra in the smaller-pore palygorskite exhibit a limited number of hydrogen bonds having relatively short bond lengths. In contrast, without the distorted silicate tetrahedra, an increased number of hydrogen bonds are observed in the larger-pore sepiolite with corresponding longer bond distances. Because there is more hydrogen bonding at the pore interface in sepiolite than in palygorskite, we expect librational modes to have higher overall frequencies (i.e., more restricted rotational motions); experimental neutron scattering data clearly illustrates this shift in spectroscopic signatures. Distortions of the silicate tetrahedra in these minerals effectively disrupts hydrogen bonding patterns at the silicate-water interface, and this has a greater impact on the dynamical behavior of nanoconfined water than the actual size of the pore or the presence of coordinatively-unsaturated magnesium edge sites.

Proposed for publication in the Journal of the American Chemical Society.

Greathouse, Jeffery A.; Cygan, Randall T.; Durkin, Justin S.; Nenoff, T.M.; Ockwig, Nathan O.

Inelastic neutron scattering, density functional theory, ab initio molecular dynamics, and classical molecular dynamics were used to examine the behavior of nanoconfined water in palygorskite and sepiolite. These complementary methods provide a strong basis to illustrate and correlate the significant differences observed in the spectroscopic signatures of water in two unique clay minerals. Distortions of silicate tetrahedra in the smaller-pore palygorskite exhibit a limited number of hydrogen bonds having relatively short bond lengths. However, without the distorted silicate tetrahedra, an increased number of hydrogen bonds are observed in the larger-pore sepiolite with corresponding longer bond distances. Because there is more hydrogen bonding at the pore interface in sepiolite than in palygorskite, we expect librational modes to have higher overall frequencies (i.e., more restricted rotational motions); experimental neutron scattering data clearly illustrates this shift in spectroscopic signatures. It follows that distortions of the silicate tetrahedra in these minerals effectively disrupt hydrogen-bonding patterns at the silicate?water interface, and this has a greater impact on the dynamical behavior of nanoconfined water than the actual size of the pore or the presence of coordinatively unsaturated magnesium edge sites.

Journal of Chemical Physics

Greathouse, Jeffery A.; Durkin, Justin S.; Cygan, Randall T.

Abstract not provided.

Proposed for publication in the Journal of Materials Chemistry.

Greathouse, Jeffery A.; Cygan, Randall T.

The micro- to nano-sized nature of layered materials, particularly characteristic of naturally occurring clay minerals, limits our ability to fully interrogate their atomic dispositions and crystal structures. The low symmetry, multicomponent compositions, defects, and disorder phenomena of clays and related phases necessitate the use of molecular models and modern simulation methods. Computational chemistry tools based on classical force fields and quantum-chemical methods of electronic structure calculations provide a practical approach to evaluate structure and dynamics of the materials on an atomic scale. Combined with classical energy minimization, molecular dynamics, and Monte Carlo techniques, quantum methods provide accurate models of layered materials such as clay minerals, layered double hydroxides, and clay-polymer nanocomposites.

Cygan, Randall T.; Jiang, Ying B.; Alam, Todd M.; Brinker, C.J.; Bunker, B.C.; Leung, Kevin L.; Nenoff, T.M.; Nyman, M.; Ockwig, Nathan O.; Orendorff, Christopher O.; Rempe, Susan R.; Singh, Seema S.; Criscenti, Louise C.; Stevens, Mark J.; Thurmer, Konrad T.; Van Swol, Frank; Varma, Sameer V.; Crozier, Paul C.; Feibelman, Peter J.; Houston, Jack E.; Huber, Dale L.

A molecular-scale interpretation of interfacial processes is often downplayed in the analysis of traditional water treatment methods. However, such an approach is critical for the development of enhanced performance in traditional desalination and water treatments. Water confined between surfaces, within channels, or in pores is ubiquitous in technology and nature. Its physical and chemical properties in such environments are unpredictably different from bulk water. As a result, advances in water desalination and purification methods may be accomplished through an improved analysis of water behavior in these challenging environments using state-of-the-art microscopy, spectroscopy, experimental, and computational methods.

Alam, Todd M.; Huang, Jian Y.; Cygan, Randall T.

Abstract not provided.

Greathouse, Jeffery A.; Cygan, Randall T.

Abstract not provided.

Greathouse, Jeffery A.; Cygan, Randall T.

Abstract not provided.

Greathouse, Jeffery A.; Cygan, Randall T.; Simmons, Blake S.; Dedrick, Daniel E.

Advances are reported in several aspects of clathrate hydrate desalination fundamentals necessary to develop an economical means to produce municipal quantities of potable water from seawater or brackish feedstock. These aspects include the following, (1) advances in defining the most promising systems design based on new types of hydrate guest molecules, (2) selection of optimal multi-phase reactors and separation arrangements, and, (3) applicability of an inert heat exchange fluid to moderate hydrate growth, control the morphology of the solid hydrate material formed, and facilitate separation of hydrate solids from concentrated brine. The rate of R141b hydrate formation was determined and found to depend only on the degree of supercooling. The rate of R141b hydrate formation in the presence of a heat exchange fluid depended on the degree of supercooling according to the same rate equation as pure R141b with secondary dependence on salinity. Experiments demonstrated that a perfluorocarbon heat exchange fluid assisted separation of R141b hydrates from brine. Preliminary experiments using the guest species, difluoromethane, showed that hydrate formation rates were substantial at temperatures up to at least 12 C and demonstrated partial separation of water from brine. We present a detailed molecular picture of the structure and dynamics of R141b guest molecules within water cages, obtained from ab initio calculations, molecular dynamics simulations, and Raman spectroscopy. Density functional theory calculations were used to provide an energetic and molecular orbital description of R141b stability in both large and small cages in a structure II hydrate. Additionally, the hydrate of an isomer, 1,2-dichloro-1-fluoroethane, does not form at ambient conditions because of extensive overlap of electron density between guest and host. Classical molecular dynamics simulations and laboratory trials support the results for the isomer hydrate. Molecular dynamics simulations show that R141b hydrate is stable at temperatures up to 265K, while the isomer hydrate is only stable up to 150K. Despite hydrogen bonding between guest and host, R141b molecules rotated freely within the water cage. The Raman spectrum of R141b in both the pure and hydrate phases was also compared with vibrational analysis from both computational methods. In particular, the frequency of the C-Cl stretch mode (585 cm{sup -1}) undergoes a shift to higher frequency in the hydrate phase. Raman spectra also indicate that this peak undergoes splitting and intensity variation as the temperature is decreased from 4 C to -4 C.

Cygan, Randall T.

Within reactive geochemical transport, several conceptual models exist for simulating sorption processes in the subsurface. Historically, the K{sub D} approach has been the method of choice due to ease of implementation within a reactive transport model and straightforward comparison with experimental data. However, for modeling complex sorption phenomenon (e.g. sorption of radionuclides onto mineral surfaces), this approach does not systematically account for variations in location, time, or chemical conditions, and more sophisticated methods such as a surface complexation model (SCM) must be utilized. It is critical to determine which conceptual model to use; that is, when the material variation becomes important to regulatory decisions. The geochemical transport tool GEOQUIMICO has been developed to assist in this decision-making process. GEOQUIMICO provides a user-friendly framework for comparing the accuracy and performance of sorption conceptual models. The model currently supports the K{sub D} and SCM conceptual models. The code is written in the object-oriented Java programming language to facilitate model development and improve code portability. The basic theory underlying geochemical transport and the sorption conceptual models noted above is presented in this report. Explanations are provided of how these physicochemical processes are instrumented in GEOQUIMICO and a brief verification study comparing GEOQUIMICO results to data found in the literature is given.

Nenoff, T.M.; Criscenti, Louise C.; Cygan, Randall T.

Abstract not provided.

Greathouse, Jeffery A.; Cygan, Randall T.

Abstract not provided.

J. Phys. Chem. C

Leung, Kevin L.; Cygan, Randall T.; Alam, Todd M.; Ockwig, Nathan O.

Abstract not provided.