SEARCH · Engineering Papers

Results for “DFT”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 703 records · Page 39

All sources NASA NTRS DOE OSTI

DOE OSTI2023

On Proton Defects and the Phase Transformation of NiO 2

One approach to high-capacity Li-ion battery cathodes is to increase the Ni content. Unfortunately, Ni-rich materials undergo a phase transition and volume collapse at a high state of charge (SOC), which degrades battery performance. At a high SOC, NiO 2 undergoes a phase transition from the O 3 to the O 1 phase. The O 1 phase appears to be more thermodynamically stable than the O 3 phase. However, density functional theory (DFT) often predicts higher stability for O 3 , depending on the DFT approximation. The energy difference between the O 1 and O 3 phases is on the order of 10 meV/atom. Furthermore, we find that proton defects are more stable in the O 1 phase than in the O 3 phase by about 90 meV per proton. Defects in NiO 2 energetically favor the O 3 -O 1 phase transition observed in a high SOC in Li $x$ NiO 2 .

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Improving and Understanding the Hydrogen Evolving Activity of a Cobalt Dithiolene Metal–Organic Framework

Despite the promising previous reports on the development of electrocatalytic dithiolene-based metal–organic frameworks (MOFs) for the hydrogen evolution reaction (HER), these materials often display poor reproducibility of the HER performance because of their poor bulk properties upon integration with electrode materials. In this work, we demonstrate an in-depth investigation of the electrocatalytic HER activity of a cobalt 2,3,6,7,10,11-triphenylenehexathiolate (CoTHT) MOF. To enhance the durability and charge transport properties of the constructed CoTHT/electrode architecture, CoTHT is deposited as an ink composite (1) composed of Nafion and carbon black. We leverage here the well-established use of catalyst inks in the literature to increase adhesion of the catalyst to the electrode surface and to improve the overall electrical conductivity of the integrated catalyst/electrode. The utilization of the composite 1 leads to a significant improvement in the overpotential (η) to reach a current density of 10 mA/cm 2 (η = 143 mV) compared to prior reports, resulting in the most active MOF-based electrocatalyst for the HER that contains only earth-abundant elements. Extensive density functional theory (DFT) calculations were applied to understand the structure of CoTHT and the mechanistic pathways of the HER. The computational results suggest that an AB stacking geometry is energetically favorable, where one layer is slipped by 1.6 Å relative to the neighboring one along the a and b vectors. Additionally, the DFT calculations indicate that the catalytic cycle likely involves a Volmer discharge step to generate a cobalt hydride, followed by a Heyrovsky step to form a cobalt–H 2 intermediate, and finally the dissociation of H 2 .

DFT calculations↗

DOE OSTI2021

Computational Investigation of the Role of Active Site Heterogeneity for aSupported Organovanadium(III) Hydrogenation Catalyst

A crucial consideration for supported heterogeneous catalysts is the nonuniformity of the active sites, particularly for supported organometallic catalysts. Standard spectroscopic techniques, such as X-ray absorption spectroscopy, reflect the nature of the most populated sites, which are often intrinsically structurally distinct from the most active catalytic sites. Additionally, with computational models, often, only a few representative structures are used to depict catalytic active sites on a surface, even though there are numerous observable factors of surface heterogeneity that contribute to the kinetically favorable active species. A previously reported study on the mechanism of a surface organovanadium-(III) catalyst [(SiO 2 )V III (Mes)(THF)] for styrene hydrogenation yielded two possible mechanisms: heterolytic cleavage and redox cycling. These two mechanistic scenarios are challenging to differentiate experimentally since the kinetic readouts of the catalyst are identical. To showcase the importance of modeling surface heterogeneity and its effect on catalytic activity, density functional theory (DFT) computational models of a series of potential active sites of [(SiO 2 )V III (Mes)(THF)] for the reaction pathways are applied in combination with kinetic Monte Carlo (kMC) simulations. Computed results were then compared to the previously reported experimental kinetic study: (1) DFT free-energy reaction pathways indicated the likely active site and pathway for styrene hydrogenation, a heterolytic cleavage pathway requiring a bare tripodal vanadium site. (2) From the kMC simulations, a mixture of different bond lengths from the support oxygen to the metal center was required to qualitatively describe the experimentally observed kinetic aspects of a supported organovanadium(III) catalyst for olefin hydrogenation. This work underscores the importance of modeling surface heterogeneity in computational catalysis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Mechanistic Understanding of Methane Combustion over Ni/CeO 2 : A Combined Experimental and Theoretical Approach

