Search NASA⌕ Search

SEARCH · Search NASA

Results for “Density Functional calculations”

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

Thermochemical Data for Furan-based Monomer Candidates for Frontal Ring-Opening Metathesis Polymerization (FROMP)

This dataset includes 471 furan-based monomer candidates for frontal ring-opening metathesis polymerization (FROMP) and relevant thermochemistry as calculated with density functional theory (DFT). The monomer candidates were combinatorically enumerated using Diels-Alder reactions of furan derivatives as dienes and four types of dienophiles (alkenes, alkynes, allenes, and benzynes). Common substituents were enumerated for the dienophile classes, and methyl substitution on the diene was explored. We used the SMILES arbitrary target specification (SMARTS) language to produce monomers and ring-opened structures from diene and dienophile precursor SMILES, and we studied the ring-opening reaction using a homodesmotic equation with ethene. RDKit conformers were initially generated from SMILES, then optimized with GFN2-xTB. The two conformers lowest in energy were then optimized with DFT using the wb97x-D3 functional, def2-TZVP basis set, and def2/J auxiliary basis set. Gibbs free energy corrections were obtained through frequency calculations. Structures with imaginary frequencies below -50 cm^{-1} were excluded from this work, and smaller imaginary modes were flipped to be positive for free energy calculations. Modes below 50 cm^{-1} were treated with the modified rigid rotor approximation, and all thermochemical values were calculated at T=200C. The CSV file contains the monomer SMILES, the free energy of reaction for Diels-Alder addition (G_DA_200), and the enthalpy of the ring-opening reaction (H_RO_200). All energies are given in kcal/mol. An interactive HTML is also included to visualize the monomers in this dataset.

Chua, Lauren↗
DOE OSTI2024

Enhancing Nitrogen Activation in Electrochemical Reduction: The Role of Rare Earth Oxide Surface Configurations

The pursuit of sustainable ammonia synthesis has prompted the exploration of ambient electrochemical nitrogen reduction reaction (e-NRR) as an alternative to the energy-intensive Haber-Bosch process. Here, this study conducted a theoretical investigation into the use of rare earth oxide materials, specifically dysprosium oxide (Dy 2 O 3 ), as potential electrocatalysts for NRR. Utilizing spin-polarized density functional theory calculations, we explored the interaction between Dy 2 O 3 surfaces and nitrogen (N 2 ) molecules, examining the capability of Dy 2 O 3 to adsorb and activate N 2 under ambient conditions. The results indicate that Dy 2 O 3 surfaces exhibit diverse configurations and bonding environments, providing a variety of reactive sites that display different behaviors in N 2 adsorption and activation. The distinctive electronic structure and surface chemistry of a particular Dy 2 O 3 surface configuration were found to significantly enhance the activation of N 2 by promoting charge transfer, which facilitates the NRR process. This research provides deep insights into the mechanistic pathways of N 2 reduction over Dy 2 O 3 , highlighting the surface properties as pivotal in catalysis. These theoretical insights serve as a foundation for the development of novel rare earth-based electrocatalytic materials for efficient ambient e-NRR, potentially transforming ammonia production into a greener and more energy-efficient process.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2021

Activation of N 2 on Manganese Nitride-Supported Ni 3 and Fe 3 Clusters and Relevance to Ammonia Formation

Dual-site models were constructed to represent manganese nitride (Mn 4 N)-supported Ni 3 and Fe 3 clusters for NH 3 synthesis. Density functional theory calculations produced an energy barrier of approximately 0.55 eV for N–N bond activation at the interfacial nitrogen vacancy sites (N v ); also, the hydrogenation and removal of interfacial N is promoted by earth-abundant Ni and Fe metals. Steady-state microkinetic modeling revealed that the turnover frequencies of NH 3 production follow an order of Fe 3 @Mn 4 N ≈ Ni 3 @Mn 4 N > Mn 4 N > Fe >> Ni. Moreover, we present clear evidence that, before NH 3 formation, NH migrates from N v onto the metallic sites. Using N binding energy (BE N ) and the transition-state energy of N 2 activation (E TS ) as descriptors, we concluded that the beneficial effects owing to interfacial N v sites are the most pronounced when BE N is either too strong or too weak while E TS is high; otherwise, excessive N v sites may hinder catalyst performance.

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

