Publications

Ab initio molecular dynamics (AIMD) simulations were carried out to investigate the equation of state of Nb2O5 and its pressure-density relationship under shock conditions. The focus of this study is on the monoclinic B−Nb2O5 (C2/c) polymorph. Enthalpy calculations from AIMD trajectories at 300 K show that the pressure-induced transformation between the thermodynamically most stable crystalline monoclinic parent phase H−Nb2O5 (P2/m) and B−Nb2O5 occurs at ∼1.9 GPa. This H→B transition is energetically more favorable than the H→L(Pmm2) pressure-induced transition recently observed at ∼5.9−9.0 GPa. The predicted shock properties of Nb2O5 polymorphs are also compared to their Nb and NbO2 counterparts to assess the impact of niobium oxidation on shock response.

Abstract not provided.

Abstract not provided.

Glycoboehmite (GB) materials are synthesized by a solvothermal reaction to form layered aluminum oxyhydroxide (boehmite) modified by intercalated butanediol molecules. These hybrid materials offer a platform to design materials with potentially novel sorption, wetting, and catalytic properties. Several synthetic methods have been used, resulting in different structural and spectroscopic properties, but atomistic detail is needed to determine the interlayer structure to explore the synthetic control of GB materials. Here, we use classical molecular dynamics (MD) simulations to compare the structural properties of GB interlayers containing chemisorbed butanediol molecules as a function of diol loading. Accompanying quantum (density functional theory, DFT) static calculations and MD simulations are used to validate the classical model and compute the infrared spectra of various models. Classical MD results reveal the existence of two unique interlayer environments at higher butanediol loading, corresponding to smaller (cross-linked) and expanded interlayers. DFT-computed infrared spectra reveal the sensitivity of the aluminol O-H stretch frequencies to the interlayer environment, consistent with the spectrum of the synthesized material. Insight from these simulations will aid in the characterization of the newly synthesized GB materials.

Chromium self-diffusion through stainless steel (SS) matrix and along grain boundaries is an important mechanism controlling SS structural materials corrosion. Cr diffusion in austenitic SS was simulated using canonical ab initio molecular dynamics with realistic models of type-316 SS bulk, with and without Cr vacancies, and a low-energy Σ3 twin boundary typically observed at active corrosion sites. Cr self-diffusion coefficients at 750 and 850 °C calculated using Einstein's diffusion equation are 4.2 × 10−6 and 8.1 × 10−6 Å2 ps−1 in pristine bulk, 3.8 × 10−3 and 5.5 × 10−3 Å2 ps−1 in bulk including Cr vacancies, and 9.5 × 10−2 and 1.0 × 10−1 Å2 ps−1 at a Σ3[1 1 1]60° twin boundary.

Intermolecular Coulombic decay (ICD) in liquid water is a relatively novel type of nonlocal electronic decay mechanism, competing with the traditional mechanism of proton transfer between neighboring water molecules. Key features of ICD are its ultrafast non-radiative decay process and ultralong-range for excess energy transfer from the excited atom/molecule to its neighbors. Since detecting unambiguous ICD signatures in bulk liquid water is technically challenging, small water clusters have often been utilized to gain insights into ICD and other ionization processes in aqueous environment. Here, we present results from quantum mechanical calculations of the electronic structures of neutral to multiply-ionized water monomer, dimer, trimer, and tetramer. Core-level electrons of water are also considered here since recent studies demonstrated that emission site and energy of the electrons released during resonant-Auger-ICD cascade can be controlled by coupling ICD to resonant core excitation. Previous studies of ICD and electronic structures of neutral and ionized small water clusters and liquid water are briefly discussed.

This report represents completion of milestone deliverable M2SF-24SN010309082 Annual Status Update for OWL due on November 30, 2023. It contains the status of fiscal year 2023 (FY2023) updates for the Online Waste Library (OWL).

