Communication: Mode specific quantum dynamics of the F + CHD3 → HF + CD3 reaction

An efficient and universal parallel algorithm for high-dimensional quantum dynamics in poly-atomic reactions

This study presents a parallel algorithm for high-dimensional quantum dynamics simulations in poly atomic reactions, integrating distributed- and shared-memory models. The distributions of the wave function and potential energy matrix across message passing interface processes are based on bundled radial and angular dimensions, with implementations featuring either two- or one-sided communication schemes. Using realistic parameters for the H + NH3 reaction, performance assessment reveals linear scalability, exceeding 90% efficiency with up to 600 processors. In addition, owing to the universal and concise structure, the algorithm demonstrates remarkable extensibility to diverse reaction systems, as demonstrated by successes with six-atom and four-atom reactions. This work establishes a robust foundation for high-dimensional dynamics studies, showcasing the algorithm's efficiency, scalability, and adaptability. The algorithm's potential as a valuable tool for unraveling quantum dynamics complexities is underscored, paving the way for future advancements in the field.

Theoretical dynamics studies of the CH3 + HBr → CH4 + Br reaction: effects of isotope substitution and vibrational excitation of CH3

The rate coefficient for two deuterium substituted isotopologues of reaction CH3 + HBr → CH4 + Br has been determined using the quasiclassical trajectory (QCT) method. We used the analytical potential energy surface (PES) fitted to high-level ab initio points in earlier work. The PES exhibits a pre-reaction van der Waals complex and a submerged potential barrier. The rate coefficients of the deuterated isotopologue reactions, similarly to the pure-protium isotopologue, show significant deviation from the Arrhenius law, namely, the activation energy is negative below about 600 K and positive above it: k[CH3 + DBr] = 1.35 × 10-11 exp(- 2472/T) + 5.85 × 10-13 exp(335/T) and k[CD3 + HBr] = 2.73 × 10-11 exp(- 2739/T) + 1.46 × 10-12 exp(363/T). The CH3 + DBr reaction is slower by a factor of 1.8, whereas CD3 + HBr isotopologue is faster by a factor of 1.4 compared to the HBr + CH3 system across a wide temperature range. The isotope effects are interpreted in terms of the properties of various regions of the PES. Quantum state-resolved simulations revealed that the reaction of CH3 with HBr becomes slower when any of the vibrational modes of the methyl radical is excited. This contradicts the assumption that vibrational excitation of methyl radicals enhances its reactivity, which is of historical importance: this assumption was used as an argument against the existence of negative activation energy in a decade-long controversy in the 1980s and 1990s.

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.

Feshbach resonances in the F + CHD3 → HF + CD3 reaction

