Search NASA⌕ Search

SEARCH · Search NASA

Results for “alloy optimization”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Development of Industrial Scale Rare Earth Master Alloys from Their Native Oxides for Magnet Production

The objective of the project is to develop an energy-efficient, reduced cost, and single-step critical metal oxide reduction and alloying methodology for the production of NdFeB and SmCo magnets to facilitate the establishment of a sustainable domestic critical materials supply chain. The method consists of an immiscible molten salt flux layer and a higher density molten metal alloy pool (FeB or Co). The rare earth (RE) oxide and a reductant are added to the molten salt layer, where the reductant first strips the oxygen from the rare earth oxide. Subsequently, the separated RE metal diffuses into the molten metal pool below, creating a RE-saturated master alloy. Towards this end, thermodynamic calculations of the reactions between RE oxides, metallic reducing agent, and molten salt bath chemistry have been performed, and the ideal feeds and conditions for extraction and diffusion to produce master alloys were established. Small-scale and scaled experimentation was performed to validate and optimize the feasibility of viable reactions, temperatures, and process conditions with regards to yield, composition, and process efficiency. Finally, full-scale experiments for the production of NdFeB and SmCo magnets were performed, and their performance characteristics were established.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Texture development in magnetostrictive Fe-Ga alloys processed by laser powder bed fusion

Iron-gallium (Fe-Ga, Galfenol) alloys are promising magnetostrictive materials for actuators, sensors, and energy harvesting, but their performance is highly sensitive to microstructure and texture. Additive manufacturing by laser powder bed fusion (LPBF) offers a pathway to engineer texture and integrate functional materials into complex geometries. Here, we fabricate Fe-Ga alloys (Fe 82.2 Ga 17.8 ) by LPBF of gas-atomized powders and systematically optimize laser power and scan speed to maximize density and control texture. Nearly full-density parts (up to 99.6 %) are achieved within a narrow processing window. Electron backscatter diffraction (EBSD) reveals a strong <100> fiber texture aligned with the build direction and columnar grains up to 1 mm long. Magnetostriction measurements show saturation magnetostriction of 190 ppm in the build direction. Correlating texture data with macroscopic magnetostriction, we estimate intrinsic magnetostriction constants (λ 100 = 228 ppm, λ 111 = 12 ppm), closely matching single crystal-derived values. These results demonstrate the critical interplay between processing, texture, and functional performance in additively manufactured Fe-Ga alloys and establish LPBF as a viable route for high-performance magnetostrictive materials.

Additive manufacturing↗
DOE OSTI2024

Resolving local ordering and structure in Mn x Ge 1- x Te alloys through thermodynamic ensembles of pair distribution functions

Characterizing local bonding environments in complex materials is essential for understanding and optimizing their properties. Equally as important is the ability to predict local motifs as a function of synthesis conditions, enhancing chemists’ ability to design properties into materials. In this study, we present an approach to leverage statistical mechanics to generate temperature- and energy-informed ensemble averaged pair distribution functions (PDFs). This method, which we have named Thermodynamic Ensemble Averages of PDFs for Ordering and Transformations (TEAPOT), utilizes density functional theory (DFT) to relax supercells while incorporating energetic penalties for local order, enabling accurate and computationally efficient analysis of local structure. We apply this method to the neutron PDF measurements of the pseudobinary MnTe–GeTe (MGT) alloy, demonstrating its capability to resolve complex local distortions and chemical ordering. Our results reveal detailed insights into phase transformations and local distortions driven by Mn substitution. For compositions that globally present as rock salt, our analysis reveals that Ge coordination geometry is heavily impacted by synthesis temperature. We propose that high temperature synthesis conditions promote a lowered Ge polyhedra distortion, promoting high charge carrier mobility due to the alignment of local and global structure. Incorporating statistical mechanics and computation into experimental analysis thus guides synthesis of tailored local structure.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Laser powder bed fusion of ODS Fe–Cr–Al (0.3Zr, 0.3Y 2 O 3 ): Unveiling processing-microstructure- mechanical property relationships

