-
Accelerating the prediction of inorganic surfaces with machine learning interatomic potentials
Authors:
Kyle Noordhoek,
Christopher J. Bartel
Abstract:
The surface properties of solid-state materials often dictate their functionality, especially for applications where nanoscale effects become important. The relevant surface(s) and their properties are determined, in large part, by the materials synthesis or operating conditions. These conditions dictate thermodynamic driving forces and kinetic rates responsible for yielding the observed surface s…
▽ More
The surface properties of solid-state materials often dictate their functionality, especially for applications where nanoscale effects become important. The relevant surface(s) and their properties are determined, in large part, by the materials synthesis or operating conditions. These conditions dictate thermodynamic driving forces and kinetic rates responsible for yielding the observed surface structure and morphology. Computational surface science methods have long been applied to connect thermochemical conditions to surface phase stability, particularly in the heterogeneous catalysis and thin film growth communities. This review provides a brief introduction to first-principles approaches to compute surface phase diagrams before introducing emerging data-driven approaches. The remainder of the review focuses on the application of machine learning, predominantly in the form of learned interatomic potentials, to study complex surfaces. As machine learning algorithms and large datasets on which to train them become more commonplace in materials science, computational methods are poised to become even more predictive and powerful for modeling the complexities of inorganic surfaces at the nanoscale.
△ Less
Submitted 18 December, 2023;
originally announced December 2023.
-
Selective formation of metastable polymorphs in solid-state synthesis
Authors:
Yan Zeng,
Nathan J. Szymanski,
Tanjin He,
KyuJung Jun,
Leighanne C. Gallington,
Haoyan Huo,
Christopher J. Bartel,
Bin Ouyang,
Gerbrand Ceder
Abstract:
Metastable polymorphs often result from the interplay between thermodynamics and kinetics. Despite advances in predictive synthesis for solution-based techniques, there remains a lack of methods to design solid-state reactions targeting metastable materials. Here, we introduce a theoretical framework to predict and control polymorph selectivity in solid-state reactions. This framework presents rea…
▽ More
Metastable polymorphs often result from the interplay between thermodynamics and kinetics. Despite advances in predictive synthesis for solution-based techniques, there remains a lack of methods to design solid-state reactions targeting metastable materials. Here, we introduce a theoretical framework to predict and control polymorph selectivity in solid-state reactions. This framework presents reaction energy as a rarely used handle for polymorph selection, which influences the role of surface energy in promoting the nucleation of metastable phases. Through in situ characterization and density functional theory calculations on two distinct synthesis pathways targeting LiTiOPO4, we demonstrate how precursor selection and its effect on reaction energy can effectively be used to control which polymorph is obtained from solid-state synthesis. A general approach is outlined to quantify the conditions under which metastable polymorphs are experimentally accessible. With comparison to historical data, this approach suggests that using appropriate precursors could enable the synthesis of many novel materials through selective polymorph nucleation.
△ Less
Submitted 11 September, 2023;
originally announced September 2023.
-
Assessing Thermodynamic Selectivity of Solid-State Reactions for the Predictive Synthesis of Inorganic Materials
Authors:
Matthew J. McDermott,
Brennan C. McBride,
Corlyn Regier,
Gia Thinh Tran,
Yu Chen,
Adam A. Corrao,
Max C. Gallant,
Gabrielle E. Kamm,
Christopher J. Bartel,
Karena W. Chapman,
Peter G. Khalifah,
Gerbrand Ceder,
James R. Neilson,
Kristin A. Persson
Abstract:
Synthesis is a major challenge in the discovery of new inorganic materials. Currently, there is limited theoretical guidance for identifying optimal solid-state synthesis procedures. We introduce two selectivity metrics, primary and secondary competition, to assess the favorability of target/impurity phase formation in solid-state reactions. We used these metrics to analyze 3,520 solid-state react…
▽ More
Synthesis is a major challenge in the discovery of new inorganic materials. Currently, there is limited theoretical guidance for identifying optimal solid-state synthesis procedures. We introduce two selectivity metrics, primary and secondary competition, to assess the favorability of target/impurity phase formation in solid-state reactions. We used these metrics to analyze 3,520 solid-state reactions in the literature, ranking existing approaches to popular target materials. Additionally, we implemented these metrics in a data-driven synthesis planning workflow and demonstrated its application in the synthesis of barium titanate (BaTiO$_3$). Using an 18-element chemical reaction network with first-principles thermodynamic data from the Materials Project, we identified 82,985 possible BaTiO$_3$ synthesis reactions and selected nine for experimental testing. Characterization of reaction pathways via synchrotron powder X-ray diffraction reveals that our selectivity metrics correlate with observed target/impurity formation. We discovered two efficient reactions using unconventional precursors (BaS/BaCl$_2$ and Na$_2$TiO$_3$) that produce BaTiO$_3$ faster and with fewer impurities than conventional methods, highlighting the importance of considering complex chemistries with additional elements during precursor selection. Our framework provides a foundation for predictive inorganic synthesis, facilitating the optimization of existing recipes and the discovery of new materials, including those not easily attainable with conventional precursors.
△ Less
Submitted 27 September, 2023; v1 submitted 22 August, 2023;
originally announced August 2023.
-
Interpretable machine learning to understand the performance of semi local density functionals for materials thermochemistry
Authors:
Santosh Adhikari,
Christopher J. Bartel,
Christopher Sutton
Abstract:
This study investigates the use of machine learning (ML) to correct the enthalpy of formation (Hf) from two separate DFT functionals, PBE and SCAN, to the experimental Hf across 1011 solid-state compounds. The ML model uses a set of 25 properties that characterize the electronic structure as calculated using PBE and SCAN. The ML model significantly decreases the error in PBE-calculated Hf values f…
▽ More
This study investigates the use of machine learning (ML) to correct the enthalpy of formation (Hf) from two separate DFT functionals, PBE and SCAN, to the experimental Hf across 1011 solid-state compounds. The ML model uses a set of 25 properties that characterize the electronic structure as calculated using PBE and SCAN. The ML model significantly decreases the error in PBE-calculated Hf values from an mean absolute error (MAE) of 195 meV/atom to an MAE = 80 meV/atom when compared to the experiment. For PBE, the PDP+GAM analysis shows compounds with a high ionicity (I), i.e., I>0.22, have errors in Hf that are twice as large as compounds having I < 0.22 (246 meV/atom compared to 113 meV/atom). Conversely, no analogous trend is observed for SCAN-calculated Hfs, which explains why the ML model for PBE can more easily correct the systematic error in calculated Hfs for PBE but not for SCAN. Although the literature suggests PBE is reliable for intermetallics but less so for oxides and halides, our analysis reveals intermetallics pose a challenge for PBE only when the charge transfer is significant (I >0.22). Meanwhile, oxides and halides may be described accurately by PBE for systems in which charge transfer is relatively low (I < 0.22).
△ Less
Submitted 14 July, 2023;
originally announced July 2023.
-
Autonomous decision making for solid-state synthesis of inorganic materials
Authors:
Nathan J. Szymanski,
Pragnay Nevatia,
Christopher J. Bartel,
Yan Zeng,
Gerbrand Ceder
Abstract:
To aid in the automation of inorganic materials synthesis, we introduce an algorithm (ARROWS3) that guides the selection of precursors used in solid-state reactions. Given a target phase, ARROWS3 iteratively proposes experiments and learns from their outcomes to identify an optimal set of precursors that leads to maximal yield of that target. Initial experiments are selected based on thermochemica…
▽ More
To aid in the automation of inorganic materials synthesis, we introduce an algorithm (ARROWS3) that guides the selection of precursors used in solid-state reactions. Given a target phase, ARROWS3 iteratively proposes experiments and learns from their outcomes to identify an optimal set of precursors that leads to maximal yield of that target. Initial experiments are selected based on thermochemical data collected from first principles calculations, which enable the identification of precursors exhibiting large thermodynamic force to form the desired target. Should the initial experiments fail, their associated reaction paths are determined by sampling a range of synthesis temperatures and identifying their products. ARROWS3 then uses this information to pinpoint which intermediate reactions consume most of the available free energy associated with the starting materials. In subsequent experimental iterations, precursors are selected to avoid such unfavorable reactions and therefore maintain a strong driving force to form the target. We validate this approach on three experimental datasets containing results from more than 200 distinct synthesis procedures. When compared to several black-box optimization algorithms, ARROWS3 identifies the most effective set of precursors for each target while requiring substantially fewer experimental iterations. These findings highlight the importance of using domain knowledge in the design of optimization algorithms for materials synthesis, which are critical for the development of fully autonomous research platforms.
△ Less
Submitted 24 May, 2023; v1 submitted 18 April, 2023;
originally announced April 2023.
-
CHGNet: Pretrained universal neural network potential for charge-informed atomistic modeling
Authors:
Bowen Deng,
Peichen Zhong,
KyuJung Jun,
Janosh Riebesell,
Kevin Han,
Christopher J. Bartel,
Gerbrand Ceder
Abstract:
The simulation of large-scale systems with complex electron interactions remains one of the greatest challenges for the atomistic modeling of materials. Although classical force fields often fail to describe the coupling between electronic states and ionic rearrangements, the more accurate \textit{ab-initio} molecular dynamics suffers from computational complexity that prevents long-time and large…
▽ More
The simulation of large-scale systems with complex electron interactions remains one of the greatest challenges for the atomistic modeling of materials. Although classical force fields often fail to describe the coupling between electronic states and ionic rearrangements, the more accurate \textit{ab-initio} molecular dynamics suffers from computational complexity that prevents long-time and large-scale simulations, which are essential to study many technologically relevant phenomena, such as reactions, ion migrations, phase transformations, and degradation.
In this work, we present the Crystal Hamiltonian Graph neural Network (CHGNet) as a novel machine-learning interatomic potential (MLIP), using a graph-neural-network-based force field to model a universal potential energy surface. CHGNet is pretrained on the energies, forces, stresses, and magnetic moments from the Materials Project Trajectory Dataset, which consists of over 10 years of density functional theory static and relaxation trajectories of $\sim 1.5$ million inorganic structures. The explicit inclusion of magnetic moments enables CHGNet to learn and accurately represent the orbital occupancy of electrons, enhancing its capability to describe both atomic and electronic degrees of freedom. We demonstrate several applications of CHGNet in solid-state materials, including charge-informed molecular dynamics in Li$_x$MnO$_2$, the finite temperature phase diagram for Li$_x$FePO$_4$ and Li diffusion in garnet conductors. We critically analyze the significance of including charge information for capturing appropriate chemistry, and we provide new insights into ionic systems with additional electronic degrees of freedom that can not be observed by previous MLIPs.
△ Less
Submitted 20 June, 2023; v1 submitted 27 February, 2023;
originally announced February 2023.
-
Precursor recommendation for inorganic synthesis by machine learning materials similarity from scientific literature
Authors:
Tanjin He,
Haoyan Huo,
Christopher J. Bartel,
Zheren Wang,
Kevin Cruse,
Gerbrand Ceder
Abstract:
Synthesis prediction is a key accelerator for the rapid design of advanced materials. However, determining synthesis variables such as the choice of precursor materials is challenging for inorganic materials because the sequence of reactions during heating is not well understood. In this work, we use a knowledge base of 29,900 solid-state synthesis recipes, text-mined from the scientific literatur…
▽ More
Synthesis prediction is a key accelerator for the rapid design of advanced materials. However, determining synthesis variables such as the choice of precursor materials is challenging for inorganic materials because the sequence of reactions during heating is not well understood. In this work, we use a knowledge base of 29,900 solid-state synthesis recipes, text-mined from the scientific literature, to automatically learn which precursors to recommend for the synthesis of a novel target material. The data-driven approach learns chemical similarity of materials and refers the synthesis of a new target to precedent synthesis procedures of similar materials, mimicking human synthesis design. When proposing five precursor sets for each of 2,654 unseen test target materials, the recommendation strategy achieves a success rate of at least 82%. Our approach captures decades of heuristic synthesis data in a mathematical form, making it accessible for use in recommendation engines and autonomous laboratories.
△ Less
Submitted 19 May, 2023; v1 submitted 4 February, 2023;
originally announced February 2023.
-
Testing the r$^2$SCAN density functional for the thermodynamic stability of solids with and without a van der Waals correction
Authors:
Manish Kothakonda,
Aaron D. Kaplan,
Eric B. Isaacs,
Christopher J. Bartel,
James W. Furness,
Jinliang Ning,
Chris Wolverton,
John P. Perdew,
Jianwei Sun
Abstract:
A central aim of materials discovery is an accurate and numerically reliable description of thermodynamic properties, such as the enthalpies of formation and decomposition. The r$^2$SCAN revision of the strongly constrained and appropriately normed (SCAN) meta-generalized gradient approximation (meta-GGA) balances numerical stability with high general accuracy. To assess the r$^2$SCAN description…
▽ More
A central aim of materials discovery is an accurate and numerically reliable description of thermodynamic properties, such as the enthalpies of formation and decomposition. The r$^2$SCAN revision of the strongly constrained and appropriately normed (SCAN) meta-generalized gradient approximation (meta-GGA) balances numerical stability with high general accuracy. To assess the r$^2$SCAN description of solid-state thermodynamics, we evaluate the formation and decomposition enthalpies, equilibrium volumes, and fundamental bandgaps of more than 1,000 solids using r$^2$SCAN, SCAN, and PBE, as well as two dispersion-corrected variants, SCAN+rVV10 and r$^2$SCAN+rVV10. We show that r$^2$SCAN achieves accuracy comparable to SCAN and often improves upon SCAN's already excellent accuracy. Whereas SCAN+rVV10 is often observed to worsen the formation enthalpies of SCAN, and makes no substantial correction to SCAN's cell volume predictions, r$^2$SCAN+rVV10 predicts marginally less-accurate formation enthalpies than r$^2$SCAN, and slightly more-accurate cell volumes than r$^2$SCAN. The average absolute errors in predicted formation enthalpies are found to decrease by a factor of 1.5 to 2.5 from the GGA level to the meta-GGA level. Smaller decreases in error are observed for decomposition enthalpies. For formation enthalpies r$^2$SCAN improves over SCAN for intermetallic systems. For a few classes of systems -- transition metals, intermetallics, weakly-bound solids, and enthalpies of decomposition into compounds -- GGAs are comparable to meta-GGAs. In total, r$^2$SCAN and r$^2$SCAN+rVV10 can be recommended as stable, general-purpose meta-GGAs for materials discovery.
△ Less
Submitted 4 August, 2022;
originally announced August 2022.
-
Machine-learning rationalization and prediction of solid-state synthesis conditions
Authors:
Haoyan Huo,
Christopher J. Bartel,
Tanjin He,
Amalie Trewartha,
Alexander Dunn,
Bin Ouyang,
Anubhav Jain,
Gerbrand Ceder
Abstract:
There currently exist no quantitative methods to determine the appropriate conditions for solid-state synthesis. This not only hinders the experimental realization of novel materials but also complicates the interpretation and understanding of solid-state reaction mechanisms. Here, we demonstrate a machine-learning approach that predicts synthesis conditions using large solid-state synthesis datas…
▽ More
There currently exist no quantitative methods to determine the appropriate conditions for solid-state synthesis. This not only hinders the experimental realization of novel materials but also complicates the interpretation and understanding of solid-state reaction mechanisms. Here, we demonstrate a machine-learning approach that predicts synthesis conditions using large solid-state synthesis datasets text-mined from scientific journal articles. Using feature importance ranking analysis, we discovered that optimal heating temperatures have strong correlations with the stability of precursor materials quantified using melting points and formation energies ($ΔG_f$, $ΔH_f$). In contrast, features derived from the thermodynamics of synthesis-related reactions did not directly correlate to the chosen heating temperatures. This correlation between optimal solid-state heating temperature and precursor stability extends Tamman's rule from intermetallics to oxide systems, suggesting the importance of reaction kinetics in determining synthesis conditions. Heating times are shown to be strongly correlated with the chosen experimental procedures and instrument setups, which may be indicative of human bias in the dataset. Using these predictive features, we constructed machine-learning models with good performance and general applicability to predict the conditions required to synthesize diverse chemical systems. Codes and data used in this work can be found at: https://github.com/CederGroupHub/s4.
△ Less
Submitted 5 August, 2022; v1 submitted 17 April, 2022;
originally announced April 2022.
-
A probabilistic deep learning approach to automate the interpretation of multi-phase diffraction spectra
Authors:
Nathan J. Szymanski,
Christopher J. Bartel,
Yan Zeng,
Qingsong Tu,
Gerbrand Ceder
Abstract:
Autonomous synthesis and characterization of inorganic materials requires the automatic and accurate analysis of X-ray diffraction spectra. For this task, we designed a probabilistic deep learning algorithm to identify complex multi-phase mixtures. At the core of this algorithm lies an ensemble convolutional neural network trained on simulated diffraction spectra, which are systematically augmente…
▽ More
Autonomous synthesis and characterization of inorganic materials requires the automatic and accurate analysis of X-ray diffraction spectra. For this task, we designed a probabilistic deep learning algorithm to identify complex multi-phase mixtures. At the core of this algorithm lies an ensemble convolutional neural network trained on simulated diffraction spectra, which are systematically augmented with physics-informed perturbations to account for artifacts that can arise during experimental sample preparation and synthesis. Larger perturbations associated with off-stoichiometry are also captured by supplementing the training set with hypothetical solid solutions. Spectra containing mixtures of materials are analyzed with a newly developed branching algorithm that utilizes the probabilistic nature of the neural network to explore suspected mixtures and identify the set of phases that maximize confidence in the prediction. Our model is benchmarked on simulated and experimentally measured diffraction spectra, showing exceptional performance with accuracies exceeding those given by previously reported methods based on profile matching and deep learning. We envision that the algorithm presented here may be integrated in experimental workflows to facilitate the high-throughput and autonomous discovery of inorganic materials.
△ Less
Submitted 30 March, 2021;
originally announced March 2021.
-
Synthetic accessibility and stability rules of NASICONs
Authors:
Bin Ouyang,
Jingyang Wang,
Tanjin He,
Christopher J. Bartel,
Haoyan Huo,
Yan Wang,
Valentina Lacivita,
Haegyeom Kim,
Gerbrand Ceder
Abstract:
In this paper we develop the stability rules for NASICON structured materials, as an example of compounds with complex bond topology and composition. By applying machine learning to the ab-initio computed phase stability of 3881 potential NASICONs we can extract a simple two-dimensional descriptor that is extremely good at separating stable from unstable NASICONS. This machine-learned "tolerance f…
▽ More
In this paper we develop the stability rules for NASICON structured materials, as an example of compounds with complex bond topology and composition. By applying machine learning to the ab-initio computed phase stability of 3881 potential NASICONs we can extract a simple two-dimensional descriptor that is extremely good at separating stable from unstable NASICONS. This machine-learned "tolerance factor" contains information on the Na content, the radii and electronegativities of the elements, and the Madelung energy. We test the predictive capability of this approach by selecting six predicted NASICON compositions. Five out of the six resulted in a phase pure NASICON while the sixth composition led to a NASICON that coexisted with other phases, validating the efficacy of this approach. This work not only provide tools to understand synthetic accessibility of NASICON-type materials, but also demonstrate an efficient paradigm for discovering new materials with complicate composition and atomic structure.
△ Less
Submitted 9 September, 2021; v1 submitted 6 February, 2021;
originally announced February 2021.
-
Alloying behavior of wide band gap alkaline-earth chalcogenides
Authors:
Samantha L. Millican,
Jacob M. Clary,
Christopher J. Bartel,
Nicholas R. Singstock,
Aaron M. Holder,
Charles B. Musgrave
Abstract:
Alloying is a powerful tool for tuning materials that facilitates the targeted design of desirable properties for a variety of applications. In this work, we provide a comprehensive investigation of the synthetic accessibility and electronic properties of nine alkaline-earth chalcogenide anion alloys (CaS$_{1-x}$O$_x$, CaS$_{1-x}$Se$_x$, CaS$_{1-x}$Te$_x$, SrS$_{1-x}$O$_x$, SrS$_{1-x}$Se$_x$, SrS…
▽ More
Alloying is a powerful tool for tuning materials that facilitates the targeted design of desirable properties for a variety of applications. In this work, we provide a comprehensive investigation of the synthetic accessibility and electronic properties of nine alkaline-earth chalcogenide anion alloys (CaS$_{1-x}$O$_x$, CaS$_{1-x}$Se$_x$, CaS$_{1-x}$Te$_x$, SrS$_{1-x}$O$_x$, SrS$_{1-x}$Se$_x$, SrS$_{1-x}$Te$_x$, MgS$_{1- x}$O$_x$, MgS$_{1-x}$Se$_x$, and MgS$_{1-x}$Te$_x$). We show that isostructural alloying within the rock salt structure is favored for all systems except MgS$_{1-x}$Te$_x$, which is predicted to be a heterostructural alloy between the rock salt and wurtzite structures. Alloys of S and Se are shown to be readily accessible for all cations with low miscibility critical temperatures, enabling continuous tuning of electronic properties across this composition space. Alloys of S and Te have higher critical temperatures but may be accessible through non-equilibrium synthesis strategies and are predicted here to have desirable electronic properties for optoelectronics with wide band gaps and lower effective masses than alloys of S and Se. Anion alloying in MgS$_{1-x}$Te$_x$ stabilizes the wurtzite structure across a significant fraction of composition space, which may make it of particular interest as a transparent conducting material due to its lower effective masses and a higher band gap than the rock salt structure. Zero-point corrected random phase approximation (RPA) energies were computed to resolve the small polymorph energy differences of the Mg compounds and are shown to be critical for accurately describing the thermodynamic properties of the corresponding alloys.
△ Less
Submitted 12 November, 2020;
originally announced November 2020.
-
Observing and modeling the sequential pairwise reactions that drive solid-state ceramic synthesis
Authors:
Akira Miura,
Christopher J. Bartel,
Yusuke Goto,
Yoshikazu Mizuguchi,
Chikako Moriyoshi,
Yoshihiro Kuroiwa,
Yongming Wang,
Toshie Yaguchi,
Manabu Shirai,
Masanori Nagao,
Nataly Carolina Rosero-Navarro,
Kiyoharu Tadanaga,
Gerbrand Ceder,
Wenhao Sun
Abstract:
Solid-state synthesis from powder precursors is the primary processing route to advanced multicomponent ceramic materials. Designing ceramic synthesis routes is usually a laborious, trial-and-error process, as heterogeneous mixtures of powder precursors often evolve through a complicated series of reaction intermediates. Here, we show that phase evolution from multiple precursors can be modeled as…
▽ More
Solid-state synthesis from powder precursors is the primary processing route to advanced multicomponent ceramic materials. Designing ceramic synthesis routes is usually a laborious, trial-and-error process, as heterogeneous mixtures of powder precursors often evolve through a complicated series of reaction intermediates. Here, we show that phase evolution from multiple precursors can be modeled as a sequence of pairwise interfacial reactions, with thermodynamic driving forces that can be efficiently calculated using ab initio methods. Using the synthesis of the classic high-temperature superconductor YBa$_2$Cu$_3$O$_{6+x}$ (YBCO) as a representative system, we rationalize how replacing the common BaCO$_3$ precursor with BaO$_2$ redirects phase evolution through a kinetically-facile pathway. Our model is validated from in situ X-ray diffraction and in situ microscopy observations, which show rapid YBCO formation from BaO$_2$ in only 30 minutes. By combining thermodynamic modeling with in situ characterization, we introduce a new computable framework to interpret and ultimately design synthesis pathways to complex ceramic materials.
△ Less
Submitted 13 January, 2021; v1 submitted 22 September, 2020;
originally announced September 2020.
-
A critical examination of compound stability predictions from machine-learned formation energies
Authors:
Christopher J. Bartel,
Amalie Trewartha,
Qi Wang,
Alexander Dunn,
Anubhav Jain,
Gerbrand Ceder
Abstract:
Machine learning has emerged as a novel tool for the efficient prediction of materials properties, and claims have been made that machine-learned models for the formation energy of compounds can approach the accuracy of Density Functional Theory (DFT). The models tested in this work include five recently published compositional models, a baseline model using stoichiometry alone, and a structural m…
▽ More
Machine learning has emerged as a novel tool for the efficient prediction of materials properties, and claims have been made that machine-learned models for the formation energy of compounds can approach the accuracy of Density Functional Theory (DFT). The models tested in this work include five recently published compositional models, a baseline model using stoichiometry alone, and a structural model. By testing seven machine learning models for formation energy on stability predictions using the Materials Project database of DFT calculations for 85,014 unique chemical compositions, we show that while formation energies can indeed be predicted well, all compositional models perform poorly on predicting the stability of compounds, making them considerably less useful than DFT for the discovery and design of new solids. Most critically, in sparse chemical spaces where few stoichiometries have stable compounds, only the structural model is capable of efficiently detecting which materials are stable. The non-incremental improvement of structural models compared with compositional models is noteworthy and encourages the use of structural models for materials discovery, with the constraint that for any new composition, the ground-state structure is not known a priori. This work demonstrates that accurate predictions of formation energy do not imply accurate predictions of stability, emphasizing the importance of assessing model performance on stability predictions, for which we provide a set of publicly available tests.
△ Less
Submitted 11 May, 2020; v1 submitted 28 January, 2020;
originally announced January 2020.
-
The role of decomposition reactions in assessing first-principles predictions of solid stability
Authors:
Christopher J. Bartel,
Alan W. Weimer,
Stephan Lany,
Charles B. Musgrave,
Aaron M. Holder
Abstract:
The performance of density functional theory (DFT) approximations for predicting materials thermodynamics is typically assessed by comparing calculated and experimentally determined enthalpies of formation from elemental phases, ΔHf. However, a compound competes thermodynamically with both other compounds and their constituent elemental forms, and thus, the enthalpies of the decomposition reaction…
▽ More
The performance of density functional theory (DFT) approximations for predicting materials thermodynamics is typically assessed by comparing calculated and experimentally determined enthalpies of formation from elemental phases, ΔHf. However, a compound competes thermodynamically with both other compounds and their constituent elemental forms, and thus, the enthalpies of the decomposition reactions to these competing phases, ΔHd, determines thermodynamic stability. We evaluated the phase diagrams for 56,791 compounds to classify decomposition reactions into three types: 1. those that produce elemental phases, 2. those that produce compounds, and 3. those that produce both. This analysis shows that the decomposition into elemental forms is rarely the competing reaction that determines compound stability and that approximately two-thirds of decomposition reactions involve no elemental phases. Using experimentally reported formation enthalpies for 1,012 solid compounds, we assess the accuracy of the generalized gradient approximation (GGA) (PBE) and meta-GGA (SCAN) density functionals for predicting compound stability. For 646 decomposition reactions that are not trivially the formation reaction, PBE (MAD = 70 meV/atom) and SCAN (MAD = 59 meV/atom) perform similarly, and commonly employed correction schemes using fitted elemental reference energies make only a negligible improvement (~2 meV/atom). Furthermore, for 231 reactions involving only compounds (Type 2), the agreement between SCAN, PBE, and experiment is within ~35 meV/atom and is thus comparable to the magnitude of experimental uncertainty.
△ Less
Submitted 7 January, 2019; v1 submitted 18 October, 2018;
originally announced October 2018.
-
Physical descriptor for the Gibbs energy of inorganic crystalline solids and temperature-dependent materials chemistry
Authors:
Christopher J. Bartel,
Samantha L. Millican,
Ann M. Deml,
John R. Rumptz,
William Tumas,
Alan W. Weimer,
Stephan Lany,
Vladan Stevanović,
Charles B. Musgrave,
Aaron M. Holder
Abstract:
The Gibbs energy, G, determines the equilibrium conditions of chemical reactions and materials stability. Despite this fundamental and ubiquitous role, G has been tabulated for only a small fraction of known inorganic compounds, impeding a comprehensive perspective on the effects of temperature and composition on materials stability and synthesizability. Here, we use the SISSO (sure independence s…
▽ More
The Gibbs energy, G, determines the equilibrium conditions of chemical reactions and materials stability. Despite this fundamental and ubiquitous role, G has been tabulated for only a small fraction of known inorganic compounds, impeding a comprehensive perspective on the effects of temperature and composition on materials stability and synthesizability. Here, we use the SISSO (sure independence screening and sparsifying operator) approach to identify a simple and accurate descriptor to predict G for stoichiometric inorganic compounds with ~50 meV/atom (~1 kcal/mol) resolution, and with minimal computational cost, for temperatures ranging from 300-1800 K. We then apply this descriptor to ~30,000 known materials curated from the Inorganic Crystal Structure Database (ICSD). Using the resulting predicted thermochemical data, we generate thousands of temperature-dependent phase diagrams to provide insights into the effects of temperature and composition on materials synthesizability and stability and to establish the temperature-dependent scale of metastability for inorganic compounds.
△ Less
Submitted 6 January, 2019; v1 submitted 21 May, 2018;
originally announced May 2018.
-
New Tolerance Factor to Predict the Stability of Perovskite Oxides and Halides
Authors:
Christopher J. Bartel,
Christopher Sutton,
Bryan R. Goldsmith,
Runhai Ouyang,
Charles B. Musgrave,
Luca M. Ghiringhelli,
Matthias Scheffler
Abstract:
Predicting the stability of the perovskite structure remains a longstanding challenge for the discovery of new functional materials for many applications including photovoltaics and electrocatalysts. We developed an accurate, physically interpretable, and one-dimensional tolerance factor, τ, that correctly predicts 92% of compounds as perovskite or nonperovskite for an experimental dataset of 576…
▽ More
Predicting the stability of the perovskite structure remains a longstanding challenge for the discovery of new functional materials for many applications including photovoltaics and electrocatalysts. We developed an accurate, physically interpretable, and one-dimensional tolerance factor, τ, that correctly predicts 92% of compounds as perovskite or nonperovskite for an experimental dataset of 576 $ABX_3$ materials ($\textit{X} =$ $O^{2-}$, $F^-$, $Cl^-$, $Br^-$, $I^-$) using a novel data analytics approach based on SISSO (sure independence screening and sparsifying operator). τ is shown to generalize outside the training set for 1,034 experimentally realized single and double perovskites (91% accuracy) and is applied to identify 23,314 new double perovskites ($A_2$$\textit{BB'}$$X_6$) ranked by their probability of being stable as perovskite. This work guides experimentalists and theorists towards which perovskites are most likely to be successfully synthesized and demonstrates an approach to descriptor identification that can be extended to arbitrary applications beyond perovskite stability predictions.
△ Less
Submitted 6 January, 2019; v1 submitted 23 January, 2018;
originally announced January 2018.