Quasi-classical treatment of the Stereodynamics of chemical reactions: k-r-k′ vector correlation for the Li+HF(v=1,j=1)→LiF+H reaction

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.

Dynamical and Mechanical Insights into the Li(2 S)+ HCl( X 1 Σ + ${X^1 {\rm{\Sigma }}^ + }$ ) Reaction: A Detailed Quantum Wavepacket Study

Quantum wave packet dynamics of the Li(² S)+HCl( X 1 Σ + ${X^1 \Sigma ^ + }$ ) reaction in its electronic ground state is studied. The initial state-selected and energy-resolved dynamical attributes such as reaction probability, integral cross section, and thermal rate constant for the Cl-abstraction and H-abstraction pathways are reported. All partial wave contributions of J up to 120 were found to be necessary for the title reaction up to the collision energy of ∼1.0 eV. The dynamical results reveal that the Cl-abstraction is more favored over the H-abstraction for the different rovibrational (v, j) excitations. Due to the existence of an early barrier in the potential energy surface, the cross sections increase with increasing collision energy. The rate constants also monotonously increase with temperature for both channels. Resonances are identified and characterized in terms of eigenfunctions and lifetimes. Nearly 120 well-resolved eigenstates are reported for the LiHCl complex, and they are categorized as van der Waals (vdW), barrier and product states according to the nodal progressions along (R, r, γ). The vdW resonances reveal a local-mode behavior of quasibound type at low energies and extended progressions at high energies. Further, the single-quantized periodic orbit type is also observed in the barrier region, which decays very fast. Finally, the lifetime analysis reveals that the vdW resonances can survive as long as ∼2.2 ps, which is much longer than the lifetime of the resonances in the barrier region.

Controlling the Spin-Orbit Branching Fraction in Molecular Collisions

The collision geometry, that is, the relative orientation of reactants before interaction, can have a large effect on how a collision or reaction proceeds. Certain geometries may prevent access to a given product channel, while others might enhance it. In this Letter, we demonstrate how the initial orientation of NO molecules relative to approaching Ar atoms determines the branching between the spin-orbit changing and the spin-orbit conserving rotational product channels. We use a recently developed quantum treatment to calculate differential and integral branching fractions, at any arbitrary orientation, from theoretical and experimental data points. Our results show that a substantial degree of control over the final spin-orbit state of the scattering products can be achieved by tuning the initial collision geometry.

Steric Effects in the Inelastic Scattering of NO(X) + Ar: Side-on Orientation

The rotationally inelastic collisions of NO(X) with Ar, in which the NO bond-axis is oriented side-on (i.e., perpendicular) to the incoming collision partner, are investigated experimentally and theoretically. The NO(X) molecules are selected in the |j = 0.5, Ω = 0.5, ε = -1, f⟩ state prior to bond-axis orientation in a static electric field. The scattered NO products are then state selectively detected using velocity-map ion imaging. The experimental bond-axis orientation resolved differential cross sections and integral steric asymmetries are compared with quantum mechanical calculations, and are shown to be in good agreement. The strength of the orientation field is shown to affect the structure observed in the differential cross sections, and to some extent also the steric preference, depending on the ratio of the initial e and f Λ-doublets in the superposition determined by the orientation field. Classical and quantum calculations are compared and used to rationalize the structures observed in the differential cross sections. It is found that these structures are due to quantum mechanical interference effects, which differ for the two possible orientations of the NO molecule due to the anisotropy of the potential energy surface probed in the side-on orientation. Side-on collisions are shown to maximize and afford a high degree of control over the scattering intensity at small scattering angles (θ < 90°), while end-on collisions are predicted to dominate in the backward scattered region (θ > 90°).