Surface Reactivity of Ag-Modified Ceria to Hydrogen: A Combined Experimental and Theoretical Investigation

We investigate the mechanism of H 2 activation on Ag-modified cerium oxide surfaces, of interest for different catalytic applications. The study is performed on thin epitaxial cerium oxide films, investigated by X-ray photoemission spectroscopy to assess the changes of both the Ag oxidation state and the concentration of Ce 3+ ions, O vacancies, and hydroxyl groups on the surface during thermal reduction cycles in vacuum and under hydrogen exposure. The results are interpreted using density functional theory calculations to model pristine and Ag-modified ceria surfaces. Although the reactivity of ceria toward H 2 oxidation improves when a fraction of Ce cations is substituted with Ag, the concentration of reduced Ce 3+ ions in Ag-modified ceria is found to be lower than in pure ceria under the same conditions. This behavior is observed even though the number of surface oxygen vacancies caused by the thermal treatment under hydrogen exposure is larger for the Ag-modified surface. These results are explained in terms of a change of the oxidation state of the surface Ag, which is able to acquire some of the extra surface electrons created by the oxygen vacancies and the adsorbed hydrogen atoms. Our findings provide new insights into the reactivity of Ag-modified ceria, which has been proposed as a promising alternative to platinum electrodes in electrochemical devices.

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

Soft x-ray spectroscopy of ammonium sulfate aerosols at varied pH conditions

Aqueous aerosols are useful model systems for understanding the physical chemistry that takes place in the Earth’s atmosphere, as well as the chemistry of other solar system objects. Ammonium sulfate aerosols are impacted by anthropogenic sources such as farming and shipping and are understood to be important seeds of cloud nucleation in the atmosphere, affecting climate. X-ray photoelectron spectroscopy and near edge x-ray absorption fine structure spectroscopy were used in tandem to probe the surface and bulk properties of aqueous (NH 4 ) 2 SO 4 , respectively. Aerosolized solutions of (NH 4 ) 2 SO 4 , some altered with NaOH to modify pH, were introduced via an aerodynamic lens to a velocity map imaging instrument for detection. The results show that as pH and sodium cation concentration increase, sulfate anion and ammonia become more prevalent at the surface, and the aerosol surface appears to become drier. However, most of the aerosol remains aqueous, with ion concentrations insensitive to the addition of NaOH. Density functional theory calculations were also performed to probe the coordination environment at the aerosol surfaces to understand the underlying phenomena. Finally, the changes seen on the surface of the aerosol underscore the importance of pH in regulating the structure and chemical properties of atmospherically relevant aerosols.

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

Modeling carrier mobility in graphene as a sensitive probe of molecular magnets

Carrier mobility in graphene on a GaAs substrate and its change due to the adsorption of molecular magnets, Mn 12 and [Mn 3 ] 2 , on the surface of graphene, is calculated from first principles. Phonon limited mobility is also calculated for comparison. For Mn 12 adsorption on graphene, the mobility is compared for different organic ligands of Mn 12 (-H, -CH 3 , and - CHCl 2 ), while for [Mn 3 ] 2 dimers, the mobility is calculated for different linkers that yield the ferromagnetic and anti-ferromagnetic configurations of [Mn 3 ] 2 , as well as for different orientations of the molecule. Significant changes in carrier mobility due to the adsorption of the molecules and due to differences in linkers suggest mobility measurement as a possible sensitive probe of magnetic molecules.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Band structure, superconductivity, and polytypism in AuSn 4

