Dynamical effects on the O(3P) + D2 reaction and its impact on the Λ-doublet population

The O(3P) + D2 → OD(2Π) + D reaction presents the peculiarity of taking place on two different potential energy surfaces (PESs) of different symmetry, 3A' and 3A'', which become degenerate for collinear configurations where the saddle-point of the reaction is located. The degeneracy is broken for non-collinear approaches with the energy on the 3A' PES rising more abruptly with the bending angle, making the frequency of this mode higher on the 3A' state. Consequently, the 3A' PES should be less reactive than the 3A'' one. Nevertheless, quantum scattering calculations show that the cross section is higher on the 3A' PES for energies close to the classical reaction threshold and rotationless reactant. It is found that the differences between the reactivity on the two PESs are greater for low values of total angular momentum, where the centrifugal barrier is lower and contribute to the higher population of the Π(A') Λ-doublet states of OD at low collision energies. At high collision energies, the Π(A') Λ-doublet state is also preferentially populated. Analysis of the differential cross sections reveals that the preponderance for the Π(A') Λ-doublet at low energies comes from backward scattering, originating from the reaction on the 3A' PES, while at high energies, it proceeds from a different mechanism that leads to sideways scattering on the 3A'' PES and that populates the Π(A') manifold.

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.

The F + HD(v = 0, 1; j = 0, 1) reactions: stereodynamical properties of orbiting resonances

The excitation functions (reaction cross-section as a function of collision energy) of the F + HD(v = 0, 1; j = 0, 1) benchmark system have been calculated in the 0.01-6 meV collision energy interval using a time-independent hyperspherical quantum dynamics methodology. Special attention has been paid to orbiting resonances, which bring about detailed information on the three-atom interaction during the reactive encounter. The location of the resonances depends on the rovibrational state of the reactants HD(v,j), but is the same for the two product channels HF + D and DF + H, as expected for these resonances that are linked to the van der Waals well at the entrance. The resonance intensities depend both on the entrance and on the exit channels. The peak intensities for the HF + D channel are systematically larger than those for DF + H. Vibrational excitation leads to an increase of the peak intensity by more than an order of magnitude, but rotational excitation has a less drastic effect. It deceases the resonance intensity of the F + HD(v = 1) reaction, but increases somewhat that of F + HD(v = 0). Polarization of the rotational angular momentum with respect to the initial velocity reveals intrinsic directional preferences in the F + HD(v = 0, 1; j = 1) reactions that are manifested in the resonance patterns. The helicities (Ω = 0, Ω = ±1) possible for j = 1 contribute to the resonances, but that from Ω± 1 is, in general, dominant and in some cases exclusive. It corresponds to a preferential alignment of the HD internuclear axis perpendicular to the initial direction of approach and, thus, to side-on collisions. This work also shows that external preparation of the reactants, following the intrinsic preferences, would allow the enhancement or reduction of specific resonance features, and would be of great help for their eventual experimental detection.

Quantum study of reaction O(3 P) + H2 (v, j) → OH + H: OH formation in strongly UV-irradiated gas

The reaction between atomic oxygen and molecular hydrogen is an important one in astrochemistry as it regulates the abundance of the hydroxyl radical and serves to open the chemistry of oxygen in diverse astronomical environments. However, the existence of a high activation barrier in the reaction with ground state oxygen atoms limits its efficiency in cold gas. In this study we calculate the dependence of the reaction rate coefficient on the rotational and vibrational state of H2 and evaluate the impact on the abundance of OH in interstellar regions strongly irradiated by far-UV photons, where H2 can be efficiently pumped to excited vibrational states. We use a recently calculated potential energy surface and carry out time-independent quantum mechanical scattering calculations to compute rate coefficients for the reaction O(3 P) + H2 (v, j) → OH + H, with H2 in vibrational states v = 0-7 and rotational states j = 0-10. We find that the reaction becomes significantly faster with increasing vibrational quantum number of H2, although even for high vibrational states of H2 (v = 4-5) for which the reaction is barrierless, the rate coefficient does not strictly attain the collision limit and still maintains a positive dependence with temperature. We implemented the calculated state-specific rate coefficients in the Meudon PDR code to model the Orion Bar PDR and evaluate the impact on the abundance of the OH radical. We find the fractional abundance of OH is enhanced by up to one order of magnitude in regions of the cloud corresponding to A V = 1.3-2.3, compared to the use of a thermal rate coefficient for O + H2, although the impact on the column density of OH is modest, of about 60%. The calculated rate coefficients will be useful to model and interpret JWST observations of OH in strongly UV-illuminated environments.

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.

