Theoretical study of the dynamics of the H+CH4 and H+C2H6 reactions using a specific-reaction-parameter semiempirical Hamiltonian

We present a theoretical study of the reactions of hydrogen atoms with methane and ethane molecules and isotopomers. High-accuracy electronic-structure calculations have been carried out to characterize representative regions of the potential-energy surface (PES) of various reaction pathways, including H abstraction and H exchange. These ab initio calculations have been subsequently employed to derive an improved set of parameters for the modified symmetrically-orthogonalized intermediate neglect of differential overlap (MSINDO) semiempirical Hamiltonian, which are specific to the H+alkane family of reactions. The specific-reaction-parameter (SRP) Hamiltonian has then been used to perform a quasiclassical-trajectory study of both the H+CH4 and H+C2H6 reactions. The calculated values of dynamics properties of the H+CH4-->H2+CH3 reaction and isotopologues, including alkyl product speed distributions, diatomic product internal-state distributions, and cross sections, are generally in good agreement with experiment and with the results provided by the ZBB3 PES [Z. Xie et al., J. Chem. Phys. 125, 133120 (2006)]. The results of trajectories propagated with the SRP Hamiltonian for the H+C2H6-->H2+C2H5 reaction also agree with experiment. The level of agreement between the results calculated with the SRP Hamiltonian and experiment in both the H+methane and H+ethane reactions indicates that semiempirical Hamiltonians can be improved for not only a specific reaction but also a family of reactions.     

Quasiclassical trajectory study of the CH3++HD-->CH2D++H2 Reaction

A full dimensional ab initio potential energy surface for the CH5+ system based on coupled cluster electronic structure calculations and capable of describing the dissociation of methonium ion into methyl cation and molecular hydrogen (J. Phys. Chem. A 2006, 110, 1569) is used in quasiclassical trajectory calculations of the reaction CH3++HD-->CH2D++H2 for low collision energies of relevance to astrochemistry. Cross sections for the exchange are obtained at several relative translational energies and a fit to the energy dependence of the cross sections is used to obtain the rate constant at temperatures between 10 and 50 K. The calculated rate constant at 10 K agrees well with the previously reported experimental value. Internal energy distributions of the products are presented and discussed in the context of zero-point energy "noncompliance".     

Quantum mechanical and quasiclassical trajectory scattering calculations for the C(1D) + H2 reaction on the second excited 1 1A" potential energy surface

Time-independent quantum mechanical (QM) and quasiclassical trajectory (QCT) scattering calculations have been carried out for the C(1D) + H2 --> CH + H reaction at a collision energy of 80 meV on a newly developed ab initio potential energy surface [B. Bussery-Honvault et al., Phys. Chem. Chem. Phys. 7, 1476 (2005)] of 1 1A" symmetry, corresponding to the second singlet state 1 1B1 of CH2. A general good agreement has been found between the QM and QCT rotational distributions and differential cross sections (DCSs). In both cases, DCSs are strongly peaked in the forward direction with a small contribution in the backward direction in contrast with those obtained on the 1 1A' surface, which are nearly symmetric. Rotational distributions obtained on the 1 1A" surface are somewhat colder than those calculated on the 1 1A' surface. The specific dynamics and the contribution of the 1 1A" surface to the overall reactivity of this system are discussed.     

Quantum dynamics study of H+NH3-->H2+NH2 reaction

We report in this paper a quantum dynamics study for the reaction H+NH3-->NH2+H2 on the potential energy surface of Corchado and Espinosa-Garcia [J. Chem. Phys. 106, 4013 (1997)]. The quantum dynamics calculation employs the semirigid vibrating rotor target model [J. Z. H. Zhang, J. Chem. Phys. 111, 3929 (1999)] and time-dependent wave packet method to propagate the wave function. Initial state-specific reaction probabilities are obtained, and an energy correction scheme is employed to account for zero point energy changes for the neglected degrees of freedom in the dynamics treatment. Tunneling effect is observed in the energy dependency of reaction probability, similar to those found in H+CH4 reaction. The influence of rovibrational excitation on reaction probability and stereodynamical effect are investigated. Reaction rate constants from the initial ground state are calculated and are compared to those from the transition state theory and experimental measurement.     

