Jiaxu Zhang

Imaging nucleophilic substitution dynamics.

Anion-molecule nucleophilic substitution (SN2) reactions are known for their rich reaction dynamics, caused by a complex potential energy surface with a submerged barrier and by weak coupling of the relevant rotational-vibrational quantum states. The dynamics of the SN2 reaction of Cl– + CH3I were uncovered in detail by using crossed molecular beam imaging. As a function of the collision energy, the transition from a complex-mediated reaction mechanism to direct backward scattering of the I– product was observed experimentally. Chemical dynamics calculations were performed that explain the observed energy transfer and reveal an indirect roundabout reaction mechanism involving CH3 rotation.

Chemical Dynamics Simulations of X– + CH3Y → XCH3 + Y– Gas-Phase SN2 Nucleophilic Substitution Reactions. Nonstatistical Dynamics and Nontraditional Reaction Mechanisms

Extensive classical chemical dynamics simulations of gas-phase X(-) + CH(3)Y → XCH(3) + Y(-) S(N)2 nucleophilic substitution reactions are reviewed and discussed and compared with experimental measurements and predictions of theoretical models. The primary emphasis is on reactions for which X and Y are halogen atoms. Both reactions with the traditional potential energy surface (PES), which include pre- and postreaction potential energy minima and a central barrier, and reactions with nontraditional PESs are considered. These S(N)2 reactions exhibit important nonstatistical atomic-level dynamics. The X(-) + CH(3)Y → X(-)---CH(3)Y association rate constant is less than the capture model as a result of inefficient energy transfer from X(-)+ CH(3)Y relative translation to CH(3)Y rotation and vibration. There is weak coupling between the low-frequency intermolecular modes of the X(-)---CH(3)Y complex and higher frequency CH(3)Y intramolecular modes, resulting in non-RRKM kinetics for X(-)---CH(3)Y unimolecular decomposition. Recrossings of the [X--CH(3)--Y](-) central barrier is important. As a result of the above dynamics, the relative translational energy and temperature dependencies of the S(N)2 rate constants are not accurately given by statistical theory. The nonstatistical dynamics results in nonstatistical partitioning of the available energy to XCH(3) +Y(-) reaction products. Besides the indirect, complex forming atomic-level mechanism for the S(N)2 reaction, direct mechanisms promoted by X(-) + CH(3)Y relative translational or CH(3)Y vibrational excitation are possible, e.g., the roundabout mechanism.

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.

Indirect Dynamics in a Highly Exoergic Substitution Reaction

The highly exoergic nucleophilic substitution reaction F(-) + CH3I shows reaction dynamics strikingly different from that of substitution reactions of larger halogen anions. Over a wide range of collision energies, a large fraction of indirect scattering via a long-lived hydrogen-bonded complex is found both in crossed-beam imaging experiments and in direct chemical dynamics simulations. Our measured differential scattering cross sections show large-angle scattering and low product velocities for all collision energies, resulting from efficient transfer of the collision energy to internal energy of the CH3F reaction product. Both findings are in strong contrast to the previously studied substitution reaction of Cl(-) + CH3I [Science 2008, 319, 183-186] at all but the lowest collision energies, a discrepancy that was not captured in a subsequent study at only a low collision energy [J. Phys. Chem. Lett. 2010, 1, 2747-2752]. Our direct chemical dynamics simulations at the DFT/B97-1 level of theory show that the reaction is dominated by three atomic-level mechanisms, an indirect reaction proceeding via an F(-)-HCH2I hydrogen-bonded complex, a direct rebound, and a direct stripping reaction. The indirect mechanism is found to contribute about one-half of the overall substitution reaction rate at both low and high collision energies. This large fraction of indirect scattering at high collision energy is particularly surprising, because the barrier for the F(-)-HCH2I complex to form products is only 0.10 eV. Overall, experiment and simulation agree very favorably in both the scattering angle and the product internal energy distributions.

The need for enzymatic steering in abietic acid biosynthesis: gas-phase chemical dynamics simulations of carbocation rearrangements on a bifurcating potential energy surface.

Abietic acid, a constituent of pine resin, is naturally derived from abietadiene --a process that requires four enzymes: one (abietadiene synthase) for conversion of the acyclic, achiral geranylgeranyl diphosphate to the polycyclic, chiral abietadiene (a complex process involving the copalyl diphosphate intermediate) and then three to oxidize a single methyl group of abietadiene to the corresponding carboxylic acid. In previous work (Nature Chem.2009, 1, 384), electronic structure calculations on carbocation rearrangements leading to abietadienyl cation revealed an interesting potential energy surface with a bifurcating reaction pathway (two transition-state structures connected directly with no intervening minimum), which links two products--one natural and one not yet isolated from Nature. Herein we describe direct dynamics simulations of the key step in the formation of abietadiene (in the gas phase and in the absence of the enzyme). The simulations reveal that abietadiene synthase must intervene in order to produce abietadiene selectively, in essence steering this reaction to avoid the generation of byproducts with different molecular architectures.

Simulation studies of the Cl− + CH3I SN2 nucleophilic substitution reaction: Comparison with ion imaging experiments

