Hydrogen-graphite interaction: Experimental evidences of an adsorption barrier

Exploration of Graphene Defect Reactivity toward a Hydrogen Radical Utilizing a Preactivated Circumcoronene Model

A preexisting chemisorbed defect is well-known to increase the reactivity of graphene which is normally chemically inert. Specifically, the presence of chemisorbed hydrogen atoms forming an sp³-hybridized C-H bond is known to increase the reactivity of neighboring carbon atoms toward additional hydrogenation with wide-ranging applications from materials science to astrochemistry. In this work, static DFT and DFT-based direct dynamics simulations are used to characterize the reactivity of a graphene sheet around an existing C-H bond defect. The spin density landscape shows how to guide subsequent H atom additions, always bonding most strongly to the carbon atom with greatest spin density. Molecular dynamics of an impinging H atom under thermal conditions with defect graphene was used to determine the statistics of probable reactions. The most frequent outcome is inelastic scattering (48%) and then Eley-Rideal (ER) abstraction of the chemisorbed H atom as vibrationally hot H₂ (40%), while the least likely, but probably most interesting, result is formation of a novel C-H bond (12%). The C-H bonds always form in the β sublattice. The carbon atom in the para position shows to be most reactive toward the incoming H atom, followed by the ortho carbon, in agreement with the spin density computed in the static calculations. Globally, the graphene energy surface is repulsive, but the defects create local channels into this energy surface through which reactants can move locally through and react with the activated surface without a barrier.

Dual-Route Hydrogenation of the Graphene/Ni Interface

Nanostructured architectures based on graphene/metal interfaces might be efficiently exploited in hydrogen storage due to the attractive capability to provide adsorption sites both at the top side of graphene and at the metal substrate after intercalation. We combined in situ high-resolution X-ray photoelectron spectroscopy and scanning tunneling microscopy with theoretical calculations to determine the arrangement of hydrogen atoms at the graphene/Ni(111) interface at room temperature. Our results show that at low coverage H atoms predominantly adsorb as monomers and that chemisorption saturates when ∼25% of the surface is hydrogenated. In parallel, with a much lower rate, H atoms intercalate below graphene and bind to Ni surface sites. Intercalation progressively destabilizes the C-H bonds and triggers the release of the hydrogen chemisorbed on graphene. Valence band and near-edge absorption spectroscopy demonstrate that the graphene layer is fully lifted when the Ni surface is saturated with H. Thermal programmed desorption was used to determine the stability of the hydrogenated interface. Whereas the H atoms chemisorbed on graphene remain unperturbed over a wide temperature range, the intercalated phase abruptly desorbs 50-100 K above room temperature.

Sticking of atomic hydrogen on graphene

Recent years have witnessed an ever growing interest in the interactions between hydrogen atoms and a graphene sheet. Largely motivated by the possibility of modulating the electric, optical and magnetic properties of graphene, a huge number of studies have appeared recently that added to and enlarged earlier investigations on graphite and other carbon materials. In this review we give a glimpse of the many facets of this adsorption process, as they emerged from these studies. The focus is on those issues that have been addressed in detail, under carefully controlled conditions, with an emphasis on the interplay between the adatom structures, their formation dynamics and the electric, magnetic and chemical properties of the carbon sheet.

Thermodynamics and kinetics of graphene chemistry: a graphene hydrogenation prototype study

The thermodynamic and kinetic controls of graphene chemistry are studied computationally using a graphene hydrogenation reaction and polyaromatic hydrocarbons to represent the graphene surface. Hydrogen atoms are concertedly chemisorped onto the surface of graphene models of different shapes (i.e., all-zigzag, all-armchair, zigzag-armchair mixed edges) and sizes (i.e., from 16-42 carbon atoms). The second-order Z-averaged perturbation theory (ZAPT2) method combined with Pople double and triple zeta basis sets are used for all calculations. It is found that both the net enthalpy change and the barrier height of graphene hydrogenation at graphene edges are lower than at their interior surfaces. While the thermodynamic product distribution is mainly determined by the remaining π-islands of functionalized graphenes (Phys. Chem. Chem. Phys., 2013, 15, 3725-3735), the kinetics of the reaction is primarily correlated with the localization of the electrostatic potential of the graphene surface.

