Dynamics of the HCl + C2H5 Multichannel Reaction on a Full-Dimensional Ab Initio Potential Energy Surface

We report a full-dimensional ab initio analytical potential energy surface (PES), which accurately describes the HCl + C2H5 multichannel reaction. The new PES is developed by iteratively adding selected configurations along HCl + C2H5 quasi-classical trajectories (QCTs), thereby improving our previous Cl(2P3/2) + C2H6 PES using the Robosurfer program package. QCT simulations for the H'Cl + C2H5 reaction reveal hydrogen-abstraction, chlorine-abstraction, and hydrogen-exchange channels leading to Cl + C2H5H', H' + C2H5Cl, and HCl + C2H4H', respectively. Hydrogen abstraction dominates in the collision energy (Ecoll) range of 1-80 kcal/mol and proceeds with indirect isotropic scattering at low Ecoll and forward-scattered direct stripping at high Ecoll. Chlorine abstraction opens around 40 kcal/mol collision energy and becomes competitive with hydrogen abstraction at Ecoll = 80 kcal/mol. A restricted opening of the cone of acceptance in the Cl-abstraction reaction is found to result in the preference for a backward-scattering direct-rebound mechanism at all energies studied. Initial attack-angle distributions show mainly side-on collision preference of C2H5 for both abstraction reactions, and in the case of the HCl reactant, H/Cl-side preference for the H/Cl abstraction. For hydrogen abstraction, the collision energy transfer into the product translational and internal energy is almost equally significant, whereas in the case of chlorine abstraction, most of the available energy goes into the internal degrees of freedom. Hydrogen exchange is a minor channel with nearly constant reactivity in the Ecoll range of 10-80 kcal/mol.

Vibrational Mode-Specific Dynamics of the OH + C2H6 Reaction

We investigate the effects of the initial vibrational excitations on the dynamics of the OH + C2H6 → H2O + C2H5 reaction using the quasi-classical trajectory method and a full-dimensional analytical ab initio potential energy surface. Excitation of the initial CH, CC, and OH stretching modes enhances, slightly inhibits, and does not affect the reactivity, respectively. Translational energy activates the early-barrier title reaction more efficiently than OH and CC stretching excitations, in accord with the Polanyi rules whereas CH stretching modes have similar or higher efficacy than translation, showing that these rules are not always valid in polyatomic processes. Scattering angle, initial attack angle, and product translational energy distributions show the dominance of direct stripping with increasing collision energy, side-on OH and isotropic C2H6 attack preferences, and substantial reactant-product translational energy transfer without any significant mode specificity. The reactant vibrational excitation energy of OH and C2H6 flows into the H2O and C2H5 product vibrations, respectively, whereas product rotations are not affected. The computed mode-specific H2O vibrational distributions show that initial OH excitation appears in the asymmetric stretching vibration of the H2O product and allow comparison with experiments.

Full-dimensional automated potential energy surface development and dynamics for the OH + C2H6 reaction

We develop a full-dimensional analytical potential energy surface (PES) for the OH + C2H6 reaction using the Robosurfer program system, which automatically (1) selects geometries from quasi-classical trajectories, (2) performs ab initio computations using a CCSD(T)-F12/triple-zeta-quality composite method, (3) fits the energies utilizing the permutationally-invariant monomial symmetrization approach, and iteratively improves the PES via steps (1)−(3). Quasi-classical trajectory simulations on the new PES reveal that hydrogen abstraction leading to H2O + C2H5 dominates in the collision energy range of 10−50 kcal/mol. The abstraction cross sections increase and the dominant mechanism shifts from rebound (small impact parameters and backward scattering) to stripping (larger impact parameters and forward scattering) with increasing collision energy as opacity functions and scattering angle distributions indicate. The abstraction reaction clearly favors side-on OH attack over O-side and the least-preferred H-side approach, whereas C2H6 behaves like a spherical object with only slight C−C-perpendicular side-on preference. Collision energy efficiently flows into the relative translation of the products, whereas product internal energy distributions show only little collision energy dependence. H2O/C2H5 vibrational distributions slightly/significantly violate zero-point energy and are nearly independent of collision energy, whereas the rotational distributions clearly blue-shift as collision energy increases.

