Publications

Gordon, Margaret E.; Bryan, Charles R.; Greathouse, Jeffery A.; Weck, Philippe F.; Kim, Eunja; Payne, Clay P.

Abstract not provided.

Greathouse, Jeffery A.; Weck, Philippe F.; Gordon, Margaret E.; Kim, Eunja; Bryan, Charles R.

Abstract not provided.

Hanz, Elizabeth K.; Gordon, Margaret E.; Weck, Philippe F.; Greathouse, Jeffery A.; P., Meserole S.; Cannady, Manuela R.; Payne, Clay P.; Bryan, Charles R.

Abstract not provided.

Physical Chemistry Chemical Physics

Weck, Philippe F.; Jove Colon, Carlos F.; Kim, Eunja

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

International High-Level Radioactive Waste Management 2019, IHLRWM 2019

Greathouse, Jeffery A.; Weck, Philippe F.; Gordon, Margaret E.; Kim, Eunja; Bryan, Charles R.

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.

International High-Level Radioactive Waste Management 2019, IHLRWM 2019

Bryan, Charles R.; Gordon, Margaret E.; Weck, Philippe F.; Greathouse, Jeffery A.; Kim, Eunja; Payne, Clay P.

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.

International High-Level Radioactive Waste Management 2019, IHLRWM 2019

Kim, Eunja; Gordon, M.E.; Weck, Philippe F.; Greathouse, Jeffery A.; Meserole, S.P.; Rodriguez, Mark A.; Payne, Clay P.; Bryan, Charles R.

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.

Surface Science

Jooya, H.Z.; McKay, K.S.; Kim, E.; Weck, Philippe F.; Pappas, D.P.; Hite, D.A.; Sadeghpour, H.R.

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.

Sassani, David C.; Brady, Patrick V.; Gelbard, Fred G.; Price, Laura L.; Prouty, Jeralyn L.; Rechard, Robert P.; Rigali, Mark J.; Rogers, Ralph D.; Sanchez, Amanda C.; Walkow, Walter M.; Weck, Philippe F.

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.

Gordon, Margaret E.; Bryan, Charles R.; Weck, Philippe F.; Greathouse, Jeffery A.; Payne, Clay P.; Kim, Eunja K.

Abstract not provided.

Dingreville, Remi P.; Hattar, Khalid M.; Boyle, Timothy J.; Monterrosa, Anthony M.; Barr, Christopher M.; Weck, Philippe F.; Juan, Pierre-Alexandre J.

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

Gordon, Margaret E.; Bryan, Charles R.; Weck, Philippe F.; Greathouse, Jeffery A.; Payne, Clay P.

Abstract not provided.

Journal of Raman Spectroscopy

Weck, Philippe F.; Gordon, Margaret E.; Greathouse, Jeffery A.; Bryan, Charles R.; Meserole, Stephen M.; Rodriguez, Mark A.; Payne, Clay P.; Kim, Eunja

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.

Physical Review B

Weck, Philippe F.; Cochrane, Kyle C.; Root, Seth R.; Lane, J.M.; Shulenburger, Luke N.; Carpenter, John H.; Sjostrom, Travis; Mattsson, Thomas M.; Vogler, Tracy V.

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.

Acta Crystallographica Section C: Structural Chemistry

Borjas Nevarez, Rosendo; Balasekaran, Samundeeswari M.; Kim, Eunja; Weck, Philippe F.; Poineau, Frederic

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.

Weck, Philippe F.; Cochrane, Kyle C.; Walker, Seth R.; Lane, James M.; Shulenburger, Luke S.; Carpenter, John H.; Mattsson, Thomas M.; Vogler, Tracy V.

Abstract not provided.

Nature

Evlyukhin, E E.; Kim, E K.; Goldberger, D G.; Cifligu, P C.; Weck, Philippe F.; Pravica, M P.

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.

AISTech - Iron and Steel Technology Conference Proceedings

Kim, Eunja; Weck, Philippe F.; Borjas, Rosendo; Balasekaran, Eswari; Poineau, Frederic

Abstract not provided.

Chemical Physics Letters

Kim, Eunja; Weck, Philippe F.; Borjas, Rosendo; Poineau, Frederic

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.

Scientific Reports

Weck, Philippe F.; Kim, Eunja; Wang, Yifeng; Kruichak, Jessica N.; Mills, Melissa M.; Matteo, Edward N.; Pellenq, Roland J.M.

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.

Sassani, David C.; Price, Laura L.; Rechard, Robert P.; Rogers, Ralph D.; Walkow, Walter M.; Johnson, Ava J.; Sanchez, Amanda C.; Mariner, Paul M.; Rigali, Mark J.; Stein, Emily S.; Weck, Philippe F.

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.

Scientific Reports

Biswas, S B.; Champion, C C.; Weck, Philippe F.; Tribedi, L C.

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.

Journal of Physical Chemistry. C

Weck, Philippe F.; Juan, Pierre-Alexandre J.; Dingreville, Remi P.; Kim, Eunja K.

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

European Physical Journal D

Quinto, Michele A.; Monti, Juan M.; Weck, Philippe F.; Fojón, Omar A.; Hanssen, Jocelyn; Rivarola, Roberto D.; Weck, Philippe F.; Champion, Christophe

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

Shulenburger, Luke N.; Cochrane, Kyle C.; Mattsson, Thomas M.; Lane, James M.; Weck, Philippe F.; Vogler, Tracy V.; Desjarlais, Michael P.

Abstract not provided.