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

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

Reaction dynamics of the methoxy anion CH3O- with methyl iodide CH3I

Studying larger nucleophiles in bimolecular nucleophilic substitution (SN2) reactions bridges the gap from simple model systems to those relevant to organic chemistry. Therefore, we investigated the reaction dynamics between the methoxy anion (CH3O-) and iodomethane (CH3I) in our crossed-beam setup combined with velocity map imaging at the four collision energies 0.4, 0.7, 1.2, and 1.6 eV. We find the two ionic products I- and CH2I-, which can be attributed to the SN2 and proton transfer channels, respectively. The proton transfer channel progresses in a previously observed fashion from indirect to direct scattering with increasing collision energy. Interestingly, the SN2 channel exhibits direct dynamics already at low collision energies. Both the direct stripping, leading to forward scattering, and the direct rebound mechanism, leading to backward scattering into high angles, are observed.

Dynamics of carbene formation in the reaction of methane with the tantalum cation in the gas phase

The controlled activation of methane has drawn significant attention throughout various disciplines over the last few decades. In gas-phase experiments, the use of model systems with reduced complexity compared to condensed-phase catalytic systems allows us to investigate the intrinsic reactivity of elementary reactions down to the atomic level. Methane is rather inert in chemical reactions, as the weakening or cleavage of a C-H bond is required to make use of methane as C1-building block. The simplest model system for transition-metal-based catalysts is a mono-atomic metal ion. Only a few atomic transition-metal cations activate methane at room temperature. One of the most efficient elements is tantalum, which forms a carbene and releases molecular hydrogen in the reaction with methane: Ta+ + CH4 → TaCH2+ + H2. The reaction takes place at room temperature due to efficient intersystem crossing from the quintet to the triplet surface, i.e., from the electronic ground state of the tantalum cation to the triplet ground state of the tantalum carbene. This multi-state reactivity is often seen for reactions involving transition-metal centres, but leads to their theoretical treatment being a challenge even today. Chemical reactions, or to be precise reactive collisions, are dynamic processes making their description even more of a challenge to experiment and theory alike. Experimental energy- and angle-differential cross sections allow us to probe the rearrangement of atoms during a reactive collision. By interpreting the scattering signatures, we gain insight into the atomistic mechanisms and can move beyond stationary descriptions. Here, we present a study combining collision energy dependent experimentally measured differential cross sections with ab initio calculations of the minimum energy pathway. Product ion velocity distributions were recorded using our crossed-beam velocity map imaging experiment dedicated to studying transition-metal ion molecule reactions. TaCH2+ velocity distributions reveal a significant degree of indirect dynamics. However, the scattering distributions also show signatures of rebound dynamics. We compare the present results to the oxygen transfer reaction between Ta+ and carbon dioxide, which we recently studied.

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.

Understanding the Formation and Growth of New Atmospheric Particles at the Molecular Level through Laboratory Molecular Beam Experiments

Atmospheric new particle formation (NPF), which exerts comprehensive implications for climate, air quality and human health, has received extensive attention. From molecule to cluster is the initial and most important stage of the nucleation process of atmospheric new particles. However, due to the complexity of the nucleation process and limitations of experimental characterization techniques, there is still a great uncertainty in understanding the nucleation mechanism at the molecular level. Laboratory-based molecular beam methods can experimentally implement the generation and growth of typical atmospheric gas-phase nucleation precursors to nanoscale clusters, characterize the key physical and chemical properties of clusters such as structure and composition, and obtain a series of their physicochemical parameters, including association rate coefficients, electron binding energy, pickup cross section and pickup probability and so on. These parameters can quantitatively illustrate the physicochemical properties of the cluster, and evaluate the effect of different gas phase nucleation precursors on the formation and growth of atmospheric new particles. We review the present literatures on atmospheric cluster formation and reaction employing the experimental method of laboratory molecular beam. The experimental apparatuses were classified and summarized from three aspects of cluster generation, growth and detection processes. Focus of this review is on the properties of nucleation clusters involving different precursor molecules of water, sulfuric acid, nitric acid and NxOy, respectively. We hope this review will provide a deep insight for effects of cluster physicochemical properties on nucleation, and reveal the formation and growth mechanism of atmospheric new particle at the molecular level.

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.

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.

