Steric Effects of Solvent Molecules on SN2 Substitution Dynamics

Benchmark Ab Initio Mapping of the F- + CH2ClI SN2 and Proton-Abstraction Reactions

The experimental and theoretical studies of gas-phase SN2 reactions have significantly broadened our understanding of the mechanisms governing even the simplest chemical processes. These investigations have not only advanced our knowledge of reaction pathways but also provided critical insights into the fundamental dynamics of chemical systems. Nevertheless, in the case of the prototypical X- + CH3Y → Y- + CH3X [X, Y = F, Cl, Br, and I] SN2 reactions, the effect of the additional halogenation of CH3Y has not been thoroughly explored. Thus, here, we perform the first high-level ab initio characterization of the F- + CH2ClI SN2 and proton-abstraction reactions utilizing the explicitly-correlated CCSD(T)-F12b method. Two possible SN2 channels leading to the Cl- + CH2FI and I- + CH2FCl products are distinguished, in which we investigate four different pathways of back-side attack Walden inversion, front-side attack, double inversion, and halogen-bonded complex formation. In order to obtain the benchmark energies of the geometries of the stationary points, determined at the CCSD(T)-F12b/aug-cc-pVTZ level of theory, additional computations are carried out considering the basis set effects, post-CCSD(T) correlations, and core corrections. Using the benchmark data, we assess the accuracy of the MP2, DF-MP2, MP2-F12, and DF-MP2-F12 methods as well. By comparing the present F- + CH2ClI system with the corresponding F- + CH3Y [Y = Cl and I] reactions, this study demonstrates that further halogenation of CH3Y significantly promotes the corresponding proton-abstraction and SN2 retention channels as well as the halogen-bonded complex formation, and as a consequence, the traditional back-side attack Walden-inversion mechanism becomes less pronounced.

Roundabout Mechanism of Ion-Molecule Nucleophilic Substitution Reactions

Roundabout (RA) is an important indirect mechanism for gas-phase X- + CH3Y → XCH3 + Y- SN2 reactions at a high collision energy. It refers to the rotation of the CH3-group by half or multiple circles upon the collision of incoming nucleophiles before substitution takes place. The RA mechanism was first discovered in the Cl- + CH3I SN2 reaction to explain the energy transfer observed in crossed molecular beam imaging experiments in 2008. Since then, the RA mechanism and its variants have been observed not only in multiple C-centered SN2 reactions, but also in N-centered SN2 reactions, proton transfer reactions, and elimination reactions. This work reviewed recent studies on the RA mechanism and summarized the characteristics of RA mechanisms in terms of variant types, product energy partitioning, and product velocity scattering angle distribution. RA mechanisms usually happen at small impact parameters and tend to couple with other mechanisms at relatively low collision energy, and the available energy of roundabout trajectories is primarily partitioned to internal energy. Factors that affect the importance of the RA mechanism were analyzed, including the type of leaving group and nucleophile, collision energy, and microsolvation. A massive leaving group and relatively high collision energy are prerequisite for the occurrence of the roundabout mechanism. Interestingly, when reacting with CH3I, the importance of RA mechanisms follows an order of Cl- > HO- > F-, and such a nucleophile dependence was attributed to the difference in proton affinity and size of the nucleophile.

Solvent-induced dual nucleophiles and the α-effect in the SN2 versus E2 competition

We have quantum chemically investigated how microsolvation affects the various E2 and SN2 pathways, their mutual competition, and the α-effect of the model reaction system HOO-(H2O)n + CH3CH2Cl, at the CCSD(T) level. Interestingly, we identify the dual nature of the α-nucleophile HOO- which, upon solvation, is in equilibrium with HO-. This solvent-induced dual appearance gives rise to a rich network of competing reaction channels. Among both nucleophiles, SN2 is always favored over E2, and this preference increases upon increasing microsolvation. Furthermore, we found a pronounced α-effect, not only for SN2 substitution but also for E2 elimination, i.e., HOO- is more reactive than HO- in both cases. Our activation strain and quantitative molecular orbital analyses reveal the physical mechanisms behind the various computed trends. In particular, we demonstrate that two recently proposed criteria, required for solvent-free nucleophiles to display the α-effect, must also be satisfied by microsolvated HOO-(H2O)n nucleophiles.

