How a Solvent Molecule Affects Competing Elimination and Substitution Dynamics. Insight into Mechanism Evolution with Increased Solvation

Competing C and N as Reactive Centers for Microsolvated Ambident Nucleophiles CN-(H2O)n=0-3: A Theoretical Study of E2/SN2 Reactions with CH3CH2X (X = Cl, Br, I)

As an ambident nucleophile, CN- has both C and N atoms that can act as the reactive center to facilitate substitution reactions. We investigate in detail the potential energy profiles of CN-(H2O)0-3 with CH3CH2X (X = Cl, Br, I) to explore the influence of solvent molecules on competition between the different nucleophilic atoms C and N involving the SN2 and E2 pathways. The energy barrier sequence for the transition states follows C@inv-SN2 < N@inv-SN2 < C@anti-E2 < N@anti-E2. When two different atoms act as nucleophilic atoms, the SN2 reaction is always preferred over the E2 reaction, and this preference increases with microsolvation. For the ambident nucleophiles CN-(H2O)0-3, C as the reactive center always has stronger nucleophilicity and basicity than N acting as the reactive center. Regarding the leaving group, the height of the energy barrier is positively correlated with the acidity of the CH3CH2X substrate for the E2 pathway and with X-heterolysis for the SN2 pathway. Furthermore, we found that in the gas phase, the energy barrier for different leaving group systems decreases gradually in the order Cl > Br > I, while in the SMD solution, the energy barrier and product energy increase slightly in the system from X = Cl to Br; this change may be due to the significantly weakened transition-state interaction for the X = Br system. Our activation strain, quantitative molecular orbital, and charge analyses reveal the physical mechanisms underlying the various computed trends. In addition, we also demonstrate the two points recently proposed by Vermeeren, P. . Chem. Eur. J. 2020, 26, 15538-15548.

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.

Computational Insights into SN 2 and Proton Transfer Reactions of CH3 O- with NH2 Y and CH3 Y

Bimolecular nucleophilic substitution (SN 2) reactions have been extensively studied in both theory and experiment. While research on C-centered SN 2 reactions (SN 2@C) has been ongoing, SN 2 reactions at neutral nitrogen (SN 2@N) have received increased attention in recent years. To recommend methods for dynamics simulations, the comparison for the properties of the geometries, vibrational frequencies, and energies is done between MP2 and six DFT functional calculations and experimental data as well as the high-level CCSD(T) method for CH3 O- +NH2 Cl/CH3 Cl reactions. The relative energy diagrams at the M06 method for CH3 O- with CH3 Y/NH2 Y reactions (Y=F, Cl, Br, I) in the gas and solution phase are explored to investigate the effects of the leaving groups, different reaction centers, and solvents. We mainly focus on the computational of inv-SN 2 and proton transfer (PT) pathways. The PT channel in the gas phase is more competitive than the SN 2 channel for N-center reactions, while the opposite is observed for C-centered reactions. Solvation completely inhibits the PT channel, making SN 2 the dominant pathway. Our study provides new insight into the SN 2 reaction mechanisms and rich the novel reaction model in gas-phase 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.

Competition between Elimination and Substitution for Ambident Nucleophiles CN- and Iodoethane Reactions in Gaseous and Aqueous Medium

Nucleophilic substitution (SN2) and elimination (E2) reactions between ambident nucleophiles have long been considered as typical reactions in organic chemistry, and exploring the competition between the two reactions is of great importance in chemical synthesis. As a nucleophile, CN- can use its C and N atoms as the reactive centers to undergo E2 and SN2 reactions, but related research is currently limited. This study uses the CCSD(T)/pp/t//MP2/ECP/d electronic structure method to perform detailed investigations on the potential energy profiles for SN2 and E2 reactions between CN- and CH3CH2I in gaseous and aqueous media. The potential energy profiles reveal that the energy barriers for SN2 and E2 reactions with the C atom as the reactive center are consistently lower than those with the N atom, indicating that the C atom has a stronger nucleophilic ability and stronger basicity. Furthermore, the potential energy profiles in both gas and aqueous environments show that the barriers of SN2 reactions are lower than those for E2 reactions with both C and N as the attacking atom. By using the frontier molecular orbital and activation strain models to explain the interesting phenomenon, the transition from the gas phase to solution was investigated, specifically in the gas-microsolvation-water transition. The results show that water molecules reduce the nucleophilicity and basicity of CN-, while strain energy (ΔEstrain) causes a greater increase in the energy barrier for E2 reactions. This study provides new insights and perspectives on the understanding of CN- as a nucleophile in SN2 reactions and serves as theoretical guidance for organic synthesis.