Quasi-classical trajectory study of the OH- + CH3I reaction: theory meets experiment

Regarding OH^- + CH₃I, several studies have focused on the dynamics of the reaction. Here, high-level quasi-classical trajectory simulations are carried out at four different collision energies on our recently developed potential energy surface. In all, more than half a million trajectories are performed, and for the first time, the detailed quasi-classical trajectory results are compared with the reanalysed crossed-beam ion imaging experiments. Concerning the previously reported direct dynamics study of OH^- + CH₃I, a better agreement can be obtained between the revised experiment and our novel theoretical results. Furthermore, in the present work, the benchmark geometries, frequencies and relative energies of the stationary points are also determined for the OH^- + CH₃I proton-abstraction channel along with the earlier characterized S_N2 channel.

Imaging the Atomistic Dynamics of Single Proton Transfer and Combined Hydrogen/Proton Transfer in the O- + CH3I Reaction

We report on reactive scattering studies of the proton transfer and combined hydrogen/proton transfer in the O^- + CH₃I reaction. We combine state-of-the-art crossed-beam velocity map imaging and quantum chemistry calculations to understand the dynamics for the formations of the CH₂I^- + OH and CHI^- + H₂O products. The experimental velocity- and angle-differential cross section show for both products and at all collision energies (between 0.3 and 2.0 eV) that the product ions are predominantly forward scattered. For the CHI^- + H₂O channel, the data show lower product velocities, indicative of higher internal excitation, than in the case of single proton transfer. Furthermore, our results suggest that the combined hydrogen/proton transfer proceeds via a two-step process: In the first step, O^- abstracts one H atom to form OH^-, and then the transient OH^- removes an additional proton from CH₂I to form the energetically stable H₂O coproduct.

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.

Reaction mechanism of an intracluster SN2 reaction induced by electron capture

Bimolecular nucleophilic substitution (S_N2) reactions have been widely investigated from both experimental and theoretical points of view because they represent one of the simplest organic reactions. Most studies on S_N2 reactions have been focused on bimolecular collision. In contrast, information on intracluster S_N2 reactions is limited. In this study, an intracluster S_N2 reaction of NF₃-CH₃Cl triggered by electron attachment was investigated using a direct ab initio molecular dynamics (AIMD) method. In the structure of NF₃-CH₃Cl, the N-F bond in NF₃ is oriented collinearly toward the carbon atom of CH₃Cl. After electron capture by NF₃-CH₃Cl, the F^- ion that is generated from the (NF₃)^- moiety collides with the carbon atom of CH₃Cl. The intracluster S_N2 reaction occurs as follows: (NF₃-CH₃Cl)^- (electron capture state) → NF₂-(F^-)-CH₃Cl (pre-reaction complex) → transition state (TS) → NF₂-CH₃F-Cl^- (post-reaction complex) → NF₂ + CH₃F + Cl^- (product state). The reaction energy is efficiently transferred to the translational mode of Cl^-, and the Cl^- ion with a high translational energy is then removed from the system. This energy is significantly larger than that of Cl^- formed in the bimolecular S_N2 reaction (F^- + CH₃Cl). The reaction mechanism is discussed based on the theoretical results.

A photoionized pulsed low-energy ion beam source for quantum state-to-state crossed ion-molecule scattering

Here, we report the design and test of a pulsed low-energy ion beam source for crossed ion-molecule scattering studies. The ions are produced by laser photoionization based methods and thus can be prepared in well-defined quantum states. By using the combination of a double Einzel lenses setup and a specially designed shielding tube, a well spatially confined ion bunch with tunable kinetic energies in the range of 1.0-5.0 eV and typical spreads of ∼150 meV (full width at half maximum) can be formed in the center of a velocity-map imaging (VMI) stack. By combining it with a recently constructed three-dimensional VMI system, the present apparatus is readily available for quantum state-to-state crossed ion-molecule scattering studies.

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.