Here, this study investigates the fabrication of oxide dispersion strengthened (ODS) Fe-Cr-Al alloys via laser powder bed fusion (LPBF) with strategic additions of 0.3 wt% Zr and 0.3 wt% Y 2 O 3 for enhanced mechanical performance in nuclear applications. Systematic processing parameter optimization yielded three distinct conditions: one low-density product with significant defects and two near-full-density materials with improved consolidation. Comprehensive characterization confirmed single-phase α-ferrite matrix formation with successful incorporation of Y-, Zr-, O-, and C-rich precipitates characteristic of ODS alloys. However, precipitate density remained low (∼10 7 cm −3 ), resulting in sink strength values substantially below optimal levels for radiation resistance. Microhardness values (mid-200s HV) correlated inversely with grain size following the Hall-Petch relationship, indicating grain boundary strengthening as the dominant mechanism rather than precipitation strengthening. The optimized processing conditions achieved excellent mechanical properties with room temperature yield strength of approximately 500 MPa and 30 % elongation, demonstrating superior strength-ductility synergy compared to other additively manufactured ODS materials and performance consistent with literature values for LPBF-processed ODS-FeCrAl alloys. This investigation reveals both the potential and limitations of LPBF processing for ODS Fe-Cr-Al alloys. While successful defect-free fabrication was achieved, results highlight the critical need for systematic optimization of processing parameters and post-processing heat treatments to enhance precipitate density for effective dispersion strengthening and radiation resistance while maintaining additive manufacturing advantages.

Additive manufacturing↗
DOE OSTI2025

Rapid Screening of Liquid Metal Wetting for a Materials Compatibility Library

Wetting behavior of molten metals on solid substrates is a critical phenomenon influencing numerous industrial applications, including welding, anti-corrosion coatings, and metal additive manufacturing (AM). In particular, molten metal jetting (MMJ), an emerging AM technology, requires that the molten metal remain pinned at the nozzle exit. Thus, each new metal requires a specific nozzle material to ensure consistent droplet ejection and deposition, making it important to rapidly identify the appropriate wetting combinations. However, traditional measurements of wetting angles require expensive equipment and only allow one combination of materials to be investigated at a time which can be time consuming. This work introduces a rapid screening method based on sessile droplet experiments to evaluate wetting profiles across multiple metal–substrate combinations simultaneously. This study investigates the wetting interactions of molten Al alloy (Al4008), Cu, and Sn on various ceramic and metal substrates to identify optimal material combinations for MMJ nozzle designs. Results demonstrate that Al4008 achieves wetting on ceramic substrates such as AlN, TiO 2 , and SiC, with varying mechanisms including chemical reactions and weak surface interactions. Additionally, theoretical predictions regarding miscibility gaps and melting point differences were verified for Cu and Sn on refractory metals like Mo and W. Findings from this study contribute to the establishment of a materials compatibility library, enabling the selection of wetting/non-wetting combinations for stable MMJ operation. This resource not only advances MMJ technologies but also provides valuable insights for broader applications such as welding, coating, and printed electronics.

Materials science↗
DOE OSTI2025

An innovative radial gradient material design using hot isostatic pressing for applications in extreme environments

Functionally graded materials (FGMs) are highly advanced continuous or discontinuous structures whose structural and material properties vary along a singular geometric dimension either in the axial or radial direction. Here, the radial gradient FGM design makes for an optimal structural design to incorporate a bi-metallic structure with a copper-based high entropy alloy (Cu-HEA) with good mechanical properties and high irradiation resistance, and Chromium (Cr) with great corrosion resistance. This study focuses on the experimental design of a metal powder loading mechanism to fabricate a bi-metallic radial gradient structure using Cu-HEA and 99.9 % pure Cr metal powders. The powder loading strategy uses custom-designed concentric cylindrical dividers to separate the individual compositions. Two benchtop trial runs were performed for design optimization. The optimized design was then implemented to eventually load the HEA and Cr powders for consolidation via powder metallurgy hot isostatic pressing (PM-HIP). The electron microscopy analysis reveals the successful fabrication of the radial gradient structure with the chemical mapping analysis, demonstrating the gradual composition shift from the HEA at the center to the pure-Cr at the periphery via a three-step gradient.

