Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Density functional”

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 1,369 records · Page 76
All sourcesNASA NTRSDOE OSTI
NASA NTRS2013

Galactic Cosmic Ray Event-Based Risk Model (GERM) Code

This software describes the transport and energy deposition of the passage of galactic cosmic rays in astronaut tissues during space travel, or heavy ion beams in patients in cancer therapy. Space radiation risk is a probability distribution, and time-dependent biological events must be accounted for physical description of space radiation transport in tissues and cells. A stochastic model can calculate the probability density directly without unverified assumptions about shape of probability density function. The prior art of transport codes calculates the average flux and dose of particles behind spacecraft and tissue shielding. Because of the signaling times for activation and relaxation in the cell and tissue, transport code must describe temporal and microspatial density of functions to correlate DNA and oxidative damage with non-targeted effects of signals, bystander, etc. These are absolutely ignored or impossible in the prior art. The GERM code provides scientists data interpretation of experiments; modeling of beam line, shielding of target samples, and sample holders; and estimation of basic physical and biological outputs of their experiments. For mono-energetic ion beams, basic physical and biological properties are calculated for a selected ion type, such as kinetic energy, mass, charge number, absorbed dose, or fluence. Evaluated quantities are linear energy transfer (LET), range (R), absorption and fragmentation cross-sections, and the probability of nuclear interactions after 1 or 5 cm of water equivalent material. In addition, a set of biophysical properties is evaluated, such as the Poisson distribution for a specified cellular area, cell survival curves, and DNA damage yields per cell. Also, the GERM code calculates the radiation transport of the beam line for either a fixed number of user-specified depths or at multiple positions along the Bragg curve of the particle in a selected material. The GERM code makes the numerical estimates of basic physical and biophysical quantities of high-energy protons and heavy ions that have been studied at the NASA Space Radiation Laboratory (NSRL) for the purpose of simulating space radiation biological effects. In the first option, properties of monoenergetic beams are treated. In the second option, the transport of beams in different materials is treated. Similar biophysical properties as in the first option are evaluated for the primary ion and its secondary particles. Additional properties related to the nuclear fragmentation of the beam are evaluated. The GERM code is a computationally efficient Monte-Carlo heavy-ion-beam model. It includes accurate models of LET, range, residual energy, and straggling, and the quantum multiple scattering fragmentation (QMSGRG) nuclear database.

Cucinotta, Francis A.↗
DOE OSTI2024

First-principles study of the structures and redox mechanisms of Ni-rich lithium nickel manganese cobalt oxides

To reduce the cobalt (Co) content in lithium-ion batteries, Ni-rich (high-Ni) lithium nickel manganese cobalt oxides (NMC) are pursued as one of the next-generation cathode materials. However, there is still debate on the crystal and electronic structures of the baseline, LiNiO 2 . Density Functional Theory (DFT) calculations were performed to provide a theoretical understanding of Ni-rich NMC. First, it was found that the commonly used $R\bar{3}m$ structure for LiNiO 2 is metallic, contrary to the experimentally reported mix-conducting behavior. Among the four different space groups, $R\bar{3}m$, C2/m, P2 1 /c, and P2/c, P2/c with charge disproportionation of Ni 2+ and Ni 4+ is the most energetically stable and semiconducting structure of LiNiO 2 . Therefore, the atomic structures of representative Ni-rich NMC were built by partially replacing Ni with Co or Mn in the P2/c LiNiO 2 to form Li x Ni y Mn z Co 1-y-z O 2 . In the fully lithiated (x=1.0) high Ni content NMC (y>0.5), the oxidation state of all Mn ions becomes 4+, while Co ions still maintain 3+, and part of the Ni ions become 3+ to compensate for the charge. Upon delithiation, the local environment shows more variation of the charge states on the transition metal (TM) ions. The average oxidation on each TM follows a sequence of losing electrons that starts from Ni 2+ to Ni 3+ , then oxidizing Ni 3+ and Co 3+ , while Mn 4+ remains electrochemically inactive till x=0. Finally, a general relationship for the oxidation state change in each TM as a function of x is derived and shows agreement with both modeling and experimental data.

25 ENERGY STORAGE↗
DOE OSTI2023

Impact of Pendent Ammonium Groups on Solubility and Cycling Charge Carrier Performance in Nonaqueous Redox Flow Batteries

