Neutral and charged gallium clusters: structures, physical properties and implications for the melting features

We report the putative Global Minimum (GM) structures and electronic properties of Ga(N)(+), Ga(N) and Ga(N)(-) clusters with N = 13-37 atoms, obtained from first-principles density functional theory structural optimizations. The calculations include spin polarization and employ an exchange-correlation functional which accounts for van der Waals dispersion interactions (vdW-DFT). We find a wide diversity of structural motifs within the located GM, including decahedral, polyicosahedral, polytetrahedral and layered structures. The GM structures are also extremely sensitive to the number of electrons in the cluster, so that the structures of neutral and charged clusters differ for most sizes. The main magic numbers (clusters with an enhanced stability) are identified and interpreted in terms of electronic and geometric shell closings. The theoretical results are consistent with experimental abundance mass spectra of Ga(N)(+) and with photoelectron spectra of Ga(N)(-). The size dependence of the latent heats of melting, the shape of the heat capacity peaks, and the temperature dependence of the collision cross-sections, all measured for Ga(N)(+) clusters, are properly interpreted in terms of the calculated cohesive energies, spectra of configurational excitations, and cluster shapes, respectively. The transition from "non-melter" to "magic-melter" behaviour, experimentally observed between Ga(30)(+) and Ga(31)(+), is traced back to a strong geometry change. Finally, the higher-than-bulk melting temperatures of gallium clusters are correlated with a more typically metallic behaviour of the clusters as compared to the bulk, contrary to previous theoretical claims.