Publications

Absolute double differential cross sections (DDCS) of electrons emitted from uracil and 5-bromouracil (BrU) in collisions with protons of energy 200 keV have been measured for various forward and backward emission angles over wide range of electron energies. The measured DDCS are compared with the continuum distorted wave-eikonal initial state (CDW-EIS) calculations. The optimized structure of the BrU was estimated along with the population analysis of all the occupied orbitals using a self-consistent field density. A comparison between the measured DDCS data for the two molecules show that the cross section of low energy electrons emitted from BrU is substantially larger than that for uracil. The BrU-to-uracil DDCS ratios obtained from the present measurements indicate an enhancement of the electron emission by a factor which is as large as 2.0 to 2.5. These electrons being the major agent for damaging the DNA/RNA of the malignant tissues, the present results are expected to provide an important input for the radiosensitization effect in hadron therapy. It is noteworthy to mention that the CDW-EIS calculations for Coulomb ionization cannot predict such enhancement. A large angular asymmetry is observed for uracil with a broad structure, which is absent in case of BrU.

Structural alloys may experience corrosion when exposed to molten chloride salts due to selective dissolution of active alloying elements. One way to prevent this is to make the molten salt reducing. For the KCl + MgCl2 eutectic salt mixture, pure Mg can be added to achieve this. However, Mg can form intermetallic compounds with nickel at high temperatures, which may cause alloy embrittlement. This study shows that an optimum level of excess Mg could be added to the molten salt which will prevent corrosion of alloys like 316 H, while not forming any detectable Ni-Mg intermetallic phases on Ni-rich alloy surfaces.

This report represents completion of milestone deliverable M2SF-22SN010309082 Annual Status Update for OWL, which is due on November 30, 2021 as part of the fiscal year 2022 (FY2022) work package SF-22SN01030908. This report provides an annual update on status of FY2021 activities for the work package “OWL - Inventory – SNL”. 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), DOE-managed spent nuclear fuel (DSNF), and other wastes that are likely candidates for deep geologic disposal. Links to the current supporting documents for the data are provided when possible; however, 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 report contains new information on sodium-bonded spent fuel waste types and wastes forms, which are included in the next release of OWL, Version 3.0, on the Sandia National Laboratories (SNL) External Collaboration Network (ECN). The report also provides an update on the effort to include information regarding the types of vessels capable of disposing of DOE-managed waste.

This project focused on providing a fundamental physico-chemical understanding of the coupling mechanisms of corrosion- and radiation-induced degradation at material-salt interfaces in Ni-based alloys operating in emulated Molten Salt Reactor(MSR) environments through the use of a unique suite of aging experiments, in-situ nanoscale characterization experiments on these materials, and multi-physics computational models. The technical basis and capabilities described in this report bring us a step closer to accelerate the deployment of MSRs by closing knowledge gaps related to materials degradation in harsh environments.

In this work, we have studied the electron emission from one of the polycyclic aromatic hydrocarbon (PAH) molecules namely, fluorene (C13H10), upon 3.5 MeV/u Si8+ ion impact. The experimentally measured absolute double differential cross sections (DDCS) are compared with the continuum distorted wave-eikonal initial state (CDW-EIS) model and the first Born approximation including correct boundary conditions (CB1). The measurements are carried out in the ejected e⁻-energy range of 1 eV–400 eV and in the angular range of 20°–160°. We have obtained the single differential and the total cross sections (TCSs) of e⁻-emission as well. The CB1 calculation largely underestimates the data. The CDW-EIS model, which is applied for the PAH molecule for the first time, provides an overall better agreement with the double differential, single differential and TCS data. The DDCS data for fluorene has also been compared with that for CH4 molecule, at a few angles. The forward–backward angular asymmetry shows a relatively flatter distribution compared to the theoretical predictions. The contribution due to the giant plasmon resonance could not be clearly observed except a mild indication in the asymmetry parameter. The angular distribution of the carbon KLL Auger electron cross section shows certain variations. The study of the KLL hyper-satellite component indicates the double K-ionization cross section is about 8.6% of the single K-ionization one.

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).

Abstract not provided.

Abstract not provided.

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).

The Online Waste Library (OWL) provides a consolidated source of information on Department of Energy-managed radioactive waste likely to require deep geologic disposal. With the release of OWL Version 1.0 in fiscal year 2019 (FY2019), much of the FY2020 work involved developing the OWL change control process and the OWL release process. These two processes (in draft form) were put into use for OWL Version 2.0, which was released in early FY2021. With the knowledge gained, the OWL team refined and documented the two processes in two separate reports. This report focuses on the change control process and discusses the following: (1) definitions and system components; (2) roles and responsibilities; (3) origin of changes; (4) the change control process including the Change List, Task List, activity categories, implementation examples, and checking and review; and (5) the role of the re lease process in ensuring changes in the Change List are incorporated into a public release.

The Online Waste Library (OWL) provides one consolidated source of information on Department of Energy-managed wastes likely to require deep geologic disposal. With the release of OWL Version 1.0 in fiscal year (FY) 2019, much of the FY2020 work involved developing the OWL change control process and the OWL release process. These two processes (in draft form) were put into use for OWL Version 2.0, which was released in early FY2021. With the knowledge gained, the OWL team refined and documented the two processes in two separate reports. This report addresses the release process starting with a definition of release management in Section 2. Section 3 describes the Information Technology Infrastructure Library (ITIL) framework, part of which includes the three different environments used for release management. Section 4 presents the OWL components existing in the different environments and provides details on the release schedule and procedures.

The synthesis, structure, and thermal stability of the periodate double perovskites A₂NaIO₆ (A= Ba, Sr, Ca) were investigated in the context of potential application for the immobilization of radioiodine. A combination of X-ray diffraction and neutron diffraction, Raman spectroscopy, and DFT simulations were applied to determine accurate crystal structures of these compounds and understand their relative stability. The compounds were found to exhibit rock-salt ordering of Na and I on the perovskite B-site; Ba₂NaIO₆ was found to adopt the Fm-3m aristotype structure, whereas Sr₂NaIO₆ and Ca₂NaIO₆ adopt the P2₁/n hettotype structure, characterized by cooperative octahedral tilting. DFT simulations determined the Fm-3m and P2₁/n structures of Ba₂NaIO₆ to be energetically degenerate at room temperature, whereas diffraction and spectroscopy data evidence only the presence of the Fm-3m phase at room temperature, which may imply an incipient phase transition for this compound. The periodate double perovskites were found to exhibit remarkable thermal stability, with Ba₂NaIO₆ only decomposing above 1050 °C in air, which is apparently the highest recorded decomposition temperature so far recorded for any iodine bearing compound. As such, these compounds offer some potential for application in the immobilization of iodine-129, from nuclear fuel reprocessing, with an iodine incorporation rate of 25–40 wt%. The synthesis of these compounds, elaborated here, is also compatible with both current conventional and future advanced processes for iodine recovery from the dissolver off-gas.