SN 2 versus E2 Competition of Cyclic Ethers

We have quantum chemically studied the influence of ring strain on the competition between the two mechanistically different SN 2 and E2 pathways using a series of archetypal ethers as substrate in combination with a diverse set of Lewis bases (F- , Cl- , Br- , HO- , H3 CO- , HS- , H3 CS- ), using relativistic density functional theory at ZORA-OLYP/QZ4P. The ring strain in the substrate is systematically increased on going from a model acyclic ether to a 6- to 5- to 4- to 3-membered ether ring. We have found that the activation energy of the SN 2 pathway sharply decreases when the ring strain of the system is increased, thus on going from large to small cyclic ethers, the SN 2 reactivity increases. In contrast, the activation energy of the E2 pathway generally rises along this same series, that is, from large to small cyclic ethers. The opposing reactivity trends induce a mechanistic switch in the preferred reaction pathway for strong Lewis bases from E2, for large cyclic substrates, to SN 2, for small cyclic substrates. Weak Lewis bases are unable to overcome the higher intrinsic distortivity of the E2 pathway and, therefore, always favor the less distortive SN 2 reaction.

Quasiclassical Trajectory Simulation as a Protocol to Build Locally Accurate Machine Learning Potentials

Direct trajectory calculations have become increasingly popular in recent computational chemistry investigations. However, the exorbitant computational cost of ab initio trajectory calculations usually limits its application in mechanistic explorations. Recently, machine learning-based potential energy surface (ML-PES) provides a powerful strategy to circumvent the heavy computational cost and meanwhile maintain the required accuracy. Despite the appealing potential, constructing a robust ML-PES is still challenging since the training set of the PES should cover a broad enough configuration space. In this work, we demonstrate that when the concerned properties could be collected by the localized sampling of the configuration space, quasiclassical trajectory (QCT) calculations can be invoked to efficiently obtain locally accurate ML-PESs. We prove our concept with two model reactions: methyl migration of i-pentane cation and dimerization of cyclopentadiene. We found that the locally accurate ML-PESs are sufficiently robust for reproducing the static and dynamic features of the reactions, including the time-resolved free energy and entropy changes, and time gaps.

Unsupervised Reaction Pathways Search for the Oxidation of Hypergolic Ionic Liquids: 1-Ethyl-3-methylimidazolium Cyanoborohydride (EMIM+/CBH-) as a Case Study

Hypergolic ionic liquids have come under increased study for having several desirable properties as a fuel source. One particular ionic liquid, 1-ethyl-3-methylimidazolium/cyanoborohydride (EMIM⁺/CBH^-), and oxidant, nitric acid (HNO₃), has been reported to be hypergolic experimentally, but its mechanism is not well-understood at a mechanistic level. In this computational study, the reaction is first probed with ab initio molecular dynamics simulations to confirm that anion-oxidant interactions likely are the first step in the mechanism. Second, the potential energy surface of the anion-oxidant system is studied with an in-depth search over possible isomerizations, and a network of possible intermediates are found. The critical point search is unsupervised and thus has the potential of identifying structures that deviate from chemical intuition. Molecular graphs are employed for analyzing 3000+ intermediates found, and nudged elastic band calculations are employed to identify transition states between them. Finally, the reactivity of the system is discussed through examination of minimal energy paths connecting the reactant to various common products from hypergolic ionic liquid oxidation. Eight products are reported for this system: NO, N₂O, NO₂, HNO, HONO, HNO₂, HCN, and H₂O. All reaction paths leading to these exothermic products have overall reaction barriers of 6-7 kcal/mol.

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.

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.

Microsolvated Ion-Molecule SN2 Reactions with Dual Nucleophiles Induced by Solvent Molecules

