Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “zirconium”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 775 records · Page 43
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Aluminum Boron Nitride Ferroelectric Field-Effect Transistors With ZnO Semiconductor Channel

The discovery of ferroelectricity in hafnium zirconium oxide (HfxZr 1−x O2) and related fluorite materials has spurred interest in ferroelectric devices suitable for integration with silicon integrated circuits (ICs), especially those that can be embedded in the back-end-of-line (BEOL) process. More recently, ferroelectricity has been found in wurtzite aluminum nitride-based materials, such as scandium and boron-doped aluminum nitride (Al 1−x ScxN and Al 1−x BxN). Although these materials currently have undesirably large coercive electric fields, and small breakdown electric field-to-coercive electric field ratio, their low processing temperatures and large remanent polarization offer intriguing possibilities for device applications. Furthermore, we report ferroelectric field-effect transistors (FeFETs) with a 15 nm thick Al0.88B0.12N layer and an 11 nm ZnO semiconductor channel, achieving a memory window >1 V and switching voltages ( V switch ) <±10 V.

42 ENGINEERING↗
DOE OSTI2022

Material Plasma Exposure eXperiment (MPEX) High Heat Flux Bellows for Component Alignment in High Microwave Environment

The Material Plasma Exposure eXperiment (MPEX) device is a linear plasma device developed to perform plasma material interaction experiments under the conditions prototypic of a fusion reactor divertor. MPEX has multiple systems that must be precisely aligned to the plasma axis, including an electron cyclotron heating system that emits up to 400 kW of microwave power into the vacuum vessel. Five distinct systems require precise alignment on the MPEX device, thus requiring four bellows, all of which are adjacent to the plasma at a relatively high heat flux of approximately 47 kW/m2 and microwave power regions. The MPEX high heat flux bellows (HHFB) is designed to deliver 6 degrees of freedom positioning. The HHFB includes titanium–zirconium–molybdenum (TZM) inserts that are brazed into a Glidcop AL-15 body using a high-temperature braze alloy, thus blocking direct line of sight to an edge-welded bellows from the plasma and microwave screen to block microwaves from the bellows. Further, a custom ConFlat knife edge is machined into the Glidcop AL-15, so the vacuum flanges do not need a braze or weld joint on the vacuum interface. Fingerstock or copper mesh is used to restrict microwave power from entering the interstitial space between the water-cooled Glidcop AL-15 body and edge-welded bellows. Glidcop AL-15 was selected as the material of choice for the water-cooled body because it can maintain mechanical integrity at elevated temperatures, and it also allows for a high-temperature braze. TZM was selected for its machinability and compatibility with vacuum and plasma requirements. A test article is also being considered for a similarly shaped component, the MPEX limiter, that will demonstrate the integrity of the braze joint under high thermal load. Results from this testing will be extrapolated to deduce the lifetime and integrity of the HHFB design.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

MPEX High Heat Flux Plasma Dump Design, Manufacture, and Articles Test

Here, the Material Plasma Exposure eXperiment (MPEX) device is a steady-state linear plasma device currently in the final design phase at the Oak Ridge National Laboratory. This device will reach ion fluences up to 10 31 m -2 and will be used to expose neutron-irradiated materials to divertor-relevant plasmas and to study the effects of plasma-material interactions. These studies will elucidate the complex effects of plasmas with divertor candidate materials capable of withstanding high heat flux and high fluences for next-generation fusion devices. Bidirectional plasma will be generated using a high-power (200 kW) helicon source. Plasma will be confined by superconducting magnets. The last plasma-facing component on the upstream side of the MPEX device is the dump, which has been designed to intercept plasma and energetic particles. The dump will have a total heat load of 9.2 kW. A copper alloy (Glidcop AL-15) was selected for use in the water-cooled flange design because of its high thermal conductivity, its retaining strength at elevated temperatures, and its ability to be used in the high-temperature braze joints used in this application. Titanium-zirconium-molybdenum (TZM) tiles are brazed to the Glidcop AL-15 flange using high-temperature braze alloy. External water-cooling channels are used on the dump flange to prevent water leakage inside the vacuum space. This article discusses the details of the high heat flux dump design, including the computational fluid dynamics (CFD) and structural analyses performed to validate the design to meet the operational requirements of the MPEX device.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

A Study on the Impact of Temperature-Dependent Ferroelectric Switching Behavior in 3D Memory Architecture

