Kinetic Monte Carlo studies of hydrogen abstraction from graphite

First-Principles Study of the Reaction Mechanism of CHO + H on Graphene Surface

Many organic molecules observed in the interstellar medium are considered to be formed on dust grains and populated into the gas phase. We analyzed the reaction of HCO + H on a graphene surface using ab initio molecular dynamics simulations as a case study of the formation and desorption of organic molecules on interstellar dust particles. During the reactions of chemisorbed CHO (chemisorbed at the C atom) with free H, CO was generated and efficiently desorbed from the surface. These results suggest that the reactions, of which the reactant forms a covalent bond with the surface while the product does not, cause efficient desorption of the product upon reaction. In such reactions a repulsive force between the product and the surface would be generated and accelerate translation of the product in a specific direction. In addition, it was also shown that the branching ratio of the reactions between radical species on the surface would be affected by the form of the adsorption on the surface, e.g., when a free H reacted with the CHO chemisorbed at the C atom, CH₂O was not generated.

Ortho and para hydrogen dimers on G/SiC(0001): combined STM and DFT study

The hydrogen (H) dimer structures formed upon room-temperature H adsorption on single layer graphene (SLG) grown on SiC(0001) are addressed using a combined theoretical-experimental approach. Our study includes density functional theory (DFT) calculations for the full (6√3 × 6√3)R30° unit cell of the SLG/SiC(0001) substrate and atomically resolved scanning tunneling microscopy images determining simultaneously the graphene lattice and the internal structure of the H adsorbates. We show that H atoms normally group in chemisorbed coupled structures of different sizes and orientations. We make an atomic scale determination of the most stable experimental geometries, the small dimers and ellipsoid-shaped features, and we assign them to hydrogen adsorbed in para dimers and ortho dimers configuration, respectively, through comparison with the theory.

Nanoconfined NaAlH4: prolific effects from increased surface area and pore volume

Nanoconfinement is a promising technique to improve the properties of nanomaterials such as the kinetics for hydrogen release and uptake and the stability during cycling. Here we present a systematic study of nanoconfined NaAlH4 in nanoporous scaffolds with increasing surface area and pore volume and almost constant pore sizes in the range of 8 to 11 nm. A resorcinol formaldehyde carbon aerogel was CO2-activated under different conditions and provided aerogels with BET surface areas of 704, 1267 and 2246 m(2) g(-1) and total pore volumes of 0.91, 1.30 and 2.21 mL g(-1), respectively. Nanoconfinement of NaAlH4 was achieved by melt infiltration and (27)Al MAS NMR reveals that the respective scaffolds incorporate 68, 82 and 91 wt% NaAlH4, for the above-mentioned samples, while the remaining fraction decomposes to metallic Al indicating that increasing CO2-activation tends to facilitate the infiltration process. The frequencies for the (23)Na and (27)Al MAS NMR centerband resonances from NaAlH4 vary systematically for the infiltrated samples and are shifted towards higher frequency and become more narrow with increasing degree of CO2 activation of the scaffolds. This new effect is attributed to increasing interactions with conduction electrons from increasingly graphite-/graphene-like scaffolds. The bulk versus nanoconfined ratio of NaAlH4 was investigated using Rietveld refinement, revealing that the majority of added NaAlH4 is confined inside the nanopores. The hydrogen desorption kinetics decreased with increasing surface area and the hydrogen storage capacity is more stable and decreases less during continuous hydrogen release and uptake cycles. In fact, the available amount of hydrogen (2.7 wt% H2) was more than doubled compared to the nanoconfinement in the non-activated carbon aerogel (1.3 wt% H2). Furthermore, it was demonstrated that Ti-functionalization of the CO2-activated aerogels combines the high storage capacity with fast hydrogen release kinetics from NaAlH4 which fully decomposes into Na3AlH6 at T ≤ 100 °C.

The formation of ice mantles on interstellar grains revisited – the effect of exothermicity

