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.