Mapping Reaction Dynamics via State-to-State Measurements: Rotations Tell the Tale

The hydrogen abstraction reaction H + C2H6 → H2(v,j) + C2H5. Part II. Theoretical kinetics and dynamics study

Two important issues motivated the present study: the role of the tunnelling contribution at low temperatures and the role of the alkyl fragment in the dynamics.

State-to-state reaction dynamics: A selective review

A selective review of state-to-state reaction dynamics experiments is presented. The review focuses on three classes of reactions that exemplify the rich history and illustrate the current state of the art in such work. These three reactions are (1) the hydrogen exchange reaction, H+H2-->H2+H and its isotopomers; (2) the H+RH-->H2+R reactions, where RH is an alkane, beginning with H+CH4-->H2+CH3 and extending to much larger alkanes; and (3) the Cl+RH-->HCl+R reactions, principally Cl+CH4-->HCl+CH3. We describe the experiments, discuss their results, present comparisons with theory, and introduce heuristic models.

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.

The roaming atom pathway in formaldehyde decomposition

We present a detailed experimental and theoretical investigation of formaldehyde photodissociation to H(2) and CO following excitation to the 2(1)4(1) and 2(1)4(3) transitions in S(1). The CO velocity distributions were obtained using dc slice imaging of single CO rotational states (v=0, j(CO)=5-45). These high-resolution measurements reveal the correlated internal state distribution in the H(2) cofragments. The results show that rotationally hot CO (j(CO) approximately 45) is produced in conjunction with vibrationally "cold" H(2) fragments (v=0-5): these products are formed through the well-known skewed transition state and described in detail in the accompanying paper. After excitation of formaldehyde above the threshold for the radical channel (H(2)CO-->H+HCO) we also find formation of rotationally cold CO (j(CO)=5-28) correlated to highly vibrationally excited H(2) (v=6-8). These products are formed through a novel mechanism that involves near dissociation followed by intramolecular H abstraction [D. Townsend et al., Science 306, 1158 (2004)], and that avoids the region of the transition state entirely. The dynamics of this "roaming" mechanism are the focus of this paper. The correlations between the vibrational states of H(2) and rotational states of CO formed following excitation on the 2(1)4(3) transition allow us to determine the relative contribution to molecular products from the roaming atom channel versus the conventional molecular channel.

Reaction Products with Internal Energy beyond the Kinematic Limit Result from Trajectories Far from the Minimum Energy Path: An Example from H + HBr → H2 + Br

The importance of reactive trajectories straying far from the minimum energy path is demonstrated for the bimolecular reaction H + HBr --> H2(v', j') + Br at 53 kcal/mol collision energy. Product quantum state distributions are measured and calculated using the quasi-classical trajectory technique, and the calculations indicate that highly internally excited H2 products result from indirect reactive trajectories with bent transition states. A general argument is made suggesting that reaction products with internal energy exceeding a kinematic constraint can, in general, be attributed to reactive collisions straying far from the minimum energy path.

Imaging the Reaction Dynamics of OH + CD4. 2. Translational Energy Dependencies

The reaction of OH + CD4 is investigated in a crossed-beam experiment over the collisional energies ranging from reaction threshold of about 5 to 16 kcal/mol. Exploiting a time-sliced ion velocity imaging detection scheme, the coincident information on the two polyatomic product pairs, HOD and CD3, is revealed in a state-correlated manner. The recently discovered vibrational mode-correlation between the two products is found to persist over the full range of collision energies of this study. In addition, the energy dependencies of the correlated cross section, state distribution, and angular distribution are elucidated, providing an unprecedented insight into this important reaction.

Imaging photon-initiated reactions: A study of the Cl(P3∕22)+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.

Quasiclassical trajectory calculations to evaluate a kinematic constraint on internal energy in suprathreshold collision energy abstraction reactions

Experimentally observed product quantum state distributions across a wide range of abstraction reactions at suprathreshold collision energies have shown a strong bias against product internal energy. Only a fraction, sometimes quite a small fraction, of the energetically accessible product quantum states are populated. Picconatto et al. [J. Chem. Phys. 114, 1663 (2001)] noted a simple mathematical relationship between the highest-energy rovibrational states observed and the kinematics of the reaction system. They proposed a reaction model based on reaction kinematics that quantitatively explains this behavior. The model is in excellent agreement with measured quantum state distributions. The assumptions of the model invoke detailed characteristics of reactive trajectories at suprathreshold collision energies. Here we test those assumptions using quasiclassical trajectory calculations for the abstraction reactions H+HCl-->H2+Cl, D+HCl-->HD+Cl, and H+DCl-->HD+Cl. Trajectories were run on a potential-energy surface calculated with a London-Eyring-Polyani-Sato function with a localized 3-center term (LEPS-3C) previously shown to accurately reproduce experimentally observed product state distributions for the H+HCl abstraction reaction. The trajectories sample collision energies near threshold and also substantially above it. Although the trajectories demonstrate some aspects of the model, they show that it is not valid. However, the inadequacy of the proposed model does not invalidate the apparent kinematic basis of the observed energy constraint. The present results show that there must be some other molecular behavior rooted in the reaction kinematics that is the explanation and the source of the constraint.

Global potential energy surfaces for the Al+(1S)+H2 system

Global, three-dimensional multireference ab initio potential energy surfaces have been calculated for the AlH2+ system for the two lowest energy singlet states and the lowest energy triplet state. These surfaces were calculated using the multireference configuration interaction level of theory with a large basis set. The accuracy of the surfaces were checked against available experimental data and previous theoretical investigations. The areas of surface crossings between the ground state singlet surface and the lowest energy triplet surface and the first excited singlet surface have been thoroughly investigated in all three dimensions and found to give rise to two regions of surface crossings--an "early" crossing (reduced H2 distance) and a "late" crossing (enlarged H2 distance). It is anticipated that both of these crossings will be important in modeling the dynamics of the system. Each of the global potential energy surfaces were fit by interpolation methodology to obtain analytic representations of the surfaces. A representative classical simulation on the ground state singlet surface was performed and discussed.

Disagreement between theory and experiment in the simplest chemical reaction: Collision energy dependent rotational distributions for H+D2→HD(ν′=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.

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.