Adsorption
is a scientifically and technologically important interfacial
phenomenon, which however presents challenges to conventional density
functional theory (DFT) due to the long-range van der Waals (vdW)
interactions. We have developed a model of long-range vdW correction
for physisorption of graphene (G) on metals with the Lifshitz–Zaremba–Kohn
second-order perturbation theory, by incorporating dipole- and quadrupole-surface
interactions and screening effects. The physisorption energies calculated
by the model between graphene and eight metal surfaces (Al, Ni, Co,
Pd, Pt, Cu, Ag, and Au), and the adsorption energies for the same
G/metal structures from self-consistent DFT PBE (Perdew–Burke–Ernzerhof)
calculations, are obtained in a range of distances between G and the
metal surfaces. The sum of these two parts is the total adsorption
energy as a function of the distance, from which the equilibrium distance
and the binding energy are determined simultaneously. The results
show high accuracy, with the mean absolute error (MAE) of binding
energy of 7 meV and the MAE of equilibrium distance of 0.2 Å,
significantly improving upon other vdW methods. The PBE + vdW binding
energy curves give better fits to the random phase approximation curves
around the equilibrium distances than do the curves of other methods
considered here. The higher-order quadrupole-surface correction is
important and accounts for about 30% of the total vdW correction.