Theoretical study of the adsorption of aromatic units on carbon allotropes including explicit (empirical) DFT dispersion corrections and implicitly dispersion-corrected functionals: the pyridine case

Abstract

The suitability of implicitly dispersion-corrected functionals, namely the M06-2X, for the determination of interaction energies and electron polarization densities in adsorption studies of aromatic molecules on carbon allotropes surfaces is analysed by comparing the results with those obtained using explicit dispersion through Grimme’s empirical corrections. Several models of increasing size for the graphene sheet together with one-dimensional curved carbon structures, (5,5), (6,6) and (7,7) armchair single-walled nanotubes, and two-dimensional curved carbon structures, C60 fullerene, have been considered as substrates in this work, whereas pyridine has been chosen as an example for the adsorbed aromatic molecule. Comparison with recent experimental estimations of the adsorption energy and calculations using periodic boundary conditions on a supercell of 72 carbon atoms indicates that a finite model containing ninety six carbon atoms (C96) approaches quite well the adsorption on a graphene sheet. Analysis of the interaction energy components reveals that the M06-2X functional accounts for most of the dispersion energy implicitly, followed far by wB97X and B3LYP, whereas B97 and BLYP do not differ too much from HF. It has been found that M06-2X corrects only the energy component associated to dispersion and leaves the rest, electrostatic, Pauli and induction “unaltered” with respect to the other DFT functionals investigated. Moreover, only the M06-2X functional reflects the effect of dispersion on the electron polarization density, whereas for the remaining functionals the polarization density does not differ too much from the HF density. This makes the former functional more suitable a priori for the calculation of electron density related properties in these adsorption complexes.