Engineering PapersSearch

SEARCH · Engineering Papers

Results for “solid solution”

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 73 records · Page 4

First-principles insights into chromium-induced oxide phases in NiO

The high-temperature oxidation of Ni–Cr alloys leads to complex oxide scales comprising Ni(Cr)O solid solutions, NiCr 2 O 4 spinel, and Cr 2 O 3 corundum within the NiO matrix. Understanding the atomic-scale mechanisms of Cr segregation and precipitate formation is crucial for enhancing oxidation resistance. Here, we employ density functional theory calculations to investigate Cr behavior on NiO(100), (110), and (111) surfaces and in the bulk. Our results reveal that isolated Cr atoms preferentially segregate to the surfaces, stabilizing Ni(Cr)O solid solutions via strong Cr–O bonding, whereas Cr pairs and clusters favor subsurface migration and bulk aggregation, promoting nucleation of NiCr 2 O 4 and Cr 2 O 3 phases. These findings elucidate a size-dependent segregation mechanism linking Cr coordination environments to oxide phase evolution. This atomic-scale insight informs strategies to tailor oxide microstructures and enhance the high-temperature oxidation resistance of Ni–Cr alloys.

Chromium Segregation

Computational Design of Interlayers for Thermally Stable Compositionally Graded Coatings on Nickel Alloys

To extend the service life of Ni-based superalloys, refractory metal coatings are often used. However, direct bonding between metals with dissimilar crystal structure promotes brittle intermetallic phase formation. This work presents a computational thermodynamic framework for high throughput design of functionally graded interlayers to suppress deleterious phases that may form at the interlayer. The Thermo-Calc software package was used to screen candidate metallic interlayer elements based on stability of solid-solution phases. Vanadium was identified as a promising interlayer due to its consistent suppression of intermetallic phases. Temperature-dependent phase diagram mapping between 600 and 1000 °C guided selection of a compositional pathway that significantly reduced intermetallic formation compared to directly joining the Ni-based and Nb refractory alloys. Time–temperature–transformation analysis was performed to assess whether equilibrium-predicted phases are kinetically accessible along regions of the graded path where non-solid-solution phases are not fully suppressed. The methodology was further applied to additional Ni-based alloy and coating systems, illustrating its transferability as an approach for rapid computational design of graded interlayers in dissimilar high-temperature materials.

36 MATERIALS SCIENCE

Fundamental Thermodynamic, Kinetic, and Mechanical Properties of Lithium and Its Alloys

Lithium alloying reactions are beneficial in promoting uniform plating and stripping of lithium metal in all-solid-state batteries. First-principles calculations are performed to predict thermodynamic, kinetic, and mechanical properties of lithium and several important Li–M alloys (M = Mg, Ag, Zn, Al, Ga, In, Sn, Sb, and Bi). While the Li–Mg binary system forms a solid solution, most other lithium–metal alloys prefer stoichiometric intermetallic compounds with common local motifs that enable fast Li diffusion. Lithium and Li-rich alloys exhibit an unusually flat energy landscape along paths that connect BCC to close-packed structures like FCC and HCP, with important implications for mechanical properties. Very low migration barriers for Li diffusion that rival those of superion conductors are predicted, both in pure Li and in Li–M intermetallics. However, vacancy concentration, which is crucial for substitutional diffusion, is predicted to be low in metallic Li and most Li–M intermetallics. Compounds such as B32 LiAl and LiGa as well as D0 3 Li 3 Sb and Li 3 Bi exhibit structural vacancies at higher ends of their voltage windows, which together with low migration barriers leads to exceptionally high Li mobilities. In the Li–Mg solid solution, the addition of Mg is found to decrease the vacancy tracer diffusion coefficient by an order of magnitude.

36 MATERIALS SCIENCE

Multi-material direct ink writing and co-sintering of gadolinium oxide – zirconium oxide components

