CHD3 dissociation on Pt(111): A comparison of the reaction dynamics based on the PBE functional and on a specific reaction parameter functional.

We present a comparison of ab initio molecular dynamics calculations for CHD3 dissociation on Pt(111) using the Perdew, Burke and Ernzerhof (PBE) functional and a specific reaction parameter (SRP) functional. Despite the two functionals predicting approximately the same activation barrier for the reaction, the calculations using the PBE functional consistently overestimate the experimentally determined dissociation probability, whereas the SRP functional reproduces the experimental values within a chemical accuracy (4.2 kJ/mol). We present evidence that suggests that this difference in reactivity can at least in part be attributed to the presence of a van der Waals well in the potential of the SRP functional which is absent from the PBE description. This leads to the CHD3 molecules being accelerated and spending less time near the surface for the trajectories run with the SRP functional, as well as more energy being transferred to the surface atoms. We suggest that both these factors reduce the reactivity observed in the SRP calculations compared to the PBE calculations.

Exploring surface landscapes with molecules: rotationally induced diffraction of H2 on LiF(001) under fast grazing incidence conditions.

Atomic diffraction by surfaces under fast grazing incidence conditions has been used for almost a decade to characterize surface properties with more accuracy than with more traditional atomic diffraction methods. From six-dimensional solutions of the time-dependent Schrödinger equation, we show that diffraction of H2 molecules under fast grazing incidence conditions could be even more informative for the characterization of ionic surfaces, due to the large anisotropic electrostatic interaction between the quadrupole moment of the molecule and the electric field created by the ionic crystal. Using the LiF(001) surface as a benchmark, we show that fast grazing incidence diffraction of H2 strongly depends on the initial rotational state of the molecule, while rotationally inelastic processes are irrelevant. We demonstrate that, as a result of the anisotropy of the impinging projectile, initial rotational excitation leads to an increase in intensity of high-order diffraction peaks at incidence directions that satisfy precise symmetry constraints, thus providing a more detailed information on the surface characteristics than that obtained from low-order atomic diffraction peaks under fast grazing incidence conditions. As quadrupole-ion surface potentials are expected to accurately represent the interaction between H2 and any surface with a marked ionic character, our analysis should be of general applicability to any of such surfaces. Finally, we show that a density functional theory description of the molecule-ion surface potential catches the main features observed experimentally.

Possible effect of static surface disorder on diffractive scattering of H2 from Ru(0001): Comparison between theory and experiment.

Specific features of diffractive scattering of H2 from metal surfaces can serve as fingerprints of the reactivity of the metal towards H2, and in principle theory-experiment comparisons for molecular diffraction can help with the validation of semi-empirical functionals fitted to experiments of sticking of H2 on metals. However, a recent comparison of calculated and Debye-Waller (DW) extrapolated experimental diffraction probabilities, in which the theory was done on the basis of a potential energy surface (PES) accurately describing sticking to Ru(0001), showed substantial discrepancies, with theoretical and experimental probabilities differing by factors of 2 and 3. We demonstrate that assuming a particular amount of random static disorder to be present in the positions of the surface atoms, which can be characterized through a single parameter, removes most of the discrepancies between experiment and theory. Further improvement might be achievable by improving the accuracy of the DW extrapolation, the model of the H2 rotational state distribution in the experimental beams, and by fine-tuning the PES. However, the question of whether the DW model is applicable to attenuation of diffractive scattering in the presence of a sizable van der Waals well (depth ≈ 50 meV) should also receive attention, in addition to the question of whether the amount of static surface disorder effectively assumed in the modeling by us could have been present in the experiments.

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