The flourishing development of neural networks that require exponentially growing amounts of data has presented an elevated demand for memory footprint. To address this, researchers have been exploring hardware accelerators with innovative memory architectures like 3D memory. These 3D memory architectures offer enhanced storage capacity and processing capabilities, at a cost of rising on-chip temperature during operation. Hafnium Zirconium Oxide (HZO) based Ferroelectric Random Access Memory (FeRAM) is a promising nonvolatile memory candidate in neural network hardware accelerators for its outstanding write performance and reliability. However, its implementation in the architecture regarding the temperature-dependent ferroelectric switching behavior has not been well studied. In this work, we study the thermal impacts on polarization switching through experimental devices and simulation results. We conduct the circuit and architecture-level simulations to showcase that one can exploit this temperature rise to reduce FeRAM's write voltage and write energy due to its unique temperature-activated polarization switching mechanisms. As the on-chip temperature increases to 351K (ambient temperature at 300K) due to neural network workloads, the access energy per bit can be reduced by 27.6% when a dynamic write voltage is applied.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Final-stage densification kinetics of direct current–sinteredZrB2

Final-stage sintering was analyzed for nominally phase pure zirconium diboride synthesized by borothermal reduction of high-purity ZrO2. Analysis was conducted on ZrB2 ceramics with relative densities greater than 90% using the Nabarro–Herring stress–directed vacancy diffusion model. Temperatures of 1900°C or above and an applied uniaxial pressure of 50 MPa were required to fully densify ZrB2 ceramics by direct current sintering. Ram travel data were collected and used to determine the relative density of the specimens during sintering. Specimens sintered between 1900 and 2100°C achieved relative densities greater than 97%, whereas specimens sintered below 1900°C failed to reach the final stage of sintering. The average grain size ranged from 1.0 to 14.7 μm. The activation energy was calculated from the slope of an Arrhenius plot that used the Kalish equation. The activation energy was 162 ± 34 kJ/mol, which is consistent with the activation energy for dislocation movement in ZrB2. The diffusion coefficients for dislocation motion that controls densification were 5.1 × 10−6 cm2/s at 1900°C and 5.1 × 10−5 cm2/s at 2100°C, as calculated from activation energy and average grain sizes. This study provides evidence that the dominant mechanism for final-stage sintering of ZrB2 ceramics is dislocation motion.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Unraveling Li growth kinetics in solid electrolytes due to electron beam charging

Revealing the local structure of solid electrolytes (SEs) with electron microscopy is critical for the fundamental understanding of the performance of solid-state batteries (SSBs). However, the intrinsic structural information in the SSB can be misleading if the sample’s interactions with the electron beams are not fully understood. In this work, we systematically investigate the effect of electron beams on Al-doped lithium lanthanum zirconium oxide (LLZO) under different imaging conditions. Li metal is observed to grow directly on the clean surface of LLZO. The Li metal growth kinetics and the morphology obtained are found to be heavily influenced by the temperature, accelerating voltage, and electron beam intensity. We prove that the lithium growth is due to the LLZO delithiation activated by a positive charging effect under electron beam emission. Our results deepen the understanding of the electron beam impact on SEs and provide guidance for battery material characterization using electron microscopy.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Emergent ferroelectricity in subnanometer binary oxide films on silicon

The critical size limit of voltage-switchable electric dipoles has extensive implications for energy-efficient electronics, underlying the importance of ferroelectric order stabilized at reduced dimensionality. Here, we report on the thickness-dependent antiferroelectric-to-ferroelectric phase transition in zirconium dioxide (ZrO 2 ) thin films on silicon. The emergent ferroelectricity and hysteretic polarization switching in ultrathin ZrO 2 , conventionally a paraelectric material, notably persists down to a film thickness of 5 angstroms, the fluorite-structure unit-cell size. This approach to exploit three-dimensional centrosymmetric materials deposited down to the two-dimensional thickness limit, particularly within this model fluorite-structure system that possesses unconventional ferroelectric size effects, offers substantial promise for electronics, demonstrated by proof-of-principle atomic-scale nonvolatile ferroelectric memory on silicon. Additionally, it is also indicative of hidden electronic phenomena that are achievable across a wide class of simple binary materials.

42 ENGINEERING↗
DOE OSTI2024

Supramolecular assembly of blue and green halide perovskites with near-unity photoluminescence

The metal-halide ionic octahedron is the optoelectronic unit for halide perovskites, and a crown ether–assisted supramolecular assembly approach can pack various ionic octahedra into tunable symmetries. In this work, we demonstrate near-unity photoluminescence quantum yield (PLQY) blue and green emission with the supramolecular assembly of hafnium (Hf) and zirconium (Zr) halide octahedral clusters. (18C6@K) 2 HfBr 6 powders showed blue emission with a near-unity PLQY (96.2%), and green emission was also achieved with (18C6@K) 2 ZrCl 4 Br 2 powders at a PLQY of 82.7%. In conclusion, these highly emissive powders feature facile low-temperature solution-based synthesis conditions and maintain high PLQY in solution-processable semiconductor inks under ambient conditions, and they were used in thin-film displays and emissive three-dimensional–printed architectures that exhibited high spatial resolution.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Shape evolution in even-mass 98-104 Zr isotopes via lifetime measurements using the $\gamma \gamma $-coincidence technique

