Water Clusters on Graphite: Methodology for Quantum Chemical A Priori Prediction of Reaction Rate Constants

Ab initio insight into the electrolysis of water on basal and edge (fullerene C20) surfaces of 4 Å single-walled carbon nanotubes

The extreme surface reactivity of 4 Å single-walled carbon nanotubes (SWCNTs) makes for a very promising catalytic material, however, controlling it experimentally has been found to be challenging. Here, we employ ab initio calculations to investigate the extent of surface reactivity and functionalization of 4 Å SWCNTs. We study the kinetics of water dissociation and adsorption on the surface of 4 Å SWCNTs with three different configurations: armchair (3,3), chiral (4,2) and zigzag (5,0). We reveal that out of three different configurations of 4 Å SWCNTs, the surface of tube (5,0) is the most reactive due to its small HOMO-LUMO gap. The dissociation of 1 H2O molecule into an OH/H pair on the surface of tube (5,0) has an adsorption energy of -0.43 eV and an activation energy barrier of 0.66 eV at 298.15 K in pure aqueous solution, which is less than 10% of the activation energy barrier of the same reaction without the catalyst present. The four steps of H+/e- transfer in the oxygen evolution reaction have also been studied on the surface of tube (5,0). The low overpotential of 0.38 V indicates that tube (5,0) has the highest potential efficiency among all studied carbon-based catalysts. We also reveal that the armchair edge of tube (5,0) is reconstructed into fullerene C20. The dangling bonds on the surface of fullerene C20 result in a more reactive surface than the basal surface of tube (5,0), however the catalytic ability was also inhibited in the later oxygen reduction processes.

Exploring the effect of oxygen-containing functional groups on the water-holding capacity of lignite

Graphene oxide with different degrees of oxidation was prepared and selected as a model compound of lignite to study quantitatively, using both experiment and theoretical calculation methods, the effect on water-holding capacity of oxygen-containing functional groups. The experimental results showed that graphite can be oxidized, and forms epoxy groups most easily, followed by hydroxyl and carboxyl groups. The prepared graphene oxide forms a membrane-state as a single layer structure, with an irregular surface. The water-holding capacity of lignite increased with the content of oxygen-containing functional groups. The influence on the configuration of water molecule clusters and binding energy of water molecules of different oxygen-containing functional groups was calculated by density functional theory. The calculation results indicated that the configuration of water molecule clusters was totally changed by oxygen-containing functional groups. The order of binding energy produced by oxygen-containing functional groups and water molecules was as follows: carboxyl > edge phenol hydroxyl >epoxy group. Finally, it can be concluded that the potential to form more hydrogen bonds is the key factor influencing the interaction energy between model compounds and water molecules.

Roles of sliding-induced defects and dissociated water molecules on low friction of graphene

Sliding contact experiments and first-principles calculations were performed to elucidate the roles of structural defects and water dissociative adsorption process on the tribo-chemical mechanisms responsible for low friction of graphene. Sliding friction tests conducted in ambient air and under a dry N2 atmosphere showed that in both cases a high running-in coefficient of friction (COF) occurred initially but a low steady-state COF was reached only when the sliding was continued in air with moisture. Density functional theory (DFT) calculations indicated that the energy barrier (E b ) for dissociative adsorption of H2O was significantly lower in case of reconstructed graphene with a monovacancy compared to pristine graphene. Cross-sectional transmission electron microscopy of graphene transferred to the counterface revealed a partly amorphous structure incorporating damaged graphene layers with d-spacings larger than that of the original layers. DFT calculations on the reconstructed bilayer AB graphene systems revealed an increase of d-spacing due to the chemisorption of H, O, and OH at the vacancy sites and a reduction in the interlayer binding energy (E B ) between the bilayer graphene interfaces compared to pristine graphene. Thus, sliding induced defects facilitated dissociative adsorption of water molecules and reduced COF of graphene for sliding tests under ambient and humid environments but not under an inert atmosphere.

Water Multilayers on TiO2 (101) Anatase Surface: Assessment of a DFTB-Based Method

A water/(101) anatase TiO₂ interface has been investigated with the DFT-based self-consistent-charge density functional tight-binding theory (SCC-DFTB). By comparison of the computed structural, energetic, and dynamical properties with standard DFT-GGA and experimental data, we assess the accuracy of SCC-DFTB for this prototypical solid–liquid interface. We tested different available SCC-DFTB parameters for Ti-containing compounds and, accordingly, combined them to improve the reliability of the method. To better describe water energetics, we have also introduced a modified hydrogen-bond-damping function (HBD). With this correction, equilibrium structures and adsorption energies of water on (101) anatase both for low (0.25 ML) and full (1 ML) coverages are in excellent agreement with those obtained with a higher level of theory (DFT-GGA). Furthermore, Born–Oppenheimer molecular dynamics (MD) simulations for mono-, bi-, and trilayers of water on the surface, as computed with SCC-DFTB, evidence similar ordering and energetics as DFT-GGA Car–Parrinello MD results. Finally, we have evaluated the energy barrier for the dissociation of a water molecule on the anatase (101) surface. Overall, the combined set of parameters with the HBD correction (SCC-DFTB+HBD) is shown to provide a description of the water/water/titania interface, which is very close to that obtained by standard DFT-GGA, with a remarkably reduced computational cost. Hence, this study opens the way to the future investigations on much more extended and realistic TiO₂/liquid water systems, which are extremely relevant for many modern technological applications.

