Publications

Abstract not provided.

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.

Abstract not provided.

Molecular simulations are used to assess the ability of metal-organic framework (MOF) materials to store and separate noble gases. Specifically, grand canonical Monte Carlo simulation techniques are used to predict noble gas adsorption isotherms at room temperature. Experimental trends of noble gas inflation curves of a Zn-based material (IRMOF-1) are matched by the simulation results. The simulations also predict that IRMOF-1 selectively adsorbs Xe atoms in Xe/Kr and Xe/Ar mixtures at total feed gas pressures of 1 bar (14.7 psia) and 10 bar (147 psia). Finally, simulations of a copper-based MOF (Cu-BTC) predict this material's ability to selectively adsorb Xe and Kr atoms when present in trace amounts in atmospheric air samples. These preliminary results suggest that Cu-BTC may be an ideal candidate for the pre-concentration of noble gases from air samples. Additional simulations and experiments are needed to determine the saturation limit of Cu-BTC for xenon, and whether any krypton atoms would remain in the Cu-BTC pores upon saturation.

Abstract not provided.

Abstract not provided.

Molecular dynamics simulations were conducted to validate a hybrid force field for metal-organic framework-5 (IRMOF-1). In this force field, only nonbonded parameters are used to describe Zn-O interactions. The CVFF force field was used with slight modifications to describe the benzene dicarboxylate linker. The force field correctly predicts a wide range of structural properties of this MOF, including a negative thermal expansion of approximately 1% at 30 and 293 K, in agreement with both theory and experiment. Compressibility results and the associated elastic moduli are in also good agreement with published density functional theory calculations and nanoindentation experiments. The force field predicts a decrease in elastic moduli as temperature increases, which would greatly affect the mechanical properties of MOFs. Calculated vibrational frequencies for Zn - O modes agree with experiment, and a low-frequency mode representing a 180° rotation of the phenyl groups is seen. This rotation becomes more prevalent as the temperature is increased from 300 to 400 K, in agreement with NMR data. Simulations were also carried out with adsorbed guests, including ethanol, cyclohexane, and several chloromethanes. It is shown that the IRMOF-1 lattice parameter depends on the nature of the guest-framework interaction; strongly hydrophilic guests, such as ethanol, cause a decrease (-0.9%) in unit cell volume, while hydrophobic guests cause an increase (0.7-1.5%) in unit cell volume. The calculated free volumes in IRMOF-1 range from 53.5% to 56.0%, in good agreement with experiment. Finally, the activation energy for benzene self-diffusion calculated at low loadings is in good agreement with previous simulations and NMR results, but the magnitude of the diffusion constant is underestimated, most likely because of deficiencies in the CVFF portion of the force field. The results demonstrate, however, that employing a rigid force field results in much poorer agreement with experimental data. Additionally, a flexible force field approach is required when simulating framework stability because of physical changes or the presence of adsorbates. The use of a general-purpose force field for the organic components allows our approach to be extended to other Zn-based frameworks. © 2008 American Chemical Society.

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.

Abstract not provided.

The structural and vibrational properties of two uncharged 2:1 phyllosilicates, talc and pyrophyllite, are investigated via ab initio and classical molecular dynamics simulations. The quantum mechanical simulations are based on plane-wave pseudopotential density functional theory (DFT), which is shown to be sufficiently accurate in predicting the clay mineral structural and vibrational properties. The classical molecular dynamics (MD) simulations, using the CLAYFF force field, faithfully reproduce the crystal structures with relatively simple analytical functions that include primarily nonbonded interactions. The adsorption properties in these clay minerals are strongly dependent upon the disposition of the hydroxyl group in the octahedral sheet. With the assistance of molecular simulation, the relationship between the hydroxyl group vibrational modes and the molecular-scale structure is explored. The talc hydroxyl groups are oriented perpendicular to the ab plane, while the presence of the dioctahedral vacancies associated with pyrophyllite significantly alters the hydroxyl group structural and vibrational character. Overall, a detailed comparison between the ab initio and the classical MD structural and vibrational properties provides guidance for future refinements to the empirical force field. © 2007 American Chemical Society.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

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.

We report on our studies of the structural properties of a hydrogen molecule dissolved in liquid water. The radial distribution function, coordination number and coordination number distribution are calculated using different representations of the interatomic forces within molecular dynamics (MD), Monte Carlo (MC) and ab initio molecular dynamics (AIMD) simulation frameworks. Although structural details differ in the radial distribution functions generated from the different force fields, all approaches agree that the average and most probable number of water molecules occupying the inner hydration sphere around hydrogen is 16. Furthermore, all results exclude the possibility of clathrate-like organization of water molecules around the hydrophobic molecular hydrogen solute.

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.

Abstract not provided.

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.

Abstract not provided.