Experimental Spin-Orbit State-Resolved Differential Cross Sections of the S(1D) + D2 → SD(2Π3/2,1/2) + D Reaction at Collision Energies of 266.2 and 206.5 cm-1

The S(1D) + D2 → SD + D reaction is a prototype insertion chemical reaction that involves spin-orbit interactions in the exit channel. In this work, we report spin-orbit state-resolved differential cross sections (DCSs) of this reaction obtained by crossed beam experiments at collision energies of 266.2 and 206.5 cm-1. The DCSs of specific rovibrational states exhibit a slight preference for forward scattering. When integrated over all rotational quantum states within each spin-orbit manifold, the total angular distributions of the two manifolds show nearly forward-backward symmetry, indicating that the deep well responsible for the long-living complex-forming mechanism predominates the entire reaction dynamics. Moreover, significant spin-orbit preference was observed at rotational quantum number N > 9 in the vibrationally ground state of SD products. It was also observed that SD products in the vibrationally excited state v' = 1 prefer to populate in the 2Π3/2 manifold, with the 2Π3/2/2Π1/2 ratio of 15.8 and 25.2 at collision energies of 266.2 and 206.5 cm-1, respectively. The experimental spin-orbit state-resolved DCSs obtained in this work will be of great importance for developing an accurate diabatic theory that includes spin-orbit interactions for this title reaction.

Scattering resonances in the rotational excitation of HDO by Ne and normal-H2: theory and experiment

The rotational excitation of a singly deuterated water molecule (HDO) by a heavy atom (Ne) and a light diatomic molecule (H2) is investigated theoretically and experimentally in the near-threshold regime. Crossed-molecular-beam measurements with a variable crossing angle are compared to close-coupling calculations based on high-accuracy potential energy surfaces. The two lowest rotational transitions, 000 → 101 and 000 → 111, are probed in detail and a good agreement between theory and experiment is observed for both transitions in the case of HDO + Ne, where scattering resonances are however blurred out experimentally. In the case of HDO + H2, the predicted theoretical overlapping resonances are faithfully reproduced by experiment for the 000 → 111 transition, while the calculated strong signal for the 000 → 101 transition is not detected. Future work is needed to reconcile this discrepancy.

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.

Universal behavior in complex-mediated reactions: Dynamics of S(1D) + o-D2 → D + SD at low collision energies

Reactive and elastic cross sections and rate coefficients have been calculated for the S(1D) + D2(v = 0, j = 0) reaction using a modified hyperspherical quantum reactive scattering method. The considered collision energy ranges from the ultracold regime, where only one partial wave is open, up to the Langevin regime, where many of them contribute. This work presents the extension of the quantum calculations, which in a previous study were compared with the experimental results, down to energies in the cold and ultracold domains. Results are analyzed and compared with the universal case of the quantum defect theory by Jachymski et al. [Phys. Rev. Lett. 110, 213202 (2013)]. State-to-state integral and differential cross sections are also shown covering the ranges of low-thermal, cold, and ultracold collision energy regimes. It is found that at E/kB < 1 K, there are substantial departures from the expected statistical behavior and that dynamical features become increasingly important with decreasing collision energy, leading to vibrational excitation.

Time-dependent wave packet dynamics study of the resonances in the H + LiH+(v = 0, j = 0) → Li+ + H2 reaction at low collision energies

The depletion process of LiH⁺ by H collision plays an important role in the evolution of the early universe and astrophysical processes, including the eventual charge-states, abundances of atomic and molecular species and ensuing astrochemistry. Here, a quantum dynamics study on the H + LiH⁺(v = 0, j = 0) → Li⁺ + H₂ reaction is performed at the low collision energy range from 0.1 meV to 10 meV using the time-dependent wave packet method. A Feshbach resonance peak is observed near 0.8 meV collision energy on the total reaction probability curves. This resonance originates from the coupling with the v = 0, j = 1 energy level of the reactant LiH⁺, and it is dominated by the contributions of J = 0-4 partial waves. Another partial wave resonance is also found on the total integral cross section at 1.2 meV, which is closely connected to the opening of the J = 7 partial wave. The opening of the J = 7 partial wave generates a notable forward scattering peak, and the Feshbach resonance can promote both the forward and backward scatterings. Moreover, the total and product vibrational state-resolved rate coefficients for the temperature range of 1-100 K are also reported.

Low-temperature reaction dynamics of paramagnetic species in the gas phase