Evaluation of the SCF Combination of KS-DFT and 3D-RISM-KH; Solvation Effect on Conformational Equilibria, Tautomerization Energies, and Activation Barriers

The effect of solvation on conformational equilibria, tautomerization energies, and activation barriers in simple SN2 reactions is reproduced by using the self-consistent field coupling of the Kohn-Sham density functional theory (KS-DFT) for electronic structure and the three-dimensional reference interaction site model with the closure approximation of Kovalenko and Hirata (3D-RISM-KH) for molecular solvation structure. These examples are used in order to validate the implementation of the 3D-RISM-KH method in the Amsterdam Density Functional package. The computations of the free energy difference in the trans/gauche conformational equilibrium for 1,2-dichloroethane in different solvents; the relative tautomerization free energy for cytosine, isocytosine, and guanine; and the free energy activation barrier for a CH3X-type (X = F, Cl, Br) SN2 reaction exhibit agreement with the experimental data. The method is also applied to the electronic and hydration structure of carbon single-wall nanotubes of different diameters, including the effect of water located in the inner space of the nanotubes. A comparison with continuum models of solvation (including COSMO) as well as with other more precise and computationally more expensive calculations is made to demonstrate the accuracy and predictive capability of the new KS-DFT/3D-RISM-KH method.

Water on BN doped benzene: a hard test for exchange-correlation functionals and the impact of exact exchange on weak binding

Density functional theory (DFT) studies of weakly interacting complexes have recently focused on the importance of van der Waals dispersion forces, whereas the role of exchange has received far less attention. Here, by exploiting the subtle binding between water and a boron and nitrogen doped benzene derivative (1,2-azaborine) we show how exact exchange can alter the binding conformation within a complex. Benchmark values have been calculated for three orientations of the water monomer on 1,2-azaborine from explicitly correlated quantum chemical methods, and we have also used diffusion quantum Monte Carlo. For a host of popular DFT exchange-correlation functionals we show that the lack of exact exchange leads to the wrong lowest energy orientation of water on 1,2-azaborine. As such, we suggest that a high proportion of exact exchange and the associated improvement in the electronic structure could be needed for the accurate prediction of physisorption sites on doped surfaces and in complex organic molecules. Meanwhile to predict correct absolute interaction energies an accurate description of exchange needs to be augmented by dispersion inclusive functionals, and certain non-local van der Waals functionals (optB88- and optB86b-vdW) perform very well for absolute interaction energies. Through a comparison with water on benzene and borazine (B3N3H6) we show that these results could have implications for the interaction of water with doped graphene surfaces, and suggest a possible way of tuning the interaction energy.

Anomalous stability of graphene containing defects covered by a water layer

Defects are inevitably present in graphene and can alter its properties and thus its applications. Interestingly, we find that commonly observed Stone-Wales and double vacancy defects do not affect graphene's hydrophilic and hydrophobic properties and that an adsorbed single water layer does not noticeably affect the defect-containing graphene's electronic properties. Our findings are based on calculations using a density functional tight-binding theory. Specifically, we observe negligible alteration in the interaction strength (less than 0.1 kcal mol(-1)) between a single water layer and graphene upon the incorporation of the various types of defects, which indicates that graphene has relatively stable hydrophilic and hydrophobic properties. The presence of a single water layer causes only negligible changes in the energy gap and a small charge transfer to the aqueous layer (less than 0.1 e). The results indicate that the electronic properties of graphene are determined mainly by its own structural characteristics and are not considerably affected by the adsorbed water layer. Further electronic structure analysis reveals that the two commonly observed defects do not change the sp(2) hybridization characteristics of the C atoms of graphene even in the water environment. Our results are significant for graphene studies and applications in areas such as life sciences and materials science where hydrophilic and hydrophobic properties and electronic properties are important.

Probing the Structure of Salt Water under Confinement with First-Principles Molecular Dynamics and Theoretical X-ray Absorption Spectroscopy