We report the measurement of the absolute double differential cross sections (DDCS) of secondary electrons emitted due to the ionization of N2 molecule in collisions with fast electrons having energies between 3 and 5 keV. The emitted electrons with energies from 1-500 eV have been measured for different forward and backward emission angles. The measured DDCS have been compared with the state-of-the-art first Born approximation with correct boundary condition (CB1) model calculations as well as with the classical trajectory Monte Carlo (CTMC) method. From the measured DDCS, the single differential cross sections (SDCS) as a function of the emission energies have been computed and eventually the total ionization cross sections (TCS) have been derived. The TCS values are also compared with a semi-empirical calculation, namely, the CSP-ic (complex scattering potential-ionization contribution) model.

Density functional perturbation theory (DFPT) calculations of the thermodynamic properties of metaschoepite, (UO2)8O2(OH)12·10H2O, are reported. Using a recently revised crystal structure of metaschoepite, the predicted molar entropy and isobaric heat capacity are overall significantly smaller than previous calculations using an earlier orthorhombic crystal structure model. The present DFPT calculations also show large differences between the thermodynamic functions of metaschoepite and schoepite, which might reflect the change in phonon properties upon removal of two H2O molecules per formula unit and alteration of the H-bonded interlayer water network from schoepite to metaschoepite.

We use ab initio spin-polarized density functional theory to study the magnetic order in a Kagomé-like 2D metamaterial consisting of pristine or substitutionally doped phenalenyl radicals polymerized into a nanoporous, graphene-like structure. In this and in a larger class of related structures, the constituent polyaromatic hydrocarbon molecules can be considered as quantum dots that may carry a net magnetic moment. The structure of this porous system and the coupling between the quantum dots may be changed significantly by applying moderate strain, thus allowing to control the magnetic order and the underlying electronic structure.

In radiobiology, predicting the evolution of irradiated biological matter is nowadays an active field of research to identify DNA lesions or to adapt the radiotherapeutic protocols in radiation oncology. In this context, the numerical methods, based on Monte Carlo track-structure simulations, represent the most suitable and powerful tools for understanding the radiobiological damages induced by ionizing particles. In the present work, we report the theoretical differential and total cross sections, computed within the quantum mechanical continuum distorted wave-eikonal initial state (CDW-EIS) approach, for ion impact on water vapor and DNA nucleobases. These cross sections have been used to build up the input database for the homemade Monte Carlo track-structure TILDA-V. A comparison between the theoretical predictions and the available experimental data is presented. Micro-dosimetry results obtained with TILDA-V are also reported.

Classical molecular dynamics (MD) simulations were performed to provide a conceptual understanding of the amorphous-crystalline interface for a candidate negative thermal expansion (NTE) material, ZrW2O8. Simulations of pressure-induced amorphization at 300 K indicate that an amorphous phase forms at pressures of 10 GPa and greater, and this phase persists when the pressure is subsequently decreased to 1 bar. However, the crystalline phase is recovered when the slightly distorted 5 GPa phase is relaxed to 1 bar. Simulations were also performed on a two-phase model consisting of the high-pressure amorphous phase in direct contact with the crystalline phase. Upon equilibration at 300 K and 1 bar, the crystalline phase remains unchanged beyond a thin layer of disrupted structure at the crystalline-amorphous interface. Differences in local atomic structure at the interface are quantified from the simulation trajectories.

This report describes the potential of a novel class of materials--a-ZrW 2 0 8 , Zr 2 WP 2 0 12 , and related compounds that contract upon amorphization as possible radionuclide waste-forms. The proposed ceramic waste-forms would consist of zoned grains, or sintered ceramics with center- loaded radionuclides and barren shells. Radiation-induced amorphization would result in core shrinkage but would not fracture the shells or overgrowths, maintaining isolation of the radionuclide. In this report, we have described synthesis techniques to produce phase-pure forms of the materials, and how to fully densify those materials. Structural models for the materials were developed and validated using DFPT approaches, and radionuclide substitution was evaluated; U(IV), Pu(IV), Tc(IV) and Tc(VII) all readily substitute into the material structures. MD modeling indicated that strain associated with radiation-induced amorphization would not affect the integrity of surrounding crystalline materials, and these results were validated via ion beam experimental studies. Finally, we have evaluated the leach rates of the barren materials, as determined by batch and flow-through reactor experiments. ZrW 2 0 8 leaches rapidly, releasing tungstate while Zr is retained as a solid oxide or hydroxide. Tungsten release rates remain elevated over time and are highly sensitive to contact times, suggesting that this material will not be an effective waste-form. Conversely, tungsten releases rates from Zr2WP2012 rapidly drop, show little dependence on short-term changes in fluid contact time, and in over time, become tied to P release rates. The results presented here suggest that this material may be a viable waste-form for some hard-to-handle radionuclides such as Pu and Tc. ACKNOWLEDGEMENTS The authors acknowledge the contributions to this report from Sandia National Laboratories researchers Steven Meserole, Mark Rodriguez, Clay Payne, Tim Boyle, Nate Padilla, Khalid Hattar, Anthony Monterrosa, Trevor Clark, and Daniel Perry.

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.

The equation of state (EOS) of bulk niobium (Nb) was investigated within the framework of density functional theory, with Mermin's generalization to finite temperatures. The shock Hugoniot for fully-dense and porous Nb was obtained from canonical ab initio molecular dynamics simulations with Erpenbeck's approach based on the Rankine-Hugoniot jump conditions. The phase space was sampled along isotherms between 300 and 4000 K, for densities ranging from ρ=5.5 to 12 g/cm3. Results from simulations compare favorably with room-temperature multianvil and diamond anvil cell data for fully-dense Nb samples and with a recent tabulated SESAME EOS. The results of this study indicate that, for the application of weak and intermediate shocks, the tabular EOS models are expected to give reliable predictions.

The phonon, infrared, and Raman spectroscopic properties of zirconium tungsten phosphate, Zr2(WO4)(PO4)2 (space group Pbcn, IT No. 60; Z = 4), have been extensively investigated using density functional perturbation theory (DFPT) calculations with the Perdew, Burke, and Ernzerhof exchange-correlation functional revised for solids (PBEsol) and validated by experimental characterization of Zr2(WO4)(PO4)2 prepared by hydrothermal synthesis. Using DFPT-simulated infrared, Raman, and phonon density-of-state spectra combined with Fourier transform infrared and Raman measurements, new comprehensive and extensive assignments have been made for the spectra of Zr2(WO4)(PO4)2, resulting in the characterization of its 29 and 34 most intense IR- and Raman-active modes, respectively. DFPT results also reveal that ν1(PO4) symmetric stretching and ν3(PO4) antisymmetric stretching bands have been interchanged in previous Raman experimental assignments. Negative thermal expansion in Zr2(WO4)(PO4)2 appears to have very limited impact on the spectral properties of this compound. This work shows the high accuracy of the PBEsol exchange-correlation functional for studying the spectroscopic properties of crystalline materials using first-principles methods.

Abstract not provided.

