摘要： The characterisation of the interaction between nano- or sub-nano- particles with a support nowadays increasingly relies on computational modelling by means of the density functional theory calculations. These provide valuable atomic-detail understanding of the structure and energetics of supported clusters, but it is still challenging to find (or design) structural models that are representative of real systems in terms of size, structure, and composition. In this study, we have applied an extensive and systematic approach combining global optimisation based on an evolutionary algorithm with atomistic ab initio thermodynamics for finding stable structures of a relevant material for catalytic methanol synthesis: Cu(111)-supported ZnyOx clusters. We identify the ZnO3 motif as the elementary building block of such clusters, on which we recently have investigated the full catalytic process for methanol synthesis. With the collection of global minima of Cu(111)-supported ZnyOx clusters resulting from this large-scale global optimisation effort, we assess the effect of size, gas-phase conditions, and support interactions on the phase diagrams, reactivity, and structural properties of the ZnyOx particles. We find moderate size-effects that are mostly related to the differences in stable Zn:O ratios of the identified global minima and to the formation of different sites in larger clusters. In contrast, large differences in the oxidation state of the clusters as defined by the gas-phase conditions significantly affect the geometry, electronic structure, and reactivity of the ZnyOx particles. This highlights the importance of thoroughly sampling structures with different stoichiometry and appropriately assessing their stability using a detailed thermodynamics analysis.