Unimolecular HBr and HF Elimination Reactions of Vibrationally Excited C2H5CH2Br and C2D5CHFBr: Identification of the 1,1-HBr Elimination Reaction from C2D5CHFBr and Search for the C2D5(F)C:HBr Adduct.

Chemical activation experiments and computational methods have been used to study the unimolecular reactions of C2H5CH2Br and C2D5CHFBr with 90 and 93 kcal mol-1 of vibrational energy, respectively. The four-centered elimination reactions of HBr and DBr are the dominant reactions; however, 2,1-DF, 1,1-HBr, and 1,1-HF reactions are also observed from C2D5CHFBr. The main focus was to search for the role of the C2D5(F)C:HBr adduct in the 1,1-HBr elimination for comparison with carbene adducts in 1,1-HX(Y) elimination from RCHXY (X,Y = Cl and F) molecules. Models of transition states and molecules from electronic structure calculations were used in statistical calculations of the rate constants to assign threshold energies for each reaction based on the experimental rate constants. The threshold energy for 2,1-HBr elimination from 1-bromopropane is 50 kcal mol-1, which is in basic agreement with thermal activation experiments. Comparison of the 2,1-DBr and 2,1-HBr rate constants permits discussion of the kinetic isotope effects and the effect of F atom substitution on the threshold energy for 2,1-HBr elimination. Although CD3CD═CDF from 1,1-HBr elimination of C2D5CHFBr followed by D atom migration is an experimentally observed product, dissociation of the C2D5(F)C:HBr adduct may be the rate-limiting step rather than crossing the barrier associated with the transition state for 1,1-HBr elimination. The calculated dissociation energies of C2H5(X)C:HF adducts are 9.9, 9.3, and 9.0 kcal mol-1 for X = F, Cl, and Br, and the values for C2H5(F)C:HX are 9.9, 6.4, and ∼4.9 kcal mol-1.

Post-transition state dynamics and product energy partitioning following thermal excitation of the F⋯HCH2CN transition state: Disagreement with experiment.

Born-Oppenheimer direct dynamics simulations were performed to study atomistic details of the F + CH3CN → HF + CH2CN H-atom abstraction reaction. The simulation trajectories were calculated with a combined M06-2X/MP2 algorithm utilizing the 6-311++G** basis set. The experiments were performed at 300 K, and assuming the accuracy of transition state theory (TST), the trajectories were initiated at the F⋯HCH2CN abstraction TS with a 300 K Boltzmann distribution of energy and directed towards products. Recrossing of the TS was negligible, confirming the accuracy of TST. HF formation was rapid, occurring within 0.014 ps of the trajectory initiation. The intrinsic reaction coordinate (IRC) for reaction involves rotation of HF about CH2CN and then trapping in the CH2CN⋯HF post-reaction potential energy well of ∼10 kcal/mol with respect to the HF + CH2CN products. In contrast to this IRC, five different trajectory types were observed: the majority proceeded by direct H-atom transfer and only 11% approximately following the IRC. The HF vibrational and rotational quantum numbers, n and J, were calculated when HF was initially formed and they increase as potential energy is released in forming the HF + CH2CN products. The population of the HF product vibrational states is only in qualitative agreement with experiment, with the simulations showing depressed and enhanced populations of the n = 1 and 2 states as compared to experiment. Simulations with an anharmonic zero-point energy constraint gave product distributions for relative translation, HF rotation, HF vibration, CH2CN rotation, and CH2CN vibration as 5%, 11%, 60%, 7%, and 16%, respectively. In contrast, the experimental energy partitioning percentages to HF rotation and vibration are 6% and 41%. Comparisons are made between the current simulation and those for other F + H-atom abstraction reactions. The simulation product energy partitioning and HF vibrational population for F + CH3CN → HF + CH2CN resemble those for other reactions. A detailed discussion is given of possible origins of the difference between the simulation and experimental energy partitioning dynamics for F + CH3CN → HF + CH2CN. The F + CH3CN reaction also forms the CH3C(F)N intermediate, in which the F-atom adds to the C≡N bond. However, this intermediate and F⋯CH3CN and CH3CN⋯F van der Waals complexes are not expected to affect the F + CH3CN → HF + CH2CN product energy partitioning.

Unimolecular reactions of chemically activated CF2BrCF2CH3 and CF2BrCF2CD3: evidence for 1,2-FBr interchange.