Fifty years of nucleophilic substitution in the gas phase

Bimolecular nucleophilic substitution ( S N 2 ) reactions have become a model system for the investigation of structure-reactivity relationships, stereochemistry, solvent influences, and detailed atomistic dynamics. In this review, the progress during five decades of experimental and theoretical research on gas phase S N 2 reactions is discussed. Many advancements of the employed methods have led to a tremendous increase in our understanding of the properties and the dynamics of these reactions. For reactions involving six atoms a quantitative agreement of the differential reactive scattering cross sections has already been achieved, in the future it is expected that even larger polyatomic reactions systems become tractable. Furthermore, studies with higher precision, improved reactant control, and a more accurate theoretical treatment of quantum effects are envisioned.

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.

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.

Photoelectron-Photofragment Coincidence Spectroscopy With Ions Prepared in a Cryogenic Octopole Accumulation Trap: Collisional Excitation and Buffer Gas Cooling

A cryogenic octopole accumulation trap (COAT) has been coupled to a photoelectron-photofragment coincidence (PPC) spectrometer allowing for improved control over anion vibrational excitation. The anions are heated and cooled via collisions with buffer gas <17 K. Shorter trapping times (500 μs) prevent thermalization and result in anions with high internal excitation while longer trapping times (80 ms) at cryogenic temperatures thermalize the ions to the temperature of the buffer gas. The capabilities of the COAT are demonstrated using PPC spectroscopy ofO 3 - at 388 nm (Ehν = 3.20 eV). Cooling the precursor anions with COAT resulted in the elimination of the autodetachment of vibrationally excitedO 2 - produced by the photodissociationO 3 - + hν → O +O 2 - (v ≥ 4). Under heating conditions, a lower limit temperature for the anions was determined to be 1,500 K through Franck-Condon simulations of the photodetachment spectrum ofO 3 - , considering a significant fraction of the ions undergo photodissociation in competition with photodetachment. The ability to cool or heat ions by varying ion injection and trapping duration in COAT provides a new flexibility for studying the spectroscopy of cold ions as well as thermally activated processes.

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.

Mass spectrometry of aerosol particle analogues in molecular beam experiments

Nanometer-size particles such as ultrafine aerosol particles, ice nanoparticles, water nanodroplets, etc, play an important, however, not yet fully understood role in the atmospheric chemistry and physics. These species are often composed of water with admixture of other atmospherically relevant molecules. To mimic and investigate such particles in laboratory experiments, mixed water clusters with atmospherically relevant molecules can be generated in molecular beams and studied by various mass spectrometric methods. The present review demonstrates that such experiments can provide unprecedented details of reaction mechanisms, and detailed insight into the photon-, electron-, and ion-induced processes relevant to the atmospheric chemistry. After a brief outline of the molecular beam preparation, cluster properties, and ionization methods, we focus on the mixed clusters with various atmospheric molecules, such as hydrated sulfuric acid and nitric acid clusters, N_x O_y and halogen-containing molecules with water. A special attention is paid to their reactivity and solvent effects of water molecules on the observed processes.

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.

Imaging Proton Transfer and Dihalide Formation Pathways in Reactions of F(-) + CH3I

Ion–molecule reactions of the type X^– + CH₃Y are commonly assumed to produce Y^– through bimolecular nucleophilic substitution (S_N2). Beyond this reaction, additional reaction products have been observed throughout the last decades and have been ascribed to different entrance channel geometries differing from the commonly assumed collinear approach. We have performed a crossed beam velocity map imaging experiment on the F^– + CH₃I reaction at different relative collision energies between 0.4 and 2.9 eV. We find three additional channels competing with nucleophilic substitution at high energies. Experimental branching ratios and angle- and energy differential cross sections are presented for each product channel. The proton transfer product CH₂I^– is the main reaction channel, which competes with nucleophilic substitution up to 2.9 eV relative collision energy. At this level, the second additional channel, the formation of IF^– via halogen abstraction, becomes more efficient. In addition, we present the first evidence for an [FHI]⁻ product ion. This [FHI]⁻ product ion is present only for a narrow range of collision energies, indicating possible dissociation at high energies. All three products show a similar trend with respect to their velocity- and scattering angle distributions, with isotropic scattering and forward scattering of the product ions occurring at low and high energies, respectively. Reactions leading to all three reaction channels present a considerable amount of energy partitioning in product internal excitation. The internally excited fraction shows a collision energy dependence only for CH₂I^–. A similar trend is observed for the isoelectronic OH^– + CH₃I system. The comparison of our experimental data at 1.55 eV collision energy with a recent theoretical calculation for the same system shows a slightly higher fraction of internal excitation than predicted, which is, however, compatible within the experimental accuracy.