Zirconium tetrachloride was synthesized from the reaction between zirconium metal and chlorine gas at 300 °C and was analyzed by electron impact mass spectrometry (EI-MS). Substantial fragmentation products of ZrCl4 were observed in the mass spectra, with ZrCl3 being the most abundant species, followed by ZrCl2, ZrCl, and Zr. The predicted geometry and kinetic stability of the fragments previously mentioned were investigated by density functional theory (DFT) calculations. Energetics of the dissociation processes support the most stable fragment to be ZrCl3 while the least abundant are ZrCl and ZrCl2.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

The relationship between the structure and thermodynamic properties of schoepite, an important uranyl phase with formula [(UO2)8O2(OH)12]·12H2O formed upon corrosion of UO2, has been investigated within the framework of density functional perturbation theory (DFPT). Experimental crystallographic lattice parameters are well reproduced in this study using standard DFT. Phonon calculations within the quasi-harmonic approximation predict standard molar entropy and isobaric heat capacity of S0 = 179.60 J mol-1 K-1 and C0P = 157.4 J mol-1 K-1 at 298.15 K, i.e., ∼6% and ∼4% larger than existing DFPT-D2 calculations. The computed variation of the standard molar isobaric heat capacity with water content from schoepite (UO3·xH2O, x = 2.25) to dehydrated schoepite (x = 1) is predicted to be essentially linear along isotherms ranging from 100 to 500 K. These findings have important implications for the dehydration of layered uranyl corrosion phases and hygroscopic materials.

We use molecular simulations to provide a conceptual understanding of a crystalline-amorphous interface for a candidate negative thermal expansion (NTE) material. Specifically, classical molecular dynamics (MD) simulations were used to investigate the temperature and pressure dependence on structural properties of ZrW2O8. Polarizability of oxygen atoms was included to better account for the electronic charge distribution within the lattice. Constant-pressure simulations of cubic crystalline ZrW2O8 at ambient pressure reveal a slight NTE behavior, characterized by a small structural rearrangement resulting in oxygen sharing between adjacent WO4 tetrahedra. Periodic quantum calculations confirm that the MD-optimized structure is lower in energy than the idealized structure obtained from neutron diffraction experiments. Additionally, simulations of pressure-induced amorphization of ZrW2O8 at 300 K indicate that an amorphous phase forms at pressures greater than 10 GPa, and this phase persists when the pressure is decreased to 1 bar. Simulations were performed on a hybrid model consisting of amorphous ZrW2O8 in direct contact with the cubic crystalline phase. Upon equilibration at 300 K and 1 bar, the crystalline phase remains unchanged beyond a thin layer of disrupted structure at the amorphous interface. Detailed analysis reveals the transition in metal coordination at the interface.

Appropriate waste-forms for radioactive materials must isolate the radionuclides from the environment for long time periods. To accomplish this typically requires low waste-form solubility, to minimize radionuclide release to the environment. However, radiation eventually damages most waste-forms, leading to expansion, crumbling, increased exposed surface area, and faster dissolution. We have evaluated the use of a novel class of materials-ZrW2O8, Zr2P2WO12 and related compounds-that contract upon amorphization. The proposed ceramic waste-forms would consist of zoned grains, or sintered ceramics with center-loaded radionuclides and barren shells. Radiation-induced amorphization would result in core shrinkage but would not fracture the shells or overgrowths, maintaining isolation of the radionuclide. We have synthesized these phases and have evaluated their leach rates. Tungsten forms stable aqueous species at neutral to basic conditions, making it a reliable indicator of phase dissolution. ZrW2O8 leaches rapidly, releasing tungstate while Zr is retained as a solid oxide or hydroxide. Tungsten release rates remain elevated over time and are highly sensitive to contact times, suggesting that this material will not be an effective waste-form. Conversely, tungsten release rates from Zr2P2WO12 rapidly drop and are tied to P release rates; we speculate that a low-solubility protective Zr-phosphate leach layer forms, slowing further dissolution.

We have investigated cubic zirconium tungstate (ZrW2O8) using density functional perturbation theory (DFPT), along with experimental characterization to assess and validate computational results. Cubic zirconium tungstate is among the few known materials exhibiting isotropic negative thermal expansion (NTE) over a broad temperature range, including room temperature where it occurs metastably. Isotropic NTE materials are important for technological applications requiring thermal-expansion compensators in composites designed to have overall zero or adjustable thermal expansion. While cubic zirconium tungstate has attracted considerable attention experimentally, a very few computational studies have been dedicated to this well-known NTE material. Therefore, spectroscopic, mechanical and thermodynamic properties have been derived from DFPT calculations. A systematic comparison of the calculated infrared, Raman, and phonon density-of-state spectra has been made with Fourier transform far-/mid-infrared and Raman data collected in this study, as well as with available inelastic neutron scattering measurements. The thermal evolution of the lattice parameter computed within the quasi-harmonic approximation exhibits negative values below the Debye temperature, consistent with the observed negative thermal expansion characteristics of cubic zirconium tungstate, α-ZrW2O8. These results show that this DFPT approach can be used for studying the spectroscopic, mechanical and thermodynamic properties of prospective NTE ceramic waste forms for encapsulation of radionuclides produced during the nuclear fuel cycle.

The variation of the work function upon carbon adsorption on the reconstructed Au(110) surface is measured experimentally and compared to density functional calculations. The adsorption dynamics is simulated with ab-initio molecular dynamics techniques. The contribution of various energetically available adsorption sites on the deposition process is analyzed, and the work function behavior with carbon coverage is explained by the resultant electron charge density distributions.

This report represents completion of milestone deliverable M2SF-18SNO10309013 "Inventory and Waste Characterization Status Report and OWL Update that reports on FY2018 activities for the work package (WP) SF-18SNO1030901. This report provides the detailed final information for completed FY2018 work activities for WP SF-18SN01030901, and a summary of priorities for FY2019. This status report on FY2018 activities includes evaluations of waste form characteristics and waste form performance models, updates to the OWL development, and descriptions of the two planned 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. There are two processes being planned in FY2018, which will be implemented in FY2019. One process covers methods for interfacing with the DOE SNF DB (DOE 2007) at INL on the numerous entries for DOE managed SNF, and the other process covers the management of updates to, and version control/archiving of, the OWL database. In FY2018, we have pursued three studies to evaluate/redefine waste form characteristics and/or performance models. First characteristic isotopic ratios for various waste forms included in postclosure performance studies are being evaluated to delineate isotope ratio tags that quantitatively identify each particular waste form. This evaluation arose due to questions regarding the relative contributions of radionuclides from disparate waste forms in GDSA results, particularly, radionuclide contributions of DOE-managed SNF vs HLW glass. In our second study we are evaluating the bases of glass waste degradation rate models to the HIP calcine waste form. The HIP calcine may likely be a ceramic matrix material, with multiple ceramic phases with/without a glass phase. The ceramic phases are likely to have different degradation performance from the glass portion. The distribution of radionuclides among those various phases may also be a factor in the radionuclide release rates. Additionally, we have an ongoing investigation of the performance behavior of TRISO particle fuels and are developing a stochastic model for the degradation of those fuels that accounts for simultaneous corrosion of the silicon carbide (SiC) layer and radionuclide diffusion through it. The detailed model of the TRISO particles themselves, will be merged with models of the degradation behavior(s) of the graphite matrix (either prismatic compacts or spherical "pebbles") containing the particles and the hexagonal graphite elements holding the compacts.