Catalytic oxidation of methane (CH 4 ) over nonprecious Ni/CeO 2 catalysts has received a lot of attention due to the large natural gas reserves found in North America and the prohibitive cost of palladium-based catalysts, commonly used for CH 4 oxidation. However, the catalytic mechanism of CH 4 oxidation over Ni/CeO 2 still remains unclear. Moreover, the parameters affecting the reaction rates, the interaction between nickel and CeO 2 , and the reaction intermediates are still not well understood. In this study, kinetic model fitting, CH 4 temperature-programmed reduction-mass spectroscopy (CH 4 TPR-MS), in situ diffuse reflectance infrared Fourier transform spectroscopy (DRIFTS), and density functional theory (DFT) calculations were combined to elucidate the mechanism of complete oxidation of CH 4 over Ni/CeO 2 . CH 4 TPR-MS showed that the complete oxidation of CH 4 over Ni/CeO 2 requires 55–120 °C lower compared to bare CeO 2 or Ni/quartz sand; complete oxidation of CH 4 took place when the surface oxygen species were abundant, while partial oxidation products (CO, H 2 ) were formed when the oxygen species were depleted. In situ DRIFTS showed that CH 3 , CH 2 , CO, and CO 2 were formed after CH 4 activation over Ni/CeO 2 , while CH 3 O species were not observed. Combining those findings with kinetic model fitting, a redox Mars–van Krevelen (MvK) mechanism showed the best description of the experimental observations. The MvK mechanism involves the reaction of dissociated oxygen species with gas-phase CH 4 while water inhibits the reaction rate by adsorbing on the oxidized sites. Moreover, CH 4 activation leads to the reduction of the active sites and oxygen vacancy formation followed by reoxidation of the active sites by gas-phase O 2 . A CH 4 oxidation reaction pathway over Ni/CeO 2 is proposed by DFT calculations. In summary, the findings shown here suggest that CH 4 oxidation over Ni/CeO 2 follows a redox MvK mechanism and provides guidance for the rational design of non-precious-metal catalysts for CH 4 oxidation reactions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Shape-Selective Synthesis of Alkylcyclopentenyl Cations in Zeolites and Spectroscopic Distinction of Constitutional Isomers

Alkylcyclopentenyl cations belong to the long-lived intermediates that make up the “hydrocarbon pool” during the catalytic conversion of methanol on zeolites, and recent works show that such cations contribute to olefin and aromatics formation. From liquid phase chemistry, two types of alkylcyclopentenyl cations are known and distinguished by the substituent at the central carbon (C-2) of the allylic system: a more stable type with a methyl group at the C-2 and a less stable type with hydrogen at the C-2. The following three linked objectives are pursued: (i) IR spectroscopic distinction of the different substitution patterns at the allylic system of alkylcyclopentenyl cations, (ii) the role of the zeolite framework in determining the substitution pattern, and (iii) identification of alkylcyclopentenyl cations from precursors relevant to methanol-to-olefins conversion. UV–vis and IR spectroscopy are applied in situ to characterize alkylcyclopentenyl cations produced by adsorption of a pentaalkylcyclopentadiene or by adsorption and thermally induced cyclization of 2,6-dimethyl-2,4,6-octatriene. Prior knowledge of electronic spectra is combined with DFT-computed and experimental IR spectra to establish the frequency of the C–H vibration of the C-2 hydrogen-substituted type, and a characteristic red shift of the asymmetric allylic stretching vibration (Δν ≈ -20 to -30 cm –1 ) after replacing hydrogen by methyl at the central carbon of the allylic system. Although DFT demonstrates that both types of ions fit into medium- and large-pore zeolites and that the C-2 methyl-substituted type is thermodynamically favored even in the pores of the considered zeolites, formation of alkylcyclopentenyl cations is found by UV–vis and IR spectra to be shape-selective. The bulkier C-2 methyl-substituted type is detected in large-pore zeolites (MOR, BEA) and in the intersections of medium-pore zeolites (MFI), whereas in channels of medium size (TON), the less bulky C-2 hydrogen-substituted type is exclusively formed. Finally, the type of ion formed and its overall size are dictated by the zeolite framework and are independent of the precursor; the same type of alkylcyclopentenyl cation as found through cyclization of dimethyloctatriene could be generated from ethene.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Correlating Experimentally Determined Hydrogen Binding Energy with Hydrogen Evolution Activity over Metal Monolayers on Molybdenum Nitride

