Publications

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

Abstract not provided.

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

Abstract not provided.

Physical Review A

Kim, E.; Safavi-Naini, A.; Hite, D.A.; McKay, K.S.; Pappas, D.P.; Weck, Philippe F.; Sadeghpour, H.R.

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.

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

Abstract not provided.

Kim, E.K.; Weck, Philippe F.; Borjas, R.B.; Poineau, F P.

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.

Kim, E.K.; Weck, Philippe F.; Poineau, F.P.; Paviet, P.P.

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.

Sneed, D.S.; Pravica, M.P.; Kim, E.K.; Weck, Philippe F.

Abstract not provided.

Sassani, David C.; Jang, Jay J.; Mariner, Paul M.; Price, Laura L.; Rechard, Robert P.; Rigali, Mark J.; Rogers, Ralph R.; Stein, Emily S.; Walkow, Walter M.; Weck, Philippe F.

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.

Journal of Physical Chemistry C

Weck, Philippe F.; Kim, Eunja

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.

Walker, Matthew W.; Walker, Matthew W.; Kruizenga, Alan M.; Kruizenga, Alan M.; Weck, Philippe F.; Weck, Philippe F.; Withey, Elizabeth A.; Withey, Elizabeth A.; Fleming, Darryn F.; Fleming, Darryn F.; Rochau, Gary E.; Rochau, Gary E.

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

Lavina, B L.; Kim, E.K.; Cynn, H.C.; Weck, Philippe F.; Siska, E S.

Abstract not provided.

Sandia journal manuscript; Not yet accepted for publication

Lavina, Barbara L.; Kim, Eunja K.; Cynn, Hyunchae C.; Weck, Philippe F.; Seaborg, Kelly S.; Siska, Emily S.; Meng, Yue M.; Evans, Williams E.

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.

Chemical Physics Letters

Weck, Philippe F.; Kim, Eunja; Wang, Yifeng

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.

Journal of Physics and Chemistry of Solids

Mast, Daniel S.; Lavina, Barbara L.; Kim, Eunja K.; Siska, Emily S.; Weck, Philippe F.; Poineau, F.P.; Czerwinski, K R.; Forster, P M.

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.

Physical Chemistry Chemical Physics

Weck, Philippe F.; Kim, Eunja

The structure-property relationships of bulk CeO2 and Ce2O3 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 CeO2 and Ce2O3 with superior accuracy, especially when +U is chosen close to its value derived by the linear-response approach. 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.

RSC Advances

Weck, Philippe F.; Kim, Eunja

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.

Radiation Physics and Chemistry

Fromm, Michel; Quinto, Michele A.; Weck, Philippe F.; Champion, Christophe

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.

Journal of Physical Chemistry A

Johnson, Timothy J.; Sweet, Lucas E.; Meier, David E.; Mausolf, Edward J.; Kim, Eunja; Weck, Philippe F.; Buck, Edgar C.; McNamara, Bruce K.

Uranyl nitrate is a key species in the nuclear fuel cycle. However, this species is known to exist in different states of hydration, including the hexahydrate ([UO2(NO3)2(H2O)6] often called UNH), the trihydrate [UO2(NO3)2(H2O)3 or UNT], and in very dry environments the dihydrate form [UO2(NO3)2(H2O)2]. Their relative stabilities depend on both water vapor pressure and temperature. In the 1950s and 1960s, the different phases were studied by infrared transmission spectroscopy but were limited both by instrumental resolution and by the ability to prepare the samples for transmission. We have revisited this problem using time-resolved reflectance spectroscopy, which requires no sample preparation and allows dynamic analysis while the sample is exposed to a flow of N2 gas. Samples of known hydration state were prepared and confirmed via X-ray diffraction patterns of known species. In reflectance mode the hexahydrate UO2(NO3)2(H2O)6 has a distinct uranyl asymmetric stretch band at 949.0 cm^-1 that shifts to shorter wavelengths and broadens as the sample desiccates and recrystallizes to the trihydrate, first as a shoulder growing in on the blue edge but ultimately results in a doublet band with reflectance peaks at 966 and 957 cm^-1. The data are consistent with transformation from UNH to UNT as UNT has two inequivalent UO2²⁺ sites. The dehydration of UO2(NO3)2(H2O)6 to UO2(NO3)2(H2O)3 is both a structural and morphological change that has the lustrous lime green UO2(NO3)2(H2O)6 crystals changing to the matte greenish yellow of the trihydrate solid. The phase transformation and crystal structures were confirmed by density functional theory calculations and optical microscopy methods, both of which showed a transformation with two distinct sites for the uranyl cation in the trihydrate, with only one in the hexahydrate.

Sandia journal manuscript; Not yet accepted for publication

Sedlacek, J.A.; Kim, E.K.; Rittenhouse, S.T.; Weck, Philippe F.; Sadeghpour, H.R.; Shaffer, J.P.

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.

Physics in Medicine and Biology

Champion, C.; Quinto, M.A.; Monti, J.M.; Galassi, M.E.; Weck, Philippe F.; Fojón, O.A.; Hanssen, J.; Rivarola, R.D.

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.