Benchmark Ab Initio Characterization of the Abstraction and Substitution Pathways of the Cl + CH3CN Reaction

We investigate the reaction pathways of the Cl + CH₃CN system: hydrogen abstraction, methyl substitution, hydrogen substitution, and cyanide substitution, leading to HCl + CH₂CN, ClCN/CNCl + CH₃, ClCH₂CN + H, and CH₃Cl + CN, respectively. Hydrogen abstraction is exothermic and has a low barrier, whereas the other channels are endothermic with high barriers. The latter two can proceed via a Walden inversion or front-side attack mechanism, and the front-side attack barriers are always higher. The C-side methyl substitution has a lower barrier and also a lower endothermicity than the N-side reaction. The computations utilize an accurate composite ab initio approach and the explicitly correlated CCSD(T)-F12b method. The benchmark classical and vibrationally adiabatic energies of the stationary points are determined with the most accurate CCSD(T)-F12b/aug-cc-pVQZ energies adding further contributions of the post-(T) and core correlation, scalar relativistic effects, spin–orbit coupling, and zero-point energy corrections. These contributions are found to be non-negligible to reach subchemical accuracy.

Rotational Mode-Specificity in the Cl + C2H6 → HCl + C2H5 Reaction

We perform rotational mode-specific quasi-classical trajectory simulations using a high-quality ab initio analytical potential energy surface for the Cl(²P_3/2) + C₂H₆ → HCl + C₂H₅ reaction. As ethane, being a prolate-type symmetric top, can be characterized by the J and K rotational quantum numbers, the excitation of two rotational modes, the tumbling (J, K = 0) and spinning (J, K = J) rotations of the reactant is carried out with J = 10, 20, 30, and 40 at a wide range of collision energies. The impacts of rotational excitation on the reactivity, the mechanism, and the post-reaction distribution of energy are investigated: (1) exciting both rotational modes enhances the reactivity with the spinning rotation being more effective due to its coupling to the C-H stretching vibrational normal modes (C-H bond elongating effect) and larger rotational energies, (2) rotational excitation increases the dominance of direct rebound over the stripping mechanism, while collision energy favors the latter, (3) investing energy in tumbling rotation excites the translational motion of the products, while the excess spinning rotational energy readily flows into the internal degrees of freedom of the ethyl radical or, less significantly, into the HCl vibration, probably due to the pronounced rovibrational coupling in this case. We also study the relative efficiency of vibrational and rotational excitation on the reactivity of the barrierless and thus translationally hindered title reaction.

Feshbach resonances in D + HD(v = 1, j = 0) reaction at low collision energies

Feshbach resonances in D + HD(v = 1, j = 0) reaction are studied by using the time-independent quantum method. The integral cross section (ICS) results present three Feshbach resonance peaks, which are different from H + HD(v = 1, j = 0) reaction dominated by only one peak. These resonances are attributed to coupling with adiabatic effective potentials of D + HD(v = 1, j = 1) reaction, and the most obvious peak is contributed by J = 1 at 83.16 cm^-1 collision energy. For J = 0 and 2, the resonances are related with the same L partial wave and present a double-peak structure in total ICS. The characteristics of product angular distribution show that the resonance of J = 1 is long-lived, while the lifetimes are relatively shorter for the resonance of J = 0 and 2.

Time-dependent quantum dynamics study of the F + C2H6 → HF + C2H5 reaction

The reaction probabilities, integral cross sections, energy efficiency and rate constants are investigated for the F + C₂H₆ reaction using the quantum reaction dynamics, wave packet method. The ground-state integral cross section calculated using a six-degree-of-freedom approach is in very good agreement with the quasi-classical trajectory results. We find that the H-CH₂CH₃ stretching motion has the largest enhancement to reactivity, followed by the H-CH₂-CH₃ bending motion. However, the stretching motion between CH₂ and CH₃ slightly hinders the reactivity. The energy-form efficacy based on an equal amount of total energy shows that translational energy is more effective in enhancing the reactivity than vibrational energy of the H-CH₂CH₃ stretching motion at a relatively lower translational energy, while the reverse is true at a relatively high translational energy. An energy-shifting method is employed to calculate the full-dimensional rate constants. The quantum rate constants agree well with one of the two main experimental measurements, and the activation energy has an excellent agreement with the one calculated using canonical variational transition-state theory.

