Quantum Dynamics of the Eley−Rideal Hydrogen Formation Reaction on Graphite at Typical Interstellar Cloud Conditions

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.

Full quantum dynamical investigation of the Eley-Rideal reaction forming H2 on a movable graphitic substrate at T = 0 K

The dynamics of the Eley-Rideal abstraction reaction of hydrogen atoms on a movable graphitic surface is investigated for the first time in a numerically exact fully quantum setting. A system-bath strategy was applied where the two recombining H atoms and a substrate C atom form a relevant subsystem, while the rest of the lattice takes the form of an independent oscillator bath. High-dimensional wavepacket simulations were performed in the collision energy range 0.2-1.0 eV with the help of the multi-layer multi-configuration time-dependent Hartree method, focusing on the collinear reaction on a zero-temperature surface. Results show that the dynamics is close to a sudden limit in which the reaction is much faster than the substrate motion. Unpuckering of the surface is fast (some tens of fs) but starts only after the formation of H₂ is completed, thereby determining a considerable substrate heating (∼0.8 eV per reactive event). Energy partitioning in the product molecule favors translational over vibrational energy, and H₂ molecules are vibrationally hot (∼1.5 eV) though to a lesser extent than previously predicted.

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.

Etching of graphene in a Hydrogen-rich Atmosphere towards the Formation of Hydrocarbons in Circumstellar Clouds

We describe a mechanism that explains the formation of hydrocarbons and hydrocarbyls from hydrogenated graphene/graphite; hard C-C bonds are weakened and broken by the synergistic effect of chemisorbed hydrogen and high temperature vibrations. Total energies, optimized structures, and transition states are obtained from Density Functional Theory simulations. These values have been used to determine the Boltzman probability for a thermal fluctuation to overcome the kinetic barriers, yielding the time scale for an event to occur. This mechanism can be used to rationalize the possible routes for the creation of small hydrocarbons and hydrocarbyls from etched graphene/graphite in stellar regions.

Adiabatic potential energy surfaces for the low-energy collisional dynamics of C(+)((2)P) ions with H2 molecules

The low-energy electronic states of the CH2(+) molecular ion are investigated with multireference configuration interaction calculations based on complete active space self-consistent field reference wave functions using a large C(6s5p4d3f)/H(8s6p3d1f) basis set. The focus is on the three lowest-lying states describing formation and destruction of the astrophysically relevant methylidine cation CH(+). Both processes are discussed in light of the topology of the relevant potential energy surfaces and their intersections.

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.

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.