Shapeshifting Nucleophiles HO-(NH3)n React with Methyl Chloride

The microsolvated anions HO-(NH3)n were found to induce new nucleophile NH2-(H2O)(NH3)n-1 via intramolecular proton transfer. Hence, the ion-molecule nucleophilic substitution (SN2) reaction between CH3Cl and these shapeshifting nucleophiles lead to both the HO- path and NH2- path, meaning that the respective attacking nucleophile is HO- or NH2-. The CCSD(T) level of calculation was performed to characterize the potential energy surfaces. Calculations indicate that the HO- species are lower in energy than the NH2- species, and the SN2 reaction barriers are lower for the HO- path than the NH2--path. Incremental solvation increases the barrier for both paths. Comparison between HO-(NH3)n and HOO-(NH3)n confirmed the existence of an α-effect under microsolvated conditions. Comparison between HO-(NH3)n and HO-(H2O)n indicated that the more polarized H2O stabilizes the nucleophiles more than NH3, and thus, the hydrated systems have higher SN2 reaction barriers. The aforementioned barrier changes can be explained by the differential stabilization of the nucleophile and HOMO levels upon solvation, thus affecting the HOMO-LUMO interaction between the nucleophile and substrate. For the same kind of nucleophilic attacking atom, O or N, the reaction barrier has a good linear correlation with the HOMO level of the nucleophiles. Hence, the HOMO level or the binding energy of microsolvated nucleophiles is a good indicator to evaluate the order of barrier heights. This work expands our understanding of the microsolvation effect on prototype SN2 reactions beyond the water solvent.

Mechanistic Insights into the Proton Transfer and Substitution Dynamics of N-Atom Center Reactions: A Study of CH3O- with NH2Cl

Bimolecular substitution reactions involving N as the central atom have continuously improved our understanding of substitution dynamics. This work used chemical dynamics simulations to investigate the dynamics of NH2Cl with N as the central atom and the multiatomic nucleophile CH3O- and compared these results with the F- + NH2Cl reaction. The most noteworthy difference is in the competition between proton transfer (PT) and the SN2 pathways. Our results demonstrate that, for the CH3O- + NH2Cl system, the PT pathway is considerably more favorable than the SN2 pathway. In contrast, no PT pathway was observed for the F- + NH2Cl system at room temperature. This can be attributed to the exothermic reaction of the PT pathway for the CH3O- + NH2Cl reaction and is coupled with a more stable transition state compared to the substitution pathway. Furthermore, the bulky nature of the CH3O- group impedes its participation in SN2 reactions, which enhances both the thermodynamic and the dynamic advantages of the PT reaction. Interestingly, the atomic mechanism reveals that the PT pathway is primarily governed by indirect mechanisms, similar to the SN2 pathway, with trajectories commonly trapped in the entrance channel being a prominent feature. These trajectories are often accompanied by prolonged and frequent proton exchange or proton abstraction processes. This current work provides insights into the dynamics of N-centered PT reactions, which are useful in gaining a comprehensive understanding of the dynamics behavior of similar reactions.

Shapeshifting Nucleophile Singly Hydrated Hydroperoxide Anion Leads to the Occurrence of the Thermodynamically Unfavored SN2 Product

Single water molecules alone may introduce unusual features into the kinetics and dynamics of chemical reactions. The singly hydrated hydroperoxide anion, HOO-(H2O), was found to be a shapeshifting nucleophile, which can be transformed to HO- solvated by hydrogen peroxide HO-(HOOH). Herein, we performed direct dynamics simulations of its reaction with methyl iodide to investigate the effect of individual water molecules. In addition to the normal SN2 product CH3OOH, the thermodynamically unfavored proton transfer-induced HO--SN2 path (produces CH3OH) was also observed, contributing ∼4%. The simulated branching ratio of the HO--SN2 path exceeded the statistical estimation (0.6%) based on the free energy barrier difference. The occurrence of the HO--SN2 path was attributed to the shallow entrance channel well before a submerged saddle point, thus providing a region for extensive proton exchange and ultimately leading to the formation of CH3OH. In comparison, changing the leaving group from Cl to I increased the overall reaction rate as well as the proportion of the HO--SN2 path because the CH3I system has a smaller internal barrier. This work elucidates the importance of the dynamic effect introduced by a single solvent molecule to alter the product channel and kinetics of typical ion-molecule SN2 reactions.

