Quasiclassical trajectory study of the Cl+CH4 reaction dynamics on a quadratic configuration interaction with single and double excitation interpolated potential energy surface

An ab initio interpolated potential energy surface (PES) for the Cl+CH(4) reactive system has been constructed using the interpolation method of Collins and co-workers [J. Chem. Phys. 102, 5647 (1995); 108, 8302 (1998); 111, 816 (1999); Theor. Chem. Acc. 108, 313 (2002)]. The ab initio calculations have been performed using quadratic configuration interaction with single and double excitation theory to build the PES. A simple scaling all correlation technique has been used to obtain a PES which yields a barrier height and reaction energy in good agreement with high level ab initio calculations and experimental measurements. Using these interpolated PESs, a detailed quasiclassical trajectory study of integral and differential cross sections, product rovibrational populations, and internal energy distributions has been carried out for the Cl+CH(4) and Cl+CD(4) reactions, and the theoretical results have been compared with the available experimental data. It has been shown that the calculated total reaction cross sections versus collision energy for the Cl+CH(4) and Cl+CD(4) reactions is very sensitive to the barrier height. Besides, due to the zero-point energy (ZPE) leakage of the CH(4) molecule to the reaction coordinate in the quasiclassical trajectory (QCT) calculations, the reaction threshold falls below the barrier height of the PES. The ZPE leakage leads to CH(3) and HCl coproducts with internal energy below its corresponding ZPEs. We have shown that a Gaussian binning (GB) analysis of the trajectories yields excitation functions in somehow better agreement with the experimental determinations. The HCl(v'=0) and DCl(v'=0) rotational distributions are as well very sensitive to the ZPE problem. The GB correction narrows and shifts the rotational distributions to lower values of the rotational quantum numbers. However, the present QCT rotational distributions are still hotter than the experimental distributions. In both reactions the angular distributions shift from backward peaked to sideways peaked as collision energy increases, as seen in the experiments and other theoretical calculations.     

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.     

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.     

Unravelling the reactivity of antisymmetric stretch-excited CH4 with Cl by-product pair-correlation measurements

The vibrational branching ratio [CH(3)(v=0) + HCl(v'=1)]/[CH(3)(v=0) + HCl(v'=0)] of two correlated product pairs from the title reaction shows a dramatic E(c)-dependence, in sharp contrast to the previously observed behavior in Cl + CHD(3)(v(1)=1), while the vibrational enhancement factors in reactivity of the two reactions are remarkably similar.     

Quantum and quasiclassical studies of the O(3P)+HCl-->OH+Cl(2P) reaction using benchmark potential surfaces

We have performed quantum mechanical (QM) dynamics calculations within the independent-state approximation with new benchmark triplet A" and A' surfaces [B. Ramachandran et al., J. Chem. Phys. 119, 9590 (2003)] for the rovibronic state-to-state measurements of the reaction O(3P)+HCl(v=2,j=1,6,9)-->OH(v'j')+Cl(2P) [Zhang et al., J. Chem. Phys. 94, 2704 (1991)]. The QM and experimental rotational distributions peak at similar OH(j') levels, but the QM distributions are significantly narrower than the measurements and previous quasiclassical dynamics studies. The OH(low j) populations observed in the measurements are nearly absent in the QM results. We have also performed quasiclassical trajectory with histogram binning (QCT-HB) calculations on these same benchmark surfaces. The QCT-HB rotational distributions, which are qualitatively consistent with measurements and classical dynamics studies using other surfaces, are much broader than the QM results. Application of a Gaussian binning correction (QCT-GB) dramatically narrows and shifts the QCT-HB rotational distributions to be in very good agreement with the QM results. The large QCT-GB correction stems from the special shape of the joint distribution of the classical rotational/vibrational action of OH products. We have also performed QM and QCT calculations for the transition, O+HCl(v=0,T=300 K)-->OH(v'j')+Cl from threshold to approximately 130 kcal mol(-1) collision energy as a guide for possible future hyperthermal O-atom measurements. We find in general a mixed energy release into translation and rotation consistent with a late barrier to reaction. Angular distributions at high collision energy are forward peaked, consistent with a stripping mechanism. Direct collisional excitation channel cross sections, O+HCl(v=0,T=300 K)-->O+HCl(v'=1), in the same energy range are large, comparable in magnitude to the reactive channel cross sections. Although the (3)A" state dominates most collision processes, above approximately 48 kcal mol(-1), the (3)A' state plays the major role in collisional excitation.     