Using quantum dynamics to study the effect of energy efficiency on the reactivity of the OH + DBr reaction

We report a time-dependent, full dimensional, wave-packet calculation for the reaction of OH + DBr to examine the effect of the energy efficiency on the reactivity. This study shows that the vibrational excitations of the OH and DBr enhance the reaction. However, the rotational excitations of OH and DBr both hinder the reaction. As a result, the vibrational energies of both the OH and DBr reactants are more efficient at promoting the reactivity than the translational energy, while the rotational energies of OH and DBr are less effective than the translational energy. By analyzing the state population of the vibrational and rotational states along the reaction pathway, we also developed an approach in order to explain the enhancement of the vibrational excitation and the hindrance of the rotational excitation of the reaction. We found that the initial-state selected vibrational excited states of OH and DBr are the dominant components, respectively, for surmounting the barrier. However, the initial-state selected rotational excited states of OH and DBr are no longer the dominant states for surmounting the transition state owing to their population changes in the van der Waals well. This quantitative analysis demonstrates the potential well in the entrance valley plays an important role in the energy efficiency with regards to the reactivity.

Vibrational mode-specific dynamics of the F(2P3/2) + C2H6 → HF + C2H5 reaction

We investigate the competing effect of vibrational and translational excitation and the validity of the Polanyi rules in the early- and negative-barrier F(²P_3/2) + C₂H₆ → HF + C₂H₅ reaction by performing quasi-classical dynamics simulations on a recently developed full-dimensional multi-reference analytical potential energy surface. The effect of five normal-mode excitations of ethane on the reactivity, the mechanism, and the post-reaction energy flow is followed through a wide range of collision energies. Promoting effects of vibrational excitations and interaction time, related to the slightly submerged barrier, are found to be suppressed by the early-barrier-induced translational enhancement, in contrast to the slightly late-barrier Cl + C₂H₆ reaction. The excess vibrational energy mostly converts into ethyl internal excitation while collision energy is transformed into product separation. The substantial reaction energy excites the HF vibration, which tends to show mode-specificity and translational energy dependence as well. With increasing collision energy, direct stripping becomes dominant over the direct rebound and indirect mechanisms, being basically independent of reactant excitation.

Vibrational mode-specificity in the dynamics of the Cl + C2H6 → HCl + C2H5 reaction

We report a detailed dynamics study on the mode-specificity of the Cl + C₂H₆ → HCl + C₂H₅ H-abstraction reaction. We perform quasi-classical trajectory simulations using a recently developed high-level ab initio full-dimensional potential energy surface by exciting five different vibrational modes of ethane at four collision energies. We find that all the studied vibrational excitations, except that of the CC-stretching mode, clearly promote the title reaction, and the vibrational enhancements are consistent with the predictions of the Sudden Vector Projection (SVP) model, with the largest effect caused by the CH-stretching excitations. Intramolecular vibrational redistribution is also monitored for the differently excited ethane molecule. Our results indicate that the mechanism of the reaction changes with increasing collision energy, with no mode-specificity at high energies. The initial translational energy mostly converts into product recoil, while a significant part of the excess vibrational energy remains in the ethyl radical. An interesting competition between translational and vibrational energies is observed for the HCl vibrational distribution: the effect of exciting the low-frequency ethane modes, having small SVP values, is suppressed by translational excitation, whereas a part of the excess vibrational energy pumped into the CH-stretching modes (larger SVP values) efficiently flows into the HCl vibration.

Benchmark ab initio stationary-point characterization of the complex potential energy surface of the multi-channel Cl + CH3NH2 reaction