The synthesis, characterization, electrochemical performance, and theoretical modeling of two base-metal charge carrier complexes incorporating a pendent quaternary ammonium group, [Ni(bppn-Me3)][BF4], 3′, and [Fe(PyTRENMe)][OTf]3, 4’, are described. Both complexes were produced in high yield and fully characterized using NMR, IR, and UV–vis spectroscopies as well as elemental analysis and single-crystal X-ray crystallography. The solubility of 3′ in acetonitrile showed a 283% improvement over its neutral precursor, whereas the solubility of complex 4’ was effectively unchanged. Cyclic voltammetry indicates an ∼0.1 V positive shift for all waves, with some changes in reversibility depending on the wave. Bulk electrochemical cycling demonstrates that both 3′ and 4’ can utilize the second more negative wave to a degree, whereas 4’ ceases to have a reversible positive wave. Flow cell testing of 3′ and 4’ with Fc as the posolyte reveals little improvement to the cycling performance of 3′ compared with its parent complex, whereas 4’ exhibits reductions in capacity decay when cycling either negative wave. Postcycling CVs indicate that crossover is the likely source of capacity loss in complexes 3, 3′, and 4’ because there is little change in the CV trace. Density functional theory calculations indicate that the ammonium group lowers the HOMO energy in 3′ and 4’, which may impart stability to cycling negative waves while making positive waves less accessible. The incorporation of a positively charged species can improve solubility, stored electron density, and capacity decay depending on the complex, features critical to high energy density redox flow battery performance.

25 ENERGY STORAGE↗
DOE OSTI2021

Self-interaction correction in water–ion clusters

We study the importance of self-interaction errors in density functional approximations for various water–ion clusters. We have employed the Fermi-Löwdin orbital self-interaction correction (FLOSIC) method in conjunction with LSDA, PBE, and SCAN to describe binding energies of hydrogen-bonded water–ion clusters, i.e., water–hydronium, water– hydroxide, water–halide, as well as non-hydrogen-bonded water–alkali clusters. In the hydrogen-bonded water–ion clusters, the building blocks are linked by hydrogen atoms, although the links are much stronger and longer-ranged than the normal hydrogen bonds between water molecules, because the monopole on the ion interacts with both permanent and induced dipoles on the water molecules. We find that self-interaction errors overbind the hydrogen-bonded water– ion clusters and that FLOSIC reduces the error and brings the binding energies into closer agreement with higher–level calculations. Here, the non-hydrogen-bonded water–alkali clusters are not significantly affected by self-interaction errors. Self-interaction corrected PBE predicts the lowest mean unsigned error in binding energies (≤ 50 meV/H 2 O) for hydrogen-bonded water–ion clusters. Self-interaction errors are also largely dependent on the cluster size, and FLOSIC does not accurately capture the subtle variation in all clusters, indicating the need for further refinement.

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

The impacts of charge transfer, localization, and metallicity on hydrogen retention and transport capacity

Solid state hydrides such as early transition metal hydrides are of inestimable importance for the future of hydrogen energy and are actively being investigated for energy conversion and storage applications such as fuel cells, solid-state batteries and neutron moderators. The retention and transport behavior of hydrogen in these hydrides has a huge role on the extended performance of components. While early transition-metal-based compounds exhibit many peculiar properties due to their unique correlated electronic signatures arising from -orbital electrons, the fundamental chemistry and transport behavior of hydrogen in such hydrides is not well understood. In the present work, using density functional theory, a highly intricate bonding feature is revealed through the theoretical investigation of the electronic structure of early transition metal hydrides YH 2 and ZrH 2 . In particular, a pronounced charge transfer from the transition element to H, results in localized electron densities at deep energy levels. The interplay between intrinsic charge transfer, charge localization, and metallicity in YH 2 and ZrH 2 leads to strong chemical bonding between metal and hydrogen atoms and large energy barriers for the migration of hydrogen vacancies. Specifically, hydrogen vacancies are found to be stable in the neutral state due to electron screening effects, accompanied by substantially high migration barriers between 0.8–1.2 eV along different crystallographic directions. In contrast, recent literature shows the migration barrier for charged H vacancies in insulating s-block metal hydrides lie between 0.1–0.4 eV, which is suitable for fast conduction applications. This pivotal electron structure difference exploited between early transition metal hydrides and alkali/alkaline earth metal hydrides determines extended hydrogen retention in these early transition metal hydrides. Finally, this work explains fundamental differences between the electronic structure of s-block and d-block metal hydrides, and its impact on the mobility of hydrogen vacancies.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Microscopic calculation of fission product yields with particle-number projection