High Entropy Alloys (HEAs)↗
DOE OSTI2024

Alloying Effects on the Transport Properties of Refractory High-entropy Alloys

Additive Manufacturing (AM) has opened new frontiers for the design of refractory high-entropy alloys (HEAs) for high-temperature applications. The thermal conductivity of the AM feedstock is among the most important thermo-physical properties that control the melting and solidification process. Despite its significance, there remains a notable gap in both computational and experimental research concerning the thermal conductivity of HEAs. Here, we use density functional theory (DFT) to systematically investigate the alloying effects on the transport properties of Ti-Cr-Mo-W-V-Nb-Ta RHEAs, including electrical and thermal conductivities and Seebeck coefficient. The relaxation time of charge carriers is a key underlying parameter determining thermal conductivity that is exceedingly challenging to predict from first principles alone, and we thus follow the approach by Mukherjee, Satsangi, and Singh [Chem Mater 32, 6507 (2022)] to optimize the relaxation time for RHEAs. Here we validated thermal conductivity predictions on elemental solids, binary and ternary alloys, and RHEAs and compared them against thermodynamic (CALPHAD) predictions and our experiments with good correlations. To understand observed trends in thermal conductivity, we assessed the phase stability, electronic structure, phonon, and intrinsic- and tensile strength of down-selected RHEAs. Our electronic structure and phonon results connect well with the observed compositional trends for thermal transport in RHEAs. Our DFT assessment and CALPHAD predictions provide a unique design guide for RHEAs with tailored thermal conductivity, a critical consideration for AM and thermal-management applications.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Surface Segregation of Liquid Metal Plasma-Facing Component Alloys: A ReaxFF Investigation

Engineering liquid metal alloys offers a transformative pathway for plasma-facing components by enabling chemically tailored surfaces that can simultaneously optimize plasma-material interactions, reduce divertor heat flux, and enhance core plasma confinement, thereby advancing the commercial viability of nuclear fusion power plants. This study, employing an atomistic simulation approach, provides direct evidence that incorporating nonmetal surface-active agents (such as O and H, or their combination) enables strong surface segregation. This capability makes tin−aluminum (Sn−Al) and tin−lithium (Sn−Li) alloys, with suitable compositions, good candidates for PFC applications. Specifically, the presence of low-Z solutes (Li, Al) leads to preferential surface enrichment, which imparts low-Z sputtering characteristics, while the Sn solvent maintains thermophysical stability. To systematically examine this behavior, we developed a ReaxFF force field spanning the full Sn/Al/Li/O/H chemistry, validated it against formation energies and elastic constants, and applied it in reactive molecular dynamics simulations at fusion-relevant temperatures. We also introduced an overlapbased segregation index that captures interfacial compositional separation directly from atomistic density distributions. Here, this metric reveals a clear hierarchy of segregation regimes and provides a unified view across all systems studied. Together, these findings establish a mechanistic link between nonmetal chemistry and interfacial structure, providing a predictive framework for designing self-adaptive, low-sputtering liquid metal alloys for fusion applications.

Alloys↗
DOE OSTI2026

Tantalum alloy–based resonators for quantum information systems

