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

Vibrational mode-specificity in the dynamics of the OH- + CH3I multi-channel reaction

We report a comprehensive characterization of the vibrational mode-specific dynamics of the OH- + CH3I reaction. Quasi-classical trajectory simulations are performed at four different collision energies on our previously-developed full-dimensional high-level ab initio potential energy surface in order to examine the impact of four different normal-mode excitations in the reactants. Considering the 11 possible pathways of OH- + CH3I, pronounced mode-specificity is observed in reactivity: In general, the excitations of the OH- stretching and CH stretching exert the greatest influence on the channels. For the SN2 and proton-abstraction products, the reactant initial attack angle and the product scattering angle distributions do not show major mode-specific features, except for SN2 at higher collision energies, where forward scattering is promoted by the CI stretching and CH stretching excitations. The post-reaction energy flow is also examined for SN2 and proton abstraction, and it is unveiled that the excess vibrational excitation energies rather transfer into the product vibrational energy because the translational and rotational energy distributions of the products do not represent significant mode-specificity. Moreover, in the course of proton abstraction, the surplus vibrational energy in the OH- reactant mostly remains in the H2O product owing to the prevailing dominance of the direct stripping mechanism.

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.

Predicting atomic-level reaction mechanisms for SN2 reactions via machine learning

Identifying atomic-level reaction mechanisms is an essential step in chemistry. In this study, we develop a joint-voting model based on three parallel machine-learning algorithms to predict atomic-level and dynamical mechanisms trained with 1700 trajectories. Three predictive experiments are carried out with the training trajectories divided into ten, seven, and five classes. The results indicate that, as the number of trajectories in each class increases from the ten- to five-class model, the five-class model converges the fastest and the prediction success rate increases. The number of trajectories in each experiment to get the predictive models converged is 100, 100, and 70, respectively. The prediction accuracy increases from 88.3% for the ten-class experiment, to 91.0% for the seven-class, and to 92.0% for the five-class. Our study demonstrates that machine learning can also be used to predict elementary dynamical processes of structural evolution along time, that is, atomic-level reaction mechanisms.

Fifty years of nucleophilic substitution in the gas phase

Bimolecular nucleophilic substitution ( S N 2 ) reactions have become a model system for the investigation of structure-reactivity relationships, stereochemistry, solvent influences, and detailed atomistic dynamics. In this review, the progress during five decades of experimental and theoretical research on gas phase S N 2 reactions is discussed. Many advancements of the employed methods have led to a tremendous increase in our understanding of the properties and the dynamics of these reactions. For reactions involving six atoms a quantitative agreement of the differential reactive scattering cross sections has already been achieved, in the future it is expected that even larger polyatomic reactions systems become tractable. Furthermore, studies with higher precision, improved reactant control, and a more accurate theoretical treatment of quantum effects are envisioned.

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.

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.

The importance of the composite mechanisms with two transition states in the F- + NH2I SN2 reaction

The dynamics of the bimolecular nucleophilic substitution (S_N2) reactions at nitrogen are less understood than those of their corresponding reactions at carbon. In this paper, we report an ab initio molecular dynamics approach to investigate the reaction mechanisms of the F^- + NH₂I S_N2 reaction at nitrogen. We found not only the one-transition-state mechanisms, but also the composite mechanisms with two and three transition states. For the two-transition-state mechanisms, the double inversion mechanism and the proton-abstraction roundabout followed by the backside-attack reaction mechanism have been reported before; but we discovered that there is also a new, front-side attack followed by the backside-attack Walden-inversion mechanism. Furthermore, a composite mechanism with three transition states also shows up in the reactive trajectories. Our results show that, as the collision energy increases, the S_N2 reactivity decreases, and the proton-abstraction reactivity increases. The two-transition-state mechanisms, especially the double-inversion mechanism, make the largest contribution to the S_N2 reactivity, followed then by the one-transition-state mechanisms, with the three-transition-state mechanism contributing the least. The potential energy profiles of the reaction mechanisms are characterized at the CCSD(T)/aug-cc-pVTZ(PP) level of theory. The analysis on stationary points shows that the proton-abstraction inversion transition state is ∼12.4 kcal mol^-1 lower than the Walden-inversion transition state in contrast to the corresponding reaction at carbon F^- + CH₃I, in which the former is ∼26.1 kcal mol^-1 higher than the latter. This might explain why the composite mechanism of the double inversion mechanism contributes the most to the S_N2 reactivity in the F^- + NH₂I reaction.

Proton transfer dynamics modified by CH-stretching excitation

Gaining insight how specific rovibrational states influence reaction kinetics and dynamics is a fundamental goal of physical chemistry. Purely statistical approaches often fail to predict the influence of a specific state on the reaction outcome, evident in a great number of both experimental and theoretical studies. Most detailed insight in atomistic reaction mechanisms is achieved using accurate collision experiments and high level dynamics calculations. For ion-molecule reactions such experiments are scarce. Here we show the influence of symmetric CH-stretching vibration on the rate and dynamics of proton transfer in the reaction of F- + CH3I. We find a pronounced shift in the reaction dynamics for excited reactions from indirect to preferred direct dynamics at higher collision energy. Moreover, excited reactions occur at larger impact parameters. Finally, we compare vibrational excitation with collision energy and find that vibration is overall more efficient in promoting reactivity, which agrees with recent theoretical calculations.

Influence of Vibrational Excitation on the Reaction of F- with CH3I: Spectator Mode Behavior, Enhancement, and Suppression

Detailed insight into chemical reaction dynamics can be obtained by probing the effect of mode-specific vibrational excitation. Suppression or enhancement of reactivity is possible as is already known from the Polanyi rules. In the reaction F- + CH3I, we found vibrational enhancement, suppression, and spectator mode dynamics in the four different reaction channels. For this system we have probed the influence of symmetric CH-stretching vibration over a collision energy range of 0.7-2.3 eV. Proton transfer is significantly enhanced, while for the nucleophilic substitution channel the spectator mode dynamics at lower collision energies unexpectedly move toward enhancement at higher collision energies. In contrast, for two halide abstraction channels, forming FI- and FHI-, we found an overall suppression, which stems mainly from a suppression of the FHI- product. We compare these results to quasiclassical trajectory calculations and with the sudden vector projection model.