Comparing the dynamical effects of symmetric and antisymmetric stretch excitation of methane in the Cl+CH4 reaction

The effects of two nearly isoenergetic C-H stretching motions on the gas-phase reaction of atomic chlorine with methane are examined. First, a 1:4:9 mixture of Cl(2), CH(4), and He is coexpanded into a vacuum chamber. Then, either the antisymmetric stretch (nu(3)=3019 cm(-1)) of CH(4) is prepared by direct infrared absorption or the infrared-inactive symmetric stretch (nu(1)=2917 cm(-1)) of CH(4) is prepared by stimulated Raman pumping. Photolysis of Cl(2) at 355 nm generates fast Cl atoms that initiate the reaction with a collision energy of 1290+/-175 cm(-1) (0.16+/-0.02 eV). Finally, the nascent HCl or CH(3) products are detected state-specifically via resonance enhanced multiphoton ionization and separated by mass in a time-of-flight spectrometer. We find that the rovibrational distributions and state-selected differential cross sections of the HCl and CH(3) products from the two vibrationally excited reactions are nearly indistinguishable. Although Yoon et al. [J. Chem. Phys. 119, 9568 (2003)] report that the reactivities of these two different types of vibrational excitation are quite different, the present results indicate that the reactions of symmetric-stretch excited or antisymmetric-stretch excited methane with atomic chlorine follow closely related product pathways. Approximately 37% of the reaction products are formed in HCl(v=1,J) states with little rotational excitation. At low J states these products are sharply forward scattered, but become almost equally forward and backward scattered at higher J states. The remaining reaction products are formed in HCl(v=0,J) and have more rotational excitation. The HCl(v=0,J) products are predominantly back and side scattered. Measurements of the CH(3) products indicate production of a non-negligible amount of umbrella bend excited methyl radicals primarily in coincidence with the HCl(v=0,J) products. The data are consistent with a model in which the impact parameter governs the scattering dynamics.     

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.     

Quasiclassical trajectory study of the C(1D) + H2 reaction and isotopomeric variants: kinetic isotope effect and CD/CH branching ratio

The recently proposed ab initio single-sheeted double many-body expansion potential energy for the methylene molecule has been used to perform quasiclassical trajectory (QCT) calculations for the title reaction. Thermal and initial state-specific (v = 0, j = 0) rate constants for the C((1)D) + H(2)/HD/D(2) reactions have been obtained over a wide range of temperatures. Cross sections for the reaction C((1)D) + H(2) and its deuterated isotopes have also been calculated, as well as the CD/CH branching ratios for the C((1)D) + HD reaction. It is found that the CD + H product channel in the C((1)D) + HD reaction is preferred relative to the CH + D channel. The estimated rate constants are predicted to be in the order k(H2) > k(HD) > k(D2) and the calculated cross sections and rate constants compared with available theoretical and experimental data.     

Bond and mode selectivity in the reaction of atomic chlorine with vibrationally excited CH2D2

The title reaction is investigated by co-expanding a mixture of Cl2 and CH2D2 into a vacuum chamber and initiating the reaction by photolyzing Cl2 with linearly polarized 355 nm light. Excitation of the first C-H overtone of CH2D2 leads to a preference for hydrogen abstraction over deuterium abstraction by at least a factor of 20, whereas excitation of the first C-D overtone of CH2D2 reverses this preference by at least a factor of 10. Reactions with CH2D2 prepared in a local mode containing two quanta in one C-H oscillator /2000>- or in a local mode containing one quantum each in two C-H oscillators /1100> lead to products with significantly different rotational, vibrational, and angular distributions, although the vibrational energy for each mode is nearly identical. The Cl+CH2D2/2000>- reaction yields methyl radical products primarily in their ground state, whereas the Cl+CH2D2/1100> reaction yields methyl radical products that are C-H stretch excited. The HCl(v=1) rotational distribution from the Cl+CH2D2/2000>- reaction is significantly hotter than the HCl(v=1) rotational distribution from the Cl+CH2D2/1100> reaction, and the HCl(v=1) differential cross-section (DCS) of the Cl+CH2D2/2000>- reaction is more broadly side scattered than the HCl(v=1) DCS of the Cl+CH2D2/1100> reaction. The results can be explained by a simple spectator model and by noting that the /2000>- mode leads to a wider cone of acceptance for the reaction than the /1100> mode. These measurements represent the first example of mode selectivity observed in a differential cross section, and they demonstrate that vibrational excitation can be used to direct the reaction pathway of the Cl+CH2D2 reaction.     

