Temperature Dependence of the Rate Constants and Branching Ratios for the Reactions of Cl-(D2O)1-3 with CH3Br and Thermal Dissociation Rates for Cl-(CH3Br)

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.

Gas-Phase Reactivity of Ozone with Lanthanide Ions (Sm+, Nd+) and Their Higher Oxides

The kinetics of SmO_n⁺ (n = 0-2) and NdO_n⁺ (n = 0-2) with O₃ are measured using a selected-ion flow tube. Reaction of Nd⁺ to yield NdO⁺ + O₂ occurs rapidly, with a rate constant near the capture-controlled limit of ∼8 × 10^-10 cm³ s^-1. NdO⁺ reacts at ∼40% of the capture limit to yield NdO₂⁺ with little temperature dependence from 200 to 400 K. NdO₂⁺ likely reacts very slowly (k ∼ 10^-13 cm³ s^-1) to yield NdO⁺ + 2O₂, does not react to yield NdO₃⁺, and associates slowly (k ∼ 10^-12 cm³ s^-1) to yield NdO₂⁺(O₃)_1-3. Reaction of Sm⁺ also yields SmO⁺ at near the capture limit at all temperatures, but a significant fraction (∼50%) of the SmO⁺ is produced in excited states that are long-lived compared to the millisecond time scale of the experiment. These states are evidently resistant to both radiative and collisional relaxation. The excited-state production is likely due to a spin-conservation constraint on the reaction, despite the large spin-orbit coupling typical for lanthanide-containing species. Ground-state SmO⁺ reacts inefficiently (k = 2 × 10^-11 (T/300)^-2.5 cm³ s^-1) to yield SmO₂⁺ + O₂, while the excited-state SmO^+* reacts at the capture limit, with branching to yield Sm⁺ + 2O₂ (ΔH_r,0K = 148.7 ± 0.4 kJ mol^-1 for ground-state SmO⁺) approximately 60% of the time, the remainder forming SmO₂⁺, which further reacts with O₃ to yield SmO⁺ at about 1% of the collisional value.

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.

Temperature and energy dependences of ion-molecule reactions: Studies inspired by Diethard Böhme

Diethard Böhme has had a long career covering many topics in ion-molecule reactivity. In this review, we describe the work done at the Air Force Research Laboratory (and its variously named preceding organizations) that was inspired by his studies. These fall into two main areas: nucleophilic displacement (S_N 2) and metal cation chemistry. In S_N 2 chemistry, we revisited many of the reactions Diethard pioneered and studied them in more detail. We found nonstatistical behavior, both competition and noncompetition between multiple channels. New channels were found as hydration occurred, with more solution-like behavior occurring as only a few ligands were added. Temperature-dependent studies revealed details that were not observable at room temperature. These and other highlights will be discussed. In metal cation reactions, Diethard's use of an inductively coupled ion source allowed him to systematically study the periodic table of elements with a number of simple neutrals. We have taken the most interesting of these and studied them in greater detail. In doing so, we were able to identify curve crossing rates, in a few instances information about product states, and the importance of multiple entrance channels.

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.

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

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

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.

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.

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.

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.

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.

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.

Microsolvated F(-)(H2O) + CH3I S(N)2 Reaction Dynamics. Insight into the Suppressed Formation of Solvated Products

Microsolvation offers a bottom-up approach to investigate details of how solute-solvent interactions affect chemical reaction dynamics. The dynamics of the microsolvated S(N)2 reaction F(-)(H2O) + CH3I are uncovered in detail by using direct chemical dynamics simulations. Direct rebound and stripping and indirect atomic-level mechanisms are observed. The indirect events comprise ∼70% of the solvated reaction and occur predominantly via a hydrogen-bonded F(-)(H2O)···HCH2I prereaction complex. The reaction dynamics show propensity for the direct three-body dissociation channel F(-)(H2O) + CH3I → CH3F + I(-) + H2O after passing the reaction's dynamical bottleneck. The water molecule leaves the reactive system before traversing the postreaction region of the PES, where water transfer toward the product species occurs. This provides an insight into the very interesting finding of strongly suppressed formation of energetically favored solvated products for almost all SN2 reactions under microsolvation.

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.