In dynamics calculations of H2 dissociating on metal surfaces often clean, high-symmetry surfaces are used. Few such dynamics studies have been performed on surfaces with pre-adsorbed molecules, especially when also the motion of the surface and the adsorbate are considered. In this study, the dissociation of H2 on a carbon monoxide-covered Ru(0001) surface is considered. Ab initio molecular dynamics (AIMD) calculations are performed on this system using the PBE-vdW-DF2 functional, which accurately describes the reaction probability for H2 on Ru(0001). Using this functional, the reaction probability of H2 on the CO-covered Ru(0001) surface is found to be too low when compared to experiments. This suggests that exchange-correlation functionals that can describe the reaction of H2 on a bare metal surface are not in general able to describe the reaction of H2 on a CO-precovered surface of the same metal, with the same accuracy. However, it cannot be ruled out that the discrepancy between theory and experiment is partly due to an inhomogeneous coverage of the surface by CO in the experiments. The incorporation of the motion of the surface has only a small effect on the reaction probability. It is found that when including surface motion for this system, the size of the simulation cell can be important. Upon collision, a considerable amount of energy is transferred to the surface, causing the adsorbed CO molecules to move apart, which opens the surface for reaction. In order to obtain converged reaction probabilities with respect to the size of the simulation cell, at least a 3 × 3 simulation cell is needed, because in the smaller cell the CO molecules cannot be pushed apart as only a single independent CO molecule is present, also leading to less energy exchange with the surface.

Electronic friction dominates hydrogen hot-atom relaxation on Pd(100).

We study the dynamics of transient hot H atoms on Pd(100) that originated from dissociative adsorption of H2. The methodology developed here, denoted AIMDEF, consists of ab initio molecular dynamics simulations that include a friction force to account for the energy transfer to the electronic system. We find that the excitation of electron-hole pairs is the main channel for energy dissipation, which happens at a rate that is five times faster than energy transfer into Pd lattice motion. Our results show that electronic excitations may constitute the dominant dissipation channel in the relaxation of hot atoms on surfaces.

The effect of the exchange-correlation functional on H2 dissociation on Ru(0001).

The specific reaction parameter (SRP) approach to density functional theory (DFT) has enabled a chemically accurate description of reactive scattering experiments for activated H2-metal systems (H2 + Cu(111) and Cu(100)), but its application has not yet resulted in a similarly accurate description of non-activated or weakly activated H2-metal systems. In this study, the effect of the choice of the exchange-correlation functional in DFT on the potential energy surface and dynamics of H2 dissociation on Ru(0001), a weakly activated system, is investigated. In total, full potential energy surfaces were calculated for over 20 different functionals. The functionals investigated include functionals incorporating an approximate description of the van der Waals dispersion in the correlation functional (vdW-DF and vdW-DF2 functionals), as well as the revTPSS meta-GGA. With two of the functionals investigated here, which include vdW-DF and vdW-DF2 correlation, it has been possible to accurately reproduce molecular beam experiments on sticking of H2 and D2, as these functionals yield a reaction probability curve with an appropriate energy width. Diffraction probabilities computed with these two functionals are however too high compared to experimental diffraction probabilities, which are extrapolated from surface temperatures (Ts) ⩾ 500 K to 0 K using a Debye-Waller model. Further research is needed to establish whether this constitutes a failure of the two candidate SRP functionals or a failure of the Debye-Waller model, the use of which can perhaps in future be avoided by performing calculations that include the effect of surface atom displacement or motion, and thereby of the experimental Ts.

Molecular dynamics simulations of CO2 formation in interstellar ices.

CO2 ice is one of the most abundant components in ice-coated interstellar ices besides H2O and CO, but the most favorable path to CO2 ice is still unclear. Molecular dynamics calculations on the ultraviolet photodissociation of different kinds of CO-H2O ice systems have been performed at 10 K in order to demonstrate that the reaction between CO and an OH molecule resulting from H2O photodissociation through the first excited state is a possible route to form CO2 ice. However, our calculations, which take into account different ice surface models, suggest that there is another product with a higher formation probability ((3.00 ± 0.07) × 10(-2)), which is the HOCO complex, whereas the formation of CO2 has a probability of only (3.6 ± 0.7) × 10(-4). The initial location of the CO is key to obtain reaction and form CO2: the CO needs to be located deep into the ice. The HOCO complex becomes trapped in the cold ice surface in the trans-HOCO minimum because it quickly loses its internal energy to the surrounding ice, preventing further reaction to H + CO2. Several laboratory experiments have been carried out recently, and they confirm that CO2 can also be formed through other, different routes. Here we compare our theoretical results with the data available from experiments studying the formation of CO2 through a similar pathway as ours, even though the initial conditions were not exactly the same. Our results also show that the HCO van der Waals complex can be formed through the interaction of CO with the H atom that is formed as a product of H2O photodissociation. Thus, the reaction of the H atom photofragment following H2O photodissociation with CO can be a possible route to form HCO ice.