Harnessing the High Interfacial Electric Fields on Water Microdroplets to Accelerate Menshutkin Reactions

Even though it is still an emerging field, the application of a high external electric field (EEF) as a green and efficient catalyst in synthetic chemistry has recently received significant attention for the ability to deliver remarkable control of reaction selectivity and acceleration of reaction rates. Here, we extend the application of the EEF to Menshutkin reactions by taking advantage of the spontaneous high electric field at the air-water interfaces of sprayed water microdroplets. Experimentally, a series of Menshutkin reactions were accelerated by 7 orders of magnitude. Theoretically, both density functional theory calculations and ab initio molecular dynamics simulations predict that the reaction barrier decreases significantly in the presence of oriented external electric fields, thereby supporting the notion that the electric fields in the water droplets are responsible for the catalysis. In addition, the ordered solvent and reactant molecules oriented by the electric field alleviate the steric effect of solvents and increase the successful collision rates, thus facilitating faster nucleophilic attack. The success of Menshutkin reactions in this study showcases the great potential of microdroplet chemistry for green synthesis.

Accurate fundamental invariant-neural network representation of ab initio potential energy surfaces

Highly accurate potential energy surfaces are critically important for chemical reaction dynamics. The large number of degrees of freedom and the intricate symmetry adaption pose a big challenge to accurately representing potential energy surfaces (PESs) for polyatomic reactions. Recently, our group has made substantial progress in this direction by developing the fundamental invariant-neural network (FI-NN) approach. Here, we review these advances, demonstrating that the FI-NN approach can represent highly accurate, global, full-dimensional PESs for reactive systems with even more than 10 atoms. These multi-channel reactions typically involve many intermediates, transition states, and products. The complexity and ruggedness of this potential energy landscape present even greater challenges for full-dimensional PES representation. These PESs exhibit a high level of complexity, molecular size, and accuracy of fit. Dynamics simulations based on these PESs have unveiled intriguing and novel reaction mechanisms, providing deep insights into the intricate dynamics involved in combustion, atmospheric, and organic chemistry.

Competitive dynamics of E2 and SN2 reaction driven by collision energy and leaving group

The competition between E2 and SN2 reactions is essential in organic chemistry. In this paper, the reaction mechanism of F- + CH3CH2Cl is investigated utilizing direct dynamics simulations, and unravel how the collision energy (Ecoll) and the leaving group affect the competition between SN2 and E2 in the F- + CH3CH2Y (Y = Cl and Br) reactions. Simulation results for F- + CH3CH2Cl reaction show that the anti-E2 channel is dominant, but with the increase of Ecoll from 0.04 to 1.9 eV the branching ratio of the anti-E2 pathway significantly decreases by 21%, and the SN2 pathway becomes more important. A transition from indirect to direct reaction has been revealed when Ecoll is increased from 0.04 to 1.90 eV. At lower Ecoll, a large ratio of indirect events occurs via a long-lived hydrogen-bonded complex, and as the collision energy is increased, the lifetimes of the hydrogen-bonded complexes are shortened, due to an initial faster relative velocity. The simulation results of F- + CH3CH2Cl are further compared with the F- + CH3CH2Br reaction at Ecoll of 0.04 eV. Changing the leaving group from Cl to Br drastically suppresses the indirect events of anti-E2 with a branching ratio decreasing from 0.46 to 0.36 due to the mass effect, and promotes direct rebound mechanism resulting from a looser transition state geometry caused by varied electronegativity.

The Influence of Hydrogen Bonds on the Roaming Reaction

