Publications

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.

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.

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.

Abstract not provided.

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.

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.