Scalable properties of metal clusters: a comparative study of modern exchange-correlation functionals

The performance of eight generalized gradient approximation exchange-correlation (xc) functionals is assessed by a series of scalar relativistic all-electron calculations on octahedral palladium model clusters Pd(n) with n = 13, 19, 38, 55, 79, 147 and the analogous clusters Au(n) (for n up through 79). For these model systems, we determined the cohesive energies and average bond lengths of the optimized octahedral structures. We extrapolate these values to the bulk limits and compare with the corresponding experimental values. While the well-established functionals BP, PBE, and PW91 are the most accurate at predicting energies, the more recent forms PBEsol, VMTsol, and VT{84}sol significantly improve the accuracy of geometries. The observed trends are largely similar for both Pd and Au. In the same spirit, we also studied the scalability of the ionization potentials and electron affinities of the Pd clusters, and extrapolated those quantities to estimates of the work function. Overall, the xc functionals can be classified into four distinct groups according to the accuracy of the computed parameters. These results allow a judicious selection of xc approximations for treating transition metal clusters.