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.     

Competitive C-H and O-D bond fission channels in the UV photodissociation of the deuterated hydroxymethyl radical CH2OD

Photodissociation studies of the CH2OD radical in the region 28,000-41,000 cm(-1) (357-244 nm), which includes excitation to the 3s, 3p(x), and 3p(z) states, are reported. H and D photofragments are monitored by using resonance-enhanced multiphoton ionization (REMPI) from the onset of H formation at approximately 30,500 cm(-1) to the origin band region of the 3pz(2A")<--1 2A" transition at 41,050 cm(-1). Kinetic energy distributions P(ET) and recoil anisotropy parameters as a function of kinetic energy, beta(eff)(ET), are determined by the core sampling technique for the channels producing H and D fragments. Two dissociation channels are identified: (I) D+CH2O and (II) H+CHOD. The contribution of channel II increases monotonically as the excitation energy is increased. Based on the calculations of Hoffmann and Yarkony [J. Chem. Phys. 116, 8300 (2002)], it is concluded that conical intersections between 3s and the ground state determine the final branching ratio even when initial excitation accesses the 3px) and 3pz states. The different beta(eff) values obtained for channels I and II (-0.7 and approximately 0.0, respectively) are attributed to the different extents of out-of-plane nuclear motions in the specific couplings between 3s and the ground state (of A' and A' symmetry, respectively) that lead to each channel. The upper limit to the dissociation energy of the C-H bond, determined from P(ET), is D0(C-H)=3.4+/-0.1 eV (79+/-2 kcal/mol). Combining this value with the known heats of formation of H and CH2OD, the heat of formation of CHOD is estimated at DeltaHf(0)(CHOD)=24+/-2 kcal/mol.     