This report describes research and development (R&D) activities conducted during Fiscal Year 2023 (FY23) in the Advanced Fuels and Advanced Reactor Waste Streams Strategies work package in the Spent Fuel Waste Science and Technology (SFWST) Campaign supported by the United States (U.S.) Department of Energy (DOE). This report is focused on evaluating and cataloguing Advanced Reactor Spent Nuclear Fuel (AR SNF) and Advanced Reactor Waste Streams (ARWS) and creating Back-end Nuclear Fuel Cycle (BENFC) strategies for their disposition. The R&D team for this report is comprised of researchers from Sandia National Laboratories and Enviro Nuclear Services, LLC.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Polymorphism and phase transitions in sodium diuranate, Na2U2O7, are investigated with density functional perturbation theory (DFPT). Thermal properties of crystalline α-, β- and γ-Na2U2O7 polymorphs are predicted from DFPT phonon calculations, i.e., the first time for the high-temperature γ-Na2U2O7 phase (R3̄m symmetry). The standard molar isochoric heat capacities predicted within the quasi-harmonic approximation are for P21/a α-Na2U2O7 and C2/m β-Na2U2O7, respectively. Gibbs free energy calculations reveal that α-Na2U2O7 (P21/a) and β-Na2U2O7 (C2/m) are almost energetically degenerate at low temperature, with β-Na2U2O7 becoming slightly more stable than α-Na2U2O7 as temperature increases. These findings are consistent with XRD data showing a mixture of α and β phases after cooling of γ-Na2U2O7 to room temperature and the observation of a sluggish α → β phase transition above ca. 600 K. A recently observed α-Na2U2O7 structure with P21 symmetry is also shown to be metastable at low temperature. Based on Gibbs free energy, no direct β → γ solid-solid phase transition is predicted at high temperature, although some experiments reported the existence of such phase transition around 1348 K. This, along with recent experiments, suggests the occurrence of a multi-step process consisting of initial β-phase decomposition, followed by recrystallization into γ-phase as temperature increases.

Abstract not provided.

Abstract not provided.

The overarching goal of the combined computational and experimental R&D activities proposed in this project is to enhance understanding of the mechanisms and thermal-mechanical-chemical (TMC) parameters controlling the instant release fraction (IRF) and matrix dissolution of high-burnup (HB; burnup) spent nuclear fuels (SNFs) and the subsequent formation, stability, and phase transformations of SNF alteration products under long-term storage and geological disposal conditions. Uranium dioxide may undergo oxidative corrosion/alteration, and the IRF may be increased for HB SNF, both of which may affect environmental systems associated with SNF long-term storage and disposal. The oxidative matrix dissolution may form various complex uranyl-based phases, including a rich variety of oxides, silicates, carbonates and other secondary minerals in varied geological environments (e.g., studtite, metastudtite, amorphous uranyl peroxide, uranium trioxide, triuranium octoxide, schoepite, dehydrated schoepite, metaschoepite, becquerelite, soddyite, rutherfordine,...). These uranyl phases generally have higher mobility UO₂⁺² species than less soluble U⁴⁺ phases. However, limited information on the thermodynamic properties and formation kinetics of these uranyl-bearing phases is available to predict explicitly paragenesis under the conditions relevant to long-term storage or disposal. The proposed project draws on complementary expertise and research backgrounds from the team members: (i) to apply a combined ab initio modeling (UNLV/UTEP and SNL) and experimental (UNLV) strategy investigating the high-temperature TMC mechanisms of alteration of SNF under α-radiolysis conditions; (ii) to investigate the mechanistic of phase transformations in UNF degradation products under various conditions expected in long-term storage systems (e.g. (UO₂)O₂(H₂O)₄ → (UO₂)O₂(H₂O)₂ → U₂O₇ → UO₃ → U₃O₈); (iii) to determine high-accuracy TMC parameters for complex uranyl-based phases formed in storage or geological disposal environments (e.g. UO₃(H₂O)₂, Ca[(UO₂)₆O₄(OH)₈]₈H₂O, (UO₂)₂(SiO₄)₃2H₂O,…). The unforeseen COVID-19 pandemic led to the laboratory/campus closure since March 2020, that resulted in a significant delay in reaching milestones in a satisfactory manner, due to (i) the statewide recommendation from stop-working to later limited work in the lab and work-from-home (WFH), (ii) no in-person interactions, and (iii) a hiring freeze at UNLV. Therefore, a no cost extension (10/01/2021- 9/30/2022) was requested to help make up the time we lost during the global pandemic in 2020-2021, leading to paradigm shifts in the focus of the project in the following three main tasks: Task 1 (Computational), Task 2 (Experimental), and Task 3 (Final report, due on 12/29/2022).

