State-to-state reaction probabilities for the H+O2(v,j)-->O+OH(v',j') reaction on three potential energy surfaces

We report state-to-state and total reaction probabilities for J=0 and total reaction probabilities for J=2 and 4 for the title reaction, both for ground-state and initially rovibrationally excited reactants. The results for three different potential energy surfaces are compared and contrasted. The potential energy surfaces employed are the DMBE IV surface by Pastrana et al. [J. Phys. Chem. 94, 8073 (1990)], the surface by Troe and Ushakov (TU) [J. Chem. Phys. 115, 3621 (2001)], and the new XXZLG ab initio surface by Xu et al. [J. Chem. Phys. 122, 244305 (2005)]. Our results show that the total reaction probabilities from both the TU and XXZLG surfaces are much smaller in magnitude for collision energies above 1.2 eV compared to the DMBE IV surface. The three surfaces also show different behavior with regards to the effect of initial state excitation. The reactivity is increased on the XXZLG and the TU surfaces and decreased on the DMBE IV surface. Vibrational and rotational product state distributions for the XXZLG and the DMBE IV surface show different behaviors for both types of distributions. Our results show that for energies above 1.25 eV the dynamics on the DMBE IV surface are not statistical. However, there is also evidence that the dynamics on the XXZLG surface are not purely statistical for energies above the onset of the first excited product vibrational state v'=1. The magnitude of the total reaction probability is decreased for J>0 for the DMBE IV and the XXZLG surfaces for ground-state reactants. However, for initially rovibrationally excited reactants, the total reaction probability does not decrease as expected for both surfaces. As a result the total cross section averaged over all Boltzmann accessible rotational states may well be larger than the cross section reported in the literature for j=1.     

State-to-state photodissociation dynamics of triatomic molecules: H2O in the B band

State-to-state photodissociation dynamics of H(2)O in its B band has been investigated quantum mechanically on a new set of non-adiabatically coupled potential energy surfaces for the lowest two (1)A' states of H(2)O, which are developed at the internally contracted multi-reference configuration interaction level with the aug-cc-pVQZ basis set. Quantum dynamical calculations carried out using the Chebyshev propagator yield absorption spectra, product state distributions, branching ratios, and differential cross sections, which are in reasonably good agreement with the latest experimental results. Particular focus is placed here on the dependence of various dynamical observables on the photon energy. Detailed analyses of the dynamics have assigned the diffuse structure in absorption spectrum to short-time recurring dynamics near the HOH conical intersection. The non-adiabatic dissociation to the ground state OH product via the HOH conical intersection is facile, direct, fast, and produces rotationally hot OH(X?) products. On the other hand, the adiabatic channel on the excited state leading to the OH(A?) product is dominated by long-lived resonances, which depend sensitively on the potential energy surfaces.     

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.     

Macroscopic evidences for non-Rice-Ramsperger-Kassel effects in the reaction between H3O+ and D2O: the occurrence of nonstatistical isotopic branching ratio

The dynamics of the isotopic scrambling in the energized and metastable complex D2O-H3O+ has been studied using classical molecular dynamics (MD) trajectories starting from regions of phase space corresponding to an already formed collisional complex. The simulations cover the range of internal energies spanned by gas phase collision experiments. Rate constants for the isotopic exchange and the complex dissociation have been computed; the isotopic branching ratio R=[HD2O+]/[H2DO+] has also been obtained from MD simulations and found to deviate substantially from an equivalent prediction based on a previously proposed kinetic scheme. This finding suggests the possibility that details of the reaction dynamics play a role in defining the isotopic branching ratio. The analysis of trajectory results indicated a relatively long lifetime for the collisional complex and the presence of multiple time scales for the exchange process, with a large fraction of the exchange events being separated only by a single oxygen-oxygen vibration or half of it. The occurrence of these fast consecutive jumps and their different probabilities as a function of the relative direction between first and second jumps suggest the presence of ballistic motion in the complex following each reactive event. This can be explained on the basis of overlapping regions in phase space and it is used to provide an explanation of the difference between kinetic and MD branching ratios.     

State-to-state quantum dynamics of the H + HBr reaction: competition between the abstraction and exchange reactions

Quantum state-to-state dynamics for the H + HBr(υ(i) = 0, j(i) =0) reaction was studied on an accurate ab intio potential energy surface for the electronic ground state of BrH(2). Both the H + HBr → H(2) + Br abstraction reaction and the H' + HBr → H'Br + H exchange reaction were investigated up to a collision energy of 2.0 eV. It was found that the abstraction channel is dominant at lower collision energies, while the exchange channel becomes dominant at higher collision energies. The total integral cross section of the abstraction reaction at a collision energy of 1.6 eV was found to be 1.37 A?(2), which is larger than a recent quantum mechanical result (1.06 A?(2)) and still significantly smaller than the experimental value (3 ± 1 A?(2)). Meanwhile, similar to the previous theoretical study, our calculations also predicted much hotter product rotational state distributions than those from the experimental study. This suggests that further experimental investigations are highly desirable to elucidate the dynamic properties of the title reactions.     

State-to-state dynamics of the hydrogen exchange reaction