Singly-hydrated HOO - anion was found to induce alternative nucleophile HO - via proton transfer from water molecule as react with CH 3 Cl recently. To investigate the generality of this effect, the competition between the solvent-induced HO - -S N 2 pathway and the normal HOO - -S N 2 pathway is studied for the microsolvated HOO - (H 2 O) n=1,2,3 + CH 3 X (X = F, Cl, Br, I) reaction by quantum chemistry calculation. Incremental hydration increases the barrier heights of both pathways and enlarges the barrier difference between them, which favors the HOO - -S N 2 pathway. Interestingly, the barrier difference is insensitive to the leaving group. Calculation shows the water induced HO - -S N 2 pathway is highly suppressed as the degree of hydration increases beyond two. The differential barrier under incremental hydration can be explained by solvent molecules stabilizing the HOMO level of HO - (HOOH)(H 2 O) n-1 nucleophile more than that of HOO - (H 2 O) n nucleophile. Comparison between these HO - -nucleophiles and HOO - -nucleophiles suggests that α-effect exists. Activation strain analysis attributes the barrier differences to the stronger distortion of the TS of HO - -S N 2 pathway than the counterparts of HOO - -S N 2 pathway. This work adds our understanding of the role of individual solvent molecules to induce new nucleophiles of the fundamental organic reaction.

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.

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.

Nucleophilic Radiofluorination Using Tri-tert-Butanol Ammonium as a Bifunctional Organocatalyst: Mechanism and Energetics