The orthorhombic compound AuSn 4 is compositionally similar to the Dirac node arc semimetal PtSn 4 . AuSn 4 is, contrary to PtSn 4 , superconducting with a critical temperature of T c = 2.35 K. Recent measurements present indications for quasi-two-dimensional superconducting behavior in AuSn 4 . Here we present measurements of the superconducting density of states and the band structure of AuSn 4 through scanning tunneling microscopy and angular resolved photoemission spectroscopy (ARPES). The superconducting gap values in different portions of the Fermi surface are spread around Δ 0 = 0.4 meV, which is close to but somewhat larger than Δ = 1.76k B T c expected from BCS theory. We observe superconducting features in the tunneling conductance at the surface up to temperatures about 20% larger than bulk T c . The band structure calculated with density functional theory follows well the results of ARPES. The crystal structure presents two possible stackings of Sn layers, giving two nearly degenerate polytypes. This makes AuSn 4 a rather unique case with a three-dimensional electronic band structure but properties ressembling those of low-dimensional layered compounds.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Enhanced Guest@MOF Interaction via Stepwise Thermal Annealing: TCNQ@Cu 3 (BTC) 2

Confinement of guest molecules in porous materials such as metal organic frameworks (MOFs) promises to deliver emergent properties separate from those of the individual components. Understanding the confinement mechanism is therefore important for the development of new synthesis routes that adjust MOF properties for specific applications. In this work, we developed a new synthetic method to confine guest molecules into MOF pores through a stepwise thermal annealing process, wherein the confinement of 7,7,8,8-tetracyanoquinodimethane (TCNQ) guest molecules into Cu 3 (BTC) 2 (BTC = benzene-1,3,5-tricarboxylic acid) MOF host is used as an example of how novel materials can be created with new physical properties. The stepwise thermal annealing process includes 1) an activation process of pristine Cu 3 (BTC) 2 MOF to maximizes the TCNQ guest loading in the MOF host by effectively removing the residual solvents and 2) post-annealing of the TCNQ infiltrated MOF to enhances the interaction of the confined guest molecules with the MOF host. Obtained experimental results based on thermogravimetric analysis, N 2 gas adsorption, electron microscopy, X-ray diffraction and infrared absorption, combined with density functional theory calculations provide evidence that the use of a stepwise thermal annealing process yields enhancements in the guest loading, packing and interaction between the TCNQ guest and the MOF host. The new hybrid TCNQ@Cu 3 (BTC) 2 system is stable and shows no significant signs of structural degradation even after submersion in water. This is due to the presence of significantly stronger interactions of TCNQ with the frame-work metal ions compared to those of the water molecules competing for the same framework binding sites. It was also found that TCNQ@Cu 3 (BTC) 2 system maintains a significant CO 2 and CH 4 adsorption potential compared to the pristine MOF. The synthetic route developed in this work yields novel guest@MOF hybrid systems that will be useful for many MOF-based applications such as gas separations and chemical sensors performed under humid conditions.

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

Influence of surface chemistry on Li nucleation energetics on graphene-based surfaces

Lithium metal is a promising high-capacity anode material for solid-state batteries, but it typically suffers from poor cyclability. Carbon scaffold hosts have the potential to improve this performance due to their high electronic conductivity and large surface area, which facilitates lithium-ion adsorption and desorption. Scaffold surface chemistry is known to significantly influence performance outcomes, but the details of these interactions are not fully understood. Here, this study employs first-principles simulations to explore lithium transport and nucleation on graphene anodes with various surface chemistries. Using enhanced sampling techniques, ab initio molecular dynamics, and density functional theory calculations, we find that although surface chemistry has a minimal impact on lithium interfacial transport, it influences surface nucleation significantly. Both heteroatom dopants and intrinsic defects lower the nucleation barrier, creating a more favorable environment for lithium nucleation compared to pristine graphene. In addition, our results reveal a complex interplay between surface lithium concentration, lithium transport, and nucleation kinetics. These findings highlight the potential of surface modifications to precisely control nucleation processes on carbon-based anodes and provide design guidance for reducing dendrite formation and improving the cycle life of solid-state batteries.

36 MATERIALS SCIENCE↗
DOE OSTI2020

High-throughput investigation of the formation of double spinels