We characterize the exothermic low/submerged-barrier hydrogen-abstraction (HCl + CH2NH2/CH3NH) as well as, for the first time, the endothermic high-barrier amino-substitution (CH3Cl + NH2), methyl-substitution (NH2Cl + CH3), and hydrogen-substitution (CH2ClNH2/CH3NHCl + H) pathways of the Cl + CH3NH2 reaction using an accurate composite ab initio approach. The computations reveal a CH3NH2Cl complex in the entrance channel, nine transition states corresponding to different abstractions, Walden-inversion substitution, and configuration-retaining front-side attack substitution pathways, as well as nine post-reaction complexes. The global minima of the electronic and vibrationally adiabatic potential energy surfaces correspond to the pre-reaction CH3NH2Cl and post-reaction CH2NH2HCl complexes, respectively. The benchmark composite energies of the stationary points are obtained by considering basis-set effects up to the correlation-consistent polarized valence quadruple-zeta basis augmented with diffuse functions (aug-cc-pVQZ) using the explicitly-correlated coupled-cluster singles, doubles, and perturbative triples CCSD(T)-F12b method, post-(T) correlation up to CCSDT(Q) including full triples and perturbative quadruples, core correlation, and scalar relativistic and spin-orbit effects, as well as harmonic zero-point energy corrections.

Rotational Mode Specificity in the F- + CH3I(v = 0, JK) SN2 and Proton-Transfer Reactions

Quasiclassical trajectory computations are performed for the F^– + CH₃I(v = 0, JK) → I^– + CH₃F (S_N2) and HF + CH₂I^– (proton-transfer) reactions considering initial rotational states characterized by J = {0, 2, 4, 6, 8, 12, and 16} and K = {0 and J} in the 1–30 kcal/mol collision energy (E_coll) range. Tumbling rotation (K = 0) counteracts orientation effects, thereby hindering the S_N2 reactivity by about 15% for J = 16 in the 1–15 kcal/mol E_coll range and has a negligible effect on proton transfer. Spinning about the C–I bond (K = J), which is 21 times faster than tumbling, makes the reactions more direct, inhibiting the S_N2 reactivity by 25% in some cases, whereas significantly enhancing the proton-transfer channel by a factor of 2 at E_coll = 15 kcal/mol due to the fact that the spinning-induced centrifugal force hinders complex formation by breaking H-bonds and activates C–H bond cleavage, thereby promoting proton abstraction on the expense of substitution. At higher E_coll, as the reactions become more direct, the rotational effects are diminishing.

Reactivity oscillation in the heavy-light-heavy Cl + CH4 reaction

It has long been predicted that oscillatory behavior exists in reactivity as a function of collision energy for heavy-light-heavy (HLH) chemical reactions in which a light atom is transferred between two heavy atoms or groups of atoms, but direct observation of such a behavior in bimolecular reactions remains a challenge. Here we report a joint theoretical and crossed-molecular-beam study on the Cl + CH₄ → HCl + CH₃ reaction. A distinctive peak at a collision energy of 0.15 eV for the CH₃(v = 0) product was experimentally detected in the backward scattering direction. Detailed quantum-dynamics calculations on a highly accurate potential energy surface revealed that this feature originates from the reactivity oscillation in this HLH polyatomic reaction. We anticipate that such reactivity oscillations exist in many HLH reactions involving polyatomic reagents.

Benchmark ab initio characterization of the abstraction and substitution pathways of the OH + CH4/C2H6 reactions

We report benchmark ab initio stationary-point properties for the hydrogen-abstraction, hydrogen-substitution, and methyl-substitution pathways of the OH + CH₄/C₂H₆ reactions.

An accurate potential energy surface and ring polymer molecular dynamics study of the Cl + CH4→ HCl + CH3reaction

Thermal rate coefficients for the Cl + CH₄/CD₄reactions were studied on a new full-dimensional accurate potential energy surface with the spin–orbit corrections considered in the entrance channel.

Detailed benchmark ab initio mapping of the potential energy surfaces of the X + C2H6 [X = F, Cl, Br, I] reactions

We provide benchmark relative energies for the stationary points of three different channels of the halogen atom + ethane reactions.

Mode-Specific Quasiclassical Dynamics of the F- + CH3I SN2 and Proton-Transfer Reactions

