|Home | About | Journals | Submit | Contact Us | Français|
Reactions occurring at a carbon atom through the Walden inversion mechanism are one of the most important and useful classes of reactions in chemistry. Here we report an accurate theoretical study of the simplest reaction of that type: the H+CH4 substitution reaction and its isotope analogues. It is found that the reaction threshold versus collision energy is considerably higher than the barrier height. The reaction exhibits a strong normal secondary isotope effect on the cross-sections measured above the reaction threshold, and a small but reverse secondary kinetic isotope effect at room temperature. Detailed analysis reveals that the reaction proceeds along a path with a higher barrier height instead of the minimum-energy path because the umbrella angle of the non-reacting methyl group cannot change synchronously with the other reaction coordinates during the reaction due to insufficient energy transfer from the translational motion to the umbrella mode.
Reactions occurring at a carbon atom in a tetrahedral environment, proceeding through the back-side attack Walden inversion mechanism leading to inversion of the chirality of a molecule, are one of the most important and useful classes of reactions in chemistry. The bimolecular nucleophilic substitution (SN2) reaction is the most common reaction of that type, in which a nucleophile (often negatively charged) approaches a saturated carbon from one side, displaces a leaving group on the opposite side of the carbon atom, resulting in inversion of the carbon centre. Strong solvent effects of these reactions in solution have prompted investigations of the gas-phase SN2 reaction to probe the intrinsic reaction mechanisms without solvent1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27. Numerous experimental and theoretical studies have revealed that an inverse secondary kinetic isotope effect (KIE), that is, kH/kD<1 (where k denotes the thermal rate constant), is characteristic for a thermal SN2 reaction when the isotopically substituted atom is not directly involved in a reaction1,2,3,4,5,6,7,8,9,10. In strong contrast, by using a guided ion beam technique, Ervin and co-workers discovered that the symmetric 37Cl−+CH335Cl→35Cl−+CH337Cl reaction exhibits a large normal secondary isotope effect on the cross-section, in addition to a reaction threshold versus collision energy substantially higher than the calculated barrier height11,14. They also observed that the reaction threshold for the endothermic Cl−+ CH3F→CH3Cl+F− reaction is considerably higher than the reaction endothermicity12. Hase and co-workers performed extensive direct molecular dynamics studies on many exothermic and a few centre barrier SN2 reactions, and uncovered the non-statistical nature and inefficiency of energy transfer between intermolecular and intramolecular modes in these reactions13,14,15,16,17,18. Furthermore, they found that the reaction thresholds for SN2 reactions, which are larger than their corresponding barrier heights, result from direct reactions15,17. Many reduced dimensionality quantum scattering studies were also carried out to investigate dynamics in SN2 reactions19,20,21,22. Recently, the molecular beam experiments in the Wester's group, in combination with theory, have revealed unprecedented dynamical details on some exothermic SN2 reactions23,24,25,26,27. However, despite so many experimental and theoretical investigations, the origin of the intriguing difference between the isotope effects on the kinetic and cross-sections as well as the high-energy threshold for centre barrier and endothermic SN2 reactions remain unclear.
The H+CH4 reaction is the simplest and most prototypical reaction occurring at a carbon atom in a tetrahedral environment. The abstraction H+CH4→H2+CH3 reaction has been the subject of numerous experimental investigations, and has also become a benchmark system for theoretical studies of polyatomic reactions28,29,30,31,32,33,34,35,36,37,38,39. In addition to the abstraction reaction, there exists a substitution reaction, H′+CH4→H+CH′H3, with a D3h transition state and a static barrier height of 1.6eV. It is the simplest reaction proceeding through the back-side attack Walden inversion mechanism, very similar to the gas-phase SN2 reactions with central barriers, except that there exist pre- and post-reaction wells in SN2 reactions stemming from strong ion–dipole interaction between reagents/products. The experimental study of the substitution reaction has a long history. Rowland and co-workers observed a threshold energy of ~1.5eV for the T+CD4→CD3T+H reaction40, and obtained a relative yield of 7.2 for the T+CH4 and T+CD4 reaction using tritium atoms with a translational energy of 2.8eV (ref. 41). Bersohn and co-worker measured the substitution cross-sections for the H+CD4→CHD3+D and H+CH3D→CH4+D reactions at a collision energy near 2eV to directly study the secondary isotope effects, and found that the cross-section was reduced by about a factor of 2 from 0.040±0.015Å2 for H+CH3D to 0.021±0.005Å2 for H+CD4 (per D atom), indicating that the secondary isotope effect is strong in the reaction, although not as prominent as observed by Rowland and co-workers42. On the basis of measured cross-sections and the velocity distribution measured for the reacting H atom and displaced D atom, they concluded that the substitution reaction takes place by an inversion mechanism.
Despite these intriguing experimental discoveries, to the best of our knowledge there is no theoretical study of the substitution reaction on a high-quality potential energy surface (PES), except Bunker and co-workers performed trajectory calculations on a model analytic PES43. Quantitative theory for such a reaction faces two difficult tasks: the construction of an accurate global PES involving six atoms and the performance of reactive scattering calculations, preferably quantum mechanical, on an accurate PES. The past decade has witnessed significant progress on accurate PES construction and quantum reactive scattering study of polyatomic reactions38,39,44,45,46,47,48,49,50. These advances have been combined with the rise in computational power to make accurate dynamics study practical for such a reaction process.
Here we report an accurate theoretical study of the H+CH4 substitution reaction. Our calculation shows the reaction has (a) a threshold energy considerably higher than the static barrier height, (b) a large normal secondary isotope effect on cross-section and (c) a small but reverse secondary KIE, in exactly the same way as have been observed in many SN2 reactions. Therefore, the study not only provides unprecedented dynamical details for this simplest Walden inversion reaction, but also sheds light on the dynamics of gas-phase SN2 reactions.
The total reaction probabilities for the substitution reactions
are presented as a function of collision energy for their corresponding ground rovibrational initial states in Fig. 1. (Note the probability for reaction (1) is for one H atom substitution). The reaction probabilities increase smoothly and quickly with the increasing of the translational energy, but the overall values for all the three reactions are small. The probability for reaction (1), which is the largest among these three reactions, only reaches 0.0124 at E=2.5eV. With the substituted atom changed from H to D, the substitution probability for reaction (2) drops substantially, by a factor of 2.6 at E=2.5eV. This large primary isotope effect, as also found in the H′+H2O→H′OH+H and H′+HOD→H′OH+D reactions51 is apparently due to the fact that the substituted H atom is faster than D atom in responding to the attack of the incoming H atom.
In contrast to the H′+H2O→H′OH+H and H′+HOD→H′OD+H reactions51, the reaction probability for reactions (1) and (3), with the same newly formed bond and cleaved bond, decreases substantially as the non-reacting group changes from CH3 to CD3, manifesting strong secondary isotope effects. As these two reactions process on the same PES with their variational bottlenecks located at the saddle point and their zero point energies (ZPEs) corrected barrier heights differ <0.01eV (1.591 versus 1.585 for reactions (1) and (3)), the substantial difference between two reactions apparently cannot be explained by the barrier height or the topography of the PES.
To explore the dynamics difference between these two isotopic reactions, we calculated the average value of the umbrella angle χ for methyl from the scattering wavefunction at the translational energy E=1.8eV. Figure 2a shows the average value of the umbrella angle, <χ>, as a function of R (the distance between incoming H and the centre of mass of CHD3) and r (the bond length of breaking CH bond) for the H+CHD3 substitution reaction, together with the potential energy contour obtained by minimizing the other degrees of freedom. The static saddle point locates at R=2.7Bohr and r=2.5Bohr with the corresponding umbrella angle χ=90°. However, the calculated <χ> equals to 108.5° at the static saddle point, close to the corresponding vibrationally averaged value for CHD3 reagent of 109.36°. With the increase of r (elongation of CH bond) or as the wavefunction proceeding to the product region from the static saddle point, <χ> decreases in an accelerated way, eventually passes through the 90° line at r~3.1Bohr, which is considerably larger than the value at the saddle point. This means that the reaction does not proceed along the minimum-energy path on which the umbrella angle <χ> changes synchronously during the reaction with the incoming of H atom and elongation of the breaking CH bond to reach a value close to 90° at the static saddle point. Instead, it essentially does not vary while the incoming H atom approaches to the C atom and the breaking bond CH bond is elongated to their corresponding desired values for reaction. Only after that, the umbrella angle χ starts to decrease, initially slowly, then increasingly fast. Therefore, in reality, the reaction proceeds not on the potential shown in Fig. 2a, but on a potential shown in Fig. 2b, which is obtained by taking the umbrella angle at the calculated average value for every combination of R and r. In Fig. 2b, the ‘saddle point' moves to R~2.85, r~3.1Bohr. This is the very dynamical saddle point for the reaction with a corresponding barrier height of ~1.93eV, higher than the original barrier height by 0.32eV, also locating much later than the original barrier in term of C–H bond length.
The above finding clearly reveals that the umbrella motion of the non-reacting CD3 group is slow in responding to the attack of the incoming H atom during the reaction. Because the umbrella motion depends on the moment of inertia (masses of atom in the methyl group), it is expected that isotope replacement of D atom to H atom can change the umbrella motion substantially. As seen from Fig. 2c for the H+CH4 reaction, the umbrella angle starts to decrease even before the incoming H atom reaches the static saddle point, reaches 90° at r~2.9Bohr in an accelerated way, and then begins to slow down with the further increase of r. As a result, the reaction has a dynamical saddle point at r~2.8Bohr with a corresponding dynamical barrier height of 1.73eV, lower than that for H+CHD3 by 0.2eV. Therefore, it is clear that the strong secondary isotope effects in the reaction system in term of the huge differences on the reaction thresholds versus collision energy and the magnitudes of reaction probabilities for these two isotopic reactions shown in Fig. 1 originate from the important effect of the umbrella motion of the methyl group on the reactions.
Figure 3a compares the reaction probabilities for a number of initial vibrational states for reaction (3) as a function of collision energy. Both the CH stretch excitation and CD3 umbrella excitation substantially enhance the reactivity and reduce the reaction threshold versus collision energy (see Supplementary Fig. 2 for more initial states and Supplementary Fig. 3 for reaction (1)). The reaction probability for the first umbrella excited (0,1) state, with an excitation energy of 0.126eV, is higher than that for the ground state by a factor of 67, 8.0 and 3.8, respectively, at E=1.5, 2.0 and 2.5eV. The reactivity for the (0,2) state, with an excitation energy of 0.252eV, is even higher than that for the first CH stretch excited state with an excitation energy of 0.369eV, in particular in low energy region, indicating that the umbrella excitation is more efficient on promoting the reaction than CH stretch excitation. This can be seen more clearly from Fig. 3b, which shows these reaction probabilities as a function of total energy measured from the ground vibrational of CHD3 (see Supplementary Fig. 4 for more initial states and Supplementary Fig. 5 for reaction (1)). Both the CH stretch excitation and CD3 umbrella excitation are more efficient than the translational energy on promoting the reaction, in particular for the umbrella motion the efficacy increases with initial excitation up to the 5th/6th overtone state (Supplementary Fig. 4). For reaction (1), the situation is similar as for reaction (3), with both the CH stretch excitation and CH3 umbrella excitation more efficient than the translational energy on promoting the reaction (Supplementary Fig. 5).
Because reaction (3) possess a late barrier, in particular a very late dynamical barrier, in term of the CH stretch as shown in Fig. 2a,b, it is expected that the CH stretch excitation enhances the reactivity substantially, according to Polanyi's rules, which state that vibrational energy is more efficient than translational energy in promoting a late barrier reaction52. To understand the pronounced effect of the umbrella excitation on the reactivity, we present in Fig. 2d the average umbrella angle for methyl group as a function of R and r for the H+CHD3(0,1) reaction at the translational energy E=1.8eV. Around the static saddle point, the average umbrella angle for the (0,1) state is very close to ground state shown in Fig. 2a. However, with the further elongation of r, the umbrella angle for the excited state decreases much faster than the ground state, resulting in a dynamical barrier height of ~1.80eV at (R=2.9, r=2.9). This dynamical barrier is lower than the ground state by ~0.19eV, larger than the (0,1) excitation energy of 0.126eV. Clearly, more vibrational energy in the umbrella mode due to initial excitation makes the umbrella angle faster to respond, giving rise to a reduced dynamical barrier height and enhancing the reactivity substantially. Therefore, the large difference on reactivity between the H+CH4/CHD3 reactions shown in Fig. 1 should also relate to the difference in ZPE, in addition to the difference in the moment of inertia discussed above. This, on the other hand, also indicates that the energy transfer from translational motion to umbrella mode is inefficient because an efficient energy transfer to the umbrella mode would enhance the efficacy of translational energy on prompting reaction.
From Fig. 3a, one can see that for the H+CHD3 reaction, the reactivity at a given collision energy above threshold, in particular at high energy, is mainly determined by that for the ground vibrational state, because even the first umbrella excited state with an excitation energy of 0.126eV has a small population at a moderate temperature (Say, at T=300K the population is ~0.75%). This means that measured at a given collision energy, the cross-sections for reaction (1) will be substantially larger than that for reaction (3), that is, one will observe a large normal secondary isotope effect on the integral cross-section.
However, the thermal rate constant for the reaction, as can be seen from Fig. 3b, is contributed by the vibrational excited states, in particular the highly umbrella excited states (Supplementary Fig. 4). With many vibrationally excited states contributed to the thermal rate constant, it is more efficient to calculate the cumulative reaction probabilities, NE(E), (the sum of the reaction probabilities for all the initial state with a fixed total energy) for the reactions (1) and (3) from which the thermal rate constants can be reliably estimated35,39 (see Supplementary Information for details). Figure 3c shows the NE(E) for reactions (1) and (3) as a function of total energy measured from their corresponding ground-state energies. Surprisingly, although the reaction probabilities for reaction (3) for the ground initial state is substantially smaller than reaction (1) as shown in Fig. 1, NE(E) for the reaction (3) is actually slightly larger than reaction (1) for the total energy up to 1.72eV. On the basis of NE(E), we obtained a thermal rate constant at 300K for reaction (1) of 3.1 × 10−36cm3s−1, for reaction (3) of 3.5 × 10−36cm3s−1, which results in kH/kD=0.89, a reverse KIE. It is worthwhile to point out here that the tiny rate constant is due to the 1.62eV static barrier height of the reaction, but the relative rate constant is still a good indicator of KIE.
The cumulative reaction probabilities shown in Fig. 3c for reactions (1) and (3) look intriguing because the reaction probability for the ground vibrational state for reaction (1) is much larger than that for reaction (3) as shown in Fig. 1. From the initial state selected reaction probability point of view, larger NE(E) for reaction (3) indicates contributions from the umbrella excited states for the reaction eventually beat those for reaction (1). While from the transition state point of view, this means in low energy the umbrella motion both for CH3 and CD3 can follow the motion of H atom to product region without much reflection of the wavefunction (or re-crossing of trajectory in classical terminology). Because the reaction (3) has a lower ZPE corrected barrier height and a relatively higher density of the umbrella states, it has a slightly larger NE(E). However, with the increase of energy, the umbrella motion of CD3 cannot act as fast as CH3 on responding to the H atom motion any more, causing sever reflection of the wavefunction that results in a smaller value of NE(E).
Finally, we depict integral cross-sections for T+CH4/CD4 substitution reactions in Fig. 4a and for the H+CD4/CH3D reaction in Fig. 4b, which have been measured experimentally. Same as the reaction probabilities for J=0 shown in Fig. 1, large but normal isotope effects can be observed for the reaction cross-sections. Agreement between theory and experiment on the integral cross-sections is both positive and negative. We obtained a relative yield of 6.5 for the T+CH4/CD4 substitution reactions for tritium atoms with a translational energy of 2.8eV, in comparison with the experimental value of 7.2 (ref. 41), despite the fact that there are some uncertainties in the comparison due to partial thermalization of the hot atoms before reaction in the experiment as discussed by Raff and co-workers53. However, for the H+CH3D and H+CD4 reactions, the relative theoretical yield at the collision energy of 2.2eV is 8.4, in comparison with experimental value of 2.0 (ref. 42). Even worse, the absolute theoretical integral cross-section for the H+CH3D→CH4+D reaction is only 0.001Å2 at that collision energy, smaller than the reported experimental value by a factor of 20 (ref. 42). Although the PES used in the present study is highly accurate, the theoretical cross-sections presented here are only for the ground rotational states, while the experimental results contain the contributions from all the rotational states populated in the experiment. Further studies should be performed to investigate the effects of initial rotational excitation of reagent on the reactions. In addition, the reduced dimensionality model used in this study by constraining the non-reacting methyl group to C3v symmetry32,34 may also introduce some errors, despite the fact that the model is at a quantitative level of accuracy for the abstraction reaction39.
We have carried out a detailed theoretical study of the dynamics of the H+CH4 substitution reaction and its isotope analogues. Our calculation reveals that the H+CH4 substitution reaction may manifest different isotope effects. In terms of the cross-section beyond the reaction threshold energy, it has a large normal secondary isotope effect because the measured cross-section at a moderate temperature is mainly contributed by the ground vibrational state. However, in terms of the thermal rate constant, it has a reverse secondary KIE due to contributions from the umbrella excited states.
All of the phenomena found from our calculation for the H+CH4 substitution reactions due to disparity between the translational motion and umbrella motion and inefficient energy transfer between them, including higher threshold energy, a reverse secondary KIE at room temperature and a large normal isotope effect for reaction cross-section, have been observed for gas-phase SN2 reactions with barriers. Extensive direct molecular dynamics simulations have revealed similar inefficient energy transfer as observed here and a non-statistical nature in SN2 reactions, despite the fact there exist pre- and post-barrier wells in gas-phase SN2 reactions13,14,15,16,17,18. Inefficient energy transfer between intermolecular and intramolecular modes render the umbrella mode inert to approach of the nucleophile until the nucleophile is very close to the C atom (the distance between the nucleophile and C atom close to that at the static saddle point), making the pre-barrier well not useful on activating the umbrella mode. Therefore, we anticipate that the mechanisms uncovered from this study may play important roles in these gas-phase SN2 reactions. More accurate dynamics studies should be carried out on accurate PESs for SN2 reactions with barriers to provide a more definitive answer.
The PES used in the calculation is the full dimensional global PES constructed in this group with neural network method based on ~46,000 UCCSD(T)-F12a/aug-cc-pVTZ ab initio energies. With a fitting error is 4meV, measured in term of root-mean-square error, the PES is highly accurate and capable of providing definite dynamical information for the reaction48.
Our quantum dynamics calculations employ the eight-dimensional model originally proposed by Palma and Clary32 by restricting the non-reacting CH3 group under C3V symmetry. Since the substitution reaction has a saddle point with D3h symmetry, the assumption should hold very well in this study and the model is expected to have a high level of accuracy as having been demonstrated in many studies on the abstraction reaction process38,39,49,50. Furthermore, because the length of the non-reacting CH bond essentially does not change during the reaction, it was fixed at the equilibrium value of 2.06Bohr for CH4, reducing the degrees of freedom included in our calculations to seven. Initial state selected wave packet method was employed to calculate total reaction probabilities for ground and some vibrationally excited initial states as a function of collision energy for some isotope combinations34. Transition state wave packet calculations were also carried out to obtain cumulative reaction probabilities from which the thermal rate constants can be evaluated35,39. Details of the calculations are provided in Supplementary Information.
The data that support the findings of this study are available from the corresponding author on request.
How to cite this article: Zhao, Z. et al. Dynamical barrier and isotope effects in the simplest substitution reaction via Walden inversion mechanism. Nat. Commun. 8, 14506 doi: 10.1038/ncomms14506 (2017).
Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Figures, Supplementary Methods and Supplementary References.
This work was supported by the National Natural Science Foundation of China (grant nos 21433009, 21688102 and 21403223), the Ministry of Science and Technology of China (2013CB834601) and the Chinese Academy of Sciences (XDB17010000).
The authors declare no competing financial interests.
Author contributions D.H.Z. and S.L. conceived and supervised the research; Z. Zhao, Z. Zhang and S.L. performed the research; Z. Zhao and D.H.Z. analysed the data; and D.H.Z. and S.L. wrote the manuscript.