Abstract not provided.

This project focused on providing a fundamental mechanistic understanding of the complex degra- dation mechanisms associated with Pellet/Clad Debonding (PCD) through the use of a unique suite of novel synthesis of surrogate spent nuclear fuel, in-situ nanoscale experiments on surrogate interfaces, multi-modeling, and characterization of decommissioned commercial spent fuel. The understanding of a broad class of metal/ceramic interfaces degradation studied within this project provided the technical basis related to the safety of high burn-up fuel, a problem of interest to the DOE.

Abstract not provided.

Cubic zirconium tungstate (α-ZrW2O8), a well-known negative thermal expansion material, has been investigated within the framework of density functional perturbation theory (DFPT), combined with experimental characterization to assess and validate computational results. Using combined Fourier transform infrared measurements and DFPT calculations, new and extensive assignments were made for the far-infrared (<400 cm−1) spectrum of α-ZrW2O8. A systematic comparison of DFPT-simulated infrared, Raman, and phonon density-of-state spectra with Fourier transform far-/mid-infrared and Raman data collected in this study, as well as with available inelastic neutron scattering measurements, shows the superior accuracy of the PBEsol exchange-correlation functional over standard PBE calculations for studying the spectroscopic properties of this material.

The shock Hugoniot for full-density and porous CeO2 was investigated in the liquid regime using ab initio molecular dynamics (AIMD) simulations with Erpenbeck's approach based on the Rankine-Hugoniot jump conditions. The phase space was sampled by carrying out NVT simulations for isotherms between 6000 and 100 000 K and densities ranging from ρ=2.5 to 20g/cm3. The impact of on-site Coulomb interaction corrections +U on the equation of state (EOS) obtained from AIMD simulations was assessed by direct comparison with results from standard density functional theory simulations. Classical molecular dynamics (CMD) simulations were also performed to model atomic-scale shock compression of larger porous CeO2 models. Results from AIMD and CMD compression simulations compare favorably with Z-machine shock data to 525 GPa and gas-gun data to 109 GPa for porous CeO2 samples. Using results from AIMD simulations, an accurate liquid-regime Mie-Grüneisen EOS was built for CeO2. In addition, a revised multiphase SESAME-Type EOS was constrained using AIMD results and experimental data generated in this work. This study demonstrates the necessity of acquiring data in the porous regime to increase the reliability of existing analytical EOS models.

Zirconium tetrachloride, ZrCl4, is a strategic material with wide-ranging applications. Until now, only one crystallographic study on ZrCl4 has been reported [Krebs (1970). Z. Anorg. Allg. Chem.378, 263-272] and that was more than 40 years ago. The compound used for the previous determination was prepared from ZrO2 and Cl2-CCl4, and single-crystal X-ray diffraction (SCXRD) studies on ZrCl4 obtained from Zr metal have not yet been reported. In this context, we prepared ZrCl4 from the reaction of Zr metal and Cl2 gas in a sealed tube and investigated its structure at 100, 150, 200, 250, and 300 K. At 300 K, the SCXRD analysis indicates that ZrCl4 crystallizes in the orthorhombic space group Pca21 [a = 6.262 (9), b = 7.402 (11), c = 12.039 (17) Å, and V = 558.0 (14) Å3] and consists of infinite zigzag chains of edge-sharing ZrCl6 octahedra. This chain motif is similar to that observed previously in ZrCl4, but the structural parameters and space group differ. In the temperature range 100-300 K, no phase transformation was identified, while elongation of intra-chain Zr...Zr [3.950 (1) Å at 100 K and 3.968 (5) Å at 300 K] and inter-chain Cl...Cl [3.630 (3) Å at 100 K and 3.687 (9) Å at 300 K] distances occurred.ZrCl4 has been prepared from the reaction of Zr metal and Cl2 gas in a sealed tube and its structure determined at 100, 150, 200, 250, and 300 K. The chain motif is similar to that observed previously in ZrCl4, but the structural parameters and space group differ.

Abstract not provided.

X-ray radiation induced damage has been known for decades and has largely been viewed as a tremendous nuisance; e.g., most X-ray-related studies of organic and inorganic materials suffer X-ray damage to varying degrees. Although, recent theoretical and experimental investigation of the response of simple chemical systems to X-rays offered better understanding of the mechanistic details of X-ray induced damage, the question about useful applicability of this technique is still unclear. Furthermore we experimentally demonstrate that by tuning pressure and X-ray energy, the radiation induced damage can be controlled and used for synthesis of novel materials.

Abstract not provided.

The crystal structure, lattice dynamics and themomechanical properties of bulk monoclinic zirconium tetrachloride (ZrCl4) have been investigated using zero-damping dispersion-corrected density functional theory [DFT-D3(zero)]. Phonon analysis reveals that ZrCl4(cr) undergoes negative thermal expansion (NTE) near T≈10 K, with a coefficient of thermal expansion of α=-1.2 ppm K−1 and a Grüneisen parameter of γ=-1.1. The bulk modulus is predicted to vary from K0=8.7 to 7.0 GPa in the temperature range 0–550 K. The isobaric molar heat capacity derived from phonon calculations within the quasi-harmonic approximation is in fair agreement with existing calorimetric data.

Molecular structures of kerogen control hydrocarbon production in unconventional reservoirs. Significant progress has been made in developing model representations of various kerogen structures. These models have been widely used for the prediction of gas adsorption and migration in shale matrix. However, using density functional perturbation theory (DFPT) calculations and vibrational spectroscopic measurements, we here show that a large gap may still remain between the existing model representations and actual kerogen structures, therefore calling for new model development. Using DFPT, we calculated Fourier transform infrared (FTIR) spectra for six most widely used kerogen structure models. The computed spectra were then systematically compared to the FTIR absorption spectra collected for kerogen samples isolated from Mancos, Woodford and Marcellus formations representing a wide range of kerogen origin and maturation conditions. Limited agreement between the model predictions and the measurements highlights that the existing kerogen models may still miss some key features in structural representation. A combination of DFPT calculations with spectroscopic measurements may provide a useful diagnostic tool for assessing the adequacy of a proposed structural model as well as for future model development. This approach may eventually help develop comprehensive infrared (IR)-fingerprints for tracing kerogen evolution.

This report provides an update to Sassani et al. (2016) and includes: (1) an updated set of inputs (Sections 2.3) on various additional waste forms (WF) covering both DOE-managed spent nuclear fuel (SNF) and DOE-managed (as) high-level waste (HLW) for use in the inventory represented in the geologic disposal safety analyses (GDSA); (2) summaries of evaluations initiated to refine specific characteristics of particular WF for future use (Section 2.4); (3) updated development status of the Online Waste Library (OWL) database (Section 3.1.2) and an updated user guide to OWL (Section 3.1.3); and (4) status updates (Section 3.2) for the OWL inventory content, data entry checking process, and external OWL BETA testing initiated in fiscal year 2017.