Vibrationally induced charge transfer in a bimolecular model complex in vacuo

We report vibrationally induced charge transfer from nitromethane anion to methyliodide in a molecular complex. Excitation of a CH stretching vibrational transition in either of the molecular constituents results in dissociative electron transfer to the CH3I molecule, resulting in I(-) product anions. Solvation of the pre-reactive complex with more than two Ar atoms leads to complete quenching of the reaction and can be used to estimate the barrier for this reaction. We discuss the results in the framework of electronic structure calculations and compare the intra-complex electron transfer with vibrationally mediated electron emission in bare nitromethane anion.

Quantum chemical calculations of the Cl- + CH3I --> CH3Cl + I- potential energy surface

Electronic structure theory calculations, using MP2 theory and the DFT functionals OPBE, OLYP, HCTH407, BhandH, and B97-1, were performed to characterize the structures, vibrational frequencies, and energies for stationary points on the Cl(-) + CH(3)I --> ClCH(3) + I(-) potential energy surface. The aug-cc-pVDZ and aug-cc-pVTZ basis sets, with an effective core potential (ECP) for iodine, were employed. Single-point CCSD(T) calculations were performed to obtain the complete basis set (CBS) limit for the reaction energies. DFT was found to give significantly longer halide ion/carbon atom bond lengths for the ion-dipole complexes and central barrier transition state than MP2. BhandH, with either the aug-cc-pVDZ or aug-cc-pVTZ basis sets, gives good agreement with the experimental structures for both CH(3)I and CH(3)Cl. The frequencies of CH(3)I and CH(3)Cl, obtained with the different levels of theory and basis sets, are in excellent agreement with experiment. The major difference between the MP2 and DFT frequencies is for the imaginary frequency of the central barrier. Using the aug-cc-pVTZ basis the MP2 value for this frequency ranges from 1.26 to 1.59 times larger than those for the DFT functionals. Thus, the MP2 and DFT theories have different PES shapes in the vicinity of the [Cl--CH(3)--I](-) central barrier. The CCSD(T)/CBS energies are in good agreement with experiments for the complexation energies and reaction exothermicity, with a small 1 kcal/mol difference for the latter. The CCSD(T)/CBS central barrier height is lower than values deduced by using statistical theoretical models to fit the Cl(-) + CH(3)I --> ClCH(3) + I(-) experimental rate constant, which is consistent with the expected nonstatistical dynamics for the reaction. The BhandH energies are in overall best agreement with the CCSD(T) values, with a largest difference of only 0.7 kcal/mol.

Interpretation of the Gas-Phase Solvent Deuterium Kinetic Isotope Effects in the SN2 Reaction Mechanism: Comparison of Theoretical and Experimental Results in the Reaction of Microsolvated Fluoride Ions with Methyl Halides

We carried out a comprehensive ab initio calculation and transition-state theory analysis of the solvent and secondary deuterium kinetic isotope effects in the SN2 reactions of microsolvated fluoride ions with methyl halides. Water, methanol, and hydrogen fluoride were used as solvents, and the results are compared with recent experiments. Kinetic isotope effects were dissected into contributions from translations, rotations, and different vibration modes, and the validity of such analysis is also discussed. Excellent agreement was found for some reactions, whereas the agreement was poor for other reactions. We showed that the deviation between theory and experiments is related to the reaction kinetics; a faster reaction produced a kinetic isotope effect that was systematically larger (less inverse) than the calculated value. In addition, we also found that the magnitude of the deviation is proportional to the reaction efficiency. We rationalize the disagreement as a failure of the transition-state theory to model barrierless reactions, and we propose a very simple scheme to interpret these findings and predict the deviation between experimental and theoretical values in those reactions.

Gas Phase SN2 Reactions of Halide Ions with Trifluoromethyl Halides: Front- and Back-Side Attack vs. Complex Formation

