Cold and ultracold chemical reactions of F+HCl and F+DCl

Time-Dependent Wave Packet Dynamics Calculations of Cross Sections for Ultracold Four-Atom Reactions

We report here the first time-dependent wave packet dynamics study for an ultracold four-atom reaction. Our calculations provide accurate integral cross sections and rate constants all the way down to the Bethe-Wigner threshold regime for the benchmark OH + H₂(v = 2, j = 0) → H₂O + H reaction, indicating that the time-dependent wave packet method is a powerful tool for studying ultracold four-atom reactions.

Controlling rotational quenching rates in cold molecular collisions

The relative orientation and alignment of colliding molecules plays a key role in determining the rates of chemical processes. Here, we examine in detail a prototypical example: rotational quenching of HD in cold collisions with H₂. We show that the rotational quenching rate from j = 2 → 0, in the v = 1 vibrational level, can be maximized by aligning the HD along the collision axis and can be minimized by aligning the HD at the so called magic angle. This follows from quite general helicity considerations and suggests that quenching rates for other similar systems can also be controlled in this manner.

Restricted basis set coupled-channel calculations on atom-molecule collisions in magnetic fields

Rigorous coupled-channel quantum scattering calculations on molecular collisions in external fields are computationally demanding due to the need to account for a large number of coupled channels and multiple total angular momenta J of the collision complex. We show that by restricting the total angular momentum basis to include only the states with helicities K ≤ K_max, it is possible to obtain accurate elastic and inelastic cross sections for low-temperature He + CaH, Li + CaH, and Li + SrOH collisions in the presence of an external magnetic field at a small fraction of the computational cost of the full coupled-channel calculations (where K is the projection of the molecular rotational angular momentum on the atom-diatom axis). The optimal size of the truncated helicity basis set depends on the mechanism of the inelastic process and on the magnitude of the external magnetic field, with the minimal basis set (K_max = 0) producing quantitatively accurate results for, e.g., ultracold Li + CaH and Li + SrOH scattering at low magnetic fields, leading to nearly 90-fold gain in computational efficiency. Larger basis sets are required to accurately describe the resonance structure in the magnetic field dependence of Li + CaH and Li + SrOH inelastic cross sections in the few partial wave-regime as well as indirect spin relaxation in He + CaH collisions. Our calculations indicate that the resonance structure is due to an interplay of the spin-rotation and Coriolis couplings between the basis states of different K and the couplings between the rotational states of the same K induced by the anisotropy of the interaction potential.

Ring-polymer molecular dynamical calculations for the F + HCl → HF + Cl reaction on the ground 12A' potential energy surface

The reaction kinetics of the heavy-light-heavy abstraction reaction F + HCl → HF + Cl on the ground electronic state potential energy surface (PES) is investigated theoretically by a recently developed ring polymer molecular dynamics (RPMD) approach. First, a new PES is developed by the permutation invariant polynomial neural network (PIP-NN) approach based on 30 620 points sampled over a large configuration space from the latest and most accurate Deskevich-Hayes-Takahashi-Skodje-Nesbitt (DHTSN) PES (J. Chem. Phys., 2006, 124, 224303). Excellent fitting performance was achieved with only 521 parameters. The PIP-NN PES is 11 times faster than the DHTSN PES. Besides, the first analytical derivatives with respect to the coordinates of the atoms have been obtained for the PIP-NN PES. The RPMD rate coefficients on the PIP-NN PES are calculated and compared with available theoretical and experimental results. It is found that the experimental rate coefficients are significantly larger than the theoretical results on the DHTSN PES, due to its overestimated reaction barrier. We conclude that a reliable PES for this important heavy-light-heavy reaction is highly desirable.

Quasi-Classical Trajectory Study of Atom-Diatomic Molecule Collisions in Symmetric Hyperspherical Coordinates: The F + HCl Reaction as a Test Case