Utilizing tantalum (Ta) in superconducting circuits has led to significant improvements, such as high qubit lifetime (T 1 ) and quality factors in both qubits and resonators, suggesting that material optimization plays an important role in the development of superconducting circuits. Thus we here explore superconducting gap engineering in Ta-based devices as a powerful strategy for expanding the range of suitable host materials. By alloying 20 atomic percent (at.%) hafnium into Ta thin films, we achieve a superconducting transition temperature (T c ) of 6.09 K as observed in direct current (DC) transport measurements, reflecting an increase in the superconducting gap. We systematically vary deposition conditions to control film orientation and transport properties of Ta-Hf alloy thin films. We then confirm the enhancement in T c via microwave measurements at millikelvin temperatures. We verify the ~40 increase in T c relative to bare Ta devices, while the loss contributions from two-level systems and quasi-particles remain unchanged in the low temperature regime. These findings emphasize the promise of material engineering in superconducting circuits and point to many potential material candidates for further exploration.

Ta-Hf alloy thin films↗
DOE OSTI2026

First-Principles Evaluation of Proton Hopping in Tetrahedral Oxide Motifs

Proton-conducting oxides (PCOs) are important materials used as ionic conductors for energy conversion technologies. Existing research efforts on PCO optimization and discovery generally focus on complex perovskite-based oxides that require doping and alloying to engineer oxygen deficiency and high proton conductivity. However, the variety of chemical compositions and coordination environments in oxides poses challenges for efficient materials design. In this computational study, we construct a database of simplified motifs to elucidate the relationship between fundamental materials chemistry and proton kinetics. Specifically, we focus on the zincblende crystal structure as a proxy for tetrahedral metal–oxide (M–O) coordination environments. We systematically quantified the effects of cation type, oxidation states, and M–O bond lengths on the proton hopping barrier, and found that strong M–O bonds and metal cations with large and variable oxidation states (e.g., Mo 6+ , V 5+ ) lead to smaller proton hopping barriers. By mapping the candidate cations and their preferred bond geometries onto materials databases such as the Inorganic Crystal Structure Database (ICSD) and Materials Project, we identified real materials containing the corresponding metal–oxide units. In general, we observed good agreement between the calculated proton hopping barriers obtained in real crystal structures and those predicted by our motif database. We also discuss the limitations of our model and possible future extensions to improve its predictive capabilities. Overall, our model provides a first step for the rational design and quick screening of energy-efficient PCOs.

organic↗
DOE OSTI2024

Laser powder bed fusion of oxide dispersion-strengthened IN718 alloys: A complementary study on microstructure and mechanical properties

In this study, two new grades of oxide dispersion strengthened (ODS) Inconel 718 (IN718) alloys were designed by the thermochemical CALPHAD method and produced by laser powder bed fusion (LPBF) technique. Alloys designated as IN718-YF and IN718-YFH, that consist Y 2 O 3 –FeO and Y 2 O 3 –FeO–Hf, respectively, were fabricated with >99.9 % densification using optimized process parameters. CALPHAD calculations were highly consistent with experimental findings, highlighting the formation of Al-containing Y–Ti–O and Y–Hf–O nano-oxides in both alloy types. Texture analyzes revealed no significant texture development in as-built (AB) or heat-treated (HT) alloys. Heat treatment was applied at 1050 °C for 1 h to enhance nano-oxide density. Further, the nano-oxide number density remained similar in IN718-YF while it decreased in IN718-YFH alloy as a result of carbide formation after the heat treatment. Besides, formation of secondary γ' particles was observed in the IN718-YFH/HT alloy. Even though the yield strengths of IN718-YF and IN718-YFH alloys in both AB and HT conditions were similar, the ductility of IN718-YFH was ~50 % less in almost all conditions compared to the ductility of IN718-YF. This has been shown to be as a result of irregular shaped micron-sized Y-Hf-O oxides, martensite formation in AB condition, increased amount of carbides and existence of secondary γ' particles in HT condition in IN718-YFH. High density of stacking faults (SF) forming at the interface of the nano-oxides have been detected in IN718-YF alloys. Besides dislocation/nanoparticle interactions, SFs which are responsible for the delocalization of the deformation improve the ductility of IN718-YF alloys. Overall, high temperature mechanical tests exhibit that both alloys have higher strength with improved ductility compared to the standard IN718 alloys, indicating the contribution of the nano-oxides.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Electronic structure prediction of medium and high entropy alloys across composition space

