Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “tantalum”

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 631 records · Page 35
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Constraints on the Decay of 180⁢𝑚 Ta

180m Ta is a rare nuclear isomer whose decay has never been observed. Its remarkably long lifetime surpasses the half-lives of all other known β and electron capture decays due to the large K-spin differences and small energy differences between the isomeric and lower energy states. Detecting its decay presents a significant experimental challenge but could shed light on neutrino-induced nucleosynthesis mechanisms, the nature of dark matter and K-spin violation. For this study, we repurposed the Majorana Demonstrator, an experimental search for the neutrinoless double-beta decay of 76 Ge using an array of high-purity germanium detectors, to search for the decay of 180m Ta. More than 17 kilograms, the largest amount of tantalum metal ever used for such a search was installed within the ultra-low background Majorana Demonstrator detector array. In this paper we present results from the first year of Ta data taking and provide an updated limit for the 180m Ta half-life on the different decay channels. With new limits up to 1.5 × 10 19 years, we improved existing limits by one to two orders of magnitude. Furthermore, this result is the most sensitive search for a single β and electron capture decay ever achieved.

150 ≤ A ≤ 189↗
DOE OSTI2020

Low friction in bcc metals via grain boundary sliding

In this work, low friction is demonstrated with pure polycrystalline tantalum sliding contacts in both molecular dynamics simulations and ultrahigh vacuum experiments. This phenomenon is shown to be correlated with deformation occurring primarily through grain boundary sliding and can be explained using a recently developed predictive model for the shear strength of metals. Specifically, low friction is associated with grain sizes at the interface being smaller than a critical, material-dependent value, where a crossover from dislocation mediated plasticity to grain-boundary sliding occurs. Low friction is therefore associated with inverse Hall-Petch behavior and softening of the interface. Direct quantitative comparisons between experiments and atomistic calculations are used to illustrate the accuracy of the predictions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Optimizing MeV photon dose generation in laser-solid interactions via stable beam propagation

Bremsstrahlung MeV x rays with high doses are essential for radiography of dense objects. Intense, short-pulse laser-driven MeV x-ray sources offer a promising solution, but understanding how the dose scales with laser intensity is crucial for their development. While theory predicts a square root scaling with laser intensity, experimental verification has been hindered by significant shot-to-shot variations, often exacerbated by long preplasma scale lengths from low laser contrast. Here, we demonstrate low shot-to-shot dose variation and confirm the theoretical scaling using high-contrast OMEGA EP laser interactions with millimeter-thick tantalum foils, achieving MeV x-ray doses up to 5.4 rad at 1 m for E photon > 500 keV. Simulations reveal that the unique focal spot of the OMEGA EP laser enables efficient x-ray generation and reduced variability by facilitating deeper laser penetration into dense plasma, even with shorter preplasma scale lengths. Published by the American Physical Society 2025

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Quantum imaging with positronium-decay-emitted gamma rays

The use of entangled gamma rays from positronium decay for quantum-enhanced imaging of dense materials is demonstrated. Quantum ghost images, where only one of the entangled 511-keV photons interacts with the object, are obtained for tantalum samples of varying density using a 210-ps time-resolution dual detector system and a Na-22 positron source. An analysis comparing both classical and quantum imaging modalities is employed to isolate true 511-keV events from background noise. Image quality is quantitatively assessed using transmission ratios and the Michelson contrast. Quantum-correlated images are found to exhibit superior (up to approximately 1.7x, from 0.49 to 0.83 in the thickest sample measured) contrast compared to classical methods and align well with theoretical expectations. These results suggest that quantum ghost imaging with positronium-based entangled gamma rays could significantly enhance noninvasive imaging of high-density objects, with potential applications in areas such as cargo inspection and security screening.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Buried Dirac Points in Quantum Spin Hall Insulators: Implications for Majorana Kramers Pair-Based Quantum Computing

