State-to-state dynamics of the Cl + CH3OH --> HCl + CH2OH reaction

Molecular chlorine, methanol, and helium are co-expanded into a vacuum chamber using a custom designed "late-mixing" nozzle. The title reaction is initiated by photolysis of Cl2 at 355 nm, which generates monoenergetic Cl atoms that react with CH3OH at a collision energy of 1960 +/- 170 cm(-1) (0.24 +/- 0.02 eV). Rovibrational state distributions of the nascent HCl products are obtained via 2 + 1 resonance enhanced multiphoton ionization, center-of-mass scattering distributions are measured by the core-extraction technique, and the average internal energy of the CH3OH co-products is deduced by measuring the spatial anisotropy of the HCl products. The majority (84 +/- 7%) of the HCl reaction products are formed in HCl(v = 0) with an average rotational energy of [Erot] = 390 +/- 70 cm(-1). The remaining 16 +/- 7% are formed in HCl(v = 1) and have an average rotational energy of [Erot] = 190 +/- 30 cm(-1). The HCl(v = 1) products are primarily forward scattered, and they are formed in coincidence with CH2OH products that have little internal energy. In contrast, the HCl(v = 0) products are formed in coincidence with CH2OH products that have significant internal energy. These results indicate that two or more different mechanisms are responsible for the dynamics in the Cl + CH3OH reaction. We suggest that (1) the HCl(v = 1) products are formed primarily from collisions at high impact parameter via a stripping mechanism in which the CH2OH co-products act as spectators, and (2) the HCl(v = 0) products are formed from collisions over a wide range of impact parameters, resulting in both a stripping mechanism and a rebound mechanism in which the CH2OH co-products are active participants. In all cases, the reaction of fast Cl atoms with CH3OH is with the hydrogen atoms on the methyl group, not the hydrogen on the hydroxyl group.     

High-temperature rate constants for CH3OH + Kr --> products, OH + CH3OH --> products, OH + (CH3)(2)CO --> CH2COCH3 + H2O, and OH + CH3 --> CH) + H2O

The reflected shock tube technique with multipass absorption spectrometric detection of OH radicals at 308 nm (corresponding to a total path length of approximately 4.9 m) has been used to study the dissociation of methanol between 1591 and 2865 K. Rate constants for two product channels [CH3OH + Kr --> CH3 + OH + Kr (1) and CH3OH + Kr --> 1CH2 + H2O + Kr (2)] were determined. During the course of the study, it was necessary to determine several other rate constants that contributed to the profile fits. These include OH + CH3OH --> products, OH + (CH3)2CO --> CH2COCH3 + H2O, and OH + CH3 --> 1,3CH2 + H2O. The derived expressions, in units of cm(3) molecule(-1) s(-1), are k(1) = 9.33 x 10(-9) exp(-30857 K/T) for 1591-2287 K, k(2) = 3.27 x 10(-10) exp(-25946 K/T) for 1734-2287 K, kOH+CH3OH = 2.96 x 10-16T1.4434 exp(-57 K/T) for 210-1710 K, k(OH+(CH3)(2)CO) = (7.3 +/- 0.7) x 10(-12) for 1178-1299 K and k(OH+CH3) = (1.3 +/- 0.2) x 10(-11) for 1000-1200 K. With these values along with other well-established rate constants, a mechanism was used to obtain profile fits that agreed with experiment to within <+/-10%. The values obtained for reactions 1 and 2 are compared with earlier determinations and also with new theoretical calculations that are presented in the preceding article in this issue. These new calculations are in good agreement with the present data for both (1) and (2) and also for OH + CH3 --> products.     

Quantum Monte Carlo study of the reaction: Cl+CH3OH-->CH2OH+HCl

