Abstraction of D chemisorbed on graphite (0001) with gaseous H atoms

Chemically-resolved determination of hydrogenated graphene-substrate interaction

Functionalization of graphene on Ir(111) is a promising route to modify graphene by chemical means in a controlled fashion at the nanoscale. Yet, the nature of such functionalized sp3 nanodots remains unknown. Density functional theory (DFT) calculations alone cannot differentiate between two plausible structures, namely true graphane and substrate stabilized graphane-like nanodots. These two structures, however, interact dramatically differently with the underlying substrate. Discriminating which type of nanodots forms on the surface is thus of paramount importance for the applications of such prepared nanostructures. By comparing X-ray standing wave measurements against theoretical model structures obtained by DFT calculations we are able to exclude the formation of true graphane nanodots and clearly show the formation graphane-like nanodots.

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.

How the alignment of adsorbed ortho H pairs determines the onset of selective carbon nanotube etching

Unlocking the enormous technological potential of carbon nanotubes strongly depends on our ability to specifically produce metallic or semiconducting tubes. While selective etching of both has already been demonstrated, the underlying reasons, however, remain elusive as yet. We here present computational and experimental evidence on the operative mechanisms at the atomic scale. We demonstrate that during the adsorption of H atoms and their coalescence, the adsorbed ortho hydrogen pairs on single-walled carbon nanotubes induce higher shear stresses than axial stresses, leading to the elongation of HC-CH bonds as a function of their alignment with the tube chirality vector, which we denote as the γ-angle. As a result, the C-C cleavage occurs more rapidly in nanotubes containing ortho H-pairs with a small γ-angle. This phenomenon can explain the selective etching of small-diameter semiconductor nanotubes with a similar curvature. Both theoretical and experimental results strongly indicate the important role of the γ-angle in the selective etching mechanisms of carbon nanotubes, in addition to the nanotube curvature and metallicity effects and lead us to clearly understand the onset of selective synthesis/removal of CNT-based materials.

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.

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.

Controlling hydrogenation of graphene on Ir(111)

Combined fast X-ray photoelectron spectroscopy and density functional theory calculations reveal the presence of two types of hydrogen adsorbate structures at the graphene/Ir(111) interface, namely, graphane-like islands and hydrogen dimer structures. While the former give rise to a periodic pattern, dimers tend to destroy the periodicity. Our data reveal distinctive growth rates and stability of both types of structures, thereby allowing one to obtain well-defined patterns of hydrogen clusters. The ability to control and manipulate the formation and size of hydrogen structures on graphene facilitates tailoring of its properties for a wide range of applications by means of covalent functionalization.

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.

Reduced density matrix quantum approach for particle trapping and sticking on corrugated moving surfaces

A short time propagation algorithm for the reduced density matrix is derived to model the interaction of a quantum particle with a moving corrugated surface. The algorithm includes dissipative terms, which can be derived directly from the full Hamiltonian. The scattering of He from a corrugated Cu surface is examined as a function of incident energy and angle and the temperature of the substrate, with a focus on the nature of trapping. It is found that corrugation can make a significant contribution to trapping, even on a metal surface. Energy exchange with the phonons is shown to significantly modify the nature of diffraction mediated selective adsorption.

Clustering of chemisorbed H(D) atoms on the graphite (0001) surface due to preferential sticking

We present scanning tunneling microscopy experiments and density functional theory calculations which reveal a unique mechanism for the formation of hydrogen adsorbate clusters on graphite surfaces. Our results show that diffusion of hydrogen atoms is largely inactive and that clustering is a consequence of preferential sticking into specific adsorbate structures. These surprising findings are caused by reduced or even vanishing adsorption barriers for hydrogen in the vicinity of already adsorbed H atoms on the surface and point to a possible novel route to interstellar H2 formation.

Classical Studies of H Atom Trapping on a Graphite Surface†

