Laser assisted associative desorption of N2 and CO from Ru(0001)

SBH17: Benchmark Database of Barrier Heights for Dissociative Chemisorption on Transition Metal Surfaces

Accurate barriers for rate controlling elementary reactions on metal surfaces are key to understanding, controlling, and predicting the rate of heterogeneously catalyzed processes. While barrier heights for gas phase reactions have been extensively benchmarked, dissociative chemisorption barriers for the reactions of molecules on metal surfaces have received much less attention. The first database called SBH10 and containing 10 entries was recently constructed based on the specific reaction parameter approach to density functional theory (SRP-DFT) and experimental results. We have now constructed a new and improved database (SBH17) containing 17 entries based on SRP-DFT and experiments. For this new SBH17 benchmark study, we have tested three algorithms (high, medium, and light) for calculating barrier heights for dissociative chemisorption on metals, which we have named for the amount of computational effort involved in their use. We test the performance of 14 density functionals at the GGA, GGA+vdW-DF, and meta-GGA rungs. Our results show that, in contrast with the previous SBH10 study where the BEEF-vdW-DF2 functional seemed to be most accurate, the workhorse functional PBE and the MS2 density functional are the most accurate of the GGA and meta-GGA functionals tested. Of the GGA+vdW functionals tested, the SRP32-vdW-DF1 functional is the most accurate. Additionally, we found that the medium algorithm is accurate enough for assessing the performance of the density functionals tested, while it avoids geometry optimizations of minimum barrier geometries for each density functional tested. The medium algorithm does require metal lattice constants and interlayer distances that are optimized separately for each functional. While these are avoided in the light algorithm, this algorithm is found not to give a reliable description of functional performance. The combination of relative ease of use and demonstrated reliability of the medium algorithm will likely pave the way for incorporation of the SBH17 database in larger databases used for testing new density functionals and electronic structure methods.

STM Visualization of N2 Dissociative Chemisorption on Ru(0001) at High Impinging Kinetic Energies

This paper examines the reactive surface dynamics of energy- and angle-selected N₂ dissociation on a clean Ru(0001) surface. Presented herein are the first STM images of highly energetic N₂ dissociation on terrace sites utilizing a novel UHV instrument that combines a supersonic molecular beam with an in situ STM that is in-line with the molecular beam. Atomically resolved visualization of individual N₂ dissociation events elucidates the fundamental reactive dynamics of the N₂/Ru(0001) system by providing a detailed understanding of the on-surface dissociation dynamics: the distance and angle between nitrogen atoms from the same dissociated N₂ molecule, site specificity and coordination of binding on terrace sites, and the local evolution of surrounding nanoscopic areas. These properties are precisely measured over a range of impinging N₂ kinetic energies and angles, revealing previously unattainable information about the energy dissipation channels that govern the reactivity of the system. The experimental results presented in this paper provide insight into the fundamental N₂ dissociation mechanism that, in conjunction with ongoing theoretical modeling, will help determine the role of dynamical processes such as energy transfer to surface phonons and nonadiabatic excitation of electron-hole pairs (ehps). These results will not only help uncover the underlying chemistry and physics that give rise to the unique behavior of this activated dissociative chemisorption system but also represent an exciting approach to studying reaction dynamics by pairing the angstrom-level spatiotemporal resolution of an in situ STM with nonequilibrium fluxes of reactive gases generated in a supersonic molecular beam to access highly activated chemical dynamics and observe the results of individual reaction events.

Neural Network Sampling of the Free Energy Landscape for Nitrogen Dissociation on Ruthenium

In heterogeneous catalysis, free energy profiles of reactions govern the mechanisms, rates, and equilibria. Energetics are conventionally computed using the harmonic approximation (HA), which requires determination of critical states a priori. Here, we use neural networks to efficiently sample and directly calculate the free energy surface (FES) of a prototypical heterogeneous catalysis reaction-the dissociation of molecular nitrogen on ruthenium-at density-functional-theory-level accuracy. We find that the vibrational entropy of surface atoms, often neglected in HA for transition metal catalysts, contributes significantly to the reaction barrier. The minimum free energy path for dissociation reveals an "on-top" adsorbed molecular state prior to the transition state. While a previously reported flat-lying molecular metastable state can be identified in the potential energy surface, it is absent in the FES at relevant reaction temperatures. These findings demonstrate the importance of identifying critical points self-consistently on the FES for reactions that involve considerable entropic effects.