Unveiling shape resonances in H + HF collisions at cold energies

Scattering resonances are pure quantum effects that appear whenever the collision energy matches the energy of a quasi-bound state of the intermolecular complex. Here we show that rotational quenching of HF(j = 1, 2) with H is strongly influenced by the presence of two resonance peaks, leading to up to a two-fold increase in the thermal rate coefficients at the low temperatures characteristic of the interstellar medium. Our results show that each resonance peak is formed by a cluster of shape resonances, each of them characterized by the same value of the orbital angular momentum but different values of the total angular momentum. The relative intensity of these resonances depends on the relative geometry of the incoming reactants, and our results predict that by changing the alignment of the HF rotational angular momentum it is possible to decompose the resonance peaks, disentangling the underlying resonance pattern and the contribution of different total angular momenta to the resonance.

Unexpected dynamical effects change the lambda-doublet propensity in the tunneling region for the O(3P) + H2 reaction

One of the most relevant features of the O(³P) + H₂ reaction is that it occurs on two different potential energy surfaces (PESs) of symmetries A' and A'' that correlate reactants and products. The respective saddle points, which correspond to a collinear arrangement, are the same for both PESs, whilst the barrier height rises more abruptly on the ³A' PES than on the ³A'' PES. Accordingly, the reactivity on the ³A'' PES should be always higher than on the ³A' PES. In this work, we present accurate quantum-scattering calculations showing that this is not always the case for rotationless reactants, where dynamical factors near the reaction threshold cause the ³A' PES to dominate at energies around the barrier. Further calculation of cross sections and Λ-doublet populations has allowed us to establish how the reaction mechanism changes from the deep tunneling regime to hyperthermal energies.

New Stress Test for Ring Polymer Molecular Dynamics: Rate Coefficients of the O(3P) + HCl Reaction and Comparison with Quantum Mechanical and Quasiclassical Trajectory Results

In the past decade, ring polymer molecular dynamics (RPMD) has emerged as a very efficient method to determine thermal rate coefficients for a great variety of chemical reactions. This work presents the application of this methodology to study the O(³P) + HCl reaction, which constitutes a stringent test for any dynamical calculation due to rich resonant structure and other dynamical features. The rate coefficients, calculated on the ³A' and ³A″ potential energy surfaces (PESs) by Ramachandran and Peterson [ J. Chem. Phys. 2003 , 119 , 9590 ], using RPMD and quasiclassical trajectories (QCT) are compared with the existing experimental and the quantum mechanical (QM) results by Xie et al. [ J. Chem. Phys. 2005 122 , 014301 ]. The agreement is very good at T > 600 K, although RPMD underestimates rate coefficients by a factor between 4 and 2 in the 200-500 K interval. The origin of these discrepancies lies in the large contribution from tunneling on the ³A″ PES, which is enhanced by resonances due to quasibound states in the van der Waals wells. Although tunneling is fairly well accounted for by RPMD even below the crossover temperature, the effect of resonances, a long-time effect, is not included in the methodology. At the highest temperatures studied in this work, 2000-3300 K, the RPMD rate coefficients are somewhat larger than the QM ones, but this is shown to be due to limitations in the QM calculations and the RPMD are believed to be more reliable.

Differential steric effects in the inelastic scattering of NO(X) + Ar: spin-orbit changing transitions

Spin-orbit changing transitions for bond-axis oriented collisions of NO(X) with Ar have been investigated with full quantum state selection via a crossed molecular beam experiment at collision energies of 532 cm^-1 and 651 cm^-1. NO(X) molecules were selected in their ground rotational state (Ω = 0.5, j = 0.5, f) before being adiabatically oriented using a static electric field, such that either the N- or O-end of the molecule was directed towards the incoming Ar atom. After collision, NO(X, Ω' = 1.5, j', e) molecules were probed quantum state specifically using velocity-map ion imaging, coupled with resonantly enhanced multi-photon ionization. Differences were observed between the experimental ion images and differential cross sections for collisions occurring at the two ends of the molecule, with results that could largely be accounted for by quantum mechanical scattering calculations. The bond-axis oriented data for the spin-orbit changing collisions are compared with similar results obtained previously for spin-orbit conserving transitions, and for field free scattering of NO(X) with Ar.

