Adsorption of Water on Pristine Graphene: A van der Waals Density Functional Study with the vdW-C09 Approach

Understanding how water interacts with graphene at the molecular level is essential for advancing nanomaterial applications in filtration, catalysis, and environmental technologies. This study establishes a quantitative baseline for assessing how structural defects, dopants, or surface functionalization may enhance water adsorption, providing insights for the rational design of graphene-based materials in water purification, sensing, and nanofluidic applications. In this work, we employed density functional theory (DFT) with the vdW-C09 functional to investigate the adsorption of a single water molecule on pristine graphene, accurately accounting for long-range dispersion forces. Three high-symmetry adsorption sites-the center of the hexagonal ring, the C-C bond, and the top site-were explored in combination with three molecular orientations: Down, H-bond, and Up configurations. The calculated adsorption energies range from -93 to -145 meV (milli-electron volts), indicating that the interaction is dominated by weak van der Waals forces characteristic of physisorption. The most stable configuration corresponds to the Down orientation above the center of the hexagonal ring, with an adsorption energy of -145 meV and an equilibrium distance of 3.27 A (angstrom), defined as the vertical separation between the oxygen atom of the water molecule and the graphene plane. These results are in close agreement with previous theoretical studies and confirm the non-reactive and hydrophobic nature of pristine graphene.

💡 Research Summary

This paper presents a systematic first‑principles investigation of a single water molecule adsorbed on pristine graphene using density functional theory (DFT) with the van der Waals‑C09 (vdW‑C09) exchange‑correlation functional. The authors aim to establish a reliable quantitative baseline for water‑graphene interactions, which is essential for designing graphene‑based membranes, filters, and catalytic platforms for water treatment and environmental applications.

Computationally, the study employs the SIESTA code (version 4.1‑b4) with norm‑conserving Kleinman‑Bylander pseudopotentials and a double‑zeta plus polarization (DZP) basis set for carbon, hydrogen, and oxygen. A 6 × 6 graphene supercell containing 144 carbon atoms is used, together with a ~15 Å vacuum slab to eliminate spurious inter‑layer interactions. Spin polarization is enabled in all calculations, although the final ground state is non‑magnetic. The real‑space grid cutoff is set to 350 Ry, structural relaxations continue until forces fall below 0.05 eV/Å, and self‑consistent field (SCF) convergence criteria are 10⁻⁵ eV for total energy and density matrix. To correct for basis‑set superposition error (BSSE), the counter‑poise method is applied to every adsorption configuration.

Three high‑symmetry adsorption sites are examined: the center of a hexagonal carbon ring, the midpoint of a C–C bond, and the top of a carbon atom. For each site, three molecular orientations are considered: “Down” (hydrogens pointing toward the surface), “H‑bond” (one hydrogen oriented toward the surface while the oxygen points away), and “Up” (oxygen facing the surface). Geometry optimizations yield adsorption energies (E_ads) ranging from –93 meV to –145 meV and equilibrium O‑graphene distances (d) between 2.6 Å and 3.3 Å. The most stable configuration is the Down orientation above the hexagon center, with E_ads = –145 meV and d = 3.27 Å. The H‑bond orientation is moderately stable (E_ads ≈ –128 to –134 meV, d ≈ 2.6–2.9 Å), while the Up orientation is the least favorable (E_ads ≈ –93 to –98 meV, d ≈ 3.07–3.30 Å).

These energy values are characteristic of physisorption, confirming that van der Waals dispersion forces dominate the water‑graphene interaction. The relatively small energy differences among the three adsorption sites (≈20–50 meV) reflect the homogeneous, non‑polar nature of pristine graphene, which offers no strongly preferred binding positions. The authors attribute the enhanced stability of the Down orientation to favorable dipole‑induced‑dipole interactions between the water dipole (hydrogens positive) and the polarizable π‑electron cloud of graphene, in addition to dispersion forces.

The computed adsorption energies and distances are in excellent agreement with prior theoretical works employing other vdW‑corrected functionals (e.g., PBE‑D2, optB88‑vdW, vdW‑DF2) that reported values in the –100 to –150 meV range and equilibrium separations around 3.0 Å. Experimental observations, such as water contact angles on graphene exceeding 85°, also support the weak, hydrophobic character inferred from the calculations.

The paper emphasizes that, while pristine graphene interacts only weakly with water, the established reference data provide a solid foundation for future studies on defected, doped, or functionalized graphene. Introducing vacancies, hetero‑atom dopants, or oxygen‑containing groups is expected to increase local polarizability, create specific adsorption sites, and thereby raise the binding energy and possibly induce partial chemisorption. Such modifications could be strategically employed to tailor graphene’s wettability, enhance water uptake, and improve performance in desalination membranes, sensing devices, or catalytic reactors.

In conclusion, the vdW‑C09 functional accurately captures long‑range dispersion in the water‑graphene system, confirming that adsorption is governed by weak van der Waals forces and that the most favorable geometry is the Down orientation over the hexagonal ring center. The study delivers a rigorously benchmarked dataset that can be used as a reference point for assessing the impact of surface engineering on water adsorption. The work is funded by Brazilian agencies (FAPESB, FAPEMIG, CAPES, CNPq, INCT) and computational resources provided by the LNCC supercomputer.