A theoretical study is reported of the Cl+CH3OH-->CH2OH+HCl reaction based on the diffusion Monte Carlo (DMC) variant of the quantum Monte Carlo method. Using a DMC trial function constructed as a product of Hartree-Fock and correlation functions, we have computed the barrier height, heat of reaction, atomization energies, and heats of formation of reagents and products. The DMC heat of reaction, atomization energies, and heats of formation are found to agree with experiment to within the error bounds of computation and experiment. M?ller-Plesset second order perturbation theory (MP2) and density functional theory, the latter in the B3LYP generalized gradient approximation, are found to overestimate the experimental heat of reaction. Intrinsic reaction coordinate calculations at the MP2 level of theory demonstrate that the reaction is predominantly direct, i.e., proceeds without formation of intermediates, which is consistent with a recent molecular beam experiment. The reaction barrier as determined from MP2 calculations is found to be 2.24 kcal/mol and by DMC it is computed to be 2.39(49) kcal/mol.     

Imaging photon-initiated reactions: a study of the Cl(2P(3/2))+CH4-->HCl+CH3 reaction

The hydrogen or deuterium atom abstraction reactions between Cl((2)P(3/2)) and methane, or its deuterated analogues CD(4) and CH(2)D(2), have been studied at mean collision energies around 0.34 eV. The experiments were performed in a coexpansion of molecular chlorine and methane in helium, with the atomic Cl reactants generated by polarized laser photodissociation of Cl(2) at 308 nm. The Cl-atom reactants and the methyl radical products were detected using (2+1) resonantly enhanced multiphoton ionization, coupled with velocity-map ion imaging. Analysis of the ion images reveals that in single-beam experiments of this type, careful consideration must be given to the spread of reagent velocities and collision energies. Using the reactions of Cl with CH(4), CD(4), and CH(2)D(2), as examples, it is shown that the data can be fitted well if the reagent motion is correctly described, and the angular scattering distributions can be obtained with confidence. New evidence is also provided that the CD(3) radicals from the Cl+CD(4) reaction possess significant rotational alignment under the conditions of the present study. The results are compared with previous experimental and theoretical works, where these are available.     

Reflected shock tube studies of high-temperature rate constants for OH + CH4 --> CH3 + H2O and CH3 + NO2 --> CH3O + NO

The reflected shock tube technique with multipass absorption spectrometric detection of OH radicals at 308 nm has been used to study the reactions OH + CH(4) --> CH(3) + H(2)O and CH(3) + NO(2) --> CH(3)O + NO. Over the temperature range 840-2025 K, the rate constants for the first reaction can be represented by the Arrhenius expression k = (9.52 +/- 1.62) x 10(-11) exp[(-4134 +/- 222 K)/T] cm(3) molecule(-1) s(-1). Since this reaction is important in both combustion and atmospheric chemistry, there have been many prior investigations with a variety of techniques. The present results extend the temperature range by 500 K and have been combined with the most accurate earlier studies to derive an evaluation over the extended temperature range 195-2025 K. A three-parameter expression describes the rate behavior over this temperature range, k = (1.66 x 10(-18))T(2.182) exp[(-1231 K)/T] cm(3) molecule(-1) s(-1). Previous theoretical studies are discussed, and the present evaluation is compared to earlier theoretical estimates. Since CH(3) radicals are a product of the reaction and could cause secondary perturbations in rate constant determinations, the second reaction was studied by OH radical production from the fast reactions CH(3)O --> CH(2)O + H and H + NO(2) --> OH + NO. The measured rate constant is 2.26 x 10(-11) cm(3) molecule(-1) s(-1) and is not dependent on temperature from 233 to 1700 K within experimental error.     

Ab initio and direct quasiclassical-trajectory study of the Cl + CH4-->HCl + CH3 reaction