Mode-specific quasiclassical trajectory computations are performed for the F^- + CH₃I( v _k = 0, 1) S_N2 and proton-transfer reactions at nine different collision energies in the range of 1.0-35.3 kcal/mol using a full-dimensional high-level ab initio analytical potential energy surface with ground-state and excited CI stretching ( v₃), CH₃ rocking ( v₆), CH₃ umbrella ( v₂), CH₃ deformation ( v₅), CH symmetric stretching ( v₁), and CH asymmetric stretching ( v₄) initial vibrational modes. Millions of trajectories provide statistically definitive mode-specific cross sections, opacity functions, scattering angle distributions, and product internal energy distributions. The excitation functions reveal slight vibrational S_N2 inversion inhibition/enhancement at low/high collision energies ( E_coll), whereas large decaying-with- E_coll vibrational enhancement effects for the S_N2 retention (double inversion) and proton-transfer channels. The most efficient vibrational enhancement is found by exciting the CI stretching (high E_coll) for S_N2 inversion and the CH stretching modes (low E_coll) for double inversion and proton transfer. Mode-specific effects do not show up in the scattering angle distributions and do blue-shift the hot/cold S_N2/proton-transfer product internal energies.

Understanding mode-specific dynamics in the local mode representation

Mode specificity is a main characteristic of transition state control of reaction dynamics. The normal mode representation has been widely employed to describe the mode specificity in elementary chemical reactions. However, spectroscopists have demonstrated that the local mode representation has advantages in analyzing the overtone and combination band spectra. In this work, the mode-specific reaction dynamics between the hydrogen atom and the molecules H2S and H2O is studied using a full-dimensional quantum scattering model in the (2 + 1) Radau-Jacobi coordinates. The mode specificities in the reactions that violates our physical intuition in the normal mode representation are well rationalized in the local mode representation. The energy flow between different XH bonds resulting from the intramolecular interaction and/or intermolecular interaction is unveiled, together with its impacts on dynamics of the abstraction and exchange reactions.

Angular Scattering Dynamics of the CH4 + Cl → CH3 + HCl Reaction Using Nearside-Farside, Local Angular Momentum, and Resummation Theories

The differential cross section (DCS) for the CH4 + Cl → CH3 + HCl reaction is studied at six total energies where all of the species are in their ground states. The scattering (S) matrix elements have been calculated by the rotating line umbrella method for a dual-level ab initio analytic potential energy surface. We make the first application to this reaction of nearside-farside (NF) and local angular momentum (LAM) techniques, including resummation orders (r) of 0, 1, 2, and 3 for the partial-wave series representation of the full scattering amplitude. We find that resummation usually cleans the NF r = 0 DCSs of unphysical oscillations, especially at small angles. This cleaning effect is typically most pronounced when changing from no resummation (r = 0) to r = 1; further resummations from r = 1 to r = 2 and from r = 2 to r = 3 have smaller effects. The NF DCS analyses show that the reaction is N-dominated at sideward and large angles, whereas at small angles there are oscillations caused by NF interference. The NF LAM analysis provides consistent and complementary information, in particular for the total angular momenta that contribute to the reaction at different scattering angles. The NF analyses also provide justification for simpler N-dominant dynamical theories such as the semiclassical optical model, which provides an explanation for the distorted mirror image effect for the moduli of the S matrix elements and the DCSs, as well as the use of a hard-sphere DCS over limited angular ranges.

Energy efficiency in surmounting the central energy barrier: a quantum dynamics study of the OH + CH3 → O + CH4 reaction

A quantum dynamics study of the OH + CH₃ with a slightly early barrier shows that vibrational energy is more effective in promoting the reactivity than translational energy, which is just opposite to the Polanyi rules.

Quantum rate coefficients and kinetic isotope effect for the reaction Cl + CH4 → HCl + CH3 from ring polymer molecular dynamics

Thermal rate coefficients and kinetic isotope effect have been calculated for prototypical heavy-light-heavy polyatomic bimolecular reactions Cl + CH4/CD4 → HCl/DCl + CH3/CD3, using a recently proposed quantum dynamics approach: ring polymer molecular dynamics (RPMD). Agreement with experimental rate coefficients, which are quite scattered, is satisfactory. However, differences up to 50% have been found between the RPMD results and those obtained from the harmonic variational transition-state theory on one of the two full-dimensional potential energy surfaces used in the calculations. Possible reasons for such discrepancy are discussed. The present work is an important step in a series of benchmark studies aimed at assessing accuracy for RPMD for chemical reaction rates, which demonstrates that this novel method is a quite reliable alternative to previously developed techniques based on transition-state theory.