Roaming bypasses the conventional transition state and is a significant reaction pathway due to the unusual energy distributions of its products; however, its reaction pathway under external environmental interactions remains unclear. Herein, we report for the first time the roaming process of nitrobenzene, which is influenced by the hydrogen bonds (H-bonds) between nitro- and phenyl radicals and water molecules in the gas phase. Notably, despite the fact that the single water structure produces a higher but narrower barrier, whereas the double water structure leads to a lower but wider barrier, the roaming reaction still occurs. The underlying mechanism responsible for these influences of H-bonds is ascribed to the dramatically changed polarization and correlation interactions between the roaming radicals. The reaction rates and thermal perturbation probabilities are also remarkably influenced due to the presence of the H-bonds, by approximately 2 orders of magnitude. It is anticipated that this work will encourage the promising feasibility of introducing environmental molecules to modulate the roaming reaction.

A facile cost-effective electrolyte-assisted approach and comparative study towards the Greener synthesis of silica nanoparticles

Nowadays, silica nanoparticles are gaining tremendous importance because of their wide applications across different domains such as drug delivery, chromatography, biosensors, and chemosensors. The synthesis of silica nanoparticles generally requires a high percentage composition of organic solvent in an alkali medium. The eco-friendly synthesis of silica nanoparticles in bulk amounts can help save the environment and is cost-effective. Herein, efforts have been made to minimize the concentration of organic solvents used during synthesis via the addition of a low concentration of electrolytes, e.g., NaCl. The effects of electrolytes and solvent concentrations on nucleation kinetics, particle growth, and particle size were investigated. Ethanol was used as a solvent in various concentrations, ranging from 60% to 30%, and to optimize and validate the reaction conditions, isopropanol and methanol were also utilized as solvents. The concentration of aqua-soluble silica was determined using the molybdate assay to establish reaction kinetics, and this approach was also utilized to quantify the relative concentration changes in particles throughout the synthesis. The prime feature of the synthesis is the reduction in organic solvent usage by up to 50% using 68 mM NaCl. The surface zeta potential was reduced after the addition of an electrolyte, which made the condensation process faster and helped reaching the critical aggregation concentration in a shorter time. The effect of temperature was also monitored, and we obtained homogeneous and uniform nanoparticles by increasing the temperature. We found that it is possible to tune the size of the nanoparticles by changing the concentration of electrolytes and the temperature of the reaction using an eco-friendly approach. The overall cost of the synthesis can also be reduced by ∼35% by adding electrolytes.

Origin of the Microsolvation Effect on the Central Barriers of SN2 Reactions

We have quantitatively analyzed the microsolvation effect on the central barriers of microsolvated bimolecular nucleophilic substitution (S_N2) reactions by means of a two-step energy decomposition procedure. According to the first energy decompositions, an obvious increase in the central barrier for a microsolvated S_N2 reaction against its unsolvated counterpart can be mainly ascribed to the fact that the interaction between the solute and the conjunct solvent becomes less attractive from the reactant complex to the transition state. On the basis of the second energy decompositions with symmetry-adapted perturbation theory, this less attractive interaction in the transition state is primarily due to the interplay of the changes in the electrostatic, exchange, and induction components. However, the contribution of the change for the dispersion component is relatively small. A distinct linear correlation has also been observed between the changes of the total interaction energies and those of the corresponding electrostatic components for the microsolvated S_N2 reactions studied in this work. Moreover, the two-step energy decomposition procedure employed in this work is expected to be extensively applied to the gas phase reactions mediated by molecules or clusters.

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.

Investigating the competing E2 and SN2 mechanisms for the microsolvated HO-(H2O)n=0-4 + CH3CH2X (X = Cl, Br, I) reactions