Inelastic scattering matrix elements for the nonadiabatic collision B(2P1/2)+H2(1Sigmag+,j)<-->B(2P3/2)+H2(1Sigmag+,j')

Inelastic scattering matrix elements for the nonadiabatic collision B(2P1/2)+H2(1Sigmag+,j)<-->B(2P3/2)+H2(1Sigmag+,j') are calculated using the time dependent channel packet method (CPM). The calculation employs 1 2A', 2 2A', and 1 2A" adiabatic electronic potential energy surfaces determined by numerical computation at the multireference configuration-interaction level [M. H. Alexander, J. Chem. Phys. 99, 6041 (1993)]. The 1 2A' and 2 2A', adiabatic electronic potential energy surfaces are transformed to yield diabatic electronic potential energy surfaces that, when combined with the total B+H2 rotational kinetic energy, yield a set of effective potential energy surfaces [M. H. Alexander et al., J. Chem. Phys. 103, 7956 (1995)]. Within the framework of the CPM, the number of effective potential energy surfaces used for the scattering matrix calculation is then determined by the size of the angular momentum basis used as a representation. Twenty basis vectors are employed for these calculations, and the corresponding effective potential energy surfaces are identified in the asymptotic limit by the H2 rotor quantum numbers j=0, 2, 4, 6 and B electronic states 2Pja, ja=1/2, 3/2. Scattering matrix elements are obtained from the Fourier transform of the correlation function between channel packets evolving in time on these effective potential energy surfaces. For these calculations the H2 bond length is constrained to a constant value of req=1.402 a.u. and state to state scattering matrix elements corresponding to a total angular momentum of J=1/2 are discussed for j=0<-->j'=0,2,4 and 2P1/2<-->2P1/2, 2P3/2 over a range of total energy between 0.0 and 0.01 a.u.     

Potential energy surface intersections in the C(1D)H2 reactive system

Potential energy surface (PES) intersection seams of two or more electronic states from the 1 1A', 2 1A', 3 1A', 1 1A", and 2 1A" states in the C(1D)H2 reactive system are investigated using the internally contracted multireference configuration interaction method and the aug-cc-pVQZ basis set. Intersection seams with energies less than 20 kcal/mol relative to the C(1D) + H2 asymptote are searched systematically, and finally several seam lines (at the linear H-C-H, linear C-H-H, and C(2v), geometries, respectively) and a seam surface (at Cs geometries) are discovered and determined. The minimum energy crossing points on these seams are reported and the influences of the PES intersections, in particular, conical intersections, on the CH2 spectroscopy and the C(1D) + H2 reaction dynamics are discussed. In addition, geometries and energies of the 1 1A2 and 1 1B2 states of methylene biradical CH2 are reported in detail for the first time.     

A new reaction path for the C + NO reaction: dynamics on the 4A' potential-energy surface

We present a new reaction path without significant barriers for the C + NO reaction, forming ground state N((4)S) and CO. Electronic structure (CASPT2) calculations have been performed for the two lowest (4)A' states of the CNO system. The lowest of these states shows no significant barriers against reaction in the C + NO or O + CN channels. This surface has been fitted to an analytical function using a many-body expansion. Using this surface, and the previously published (2)A' and (2)A' surfaces [Andersson et al., Phys. Chem. Chem. Phys., 2000, 2, 613; Andersson et al., Chem. Phys., 2000, 259, 99], we have performed quasiclassical trajectory (QCT) calculations, obtaining rate coefficients for the C((3)P) + NO(X(2)Pi) --> CO(X(1)Sigma(+)) + N((4)S,(2)D) and C((3)P) + NO(X(2)Pi) --> O((3)P) + CN(X(2)Sigma(+)) reactions. We have also simulated the crossed molecular beam experiments of Naulin et al. [Chem. Phys., 1991, 153, 519] The inclusion of the (4)A' surface in the QCT calculations gives excellent agreement with experiments. This is the first time an adiabatic pathway from C((3)P) + NO(X(2)Pi) to CO(X(1)Sigma(+))+N((4)S) has been reported.     

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.     

Cl-loss and H-loss dissociations in low-lying electronic states of the CH3Cl+ ion studied using multiconfiguration second-order perturbation theory

To examine the experimentally suggested scheme of the pathways for Cl- and H-loss dissociations of the CH(3)Cl(+) ion in the X(2)E (1(2)A', 1(2)A' '), A(2)A(1) (2(2)A'), and B(2)E (3(2)A', 2(2)A") states, the complete active space-self-consistent field (CASSCF) and multiconfiguration second-order perturbation theory (CASPT2) calculations with an atomic natural orbital (ANO) basis were performed for the 1(2)A' (X(2)A'), 1(2)A", 2(2)A', and 2(2)A'" states. The potential energy curves describing dissociation from the four C(s) states were obtained on the basis of the CASSCF partial geometry optimization calculations at fixed C-Cl or C-H distance values, followed by the CASPT2 energy calculations. The electronic states of the CH3(+) and CH(2)Cl(+) ions produced by Cl-loss and H-loss dissociation, respectively, were carefully determined. Our calculations confirm the following experimental facts: Cl-loss dissociation occurs from the 1(2)A' (X(2)A'), 1(2)A", and 2(2)A' states (all leading to CH3(+) (X(1)A(1)') + Cl), and H-loss dissociation does not occur from 2(2)A'. The calculations indicate that H-loss dissociation occurs from the 1(2)A' and 1(2)A' ' states (leading to CH(2)Cl(+) (X(1)A(1)) + H and CH(2)Cl(+) (1(3)A") + H, respectively). The calculations also indicate that H-loss dissociation occurs (with a barrier) from the 2(2)A" state (leading to CH(2)Cl(+) (1(1)A") + H), supporting the observation of direct dissociation from the B state to CH(2)Cl(+) and that Cl-loss dissociation occurs from the 2(2)A" state (leading to CH3(+) (1(3)A") + Cl), not supporting the previously proposed Cl-loss dissociation of the B state via internal conversion of B to A. The predicted appearance potential values for CH3(+) (X(1)A(1)') and CH(2)Cl(+) (X(1)A(1)) are in good agreement with the experimental values.     

Quasiclassical trajectory calculations comparing the reactivity and dynamics of symmetric and asymmetric stretch and the role of the bending mode excitations of methane in the Cl + CH4 reaction

To analyze the effects of the symmetric (nu(1)) and asymmetric (nu(3)) stretch mode excitations and the role played by the "umbrella" bending (nu(4)) mode excitation in the reactivity and the dynamics of the gas-phase Cl+CH(4) reaction, an exhaustive dynamics study was performed. Quasiclassical trajectory (QCT) calculations, including corrections to avoid zero-point energy leakage along the trajectories, were used in this work on an analytical potential energy surface previously developed by Espinosa-Garcia et al. [J. Chem. Phys. (to be published)]. First, with respect to the reactivity, we found that the nu(1) mode excitation is more reactive than the nu(3) mode by a factor of 1.20, in agreement with the experimental tendency between these modes. The inclusion of the nu(4) bending mode practically does not affect this relative reactivity, (nu(1+)nu(4))(nu(3+)nu(4)) = 1.16. Second, with respect to the dynamics (rotovibrational and angular distributions of the products), the two stretch modes, nu(1) and nu(3), give very similar pictures, reproducing the experimental behavior, and the nu(4) "umbrella" mode does not affect the dynamics. The satisfactory reproduction (always qualitatively acceptable and sometimes even quantitatively) of a great variety of experimental data by the QCT study presented here lends confidence to the potential energy surface constructed by Espinosa-Garcia et al. [J. Chem. Phys. (to be published)].     

Experimental and theoretical differential cross sections for the N(2D) + H2 reaction

In this paper, we report a combined experimental and theoretical study on the dynamics of the N(2D) + H2 insertion reaction at a collision energy of 15.9 kJ mol(-1). Product angular and velocity distributions have been obtained in crossed beam experiments and simulated by using the results of quantum mechanical (QM) scattering calculations on the accurate ab initio potential energy surface (PES) of Pederson et al. (J. Chem. Phys. 1999, 110, 9091). Since the QM calculations indicate that there is a significant coupling between the product angular and translational energy distributions, such a coupling has been explicitly included in the simulation of the experimental results. The very good agreement between experiment and QM calculations sustains the accuracy of the NH2 ab initio ground state PES. We also take the opportunity to compare the accurate QM differential cross sections with those obtained by two approximate methods, namely, the widely used quasiclassical trajectory calculations and a rigorous statistical method based on the coupled-channel theory.     

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".     

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.     

Quasiclassical trajectory study of the reaction H+CH4(nu3 = 0,1)-->CH3+H2 using a new ab initio potential energy surface

Detailed quasiclassical trajectory calculations of the reaction H+CH4(nu3 = 0,1)-->CH3 + H2 using a slightly updated version of a recent ab initio-based CH5 potential energy surface [X. Zhang et al., J. Chem. Phys. 124, 021104 (2006)] are reported. The reaction cross sections are calculated at initial relative translational energies of 1.52, 1.85, and 2.20 eV in order to make direct comparison with experiment. The relative reaction cross section enhancement ratio due to the excitation of the C-H antisymmetric stretch varies from 2.2 to 3.0 over this energy range, in good agreement with the experimental result of 3.0 +/- 1.5 [J. P. Camden et al., J. Chem. Phys. 123, 134301 (2005)]. The laboratory-frame speed and center-of-mass angular distributions of CH3 are calculated as are the vibrational and rotational distributions of H2 and CH3. We confirm that this reaction occurs with a combination of stripping and rebound mechanisms by presenting the impact parameter dependence of these distributions and also by direct examination of trajectories.     

Nonadiabatic quantum dynamics in O(3P)+H2→OH+H: a revisited study

To investigate the extent of nonadiabatic effects in the title reaction, quasi-classical trajectory and nonadiabatic quantum scattering as well as the nonadiabatic quantum-classical trajectory calculations were performed on the accurate ab initio benchmark potential energy surfaces of the lowest (3)A' and (3)A" electronic states [Rogers et al., J Phys Chem A 2000, 104, 2308], together with the spin-orbit coupling matrix [Maiti and Schatz, J Chem Phys 2003, 119, 12360] and the lowest singlet (1) A' potential energy surface [Dobby and Knowles, Faraday Discuss 1998, 110, 247]. Comparison of the calculated total cross sections from both adiabatic and nonadiabatic calculations has demonstrated that for adiabatic channels including (3)A'→(3)A' and (3)A"→(3)A", difference does exist between the two kinds of adiabatic and nonadiabatic calculations, showing nonadiabatic effects to some extent. Such nonadiabatic effects tend to become more conspicuous at high collision energies and are found to be more pronounced with trajectories/quantum wave packet initiated on (3)A' than on (3)A". Furthermore, the present study also showed that nonadiabatic effects can bring the component of forward-scattering in the product angular distributions.     

Cumulative reaction probabilities: A comparison between quasiclassical and quantum mechanical results

This article presents a quasiclassical trajectory (QCT) method for determining the cumulative reaction probability (CRP) as a function of the total energy. The method proposed is based on a discrete sampling using integer values of the total and orbital angular momentum quantum numbers for each trajectory and on the development of equations that have a clear counterpart in the quantum mechanical (QM) case. The calculations comprise cumulative reaction probabilities at a given total angular momentum J, as well as those summed over J. The latter are used to compute QCT rate constants. The method is illustrated by comparing QCT and exact QM results for the H+H2, H+D2, D+H2, and H+HD reactions. The agreement between QCT and QM results is very good, with small discrepancies between the two data sets indicating some genuine quantum effects. The most important of these involves the value of the CRP at low energies which, due to the absence of tunneling, is lower in the QCT calculations, causing the corresponding rate constants to be smaller. The second is the steplike structure that is clearly displayed in the QM CRP for J = 0, which is much smoother in the corresponding QCT results. However, when the QCT density of reactive states, i.e., the derivatives of the QCT CRP with respect to the energy, is calculated, a succession of maxima and minima is obtained which roughly resembles those found in the QM calculations, although the latter are considerably sharper. The analysis of the broad peaks in the QCT density of reactive states indicates that the distributions of collision times associated with the maxima are somewhat broader, with a tail extending to larger collision times, than those associated with the minima. In addition, the QM and QCT dynamics of the isotopic variants mentioned above are compared in the light of their CRPs. Issues such as the compliance of the QCT CRP with the law of microscopic reversibility, as well as the similarity between the CRPs for ortho and para species in the QM and QCT cases, are also addressed.     

Theoretical reinvestigation of the O(3P) + C6H6 reaction: quantum chemical and statistical rate calculations

The lowest-lying triplet and singlet potential energy surfaces for the O(3P) + C6H6 reaction were theoretically characterized using the "complete basis set" CBS-QB3 model chemistry. The primary product distributions for the multistate multiwell reactions on the individual surfaces were then determined by RRKM statistical rate theory and weak-collision master equation analysis using the exact stochastic simulation method. It is newly found that electrophilic O-addition onto a carbon atom in benzene can occur in parallel on two triplet surfaces, 3A' and 3A' '; the results predict O-addition to be dominant up to combustion temperatures. Major expected end-products of the addition routes include phenoxy radical + H*, phenol and/or benzene oxide/oxepin, in agreement with the experimental evidence. While c-C6H5O* + H* are nearly exclusively formed via a spin-conservation mechanism on the lowest-lying triplet surface, phenol and/or benzene oxide/oxepin are mainly generated from the lowest-lying singlet surface after inter-system crossing from the initial triplet surface. CO + c-C5H6 are predicted to be minor products in flame conditions, with a yield < or = 5%. The O + C6H6 --> c-C5H5* + *CHO channel is found to be unimportant under all relevant combustion conditions, in contrast with previous theoretical conclusions (J. Phys. Chem. A 2001, 105, 4316). Efficient H-abstraction pathways are newly identified, occurring on two different electronic state surfaces, 3B1 and 3B2, resulting in hydroxyl plus phenyl radicals; they are predicted to play an important role at higher temperatures in hydrocarbon combustion, with estimated contributions of ca. 50% at 2000 K. The overall thermal rate coefficient k(O + C6H6) at 300-800 K was computed using multistate transition state theory: k(T) = 3.7 x 10-16 x T 1.66 x exp(-1830 K/T) cm(3) molecule(-1) s(-1), in good agreement with the experimental data available.     

Collision energy dependence of the O(1D) + HCl --> OH + Cl(2P) reaction studied by crossed beam scattering and quasiclassical trajectory calculations on ab initio potential energy surfaces

The dynamics of the O(1D) + HCl --> OH + Cl(2P) reaction are investigated by a crossed molecular beam ion-imaging method and quasiclassical trajectory calculations on the three ab initio potential energy surfaces, the ground 1(1)A' and two excited (1(1)A' and 2(1)A') states. The scattering experiment was carried out at collision energies of 4.2, 4.5, and 6.4 kcal/mol. The observed doubly differential cross sections (DCSs) for the Cl(2P) product exhibit almost no collision energy dependence over this inspected energy range. The nearly forward-backward symmetric DCS indicates that the reaction proceeds predominantly on the ground-state potential energy surface at these energies. Variation of the forward-backward asymmetry with collision energy is interpreted using an osculating complex model. Although the potential energy surfaces obtained by CASSCF-MRCI ab initio calculations exhibit relatively low potential barriers of 1.6 and 6.5 kcal/mol for 1(1)A' and 2(1)A', respectively, the dynamics calculations indicate that contributions of these excited states are small at the collision energies lower than 15.0 kcal/mol. Theoretical DCSs calculated for the ground-state reaction pathway agree well with the observed ones. These experimental and theoretical results suggest that the titled reaction at collision energies less than 6.5 kcal/mol is predominantly via the ground electronic state.     

F-loss and H-loss dissociations in low-lying electronic states of the CH3F+ ion studied using multiconfiguration second-order perturbation theory

Complete active space self-consistent field (CASSCF) and multiconfiguration second-order perturbation theory (CASPT2) calculations with an atomic natural orbital basis were performed for the 1(2)A', 1(2)A', 2(2)A', 2(2)A', and 3(2)A' (X2E, A2A1, and B2E) states of the CH3F+ ion. The 1(2)A' state is predicted to be the ground state, and the C(s)-state energy levels are different from those of the CH3Cl+ ion. The 2(2)A' (A2A1) state is predicted to be repulsive, and the calculated adiabatic excitation energies for 2(2)A' and 3(2)A' are very close to the experimental value for the B state. The CASPT2//CASSCF potential energy curves (PECs) were calculated for F-loss dissociation from the five C(s) states and H-loss dissociation from the 1(2)A', 1(2)A', and 2(2)A' states. The electronic states of the CH3+ and CH2F+ ions as the dissociation products were carefully determined by checking the energies and geometries of the asymptote products, and appearance potentials for the two ions in different states are predicted. The F-loss PEC calculations for CH3F+ indicate that F-loss dissociation occurs from the 1(2)A', 1(2)A', and 2(2)A' states [all correlating with CH3+(X1A1')], which supports the experimental observations of direct dissociation from the X and A states, and that direct F-loss dissociation can occur from the two Jahn-Teller component states of B2E, 2(2)A' and 3(2)A' [correlating with CH3+(1(3)A') and CH3+(1(3)A'), respectively]. Some aspects of the 3(2)A' Cl-loss PEC of the CH3Cl+ ion are inferred on the basis of the calculation results for CH3F+. The H-loss PEC calculations for CH3F+ indicate that H-loss dissociation occurs from the 1(2)A', 1(2)A', and 2(2)A' states [correlating with CH2F+(1(3)A'), CH2F+(X1A1), and CH2F+(1(1)A'), respectively], which supports the observations of direct dissociation from the X and B states. As the 2(2)A' H-loss PEC of CH3Cl+, the 2(2)A' H-loss PEC of CH3F+ does not lead to H + CH2X+, but the PECs of the two ions represent different types of reactions.     

Photodissociation dynamics of the ethoxy radical investigated by photofragment coincidence imaging

The photodissociation dynamics of the ethoxy radical (CH3CH2O) have been studied at energies from 5.17 to 5.96 eV using photofragment coincidence imaging. The upper state of the electronic transition excited at these energies is assigned to the C2A'state on the basis of electronic structure calculations. Fragment mass distributions show two photodissociation channels, OH + C2H4 and CH3 + CH2O. The presence of an additional photodissociation channel, identified as D + C2D4O, is revealed in time-of-flight distributions from the photodissociation of CD3CD2O. The product branching ratios and fragment translational energy distributions for all of the observed mass channels are nonstatistical. Moreover, the significant yield of OH + C2H4 product suggests that the mechanism for this channel involves isomerization on the excited-state surface. Photodissociation at a much lower yield is seen following excitation at 3.91 eV, corresponding to a vibronic band of the B2A' <-- X2A' transition.     

Dynamics of the D+ + H2 and H+ + D2 reactions: a detailed comparison between theory and experiment

An extensive set of experimental measurements on the dynamics of the H(+) + D(2) and D(+) + H(2) ion-molecule reactions is compared with the results of quantum mechanical (QM), quasiclassical trajectory (QCT), and statistical quasiclassical trajectory (SQCT) calculations. The dynamical observables considered include specific rate coefficients as a function of the translational energy, E(T), thermal rate coefficients in the 100-500 K temperature range. In addition, kinetic energy spectra (KES) of the D(+) ions reactively scattered in H(+) + D(2) collisions are also presented for translational energies between 0.4 eV and 2.0 eV. For the two reactions, the best global agreement between experiment and theory over the whole energy range corresponds to the QCT calculations using a gaussian binning (GB) procedure, which gives more weight to trajectories whose product vibrational action is closer to the actual integer QM values. The QM calculations also perform well, although somewhat worse over the more limited range of translational energies where they are available (E(T) < 0.6 eV and E(T) < 0.2 eV for the H(+) + D(2) and D(+) + H(2) reactions, respectively). The worst agreement is obtained with the SQCT method, which is only adequate for low translational energies. The comparison between theory and experiment also suggests that the most reliable rate coefficient measurements are those obtained with the merged beams technique. It is worth noting that none of the theoretical approaches can account satisfactorily for the experimental specific rate coefficients of H(+) + D(2) for E(T)≤ 0.2 eV although there is a considerable scatter in the existing measurements. On the whole, the best agreement with the experimental laboratory KES is obtained with the simulations carried out using the state resolved differential cross sections (DCSs) calculated with the QCT-GB method, which seems to account for most of the observed features. In contrast, the simulations with the SQCT data predict kinetic energy spectra (KES) considerably cooler than those experimentally determined.