Buffer-Gas Cooling of a Single Ion in a Multipole Radio Frequency Trap Beyond the Critical Mass Ratio

We theoretically investigate the dynamics of a trapped ion immersed in a spatially localized buffer gas. For a homogeneous buffer gas, the ion's energy distribution reaches a stable equilibrium only if the mass of the buffer gas atoms is below a critical value. This limitation can be overcome by using multipole traps in combination with a spatially confined buffer gas. Using a generalized model for elastic collisions of the ion with the buffer-gas atoms, the ion's energy distribution is numerically determined for arbitrary buffer-gas distributions and trap parameters. Three regimes characterized by the respective analytic form of the ion's equilibrium energy distribution are found. Final ion temperatures down to the millikelvin regime can be achieved by adiabatically decreasing the spatial extension of the buffer gas and the effective ion trap depth (forced sympathetic cooling).

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

The potential energy profile of microhydrated fluorine ion reaction with methyl iodine has been characterized by extensive electronic structure calculations. Both hydrogen-bonded F(-)(H2O)---HCH2I and ion-dipole F(-)(H2O)---CH3I complexes are formed for the reaction entrance and the PES in vicinity of these complexes is very flat, which may have important implications for the reaction dynamics. The water molecule remains on the fluorine side until the reactive system goes to the SN2 saddle point. It can easily move to the iodine side with little barrier, but in a nonsynchronous reaction path after the dynamical bottleneck to the reaction, which supports the previous prediction for microsolvated SN2 systems. The influence of solvating water molecule on the reaction mechanism is probed by comparing with the influence of the nonsolvated analogue and other microsolvated SN2 systems. Taking the CCSD(T) single-point calculations based on MP2-optimized geometries as benchmark, the DFT functionals B97-1 and B3LYP are found to better characterize the potential energy profile for the title reaction and are recommended as the preferred methods for the direct dynamics simulations to uncover the dynamic behaviors.

Chemical dynamics simulations of the monohydrated OH(-)(H2O) + CH3I reaction. Atomic-level mechanisms and comparison with experiment

Direct dynamics simulations, with B97-1/ECP/d theory, were performed to study the role of microsolvation for the OH(-)(H2O) + CH3I reaction. The SN2 reaction dominates at all reactant collision energies, but at higher collision energies proton transfer to form CH2I(-), and to a lesser extent CH2I(-) (H2O), becomes important. The SN2 reaction occurs by direct rebound and stripping mechanisms, and 28 different indirect atomistic mechanisms, with the latter dominating. Important components of the indirect mechanisms are the roundabout and formation of SN2 and proton transfer pre-reaction complexes and intermediates, including [CH3--I--OH](-). In contrast, for the unsolvated OH(-) + CH3I SN2 reaction, there are only seven indirect atomistic mechanisms and the direct mechanisms dominate. Overall, the simulation results for the OH(-)(H2O) + CH3I SN2 reaction are in good agreement with experiment with respect to reaction rate constant, product branching ratio, etc. Differences between simulation and experiment are present for the SN2 velocity scattering angle at high collision energies and the proton transfer probability at low collision energies. Equilibrium solvation by the H2O molecule is unimportant. The SN2 reaction is dominated by events in which H2O leaves the reactive system as CH3OH is formed or before CH3OH formation. Formation of solvated products is unimportant and participation of the (H2O)CH3OH---I(-) post-reaction complex for the SN2 reaction is negligible.