Quantum dynamical investigation of the isotope effect in H2 formation on graphite at cold collision energies

The Eley-Rideal abstraction of hydrogen atoms on graphitic surfaces at cold collision energies was investigated using a time-dependent wave packet method within the rigid-flat surface approximation, with a focus on hydrogen-deuterium isotopic substitutions. It is found that the marked isotope effect of collinear collisions disappears when the full dimensionality of the problem is taken into account, thereby suggesting that abstraction is less direct than commonly believed and proceeds through glancing rather than head-on collisions. In contrast, a clear isotope effect is observed for "hot-atom" formation, which appears to be strongly favored for heavy projectiles because of their higher density of physisorbed states. Overall, the dynamics is essentially classical and reasonably well described by quasi-classical trajectory methods at all but the lowest energies (≲10 meV). A comparison of the results obtained in the (substrate) adiabatic and diabatic limits suggests that the reaction is only marginally affected by the lattice dynamics, but highlights the importance of including energy dissipation processes in order to accurately describe the internal excitation of the product molecules.

Quantum dynamics of hydrogen atoms on graphene. II. Sticking

Following our recent system-bath modeling of the interaction between a hydrogen atom and a graphene surface [Bonfanti et al., J. Chem. Phys. 143, 124703 (2015)], we present the results of converged quantum scattering calculations on the activated sticking dynamics. The focus of this study is the collinear scattering on a surface at zero temperature, which is treated with high-dimensional wavepacket propagations with the multi-configuration time-dependent Hartree method. At low collision energies, barrier-crossing dominates the sticking and any projectile that overcomes the barrier gets trapped in the chemisorption well. However, at high collision energies, energy transfer to the surface is a limiting factor, and fast H atoms hardly dissipate their excess energy and stick on the surface. As a consequence, the sticking coefficient is maximum (∼0.65) at an energy which is about one and half larger than the barrier height. Comparison of the results with classical and quasi-classical calculations shows that quantum fluctuations of the lattice play a primary role in the dynamics. A simple impulsive model describing the collision of a classical projectile with a quantum surface is developed which reproduces the quantum results remarkably well for all but the lowest energies, thereby capturing the essential physics of the activated sticking dynamics investigated.

Direct observation of ordered configurations of hydrogen adatoms on graphene

Ordered configurations of hydrogen adatoms on graphene have long been proposed, calculated, and searched for. Here, we report direct observation of several ordered configurations of H adatoms on graphene by scanning tunneling microscopy. On the top side of the graphene plane, H atoms in the configurations appear to stick to carbon atoms in the same sublattice. Scanning tunneling spectroscopy measurements revealed a substantial gap in the local density of states in H-contained regions as well as in-gap states below the conduction band due to the incompleteness of H ordering. These findings can be well explained by density functional theory calculations based on double-sided H configurations. In addition, factors that may influence H ordering are discussed.

Cooperative interplay of van der Waals forces and quantum nuclear effects on adsorption: H at graphene and at coronene

The energetic barriers that atoms and molecules often experience when binding to surfaces are incredibly important to a myriad of chemical and physical processes. However, these barriers are difficult to describe accurately with current computer simulation approaches. Two prominent contemporary challenges faced by simulation are the role of van der Waals forces and nuclear quantum effects. Here we examine the widely studied model systems of hydrogen on graphene and coronene using a van der Waals inclusive density functional theory approach together with path integral molecular dynamics at 50 K. We find that both van der Waals and quantum nuclear effects work together in a cooperative manner to dramatically reduce the barriers for hydrogen atoms to adsorb. This suggests that the low temperature hydrogenation of graphene is easier than previously thought and in more general terms that the combined roles of van der Waals and quantum tunnelling can lead to qualitative changes in adsorption.