We present a quantum chemical analysis of the 18F-fluorination of 1,3-ditosylpropane, promoted by a quaternary ammonium salt (tri-(tert-butanol)-methylammonium iodide (TBMA-I) with moderate to good radiochemical yields (RCYs), experimentally observed by Shinde et al. We obtained the mechanism of the SN2 process, focusing on the role of the -OH functional groups facilitating the reactions. We found that the counter-cation TBMA+ acts as a bifunctional promoter: the -OH groups function as a bidentate 'anchor' bridging the nucleophile [18F]F- and the -OTs leaving group or the third -OH. These electrostatic interactions cooperate for the formation of the transition states of a very compact configuration for facile SN2 18F-fluorination.

How Solvation Influences the SN2 versus E2 Competition

We have quantum chemically investigated how solvation influences the competition between the S_N2 and E2 pathways of the model F^– + C₂H₅Cl reaction. The system is solvated in a stepwise manner by going from the gas phase, then via microsolvation of one to three explicit solvent molecules, then last to bulk solvation using relativistic density functional theory at (COSMO)-ZORA-OLYP/QZ4P. We explain how and why the mechanistic pathway of the system shifts from E2 in the gas phase to S_N2 upon strong solvation of the Lewis base (i.e., nucleophile/protophile). The E2 pathway is preferred under weak solvation of the system by dichloromethane, whereas a switch in reactivity from E2 to S_N2 is observed under strong solvation by water. Our activation strain and Kohn–Sham molecular orbital analyses reveal that solvation of the Lewis base has a significant impact on the strength of the Lewis base. We show how strong solvation furnishes a weaker Lewis base that is unable to overcome the high characteristic distortivity associated with the E2 pathway, and thus the S_N2 pathway becomes viable.

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.

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.

Single Solvent Molecules Induce Dual Nucleophiles in Gas-Phase Ion-Molecule Nucleophilic Substitution Reactions

Direct dynamics simulation of singly hydrated peroxide ion reacting with CH₃Cl reveals a new product channel that forms CH₃OH + Cl^- + HOOH, besides the traditional channel that forms CH₃OOH + Cl^- + H₂O. This finding shows that singly hydrated peroxide ion behaves as a dual nucleophile through proton transfer between HOO^-(H₂O) and HO^-(HOOH). Trajectory analysis attributes the occurrence of the thermodynamically and kinetically unfavored HO^--induced pathway to the entrance channel dynamics, where extensive proton transfer occurs within the deep well of the prereaction complex. This study represents the first example of a single solvent molecule altering the nucleophile in a gas-phase ion-molecule nucleophilic substitution reaction, in addition to reducing the reactivity and affecting the dynamics, signifying the importance of dynamical effects of solvent molecules.

Resolving Entangled Reactivity Modes through External Electric Fields and Substitution: Application to E2/SN2 Reactions

In this study, we explore strategies to resolve entangled reactivity modes. More specifically, we consider the competition between S_N2 and E₂ reaction pathways for alkyl halides and nucleophiles/bases. We first demonstrate that the emergence of an E₂-preference is associated with an enhancement of the magnitude of the resonance stabilization in the transition-state (TS) region, resulting from the improved mixing of electrostatically stabilized valence bond structures into the TS wavefunction. Subsequently, we show that the TS resonance energy can be tuned selectively and rationally either through the application of an oriented external electric field directed along the C-C axis of the alkyl halide or through a regular substitution approach of the C-C moiety. We end our study by demonstrating that the insights gained from our analysis enable one to rationalize the main reactivity trends emerging from a recently published large database of competing S_N2 and E₂ reaction pathways.

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.

Ab Initio Molecular Dynamics Simulations of the SN1/SN2 Mechanistic Continuum in Glycosylation Reactions

We report a computational approach to evaluate the reaction mechanisms of glycosylation using ab initio molecular dynamics (AIMD) simulations in explicit solvent. The reaction pathways are simulated via free energy calculations based on metadynamics and trajectory simulations using Born-Oppenheimer molecular dynamics. We applied this approach to investigate the mechanisms of the glycosylation of glucosyl α-trichloroacetimidate with three acceptors (EtOH, i-PrOH, and t-BuOH) in three solvents (ACN, DCM, and MTBE). The reactants and the solvents are treated explicitly using density functional theory. We show that the profile of the free energy surface, the synchronicity of the transition state structure, and the time gap between leaving group dissociation and nucleophile association can be used as three complementary indicators to describe the glycosylation mechanism within the SN1/SN2 continuum for a given reaction. This approach provides a reliable means to rationalize and predict reaction mechanisms and to estimate lifetimes of oxocarbenium intermediates and their dependence on the glycosyl donor, acceptor, and solvent environment.

Investigating the role of halogen-bonded complexes in microsolvated Y−(H2O)n + CH3I SN2 reactions

A halogen-bonded complex pathway is computed for Y⁻(H₂O)_n + CH₃I (Y = HO, F, Cl, Br, and I) ion–molecule nucleophilic substitution reactions and is compared with back-side and front-side attack pathways.

Tackling Solvent Effects by Coupling Electronic and Molecular Density Functional Theory

Solvation effects can have a tremendous influence on chemical reactions. However, precise quantum chemistry calculations are most often done either in vacuum neglecting the role of the solvent or using continuum solvent model ignoring its molecular nature. We propose a new method coupling a quantum description of the solute using electronic density functional theory with a classical grand-canonical treatment of the solvent using molecular density functional theory. Unlike a previous work, both densities are minimized self-consistently, accounting for mutual polarization of the molecular solvent and the solute. The electrostatic interaction is accounted using the full electron density of the solute rather than fitted point charges. The introduced methodology represents a good compromise between the two main strategies to tackle solvation effects in quantum calculation. It is computationally more effective than a direct quantum mechanics/molecular mechanics coupling, requiring the exploration of many solvent configurations. Compared to continuum methods, it retains the full molecular-level description of the solvent. We validate this new framework onto two usual benchmark systems: a water solvated in water and the symmetrical nucleophilic substitution between chloromethane and chloride in water. The prediction for the free energy profiles are not yet fully quantitative compared to experimental data, but the most important features are qualitatively recovered. The method provides a detailed molecular picture of the evolution of the solvent structure along the reaction pathway.

Role of hydrogen bond capacity of solvents in reactions of amines with CO2: A computational study

Various computational methods were employed to investigate the zwitterion formation, a critical step for the reaction of monoethanolamine with CO₂, in five solvents (water, monoethanolamine, propylamine, methanol and chloroform) to probe the effect of hydrogen bond capacity of solvents on the reaction of amine with CO₂ occurring in the amine-based CO₂ capture process. The results indicate that the zwitterion can be formed in all considered solvents except chloroform. For two pairs of solvents (methanol and monoethanolamine, propylamine and chloroform) with similar dielectric constant but different hydrogen bond capacity, the solvents with higher hydrogen bond capacity (monoethanolamine and propylamine) facilitate the zwitterion formation. More importantly, kinetics parameters such as activation free energy for the zwitterion formation are more relevant to the hydrogen bond capacity than to dielectric constant of the considered solvents, clarifying the hydrogen bond capacity could be more important than dielectric constant in determining the kinetics of monoethanolamine with CO₂.

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.

Synthesis of bodinieric acids A and B, both C-18 and C-19-functionalized abietane diterpenoids: DFT study of the key aldol reaction

The first synthesis of C-18- and C-19-bifunctionalized abietane diterpenoids, bodinieric (or callicapoic) acids, via an aldol reaction has been developed. This key aldol reaction was very sensitive to steric hindrance. This fact has been studied by deuterium exchange experiments and DFT methods. Optimization of this reaction led to the synthesis of anti-inflammatory bodinieric acids A and B, starting from abietic acid.

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.

Structural, kinetics and mechanistic studies of transfer hydrogenation of ketones catalyzed by chiral (pyridyl)imine nickel(ii) complexes

The chiral synthons (S-)-1-phenyl-N-(pyridine-2-yl)ethylidine)ethanamine (L1), (R-)-1phenyl-N-(pyridine-2-yl)ethylidine))ethanamine (L2) (S)-1-phenyl-N-(pyridine-2-yl methylene) ethanamine (L3), and (R)-1-phenyl-N-(pyridine-2-yl methylene) ethanamine (L4) were synthesized in good yields. Treatments of L1-L4 with NiBr₂(DME) and NiCl₂ precursor afforded dinuclear complexes [Ni₂(L1)₄-μ-Br₂]NiBr₄ (Ni1), [Ni₂(L2)₄-μ-Br₂]NiBr₄ (Ni2), [Ni₂(L3)₄-μBr₂]Br₂ (Ni3), [Ni₂(L4)₄-μ-Br₂]NiBr₄ (Ni4) and [Ni(L4)₂Cl₂] (Ni5). The identities of the compounds were established using NMR, FT-IR and EPR spectroscopy, mass spectrometry, magnetic moments, elemental analysis and single crystal X-ray crystallography. The dinuclear dibromide nickel complexes dissociate into mononuclear species in the presence of strongly coordinating solvents. Compounds Ni1-Ni5 displayed moderate catalytic activities in the asymmetric transfer hydrogenation (ATH) of ketones, but with low enantiomeric excess (ee%). Both mercury and substoichiometric poisoning tests pointed to the homogeneous nature of the active species with the partial formation of catalytically active Ni(0) nanoparticles. Low resolution mass spectrometry analyses of the intermediates supported a dihydride mechanistic pathway for the transfer of hydrogenation reactions.