Methods for fabrication of multi-material or functionally graded ceramic composite architectures are of interest for numerous applications. However, conventional co-sintering of multi-material ceramic parts is a challenge because differences in the sintering behavior of the two materials leads to interfacial strain and, ultimately, component failure. Direct ink writing (DIW) is an extrusion-based additive manufacturing process that excels at multi-material printing because multiple extrusion nozzles can be installed on the same gantry system. Furthermore, the use of DIW as a method to fabricate multi-material ceramic green bodies offers an additional variable for controlling and potentially matching sintering kinetics in the slurry formulation used for two dissimilar feedstocks. In the work documented in this manuscript, we explored two strategies to successfully co-sinter multi-material ceramic oxides: slurry optimization to match sintering kinetics and material gradients to step from one material to another. This manuscript also quantifies the allowable mismatch that avoids part cracking in solid solution forming multi-material systems and discusses best strategies to reduce mismatch during co-sintering. Inks composed of gadolinium oxide (Gd 2 O 3 ) and zirconium oxide (ZrO 2 ), a surrogate for uranium oxide (UO 2 ), were thermally matched, which resulted in a sintering mismatch reduction of over 10%. It was found that ~1% mismatch is tolerable during debind cycles and that ~5% mismatch is manageable during sintering cycles after slurry formulations are optimized to match the sintering behavior. Use of continuous gradients is shown to reduce sintering mismatch, although geometric resolution may be lost due to solid solution formation.

