Publications

Abstract not provided.

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.

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.

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.

Abstract not provided.

Abstract not provided.

Abstract not provided.

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

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.

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 structure–property relationships of bulk CeO₂ and Ce₂O₃ 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 CeO₂ and Ce₂O₃ with superior accuracy, especially when +U is chosen close to its value derived by the linear-response approach. Lastly, 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 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.

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.

The supercritical carbon dioxide (S-CO2) Brayton Cycle has gained significant attention in the last decade as an advanced power cycle capable 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 ready 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₂, and (5) Identifying the influence of gas chemistry on alloy corrosion. Progress in each of these areas is provided in this report.

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

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.

A C++ library (called Eshelby) was implemented in fiscal year 2015 based upon the formulas documented in this report. The library implements a generalized version of Eshelby's inclusion problem. The library was written as a set of functions which can be called from another program; the principle intended use cases are kinetic models of precipitate formation in zirconium claddings where use of the Eshelby library provides needed elastic energy density calculations, as well as calculations of stress and strain in and around precipitates; it is intended that the library will be made open source. For isotropic inclusions in the form of oblate and prolate ellipsoids, the Eshelby library can be used for nearly any relevant/appropriate shape parameters to calculate strains, stresses and energy density at interior and exterior points. The Eshelby library uses a combination of analytical formulas and numerical routines making it very extensible. For example, the library can can easily be extended to include inclusions such as spheres since analytical expressions exist for the required elliptic integrals; similarly, general ellipsoids do not have analytical results for the required elliptic integrals but those integrals can be numerically evaluated and thus fit naturally into the Eshelby library. This report documents all formulas implemented in the Eshelby library and presents some demonstration calculations relevant to the intended application.

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.