Assessment of Density Functional Theory in Predicting Interaction Energies Between Water and Polycyclic Aromatic Hydrocarbons: From Water on Benzene to Water on Graphene

27 November 2018, Version 1
This content is a preprint and has not undergone peer review at the time of posting.


The interaction of water with polycyclic aromatic hydrocarbons, from benzene to graphene, is investigated using various exchange-correlation functionals selected across generalized gradient approximation (GGA), meta-GGA, and hybrid families within the density functional theory (DFT) hierarchy. The accuracy of the different functionals is assessed through comparisons with high-level electronic structure methods, including random phase approximation (RPA), diffusion Monte Carlo (DMC), and coupled-cluster with single, double, and perturbative triple excitations (CCSD(T)). Relatively large variations are found in the interaction energies predicted by different DFT models, with GGA functionals underestimating the interaction strength for configurations with the water oxygen pointing toward the aromatic molecules, and the meta-GGA B97M-rV and hybrid ωB97M-V functionals providing nearly quantitative agreement with CCSD(T) values available for the water-benzene, water-coronene, and water-circumcoronene dimers, which, in turn, are within ∼1 kcal/mol of the corresponding RPA and DMC results. Similar trends among GGA, meta-GGA, and hybrid functionals are observed for the larger polycyclic aromatic hydrocarbon molecules considered in this analysis (up to C216H36). By performing absolutely localized molecular orbital energy decomposition analyses (ALMO-EDA) of the DFT results, it is found that, independently of the number of carbon atoms and exchange-correlation functional, the dominant contributions to the interaction energies between water and polycyclic aromatic hydrocarbon molecules are the electrostatic and dispersion terms while polarization and charge transfer effects are negligibly small. Calculations carried out with GGA and meta-GGA functionals indicate that, as the number of carbon atoms increases, the interaction energies slowly converge to the corresponding values obtained for an infinite graphene sheet. Importantly, water-graphene interaction energies calculated with the B97M-rV functional appear to deviate by more than 1 kcal/mol from the available RPA and DMC values.


polyaromatic hydrocarbons
molecular interactions
density functional theory
random phase approximation
coupled cluster theory


Comments are not moderated before they are posted, but they can be removed by the site moderators if they are found to be in contravention of our Commenting Policy [opens in a new tab] - please read this policy before you post. Comments should be used for scholarly discussion of the content in question. You can find more information about how to use the commenting feature here [opens in a new tab] .
This site is protected by reCAPTCHA and the Google Privacy Policy [opens in a new tab] and Terms of Service [opens in a new tab] apply.