Fission fragments' charge and mass distribution is an important input to applications ranging from basic science to energy production or nuclear nonproliferation. In simulations of nucleosynthesis or calculations of superheavy elements, these quantities must be computed from models, as they are needed in nuclei where no experimental information is available. Until now, standard techniques to estimate these distributions were not capable of accounting for fine-structure effects, such as the odd-even staggering of the charge distributions. In this work, we combine a fully microscopic collective model of fission dynamics with a recent extension of the particle number projection formalism to provide the highest-fidelity prediction of the primary fission fragment distributions for the neutron-induced fission of 235 U and 239 Pu. Here, we show that particle-number projection is an essential ingredient to reproduce odd-even staggering in the charge yields and benchmark the performance of various empirical probability laws that could simulate its effect. This new approach also enables for the first time the realistic determination of two-dimensional isotopic yields within nuclear density functional theory.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

P–V–T Equation of State of Iridium Up to 80 GPa and 3100 K

In the present study, the high-pressure high-temperature equation of the state of iridium has been determined through a combination of in situ synchrotron X-ray diffraction experiments using laser-heating diamond-anvil cells (up to 48 GPa and 3100 K) and density-functional theory calculations (up to 80 GPa and 3000 K). The melting temperature of iridium at 40 GPa was also determined experimentally as being 4260 (200) K. The results obtained with the two different methods are fully consistent and agree with previous thermal expansion studies performed at ambient pressure. The resulting thermal equation of state can be described using a third-order Birch–Murnaghan formalism with a Berman thermal-expansion model. The present equation of the state of iridium can be used as a reliable primary pressure standard for static experiments up to 80 GPa and 3100 K. A comparison with gold, copper, platinum, niobium, rhenium, tantalum, and osmium is also presented. On top of that, the radial-distribution function of liquid iridium has been determined from experiments and calculations.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Be Surface Structures on W(110) and W(211): A DFT Study

Beryllium and Tungsten are promising candidates for use as plasma facing materials (PFMs) in upcoming fusion reactors. Many complex and competing phenomena however complicate the understanding and development of materials that operate in the harsh conditions of the reactor. In particular, redeposition of Be onto the W divertor must be considered due to the expected erosion of the Be first wall under exposure to the plasma. It is known that a build up of Be on W allows for the formation of BeW alloys which can harm the longevity and performance of the fusion divertor. In an effort to understand the interaction of Be with W surfaces, a study of Be structures on W(110) and W(211) as a function of Be coverage has been carried out using Density Functional Theory. We have found that both W surfaces develop monolayers of Be characterized by a densely packed hexagonal structure. Below this monolayer coverage, Be structures are found to have two motifs of coverage, a densely packed hexagonal pattern punctuated by areas of low-density coverage. The structures found here produce work function values and trends in good agreement with experimental measurements.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Unveiling the ambipolar carrier transport property of SnO 2–X for multiple-functional interlayers in perovskite solar cells

Tin oxides are the most promising electron transport layers in perovskite solar cells. An ambipolar carrier transport property has been recently demonstrated which enables a simple interconnection structure for all-perovskite tandem solar cells. However, the underlying mechanism for its ambipolar behavior is unclear which cannot be explained by the intrinsic defects in SnO 2-x . Here by using density functional theory calculations, we unveil the origin of the ambipolar carrier transport of non-stoichiometry SnO 2-x with a structure of SnO embedded in SnO 2 matrix. The hybridization of O 2p and Sn 5s orbitals of SnO introduces mid-gap states in the bandgap of SnO 2 , enabling hole transport property for SnO 2-x when x is > 0.2. Increasing the percentage of SnO in SnO 2-x significantly enhances the hole transport capability of SnO 2-x due to the enlarged Sn-O-Sn angles that increases orbital overlapping between O and Sn atoms, thus providing strategies for the further tuning of the carrier transport properties of SnO 2-x by compositional and structural designs.

14 SOLAR ENERGY↗
DOE OSTI2026

Local Spin Density Approximation Strongly Improved by a Better-Informed Local Scaling of Its Self-Interaction Correction

