Benchmark ab initio characterization of the multi-channel Cl + CH3X [X = F, Cl, Br, I] reactive potential energy surfaces

We determine benchmark geometries and relative energies for the stationary points of the Cl + CH3X [X = F, Cl, Br, I] reactions. We consider four possible reaction pathways: hydrogen abstraction, hydrogen substitution, halogen abstraction, and halogen substitution, where the substitution processes can proceed via either Walden inversion or front-side attack. We perform geometry optimizations and obtain harmonic vibrational frequencies at the explicitly-correlated UCCSD(T)-F12b/aug-cc-pVTZ level of theory, followed by UCCSD(T)-F12b/aug-cc-pVQZ single-point computations to make finite-basis-set error negligible. To reach chemical (<1 kcal mol-1), or even subchemical (<0.5 kcal mol-1) accuracy, we include core-correlation, scalar relativistic, post-(T), spin-orbit-splitting and zero-point-energy contributions, as well, in the relative energies of all the stationary points. Our benchmark 0 K reaction enthalpies are compared to available experimental results and show good agreement. The stationary-point structures and energetics are interpreted in terms of Hammond's postulate and used to make predictions related to the dynamical behavior of these reactive systems.

Benchmark ab initio characterization of the complex potential energy surfaces of the HOO- + CH3Y [Y = F, Cl, Br, I] reactions

The α-effect is a well-known phenomenon in organic chemistry, and is related to the enhanced reactivity of nucleophiles involving one or more lone-pair electrons adjacent to the nucleophilic center. The gas-phase bimolecular nucleophilic substitution (SN2) reactions of α-nucleophile HOO- with methyl halides have been thoroughly investigated experimentally and theoretically; however, these investigations have mainly focused on identifying and characterizing the α-effect of HOO-. Here, we perform the first comprehensive high-level ab initio mapping for the HOO- + CH3Y [Y = F, Cl, Br and I] reactions utilizing the modern explicitly-correlated CCSD(T)-F12b method with the aug-cc-pVnZ [n = 2-4] basis sets. The present ab initio characterization considers five distinct product channels of SN2: (CH3OOH + Y-), proton abstraction (CH2Y- + H2O2), peroxide ion substitution (CH3OO- + HY), SN2-induced elimination (CH2O + HY + HO-) and SN2-induced rearrangement (CH2(OH)O- + HY). Moreover, besides the traditional back-side attack Walden inversion, the pathways of front-side attack, double inversion and halogen-bond complex formation have also been explored for SN2. With regard to the Walden inversion of HOO- + CH3Cl, the previously unaddressed discrepancies concerning the geometry of the corresponding transition state are clarified. For the HOO- + CH3F reaction, the recently identified SN2-induced elimination is found to be more exothermic than the SN2 channel, submerged by ∼36 kcal mol-1. The accuracy of our high-level ab initio calculations performed in the present study is validated by the fact that our new benchmark 0 K reaction enthalpies show excellent agreement with the experimental data in nearly all cases.

Phosphorus-centered ion-molecule reactions: benchmark ab initio characterization of the potential energy surfaces of the X- + PH2Y [X, Y = F, Cl, Br, I] systems

In the present work we determine the benchmark relative energies and geometries of all the relevant stationary points of the X- + PH2Y [X, Y = F, Cl, Br, I] identity and non-identity reactions using state-of-the-art electronic-structure methods. These phosphorus-centered ion-molecule reactions follow two main reaction routes: bimolecular nucleophilic substitution (SN2), leading to Y- + PH2X, and proton transfer, resulting in HX + PHY- products. The SN2 route can proceed through Walden-inversion, front-side-attack retention, and double-/multiple-inversion pathways. In addition, we also identify the following product channels: H--formation, PH2-- and PH2-formation, 1PH- and 3PH-formation, H2-formation and HY + PHX- formation. The benchmark classical relative energies are obtained by taking into account the core-correlation, scalar relativistic, and post-(T) corrections, which turn out to be necessary to reach subchemical (<1 kcal mol-1) accuracy of the results. Classical relative energies are augmented with zero-point-energy contributions to gain the benchmark adiabatic energies.

Reaction mechanism conversion induced by the contest of nucleophile and leaving group

