Using molecular dynamics simulations, we investigate the molecular scale origin of crystal face selectivity when one gibbsite particle attaches to another in water. A comparison of the free energy per unit surface area of particle–particle attachment indicates that particle attachment through edge surfaces, where the edge surfaces are either (1 0 0) or (1 1 0) crystal faces, is more energetically favorable compared to attachment between two basal surfaces (i.e., (0 0 1) crystal faces) or between the basal surface of one particle and the edge surface of another. This result suggests that gibbsite crystals with low basal/edge surface area ratio will preferentially attach through edge surfaces, potentially helping the crystals grow laterally. However, for larger gibbsite particles (high basal/edge surface area ratio) the total free energy, not normalized by surface area, of particle attachment through the basal surfaces is lower (more negative) than attachment through the edge surfaces, indicating that larger gibbsite particles will preferentially aggregate through basal surface attachments. The short-range electrostatic interactions including the interparticle hydrogen bonds from surface –OH groups drive particle attachment, and the dominant contribution to the free energy minimum is enthalpic rather than entropic. However, the enthalpy of basal-edge attachment is significantly offset by the entropy leading to a higher free energy (less negative) compared to that of basal-basal attachment. Study of the free energy for a few imperfect attachments of two particles indicates a higher free energy (i.e., less negative, less stable), compared to a perfect attachment
We present a combined molecular dynamics (MD) simulation and X-ray absorption fine structure (XAFS) spectroscopic investigation of aqueous iron adsorption on nanoconfined amorphous silica surfaces. The simulation models examine the effects of pore size, pH (surface charge), iron valency, and counter-ion (chloride or hydroxide). The simulation methods were validated by comparing the coordination environment of adsorbed iron with coordination numbers and bond lengths derived from XAFS. In the MD models, nanoconfinement effects on local iron coordination were investigated by comparing results for unconfined silica surfaces and in confined domains within 2 nm, 4 nm, and 8 nm pores. Experimentally, coordination environments of iron adsorbed onto mesoporous silica with 4 nm and 8 nm pores at pH 7.5 were investigated. The effect of pH in the MD models was included by simulating Fe(ii) adsorption onto negatively charged SiO2surfaces and Fe(iii) adsorption on neutral surfaces. The simulation results show that iron adsorption depends significantly on silica surface charge, as expected based on electrostatic interactions. Adsorption on a negatively charged surface is an order of magnitude greater than on the neutral surface, and simulated surface coverages are consistent with experimental results. Pore size effects from the MD simulations were most notable in the adsorption of Fe(ii) at deprotonated surface sites (SiO−), but adsorption trends varied with concentration and aqueous Fe speciation. The coordination environment of adsorbed iron varied significantly with the type of anion. Considerable ion pairing with hydroxide anions led to the formation of oligomeric surface complexes and aqueous species, resulting in larger iron hydroxide clusters at higher surface loadings.
This report represents the milestone deliverable M4SF-21SN010309021 “Modeling Activities Related to Waste Form Degradation: Progress Report” that describes the progress of R&D activities of ongoing modeling investigations specifically on nuclear waste glass degradation, Density Functional Theory (DFT) studies on clarkeite structure and stability, and electrochemical modeling of spent nuclear fuel (SNF). These activities are part of the newly-created Waste form Testing, Modeling, and Performance work package at Sandia National Laboratories (SNL). This work package is part of the “Inventory and Waste Form Characteristics and Performance” control account that includes various experimental and modeling activities on nuclear waste degradation conducted at Oak Ridge National Laboratory (ORNL), SNL, Argonne National Laboratory (ANL), and Pacific Northwest National Laboratory (PNNL).
This report represents completion of milestone deliverable M2SF-21SN010309012 “Annual Status Update for OWL and Waste Form Characteristics” that provides an annual update on status of fiscal year (FY 2020) activities for the work package SF-20SN01030901 and is due on January 29, 2021. The Online Waste Library (OWL) has been designed to contain information regarding United States (U.S.) Department of Energy (DOE)-managed (as) high-level waste (DHLW), spent nuclear fuel (SNF), and other wastes that are likely candidates for deep geologic disposal, with links to the current supporting documents for the data (when possible; note that no classified or official-use-only (OUO) data are planned to be included in OWL). There may be up to several hundred different DOE-managed wastes that are likely to require deep geologic disposal. This draft report contains versions of the OWL model architecture for vessel information (Appendix A) and an excerpt from the OWL User’s Guide (Appendix B and SNL 2020), which are for the current OWL Version 2.0 on the Sandia External Collaboration Network (ECN).
This report describes the high-level accomplishments from the Plasma Science and Engineering Grand Challenge LDRD at Sandia National Laboratories. The Laboratory has a need to demonstrate predictive capabilities to model plasma phenomena in order to rapidly accelerate engineering development in several mission areas. The purpose of this Grand Challenge LDRD was to advance the fundamental models, methods, and algorithms along with supporting electrode science foundation to enable a revolutionary shift towards predictive plasma engineering design principles. This project integrated the SNL knowledge base in computer science, plasma physics, materials science, applied mathematics, and relevant application engineering to establish new cross-laboratory collaborations on these topics. As an initial exemplar, this project focused efforts on improving multi-scale modeling capabilities that are utilized to predict the electrical power delivery on large-scale pulsed power accelerators. Specifically, this LDRD was structured into three primary research thrusts that, when integrated, enable complex simulations of these devices: (1) the exploration of multi-scale models describing the desorption of contaminants from pulsed power electrodes, (2) the development of improved algorithms and code technologies to treat the multi-physics phenomena required to predict device performance, and (3) the creation of a rigorous verification and validation infrastructure to evaluate the codes and models across a range of challenge problems. These components were integrated into initial demonstrations of the largest simulations of multi-level vacuum power flow completed to-date, executed on the leading HPC computing machines available in the NNSA complex today. These preliminary studies indicate relevant pulsed power engineering design simulations can now be completed in (of order) several days, a significant improvement over pre-LDRD levels of performance.
This report represents completion of milestone deliverable M2SF-19SNO10309013 "Online Waste Library (OWL) and Waste Forms Characteristics Annual Report" that reports annual status on fiscal year (FY) 2019 activities for the work package SF-19SN01030901 and is due on August 2, 2019. The online waste library (OWL) has been designed to contain information regarding United States (U.S.) Department of Energy (DOE)-managed (as) high-level waste (DHLW), spent nuclear fuel (SNF), and other wastes that are likely candidates for deep geologic disposal, with links to the current supporting documents for the data (when possible; note that no classified or official-use-only (OUO) data are planned to be included in OWL). There may be up to several hundred different DOE-managed wastes that are likely to require deep geologic disposal. This annual report on FY2019 activities includes evaluations of waste form characteristics and waste form performance models, updates to the OWL development, and descriptions of the management processes for the OWL. Updates to the OWL include an updated user's guide, additions to the OWL database content for wastes and waste forms, results of the beta testing and changes implemented from it. Also added are descriptions of the management/control processes for the OWL development, version control, and archiving. These processes have been implemented as part of the full production release of OWL (i.e., OWL Version 1.0), which has been developed on, and will be hosted and managed on, Sandia National Laboratories (SNL) systems. The version control/update processes will be implemented for updates to the OWL in the future. Additionally, another process covering methods for interfacing with the DOE SNF Database (DOE 2007) at Idaho National Laboratory on the numerous entries for DOE-managed SNF (DSNF) has been pushed forward by defining data exchanges and is planned to be implemented sometime in FY2020. The INL database is also sometimes referred to as the Spent Fuel Database or the SFDB, which is the acronym that will be used in this report. Once fully implemented, this integration effort will serve as a template for interfacing with additional databases throughout the DOE complex.
A molecular-scale understanding of the transition between hydration states in clay minerals remains a challenging problem because of the very fast stepwise swelling process observed from X-ray diffraction (XRD) experiments. XRD profile modeling assumes the coexistence of multiple hydration states in a clay sample to fit the experimental XRD pattern obtained under humid conditions. While XRD profile modeling provides a macroscopic understanding of the heterogeneous hydration structure of clay minerals, a microscopic model of the transition between hydration states is still missing. Here, for the first time, we use molecular dynamics simulation to investigate the transition states between a dry interlayer, one-layer hydrate, and two-layer hydrate. We find that the hydrogen bonds that form across the interlayer at the clay particle edge make an important contribution to the energy barrier to interlayer hydration, especially for initial hydration.
Glassy silicates are substantially weaker when in contact with aqueous electrolyte solutions than in vacuum due to chemical interactions with preexisting cracks. To investigate this silicate weakening phenomenon, classical molecular dynamics (MD) simulations of silica fracture were performed using the bond-order based, reactive force field ReaxFF. Four different environmental conditions were investigated: vacuum, water, and two salt solutions (1M NaCl, 1M NaOH) that form relatively acidic and basic solutions, respectively. Any aqueous environment weakens the silica, with NaOH additions resulting in the largest decreases in the effective fracture toughness (eKIC) of silica or the loading rate at which the fracture begins to propagate. The basic solution leads to higher surface deprotonation, narrower radius of curvature of the crack tip, and greater weakening of the silica, compared with the more acidic environment. The results from the two different electrolyte solutions correspond to phenomena observed in experiments and provide a unique atomistic insight into how anions alter the chemical-mechanical fracture response of silica.
Observation of vibrational properties of phyllosilicate edges via a combined molecular modeling and experimental approach was performed. Deuterium exchange was utilized to isolate edge vibrational modes from their internal counterparts. The appearance of a specific peak within the broader D2O band indicates the presence of deuteration on the edge surface, and this peak is confirmed with the simulated spectra. These results are the first to unambiguously identify spectroscopic features of phyllosilicate edge sites.
In this work we investigate the Orowan hypothesis, that decreases in surface energy due to surface adsorbates lead directly to lowered fracture toughness, at an atomic/molecular level. We employ a Lennard-Jones system with a slit crack and an infiltrating fluid, nominally with gold-water properties, and explore steric effects by varying the soft radius of fluid particles and the influence of surface energy/hydrophobicity via the solid–fluid binding energy. Using previously developed methods, we employ the J-integral to quantify the sensitivity of fracture toughness to the influence of the fluid on the crack tip, and exploit dimensionless scaling to discover universal trends in behavior.
A fluid flow in a nanochannel highly depends on the wettability of the channel surface to the fluid. The permeability of the nanochannel is usually very low, largely due to the adhesion of fluid at the solid interfaces. Using molecular dynamics (MD) simulations, we demonstrate that the flow of water in a nanochannel with rough hydrophilic surfaces can be significantly enhanced by the presence of a thin layer of supercritical carbon dioxide (scCO2) at the water-solid interfaces. The thin scCO2 layer acts like an atomistic lubricant that transforms a hydrophilic interface into a super-hydrophobic one and triggers a transition from a stick- to- a slip boundary condition for a nanoscale flow. This work provides an atomistic insight into multicomponent interactions in nanochannels and illustrates that such interactions can be manipulated, if needed, to increase the throughput and energy efficiency of nanofluidic systems.
Fracture toughness of silicates is reduced in aqueous environments due to water-silica interactions at the crack tip. To investigate this effect, classical molecular dynamics simulations using the bond-order-based reactive force field (ReaxFF) were used to simulate silica fracture. The chemical and mechanical aspects were separated by simulating fracture in (a) a vacuum with dynamic loading, (b) an aqueous environment with dynamic loading, and (c) an aqueous environment with static subcritical mechanical loading to track silica dissolution. The addition of water to silica fracture reduced the silica fracture toughness by ~25%, a trend consistent with experimentally reported results. Analysis of Si─O bonds in the process zone and calculations of dissipation energy associated with fracture indicated that water relaxes the entire process zone and not just the surface. Additionally, the crack tip sharpens during fracture in water and an increased number of microscopic propagation events occur. This results in earlier fracture in systems with increasing mechanical loading in aqueous conditions, despite the lack of significant silica dissolution. Therefore, the threshold for Si─O bond breakage has been lowered in the presence of water and the reduction in fracture toughness is due to structural and energetic changes in the silica, rather than specific dissolution events.
Predicting chemical-mechanical fracture initiation and propagation in materials is a critical problem, with broad relevance to a host of geoscience applications including subsurface storage and waste disposal, geothermal energy development, and oil and gas extraction. In this project, we have developed molecular simulation and coarse- graining techniques to obtain an atomistic-level understanding of the chemical- mechanical mechanisms that control subcritical crack propagation in materials under tension and impact the fracture toughness. We have applied these techniques to the fracture of fused quartz in vacuum, in distilled water, and in two salt solutions - 1M NaC1, 1M NaOH - that form relatively acidic and basic solutions respectively. We have also established the capability to conduct double-compression double-cleavage experiments in an environmental chamber to observe material fracture in aqueous solution. Both simulations and experiments indicate that fractures propagate fastest in NaC1 solutions, slower in distilled water, and even slower in air.
The adsorption equilibrium constants of monovalent and divalent cations to material surfaces in aqueous media are central to many technological, natural, and geochemical processes. Cation adsorption-desorption is often proposed to occur in concert with proton transfer on hydroxyl-covered mineral surfaces, but to date this cooperative effect has been inferred indirectly. This work applies density functional theory-based molecular dynamics simulations of explicit liquid water/mineral interfaces to calculate metal ion desorption free energies. Monodentate adsorption of Na+, Mg2+, and Cu2+ on partially deprotonated silica surfaces are considered. Na+ is predicted to be unbound, while Cu2+ exhibits binding free energies to surface SiO- groups that are larger than those of Mg2+. The predicted trends agree with competitive adsorption measurements on fumed silica surfaces. As desorption proceeds, Cu2+ dissociates one of the H2O molecules in its first solvation shell, turning into Cu2+(OH-)(H2O)3, while Mg remains Mg2+(H2O)6. The protonation state of the SiO- group at the initial binding site does not vary monotonically with cation desorption.
Using an atomistic technique consistent with continuum balance laws and drawing on classical fracture mechanics theory, we estimate the resistance to fracture propagation of amorphous silica. We discuss correspondence and deviations from classical linear elastic fracture mechanics theory including size dependence, rigid/floppy modes of deformation, and the effects of surface energy and stress.
Methane (CH4) and carbon dioxide (CO2), the two major components generated from kerogen maturation, are stored dominantly in nanometer-sized pores in shale matrix as (1) a compressed gas, (2) an adsorbed surface species and/or (3) a species dissolved in pore water (H2O). In addition, supercritical CO2 has been proposed as a fracturing fluid for simultaneous enhanced oil/gas recovery (EOR) and carbon sequestration. A mechanistic understanding of CH4-CO2-H2O interactions in shale nanopores is critical for designing effective operational processes. Using molecular simulations, we show that kerogen preferentially retains CO2 over CH4 and that the majority of CO2 either generated during kerogen maturation or injected in EOR will remain trapped in the kerogen matrix. The trapped CO2 may be released only if the reservoir pressure drops below the supercritical CO2 pressure. When water is present in the kerogen matrix, it may block CH4 release. However, the addition of CO2 may enhance CH4 release because CO2 can diffuse through water and exchange for adsorbed methane in the kerogen nanopores.
Hypothesis: Sodium adsorption on silica surfaces depends on the solution counter-ion. Here, we use NaOH solutions to investigate basic environments. Simulations: Sodium adsorption on hydroxylated silica surfaces from NaOH solutions were investigated through molecular dynamics with a dissociative force field, allowing for the development of secondary molecular species. Findings: Across the NaOH concentrations (0.01 M − 1.0 M), ∼50% of the Na+ ions were concentrated in the surface region, developing silica surface charges between − 0.01 C/m2 (0.01 M NaOH) and − 0.76 C/m2 (1.0 M NaOH) due to surface site deprotonation. Five inner-sphere adsorption complexes were identified, including monodentate, bidentate, and tridentate configurations and two additional structures, with Na+ ions coordinated by bridging oxygen and hydroxyl groups or water molecules. Coordination of Na+ ions by bridging oxygen atoms indicates partial or complete incorporation of Na+ ions into the silica surface. Residence time analysis identified that Na+ ions coordinated by bridging oxygen atoms stayed adsorbed onto the surface four times longer than the mono/bi/tridentate species, indicating formation of relatively stable and persistent Na+ ion adsorption structures. Such inner-sphere complexes form only at NaOH concentrations of > 0.5 M. Na+ adsorption and lifetimes have implications for the stability of silica surfaces.
Mechanistic insight into the process of crack growth can be obtained through molecular dynamics (MD) simulations. In this investigation of fracture propagation, a slit crack was introduced into an atomistic amorphous silica model and mode I stress was applied through far-field loading until the crack propagates. Atomic displacements and forces and an Irving–Kirkwood method with a Lagrangian kernel estimator were used to calculate the J-integral of classical fracture mechanics around the crack tip. The resulting fracture toughness (KIC), 0.76 ± 0.16 MPa√m, agrees with experimental values. In addition, the stress fields and dissipation energies around the slit crack indicate the development of an inelastic region ~30Å in diameter. This is one of the first reports of KIC values obtained from up-scaled atomic-level energies and stresses through the J-integral. The application of the ReaxFF classical MD force field in this study provides the basis for future research into crack growth in multicomponent oxides in a variety of environmental conditions.
Kerogen plays a central role in hydrocarbon generation in an oil/gas reservoir. In a subsurface environment, kerogen is constantly subjected to stress confinement or relaxation. The interplay between mechanical deformation and gas adsorption of the materials could be an important process for shale gas production but unfortunately is poorly understood. Using a hybrid Monte Carlo/molecular dynamics simulation, we show here that a strong chemo-mechanical coupling may exist between gas adsorption and mechanical strain of a kerogen matrix. The results indicate that the kerogen volume can expand by up to 5.4% and 11% upon CH4 and CO2 adsorption at 192 atm, respectively. The kerogen volume increases with gas pressure and eventually approaches a plateau as the kerogen becomes saturated. The volume expansion appears to quadratically increase with the amount of gas adsorbed, indicating a critical role of the surface layer of gas adsorbed in the bulk strain of the material. Furthermore, gas uptake is greatly enhanced by kerogen swelling. Swelling also increases the surface area, porosity, and pore size of kerogen. Our results illustrate the dynamic nature of kerogen, thus questioning the validity of the current assumption of a rigid kerogen molecular structure in the estimation of gas-in-place for a shale gas reservoir or gas storage capacity for subsurface carbon sequestration. The coupling between gas adsorption and kerogen matrix deformation should be taken into consideration.
The porosity of clay aggregates is an important property governing chemical reactions and fluid flow in low-permeability geologic formations and clay-based engineered barrier systems. Pore spaces in clays include interlayer and interparticle pores. Under compaction and dewatering, the size and geometry of such pore spaces may vary significantly (sub-nanometer to microns) depending on ambient physical and chemical conditions. Here we report a molecular dynamics simulation method to construct a complex and realistic clay-like nanoparticle aggregate with interparticle pores and grain boundaries. The model structure is then used to investigate the effect of dewatering and water content on micro-porosity of the aggregates. The results suggest that slow dewatering would create more compact aggregates compared to fast dewatering. Furthermore, the amount of water present in the aggregates strongly affects the particle-particle interactions and hence the aggregate structure. Detailed analyses of particle-particle and water-particle interactions provide a molecular-scale view of porosity and texture development of the aggregates. The simulation method developed here may also aid in modeling the synthesis of nanostructured materials through self-assembly of nanoparticles.
The complexation of toxic and/or radioactive ions on to mineral surfaces is an important topic in geochemistry. We apply periodic-boundary-conditions density functional theory (DFT) molecular dynamics simulations to examine the coordination of Pb(II), SeO2-3, and their contact ion pairs to goethite (1 0 1) and (2 1 0) surfaces. The multitude of Pb(II) adsorption sites and possibility of Pb(II)-induced FeOH deprotonation make this a complex problem. At surface sites where Pb(II) is coordinated to three FeO and/or FeOH groups, and with judicious choices of FeOH surface group protonation states, the predicted Fe-Pb distances are in good agreement with EXAFS measurements. Trajectories where Pb(II) is in part coordinated to only two surface Fe-O groups exhibit larger fluctuations in Pb-O distances. Pb(II)/ SeO2-3 contact ion pairs are at least metastable on goethite (2 1 0) surfaces if the SeO2-3has a monodentate Se-O-Fe bond. Our DFT-based molecular dynamics calculations are a prerequisite for calculations of finite temperature equilibrium binding constants of Pb(II) and Pb(II)/ ion pairs to goethite adsorption sites.
Surface energies of silicates influence crack propagation during brittle fracture and decrease with surface relaxation caused by annealing and hydroxylation. Molecular-level simulations are particularly suited for the investigation of surface processes. In this work, classical MD simulations of silica surfaces are performed with two force fields (ClayFF and ReaxFF) to investigate the effect of force field reactivity on surface structure and energy as a function of surface hydroxylation. An unhydroxylated fracture surface energy of 5.1 J/m2 is calculated with the ClayFF force field, and 2.0 J/m2 is calculated for the ReaxFF force field. The ClayFF surface energies are consistent with the experimental results from double cantilever beam fracture tests (4.5 J/m2), whereas ReaxFF underestimated these surface energies. Surface relaxation via annealing and hydroxylation was performed by creating a low-energy equilibrium surface. Annealing condensed neighboring siloxane bonds increased the surface connectivity, and decreased the surface energies by 0.2 J/m2 for ClayFF and 0.8 J/m2 for ReaxFF. Posthydroxylation surface energies decreased further to 4.6 J/m2 with the ClayFF force field and to 0.2 J/m2 with the ReaxFF force field. Experimental equilibrium surface energies are ∼0.35 J/m2, consistent with the ReaxFF force field. Although neither force field was capable of replicating both the fracture and equilibrium surface energies reported from experiment, each was consistent with one of these conditions. Therefore, future computational investigations that rely on accurate surface energy values should consider the surface state of the system and select the appropriate force field.
Ilgen, Anastasia G.; Burgos, William B.; Criscenti, Louise C.; Furukawa, Yoko F.; Jun, Young-Shin J.; Kerisit, Sebastien K.; Kubicki, James K.; Lee, Sang H.; Smith, Frances S.; Stack, Andrew S.; Stillings, Lisa S.
This report describes the progress in fiscal year 2010 in developing the Waste Integrated Performance and Safety Codes (IPSC) in support of the U.S. Department of Energy (DOE) Office of Nuclear Energy Advanced Modeling and Simulation (NEAMS) Campaign. The goal of the Waste IPSC is to develop an integrated suite of computational modeling and simulation capabilities to quantitatively assess the long-term performance of waste forms in the engineered and geologic environments of a radioactive waste storage or disposal system. The Waste IPSC will provide this simulation capability (1) for a range of disposal concepts, waste form types, engineered repository designs, and geologic settings, (2) for a range of time scales and distances, (3) with appropriate consideration of the inherent uncertainties, and (4) in accordance with robust verification, validation, and software quality requirements. Waste IPSC activities in fiscal year 2010 focused on specifying a challenge problem to demonstrate proof of concept, developing a verification and validation plan, and performing an initial gap analyses to identify candidate codes and tools to support the development and integration of the Waste IPSC. The current Waste IPSC strategy is to acquire and integrate the necessary Waste IPSC capabilities wherever feasible, and develop only those capabilities that cannot be acquired or suitably integrated, verified, or validated. This year-end progress report documents the FY10 status of acquisition, development, and integration of thermal-hydrologic-chemical-mechanical (THCM) code capabilities, frameworks, and enabling tools and infrastructure.
Classical molecular dynamics simulations were used to investigate the formation of water droplets on two kaolinite surfaces: the gibbsite-like surface which is hydrophilic and the silica surface which is hydrophobic. Two methods for calculating contact angles were investigated in detail. The method of Giovambattista et al. was successful in calculating contact angles on both surfaces that compare well to the experimental data available. This is the first time that contact angles have been calculated for kaolinite surfaces from molecular simulations. This preliminary study provides the groundwork for investigating contact angles for more complex systems involving multiple fluids (water, CO{sub 2}, oil) in contact with different minerals in the subsurface environment.
Bacteria, algae and plants produce metal-specific chelators to capture required nutrient or toxic trace metals. Biological systems are thought to be very efficient, honed by evolutionary forces over time. Understanding the approaches used by living organisms to select for specific metals in the environment may lead to design of cheaper and more effective approaches for metal recovery and contaminant-metal remediation. In this study, the binding of a common siderophore, desferrioxamine B (DFO-B), to three aqueous metal cations, Fe(II), Fe(III), and UO{sub 2}(VI) was investigated using classical molecular dynamics. DFO-B has three acetohydroxamate groups and a terminal amine group that all deprotonate with increasing pH. For all three metals, complexes with DFO-B (-2) are the most stable and favored under alkaline conditions. Under more acidic conditions, the metal-DFO complexes involve chelation with both acetohydroxamate and acetylamine groups. The approach taken here allows for detailed investigation of metal binding to biologically-designed organic ligands.
This one-year exploratory LDRD aims to provide fundamental understanding of the mechanism of CO2 scrubbing platforms that will reduce green house gas emission and mitigate the effect of climate change. The project builds on the team members expertise developed in previous LDRD projects to study the capture or preferential retention of CO2 in nanoporous membranes and on metal oxide surfaces. We apply Density Functional Theory and ab initio molecular dynamics techniques to model the binding of CO2 on MgO and CaO (100) surfaces and inside water-filled, amine group functionalized silica nanopores. The results elucidate the mechanisms of CO2 trapping and clarify some confusion in the literature. Our work identifies key future calculations that will have the greatest impact on CO2 capture technologies, and provides guidance to science-based design of platforms that can separate the green house gas CO2 from power plant exhaust or even from the atmosphere. Experimentally, we modify commercial MFI zeolite membranes and find that they preferentially transmit H2 over CO2 by a factor of 34. Since zeolite has potential catalytic capability to crack hydrocarbons into CO2 and H2, this finding paves the way for zeolite membranes that can convert biofuel into H2 and separate the products all in one step.
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.
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.
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.