Interaction between polycyclic aromatic hydrocarbon (PAH) molecule and energetic ion is a subject of interest in different areas of modern physics. Here, we present measurements of energy and angular distributions of absolute double differential electron emission cross section for coronene (C₂₄H₁₂) and fluorene (C₁₃H₁₀) molecules under fast bare oxygen ion impact. For coronene, the angular distributions of the low energy electrons are quite different from that of simpler targets like Ne or CH₄, which is not the case for fluorene. The behaviour of the higher electron energy distributions for both the targets are similar to that for simple targets. In case of coronene, a clear signature of plasmon resonance is observed in the analysis of forward-backward angular asymmetry of low energy electron emission. For fluorene, such signature is not identified probably due to lower oscillator strength of plasmon compared to the coronene. The theoretical calculation based on the first-order Born approximation with correct boundary conditions (CB1), in general, reproduced the experimental observations qualitatively, for both the molecules, except in the low energy region for coronene, which again indicates the role of collective excitation. Single differential and total cross sections are also deduced. An overall comparative study is presented.

The structures and properties of Ce_1–xZr_xO₂ (x = 0–1) solid solutions, selected Ce_1–xZr_xO₂ surfaces, and Ce_1–xZr_xO₂/CeO₂ interfaces were computed within the framework of density functional theory corrected for strong electron correlation (DFT+U). The calculated Debye temperature increases steadily with Zr content in (Ce, Zr)O₂ phases, indicating a significant rise in microhardness from CeO₂ to ZrO₂, without appreciable loss in ductility as the interfacial stoichiometry changes. Surface energy calculations for the low-index CeO₂(111) and (110) surfaces show limited sensitivity to strong 4f-electron correlation. The fracture energy of Ce_1–xZr_xO₂(111)/CeO₂(111) increases markedly with Zr content, with a significant decrease in energy for thicker Ce_1–xZr_xO₂ films. These findings suggest the crucial role of Zr acting as a binder at the Ce_1–xZr_xO₂/CeO₂ interfaces, due to the more covalent character of Zr–O bonds compared to Ce–O. Finally, the impact of surface relaxation upon interface cracking was assessed and found to reach a maximum for Ce_0.25Zr_0.75O₂/CeO₂ interfaces.

Abstract: Understanding the radiation-induced effects at the cellular and subcellular levels remains crucial for predicting the evolution of irradiated biological matter. In this context, Monte Carlo track-structure simulations have rapidly emerged among the most suitable and powerful tools. However, most existing Monte Carlo track-structure codes rely heavily on the use of semi-empirical cross sections as well as water as a surrogate for biological matter. In the current work, we report on the up-to-date version of our homemade Monte Carlo code TILDA-V – devoted to the modeling of the slowing-down of 10 keV–100 MeV protons in both water and DNA – where the main collisional processes are described by means of an extensive set of ab initio differential and total cross sections. Graphical abstract: [Figure not available: see fulltext.].

Abstract not provided.

Abstract not provided.

Abstract not provided.

The decoherence of trapped-ion quantum gates due to heating of their motional modes is a fundamental science and engineering problem. This heating is attributed to electric-field noise arising from the trap-electrode surfaces. In this work, we investigate the source of this noise by focusing on the diffusion of carbon-containing adsorbates on the surface of Au(110). We show by density functional theory, based on detailed scanning probe microscopy, how the carbon adatom diffusion on the gold surface changes the energy landscape and how the adatom dipole moment varies with the diffusive motion. A simple model for the diffusion noise, which varies quadratically with the variation of the dipole moment, predicts a noise spectrum, in accordance with the measured values.

Abstract not provided.

Here, the mechanisms and energetics of Zr(0001) surface chlorination by dissociative adsorption of gaseous Cl₂, and associated speciation and surface degradation processes, have been investigated within the framework of density functional theory. Chlorination of Zr(0001) is predicted to be exothermic by ~3 eV/Cl for dissociative adsorption of a single Cl₂ molecule, followed by exothermic chlorination to 1ML and 2 ML under Cl-rich conditions, with respective energy gains of 1.93 and 2.79 eV/Cl. Calculations also show that exfoliation of the top Cl-Zr-Cl sandwich layers is exothermic and most energetically favorable, and can thus be considered as a leading mechanism for Zr(0001) surface dissolution. Consistent with experimental findings, formation of ZrCl₄ molecular products is also found to be dominant during Zr(0001) chlorination.

The mechanisms and energetics of Zr(0001) surface chlorination by dissociative adsorption of gaseous Cl₂, and associated speciation and surface degradation processes, have been investigated within the framework of density functional theory. Chlorination of Zr(0001) is predicted to be exothermic by 3 eV/Cl for dissociative adsorption of a single Cl₂ molecule, followed by exothermic chlorination to 1ML and 2 ML under Cl-rich conditions, with respective energy gains of 1.93 and 2.79 eV/Cl. Calculations also show that exfoliation of the top Cl-Zr-Cl sandwich layers is exothermic and most energetically favorable, and can thus be considered as a leading mechanism for Zr(0001) surface dissolution. Finally, consistent with experimental findings, formation of ZrCl₄ molecular products is also found to be dominant during Zr(0001) chlorination.

Abstract not provided.

The Waste Form Disposal Options Evaluation Report (SNL 2014) evaluated disposal of both Commercial Spent Nuclear Fuel (CSNF) and DOE-managed HLW and Spent Nuclear Fuel (DHLW and DSNF) in the variety of disposal concepts being evaluated within the Used Fuel Disposition Campaign. That work covered a comprehensive inventory and a wide range of disposal concepts. The primary goal of this work is to evaluate the information needs for analyzing disposal solely of a subset of those wastes in a Defense Repository (DRep; i.e., those wastes that are either defense related, or managed by DOE but are not commercial in origin). A potential DRep also appears to be safe in the range of geologic mined repository concepts, but may have different concepts and features because of the very different inventory of waste that would be included. The focus of this status report is to cover the progress made in FY16 toward: (1) developing a preliminary DRep included inventory for engineering/design analyses; (2) assessing the major differences of this included inventory relative to that in other analyzed repository systems and the potential impacts to disposal concepts; (3) designing and developing an on-line waste library (OWL) to manage the information of all those wastes and their waste forms (including CSNF if needed); and (4) constraining post-closure waste form degradation performance for safety assessments of a DRep. In addition, some continuing work is reported on identifying potential candidate waste types/forms to be added to the full list from SNL (2014 – see Table C-1) which also may be added to the OWL in the future. The status for each of these aspects is reported herein.

The interplay between thermodynamics and mechanical properties in the transformation of studtite, (UO2)(O2)(H2O)2·2H2O, into metastudtite, (UO2)(O2)(H2O)2, two important corrosion phases observed on the surface of uranium dioxide exposed to water, is revealed using density functional perturbation theory. Phonon calculations within the quasi-harmonic approximation predict that the standard entropy change for the (UO2)(O2)(H2O)2·2H2O → (UO2)(O2)(H2O)2 + 2H2O reaction is ΔS0 = +80 J·mol-1·K-1 for the production of water in the liquid state and +389 J·mol-1·K-1 for water vapor. Similar to bulk H2O(l), the bulk modulus of (UO2)(O2)(H2O)2·2H2O increases with temperature, contrasting with (UO2)(O2)(H2O)2 which features the typical Anderson-Gruneisen temperature dependence of oxide solids. Upon removal of interstitial H2O in studtite, the most important changes in the shear modulus, the parameter limiting the mechanical stability, arise in the planes normal to chain propagation directions. The present findings have important implications for the dehydration of other hygroscopic materials.