Quantum spin Hall insulators (QSHIs) host helical electronic edge states that are protected from backscattering due to time-reversal symmetry (TRS). Despite considerable work investigating QSHI edge states, there is still an open question about their unexpected resilience to large magnetic fields where TRS is undoubtedly broken. In this work, we investigate the transport properties of helical edge states in a QSHI-superconductor (QSHI-SC) junction formed by a In⁢As(15 nm)/Ga⁢Sb(5 nm) double quantum well and a superconducting tantalum (Ta) constriction. We observe a robust conductance plateau up to 2 T, signaling resilient edge-state transport. Using a modified Landauer-Büttiker analysis, we find that the zero-field conductance is consistent with 98% Andreev-reflection probability owing to the high transparency of the (In⁢As/Ga⁢Sb)-Ta interface. Such resilience is consistent with the Dirac point for the edge states being buried in the bulk valence band. We further theoretically show that a buried Dirac point does not affect the robustness of the quasi-one-dimensional topological superconducting phase. We find that a buried Dirac point favors the hybridization of Majorana Kramers pairs (MKPs)—predicted to exist in a QSHI-SC constriction—and fermionic modes in the QSHI vacuum edge resulting in extended MKP states, highlighting the subtle role of buried Dirac points in probing MKPs.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Ta Based Damascene Resonators

Transmon qubits have become a leading technology in quantum information science (QIS), with fabrication techniques widely accessible. Historically, devices were predominantly made from aluminum (Al) and niobium (Nb), but state-of-the-art transmons now use tantalum (Ta) achieving coherence times up to and exceeding 0.5 milliseconds [1-3]. The differentiation in performance between Ta and Nb-based devices is attributed to the surface oxidation states affecting Two-Level-System (TLS) noise. One approach to quickly evaluating material and device performance consists of fabricating only the readout resonator of a qubit and evaluating the quality factor which is the strategy employed in this project. Our research focuses on Ta coplanar waveguide (CPW) resonators fabricated using a damascene technique producing pristine metal structures within silicon substrates. This approach increases the dielectric substrate's participation ratio while keeping device edges oxide-free due to entrenchment. This collaborative effort, involving Pacific Northwest National Laboratory (PNNL), NY CREATES/SUNY Poly, and Brookhaven National Laboratory (BNL), has resulted in resonators with Q factors as high as 104 at low power. Ongoing studies aim to further understand and improve device performance by systematic analysis of oxygen trapped between device layers during processing. This research showcases the potential for significant improvements in superconducting qubit performance via the damascene process, heralding new possibilities for the advancement of scalable QIS technology.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2026

Aluminum-Based Superconducting Tunnel Junction Sensors for Nuclear Recoil Spectroscopy

The BeEST experiment is searching for sub-MeV sterile neutrinos by measuring nuclear recoil energies from the decay of 7 Be implanted into superconducting tunnel junction (STJ) sensors. The recoil spectra are affected by interactions between the radioactive implants and the sensor materials. We are therefore developing aluminum-based STJs (Al-STJs) as an alternative to existing tantalum devices (Ta-STJs) to investigate how to separate material effects in the recoil spectrum from potential signatures of physics beyond the Standard Model. Three iterations of Al-STJs were fabricated. The first had electrode thicknesses similar to existing Ta-STJs. They had low responsivity and reduced resolution, but were used successfully to measure 7 Be nuclear recoil spectra. The second iteration had STJs suspended on thin SiN membranes by backside etching. These devices had low leakage current, but also low yield. The final iteration was not backside etched, and the Al-STJs had thinner electrodes and thinner tunnel barriers to increase signal amplitudes. These devices achieved 2.96 eV FWHM energy resolution at 50 eV using a pulsed 355 nm (~3.5 eV) laser. These results establish Al-STJs as viable detectors for systematic material studies in the BeEST experiment.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Zone plate-based extreme ultraviolet mask microscope with through-pellicle imaging capability