Radicals are abundant in a range of important gas-phase environments. They are prevalent in the atmosphere, in interstellar space, and in combustion processes. As such, understanding how radicals react is essential for the development of accurate models of the complex chemistry occurring in these gas-phase environments. By controlling the properties of the colliding reactants, we can also gain insights into how radical reactions occur on a fundamental level. Recent years have seen remarkable advances in the breadth of experimental methods successfully applied to the study of reaction dynamics involving paramagnetic species-from improvements to the well-known crossed molecular beams approach to newer techniques involving magnetically guided and decelerated beams. Coupled with ever-improving theoretical methods, quantum features are being observed and interesting insights into reaction dynamics are being uncovered in an increasingly diverse range of systems. In this highlight article, we explore some of the exciting recent developments in the study of chemical dynamics involving paramagnetic species. We focus on low-energy reactive collisions involving neutral radical species, where the reaction parameters are controlled. We conclude by identifying some of the limitations of current methods and exploring possible new directions for the field.

Signature of shape resonances on the differential cross sections of the S(1D)+H2 reaction

Shape resonances appear when the system is trapped in an internuclear potential well after tunneling through a barrier. They manifest as peaks in the collision energy dependence of the cross section (excitation function), and in many cases, their presence can be observed experimentally. High-resolution crossed-beam experiments on the S(¹D) + H₂(j = 0) reaction in the 0.81-8.5 meV collision energy range reaction revealed non-monotonic behavior and the presence of oscillations in the reaction cross section as a function of the collision energy, as predicted by quantum mechanical (QM) calculations. In this work, we have analyzed the effect of shape resonances on the differential cross sections for this insertion reaction by performing additional QM calculations. We have found that, in some cases, the resonance gives rise to a large enhancement of extreme backward scattering for specific final states. Our results also show that, in order to yield a significant change in the state-resolved differential cross section, the resonance has to be associated with constructive interference between groups of partial waves, which requires not getting blurred by the participation of many product helicity states.

Interaction-Asymptotic Region Decomposition Method for an Insertion Reaction: Application to the S(1D) + H2 Reaction

With adjusting principal axes hyperspherical (APH) coordinate in the interaction region, and the Jacobi coordinates in the asymptotic regions, an efficient multidomain interaction-asymptotic region decomposition (IARD) method has been developed to solve the "coordinate problem" in a product-state-resolved reactive scattering calculation using the quantum wave packet method. Although the APH coordinate treats with all three channels equally, and is efficient for describing the interaction region for some direct reactions, it is inefficient for describing the insertion-type reaction due to the singularity problem, such as the S(¹D) + H₂ reaction. To deal with this issue, in this work, the channel-dependent Delves hyperspherical (DH) coordinate is proposed to describe the interaction region using the IARD method. The proposed DH-IARD method was applied to calculate the product-state-resolved reaction probabilities of the H + HD reaction, and the differential and integral cross sections of the typical insertion reaction S(¹D) + H₂. It is found that the new DH-IARD method is much more efficient than the previous APH-IARD method for dealing with insertion reactions. The partial wave resonance structures were observed in the integral cross section. It is found that at a low collision energy, the position of the initial wave packet has to be put far away. Otherwise, the partial wave resonance structures could not be correctly reproduced due to the reef well arising with a large total angular momentum J.

Dynamical resonances in chemical reactions

The transition state is a key concept in the field of chemistry and is important in the study of chemical kinetics and reaction dynamics.

A review of dynamical resonances in A + BC chemical reactions

The concept of the transition state has played an important role in the field of chemical kinetics and reaction dynamics. Reactive resonances in the transition-state region can dramatically enhance the reaction probability; thus investigation of the reactive resonances has attracted great attention from chemical physicists for many decades. In this review, we mainly focus on the recent progress made in probing the elusive resonance phenomenon in the simple A  +  BC reaction and understanding its nature, especially in the benchmark F/Cl  +  H₂ and their isotopic variants. The signatures of reactive resonances in the integral cross section, differential cross section (DCS), forward- and backward-scattered DCS, and anion photodetachment spectroscopy are comprehensively presented in individual prototype reactions. The dynamical origins of reactive resonances are also discussed in this review, based on information on the wave function in the transition-state region obtained by time-dependent quantum wave-packet calculations.

An Empirical Dynamical Barrier for Statistical Theory of Low-Energy Reactive S(1D) + HD(j = 0), H2(j = 0) Collisions

A simple model potential is proposed to describe the dynamical barrier in the mean interaction potential at small distances between the reactants in S(¹D) + HD(¹Σ, v = 0, j = 0) reaction. The statistical theory of collision complex formation and complex decay is applied to calculate the total reaction cross sections and the cross sections for SH and SD productions in the range of low collision energies E_c = (0.4-60) meV. The results are compared with measured cross sections and results of hyperspherical close coupling calculations. As a check of consistency the same comparisons are presented for the case of S(¹D) + H₂(¹Σ, v = 0, j = 0) reaction.