We characterized the anti-E2, syn-E2, inv-S_N2, and ret-S_N2 reaction channels for the reaction of microsolvated HO^-(H₂O)_n anions with CH₃CH₂X (X = Cl, Br, I), using the CCSD(T)/PP/t//MP2/ECP/d level method, to understand how a solvent influences the competing E2 and S_N2 reactions. The calculated sequence of barrier for the four channels is ret-S_N2 > syn-E2 > anti-E2 > inv-S_N2. The barrier heights increase with incremental hydration as the system transfers from the gas phase to microsolvation, and to bulk solvation (using the PCM implicit solvent model). As the degree of hydration n increases, good correlations have been found between barrier heights and several thermodynamic, geometric and charge parameters, including the reaction enthalpy, proton/ethyl-cation affinity of the hydrated nucleophile, geometric looseness (%L^‡) and asymmetry (%AS^‡) and charge asymmetry (Δq(X-O)) of the transition structures. Under a molecular orbital scheme, the HOMOs of nucleophiles are stabilized by stepwise hydration, explaining the rise in the barriers. Considering the effect of the leaving group, the barrier heights exhibit linear correlation with the halogen electronegativity and H-acidity of substrate CH₃CH₂X. In terms of E2/S_N2 competition, the barrier difference, , first increases then decreases as the number of explicit water molecules increases, under both microsolvation and bulk solvation conditions, but the inv-S_N2 pathway is always favored over the anti-E2 pathway. Energy decomposition analysis attributes the increase of barrier difference to the greater geometric distortion in the anti-E2 transition structure.

Free-Energy Landscape of the SN2 Reaction CH3Br + Cl- → CH3Cl + Br- in Different Liquid Environments

This work describes in detail the reaction path of the well-known S_N2 reaction CH₃Br + Cl^- → CH₃Cl + Br^-, whose reaction rate has a huge variation with the solvent in the gas phase and in protic and aprotic liquid environments. We employed the ASEC-FEG method to optimize for minima (reactants and products) and saddle points (transition states) in the in-solution free-energy hypersurface. The method takes atomistic details of the solvent into account. A polarizable continuum model (PCM) has also been employed for comparison. The most perceptive structural changes are noted in aqueous solution by using the ASEC-FEG approach. The activation energies in all solvents, estimated by means of free-energy perturbation calculations, are in good agreement with the experimental data. The total solute-solvent hydrogen bonds play an important role in the increased barrier height observed in water and are therefore crucial to explain the huge decrease in the kinetic constant. It is also found that the hydration shell around the ions breaks itself spontaneously to accommodate the molecule, thus forming minimum energy complexes.

Proton transfer-induced competing product channels of microsolvated Y-(H2O)n + CH3I (Y = F, Cl, Br, I) reactions

The potential energy profiles of three proton transfer-involved product channels for the reactions of Y^-(H₂O)_1,2 + CH₃I (Y = F, Cl, Br, I) were characterized using the B97-1/ECP/d method. These three channels include the (1) PT_CH3 product channel that transfers a proton from methyl to nucleophile, (2) HO^--induced nucleophilic substitution (HO^--S_N2) product channel, and (3) oxide ion substitution (OIS) product channel that gives CH₃O^- and HY products. The reaction enthalpies and barrier heights follow the order OIS > PT_CH3 > HO^--S_N2 > Y^--S_N2, and thus HO^--S_N2 can compete with the most favored Y^--S_N2 product channel under singly-/doubly-hydrated conditions, while the PT_CH3 channel only occurs under high collision energy and the OIS channel is the least probable. All product channels share the same pre-reaction complex, Y^-(H₂O)_n-CH₃I, in the entrance of the potential energy profile, signifying the importance of the pre-reaction complex. For HO^-/Y^--S_N2 channels, we considered front-side attack, back-side attack, and halogen-bonded complex mechanisms. Incremental hydration increases the barriers of both HO^-/Y^--S_N2 channels as well as their barrier difference, implying that the HO^--S_N2 channel becomes less important when further hydrated. Varying the nucleophile Y^- from F^- to I^- also increases the barrier heights and barrier difference, which correlates with the proton affinity of the nucleophiles. Energy decomposition analyses show that both the orbital interaction energy and structural deformation energy of the transition states determine the S_N2 barrier change trend with incremental hydration and varying Y^-. In brief, this work computes the comprehensive potential energy surfaces of the HO^--S_N2 and PT_CH3 channels and shows how proton transfer affects the microsolvated Y^-(H₂O)_1,2 + CH₃I reaction by competing with the traditional Y^--S_N2 channel.

Effect of microsolvation on the mode specificity of the OH˙(H2O) + HCl reaction

