Adsorption, diffusion, and recombination of hydrogen on pure and boron-doped graphite surfaces

Stability of Single Gold Atoms on Defective and Doped Diamond Surfaces

Polycrystalline boron-doped diamond (BDD) is widely used as a working electrode material in electrochemistry, and its properties, such as its stability, make it an appealing support material for nanostructures in electrocatalytic applications. Recent experiments have shown that electrodeposition can lead to the creation of stable small nanoclusters and even single gold adatoms on the BDD surfaces. We investigate the adsorption energy and kinetic stability of single gold atoms adsorbed onto an atomistic model of BDD surfaces by using density functional theory. The surface model is constructed using hybrid quantum mechanics/molecular mechanics embedding techniques and is based on an oxygen-terminated diamond (110) surface. We use the hybrid quantum mechanics/molecular mechanics method to assess the ability of different density functional approximations to predict the adsorption structure, energy, and barrier for diffusion on pristine and defective surfaces. We find that surface defects (vacancies and surface dopants) strongly anchor adatoms on vacancy sites. We further investigated the thermal stability of gold adatoms, which reveals high barriers associated with lateral diffusion away from the vacancy site. The result provides an explanation for the high stability of experimentally imaged single gold adatoms on BDD and a starting point to investigate the early stages of nucleation during metal surface deposition.

Hydrogen Storage in Pure and Boron-Substituted Nanoporous Carbons-Numerical and Experimental Perspective

Nanoporous carbons remain the most promising candidates for effective hydrogen storage by physisorption in currently foreseen hydrogen-based scenarios of the world’s energy future. An optimal sorbent meeting the current technological requirement has not been developed yet. Here we first review the storage limitations of currently available nanoporous carbons, then we discuss possible ways to improve their storage performance. We focus on two fundamental parameters determining the storage (the surface accessible for adsorption and hydrogen adsorption energy). We define numerically the values nanoporous carbons have to show to satisfy mobile application requirements at pressures lower than 120 bar. Possible necessary modifications of the topology and chemical compositions of carbon nanostructures are proposed and discussed. We indicate that pore wall fragmentation (nano-size graphene scaffolds) is a partial solution only, and chemical modifications of the carbon pore walls are required. The positive effects (and their limits) of the carbon substitutions by B and Be atoms are described. The experimental ‘proof of concept’ of the proposed strategies is also presented. We show that boron substituted nanoporous carbons prepared by a simple arc-discharge technique show a hydrogen adsorption energy twice as high as their pure carbon analogs. These preliminary results justify the continuation of the joint experimental and numerical research effort in this field.

Interaction of Hydrogen with MB6 (M = Ba, Ca, La, and Sr) Surfaces from First Principles

We show results of basic energetics and interacting behavior of hydrogen with metal hexaboride surfaces using a combination of self-consistent density functional calculations and dynamics based on the Car-Parrinello method. Our results show that hydrogen is strongly attracted to localized exposed boron atoms and interactions with the terminal cations are strictly repulsive. From these, preliminary local adsorption energy calculations suggest that a single hydrogen molecule per surface unit-cell is possible (one ML). Strongest bonds are found when hydrogen is above the terminal boron atoms affected by reduced coordination and dangling bonds. This location serves to restore the hexaboride unit to a more stable structure by providing electronic density to the deficient surface octahedra. Additionally, trajectories from dynamic simulations provide insight into how hydrogen recombination reactions occur on the surface through dissociative adsorption and the method of travel prior to recombination to be along the octahedral face and bridging sites connecting separate unit cells on the surface. Upon adsorption, a single hydrogen atom becomes localized at the dangling bond site while the second interacts with the surface along a weaker potential energy path. Desorption at lower temperatures occurs when migrating atoms from separate adsorption sites intersect to form a new pair.

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.

Site specificity in femtosecond laser desorption of neutral H atoms from graphite(0001)

Femtosecond laser excitation and density functional theory reveal site and vibrational state specificity in neutral atomic hydrogen desorption from graphite induced by multiple electronic transitions. Multimodal velocity distributions witness the participation of ortho and para pair states of chemisorbed hydrogen in the desorption process. Very slow velocities of 700 and 400  ms^{-1} for H and D atoms are associated with the desorption out of the highest vibrational state of a barrierless potential.

Density Functional Study of Methyl Chemisorption on Polycyclic Aromatic Hydrocarbons

The reactions of methyl radicals with large (up to C(96)H(24)) polycyclic aromatic hydrocarbons (PAHs) are studied by density functional calculations to shed light on the experimentally observed deposition of carbon on highly oriented pyrolytic graphite (HOPG), which occurs when hot HOPG (decorated by nanometre-sized defects) is exposed to methyl radicals. The equilibrium structures of the reaction products, together with transition structures for PAHs up to the size of phenanthroperylene, are determined using the density functionals B3LYP, TPSSh, BP86 and TPSS. The structures are analysed by computing the pi orbital axis vector (POAV) and the altitude of the reactive carbon above the molecular plane of the PAH. The strongest C-CH(3) bonds are found at the edges of the PAHs, where the s character of the C orbital involved in the bond is roughly 25 % (sp(3) hybrid orbital). Carbon atoms inside the PAH form bonds with the methyl radical through atomic orbitals with about 16 % s character in the POAV analysis. These bonds are much weaker than those at the edges of the PAH, while the reactive carbon has moved about 40 pm above the molecular plane. At the edges, the PAH carbon atoms do not leave the molecular plane to this extent. The computed barrier heights and geometrical parameters of the transition structures are in agreement with Hammond's postulate, and the relative energies of all of the equilibrium structures can be rationalized by Hückel molecular orbital (HMO) theory.