Ab initio study of water anchored in graphene pristine and vacancy-type defects.

Abstract

In this paper, we have addressed two issues that are relevant to the interaction of water in pristine and vacant graphene through first-principles calculations based on the Density Functional Theory (DFT). The results showed that for the interaction of pristine graphene with water, the DOWN configuration (with the hydrogen atoms facing downwards) was the most stable, presenting binding energies in the order of -13.62 kJ/mol at a distance of 2.375 Å in the TOP position. We also evaluated the interaction of water with two vacancy models, removing one carbon atom (Vac-1C) and four atoms (Vac-4C). In the Vac-1C system, the most favourable system was the DOWN configuration, with binding energies ranging from -20.60 kJ/mol to -18.41 kJ/mol in the TOP and UP positions, respectively. A different behaviour was observed for the interaction of water with Vac-4C; regardless of the configuration of the water, it is always more favourable for the interaction to occur through the vacancy centre, with binding energies ​​between -13.28 kJ/mol and -20.49 kJ/mol. Thus, the results presented open perspectives for the technological development of nanomembranes as well as providing a better understanding of the wettability effects of graphene sheets, whether pristine or with defects. We evaluated the interaction of pristine and vacant graphene with the water molecule, through calculations based on Density Functional Theory (DFT); implemented by the SIESTA program. The electronic, energetic, and structural properties were analyzed by solving self-consistent Kohn-Sham equations. In all calculations, a double ζ plus a polarized function (DZP) was used for the numerical baise set. Local Density Approximation (LDA) with the Perdew and Zunger (PZ) parameterisation along with a basis set superposition error (BSSE) correction were used to describe the exchange and correlation potential (Vxc). The water and isolated graphene structures were relaxed until the residual forces were less than 0.05 eV/Å-1 in all atomic coordinates.