Nascent product quantum state distributions have been measured for reactive, H + D₂ → HD + D, and inelastic, H + D₂ → D₂? + H, collisions of H with D₂, at collision energies of 0.98 eV, 1.1 eV, and 1.3 eV. These distributions are extracted from coherent anti-Stokes Raman scattering (CARS) spectra of HD and D₂ (recorded under single-collision conditions) following laser photodissociation of HI to generate translationally hot H atoms in a D₂/HI gas mixture. Variation of the photolysis wavelength allows selection of the H atom translational energy. These rotational and vibrational state distributions are compared with those obtained from quasiclassical trajectory calculations on an ab initio potential energy surface. The agreement between the calculated and measured distributions is extremely good. TheHD and D₂ product quantum state distributions are well represented by simple linear surprisal functions, with large positive vibrational and rotational surprisal.     

State-to-state quantum reactive scattering for four-atom chemical reactions: differential cross section for the H+H2O-->H2+OH abstraction reaction

The time-dependent wave packet method was extended to calculate the state-to-state differential cross section for the title four-atom abstraction reaction with H2O in the ground rovibrational state. One spectator OH bond length was fixed in the study, but the remaining five degrees of freedom were treated exactly. It was found that (a) the differential cross section changes from being strongly backward peaked at low collision energy to sideward scattering at E = 1.4 eV, and (b) the rotational state-resolved differential cross section for H2 differs substantially from that for OH.     

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.     

State-to-state rotational transitions in H2+H2 collisions at low temperatures

We present quantum mechanical close-coupling calculations of collisions between two hydrogen molecules over a wide range of energies, extending from the ultracold limit to the superthermal region. The two most recently published potential energy surfaces for the H(2)-H(2) complex, the so-called Diep-Johnson (DJ) [J. Chem. Phys. 112, 4465 (2000); 113, 3480 (2000)] and Boothroyd-Martin-Keogh-Peterson (BMKP) [J. Chem. Phys. 116, 666 (2002)] surfaces, are quantitatively evaluated and compared through the investigation of rotational transitions in H(2)+H(2) collisions within rigid rotor approximation. The BMKP surface is expected to be an improvement, approaching chemical accuracy, over all conformations of the potential energy surface compared to previous calculations of H(2)-H(2) interaction. We found significant differences in rotational excitation/deexcitation cross sections computed on the two surfaces in collisions between two para-H(2) molecules. The discrepancy persists over a large range of energies from the ultracold regime to thermal energies and occurs for several low-lying initial rotational levels. Good agreement is found with experiment B. Mate et al., [J. Chem. Phys. 122, 064313 (2005)] for the lowest rotational excitation process, but only with the use of the DJ potential. Rate coefficients computed with the BMKP potential are an order of magnitude smaller.     

Direct dynamics study on the hydrogen abstraction reaction CH2O + HO2 --> CHO + H2O2

We present a direct ab initio dynamics study on the hydrogen abstraction reaction CH2O + HO2 --> CHO + H2O2, which is predicted to have four possible reaction channels caused by different attacking orientations of HO2 radical to CH2O. The structures and frequencies at the stationary points and the points along the minimum energy paths (MEPs) of the four reaction channels are calculated at the B3LYP/cc-pVTZ level of theory. Energetic information of stationary points and the points along the MEPs is further refined by means of some single-point multilevel energy calculations (HL). The rate constants of these channels are calculated using the improved canonical variational transition-state theory with the small-curvature tunneling correction (ICVT/SCT) method. The calculated results show that, in the whole temperature range, the more favorable reaction channels are Channels 1 and 3. The total ICVT/SCT rate constants of the four channels at the HL//B3LYP/cc-pVTZ level of theory are in good agreement with the available experiment data over the measured temperature ranges, and the corresponding three-parameter expression is k(ICVT/SCT) = 3.13 x 10(-20) T(2.70) exp(-11.52/RT) cm3 mole(-1) s(-1) in the temperature range of 250-3000 K. Additionally, the flexibility of the dihedral angle of H2O2 is also discussed to explain the different experimental values.     

Product spin-orbit state resolved dynamics of the H+H2O and H+D2O abstraction reactions

