P. Larrégaray

On the statistical behavior of the O+OH→H+O2 reaction: A comparison between quasiclassical trajectory, quantum scattering, and statistical calculations

The dynamics of the O + OH reaction on the ground state potential energy surface (PES) is investigated by means of the quasiclassical trajectory method and two statistical methods: phase space theory and statistical quantum method. Preliminary calculations with an exact quantum method are also reported. The quasiclassical trajectory calculations show evidence for a phase space bottleneck inhibiting the intramolecular energy transfer between the O-H and O-O bonds. As a result, the probability of the intermediate complex dissociating back toward the reactants is high, thereby yielding a reaction probability significantly lower than expected for a barrierless and exothermic reaction. The features of the PES, which are the cause of this dynamical effect, are identified. This is essentially the conservation of the equilibrium distance of the O-H bond, hardly changed by a close encounter with an oxygen atom. The statistical calculations, which do not take into account the PES in the complex region, yield a high reaction probability, much larger than the probability calculated from the dynamical methods, both classical and quantum. If the statistical cross sections are corrected by a scaling factor, which corresponds actually to scaling the capture probability, then a good agreement is observed between dynamical and statistical calculations of the product state distributions. The differential cross sections calculated with all the methods show a backward-forward symmetry, with sharp polarization peaks. The complex lifetime is divided into two parts by the bottleneck. During the first part, the system remains trapped in a small region of the phase space and has a high probability to dissociate back toward the reactants. This is a nonstatistical effect due to the PES shape. During the second part, fast intramolecular vibrational energy redistribution takes place, leading to a statistical distribution of energy on the rovibrational states of the products. These findings indicate that the O + OH reaction has mixed dynamics, both with statistical and nonstatistical aspects.

Dynamics simulation of N2 scattering onto W(100,110) surfaces: A stringent test for the recently developed flexible periodic London-Eyring-Polanyi-Sato potential energy surface

An efficient method to construct the six dimensional global potential energy surface (PES) for two atoms interacting with a periodic rigid surface, the flexible periodic London-Eyring-Polanyi-Sato model, has been proposed recently. The main advantages of this model, compared to state-of-the-art interpolated ab initio PESs developed in the past, reside in its global nature along with the small number of electronic structure calculations required for its construction. In this work, we investigate to which extent this global representation is able to reproduce the fine details of the scattering dynamics of N(2) onto W(100,110) surfaces reported in previous dynamics simulations based on locally interpolated PESs. The N(2)/W(100) and N(2)/W(110) systems are chosen as benchmarks as they exhibit very unusual and distinct dissociative adsorption dynamics although chemically similar. The reaction pathways as well as the role of dynamic trapping are scrutinized. Besides, elastic/inelastic scattering dynamics including internal state and angular distributions of reflected molecules are also investigated. The results are shown to be in fair agreement with previous theoretical predictions.

Classical treatment of molecular collisions: striking improvement of the description of recoil energy distributions using Gaussian weighted trajectories

The Gaussian weighting (GW) procedure, recently used in the classical treatment of molecular collisions, is a practical way of taking into account quantization of product vibrational actions. The goal of this brief communication is to show that the GW procedure may drastically improve the predictions of the recoil energy distribution between final fragments, an observable frequently measured in molecular beam experiments.

Mean potential phase space theory of chemical reactions.

A nonconventional application of phase space theory to the insertion reactions A+H(2), with A=C((1)D) and S((1)D), is presented. Instead of approximating the potential energies of interaction between separated fragments by their isotropic long-range contributions, as in the original theory, the latter are replaced by the accurate potential energies averaged with respect to Jacobi angles. The integral and differential cross sections obtained from this mean potential phase space theory (MPPST) turn out to be in very satisfying agreement with the benchmark predictions of the time-independent and time-dependent statistical quantum methods. The formal and numerical simplicity of MPPST with respect to any approach combining statistical assumptions and dynamical calculations makes it a promising tool for studying indirect polyatomic reactions.

On product state distributions in triatomic unimolecular reactions: Beyond phase space theory and the adiabatic assumption

Our goal is to derive a simple dynamically corrected statistical treatment of state distributions in the products of triatomic unimolecular reactions involving efficient energy transfers between rotation and translation motions en route to products. For, phase space theory or the adiabatic channel model—the only statistical approaches of final state distributions—may not be applicable to such processes. We thus analyze and model how the departing atom perturbs the rotation of the diatomic molecule. Since the general problem is intricate, we limit our study to the basic, yet realistic and instructive case where the bending force is harmonic and the total angular momentum is zero. A remarkable fact is the one-to-one relation between the perturbation and a coefficient χ, related in a simple way to the mechanical parameters of the system. Transition state theory combined with our model—the linear transformation model—leads to the desired treatment of product state distributions. Its predictions are in good ag...

Surface temperature effects on the dynamics of N2 Eley-Rideal recombination on W(100)