The Zirconium (Z = 40) isotopic chain has attracted interest for more than four decades. The abrupt lowering of the energy of the first 2 + state and the increase in the transition strength B(E2; 2$^+_1$ →0$^+_1$) going from 98 Zr to 100 Zr has been the first example of “quantum phase transition” in nuclear shapes, which has few equivalents in the nuclear chart. Although a multitude of experiments have been performed to measure nuclear properties related to nuclear shapes and collectivity in the region, none of the measured lifetimes were obtained using the Recoil Distance Doppler Shift method in the γγ-coincidence mode where a gate on the direct feeding transition of the state of interest allows a strict control of systematical errors. Here this work reports the results of lifetime measurements for the first yrast excited states in 98-104 Zr carried out to extract reduced transition probabilities. The new lifetime values in γγ-coincidence and γ-single mode are compared with the results of former experiments. Recent predictions of the Interacting Boson Model with Configuration Mixing, the Symmetry Conserving Configuration Mixing model based on the Hartree–Fock–Bogoliubov approach and the Monte Carlo Shell Model are presented and compared with the experimental data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Lowering the Impedance of Lanthanum Strontium Manganite-Based Electrodes with Lanthanum Oxychloride and Lanthanum Scavenging Chloride Salts

The impact of infiltrating chloride salts on the electrochemical behavior of lanthanum strontium manganite-yttria stabilized zirconia (LSM-YSZ) cathodes was investigated under solid oxide fuel cell operation. Infiltrating a lanthanum chloride solution resulted in the formation of a lanthanum oxychloride (LaOCl) phase. A LaOCl phase also formed by infiltrating an ammonium chloride solution; however, lanthanum was scavenged from the LSM phase to form LaOCl. The third infiltrating solution, a combination of zirconium chloride and yttrium nitrate, formed LaOCl by scavenging lanthanum from LSM and produced YSZ nanoparticles. Electrochemical impedance spectroscopy results suggest that LaOCl improves oxygen adsorption kinetics compared to a baseline LSM-YSZ cathode, reducing the low frequency impedance by 30%. In addition, scavenging lanthanum from LSM improved oxygen ion diffusion polarization as indicated by the observed 40% reduction in high frequency impedance and improved serial ohmic resistance by 19%. Finally, YSZ nanoparticles further reduced the high frequency impedance and ohmic resistance by 45% and 23%, respectively. The findings reveal new strategies for lowering the impedance of LSM-YSZ cathodes.

Taylor, Thomas H.↗
DOE OSTI2025

Behavior of the Cadmium Pool in the INL Mk-IV Electrorefiner

The Mk-IV electrorefiner has been in service since 1996 recovering uranium from spent fuels from the Experimental Breeder Rector II and Fast Flux Test Facility sodium-cooled fast reactors. The electrorefiner vessel includes a cadmium pool beneath the salt. The voltage signal between the cadmium pool and reference electrode (Ag/AgCl type) provides information about the chemical condition of the cadmium pool with respect to uranium and zirconium saturation. This information is used to guide process control and sampling decisions. This paper describes how the voltage signal is related to chemistry by the analyses of data from Mk-IV operations and confirmatory experimentation.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Viscoplastic selfconsistent (VPSC) code (version 8, open source)

VPSC is a mean-field polycrystal plasticity code for the prediction of the mechanical response and microstructure evolution of polycrystalline aggregates. The VPSC code is the computational realization of the visco-plastic self-consistent theory of polycrystal deformation. Both the theory and details of its numerical implementation were originally reported by R.A. Lebensohn and C.N. Tomé: "A self-consistent approach for the simulation of plastic deformation and texture development of polycrystals: application to Zirconium alloys", Acta Metallurgica et Materialia 41, 2611 (1993). Since then, the VPSC code has experienced several improvements and extensions and it is nowadays extensively used to simulate plastic deformation of polycrystalline aggregates and to interpret experimental evidence on metals, minerals and polymers.

Lebensohn, Ricardo↗
DOE OSTI2024

NbZr_BCC_SolidSolution_128atoms_VASP6