Mirror-based and zone plate-based imaging systems are being used in actinic extreme ultraviolet (EUV) reticle review tools. With regard to zone plates, a short working distance is advantageous in terms of the required spectral bandwidth, manufacturability, and potential throughput and imaging performance. Zone plates therefore typically have a short working distance. The industry has adopted the use of an EUV pellicle to protect the photomask. Imaging photomask through-pellicle requires a working distance larger than 2.5 mm. A zone-plate-based EUV mask microscope with a 3-mm working distance has been commissioned at beamline 11.3.2 of the Advanced Light Source. Through-pellicle imaging at an exposure time of 2 s is demonstrated. The instrument achieves an image contrast of 95% on large features on a photomask with a tantalum-based absorber. Imaging down to 45-nm half pitch (mask scale) is demonstrated. A NILS of 2.55 is achieved on 60-nm half-pitch (mask scale) lines and spaces. These results demonstrate that zone-plate-based imaging systems can meet the requirements of an actinic EUV mask review tool in terms of imaging performance and throughput in an instrument compatible with EUV pellicles.

Extreme ultraviolet↗
DOE OSTI2022

Photon Doppler velocimetry surface return study for common surface preparations

Photon Doppler velocimetry (PDV) is a ubiquitous diagnostic method in dynamic compression experiments. It typically requires the roughening of the test surface to ensure scattering of the probe light back into the PDV probe at non-normal angles of incidence. Although surface scattering plays an essential role in the quality of PDV data, surface scattering properties for common materials and preparations have not previously been studied with specific relevance to PDV. Here, we present detailed measurements of the bidirectional reflectance distribution function for four different surface preparations (Scotch-Brite hand and drill roughened, ball-rolled, and grit blasted) and four different materials (copper, aluminum, stainless steel, and tantalum). These measurements employ a conventional PDV probe and obtain a diffraction limited angular resolution with 10 pW accuracy and 2.6% repeatability. In addition to scattering data, we employ scattering theory and simulations to accurately emulate the measured data. We also present a straightforward method to derive the average scattering distribution from surface profilometry and observe several qualitative aspects of the scattering data that may help to optimize PDV signals.

42 ENGINEERING↗
DOE OSTI2023

Origin of age softening in the refractory high-entropy alloys

Refractory high-entropy alloys (RHEAs) are emerging materials with potential for use under extreme conditions. As a newly developed material system, a comprehensive understanding of their long-term stability under potential service temperatures remains to be established. This study examined a titanium-vanadium-niobium-tantalum alloy, a promising RHEA known for its superior high-temperature strength and room-temperature ductility. Using a combination of advanced analytical microscopies, Calculation of Phase Diagrams (CALPHAD) software, and nanoindentation, we investigated the evolution of its microstructure and mechanical properties upon aging at 700°C. Trace interstitials such as oxygen and nitrogen, initially contributing to solid solution strengthening, promote phase segregation during thermal aging. As a result of the depletion of solute interstitials within the metal matrix, a progressive softening is observed in the alloy as a function of aging time. This study, therefore, underscores the need for a better control of impurities in future development and application of RHEAs.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Electronic-grade epitaxial (111) KTaO 3 heterostructures

KTaO 3 heterostructures have recently attracted attention as model systems to study the interplay of quantum paraelectricity, spin-orbit coupling, and superconductivity. However, the high and low vapor pressures of potassium and tantalum present processing challenges to creating heterostructure interfaces clean enough to reveal the intrinsic quantum properties. Here, we report superconducting heterostructures based on high-quality epitaxial (111) KTaO 3 thin films using an adsorption-controlled hybrid PLD to overcome the vapor pressure mismatch. Electrical and structural characterizations reveal that the higher-quality heterostructure interface between amorphous LaAlO 3 and KTaO 3 thin films supports a two-dimensional electron gas with substantially higher electron mobility, superconducting transition temperature, and critical current density than that in bulk single-crystal KTaO 3 -based heterostructures. Our hybrid approach may enable epitaxial growth of other alkali metal–based oxides that lie beyond the capabilities of conventional methods.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Dealing with Measurement Uncertainties as Nuisance Parameters in Bayesian Model Calibration

