Abstract
The ground state structures of copper clusters with different sizes along with their aggregation have been systematic investigated using Amsterdam Density Functional (ADF) and Atomistix ToolKit (ATK) programs. On the basis of geometry optimization, some Cu clusters with more stable structures which were not reported previously have been revealed. In most cases, these Cu clusters prefer to adopt icosahedral structures which originate from the 13-atom icosahedron. It has also been demonstrated that the interaction between two Cu clusters is anisotropic, which is attributed to their charge distribution, especially the highest occupied molecular orbital (HOMO) and lowest unoccupied molecular orbital (LUMO) of Cu clusters. Moreover, we have carried out the simulation of Cu clusters aggregation on the silicone oil substrate by means of Monte Carlo (MC) method, which shows good consistence with our previous experimental studies.