How to predict a new material: a marriage of machine learning and quantum mechanics ✛ ✦

Original: "Data-Driven Prediction of NaCl-Type Entropy-Stabilized Oxide Compositions from First-Principles and Supervised Learning"
arXiv:2607.04502 · 2026-07-05 · CC BY 4.0 · 5 min
High-throughput screening of entropy-stabilized oxides using DFT and a neural network paves the way to accelerated design of functional materials.

Context

The search for new inorganic materials with unique properties is a central challenge of modern solid-state chemistry. The traditional "synthesize–test" approach is extremely labor-intensive, especially when it comes to multicomponent systems: the space of possible elemental combinations is vast, like finding a needle in a haystack. A breakthrough came with the discovery of high-entropy oxides (HEOs)—solid solutions of five or more metals in equimolar ratios, stabilized precisely by high configurational entropy of mixing. As Ludwig Boltzmann once foresaw, disorder can be a driving force for stability—the second law of thermodynamics kicks in, making the high-entropy state favorable at elevated temperatures. Until now, however, there were no reliable methods for rapidly screening thousands of possible compositions to determine whether they form a single-phase ESO or decompose into a mixture of simpler oxides. This limitation held back practical application of the entropy-stabilization concept.

Methods

The authors built a multi-layered virtual pipeline. In the first stage, using numerical simulation within density functional theory (DFT), they calculated the formation enthalpies of all possible ordered oxides—binary and ternary—for 16 selected cations: from alkaline-earth to transition metals (including iron, cobalt, nickel, copper, zinc) and post-transition elements. Special attention was paid to elements with high metallicity—that is, metals capable of forming oxides with variable valence. Disorder in five-component solid solutions was modeled using special quasi-random supercells (SQS) of up to 60 atoms, which reproduce a statistically random distribution of cations on the sublattice against a rigid anion lattice of oxygen. To train the models, full DFT calculations were performed on 416 five-component configurations (about one million core-hours). Then machine learning entered the scene: three models (linear regression, random forest, and a perceptron neural network) were trained to predict the key parameter—the distance to the convex stability hull in multidimensional composition space and the derived stabilization temperature Tstab. The descriptors used were atomic number, radius, electronegativity, and number of valence electrons. The thermodynamic balance was described via quantum thermodynamics: the free energy of the mixture included the DFT enthalpy and the ideal configurational entropy obeying the Boltzmann distribution.

Results

The best performance was delivered by a multilayer perceptron (MLP): on a test set of 104 systems, the error was only 4.2 kJ/mol (vs. 6.1 for the linear model). This means the neural network captures subtle nonlinear relationships between the chemical nature of cations and oxide stability. Applying the MLP to all 4368 possible equimolar combinations, the researchers obtained a ranked list of candidates. Already known ESOs, such as (Co,Cu,Mg,Ni,Zn)O, were among the leaders: after correction with the more accurate r2SCAN functional, its Tstab dropped from 2745 K to 1249 K, approaching the experimental value of ~1100 K. New promising compositions included (Ca,Cu,Fe,Mg,Zn)O with a Tstab around 1676 K and (Ca,Cu,In,Sr,Zn)O at 1750 K. Experimental verification of three selected compositions confirmed the main trends: in some cases a rock-salt phase appeared (though with partial calcium segregation), in others melting occurred before the calculated stabilization, indicating the need to account for the liquid phase. The model correctly predicted the dominant competing phases.

Implications

The developed approach changes the paradigm for discovering new entropy-stabilized oxides. Instead of blind trial and error, experimenters now have a map of prospective compositions based on rigorous physical principles and supported by quantum calculations. This not only saves years of lab work but also enables rational interpretation of synthesis failures—for instance, revealing which competing phases hinder the formation of the target material. From a fundamental standpoint, the work confirms that the stability of high-entropy oxides is not solely determined by configurational entropy: enthalpic factors (charge transfer, oxidation state compatibility, local distortions) play an equally important role, and it is precisely their nonlinear combinations that the neural network successfully learns.

Future development

The next generation of models should include thermal effects beyond ideal configurational entropy: vibrational entropy (phonon calculations), magnetic contributions, and the possibility of non-equimolar compositions. It is also important to integrate the thermodynamics of the liquid phase and more accurate exchange-correlation functionals (hybrid, with strong correlations) to reduce systematic errors in absolute Tstab values. In the future, the method could be extended to other structure types (fluorite, perovskite, spinel) and even more multicomponent systems.

Impact