Ab initio potential energy surface, variational transition state theory, and quasiclassical trajectory studies of the F+CH4-->HF+CH3 reaction

In this work we present a study of the F+CH(4)-->HF+CH(3) reaction (DeltaHdegrees(298 K)=-32.0 kcal mol(-1)) using different methods of the chemical reaction theory. The ground potential energy surface (PES) is characterized using several ab initio methods. Full-dimensional rate constants have been calculated employing the variational transition state theory and using directly ab initio data. A triatomic analytical representation of the ground PES was derived from ab initio points calculated at the second- and fourth-order M?ller-Plesset levels with the 6-311+G(2df,2pd) basis set, assuming the CH(3) fragment to be a 15 a.m.u. pseudoatom in the fitting process. This is suggested from experiments that indicate that the methyl group is uncoupled to the reaction coordinate. A dynamics study by means of the quasiclassical trajectory (QCT) method and employing this analytical surface was also carried out. The experimental data available on the HF internal states distributions are reproduced by the QCT results. Very recent experimental information about the reaction stereodynamics is also borne out by our QCT calculations. Comparisons with the benchmark F+H(2) and analogous Cl+CH(4) reactions are established throughout.     

Imaging the nature of the mode-specific chemistry in the reaction of Cl atom with antisymmetric stretch-excited CH4