This report represents completion of milestone deliverable M2SF-23SN010309082 Annual Status Update for OWL due on November 30, 2022. It provides the status of fiscal year 2022 (FY2022) updates for the Online Waste Library (OWL).

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.

This report represents the milestone deliverable M4SF-22SN010309092 “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 model development of spent nuclear fuel (SNF). These activities are part of the 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-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.

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 + MgCl₂ 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 work 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 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.

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

The equation of state (EOS) and shock compression of bulk vanadium were investigated using canonical ab initio molecular dynamic simulations, with experimental validation to 865 GPa from shock data collected at Sandia's Z Pulsed Power Facility. In simulations the phase space was sampled along isotherms ranging from 3000 K to 50000 K, for densities between -ü=3 and 15g/cm3, with a focus on the liquid regime and the body-centered-cubic phase in the vicinity of the melting limit. The principal Hugoniot predicted from first principles is overall consistent with shock data, while it showed that current multiphase SESAME-type EOS for vanadium needed revision in the liquid regime. A more accurate SESAME EOS was developed using constraints from experiments and simulations. This work emphasizes the need to use a combined theoretical and experimental approach to develop high-fidelity EOS models for extreme conditions.

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.

The DOE R&D program under the Spent Fuel Waste Science Technology (SFWST) campaign has made key progress in modeling and experimental approaches towards the characterization of chemical and physical phenomena that could impact the long-term safety assessment of heat-generating nuclear waste disposition in deep clay/shale/argillaceous rock. International collaboration activities such as heater tests and postmortem analysis of samples recovered from these have elucidated key information regarding changes in the engineered barrier system (EBS) material exposed to years of thermal loads. Chemical and structural analyses of sampled bentonite material from such tests has as well as experiments conducted on these are key to the characterization of thermal effects affecting bentonite clay barrier performance and the extent of sacrificial zones in the EBS during the thermal period. Thermal, hydrologic, and chemical data collected from heater tests and laboratory experiments has been used in the development, validation, and calibration of THMC simulators to model near-field coupled processes. This information leads to the development of simulation approaches (e.g., continuum vs. discrete) to tackle issues related to flow and transport at various scales of the host-rock and EBS design concept. Consideration of direct disposal of large capacity dual-purpose canisters (DPCs) as part of the back-end SNF waste disposition strategy has generated interest in improving our understanding of the effects of elevated temperatures on the EBS design. This is particularly important for backfilled repository concepts where temperature plays a key role in the EBS behavior and long-term performance. This report describes multiple R&D efforts on disposal in argillaceous geologic media through development and application of coupled THMC process models, experimental studies on clay/metal/cement barrier and host-rock (argillite) material interactions, molecular dynamic (MD) simulations of water transport during (swelling) clay dehydration, first-principles studies of metaschoepite (UO₂ corrosion product) stability, and advances in thermodynamic plus surface complexation database development. Drift-scale URL experiments provides key data for testing hydrological-chemical (HC) model involving strong couplings of fluid mixing and barrier material chemical interactions. The THM modeling focuses on heater test experiments in argillite rock and gas migration in bentonite as part of international collaboration activities at underground research laboratories (URLs). In addition, field testing at an URL involves in situ analysis of fault slip behavior and fault permeability. Pore-scale modeling of gas bubble migration is also being investigated within the gas migration modeling effort. Interaction experiments on bentonite samples from heater test under ambient and elevated temperatures permit the evaluation of ion exchange, phase stability, and mineral transformation changes that could impact clay swelling. Advances in the development, testing, and implementation of a spent nuclear fuel (SNF) degradation model coupled with canister corrosion focus on the effects of hydrogen gas generation and its integration with Geologic Disposal Safety Assessment (GDSA). GDSA integration activities includes evaluation of groundwater chemistries in shale formations.