The supercritical carbon dioxide (S - CO2) Brayton Cycle has gained significant attention in the last decade as an advanced power cycle capab le of achieving high efficiency power conversion. Sandia National Laboratories, with support from the U.S. Department of Energy Office of Nuclear Energy (US DOE - NE), has been conducting research and development in order to deliver a technology that is rea dy for commercialization. There are a wide range of materials related challenges that must be overcome for the success of this technology. At Sandia, recent work has focused on the following main areas: (1) Investigating the potential for system cost re duction through the introduction of low cost alloys in low temperature loop sections, (2) Identifying material options for 10MW RCBC systems, (3) Understanding and resolving turbine degradation, (4) Identifying gas foil bearing behavior in CO 2 , and (5) Ide ntifying the influence of gas chemistry on alloy corrosion. Progress in each of these areas is provided in this report.

Abstract not provided.

We report on the experimental and theoretical characterization of a novel GaP polymorph formed by laser heating of a single crystal of GaP-II in its stable region near 43 GPa. Thereby formed unstrained multigrain sample at 43 GPa and 1300 K, allowed high-resolution crystallographic analysis. We find an oS24 as an energetically optimized crystal structure contrary to oS8 reported by Nelmes et al. (1997). Our DFT calculation confirms a stable existence of oS24 between 18 – 50 GPa. The emergence of the oS24 structure is related to the differentiation of phosphorous atoms between those forming P-P dimers and those forming P-Ga bonds only. Bonding anisotropy explains the symmetry lowering with respect to what is generally expected for semiconductors high-pressure polymorphs. The metallization of GaP does not occur through a uniform change of the nature of its bonds but through the formation of an anisotropic phase containing different bond types.

Interactions between CH4, COOH, NH3, OH, SH and armchair (n,n) (n=4,7,14) and zigzag (n,0) (n=7,12,25) single-walled carbon nanotubes (SWCNTs) have been systematically investigated within the framework of dispersion-corrected density functional theory (DFT-D2). Endohedral and exohedral molecular adsorption on SWCNT walls is energetically unfavorable or weak, despite the use of C6/r6 pairwise London-dispersion corrections. The effects of pore size and chirality on the molecule/SWCNTs interaction were also assessed. Chemisorption of COOH, NH3, OH and SH at SWCNT edge sites was examined using a H-capped (7,0) SWCNT fragment and its impact on electrophilic, nucleophilic and radical attacks was predicted by means of Fukui functions.

Here, the ambient temperature equation of state (EoS) of technetium metal has been measured by X-ray diffraction. The metal was compressed using a diamond anvil cell and using a 4:1 methanol-ethanol pressure transmitting medium. The maximum pressure achieved, as determined from the gold pressureEquation of state for technetium from X-ray diffraction and first-principle calculations scale, was 67 GPa. The compression data shows that the HCP phase of technetium is stable up to 67 GPa. The compression curve of technetium was also calculated using first-principles total-energy calculations. Utilizing a number of fitting strategies to compare the experimental and theoretical data it is determined that the Vinet equation of state with an ambient isothermal bulk modulus of B_0T=288 GPa and a first pressure derivative of B'=5.9(2) best represent the compression behavior of technetium metal.

The structure-property relationships of bulk CeO2 and Ce2O3 have been investigated using AM05 and PBEsol exchange-correlation functionals within the frameworks of Hubbard-corrected density functional theory (DFT+U) and density functional perturbation theory (DFPT+U). Compared with conventional PBE+U, RPBE+U, PW91+U and LDA+U functionals, AM05+U and PBEsol+U describe experimental crystalline parameters and properties of CeO2 and Ce2O3 with superior accuracy, especially when +U is chosen close to its value derived by the linear-response approach. The present findings call for a reexamination of some of the problematic oxide materials featuring strong f- and d-electron correlation using AM05+U and PBEsol+U.

The energetics of Sn2+ substitution into the Ca2+ sublattice of hydroxylapatite (HA), Ca10(PO4)6(OH)2, has been investigated within the framework of density functional theory. Calculations reveal that Sn2+ incorporation via coupled substitutions at Ca(ii) sites is energetically favourable up to a composition of Sn6Ca4(PO4)6(OH)2, and further substitutions at Ca(i) sites proceed once full occupancy of Ca(ii) sites by Sn2+ is achieved. Compositions of SnxCa10−x(PO4)6(OH)2 (x = 4-9) are predominant, with an optimal stoichiometry of Sn8Ca2(PO4)6(OH)2, and Sn-substituted HA follows approximately Vegard's law across the entire composition range.

The current work aims at providing an accurate description of the ion track-structure in poly-allyl dyglycol carbonate (PADC) by using an up-to-date Monte-Carlo code-called TILDA-V (a French acronym for Transport d'Ions Lourds Dans l'Aqua & Vivo). In this simulation the ion track-structure in PADC is mainly described in terms of ejected electrons with a particular attention done to the Low Energy Electrons (LEEs). After a brief reminder of the most important channels through which LEEs are prone to break a chemical bond, we will report on the simulated energetic distributions of LEEs along an ion track in PADC for particular incident energies located on both sides of the Bragg-peak position. Finally, based on the rare data dealing with LEEs interaction with polymers or organic molecules, we will emphasise the role played by the LEEs in the formation of a latent track in PADC, and more particularly the one played by the sub-ionization electrons.

Uranyl nitrate is a key species in the nuclear fuel cycle. However, this species is known to exist in different states of hydration, including the hexahydrate ([UO2(NO3)2(H2O)6] often called UNH), the trihydrate [UO2(NO3)2(H2O)3 or UNT], and in very dry environments the dihydrate form [UO2(NO3)2(H2O)2]. Their relative stabilities depend on both water vapor pressure and temperature. In the 1950s and 1960s, the different phases were studied by infrared transmission spectroscopy but were limited both by instrumental resolution and by the ability to prepare the samples for transmission. We have revisited this problem using time-resolved reflectance spectroscopy, which requires no sample preparation and allows dynamic analysis while the sample is exposed to a flow of N2 gas. Samples of known hydration state were prepared and confirmed via X-ray diffraction patterns of known species. In reflectance mode the hexahydrate UO2(NO3)2(H2O)6 has a distinct uranyl asymmetric stretch band at 949.0 cm^-1 that shifts to shorter wavelengths and broadens as the sample desiccates and recrystallizes to the trihydrate, first as a shoulder growing in on the blue edge but ultimately results in a doublet band with reflectance peaks at 966 and 957 cm^-1. The data are consistent with transformation from UNH to UNT as UNT has two inequivalent UO2²⁺ sites. The dehydration of UO2(NO3)2(H2O)6 to UO2(NO3)2(H2O)3 is both a structural and morphological change that has the lustrous lime green UO2(NO3)2(H2O)6 crystals changing to the matte greenish yellow of the trihydrate solid. The phase transformation and crystal structures were confirmed by density functional theory calculations and optical microscopy methods, both of which showed a transformation with two distinct sites for the uranyl cation in the trihydrate, with only one in the hexahydrate.

