Numerical separation of the front-side attack and double-inversion retention pathways of SN2 reactions

Automated potential energy surface development and quasi-classical dynamics for the F- + SiH3I system

We report a potential energy surface (PES) development for the F- + SiH3I system to study its gas-phase reactions through quasi-classical dynamics simulations. The PES is represented by a full-dimensional permutationally invariant polynomial fitted to composite coupled cluster energy points obtained at the ManyHF-[CCSD-F12b + BCCD(T) - BCCD]/aug-cc-pVTZ(-PP) level of theory. The development was automated by Robosurfer, which samples the configurational space, manages ab initio calculations, and iteratively extends the fitting set. When selecting the ab initio method, we address two types of electronic structure calculation issues: first, the gold standard CCSD(T)-F12b is prone to occasional breakdown due to the perturbative (T) contribution, whereas CCSD-F12b + BCCD(T) - BCCD, with the Brueckner (T) term, is more robust; second, the underlying Hartree-Fock calculation may not always converge to the global minimum, resulting in highly erroneous energies. To mitigate this, we employed ManyHF, configuring the Hartree-Fock calculations with multiple initial guess orbitals and selecting the solution with the lowest energy. According to the simulations, the title system exhibits exceptionally high and diverse reactivity. We observe two dominant product formations: SN2 and proton abstraction. Moreover, SiH2F- + HI, SiHFI- + H2, SiH2FI + H-, SiH2 + FHI-, SiH2 + HF + I-, and SiHF + H2 + I- formations are found at lower probabilities. We differentiated inversion and retention for SN2, both being significant throughout the entire collision energy range. Opacity- and excitation functions are reported, and the details of the atomistic dynamics are visually examined via trajectory animations.

Backside versus Frontside SN2 Reactions of Alkyl Triflates and Alcohols

Nucleophilic substitution reactions are elementary reactions in organic chemistry that are used in many synthetic routes. By quantum chemical methods, we have investigated the intrinsic competition between the backside SN2 (SN2-b) and frontside SN2 (SN2-f) pathways using a set of simple alkyl triflates as the electrophile in combination with a systematic series of phenols and partially fluorinated ethanol nucleophiles. It is revealed how and why the well-established mechanistic preference for the SN2-b pathway slowly erodes and can even be overruled by the unusual SN2-f substitution mechanism going from strong to weak alcohol nucleophiles. Activation strain analyses disclose that the SN2-b pathway is favored for strong alcohol nucleophiles because of the well-known intrinsically more efficient approach to the electrophile resulting in a more stabilizing nucleophile-electrophile interaction. In contrast, the preference of weaker alcohol nucleophiles shifts to the SN2-f pathway, benefiting from a stabilizing hydrogen bond interaction between the incoming alcohol and the leaving group. This hydrogen bond interaction is strengthened by the increased acidity of the weaker alcohol nucleophiles, thereby steering the mechanistic preference toward the frontside SN2 pathway.

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.

Exploring the versatile reactivity of the F- + SiH3Cl system on a full-dimensional coupled-cluster potential energy surface

We develop a full-dimensional analytical potential energy surface (PES) for the F- + SiH3Cl reaction using Robosurfer for automatically sampling the configuration space, the robust [CCSD-F12b + BCCD(T) - BCCD]/aug-cc-pVTZ composite level of theory for computing the energy points, and the permutationally invariant polynomial method for fitting. Evolution of the fitting error and the percentage of the unphysical trajectories are monitored as a function of the iteration steps/number of energy points and polynomial order. Quasi-classical trajectory simulations on the new PES reveal rich dynamics resulting in high-probability SN2 (SiH3F + Cl-) and proton-transfer (SiH2Cl- + HF) products as well as several lower-probability channels, such as SiH2F- + HCl, SiH2FCl + H-, SiH2 + FHCl-, SiHFCl- + H2, SiHF + H2 + Cl-, and SiH2 + HF + Cl-. The Walden-inversion and front-side-attack-retention SN2 pathways are found to be competitive, producing nearly racemic products at high collision energies. The detailed atomic-level mechanisms of the various reaction pathways and channels as well as the accuracy of the analytical PES are analyzed along representative trajectories.

Spectroscopic Study of the Br- +CH3 I→I- +CH3 Br SN 2 Reaction

Mass spectrometry and anion photoelectron spectroscopy have been used to study the gas-phaseS N 2 ${{{\rm S}}_{{\rm N}}2}$ reaction involvingB r - ${{{\rm B}{\rm r}}^{-}}$ andC H 3 I ${{{\rm C}{\rm H}}_{3}{\rm I}}$ . The anion photoelectron spectra associated with the reaction intermediates of thisS N 2 ${{{\rm S}}_{{\rm N}}2}$ reaction are presented. High-level CCSD(T) calculations have been utilised to investigate the reaction intermediates that may form as a result of theS N 2 ${{{\rm S}}_{{\rm N}}2}$ reaction along various different reaction pathways, including back-side attack and front-side attack. In addition, simulated vertical detachment energies of each reaction intermediate have been calculated to rationalise the photoelectron spectra.

Dynamical Effects of SN2 Reactivity Suppression by Microsolvation: Dynamics Simulations of the F-(H2O) + CH3I Reaction on a 21-Dimensional Potential Energy Surface