Direct dynamic simulations have been employed to investigate the OH^- + CH₃Cl reaction with the chosen B3LYP/aug-cc-pVDZ method. The calculated rate coefficient for the bimolecular nucleophilic substitution reaction (S_N2), 1.0 × 10^-9 cm³ mol^-1 s^-1 at 300 K, agrees well with the experimental result of (1.3-1.6) × 10^-9 cm³ mol^-1 s^-1. The simulations reveal that the majority of the S_N2 reactions are temporarily trapped in the hydrogen-bonded complex at E_coll = 0.89 kcal mol^-1. Importantly, the influences of the leaving group and nucleophile have been discussed by comparisons of X^- + CH₃Y (X = F, OH; Y = Cl, I) reactions. For the X = F^- reactions, the reaction probability of S_N2 increases along the increased leaving group ability Cl < I, suggesting that the thermodynamic factor plays a key role. The indirect mechanisms were found to be dominant for both reactions. In contrast, for X = OH^-, the fraction of S_N2 drops with the enhanced leaving group ability. In particular, a dramatic transition occurs for the dominant atomic reaction mechanisms, i.e., from complex-mediated indirect to direct, implying an interesting contest between the leaving group and the nucleophile and the importance of the dynamic factors, i.e., the dipole moment, steric hindrance, and electronegativity.

Current Status of the X + C2H6 [X ≡ H, F(2P), Cl(2P), O(3P), OH] Hydrogen Abstraction Reactions: A Theoretical Review

This paper is a detailed review of the chemistry of medium-size reactive systems using the following hydrogen abstraction reactions with ethane, X + C₂H₆ → HX + C₂H₅; X ≡ H, F(²P), Cl(²P), O(³P) and OH, and focusing attention mainly on the theoretical developments. These bimolecular reactions range from exothermic to endothermic systems and from barrierless to high classical barriers of activation. Thus, the topography of the reactive systems changes from reaction to reaction with the presence or not of stabilized intermediate complexes in the entrance and exit channels. The review begins with some reflections on the inherent problems in the theory/experiment comparison. When one compares kinetics or dynamics theoretical results with experimental measures, one is testing both the potential energy surface describing the nuclei motion and the kinetics or dynamics method used. Discrepancies in the comparison may be due to inaccuracies of the surface, limitations of the kinetics or dynamics methods, and experimental uncertainties that also cannot be ruled out. The paper continues with a detailed review of some bimolecular reactions with ethane, beginning with the reactions with hydrogen atoms. The reactions with halogens present a challenge owing to the presence of stabilized intermediate complexes in the entrance and exit channels and the influence of the spin-orbit states on reactivity. Reactions with O(³P) atoms lead to three surfaces, which is an additional difficulty in the theoretical study. Finally, the reactions with the hydroxyl radical correspond to a reactive system with ten atoms and twenty-four degrees of freedom. Throughout this review, different strategies in the development of analytical potential energy surfaces describing these bimolecular reactions have been critically analyzed, showing their advantages and limitations. These surfaces are fitted to a large number of ab initio calculations, and we found that a huge number of calculations leads to accurate surfaces, but this information does not guarantee that the kinetics and dynamics results match the experimental measurements.

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.

High-Level Systematic Ab Initio Comparison of Carbon- and Silicon-Centered SN2 Reactions

We characterize the stationary points along the Walden inversion, front-side attack, and double-inversion pathways of the X^– + CH₃Y and X^– + SiH₃Y [X, Y = F, Cl, Br, I] S_N2 reactions using chemically accurate CCSD(T)-F12b/aug-cc-pVnZ [n = D, T, Q] levels of theory. At the carbon center, Walden inversion dominates and proceeds via prereaction (X^–···H₃CY) and postreaction (XCH₃···Y^–) ion-dipole wells separated by a usually submerged transition state (X–H₃C–Y)⁻, front-side attack occurs over high barriers, double inversion is the lowest-energy retention pathway for X = F, and hydrogen- (F^–···HCH₂Y) and halogen-bonded (X^–···YCH₃) complexes exist in the entrance channel. At the silicon center, Walden inversion proceeds through a single minimum (X–SiH₃–Y)⁻, the front-side attack is competitive via a usually submerged transition state separating pre- and postreaction minima having X–Si–Y angles close to 90°, double inversion occurs over positive, often high barriers, and hydrogen- and halogen-bonded complexes are not found. In addition to the S_N2 channels (Y^– + CH₃X/SiH₃X), we report reaction enthalpies for proton abstraction (HX + CH₂Y^–/SiH₂Y^–), hydride substitution (H^– + CH₂XY/SiH₂XY), XH···Y^– complex formation (XH···Y^– + ¹CH₂/¹SiH₂), and halogen abstraction (XY + CH₃^–/SiH₃^– and XY^– + CH₃/SiH₃).

