We identify candidate low-energy cluster structures using one of the four methods listed above. DFT calculations were performed on these clusters structures before adding them to the database. An outline of the high-throughput workflow used in these DFT calculations is provided in Figure . In the first step we check whether the ionic relaxation has converged. If it has not, we rerun the relaxation with modified VASP settings such as increased wall time or a different optimization algorithm.