Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “extreme value”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Global Warming Amplifies Outdoor Extreme Moist Heat During the Indian Summer Monsoon

Because of the climatological prevalence of hot, humid conditions, moist heat extremes are a significant challenge to the health and wellbeing of the people in India. While research has demonstrated the importance of summer monsoon to moist heat in India, impact of monsoon–break and warm spells in modulating extreme moist heat regionally has not been fully investigated. Here we investigate moist heat extremes, as measured by the Wet–Bulb Globe Temperature (WBGT) metric, specifically during monsoon and monsoon–break periods and find that they pose a major threat to physical labor and health relative to other seasons. During the 1951–2020 break period, an increase in area exposed (~42.76 million km 2 ), representing at least 670 million people, to extreme and detrimental WBGT values >31°C occur. Our results imply that future studies on extreme moist heat must pay close attention to the variation of weather systems on synoptic to subseasonal time scales that are superimposed on the seasonal monsoon migration.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Empirical Estimates of Transmission Value using Locational Marginal Prices [Slides]

This study focuses on one potential benefit of transmission infrastructure—congestion relief. It explores historical grid conditions from 2012 through the first half of 2022, and evaluates the marginal value of transmission in facilitating trade within and across regional boundaries by calculating differences in observed nodal wholesale power prices. The study finds that wholesale power prices exhibit stark geographic differences that, in many cases, are stable over time. Many regional and interregional transmission links have significant potential economic value from reducing congestion and expanding opportunities for trade. In fact, many links have hourly average pricing differences in 2021 that exceeded $\$15$/MWh—equivalent to $130 million per year for a 1000 MW link. The value of transmission is correlated with overall energy prices and varies by region and year. Critically, extreme conditions and high-value periods play an outsized role in the value of transmission, with 50% of transmission’s congestion value coming from only 5% of hours. Transmission planners run the risk of understating the benefits of regional and interregional transmission if extreme conditions and high-value periods are not adequately considered. These periods are natural features of actual market operations. As such, the study highlights the need for planners to more-comprehensively assess the value of transmission under both normal and extreme conditions.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2024

Three-dimensional superconductivity induced by an extremely small amount of Li in Li x ⁢SnSe 2

Unconventional superconductivity occurs often in materials with low dimensionality. Furthermore, we report superconductivity observed in layered Li x SnSe 2 with the superconducting transition temperature T c ~ 6 K. Through L⁢i + intercalation in semiconducting SnSe 2 via electrochemical process, Li x ⁢SnSe 2 is formed with an extremely small x value as estimated from the c-axis lattice parameter, carrier concentration, and first-principles calculations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Accelerated Depth Computation for Surface Boxplots with Deep Learning

Functional depth is a well-known technique used to derive descriptive statistics (e.g., median, quartiles, and outliers) for 1D data. Surface boxplots extend this concept to ensembles of images, helping scientists and users identify representative and outlier images. However, the computational time for surface boxplots increases cubically with the number of ensemble members, making it impractical for integration into visualization tools. In this paper, we propose a deep-learning solution for efficient depth prediction and computation of surface boxplots for time-varying ensemble data. Our deep learning framework accurately predicts member depths in a surface boxplot, achieving average speedups of 6X on a CPU and 15X on a GPU for the 2D Red Sea dataset with 50 ensemble members compared to the traditional depth computation algorithm. Our approach achieves at least a 99% level of rank preservation, with order flipping occurring only at pairs with extremely similar depth values that pose no statistical differences. This local flipping does not significantly impact the overall depth order of the ensemble members.

Han, Mengjiao↗
DOE OSTI2026

Heat exposure and maternal stress: evidence from the GRAPHS pregnancy cohort in Ghana