In the presence of model discrepancy, the calibration of physics-based models for physical parameter inference is a challenging problem. Lack of identifiability between calibration parameters and model discrepancy requires additional identifiability constraints to be placed on the model discrepancy to obtain unique physical parameter estimates. If these assumptions are violated, the inference for the calibration parameters can be systematically biased. In many applications, such as in dynamic material property experiments, many of the calibration inputs refer to measurement uncertainties. Here, we develop a metric for identifying overfitting of these measurement uncertainties, propose a prior capable of reducing this overfitting, and show how this leads to a diagnostic tool for validation of physical parameter inference. The approach is demonstrated for a benchmark example and applied for a material property application to perform inference on the equation of state parameters of tantalum.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Modeling Electrokinetics of Oxygen Electrodes in Solid Oxide Electrolyzer Cells

A microscale model is presented in this study to simulate electrode kinetics of the oxygen electrode in a solid oxide electrolyzer cell (SOEC). Two mixed ionic/electronic conducting structures are examined for the oxygen producing electrode in this work: single layer porous lanthanum strontium cobalt ferrite (LSCF), and bilayer LSCF/SCT (strontium cobalt tantalum oxide) structures. A yttrium-stabilized zirconia (YSZ) electrolyte separates the hydrogen and oxygen electrodes, as well as a gadolinium doped-ceria (GDC) buffer layer on the oxygen electrode side. Electrochemical reactions occurring at the two-phase boundaries (2PBs) and three-phase boundaries (3PBs) of single-layer LSCF and bilayer LSCF/SCT oxygen electrodes are modeled under various SOEC voltages with lattice oxygen stoichiometry as the key output. The results reveal that there exists a competition in electrode kinetics between 2PBs and 3PBs, but 3PBs are the primary reactive sites for single-layer LSCF oxygen electrode under high voltages. These locations experience the greatest oxygen stoichiometry variations and are therefore the most likely locations for dimensional changes. By applying an active SCT layer over LSCF, the 2PBs become activated to compete with the 3PBs, thus alleviating oxygen stoichiometry variations and reducing the likelihood of dimensional change. This strategy could reduce lattice structural expansion, proving to be valuable for electrode-electrolyte delamination prevention and will be the focus of future work.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2023

Crack Growth Rate at Oxygen Electrode/Electrolyte Interface in Solid Oxide Electrolysis Cells Predicted by Experiment Coupled Multiphysics Modeling

Solid oxide electrolysis cell (SOEC) is a very efficient hydrogen production technology, but the cell degradation is a serious limiting factor for its long-term implementation. Oxygen electrode (OE) delamination is reported to be the critical degradation mechanism. In this study, we present a methodology to understand the delamination failure of the OE due to chemical stress in a better perspective. Several OE configurations were tested: baseline strontium-doped lanthanum cobalt iron oxide (LSCF) single layer design and tantalum-doped strontium cobalt oxide (SCT) - LSCF bilayer designs with different SCT loadings. An electro-chemo-mechanical model is developed to associate the electrochemical behavior of the cell with solid mechanics for calculating crack growth of the cell during long term test. The bilayer configuration with SCT 20 wt% has better performance as it survived in the long-term life test with the least crack length. This study implies that an additional nano-coating of SCT over the OE have improved the species transport and oxygen evolution with reduced chemical stress. As the operating current density decreases, it takes longer time for the cell to reach the delamination with the same critical crack length of 6.5 μ m (∼93% of the electrode/electrolyte interface length). Finally, it was concluded that chemical stress plays a significant role in interface delamination failure, however it may not be the only source of stresses at the interface.

08 HYDROGEN↗
DOE OSTI2023

Sensitive Detection of an Erbium Isotope in an Atomic Beam Using Cavity Ring-Down Spectroscopy