Orbital-Dependent Electronic Friction Significantly Affects the Description of Reactive Scattering of N2 from Ru(0001)

Electron-hole pair (ehp) excitation is thought to substantially affect the dynamics of molecules on metal surfaces, but it is not clear whether this can be better addressed by orbital-dependent friction (ODF) or the local density friction approximation (LDFA). We investigate the effect of ehp excitation on the dissociative chemisorption of N2 on and its inelastic scattering from Ru(0001), which is the benchmark system of highly activated dissociation, with these two different models. ODF is in better agreement with the best experimental estimates for the reaction probabilities than LDFA, yields results for vibrational excitation in better agreement with experiment, but slightly overestimates the translational energy loss during scattering. N2 on Ru(0001) is thus the first system for which the ODF and LDFA approaches are shown to yield substantially different results for easily accessible experimental observables, including reaction probabilities.

An ab initio molecular dynamics study of D2 dissociation on CO-precovered Ru(0001)

In dynamics studies of hydrogen dissociation on CO pre-covered Ru(0001) the simulation cell size is important for describing energy exchange.

Challenges in reduction of dinitrogen by proton and electron transfer

Catalytic reduction of dinitrogen with protons and electrons is a very challenging alternative to the energy expensive Haber–Bosch reaction.

Reactive and nonreactive scattering of N2 from Ru(0001): a six-dimensional adiabatic study

We have studied the dissociative chemisorption and scattering of N(2) on and from Ru(0001), using a six-dimensional quasiclassical trajectory method. The potential energy surface, which depends on all the molecular degrees of freedom, has been built applying a modified Shepard interpolation method to a data set of results from density functional theory, employing the RPBE generalized gradient approximation. The frozen surface and Born-Oppenheimer [Ann. Phys. (Leipzig) 84, 457 (1927)] approximations were used, neglecting phonons and electron-hole pair excitations. Dissociative chemisorption probabilities are found to be very small even for translational energies much higher than the minimum reaction barrier, in good agreement with experiment. A comparison to previous low dimensional calculations shows the importance of taking into account the multidimensional effects of N(2) rotation and translation parallel to the surface. The new calculations strongly suggest a much smaller role of nonadiabatic effects than previously assumed on the basis of a comparison between low dimensional results and experiments [J. Chem. Phys. 115, 9028 (2001)]. Also in agreement with experiment, our theoretical results show a strong dependence of reaction on the initial vibrational state. Computed angular scattering distributions and parallel translation energy distributions are in good agreement with experiments on scattering, but the theory overestimates vibrational and rotational excitations in scattering.

High translational energy release in H2 (D2) associative desorption from H (D) chemisorbed on C(0001)

Highly energetic translational energy distributions are reported for hydrogen and deuterium molecules desorbing associatively from the atomic chemisorption states on highly oriented pyrolytic graphite (HOPG). Laser assisted associative desorption is used to measure the time of flight of molecules desorbing from a hydrogen (deuterium) saturated HOPG surface produced by atomic exposure from a thermal atom source at around 2100 K. The translational energy distributions normal to the surface are very broad, from approximately 0.5 to approximately 3 eV, with a peak at approximately 1.3 eV. The highest translational energy measured is close to the theoretically predicted barrier height. The angular distribution of the desorbing molecules is sharply peaked along the surface normal and is consistent with thermal broadening contributing to energy release parallel to the surface. All results are in qualitative agreement with recent density functional theory calculations suggesting a lowest energy para-type dimer recombination path.

Microcanonical unimolecular rate theory at surfaces. III. Thermal dissociative chemisorption of methane on Pt(111) and detailed balance