It is well established that hydrogen binding energy (HBE) is a key descriptor for hydrogen evolution reaction (HER) activity, and such a relationship is a useful tool for searching efficient and cost-effective HER catalysts. However, in almost all cases, the HBE values are obtained from density functional theory (DFT) calculations. Here in this study, temperature programmed desorption (TPD) was used to experimentally determine the HBE values of metal monolayers supported on molybdenum nitride (Mo 2 N), and electrochemical measurements were performed on the same surfaces. Combined DFT and kinetic Monte Carlo (kMC) simulations were used to validate the trend observed with TPD and the electrochemical HER activity. Depositing one monolayer of Pt on Mo 2 N led to similar HBE values seen for bulk Pt, and electrochemical measurements showed that monolayer Pt on Mo 2 N had similar HER activity to that of bulk Pt. Similar studies were also performed for monolayer Pd on Mo 2 N to confirm the correlation. This work also demonstrates that Mo 2 N is a promising support to reduce precious metal loading in HER catalysts.

08 HYDROGEN↗

DOE OSTI2025

From Micro-environments to Macroscopic Effects: How the Alkaline Hydrogen Evolution Reaction Drives Cu Cathodic Corrosion

Cathodic corrosion of copper (Cu) has posed a significant challenge for over a century, impeding various technological progresses such as electrochemical conversion of CO 2 (eCO 2 RR) into fuels and other value-added carbon products. Here, in this study, employing a combined Density Functional Theory (DFT) and kinetic Monte Carlo (kMC) simulation approach, we delve into the atomistic level mechanism driving this phenomenon in Cu. Our hypothesis posits the pivotal role of alkaline hydrogen evolution reaction (HER) in facilitating cathodic corrosion in Cu. We rigorously develop a pH-dependent hydroxide (OH) adsorption mechanism and calculate the equilibrium OH coverage (𝜃 OH ) at varying pH levels, the thermodynamic stability of subsurface oxygen (O sub ), as well as the Cu-vacancy mediated diffusion of subsurface oxygens (O sub ). Through comprehensive analysis, we establish correlation among various microenvironments, including oxygen diffusion in subsurface layers, pH-dependent OH adsorption, and Cu dissolution into the electrolyte as (Cu-OH) complexes. Furthermore, our investigation explores the correlation between surface coordination environment of active sites and cathodic corrosion of Cu. Finally, by integrating DFT-derived thermodynamic data into a kMC model, we successfully predict the formation of experimentally observed corrosion pits on Cu-surfaces. This combined approach not only advances our fundamental understanding of Cu cathodic corrosion but also offers insights crucial for developing effective corrosion mitigation strategies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2026

Direct Ab Initio Simulation of the Synthesis of BaZrO 3 and the Microstructure Impacts on Proton Transport

Controlling and predicting the processing-structure-performance relationship in functional materials is a grand challenge in materials science, with important implications for a wide range of emerging applications; a high fidelity understanding of the performance impact of microstructures formed under synthesis conditions is required to develop advanced materials, such as solid-state fuel cells and electrolyzers. Using the ceramic BaZrO 3 as a case study, we directly simulate the synthesis and investigate how proton transport is dictated by microstructures. We develop a framework that couples density functional theory (DFT), machine-learning interatomic potential (MLIP) driven molecular dynamics, and grand canonical Monte Carlo to perform large-scale, microstructure-resolved, atomistic simulations of proton transport in experimentally representative polycrystalline structures. Our fully ab initio approach, using a MLIP as a proxy for DFT, allows us to quantify the competition between two distinct diffusion mechanisms: one associated with grain-boundary regions and another within grains. When the impacts of grain boundaries are taken into account, proton transport exhibits substantial deviation from the bulk oxide limit. This addresses long-standing discrepancies between theory and experiments. Our integrated approach provides atomistic insight into microstructure-dependent proton pathways in BaZrO 3 and establishes a general protocol for predicting processing-structure-performance relationships.