Isotope effects on the photodesorption processes of X2O (X = H,D) and HOD ice.

To investigate the isotope effects on the photodesorption processes of X2O (X = H,D) ice, molecular dynamics calculations have been performed on the ultraviolet photodissociation of an H2O or a D2O molecule in an H2O or a D2O amorphous ice surface, and on HOD photodissociation in an H2O amorphous ice surface, where the photodissociated molecules were located in the top four or five monolayers at ice temperatures of 10, 20, 30, 60, and 90 K. Three photodesorption processes can occur upon X2O photodissociation: X atom photodesorption, OX radical photodesorption, and X2O (or HOD) molecule photodesorption. X2O (or HOD) photodesorption can occur after recombination of X and OX, or after an energetic X atom photofragment kicks a surrounding X2O molecule from the ice surface. Isotope effects are observed for the X atom and the OX radical photodesorption as well as for the kick-out photodesorption. However, no isotope effects were noticeable for the photodesorption of recombined X2O molecules. The average D atom photodesorption probabilities are about a factor 0.9 smaller than those for the H atom, regardless of the isotope of the surrounding ice system. Also, the kick-out mechanism is more likely to occur if a D photofragment is created upon dissociation than if an H atom is created. These observations can be explained by more efficient energy transfer from the D atom to water molecules than from the H atom. Reasoning based on the X2O phonon frequencies associated with the librational modes and energy transfer efficiencies explain why the OX radical photodesorption probabilities are noticeably larger if the OX radical desorbs from a D2O ice system than from an H2O ice system. Also, the OX radical photodesorption is more probable upon dissociation of DOX (X = H,D) than upon dissociation of HOX (X = H,D), because the initial kinetic energy of the OX radical is larger if the dissociation products are D + OX than H + OX. The branching ratio of OD/OH desorption following photodissociation of an HOD molecule in ice (about 1.0) is much lower than the OD/OH branching ratio in gas-phase HOD photodissociation. This may lead to differences in isotope fractionation in OH(g) formation in dense and diffuse clouds in the interstellar medium.

Towards a specific reaction parameter density functional for reactive scattering of H2 from Pd(111).

Recently, an implementation of the specific reaction parameter (SRP) approach to density functional theory (DFT) was used to study several reactive scattering experiments of H2 on Cu(111). It was possible to obtain chemical accuracy (1 kcal/mol ≈ 4.2 kJ/mol), and therefore, accurately model this paradigmatic example of activated H2 dissociation on a metal surface. In this work, the SRP-DFT methodology is applied to the dissociation of hydrogen on a Pd(111) surface, in order to test whether the SRP-DFT approach is also applicable to non-activated H2-metal systems. In the calculations, the Born-Oppenheimer static surface approximations are used. A comparison to molecular beam sticking experiments, performed at incidence energies ≥125 meV, on H2 + Pd(111) suggested the PBE-vdW [where the Perdew, Burke, and Ernzerhof (PBE) correlation is replaced by van der Waals correlation] functional as a candidate SRP density functional describing the reactive scattering of H2 on Pd(111). Unfortunately, quantum dynamics calculations are not able to reproduce the molecular beam sticking results for incidence energies <125 meV. From a comparison to initial state-resolved (degeneracy averaged) sticking probabilities it seems clear that for H2 + Pd(111) dynamic trapping and steering effects are important, and that these effects are not yet well modeled with the potential energy surfaces considered here. Applying the SRP-DFT method to systems where H2 dissociation is non-activated remains difficult. It is suggested that a density functional that yields a broader barrier distribution and has more non-activated pathways than PBE-vdW (i.e., non-activated dissociation at some sites but similarly high barriers at the high energy end of the spectrum) should allow a more accurate description of the available experiments. Finally, it is suggested that new and better characterized molecular beam sticking experiments be done on H2 + Pd(111), to facilitate the development of a more accurate theoretical description of this system.