Comparison of direct dynamics simulations with different electronic structure methods. F(-) + CH3I with MP2 and DFT/B97-1

In previous work, ion imaging experiments and direct chemical dynamics simulations with DFT/B97-1 were performed to study the atomic-level dynamics of the F(-) + CH3I → FCH3 + I(-) SN2 nucleophilic substitution reaction at different collision energies. Overall, the simulations are in quite good agreement with experiment at the low collision energy of 0.32 eV, however there are differences between experiment and simulation at the high collision energy of 1.53 eV. A recent CCSD(T) study of the potential energy surface for the F(-) + CH3I → FCH3 + I(-) SN2 reaction shows that it has both a traditional C3v and a hydrogen-bond entrance channel. They are represented by MP2 but not by B97-1, which has only the latter channel. On the other hand, B97-1 gives a reaction exothermicity in excellent agreement with experiment, while MP2 is in error by 24.3 kJ mol(-1). In the work presented here, direct dynamics simulations using MP2/aug-cc-pvdz/ECP/d were performed for the F(-) + CH3I → FCH3 + I(-) reaction at a 1.53 eV collision energy. The same direct rebound and stripping and indirect atomistic reaction mechanisms are found in the B97-1 and MP2 simulations. Both the B97-1 and MP2 simulations agree with the experimental fraction of the available product energy partitioned to CH3F internal energy, i.e. fint = 0.59 ± 0.08. However, the MP2 fint distribution is broader and in better agreement with experiment than B97-1. The two simulations methods give the same product energy partitioning for the stripping mechanism, but different partitionings for the rebound and indirect mechanisms. Compared to experiment, the principal difference between the B97-1 and MP2 results is the differential cross section which is nearly isotropic for B97-1. For MP2 backward scattering is more important than forward, as found in the experiments. Though there is no overall barrier for the reaction, B97-1 gives a reaction cross section appreciably larger than that for MP2, i.e. 8.6 ± 2.2 Å(2)versus 1.8 ± 0.3 Å(2). For B97-1 59% of the reaction consists of indirect mechanisms, while for MP2 the indirect mechanisms only comprise 11% of the reaction. The experimental differential cross section is more consistent with the atomistic mechanisms for MP2 than for B97-1.

Identification of atomic-level mechanisms for gas-phase X- + CH3Y SN2 reactions by combined experiments and simulations

For the traditional model of gas-phase X(-) + CH3Y SN2 reactions, C3v ion-dipole pre- and postreaction complexes X(-)---CH3Y and XCH3---Y(-), separated by a central barrier, are formed. Statistical intramolecular dynamics are assumed for these complexes, so that their unimolecular rate constants are given by RRKM theory. Both previous simulations and experiments have shown that the dynamics of these complexes are not statistical and of interest is how these nonstatistical dynamics affect the SN2 rate constant. This work also found there was a transition from an indirect, nonstatistical, complex forming mechanism, to a direct mechanism, as either the vibrational and/or relative translational energy of the reactants was increased. The current Account reviews recent collaborative studies involving molecular beam ion-imaging experiments and direct (on-the-fly) dynamics simulations of the SN2 reactions for which Cl(-), F(-), and OH(-) react with CH3I. Also considered are reactions of the microsolvated anions OH(-)(H2O) and OH(-)(H2O)2 with CH3I. These studies have provided a detailed understanding of the atomistic mechanisms for these SN2 reactions. Overall, the atomistic dynamics for the Cl(-) + CH3I SN2 reaction follows those found in previous studies. The reaction is indirect, complex forming at low reactant collision energies, and then there is a transition to direct reaction between 0.2 and 0.4 eV. The direct reaction may occur by rebound mechanism, in which the ClCH3 product rebounds backward from the I(-) product or a stripping mechanism in which Cl(-) strips CH3 from the I atom and scatters in the forward direction. A similar indirect to direct mechanistic transition was observed in previous work for the Cl(-) + CH3Cl and Cl(-) + CH3Br SN2 reactions. At the high collision energy of 1.9 eV, a new indirect mechanism, called the roundabout, was discovered. For the F(-) + CH3I reaction, there is not a transition from indirect to direct reaction as Erel is increased. The indirect mechanism, with prereaction complex formation, is important at all the Erel investigated, contributing up ∼60% of the reaction. The remaining direct reaction occurs by the rebound and stripping mechanisms. Though the potential energy curve for the OH(-) + CH3I reaction is similar to that for F(-) + CH3I, the two reactions have different dynamics. They are akin, in that for both there is not a transition from an indirect to direct reaction. However, for F(-) + CH3I indirect reaction dominates at all Erel, but it is less important for OH(-) + CH3I and becomes negligible as Erel is increased. Stripping is a minor channel for F(-) + CH3I, but accounts for more than 60% of the OH(-) + CH3I reaction at high Erel. Adding one or two H2O molecules to OH(-) alters the reaction dynamics from that for unsolvated OH(-). Adding one H2O molecule enhances indirect reaction at low Erel, and changes the reaction mechanism from primarily stripping to rebound at high Erel. With two H2O molecules the dynamics is indirect and isotropic at all collision energies.