We performed density functional theory (DFT) calculations for body-centered-cubic (BCC) structures with 128 lattices sites of solid solution binary alloys niobium-zirconium (Nb-Zr). The electronic structures of alloys have been calculated using Vienna Ab initio Simulation Package (VASP). Within this package the DFT approach is used to reduce many-body Schrodinger equation to set of single particle Kohn-Sham (KS) equations. The generalized electronic exchange-correlation functional is described by generalized gradient approximation with the Perdew-Burke-Ernzerhof parametrization. The electron-ion interactions is described by pseudopotentials developed within the plane-wave basis projector augmented-wave (PAW) approach \cite{PAW}. These pseudopotentials are available at the VASP portal (http://cms.mpi.univie.ac.at/vasp/). Our calculations have been run with the pseudopotentials treating s and p semi-core states as valence in case for the elements Nb and Zr. The electronic densities and potentials are expanded over plane-waves with energy cutoff of 350 eV. 2x2x2 k-mesh and normal precision were used. The alloys were modeled by supercell containing 128 randomly distributed atoms. At initial step the atoms occupy perfect bcc lattice cites. This initial structure was optimized until energy changes less than 1e-6 eV, while forces acting on atoms don't exceed 1e-2 eV/angstrom. The electron-ion interaction is described by PAW pseudopotentials. The calculations have been collected by sampling chemical compositions across the entire compositional range. The chemical compositions have been sampled by progressively changing the number of atoms per constituent by 4. For each chemical composition of binaries and ternaries, the first-principle calculations have been run for 100 randomized arrangements of the constituents on the BCC lattice sites. We collected data for a total of 3,100 randomized atomic structures over 31 chemical compositions. The calculations have been collected on NERSC-Perlmutter and OLCF-Summit using the VASP 6.3.2. The VASP calculations for every atomic structure have been performed in 2 main steps: 1. Starting from an ideal body-centered-cubic (BCC) structure, geometry optimization with low precision has been executed to perform a preliminary optimization of the atomic structure. The output for this calculations is available in the files 0.CONTCAR, 0.OUTCAR, rlx1.out. 2. Using the atomic structure resulting from the preliminary geometry optimization, a second geometry optimization has been performed using normal precision. The output for this calculations is available in the files CONTCAR, OUTCAR, rlx2.out, vaspout.h5, and vasprun.xml. Cases 1-10 have been run without generating the file 'vaspout.h5'. Every chemical composition sampled across the composition range in the dataset has its own directory. The convention used to name the directories for binary alloys is AXBY, where A and B refer to the constituents, whereas X and Y are positive integers that represent the number of atoms for each constituent and their values still sum up to 128. Each atomic structure associated with a specific chemical composition has its own sub-directory within the directory of the corresponding chemical composition. The sub-directories for each atomic structure for each chemical composition are named 'case-*', where * is a positive integer that spans all the values from 1 through 100, extremes included. The files contained in each sub-directory 'case-*' for each atomic structure are as follows: FILES contained in each subdirectory with name "case-N" where N ranges between 11 and 100, extremes included: 1. INCAR: input file that contains various parameters and settings for controlling the behavior of the electronic structure calculations 2. KPOINTS: input file that specifies the Bloch vectors (k points) used to sample the Brillouin zone 3. 0.POSCAR: input file that defines the atomic structure of a system 4. 0.CONTCAR: output file that provides the atomic positions and cell parameters after the first geometry optimization has been run with the precision variable set to PREC=Low in the INCAR file 5. 0.OUTCAR: output file that contains detailed information about the progress of a calculation after the first geometry optimization has been run with the precision variable set to PREC=Low in the INCAR file 6. rlx1.out: file with diagnostic information about the execution of the first geometry optimization with precision variable set to PREC=Low in the INCAR file 7. POSCAR: input file that defines the atomic structure of a system after the first geometry optimization has been run at low precision. This represents the input for the second geometry optimization run with the precision variable set to PREC=Normal in the INCAR file 8. CONTCAR: output file that provides the atomic positions and cell parameters after the second geometry optimization has been run with the precision variable set to PREC=Normal in the INCAR file 9. OUTCAR: output file that contains detailed information about the progress of a calculation after the second geometry optimization has been run with the precision variable set to PREC=Normal in the INCAR file 10. rlx2.out: file with diagnostic information about the execution of the second geometry optimization with precision variable set to PREC=Normal in the INCAR file 11. vaspout.h5: hierarchical HDF5 file containing the inputs and outputs of a VASP calculation. To analyze the data in this file we recommend using py4vasp. This file is only produced if the VASP version used is compiled with HDF5 support 12. vasprun.xml: contains similar information to OUTCAR, but in an xml format. 13. CHGCAR: contains the charge density data of a system. This data is crucial for analyzing electronic structures, calculating electrostatic potential, and studying the distribution of charge in a crystal or molecular system FILES contained in each subdirectory with name "case-N" where N ranges between 1 and 10, extremes included: 1. INCAR: input file that contains various parameters and settings for controlling the behavior of the electronic structure calculations 2. KPOINTS: input file that specifies the Bloch vectors (k points) used to sample the Brillouin zone 3. {ID}.POSCAR: input file that defines the atomic structure of a system at the beginning of ID execution of geometry optimization with PREC=LOW 4. {ID}.CONTCAR: output file that provides the atomic positions and cell parameters at the end of ID execution of geometry optimization with PREC=LOW in the INCAR file 5. {ID}.OUTCAR: output file that contains detailed information about the progress of a calculation after the ID execution of geometry optimization that has been run with the precision variable set to PREC=Low in the INCAR file 6. rlx1.{ID}.out: file with diagnostic information about the execution of the ID execution of the geometry optimization with precision variable set to PREC=Low in the INCAR file 7. N{ID}.POSCAR: input file that defines the atomic structure of a system after the geometry optimization run at low precision. This represents the input for the ID execution of the second geometry optimization run with the precision variable set to PREC=Normal in the INCAR file 8. N{ID}.CONTCAR: output file that provides the atomic positions and cell parameters after the ID execution of the geometry optimization run with the precision variable set to PREC=Normal in the INCAR file 9. N{ID}.OUTCAR: output file that contains detailed information about the progress of a calculation after the ID execution of the geometry optimization run with the precision variable set to PREC=Normal in the INCAR file 10. rlx2.{ID}.out: file with diagnostic information about the ID execution of geometry optimization with precision variable set to PREC=Normal in the INCAR file 11. vaspout.h5: hierarchical HDF5 file containing the inputs and outputs of a VASP calculation. To analyze the data in this file we recommend using py4vasp. This file is only produced if the VASP version used is compiled with HDF5 support 12. vasprun.xml: contains similar information to OUTCAR, but in an xml format. 13. CHGCAR: contains the charge density data of a system. This data is crucial for analyzing electronic structures, calculating electrostatic potential, and studying the distribution of charge in a crystal or molecular system This research is sponsored by the Artificial Intelligence Initiative as part of the Laboratory Directed Research and Development (LDRD) Program of Oak Ridge National Laboratory, managed by UT-Battelle, LLC, for the US Department of Energy under contract DE-AC05-00OR22725. This work used resources of the Oak Ridge Leadership Computing Facility, which is supported by the Office of Science of the U.S. Department of Energy under Contract No. DE-AC05-00OR22725, under Directorate Discretionary awards MAT025 (Materials Science) and LRN026 (Machine Learning), and INCITE award MAT201. This work also used resources of the National Energy Research Scientific Computing Center, which is supported by the Office of Science of the U.S. Department of Energy under Contract No. DE-AC02-05CH11231, under award ERCAP0025216. REFERENCES (1) Kresse, G. & Hafner, J. Ab initio molecular dynamics for liquid metals. Phys. review B 47, 558 (1993). (2) Kresse, G. & Hafner, J. Ab initio molecular-dynamics simulation of the liquid-metal–amorphous-semiconductor transition in germanium. Phys. Rev. B 49, 14251 (1994) (3) Kresse, G. & Furthmüller, J. Efficiency of ab-initio total energy calculations for metals and semiconductors using a plane-wave basis set. Comput. materials science 6, 15–50 (1996) (4) Kresse, G. & Furthmüller, J. Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set. Phys. review B 54, 11169 (1996) (5) Kresse, G. & Joubert, D. From ultrasoft pseudopotentials to the projector augmented-wave method. Phys. review b 59, 1758 (1999)

36 MATERIALS SCIENCE↗
DOE OSTI2024

TaZr_BCC_SolidSolution_128atoms_VASP6

We performed density functional theory (DFT) calculations for body-centered-cubic (BCC) structures with 128 lattices sites of solid solution binary alloys tantalum-zirconium (Ta-Zr). The electronic structures of alloys have been calculated using Vienna Ab initio Simulation Package (VASP). Within this package the DFT approach is used to reduce many-body Schrodinger equation to set of single particle Kohn-Sham (KS) equations. The generalized electronic exchange-correlation functional is described by generalized gradient approximation with the Perdew-Burke-Ernzerhof parametrization. The electron-ion interactions is described by pseudopotentials developed within the plane-wave basis projector augmented-wave (PAW) approach \cite{PAW}. These pseudopotentials are available at the VASP portal (http://cms.mpi.univie.ac.at/vasp/). Our calculations have been run with the pseudopotentials treating s and p semi-core states as valence in case for the elements Ta and Zr. The electronic densities and potentials are expanded over plane-waves with energy cutoff of 350 eV. 2x2x2 k-mesh and normal precision were used. The alloys were modeled by supercell containing 128 randomly distributed atoms. At initial step the atoms occupy perfect bcc lattice cites. This initial structure was optimized until energy changes less than 1e-6 eV, while forces acting on atoms don't exceed 1e-2 eV/angstrom. The electron-ion interaction is described by PAW pseudopotentials. The calculations have been collected by sampling chemical compositions across the entire compositional range. The chemical compositions have been sampled by progressively changing the number of atoms per constituent by 4. For each chemical composition of binaries and ternaries, the first-principle calculations have been run for 100 randomized arrangements of the constituents on the BCC lattice sites. We collected data for a total of 3,100 randomized atomic structures over 31 chemical compositions. The calculations have been collected on NERSC-Perlmutter and OLCF-Summit using the VASP 6.3.2. The VASP calculations for every atomic structure have been performed in 2 main steps: 1. Starting from an ideal body-centered-cubic (BCC) structure, geometry optimization with low precision has been executed to perform a preliminary optimization of the atomic structure. The output for this calculations is available in the files 0.CONTCAR, 0.OUTCAR, rlx1.out. 2. Using the atomic structure resulting from the preliminary geometry optimization, a second geometry optimization has been performed using normal precision. The output for this calculations is available in the files CONTCAR, OUTCAR, rlx2.out, vaspout.h5, and vasprun.xml. Cases 1-10 have been run without generating the file 'vaspout.h5'. Every chemical composition sampled across the composition range in the dataset has its own directory. The convention used to name the directories for binary alloys is AXBY, where A and B refer to the constituents, whereas X and Y are positive integers that represent the number of atoms for each constituent and their values still sum up to 128. Each atomic structure associated with a specific chemical composition has its own sub-directory within the directory of the corresponding chemical composition. The sub-directories for each atomic structure for each chemical composition are named 'case-*', where * is a positive integer that spans all the values from 1 through 100, extremes included. The files contained in each sub-directory 'case-*' for each atomic structure are as follows: FILES contained in each subdirectory with name "case-N" where N ranges between 11 and 80, extremes included: 1. INCAR: input file that contains various parameters and settings for controlling the behavior of the electronic structure calculations 2. KPOINTS: input file that specifies the Bloch vectors (k points) used to sample the Brillouin zone 3. 0.POSCAR: input file that defines the atomic structure of a system 4. 0.CONTCAR: output file that provides the atomic positions and cell parameters after the first geometry optimization has been run with the precision variable set to PREC=Low in the INCAR file 5. 0.OUTCAR: output file that contains detailed information about the progress of a calculation after the first geometry optimization has been run with the precision variable set to PREC=Low in the INCAR file 6. rlx1.out: file with diagnostic information about the execution of the first geometry optimization with precision variable set to PREC=Low in the INCAR file 7. POSCAR: input file that defines the atomic structure of a system after the first geometry optimization has been run at low precision. This represents the input for the second geometry optimization run with the precision variable set to PREC=Normal in the INCAR file 8. CONTCAR: output file that provides the atomic positions and cell parameters after the second geometry optimization has been run with the precision variable set to PREC=Normal in the INCAR file 9. OUTCAR: output file that contains detailed information about the progress of a calculation after the second geometry optimization has been run with the precision variable set to PREC=Normal in the INCAR file 10. rlx2.out: file with diagnostic information about the execution of the second geometry optimization with precision variable set to PREC=Normal in the INCAR file 11. vaspout.h5: hierarchical HDF5 file containing the inputs and outputs of a VASP calculation. To analyze the data in this file we recommend using py4vasp. This file is only produced if the VASP version used is compiled with HDF5 support 12. vasprun.xml: contains similar information to OUTCAR, but in an xml format. 13. CHGCAR: contains the charge density data of a system. This data is crucial for analyzing electronic structures, calculating electrostatic potential, and studying the distribution of charge in a crystal or molecular system FILES contained in each subdirectory with name "case-N" where N ranges between 1 and 10 and between 81 and 100, extremes included: 1. INCAR: input file that contains various parameters and settings for controlling the behavior of the electronic structure calculations 2. KPOINTS: input file that specifies the Bloch vectors (k points) used to sample the Brillouin zone 3. {ID}.POSCAR: input file that defines the atomic structure of a system at the beginning of ID execution of geometry optimization with PREC=LOW 4. {ID}.CONTCAR: output file that provides the atomic positions and cell parameters at the end of ID execution of geometry optimization with PREC=LOW in the INCAR file 5. {ID}.OUTCAR: output file that contains detailed information about the progress of a calculation after the ID execution of geometry optimization that has been run with the precision variable set to PREC=Low in the INCAR file 6. rlx1.{ID}.out: file with diagnostic information about the execution of the ID execution of the geometry optimization with precision variable set to PREC=Low in the INCAR file 7. N{ID}.POSCAR: input file that defines the atomic structure of a system after the geometry optimization run at low precision. This represents the input for the ID execution of the second geometry optimization run with the precision variable set to PREC=Normal in the INCAR file 8. N{ID}.CONTCAR: output file that provides the atomic positions and cell parameters after the ID execution of the geometry optimization run with the precision variable set to PREC=Normal in the INCAR file 9. N{ID}.OUTCAR: output file that contains detailed information about the progress of a calculation after the ID execution of the geometry optimization run with the precision variable set to PREC=Normal in the INCAR file 10. rlx2.{ID}.out: file with diagnostic information about the ID execution of geometry optimization with precision variable set to PREC=Normal in the INCAR file 11. vaspout.h5: hierarchical HDF5 file containing the inputs and outputs of a VASP calculation. To analyze the data in this file we recommend using py4vasp. This file is only produced if the VASP version used is compiled with HDF5 support 12. vasprun.xml: contains similar information to OUTCAR, but in an xml format. 13. CHGCAR: contains the charge density data of a system. This data is crucial for analyzing electronic structures, calculating electrostatic potential, and studying the distribution of charge in a crystal or molecular system This research is sponsored by the Artificial Intelligence Initiative as part of the Laboratory Directed Research and Development (LDRD) Program of Oak Ridge National Laboratory, managed by UT-Battelle, LLC, for the US Department of Energy under contract DE-AC05-00OR22725. This work used resources of the Oak Ridge Leadership Computing Facility, which is supported by the Office of Science of the U.S. Department of Energy under Contract No. DE-AC05-00OR22725, under Directorate Discretionary awards MAT025 (Materials Science) and LRN026 (Machine Learning), and INCITE award MAT201. This work also used resources of the National Energy Research Scientific Computing Center, which is supported by the Office of Science of the U.S. Department of Energy under Contract No. DE-AC02-05CH11231, under award ERCAP0025216. REFERENCES (1) Kresse, G. & Hafner, J. Ab initio molecular dynamics for liquid metals. Phys. review B 47, 558 (1993). (2) Kresse, G. & Hafner, J. Ab initio molecular-dynamics simulation of the liquid-metal–amorphous-semiconductor transition in germanium. Phys. Rev. B 49, 14251 (1994) (3) Kresse, G. & Furthmüller, J. Efficiency of ab-initio total energy calculations for metals and semiconductors using a plane-wave basis set. Comput. materials science 6, 15–50 (1996) (4) Kresse, G. & Furthmüller, J. Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set. Phys. review B 54, 11169 (1996) (5) Kresse, G. & Joubert, D. From ultrasoft pseudopotentials to the projector augmented-wave method. Phys. review b 59, 1758 (1999)

36 MATERIALS SCIENCE↗
DOE OSTI2024

VZr_BCC_SolidSolution_128atoms_VASP6

We performed density functional theory (DFT) calculations for body-centered-cubic (BCC) structures with 128 lattices sites of solid solution binary alloys vanadium-zirconium (V-Zr). The electronic structures of alloys have been calculated using Vienna Ab initio Simulation Package (VASP). Within this package the DFT approach is used to reduce many-body Schrodinger equation to set of single particle Kohn-Sham (KS) equations. The generalized electronic exchange-correlation functional is described by generalized gradient approximation with the Perdew-Burke-Ernzerhof parametrization. The electron-ion interactions is described by pseudopotentials developed within the plane-wave basis projector augmented-wave (PAW) approach \cite{PAW}. These pseudopotentials are available at the VASP portal (http://cms.mpi.univie.ac.at/vasp/). Our calculations have been run with the pseudopotentials treating s and p semi-core states as valence in case for the elements V and Zr. The electronic densities and potentials are expanded over plane-waves with energy cutoff of 350 eV. 2x2x2 k-mesh and normal precision were used. The alloys were modeled by supercell containing 128 randomly distributed atoms. At initial step the atoms occupy perfect bcc lattice cites. This initial structure was optimized until energy changes less than 1e-6 eV, while forces acting on atoms don't exceed 1e-2 eV/angstrom. The electron-ion interaction is described by PAW pseudopotentials. The calculations have been collected by sampling chemical compositions across the entire compositional range. The chemical compositions have been sampled by progressively changing the number of atoms per constituent by 4. For each chemical composition of binaries and ternaries, the first-principle calculations have been run for 100 randomized arrangements of the constituents on the BCC lattice sites. We collected data for a total of 3,100 randomized atomic structures over 31 chemical compositions.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Designing an impedance-based technique for studying corrosion on cladding materials

In nuclear cladding materials, changes in structure and chemistry over the lifetime of claddings can have a significant impact on fuel and reactor performance. Currently, chemical evolution over the lifetime of the fuel and cladding materials can only be inferred using post-irradiation examination. As an alternative to post-irradiation measurements, electrochemical impedance spectroscopy (EIS) can provide an in-pile measurement of changes in chemistry in oxides and hydrides. The objective of this study is to develop an impedance based sensing technology for measuring spatial- and time-resolved in-pile changes in cladding chemistry in order to understand the mechanisms of hydride formation, hydride dealloying and associated cladding corrosion so that a mitigation strategy and technology can be developed for reactor safe operation. Two prototype impedance-based sensors were designed and assembled using Zr-4 as working electrode and Pt as counter and reference electrode. One sensor is tubular type and the other one is plate-type. These sensors will be tested in static autoclave to investigate the effects of various environmental conditions (including temperature, pressure, and water chemistry) on cladding material corrosion mechanisms. These in-situ EIS sensors would provide corrosion-layer parameters in high-temperature water, such as (1) total layer thickness, (2) the evolution of oxide structures as uniform layer dielectric or sublayer(s), and (3) the possible electrochemistry of the volume corrosion process and the solid/liquid interface, such as the cyclic nature of corrosion of Zr alloys in high-temperature water. This work was focused on electrochemical sensor development, targeting cladding oxides/hydrides under thermodynamic equilibria; characterization of samples using Kelvin probe force microscopy (SKPFM), Raman spectroscopy, and scanning electron microscopy (SEM) with energy-dispersive X-ray spectroscopy (EDS) to provide high-resolution co-localized characterization of zirconium oxide and metal/oxide interface.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

A review of cladding failure thresholds in RIA conditions based on transient reactor test data and the need for continued testing

Transient reactor experiments on light water reactor (LWR) fuel pins had been conducted since the beginning of the nuclear era to help determine core coolability and cladding failure thresholds. During one such test in November of 1993 at the CABRI transient test reactor on a test involving a high burnup fuel rod with a corroded Zircaloy-4 cladding it was first observed that cladding failures could occur prior to a departure from nucleate boiling (pre-DNB) at lower-than-expected peak radial average enthalpies. This paper will present an independent review of the publicly available transient reactor test database on higher burnup LWR pins conducted at the CABRI and NSRR reactors as well as review of a selection of published out of pile mechanical testing methods. The purpose of the review is to determine how well the new regulatory limits are supported by experimental data. The review will identify if additional transient reactor tests could provide additional support for the NRC guidance or identify the need for revisions. The evaluation will consider how far the existing database can be extrapolated when considering low hydrogen zirconium alloy claddings (with and without protective coatings) containing very high burnup (> 70 MWd/kgU) UO2 fuel pellets. Finally, the authors will suggest how out of pile mechanical tests can be used in conjunction with a limited number of transient reactor tests to develop cladding specific failure thresholds in RIA type transients.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Effect of surface oxides on tritium entrance and permeation in FeCrAl alloys for nuclear fuel cladding: a review

Abstract Iron-chromium-aluminum (FeCrAl) alloys are being considered for the cladding of uranium dioxide fuel in light water reactors (LWRs). FeCrAl alloys have good mechanical properties at temperatures of 300 °C and higher, and have superlative resistance to attack by steam at temperatures of up to 1000 °C and higher. A concern has been raised that the use of FeCrAl for cladding would result in a higher content of tritium in the reactor coolant as compared with the current system where the cladding is a zirconium based alloy. This review shows that the flux of tritium from the fuel rod cavities to the coolant across the fuel cladding wall will be greatly reduced by the presence of oxides on the surface of the cladding. The review of current literature and permeation data show that (a) protective oxides are expected to be present on both sides of the FeCrAl cladding, and (b) depending on the characteristics of these oxide layers it is reasonable to expect about two–three orders of magnitude reduction in tritium permeation, relative to the permeation response in clean, unoxidized condition for FeCrAl steels of interest, around 277 °C–377 °C temperatures.

Electrochemistry↗