The present study investigates the mode specificity in the microsolvated OH˙(H₂O) + HCl reaction using on-the-fly direct dynamics simulation. To the best of our knowledge, this is the first study which aims to gain insights into the effect of microsolvation on the mode selectivity. Our investigation reveals that, similar to the gas phase OH˙ + HCl reaction, the microsolvated reaction is also predominantly affected by the vibrational excitation of the HCl mode, whereas the OH vibrational mode behaves as a spectator. Interestingly, in contrast to the behavior of the bare reaction, the integral cross section at the ground state of the microsolvated reaction decreases with an increase in translational energy. However, for the vibrational excited states, the reactivity of the microsolvated reaction is found to be higher than that of the bare reaction within the selected range of translational energies.

Nonstatistical Reaction Dynamics

Nonstatistical dynamics is important for many chemical reactions. The Rice-Ramsperger-Kassel-Marcus (RRKM) theory of unimolecular kinetics assumes a reactant molecule maintains a statistical microcanonical ensemble of vibrational states during its dissociation so that its unimolecular dynamics are time independent. Such dynamics results when the reactant's atomic motion is chaotic or irregular. Intrinsic non-RRKM dynamics occurs when part of the reactant's phase space consists of quasiperiodic/regular motion and a bottleneck exists, so that the unimolecular rate constant is time dependent. Nonrandom excitation of a molecule may result in short-time apparent non-RRKM dynamics. For rotational activation, the 2J + 1 K levels for a particular J may be highly mixed, making K an active degree of freedom, or K may be a good quantum number and an adiabatic degree of freedom. Nonstatistical dynamics is often important for bimolecular reactions and their intermediates and for product-energy partitioning of bimolecular and unimolecular reactions. Post–transition state dynamics is often highly complex and nonstatistical.

Imaging Reaction Dynamics of F-(H2O) and Cl-(H2O) with CH3I

The dynamics of microhydrated nucleophilic substitution reactions have been studied using crossed beam velocity map imaging experiments and quasiclassical trajectory simulations at different collision energies between 0.3 and 2.6 eV. For F^–(H₂O) reacting with CH₃I, a small fraction of hydrated product ions I^–(H₂O) is observed at low collision energies. This product, as well as the dominant I^–, is formed predominantly through indirect reaction mechanisms. In contrast, a much smaller indirect fraction is determined for the unsolvated reaction. At the largest studied collision energies, the solvated reaction is found to also occur via a direct rebound mechanism. The measured product angular distributions exhibit an overall good agreement with the simulated angular distributions. Besides nucleophilic substitution, also ligand exchange reactions forming F^–(CH₃I) and, at high collision energies, proton transfer reactions are detected. The differential scattering images reveal that the Cl^–(H₂O) + CH₃I reaction also proceeds predominantly via indirect reaction mechanisms.

High-Level-Optimized Stationary Points for the F-(H2O) + CH3I System: Proposing a New Water-Induced Double-Inversion Pathway

We report 29 stationary points for the F^-(H₂O) + CH₃I reaction obtained by using the high-level explicitly correlated CCSD(T)-F12b method with the aug-cc-pVDZ basis set for the determination of the benchmark structures and frequencies and the aug-cc-pVQZ basis for energy computations. The stationary points characterize the monohydrated F^-- and OH^--induced Walden-inversion pathways and, for the first time, the front-side attack and F^--induced double-inversion mechanisms leading to CH₃F with retention as well as the novel H₂O-induced double-inversion retention pathway producing CH₃OH. Hydration effectively increases the relative energies of the stationary points, but the monohydrated inversion pathways are still barrierless, whereas the front-side attack and double-inversion barrier heights are around 30 and 20 kcal/mol, respectively.

Competition of F/OH-Induced SN2 and Proton-Transfer Reactions with Increased Solvation