Following up our preliminary communication [Kawamata et al., Phys. Chem. Chem. Phys. 10, 4378 (2008)], the effects of the antisymmetric-stretching excitation of methane on the Cl((2)P(3/2))+CH(4) reaction are examined here over a wide range of initial collision energy in a crossed molecular beam imaging experiment. The antisymmetric stretch of CH(4) is prepared in a single rovibrational state of (v(3)=1, j=2) by direct infrared absorption, and the major product states of CH(3)(v=0) are probed by a time-sliced velocity-map imaging method. We find that at fixed collision energies, the stretching excitation promotes reaction rate. Compared to the ground-state reaction, this vibrational enhancement factor is, however, no more effective than the translational enhancement. The correlated HCl(v'=1) vibrational branching fraction shows a striking dependence on collision energies, varying from 0.7 at E(c)=2 kcal mol(-1) to about 0.2 at 13 kcal mol(-1). This behavior resembles the previously studied Cl+CH(2)D(2)(v(6)=1), but is in sharp contrast to the Cl+CHD(3)(v(1)=1) and CH(2)D(2)(v(1)=1) reactions. Dependences of experimental results on the probed rotational states of CH(3)(v=0) are also elucidated. We qualitatively interpret those experimental observations based on a conceptual framework proposed recently.     

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.     

The dynamics of the Cl+n-C4H10-->HCl (v',j') + C4H9 reaction at 0.32 eV

Rotational state resolved center-of-mass angular scattering and kinetic energy release distributions have been determined for the HCl (v' = 0, j' = 0-6) products of the reaction of chlorine with n-butane using the photon-initiated reaction technique, coupled with velocity-map ion imaging. The angular and kinetic energy release distributions derived from the ion images are very similar to those obtained previously for the Cl plus ethane reaction. The angular distributions are found to shift from forward scattering to more isotropic scattering with increasing HCl rotational excitation. The kinetic energy release distributions indicate that around 30% of the available energy is channeled into internal excitation of the butyl radical products. The data analysis also suggests that H-atom abstraction takes place from both primary and secondary carbon atom sites, with the primary site producing rotationally cold, forward scattered HCl (v' = 0) products, and the secondary site yielding more isotropically scattered HCl (v' = 0) possessing higher rotational excitation. The mechanisms leading to these two product channels are discussed in the light of the present findings, and in comparison with studies of other Cl plus alkane reactions.     

Reaction dynamics of Cl+CH3SH: rotational and vibrational distributions of HCl probed with time-resolved Fourier-transform spectroscopy

Rotationally resolved infrared emission spectra of HCl(v=1-3) in the reaction of Cl+CH3SH, initiated with radiation from a laser at 308 nm, are detected with a step-scan Fourier-transform spectrometer. Observed rotational temperature of HCl(v=1-3) decreases with duration of reaction due to collisional quenching; a short extrapolation to time zero based on data in the range 0.25-4.25 micros yields a nascent rotational temperature of 1150+/-80 K. The rotational energy averaged for HCl(v=1-3) is 8.2+/-0.9 kJ mol(-1), yielding a fraction of available energy going into rotation of HCl, fr=0.10+/-0.01, nearly identical to that of the reaction Cl+H(2)S. Observed temporal profiles of the vibrational population of HCl(v=1-3) are fitted with a kinetic model of formation and quenching of HCl(v=1-3) to yield a branching ratio (68+/-5):(25+/-4):(7+/-1) for formation of HCl(v=1):(v=2):(v=3) from the title reaction and its thermal rate coefficient k(2a)=(2.9+/-0.7)x10(-10) cm(3) molecule(-1) s(-1). Considering possible estimates of the vibrational population of HCl(v=0) based on various surprisal analyses, we report an average vibrational energy 36+/-6 kJ mol(-1) for HCl. The fraction of available energy going into vibration of HCl is f(v)=0.45+/-0.08, significantly greater than a value fv=0.33+/-0.06 determined previously for Cl+H2S. Reaction dynamics of Cl+H(2)S and Cl+CH3SH are compared; the adduct CH3S(Cl)H is likely more transitory than the adduct H(2)SCl.     

Quantum and classical studies of the O(3P) + H2(v = 0-3,j = 0) --> OH + H reaction using benchmark potential surfaces

We present results of time-dependent quantum mechanics (TDQM) and quasiclassical trajectory (QCT) studies of the excitation function for O(3P) + H2(v = 0-3,j = 0) --> OH + H from threshold to 30 kcal/mol collision energy using benchmark potential energy surfaces [Rogers et al., J. Phys. Chem. A 104, 2308 (2000)]. For H2(v = 0) there is excellent agreement between quantum and classical results. The TDQM results show that the reactive threshold drops from 10 kcal/mol for v = 0 to 6 for v = 1, 5 for v = 2 and 4 for v = 3, suggesting a much slower increase in rate constant with vibrational excitation above v = 1 than below. For H2(v > 0), the classical results are larger than the quantum results by a factor approximately 2 near threshold, but the agreement monotonically improves until they are within approximately 10% near 30 kcal/mol collision energy. We believe these differences arise from stronger vibrational adiabaticity in the quantum dynamics, an effect examined before for this system at lower energies. We have also computed QCT OH(v',j') state-resolved cross sections and angular distributions. The QCT state-resolved OH(v') cross sections peak at the same vibrational quantum number as the H2 reagent. The OH rotational distributions are also quite hot and tend to cluster around high rotational quantum numbers. However, the dynamics seem to dictate a cutoff in the energy going into OH rotation indicating an angular momentum constraint. The state-resolved OH distributions were fit to probability functions based on conventional information theory extended to include an energy gap law for product vibrations.     

Quasi-classical trajectory calculations analyzing the role of bending mode excitations of methane in the Cl + CH4 reaction

The effects of the methane torsional (nu(2)), umbrella (nu(4)), and the combination nu(2)+nu(4) bending mode excitations on the reactivity and dynamics of the gas-phase Cl + CH(4) --> HCl + CH(3) reaction were analyzed. Quasi-classical trajectory (QCT) calculations, including corrections to avoid zero-point energy leakage along the trajectories, were used on an analytical potential energy surface previously developed by our group. With respect to the reactivity, we found that excitation of either bending mode independently gave similar increases in the reactivity, while the increase observed upon excitation of both modes was larger than the sum of the effect of exciting them independently. Both results agree with recent experimental measures. With respect to the dynamics (rotovibrational and angular distributions of the products), the two bending modes and their combination gave very similar pictures, reproducing the experimental behavior. The satisfactory agreement obtained with a great variety of experimental data (always qualitatively acceptable and sometimes even quantitatively) of the present QCT study lends confidence to the potential energy surface constructed by our group.     

State-resolved dynamics of the CH(A2Delta) channels from single and multiple photon dissociation of bromoform in the 10-20 eV energy range

Single photon dissociation of bromoform using synchrotron radiation has been investigated by Fourier transform visible fluorescence spectroscopy (FTVIS). The photodissociation of bromoform in the 12-18 eV energy range results in several products, among which are the CH(A2Delta) and CH(B2Sigma) radicals. Vibrational and rotational state distributions of the CH(A2Delta) are determined from their fluorescence spectra. From the threshold photon energy above which emission from the CH(A2Delta) radicals is observed, the most likely process leading to CH(A) formation is CHBr3 --> CH + 3Br rather than CHBr3 --> CH + Br + Br2. The rotational Boltzmann temperatures in the CH(A --> X) emission spectra for v' = 0 and v' = 1 range between 1570 and 3650 K, depending on the excitation photon energy. From the high rotational excitation, the results suggest that the mechanism for the loss of three bromine atoms is most likely sequential. A small negative emission anisotropy of the CH(A) radicals [(Ipar - Iper)/(Ipar + 2Iper) = -0.024 +/- 0.005] is constant across the action spectrum; a small net absorption dipole of CHBr3 in the vacuum ultraviolet is parallel to the 3-fold symmetry axis of the CHBr3 molecule. The state distributions of the CH(A2Delta) radicals from multiphoton dissociation of bromoform using the 266 nm output (three photons) of a femtosecond laser (Boltzmann temperatures: T(v'=0)(rot)= 4250 +/- 300 K; T(v'=1)(rot)= 3100 +/- 550 K) are compared to those from the single photon dissociation results (Boltzmann temperatures: T(v'=0)(rot)= 3650 +/- 150 K; T(v'=1)(rot)= 2400 +/- 200 K) at the same total excitation energy under collision free conditions. The analysis of the CH(A) rotational populations shows hotter rotational populations for the femtosecond experiment, also suggesting sequential dissociation of the bromoform in the femtosecond experiment. The duration of the femtosecond laser pulse is approximately 180 fs, setting a limit on the time scales for the multiple dissociations.     

Quasi-classical trajectory calculations analyzing the reactivity and dynamics of asymmetric stretch mode excitations of methane in the H + CH4 reaction

An exhaustive dynamics study was performed at two collision energies, 1.52 and 2.20 eV, analyzing the effects of the asymmetric (nu3) stretch mode excitation in the reactivity and dynamics of the gas-phase H + CH4 reaction. Quasi-classical trajectory (QCT) calculations, including corrections to avoid zero-point energy leakage along the trajectories, were performed on an analytical potential energy surface previously developed by our group. First, strong coupling between different vibrational modes in the entry channel was observed, indicating that energy can flow between these modes, and therefore that they do not preserve their adiabatic character along the reaction path; i.e., the reaction is nonadiabatic. Second, we found that the reactant vibrational excitation has a significant influence on the vibrational and rotational product distributions. With respect to the vibrational distribution, our results confirm the purely qualitative experimental evidence, although the theoretical results presented here are also quantitative. The rotational distributions are predictive, because no experimental data have been reported. Third, with respect to the reactivity, we found that the nu3 mode excitation by one quantum is more reactive than the ground state by a factor of about 2, independently of the collision energy, and in agreement with the experimental measurement of 3.0 +/- 1.5. Fourth, the state-to-state angular distributions of the products reproduce the experimental behavior at 1.52 eV, where the CH3 products scatter sideways and backward. At 2.20 eV this experimental information is not available, and therefore the results reported here are again predictive. The satisfactory reproduction of a great variety of experimental data by the present QCT study lends confidence to the potential energy surface constructed by our group and to those results whose accuracy cannot be checked by comparison with experiment.     

Vibrational and rotational distributions of the CH(A2Delta) product of the C(2)H + O(3P) reaction studied by fourier transform visible (FTVIS) emission spectroscopy

The C(2)H + O((3)P) --> CH(A) + CO reaction is investigated using Fourier transform visible emission spectroscopy. The O((3)P) and C(2)H radicals are produced by simultaneous 193 nm photolysis of SO(2) and C(2)H(2) precursors, respectively. The nascent vibrational and rotational distributions of the CH(A) product are obtained under time-resolved, but quasi-steady-state, conditions facilitated by the short lifetime of the CH(A) emission. The vibrational temperature of the CH(A) product is found to be appreciably hotter (2800 +/- 100 K) than the rotational distributions in the v' = 0 (1400 +/- 100 K) and v' = 1 (1250 +/- 250 K) levels. The results suggest that the reaction may proceed through an electronically excited HCCO() intermediate; moreover, the vibrational excitation compared to rotational excitation is higher than expected based on a statistical distribution of energy and may be the result of geometrical changes in the transition state. The CH(A) emission is also observed in a C(2)H(2)/O/H reaction mixture using a microwave discharge apparatus to form O atoms, with subsequent H atom production. The nascent rotational and vibrational distributions of the CH(A) determined by the microwave discharge apparatus are very similar to the CH(A) distributions obtained in the photodissociation experiment. The results support the idea that the C(2)H + O((3)P) reaction may play a role in low-pressure C(2)H(2)/O/H flames, as previously concluded.     

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.