Heat exposure has been linked to psychosocial stress, an established antecedent of perinatal depression; however, evidence on heat-related stress during pregnancy in sub-Saharan Africa remains limited. We analyzed psychosocial stress scores and covariate data from the Ghana Randomized Air Pollution and Health Study, linking daily maximum and minimum shaded wet bulb globe temperature (WBGT) metrics to participants’ stress scores derived from the Crisis in Family Systems-Revised Life Events Questionnaire. We evaluated associations using ordinal logistic regression of pregnancy-average and trimester-average exposures and distributed lag non-linear models (DLNMs) to assess time-varying associations across gestation. Higher average maximum WBGT exposure across pregnancy was associated with increased odds of higher psychosocial stress; each 1 °C increase in maximum WBGT was associated with 64% higher odds of belonging to a higher stress category (OR = 1.64; 95% CI = 1.17–2.31; p = 0.0040). In trimester-average models, higher first-trimester maximum WBGT was also associated with higher stress (OR = 1.44; 95% CI = 1.15–1.81; p = 0.0014). DLNMs suggested that relatively cooler daily maximum WBGT values (25th percentile) were associated with decreased odds of stress in early pregnancy, whereas extreme daily maximum WBGT values (99th percentile) showed a pattern consistent with increased odds of stress in mid-to-late gestation. These findings highlight gestational windows in which heat exposure may influence stress, emphasizing the need for further research into underlying mechanisms and effective interventions to protect maternal mental health in heat-vulnerable settings.