Tikare, Veena T.; Weck, Philippe F.; Mitchell, John A.

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.

RSC Advances

Weck, Philippe F.; Kim, Eunja; Buck, Edgar C.

The mechanical properties and stability of studtite, (UO2)(O2)(H2O)2·2H2O, and metastudtite, (UO2)(O2)(H2O)2, two important corrosion phases observed on spent nuclear fuel exposed to water, have been investigated using density functional perturbation theory. While (UO2)(O2)(H2O)2 satisfies the necessary and sufficient Born criteria for mechanical stability, (UO2)(O2)(H2O)2·2H2O is found to be mechanically metastable, which might be the underlying cause of the irreversibility of the studtite to metastudtite transformation. According to Pugh's and Poisson's ratios and the Cauchy pressure, both phases are considered ductile and shear modulus is the parameter limiting their mechanical stability. Debye temperatures of 294 and 271 K are predicted for polycrystalline (UO2)(O2)(H2O)2·2H2O and (UO2)(O2)(H2O)2, suggesting a lower micro-hardness of metastudtite.

Jove Colon, Carlos F.; Weck, Philippe F.; Hammond, Glenn E.; Kuhlman, Kristopher L.; Zheng, L.Z.; Rutqvist, J.R.; Kim, K.W.; Houseworth, J.H.; Caporuscio, F.A.C.; Cheshire, M.C.; Palaich, S.P.; Norskog, K.E.; Zavarin, M.Z.; Wolery, T.J.W.; Jerden, J.J.; Copple, J.M.C.; Cruse, T.C.; Ebert, W.E.; Jove Colon, Carlos F.

Deep geological disposal of nuclear waste in clay/shale/argillaceous rock formations has received much consideration given its desirable attributes such as isolation properties (low permeability), geochemically reduced conditions, slow diffusion, sorbtive mineralogy, and geologically widespread (Jové Colón et al., 2014). There is a wealth of gained scientific expertise on the behavior of clay/shale/ argillaceous rock given its focus in international nuclear waste repository programs that includes underground research laboratories (URLs) in Switzerland, France, Belgium, and Japan. Jové Colón et al. (2014) have described some of these investigative efforts in clay rock ranging from site characterization to research on the engineered barrier system (EBS). Evaluations of disposal options that include nuclear waste disposition in clay/shale/argillaceous rock have determined that this host media can accommodate a wide range of waste types. R&D work within the Used Fuel Disposition Campaign (UFDC) assessing thermal effects and fluid-mineral interactions for the disposition of heat-generating waste have so far demonstrated the feasibility for the EBS and clay host rock to withstand high thermal loads. This report represents the continuation of disposal R&D efforts on the advancement and refinement of coupled Thermal-Hydrological-Mechanical-Chemical (THMC), hydrothermal experiments on clay interactions, used fuel degradation (source term), and thermodynamic modeling and database development. The development and implementation of a clay/shale/argillite reference case described in Jové Colón et al. (2014) for FY15 will be documented in another report (Mariner et al. 2015) – only a brief description will be given here. This clay reference case implementation is the result of integration efforts between the GDSA PA and disposal in argillite work packages. The assessment of sacrificial zones in the EBS is being addressed through experimental work along with 1D reactive-transport and reaction path modeling. The focus of these investigations into the nature of sacrificial zones is to evaluate the chemical effects of heterogeneous chemical reactions at EBS interfaces. The difference in barrier material types and the extent of chemical reactions within these interfacial domains generates changes in mineral abundances. These mineralogical alterations also result in volume changes that, although small, could affect the interface bulk porosity. As in previous deliverables, this report is structured according to various national laboratory contributions describing R&D activities applicable to clay/shale/argillite media.

Dalton Transactions

Weck, Philippe F.; Kim, Eunja

The structure, lattice dynamics and thermodynamic properties of bulk technetium were investigated within the framework of density functional theory. The phonon density of states spectrum computed with density functional perturbation theory closely matches inelastic coherent neutron scattering measurements. The thermal properties of technetium were derived from phonon frequencies calculated within the quasi-harmonic approximation (QHA), which introduces a volume dependence of phonon frequencies as a part of the anharmonic effect. The predicted thermal expansion and isobaric heat capacity of technetium are in excellent agreement with available experimental data for temperatures up to ∼1600 K.

Dalton Transactions

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

The structural, mechanical and thermodynamic properties of 1: 1 layered dioctahedral kaolinite clay, with ideal Al2Si2O5(OH)4 stoichiometry, were investigated using density functional theory corrected for dispersion interactions (DFT-D2). The bulk moduli of 56.2 and 56.0 GPa predicted at 298 K using the Vinet and Birch-Murnaghan equations of state, respectively, are in good agreement with the recent experimental value of 59.7 GPa reported for well-crystallized samples. The isobaric heat capacity computed for uniaxial deformation of kaolinite along the stacking direction reproduces calorimetric data within 0.7-3.0% from room temperature up to its thermal stability limit.