Abstract
Experimental and theoretical works have, to date, been unable to uncover the ground state configuration of prominent solid electrolyte candidate cubic Li$_7$La$_3$Zr$_2$O$_{12}$ (c-LLZO). Computational studies rely on an initial low-energy structure as a reference point. In this study, we present a methodology to identify energetically favourable configurations of c-LLZO, enabling the isolation of low-energy structures, for a crystallographically predicted structure. We begin by eliminating structures that involve overlapping Li atoms based on nearest neighbour counts. We further reduce the configuration space by eliminating symmetry images from all remaining structures. This is followed up with a machine learning-based energetic ordering of all remaining structures. By considering the geometrical constraints that emerge from this methodology we determine that a large portion of previously reported structures may not be feasible or stable. The method developed here could be extended to other ion conductors and partially occupied crystals. Furthermore, we provide all structures generated in a freely accessible database with the aim to improve accuracy and reproducibility in future c-LLZO research.