Reactive scattering of H2 from Cu(100): comparison of dynamics calculations based on the specific reaction parameter approach to density functional theory with experiment.

We present new experimental and theoretical results for reactive scattering of dihydrogen from Cu(100). In the new experiments, the associative desorption of H(2) is studied in a velocity resolved and final rovibrational state selected manner, using time-of-flight techniques in combination with resonance-enhanced multi-photon ionization laser detection. Average desorption energies and rotational quadrupole alignment parameters were obtained in this way for a number of (v = 0, 1) rotational states, v being the vibrational quantum number. Results of quantum dynamics calculations based on a potential energy surface computed with a specific reaction parameter (SRP) density functional, which was derived earlier for dihydrogen interacting with Cu(111), are compared with the results of the new experiments and with the results of previous molecular beam experiments on sticking of H(2) and on rovibrationally elastic and inelastic scattering of H(2) and D(2) from Cu(100). The calculations use the Born-Oppenheimer and static surface approximations. With the functional derived semi-empirically for dihydrogen + Cu(111), a chemically accurate description is obtained of the molecular beam experiments on sticking of H(2) on Cu(100), and a highly accurate description is obtained of rovibrationally elastic and inelastic scattering of D(2) from Cu(100) and of the orientational dependence of the reaction of (v = 1, j = 2 - 4) H(2) on Cu(100). This suggests that a SRP density functional derived for H(2) interacting with a specific low index face of a metal will yield accurate results for H(2) reactively scattering from another low index face of the same metal, and that it may also yield accurate results for H(2) interacting with a defected (e.g., stepped) surface of that same metal, in a system of catalytic interest. However, the description that was obtained of the average desorption energies, of rovibrationally elastic and inelastic scattering of H(2) from Cu(100), and of the orientational dependence of reaction of (v = 0, j = 3 - 5, 8) H(2) on Cu(100) compares less well with the available experiments. More research is needed to establish whether more accurate SRP-density functional theory dynamics results can be obtained for these observables if surface atom motion is added to the dynamical model. The experimentally and theoretically found dependence of the rotational quadrupole alignment parameter on the rotational quantum number provides evidence for rotational enhancement of reaction at low translational energies.

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.

Vibrational deexcitation and rotational excitation of H2 and D2 scattered from Cu(111): adiabatic versus non-adiabatic dynamics.

We have studied survival and rotational excitation probabilities of H(2)(v(i) = 1, J(i) = 1) and D(2)(v(i) = 1, J(i) = 2) upon scattering from Cu(111) using six-dimensional (6D) adiabatic (quantum and quasi-classical) and non-adiabatic (quasi-classical) dynamics. Non-adiabatic dynamics, based on a friction model, has been used to analyze the role of electron-hole pair excitations. Comparison between adiabatic and non-adiabatic calculations reveals a smaller influence of non-adiabatic effects on the energy dependence of the vibrational deexcitation mechanism than previously suggested by low-dimensional dynamics calculations. Specifically, we show that 6D adiabatic dynamics can account for the increase of vibrational deexcitation as a function of the incidence energy, as well as for the isotope effect observed experimentally in the energy dependence for H(2)(D(2))/Cu(100). Furthermore, a detailed analysis, based on classical trajectories, reveals that in trajectories leading to vibrational deexcitation, the minimum classical turning point is close to the top site, reflecting the multidimensionally of this mechanism. On this site, the reaction path curvature favors vibrational inelastic scattering. Finally, we show that the probability for a molecule to get close to the top site is higher for H(2) than for D(2), which explains the isotope effect found experimentally.

Diffractive and reactive scattering of H2 from Ru(0001): experimental and theoretical study.

