Broken-Symmetry Unrestricted Hybrid Density Functional Calculations on Nickel Dimer and Nickel Hydride

In the present work we investigate the adequacy of broken-symmetry unrestricted density functional theory (DFT) for constructing the potential energy curve of nickel dimer and nickel hydride, as a model for larger bare and hydrogenated nickel cluster calculations. We use three hybrid functionals: the popular B3LYP, Becke's newest optimized functional Becke98, and the simple FSLYP functional (50% Hartree-Fock and 50% Slater exchange and LYP gradient-corrected correlation functional) with two basis sets: all-electron (AE) Wachters+f basis set and Stuttgart RSC effective core potential (ECP) and basis set. We find that, overall, the best agreement with experiment, comparable to that of the high-level CASPT2, is obtained with B3LYP/AE, closely followed by Becke98/AE and Becke98/ECP. FSLYP/AE and B3LYP/ECP give slightly worse agreement with experiment, and FSLYP/ECP is the only method among the ones we studied that gives an unaceptably large error, underestimating the dissociation energy of nickel dimer by 28%, and being in the largest disagreement with the experiment and the other theoretical predictions.