Publications

Abstract not provided.

Abstract not provided.

Abstract not provided.

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.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Magnesium oxide (MgO)-engineered barriers used in subsurface applications will be exposed to high concentration brine environments and may form stable intermediate phases that can alter the effectiveness of the barrier. To explore the formation of these secondary intermediate phases, MgO was aged in water and three different brine solutions and characterized with X-ray diffraction (XRD) and 1H magic angle spinning (MAS) nuclear magnetic resonance (NMR) spectroscopy. After aging, there is ∼4% molar equivalent of a hydrogen-containing species formed. The 1H MAS NMR spectra resolved multiple minor phases not visible in XRD, indicating that diverse disordered proton-containing environments are present in addition to crystalline Mg(OH)2 brucite. Density functional theory (DFT) simulations for the proposed Mg-O-H-, Mg-Cl-O-H-, and Na-O-H-containing phases were performed to index resonances observed in the experimental 1H MAS NMR spectra. Although the intermediate crystal structures exhibited overlapping 1H NMR resonances in the spectra, Mg-O-H intermediates were attributed to the growth of resonances in the δ +1.0 to 0.0 ppm region, and Mg-Cl-O-H structures produced the increasing contributions of the δ = +2.5 to 5.0 ppm resonances in the chloride-containing brines. Overall, 1H NMR analysis of aged MgO indicates the formation of a wide range of possible intermediate structures that cannot be observed or resolved in the XRD analysis.

Abstract not provided.

Here, we apply density functional theory (DFT) to investigate rare-earth metal organic frameworks (RE-MOFs), RE12(μ3-OH)16(C8O6H4)8(C8O6H5)4 (RE = Y, Eu, Tb, Yb), and characterize the level of theory needed to accurately predict structural and electronic properties in MOF materials with 4f-electrons. A two-step calculation approach of geometry optimization with spin-restricted DFT and large core potential (LCPs), and detailed electronic structures with spin-unrestricted DFT with a full valence potential + Hubbard U correction is investigated. Spin-restricted DFT with LCPs resulted in good agreement between experimental lattice parameters and optimized geometries, while a full valence potential is necessary for accurate representation of the electronic structure. The electronic structure of Eu-DOBDC MOF indicated a strong dependence on the treatment of highly localized 4f-electrons and spin polarization, as well as variation within a range of Hubbard corrections (U = 1-9 eV). For Hubbard corrected spin-unrestricted calculations, a U value of 1-4 eV maintains the non-metallic character of the band gap with slight deviations in f-orbital energetics. When compared with experimentally reported results, the importance of the full valence calculation and the Hubbard correction in correctly predicting the electronic structure is highlighted.

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.

Abstract not provided.

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.

Abstract not provided.

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.

This report summarizes the 2018 fiscal year (FY18) field, laboratory, and modeling work funded by the US Department of Energy Office of Nuclear Energy (DOE-NE) Spent Fuel and Waste Science & Technology (SFWST) campaign as part of the Sandia National Laboratories Salt Research and Development (R&D) and Salt International work packages. This report satisfies level-two milestone M2SF-18SNO10303031and comprises three related but stand-alone sections. The first section summarizes the programmatic progress made to date in the DOE-NE salt program and its goals going forward. The second section presents brine composition modeling and laboratory activities related to salt evaporation experiments, which will be used to interpret data collected during the heater test. The third section presents theoretical and numerical modeling work done to investigate the effects brine composition have on dihedral angle and the permeability of salt.

The impact on the morphology nanoceramic materials generated from group 4 metal alkoxides ([M(OR)4]) and the same precursors modified by 6,6'-(((2-hydroxyethyl)azanediyl)bis(methylene))bis(2,4-di- tert-butylphenol) (referred to as H3-AM-DBP2 (1)) was explored. The products isolated from the 1:1 stoichiometric reaction of a series of [M(OR)4] where M = Ti, Zr, or Hf; OR = OCH(CH3)2(OPr i); OC(CH3)3(OBu t); OCH2C(CH3)3(ONep) with H3-AM-DBP2 proved, by single crystal X-ray diffraction, to be [(ONep)Ti( k4( O,O',O'',N)-AM-DBP2)] (2), [(OR)M(μ( O)- k3( O',O'',N)-AM-DBP2)]2 [M = Zr: OR = OPr i, 3·tol; OBu t, 4·tol; ONep, 5·tol; M = Hf: OR = OBu t, 6·tol; ONep, 7·tol]. The product from each system led to a tetradentate AM-DBP2 ligand and retention of a parent alkoxide ligand. For the monomeric Ti derivative (2), the metal was solved in a trigonal bipyramidal geometry, whereas for the Zr (3-5) and Hf (6, 7) derivatives a symmetric dinuclear complex was formed where the ethoxide moiety of the AM-DBP2 ligand bridges to the other metal center, generating an octahedral geometry. High quality density functional theory level gas-phase electronic structure calculations on compounds 2-7 using Gaussian 09 were used for meaningful time dependent density functional theory calculations in the interpretation of the UV-vis absorbance spectral data on 2-7. Nanoparticles generated from the solvothermal treatment of the ONep/AM-DBP2 modified compounds (2, 5, 7) in comparison to their parent [M(ONep)4] were larger and had improved regularity and dispersion of the final ceramic nanomaterials.

Abstract not provided.

Abstract not provided.

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.

Abstract not provided.

Abstract not provided.

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.