Density functional theory computations and pulsed-ionization high-pressure mass spectrometry experiments have been used to explore the potential energy surfaces for gas-phase S(N)2 reactions between halide ions and trifluoromethyl halides, X(-) + CF(3)Y --> Y(-) + CF(3)X. Structures of neutrals, ion-molecule complexes, and transition states show the possibility of two mechanisms: back- and front-side attack. From pulsed-ionization high-pressure mass spectrometry, enthalpy and entropy changes for the equilibrium clustering reactions for the formation of Cl(-)(BrCF(3)) (-16.5 +/- 0.2 kcal mol(-1) and -24.5 +/- 1 cal mol(-1) K(-1)), Cl(-)(ICF(3)) (-23.6 +/- 0.2 kcal mol(-1)), and Br(-)(BrCF(3)) (-13.9 +/- 0.2 kcal mol(-1) and -22.2 +/- 1 cal mol(-1) K(-1)) have been determined. These are in good to excellent agreement with computations at the B3LYP/6-311+G(3df)//B3LYP/6-311+G(d) level of theory. It is shown that complex formation takes place by a front-side attack complex, while the lowest energy S(N)2 reaction proceeds through a back-side attack transition state. This latter mechanism involves a potential energy profile which closely resembles a condensed phase S(N)2 reaction energy profile. It is also shown that the Cl(-) + CF(3)Br --> Br(-) + CF(3)Cl S(N)2 reaction can be interpreted using Marcus theory, in which case the reaction is described as being initiated by electron transfer. A potential energy surface at the B3LYP/6-311+G(d) level of theory confirms that the F(-) + CF(3)Br --> Br(-) + CF(4) S(N)2 reaction proceeds through a Walden inversion transition state.

Solvent Effects on the SN2 Reaction: Application of the Density Functional Theory-Based Effective Fragment Potential Method

The performance of the density functional theory (DFT)-based effective fragment potential (EFP) method is assessed using the S(N)2 reaction: Cl- + nH2O + CH3Br = CH3Cl + Br- + nH2O. The effect of the systematic addition of water molecules on the structures and relative energies of all species involved in the reaction has been studied. The EFP1 method is compared with second-order perturbation theory (MP2) and DFT results for n = 1, 2, and 3, and EFP1 results are also presented for four water molecules. The incremental hydration effects on the barrier height are the same for all methods. However, only full MP2 or MP2 with EFP1 solvent molecules are able to provide an accurate treatment of the transition state (TS) and hence the central barriers. Full DFT and DFT with EFP1 solvent molecules both predict central barriers that are too small. The results illustrate that the EFP1-based DFT method gives reliable results when combined with an accurate quantum mechanical (QM) method, so it may be used as an efficient alternative to fully QM methods in the treatment of larger microsolvated systems.

Trajectory studies of S(N)2 nucleophilic substitution. 8. Central barrier dynamics for gas phase Cl(-) + CH(3)Cl

Quasiclassical direct dynamics trajectories, calculated at the MP2/6-31G level of theory, are used to study the central barrier dynamics for the C1(-) + CH(3)Cl S(N)2 reaction. Extensive recrossings of the central barrier are observed in the trajectories. The dynamics of the Cl(-)-CH(3)Cl complex is non-RRKM and transition state theory (TST) is predicted to be an inaccurate model for calculating the Cl(-) + CH(3)Cl S(N)2 rate constant. Direct dynamics trajectories also show that Cl(-) + CH(3)Cl trajectories, which collide backside along the S(N)2 reaction path, do not form the Cl(-)-CH(3)Cl complex. This arises from weak coupling between the Cl(-)-CH(3)Cl intermolecular and CH(3)Cl intramolecular modes. The trajectory results are very similar to those of a previous trajectory study, based on a HF/6-31G* analytic potential energy function, which gives a less accurate representation of the central barrier region of the Cl(-) + CH(3)Cl reaction than does the MP2/6-31G* level of theory used here. Experiments are suggested for investigating the non-RRKM and non-TST dynamics predicted by the trajectories.