Snarr, Patrick L. [Oak Ridge National Laboratory (

Developing reliable machine learning interatomic potential for Fe–Cr–Ni austenitic alloys

Gaining atomistic understanding of mechanical behavior of heat-resistant structural materials such as Fe–Cr–Ni-based alloys requires an approach with an accuracy close to density functional theory (DFT) that considers the intrinsic properties of the bulk lattice and important defects such as stacking faults, grain boundaries, and surfaces. This work aims to develop reliable machine learning interatomic potential (MLIAP) at cross-scale for Fe–Cr–Ni ternary alloys with a focus on the face-centered-cubic (fcc) solid solution structure. Leveraging the advantages of moment tensor potentials, which typically necessitate a relatively small training dataset and enable rapid calculations using the large-scale atomic/molecular massively parallel simulator package, we ensure the stability and accuracy of the trained potentials. Important defects such as stacking faults, grain boundaries, and surfaces for wide-range compositions are investigated. Structural, thermal, elastic, and defect properties are determined from molecular dynamics simulations comprising several thousand atoms, generated via canonical Monte Carlo simulations guided by the trained potential. The trained potential allows efficient atomic simulations of structural, thermal, and mechanical properties of fcc Fe–Cr–Ni solid solution alloys as a function of composition and temperature. Therefore, the MLIAP approach represents a major advancement from DFT calculations that are limited to small simulation sizes and traditional molecular dynamics simulations using relatively low accuracy potentials. Furthermore, this work outlines a practical foundation for further investigating the structural evolution and mechanical behavior of austenitic stainless steel and nickel-based alloys in a wide array of applications in extreme environments.

Crystal structure

Short-range order and longer-range disorder revealed in germanium–tin alloy thin films by extended x-ray absorption fine structure analysis

Short-range order (SRO) in semiconductor alloys, a relatively under-studied structural phenomenon in which local atomic arrangements differ from those of a random solid solution, is investigated in molecular beam epitaxy (MBE)-grown GeSn thin films. A novel preparation technique is used to pattern these films into microscale ribbons that are released from the substrate for extended x-ray absorption fine structure (EXAFS) analysis. The results indicate a strong SRO in which the first shell around Sn atoms is greatly denuded of Sn atoms relative to the nominal atomic composition of the alloy. This effect is more pronounced than that observed recently in GeSn nanowires grown by chemical vapor deposition. Additionally, the presence of a longer-range disorder detected by EXAFS analysis in the shells of atoms more distant from the absorbers is indicative of the defects and inhomogeneous strain present in the MBE-grown films. The evident existence of the SRO in GeSn alloys deposited by different growth methods and in different strain states suggests that SRO is a general phenomenon in the thin films of this metastable solid solution.

74 ATOMIC AND MOLECULAR PHYSICS

New precious metal containing normal spinels: LiRhRu 1-x Ir x O 4 , LiFeIr 1-x Ru x O 4 , and LiCoIr 1-x Ru x O 4

Large spin-orbit-coupled cations in geometrically frustrated crystal structures have the most suitable setting for exploring novel exotic states of matter. Spinel oxides (AM 2 O 4 ) are well-known examples of geometrically frustrated systems. In this study, we report for the first time the synthesis of compositions LiRhRu 1-x Ir x O 4 (x = 0–0.5), LiFeIr 1-x Ru x O 4 (x = 0–0.5), and LiCoIr 1-x Ru x O 4 (x = 0–0.3) containing precious metal cations on edge-sharing octahedral M-sites, and systematically investigate their magnetic and electrical properties. 57 Fe Mössbauer spectroscopy revealed that iron is trivalent in all LiFeIr 1-x Ru x O 4 solid solutions. Magnetic measurements indicate deviations from theoretical spin-only magnetic moment values, indicating the influence of spin-orbit coupling owing to the presence of 4d and 5d block elements. The LiFeIr 1-x Ru x O 4 series shows spin-glass-like freezing behavior with T g ≈ 20 K, and a small frustration index (f ≈ 1-2), indicating that the frustration originates from site disorder. LiRhRu 1-x Ir x O 4 and LiCoIr 1-x Ru x O 4 exhibit strongly geometrically frustrated magnetism. Electrical resistivity measurements as a function of temperature indicate that all phases are semiconducting. Seebeck coefficient measurements show that LiRhRu 1-x Ir x O 4 and LiFeIr 1-x Ru x O 4 are p-type semiconductors with holes as the major charge carriers. A sign reversal of the Seebeck coefficient indicates both holes and electrons as carriers for LiCoIr 1-x RuxO 4 (x = 0–0.2), but only holes as major carriers for x = 0.3. Here, the Seebeck coefficient and power factor increase drastically in the LiRhRu 1-x Ir x O 4 solid solution with Ir substitution, reaching a maximum of ≈ +125 μV/K and ≈2.3×10 -6 W/mK 2 at ∼650 K for x = 0.5.

Charge carriers

Grain Boundary Segregation Suppresses Local Short‐Range Ordering in Nanocrystalline High‐Entropy Alloys

Multi-principal-element alloys like high-entropy alloys (HEAs) have potential applications in many engineering fields due to their unique mechanical/functional properties. While HEAs are generally considered random solid solutions, recent studies revealed that they are prone to short-range-ordering (SRO) due to the complex multi-pair-wise interactions among the constituent elements. Meanwhile, SROs' evolution can sometimes be deleterious, and it is necessary to have control over their evolution. Examining the AlCoCrFe-Zr model alloy, long-range ordering occurs following the expectation of enthalpic predictions. Advanced characterization techniques—transmission electron microscopy, high-energy synchrotron X-ray diffraction/pair distribution function, and atom probe tomography, reveal that SRO is suppressed in as-milled and GB-decorated NC-(AlCoCrFe)100-xZrx (x = 0–1.5 atomic %). Warren-Cowley coefficient calculations are further used to validate the suppression of SRO. Besides the low segregation enthalpies of Cr, Fe, and Zr, and the high-mixing enthalpy of Cr and Fe, the short diffusion path to GBs due to high-GB density in the NC-HEAs and the higher energy state of the GBs than the matrix promotes GB-segregation that further alters the matrix chemistry and consequently disfavors SRO formation within the matrix. Despite the GB-segregation of Cr, Fe, and Zr, the matrices and GBs remain in a random solid solution.

36 MATERIALS SCIENCE

UO2 microstructural evolutions induced by Ni, Mo, and W dopants for intentional forensics

The concept of tagging nuclear fuel with a chemical barcode to enable forensics analysis across the nuclear fuel cycle is an area of active investigation, particularly to ensure fabrication viability without disrupting current fuel performance. This study explored the feasibility of using Ni, Mo, and W isotopic double-spikes as dopants in UO2 fuel from the perspective of fuel fabrication. Doped UO2 pellets were produced using conventional fuel fabrication processes, including powder mixing, sieving, pressing, and sintering in a reductive atmosphere. Two composition levels, 100 and 1000 ppm, were evaluated for each dopant element with isotopic double-spike configurations. For the Ni system, additional dopant concentrations of 250 and 500 ppm were produced with nonperturbed isotopic ratios. The results demonstrated that successful incorporation of Ni, Mo, and W double-spikes into UO2 pellets occurred with minimal shift in final density or dopant loss during pellet fabrication. Isotopic analysis confirmed the presence of the double-spike signature even when diluted with natural isotopic material in ratio of 1:5 in the fabrication process. Microstructural examinations revealed different impacts on grain size compared with undoped UO2. This study showed that Ni incorporation up to ∼500 ppm promoted moderate grain growth, whereas the Mo and W systems caused grain size reduction at all concentrations. Changes in the UO2 lattice parameter as a function of composition were detected exclusively for Ni up to 500 ppm, indicating that the Ni solid solution was the main factor for the observed grain growth. Insoluble (Mo and W) or supersaturated (Ni > 500 ppm) conditions produced grain size reduction. The Ni-doped pellets in the solution range resulted in a final microstructure within fuel specifications, demonstrating its potential benefits of employing complex dopant systems for potential nuclear forensic applications.

36 MATERIALS SCIENCE

HfTi_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-titanium (Hf-Ti). 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 Ti. 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 Air Force HPC11 cluster using the VASP 6.5.1. Additional methodology and file structure information is available in the dataset README.txt file.

36 MATERIALS SCIENCE

Utilizing integrated neutron diffraction and elastoplastic self-consistent crystal plasticity model to quantitatively assess the strengthening mechanism in Al–12.5Ce and Al–12.5Ce–0.4Mg alloys

An integrated in-situ neutron diffraction and elastic plastic self-consistent crystal plasticity (EPSC-CP) modeling scheme is performed on a binary Al–12Ce alloy and a ternary Al–12Ce–0.4Mg alloys. Using this scheme, the constitutive parameters, i.e. elastic constants and slip system parameters of individual phases can be calibrated which can be used in microstructure-based CP models to predict materials performance. From this study, it is shown that the elastic constants of Al 11 Ce 3 intermetallics calculated from density function theory calculation in the literature are rather accurate. When applied to the EPSC-CP model, the lattice strains of both the binary and ternary alloys are correctly predicted as compared with experiments, and large lattice strain differences between Al (100) plane and Al 11 Ce 3 (010) plane are demonstrated. The slip system parameters calibrated by the scheme shows that the addition of 0.4 wt% Mg in the alloy has little influence on the critical resolved shear stress of initial dislocation glide in the Al matrix which caused plastic yield in the material. This can be explained by the very dilute Mg solute content in the Al solid solution, causing large spacing of Al–Mg lattice misfit sites and little impact on resistance of dislocation glide at initial yield. The 0.4 wt% Mg addition, on the other hand, has a large influence on the hardening term in the slip system parameters, indicating those Al–Mg misfit sites do help dislocation accumulation during the deformation. The impact of dilute Mg addition on the Al slip system parameters is also reflected in the flow behavior of the ternary alloy: little impact on the yield stress, but a large impact on working hardening and tensile strength of the materials which is consistent with the literature.

36 MATERIALS SCIENCE

LABQ3: Bayesian method for quantification of mineral compositions and nano-scale elemental mapping of 3D synchrotron XCT data

Quantitative analysis of mineral compositions is essential in understanding geochemical, mineralogical and environmental processes. Fine-resolution 3D imaging is widely done using synchrotron X-ray computed tomography (XCT), but existing analyses are limited to visualization and segmentation. This paper presents a new method, Linear Attenuation Bayesian Quantitative 3D-mapper (LABQ3), based on the linearity of X-ray attenuation with respect to elemental concentrations. To address the random variability in attenuation measurements, LABQ3 employs Bayesian decision theory to minimize classification error, using reference attenuation distributions from scans of pure mineral standards. To demonstrate LABQ3 and test its performance, we studied precipitated carbonate samples. XCT scans were done at multiple energies using the transmission X-ray microscope (TXM) at beamline 32-ID-C of the Advanced Photon Source at Argonne National Laboratory. The reconstructed 3D images have a voxel size of 20 nm. Analyses revealed rich nano-scale compositional heterogeneity within individual particles. A mixture of calcium and cadmium produced an overall stoichiometric composition of (Ca 0.78 ,Cd 0.22 )CO 3 , with some voxels containing nearly pure CdCO 3 . The addition of zinc led to an overall stoichiometric composition of 33% Ca, 28% Cd, 39% Zn, with a nearly pure CaCO 3 core and compositional zonation through the rim. These compositional gradients are related to temporal sequences of carbonate mineral formation where Cd precipitated at the beginning in (Ca,Cd)CO 3 , while Cd and Zn precipitated at the end in (Ca, Cd,Zn)CO 3 . Results differ from bulk analyses using Inductively Coupled Plasma-Mass Spectrometry (ICP-MS), showing that LABQ3 provides particle-specific insights. LABQ3 distinguishes itself by quantifying chemical compositions along a continuum, making it different from XCT analyses based on segmentation. LABQ3 allows simultaneous acquisition of morphology and chemical composition in 3D, facilitating the interpretation of chemical gradients of trace elements, quantification of solid solution compositions, inferences about temporal sequences of mineral precipitation, and addressing other concerns about solid-phase chemistry.

58 GEOSCIENCES

ThO 2 and Th 1– x U x O 2 Nanoscale Materials and Thin Films for Nuclear Science Applications

This study investigates the dynamics and mechanisms of solution combustion synthesis (SCS) for the preparation of nanoscale ThO 2 and Th 1–x U x O 2 materials, utilizing metal nitrates (Th(NO 3 ) 4 and UO 2 (NO 3 ) 2 ) and acetylacetone (C 5 H 8 O 2 ) as reactants dissolved in a 2-methoxyethanol (C 3 H 8 O 2 ) solvent. By combining thermodynamic calculations, dynamic time–temperature profile measurements with differential scanning calorimetry (DSC) and thermogravimetric analysis (TGA), this research reveals how variations in acetylacetone concentration and uranium content influence the structural parameters of the synthesized oxides. The time–temperature measurements show that the heating rate and maximum combustion temperatures are sensitive to acetylacetone concentration. DSC-TGA results indicate shifts in exothermic peak temperatures as the uranium content changes. The complexation between thorium and acetylacetone emerges as a critical factor, impacting combustion parameters and the structural characteristics of the final products. The uniform distribution of Th and U in the Th 1–x U x O 2 solid solution and the formation of nanoscale particles with strained crystallites are considered essential for the low-temperature densification of these materials for nuclear fuel pellet applications. Additionally, high-quality ThO 2 and Th 1–x U x O 2 thin (100–150 nm) films are successfully synthesized via electrospray deposition of combustible solutions followed by a brief period of heat treatment. Furthermore, these films exhibit excellent structural and morphological uniformity, making them ideal candidates for nuclear measurements, irradiation damage studies, and investigations into the physical properties of both pure and mixed oxides.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH

Microsphere LiMn 0.6 Fe 0.4 PO 4 /C cathode with unique rod-like secondary architecture for high energy lithium ion batteries

LiMn x Fe 1-x PO 4 /C is considered a promising next-generation cathode material with significant commercial potential, inheriting the safety of LiFePO 4 while offering higher energy densities. However, the extremely low conductivity and the Jahn-Teller effect induced by Mn 3+ limit its practical capacity and rate performance. Effective modifications can be achieved through particle nanonization and uniform carbon coating. Here, in this work, we synthesized microspherical LiMn 0.6 Fe 0.4 PO 4 /C cathode materials using a hydrothermal method combined with spray drying carbon coating. The cathode material exhibits a microsphere structure composed of aggregated nanorods with a uniform 3 nm carbon coating, showing good dispersibility, small specific surface area and high tap density. In-situ diffraction analysis showed that expanding the single-phase solid solution region during (de)lithiation can reduce the energy barrier for electron transport, improve the kinetics of the (dis)charge process, and enhance both cycling and rate performance. The initial capacity at 0.1C can reach 155 mAh/g, and the capacity remains at 133.5 mAh/g with a retention rate of 97.1 % after 300 cycles. The synergistic effect of particle nanonization and uniform carbon coating endows the LiMn x Fe 1-x PO 4 /C material with excellent electrochemical performance.

25 ENERGY STORAGE

Helimagnetism in CrBr 2 and CrIBr

In Cr⁢𝑋 2 (𝑋 = Br, I), a Jahn-Teller effect distorts the octahedral configuration of the anions about Cr, resulting in a ribbon chain structure. We previously observed helimagnetism in CrI 2 propagating along the ribbon chains, with an 80–90° rotation per Cr ion. Via neutron scattering, we report that CrBr 2 and the solid solution CrIBr are also helimagnetic, with Néel temperatures of 17 and 12 K, respectively, and signs of intermediate magnetic transitions in both compounds. The helical angle between spins on consecutive Cr ions along the ribbon chains increases with increasing Br substitution, from 89.7° (CrI 2 ) to 116° (CrIBr) to 147° (CrBr 2 ), possibly as a result of decreasing longer-range intrachain spin interactions with the increased electron localization of the lighter anions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND

VZr_BCC_SolidSolution_128atoms_VASP6

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

36 MATERIALS SCIENCE

TaTi_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-titanium (Ta-Ti). 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 Nb. 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 Air Force HPC11 cluster using the VASP 6.5.1. Additional methodology and file structure information is available in the dataset README.txt file.

36 MATERIALS SCIENCE