Using sound physical principles we modify the DFT-D2 atom pairwise semiempirical dispersion correction to density functional theory to work for metallic systems and in particular self-assembled monolayers of thiols on gold surfaces. We test our approximation for two functionals PBE-D and revPBE-D for lattice parameters and cohesive energies for Ni, Pd, Pt, Cu, Ag, and Au, adsorption energies of CO on (111) surfaces of Pd, Pt, Cu, Ag, and Au, and adsorption energy of benzene on Ag(111) and Au(111). Agreement with experimental data is substantially improved. We apply the method to self-assembled monolayers of alkanethiols on Au(111) and find reasonable agreement for PBE-D and revPBE-D for both physisorption of n-alkanethiols as well as dissociative chemisorption of dimethyl disulfide as an Au-adatom-dithiolate complex. By modifying the C6 coefficient for Au, we obtain quantitative agreement for physisorption and chemisorption for both PBE-D and revPBE-D using the same set of parameters. Our results confirm that inclusion of dispersion forces is crucial for any quantitative analysis of the thiol and thiolate bonds to the gold surface using quantum chemical calculations.
Modelling excited states of weakly bound complexes with density functional theory
