Electronic Structure Theory Study of the Microsolvated F(-)(H2O) + CH3I SN2 Reaction

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.

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.

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.

Characterization of a trans-trans Carbonic Acid-Fluoride Complex by Infrared Action Spectroscopy in Helium Nanodroplets

The high Lewis basicity and small ionic radius of fluoride promote the formation of strong ionic hydrogen bonds in the complexation of fluoride with protic molecules. Herein, we report that carbonic acid, a thermodynamically disfavored species that is challenging to investigate experimentally, forms a complex with fluoride in the gas phase. Intriguingly, this complex is highly stable and is observed in abundance upon nanoelectrospray ionization of an aqueous sodium fluoride solution in the presence of gas-phase carbon dioxide. We characterize the structure and properties of the carbonic acid–fluoride complex, F^–(H₂CO₃), and its deuterated isotopologue, F^–(D₂CO₃), by helium nanodroplet infrared action spectroscopy in the photon energy range of 390–2800 cm^–1. The complex adopts a C_2v symmetry structure with the carbonic acid in a planar trans–trans conformation and both OH groups forming ionic hydrogen bonds with the fluoride. Substantial vibrational anharmonic effects are observed in the infrared spectra, most notably a strong blue shift of the symmetric hydrogen stretching fundamental relative to predictions from the harmonic approximation or vibrational second-order perturbation theory. Ab initio thermostated ring-polymer molecular dynamics simulations indicate that this blue shift originates from strong coupling between the hydrogen stretching and bending vibrations, resulting in an effective weakening of the OH···F^– ionic hydrogen bonds.

Dynamics of Cl-(H2O) + CH3I Substitution Reaction: The Influences of Solvent and Nucleophile

The study of microsolvation provides a deeper understanding of solvent effects on reaction dynamics. Here, the properties of the S_N2 reaction of hydrated chloride with methyl iodide are investigated by direct dynamics simulations, and how the solute-solvent interactions and the basicity of nucleophiles can profoundly affect the atomic level dynamics is discussed in detail. The results show that the direct-rebound mechanism dominates the substitution reaction, and the roundabout mechanism, which prevails in the indirect unsolvated counterpart reaction, still accounts for a high proportion of the indirect mechanisms. The involvement of a solvent water molecule does not significantly reduce the cross section and rate constant compared to the unhydrated reaction at high collision energy. By varying solvated Cl^- to F^-, the dominant mechanisms are totally different and in contrast, the dynamics of water does not show much difference, and the departure of H₂O tends to occur prior to the substitution reaction because of the facile breakage of the hydrogen bond at high collision energy.

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.

Steric Effects of Solvent Molecules on SN2 Substitution Dynamics

Influences of solvent molecules on S_N2 reaction dynamics of microsolvated F^-(H₂O)_n with CH₃I, for n = 0-3, are uncovered by direct chemical dynamics simulations. The direct substitution mechanism, which is important without microsolvation, is quenched dramatically upon increasing hydration. The water molecules tend to force reactive encounters to proceed through the prereaction collision complex leading to indirect reaction. In contrast to F^-(H₂O), reaction with higher hydrated ions shows a strong propensity for ion desolvation in the entrance channel, diminishing steric hindrance for nucleophilic attack. Thus, nucleophilic substitution avoids the potential energy barrier with all of the solvent molecules intact and instead occurs through the less solvated barrier, which is energetically unexpected because the former barrier has a lower energy. The work presented here reveals a trade-off between reaction energetics and steric effects, with the latter found to be crucial in understanding how hydration influences microsolvated S_N2 dynamics.

Effects of microsolvation on a SN2 reaction: indirect atomistic dynamics and weakened suppression of reactivity

Systematic studies of microsolvation in the gas phase have enriched our knowledge of solvent effects. Here, the dynamics of a prototype S_N2 reaction of a hydrated fluoride ion with methyl iodide is uncovered employing direct dynamics simulations that show strikingly distinct features from those determined for an unsolvated system. An indirect scattering is found to prevail, which occurs dominantly by forming hydrated F^-(H₂O)-HCH₂I and F^-(H₂O)-CH₃I pre-reaction complexes at low energies, but proceeds through their water-free counterparts at higher energies. This finding is in strong contrast to a general evolution from indirect to direct dynamics with enhancing energy for the unsolvated substitution reactions, and this discrepancy is understood by the substantial steric hindrance introduced by a water molecule. As established in experiments, solvation suppresses the reactivity, whereas we find that this depression is remarkably frustrated upon raising the energy given that collision-induced dehydration essentially diminishes the water block for reactive collisions. The present study sheds light on how solute-solvent interactions affect the underlying dynamics at a deeper atomic level, thereby promoting our understanding of the fundamental solvent effects on chemical reactions in solution.

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.