We present a combined experimental and theoretical study of the diffraction of H(2) from Ru(0001) in the incident energy range 78-150 meV, and a theoretical study of dissociative chemisorption of H(2) in the same system. Pronounced out-of-plane diffraction was observed in the whole energy range studied. The energy dependence of the elastic diffraction intensities was measured along the two main symmetry directions for a fixed parallel translational energy. The data were compared with quantum dynamics calculations performed by using DFT-based, six-dimensional potential energy surfaces calculated with both the PW91 and RPBE functionals, as well as with a functional obtained from a weighted average of both (the MIX functional, which was earlier shown to perform quite well for H(2) + Cu(111)). Our results show that the PW91 functional describes the H(2) diffraction intensities more accurately than the RPBE and the MIX functionals, although the absolute values of these intensities are overestimated in the calculations. For the reaction probabilities a preference for one or the other functional cannot be given over the entire energy range probed by the sticking experiments. The PW91 functional yields too high reaction probabilities over the entire investigated energy range, but is better than RPBE at low collision energies (<0.1 eV). The RPBE functional gives too low reaction probabilities at low energy and somewhat too high reaction probabilities at high energy, but agrees better with experiment than PW91 for energies >0.1 eV. The results suggest that, in order to get a better description of both H(2) diffraction and dissociative chemisorption for this system, a specific reaction parameter functional for H(2) + Ru(0001) is needed that is a weighted average of functionals other than PW91 and RPBE. We speculate that differences between the H(2) + Ru(0001) system (early and low reaction barrier) and H(2) + Cu(111) (late and high reaction barrier) may well lead to fundamentally different specific reaction parameter functionals, and that including a reasonable accurate description of the van der Waals interaction might be important for H(2) + Ru(0001) which has barriers localised far away from the surface. Based on our results we advocate new, systematic combined theoretical and experimental studies of H(2) interacting with transition metals in early and late barrier systems, with the aim of determining whether specific reaction parameter functionals for these systems might differ in a systematic way.

Molecular dynamics simulations of D2O ice photodesorption.

Molecular dynamics (MD) calculations have been performed to study the ultraviolet (UV) photodissociation of D(2)O in an amorphous D(2)O ice surface at 10, 20, 60, and 90 K, in order to investigate the influence of isotope effects on the photodesorption processes. As for H(2)O, the main processes after UV photodissociation are trapping and desorption of either fragments or D(2)O molecules. Trapping mainly takes place in the deeper monolayers of the ice, whereas desorption occurs in the uppermost layers. There are three desorption processes: D atom, OD radical, and D(2)O molecule photodesorption. D(2)O desorption takes places either by direct desorption of a recombined D(2)O molecule, or when an energetic D atom produced by photodissociation kicks a surrounding D(2)O molecule out of the surface by transferring part of its momentum. Desorption probabilities are calculated for photoexcitation of D(2)O in the top four monolayers and are compared quantitatively with those for H(2)O obtained from previous MD simulations of UV photodissociation of amorphous water ice at different ice temperatures [Arasa et al., J. Chem. Phys. 132, 184510 (2010)]. The main conclusions are the same, but the average D atom photodesorption probability is smaller than that of the H atom (by about a factor of 0.9) because D has lower kinetic energy than H, whereas the average OD radical photodesorption probability is larger than that of OH (by about a factor of 2.5-2.9 depending on ice temperature) because OD has higher translational energy than OH for every ice temperature studied. The average D(2)O photodesorption probability is larger than that of H(2)O (by about a factor of 1.4-2.3 depending on ice temperature), and this is entirely due to a larger contribution of the D(2)O kick-out mechanism. This is an isotope effect: the kick-out mechanism is more efficient for D(2)O ice, because the D atom formed after D(2)O photodissociation has a larger momentum than photogenerated H atoms from H(2)O, and D transfers momentum more easily to D(2)O than H to H(2)O. The total (OD + D(2)O) yield has been compared with experiments and the total (OH + H(2)O) yield from previous simulations. We find better agreement when we compare experimental yields with calculated yields for D(2)O ice than when we compare with calculated yields for H(2)O ice.

Six-dimensional dynamics study of reactive and non reactive scattering of H(2) from Cu(111) using a chemically accurate potential energy surface.