Spinel compounds, with the general chemical formula AB 2 O 4 , are a wide class of materials, where A and B can be a variety of cations, providing this structure with a great deal of functional flexibility and giving rise to its considerable scientific interest. Recently, a spinel with the general formula ABB'O 4 has been predicted, increasing the possible usability of the spinels due to the higher cation diversity in the so-called double spinel structure. In this study, we use density functional theory calculations to predict if double spinels can be formed between experimentally synthesized normal and inverse single spinels. Our computations reveal that 49 double spinels have negative mixing enthalpies and are thus thermodynamically stable, with most of the stable compounds being formed from one of two distinct cation orderings. We show that the 17 different cations that form the different double spinels have a preferred site, tetrahedral or octahedral, except for Mn, Fe and Co which can occupy both sites interchangeably. We also study the relation between mixing enthalpies and cation-specific properties, as well as ways to classify the double spinels into distinct types and spinel groups depending on the cation ordering and cation oxidation states, respectively. Furthermore, by applying the Sure Independence Screening and Sparsifying Operator (SISSO) approach on the coordination-dependent ionic radii of the elemental constituents, we show that an interplay of local strain and electrostatic dominated terms can be used to separate the double spinels into distinct structural types depending on the cation order and their oxidation states.

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

Low-frequency Raman spectrum of 2D layered perovskites: Local atomistic motion or superlattice modes?

We report the low-frequency Raman spectrum (ω = 10 cm -1 –150 cm -1 ) of a wide variety of alkylammonium iodide based 2D lead halide perovskites (2D LHPs) as a function of A-site cation (MA = methylammonium and FA = formamidinium), octahedral layer thickness (n = 2–4), organic spacer chain length (butyl-, pentyl-, hexyl-), and sample temperature (T = 77 K–293 K). Using density functional theory calculations under the harmonic approximation for n = 2 BA:MAPbI, we assign several longitudinal/transverse optical phonon modes between 30 cm -1 and 100 cm -1 , the eigen displacements of which are analogous to that observed previously for octahedral twists/distortions in bulk MAPbI. Additionally, we propose an alternative assignment for low-frequency modes below this band (<30 cm -1 ) as zone-folded longitudinal acoustic phonons corresponding to the periodicity of the entire layered structure. Furthermore, we compare measured spectra to predictions of the Rytov elastic continuum model for zone-folded dispersion in layered structures. Our results are consistent across the various 2D LHPs studied herein, with energetic shifts of optical phonons corresponding to microscopic structural differences between materials and energetic shifts of acoustic phonons according to changes in the periodicity and elastic properties of the perovskite/organic subphases. This study highlights the importance of both the local atomic order and the superlattice structure on the vibrational properties of layered 2D materials.

2D materials↗
DOE OSTI2024

Coexistence of Interacting Charge Density Waves in a Layered Semiconductor

Coexisting orders are key features of strongly correlated materials and underlie many intriguing phenomena from unconventional superconductivity to topological orders. Here, we report the coexistence of two interacting charge-density-wave (CDW) orders in EuTe 4 , a layered crystal that has drawn considerable attention owing to its anomalous thermal hysteresis and a semiconducting CDW state despite the absence of perfect FS nesting. By accessing unoccupied conduction bands with time- and angle-resolved photoemission measurements, we find that mono- and bi-layers of Te in the unit cell host different CDWs that are associated with distinct energy gaps. The two gaps display dichotomous evolutions following photoexcitation, where the larger bilayer CDW gap exhibits less renormalization and faster recovery. Surprisingly, the CDW in the Te monolayer displays an additional momentum-dependent gap renormalization that cannot be captured by density-functional theory calculations. This phenomenon is attributed to interlayer interactions between the two CDW orders, which account for the semiconducting nature of the equilibrium state. Finally, our findings not only offer microscopic insights into the correlated ground state of EuTe 4 but also provide a general non-equilibrium approach to understand coexisting, layer-dependent orders in a complex system.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Superstructures and magnetic order in heavily Cu-substituted ( Fe 1 – x Cu x ) 1 + y Te