A comparison of atomistic dynamics between microsolvated and unsolvated reactions can expose the precise role of solvent molecules and thus provide deep insight into how solvation influences chemical reactions. Here we developed the first full-dimensional analytical potential energy surface of the F-(H2O) + CH3I reaction, which facilitates the efficient dynamics simulations on a quantitatively accurate level. The computed SN2 reactivity suppression ratio of the monosolvated F-(H2O) + CH3I reaction relative to the unsolvated F- + CH3I reaction as a function of collision energy first increases and then decreases steadily, forming an inverted-V shape, due to the combined dynamical effects of interaction time, steric hindrance, and collision-induced dehydration. Moreover, further analysis reveals that the steric effect of the F-(H2O) + CH3I reaction resulting from the single water molecule is manifested mainly in dragging the F- anion away from the central C atom, rather than shielding F- from C. Our study shows there is great potential in rigorously studying the role of the solvent in more complicated reactions.

Unconventional SN2 retention pathways induced by complex formation: High-level dynamics investigation of the NH2 - + CH3I polyatomic reaction

Investigations on the dynamics of chemical reactions have been a hot topic for experimental and theoretical studies over the last few decades. Here, we carry out the first high-level dynamical characterization for the polyatom-polyatom reaction between NH₂ ^- and CH₃I. A global analytical potential energy surface is developed to describe the possible pathways with the quasi-classical trajectory method at several collision energies. In addition to S_N2 and proton abstraction, a significant iodine abstraction is identified, leading to the CH₃ + [NH₂⋯I]^- products. For S_N2, our computations reveal an indirect character as well, promoting the formation of [CH₃⋯NH₂] complexes. Two novel dominant S_N2 retention pathways are uncovered induced by the rotation of the CH₃ fragment in these latter [CH₃⋯NH₂] complexes. Moreover, these uncommon routes turn out to be the most dominant retention paths for the NH₂ ^- + CH₃I S_N2 reaction.

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.

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₃).

Uncovering an oxide ion substitution for the OH- + CH3F reaction

Theoretical investigations on chemical reactions allow us to understand the dynamics of the possible pathways and identify new unexpected routes. Here, we develop a global analytical potential energy surface (PES) for the OH⁻ + CH₃F reaction in order to perform high-level dynamics simulations. Besides bimolecular nucleophilic substitution (S_N2) and proton abstraction, our quasi-classical trajectory computations reveal a novel oxide ion substitution leading to the HF + CH₃O⁻ products. This exothermic reaction pathway occurs via the CH₃OH⋯F⁻ deep potential well of the S_N2 product channel as a result of a proton abstraction from the hydroxyl group by the fluoride ion. The present detailed dynamics study of the OH⁻ + CH₃F reaction focusing on the surprising oxide ion substitution demonstrates how incomplete our knowledge is of fundamental chemical reactions.

Reaction dynamics simulations on a high-level ab initio analytical potential energy surface reveal a novel oxide ion substitution channel for the OH⁻ + CH₃F reaction.

Detailed quasiclassical dynamics of the F- + CH3Br reaction on an ab initio analytical potential energy surface

Dynamics and mechanisms of the F^- + CH₃Br(v = 0) → Br^- + CH₃F (S_N2 via Walden inversion, front-side attack, and double inversion), F^- + inverted-CH₃Br (induced inversion), HF + CH₂Br^- (proton abstraction), and FH⋯Br^- + ¹CH₂ reactions are investigated using a high-level global ab initio potential energy surface, the quasiclassical trajectory method, as well as non-standard configuration- and mode-specific analysis techniques. A vector-projection method is used to identify inversion and retention trajectories; then, a transition-state-attack-angle-based approach unambiguously separates the front-side attack and the double-inversion retention pathways. The Walden-inversion S_N2 channel becomes direct rebound dominated with increasing collision energy as indicated by backward scattering, initial back-side attack preference, and the redshifting of product internal energy peaks in accord with CF stretching populations. In the minor retention and induced-inversion pathways, almost the entire available energy transfers into product rotation-vibration, and retention mainly proceeds with indirect, slow double inversion following induced inversion with about 50% probability. Proton abstraction is dominated by direct stripping (evidenced by forward scattering) with CH₃-side initial attack preference, providing mainly vibrationally ground state products with significant zero-point energy violation.

Facilitated inversion complicates the stereodynamics of an SN2 reaction at nitrogen center

Bimolecular nucleophilic substitution (S_N2) reactions at carbon center are well known to proceed with the stereospecific Walden-inversion mechanism. Reaction dynamics simulations on a newly developed high-level ab initio analytical potential energy surface for the F⁻ + NH₂Cl nitrogen-centered S_N2 and proton-transfer reactions reveal a hydrogen-bond-formation-induced multiple-inversion mechanism undermining the stereospecificity of the N-centered S_N2 channel. Unlike the analogous F⁻ + CH₃Cl S_N2 reaction, F⁻ + NH₂Cl → Cl⁻ + NH₂F is indirect, producing a significant amount of NH₂F with retention, as well as inverted NH₂Cl during the timescale within the unperturbed NH₂Cl molecule gets inverted with only low probability, showing the important role of facilitated inversions via an FH…NHCl⁻-like transition state. Proton transfer leading to HF + NHCl⁻ is more direct and becomes the dominant product channel at higher collision energies.

Multiple-inversion, the analogue of the double-inversion pathway recently revealed for S_N2@C, is the key mechanism in S_N2 at N center undermining stereospecificity.