Publications

Proposed for publication in the Materials Research Society Bulletin.

Cygan, Randall T.; Bunker, B.C.; Brinker, C.J.; Orendorff, Christopher O.; Nyman, M.; Leung, Kevin L.; Rempe, Susan R.; Huber, Dale L.

Abstract not provided.

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

Abstract not provided.

Journal of Physical Chemistry C

Vasconcelos, Igor F.; Bunker, Bruce A.; Cygan, Randall T.

Molecular dynamics simulation is used to study the mechanisms involved in the adsorption of various ions to the basal surfaces of kaolinite. Analysis of simulation data indicates that cations and anions adsorb preferably on the siloxane and gibbsite surfaces of kaolinite, respectively. Strong inner-sphere adsorption of chlorine at aluminum vacancies on the gibbsite surface and the occurrence of chlorine-driven inner-sphere adsorption of cesium and sodium on the gibbsite surface for high ionic strengths are observed. Cesium ions form strong inner-sphere complexes at ditrigonal cavities on the siloxane surface. Outer-sphere cesium is highly mobile and only weak adsorption may occur. A small amount of sodium adsorbs on the siloxane surface as innersphere complexes at less clearly defined sites. Like cesium, sodium only forms very weak outer-sphere complexes on this surface. Inner-sphere complexes of cadmium and lead do not occur on either surface. Relatively strong outer-sphere cadmium and lead complexes are present on the siloxane surface at ditrigonal cavities. © 2007 American Chemical Society.

Physical Chemistry Chemical Physics

Larentzos, James P.; Greathouse, Jeffery A.; Cygan, Randall T.

Abstract not provided.

Sabo, Dubravko S.; Greathouse, Jeffery A.; Leung, Kevin L.; Cygan, Randall T.; Alam, Todd M.; Varma, Sameer V.; Martin, Marcus G.

Hydrogen storage is a key enabling technology required for attaining a hydrogen-based economy. Fundamental research can reveal the underlying principles controlling hydrogen uptake and release by storage materials, and also aid in characterizing and designing novel storage materials. New ideas for hydrogen storage materials come from exploiting the properties of hydrophobic hydration, which refers to water s ability to stabilize, by its mode of association, specific structures under specific conditions. Although hydrogen was always considered too small to support the formation of solid clathrate hydrate structures, exciting new experiments show that water traps hydrogen molecules at conditions of low temperatures and moderate pressures. Hydrogen release is accomplished by simple warming. While these experiments lend credibility to the idea that water could form an environmentally attractive alternative storage compound for hydrogen fuel, which would advance our nation s goals of attaining a hydrogen-based economy, much work is yet required to understand and realize the full potential of clathrate hydrates for hydrogen storage. Here we undertake theoretical studies of hydrogen in water to establish a firm foundation for predictive work on clathrate hydrate H{sub 2} storage capabilities. Using molecular simulation and statistical mechanical theories based in part on quantum mechanical descriptions of molecular interactions, we characterize the interactions between hydrogen and liquid water in terms of structural and thermodynamic properties. In the process we validate classical force field models of hydrogen in water and discover new features of hydrophobic hydration that impact problems in both energy technology and biology. Finally, we predict hydrogen occupancy in the small and large cages of hydrogen clathrate hydrates, a property unresolved by previous experimental and theoretical work.

Proposed for publication in the Vadose Zone Journal.

Cygan, Randall T.; Siegel, Malcolm D.

Abstract not provided.

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

Abstract not provided.

Journal of Physical Chemistry B

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

Molecular dynamics simulations were performed on methane clathrate hydrates at ambient conditions. Thermal expansion results over the temperature range 60-300 K show that the unit cell volume increases with temperature in agreement with experiment. Power spectra were obtained at 273 K from velocity autocorrelation functions for selected atoms, and normal modes were assigned. The spectra were further classified according to individual atom types, allowing the assignment of contributions from methane molecules located in small and large cages within the structure I unit cell. The symmetric C-H stretch of methane in the small cages occurs at a higher frequency than for methane located in the large cages, with a peak separation of 14 cm-1. Additionally, we determined that the symmetric C-H stretch in methane gas occurs at the same frequency as methane in the large cages. Results of molecular dynamics simulations indicate the use of power spectra obtained from the velocity autocorrelation function is a reliable method to investigate the vibrational behavior of guest molecules in clathrate hydrates. © 2006 American Chemical Society.

