From clusters to bulk: A relativistic density functional investigation on a series of gold clusters Aun, n=6,…,147

A series of gold clusters spanning the size range from Au6 through Au147 (with diameters from 0.7 to 1.7 nm) in icosahedral, octahedral, and cuboctahedral structure has been theoretically investigated by means of a scalar relativistic all-electron density functional method. One of the main objectives of this work was to analyze the convergence of cluster properties toward the corresponding bulk metal values and to compare the results obtained for the local density approximation (LDA) to those for a generalized gradient approximation (GGA) to the exchange-correlation functional. The average gold–gold distance in the clusters increases with their nuclearity and correlates essentially linearly with the average coordination number in the clusters. An extrapolation to the bulk coordination of 12 yields a gold–gold distance of 289 pm in LDA, very close to the experimental bulk value of 288 pm, while the extrapolated GGA gold–gold distance is 297 pm. The cluster cohesive energy varies linearly with the inverse of the calculated cluster radius, indicating that the surface-to-volume ratio is the primary determinant of the convergence of this quantity toward bulk. The extrapolated LDA binding energy per atom, 4.7 eV, overestimates the experimental bulk value of 3.8 eV, while the GGA value, 3.2 eV, underestimates the experiment by almost the same amount. The calculated ionization potentials and electron affinities of the clusters may be related to the metallic droplet model, although deviations due to the electronic shell structure are noticeable. The GGA extrapolation to bulk values yields 4.8 and 4.9 eV for the ionization potential and the electron affinity, respectively, remarkably close to the experimental polycrystalline work function of bulk gold, 5.1 eV. Gold 4f core level binding energies were calculated for sites with bulk coordination and for different surface sites. The core level shifts for the surface sites are all positive and distinguish among the corner, edge, and face-centered sites; sites in the first subsurface layer show still small positive shifts.