Visualization of reaction route map and dynamical trajectory in reduced dimension

In the quantum chemical approach, chemical reaction mechanisms are investigated based on a potential energy surface (PES). Automated reaction path search methods enable us to construct a global reaction route map containing multiple reaction paths corresponding to a series of elementary reaction processes. The on-the-fly molecular dynamics (MD) method provides a classical trajectory exploring the full-dimensional PES based on electronic structure calculations. We have developed two reaction analysis methods, the on-the-fly trajectory mapping method and the reaction space projector (ReSPer) method, by introducing a structural similarity to a pair of geometric structures and revealed dynamic aspects affecting chemical reaction mechanisms. In this review, we will present the details of these analysis methods and discuss the dynamics effects of reaction path curvature and reaction path bifurcation with applications to the CH₃OH + OH^- collision reaction and the Au₅ cluster branching and isomerization reactions.

Atomistic dynamics of elimination and nucleophilic substitution disentangled for the F- + CH3CH2Cl reaction

Chemical reaction dynamics are studied to monitor and understand the concerted motion of several atoms while they rearrange from reactants to products. When the number of atoms involved increases, the number of pathways, transition states and product channels also increases and rapidly presents a challenge to experiment and theory. Here we disentangle the dynamics of the competition between bimolecular nucleophilic substitution (S_N2) and base-induced elimination (E2) in the polyatomic reaction F^- + CH₃CH₂Cl. We find quantitative agreement for the energy- and angle-differential reactive scattering cross-sections between ion-imaging experiments and quasi-classical trajectory simulations on a 21-dimensional potential energy hypersurface. The anti-E2 pathway is most important, but the S_N2 pathway becomes more relevant as the collision energy is increased. In both cases the reaction is dominated by direct dynamics. Our study presents atomic-level dynamics of a major benchmark reaction in physical organic chemistry, thereby pushing the number of atoms for detailed reaction dynamics studies to a size that allows applications in many areas of complex chemical networks and environments.

A benchmark ab initio study of the complex potential energy surfaces of the OH- + CH3CH2Y [Y = F, Cl, Br, I] reactions

We provide the first benchmark characterization of the OH^- + CH₃CH₂Y [Y = F, Cl, Br, I] reactions utilizing the high-level explicitly-correlated CCSD(T)-F12b method with the aug-cc-pVnZ [n = 2(D), 3(T), 4(Q)] basis sets. We explore and analyze the stationary points of the elimination (E2) and substitution (S_N2) reactions, including anti-E2, syn-E2, back-side attack, front-side attack, and double inversion. In all cases, S_N2 is thermodynamically more preferred than E2. In the entrance channel of S_N2 a significant front-side complex formation is revealed, and in the product channel the global minimum of the title reactions is obtained at the hydrogen-bonded CH₃CH₂OHY^- complex. Similar to the OH^- + CH₃Y reactions, double inversion can proceed via a notably lower-energy pathway than front-side attack, moreover, for Y = I double inversion becomes barrier-less. For the transition state of the anti-E2, a prominent ZPE effect emerges, giving an opportunity for a kinetically more favored pathway than back-side attack. In addition to S_N2 and E2, other possible product channels are considered, and in most cases, the benchmark reaction enthalpies are in excellent agreement with the experimental data.

First-Principles Reaction Dynamics beyond Six-Atom Systems

Moving beyond the six-atomic benchmark systems, we discuss the new age and future of first-principles reaction dynamics, which investigates complex, multichannel chemical reactions. We describe the methodology starting from the benchmark ab initio characterization of the stationary points, followed by full-dimensional potential energy surface (PES) developments and reaction dynamics computations. We highlight our composite ab initio approach providing benchmark stationary-point properties with subchemical accuracy, the Robosurfer program system enabling automatic PES development, and applications for the Cl + C2H6, F + C2H6, and OH- + CH3I post-six-atom reactions focusing on ab initio issues and their solutions as well as showing the excellent agreement between theory and experiment.

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.

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.