We present an electronic structure and dynamics study of the Cl + CH(4)--> HCl + CH(3) reaction. We have characterized the stationary points of the ground-state potential-energy surface using various electronic structure methods and basis sets. Our best calculations, CCSD(T) extrapolated to the complete basis-set limit based on geometries and harmonic frequencies obtained at the CCSD(T)/aug-cc-pvtz level, are in agreement with the experimental reaction energy and indirect measurements of the barrier height. Using ab initio information, we have reparametrized a semiempirical Hamiltonian so that the predictions of the improved Hamiltonian agree with the higher-level calculations in various regions of the potential-energy surface. This improved semiempirical Hamiltonian is then used to propagate quasiclassical trajectories and characterize the reaction dynamics. The good agreement of the calculated HCl rotational and angular distributions with the experiment indicates that reparametrizing semiempirical Hamiltonians is a promising approach to derive accurate potential-energy surfaces for polyatomic reactions. However, excessive energy leakage from the initial vibrational energy of the CH(4) molecule to the reaction coordinate in the trajectory calculations calls into question the suitability of the standard quasiclassical-trajectory method to describe energy partitioning in polyatomic reactions.     

Quantum calculations of the O(3P)+H2-->OH+H reaction

Quantum scattering calculations are reported for the O(3P)+H2(v=0,1) reaction using chemically accurate potential energy surfaces of 3A' and 3A" symmetry. We present state-to-state reaction cross sections and rate coefficients as well as thermal rate coefficients for the title reaction using accurate quantum calculations. Our calculations yield reaction cross sections that are in quantitative accord with results of recent crossed molecular beam experiments. Comparisons with results obtained using the J-shifting calculations show that the J-shifting approximation is quite reliable for this system. Thermal rate coefficients from the exact calculations and the J-shifting approximation agree remarkably well with experimental results. Our calculations also reproduce the markedly different OH(v'=0)/OH(v'=1) branching in O(3P)+H2(v=1) reaction, observed in experiments that use different O(3P) atom sources. In particular, we show that the branching ratio is a strong function of the kinetic energy of the O(3P) atom.     

Vibrational spectroscopic study of the electron transfer for the Cl- + CH3Cl S(N)2 reaction in the gas phase

The vibrational spectra study of electron transfer of Cl- + CH3Cl --> ClCH3 + Cl- reaction has been examined by density function theory (DFT) calculations at 6-311 + +G** level in this paper. This reaction includes old bond rupture, new bond formation and electronic transfer in the intermolecular. The vibrational frequencies and vibrational modes of reactant, precursor complex, transition state, successor complex and product are analyzed. The relationship and the change among them can confirm the rupture of bond, the formation of bond and the process of electron transfer.     

Effect of bending and torsional mode excitation on the reaction Cl+CH4-->HCl+CH3

A beam containing CH(4), Cl(2), and He is expanded into a vacuum chamber where CH(4) is prepared via infrared excitation in a combination band consisting of one quantum of excitation each in the bending and torsional modes (nu(2)+nu(4)). The reaction is initiated by fast Cl atoms generated by photolysis of Cl(2) at 355 nm, and the resulting CH(3) and HCl products are detected in a state-specific manner using resonance-enhanced multiphoton ionization (REMPI). By comparing the relative amplitudes of the action spectra of Cl+CH(4)(nu(2)+nu(4)) and Cl+CH(4)(nu(3)) reactions, we determine that the nu(2)+nu(4) mode-driven reaction is at least 15% as reactive as the nu(3) (antisymmetric stretch) mode-driven reaction. The REMPI spectrum of the CH(3) products shows no propensity toward the formation of umbrella bend mode excited methyl radical, CH(3)(nu(2)=1), which is in sharp distinction to the theoretical expectation based on adiabatic correlations between CH(4) and CH(3). The rotational distribution of HCl(v=1) products from the Cl+CH(4)(nu(2)+nu(4)) reaction is hotter than the corresponding distribution from the Cl+CH(4)(nu(3)) reaction, even though the total energies of the two reactions are the same within 4%. An explanation for this enhanced rotational excitation of the HCl product from the Cl+CH(4)(nu(2)+nu(4)) reaction is offered in terms of the projection of the bending motion of the CH(4) reagent onto the rotational motion of the HCl product. The angular distributions of the HCl(nu=0) products from the Cl+CH(4)(nu(2)+nu(4)) reaction are backward scattered, which is in qualitative agreement with theoretical calculation. Overall, nonadiabatic product vibrational correlation and mode specificity of the reaction indicate that either the bending mode or the torsional mode or both modes are strongly coupled to the reaction coordinate.     

Quasiclassical trajectory study of the O(3P) + CH4 --> OH + CH3 reaction with a specific reaction parameters semiempirical Hamiltonian

We present a theoretical study of the O(3P) + CH4 --> OH + CH3 reaction using electronic structure, kinetics, and dynamics calculations. We calculate a grid of ab initio points at the PMP2/AUG-cc-pVDZ level to characterize the potential energy surface in regions of up to 1.3 eV above reagents. This grid of ab initio points is used to derive a set of specific reaction parameters (SRP) for the MSINDO semiempirical Hamiltonian. The resulting SRP-MSINDO Hamiltonian improves the quality of the standard Hamiltonian, particularly in regions of the potential energy surface beyond the minimum-energy reaction path. Quasiclassical-trajectory calculations are used to study the reaction dynamics with the original and the improved MSINDO semiempirical Hamiltonians, and a prior surface. The SRP-MSINDO semiempirical Hamiltonian yields OH rotational distributions in agreement with experimental results, improving over the results of the other surfaces. Thermal rate constants estimated with Variational Transition State Theory using the SRP-MSINDO Hamiltonian are also in agreement with experiments. Our results indicate that reparametrized semiempirical Hamiltonians are a good alternative to generating potential energy surfaces for accurate dynamics studies of polyatomic reactions.     

Dependence of the rate constants on the treatment of internal rotation modes: the reaction OH + CH3SH --> CH3S + H2O as an example

The title reaction has been used as an example to test the importance of using a hindered rotor treatment instead of a harmonic oscillator model for calculating vibrational partition functions corresponding to low-frequency internal rotation modes. First, a normal-mode analysis according to the Ayala and Schlegel's algorithm has been used to identify the internal rotation modes of methanethiol and the transition state structure. Then, after calculation of the energy barrier for each internal rotation, the corresponding hindered rotor partition functions have been calculated following the CW scheme of Chuang and Truhlar. The results show that the anharmonic treatment produces a rather modest improvement of the rate constants at room temperature or below.     

Potential energy surface, kinetics, and dynamics study of the Cl+CH4-->HCl+CH3 reaction

A modified and recalibrated potential energy surface for the gas-phase Cl+CH4-->HCl+CH3 reaction is reported and tested. It is completely symmetric with respect to the permutation of the four methane hydrogen atoms and is calibrated with respect to updated experimental and theoretical stationary point properties and experimental forward thermal rate constants. From the kinetics point of view, the forward and reverse thermal rate constants and the activation energies were calculated using the variational transition-state theory with semiclassical transmission coefficients over a wide temperature range of 150-2500 K. The theoretical results reproduce the available experimental data, with a small curvature of the Arrhenius plot which indicates the role of tunneling in this hydrogen abstraction reaction. A dynamics study was also performed on this PES using quasiclassical trajectory (QCT) calculations, including corrections to avoid zero-point energy leakage along the trajectories. First, we found a noticeable internal energy in the coproduct methyl radical, both in the ground-state [CH4 (v=0)] and vibrationally excited [CH4 (v=1)] reactions. This CH3 internal energy was directly precluded in some experiments or oversimplified in previous theoretical studies using pseudotriatomic models. Second, our QCT calculations give HCl rotational distributions slightly hotter than those in experiment, but correctly describing the experimental trend of decreasing the HCl product rotation excitation in going from HCl (v'=0) to HCl (v'=1) for the CH4 (v=1) reaction. Third, the state specific scattering distributions present a reasonable agreement with experiment, although they tend to make the reaction more forward and backward scattered than found experimentally probably because of the hotter rotational distribution and the deficiencies of the QCT methods.     

Searching for resonances in the reaction Cl+CH4-->HCl+CH3: quantum versus quasiclassical dynamics and comparison with experiments

A quantum-mechanical (QM) and quasiclassical trajectory (QCT) study was performed on the title reaction, using a pseudotriatomic ab initio based surface. Probabilities and integral cross sections present some clear peaks versus the collision energy E(col), which we assign to Feshbach resonances of the transition state, where the light H atom oscillates between the heavy Cl and CH(3) groups. For ground-state reactants, reactivity is essentially of quantum origin (QCT observables and oscillations are smaller, or much smaller, than QM ones), and the calculated integral cross section and product distributions are in reasonable agreement with the experiment. The reaction occurs through an abstraction mechanism, following both a direct and an indirect mechanism. The quasiclassical trajectory calculations show the participation of a short-lived collision complex in the microscopic reaction mechanism. Finally, QCT differential cross sections of Cl+CH(4)-->HCl (nu(')=0 and 1)+CH(3) oscillate versus E(col), whereas experimentally this only occurs for HCl (nu(')=1). This theoretical result and other oscillating properties found here could, however, be related to the existence of a Feshbach resonance for the production of HCl (nu(')=1), as suggested by experimentalists.     

Differential and integral cross sections for the H + O2 --> OH + O combustion reaction

We present accurate differential and integral cross sections for the H + O2 --> OH + O reaction obtained on a newly developed ab initio potential energy surface using time-independent and time-dependent quantum mechanical methods. The product angular distributions near the reaction threshold show pronounced forward and backward peaks, reflecting the complex-forming mechanism. However, the asymmetry of these peaks suggests certain nonstatistical behaviors, presumably due to some relatively short-lived resonances. The integral cross section increases monotonically with the collision energy above a reaction threshold.     

Theoretical study of the complex-forming CH + H2 --> CH2 + H reaction

The complex-forming CH + H2 --> CH2 + H reaction is studied employing a recently developed global potential energy function. The reaction probability in the total angular momentum J = 0 limit is estimated with a four-atom quantum wave packet method and compared with classical trajectory and statistical theory results. The formation of complexes from different reactant internal states is also determined with wave packet calculations. While there is no barrier to reaction along the minimum energy path, we find that there are angular constraints to complex formation. Trajectory-based estimates of the low-pressure rate constants are made and compared with experimental results. We find that zero-point energy violation in the trajectories is a particularly severe problem for this reaction.     

Quantum mechanical rate constants for H + O2 <--> O + OH and H + O2 --> HO2 reactions

Canonical rate constants for both the forward and reverse H + O(2) <--> O + OH reactions were calculated using a quantum wave packet-based statistical model on the DMBE IV potential energy surface of Varandas and co-workers. For these bimolecular reactions, the results show reasonably good agreement with available experimental and theoretical data up to 1500 K. In addition, the capture rate for the H + O(2) --> HO(2) addition reaction at the high-pressure limit was obtained on the same potential using a time-independent quantum capture method. Excellent agreement with experimental and quasi-classical trajectory results was obtained except for at very low temperatures, where a reaction threshold was found and attributed to the centrifugal barrier of the orbital motion.     

Kinetics study of the OH + alkene --> H2O + alkenyl reaction class

In this paper we report on the kinetics of hydrogen abstraction for the OH + alkene reaction class, using the reaction class transition state theory (RC-TST) combined with the linear energy relationship (LER) and the barrier height grouping (BHG) approaches. Parameters for the RC-TST were derived from theoretical calculations using a set of 15 reactions representing the hydrogen abstractions from the terminal and nonterminal carbon sites of the double bond of alkene compounds. Both the RC-TST/LER, where only reaction energy is needed at either density functional theory BH&HLYP or semiempirical AM1 levels, and RC-TST/BHG, where no additional information is required, are found to be promising methods for predicting rate constants for a large number of reactions in this reaction class. Detailed error analyses show that, when compared to explicit theoretical calculations, the averaged systematic errors in the calculated rate constants using both the RC-TST/LER and RC-TST/BHG methods are less than 25% in the temperature range 300-3000 K. The estimated rate constants using these approaches are in good agreement with available data in the literature.