Vibrationally excited CF2BrCF2CH3 and CF2BrCF2CD3 molecules were prepared with 96 kcal mol-1 energy at room temperature by the recombination of CF2BrCF2 and CH3 (CD3) radicals. The observed unimolecular reactions are 1,2-BrF interchange to give CF3CFBrCH3 (CD3) molecules and 2,3-FH (FD) elimination; the rate constants are 2.2 x 10(5) (1.5 x 10(5)) s(-1) and 2.0 x 105 (0.75 x 10(5)) s(-1), respectively. The CF3CFBrCH3 (CD3) molecules rapidly, relative to the reverse reaction, eliminate HBr or DBr to give the observed product CF3CF=CH2 (CD2). Density functional theory at the B3PW91/6-311+G(2d,p) level was used to obtain vibrational frequencies and moments of inertia of the molecule and transition states for subsequent calculations of statistical rate constants for CF2BrCF2CH3 and CF2BrCF2CD3. Matching experimental and calculated rate constants gave threshold energies of 62 and 66 kcal mol-1 for 1,2-BrF interchange and 2,3-FH elimination, respectively. The BrF interchange reaction is compared to ClF interchange from CF2ClCF2CH3 and CF2ClCHFCH3.

Use of a single trajectory to study product energy partitioning in unimolecular dissociation: mass effects for halogenated alkanes.

A single trajectory (ST) direct dynamics approach is compared with quasiclassical trajectory (QCT) direct dynamics calculations for determining product energy partitioning in unimolecular dissociation. Three comparisons are made by simulating C(2)H(5)F-->HF + C(2)H(4) product energy partitioning for the MP26-31G(*) and MP26-311 + + G(**) potential energy surfaces (PESs) and using the MP26-31G(*) PES for C(2)H(5)F dissociation as a model to simulate CHCl(2)CCl(3)-->HCl + C(2)Cl(4) dissociation and its product energy partitioning. The trajectories are initiated at the transition state with fixed energy in reaction-coordinate translation E(t) (double dagger). The QCT simulations have zero-point energy (ZPE) in the vibrational modes orthogonal to the reaction coordinate, while there is no ZPE for the STs. A semiquantitative agreement is obtained between the ST and QCT average percent product energy partitionings. The ST approach is used to study mass effects for product energy partitioning in HX(X = F or Cl) elimination from halogenated alkanes by using the MP26-31G(*) PES for C(2)H(5)F dissociation and varying the masses of the C, H, and F atoms. There is, at most, only a small mass effect for partitioning of energy to HX vibration and rotation. In contrast, there are substantial mass effects for partitioning to relative translation and the polyatomic product's vibration and rotation. If the center of mass of the polyatomic product is located away from the C atom from which HX recoils, the polyatomic has substantial rotation energy. Polyatomic products, with heavy atoms such as Cl atoms replacing the H atoms, receive substantial vibration energy that is primarily transferred to the wag-bend motions. For E(t) (double dagger) of 1.0 kcalmol, the ST calculations give average percent partitionings to relative translation, polyatomic vibration, polyatomic rotation, HX vibration, and HX rotation of 74.9%, 6.8%, 1.5%, 14.4%, and 2.4% for C(2)H(5)F dissociation and 39.7%, 38.1%, 0.2%, 16.1%, and 5.9% for a model of CHCl(2)CCl(3) dissociation.

Comparison of levels of electronic structure theory in direct dynamics simulations of C2H5F --> HF + C2H4 product energy partitioning.

Direct dynamics simulations at the MP2/6-311++G** level of theory were performed to study C(2)H(5)F --> HF + C(2)H(4) product energy partitioning. The simulation results are compared with experiment and a previous MP2/6-31G* simulation. The current simulation with the larger basis set releases more energy to HF vibration and less to HF + C(2)H(4) relative translation as compared to the previous simulation with the 6-31G* basis set. The HF rotation and vibration energy distributions determined from the current simulation are in overall very good agreement with previous experimental studies of C(2)H(5)F dissociation by chemical activation and IRMPA. A comparison of the simulations with experiments suggests there may be important mass effects for energy partitioning in HX elimination from haloalkanes. The transition state (TS) structures and energies calculated with MP2 and the 6-31G* and 6-311++G** basis sets are compared with those calculated using CCD, CCSD, CCSD(T), and the 6-311++G** basis set.