Tim
Michaelsen
,
Björn
Bastian
,
Patrick
Strübin
,
Jennifer
Meyer
and
Roland
Wester
*
Institut für Ionenphysik und Angewandte Physik, Universität Innsbruck, Technikerstraße 25/3, 6020 Innsbruck, Austria. E-mail: roland.wester@uibk.ac.at; Tel: +43 512 507 52620
First published on 22nd April 2020
Gaining insight how specific rovibrational states influence reaction kinetics and dynamics is a fundamental goal of physical chemistry. Purely statistical approaches often fail to predict the influence of a specific state on the reaction outcome, evident in a great number of both experimental and theoretical studies. Most detailed insight in atomistic reaction mechanisms is achieved using accurate collision experiments and high level dynamics calculations. For ion–molecule reactions such experiments are scarce. Here we show the influence of symmetric CH-stretching vibration on the rate and dynamics of proton transfer in the reaction of F− + CH3I. We find a pronounced shift in the reaction dynamics for excited reactions from indirect to preferred direct dynamics at higher collision energy. Moreover, excited reactions occur at larger impact parameters. Finally, we compare vibrational excitation with collision energy and find that vibration is overall more efficient in promoting reactivity, which agrees with recent theoretical calculations.
In order to extend the Polanyi rules to poly-atomic systems and predict the influence of vibrations, the sudden vector projection model (SVP) has been developed by Guo and coworkers.24,25 It projects the different internal modes and translational coordinates onto the intrinsic reaction coordinate at the transition state and calculates the individual overlap. From these overlap integrals a qualitative ranking is obtained of how effective each mode is at promoting reactivity. The SVP method has been successfully applied to explain many different experimental observations.25
Here we present an experimental study about the influence of symmetric CH-stretching excitation on the proton transfer channel of F− reacting with CH3I. The dominant channel in this reaction is the bimolecular nucleophilic substitution (SN2), but three additional endothermic pathways exist (reaction enthalpies taken from ref. 26 and ref. 27).
F− + CH3I → I− + CH3F ΔH = −1.96 eV (SN2) | (1) |
→ CH2I− + HF ΔH = 0.69 eV (proton tr.) | (2) |
→ FI− + CH3 ΔH = 0.7 eV (halide abstr.) | (3) |
→ FHI− + CH2 ΔH = 0.9 eV (halide abstr.) | (4) |
We have recently shown that for the SN2 channel the symmetric stretching mode acts as a spectator, which is confirmed by SVP and QCT calculations.28 For the proton transfer reaction the situation is different and more similar to neutral halide reactions with CHD3,19 as we excite the bond that is broken during the reaction. This leads to interesting questions: Is vibration more efficient in promoting the reaction compared to collision energy and if so, does the enhancement depend on the collision energy? Is there a difference in the reaction dynamics for ground and excited state? Are stereodynamic effects important? To answer these questions we have obtained differential cross sections with and without vibrational excitation and compare our results to recently published QCT calculations by the Czakó group.26
We obtain the fraction of excited molecules experimentally using a combination of spatial map imaging (SMI),31 electron photo-detachment of the F− in the ion beam and a photo-dissociation/REMPI scheme for CH3I adapted from Hu et al.32 This requires knowledge of the overlap of the two reactant beams and the IR laser. Fig. 1a shows a sketch of this overlap inside the bottom two VMI electrodes (beams scaled up for better visibility). The size of the ion and neutral beam is measured using SMI, photo-detachement and photo-dissociation/REMPI. The neutral beam was considerably smaller compared to the ion beam in all measurements, which reduces the overlap problem to the IR laser beam and the neutral beam. The excitation can be directly mapped using REMPI of the stretch-excited state (v = 1) of the CH3 fragment. By moving the spherical lens that focuses the photo-dissociation and REMPI UV lasers vertically through the neutral beam, as shown by the black arrow in the sketch in Fig. 1a, we obtain the spatial profile shown in green without and in red with the IR laser present in Fig. 1b. To quantify the excitation we use the depletion of the ground state at the center of the neutral beam (marked by shaded area in Fig. 1b) and extrapolate this result to the entire beam. The excited fraction was determined regularly and the mean value is used to scale our results to a fully excited beam. It is important to note that the vibrational mode stays for relatively long times (> 4 μs) without intramolecular vibrational redistribution (IVR), as tested using delayed REMPI.
We record the velocity and arrival time for each product ion, which allows us to obtain branching ratios between competing product channels and differential cross sections for each product in a single measurement. To obtain sufficient statistics (roughly 1 million counts in the SN2 channel) at the three investigated collision energies a total measurement time of approx. 650 hours was necessary. Regular systematic checks of the collision energy (approx. every 6 h), the excited fraction and the beam overlaps (approx. every 30 h) were performed to ensure stable experimental conditions. The product velocity distributions are represented in the center of mass frame, where we choose the neutral beam direction as forward (see illustration at the top of Fig. 2). We project the 3D velocity distribution onto a 2D representation of the velocity along the collision axis vx and the radial velocity vr perpendicular to the collision axis.
The velocity distribution without vibrational excitation reveals two mechanisms leading to proton transfer (see Fig. 2a). A direct, forward scattered mechanism, often referred to as stripping and an indirect, isotropically scattered mechanism around zero product velocity. We know from previous studies,27 that the general trend with additional collision energy in this reaction leads to more forward, direct scattering. For the vibrationally excited reaction, the additional energy instead leads only to an indirect scattering signature (see Fig. 2b). The direct mechanism is suppressed. The indirect mechanism shows high internal excitation, in fact the entire additional energy seems to be efficiently transferred into internal excitation, which is evident in the Eint distribution shown in red in Fig. 2c. Taking a closer look at the angular distribution of the indirect mechanism for vibrationally excited reactions shown in red in Fig. 2d reveals forward–backward symmetric scattering. Such features are typically attributed to reactive collisions that occur at large impact parameters,33–36 where conservation of angular momentum constrains the leaving products into directions along the collision axis. In contrast, the indirect mechanism for the ground state reaction, marked in green in Fig. 2c and d, shows completely isotropic scattering (the maximum velocity chosen for this mechanism is marked in green in Fig. 2a). For isotropic scattering the mechanical rotation of the molecular reactant contributes significantly to the total angular momentum. Assuming the coupling of angular momentum does not change, we can attribute the shift from isotropic scattering to a forward–backward symmetric 1/cos(θ) distribution to an increase in initial orbital angular momentum and therefore larger reactive impact parameters.
At the intermediate collision energy of 1.2 eV, we can no longer clearly separate the two mechanisms in the ground state velocity distribution shown in Fig. 2e. The forward and isotropic parts overlap to a combined distribution. Interestingly, in the excited state image in Fig. 2f two contributions might be separable into a direct forward and an indirect, isotropic mechanism. This is also evident in the internal energy distribution (Fig. 2g), which is broader compared to the laser off distribution. The FWHM are 0.44 and 0.81 eV respectively and the difference is therefore very close the vibrational energy (0.37 eV). Instead of the forward–backward symmetric angular distribution at 0.7 eV, the excited reactions lead to more forward, direct scattering compared to the ground state (see Fig. 2h). This trend again fits the observations of larger impact parameters for excited reactions, as more forward scattering means less deflection due to short range collisions. This shows that there is a change from the observed dominantly indirect scattering at 0.7 eV to a mix of both indirect and direct dynamics at 1.2 eV for reactions of vibrationally excited reactants.
At 2.3 eV, the highest investigated collision energy, the indirect isotropically scattered mechanism is almost completely absent. Instead, direct forward stripping dominates (Fig. 2i). There is an additional small backward-scattered component. The forward scattered distribution shows lower velocities for larger angles. Lower velocities are equivalent to higher internal excitation and large scattering angles generally correspond to small impact parameters. Thus small impact parameter collisions lead to higher internal excitation. We have seen a similar effect in cation-neutral charge transfer reactions.37,38 At small impact parameters the short range repulsive part of the interaction potential is explored, which leads to more efficient coupling to vibration and thus to the observed higher internal excitation. The backward contribution is known from theory as a direct rebound mechanism,39 similar to the one known from SN2 reactions.
In contrast to the lower collision energies vibrationally excited reactions show more direct dynamics, compared to the ground state, at the highest collision energy. This can best be identified in the internal energy distribution (Fig. 2k). We find a shift towards lower internal energies, indicative of more direct dynamics. Furthermore the angular distribution is dominated almost exclusively by direct forward scattering. This observation is again interpreted as a larger contribution of large impact parameter collisions leading to stripping with very little angular deflection.
Larger impact parameters for vibrationally excited reactants have been seen before in the reactions of O + CHD3 and O + HCl.40,41 The authors argue that it is evidence for an enlarged cone of acceptance. For the reaction with CHD3, Wang et al. find a change from dominant backward scattering for the ground state to forward and sideways scattering with the CH-stretching vibration excited. They attribute this change to larger impact parameters and argue that the change of dipole moment due to the vibration modifies the long range interaction leading to an attraction at larger distances. This explanation is supported by quasi-classical trajectory calculations on an accurate potential energy surface.20,22,42 A similar steric effect can also explain the suppression of reactivity in the reaction of F + CHD3.21,43 Liu and coworkers refer to these effects as being analogous to an optical lens for approaching reactants steering them towards or away from the preferred transition state geometry.19 It is surprising that we find a similar effect despite the very different strength of the long range interaction due to the ion–dipole interaction in the title reaction.
Collision energy | PT counts IR off | PT counts IR on | Change (%) | Excited fraction (%) | Scaled factor | Theory non/soft/hard ZPE constraint |
---|---|---|---|---|---|---|
0.7 eV | 5673 | 7156 | +26 ± 2.0 | 4.5 ± 0.5 | 6.9 ± 0.8 | 4.2/870/∞ |
1.2 eV | 78713 | 83224 | +5.7 ± 0.5 | 2.7 ± 0.4 | 3.0 ± 0.4 | 2.3/4.7/6.7 |
2.3 eV | 88551 | 89907 | +1.5 ± 0.5 | 2.9 ± 0.3 | 1.5 ± 0.2 | — |
To gain insight how the enhancement depends on the collision energy, we have calculated the branching ratios for proton transfer relative to SN2 products at similar total energies, i.e., collision energy plus energy in the vibrational mode. Assuming no change in the SN2 channel due to vibrational excitation,28 the branching ratios yield insight whether translation or vibration is more efficient in promoting reactivity. At the highest collision energy the assumption of spectator mode SN2 dynamics is not strictly correct, as we recently found a small but significant enhancement of the SN2 channel, which will be discussed in a separate publication. Here, we have calculated the branching ratio relative to the ground state SN2 counts. The results are shown in Fig. 3. The results from the current study are shown in blue for ground state and red for vibrationally excited reactions. The same symbol is chosen for identical collision energies. The grey symbols connected by a line are results from a previous study over a large energy range on the same system by our group27 to illustrate the trend over the investigated energy range. Two additional points taken from ref. 28 are shown in light blue for the ground state and pink with vibrational excitation.
Fig. 3 Proton transfer to SN2 branching ratios versus total energy (Etotal = Ecoll + Evib). The points for excited reactions are scaled to 100% excitation in the neutral beam. The grey points27 and the light blue (ground state) and pink points (vibrationally excited)28 are from previous studies by our group. Points shown in blue are for ground state and in red for vibrationally excited reactions. The same symbol marks the same collision energy. Etotal uncertainties are calculated from the relative energy distributions. The branching ratio uncertainty is dominated by the error on the excited fraction for the excited points and is otherwise smaller than the symbols. |
Our ground state results are in good agreement with the branching ratios of the previous studies. We can now use points of similar total energy to compare whether vibration or translation is more efficient in promoting reactivity. It appears that there is a clear difference between the lowest collision energy and the two higher energy points. Comparing the two excited branching ratios (red and pink) at 1.07 eV with the ground state results at 1.2 eV (grey and blue) suggests that collision energy is more effective in enhancing the proton transfer. At the two higher collision energies we find the opposite effect. The effect that at lower collision energies translation is more efficient in promoting reactivity has been mentioned previously.19 It even occurs in ‘late’ barrier reactions.16,17 Guo and coworkers propose that this occurs, because a minimum of momentum transfer is required for the reaction to proceed, which leads to a preference for translation at low collision energies.24
Comparing the branching ratio without infrared excitation at 1.2 and 1.5 eV collision energy (grey points), we find that the additional collision energy leads to an enhancement by a factor of 2 due to 0.3 eV of additional collision energy. In contrast, vibrational excitation by a similar amount of energy enhances the branching ratio by a factor of 3 (red circle at 1.5 eV total energy). This amounts to a relative vibrational enhancement over translation by a factor of 1.5.
In the collision energy range above 1.5 eV, the branching ratios decrease slightly for the ground state reaction. This coincides with the opening of two halide abstraction channels.27 Vibrational excitation again leads to an increase in the proton transfer to SN2 ratio at 2.5 eV total energy, although this increase is not as strong compared to the data at 1.2 eV collision energy. We still observe a significant increase by a factor of 1.5. The ground and vibrationally excited state reactions therefore show opposite behaviour at the highest investigated collision energy. As the data at the appearance threshold is hardest to evaluate we can conclude from the other two, that the CH-stretching mode is indeed more effective in promoting reactivity towards the proton transfer in this system.
At the intermediate collision energy, theory and experiment should compare better, as above threshold ZPE effects and the influence of experimental uncertainties have less relative contribution. Once again the angular and internal energy distribution agree well for the ground state, but the change in dynamics, evident most prominently in the experimental internal energy distribution (see Fig. 2g), is not reproduced. In the calculations, the Eint distribution is only shifted but is not significantly broadened.
Depending on the chosen ZPE constraint the trajectory calculations give an enhancement factor of the cross section at 1.2 eV collision energy due to symmetric stretching vibration of between 2 and 7 (Table 1). They also provide data at both 1.2 and 1.5 eV, which can be compared to our ground state results. Depending on the chosen ZPE constraint they find an increase of by a factor of 2–3 due to the additional collision energy. The trend that vibrational excitation is more efficient in promoting reactivity is therefore captured both in the quasi-classical trajectory calculations and our experiment. The experimental enhancement factors are slightly larger compared to the non-ZPE values (Table 1), but smaller than the soft-constraint values. It could be interesting to find a constraint that matches the experimental values. In both theory and experiment the vibrational enhancement is collision energy dependent and becomes smaller at high collision energies (Fig. 2 in ref. 26).
The preference for vibration is also seen in calculations using the sudden vector projection model by Guo and coworkers28 for this system. For the symmetric stretching mode they find a value of 0.30 and only a value of 0.03 for translation and therefore a clear preference for stretching vibration over translation, in qualitative agreement with our experimental results.
This study shows that the atomistic reaction dynamics of polyatomic reactions are still far from being understood and that further theoretical work is needed. Our study also represents a step towards manipulating product branching ratios by quantum state preparation in halide anion reactions with methyl halides. The CH-bond cleavage necessary for proton transfer is also an essential step in elimination reactions.44 A similar enhancement by the CH-stretching mode could enable control of the competition between the elimination and nucleophilic substitution pathways.
This journal is © the Owner Societies 2020 |