organic↗

DOE OSTI2024

Kinetic Study of Polyvinyl Chloride Pyrolysis with Characterization of Dehydrochlorinated PVC

In this paper, we study the kinetics of polyvinyl chloride (PVC) decomposition using a combination of experimental and computational approaches. Here, we develop a simplified kinetic model that contains only two steps: dehydrochlorination of PVC and further decomposition of the PVC residue. The model is consistent with density functional theory (DFT) calculations and experimental data. Dehydrochlorination is an autocatalytic reaction that starts with a tertiary chloride (Cl) and generates hydrogen chloride (HCl) and benzene as the main products. Benzene and HCl formation rates showed similar trends, indicating that HCl likely catalyzes a homolytic carbon-carbon (C-C) bond cleavage, which gives rise to benzene and an aliphatic fragment. We characterized the structure of dehydrochlorinated PVC (PVC residue) by using thermal gravimetric analysis (TGA), Fourier-transform infrared spectroscopy (FTIR), and nuclear magnetic resonance spectroscopy (NMR). FTIR and NMR results indicate that the PVC residue contains 20% quaternary carbon content, indicating a high concentration of cross-linked molecules. We predict that the most probable structure in the cross-linked centers of the PVC residue is cyclohexadiene, which is supported by DFT calculations, FTIR, and NMR.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Recreating Giants Impacts in the Laboratory: Shock Compression of MgSiO 3 Bridgmanite to 14 Mbar

Understanding giant impacts requires accurate description of how extreme pressures and temperatures affect the physical properties of the constituent materials. Here, we report shock experiments on two polymorphs of MgSiO 3 : enstatite and bridgmanite (perovskite) crystals.We obtain pressure-density shock equation of state to 14 Mbar and more than 9 g/cm 3 , a 40% increase in density from previous data on MgSiO 3 . Density-functional-theory molecular dynamics (DFT-MD) simulations provide predictions for the shock Hugoniot curves for bridgmanite and enstatite and suggest that the Grüneisen parameter decreases with increasing density. The good agreement between the simulations and the experimental data, including for the shock temperature along the enstatite Hugoniot reveals that DFT-MD simulations reproduce well the behavior of dense fluid MgSiO 3 .We also reveal a high optical reflectance indicative of a metal-like electrical conductivity which supports the hypothesis that magma oceans may contribute to planetary magnetic field generation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2022

Topological electronic structure of YbMg 2 Bi 2 and CaMg 2 Bi 2

Zintl compounds have been extensively studied for their outstanding thermoelectric properties, but their electronic structure remains largely unexplored. Here, we present a detailed investigation of the electronic structure of the isostructural thermopower materials YbMg 2 Bi 2 and CaMg 2 Bi 2 using angle-resolved photoemission spectroscopy (ARPES) and density functional theory (DFT). The ARPES results show a significantly smaller Fermi surface and Fermi velocity in CaMg 2 Bi 2 than in YbMg 2 Bi 2 . Our ARPES results also reveal that in the case of YbMg 2 Bi 2 , Yb- 4 f states reside well below the Fermi level and likely have a negligible impact on transport properties. To properly model the position of 4f-states, as well as the overall electronic structure, a Hubbard U at the Yb sites and spin-orbit coupling (SOC) have to be included in the DFT calculations. The theoretical results reveal that both materials belong to a Z 2 topological class and host topological surface states around EF. Due to the intrinsic hole doping, the topological states reside above the Fermi level, inaccessible by ARPES. Our results also suggest that in addition to SOC, vacancies and the resulting hole doping play an important role in the transport properties of these materials.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Bonding and reactivity of isostructural uranyl and neptunyl peroxide phases

