These are preliminary reports that have not been peer-reviewed. They should not be regarded as conclusive, guide clinical practice/health-related behavior, or be reported in news media as established information. For more information, please see our FAQs.
water_on_graphene.pdf (3.17 MB)

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

submitted on 26.11.2018 and posted on 27.11.2018 by Adeayo Ajala, Vamsee K. Voora, Narbe Mardirossian, Filipp Furche, Francesco Paesani

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.


National Science Foundation

Air Force Office of Scientific Research


Email Address of Submitting Author


University of California, San Diego


United States

ORCID For Submitting Author


Declaration of Conflict of Interest

No conflict of interests