We investigate the reactive dynamics of the triatomic system F + HCl → HF + Cl for total angular momentum equal zero and for different low-lying rovibrational states of the diatomic molecule. For each of the initial vibrational quantum numbers, the time evolution of the atom-diatom collision process is investigated for a wide range of impact angles and collision energies. To this purpose, the Quasi-Classical Trajectories (QCT) method was implemented in a hyperspherical configuration space. The Hamilton equations of motion are solved numerically in an intermediate effective Cartesian space to exploit the relative simplicity of this intermediate representation. Interatomic interactions are described by a London-Eyring-Polanyi-Sato potential energy surface, specifically developed for the title reaction, and the results of the QCT simulations are discussed in terms of the time-evolution of the hyperangles. The analysis of the collision dynamics using symmetric hyperspherical coordinates provides, in addition to the description in terms of a natural reaction coordinate (the hyperradius), a more striking representation of the exchange dynamics, in terms of the time-dependent probability distribution along the kinematic rotation hyperangle, and a precise distinction between direct and indirect mechanisms of the reaction.

Product rotational angular momentum polarization of H+FCl (v=0-5; j=0, 3, 6, 9) → HF+Cl and HCl+F at Erel=0.5-20 kcal mol(-1)

The rotational angular momentum polarizations of product molecules of the title reactions on the ground potential energy surface 1 (2)A' of DHTSN [Deskevic et al. J Chem Phys 2006, 124, 224303] have been studied using the quasi-classical trajectory method. Reaction dynamic results of the HF product channel comparing with another channel of HCl with 100,000 trajectories can be accurately resolved. We show the value of the polar p(ϑr) in the range of 0° ≤ ϑr ≤ 180(°), azimuthal p(φr) in the range of 0° ≤ φr ≤ 360(°), and dihedral p(ϑr, φr) in the ranges of 0(°) ≤ ϑr ≤ 180(°) and 0(°) ≤ φr ≤ 360(°); the angular distributions of the product molecules HF and HCl at relative Erel = 0.5, 1, 2, 5, 10, 15, and 20 kcal mol(-1); and four polarization-dependent differential cross sections (PDDCSs) of HF and HCl at Erel = 0.5, 1, 2, 5, 10, and 15 kcal mol(-1). p(φr) distributions at v = 0-5, and j = 0, 3, 6, 9 at every Erel are plotted cylindrically together. The stereo dynamic transformation reaction dependent upon the rovibrational states of the reactant molecule FCl and its relative translational energies around 0.5-5 kcal mol(-1) can be significantly differentiated. Translational and rovibrational enhancements of the title reactions on both early barrier potential energy surfaces have been shown in great detail and clarified. Reaction mechanisms of forward and backward scattering of the product molecules HF and HCl, respectively, have been obtained. Graphical Abstract H + FCl → either HF + Cl (left) or HCl + F (right) is moving along a trajectory on the respective PES.

Quantum mechanical calculations of state-to-state cross sections and rate constants for the F + DCl → Cl + DF reaction

We present accurate state-to-state quantum wave packet calculations of integral cross sections and rate constants for the title reaction. Calculations are carried out on the best available ground 1(2)A' global adiabatic potential energy surface of Deskevich et al. [J. Chem. Phys. 124, 224303 (2006)]. Converged state-to-state reaction cross sections have been calculated for collision energies up to 0.5 eV and different initial rotational and vibrational excitations, DCl(v = 0, j = 0 - 1; v = 1, j = 0). Also, initial-state resolved rate constants of the title reaction have been calculated in a temperature range of 100-400 K. It is found that the initial rotational excitation of the DCl molecule does not enhance reactivity, in contract to the reaction with the isotopologue HCl in which initial rotational excitation produces an important enhancement. These differences between the isotopologue reactions are analyzed in detail and attributed to the presence of resonances for HCl(v = 0, j), absent in the case of DCl(v = 0, j). For vibrational excited DCl(v = 1, j), however, the reaction cross section increases noticeably, what is also explained by another resonance.

Time-dependent wave packet quantum scattering and quasi-classical trajectory calculations of the H + FCl(v=0,j=0) → HF + Cl/HCl + F reaction

The dynamics of the title reaction are investigated using both time-dependent wave packet quantum scattering and quasi-classical trajectory (QCT) methods on adiabatic ground 1(2)A' potential energy surface (PES). Compared with the quantum results of reaction probabilities of H + FCl(J=0) → HF + Cl/HCl + F, the QCT method is proven feasible and further employed to produce integral cross sections and rate constants. Significant resonance structures are observed in the reaction probabilities using the quantum method; however, there are some undulations in the calculated QCT integral cross sections for both product channels. A comparison between the quantum mechanical coupled-channel (CC) calculation and centrifugal sudden approximation calculation reveals the very important role of Coriolis coupling effects in the quantum calculation. Comparisons between the calculated thermal rate constants for both reactions and the previous theoretical and experimental results have been done. HCl product formation is favored over the HF product in the reactive system. Finally, the HF products are found to be mainly forward scattering, and the HCl products are mainly backward scattering.