We investigate the (0001) surface of single crystal quartz with a submonolayer of Rb adsorbates. Using Rydberg atom electromagnetically induced transparency, we investigate the electric elds resulting from Rb adsorbed on the quartz surface, and measure the activation energy of the Rb adsorbates. We show that the Rb induces a negative electron affnity (NEA) on the quartz surface. The NEA surface allows for low energy electrons to bind to the surface and cancel the electric eld from the Rb adsorbates. Our results have implications for integrating Rydberg atoms into hybrid quantum systems and the fundamental study of atom-surface interactions, as well as applications for electrons bound to a 2D surface.

Water is a common surrogate of DNA for modelling the charged particle-induced ionizing processes in living tissue exposed to radiations. The present study aims at scrutinizing the validity of this approximation and then revealing new insights into proton-induced energy transfers by a comparative analysis between water and realistic biological medium. In this context, a self-consistent quantum mechanical modelling of the ionization and electron capture processes is reported within the continuum distorted wave-eikonal initial state framework for both isolated water molecules and DNA components impacted by proton beams. Their respective probability of occurrence - expressed in terms of total cross sections - as well as their energetic signature (potential and kinetic) are assessed in order to clearly emphasize the differences existing between realistic building blocks of living matter and the controverted water-medium surrogate. Consequences in radiobiology and radiotherapy will be discussed in particular in view of treatment planning refinement aiming at better radiotherapy strategies.

In this report, we present a thermodynamic-­based model of hydride precipitation in Zr-based claddings. The model considers the state of the cladding immediately following drying, after removal from cooling-pools, and presents the evolution of precipitate formation upon cooling as follows: The pilgering process used to form Zr-based cladding imparts strong crystallographic and grain shape texture, with the basal plane of the hexagonal α-Zr grains being strongly aligned in the rolling-­direction and the grains are elongated with grain size being approximately twice as long parallel to the rolling direction, which is also the long axis of the tubular cladding, as it is in the orthogonal directions.

The mechanical properties and stability of studtite, (UO2)(O2)(H2O)2·2H2O, and metastudtite, (UO2)(O2)(H2O)2, two important corrosion phases observed on spent nuclear fuel exposed to water, have been investigated using density functional perturbation theory. While (UO2)(O2)(H2O)2 satisfies the necessary and sufficient Born criteria for mechanical stability, (UO2)(O2)(H2O)2·2H2O is found to be mechanically metastable, which might be the underlying cause of the irreversibility of the studtite to metastudtite transformation. According to Pugh's and Poisson's ratios and the Cauchy pressure, both phases are considered ductile and shear modulus is the parameter limiting their mechanical stability. Debye temperatures of 294 and 271 K are predicted for polycrystalline (UO2)(O2)(H2O)2·2H2O and (UO2)(O2)(H2O)2, suggesting a lower micro-hardness of metastudtite.

Deep geological disposal of nuclear waste in clay/shale/argillaceous rock formations has received much consideration given its desirable attributes such as isolation properties (low permeability), geochemically reduced conditions, slow diffusion, sorbtive mineralogy, and geologically widespread (Jové Colón et al., 2014). There is a wealth of gained scientific expertise on the behavior of clay/shale/ argillaceous rock given its focus in international nuclear waste repository programs that includes underground research laboratories (URLs) in Switzerland, France, Belgium, and Japan. Jové Colón et al. (2014) have described some of these investigative efforts in clay rock ranging from site characterization to research on the engineered barrier system (EBS). Evaluations of disposal options that include nuclear waste disposition in clay/shale/argillaceous rock have determined that this host media can accommodate a wide range of waste types. R&D work within the Used Fuel Disposition Campaign (UFDC) assessing thermal effects and fluid-mineral interactions for the disposition of heat-generating waste have so far demonstrated the feasibility for the EBS and clay host rock to withstand high thermal loads. This report represents the continuation of disposal R&D efforts on the advancement and refinement of coupled Thermal-Hydrological-Mechanical-Chemical (THMC), hydrothermal experiments on clay interactions, used fuel degradation (source term), and thermodynamic modeling and database development. The development and implementation of a clay/shale/argillite reference case described in Jové Colón et al. (2014) for FY15 will be documented in another report (Mariner et al. 2015) – only a brief description will be given here. This clay reference case implementation is the result of integration efforts between the GDSA PA and disposal in argillite work packages. The assessment of sacrificial zones in the EBS is being addressed through experimental work along with 1D reactive-transport and reaction path modeling. The focus of these investigations into the nature of sacrificial zones is to evaluate the chemical effects of heterogeneous chemical reactions at EBS interfaces. The difference in barrier material types and the extent of chemical reactions within these interfacial domains generates changes in mineral abundances. These mineralogical alterations also result in volume changes that, although small, could affect the interface bulk porosity. As in previous deliverables, this report is structured according to various national laboratory contributions describing R&D activities applicable to clay/shale/argillite media.

The structure, lattice dynamics and thermodynamic properties of bulk technetium were investigated within the framework of density functional theory. The phonon density of states spectrum computed with density functional perturbation theory closely matches inelastic coherent neutron scattering measurements. The thermal properties of technetium were derived from phonon frequencies calculated within the quasi-harmonic approximation (QHA), which introduces a volume dependence of phonon frequencies as a part of the anharmonic effect. The predicted thermal expansion and isobaric heat capacity of technetium are in excellent agreement with available experimental data for temperatures up to ∼1600 K.

The structural, mechanical and thermodynamic properties of 1: 1 layered dioctahedral kaolinite clay, with ideal Al2Si2O5(OH)4 stoichiometry, were investigated using density functional theory corrected for dispersion interactions (DFT-D2). The bulk moduli of 56.2 and 56.0 GPa predicted at 298 K using the Vinet and Birch-Murnaghan equations of state, respectively, are in good agreement with the recent experimental value of 59.7 GPa reported for well-crystallized samples. The isobaric heat capacity computed for uniaxial deformation of kaolinite along the stacking direction reproduces calorimetric data within 0.7-3.0% from room temperature up to its thermal stability limit.

Abstract not provided.

Abstract not provided.

Abstract not provided.

This present work describes a quantum-mechanically based model of the electron- and proton-induced ionization of isolated pyrimidine molecules. The impact energies range from the target ionization threshold up to ~1 keV for electrons and from 10 keV up to 10 MeV for protons. The cross-section calculations are performed within the 1st Born approximation in which the ejected electron is described by a Coulomb wave whereas the incident and the scattered projectiles are both described by plane waves. The pyrimidine target is described using the Gaussian 09 software package. Furthermore, our theoretical predictions obtained are in good agreement with experimental absolute total cross sections, while large discrepancies are observed between existing semi-empirical models and the present calculations.

