Search NASA⌕ Search

SEARCH · Search NASA

Results for “ZIRCONIUM ALLOY”

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 109 records · Page 6

Microstructure, Thermal, and Mechanical Properties Relationships in U and UZr Alloys (Final Report)

Uranium-zirconium (U-Zr) alloys are candidate fuel systems for transmutation based reactors that can be used to burn long-lived minor actinides and fission products in fast spectrum reactors. Metallic fuels have also been gaining more recent attention for applications as accident tolerant fuels, as well as for use in small modular reactors. This research focused on a “science-based” approach to capture the connections between U and U-Zr alloys’ three-dimensional (3-D) microstructure, thermal properties, and mechanical properties through closely coordinated experiments and modeling efforts from the unirradiated to the irradiated fuels. Advanced characterization and modeling techniques were used to understand irradiation-induced microstructural evolution and its direct impact on the thermal and mechanical properties of U and U-Zr fuel. Closely coordinated experiments and modeling were performed to provide crucial data that does not currently exist. Overall, this research spanned multiple length and time scales within the models and experiments. The scope of the research encompassed the understanding of the irradiation effects in U and various U-Zr alloys with particular attention paid to three task areas: (1) microstructural evolution, (2) in-situ/ex-situ thermal and mechanical properties, and (3) multiscale modeling connections to microstructure, thermal, and mechanical properties. This research resulted in (1) the 3-D characterization of neutron irradiated U-Zr fuel in multiple phase regions to better understand fission gas swelling and constituent redistribution, (2) development of a microstructural model linking thermal and mechanical properties via in situ Raman and nanoindentation, (3) and mesoscale phase field modeling was coupled with the AEH method in the MOOSE framework was used to calculate the effective thermal conductivities of U-Zr fuels consisting of α-U and δ-UZr 2 heterogeneous microstructures.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Development of Innovative Non-Concentrated Alloy (NCA) Fuel Cladding for Advanced Nuclear Reactors

As the nuclear energy sector advances toward next-generation reactors, the need for high-performance fuel cladding materials has become increasingly urgent. Traditional alloys like zirconium and stainless steel are reaching their performance limits under higher temperatures, more corrosive coolants, and extended irradiation. This report presents the development of a new class of fuel cladding materials based on Non-Concentrated Alloys (NCAs) ? multi-element systems designed to deliver enhanced mechanical strength, corrosion resistance, and radiation tolerance. Through a combination of computational modeling (Computer Coupling of Phase Diagrams and Thermochemistry: CALPHAD), simulation-guided alloy selection, and experimental fabrication via arc melting and spark plasma sintering (SPS), three strategic alloy design paths were explored: (1) FeCrAl-based NCAs, (2) refractory-lean neutron-efficient alloys, and (3) equimolar high-entropy compositions. Microstructural analysis confirmed the formation of stable body-centered cubic BCC_A2 phases, while mechanical testing demonstrated hardness values significantly exceeding those of conventional cladding materials. The results highlight the tunability of NCA systems and their potential for balancing strength and ductility ? a critical consideration for in-reactor performance. Looking forward, future work will focus on thermomechanical optimization, CALPHAD refinement, and benchmarking against industry standards to enable scalable deployment. This work not only advances the science of nuclear materials but also supports broader goals in nuclear safety, performance, and nuclear energy innovation.

36 - MATERIALS SCIENCE↗

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↗

Impact of Coating Defects on Performance of Coated Zirconium Cladding

Research on accident tolerant fuels (ATF) has started after the Fukushima accident [1–3]. While efforts have been expended on both fuel and cladding ATF concepts, the bulk of work has been devoted to improved cladding. The overarching goal of these approaches is to extend the coping time available during a severe accident before any event would result in release of radioactivity to the public. The most basic ATF cladding concept is obtained by applying a thin coating of highly corrosion-resistant material on the surface of a licensed zirconium cladding alloy. This thin coating is intended to not interfere with the neutronic or mechanical performance of the base cladding under normal operating conditions. Different coating materials, thicknesses, coating processes, process parameters, and testing methods have an impact on the microstructure and mechanical properties and therefore on the results of the applied investigation methods. These challenges have motivated an initial focus on demonstrating that the presence of coatings do not perturb the critical performance benchmarks of uncoated material. Ongoing lead test assembly irradiations of coated zirconium concepts in commercial reactions is intended to establish baseline performance in this regard in the coming years.

36 MATERIALS SCIENCE↗

DISSOLUTION OF SURROGATE U-Zr FUEL USING ALNIFLEX CONDITIONS