INVESTIGATION OF THE EXCHANGE REACTION H + H′S → HS + H′ ON THE 1A′ STATE POTENTIAL ENERGY SURFACE

The quantum time dependent wave packet (TDWP) and quasiclassical trajectory (QCT) calculations were carried out to study the exchange reaction H(2S) + H′S(2Π) → HS(2Π) + H′(2S) on the 1A′ potential energy surface (PES). The integral cross sections of the H + H′S (v = j = 0) → HS + H′ reaction calculated by the two methods were presented. The results reveal that the integral cross sections (ICS) decrease with the collision energy increasing. The result of the QCT calculations is reasonably consistent with the time-dependent wave packet. Moreover, the differential cross sections (DCS) were calculated by the QCT method at the four different collision energies, which display a forward–backward symmetry. A long-lifetime H2S intermediate complex of the exchange reaction was found according to the trajectories. In the stereodynamics investigation, the polar and dihedral angle distribution functions were calculated, which have the distinct oscillations. The oscillations could be attributed to the deep well on the 1A′ PES. However, based on the polar-angle and dihedral angle distribution functions, it could be predicted that the main product rotational angular momentum preferentially point to the positive or negative direction of y-axes.

Rate coefficients from quantum and quasi-classical cumulative reaction probabilities for the S(1D) + H2 reaction

Cumulative reaction probabilities (CRPs) at various total angular momenta have been calculated for the barrierless reaction S((1)D) + H(2) → SH + H at total energies up to 1.2 eV using three different theoretical approaches: time-independent quantum mechanics (QM), quasiclassical trajectories (QCT), and statistical quasiclassical trajectories (SQCT). The calculations have been carried out on the widely used potential energy surface (PES) by Ho et al. [J. Chem. Phys. 116, 4124 (2002)] as well as on the recent PES developed by Song et al. [J. Phys. Chem. A 113, 9213 (2009)]. The results show that the differences between these two PES are relatively minor and mostly related to the different topologies of the well. In addition, the agreement between the three theoretical methodologies is good, even for the highest total angular momenta and energies. In particular, the good accordance between the CRPs obtained with dynamical methods (QM and QCT) and the statistical model (SQCT) indicates that the reaction can be considered statistical in the whole range of energies in contrast with the findings for other prototypical barrierless reactions. In addition, total CRPs and rate coefficients in the range of 20-1000 K have been calculated using the QCT and SQCT methods and have been found somewhat smaller than the experimental total removal rates of S((1)D).

Dynamics of the S(1D2)+HD(j=0) reaction at collision energies approaching the cold regime: a stringent test for theory

We report integral cross sections for the S(1D2)+HD(j=0)→DS+H and HS+D reaction channels obtained through crossed-beam experiments reaching collision energies as low as 0.46 meV and from adiabatic time-independent quantum-mechanical calculations. While good overall agreement with experiment at energies above 10 meV is observed, neither the product channel branching ratio nor the low-energy resonancelike features in the HS+D channel can be theoretically reproduced. A nonadiabatic treatment employing highly accurate singlet and triplet potential energy surfaces is clearly needed to resolve the complex nature of the reaction dynamics.

Low temperature kinetics of unstable radical reactions

Recent advances in Earth and satellite based observations of molecules in interstellar environments and in planetary atmospheres have highlighted the lack of information regarding many important gas-phase formation mechanisms involving neutral species at low temperatures. Whilst significant progress has been made towards a better understanding of radical-molecule reactions in these regions, the inherent difficulties involved in the investigation of reactions between two unstable radical species have hindered progress in this area. This perspective article provides a brief review of the most common techniques applied to study radical-radical reactions below room temperature, before outlining the developments in our laboratory that have allowed us to extend such measurements to temperatures relevant to astrochemical environments. These developments will be discussed with particular emphasis on our recent investigations of the reactions of atomic nitrogen with diatomic radicals.

Appearance of low energy resonances in CO-para-H2 inelastic collisions

We report on crossed-beam experiments and quantum-mechanical calculations performed on the CO(j=0) + H2(j=0) → CO(j=1) + H2(j=0) system. The experimental cross sections determined in the threshold region of the CO(j=0 → j=1) transition at 3.85  cm(-1) show resonance structures in good qualitative agreement with the theoretical ones. These results suggest that the potential energy surface which describes the CO-H2 van der Waals interaction should be reinvestigated for good quantitative agreement.