Here, we have applied cavity ring-down spectroscopy (CRDS) to the study of the 166Er isotope in an atomic beam. These measurements were realized with an external cavity diode laser tuned to the 400.9 nm atomic transition of erbium and a customized high-finesse ring-down optical cavity under vacuum. Erbium atomic beams of different number densities were generated in a tantalum foil micro-crucible within the cavity. Absorbance values of the 166Er isotope between 3 × 10–6 and 7 × 10–5 were measured with a best-case precision on the order of 10–6, which is remarkable when considering the extreme temperatures at which the measurements were conducted, and the short detection path which is characteristic of collimated atomic beams. Number densities of erbium atoms were inferred to be between 2 × 106 and 6 × 107 cm–3. This work demonstrates for the first time the ability of studying dilute atomic beams of refractory materials with high accuracy utilizing CRDS. In these initial studies, we used erbium as a model system, but we expect to extend the proposed approach to the measurement of isotopes of uranium and plutonium for nuclear non-proliferation applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

NbTaV_BCC_SolidSolution_128atoms_VASP6

We performed density functional theory (DFT) calculations for body-centered-cubic (BCC) structures with 128 lattices sites of solid solution ternary alloys niobium-tantalum-vanadium (Nb-Ta-V). The first-principle code that has been used to run the calculations is the closed-source Vienna Ab-Initio Simulation Package (VASP). 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 8. 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 10,500 randomized atomic structures over 105 chemical compositions. The calculations have been collected on NERSC-Perlmutter 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. Every chemical composition sampled across the composition range in the dataset has its own directory. The convention used to name the directories for ternary alloys is AXBYCZ, where A, B, and C refer to the constituents, and X, Y, and Z 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: 1. INCAR: input file that contains various parameters and settings for controlling the behavior of the electronic structure calculations 2. 0.POSCAR: input file that defines the atomic structure of a system 3. 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 4. 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 5. 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 6. 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 7. 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 8. OUTCARL: 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 9. 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 10. 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 11. vasprun.xml: contains similar information to OUTCAR, but in an xml format. 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.

36 MATERIALS SCIENCE↗
DOE OSTI2023

NbTa_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-tantalum (Nb-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 \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. For Ta, p semi-core states as valence were chosen. 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 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 ternary alloys is AXBYCZ, where A, B, and C refer to the constituents, and X, Y, and Z 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. Subdirectories with name case-N, where N ranges between 1 and 10 (extremes included) contain all the files listed above except 'vaspout.h5'. Subdirectories with name case-N, where N ranges between 41 and 60 (extremes included), contain a duplicate copy of the files listed above except for KPOINTS. The names of the duplicate files end with -bis, and correspond to a second VASP calculation that has converged to a different optimized geometry. 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. and Hafner, J. Ab initio molecular dynamics for liquid metals. Phys. review B 47, 558 (1993). (2) Kresse, G. and 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. and 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. and 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. and Joubert, D. From ultrasoft pseudopotentials to the projector augmented-wave method. Phys. review b 59, 1758 (1999)

36 MATERIALS SCIENCE↗
DOE OSTI2023

TaV_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-vanadium (Tanu-V). 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 element V. For Ta, p semi-core states as valence were chosen. 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. Subdirectories with name case-N, where N ranges between 1 and 10 (extremes included) contain all the files listed above except 'vaspout.h5'. Subdirectories with name case-N, where N ranges between 41 and 60 (extremes included), contain a duplicate copy of the files listed above except for KPOINTS. The names of the duplicate files end with -bis, and correspond to a second VASP calculation that has converged to a different optimized geometry. 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. and Hafner, J. Ab initio molecular dynamics for liquid metals. Phys. review B 47, 558 (1993). (2) Kresse, G. and 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. and Furthmuller, 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. and Furthmuller, 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. and Joubert, D. From ultrasoft pseudopotentials to the projector augmented-wave method. Phys. review b 59, 1758 (1999)

36 MATERIALS SCIENCE↗