Cygan, Randall T.

Cation binding by polysaccharides is observed in many environments and is important for predictive environmental modeling, and numerous industrial and food technology applications. The complexities of these organo-cation interactions are well suited to predictive molecular modeling studies for investigating the roles of conformation and configuration of polysaccharides on cation binding. In this study, alginic acid was chosen as a model polymer and representative disaccharide and polysaccharide subunits were modeled. The ability of disaccharide subunits to bind calcium and to associate with the surface of calcite was investigated. The findings were extended to modeling polymer interactions with calcium ions.

Cygan, Randall T.; Altman, Susan J.

Abstract not provided.

Cygan, Randall T.; Greathouse, Jeffery A.

Water is the critical natural resource of the new century. Significant improvements in traditional water treatment processes require novel approaches based on a fundamental understanding of nanoscale and atomic interactions at interfaces between aqueous solution and materials. To better understand these critical issues and to promote an open dialog among leading international experts in water-related specialties, Sandia National Laboratories sponsored a workshop on April 24-26, 2005 in Santa Fe, New Mexico. The ''Frontiers of Interfacial Water Research Workshop'' provided attendees with a critical review of water technologies and emphasized the new advances in surface and interfacial microscopy, spectroscopy, diffraction, and computer simulation needed for the development of new materials for water treatment.

Journal of Physical Chemistry B

Wang, Jianwei; Kalinichev, Andrey G.; Kirkpatrick, R.J.; Cygan, Randall T.

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 H 2O 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 H 2O 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 1 and earlier Monte Carlo simulations. 2 Instead, a more structurally and orientationally restricted behavior of surface H 2O 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 H 2O 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. 3 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 H 2O molecules along the surface. Calculated energies of water adsorption are in good agreement with available experimental data. 4 © 2005 American Chemical Society.

Greathouse, Jeffery A.; Cygan, Randall T.

Abstract not provided.

Proposed for publication in Physical Chemistry Chemical Physics.

Greathouse, Jeffery A.; Cygan, Randall T.

We used molecular dynamics simulations to study the adsorption of aqueous uranyl species (UO22+) onto clay mineral surfaces in the presence of sodium counterions and carbonato ligands. The large system size (10 000 atoms) and long simulation times (10 ns) allowed us to investigate the thermodynamics of ion adsorption, and the atomistic detail provided clues for the observed adsorption behavior. The model system consisted of the basal surface of a low-charge Na-montmorillonite clay in contact with aqueous uranyl carbonate solutions with concentrations of 0.027 M, 0.081 M, and 0.162 M. Periodic boundary conditions were used in the simulations to better represent an aqueous solution interacting with an external clay surface. Uranyl adsorption tendency was found to decrease as the aqueous uranyl carbonate concentration was increased, while sodium adsorption remained constant. The observed behavior is explained by physical and chemical effects. As the ionic strength of the aqueous solution was increased, electrostatic factors prevented further uranyl adsorption once the surface charge had been neutralized. Additionally, the formation of aqueous uranyl carbonate complexes, including uranyl carbonato oligomers, contributed to the decreased uranyl adsorption tendency.

Voigt, James A.; Hsu, Julia W.; Criscenti, Louise C.; Cygan, Randall T.; McKenzie, Bonnie B.

Abstract not provided.

Environmental Science and Technology

Altman, Susan J.; Rivers, Mark L.; Reno, Marissa D.; Cygan, Randall T.; Mclain, Angela A.

Synchrotron-source X-ray computerized microtomography (CMT) was used to evaluate the adsorptive properties of aggregate soil samples. A linear relationship between measured mean mass attenuation coefficient (σ) and mass fraction iron was generated by imaging mineral standards with known iron contents. On the basis of reported stoichiometries of the clay minerals and identifications of iron oxyhydroxides (1), we calculated the mass fraction iron and iron oxyhydroxide in the intergranular material. The mass fractions of iron were estimated to range from 0.17 to 0.22 for measurements made at 18 keV and from 0.18 to 0.21 for measurements made at 26 keV. One aggregate sample also contained regions within the intergranular material with mass fraction iron ranging from 0.29 to 0.31 and from 0.33 to 0.36 for the 18 and 26 keV measurements, respectively. The mass fraction iron oxyhydroxide ranged from 0.18 to 0.35 for the low-iron intergranular material and from 0.40 to 0.59 for the high-iron intergranular material. Using absorption edge difference imaging with CMT, we visualized cesium on the intergranular material, presumably because of adsorption and possible exchange reactions. By characterizing the mass fraction iron, the mass fraction iron oxyhydroxide, and the adsorptive capacity of these soil mineral aggregates, we provide information useful for conceptualization, development, and parametrization of transport models. © 2005 American Chemical Society.

