Reaction Cross Section and Rate Constant Calculations for the D + H2(v=0,1) → HD + H Reaction on Three ab Initio Potential Energy Surfaces. A Quasiclassical Trajectory Study

Li + HF and Li + HCl Reactions Revisited I: QCT Calculations and Simulation of Experimental Results

The Li + HF and Li + HCl reactions share some common features. They have the same kinematics, relatively small barrier heights, bent transition states, and are both exothermic when the zero point energy is considered. Nevertheless, the pioneering crossed beam experiments by Lee and co-workers in the 80s (Becker et al., J. Chem. Phys. 1980, 73, 2833) revealed that the dynamics of the two reactions differ significantly, especially at low collision energies. In this work, we present theoretical simulations of their results in the laboratory frame (LAB), based on quasiclassical trajectories and obtained using accurate potential energy surfaces. The calculated LAB angular distributions and time-of-flight spectra agree well with the raw experimental data, although our simulations do not reproduce the experimentally derived center-of-mass (CM) differential cross section and velocity distributions. The latter were derived by forward convolution fitting under the questionable assumption that the CM recoil velocity and scattering angle distribution were uncoupled, while our results show that the coupling between them is relevant. Some important insights into the reaction mechanism discussed in the article by Becker et al. had not been contrasted with those that can be extracted from the theoretical results. Among them, the correlation between the angular momenta involved in the reactions has also been examined. Given the kinematics of both systems, the reagent orbital angular momentum, l , is almost completely transformed into the rotation of the product diatom, j'. However, contrary to the coplanar mechanism proposed in the original paper, we find that the initial and final relative orbital angular momenta are not necessarily parallel. Both reactions are found to be essentially direct, although about 15% of the LiFH complexes live longer than 200 fs.

Prediction of a Feshbach Resonance in the Below-the-Barrier Reactive Scattering of Vibrationally Excited HD with H

Quantum reactive scattering calculations on the vibrational quenching of HD due to collisions with H were carried out employing an accurate potential energy surface. The state-to-state cross sections for the chemical reaction HD(v = 1, j = 0) + H → D + H₂(v' = 0, j') at collision energies between 1 and 10 000 cm^-1 are presented, and a Feshbach resonance in the low-energy regime, below the reaction barrier, is observed for the first time. The resonance is attributed to coupling with the vibrationally adiabatic potential correlating to the v = 1, j = 1 level of the HD molecule, and it is dominated by the contribution from a single partial wave. The properties of the resonance, such as its dynamic behavior, phase behavior, and lifetime, are discussed.

Experimental and theoretical studies of the Xe-OH(A/X) quenching system

New multi-reference, global ab initio potential energy surfaces (PESs) are reported for the interaction of Xe atoms with OH radicals in their ground X²Π and excited A²Σ⁺ states, together with the non-adiabatic couplings between them. The 2A' excited potential features a very deep well at the collinear Xe-OH configuration whose minimum corresponds to the avoided crossing with the 1A' PES. It is therefore expected that, as with collisions of Kr + OH(A), electronic quenching will play a major role in the dynamics, competing favorably with rotational energy transfer within the 2A' state. The surfaces and couplings are used in full three-state surface-hopping trajectory calculations, including roto-electronic couplings, to calculate integral cross sections for electronic quenching and collisional removal. Experimental cross sections, measured using Zeeman quantum beat spectroscopy, are also presented here for comparison with these calculations. Unlike similar previous work on the collisions of OH(A) with Kr, the surface-hopping calculations are only able to account qualitatively for the experimentally observed electronic quenching cross sections, with those calculated being around a factor of two smaller than the experimental ones. However, the predicted total depopulation of the initial rovibrational state of OH(A) (quenching plus rotational energy transfer) agrees well with the experimental results. Possible reasons for the discrepancies are discussed in detail.

Influence of the Reactants Rotational Excitation on the H + D2(v = 0, j) Reactivity

We have analyzed the influence of the rotational excitation on the H + D2(v = 0, j) reaction through quantum mechanical (QM) and quasiclassical trajectories (QCT) calculations at a wide range of total energies. The agreement between both types of calculations is excellent. We have found that the rotational excitation largely increases the reactivity at large values of the total energy. Such an increase cannot be attributed to a stereodynamical effect but to the existence of recrossing trajectories that become reactive as the target molecule gets rotationally excited. At low total energies, however, recrossing is not significant and the reactivity evolution is dominated by changes in the collision energy; the reactivity decreases with the collision energy as it shrinks the acceptance cone. When state-to-state results are considered, rotational excitation leads to cold product's rovibrational distributions, so that most of the energy is released as recoil energy.

Vibrational excitation influence on the H + BrO → HBr + O reaction: a quasi-classical trajectory investigation

Quasi-classical trajectory calculations are employed to investigate the vibrational excitation effect on the scalar and vector properties of the H + BrO → HBr + O reaction using a X1A′ state ab initio potential energy surface (J. Chem. Phys. 2000, 113, 4598). The reaction probability, cross section, and rate constant are carried out with the effect of the collision energy (Ecol = 0.1–6 kcal/mol) and vibrational levels (v = 0–3). A significant vibrational dependency has been observed in the reaction probability and cross section at a relatively low collision energy area and has also been found in a low-temperature (T < 150 K) region of the rate constant. In addition, two product angular distributions, P(θr) and P(ϕr), and two generalized polarization-dependent differential cross sections, PDDCS00 and PDDCS20, are calculated as well. All of these scalar and vector properties have shown sensitive behaviors to the vibrational levels.

The Cl + O3 reaction: a detailed QCT simulation of molecular beam experiments

QCT calculations have been carried out to determine angle–velocity differential cross-sections to simulate the results of molecular beam experiments.

Cumulative reaction probabilities: A comparison between quasiclassical and quantum mechanical results

This article presents a quasiclassical trajectory (QCT) method for determining the cumulative reaction probability (CRP) as a function of the total energy. The method proposed is based on a discrete sampling using integer values of the total and orbital angular momentum quantum numbers for each trajectory and on the development of equations that have a clear counterpart in the quantum mechanical (QM) case. The calculations comprise cumulative reaction probabilities at a given total angular momentum J, as well as those summed over J. The latter are used to compute QCT rate constants. The method is illustrated by comparing QCT and exact QM results for the H+H2, H+D2, D+H2, and H+HD reactions. The agreement between QCT and QM results is very good, with small discrepancies between the two data sets indicating some genuine quantum effects. The most important of these involves the value of the CRP at low energies which, due to the absence of tunneling, is lower in the QCT calculations, causing the corresponding rate constants to be smaller. The second is the steplike structure that is clearly displayed in the QM CRP for J = 0, which is much smoother in the corresponding QCT results. However, when the QCT density of reactive states, i.e., the derivatives of the QCT CRP with respect to the energy, is calculated, a succession of maxima and minima is obtained which roughly resembles those found in the QM calculations, although the latter are considerably sharper. The analysis of the broad peaks in the QCT density of reactive states indicates that the distributions of collision times associated with the maxima are somewhat broader, with a tail extending to larger collision times, than those associated with the minima. In addition, the QM and QCT dynamics of the isotopic variants mentioned above are compared in the light of their CRPs. Issues such as the compliance of the QCT CRP with the law of microscopic reversibility, as well as the similarity between the CRPs for ortho and para species in the QM and QCT cases, are also addressed.