White, Lewis [Columbia University] (ORCID:00090005↗
DOE OSTI2025

Probabilistic Hazard Assessment for Tornadoes, Straight-Line Wind, and Extreme Precipitation at the Savannah River Site

Recent data sets for three meteorological phenomena with the potential to inflict damage on SRS facilities – tornadoes, straight-line winds, and heavy precipitation – are analyzed using appropriate statistical techniques to estimate the occurrence probabilities for these events in the future. Summaries of the results for DOE-mandated return periods and comparisons to similar calculations performed in 2013 by Werth et al. (W2013) are given. Using tornado statistics for i) the combined states of Georgia and South Carolina, and ii) a 2⁰ square area surrounding SRS, we calculated the probability per year of any location at SRS being struck by a tornado (the ‘strike’ probability) and the probability that any point will experience winds above set thresholds. The strike probability was calculated to be 7.04E-4 (1 chance in 1420) per year and tornadic wind speeds for DOE mandated return periods of 50,000 years (corresponding to wind design category 3 (WDC-3), and 125,000 years (meeting WDC-4) (USDOE, 2016) were estimated to be 132 mph and 147 mph, respectively. By contrast, default tornado wind speeds taken from ANSI/ANS-2.3-2011 are somewhat higher: 161 mph for return periods of 50,000 years and 173 mph every 125,000 years (ANS, 2011). Although the ANS and the SRS evaluation used the same basic model (Ramsdell and Rishel, 2007), the region defined in ANS 2.3 that encompasses the SRS also includes areas of the Great Plains and lower Midwest, regions with much higher occurrence frequencies of strong tornadoes. The SRS straight-line wind values associated with various return periods were calculated by fitting existing wind data to a GEV1 distribution and extrapolating the values for any return period from the tail of that function. For the DOE mandated return periods, we expect straight-line winds of 117 mph every 2500 years (the required WDC-3 standard) and 125 mph every 6250 years (WDC-4) at any point within the SRS. These values are similar to those from the ANS-2.3-2011 report, which has wind speeds of 125mph and 133 mph for return periods of 2500 years and 6250 years, respectively. For extreme precipitation, we compared the fits of two different theoretical extreme-value distributions and applied the one that fit the data best for each of several accumulation periods. The DOE mandated 6-hr accumulated rainfall for return periods of 10,000 years (corresponding to precipitation design category 3 (PDC-3) and 25,000 years (PDC-4) were estimated as 9.1 inches and 10.1 inches, respectively. For the 24-hr rainfall return periods of 10,000 years and 25,000 years, total rainfall estimates were 12.02 inches and 13.17 inches, respectively, higher than comparable values provided in the W2013 report.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

The CORDEX-CORE EXP-I Initiative: Description and Highlight Results from the Initial Analysis

Here we describe the first effort within the Coordinated Regional Climate Downscaling Experiment–Coordinated Output for Regional Evaluation, or CORDEX-CORE EXP-I. It consists of a set of twenty-first-century projections with two regional climate models (RCMs) downscaling three global climate model (GCM) simulations from the CMIP5 program, for two greenhouse gas concentration pathways (RCP8.5 and RCP2.6), over nine CORDEX domains at ~25-km grid spacing. Illustrative examples from the initial analysis of this ensemble are presented, covering a wide range of topics, such as added value of RCM nesting, extreme indices, tropical and extratropical storms, monsoons, ENSO, severe storm environments, emergence of change signals, and energy production. They show that the CORDEX-CORE EXP-I ensemble can provide downscaled information of unprecedented comprehensiveness to increase understanding of processes relevant for regional climate change and impacts, and to assess the added value of RCMs. The CORDEX-CORE EXP-I dataset, which will be incrementally augmented with new simulations, is intended to be a public resource available to the scientific and end-user communities for application to process studies, impacts on different socioeconomic sectors, and climate service activities. The future of the CORDEX-CORE initiative is also discussed.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Mechanical behavior and shear band of a powder-metallurgy-fabricated CoCrFeMnNi high-entropy alloy during high strain-rate deformation

High-entropy alloys (HEAs) fabricated by powder metallurgy are considered to exhibit extremely-great potential application values in the field for high-speed impact deformation due to their uniform microstructures and compositions. The mechanical behavior of a powder-metallurgy-fabricated (PM) CoCrFeMnNi HEA was investigated at room temperature with strain rates ranging from 1200 s -1 to 2800 s -1 . It was found that the serration behavior and shear localization of the alloy under dynamic conditions occurred. Results also indicated that the serrated flow took place in the PM CoCrFeMnNi HEA during mechanical deformation at high strain rates (>1000 s -1 ) and became more pronounced as the strain rate increased. Furthermore, it was determined that the yield-strength values of the PM CoCrFeMnNi HEA were sensitive to strain rates, where the values increased from 500 MPa to 700 MPa when the strain rates increased from 1200 s -1 to 2800 s -1 . Shear bands with a width of about 20 μm were generated in the PM CoCrFeMnNi HEA after the shear stress reached a maximum value of 635 MPa, which corresponded to a nominal strain of about 8.71. Nanotwins and ultrafine-equiaxed grains with a diameter of about 150 nm also developed at the core of the shear band during deformation. The results revealed that the strength of the shear band was much higher than that of the matrix grain. Nanotwins in the shear band were composed of a single face-centered-cubic (FCC) structure with the twin plane of ($\bar{1}$1$\bar1$), which could be formed via a rotational dynamic recrystallization (RDR).

36 MATERIALS SCIENCE↗
DOE OSTI2023

Interplay between Facets and Defects during the Dissociative and Molecular Adsorption of Water on Metal Oxide Surfaces

Both surface terminations and defects play a central role in determining how water interacts with metal oxides, thereby setting important properties of the interface that govern reactivity such as the type and distribution of hydroxyl groups. However, the interconnections between facets and defects remain poorly understood, limiting the usefulness of conventional notions such as that hydroxylation is controlled by metal cation exposure at the surface. Here, using hematite (α-Fe 2 O 3 ) as a model system, we show how oxygen vacancies overwhelm surface cation-dependent hydroxylation behavior. Synchrotron-based ambient-pressure X-ray photoelectron spectroscopy was used to monitor the adsorption of molecular water and its dissociation to form hydroxyl groups in situ on (001), (012), or (104) facet-engineered hematite nanoparticles. Supported by density functional theory calculations of the respective surface energies and oxygen vacancy formation energies, the findings show how oxygen vacancies are more prone to form on higher energy facets and induce surface hydroxylation at extremely low relative humidity values of 5 x10 - 5 %. Further, when these vacancies are eliminated, the extent of surface hydroxylation across the facets is as expected from the areal density of exposed iron cations at the surface. These findings help answer fundamental questions about the nature of reducible metal oxide-water interfaces in natural and technological settings and lay the groundwork for rational design of improved oxide-based catalysts.

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

Public values failure associated with Hurricane Ian power outages

Power outages from extreme weather events can diminish community resilience, making it difficult for the areas impacted to bounce back after such events. For socially vulnerable populations, the frequency and duration of power outages can be even more severe. Governments have an obligation to protect public values, or those values that are most fundamental to society, which includes equitable resilience. Using Jørgensen and Bozeman's inventory of public values, this manuscript explores how power outages from extreme weather events create public values failures. More specifically, the manuscript evaluates intraorganizational aspects of public administration during power outages in Florida during Hurricane Ian in 2022. Framing power outages as a public values failure may motivate greater time and effort toward improving equitable access to more resilient power systems.

Entress, Rebecca M.↗
DOE OSTI2026

Volatilized Ammonia Supports Subterranean Ecosystems with Unusual Nitrogen Isotopic Signatures

The Frasassi cave system hosts a robust subterranean ecosystem based on microbial lithoautotrophy. Curiously, acidic biofilms forming above degassing sulfidic cave streams, and the invertebrates that feed on them, are extremely depleted in nitrogen-15 (δ 15 N values less than −20‰). In this study, we tested the hypothesis that these low δ 15 N values result from the volatilization, trapping, and uptake of ammonia degassed from circumneutral streams. We found that dissolved ammonium in the streams had δ 15 N values near +3‰, whereas NH 3(g) in the cave atmosphere above streams exhibited δ 15 N values as low as −27‰, consistent with fractionation by NH 3 volatilization. Extremely acidic condensation droplets on cave walls efficiently trapped airborne NH 3 , accumulating up to 4 mM NH 4 + with δ 15 N values as low as −29‰, thereby confirming volatilized and trapped ammonia as the primary N source to cave wall biofilms and the extensive subsurface ecosystem they support. Airborne ammonia trapping represents a novel mechanism for biological N acquisition and provides abundant N for growth in an extreme subsurface environment that otherwise receives very limited nutrient input.

Nitrogen↗
DOE OSTI2024

A valuation framework for customers impacted by extreme temperature-related outages

Extreme temperature outages can lead to not just economic losses but also various non-energy impacts (NEI), such as increased mortality rates, property damage, and reduced productivity, due to significant degradation of indoor operating conditions caused by service disruptions. However, existing resilience assessment approaches lack specificity for extreme temperature conditions. They often overlook temperature-related mortality and neglect the customer characteristics and grid response in the calculation, despite the significant influence of these factors on NEI-related economic losses. This paper aims to address these gaps by introducing a comprehensive framework to estimate the impact of resilience enhancement not only on the direct economic losses incurred by customers but also on potential NEI, including mortality and the value of statistical life during extreme temperature-related outages. The proposed resilience valuation integrates customer characteristics and grid response variables based on a scalable grid simulation environment. This study adopts a holistic approach to quantify customer-oriented economic impacts, utilizing probabilistic loss scenarios that incorporate health-related factors and damage/loss models as a function of exposure for valuation. The proposed methodology is demonstrated through comparative resilient outage planning, using grid response models emulating a Texas weather zone during the 2021 winter storm Uri. The case study results show that enhanced outage planning with hardened infrastructure can improve the system resilience and thereby reduce the relative risk of mortality by 16% and save the total costs related to non-energy impacts by 74%. In conclusion, these findings underscore the efficacy of the framework by assessing the financial implications of each case, providing valuable insights for decision-makers and stakeholders involved in extreme-weather related resilience planning for risk management and mitigation strategies.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2022

Effects of Salt Concentration, Ionic Strength, and Water Activity on the Growth of a WIPP Archaeal Isolate, Halobacterium sp.

An extremely halophilic archaeon (Halobacterium sp., WIPP strain), isolated from WIPP halite, was tested for its ability to grow in media containing various combinations of sodium and magnesium chloride that influenced overall ionic strength and water activity. Halobacterium sp. grew well in most tested salt combinations, and its growth appeared to depend more upon the relative contribution of sodium to overall ionic strength than on other parameters. Growth rates were highest when media consisted of at least 2 M NaCl and lowest when NaCl dropped below 20% of the ionic strength (usually less than 0.5 M). High magnesium concentrations (1 and 1.5 M MgCl 2 ) did not prevent Halobacterium sp. growth, as long as the sodium contribution was sufficient. Water activity was not a definitive factor in determining growth, because this value did not reach an extreme lower limit in the experimental media.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Effects of Salt Concentration, Ionic Strength, and Water Activity on the Growth of a WIPP Archaeal Isolate, Halobacterium sp.

An extremely halophilic archaeon ( Halobacterium sp., WIPP strain), isolated from WIPP halite, was tested for its ability to grow in media containing various combinations of sodium and magnesium chloride that influenced overall ionic strength and water activity. Halobacterium sp. grew well in most tested salt combinations, and its growth appeared to depend more upon the relative contribution of sodium to overall ionic strength than on other parameters. Growth rates were highest when media consisted of at least 2 M NaCl and lowest when NaCl dropped below 20% of the ionic strength (usually less than 0.5 M). High magnesium concentrations ( 1 and 1.5 M MgCh 2 ) did not prevent Halobacterium sp. growth, as long as the sodium contribution was sufficient. Water activity was not a definitive factor in determining growth, because this value did not reach an extreme lower limit in the experimental media.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2023

Bayesian weighting of climate models based on climate sensitivity

Using climate model ensembles containing members that exhibit very high climate sensitivities to increasing CO 2 concentrations can result in biased projections. Various methods have been proposed to ameliorate this ‘hot model’ problem, such as model emulators or model culling. Here, we utilize Bayesian Model Averaging as a framework to address this problem without resorting to outright rejection of models from the ensemble. Taking advantage of multiple lines of evidence used to construct the best estimate of the earth’s climate sensitivity, the Bayesian Model Averaging framework produces an unbiased posterior probability distribution of model weights. The updated multi-model ensemble projects end-of-century global mean surface temperature increases of 2 °C for a low emissions scenario (SSP1-2.6) and 5 °C for a high emissions scenario (SSP5-8.5). These estimates are lower than those produced using a simple multi-model mean for the CMIP6 ensemble. The results are also similar to results from a model culling approach, but retain some weight on low-probability models, allowing for consideration of the possibility that the true value could lie at the extremes of the assessed distribution. Our results showcase Bayesian Model Averaging as a path forward to project future climate change that is commensurate with the available scientific evidence.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Particle hit clustering and identification using point set transformers in liquid argon time projection chambers

Liquid argon time projection chambers are often used in neutrino physics and dark-matter searches because of their high spatial resolution. The images generated by these detectors are extremely sparse, as the energy values detected by most of the detector are equal to 0, meaning that despite their high resolution, most of the detector is unused in a particular interaction. Instead of representing all of the empty detections, the interaction is usually stored as a sparse matrix, a list of detection locations paired with their energy values. Traditional machine learning methods that have been applied to particle reconstruction such as convolutional neural networks (CNNs), however, cannot operate over data stored in this way and therefore must have the matrix fully instantiated as a dense matrix. Operating on dense matrices requires a lot of memory and computation time, in contrast to directly operating on the sparse matrix. We propose a machine learning model using a point set neural network that operates over a sparse matrix, greatly improving both processing speed and accuracy over methods that instantiate the dense matrix, as well as over other methods that operate over sparse matrices. Compared to competing state-of-the-art methods, our method improves classification performance by 14%, segmentation performance by more than 22%, while taking 80% less time and using 66% less memory. Compared to state-of-the-art CNN methods, our method improves classification performance by more than 86%, segmentation performance by more than 71%, while reducing runtime by 91% and reducing memory usage by 61%.

calibration and fitting methods↗
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↗