The signature of dynamics resonances was observed in the benchmark polyatomic F + CH₄/CHD₃ reactions more than a decade ago; however, the dynamical origin of the resonances is still not clear due to the lack of reliable quantum dynamics studies on accurate potential energy surfaces. Here, we report a six-dimensional state-to-state quantum dynamics study on the F + CHD₃ → HF + CD₃ reaction on a highly accurate potential energy surface. Pronounced oscillatory structures are observed in the total and product rovibrational-state-resolved reaction probabilities. Detailed analysis reveals that these oscillating features originate from the Feshbach resonance states trapped in the peculiar well on the HF(v' = 3)-CD₃ vibrationally adiabatic potential caused by HF chemical bond softening. Most of the resonance structures on the reaction probabilities are washed out in the well converged integral cross sections (ICS), leaving only one distinct peak at low collision energy. The calculated HF vibrational state-resolved ICS for CD₃(v = 0) agrees quantitatively with the experimental results, especially the branching ratio, but the theoretical CD₃ umbrella vibration state distribution is found to be much hotter than the experiment.

Vibrational mode-specific dynamics of the F- + CH3CH2Cl multi-channel reaction

We investigate the mode-specific dynamics of the ground-state, C-Cl stretching (v₁₀), CH₂ wagging (v₇), sym-CH₂ stretching (v₁), and sym-CH₃ stretching (v₃) excited F^- + CH₃CH₂Cl(v_k = 0, 1) [k = 10, 7, 1, 3] → Cl^- + CH₃CH₂F (S_N2), HF + CH₃CHCl^-, FH⋯Cl^- + C₂H₄, and Cl^- + HF + C₂H₄ (E2) reactions using a full-dimensional high-level analytical global potential energy surface and the quasi-classical trajectory method. Excitation of the C-Cl stretching, CH₂ stretching, and CH₂/CH₃ stretching modes enhances the S_N2, proton abstraction, and FH⋯Cl^- and E2 channels, respectively. Anti-E2 dominates over syn-E2 (kinetic anti-E2 preference) and the thermodynamically-favored S_N2 (wider reactive anti-E2 attack angle range). The direct (a) S_N2, (b) proton abstraction, (c) FH⋯Cl^- + C₂H₄, (d) syn-E2, and (e) anti-E2 channels proceed with (a) back-side/backward, (b) isotropic/forward, (c) side-on/forward, (d) front-side/forward, and (e) back-side/forward attack/scattering, respectively. The HF products are vibrationally cold, especially for proton abstraction, and their rotational excitation increases for proton abstraction, anti-E2, and syn-E2, in order. Product internal-energy and mode-specific vibrational distributions show that CH₃CH₂F is internally hot with significant C-F stretching and CH₂ wagging excitations, whereas C₂H₄ is colder. One-dimensional Gaussian binning technique is proved to solve the normal mode analysis failure caused by methyl internal rotation.

Vibrational control of the reaction pathway in the H + CHD3 → H2 + CD3 reaction

An accurate full-dimensional quantum state-to-state simulation of the six-atom title reaction based on first-principles theory is reported. Counterintuitive effects are found: Increasing the energy in the reactant's CD₃ umbrella vibration reduces the energy in the corresponding product vibration. An in-depth analysis reveals the crucial role of the effective dynamical transition state: Its geometry is controlled by the vibrational states of the reactants and subsequently controls the quantum state distribution of the products. This finding enables generalizing the concept of transition state control of chemical reactions to the quantum state-specific level.

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.

The symmetric C-D stretching spectator mode in the H + CHD3 → H2 + CD3 reaction and its effect on dynamical modeling

The symmetric C-D stretching mode is a spectator mode in the H + CHD3 → H2 + CD3 reaction. Effects of multiple vibrational excitations of the CHD3 reactant are studied with the quantum transition-state (QTS) framework and an eight-dimensional (8D) model Hamiltonian developed by Palma and Clary. By including many thermal flux eigenstates, results have been obtained up to high energies, allowing the study of the symmetric C-D stretching spectator mode. A new concept of a state-specific thermal flux operator is proposed to analyze the C-D stretching spectator mode in detail, providing a new and insightful venue for studying transition-state control of chemical reactions. Furthermore, as a spectator mode, whether the C-D stretching motion can be excluded in a seven-dimensional (7D) model has not been fully interrogated, although the 7D model is a reasonable approximation and has provided accurate theoretical predictions. By comparing with available results of full-dimensional calculations, both the 7D and 8D models predict reasonably accurate results. However, the 7D model underestimates the mixing of two vibrational states that are in Fermi resonance. Despite its spectator nature, the C-D stretch is important in the dynamical modeling of chemical reaction systems affected by state mixing.

Full-dimensional quantum dynamics study of isotope effects for the H2 + NH2/ND2/NHD and H2/D2/HD + NH2 reactions

A full-dimensional quantum dynamical study for the bimolecular reactions of hydrogen molecules with amino radicals for different isotopologues is reported. The nonreactive amino radical is described by two Radau vectors that are very close to the valence bond coordinates. Potential-optimized discrete variable representation basis is used for the vibrational coordinates of the amino radical. Starting from the reaction H₂ + NH₂, we study the isotope effects for the two reagents separately, i.e., H₂ + NH₂/ND₂/NHD and H₂/D₂/HD + NH₂. The effects of different vibrational mode excitations of the reagents on the reactivities are studied. Physical explanations about the isotope effects are also provided thoroughly including the influence of vibrational energy differences between the different isotopologues and the impact of the tunneling effect.

Time-dependent quantum mechanical wave packet dynamics

Starting from a model study of the collinear (H, H₂) exchange reaction in 1959, the time-dependent quantum mechanical wave packet (TDQMWP) method has come a long way in dealing with systems as large as Cl + CH₄. The fast Fourier transform method for evaluating the second order spatial derivative of the wave function and split-operator method or Chebyshev polynomial expansion for determining the time evolution of the wave function for the system have made the approach highly accurate from a practical point of view. The TDQMWP methodology has been able to predict state-to-state differential and integral reaction cross sections accurately, in agreement with available experimental results for three dimensional (H, H₂) collisions, and identify reactive scattering resonances too. It has become a practical computational tool in predicting the observables for many A + BC exchange reactions in three dimensions and a number of larger systems. It is equally amenable to determining the bound and quasi-bound states for a variety of molecular systems. Just as it is able to deal with dissociative processes (without involving basis set expansion), it is able to deal with multi-mode nonadiabatic dynamics in multiple electronic states with equal ease. We present an overview of the method and its strength and limitations, citing examples largely from our own research groups.

Advances and New Challenges to Bimolecular Reaction Dynamics Theory

Dynamics of bimolecular reactions in the gas phase are of foundational importance in combustion, atmospheric chemistry, interstellar chemistry, and plasma chemistry. These collision-induced chemical transformations are a sensitive probe of the underlying potential energy surface(s). Despite tremendous progress in past decades, our understanding is still not complete. In this Perspective, we survey the recent advances in theoretical characterization of bimolecular reaction dynamics, stimulated by new experimental observations, and identify key new challenges.

Full-dimensional quantum mechanical calculations of the reaction probability of the H + CH4 reaction based on a mixed Jacobi and Radau description

A full-dimensional time-dependent wave packet study using mixed polyspherical Jacobi and Radau coordinates for the title reaction has been reported. The non-reactive moiety CH₃ has been described using three Radau vectors, whereas two Jacobi vectors have been used for the bond breaking/formation process. A potential-optimized discrete variable representation basis has been employed to describe the vibrational coordinates of the reagent CH₄. About one hundred billion basis functions have been necessary to achieve converged results. The reaction probabilities for some initial vibrational states are given. A comparison between the present approach and other methods, including reduced and full-dimensional ones, is also presented.

Imaging pair-correlated reaction cross sections in F + CH3D(νb = 0, 1) → CH2D(ν4 = 1) + HF(ν)

The title reactions were studied in a crossed-beam experiment at collisional energies (Ec) from 0.5 to 4.7 kcal mol-1. The νb (ν4) vibrational mode denotes the bending (umbrella) motion of the CH3D reactant (CH2D product). Using a time-sliced, velocity-map imaging technique, we extracted the state-specific, pair-correlated integral and differential cross sections. As with other isotopically analogous ground-state reactions, an inverted vibrational population of the HF coproduct was observed. Both the step-like excitation function near the threshold and the oscillatory forward-backward peakings in the Ec-evolution of the two dominant pair-correlated angular distributions at lower Ec suggest a resonance-mediated, time-delay mechanism. As Ec increases, the angular distribution of the HF(ν = 2) product evolves into a smooth and broad swath in the backward hemisphere, indicative of a direct rebound mechanism. One quantum excitation of the bending modes of CH3D(νb = 1) promotes the reaction rate by two- to three-fold up to Ec = 2.1 kcal mol-1. Broadly speaking, all major findings are qualitatively in line with previous results in the reactions of the F atom with other isotopologues. However, the rainbow feature recently observed in the CH2D(00) + HF(ν = 3) product channel is entirely absent. A possible rationale is put forward, which reinforces the previous reactive rainbow conjecture and calls for future theoretical investigations.

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.

Natural reaction channels in H + CHD3 → H2 + CD3

Natural reaction channels control the mode-specific chemistry of methane and its isotopomeres.

A neural network potential energy surface for the F + CH4reaction including multiple channels based on coupled cluster theory

We report here a new global and full dimensional potential energy surface (PES) for the F + CH₄reaction.

Dynamics and Novel Mechanisms of SN2 Reactions on ab Initio Analytical Potential Energy Surfaces

We describe a novel theoretical approach to the bimolecular nucleophilic substitution (S_N2) reactions that is based on analytical potential energy surfaces (PESs) obtained by fitting a few tens of thousands high-level ab initio energy points. These PESs allow computing millions of quasi-classical trajectories thereby providing unprecedented statistical accuracy for S_N2 reactions, as well as performing high-dimensional quantum dynamics computations. We developed full-dimensional ab initio PESs for the F^- + CH₃Y [Y = F, Cl, I] systems, which describe the direct and indirect, complex-forming Walden-inversion, the frontside attack, and the new double-inversion pathways as well as the proton-transfer channels. Reaction dynamics simulations on the new PESs revealed (a) a novel double-inversion S_N2 mechanism, (b) frontside complex formation,

Dynamics of transient speciesviaanion photodetachment

Recent experimental and theoretical advances in transient reaction dynamics probed by photodetachment of polyatomic anions are reviewed.

Recent advances in quantum scattering calculations on polyatomic bimolecular reactions

This review surveys quantum scattering calculations on chemical reactions of polyatomic molecules in the gas phase published in the last ten years.

CH stretching excitation promotes its cleavage in the F + CHD3(ν1 = 1) → HF + CD3 reaction at low collision energies

The effects of CH stretching excitation on the reactivity of the F + CHD₃ → HF + CD₃ reaction were experimentally studied.

Mode-Specific SN2 Reaction Dynamics

Despite its importance in chemistry, the microscopic dynamics of bimolecular nucleophilic substitution (SN2) reactions is still not completely elucidated. In this publication, the dynamics of a prototypical SN2 reaction (F(-) + CH3Cl → CH3F + Cl(-)) is investigated using a high-dimensional quantum mechanical model on an accurate potential energy surface (PES) and further analyzed by quasi-classical trajectories on the same PES. While the indirect mechanism dominates at low collision energies, the direct mechanism makes a significant contribution. The reactivity is found to depend on the specific reactant vibrational mode excitation. The mode specificity, which is more prevalent in the direct reaction, is rationalized by a transition-state-based model.