Diffusion, adsorption, and desorption of molecular hydrogen on graphene and in graphite

The diffusion of molecular hydrogen (H2) on a layer of graphene and in the interlayer space between the layers of graphite is studied using molecular dynamics computer simulations. The interatomic interactions were modeled by an Adaptive Intermolecular Reactive Empirical Bond Order (AIREBO) potential. Molecular statics calculations of H2 on graphene indicate binding energies ranging from 41 meV to 54 meV and migration barriers ranging from 3 meV to 12 meV. The potential energy surface of an H2 molecule on graphene, with the full relaxations of molecular hydrogen and carbon atoms is calculated. Barriers for the formation of H2 through the Langmuir-Hinshelwood mechanism are calculated. Molecular dynamics calculations of mean square displacements and average surface lifetimes of H2 on graphene at various temperatures indicate a diffusion barrier of 9.8 meV and a desorption barrier of 28.7 meV. Similar calculations for the diffusion of H2 in the interlayer space between the graphite sheets indicate high and low temperature regimes for the diffusion with barriers of 51.2 meV and 11.5 meV. Our results are compared with those of first principles.

Insights into H2 formation in space from ab initio molecular dynamics

Hydrogen formation is a key process for the physics and the chemistry of interstellar clouds. Molecular hydrogen is believed to form on the carbonaceous surface of dust grains, and several mechanisms have been invoked to explain its abundance in different regions of space, from cold interstellar clouds to warm photon-dominated regions. Here, we investigate direct (Eley-Rideal) recombination including lattice dynamics, surface corrugation, and competing H-dimers formation by means of ab initio molecular dynamics. We find that Eley-Rideal reaction dominates at energies relevant for the interstellar medium and alone may explain observations if the possibility of facile sticking at special sites (edges, point defects, etc.) on the surface of the dust grains is taken into account.

Reaction of the Basal Plane of Graphite with the Methyl Radical

The reaction of methyl radicals with the basal plane of graphite has been observed to occur with an activation energy of less than 0.3 eV. This reaction is initiated by Li-induced CH3Cl dissociation to produce CH3 radicals on the graphite surface. It is found that ∼3/4 of the methyl radicals remain on the graphite surface up to 700 K at puckered sp(3) carbon sites, while 1/4 of the CH3 radicals participate in CH4 formation and small amounts of C2 and C3 hydrocarbon formation. CH3 radicals become mobile over an activation energy barrier of ∼0.7 eV.

On the PES for the interaction of an H atom with an H chemisorbate on a graphenic platelet

Motivated by the problem of H(2) formation in diffuse clouds of the interstellar medium (ISM), we study the effect of including van der Waals-type corrections in DFT calculations on the entrance PES of the Eley-Rideal reaction H(b) + H(a)-GR → H(b)-H(a) + GR for a graphenic surface GR. The present calculations make use of the PBE-D3 dispersion corrected functional of Grimme et al. (2010) and are carried out on cluster models of graphenic surfaces: C(24)H(12) and C(54)H(18). To assess the soundness of the chosen functional we start by revisiting the H-GR adsorption potential. We find a satisfactory on top physisorption well (43-48 meV) correctly located at an H-GR distance of 3 Å. We then revisit the H(b)-H(a)-GR system using both the PW91 and PBE functionals. Our calculations do not reproduce the tiny potential barrier reported earlier for large H(b)distances from the surface. The barrier in the calculations of Sidis et al. (2000) and Morisset et al. (2003, 2004) has been traced to their previous use of an LSDA + POSTSCF PW91 procedure rather than the genuine PW91 one. The new PBE-D3 PES for the H(b)-H(a)-GR system is reported as a function of the H(b) distance to the surface and its impact parameter relative to the H(a) chemisorbate for the so-called "fixed puckered" ("diabatic" or "sudden") approach. The results are discussed in relation to recent experimental and theoretical work.