To understand the interaction of metal on barium titanate surfaces and its structural, energetic, and electronic properties, the lowest-energy configuration for Cu/BaTiO3 system is crucial. The shapes of gas phase Cu12 clusters are largely dictated by relativistic effects. Some computational studies mention Cu12 in a crown-like geometry with C2v symmetry. Other literature use the Cu12 cluster. However, Cu12 nanoparticles have different shapes once adsorbed on surfaces. Figure shows that simply putting the C12 crown-like or cluster does not bring the most energetically favorable structure. Global optimization identifies the position and the shape of the copper cluster that is even 0.6 eV lower than the crown-like cluster. Surfaces have a significant role in the stability, organization, and chemical characteristics of metal clusters, particularly when chemical bonding or charge transfer occur. pyGlobOpt can be an effective tool for obtaining trustworthy structures of surface-supported clusters that display electrical and structural diversity. It is also remarkable that our method completed the whole process, except for DFT calculations, within minutes on a single CPU node.