The Perdew−Zunger self-interaction correction (PZSIC) makes density functional approximations (DFAs) exact for all one-electron densities. However, it overcorrects in manyelectron regions, introducing errors for the uniform-density limit, where uncorrected DFAs are exact. The locally scaled PZSIC (LSIC), based on the iso-orbital indicator zσ [which distinguishes single-orbital and slowly varying density regions and is used with the local spin density approximation (LSDA)], restores the uniform-density limit and significantly improves results for many properties, including chemical reaction barrier heights, atomization energies, and ionization potentials. Yet, LSIC performs poorly for weakly bonded systems, leaving many unbound, due to limitations of its iso-orbital indicator. To correct this, in this work we propose a new local scaling, LSIC-α, based on the iso-orbital indicator ασ (which additionally identifies regions of overlapping density tails). A two-parameter scaling function of ασ is fitted to a subset of the nonbonded appropriate norms for the SCAN and r2SCAN meta- GGAs, and tested on many properties of main-group atoms, molecules, and molecular complexes. LSIC-α greatly improves the interaction energies of weakly bonded systems in the S22 data set while retaining LSIC’s accuracy for other properties. This work shows that the errors of LSDA (and presumably of higher-level DFAs) can be largely but not entirely repaired by a proper “do no harm” self-interaction correction.

Approximation↗
DOE OSTI2024

Solar enhanced oxygen evolution reaction with transition metal telluride

The photo-enhanced electrocatalytic method of oxygen evolution reaction (OER) shows promise for enhancing the effectiveness of clear energy generation through water splitting by using renewable and sustainable source of energy. However, despite benefits of photoelectrocatalytic (PEC) water splitting, its uses are constrained by its low efficiency as a result of charge carrier recombination, a large overpotential, and sluggish reaction kinetics. Here, we illustrate that Nickel telluride (NiTe) synthesized by hydrothermal methods can function as an extremely effective photo-coupled electrochemical oxygen evolution reaction (POER) catalyst. In this study, NiTe was synthesized by hydrothermal method at 145°C within just an hour of reaction time. In dark conditions, the NiTe deposited on carbon cloth substrate shows a small oxygen evolution reaction overpotential (261 mV) at a current density of 10 mA cm –2 , a reduced Tafel slope (65.4 mV dec -1 ), and negligible activity decay after 12 h of chronoamperometry. By virtue of its enhanced photo response, excellent light harvesting ability, and increased interfacial kinetics of charge separation, the NiTe electrode under simulated solar illumination displays exceptional photoelectrochemical performance exhibiting overpotential of 165 mV at current density of 10 mA cm -2 , which is about 96 mV less than on dark conditions. In addition, Density Functional Theory investigations have been carried out on the NiTe surface, the results of which demonstrated a greater adsorption energy for intermediate -OH on the catalyst site. Since the -OH adsorption on the catalyst site correlates to catalyst activation, it indicates the facile electrocatalytic activity of NiTe owing to favorable catalyst activation. DFT calculations also revealed the facile charge density redistribution following intermediate -OH adsorption on the NiTe surface. This work demonstrates that arrays of NiTe elongated nanostructure are a promising option for both electrochemical and photoelectrocatalytic water oxidation and offers broad suggestions for developing effective PEC devices.

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

Thermodynamics of the Metal Carbonates and Bicarbonates of Mn, Co, Ni, Cu, and Zn Relevant to Mineral Energetics

The gas phase heats of formation of ground-state MCO 3 , M(HCO 3 ) 2 , and M(HCO 3 )(OH), where M = Mn, Co, Ni, Cu, and Zn, have been predicted using the correlated molecular orbital theory at the CCSD(T) level extrapolated to the complete basis set limit using the Feller–Peterson–Dixon (FPD) approach. Cohesive energies of the carbonates were predicted based on the calculated gas phase and experimental solid heats of formation. Coulombic dissociation energies (CDEs) between metal cations and anions show a near-linear correlation with Shannon metal cation atomic radii, yet no correlation is found with the hardness of these cations. The total reaction dissociation energies (TRDEs) of transition metals are higher than their CDEs for the di-bicarbonates, in contrast to those for Mg and Ca based on our prior work. In addition to differences in the energies needed to prepare the transition metal dications, electron donation from the ligands to the 3d orbitals of open-shell transition metal dications from lone pairs of adjacent O atoms also plays a role. No electron donation from the ligands to the fully occupied 3d orbitals of Zn and Cd was found. Decomposition energies for generating MO, CO 2 , and/or H 2 O were calculated. Gas phase metal exchange energies only partially correlate with the electrochemical series for M(s) → M 2+( aq). The FPD heats of formation were used to benchmark a range of density functional theory exchange–correlation functionals, including those commonly used in solid-state mineral calculations. None of the functionals provided chemical accuracy agreement (±1 kcal/mol) with the FPD results. So, the best agreement with the FPD results is predicted for the $\tau$-HCTH functional with an average unsigned error of 8.3 kcal/mol.

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

