[Paper Review] Interaction between a Water Molecule and a Graphite Surface

D. J. W. Geldart, I Wayan Sudiarta+2 more|ArXiv.org|Jan 18, 2008

Advanced Chemical Physics Studies3 citations

TL;DR

This study uses ab initio computational methods to decompose the interaction energy between a water molecule and a graphite surface into electrostatic, induction, Pauli repulsion, and correlation components. It identifies perimeter-specific contributions that can be isolated and discarded, enabling accurate extrapolation to the bulk graphite limit and yielding an explicit analytical potential energy surface with a minimum energy configuration where both H atoms of water point symmetrically away from the surface, resulting in an electronic interaction energy of -16.8 ± 1.7 kJ/mol.

ABSTRACT

The interaction energy between a water molecule and graphitic structured clusters terminated by hydrogen atoms is analyzed by ab initio methods and decomposed into electrostatic, induction, Pauli repulsion, and correlation energy contributions. Contributions to the energy which are due solely to the perimeter of the clusters are identified. These can be isolated and discarded which greatly simplifies the problem of extrapolation to the large cluster limit. The remaining terms are intrinsic to the interaction of a water molecule with real graphitic layers and an explicit analytical form is given for the potential energy surface. The minimum energy configuration is found to have both hydrogen atoms of the water molecule pointing symmetrically away from the graphitic plane. The electronic interaction in this mode is -16.8 +/- 1.7 kJ/mol for water-graphite and the zero point energy is estimated as 1.3 kJ/mol.

Motivation & Objective

• To understand the quantum mechanical interaction between a single water molecule and a graphitic surface.
• To decompose the total interaction energy into physically meaningful components: electrostatic, induction, Pauli repulsion, and correlation.
• To identify and isolate contributions arising solely from the cluster perimeter, which distort convergence toward the bulk limit.
• To derive an explicit analytical form for the potential energy surface of water on graphite.
• To determine the most stable molecular orientation and binding energy of water on graphite.

Proposed method

• Ab initio quantum chemistry calculations were performed on hydrogen-terminated graphitic clusters of increasing size.
• Energy decomposition analysis (EDA) was applied to separate the total interaction energy into electrostatic, induction, Pauli repulsion, and correlation contributions.
• Perimeter-specific energy contributions were identified by comparing clusters of different sizes and shapes.
• These perimeter terms were isolated and removed from the energy data to enable reliable extrapolation to the infinite graphite surface limit.
• The remaining intrinsic interaction energy was fitted to derive an analytical potential energy surface.
• Zero-point energy corrections were estimated to refine the electronic binding energy.

Experimental results

Research questions

• RQ1 What are the relative contributions of electrostatic, induction, Pauli repulsion, and correlation to the water-graphite interaction energy?
• RQ2 How do edge and perimeter effects in finite graphite clusters influence the convergence of interaction energy toward the bulk graphite limit?
• RQ3 What is the most stable orientation of a water molecule on a graphite surface, and what is its binding energy?
• RQ4 Can perimeter contributions be systematically isolated and removed to improve extrapolation accuracy for large-scale surface interactions?
• RQ5 What analytical form best describes the potential energy surface of a water molecule on a graphite surface?

Key findings

• The minimum energy configuration of the water molecule on graphite features both hydrogen atoms oriented symmetrically away from the surface plane.
• The electronic interaction energy in this configuration is -16.8 ± 1.7 kJ/mol, indicating a favorable physisorption interaction.
• Perimeter contributions to the interaction energy were identified and isolated, enabling their removal to improve convergence toward the bulk limit.
• The remaining intrinsic interaction energy is dominated by dispersion and induction effects, with significant electrostatic contributions.
• The zero-point energy correction was estimated at 1.3 kJ/mol, slightly reducing the magnitude of the total binding energy.
• An explicit analytical form for the potential energy surface was derived, enabling accurate modeling of water adsorption on graphite.