The trapping and sticking of H and D atoms on the graphite (0001) surface is examined over the energy range 0.1-0.9 eV. Total electronic energy calculations based on density functional theory are used to develop a potential energy surface that allows for the full three-dimensional motion of the incident atom and the reconstruction of the bonding carbon atom, which must pucker out of the surface to form a stable bond. Classical methods are used to compute trapping cross sections as a function of incident energy. The C-H bond, once formed, rapidly dissociates without a mechanism to dissipate its excess energy. However, a number of long-lived trapping resonances exist, and for impact parameters below 1 A or so, several percent of the incident H atoms can remain trapped for 1 ps or more. This long-time trapping probability increases significantly when additional lattice degrees of freedom are added to carry energy away from the C-H stretch. Trapping can also increase with an increasing collision impact parameter, as H vibrations parallel to the surface become excited, leaving less energy in the C-H stretch. The trapping cross section at 1 ps reaches a maximum of 0.2 A2 for an H atom energy of 0.3 eV. Assuming that any atoms remaining trapped after 1 ps fully relax and stick, we estimate a lower bound for the sticking probability of H and D to be 0.024 and 0.050, respectively, about an order of magnitude below the experimental values.

Metastable structures and recombination pathways for atomic hydrogen on the graphite (0001) surface

We present scanning tunneling microscopy results which reveal the existence of two distinct hydrogen dimer states on graphite basal planes. Density functional theory calculations allow us to identify the atomic structure of these states and to determine their recombination and desorption pathways. Direct recombination is only possible from one of the two dimer states. This results in increased stability of one dimer species and explains the puzzling double peak structure observed in temperature programmed desorption spectra for hydrogen on graphite.

Quantum study of Eley-Rideal reaction and collision induced desorption of hydrogen atoms on a graphite surface. II. H-physisorbed case

Following previous investigation of collision induced (CI) processes involving hydrogen atoms chemisorbed on graphite [R. Martinazzo and G. F. Tantardini, J. Chem. Phys. 124, 124702 (2006)], the case in which the target hydrogen atom is initially physisorbed on the surface is considered here. Several adsorbate-substrate initial states of the target H atom in the physisorption well are considered, and CI processes are studied for projectile energies up to 1 eV. Results show that (i) Eley-Rideal cross sections at low collision energies may be larger than those found in the H-chemisorbed case but they rapidly decrease as the collision energy increases; (ii) product hydrogen molecules are vibrationally very excited; (iii) collision induced desorption cross sections rapidly increase, reaching saturation values greater than 10 A2; (iv) trapping of the incident atoms is found to be as efficient as the Eley-Rideal reaction at low energies and remains sizable (3-4 A2) at high energies. The latter adsorbate-induced trapping results mainly in formation of metastable hot hydrogen atoms, i.e., atoms with an excess energy channeled in the motion parallel to the surface. These atoms might contribute in explaining hydrogen formation on graphite.

The rovibrational distribution of H2 and HD formed on a graphite surface at 15–50 K

The rotational distributions of H2 and HD formed on a highly oriented pyrolitic graphite surface at temperatures of 15-50 K have been measured using laser spectroscopy. The population of the rovibrational levels nu=1, J=0-4 and nu=2, J=0-4 has been observed and the average rotational temperatures of the nascent H2 and HD molecules have been determined. We find that the average rotational temperature of the newly formed molecules is much higher than the surface temperature on which they have formed. We compare our results with other recent experimental data and theoretical calculations.

Quantum studies of H atom trapping on a graphite surface

The trapping and sticking of H and D atoms on the graphite (0001) surface is examined, over the energy range of 0.1-0.9 eV. For hydrogen to chemisorb onto graphite, the bonding carbon must pucker out of the surface plane by several tenths of an angstrom. A quantum approach in which both the hydrogen and the bonding carbon atoms can move is used to model the trapping, and a potential energy surface based on density functional theory calculations is employed. It is found, for energies not too far above the 0.2 eV barrier to chemisorption that a significant fraction of the incident H or D atoms can trap. The forces on the bonding carbon are large, and it can reconstruct within 50 fs or so. After about 100 fs, most of the trapped H atoms scatter back into the gas phase, but the 5%-10% that remain can have lifetimes on the order of a picosecond or more. Calculations of the resonance eigenstates and lifetimes confirm this. An additional lattice degree of freedom is included quantum mechanically and is shown to significantly increase the amount of H that remains trapped after 1 ps. Further increasing the incident energy destabilizes the trapped state, leading to less H remaining trapped at long times. We estimate that for a full dissipative bath, the sticking probabilities should be on the order of 0.1.