We propose machine learning (ML) models to predict the electron density — the fundamental unknown of a material’s ground state — across the composition space of concentrated alloys. From this, other physical properties can be inferred, enabling accelerated exploration. A significant challenge is that the number of descriptors and sampled compositions required for accurate prediction grows rapidly with species. To address this, we employ Bayesian Active Learning (AL), which minimizes training data requirements by leveraging uncertainty quantification capabilities of Bayesian Neural Networks. Compared to the strategic tessellation of the composition space, Bayesian-AL reduces the number of training data points by a factor of 2.5 for ternary (SiGeSn) and 1.7 for quaternary (CrFeCoNi) systems. We also introduce easy-to-optimize, body-attached-frame descriptors, which respect physical symmetries while keeping descriptor-vector size nearly constant as alloy complexity increases. Our ML models demonstrate high accuracy and generalizability in predicting both electron density and energy across composition space.

materials science↗
DOE OSTI2025

Nickel promotes selective ethylene epoxidation on silver

Over the last 80 years, chlorine (Cl) has been the primary promoter of the ethylene epoxidation reaction valued at ~40 billion USD per year, providing a ~25% selectivity increase over unpromoted silver (Ag) (~55%). Promoters such as cesium, rhenium, and molybdenum each add a few percent of selectivity enhancements to achieve 90% overall, but their codependence on Cl makes optimizing and understanding their function complex. Here, we took a theory-guided, single-atom alloy approach to identify nickel (Ni) as a dopant in Ag that can facilitate selective oxidation by activating molecular oxygen (O 2 ) without binding oxygen (O) too strongly. Surface science experiments confirmed the facile adsorption/desorption of O 2 on NiAg, as well as demonstrating that Ni serves to stabilize unselective nucleophilic oxygen. Supported Ag catalyst studies revealed that the addition of Ni in a 1:200 Ni to Ag atomic ratio provides a ~25% selectivity increase without the need for Cl co-flow and acts cooperatively with Cl, resulting in a further 10% initial increase in selectivity.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Stress-dependent χ phase transformation in a Ni-based superalloy

The ongoing push to elevate operating temperatures in aerospace gas turbine engines – driven by goals of enhanced fuel efficiency and reduced CO 2 emissions – mandates advancements in the creep resistance of Ni- and Co-based superalloys, which are integral for critical engine components. This study elucidates the role of stress assisted localized phase transformations in the creep properties of these alloys. By leveraging chemo-mechanical coupling, self-healing γ′ precipitates are designed to immobilize planar defects, thereby increasing creep resistance. Employing advanced characterization techniques such as high-resolution Scanning Transmission Electron Microscopy (HR-STEM), in conjunction with atomistic simulations and thermodynamic calculations, novel deformation pathways facilitated by χ local phase transformation (LPT) strengthening have been uncovered; notably, the formation of χ nano-laths through microtwinning and superlattice intrinsic stacking fault (SISF) shearing. This study highlights critical insights into the compositional boundaries necessary for optimizing LPT strengthening while avoiding deleterious bulk formation of η/χ phases. These advancements will guide the design of new alloys maximizing high-temperature creep strength for advanced aerospace applications.

Egan, Ashton J. [Friedrich-Alexander University Er↗
DOE OSTI2025

HfZr_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 hafnium-zirconium (Hf-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. 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 Hf 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 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 OSTI2025

HfTa_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 hafnium-tantalum (Hf-Ta). 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. 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 Hf and Ta. 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 sites. 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,040 randomized atomic structures over 31 chemical compositions. Further methodological and structural information is contained in the dataset README.txt file.

36 MATERIALS SCIENCE↗