Understanding the reactivity of actinide peroxides is critical for predicting the behavior of spent nuclear fuel in radiolytic environments. Herein, we report the synthesis and characterization of a lithium neptunyl(VI) hydroxo peroxo compound (LiNp), which is isostructural to the uranyl analogue (LiU). Single-crystal X-ray diffraction reveals that LiNp contains both [NpO 2 (O 2 ) 3 ] 4- and [NpO 2 (OH) 4 ] 2- units stabilized by Li + and hydrogen bonding and Raman spectroscopy shows systematic redshifts in neptunyl vibrational modes relative to uranyl. DFT calculations highlight the importance of secondary coordination in reproducing vibrational and structural features, but challenges remain with correctly modeling strong sigma donors. Solid-state EPR spectroscopy and DFT confirm superoxide stabilization within LiU and calculations suggest favorability of the analogous radical species in LiNp. Solution state EPR spectroscopy with the radical spin trap 5-tert-butoxycarbonyl-5-methyl-1-pyrroline N-oxide (BMPO) reveal evidence of superoxide in the LiU and LiNp phases and suggests stabilization of superoxide within actinyl triperoxide complexes, forming [AnO 2 (O 2 ) 2 (O) 2 • ] 3- .

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2025

Ba 1−x Sr x FeO 3−δ as an improved oxygen storage material for chemical looping air separation: a computational and experimental study

Chemical looping air separation (CLAS) is a promising technology to generate oxygen-rich gas streams to enable efficient carbon dioxide capture during fossil fuel combustion or gasification. CLAS relies on the capture and release of oxygen from the atmosphere using the redox properties of an oxygen-selective solid oxide carrier. This study investigates the redox characteristics of Ba 1−x Sr x FeO 3−δ (0.0 ≤ x ≤ 0.417, 0.0 ≤ δ ≤ 0.5) using a combination of density functional theory (DFT) calculations and experimental verification using X-ray diffraction, thermogravimetric analysis, and oxygen-temperature-programmed desorption. The DFT computed energies of the Ba 1−x Sr x FeO 3−δ perovskites reveal a composition-dependent transition from hexagonal to cubic phases as the Sr-concentration or oxygen vacancy concentration increases. Oxygen vacancy formation energies of the cubic perovskites are found to be lower than those of their hexagonal counterparts. A low oxygen diffusion barrier of ∼1 eV combined with the thermodynamic preference of Ba 1−x Sr x FeO 3−δ compositions that form in a cubic phase suggests them as promising candidates for oxygen storage applications. The experimental results corroborate this finding by identifying Ba 0.75 Sr 0.25 FeO 3−δ in the cubic phase as an optimal composition offering low-temperature oxygen storage capacities comparable to that of the state-of-the-art Sr 0.75 Ca 0.25 FeO 3−δ perovskite oxygen storage material at 325 °C and 350 °C.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Magnetocaloric effect observations near room temperature in few-layered chromium telluride (Cr 2 Te 3 )

Transition metal telluride compositions are explored extensively for their unique magnetic behavior. Few-layered chromium telluride (Cr 2 T e3 ) exhibits a near-room-temperature phase transition, where the material can be effectively used in applications such as magnetic refrigeration. Compared to existing magnetocaloric materials, Heusler alloys, and rare-earth-based alloys, the large-scale synthesis of mechanically exfoliated Cr 2 Te 3 involves less complexity, resulting in a stable composition. Compared to existing tellurides, Cr 2 Te 3 exhibited a large change in magnetic entropy (|ΔS M |) of 1.88 J k g−1 K −1 at a magnetic field of 4 T. A refrigeration capacity (RC) of ∼82 J kg −1 was determined from the change in magnetic entropy versus temperature curve. The results were comparable with those for existing Cr-based compounds. First-principles density functional theory (DFT) confirmed the magnetic properties of Cr 2 Te 3 , including a near-room-temperature Curie temperature, T C , consistent with experimental results. Here, structural transition was also observed using first-principles DFT, which is responsible for the magnetic behavior.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Density functional theory calculation of helium segregation and decohesion effect in W110/Ni111 interphase boundary