The α-effect and competing mechanisms: the gas-phase reactions of microsolvated anions with methyl formate

The enhanced reactivity of α-nucleophiles, which contain an electron lone pair adjacent to the reactive site, has been demonstrated in solution and in the gas phase and, recently, for the gas-phase S(N)2 reactions of the microsolvated HOO(-)(H2O) ion with methyl chloride. In the present work, we continue to explore the significance of microsolvation on the α-effect as we compare the gas-phase reactivity of the microsolvated α-nucleophile HOO(-)(H2O) with that of microsolvated normal alkoxy nucleophiles, RO(-)(H2O), in reactions with methyl formate, where three competing reactions are possible. The results reveal enhanced reactivity of HOO(-)(H2O) towards methyl formate, and clearly demonstrate the presence of an overall α-effect for the reactions of the microsolvated α-nucleophile. The association of the nucleophiles with a single water molecule significantly lowers the degree of proton abstraction and increases the S(N)2 and B(AC)2 reactivity compared with the unsolvated analogs. HOO(-)(H2O) reacts with methyl formate exclusively via the B(AC)2 channel. While microsolvation lowers the overall reaction efficiency, it enhances the B(AC)2 reaction efficiency for all anions compared with the unsolvated analogs. This may be explained by participation of the solvent water molecule in the B(AC)2 reaction in a way that continuously stabilizes the negative charge throughout the reaction.

Investigating the α-effect in gas-phase S(N)2 reactions of microsolvated anions

The α-effect-enhanced reactivity of nucleophiles with a lone-pair adjacent to the attacking center-was recently demonstrated for gas-phase S(N)2 reactions of HOO(-), supporting an intrinsic component of the α-effect. In the present work we explore the gas-phase reactivity of microsolvated nucleophiles in order to investigate in detail how the α-effect is influenced by solvent. We compare the gas-phase reactivity of the microsolvated α-nucleophile HOO(-)(H2O) to that of microsolvated normal alkoxy nucleophiles, RO(-)(H2O), in reaction with CH3Cl using a flowing afterglow-selected ion flow tube instrument. The results reveal enhanced reactivity of HOO(-)(H2O) and clearly demonstrate the presence of an α-effect for the microsolvated α-nucleophile. The association of the nucleophile with a single water molecule results in a larger Brønsted βnuc value than is the case for the unsolvated nucleophiles. Accordingly, the reactions of the microsolvated nucleophiles proceed through later transition states in which bond formation has progressed further. Calculations show a significant difference in solvent interaction for HOO(-) relative to the normal nucleophiles at the transition states, indicating that differential solvation may well contribute to the α-effect. The reactions of the microsolvated anions with CH3Cl can lead to formation of either the bare Cl(-) anion or the Cl(-)(H2O) cluster. The product distributions show preferential formation of the Cl(-) anion even though the formation of Cl(-)(H2O) would be favored thermodynamically. Although the structure of the HOO(-)(H2O) cluster resembles HO(-)(HOOH), we demonstrate that HOO(-) is the active nucleophile when the cluster reacts.