Non-aluminum clad spent nuclear fuel (NASNF) stored in L Basin at the Savannah River Site (SRS) is widely varied in fuel composition, design, packaging, and physical condition. The complexity of the NASNF inventory presents significant challenges, and technology development is necessary for successful disposition. One such fuel in the inventory is metallic uranium-zirconium (U-Zr) alloy fuel, the focus of this study. Electrolytic or nitric acid only dissolution of metallic U-Zr alloy can form insoluble zirconium oxide, which results in up to 52% loss of U to insoluble solids, and can be subject to potentially uncontrolled oxidation reactions [1, 2]. The AlNiflex process was determined to be a viable dissolution flowsheet for the U-Zr alloy fuel. Under a narrow set of solution concentrations, a combination of hydrofluoric acid (HF), nitric acid (HNO3), aluminum nitrate (Al(NO3)3), and hexavalent chromium can safely dissolve U-Zr intermetallic alloys, keep Zr soluble, and not significantly corrode stainless steel (SS) vessels [3, 4.

Gogolski, Jarrod M. [Savannah River National Labor↗

Hydrothermal corrosion and steam oxidation behavior comparison of UAM and conventional Zry-4

Here, hydrothermal corrosion and steam exposures on ultrasonic additively manufactured (UAM) Zircaloy-4 (Zry-4) specimens were conducted to compare to wrought Zry-4 performance in both simulated pressurized water reactor (PWR) operation (330°C, 15 MPa, <50 ppb oxygen, 500h) and representative accident conditions (100% H 2 O(g) at 700, 900, and 1100°C for up to 1h). UAM Zry-4 specimens showed 50% higher mass gain in PWR conditions and 50-80% more mass gain in steam oxidation due to internal oxidation between build layers in both cases. UAM Zry-4 in 900°C steam showed breakaway oxidation leading to sample deformation. Titanium used as an interlayer in the UAM process was detected in the substrate after the 1100°C exposure. Additional experiments to better link as-built defects to oxidation behavior are needed.

36 MATERIALS SCIENCE↗

Structure evolution and tin redistribution during oxidation of Zircaloy-4 at 500°C

Zirconium (Zr) alloys are widely used as fuel cladding in nuclear power reactors due to their thermal stability, mechanical durability, corrosion resistance, and low neutron absorption cross-section. However, their performance is challenged by oxidation in reactor environments, making the study of Zr alloy corrosion behavior crucial for ensuring the safety, longevity, and economic viability of nuclear power systems. While the oxidation behavior of Zr-based cladding materials has been extensively studied since the 1950s, a mechanistic understanding into the relationship between structure evolution, solute element redistribution, and properties remains elusive. Valuable insights may be obtained through advanced experimental methods, such as in-situ and high resolution microscopy techniques. Here, in this study, the oxidation behavior of Zircaloy-4 at 500 °C in O 2 is characterized using a multimodal advanced characterization approach. Using in-situ X-ray diffraction, the phase evolution from metastable to stable oxides is tracked in real time. Complementary high-resolution techniques, including electron microscopy and atom probe tomography, reveal nanoscale insights into the microstructural changes and solute redistribution across the oxide/metal interface. Nanohardness mapping across the oxide/metal interface highlights localized mechanical property variations that may be linked to changes in microstructure and crystal structure within the oxide layer. These findings offer valuable insights into the microstructure and property evolution of Zircaloy-4 during oxidation, contributing to a better understanding of microstructural changes in Zr-based alloys under oxidative environments.

APT↗

Analysis of orientation-dependent deformation mechanisms in additively manufactured Zr using in-situ micromechanical testing: Twinning and orientation gradient

Here, orientation-dependent plasticity in zirconium (Zr)-alloy, which is produced through ultrasonic additive manufacturing (UAM) with subsequent hot-isostatic pressing (HIP), was analyzed by recording electron backscattered diffraction (EBSD) data during in-situ micromechanical testing. Three (0001) plane orientations of a hexagonal close-packed (HCP) structure were analyzed parallel to (1) the rolling direction [X||tensile direction (TD)], (2) the build direction [Z], and (3) the transverse direction [Y]. The analysis revealed that the grains with ~<0001>||TD show twin-dominant plasticity with three variants from {10 $\overline{1}$ 2}<$\overline{1}$ 011>; minor slipping with ($10\overline{1}$ 1)[$\overline{2}$ 113], and ($11\overline{2}$ 2)[$\overline{1} \overline{1}$ 23] pyramidal slip has also been observed. However, grains with orientation ~<0001>⊥TD are mainly sensitive to dislocations assisted plasticity-leading to orientation gradients formation. Furthermore, the neck formation was identified as originating from higher populated micro-crack locations and their association with localized plasticity at defect points in the UAM material. These results demonstrate that minimization of impurities to enable grain growth across prior foil interfaces makes HIP an effective methodology for producing Zr plate, with deformation characteristics expected for conventionally manufactured Zr.

36 MATERIALS SCIENCE↗

Ion irradiation effects on Cr-coated zircaloy-4 surface wettability and pool boiling critical heat flux.

The concept of coating the currently used nuclear fuel cladding (zirconium-based alloy, typically Zircaloy-4 or Zirc-4) with an oxidation preventive layer is a progressing Accident tolerant Fuel (ATF) candidate alloys. The coated Zirc-4-based alloys could be a solution to suppress undesirable fast reaction kinetics with high-temperature steam. Zirc-4 has been the most preferred cladding material in pressurized water reactors (PWRs). Chromium (Cr) based alloys as a coating material provides excellent corrosion protection and good strength and wear resistance. In this paper, we present the surface wettability measurements and pool boiling Critical Heat Flux (CHF) for Cr-coated Zirc-4 claddings pre- and post-exposure to an ion irradiation environment. The wettability measurements, including static contact angle (contact angle, θ) and average surface roughness (surface roughness, Ra), are introduced for samples of different coating thicknesses (5–30 μm thick). The coatings fabricated by the cold spray of Cr-Al particles to 10 mm × 10 mm × 1.95 mm Zirc-4 substrates. Post fabrication, a Pilgering (cold rolling) process, was applied to finalize the coating thickness and resulted in a significant reduction in surface roughness of initially fabricated rough surfaces. The process produced three distinguished samples 5-μm unpolished (as machined), 5-μm, and 30-μm polished (cold rolled). The measurements are presented for the three surfaces and bare Zirc-4 as a baseline surface. The contact angle analyses were implemented in theoretical models from the literature to predict pool boiling CHF. Pool boiling experiments were conducted to measure the pool boiling CHF values and compare them to the predicted values. Scanning Electron Microscope (SEM) images and Energy Dispersive X-ray Spectroscopy (EDS) analysis was performed to characterize the surfaces for better understanding and interpreting the results. The SEM images showed localized surface damage due to ion irradiation. No recognized change in the measured surface roughness due to ion irradiation. The contact angles of irradiated Cr-coated surfaces are consistently higher (10°) than pre-irradiated surfaces. Decreasing the Cr-coating layer thickness resulted in lower contact angle pre- and post- ion irradiation. The predicted pool boiling CHF using the Kandlikar model is in good agreement with the experimentally measured CHF values within ±12% for all samples.

36 MATERIALS SCIENCE↗

Identification and time evolution of thionyl chloride (SOCl 2 ) radiolysis products

In this study, innovative solutions are needed to reduce the amount of high-level waste generated by used nuclear fuel recycling strategies to support the widespread adoption of sustainable nuclear fission energy technologies. To this end, a new sulfur chloride-based process has been developed to recycle zirconium alloy-based materials, which make up a significant fraction of high-level radioactive waste. To support the continued development of this process, we present new data on the potential reaction pathways over time of the products arising from the gamma and electron beam radiolysis of neat thionyl chloride (SOCl 2 ). Interrogation of the gamma irradiated liquid by Raman spectroscopy provided more conclusive identification of the SOCl 2 degradation products, specifically sulfur dichloride (SCl 2 ), molecular chlorine (Cl 2 ), sulfur dioxide (SO 2 ), and sulfuryl chloride (SO 2 Cl 2 ). In comparison, the high dose rate (~10 7 Gy s -1 ) electron beam irradiations formed significantly more degradation products. For both cobalt-60 gamma and electron beam irradiations, the observed degradation products were found to evolve as a function of time post-irradiation via the same reaction pathways, with indication of a solvent regeneration mechanism. These findings are fortuitous for process development, as such a mechanism would be beneficial for process longevity and cost effectiveness.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

High Level Gap Analysis for Accident Tolerant and Advanced Fuels for Storage and Transportation

This initial gap analysis considers proposed accident tolerant fuel (ATF) options currently being irradiated in commercial reactors, since these are most likely for future batch implementation. Also, advanced fuel (AF) options that may be likely for use in advanced reactors are considered. The cladding technologies considered were chromium-coated zirconium-based alloys, FeCrAl, and both monolithic and matrix composite Silicide carbide (SiC). The fuel technologies considered were chromium-doped uranium dioxide fuel, uranium alloys, uranium nitride, and uranium silicide. Numerous national labs, industry, and countries are performing significant testing and modeling on these proposed technologies to establish performance, but at this time none of the prototypes being irradiated have achieved end-of-life (EOL) burnup. There are some testing results after one burnup cycle to verify in-reactor performance, but little data beyond that. As the ATF prototypes acquire more burnup, data will be produced that is relevant to storage and transportation. The DOE:NE Spent Fuel and Waste Science and Technology (SWFST) Storage and Transportation (ST) Control Account will evaluate the performance data as it becomes available for application to the identified gaps for ST.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Fundamental Studies of Tritium Diffusivity in Pure and Defective Zircaloy-4 Getters (Tritium Science Program FY2020 Report)

Zirconium (Zr) alloys are extensively used as tritium (T) getters in nuclear reactors due to their low absorption cross section to thermal neutron, good mechanical properties at high temperature and resistant to corrosion in different environmental conditions. Zr alloys have better thermal properties than many other refractory alloys including stainless steel. The nuclear characteristic is simulated, and T is produced in tritium-producing burnable absorber rods (TPBAR) by the irradiation of neutron flux in pressurized water reactor (PWR). T thus produced diffuses through the pellets and is captured by the getter where it chemically reacts with Zr to form metal hydride (ZrTx). These hydrides are brittle and affect adversely the mechanical properties of alloys. In addition to that, mismatch in the lattice structures after hydride formation creates a stress which needs to be considered in component design and their life evaluation. Therefore, understanding the behavior of T and its species becomes significant as fuel burnup is increased, which leads to increase in hydrogen pickup and oxide formation. The transport of T and its species in pure Zr alloys, and in alloys with certain alloying elements and defects are important to understand in order to enhance the performance of the materials.

36 MATERIALS SCIENCE↗

Report on Thermal Cycling Testing and Bimaterial Weld Development for a Passively Actuated Materials Surveillance Test Article

This report describes the continuation and conclusion of a thermal cycling test program designed to evaluate a family of passively actuated material surveillance test articles. Additionally, the report details a preliminary study of bimetallic welding between 316 stainless steel and TZM (tungsten-zirconium-molybdenum) alloy. The test program was ultimately successful in demonstrating the basic viability of the passively actuated test article concept as well as specifically validating the procedures developed to size specimens to mimic the mechanical response of operating components. The report discusses and summarizes the results of the complete test program which tested samples for more than 1,000 controlled thermal cycles. The weld development program included both stir friction and electron beam joining. The stir friction welding was ultimately successful, while the electron beam welding would require additional development. Ultimately, these test programs pave the way for implementing a materials surveillance program in future molten salt reactors (MSRs). The report conclusions outline a plan for transforming the research and developed described here into a definite set of practices and procedures for a surveillance program at a future operating reactor.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Fuel Performance Analysis of Chromium-Coated Cladding under Burst Conditions

To reduce the oxidation of zirconium-based alloy cladding at high temperatures, accident tolerant fuel systems have been proposed. Of the concepts identified, chromium-coated cladding has been shown to slow oxidation without greatly impacting the fuel system geometry or neutronic performance. To determine how coated-cladding tubes will perform under high-temperature accident conditions, pressurized-tube burst tests have been performed using the Severe Accident Test Station at Oak Ridge National Laboratory. To begin modeling these tubes to better understand how the coating will impact cladding behavior, these burst tests were simulated with the BISON fuel performance code. Cladding tube surface temperatures for the burst test were developed by fitting thermocouple data into axial and azimuthal profiles, while pressure data were compared until cladding failure. The temperatures at failure and the pressure evolution show relatively good agreement between the simulation and experiment results. This is the first step of a larger effort to simulate the cladding deformation process under high-temperature transient conditions and assess the cladding margin to failure more accurately.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Fundamental Studies of Tritium Formation and Diffusivity in Pure and Defective Zircaloy-4 Getters

Zirconium (Zr) alloys have excellent mechanical properties at high temperature and are resistant to corrosion in different environmental conditions. These alloys are used as tritium (3H or T) getters in nuclear reactors due to their low absorption cross section to thermal neutrons and excellent thermo-mechanical properties. Zr alloys have better mechanical and thermal properties than many other refractory alloys including stainless steel. After the irradiation of TPBAR by neutron flux, 3H produced during irradiation diffuses though the pellets and is captured by the getter. 3H chemically reacts with Zr metal to form metal hydride (ZrTx). Hydrides thus formed are brittle and adversely affect the mechanical and thermal properties of alloys. The lattice mismatch at the interface of alloy and metal hydrides also creates a stress that significantly reduces the alloy’s performance and operation life. Therefore, understanding the behavior of 3H and its species becomes significant as fuel burnup is increased, which leads to increase in hydrogen pickup and oxide formation.

Duan, Yuhua↗

Identifying Challenges in Safeguards for Metallic Fuel Fabrication Facilities

As new advanced reactors gain popularity, there is an increasing interest in metallic fuel fabrication for fast reactors. While metallic fuels themselves are not a new idea, as many of the first reactors employed metallic fuels, new designs, compositions, and fabrication methods are appearing throughout the nuclear community. As the interest grows and facilities are constructed, both domestic and international safeguards will need to be heavily involved to support safeguards-by-design (SBD) measures from the start. This work compiles a review of historical and modern fuel types and fabrication methods, fabrication processes, safeguards gaps, and potential safeguards solutions. Metallic nuclear fuel types have been around for many decades and were included in some of the first reactors including the Experimental Breeder Reactor (EBR)-I and -II, the Fermi 1 reactor, the Integral Fast Reactor (IFR), and the Dounreay Fast Reactor (DFR). These reactors used various compositions including pure uranium (U) metal, U-zirconium (Zr) alloys, plutonium (Pu)-aluminum (Al) alloys, U-fissium (Fs) alloys, U-Pu-Zr alloys, and U-molybdenum (Mo) alloys [1, 2, 3, 4, 5]. These small alloying additions are included to improve the material properties of the pure U metal. The alpha-phase U (stable below 661C) suffers elongation in one direction causing grain boundary cracking and increasing creep rate due to irradiation growth, thermal cycling, and preferential crystal orientation. It is ideal to utilize the gamma-phase U (typically stable above 769C) by adding small amounts of alloying elements such as Zr or Mo to stabilize this phase down to room temperature [3]. Additionally, some research has been focused on U with transuranic (TRU) elements present, typically coming from the used fuel recycling process. Including these elements in fast reactor fuel can aid in the reduction of nuclear waste by burning minor long-lived actinides. However, the additions of TRU elements can cause concerns to arise when trying to fabrication or safeguard metallic fuels. A typical metallic fuel element is shown in Figure 1. Sodium is added into the cladding to create a thermal bond between the fuel slug and cladding wall. The fuel slug is then inserted and the end plug is welded on to the top of the fuel element. A gas plenum is left to create a headspace for gaseous fission products to escape rather than continue to build in the fuel itself [1, 5]. Other fuel element geometries exist as well, such as the Lightbridge twisted cruciform geometry shown in Figure 2 [6]. This design allows for better cooling performance and provides room for fuel rod swelling without impacting the fuel rod diameter. There are many different fabrication methods for metallic fuels, which is one of the many benefits of these fuel types. Many of these fabrication methods are relatively easy and cost-efficient. The most popular fabrication method is injection casting, sometimes called vacuum induction melting (VIM), shown in Figure 3 [4, 8, 9, 7, 10]. This method was largely used for EBR-II fuel fabrication. The injection casting system is contained inside of a vessel consisting of a Y2O3-coated graphite crucible surrounded by an induction coil with ZrO2-coated quartz molds suspended above the crucible. The fuel feedstock is placed inside of the graphite crucible and melted using the induction furnace. The induction furnace utilizes a dual frequency with the high frequency melting the feedstock and the low frequency causing stirring of the melted feedstock to form a homogeneous mixture. The mixture is heated to approximately 1600C in an argon environment. The vessel is evacuated and then the quartz molds are lowered into the graphite crucible containing the molten metal and the vessel is repressurized to inject the metal fuel upwards into the molds. The molds are removed and then shattered to release the fuel slugs. This fabrication method was used to fabricate 39,000 metallic fuel pins for EBR-II. While injection casting has been the most common metallic fuel fabrication method throughout the decades, many other methods have been explored including low-pressure gravity casting, microwave casting, continuous casting, centrifugal casting, coextrusion, and many others [11, 12, 8, 13, 14, 15]. Some of these methods aim to mitigate challenges that arise with americium (Am) volatilization during the casting process for TRU-containing fuel feedstocks, an issue with injection casting. Coextrusion is one of the methods explored at the Idaho National Laboratory (INL) and has been utilized for the initial fabrication tests of Lightbridge's unique fuels, as well as other metallic fuels with cladding coextruded. In this process, large billets are formed and machined and then inserted into a molten salt bath for approximately 30 minutes. The billets are then loaded into the extrusion press and extruded. This process can be seen in Figure 4 [15].

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

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↗

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↗