How reactant polarization can be used to change the effect of interference on reactive collisions

It is common knowledge that integral and differential cross sections (DCSs) are strongly dependent on the spatial distribution of the molecular axis of the reactants. Hence, by controlling the axis distribution, it is possible to either promote or hinder the yield of products into specific final states or scattering angles. This idea has been successfully implemented in experiments by polarizing the internuclear axis before the reaction takes place, either by manipulating the rotational angular distribution or by the Stark effect in the presence of an orienting field. When there is a dominant reaction mechanism, characterized by a set of impact parameters and angles of attack, it is expected that a preparation that helps the system to reach the transition state associated with that mechanism will promote the reaction, whilst a different preparation would generally impair the reaction. However, when two or more competing mechanisms via interference contribute to the reaction into specific scattering angles and final states, it is not evident which would be the effect of changing the axis preparation. To address this problem, throughout this article we have simulated the effect that different experimental preparations have on the DCSs for the H + D2 reaction at relatively high energies, for which it has been shown that several competing mechanisms give rise to interference that shapes the DCS. To this aim, we have extended the formulation of the polarization dependent DCS to calculate polarization dependent generalized deflection functions of ranks greater than zero. Our results show that interference is very sensitive to changes in the internuclear axis preparation, and that the shape of the DCS can be controlled exquisitely.

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.

Angular momentum-scattering angle quantum correlation: a generalized deflection function

A natural generalization of the classical deflection function, the functional dependence of the deflection angle on the angular momentum (or the impact parameter), is the joint probability density function of these two quantities, revealing the correlation between them. It provides, at a glance, detailed information about the reaction mechanisms and how changes in the impact parameter affect the product angular distribution. It is also useful to predict the presence of quantum phenomena such as interference. However, the classical angular momentum-scattering angle correlation function has a limited use whenever quantum effects become important. Rigorously speaking, there is not a quantum equivalent of the classical joint distribution, as the differential cross section depends on the coherences between the different values of J caused by the cross terms in the expansion of partial waves. In this article, we present a simple method to calculate a quantum analog of this correlation, a generalized deflection function that can shed light onto the reaction mechanism using just quantum mechanical results. Our results show that there is a very good agreement between the quantum and classical correlation functions as long as quantum effects are not all relevant. When this is not the case, it will also be shown that the quantum correlation function is most useful to observe the extent of quantum effects such as interference among different reaction mechanisms.

An experimental study of OH(A2Σ+) + H2: Electronic quenching, rotational energy transfer, and collisional depolarization

Zeeman quantum beat spectroscopy has been used to determine the thermal (300 K) rate constants for electronic quenching, rotational energy transfer, and collisional depolarization of OH(A²Σ⁺) by H₂. Cross sections for both the collisional disorientation and collisional disalignment of the angular momentum in the OH(A²Σ⁺) radical are reported. The experimental results for OH(A²Σ⁺) + H₂ are compared to previous work on the OH(A²Σ⁺) + He and Ar systems. Further comparisons are also made to the OH(A²Σ⁺) + Kr system, which has been shown to display significant non-adiabatic dynamics. The OH(A²Σ⁺) + H₂ experimental data reveal that collisions that survive the electronic quenching process are highly depolarizing, reflecting the deep potential energy wells that exist on the excited electronic state surface.

The dynamics of the Hg + Br2 reaction: elucidation of the reaction mechanism for the Br exchange reaction

In spite of its importance in the Hg atmospheric chemistry, the dynamics of the Hg + Br₂ → HgBr + Br reaction is poorly understood. In this article, we have carried out a comprehensive study of the reaction mechanism of this reaction by means of quasiclassical trajectories (QCTs) on an existing ab initio potential energy surface (PES). The reaction has a non trivial dynamics, as a consequence of its large endothermicity, the presence of a deep potential well, and the competition between the Br exchange and the collision induced dissociation processes. Our calculations demonstrate that insertion is only relevant at energies just above the reaction threshold and that, at energies above 2.3 eV, HgBr formation typically takes place via a sort of frustrated dissociation. In order to compare directly with the results obtained in extensive cross molecular beam experiments for the homologous reaction with I₂, angular distributions in the laboratory frame for Hg + Br₂ have been simulated under similar experimental conditions. The lack of agreement at the highest energies considered suggests that either the two reactions have substantially different mechanisms or that calculations on a single PES cannot account for the dynamics at those energies.