State-to-state reactive differential cross sections for the H+H2-->H2+H reaction on five different potential energy surfaces employing a new quantum wavepacket computer code: DIFFREALWAVE

State-to-state differential cross sections have been calculated for the hydrogen exchange reaction, H+H2-->H2+H, using five different high quality potential energy surfaces with the objective of examining the sensitivity of these detailed cross sections to the underlying potential energy surfaces. The calculations were performed using a new parallel computer code, DIFFREALWAVE. The code is based on the real wavepacket approach of Gray and Balint-Kurti [J. Chem. Phys. 108, 950 (1998)]. The calculations are parallelized over the helicity quantum number Omega' (i.e., the quantum number for the body-fixed z component of the total angular momentum) and wavepackets for each J,Omega' set are assigned to different processors, similar in spirit to the Coriolis-coupled processors approach of Goldfield and Gray [Comput. Phys. Commun. 84, 1 (1996)]. Calculations for J=0-24 have been performed to obtain converged state-to-state differential cross sections in the energy range from 0.4 to 1.2 eV. The calculations employ five different potential energy surfaces, the BKMP2 surface and a hierarchical family of four new ab initio surfaces [S. L. Mielke, et al., J. Chem. Phys. 116, 4142 (2002)]. This family of four surfaces has been calculated using three different hierarchical sets of basis functions and also an extrapolation to the complete basis set limit, the so called CCI surface. The CCI surface is the most accurate surface for the H3 system reported to date. Our calculations of differential cross sections are the first to be reported for the A2, A3, A4, and CCI surfaces. They show that there are some small differences in the cross sections obtained from the five different surfaces, particularly at higher energies. The calculations also show that the BKMP2 performs well and gives cross sections in very good agreement with the results from the CCI surface, displaying only small divergences at higher energies.     

Study of the H+O2 reaction by means of quantum mechanical and statistical approaches: the dynamics on two different potential energy surfaces

The possible existence of a complex-forming pathway for the H+O(2) reaction has been investigated by means of both quantum mechanical and statistical techniques. Reaction probabilities, integral cross sections, and differential cross sections have been obtained with a statistical quantum method and the mean potential phase space theory. The statistical predictions are compared to exact results calculated by means of time dependent wave packet methods and a previously reported time independent exact quantum mechanical approach using the double many-body expansion (DMBE IV) potential energy surface (PES) [Pastrana et al., J. Phys. Chem. 94, 8073 (1990)] and the recently developed surface (denoted XXZLG) by Xu et al. [J. Chem. Phys. 122, 244305 (2005)]. The statistical approaches are found to reproduce only some of the exact total reaction probabilities for low total angular momenta obtained with the DMBE IV PES and some of the cross sections calculated at energy values close to the reaction threshold for the XXZLG surface. Serious discrepancies with the exact integral cross sections at higher energy put into question the possible statistical nature of the title reaction. However, at a collision energy of 1.6 eV, statistical rotationally resolved cross sections managed to reproduce the experimental cross sections for the H+O(2)(v=0,j=1)-->OH(v(')=1,j('))+O process reasonably well.     

Disagreement between theory and experiment in the simplest chemical reaction: collision energy dependent rotational distributions for H + D2 --> HD(nu' = 3,j') + D