Most iron-based superconductors exhibit stripe-type magnetism, characterized by the ordering vector Q = ($\frac{1}{2},\frac{1}{2}$). In contrast, Fe 1+y Te, the parent compound of the Fe 1+y Te 1–x Se x superconductors, exhibits double-stripe magnetic order associated with the ordering vector Q = ($\frac{1}{2},0$). Here, we use elastic neutron scattering to investigate heavily Cu-substituted (Fe 1–x Cu x ) 1+y Te compounds and reveal that (1) for x ≳ 0.4, short-range magnetic order emerges around the stripe-type vector at Q = ($\frac{1}{2}$ ± δ, $\frac{1}{2}$ ± δ, $\frac{1}{2}$) with δ ≈ 0.05; (2) the short-range magnetic order is associated with a superstructure modulation at Q = ($\frac{1}{3},\frac{1}{3},\frac{1}{2}$), with the magnetic correlation length shorter than that for the superstructure; and (3) for x ≳ 0.55, we observe an additional intergrown phase with higher Cu content, characterized by a superstructure modulation vector Q = ($\frac{1}{3},\frac{1}{3},0$) and magnetic peaks at Q = ($\frac{2}{3},\frac{1}{3},\frac{1}{2}$)/($\frac{1}{3},\frac{2}{3},\frac{1}{2}$). The positions of superstructure peaks suggest that relative to the tetragonal unit cell of Fe 1+y Te, heavy Cu substitution leads to Fe-Cu orderings that expand the unit cell by $\sqrt{2}$ × 3$\sqrt{2}$ times in the ab plane, corroborated by first-principles calculations that suggest the formation of spin chains and spin ladders. Finally, our findings show that stripe-type magnetism is common in magnetically diluted iron pnictides and chalcogenides, despite the varying associated atomic orderings

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Magnetic multilayer edges in Bernal-stacked hexagonal boron nitride

Single-layer h-BN is known to have edges with unique magnetism; however, in the commonly fabricated multilayer AA'-h-BN, edge relaxations occur that create interlayer bonds and eliminate the unpaired electrons at the edge. Recently, a robust method of growing the unconventional Bernal-stacked h-BN (AB-h-BN) has been reported. Here, we use theoretical approaches to investigate the nitrogen-terminated zigzag edges in AB-h-BN that can be formed in a controlled fashion using a high-energy electron beam. We find that these "open"edges remain intact in bilayer and multilayer AB-h-BN, enabling researchers potentially to investigate these edge states experimentally. We also investigate the thermodynamics of the spin configurations at the edge by constructing a lattice model that is based on parameters extracted from a set of first-principles calculations. We find that the edge spins in neighboring layers interact very weakly, resulting in a sequence of independent spin chains in multilayer samples. By solving this model using Monte Carlo simulations, we can determine nm-scale correlation lengths at liquid-N 2 temperatures and lower. At low temperatures, these edges may be utilized in magnetoresistance and spintronics applications.

1-dimensional spin chains↗
DOE OSTI2023

High-frequency phonons drive large phonon-drag thermopower in semiconductors at high carrier density

It has been well established that (i) the thermopower of semiconductors can be enhanced through a phe- nomenon known as the drag effect, and (ii) the drag enhancement involves only low-frequency acoustic phonons and benefits from low electron densities and low temperatures. Using first-principles calculations we show that large drag enhancements to the thermopower are possible at high carrier density even at room temperature and arise from high-frequency acoustic phonons. A fascinating example is cubic boron arsenide (BAs) for which the calculated room temperature drag enhancement of the thermopower exceeds an order of magnitude at a high hole density of 10 21 cm –3 . This remarkable behavior stems from the simultaneously weak phonon-phonon and phonon-hole scattering of the high-frequency phonons in BAs that become drag active at high carrier densities through electron-phonon interactions. Furthermore, this work advances our understanding of coupled electron-phonon nanoscale transport and introduces an unexpected paradigm for achieving large thermopowers.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Deep Understanding of Strong Metal Interface Confinement: A Journey of Pd/FeO x Catalysts

