Dynamics of high-barrier surface reactions: laser-assisted associative desorption of N2 from Ru(0001)

First-Principles Study of Stability and N2 Activation on the Octahedron RuRh Clusters

The geometric and electronic structures of different octahedron RuRh clusters are studied using density functional theory calculations. The binding energy, electronic structure, and energy gap of the clusters have been obtained to determine the possible stable structures. The results show that the Ru4Rh2 cluster is the most stable structure which has D4h symmetry with the largest ionization potential, smallest affinity energy and larger energy gap. Furthermore, the information on adsorption and dissociation of multiple nitrogen molecules and the density of state for the octahedral Ru4Rh2 cluster is analyzed. The dissociation barrier of three nitrogen molecules further decreases to 1.18 eV with an increase in the number of N2 molecules. The co-adsorption of multiple N2 molecules facilitates the dissociation of N2 on the Ru4Rh2 cluster. The strong interaction between the antibonding orbital of N2 and the d orbital of the Ru4Rh2 cluster is illustrated by calculating and analyzing the results of PDOS, which stretches the N−N bond length and reduces the activation energy to dissociation. The antibonding orbital of the nitrogen molecule shows distinct and unique catalytic activity for the dissociation of the adsorbed nitrogen molecule on the octahedral Ru4Rh2 cluster.

Acid-durable electride with layered ruthenium for ammonia synthesis: boosting the activity via selective etching

We reported LnRuSi as a B₅-site-free Ru catalyst for ammonia synthesis, and its activity enhanced 2–4-fold by selective etching with EDTA-2Na.

Ruthenium (Ru) loaded catalysts are of significant interest for ammonia synthesis under mild reaction conditions. The B₅ sites have been reported as the active sites for ammonia formation, i.e., Ru with other coordinations were inactive, which has limited the utilization efficiency of Ru metal. The implantation of Ru into intermetallic compounds is considered to be a promising approach to tune the catalytic activity and utilization efficiency of Ru. Here we report an acid-durable electride, LnRuSi (Ln = La, Ce, Pr and Nd), as a B₅-site-free Ru catalyst. The active Ru plane with a negative charge is selectively exposed by chemical etching using disodium dihydrogen ethylenediaminetetraacetate (EDTA-2Na) acid, which leads to 2–4-fold enhancement in the ammonia formation rate compared with that of the original catalyst. The turnover frequency (TOF) of LnRuSi is estimated to be approximately 0.06 s^–1, which is 600 times higher than that of pure Ru powder. Density functional theory (DFT) calculations revealed that the dissociation of N₂ occurs easily on the exposed Ru plane of LaRuSi. This systematic study provides firm evidence that layered Ru with a negative charge in LnRuSi is a new type of active site that differs significantly from B₅ sites.

Mechanism and kinetics for both thermal and electrochemical reduction of N2 catalysed by Ru(0001) based on quantum mechanics

QM calculations were used to predict the free energy surfaces for N₂ thermal and electrochemical reduction (N₂TR and N₂ER) on Ru(0001), to find the detailed atomistic mechanism and kinetics, and provide the basis for improving the efficiency of N₂ER.

Analysis of Energy Dissipation Channels in a Benchmark System of Activated Dissociation: N2 on Ru(0001)

The excitation of electron-hole pairs in reactive scattering of molecules at metal surfaces often affects the physical and dynamical observables of interest, including the reaction probability. Here, we study the influence of electron-hole pair excitation on the dissociative chemisorption of N2 on Ru(0001) using the local density friction approximation method. The effect of surface atom motion has also been taken into account by a high-dimensional neural network potential. Our nonadiabatic molecular dynamics simulations with electronic friction show that the reaction of N2 is more strongly affected by the energy transfer to surface phonons than by the energy loss to electron-hole pairs. The discrepancy between the computed reaction probabilities and experimental results is within the experimental error both with and without friction; however, the incorporation of electron-hole pairs yields somewhat better agreement with experiments, especially at high collision energies. We also calculate the vibrational efficacy for the N2 + Ru(0001) reaction and demonstrate that the N2 reaction is more enhanced by exciting the molecular vibrations than by adding an equivalent amount of energy into translation.

Accurate Neural Network Description of Surface Phonons in Reactive Gas-Surface Dynamics: N2 + Ru(0001)

Ab initio molecular dynamics (AIMD) simulations enable the accurate description of reactive molecule-surface scattering especially if energy transfer involving surface phonons is important. However, presently, the computational expense of AIMD rules out its application to systems where reaction probabilities are smaller than about 1%. Here we show that this problem can be overcome by a high-dimensional neural network fit of the molecule-surface interaction potential, which also incorporates the dependence on phonons by taking into account all degrees of freedom of the surface explicitly. As shown for N₂ + Ru(0001), which is a prototypical case for highly activated dissociative chemisorption, the method allows an accurate description of the coupling of molecular and surface atom motion and accurately accounts for vibrational properties of the employed slab model of Ru(0001). The neural network potential allows reaction probabilities as low as 10^-5 to be computed, showing good agreement with experimental results.

Observation of Tunneling in the Hydrogenation of Atomic Nitrogen on the Ru(001) Surface to Form NH

The kinetics of NH and ND formation and dissociation reactions on Ru(001) were studied using time-dependent reflection absorption infrared spectroscopy (RAIRS). Our results indicate that NH and ND formation and dissociation on Ru(001) follow first-order kinetics. In our reaction temperature range (320-390 K for NH and 340-390 K for ND), the apparent activation energies for NH and ND formation were found to be 72.2 ± 1.9 and 87.1 ± 1.8 kJ/mol, respectively, while NH and ND dissociation reactions between 370 and 400 K have apparent activation barriers of 106.9 ± 4.1 and 101.8 ± 4.8 kJ/mol, respectively. The lower apparent activation energy for NH formation than that for ND as well as the comparison between experimentally measured isotope effects with theoretical results strongly indicates that tunneling already starts to play a role in this reaction at a temperature as high as 340 K.

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.

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.

First-principles study of some factors controlling the rate of ammonia decomposition on Ni and Pd surfaces

Using the plane-wave pseudopotential method within the density-functional theory with the generalized gradient approximation for exchange and correlation potential, we have calculated adsorption energies (E(ad)), diffusion barrier, and the first dissociation barrier (E(1)) for NH(3) on Ni and Pd surfaces. While the top site is found to be preferred for NH(3) adsorption on both Ni(111) and Pd(111), its calculated diffusion barrier is substantially higher for Pd(111) than for Ni(111). We also find that during the first dissociation step (NH(3)-->NH(2)+H), NH(2) moves from the top site to the nearest hollow site on Ni(111) and Pd(111) and on the stepped surfaces, Ni(211) and Pd(211), it moves from the initial top site at the step edge to the bridge site in the same atomic chain. Meanwhile H is found to occupy the hollow sites on all four surfaces. On Ni(111), E(1) is found to be 0.23 eV higher than E(ad), while at the step of Ni(211), E(1) and E(ad) are almost equal, suggesting that the probability for the molecule to dissociate is much on the step of Ni(211). In the case of Pd(211), however, we find that the dissociation barrier is much higher than E(ad). These trends are in qualitative agreement with the experimental finding that ammonia decomposition rate is much lower on Pd than on Ni.