Quasiclassical trajectories simulations are performed to study the influence of surface temperature on the dynamics of a N atom colliding a N-preadsorbed W(100) surface under normal incidence. A generalized Langevin surface oscillator scheme is used to allow energy transfer between the nitrogen atoms and the surface. The influence of the surface temperature on the N(2) formed molecules via Eley-Rideal recombination is analyzed at T = 300, 800, and 1500 K. Ro-vibrational distributions of the N(2) molecules are only slightly affected by the presence of the thermal bath whereas kinetic energy is rather strongly decreased when going from a static surface model to a moving surface one. In terms of reactivity, the moving surface model leads to an increase of atomic trapping cross section yielding to an increase of the so-called hot atoms population and a decrease of the direct Eley-Rideal cross section. The energy exchange between the surface and the nitrogen atoms is semi-quantitatively interpreted by a simple binary collision model.

Dynamical reaction pathways in Eley-Rideal recombination of nitrogen from W(100)

The scattering of atomic nitrogen over a N-pre-adsorbed W(100) surface is theoretically described in the case of normal incidence off a single adsorbate. Dynamical reaction mechanisms, in particular Eley-Rideal (ER) abstraction, are scrutinized in the 0.1-3.0 eV collision energy range and the influence of temperature on reactivity is considered between 300  and 1500 K. Dynamics simulations suggest that, though non-activated reaction pathways exist, the abstraction process exhibits a significant collision energy threshold (0.5 eV). Such a feature, which has not been reported so far in the literature, is the consequence of a repulsive interaction between the impinging and the pre-adsorbed nitrogens along with a strong attraction towards the tungsten atoms. Above threshold, the cross section for ER reaction is found one order of magnitude lower than the one for hot-atoms formation. The abstraction process involves the collision of the impinging atom with the surface prior to reaction but temperature effects, when modeled via a generalized Langevin oscillator model, do not affect significantly reactivity.

Cross sections and low temperature rate coefficients for the H + CH+ reaction: a quasiclassical trajectory study

The H + CH(+) reaction is studied by quasiclassical trajectory (QCT) calculations, along with phase space theory (PST) and quantum rigid rotor calculations, employing a global single-valued potential energy surface recently derived by our group. We report QCT total cross sections for each of the three channels, for low collision energies and different reactant rotational quantum numbers. At the lowest collision energies, all cross sections exhibit a capture-like behaviour, as expected from a barrierless reaction. At higher energies, there are important dynamical effects coming from the opening of new channels in the inelastic and reactive exchange collisions. The inelastic cross sections turn out to largely increase, while the reactive abstraction cross sections are declining faster than predicted by the capture theory. A large value of the reactant rotational quantum number tends to suppress these dynamical effects. The QCT rate coefficients are reported for a temperature range from 1-700 K. Below 20 K, the abstraction and exchange QCT rate coefficients are almost constant, as predicted by the capture theory. Above this temperature, the abstraction rate coefficient declines, while the exchange and inelastic rate coefficients are increasing, due to the opening of new channels. A good agreement is observed between the experimental abstraction rate coefficient and the QCT and PST ones. The QCT inelastic results are also compared with those obtained from rigid rotor close coupling (CCRR) calculations in order to check the ability of this approach to provide a reliable estimate of the inelastic rate coefficients for a reactive system without a barrier. The laws of variation as a function of temperature are found to be very similar and the curves are parallel above 20 K. However, reaction is not allowed in the rigid rotor approximation, therefore the CCRR results are about twice as large as their QCT counterparts.

High-harmonic transient grating spectroscopy of NO2 electronic relaxation.

We study theoretically and experimentally the electronic relaxation of NO(2) molecules excited by absorption of one ∼400 nm pump photon. Semiclassical simulations based on trajectory surface hopping calculations are performed. They predict fast oscillations of the electronic character around the intersection of the ground and first excited diabatic states. An experiment based on high-order harmonic transient grating spectroscopy reveals dynamics occurring on the same time scale. A systematic study of the detected transient is conducted to investigate the possible influence of the pump intensity, pump wavelength, and rotational temperature of the molecules. The quantitative agreement between measured and predicted dynamics shows that, in NO(2), high harmonic transient grating spectroscopy encodes vibrational dynamics underlying the electronic relaxation.

Recombination and chemical energy accommodation coefficients from chemical dynamics simulations: O/O2 mixtures reacting over a β-cristobalite (001) surface

A microkinetic model is developed to study the reactivity of an O/O(2) gas mixture over a β-cristobalite (001) surface. The thermal rate constants for the relevant elementary processes are either inferred from quasiclassical trajectory calculations or using some statistical approaches, resting on a recently developed interpolated multidimensional potential energy surface based on density functional theory. The kinetic model predicts a large molecular coverage at temperatures lower than 1000 K, in contrary to a large atomic coverage at higher temperatures. The computed atomic oxygen recombination coefficient, mainly involving atomic adsorption and Eley-Rideal recombination, is small and increases with temperature in the 700-1700 K range (0.01 < γ(O) < 0.02) in good agreement with experiments. In the same temperature range, the estimated chemical energy accommodation coefficient, the main contribution to which is the atomic adsorption process is almost constant and differs from unity (0.75 < β(O) < 0.80).