Tuning the atomic interface configuration of noble metals (NMs) and transition-metal oxides is an effective straightforward yet challenging strategy to modulate the activity and stability of heterogeneous catalysts. Herein, Pd supported on mesoporous Fe 2 O 3 with a high specific surface area was rationally designed and chosen to construct the Pd/iron oxide interface. As a versatile model, the physicochemical environments of Pd nanoparticles (NPs) could be precisely controlled by taming the reduction temperature. The experimental and density functional theory calculation results unveiled that the catalyst in the support–metal interface confinement (SMIC) state showed significantly enhanced catalytic activity and sintering resistance for CO oxidation. The constructed Fe sites at the interfaces between FeO x overlayers and Pd NPs not only provided additional coordinative unsaturated ferrous sites for the adsorption and activation of O 2 , thereby facilitating the activation efficiency of O 2 , but also impressively changed the reaction pathway of CO oxidation. As a result, the catalyst followed the Pd/Fe dual-site mechanism instead of the classical Mars–van Krevelen mechanism. For the catalyst in the strong metal–support interaction (SMSI) state, its catalytic activity was seriously suppressed because of the excessive encapsulation of the active Pd sites by FeO x overlayers. Thus, the present study therefore provides detailed insights into the SMIC and SMSI in ferric oxide-supported Pd catalysts, which could guide the preparation of highly efficient supported catalysts for practical applications.

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

Real-space solution to the electronic structure problem for nearly a million electrons

We report a Kohn–Sham density functional theory calculation of a system with more than 200 000 atoms and 800 000 electrons using a real-space high-order finite-difference method to investigate the electronic structure of large spherical silicon nanoclusters. Our system of choice was a 20 nm large spherical nanocluster with 202 617 silicon atoms and 13 836 hydrogen atoms used to passivate the dangling surface bonds. To speed up the convergence of the eigenspace, we utilized Chebyshev-filtered subspace iteration, and for sparse matrix–vector multiplications, we used blockwise Hilbert space-filling curves, implemented in the PARSEC code. For this calculation, we also replaced our orthonormalization + Rayleigh–Ritz step with a generalized eigenvalue problem step. We utilized all of the 8192 nodes (458 752 processors) on the Frontera machine at the Texas Advanced Computing Center. We achieved two Chebyshev-filtered subspace iterations, yielding a good approximation of the electronic density of states. Our work pushes the limits on the capabilities of the current electronic structure solvers to nearly 106 electrons and demonstrates the potential of the real-space approach to efficiently parallelize large calculations on modern high-performance computing platforms.

Chemistry↗
DOE OSTI2024

Accelerating the Discovery of New, Single Phase High Entropy Ceramics via Active Learning

High-entropy ceramics have garnered interest due to their remarkable hardness, compressive strength, thermal stability, and fracture toughness; yet the discovery of new high-entropy ceramics (out of a tremendous number of possible elemental permutations) still largely requires costly, inefficient, trial-and-error experimental and computational approaches. The entropy forming ability (EFA) factor was recently proposed as a computational descriptor that positively correlates with the likelihood that a 5-metal high-entropy carbide (HECs) will form the desired single phase, homogeneous solid solution; however, discovery of new compositions is computationally expensive. If you consider 8 candidate metals, the HEC EFA approach uses 49 optimizations for each of the 56 unique 5-metal carbides, requiring a total of 2744 costly density functional theory calculations. Here, we describe an orders-of-magnitude more efficient active learning (AL) approach for identifying novel HECs. To begin, we compared numerous methods for generating composition-based feature vectors (e.g., magpie and mat2vec), deployed an ensemble of machine learning (ML) models to generate an average and distribution of predictions, and then utilized the distribution as an uncertainty. Here we then deployed an AL approach to extract new training data points where the ensemble of ML models predicted a high EFA value or was uncertain of the prediction. Our approach has the combined benefit of decreasing the amount of training data required to reach acceptable prediction qualities and biases the predictions toward identifying HECs with the desired high EFA values, which are tentatively correlated with the formation of single phase HECs. Using this approach, we increased the number of 5-metal carbides screened from 56 to 15,504, revealing 4 compositions with record-high EFA values that were previously unreported in the literature. Our AL framework is also generalizable and could be modified to rationally predict optimized candidate materials/combinations with a wide range of desired properties (e.g., mechanical stability, thermal conductivity).

36 MATERIALS SCIENCE↗