We have studied the interaction of H(2) on Cu(111) using quasi-classical and quantum dynamics, and a chemically accurate six-dimensional potential energy surface (PES). The PES was computed using the specific reaction parameter (SRP) approach to density functional theory (DFT), in an implementation adapted to molecules interacting with metal surfaces. To perform this study we have applied the Born-Oppenheimer static surface (BOSS) approximation, i.e., we used both the Born-Oppenheimer (BO) and the static surface (SS) approximations. We show that our theoretical approach accurately describes experiments on dissociative adsorption, the effect of molecular vibrational and rotational motion on dissociative (associative) adsorption (desorption), and rotational excitation upon scattering. More specifically, dynamics calculations on reactive scattering of H(2) reproduce reaction probabilities measured in molecular beam experiments, effective barrier heights describing the dependence of reaction on the initial rovibrational state, and data on rotationally inelastic scattering with chemical accuracy (i.e., within 1 kcal mol(-1) approximately 4.2 kJ mol(-1)). These processes are not affected much by surface motion, either because they were measured using a low surface temperature, T(s), or because the computed observable is independent of T(s). However, we show that to account for the dependence of molecular orientation on a reaction the inclusion of surface motion is required. We have also found that vibrational excitation is poorly described within the BOSS approximation, suggesting a breakdown of this approximation.

Hydrogen dissociation on small aluminum clusters.

Transition states and reaction paths for a hydrogen molecule dissociating on small aluminum clusters have been calculated using density functional theory. The two lowest spin states have been taken into account for all the Al(n) clusters considered, with n=2-6. The aluminum dimer, which shows a (3)Π(u) electronic ground state, has also been studied at the coupled cluster and configuration interaction level for comparison and to check the accuracy of single determinant calculations in this special case, where two degenerate configurations should be taken into account. The calculated reaction barriers give an explanation of the experimentally observed reactivity of hydrogen on Al clusters of different size [Cox et al., J. Chem. Phys. 84, 4651 (1986)] and reproduce the high observed reactivity of the Al(6) cluster. The electronic structure of the Al(n)-H(2) systems was also systematically investigated in order to determine the role played by interactions of specific molecular orbitals for different nuclear arrangements. Singlet Al(n) clusters (with n even) exhibit the lowest barriers to H(2) dissociation because their highest doubly occupied molecular orbitals allow for a more favorable interaction with the antibonding σ(u) molecular orbital of H(2).

A theoretical study of H(2) dissociation on (sq.rt(3) x sq.rt(3))R30 degrees CO/Ru(0001).

We have studied the influence of preadsorbed CO on the dissociative adsorption of H(2) on Ru(0001) with density functional theory calculations. For a coverage of 1/3 ML CO, we investigated different possible reaction paths for hydrogen dissociation using nudged elastic band and adaptive nudged elastic band calculations. One reaction path was studied in detail through an energy decomposition and molecular orbital type of analysis. The minimum barrier for H(2) dissociation is found to be 0.29 eV. At the barrier the H-H bond is hardly stretched. Behind this barrier a molecular chemisorption minimum is present. Next, the molecule overcomes a second barrier, with a second local chemisorption minimum behind it. To finally dissociate to chemisorbed atoms, the molecule has to overcome a third barrier. To move along the reaction path from reactants to products, the hydrogen molecule needs to rotate, and to significantly change its center-of-mass position. The procedure of mapping out reaction paths for H(2) reacting on low-index surfaces of bare metals (computing two-dimensional elbow plots for fixed impact high-symmetry sites and H(2) orientations parallel to the surface) does not work for H(2)+CO/Ru. The first barrier in the path is recovered, but the features of the subsequent stretch to the dissociative chemisorption minimum are not captured, because the molecule is not allowed to change its center-of-mass position or to rotate. The dissociative chemisorption of H(2) on CO/Ru(0001) is endoergic, in contrast to the case of H(2) on bare Ru(0001). The zero-point energy corrected energies of molecularly and dissociatively chemisorbed H(2) are very close, suggesting that it may be possible to detect molecularly chemisorbed H(2) on (sq.rt(3) x sq.rt(3))R30 degrees CO/Ru(0001). The presence of CO on the surface increases the barrier height to dissociation compared with bare Ru(0001). Based on an energy decomposition and molecular orbital analysis we attribute the increase in the barrier height mainly to an occupied-occupied interaction between the bonding H(2) sigma(g) orbital and the (surface-hybridized) CO 1pi orbitals, i.e., to site blocking. There is a small repulsive contribution to the barrier from the interaction between the H(2) molecule and the Ru part of the CO covered Ru surface, but it is smaller than one might expect based on the calculations of H(2) interacting with a clean Ru surface, and on calculations of H(2) interacting with the CO overlayer only. Actually, the analysis suggests that the Ru surface as a subsystem is (slightly) more reactive for the reaction path studied with CO preadsorbed on it than without it. Thus, the results indicate that the influence of CO on H(2) dissociation on Ru is not only a simple site-blocking effect, the electronic structure of the underlying Ru is changed.

