DISTRIBUTION OF INTERATOMIC DISTANCES IN LARGE METALLIC CLUSTERS

Spherically averaged pseudopotential (SAPS) calculations have been done for Mg(n) clusters, with n up to 250 within the framework of density functional theory. The electronic structure is computed resorting to the Thomas-Fermi-Dirac-Weizsacker (TFDW) approximation for the kinetic energy. The equilibrium geometries have been obtained by minimizing the total cluster energy with respect to the atomic positions using the steepest-descent method. The ground state geometries obtained in this way are formed by spherical atomic shells, the number of them increasing with cluster size, up to a number of four for the biggest sizes considered here. An analysis of the distribution of the interatomic distances shows that the more internal is the shell, the more contracted are the interatomic distances. This effect diminishes progressively with increasing cluster size. For the purpose of comparison, similar calculations have been done with Cs(n) clusters in the same size range, allowing us to reproduce previous results obtained using a more elaborated density functional technique (Kohn-Sham method). The inhomogeneous contraction of interatomic distances then appears as a general fact for simple metallic clusters and not only for alkaline ones.