High-entropy perovskite oxides (HEPOs) have emerged as a promising class of materials for use in solid oxide electrolyzer cells (SOECs), particularly as oxygen electrodes, due to their excellent thermal and chemical stability. The concentration of oxygen vacancies in HEPOs directly impacts critical properties such as ionic conductivity, thermal expansion, and overall electrochemical behavior. Despite their importance, predicting how the oxygen vacancy concentration in HEPOs varies with five or more metal cation dopants and temperature remains elusive, as the net energy needed to remove a lattice oxygen is significantly complicated by variations from different metal cations and configurational entropy.In this work, we investigate how oxygen vacancy concentrations are affected by low-entropy (LE, 2–3 cations) and high-entropy (HE, 5+ cations) A-site compositions in perovskite oxides. We synthesized 14 compositions, inspired by the state-of-the-art SOEC material lanthanum strontium cobalt ferrite (LSCF), using solution combustion synthesis and measured changes in oxygen vacancy concentration between 500–1000°C via thermogravimetric analysis. Using established defect models, we estimated and compared the oxygen vacancy formation enthalpy (ΔH) and entropy (ΔS).Our results confirm the expected role of divalent A-site cations in facilitating oxygen vacancy formation. However, substantial variation was observed among HE samples with identical divalent cation concentrations, suggesting additional influencing factors. We identify the variance in A-site ionic radii as a key compositional descriptor that strongly correlates with oxygen vacancy thermodynamics. HE samples with larger A-site size variance consistently exhibited lower ΔH and ΔS of oxygen vacancy formation, as well as more linear oxygen vacancy growth with temperature compared to their LE counterparts.To investigate the origin of this behavior, we conducted high-throughput atomistic simulations using a machine-learned universal interatomic potential. We computed the energy required to form an oxygen vacancy across hundreds of oxygen sites in large supercells (1280 atoms) for each material. While the average vacancy formation energies were comparable between LE and HE samples, HE compositions exhibited a broader distribution of vacancy formation energies due to the local strain induced by size mismatches between A-site cations.We developed a statistical thermodynamic model to analyze how variations in formation energy impact a material's equilibrium oxygen vacancy concentration. By treating oxygen sites as a distribution of energy states, we derived analytical expressions for formation ΔH and ΔS as a function of vacancy energy variance. The results show that broadening the vacancy energy distribution leads to a decrease in both formation enthalpy and configurational entropy, consistent with experimental data. Notably, the model predicts a scaling relation between ΔHf and ΔSf, explaining the experimentally observed correlation between LE and HE samples.In summary, our combined experimental, computational, and theoretical approach reveals that A-site cation disorder—quantified by size variance—plays a crucial role in determining the thermodynamics of oxygen vacancy formation in HEPOs. This disorder induces variation in oxygen vacancy formation energies, which lowers the formation enthalpy (ΔH) and entropy (ΔS), leading to greater oxygen vacancy populations at lower temperatures. The findings establish A-site size variance as a practical and predictive design parameter for tuning defect behavior in complex oxides. This framework advances our understanding of defect chemistry in disordered systems and provides guidance for engineering functional materials for SOECs, catalysis, and beyond. Figure 1
Read more