EVB and polarizable MM study of energy relaxation in fluorine–acetonitrile reactions

Many-body effects can impact on rates of energy transfer from a ‘hot’ DF solute to acetonitrile solvent.

A Comparison of QM/MM Simulations with and without the Drude Oscillator Model Based on Hydration Free Energies of Simple Solutes

Maintaining a proper balance between specific intermolecular interactions and non-specific solvent interactions is of critical importance in molecular simulations, especially when predicting binding affinities or reaction rates in the condensed phase. The most rigorous metric for characterizing solvent affinity are solvation free energies, which correspond to a transfer from the gas phase into solution. Due to the drastic change of the electrostatic environment during this process, it is also a stringent test of polarization response in the model. Here, we employ both the CHARMM fixed charge and polarizable force fields to predict hydration free energies of twelve simple solutes. The resulting classical ensembles are then reweighted to obtain QM/MM hydration free energies using a variety of QM methods, including MP2, Hartree⁻Fock, density functional methods (BLYP, B3LYP, M06-2X) and semi-empirical methods (OM2 and AM1 ). Our simulations test the compatibility of quantum-mechanical methods with molecular-mechanical water models and solute Lennard⁻Jones parameters. In all cases, the resulting QM/MM hydration free energies were inferior to purely classical results, with the QM/MM Drude force field predictions being only marginally better than the QM/MM fixed charge results. In addition, the QM/MM results for different quantum methods are highly divergent, with almost inverted trends for polarizable and fixed charge water models. While this does not necessarily imply deficiencies in the QM models themselves, it underscores the need to develop consistent and balanced QM/MM interactions. Both the QM and the MM component of a QM/MM simulation have to match, in order to avoid artifacts due to biased solute⁻solvent interactions. Finally, we discuss strategies to improve the convergence and efficiency of multi-scale free energy simulations by automatically adapting the molecular-mechanics force field to the target quantum method.