Density functional theory (DFT) is employed to study the preferential distribution and decohesion effect of He in a W-NiFe composite consisting of W particles embedded in an Ni-based solid solution matrix. A slab containing {110}<100>W//{111}<110>Ni interface is employed as a surrogate model for the W-NiFe system. Firstly, the fracture energy of the W/Ni interphase boundary (IB) (4.37 J/m 2 ) is higher than the cleavage energy of Ni{111} (3.82 J/m^2) and lower than the cleavage energy of W{110} (6.60 J/m 2 ). The comparison shows that the cohesion of the IB is stronger than the Ni{111} planes that are away from the IB. However, the cohesion between the Ni{111} planes adjacent to the IB is found to be the weakest in this system, with a cleavage energy of 3.11 J/m 2 . Subsequently, formation energy of He is calculated in the Ni slab, W slab, and various interstitial sites in the IB. The DFT calculations show that He is significantly more stable in Ni than in W, by about 1.75 eV. Interestingly, He does not prefer to segregate at the IB as compared to bulk Ni. Nevertheless, it prefers to segregate to the region between the Ni{111} planes adjacent to the IB, and decreases the cohesion of the already weakest region. Based on an estimated amount of He gas production in 5 years under first wall neutron irradiation (neutron flux of 1.04x10 15 n/cm 2 /s), the He will decrease the cleavage energy of the weakest region by 21.2%, 15.4%, and 12.2% at 800, 1000, and 1200 C, respectively.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2025

Temperature-dependent lattice dynamics in polycrystalline ZrN

Zirconium nitride holds promise for improving nuclear reactor components, where thermal conductivity is important, highlighting the need for a detailed understanding of its phonons. However, the phonon dispersion derived from previous inelastic neutron scattering data is in disagreement with modern theoretical calculations, pointing to the necessity for further experimental validation. Here, in this study, we use inelastic neutron scattering experiments on polycrystalline ZrN at 6, 100, 300, 400, and 600 K, and density functional theory (DFT) calculations to study the momentum-resolved scattering and the phonon density of states. Both experimental and theoretical results reveal a phonon cutoff energy near 65 meV, in stark disagreement with previous single-crystal neutron scattering measurements indicating a cutoff near 75 meV. The 65 meV cutoff energy is, however, in reasonable agreement with previous measurements on ZrN 0.9 . Consistency between our theoretical and experimental results confirms that ZrN is a wide phonon band gap system and that DFT provides a reliable description of its phonons.

DFT↗

DOE OSTI2021

Molecular Modes of Attosecond Charge Migration

Here, first-principles calculations are employed to elucidate the modes of attosecond charge migration (CM) in halogenated hydrocarbon chains. We use constrained density functional theory (DFT) to emulate the creation of a localized hole on the halogen and follow the subsequent dynamics via time-dependent DFT. We find low-frequency CM modes (~ 1 eV) that propagate across the molecule and study their dependence on length, bond order, and halogenation. We observe that the CM speed (~ 4Å/ fs) is largely independent of molecule length, but is lower for triple-bonded versus double-bonded molecules. Additionally, as the halogen mass increases, the hole travels in a more particlelike manner as it moves across the molecule. These heuristics will be useful in identifying molecules and optimal CM detection methods for future experiments, especially for halogenated hydrocarbons which are promising targets for ionization-triggered CM.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2022

Solutes that reduce yield strength anisotropies in magnesium from first principles

Using Labusch-type solid solution strengthening models parameterized with DFT-computed solute-dislocation interaction energies, we perform a computational search for 63 solutes across the periodic table to find those that lower anisotropy ratios (non-basal to basal CRSS) of magnesium potentially increasing its ductility per the von Mises criterion. For this purpose, we compute changes in strength for solutes as a function of composition and temperature, and compute anisotropy ratios for solutes that include both rare earth and non-rare earth elements. Here we specifically focus on solute-dislocation interaction energies in the following DFT-optimized dislocations as representative of three non-basal plastic deformation modes: $\langle c + a \rangle$ edge, (10$\bar1$2) tension twinning edge, and the (10$\bar{1}$1) compression twinning edge. We find that solute-induced changes in non-basal deformation modes can be approximated using a second-order polynomial in the size misfit of the solutes, which permits rapid screening of solutes. Our approach to identify solutes known to improve strengthening incorporates solute solubility, and suggests other solutes that not have been previously explored for strengthening. The 8 rare-earth solutes that our method suggests as the best, ordered by increasing anisotropy ratios at their optimal concentrations, are: Gd, Tb, Dy, Nd, Ho, Er, Tm, and Yb. The 12 non-rare-earth solutes that our method suggests as the best, ordered by increasing anisotropy ratios, are: Y, Mn, Sc, Pb, Ca, Ag, Bi, Tl, Zn, Li, Ga, and Al. Of these, Gd, Nd, Er, Yb, Y, Mn, Ca, Zn, Li, and Al are used in commercial Mg alloys.

36 MATERIALS SCIENCE↗