A new ab initio potential energy surface of LiClH (1A') system and quantum dynamics calculation for Li + HCl (v = 0, j = 0-2) → LiCl + H reaction

A new ab initio potential energy surface (PES) for the ground state of Li + HCl reactive system has been constructed by three-dimensional cubic spline interpolation of 36 654 ab initio points computed at the MRCI+Q/aug-cc-pV5Z level of theory. The title reaction is found to be exothermic by 5.63 kcal/mol (9 kcal/mol with zero point energy corrections), which is very close to the experimental data. The barrier height, which is 2.99 kcal/mol (0.93 kcal/mol for the vibrationally adiabatic barrier height), and the depth of van der Waals minimum located near the entrance channel are also in excellent agreement with the experimental findings. This study also identified two more van der Waals minima. The integral cross sections, rate constants, and their dependence on initial rotational states are calculated using an exact quantum wave packet method on the new PES. They are also in excellent agreement with the experimental measurements.

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.

A new potential energy surface of LiHCl system and dynamic studies for the Li(2S) + HCl(X1Σ+) → LiCl(X1Σ+) + H(2S) reaction

A new global potential energy surface (PES) is constructed for the ground state of LiHCl system based on high-quality ab initio energy points calculated using multi-reference configuration interaction calculations with the Davidson correction. The AVQZ and WCVQZ basis sets are employed for H and Li atoms, respectively. To compensate the relativistic effects of heavy element, the AWCVQZ-DK basis set is employed for Cl atom. The neural network method is used for fitting the PES, and the root mean square error is small (1.36 × 10^-2 eV). The spectroscopic constants of the diatoms obtained from the new PES agree well with experimental data. The geometric characteristics of the transition state and the complex are examined and compared with the previous theoretical values. To study the reaction dynamics of the Li(²S) + HCl(X¹Σ⁺) → LiCl(X¹Σ⁺) + H(²S) reaction, quantum reactive scattering dynamics calculations using collection reactant-coordinate-based wave packet method are conducted based on the new PES. The results of the reaction probabilities indicate that a small barrier exists along the reaction path as observed from the PES. The integral cross section curves reveal that the product molecule LiCl is easily excited. In addition, the reaction is dominated by forward scattering, and similar pattern is observed from Becker's experiment.

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.

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.

EFFECTS ON THE CHEMICAL STEREODYNAMICS OF THE INITIAL VIBRATIONAL EXCITATION IN THE F + LiH (v = 0-2, j = 0) → LiF + H REACTION

In this work, the product polarization characteristics are reported for the reaction F + LiH (v = 0 - 2, j = 0) → LiF + H at the collision energy of 35 kcal/mol, by using the QCT method on the Aguado–Paniagua-potential energy surface (see Aguado et al.). The distribution of P(θr) which represents the K (reagent relative velocity vector) and J′ (product rotational angular momentum vector) correlation, the dihedral angle distribution of K-K′ (product relative velocity vector)-J′ P(φr), the angular distribution P(θr, φr) and the four polarization-dependent differential cross sections (2π/σ)(dσ00/dωt), (2π/σ)(dσ20/dωt), (2π/σ)(dσ22+/dωt), (2π/σ)(dσ21-/dωt) in each initial state are presented and discussed.

INFLUENCE OF THE COLLISION ENERGY ON STEREODYNAMICS OF THE F + LiH (v = 0, j = 0) → LiF + H REACTION

The stereodynamics of the title reaction based on the ground 2A′ potential energy surface (PES) has been investigated using the method of the quasi-classical trajectory (QCT) at different collision energies (23 kcal/mol, 35 kcal/mol and 46 kcal/mol). The vector properties of the angular momentum (described by the distribution of K - J′P(θr), the dihedral angle distribution of K - K′ - J′P(φr) and the angular distribution P(θr, ϕr)) and the four PDDCSs [(2π/σ)(dσ00/dωt), (2π/σ)(dσ20/dωt), (2π/σ)(dσ22+/dωt), (2π/σ)(dσ21-/dωt)] of the product LiF at each collision energy have been presented, respectively. Further, the collision energy effects on the behavior of the product LiF have been discussed and studied.

Stereodynamics study of the O+ + D2 (υ = 0, j = 0–2) → OD+ + D reaction

The vector correlations between products and reagents for the ion–molecule reaction O+ + D2 → OD+ + D with different rotational quantum numbers (j = 0, 1, or 2) were explored theoretically using the quasi-classical trajectory method (QCT) on a Martìnez–Millán–González (MMG) surface. The three angular distributions P(θr), P([Formula: see text]), and P(θr,[Formula: see text]), as well as four polarization-dependent differential cross sections (PDDCSs) were calculated. The results indicate that a reagent’s rotational excitation greatly influences both the vector correlations of k–k′, k–j′, and k–k′–j′ and the PDDCSs of the title reaction, which means the reactivity is very sensitive to the rotational quantum number.

A quasi-classical trajectory (QCT) study of the H + OF reaction stereodynamics

Based on the global three-dimensional adiabatic potential surface of the 13A′ excited electronic state (J. Chem. Phys. 2005, 123, 114310) of the OHF system, we investigated the H + OF → OH + F/HF + O reaction stereodynamics by using the quasi-classical trajectory (QCT) method. The four polarization-dependent differential cross sections (PDDCSs) and the three angular distributions P(θr), P([Formula: see text]), P(θr, [Formula: see text]) were calculated at a low collision energy of 0.48 eV for both product channels. The results indicated that the products are backward-scattering on the triplet state, and the product rotational angular momenta are aligned or oriented. Moreover, there is a remarkable difference between the polarization behaviors of the two product channels. Product orientation exhibited by the OH + F channel is found to be absent in the HF + O channel at this collision energy albeit the latter shows stronger alignment along the direction perpendicular to the reagent relative velocity k than OH + F.

How Reactants Polarization Can Be Used to Change and Unravel Chemical Reactivity

This article presents theoretical methods for the description of the directional effect of reactant rotation on the reactivity of atom-diatom systems and suggests an experiment that could be used to test theoretical predictions. The theory can be used in conjunction with both quantum reactive scattering and quasiclassical trajectory calculations, and is stated in general terms, which allows it to deal with arbitrary reactant polarizations. The illustrative results obtained for the benchmark H + D2 reaction are also presented and show that under experimentally achievable conditions one can largely control reactive cross sections and product state distributions, while at the same time gaining valuable and at times surprising information on the reaction mechanism.

Heavy atom tunneling in chemical reactions: Study of H+LiF collisions

The H+LiF(X (1)sigma(+),upsilon=0-2,j=0)-->HF(X (1)sigma(+),upsilon',j')+Li(2S) bimolecular process is investigated by means of quantum scattering calculations on the chemically accurate X 2A' LiHF potential energy surface of Aguado et al. [A. Aguado, M. Paniagua, C. Sanz, and J. Roncero, J. Chem. Phys. 119, 10088 (2003)]. Calculations have been performed for zero total angular momentum for translational energies from 10(-7) to 10(-1) eV. Initial-state selected reaction probabilities and cross sections are characterized by resonances originating from the decay of metastable states of the H...F-Li and Li...F-H van der Waals complexes. Extensive assignment of the resonances has been carried out by performing quasibound states calculations in the entrance and exit channel wells. Chemical reactivity is found to be significantly enhanced by vibrational excitation at low temperatures, although reactivity appears much less favorable than nonreactive processes due to the inefficient tunneling of the relatively heavy fluorine atom strongly bound in van der Waals complexes.

Quantum dynamics of the Li+HF→H+LiF reaction at ultralow temperatures

Quantum-mechanical calculations are reported for the Li+HF(v=0,1,j=0)-->H+LiF(v',j') bimolecular scattering process at low and ultralow temperatures. Calculations have been performed for zero total angular momentum using a recent high-accuracy potential-energy surface for the X2A' electronic ground state. For Li+HF(v=0,j=0), the reaction is dominated by resonances due to the decay of metastable states of the Li cdots,...F-H van der Waals complex. Assignment of these resonances has been carried out by calculating the eigenenergies of the quasibound states. We also find that while chemical reactivity is greatly enhanced by vibrational excitation, the resonances get mostly washed out in the reaction of vibrationally excited HF with Li atoms. In addition, we find that at low energies, the reaction is significantly suppressed due to the less-efficient tunneling of the relatively heavy fluorine atom.