The reaction of Tc(+7) with H2O2 has been studied in H2SO4 and the speciation of technetium performed by UV–visible and 99-Tc NMR spectroscopy. UV–visible measurements show that for H2SO4 ≥ 9 M and H2O2 = 0.17 M, TcO3(OH)(H2O)2 reacts immediately and blue solutions are obtained, while no reaction occurs for H2SO4 < 9 M. The spectra of the blue solutions exhibit bands centered around 520 and 650 nm which are attributed to Tc(+7) peroxo species. Studies in 6 M H2SO4 show that TcO4− begins to react for H2O2 = 2.12 M and red solutions are obtained. The UV–visible spectra of the red species are identical to the one obtained from the reaction of TcO4− with H2O2 in HNO3 and consistent with the presence of TcO(O2)2(H2O)(OH). The 99-Tc NMR spectrum of the red solution exhibits a broad signal centered at +5.5 ppm vs TcO4− and is consistent with the presence of a low symmetry Tc(+7) molecule.

Several numerical codes for proton and electron transport in water - a commonly used surrogate of the living matter - have been reported in the literature. In the current work, we report on a home-made step-by-step Monte Carlo code, called TILDA-V, based on a complete set of multiple-differential and total cross sections for describing all the inelastic processes occurring throughout the slowing-down of protons in water and DNA.

The interaction of Cs adatoms with mono- or bi-layered graphene (MLG and BLG), either free-standing or on a SiO2 substrate, was investigated using density functional theory. The most stable adsorption sites for Cs are found to be hollow sites on both graphene sheets and graphene-veiled SiO2(0001). Larger dipole moments are created when a MLG-veiled SiO2(0001) substrate is used for adsorption of Cs atoms compared to the adsorption on free-standing MLG, due to charge transfer occurring between the MLG and the SiO2 substrate. For the adsorption of Cs on BLG-veiled SiO2(0001) substrate, these differences are smoothed out and the binding energies corresponding to different sites are nearly degenerate; smaller dipole moments created by the Cs adatoms on BLG compared to MLG are also predicted.

Uranyl nitrate is a key species in the nuclear fuel cycle, but is known to exist in different states of hydration, including the hexahydrate [UO2(NO3)2(H2O)6] (UNH) and the trihydrate [UO2(NO3)2(H2O)3] (UNT) forms. Their stabilities depend on both relative humidity and temperature. Both phases have previously been studied by infrared transmission spectroscopy, but the data were limited by both instrumental resolution and the ability to prepare the samples as pellets without desiccating it. We report time-resolved infrared (IR) measurements using an integrating sphere that allow us to observe the transformation from the hexahydrate to the trihydrate simply by flowing dry nitrogen gas over the sample. Hexahydrate samples were prepared and confirmed via known XRD patterns, then measured in reflectance mode. The hexahydrate has a distinct uranyl asymmetric stretch band at 949.0 cm-1 that shifts to shorter wavelengths and broadens as the sample dehydrates and recrystallizes to the trihydrate, first as a blue edge shoulder but ultimately resulting in a doublet band with reflectance peaks at 966 and 957 cm-1. The data are consistent with transformation from UNH to UNT since UNT has two non-equivalent UO22+ sites. The dehydration of UO2(NO3)2(H2O)6 to UO2(NO3)2(H2O)3 is both a morphological and structural change that has the lustrous lime green crystals changing to the dull greenish yellow of the trihydrate. Crystal structures and phase transformation were confirmed theoretically using DFT calculations and experimentally via microscopy methods. Both methods showed a transformation with two distinct sites for the uranyl cation in the trihydrate, as opposed to a single crystallographic site in the hexahydrate.

The elastic properties and mechanical stability of zirconium alloys and zirconium hydrides have been investigated within the framework of density functional perturbation theory. Results show that the lowest-energy cubic Pn3m polymorph of δ-ZrH1.5 does not satisfy all the Born requirements for mechanical stability, unlike its nearly degenerate tetragonal P42/mcm polymorph. Elastic moduli predicted with the Voigt-Reuss-Hill approximations suggest that mechanical stability of α-Zr, Zr-alloy and Zr-hydride polycrystalline aggregates is limited by the shear modulus. According to both Pugh's and Poisson's ratios, α-Zr, Zr-alloy and Zr-hydride polycrystalline aggregates can be considered ductile. The Debye temperatures predicted for γ-ZrH, δ-ZrH1.5 and ε-ZrH2 are D = 299.7, 415.6 and 356.9 K, respectively, while D = 273.6, 284.2, 264.1 and 257.1 K for the α-Zr, Zry-4, ZIRLO and M5 matrices, i.e. suggesting that Zry-4 possesses the highest micro-hardness among Zr matrices.

Abstract not provided.

Abstract not provided.

Electron-induced ionization of the tetrahydrofuran molecule, the commonly used surrogate of the DNA sugar-phosphate backbone, is theoretically described in this study within the 1st Born approximation. Comparisons between theory and recent experiments are reported in terms of doubly and singly differential cross sections.

This report documents the development, demonstration and validation of a mesoscale, microstructural evolution model for simulation of zirconium hydride {delta}-ZrH{sub 1.5} precipitation in the cladding of used nuclear fuels that may occur during long-term dry storage. While the Zr-based claddings are manufactured free of any hydrogen, they absorb hydrogen during service, in the reactor by a process commonly termed ‘hydrogen pick-up’. The precipitation and growth of zirconium hydrides during dry storage is one of the most likely fuel rod integrity failure mechanisms either by embrittlement or delayed hydride cracking of the cladding. While the phenomenon is well documented and identified as a potential key failure mechanism during long-term dry storage (NUREG/CR-7116), the ability to actually predict the formation of hydrides is poor. The model being documented in this work is a computational capability for the prediction of hydride formation in different claddings of used nuclear fuels. This work supports the Used Fuel Disposition Research and Development Campaign in assessing the structural engineering performance of the cladding during and after long-term dry storage. This document demonstrates a basic hydride precipitation model that is built on a recently developed hybrid Potts-phase field model that combines elements of Potts-Monte Carlo and the phase-field models. The model capabilities are demonstrated along with the incorporation of the starting microstructure, thermodynamics of the Zr-H system and the hydride formation mechanism.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Report RWEV-REP-001, Analysis of Postclosure Groundwater Impacts for a Geologic Repository for the Disposal of Spent Nuclear Fuel and High Level Radioactive Waste at Yucca Mountain, Nye County, Nevada was issued by the DOE in 2009 and is currently being updated. Sandia National Laboratories (SNL) provided support for the original document, performing calculations and extracting data from the Yucca Mountain Performance Assessment Model that were used as inputs to the contaminant transport and dose calculations by Jason Associates Corporation, the primary developers of the DOE report. The inputs from SNL were documented in LSA-AR-037, Inputs to Jason Associates Corporation in Support of the Postclosure Repository Supplemental Environmental Impact Statement. To support the updating of the original Groundwater Impacts document, SNL has reviewed the inputs provided in LSA-AR-037 to verify that they are current and appropriate for use. The results of that assessment are documented here.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.