Accurate quantum wave packet calculations for the F + HCl → Cl + HF reaction on the ground 1(2)A' potential energy surface

We present converged exact quantum wave packet calculations of reaction probabilities, integral cross sections, and thermal rate coefficients for the title reaction. Calculations have been carried out on the ground 1(2)A' global adiabatic potential energy surface of Deskevich et al. [J. Chem. Phys. 124, 224303 (2006)]. Converged wave packet reaction probabilities at selected values of the total angular momentum up to a partial wave of J = 140 with the HCl reagent initially selected in the v = 0, j = 0-16 rovibrational states have been obtained for the collision energy range from threshold up to 0.8 eV. The present calculations confirm an important enhancement of reactivity with rotational excitation of the HCl molecule. First, accurate integral cross sections and rate constants have been calculated and compared with the available experimental data.

PRODUCT ROTATIONAL ANGULAR MOMENTUM POLARIZATION IN REACTIONS D + FCl(v = 0,j = 0) → DCl + F, DF + Cl

A comparative quasi-classical trajectory study on the product alignment and orientation of reactions D + FCl (v,j) → DCl + F , DF + Cl are carried out on a recently computed 12A′ ground-state surface reported by Deskevich et al. The reaction probabilities, integral and differential cross-sections as a function of collision energy are presented. The differential cross-sections and 〈P2(j′ ⋅ k)〉 are governed by different reaction mechanisms for these two reactions. P(θr) and P(ϕr) distributions of products DCl and DF at four selected collision energies of 5, 10, 16 and 30 kcal/mol indicate that DCl and DF product molecules are not only aligned, but also oriented along y-axis. By comparing the P(θr) and P(ϕr) distributions of reactions H + FCl (v = 0,j = 0) → HCl + F and H + FCl (v = 0,j = 0) → HF + Cl , mass factor has some influence on stereodynamics of the title reactions, but the features of PES should have major influence on the difference between dynamics of reactions D + FCl (v = 0,j = 0) → DCl + F and D + FCl (v = 0,j = 0) → DF + Cl .

Product rotational angular momentum polarization in the H+FCl(v=0-5, j=0, 3, 6, 9)→HF+Cl reaction

The product alignment and orientation of the title reaction on the ground potential energy surface of 1 (2)A' have been studied using the quasi-classical trajectory method. The calculations were carried out for case (a) at collision energies of 0.5-20 kcal mol(-1) with the initially rovibrational state of the reagent FCl molecule being at the v = 0 and j = 0 level to especially reveal in detail the dependence of the product integral cross section on collision energy. Further calculations at the collision energy of 15 kcal mol(-1) for case (b) at v = 0-5, and j = 0, and (c) at v = 0, and j = 3, 6, 9 initial states were carried out to reveal the effect of initially vibrational and rotational excitations on stereodynamics, respectively. Possessing final relative velocity k' (defined as a vector in the xz-plane), product alignment perpendicular to the reagent relative velocity vector k (defined as z- or parallel to the z-axis), for case (a) is found to be weaker at all collision energies, for case (b) is found to be vibrationally enhanced by the reactant molecule FCl, but for case (c), rather insensitive to initially rotational excitation. The rotational vector of product molecular orientation pointing to either negative or positive direction of the y-axis in the center of mass frame, e.g. origin of the coordinate system, is enhanced by collision energies regarding to 0.5-20 kcal mol(-1), while it becomes weaker at higher vibrational (v = 0-5) or rotational (j = 0, 3, 6, 9) excitation levels. Effects of collision energies and of rotational excitation at these collision energies, with 15 kcal mol(-1) as an example on the calculated PDDCSs are also shown and discussed. Detailed plots P(φ(r)) in the range of 0 ≤φ(r)≤ 360(o), and P(θ(r), φ(r)) in the ranges of 0 ≤θ(r)≤ 180° and 0 ≤φ(r)≤ 360° at collision energies 0.5-20 kcal mol(-1) have been presented. Overall, results of PDDCSs of the product alignment and product orientation at these collision energies in the title reaction are not very strongly distinguishable.