J. Phys. Chem. A (2001), 105 (33), 7905-7916.

DFT Calculations for Cu-, Ag- and Au-containing Molecules

F. Sue Legge, Graeme L. Nyberg* and J. Barrie Peel

Department of Chemistry, La Trobe University, Victoria, 3086, Australia

Received: January 17, 2001

The physico-chemical properties of a set of molecules containing Cu, Ag or Au atoms were calculated using the GAUSSIAN program suite, with the purpose of investigating the various density functional theory (DFT) approaches for subsequent application in cluster calculations. The test set comprised the copper-based molecules CuH, CuO, CuS, Cu₂, CuCl₂^-, CuCH₃, CuC₂H₂, Cu₂(HCO₂)₄(H₂O)₂ and Cu₆H₆(PH₃)₆ and the silver and gold diatomics AgH, AgO, AgS, Ag₂, AgCu, AuH, AuO, AuS, Au₂, AuCu and AuAg. The DFT methods used were SVWN, BLYP, and BPW91, together with the DFT-hybrids B3LYP and B3PW91. The calculations were carried out with the basis sets LANL2MB, LANL2DZ, 3-21G and 6-311G (when available). The properties calculated were the molecular geometry, vibrational frequencies and dissociation energies. It was found that all the DFT-based methods, particularly when allied with the LANL2DZ basis set, produced results which are significantly closer to experimental values than were those of the traditional Moller-Plesset (MP2) electron correlation method, and were also applicable to considerably larger molecules. Over the whole test set of molecules the RMS errors of the four 'BX' methods, in conjunction with LANL2DZ, were typically 3-4% for geometries, 6-8% for frequencies, and 10-16% for dissociation energies, with BPW91 and the popular B3LYP respectively at the lower and upper ends of these ranges, and with the errors being respectively over(+), under(-) and under(-)-estimates. The corresponding values for SVWN & MP2 were 2&6%, 12&12% and 33&42%, with these errors typically being +/- &+, +&- and +&-. While the best bond lengths are produced by SVWN (a local spin density approximation), which is not uncommmon, this advantage over the gradient corrected (BX) methods is only slight, and the latter are markedly superior for frequencies and especially dissociation energies. Not supported by the present results are the notions that (all) pure DFT methods underestimate metal-ligand bond lengths and overestimate bond strengths, and that hybrid methods are superior (and neither that DFT methods are overcorrelated). Testing on a subset of molecules with BPW91/LANL2DZ revealed no benefit in supplementing this basis set by the addition of diffuse functions, nor of the counterpoise correction. There appear to be specific incompatibilities with some method/basis-set combinations, and even the continuing availability of 3-21G for these metals is questionable. Because of its accuracy and reliability, the combination BPW91/LANL2DZ is recommended for these noble-metal systems (and to extensions such as the cluster-model approach to adsorption of a molecule on a metal surface).

Molecular Simulation (in press)

DFT modelling of Hydrogen on Cu(110) and (111)-type Clusters

Michelle J. S. Spencer and Graeme L. Nyberg

Department of Chemistry, La Trobe University, Victoria, 3086, Australia

Abstract

Density Functional Theory (DFT) calculations using Gaussian 98 have been performed on hydrogen adsorbed on clusters representing the (110) and (111) surfaces of Cu. Clusters were constructed to model different adsorption sites, and at least two different size clusters were used for each site. On the (111) surface, hydrogen prefers to adsorb in a hollow site, though with the hcp variant being favoured by the adsorption energy, and the fcc alternative by the vibrational frequencies. On the (110) surface, the 'fcc' site on a (1x2) reconstructed surface is preferred.

Keywords: Chemisorption; Cluster; Copper; Density Functional Theory calculations; Hydrogen