The high-pressure response of titanium dioxide (TiO2) is of interest because of its numerous industrial applications and its structural similarities to silica (SiO2). We used three platforms - Sandia's Z machine, Omega Laser Facility, and density-functional theory-based quantum molecular dynamics (QMD) simulations - to study the equation of state (EOS) of TiO2 at extreme conditions. We used magnetically accelerated flyer plates at Sandia to measure Hugoniot of TiO2 up to pressures of 855 GPa. We used a laser-driven shock wave at Omega to measure the shock temperature in TiO2. Our Z data show that rutile TiO2 reaches 2.2-fold compression at a pressure of 855 GPa and Omega data show that TiO2 is a reflecting liquid above 230 GPa. The QMD simulations are in excellent agreement with the experimental Hugoniot in both pressure and temperature. A melt curve for TiO2 is also proposed based on the QMD simulations. The combined experimental results show TiO2 is in a liquid at these explored pressure ranges and is not highly incompressible as suggested by a previous study.

The high-pressure response of titanium dioxide (TiO2) is of interest because of its numerous industrial applications and its structural similarities to silica (SiO2). We used three platforms - Sandia's Z machine, Omega Laser Facility, and density-functional theory-based quantum molecular dynamics (QMD) simulations - to study the equation of state (EOS) of TiO2 at extreme conditions. We used magnetically accelerated flyer plates at Sandia to measure Hugoniot of TiO2 up to pressures of 855 GPa. We used a laser-driven shock wave at Omega to measure the shock temperature in TiO2. Our Z data show that rutile TiO2 reaches 2.2-fold compression at a pressure of 855 GPa and Omega data show that TiO2 is a reflecting liquid above 230 GPa. The QMD simulations are in excellent agreement with the experimental Hugoniot in both pressure and temperature. A melt curve for TiO2 is also proposed based on the QMD simulations. The combined experimental results show TiO2 is in a liquid at these explored pressure ranges and is not highly incompressible as suggested by a previous study.

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.

The incorporation of uranium, plutonium and technetium in the negative thermal expansion (NTE) α-Zr(WO4)2 has been investigated within the framework of density functional theory (DFT). It is found that the vacancy formation energies of the charged vacancies are overall larger than that of its counterpart neutral Frenkel defects and Schottky defects. DFT calculations suggest that U and Pu substitutions for the Zr site are preferred in α-Zr(WO4)2. In case of Tc substitution, both Tc(IV) for the Zr site and Tc(VII) for the W site are considered under oxygen-poor and oxygen-rich conditions, while Tc(VII) substitution can be improved significantly by including Y2O3 (charge compensation).

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.

The relationship between the structure and thermodynamic properties of schoepite, an important uranyl phase with formula [(UO₂)₈O₂(OH)₁₂]·12H₂O formed upon corrosion of UO₂, 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 S⁰ = 179.60 J mol^-1 K^-1 and C⁰_P = 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 (UO₃·xH₂O, x = 2.25) to dehydrated schoepite (x = 1) is predicted to be essentially linear along isotherms ranging from 100 to 500 K. Finally, these findings have important implications for the dehydration of layered uranyl corrosion phases and hygroscopic materials.

This report describes the potential of a novel class of materials—α-ZrW₂O₈, Zr₂WP₂O₁₂, 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₂O₈ 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 Zr₂WP₂O₁₂ 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.

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.