We investigated the structure of liquid water around cations (Na(+)) and anions (Cl(-)) confined inside of a (19,0) carbon nanotube with first-principles molecular dynamics and theoretical X-ray absorption spectroscopy (XAS). We found that the ions preferentially reside near the interface between the nanotube and the liquid. Upon confinement, the XAS signal of water molecules surrounding Na(+) exhibits enhanced pre-edge and reduced post-edge features with respect to that of pure water, at variance with the solvation shell of Na(+) in bulk water. Conversely, the first solvation shell of confined Cl(-) has a main-edge intensity comparable to that of bulk solvated Cl(-), likely as a result of a high number of acceptor hydrogen bonds in the first solvation shell. Confined nonsolvating water molecules exhibit bulk-like or water-monomer-like properties, depending on whether they belong to core or interfacial layers, respectively.

Vibrational spectroscopy and molecular dynamics of water monomers and dimers adsorbed on polycyclic aromatic hydrocarbons

This paper reports structures, energetics, dynamics and spectroscopy of H2O and (H2O)2 systems adsorbed on coronene (C24H12), a compact polycyclic aromatic hydrocarbon (PAH). On-the-fly Born-Oppenheimer molecular dynamics simulations are performed for temperatures T varying from 10 to 300 K, on a potential energy surface obtained within the self-consistent-charge density-functional based tight-binding (SCC-DFTB) approach. Anharmonic infrared (IR) spectra are extracted from these simulations. We first benchmark the SCC-DFTB semi-empirical hamiltonian vs. DFT (Density Functional Theory) calculations that include dispersion, on (C6H6)(H2O)1,2 small complexes. We find that charge corrections and inclusion of dispersion contributions in DFTB are necessary to obtain consistent structures, energetics and IR spectra. Using this Hamiltonian, the structures, energetics and IR features of the low-energy isomers of (C24H12)(H2O)1,2 are found to be similar to the DFT ones, with evidence for a stabilizing edge-coordination. The temperature dependence of the motions of H2O and (H2O)2 on the surface of C24H12 is analysed, revealing ultra-fast periodic motion. The water dimer starts diffusing at a higher temperature than the water monomer (150 K vs. 10 K respectively), which appears to be consistent with the binding energies. Qualitative and quantitative analyses of the effects of T on the IR spectra are performed. Anharmonic factors in particular are derived and it is shown that they can be used as signatures for the presence of PAH-water complexes. Finally, this paper lays the foundations for the studies of larger (PAH)m(H2O)n clusters, that can be treated with the efficient computational approach benchmarked in this paper.

Quantum Chemical Prediction of Reaction Pathways and Rate Constants for Dissociative Adsorption of COx and NOx on the Graphite (0001) Surface

We present predictions of reaction rate constants for dissociative adsorption reactions of CO(x) (x = 1, 2) and NO(x) (x = 1, 2) molecules on the basal graphite (0001) surface based on potential energy surfaces (PES) obtained by the integrated ONIOM(B3LYP:DFTB-D) quantum chemical hybrid approach with dispersion-augmented density functional tight binding (DFTB-D) as low level method. Following an a priori methodology developed in a previous investigation of water dissociative adsorption reactions on graphite, we used a C(94)H(24) dicircumcoronene graphene slab as model system for the graphite surface in finite-size molecular structure investigations, and single adsorbate molecules reacting with the pristine graphene sheet. By employing the ONIOM PES information in RRKM theory we predict reaction rate constants in the temperature range between 1,000 and 5,000 K. We find that among CO(x) and NO(x) adsorbate species, the dissociative adsorption reactions of CO(2) and both radical species NO and NO(2) are likely candidates as a cause for high temperature oxidation and erosion of graphite (0001) surfaces, whereas reaction with CO is not likely to lead to long-lived surface defects. High temperature quantum chemical molecular dynamics simulations (QM/MD) at T = 5,000 K using on-the-fly DFTB-D energies and gradients confirm the results of our PES study.

Global Potential Energy Minima of C60(H2O)n Clusters

Likely candidates for the global potential energy minima of C60(H2O)n clusters with n < or = 21 are found using basin-hopping global optimization. The potential energy surfaces are constructed using the TIP4P intermolecular potential for the water molecules, a Lennard-Jones water-fullerene potential, and a water-fullerene polarization potential, which depends on the first few nonvanishing C60 multipole polarizabilities. This combination produces a rather hydrophobic water-fullerene interaction. As a consequence, the water component of the lowest C60(H2O)n minima is quite closely related to low-lying minima of the corresponding TIP4P (H2O)n clusters. In most cases, the geometrical substructure of the water molecules in the C60(H2O)n global minimum coincides with that of the corresponding free water cluster. Exceptions occur when the interaction with C60 induces a change in geometry. This qualitative picture does not change significantly if we use the TIP3P model for the water-water interaction. Structures such as C60@(H2O)60, in which the water molecules surround the C60 fullerene, correspond to local minima with much higher potential energies. For such a structure to become the global minimum, the magnitude of the water-fullerene interaction must be increased to an unphysical value.