Predicting challenging phase transitions with Bayesian active learning

Materials underpin modern technologies, from energy harvesting, storage, and conversion to information and communication technologies. Their functionality is often governed by the interplay between competing phases, as thermodynamic behavior shapes microscopic properties and ultimately determines technological performance; for instance, the light absorption of inorganic metal-halide perovskites in solar cells. Accurately predicting crystal thermodynamics, however, remains a major challenge for computational approaches because strong anharmonic effects require extensive sampling of the potential energy surface. Here, we present an on-the-fly Bayesian framework, combined with the stochastic self-consistent harmonic approximation, for learning first-principles interatomic potentials. This approach enables the prediction of thermodynamic properties over a broad temperature range with first-principles accuracy while requiring training on only a few tens to a few hundreds of atomic configurations. To demonstrate its power, we investigate the thermodynamic and dynamical properties of Li₂O, α-CsPbI₃, and δ-CsPbI₃, requiring only 44, 256, and 50 total-energy calculations, respectively. Notably, we show that this framework accurately captures the phase diagram of CsPbI₃, which explains its spontaneous degradation into the non-absorbing yellow phase, predicting the transition temperature with remarkable accuracy and efficiency. More broadly, the method presented opens a novel route toward accelerated materials engineering under realistic conditions for a wide range of technologically relevant applications, including solid-state batteries, optoelectronic devices, and memristors.