We present experimental rotational distributions for the reaction H + D2 --> HD(nu' = 3,j') + D at eight different collision energies between 1.49 and 1.85 eV. We combine a previous measurement of the state-resolved excitation function for this reaction [Ayers et al., J. Chem. Phys. 119, 4662 (2003)] with the current data to produce a map of the relative reactive cross section as a function of both collision energy and rotational quantum number (an E-j' plot). To compare with the experimental data, we also present E-j' plots resulting from both time-dependent and time-independent quantum mechanical calculations carried out on the BKMP2 surface. The two calculations agree well with each other, but they produce rotational distributions significantly colder than the experiment, with the difference being more pronounced at higher collision energies. Disagreement between theory and experiment might be regarded as surprising considering the simplicity of this system; potential causes of this discrepancy are discussed.     

Real wave packet and quasiclassical trajectory studies of the H+ + LiH reaction

Time-dependent real wave packet (RWP) and quasiclassical trajectory (QCT) calculations have been carried out to study the H(+) + LiH reaction on the ab initio potential-energy surface of Martinazzo et al. [J. Chem. Phys., 2003, 119, 11241]. Total initial state-selected and final state-resolved reaction probabilities for the two possible reaction channels, H(2)(+) + Li and LiH + H(+), have been calculated for total angular momentum J=0 at a broad range of collision energies. Integral cross sections and thermal rate coefficients have been calculated using the QCT method and from the corresponding J=0 RWP reaction probabilities by means of a capture model. The calculated thermal rate coefficients are found to be nearly independent of temperature in the 100-500 K interval with a value of approximately 10(-9) cm(3) s(-1), which is in good agreement with estimates used in evolutionary models of early-Universe lithium chemistry. The RWP results are found to be in good agreement overall with the corresponding QCT calculations.     

The role of the excited electronic states in the C+ + H2O reaction

The electronic excited states of the [COH2]+ system have been studied in order to establish their role in the dynamics of the C+ + H2O-->[COH]+ +H reaction, which is a prototypical ion-molecule reaction. The most relevant minima and saddle points of the lowest excited state have been determined and energy profiles for the lowest excited doublet and quartet electronic states have been computed along the fragmentation and isomerization coordinates. Also, nonadiabatic coupling strengths between the ground and the first excited state have been computed where they can be large. Our analysis suggests that the first excited state could play an important role in the generation of the formyl isomer, which has been detected in crossed beam experiments [D. M. Sonnenfroh et al., J. Chem. Phys. 83, 3985 (1985)], but could not be explained in quasiclassical trajectory computations [Y. Ishikawa et al., Chem. Phys. Lett. 370, 490 (2003); J. R. Flores, J. Chem. Phys. 125, 164309 (2006)].     

Shock tube investigation of CH3 + CH3OCH3

The title reaction has been investigated in a diaphragmless shock tube by laser schlieren densitometry over the temperature range 1163-1629 K and pressures of 60, 120, and 240 Torr. Methyl radicals were produced by dissociation of 2,3-butanedione in the presence of an excess of dimethyl ether. Rate coefficients for CH(3) + CH(3)OCH(3) were obtained from simulations of the experimental data yielding the following expression which is valid over the range 1100-1700 K: k = (10.19 ± 3.0)T(3.78)?exp((-4878/T)) cm(3) mol(-1)s(-1). The experimental results are in good agreement with estimates by Curran and co-workers [Fischer, S. L.; Dryer, F. L.; Curran, H. J. Int. J. Chem. Kinet.2000, 32 (12), 713-740. Curran, H. J.; Fischer, S. L.; Dryer, F. L. Int. J. Chem. Kinet.2000, 32 (12), 741-759] but about a factor of 2.6 lower than those of Zhao et al. [Zhao, Z.; Chaos, M.; Kazakov, A.; Dryer, F. L. Int. J. Chem. Kinet.2008, 40 (1), 1-18].     

Quasi-classical trajectory study of the dynamics of the Cl + CH4→ HCl + CH3 reaction

We present an on-the-fly classical trajectory study of the Cl + CH(4)→ HCl + CH(3) reaction using a specific reaction parameter (SRP) AM1 Hamiltonian that was previously optimized for the Cl + ethane reaction [S. J. Greaves et al., J. Phys Chem A, 2008, 112, 9387]. The SRP-AM1 Hamiltonian is shown to be a good model for the potential energy surface of the title reaction. Calculated differential cross sections, obtained from trajectories propagated with the SRP-AM1 Hamiltonian compare favourably with experimental results for this system. Analysis of the vibrational modes of the methyl radical shows different scattering distributions for ground and vibrationally excited products.     

Ab initio investigations of the radical-radical reaction of O(3P) + C3H3

We present ab initio calculations of the reaction of ground-state atomic oxygen [O((3)P)] with a propargyl (C(3)H(3)) radical based on the application of the density-functional method and the complete basis-set model. It has been predicted that the barrierless addition of O((3)P) to C(3)H(3) on the lowest doublet potential-energy surface produces several energy-rich intermediates, which undergo subsequent isomerization and decomposition steps to generate various exothermic reaction products: C(2)H(3)+CO, C(3)H(2)O+H, C(3)H(2)+OH, C(2)H(2)+CHO, C(2)H(2)O+CH, C(2)HO+CH(2), and CH(2)O+C(2)H. The respective reaction pathways are examined extensively with the aid of statistical Rice-Ramsperger-Kassel-Marcus calculations, suggesting that the primary reaction channel is the formation of propynal (CHCCHO)+H. For the minor C(3)H(2)+OH channel, which has been reported in recent gas-phase crossed-beam experiments [H. Lee et al., J. Chem. Phys. 119, 9337 (2003); 120, 2215 (2004)], a comparison on the basis of prior statistical calculations is made with the nascent rotational state distributions of the OH products to elucidate the mechanistic and dynamic characteristics at the molecular level.     

Trajectory dynamics study of the Ar + CH4 dissociation reaction at high temperatures: the importance of zero-point-energy effects

Large-scale classical trajectory calculations have been performed to study the reaction Ar + CH4--> CH3 +H + Ar in the temperature range 2500 < or = T/K < or = 4500. The potential energy surface used for ArCH4 is the sum of the nonbonding pairwise potentials of Hase and collaborators (J. Chem. Phys. 2001, 114, 535) that models the intermolecular interaction and the CH4 intramolecular potential of Duchovic et al. (J. Phys. Chem. 1984, 88, 1339), which has been modified to account for the H-H repulsion at small bending angles. The thermal rate coefficient has been calculated, and the zero-point energy (ZPE) of the CH3 product molecule has been taken into account in the analysis of the results; also, two approaches have been applied for discarding predissociative trajectories. In both cases, good agreement is observed between the experimental and trajectory results after imposing the ZPE of CH3. The energy-transfer parameters have also been obtained from trajectory calculations and compared with available values estimated from experiment using the master equation formalism; in general, the agreement is good.     

Vibrationally inelastic collisions in H+ +CO system: comparing quantum calculations with experiments

State-resolved cross beam experiments [H. Udseth et al., J. Chem. Phys. 60, 3051 (1974); J. Krutein and F. Linder, J. Chem. Phys. 71, 599 (1979); G. Niedner-Schatteburg and J. P. Toennies, Adv. Chem. Phys. LXXXII, 553 (1992)], coupled with proton energy loss spectroscopy for the inelastic scattering of H(+) from CO in the collision range of 10-30 eV show very low vibrational excitation of the target molecule. Stimulated by the experimentally observed low vibrational inelasticity in the system the ground and the first two low-lying excited electronic potential-energy surfaces have been computed using the ab initio multireference configuration interaction method. Quantum dynamics has been performed on the ground potential energy surface in the framework of vibrational close-coupling rotational infinite-order sudden approximation. The various computed dynamical attributes such as differential and integral cross sections, and average vibrational energy transfer are analyzed in detail, and compared successfully with the available experimental results.     

Cross sections of the O+ + H2 --> OH+ + H ion-molecule reaction and isotopic variants (D2, HD): quasiclassical trajectory study and comparison with experiments

A dynamics study [cross section and microscopic mechanism versus collision energy (E(T))] of the reaction O+ + H2 --> OH+ + H, which plays an important role in Earth's ionosphere and interstellar chemistry, was conducted using the quasiclassical trajectory method, employing an analytical potential energy surface (PES) recently derived by our group [R. Martinez et al., J. Chem. Phys. 120, 4705 (2004)]. Experimental excitation functions for the title reaction, as well as its isotopic variants with D2 and HD, were near-quantitatively reproduced in the calculations in the very broad collision energy range explored (E(T) = 0.01-6.0 eV). Intramolecular and intermolecular isotopic effects were also examined, yielding data in good agreement with experimental results. The reaction occurs via two microscopic mechanisms (direct and nondirect abstraction). The results were satisfactorily interpreted based on the reaction probability and the maximum impact parameter dependences with E(T), and considering the influence of the collinear [OHH]+ absolute minimum of the PES on the evolution from reactants to products. The agreement between theory and experiment suggests that the reaction mainly occurs through the lowest energy PES and nonadiabatic processes are not very important in the wide collision energy range analyzed. Hence, the PES used to describe this reaction is suitable for both kinetics and dynamics studies.     

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.     

Photodissociation of CH3I: a full-dimensional (9D) quantum dynamics study

The photodissociation of methyl iodide in the A band is studied by full-dimensional (9D) wave packet dynamics calculations using the multiconfigurational time-dependent Hartree approach. The potential energy surfaces employed are based on the diabatic potentials of Xie et al. [J. Phys. Chem. A 2000, 104, 1009] and the vertical excitation energy is taken from recent ab initio calculations [Alekseyev et al. J. Chem. Phys.2007, 126, 234102]. The absorption spectrum calculated for exclusively parallel excitation agrees well with the experimental spectrum of the A band. The electronic population dynamics is found to be strongly dependent on the motion in the torsional coordinate related to the H(3)-C-I bend, which presumably is an artifact of the diabatic model employed. The calculated fully product state-selected partial spectra can be interpreted based on the reflection principle and suggests strong coupling between the C-I stretching and the H(3)-C-I bending motions during the dissociation process. The computed rotational and vibrational product distributions typically reproduce the trends seen in the experiment. In agreement with experiment, a small but significant excitation of the total symmetric stretching and the asymmetric bending modes of the methyl fragment can be seen. In contrast, the umbrella mode of the methyl is found to be too highly excited in the calculated distributions.     

Trajectory dynamics study of collision-induced dissociation of the Ar + CH4 reaction at hyperthermal conditions: vibrational excitation and isotope substitution

We investigate the role of vibrational energy excitation of methane and two deuterated species (CD(4) and CH(2)D(2)) in the collision-induced dissociation (CID) process with argon at hyperthermal energies. The quasi-classical trajectory method has been applied, and the reactive Ar + CH(4) system has been modeled by using a modified version of the CH(4) potential energy surface of Duchovic et al. (J. Phys. Chem. 1984, 88, 1339) and the Ar-CH(4) intermolecular potential function obtained by Troya (J. Phys. Chem. A 2005, 109, 5814). This study clearly shows that CID is markedly enhanced with vibrational excitation and, to a lesser degree, with collision energy. In general, CID increases by exciting stretch vibrational modes of the reactant molecule. For the direct dissociation of CH(4), however, the CID cross sections appear to be essentially independent of which vibrational mode is initially excited. In all situations studied, the CID cross sections are always greater for the Ar + CD(4) reaction than for the Ar + CH(4) one, the Ar + CH(2)D(2) being an intermediate situation. A detailed analysis of the energy transfer processes, including their relation with CID, is also presented.     

Isotope branching and tunneling in O(3P)+HD-->OH+D; OD+H reactions

The O((3)P)+HD and O((3)P)+D(2) reactions are studied using quantum scattering calculations and chemically accurate potential energy surfaces developed for the O((3)P)+H(2) system by Rogers et al. [J. Phys. Chem. A 104, 2308 (2000)]. Cross sections and rate coefficients for OH and OD products are calculated using accurate quantum methods as well as the J-shifting approximation. The J-shifting approach is found to work remarkably well for both O+HD and O+D(2) collisions. The reactions are dominated by tunneling at low temperatures and for the O+HD reaction the hydrogen atom transfer leading to the OH product dominates at low temperatures. Our result for the OH/OD branching ratio is in close agreement with previous calculations over a wide range of temperatures. The computed OH/OD branching ratios are also in close agreement with experimental results of Robie et al. [Chem. Phys. Lett. 134, 579 (1987)] at temperatures above 400 K but the theoretical results do not reproduce the rapid rise in the experimental values of the branching ratio for temperatures lower than 350 K. We believe that new measurements could resolve the long-standing discrepancy between experiment and theory for this benchmark reaction.