The results will directly influence materials design in several cutting-edge areas: in energy—the search for new cathodes for lithium-ion batteries with high ionic conductivity; in catalysis—the creation of high-temperature stable catalysts; in electronics—multifunctional oxides with tunable dielectric and magnetic properties.

Next steps

The immediate next steps will be targeted synthesis and characterization of the predicted top candidates, along with phonon calculations to refine the vibrational entropy of the selected compositions.

Key open problems

The article is directly linked to the unsolved problem of predicting phase stability in multicomponent disordered solids. An open question remains: what is the true role of configurational entropy compared to enthalpic and kinetic factors? Can universal descriptors be found that a priori assess the propensity of a set of elements to form a single-phase high-entropy oxide? These questions lie at the heart of rational materials design—one of the main goals of modern solid-state physics.

The knowledge machine's view

This section is addressed to the author of the paper and their colleagues. Neighbouring works were found by meaning-based search across our archive, not by word overlap. We do not review: we show the field around the work.

What the work does

The authors built a virtual pipeline to predict the stability of high-entropy oxides: DFT calculations of enthalpies for 16 cations, disorder modeling via SQS, and training of three ML models on 416 five-component configurations. The best model (MLP) achieved an error of 4.2 kJ/mol in predicting the distance to the convex hull and the stabilization temperature.

Where its strength is

Combining DFT calculations with machine learning on descriptors derived from fundamental elemental properties allowed rapid stability assessment across a vast compositional space. This provides a practical tool for screening new ceramics without exhaustive enumeration.

Where it could go next

Worth trying

  • Try active learning: iteratively select compositions with the highest model uncertainty for additional DFT calculations.
  • Expand the descriptor set by adding, for example, ionization energies or crystal-field parameters, as done in related works.
  • Test the transferability of the trained model to other classes of high-entropy materials, such as nitrides or carbides.

Why this may matter

Fast and accurate screening of high-entropy oxides could accelerate the discovery of materials for energy and high-temperature applications, shortening the cycle from idea to synthesis.

Original abstract

as written by the authors arXiv:2607.04502

Entropy-stabilized oxides (ESOs) open access to vast multicomponent compositional spaces, but identifying promising candidates remains challenging because of the large number of possible mixtures and the need to assess their stability against competing phases. In this work, we develop a high-throughput computational framework to screen equimolar quinary ESOs in the NaCl structure type by combining density functional theory (DFT), special quasirandom structures (SQS), convex-hull thermodynamics, and supervised machine learning. A consistent reference database of binary and ternary ordered oxides, including disordered phases such as all binary cation combinations in the NaCl-type oxide, is first constructed using GGA and meta-GGA calculations. Quinary disordered phases are then described by SQS supercells and used to train machine-learning models that predict the distance to the convex hull and the corresponding stabilization temperature over the full set of 4368 possible equimolar quinary compositions generated from 16 cation species. Among the tested models, an optimized multilayer perceptron provides the best predictive performance, with a test error of about 4 kJ/mol, while requiring explicit DFT calculations for only about 10% of the quinary systems. Comparison with experimental synthesis tests and computed decomposition paths further shows that the approach captures the main stability trends and the dominant competing phases, although absolute stabilization temperatures remain affected by systematic thermodynamic approximations. These results establish an efficient route for the data-driven exploration of multicomponent oxides and provide practical guidance for the experimental search for new ESOs.

If you decided to test all 4368 five-component combinations in the lab, spending just one week on synthesis and analysis of each, you would need about 84 years of continuous work. The neural network, after training, accomplishes the task in a fraction of a second.

\Delta S_{\text{conf}} = -R \sum_{i=1}^{n} y_i \ln y_i
Configurational entropy of mixing for n components with mole fractions y_i and the universal gas constant R.
\Delta_{\text{hull}} H(\mathbf{x}_\varphi) = \Delta_f H(\mathbf{x}_\varphi) - E_{\text{GS}}(\mathbf{x}_\varphi)
Difference between the formation enthalpy of a compound and the ground-state energy (convex hull) at the same composition; if >0, the compound is unstable at 0 K.

Key numbers

  • number of studied five-component systems: 4368
  • root-mean-square error of MLP model: 4.2 kJ/mol
  • number of cations in the library: 16
  • supercell size for modeling disorder: 60 atoms
  • stabilization temperature of (Co,Cu,Mg,Ni,Zn)O (after correction): 1249 K
Scientists
Ludwig BoltzmannRudolf ClausiusJacob BekensteinStephen HawkingJames Clerk MaxwellEdward Witten
Concepts
ENSO Multivariate ENSO Index