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.

Oxygen molecule dissociation on the Al(111) surface

The dissociative adsorption of the O2 on the Al(111) surface is studied using ab initio calculations based on density-functional theory. In the sticking probability experiments an activation barrier for the O2 dissociation exists, but our calculations predict a barrier only for one trajectory. Also our results do not support the model of charge transfer from the surface to the molecule as a bond breaking mechanism. Instead, the increasing hybridization between O2 orbitals and Al states, when the adsorbate approaches the slab, is significant for the dissociation.

Quantum well states in two-dimensional gold clusters on MgO thin films.

The electronic structure of ultrasmall Au clusters on thin MgO/Ag(001) films has been analyzed by scanning tunneling spectroscopy and density functional theory. The clusters exhibit two-dimensional quantum well states, whose shapes resemble the eigenstates of a 2D electron gas confined in a parabolic potential. From the symmetry of the highest occupied (HOMO) and lowest unoccupied molecular orbital (LUMO) of a particular cluster, its electron filling and charge state is determined. In accordance with a Bader charge analysis, aggregates containing up to 20 atoms accumulate one to four extra electrons due to a charge transfer from the MgO/Ag interface. The HOMO-LUMO gap is found to close for clusters containing between 70 and 100 atoms.

Including lateral interactions into microkinetic models of catalytic reactions.

In many catalytic reactions lateral interactions between adsorbates are believed to have a strong influence on the reaction rates. We apply a microkinetic model to explore the effect of lateral interactions and how to efficiently take them into account in a simple catalytic reaction. Three different approximations are investigated: site, mean-field, and quasichemical approximations. The obtained results are compared to accurate Monte Carlo numbers. In the end, we apply the approximations to a real catalytic reaction, namely, ammonia synthesis.

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.

Multidimensional effects on dissociation of N2 on Ru(0001).

The applicability of the Born-Oppenheimer approximation to molecule-metal surface reactions is presently a topic of intense debate. We have performed classical trajectory calculations on a prototype activated dissociation reaction, of N2 on Ru(0001), using a potential energy surface based on density functional theory. The computed reaction probabilities are in good agreement with molecular beam experiments. Comparison to previous calculations shows that the rotation of N2 and its motion along the surface affect the reactivity of N2 much more than nonadiabatic effects.

Ammonia synthesis from first-principles calculations.

The rate of ammonia synthesis over a nanoparticle ruthenium catalyst can be calculated directly on the basis of a quantum chemical treatment of the problem using density functional theory. We compared the results to measured rates over a ruthenium catalyst supported on magnesium aluminum spinel. When the size distribution of ruthenium particles measured by transmission electron microscopy was used as the link between the catalyst material and the theoretical treatment, the calculated rate was within a factor of 3 to 20 of the experimental rate. This offers hope for computer-based methods in the search for catalysts.

Coadsorption of CO and NO on the Pd(111) surface: combined ab initio and Monte Carlo study.

Adsorption and coadsorption of nitric oxide (NO) and carbon monoxide (CO) on the Pd(111) surface are studied by combining first principles (FP) calculations and Monte Carlo (MC) simulations. From FP adsorption energies and molecule-molecule interactions we construct an interaction model, which is used in MC. We do several simulations with different coverages and CO/NO ratios. The simulations provide 0.75 monolayer (ML) for a saturation coverage, which is in excellent agreement with experiments. The results indicate that at 0.75 ML coverage, NO molecules take over the hollow sites and push CO molecules mainly onto bridge sites.

Spectroscopic link between adsorption site occupation and local surface chemical reactivity.

In this Letter we show that sequences of adsorbate-induced shifts of surface core level (SCL) x-ray photoelectron spectra contain profound information on surface changes of electronic structure and reactivity. Energy shifts and intensity changes of time-lapsed spectral components follow simple rules, from which adsorption sites are directly determined. Theoretical calculations rationalize the results for transition metal surfaces in terms of the energy shift of the d-band center of mass and this proves that adsorbate-induced SCL shifts provide a spectroscopic measure of local surface reactivity.