Proposed for publication in Geochemical News.

Cygan, Randall T.; Brady, Patrick V.; Siegel, Malcolm D.; Jove Colon, Carlos F.; Bryan, Charles R.

Abstract not provided.

Criscenti, Louise C.; Spoerke, Erik D.; McKenzie, Bonnie B.; Cygan, Randall T.; Voigt, James A.

Solution-based synthesis is a powerful approach for creating nano-structured materials. Although there have been significant recent successes in its application to fabricating nanomaterials, the general principles that control solution synthesis are not well understood. The purpose of this LDRD project was to develop the scientific principles required to design and build unique nanostructures in crystalline oxides and II/VI semiconductors using solution-based molecular self-assembly techniques. The ability to synthesize these materials in a range of different nano-architectures (from controlled morphology nanocrystals to surface templated 3-D structures) has provided the foundation for new opportunities in such areas as interactive interfaces for optics, electronics, and sensors. The homogeneous precipitation of ZnO in aqueous solution was used primarily as the model system for the project. We developed a low temperature, aqueous solution synthesis route for preparation of large arrays of oriented ZnO nanostructures. Through control of heterogeneous nucleation and growth, methods to predicatively alter the ZnO microstructures by tailoring the surface chemistry of the crystals were established. Molecular mechanics simulations, involving single point energy calculations and full geometry optimizations, were developed to assist in selecting appropriate chemical systems and understanding physical adsorption and ultimately growth mechanisms in the design of oxide nanoarrays. The versatility of peptide chemistry in controlling the formation of cadmium sulfide nanoparticles and zinc oxide/cadmium sulfide heterostructures was also demonstrated.

Acta Crystallographica Section B

Cygan, Randall T.

Abstract not provided.

Proposed for publication in Physics and Chemistry of Minerals.

Cygan, Randall T.

An empirical model for investigating the behavior of CaCO{sub 3} polymorphs incorporating a shell model for oxygen has been created. The model was constructed by fitting to: the structure of aragonite and calcite; their elastic, static and high-frequency dielectric constants; phonon frequencies at the wave vectors [1/2 0 2] and [0 0 0] of calcite; and vibrational frequencies of the carbonate deformation modes of calcite. The high-pressure phase transition between calcite I and II is observed. The potentials for the CO{sub 3} group were transferred to other carbonates, by refitting the interaction between CO{sub 3} and the cation to both the experimental structures and their bulk modulus, creating a set of potentials for calculating the properties of a wide range of carbonate materials. Defect energies of substitutional cation defects were analyzed for calcite and aragonite phases. The results were rationalized by studying the structure of calcite and aragonite in greater detail.

Cygan, Randall T.

Bulk and surface energies are calculated for endmembers of the isostructural rhombohedral carbonate mineral family, including Ca, Cd, Co, Fe, Mg, Mn, Ni, and Zn compositions. The calculations for the bulk agree with the densities, bond distances, bond angles, and lattice enthalpies reported in the literature. The calculated energies also correlate with measured dissolution rates: the lattice energies show a log-linear relationship to the macroscopic dissolution rates at circumneutral pH. Moreover, the energies of ion pairs translated along surface steps are calculated and found to predict experimentally observed microscopic step retreat velocities. Finally, pit formation excess energies decrease with increasing pit size, which is consistent with the nonlinear dissolution kinetics hypothesized for the initial stages of pit formation.

Langmuir

Duckworth, Owen W.; Cygan, Randall T.; Martin, Scot T.