A local hot spot model of gas-surface reactivity is used to investigate the state-resolved dynamics of methane dissociative chemisorption on Pt(111) under thermal equilibrium conditions. Three Pt surface oscillators, and the molecular vibrations, rotations, and the translational energy directed along the surface normal are treated as active degrees of freedom in the 16-dimensional microcanonical kinetics. Several energy transfer models for coupling a local hot spot to the surrounding substrate are developed and evaluated within the context of a master equation kinetics approach. Bounds on the thermal dissociative sticking coefficient based on limiting energy transfer models are derived. The three-parameter physisorbed complex microcanonical unimolecular rate theory (PC-MURT) is shown to closely approximate the thermal sticking under any realistic energy transfer model. Assuming an apparent threshold energy for CH(4) dissociative chemisorption of E(0)=0.61 eV on clean Pt(111), the PC-MURT is used to predict angle-resolved yield, translational, vibrational, and rotational distributions for the reactive methane flux at thermal equilibrium at 500 K. By detailed balance, these same distributions should be observed for the methane product from methyl radical hydrogenation at 500 K in the zero coverage limit if the methyl radicals are not subject to side reactions. Given that methyl radical hydrogenation can only be experimentally observed when the CH(3) radicals are kinetically stabilized against decomposition by coadsorbed H, the PC-MURT was used to evaluate E(0) in the high coverage limit. A high coverage value of E(0)=2.3 eV adequately reproduced the experimentally observed methane angular and translational energy distributions from thermal hydrogenation of methyl radicals. Although rigorous application of detailed balance arguments to this reactive system cannot be made because thermal decomposition of the methyl radicals competes with hydrogenation, approximate applicability of detailed balance would argue for a strong coverage dependence of E(0) with H coverage--a dependence not seen for methyl radical hydrogenation on Ru(0001), but not yet experimentally explored on Pt(111).

The reaction rate for dissociative adsorption of N2 on stepped Ru(0001): Six-dimensional quantum calculations

Quantum-mechanical calculations of the reaction rate for dissociative adsorption of N2 on stepped Ru(0001) are presented. Converged six-dimensional quantum calculations for this heavy-atom reaction have been performed using the multiconfiguration time-dependent Hartree method. A potential-energy surface for the transition-state region is constructed from density-functional theory calculations using Shepard interpolation. The quantum results are in very good agreement with the results of the harmonic transition-state theory. In contrast to the findings of previous model calculations on similar systems, the tunneling effect is found to be small.

Influence of surface morphology on D2 desorption kinetics from amorphous solid water

The influence of surface morphology/porosity on the desorption kinetics of weakly bound species was investigated by depositing D2 on amorphous solid water (ASW) films grown by low temperature vapor deposition under various conditions and with differing thermal histories. A broad distribution of binding energies of the D2 monolayer on nonporous and porous ASW was measured experimentally and correlated by theoretical calculations to differences in the degree of coordination of the adsorbed H2 (D2) to H2O molecules in the ASW depending on the nature of the adsorption site, i.e., surface valleys vs surface peaks in a nanoscale rough film surface. For porous films, the effect of porosity on the desorption kinetics was observed to be a reduction in the desorption rate with film thickness and a change in peak shape. This can be partly explained by fast diffusion into the ASW pore structure via a simple one-dimensional diffusion model and by a change in binding energy statistics with increasing total effective surface area. Furthermore, the D2 desorption kinetics on thermally annealed ASW films were investigated. The main effect was seen to be a reduction in porosity and in the number of highly coordinated binding sites with anneal temperature due to ASW restructuring and pore collapse. These results contribute to the understanding of desorption from porous materials and to the development of correct models for desorption from and catalytic processes on dust grain surfaces in the interstellar medium.

Importance of Surface Morphology in Interstellar H 2 Formation

Detailed laboratory experiments on the formation of HD from atom recombination on amorphous solid water films show that this process is extremely efficient in a temperature range of 8 to 20 kelvin, temperatures relevant for H2 formation on dust grain surfaces in the interstellar medium (ISM). The fate of the 4.5 electron volt recombination energy is highly dependent on film morphology. These results suggest that grain morphology, rather than the detailed chemical nature of the grain surface, is most important in determining the energy content of the H2 as it is released from the grain into the ISM.