The potential energy profiles of F/OH-induced nucleophilic substitution (S_N2) and proton-transfer (PT) channels evolving with solvation for reactions of F^-(H₂O) _n=1-2 + CH₃I were characterized using B3LYP/ECP/d method. The hydrogen-bonded F^-(H₂O) _n---HCH₂I prereaction complex at the entrance of potential energy surface (PES) has a significant role on the reaction dynamics for each channel. Among the above three channels, the F-S_N2 channel is the most preferred and OH-S_N2 could be competitive. In contrast, the PT channel will occur at much higher collision energy. Importantly, for each channel, the central barrier is gradually increased with the addition of water molecules. This phenomenon indicates that the reactivity will decrease with degrees of solvation and this has been confirmed by experiment and direct dynamics simulations. Moreover, compared with the previous trajectory simulations, a non-IRC behavior has been uncovered. The water delivering process from fluorine to iodine side as illustrated on PES is barely observed, and instead, the reaction tends to dehydrate before passing through the S_N2 barrier and proceeds with the less hydrated pathway in order to weaken the steric effect. The work presented here shows the comprehensive potential energy surfaces and structures information on the F-S_N2, PT, and OH-S_N2 channels, and predict their competitive relationship, which would be helpful for better understanding the dynamics behavior of the title and analogous reactions.

Nucleophilic Substitution (SN 2): Dependence on Nucleophile, Leaving Group, Central Atom, Substituents, and Solvent

The reaction potential energy surface (PES), and thus the mechanism of bimolecular nucleophilic substitution (SN 2), depends profoundly on the nature of the nucleophile and leaving group, but also on the central, electrophilic atom, its substituents, as well as on the medium in which the reaction takes place. Here, we provide an overview of recent studies and demonstrate how changes in any one of the aforementioned factors affect the SN 2 mechanism. One of the most striking effects is the transition from a double-well to a single-well PES when the central atom is changed from a second-period (e. g. carbon) to a higher-period element (e.g, silicon, germanium). Variations in nucleophilicity, leaving group ability, and bulky substituents around a second-row element central atom can then be exploited to change the single-well PES back into a double-well. Reversely, these variations can also be used to produce a single-well PES for second-period elements, for example, a stable pentavalent carbon species.

Taking the plunge: chemical reaction dynamics in liquids

The dynamics of chemical reactions in liquid solutions are now amenable to direct study using ultrafast laser spectroscopy techniques and advances in computer simulation methods. The surrounding solvent affects the chemical reaction dynamics in numerous ways, which include: (i) formation of complexes between reactants and solvent molecules; (ii) modifications to transition state energies and structures relative to the reactants and products; (iii) coupling between the motions of the reacting molecules and the solvent modes, and exchange of energy; (iv) solvent caging of reactants and products; and (v) structural changes to the solvation shells in response to the changing chemical identity of the solutes, on timescales which may be slower than the reactive events. This article reviews progress in the study of bimolecular chemical reaction dynamics in solution, concentrating on reactions which occur on ground electronic states. It illustrates this progress with reference to recent experimental and computational studies, and considers how the various ways in which a solvent affects the chemical reaction dynamics can be unravelled. Implications are considered for research in fields such as mechanistic synthetic chemistry.

Dynamic exit-channel pathways of the microsolvated HOO-(H2O) + CH3Cl SN2 reaction: Reaction mechanisms at the atomic level from direct chemical dynamics simulations

Microsolvated bimolecular nucleophilic substitution (S_N2) reaction of monohydrated hydrogen peroxide anion [HOO^-(H₂O)] with methyl chloride (CH₃Cl) has been investigated with direct chemical dynamics simulations at the M06-2X/6-31+G(d,p) level of theory. Dynamic exit-channel pathways and corresponding reaction mechanisms at the atomic level are revealed in detail. Accordingly, a product distribution of 0.85:0.15 is obtained for Cl^-:Cl^-(H₂O), which is consistent with a previous experiment [D. L. Thomsen et al. J. Am. Chem. Soc. 135, 15508 (2013)]. Compared with the HOO^- + CH₃Cl S_N2 reaction, indirect dynamic reaction mechanisms are enhanced by microsolvation for the HOO^-(H₂O) + CH₃Cl S_N2 reaction. On the basis of our simulations, further crossed molecular beam imaging experiments are highly suggested for the S_N2 reactions of HOO^- + CH₃Cl and HOO^-(H₂O) + CH₃Cl.

Imaging the dynamics of ion–molecule reactions

A range of ion–molecule reactions have been studied in the last years using the crossed-beam ion imaging technique, from charge transfer and proton transfer to nucleophilic substitution and elimination.