Bulk and surface energies are calculated for endmembers of the isostructural rhombohedral carbonate mineral family, including Ca, Cd, Co, Fe, Mg, Mn, Ni, and Zn compositions. The calculations for the bulk agree with the densities, bond distances, bond angles, and lattice enthalpies reported in the literature. The calculated energies also correlate with measured dissolution rates: the lattice energies show a log-linear relationship to the macroscopic dissolution rates at circumneutral pH. Moreover, the energies of ion pairs translated along surface steps are calculated and found to predict experimentally observed microscopic step retreat velocities. Finally, pit formation excess energies decrease with increasing pit size, which is consistent with the nonlinear dissolution kinetics hypothesized for the initial stages of pit formation.

Proposed for publication in Chemistry of Materials.

Thoma, Steven T.; Thoma, Steven T.; Bonhomme, F.; Cygan, Randall T.

A novel layered manganese(II) phosphate, Mn{sub 3}(PO{sub 4}){sub 4} {center_dot} 2(H{sub 3}NCH{sub 2}CH{sub 2}){sub 3}N {center_dot} 6(H{sub 2}O), has been synthesized solvothermally using tris(2-aminoethyl)amine (TREN) as a template. The structure was solved ab initio using X-ray powder diffraction data and confirmed by molecular modeling. The compound was further characterized by SEM, IR spectroscopy, photoluminescence, and elemental and thermal analysis. The compound crystallizes in the trigonal space group P{sub 3}c1 with a = 8.8706(4) {angstrom}, c = 26.158(2) {angstrom}, and V = 1782.6(2) {angstrom}{sup 3}. The structure consists of layers of corner sharing Mn(II)O{sub 4} and PO{sub 4} tetrahedra forming infinite [Mn{sub 3}(PO{sub 4}){sub 4}]{sup 6-} macroanions with 4.6 net topology, sandwiched by layers of TREN and water molecules. The protonated TREN molecules provide charge balancing for the inorganic sheets; the interlayer stability is accomplished mainly by a network of hydrogen bonds between water molecules and the inorganic macroanions. This hybrid organic/inorganic layered material can be reversibly dehydrated.

Proposed for publication in the Journal of Solid-State Chemistry.

Nyman, M.; Nyman, M.; Criscenti, Louise C.; Bonhomme, F.; Rodriguez, Marko A.; Cygan, Randall T.

Polyoxoniobate chemistry, both in the solid state and in solution is dominated by [Nb{sub 6}O{sub 19}]{sup 8-}, the Lindquist ion. Recently, we have expanded this chemistry through use of hydrothermal synthesis. The current publication illustrates how use of heteroatoms is another means of diversifying polyoxoniobate chemistry. Here we report the synthesis of Na{sub 8}[Nb{sub 8}Ti{sub 2}O{sub 28}] {center_dot} 34H{sub 2}O [{bar 1}] and its structural characterization from single-crystal X-ray data. This salt crystallizes in the P-1 space group (a = 11.829(4) {angstrom}, b = 12.205(4) {angstrom}, c = 12.532(4) {angstrom}, {alpha} = 97.666(5){sup o}, {beta} = 113.840(4){sup o}, {gamma} = 110.809(4){sup o}), and the decameric anionic cluster [Nb{sub 8}Ti{sub 2}O{sub 28}]{sup 8-} has the same cluster geometry as the previously reported [Nb{sub 10}O{sub 28}]{sup 6-} and [V{sub 10}O{sub 28}]{sup 6-}. Molecular modeling studies of [Nb{sub 10}O{sub 28}]{sup 6-} and all possible isomers of [Nb{sub 8}Ti{sub 2}O{sub 28}]{sup 8-} suggest that this cluster geometry is stabilized by incorporating the Ti{sup 4+} into cluster positions in which edge-sharing is maximized. In this manner, the overall repulsion between edge-sharing octahedra within the cluster is minimized, as Ti{sup 4+} is both slightly smaller and of lower charge than Nb{sup 5+}. Synthetic studies also show that while the [Nb{sub 10}O{sub 28}]{sup 6-} cluster is difficult to obtain, the [Nb{sub 8}Ti{sub 2}O{sub 28}]{sup 8-} cluster can be synthesized reproducibly and is stable in neutral to basic solutions, as well.

Proposed for publication in Molecular Simulation.

Criscenti, Louise C.; Criscenti, Louise C.; Cygan, Randall T.; Voigt, James A.

Abstract not provided.