CsV 3 Sb 5 : A $\mathbb{Z}_2$ Topological Kagome Metal with a Superconducting Ground State

Recently discovered alongside its sister compounds KV 3 Sb 5 and RbV 3 Sb 5 , CsV 3 Sb 5 crystallizes with an ideal kagome network of vanadium and antimonene layers separated by alkali metal ions. This work presents the electronic properties of CsV 3 Sb 5 , demonstrating bulk superconductivity in single crystals with a T c = 2.5 K. The normal state electronic structure is studied via angle-resolved photoemission spectroscopy and density-functional theory, which categorize CsV 3 Sb 5 as a Z 2 topological metal. Multiple protected Dirac crossings are predicted in close proximity to the Fermi level (E F ), and signatures of normal state correlation effects are also suggested by a high-temperature charge density wavelike instability. Finally, the implications for the formation of unconventional superconductivity in this material are discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Ordering in GaAs co-doped with Bi and N

Introducing only a few atomic percent of Bi or N in GaAs has a large effect on the band gap of the material. Specifically Bi doped GaAs shows potential for local band gap engineering in optoelectronic applications. The incorporation of Bi and N into GaAs is difficult due to strain effects. In this work we study the ordering of these dopants at the atomic scale in order to get a better understanding of the behavior of these dopants in the host lattice. Cross-sectional scanning tunneling microscopy is used to find the exact position of Bi and N dopants in the GaAs matrix, allowing us to study both their nearest neighbor pair occurrences and pair correlation functions. An attractive interaction between Bi dopants at short ranges (1-2 nm) is found and a similar effect is observed between N dopants. Here, we find a repulsive interaction with a similar length scale between Bi and N dopants. A similar repulsion is found in the Bi-N nearest neighbor pairs. Density functional theory is used to calculate the different nearest neighbor pair energies and test these results to the experimental pair occurrences. It is concluded from the experimental and theoretical results that the growth conditions and N inclusion greatly affects the Bi distribution in GaAs.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Be₂C as a neutron moderator for molten salt reactors: A first-principles study of structural, electronic, and defect properties

Beryllium carbide (Be 2 C), valued for its high neutron moderation efficiency and low absorption cross section, is a promising high-temperature neutron moderator for molten salt reactors. Its practical adoption, however, demands significant technological maturation, requiring comprehensive theoretical and experimental studies of its response to different conditions, including high temperature and irradiation. Here, we report initial results on the fundamental properties and radiation-induced defects of Be 2 C, focusing on antisites, vacancies, interstitial atoms, and Frenkel pairs in the Be and C sublattices. Using density functional theory (DFT) and ab initio molecular dynamics (AIMD), we calculate the defects formation and binding energies, evaluating their dependence on the supercell size, charge states, and chemical environment. In general, carbon defects exhibit higher formation energies, greater sensitivity to cell size, and stronger impacts on the density of states compared to beryllium defects, with charged state the effects being more pronounced. Static DFT reveals multiple metastable interstitial configurations, while AIMD identifies ground states as C-C <100> dumbbells and octahedral Be interstitials. In conclusion, the diversity of metastable configurations and defect states complicates the diffusion mechanisms, requiring further molecular dynamics analysis to elucidate the mechanisms and rates of radiation-induced atomic transport, as well as the structural stability of Be2C.

Ab initio modeling↗
DOE OSTI2023

Fluorinated High-Voltage Electrolytes To Stabilize Nickel-Rich Lithium Batteries

