Publications

The negative thermal expansion (NTE) material Zr2(WO4)(PO4)2 has been investigated for the first time within the framework of the density functional perturbation theory (DFPT). The structural, mechanical, and thermodynamic properties of this material have been predicted using the Perdew, Burke and Ernzerhof for solid (PBEsol) exchange-correlation functional, which showed superior accuracy over standard functionals in previous computational studies of the NTE material α-ZrW2O8. The bulk modulus calculated for Zr2(WO4)(PO4)2 using the Vinet equation of state at room temperature is K0 = 63.6 GPa, which is in close agreement with the experimental estimate of 61.3(8) at T = 296 K. The computed mean linear coefficient of thermal expansion is -3.1 × 10-6 K-1 in the temperature range ∼0-70 K, in line with the X-ray diffraction measurements. The mean Grüneisen parameter controlling the thermal expansion of Zr2(WO4)(PO4)2 is negative below 205 K, with a minimum of -2.1 at 10 K. The calculated standard molar heat capacity and entropy are CP0 = 287.6 and S0 = 321.9 J·mol-1·K-1, respectively. The results reported in this study demonstrate the accuracy of DFPT/PBEsol for assessing or predicting the relationship between structural and thermomechanical properties of NTE materials.

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.

Abstract not provided.

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

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.

Elastic and thermodynamic properties of negative thermal expansion (NTE) αα-ZrW2O8 have been calculated using PBEsol and PBE exchange-correlation functionals within the framework of density functional perturbation theory (DFPT). Measured elastic constants are reproduced within ~2% with PBEsol and 6% with PBE. The thermal evolution of the Grüneisen parameter computed within the quasi-harmonic approximation exhibits negative values below the Debye temperature, consistent with observation. The standard molar heat capacity is predicted to be C $O\atop{P}$=192.2 and 193.8 J mol^-1K^-1 with PBEsol and PBE, respectively. These results suggest superior accuracy of DFPT/PBEsol for studying the lattice dynamics, elasticity and thermodynamics of NTE materials.

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.

Using combined experimental and computational approaches, we show that at 43 GPa and 1300 K gallium phosphide adopts the super-Cmcm structure, here indicated with its Pearson notation oS24. First-principles enthalpy calculations demonstrate that this structure is more thermodynamically stable above ∼20 GPa than previously proposed polymorphs. In contrast to other polymorphs, the oS24 phase shows a strong bonding differentiation and distorted fivefold coordination geometries of both P atoms. The shortest bond of the phase is a single covalent P-P bond measuring 2.171(11) Å at synthesis pressure. Phosphorus dimerization in GaP sheds light on the nature of the super-Cmcm phase and provides critical new insights into the high-pressure polymorphism of octet semiconductors. Bond directionality and anisotropy explain the relatively low symmetry of this high-pressure phase.

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.

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.

X-ray induced damage has been known for decades and has largely been viewed as a tremendous nuisance. We, on the other hand, harness the highly ionizing and penetrating properties of hard X-rays to initiate novel decomposition and synthetic chemistry. Here, we show that powdered cesium oxalate monohydrate pressurized to ≤0.5 GPa and irradiated with X-rays of energies near the cesium K-edge undergoes molecular and structural transformations with one of the final products exhibiting a new type of bcc crystal structure that has previously not been observed. Additionally, based on cascades of ultrafast electronic relaxation steps triggered by the absorption of one X-ray photon, we propose a model explaining the X-ray induced damage of multitype bounded matter. As X-rays are ubiquitous, these results show promise in the preparation of novel compounds and novel structures that are inaccessible via conventional methods. They may offer insight into the formation of complex organic compounds in outer space.

Abstract not provided.

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.

Cubic zirconium tungstate (α-ZrW₂O₈), a notorious negative thermal expansion (NTE) material, has been investigated within the framework of density functional perturbation theory (DFPT), combined with experimental characterization to assess and validate computational results. Spectroscopic, mechanical and thermodynamic properties have been derived from DFPT calculations. 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 supe-rior accuracy of the PBEsol exchange-correlation functional over standard PBE calculations. The thermal evolution of the Grüneisen parameter computed within the quasi-harmonic approximation exhibits negative values below the Debye temperature, consistent with the observed NTE characteristics of α-ZrW₂O₈. The standard molar heat capacity is predicted to be C$0\atop{P}$=193.8 and 192.2 J.mol^-1.K^-1 with PBE and PBEsol, respectively, ca. 7% lower than calorimetric data. In conclusion, these results demonstrate the accuracy of the DFPT/PBEsol approach for studying the spectroscopic, mechanical and thermodynamic properties of materials with anomalous thermal expansion.

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 Ce1-xZrxO2 (x = 0-1) solid solutions, selected Ce1-xZrxO2 surfaces, and Ce1-xZrxO2/CeO2 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)O2 phases, indicating a significant rise in microhardness from CeO2 to ZrO2, without appreciable loss in ductility as the interfacial stoichiometry changes. Surface energy calculations for the low-index CeO2(111) and (110) surfaces show limited sensitivity to strong 4f-electron correlation. The fracture energy of Ce1-xZrxO2(111)/CeO2(111) increases markedly with Zr content, with a significant decrease in energy for thicker Ce1-xZrxO2 films. These findings suggest the crucial role of Zr acting as a binder at the Ce1-xZrxO2/CeO2 interfaces, due to the more covalent character of Zr-O bonds compared to Ce-O. The impact of surface relaxation upon interface cracking was assessed and found to reach a maximum for Ce0.25Zr0.75O2/CeO2 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.