Quantum dynamics of dissociative chemisorption of CH(4) on Ni(111): Influence of the bending vibration.

Two-dimensional, three-dimensional, and four-dimensional quantum dynamic calculations are performed on the dissociative chemisorption of CH(4) on Ni(111) using the multiconfiguration time-dependent Hartree (MCTDH) method. The potential energy surface used for these calculations is 15-dimensional (15D) and was obtained with density functional theory for points which are concentrated in the region that is dynamically relevant to reaction. Many reduced dimensionality calculations were already performed on this system, but the molecule was generally treated as pseudodiatomic. The main improvement of our model is that we try to describe CH(4) as a polyatomic molecule by including a degree of freedom describing a bending vibration in our three-dimensional and four-dimensional models. Using a polyspherical coordinate system, a general expression for the 15D kinetic energy operator is derived, which discards all the singularities in the operator and includes rotational and Coriolis coupling. We use seven rigid constraints to fix the CH(3) umbrella of the molecule to its gas phase equilibrium geometry and to derive two-dimensional, three-dimensional, and four-dimensional Hamiltonians, which were used in the MCTDH method. Only four degrees of freedom evolve strongly along the 15D minimum energy path: the distance of the center of mass of the molecule to the surface, the dissociative C[Single Bond]H bond distance, the polar orientation of the molecule, and the bending angle between the dissociative C[Single Bond]H bond and the umbrella. A selection of these coordinates is included in each of our models. The polar rotation is found to be important in determining the mode selective behavior of the reaction. Furthermore, our calculations are in good agreement with the finding of Xiang et al. [J. Chem. Phys. 117, 7698 (2002)] in their reduced dimensional calculation that the helicopter motion of the umbrella symmetry axis is less efficient than its cartwheel motion for promoting the reaction. The effect of pre-exciting the bend modes is qualitatively incorrect at higher energies, suggesting the necessity of including additional rotational and vibrational degrees of freedom in the model.

First principles study of the photo-oxidation of water on tungsten trioxide (WO3).

The photo-oxidation of water on the monoclinic P2(1)/nWO(3) (200, 020, and 002) surfaces is investigated using density functional theory calculations, employing the PW91-generalized gradient approximation, and the method developed by Norskov et al. [J. Phys. Chem. B 108, 17886 (2004)] based on the free energy differences between the reaction intermediates. We first relax the bulk material unit cell and then investigate the relative stability of different surface terminations of WO(3) and analyze the overpotential needed for the photoelectrolysis of water. We found that the rate limiting step is the transfer of a proton from the surface adsorbed OH to the electrolyte, and that the computed overpotential for O(2) evolution (1.04 V) is available upon illumination of the surface with visible light.

Predicting catalysis: understanding ammonia synthesis from first-principles calculations.

Here, we give a full account of a large collaborative effort toward an atomic-scale understanding of modern industrial ammonia production over ruthenium catalysts. We show that overall rates of ammonia production can be determined by applying various levels of theory (including transition state theory with or without tunneling corrections, and quantum dynamics) to a range of relevant elementary reaction steps, such as N(2) dissociation, H(2) dissociation, and hydrogenation of the intermediate reactants. A complete kinetic model based on the most relevant elementary steps can be established for any given point along an industrial reactor, and the kinetic results can be integrated over the catalyst bed to determine the industrial reactor yield. We find that, given the present uncertainties, the rate of ammonia production is well-determined directly from our atomic-scale calculations. Furthermore, our studies provide new insight into several related fields, for instance, gas-phase and electrochemical ammonia synthesis. The success of predicting the outcome of a catalytic reaction from first-principles calculations supports our point of view that, in the future, theory will be a fully integrated tool in the search for the next generation of catalysts.