In the previous work of Mikosch et al. [Science 319, 183 (2008)], ion imaging experiments were used to study the Cl(-) + CH3I → ClCH3 + I(-) reaction at collision energies E(rel) of 0.39, 0.76, 1.07, and 1.9 eV. For the work reported here MP2(fc)/ECP/d direct dynamics simulations were performed to obtain an atomistic understanding of the experiments. There is good agreement with the experimental product energy and scattering angle distributions for the highest three E(rel), and at these energies 80% or more of the reaction is direct, primarily occurring by a rebound mechanism with backward scattering. At 0.76 eV there is a small indirect component, with isotropic scattering, involving formation of the pre- and post-reaction complexes. All of the reaction is direct at 1.07 eV. Increasing E(rel) to 1.9 eV opens up a new indirect pathway, the roundabout mechanism. The product energy is primarily partitioned into relative translation for the direct reactions, but to CH3Cl internal energy for the indirect reactions. The roundabout mechanism transfers substantial energy to CH3Cl rotation. At E(rel) = 0.39 eV both the experimental product energy partitioning and scattering are statistical, suggesting the reaction is primarily indirect with formation of the pre- and post-reaction complexes. However, neither MP2 nor BhandH/ECP/d simulations agree with experiment and, instead, give reaction dominated by direct processes as found for the higher collision energies. Decreasing the simulation E(rel) to 0.20 eV results in product energy partitioning and scattering which agree with the 0.39 eV experiment. The sharp transition from a dominant direct to indirect reaction as E(rel) is lowered from 0.39 to 0.20 eV is striking. The lack of agreement between the simulations and experiment for E(rel) = 0.39 eV may result from a distribution of collision energies in the experiment and/or a shortcoming in both the MP2 and BhandH simulations. Increasing the reactant rotational temperature from 75 to 300 K for the 1.9 eV collisions, results in more rotational energy in the CH3Cl product and a larger fraction of roundabout trajectories. Even though a ClCH3-I(-) post-reaction complex is not formed and the mechanistic dynamics are not statistical, the roundabout mechanism gives product energy partitioning in approximate agreement with phase space theory.

Electronic Structure Theory Study of the F− + CH3I → FCH3 + I− Potential Energy Surface†

MP2 and DFT electronic structure theories, with the functionals OPBE, OLYP, HCTH407, BhandH, and B97-1 for the latter, were used to investigate stationary point properties on the F(-) + CH(3)I → FCH(3) + I(-) potential energy surface (PES). The aug-cc-pVDZ and aug-cc-pVTZ basis sets for C, H, and F, with Wadt and Hays 3s3p valence basis and an effective core potential (ECP) for iodine, were employed for both MP2 and DFT. Single-point CCSD(T) calculations were also performed to obtain the complete basis set (CBS) limit for the stationary point energies. The CCSD(T)/CBS reaction exothermicity is only 5.0 kJ/mol different than the experimental value. MP2 and DFT do not predict the same stationary points on the PES. MP2 predicts the C(3v) F(-)-CH(3)I and FCH(3)-I(-) ion-dipole complexes and traditional [F-CH(3)-I](-) central barrier as stationary points, as well as a C(s) hydrogen-bonded F(-)-HCH(2)I complex and a [F-HCH(2)-I](-) transition state connecting this latter complex to the F(-)-CH(3)I complex. A CCSD(T)/CBS relaxed potential energy curve, calculated for the MP2 structures, shows that going from the F(-)-CH(3)I complex to the [F-CH(3)-I](-) TS is a barrierless process, indicating these two structures are not stationary points. This is also suggested by the DFT calculations. The structures and frequencies for CH(3)I and CH(3)Cl given by MP2 and DFT are in overall good agreement with experiment. The calculations reported here indicate that the DFT/B97-1 functional gives the overall best agreement with the CCSD(T) energies, with a largest difference of only 7.5 kJ/mol for the FCH(3)-I(-) complex.

Dynamics of the F(-) + CH3I → HF + CH2I(-) Proton Transfer Reaction.

Direct chemical dynamics simulations, at collision energies Erel of 0.32 and 1.53 eV, were performed to obtain an atomistic understanding of the F(-) + CH3I reaction dynamics. There is only the F(-) + CH3I → CH3F + I(-) bimolecular nucleophilic substitution SN2 product channel at 0.32 eV. Increasing Erel to 1.53 eV opens the endothermic F(-) + CH3I → HF + CH2I(-) proton transfer reaction, which is less competitive than the SN2 reaction. The simulations reveal proton transfer occurs by two direct atomic-level mechanisms, rebound and stripping, and indirect mechanisms, involving formation of the F(-)···HCH2I complex and the roundabout. For the indirect trajectories all of the CH2I(-) is formed with zero-point energy (ZPE), while for the direct trajectories 50% form CH2I(-) without ZPE. Without a ZPE constraint for CH2I(-), the reaction cross sections for the rebound, stripping, and indirect mechanisms are 0.2 ± 0.1, 1.2 ± 0.4, and 0.7 ± 0.2 Å(2), respectively. Discarding trajectories that do not form CH2I(-) with ZPE reduces the rebound and stripping cross sections to 0.1 ± 0.1 and 0.7 ± 0.5 Å(2). The HF product is formed rotationally and vibrationally unexcited. The average value of J is 2.6 and with histogram binning n = 0. CH2I(-) is formed rotationally excited. The partitioning between CH2I(-) vibration and HF + CH2I(-) relative translation energy depends on the treatment of CH2I(-) ZPE. Without a CH2I(-) ZPE constraint the energy partitioning is primarily to relative translation with little CH2I(-) vibration. With a ZPE constraint, energy partitioning to CH2I(-) rotation, CH2I(-) vibration, and relative translation are statistically the same. The overall F(-) + CH3I rate constant at Erel of both 0.32 and 1.53 eV is in good agreement with experiment and negligibly affected by the treatment of CH2I(-) ZPE, since the SN2 reaction is the major contributor to the total reaction rate constant. The potential energy surface and reaction dynamics for F(-) + CH3I proton transfer are compared with those reported previously (J. Phys. Chem. A 2013, 117, 7162-7178) for the isoelectronic OH(-) + CH3I reaction.

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.

Microsolvated F-(H2O) + CH3I SN2 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 reactions 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.