Due to the need to evaluate the energy of a large number of structures during the automated search for a global minimum, the DFT calculations during the GOFEE runs have been carried out using low-precision parameters and a simplified model of the CeO2(111) surface. The best structures for each stoichiometry have been subsequently recalculated at increasingly computationally demanding (and more precise) levels of theory (see next subsection for more details).