As state-of-the-art (SOA) lithium-ion (Li-ion) batteries approach their specific energy limit (~250 Wh kg -1 ), layer-structured, nickel-rich (Ni-rich) lithium transition metal oxide-based cathode materials, e. g., LiNi 0.8 Mn 0.1 Co 0.1 O 2 (NMC811), have attracted great interest owing to their practical high specific capacities (>200 mAhg -1 ). Coupled with their high average discharge voltages (~4 V vs Li/Li + ), Ni-rich cathode-based lithium batteries possess a great potential to achieve much higher specific energies (>350 Wh kg -1 at the cell level) than the SOA Li-ion counterparts. In addition, Ni-rich oxides are low-cost battery cathode materials due to their low cobalt contents. However, Ni-rich cathode-based lithium batteries suffer quick capacity degradations upon cycling, particularly at high upper cutoff voltages (e.g., ≥= 4.5 V vs Li/Li + ), due to crystal structure changes of the active cathode materials and parasitic side reactions at the electrolyte/electrode interfaces. In this study, a fluorinated-solvent-based, high-voltage stable electrolyte (HVE), i.e., 1 M Li bis(trifluoromethanesulfonyl)imide (LiTFSI) in fluoroethylene carbonate (FEC), bis(2,2,2-trifluoroethyl) carbonate (FDEC), and methyl (2,2,2-trifluoroethyl) carbonate (FEMC) with Li difluoro(oxalate)borate (LiDFOB) additive, was formulated and evaluated in Li/NMC811 battery cells. Further, to the best of our knowledge, this class of electrolyte has not been investigated for Ni-rich cathode-based lithium batteries. Li/NMC811 cells with HVE exhibited a superior long-term cycle performance stability, maintaining ~80% capacity after ~500 cycles at a high cutoff voltage of 4.5 V (vs Li/Li + ) than a baseline carbonate-solvent-based electrolyte (BE). The superior cycle stability of the Li/NMC811 cells is attributed to the inherently high-voltage stability of HVE, supported by the physical and electrochemical analyses. This conclusion is supported by our density functional theory (DFT) modeling where HVE shows a less tendency of deprotonation/oxidation than BE, leading to the observed cycle stability. The findings in this study are important to help tackle the technical challenges facing Ni-rich cathode-based lithium batteries to realize their high energy density potentials with a long cycle life.

25 ENERGY STORAGE↗
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 OSTI2021

Navigating Transition-Metal Chemical Space: Artificial Intelligence for First-Principles Design

Conspectus The variability of chemical bonding in open-shell transition-metal complexes not only motivates their study as functional materials and catalysts but also challenges conventional computational modeling tools. Here, tailoring ligand chemistry can alter preferred spin or oxidation states as well as electronic structure properties and reactivity, creating vast regions of chemical space to explore when designing new materials atom by atom. Although first-principles density functional theory (DFT) remains the workhorse of computational chemistry in mechanism deduction and property prediction, it is of limited use here. DFT is both far too computationally costly for widespread exploration of transition-metal chemical space and also prone to inaccuracies that limit its predictive performance for localized d electrons in transition-metal complexes. These challenges starkly contrast with the well-trodden regions of small-organic-molecule chemical space, where the analytical forms of molecular mechanics force fields and semiempirical theories have for decades accelerated the discovery of new molecules, accurate DFT functional performance has been demonstrated, and gold-standard methods from correlated wavefunction theory can predict experimental results to chemical accuracy. The combined promise of transition-metal chemical space exploration and lack of established tools has mandated a distinct approach. In this Account, we outline the path we charted in exploration of transition-metal chemical space starting from the first machine learning (ML) models (i.e., artificial neural network and kernel ridge regression) and representations for the prediction of open-shell transition-metal complex properties. The distinct importance of the immediate coordination environment of the metal center as well as the lack of low-level methods to accurately predict structural properties in this coordination environment first motivated and then benefited from these ML models and representations. Once developed, the recipe for prediction of geometric, spin state, and redox potential properties was straightforwardly extended to a diverse range of other properties, including in catalysis, computational “feasibility”, and the gas separation properties of periodic metal–organic frameworks. Interpretation of selected features most important for model prediction revealed new ways to encapsulate design rules and confirmed that models were robustly mapping essential structure–property relationships. Encountering the special challenge of ensuring that good model performance could generalize to new discovery targets motivated investigation of how to best carry out model uncertainty quantification. Distance-based approaches, whether in model latent space or in carefully engineered feature space, provided intuitive measures of the domain of applicability. With all of these pieces together, ML can be harnessed as an engine to tackle the large-scale exploration of transition-metal chemical space needed to satisfy multiple objectives using efficient global optimization methods. In practical terms, bringing these artificial intelligence tools to bear on the problems of transition-metal chemical space exploration has resulted in ML-model assessments of large, multimillion compound spaces in minutes and validated new design leads in weeks instead of decades.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