Modelling of grain surface chemistry generally deals with the simulation of rare events. Usually deterministic methods or statistical approaches such as the kinetic Monte Carlo technique are applied for these simulations. All assume that the surface processes are memoryless, the Markov chain assumption, and usually also that their rates are time independent. In this paper we investigate surface reactions for which these assumptions are not valid, and discuss what the effect is on the formation of water on interstellar grains. We will particularly focus on the formation of two OH radicals by the reaction H + HO₂. Two reaction products are formed in this exothermic reaction and the resulting momentum gained causes them to move away from each other. What makes this reaction special is that the two products can undergo a follow-up reaction to form H₂O₂. Experimentally, OH has been observed, which means that the follow-up reaction does not proceed with 100% efficiency, even though the two OH radicals are formed in each other's vicinity in the same reaction. This can be explained by a combined effect of the directionality of the OH radical movement together with energy dissipation. Both effects are constrained by comparison with experiments, and the resulting parametrised mechanism is applied to simulations of the formation of water ice under interstellar conditions.

Polycyclic aromatic hydrocarbons – catalysts for molecular hydrogen formation

Polycyclic aromatic hydrocarbons (PAHs) have been shown to catalyse molecular hydrogen formation. The process occurs via atomic hydrogen addition reactions leading to the formation of super-hydrogenated PAH species, followed by molecular hydrogen forming abstraction reactions. Here, we combine quadrupole mass spectrometry data with kinetic simulations to follow the addition of deuterium atoms to the PAH molecule coronene. When exposed to sufficiently large D atom fluences, coronene is observed to be driven towards the completely deuterated state (C₂₄D₃₆) with the mass distribution peaking at 358 amu, just below the peak mass of 360 amu. Kinetic models reproduce the experimental observations for an abstraction cross-section of σ_abs = 0.01 Ų per excess H/D atom, and addition cross-sections in the range of σ_add = 0.55–2.0 Ų for all degrees of hydrogenation. These findings indicate that the cross-section for addition does not scale with the number of sites available for addition on the molecule, but rather has a fairly constant value over a large interval of super-hydrogenation levels.

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.

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.

Quantum dynamics of hydrogen interacting with single-walled carbon nanotubes: multiple H-atom adsorbates

In a previous paper [J. L. McAfee and B. Poirier, J. Chem. Phys. 130, 064701 (2009)], using spin-polarized density functional theory (DFT), the authors reported a binding energy of 0.755 eV, for a single hydrogen atom adsorbed on a pristine (unrelaxed) (5,5) single-walled carbon nanotube (SWNT) substrate. A full three-dimensional (3D) potential energy surface (PES) for the SWNT-H system was also developed, and used in a quantum dynamics calculation to compute all rovibrational bound states, and associated equatorial and longitudinal adsorbate migration rates. A highly pronounced preference for the latter migration pathway at ambient temperatures was observed. In this work, we extend the aforementioned study to include multiple H-atom adsorbates. Extensive DFT calculations are performed, in order to ascertain the most relevant dynamical pathways. For two adsorbates, the SWNT-H-H system is found to exhibit highly site-specific binding, as well as long-range correlation and pronounced binding energy enhancement. The latter effect is even more pronounced in the full-hydrogenation limit, increasing the per-adsorbate binding energy to 2.6 eV. To study migration dynamics, a single-hole model is developed, for which the binding energy drops to 2.11 eV. A global 3D PES is developed for the hole migration model, using 40 radial × 18 cylindrical ab initio geometries, fit to a Fourier basis with radially dependent expansion coefficients (rms error 4.9 meV). As compared with the single-adsorbate case, the hole migration PES does not exhibit separate chemisorption and physisorption wells. The barrier to longitudinal migration is also found to be much lower. Quantum dynamics calculations for all rovibrational states are then performed (using a mixed spectral basis/phase-space optimized discrete variable representation), and used to compute longitudinal migration rates. Ramifications for the use of SWNTs as potential hydrogen storage materials are discussed.