The product state-resolved dynamics of the reactions H+H(2)O/D(2)O-->OH/OD((2)Pi(Omega);v',N',f )+H(2)/HD have been explored at center-of-mass collision energies around 1.2, 1.4, and 2.5 eV. The experiments employ pulsed laser photolysis coupled with polarized Doppler-resolved laser induced fluorescence detection of the OH/OD radical products. The populations in the OH spin-orbit states at a collision energy of 1.2 eV have been determined for the H+H(2)O reaction, and for low rotational levels they are shown to deviate from the statistical limit. For the H+D(2)O reaction at the highest collision energy studied the OD((2)Pi(3/2),v'=0,N'=1,A') angular distributions show scattering over a wide range of angles with a preference towards the forward direction. The kinetic energy release distributions obtained at 2.5 eV also indicate that the HD coproducts are born with significantly more internal excitation than at 1.4 eV. The OD((2)Pi(3/2),v'=0,N'=1,A') angular and kinetic energy release distributions are almost identical to those of their spin-orbit excited OD((2)Pi(1/2),v'=0,N'=1,A') counterpart. The data are compared with previous experimental measurements at similar collision energies, and with the results of previously published quasiclassical trajectory and quantum mechanical calculations employing the most recently developed potential energy surface. Product OH/OD spin-orbit effects in the reaction are discussed with reference to simple models.     

A full dimensional, nine-degree-of-freedom, time-dependent quantum dynamics study for the H2+C2H reaction

A full dimensional, nine-degree-of-freedom (9DOF), time-dependent quantum dynamics wave packet approach is presented for the study of the H2+C2H-->H+C2H2 reaction system. This is the first full dimensional quantum dynamics study for a diatom-triatom reaction system. The effects of the initial vibrational and rotational excitations of the reactants on the reactivity of this reaction are investigated. This study shows that vibrational excitations of H2 enhance the reactivity; whereas, the vibrational excitations of C2H only have a small effect on the reaction probability. In addition, the bending excitations of C2H, compared to the ground state reaction probability, hinder the reactivity. Comparison of the ground state reaction probabilities of the 9DOF and 8DOF shows the reaction probability from the full dimensional calculation is larger, with more prominent resonance features.     

Transition state dynamics of the OH + H2O hydrogen exchange reaction studied by dissociative photodetachment of H3O2-

Dynamics in the transition state region of the bimolecular OH + H2O-->H2O + OH hydrogen exchange reaction have been studied by photoelectron-photofragment coincidence spectroscopy of the H3O2- negative ion and its deuterated analog D3O2-. The data reveal vibrationally resolved product translational energy distributions. The total translational energy distribution shows a vibrational progression indicating excitation of the antisymmetric stretch of the water product. Electronic structure calculations at the QCISD level of theory support this analysis. Examination of the translational energy release between the neutral products reveals a dependence on the product vibrational state. These data should provide a critical test of ab initio potential energy surfaces and dynamics calculations.     

State-to-state dynamics analysis of the F + CHD3 reaction: a quasiclassical trajectory study

An exhaustive state-to-state dynamics study was performed to analyze the F + CHD3 --> FD(nu', j') + CHD2(nu) gas-phase abstraction reaction. Quasiclassical trajectory (QCT) calculations, including corrections to avoid zero-point energy leakage along the trajectories, were performed at different collision energies on an analytical potential energy surface (PES-2006) recently developed by our group. Whereas the CHD2 coproduct appears vibrationally and rotationally cold, most of the available energy appears as FD(nu') product vibrational energy, peaking at nu' = 2 and nu' = 3, with the population in the latter level growing as the energy increases. The excitation function rises from the threshold of the reaction and then levels off at higher energies, with the maximum contribution from the FD(nu' = 3) level. The state-specific FD(nu') scattering distributions correlated with the coproduct CHD2 in the nu4 = 2 and nu3 = 1 states, at different collision energies, show a steady change from backward to forward scattering as the energy increases. This similar behavior for the two coproduct vibrational states, nu4 = 2 and nu3 = 1, agrees qualitatively with the experimental measurements. Comparison with theoretical and experimental results for the isotopic analogues, F + CH4 and F + CD4, shows that the title reaction presents a direct mechanism, similar to the perdeuterated reaction, but contrasts with that of the F + CH4 reaction. These results for the dynamics of different isotopic variants, always in qualitative and sometimes in quantitative agreement with experiment, show the capacity of the PES-2006 surface to correctly describe the title reaction, even though there are differences that could be due to deficiencies of the PES but also to the known limitations of the classical treatment in the QCT method.     

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.     

反应N+NH→N2+H的态-态量子动力学研究

采用在MRCI/aug-cc-pVQZ水平上构建的N₂H基态势能面, 并运用Chebyshev实波包法研究了N + NH→N₂ + H反应的量子动力学, 如反应几率、 积分截面以及产物振转态分布等. 在50~500 K温度范围内, 该反应的速率常数随着温度升高而递增, 与基于其它势能面的理论结果吻合. 然而, 在室温条件下, 所有理论计算的速率常数均显著大于实验值.     

Quantum dynamics of the H + O2 --> O + OH reaction on an accurate ab initio potential energy surface

We report exact time-dependent and time-independent quantum mechanical studies of the title reaction on an accurate ab initio potential energy surface of Xu et al. (J. Chem. Phys. 2005, 122, 24305). The J = 0 reaction probabilities for several reactant states show sharp resonance structures superimposed on relatively low backgrounds, and they are remarkably different from existing quantum results on an earlier potential energy surface (DMBE-IV). The new findings reported here suggest that our current understanding of this important reaction might require significant revision.     

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.     

·C2H3+O2→HC·O+H2CO 的密度泛函理论研究

应用密度泛函理论研究了@C2H3+O2→HC@O+H2CO的反应机理.在DFT(B3LYP/6-31G*)水平上对反应过程中所有反应物、中间体、过渡态和产物的几何构型进行优化,通过频率振动分析确认中间体和过渡态.计算IRC反应路径的能量,分析了中间体的异构化过程和各主要原子的自旋密度.