Angular distributions for the inelastic scattering of NO(X2Π) with O2(X3Σg-)

The inelastic scattering of NO(X²Π) by O₂(X³Σ_g^-) was studied at a mean collision energy of 550 cm^-1 using velocity-map ion imaging. The initial quantum state of the NO(X²Π, v = 0, j = 0.5, Ω=0.5, 𝜖 = -1, f) molecule was selected using a hexapole electric field, and specific Λ-doublet levels of scattered NO were probed using (1+1^') resonantly enhanced multiphoton ionization. A modified "onion-peeling" algorithm was employed to extract angular scattering information from the series of "pancaked," nested Newton spheres arising as a consequence of the rotational excitation of the molecular oxygen collision partner. The extracted differential cross sections for NO(X) f→f and f→e Λ-doublet resolved, spin-orbit conserving transitions, partially resolved in the oxygen co-product rotational quantum state, are reported, along with O₂ fragment pair-correlated rotational state population. The inelastic scattering of NO with O₂ is shown to share many similarities with the scattering of NO(X) with the rare gases. However, subtle differences in the angular distributions between the two collision partners are observed.

Stereodynamics in NO(X) + Ar inelastic collisions

The effect of orientation of the NO(X) bond axis prior to rotationally inelastic collisions with Ar has been investigated experimentally and theoretically. A modification to conventional velocity-map imaging ion optics is described, which allows the orientation of hexapole state-selected NO(X) using a static electric field, followed by velocity map imaging of the resonantly ionized scattered products. Bond orientation resolved differential cross sections are measured experimentally for a series of spin-orbit conserving transitions and compared with quantum mechanical calculations. The agreement between experimental results and those from quantum mechanical calculations is generally good. Parity pairs, which have previously been observed in collisions of unpolarized NO with various rare gases, are not observed due to the coherent superposition of the two j = 1/2, Ω = 1/2 Λ-doublet levels in the orienting field. The normalized difference differential cross sections are found to depend predominantly on the final rotational state, and are not very sensitive to the final Λ-doublet level. The differential steric effect has also been investigated theoretically, by means of quantum mechanical and classical calculations. Classically, the differential steric effect can be understood by considering the steric requirement for different types of trajectories that contribute to different regions of the differential cross section. However, classical effects cannot account quantitatively for the differential steric asymmetry observed in NO(X) + Ar collisions, which reflects quantum interference from scattering at either end of the molecule. This quantum interference effect is dominated by the repulsive region of the potential.

Integral steric asymmetry in the inelastic scattering of NO(X2Π)

The integral steric asymmetry for the inelastic scattering of NO(X) by a variety of collision partners was recorded using a crossed molecular beam apparatus. The initial state of the NO(X, v = 0, j = 1/2, Ω=1/2, ϵ=-1,f) molecule was selected using a hexapole electric field, before the NO bond axis was oriented in a static electric field, allowing probing of the scattering of the collision partner at either the N- or O-end of the molecule. Scattered NO molecules were state selectively probed using (1 + 1') resonantly enhanced multiphoton ionisation, coupled with velocity-map ion imaging. Experimental integral steric asymmetries are presented for NO(X) + Ar, for both spin-orbit manifolds, and Kr, for the spin-orbit conserving manifold. The integral steric asymmetry for spin-orbit conserving and changing transitions of the NO(X) + O₂ system is also presented. Close-coupled quantum mechanical scattering calculations employing well-tested ab initio potential energy surfaces were able to reproduce the steric asymmetry observed for the NO-rare gas systems. Quantum mechanical scattering and quasi-classical trajectory calculations were further used to help interpret the integral steric asymmetry for NO + O₂. Whilst the main features of the integral steric asymmetry of NO with the rare gases are also observed for the O₂ collision partner, some subtle differences provide insight into the form of the underlying potentials for the more complex system.

Product lambda-doublet ratios as an imprint of chemical reaction mechanism

In the last decade, the development of theoretical methods has allowed chemists to reproduce and explain almost all of the experimental data associated with elementary atom plus diatom collisions. However, there are still a few examples where theory cannot account yet for experimental results. This is the case for the preferential population of one of the Λ-doublet states produced by chemical reactions. In particular, recent measurements of the OD(²Π) product of the O(³P)+D₂ reaction have shown a clear preference for the Π(A′) Λ-doublet states, in apparent contradiction with ab initio calculations, which predict a larger reactivity on the A′′ potential energy surface. Here we present a method to calculate the Λ-doublet ratio when concurrent potential energy surfaces participate in the reaction. It accounts for the experimental Λ-doublet populations via explicit consideration of the stereodynamics of the process. Furthermore, our results demonstrate that the propensity of the Π(A′) state is a consequence of the different mechanisms of the reaction on the two concurrent potential energy surfaces

Propensity for a given Λ-doublet level is a common feature in many chemical reactions, but has so far remained unexplained. Here, the authors show how to predict computationally those propensities and relate them to the reaction mechanism on concurrent potential energy surfaces.

Influence of vibration in the reactive scattering of D + MuH: the effect of dynamical bonding

The dynamics of the D + MuH(v = 1) reaction has been investigated using time-independent quantum mechanical calculations. The total reaction cross sections and rate coefficients have been calculated for the two exit channels of the reaction leading, respectively, to DMu + H and DH + Mu. Over the 100-1000 K temperature range investigated the rate coefficients for the DMu + H channel are of the order of 10(-10) cm(3) s(-1) and those for the DH + Mu channel vary between 1 × 10(-12) and 8 × 10(-11) cm(3) s(-1). These results point to a virtually barrierless reaction for the DMu + H channel and to the presence of a comparatively small barrier for the DH + Mu channel and are consistent with the profiles of their respective collinear vibrationally adiabatic potentials (VAPs). The effective barrier in the VAP of the DH + Mu channel is located in the reactant valley and, consequently, translation is found to be more efficient than vibration for the promotion of the reaction over a large energy interval in the post threshold region. Below this barrier, the DH + Mu channel can be accessible through an indirect mechanism implying crossing from the DMu + H pathway. The most salient feature found in the present study is revealed in the total reaction cross section for the DMu + H channel, which shows a sharp resonance caused by the presence of a deep well in the vibrationally adiabatic potential. This well has a dynamical origin, reminiscent of that found recently in the vibrationally bonded BrMuBr complex [Fleming, et al., Angew. Chem., Int. Ed., 2014, 53, 1], and is due to the stabilizing effect of the light Mu atom oscillating between the heavier H and D isotopes and to the bond softening associated with vibrational excitation of MuH.

Correction: Effects of reagent rotation on interferences in the product angular distributions of chemical reactions

Correction for ‘Effects of reagent rotation on interferences in the product angular distributions of chemical reactions’ by P. G. Jambrina et al., Chem. Sci., 2016, 7, 642–649.

Effects of reagent rotation on interferences in the product angular distributions of chemical reactions

Differential cross sections (DSCs) of the HD(v', j') product for the reaction of H atoms with supersonically cooled D2 molecules in a small number of initial rotational states have been measured at a collision energy of 1.97 eV. These DCSs show an oscillatory pattern that results from interferences caused by different dynamical scattering mechanisms leading to products scattered into the same solid angle. The interferences depend on the initial rotational state j of the D2(v = 0, j) reagent and diminish in strength with increasing rotation. We present here a detailed explanation for this behavior and how each dynamical scattering mechanism has a dependence on the helicity Ω, the projection of the initial rotational angular momentum j of the D2 reagent on the approach direction. Each helicity corresponds to a different internuclear axis distribution, with the consequence that the dependence on Ω reveals the preference of the different quasiclassical mechanisms as a function of approach direction. We believe that these results are general and will appear in any reaction for which several mechanisms are operative.

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.

Rotational Orientation Effects in NO(X) + Ar Inelastic Collisions

Rotational angular momentum orientation effects in the rotationally inelastic collisions of NO(X) with Ar have been investigated both experimentally and theoretically at a collision energy of 530 cm(-1). The collision-induced orientation has been determined experimentally using a hexapole electric field to select the ϵ = -1 Λ-doublet level of the NO(X) j = 1/2 initial state. Fully quantum state resolved polarization-dependent differential cross sections were recorded experimentally using a crossed molecular beam apparatus coupled with a (1 + 1') resonance-enhanced multiphoton ionization detection scheme and subsequent velocity-map imaging. To determine the NO sense of rotation, the probe radiation was circularly polarized. Experimental orientation polarization-dependent differential cross sections are compared with those obtained from quantum mechanical scattering calculations and are found to be in good agreement. The origin of the collision-induced orientation has been investigated by means of close-coupled quantum mechanical, quantum mechanical hard shell, quasi-classical trajectory (QCT), and classical hard shell calculations at the same collision energy. Although there is evidence for the operation of limiting classical mechanisms, the rotational orientation cannot be accounted for by QCT calculations and is found to be strongly influenced by quantum mechanical effects.

Steric effects and quantum interference in the inelastic scattering of NO(X) + Ar

New measurements of the differential steric effect for NO + Ar inelastic scattering highlight the importance of quantum interference.

Rotationally inelastic collisions of NO(X) with Ar are investigated in unprecedented detail using state-to-state, crossed molecular beam experiments. The NO(X) molecules are selected in the Ω = 0.5, j = 0.5, f state and then oriented such that either the ‘N’ or ‘O’ end of the molecule is directed towards the incoming Ar atom. Velocity map ion imaging is then used to probe the scattered NO molecules in well-defined quantum states. We show that the fully quantum state-resolved differential steric asymmetry, which quantifies how the relative efficiency for scattering off the ‘O’ and the ‘N’ ends of the molecule varies with scattering angle, is strongly affected by quantum interference. Significant changes in both integral and differential cross sections are found depending on whether collisions occur with the N or O ends of the molecule. The results are well accounted for by rigorous quantum mechanical calculations, in contrast to both classical trajectory calculations and more simplistic models that provide, at best, an incomplete picture of the dynamics.

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.

A new perspective: imaging the stereochemistry of molecular collisions

The concept of the steric effect plays a central role in chemistry. This Perspective describes how the polarization of reactant molecules in space can be used to probe directly the steric effect, and highlights some of the new measurements that are made possible by coupling reactant orientation and alignment with ion imaging techniques.

The effect of the reactant internal excitation on the dynamics of the C(+) + H2 reaction

We have performed a dynamical study of the endothermic and barrierless C(+) + H2((1)Σg(+)) → CH(+)((1)Σg(+)) + H reaction for different initial rotational states of the H2(v = 0) and H2(v = 1) manifolds. The calculations have been carried out using quasiclassical trajectories and the Gaussian binning methodology on a recent potential energy surface [R. Warmbier and R. Schneider, Phys. Chem. Chem. Phys., 2011, 13, 10285]. Both state-selected integral cross sections as a function of the collision energy and rate coefficients, kv,j(T), have been determined. We show that rotational excitation of the reactants is as effective as vibrational excitation when it comes to increasing the reactivity, and that both types of excitation could contribute to explain the unexpectedly high abundance of CH(+) in the interstellar media. Such an increase in reactivity takes place by suppressing the reaction threshold when the internal energy is sufficient to overcome the endothermicity. Whenever this is the case, the excitation functions at collision energies Ecoll ≤ 0.1 eV display a ∝E(-1/2)coll dependence. However, the absolute values of the state selected kv=1(T) are one order of magnitude below the Langevin model predictions. The disagreement between the approximately derived experimental rate coefficients for v = 1 and those calculated by this and previous theoretical treatments is due to the neglect of the effect of the rotational excitation in the derivation of the former. In spite of the deep well present in the potential energy surface, the reaction does not show a statistical behaviour.

Stress Test for Quantum Dynamics Approximations: Deep Tunneling in the Muonium Exchange Reaction D + HMu → DMu + H

Quantum effects play a crucial role in chemical reactions involving light atoms at low temperatures, especially when a light particle is exchanged between two heavier partners. Different theoretical methodologies have been developed in the last decades attempting to describe zero-point energy and tunneling effects without abandoning a classical or semiclassical framework. In this work, we have chosen the D + HMu → DMu + H reaction as a stress test system for three well-established methods: two representative versions of transition state theory (TST), canonical variational theory and semiclassical instanton, and ring polymer molecular dynamics (RPMD). These calculations will be compared with accurate quantum mechanical results. Despite its apparent simplicity, the exchange of the extremely light muonium atom (0.114 u) becomes a most challenging reaction for conventional methods. The main result of this work is that RPMD provides an overall better performance than TST-based methods for such a demanding reaction. RPMD might well turn out to be a useful tool beyond TST applicability.

Origin of collision-induced molecular orientation

Collision-induced rotational angular momentum orientation is a fundamental property of molecular scattering, which is sensitive to the balance between attractive and repulsive forces at play during collision. Here, we quantify a new mechanism leading to orientation, which is purely quantum mechanical in origin. Although the new mechanism is quite general, and will operate more widely in atomic and molecular scattering, it is observed here for impulsive hard shell collisions, for which the orientation vanishes classically. The quantum mechanism can thus be studied in isolation from other processes. The orientation is proposed to originate from the nonlocal nature of the quantum mechanical collision encounter.

The reactive collision mechanism evinced: stereodynamical control of the elementary Br + H2 → H + HBr reaction

From a kinetics standpoint, reactive molecular collisions are the building blocks of the mechanisms of chemical reactions. In contrast, a dynamics standpoint reveals molecular collisions to have their own internal mechanisms, which are not mere theoretical abstractions: through suitable preparation of the reactants internal and stereochemical states, features of the mechanisms of a reactive molecular collision can be made evident and used as "handles" to control the reaction outcome. Using time-independent quantum dynamical calculations, we demonstrate this for the Br + H2(v = 0-1, j = 2) → H + HBr reaction in the 1.0-1.6 eV range of total energies. Despite its pronounced effect on reactivity, which is in agreement with the predictions from Polanyi rules, reactant vibration is found to have little effect on the mechanism of this endoergic, late-barrier reaction. Analysis of the correlations between directional reaction properties shows that the collision stereochemistry strongly depends on the total energy, but not on how this energy is partitioned between reactant translation and vibration. The stereodynamical preferences implied by the collision mechanisms determine how and to what extent one can control the reaction. Regarding the overall reaction, the extent of control is found to be large near the reaction threshold but not when the total energy is high. Regarding state-to-state reactions, the effect of reactant stereochemistry on the product rotational state distribution is found to be nontrivial and energy dependent.

Reaction dynamics and mechanism of the Cl + HD(v = 1) reaction: a quantum mechanical study

Time-independent quantum mechanical calculations have been performed in order to characterize the dynamics and stereodynamics of Cl + HD reactive collisions. Calculations have been carried out at two different total energy values and for various initial states using the adiabatic potential energy surface by Bian and Werner [J. Chem. Phys. 2000, 112, 220]. Special attention has been paid to the reaction with HD(v = 1) for which integral and differential cross-sections have been calculated and the effect of vibrational vs translational energy on the reactivity has been examined. In addition, the reactant polarization parameters and polarization-dependent differential cross-sections have been determined. From these results, the spatial preferences of the reaction and the extent of the control of the cross sections achievable through a suitable preparation of the reactants have been also studied. The directional requirements are tighter for the HCl channel than for the DCl one. Formation of the products takes place preferentially when the rotational angular momentum of the HD molecule is perpendicular to the reactants approach direction. Cross-sections and polarization moments computed from the scattering calculations have been compared with experimental results by Kandel et al. [J. Chem. Phys. 2000, 112, 670] for the reaction with HD(v = 1) produced by stimulated Raman pumping. The agreement so obtained is good, and it improves the accordance found in previous calculations with other methodologies and potential energy surfaces.

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).

Rotational alignment effects in NO(X) + Ar inelastic collisions: an experimental study

Rotational angular momentum alignment effects in the rotationally inelastic collisions of NO(X) with Ar have been investigated at a collision energy of 66 meV by means of hexapole electric field initial state selection coupled with velocity-map ion imaging final state detection. The fully quantum state resolved second rank renormalized polarization dependent differential cross sections determined experimentally are reported for a selection of spin-orbit conserving and changing transitions for the first time. The results are compared with the findings of previous theoretical investigations, and in particular with the results of exact quantum mechanical scattering calculations. The agreement between experiment and theory is generally found to be good throughout the entire scattering angle range. The results reveal that the hard shell nature of the interaction potential is predominantly responsible for the rotational alignment of the NO(X) upon collision with Ar.

H + D2 Reaction Dynamics in the Limit of Low Product Recoil Energy

Both experiment and theory recently showed that the H + D2(v = 0, j = 0) → HD(v' = 4, j') + D reactions at a collision energy of 1.97 eV display a seemingly anomalous HD product angular distribution that moves in the backward direction as the value of j' increases and the corresponding energy available for product recoil decreases. This behavior was attributed to the presence of a centrifugal barrier along the reaction path. Here, we show, using fully quantum mechanical calculations, that for low recoil energies, the collision mechanism is nearly independent of the HD internal state and the HD product becomes aligned, with its rotational angular momentum j' pointing perpendicular to the recoil momentum k'. As the kinetic energy to overcome this barrier becomes limited, the three atoms adopt a nearly collinear configuration in the transition-state region to permit reaction, which strongly polarizes the resulting HD product. These results are expected to be general for any chemical reaction in the low recoil energy limit.

Dynamics of the reactions of muonium and deuterium atoms with vibrationally excited hydrogen molecules: tunneling and vibrational adiabaticity

Quantum mechanical (QM) and quasiclassical trajectory (QCT) calculations have been carried out for the exchange reactions of D and Mu (Mu = muonium) with hydrogen molecules in their ground and first vibrational states. In all the cases considered, the QM rate coefficients, k(T), are in very good agreement with the available experimental results. In particular, QM calculations on the most accurate potential energy surfaces (PESs) predict a rate coefficient for the Mu + H(2) (ν = 1) reaction which is very close to the preliminary estimate of its experimental value at 300 K. In contrast to the D + H(2) (ν = 0,1) and the Mu + H(2) (ν = 0) reactions, the QCT calculations for Mu + H(2) (ν = 1) predict a much smaller k(T) than that obtained with the accurate QM method. This behaviour is indicative of tunneling. The QM reaction probabilities and total reactive cross sections show that the total energy thresholds for the reactions of Mu with H(2) in ν = 0 and ν = 1 are very similar, whereas for the corresponding reaction with D the ν = 0 total energy threshold is about 0.3 eV lower than that for ν = 1. The results just mentioned can be explained by considering the vibrational adiabatic potentials along the minimum energy path. The threshold for the reaction of Mu with H(2) in both ν = 0 and ν = 1 states is the same and is given by the height of the ground vibrational adiabatic collinear potential, whereas for the D + H(2) reaction the adiabaticity is preserved and the threshold for the reaction in ν = 1 is very close to the height of the ν = 1 adiabatic collinear barrier. For Mu + H(2) (ν = 1) the reaction takes place by crossing from the ν = 1 to the ν = 0 adiabat, since the exit channel leading to MuH (ν = 1) is not energetically accessible. At the lowest possible energies, the non-adiabatic vibrational crossing implies a strong tunneling effect through the ν = 1 adiabatic barrier. Absence of tunneling in the classical calculations results in a threshold that coincides with the height of the ν = 1 adiabatic barrier. Most interestingly, the expected tunneling effect in the reaction of Mu with hydrogen molecules occurs for H(2) (ν = 1) but not for H(2) (ν = 0) where zero-point-energy effects clearly dominate.

Three-vector correlation in statistical reactions: the role of the triatomic parity

This article presents a methodology for the determination of the k-j-k' three-vector correlation assuming a statistical model for atom-diatom reactions; k and k' are the reagent-approach and product-recoil directions, respectively, and j is the rotational angular momentum of the reagent diatomic. Although the polarization of reagent angular momentum is in most cases negligible, conservation of the triatomic parity imposes a certain polarization for some combinations involving low reagent and product rotational states. Statistical and quantum-mechanical polarization-dependent differential cross sections were calculated for the barrierless D(+) + H(2)(v = 0,j) → HD(v' = 0,j') + H(+) reaction. The agreement between the two is in most cases excellent, confirming the statistical character of the reaction at low and moderate collision energies.