Article about Machine Learning for Exploring Small Polaron Configurational Space Published in npj Computational Materials

Taming complexity in polaron configurational space: TACO scientists present a new machine-learning approach to accelerate polaron distribution calculation in recently published npj Computational Materials publication.

Polarons are quasiparticles that arise from the interaction of charges (electrons or holes) with the lattice vibrations (phonons) in a crystal. Polarons are a pervasive phenomenon in transition-metal oxides. They have a strong impact on the properties and functionality of materials that are important in many technical applications. For many years, polarons have been at the focus of theoretical research to understand better and predict those material properties. Despite the huge importance and considerable efforts, detailed predictions have been stalled simply by the overwhelmingly large number of possible configurations that polarons can occupy in a crystal lattice.

A team of TACO scientists has now accepted the challenge to tame this complexity: In a just-published paper, the team led by Ph.D. students Viktor Birschitzky and Florian Ellinger, under the guidance of the TACO Co-PI Michele Reticcioli and PIs Cesare Franchini and Ulrike Diebold, present a new machine-learning (ML) approach to accelerating significantly the previously prohibitively demanding calculations. The scientists trained the ML algorithm on a small sample of detailed yet computationally expensive DFT first-principles calculations (molecular dynamics and random sampling) of polaron configurations. Later, the algorithm explored the vast configurational space to identify energetically favourable and thus more likely polaron arrangements at variable polaron concentrations. Where less than 500 polaron configurations could be calculated with previous first-principle approaches, millions of structures can now be scanned and explored with the new ML approach.

The researchers validated the approach on two well-known material surface systems: rutile TiO2 (110), and Nb-doped perovskite SrTiO3 (001). Particularly for this latter surface system, the ML-aided procedure achieved a systematic improvement of the polaronic ground state configurations compared to the more traditional approaches. The researchers are optimistic that the algorithm can be adapted relatively easily to consider other types of crystal defects, such as adatoms. Also, the scalability of the code for direct comparison of the simulation results with experimental data is supposed to be good. This would open the path for applying the code to many more material systems.

The paper was published open-access on June 6 in npj Computational Materials, a member of the Nature family of journals.

Abstract: Polaron defects are ubiquitous in materials and play an important role in many processes involving carrier mobility, charge transfer and surface reactivity. Determining small polarons’ spatial distributions is essential to understand materials properties and functionalities. However, the required exploration of the configurational space is computationally demanding when using first principles methods. Here, we propose a machine-learning (ML) accelerated search that determines the ground state polaronic configuration. The ML model is trained on databases of polaron configurations generated by density functional theory (DFT) via molecular dynamics or random sampling. To establish a mapping between configurations and their stability, we designed descriptors modelling the interactions among polarons and charged point defects. We used the DFT+ML protocol to explore the polaron configurational space for two surface-systems, reduced rutile TiO2(110) and Nb-doped SrTiO3(001). The ML-aided search proposes additional polaronic configurations and can be utilized to determine optimal polaron distributions at any charge concentration.
The full article can be found here (or directly as PDF).

Authors:
Viktor Birschitzky, Florian Ellinger, Ulrike Diebold, Michele Reticcioli, and Cesare Franchini

Subprojects:
P07 – Polaron pattern recognition in correlated oxide surfaces
P02 – Surface structure and reactivity of multi-component oxides at the atomic scale