Potential surfaces and dynamics of the O(3P)+H2O(X1A1)-->OH(X2pi)+OH(X2pi) reaction

We present global potential energy surfaces for the three lowest triplet states in O(3P)+H2O(X1A1) collisions and present results of classical dynamics calculations on the O(3P)+H2O(X1A1)-->OH(X2pi)+OH(X2pi) reaction using these surfaces. The surfaces are spline-based fits of approximately 20,000 fixed geometry ab initio calculations at the complete-active-space self-consistent field+second-order perturbation theory (CASSCF+MP2) level with a O(4s3p2d1f)/H(3s2p) one electron basis set. Computed rate constants compare well to measurements in the 1000-2500 K range using these surfaces. We also compute the total, rovibrationally resolved, and differential angular cross sections at fixed collision velocities from near threshold at approximately 4 km s(-1) (16.9 kcal mol(-1) collision energy) to 11 km s(-1) (122.5 kcal mol(-1) collision energy), and we compare these computed cross sections to available space-based and laboratory data. A major finding of the present work is that above approximately 40 kcal mol(-1) collision energy rovibrationally excited OH(X2pi) products are a significant and perhaps dominant contributor to the observed 1-5 micro spectral emission from O(3P)+H2O(X1A1) collisions. Another important result is that OH(X2pi) products are formed in two distinct rovibrational distributions. The "active" OH products are formed with the reagent O atom, and their rovibrational distributions are extremely hot. The remaining "spectator" OH is relatively rovibrationally cold. For the active OH, rotational energy is dominant at all collision velocities, but the opposite holds for the spectator OH. Summed over both OH products, below approximately 50 kcal mol(-1) collision energy, vibration dominates the OH internal energy, and above approximately 50 kcal mol(-1) rotation is greater than vibrational energy. As the collision energy increases, energy is diverted from vibration to mostly translational energy. We note that the present fitted surfaces can also be used to investigate direct collisional excitation of H2O(X1A1) by O(3P) and also OH(X2pi)+OH(X2pi) collisions.     

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

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

Quasiclassical trajectory study of the vibrational quenching of hydroxyl radicals through collision with O atoms

The collisional removal of vibrationally excited OH radicals by O atoms is studied by the quasiclassical trajectory method. To evaluate the effect of different topological features on the scattering processes two different global potential energy surfaces, DMBE IV and TU, are used. Results for reactive, exchange, and inelastic scattering probabilities are reported for central collisions (with zero total angular momentum) with a fixed relative translational energy for vibrational levels of OH ranging from nu=1 to v=8. Vibrational state distributions of product molecules are also compared on the two potential energy surfaces. Both surfaces predict higher probabilities for reaction than for exchange or inelastic scattering. The vibrational state distributions of the product diatomic molecules are different on the two surfaces. In particular, the two surfaces give substantially different probabilities for multiquantum OH vibrational relaxation transitions OH(v)+O-->OH(v')+O.     

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.     

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.     

Direct dynamics simulations of O(3P) + HCl at hyperthermal collision energies

The dynamics of the O(3P) + HCl reaction at hyperthermal collision energies were investigated using the quasiclassical trajectory method. Stationary points on the OClH 3A" and 3A' potential energy surfaces (PESs) were also examined. The lowest transition state leading to OCl + H on the 3A" surface is 2.26 eV above the reagents at the CCSD(T)/cc-pVTZ level of theory. This saddle point is bent and product-like. Direct dynamics calculations at the MP2/cc-pVTZ level of theory were used to investigate the excitation functions for OH + Cl, OCl + H, and O + H + Cl formation. OCl is formed mainly from small-impact-parameter collisions, and the OCl + H excitation function peaks around 5 eV, where it is similar in magnitude to the OH + Cl excitation function. The shape of the OCl + H excitation function is discussed, and features are identified that should be general to hyperthermal collision dynamics.     

N(4S)+ NO(X2Π)→N2(X3Σg-)+O(3P) 反应的立体动力学研究

采用准经典轨线方法研究了在不同碰撞能下，碰撞反应N(4S)+NO(X2Π)→ N2(X3Σg- )+O(3P)在两个最低势能面3A 和 3A'上产物与反应物之间的矢量相关. 结果表明，对于不同的碰撞能，在两个势能面上反应产物的转动取向展示了不同的特征和趋势. 随着碰撞能的增加，发生在3A 势能面上的反应主要受外平面机理支配，而发生在 3A' 势能面上的反应倾向于受内平面机理支配. 这些差异来自于两个势能面的不同构型.     

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.     

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.     

State-to-state inelastic scattering of OH by HI: a comparison with OH-HCl and OH-HBr

Relative state-to-state cross sections and steric asymmetries have been measured for the scattering process: OH (X (2)Pi(32),v=0,J=32,M(J)=32,f)+HI ((1)Sigma,v=0,J<4)-->OH (X (2)Pi,v=0,Omega=12,J=12-52 and Omega=32,J=32-92,ef)+HI, at 690 cm(-1) collision energy. Comparison with the previously studied systems OH-HCl and OH-HBr reveals relevant features of the potential energy surfaces of these molecular systems. Some measured differences concerning the internal energy distribution after collision and the propensities for the impact with one or the other side of the OH molecule in scattering by HCl, HBr, and HI molecules are discussed.     

Coupled-states statistical investigation of vibrational and rotational relaxation of OH(2pi) by collisions with atomic hydrogen

We report state-to-state cross sections and thermal rate constants for vibrational and rotational relaxation of OH(2pi) by collision with H atoms. The cross sections are calculated by the coupled-states (CS) statistical method including the full open-shell character of the OH + H system. Four potential energy surfaces (PESs) ((1,3)A' and (1,3)A') describe the interaction of OH(X2pi) with H atoms. Of these, three are repulsive, and one (1A') correlates with the deep H2O well. Consequently, rotationally and ro-vibrationally inelastic scattering of OH in collisions with H can occur by scattering on the repulsive PESs, in a manner similar to the inelastic scattering of OH by noble gas atoms, or by collisions which enter the H2O well and then reemerge. At 300 K, we predict large (approximately 1 x 10(-10) cm3 molecule(-1) s(-1)) vibrational relaxation rates out of both v = 2 and v = 1, comparable to earlier experimental observations. This anomalously fast relaxation results from capture into the H2O complex. There exists a significant propensity toward formation of OH in the pi(A') lambda-doublet level. We also report state-resolved cross sections and rate constants for rotational excitation within the OH v = 0 manifold. Collisional excitation from the F1 to the F2 spin-orbit manifold leads to an inverted lambda-doublet population.     

Quasi-classical trajectory-gaussian binning study of the OH + D2 → HOD(v1',v2',v3') + D angle-velocity and vibrational distributions at a collision energy of 0.28 eV

The angle-velocity and product vibrational state distributions for the OH + D(2) reaction at a collision energy of 0.28 eV have been calculated using the quasi-classical trajectory-gaussian binning (QCT-GB) method and the Wu-Schatz-Lendvay-Fang-Harding (WSLFH) analytical potential energy surface. Comparison with high resolution molecular beam experiments shows that, differing from what happens when using the standard QCT method (i.e., histogram binning), very good results are obtained for both distributions. Hence, the strong differences previously observed between QCT and experimental results mainly come from an inadequate pseudoquantization of HOD rather than from other quantum effects. This is probably the first time that such a high level of agreement between theory and high resolution experimental data has been found in polyatomic reaction dynamics.     

Reaction pathway and potential barrier for the CaH product in the reaction of Ca(4s4p1P1) + H2 --> CaH(X2Sigma+) + H

The nascent CaH product in the reaction Ca(4s4p1P1) + H2 --> CaH(X2Sigma+) + H is obtained using a pump-probe technique. The CaH(v = 0,1) distributions, with a population ratio of CaH(v = 0)/CaH(v = 1) = 2.7+/-0.2, may be characterized by low Boltzmann rotational temperature. According to Arrhenius theory, the temperature dependence measurement yields a potential barrier of 3820+/-480 cm(-1) for the current reaction. As a result of the potential energy surfaces (PES) calculations, the reaction pathway favors a Ca insertion into the H2 bond along a (near) C2v geometric approach. As the H2 bond is elongated, the configurational mixing between the orbital components of the 4p and nearby low-lying 3d state with the same symmetry makes significant the nonadiabatic transition between the 5A' and 2A' surface in the repulsive limbs. Therefore, the collision species are anticipated to track along the 5A' surface, then undergo nonadiabatic transition to the inner limb of the 2A' surface, and finally cross to the reactive 1A' surface. The observed energy barrier probably accounts for the energy requirement to surmount the repulsive hill in the entrance. The findings of the nascent CaH product distributions may be reasonably interpreted from the nature of the intermediate structure and lifetime after the 2A'-1A' surface transition. The distinct product distributions between the Ca(4 1P1) and Mg(3 1P1) reactions with H2 may also be realized with the aid of the PES calculations.     

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.     

A quasiclassical trajectory study for the N(4S)+O2(X 3Sigmag-)-->NO(X 2Pi)+O(3P) reaction on the new 2A' and 4A' potential-energy surfaces

A quasiclassical trajectory study with the sixth-order explicit symplectic algorithm of the N(4S)+O2(X 3Sigmag-)-->NO(X 2Pi)+O(3P) atmospheric reaction has been performed by employing the new 2A' and 4A' potential-energy surfaces reported by Sayos et al. [J. Chem. Phys. 117, 670 (2002)]. For the translational temperature considered up to 10,000 K, the larger relative translational energy and the higher rovibrational levels of O2 molecule with respect to the previous works have been taken into account, and a clearer database about the character of the total reaction cross section has been presented in this work. The dependence of microscopic rate constants on the vibrational level of O2 molecule at T=3000, 5000, and 10,000 K has been exhibited, and we can see that the values of log10 k(T,v,J) vary almost linearly with the vibrational level of O2 molecule. The thermal rate constants at the translational temperature between 300 and 10,000 K have been evaluated and compared with the experimental and previous theoretical results. It is found that the thermal rate constants determined in this work have a better agreement with the experimental data and can provide a more valid theoretical reference at the translational temperature considered for the title reaction.     

Reaction dynamics of Y + O2--> YO(X,A',A)+ O(3P(J)) studied by the crossed-beam technique

The dynamics of the reaction, Y + O2--> YO + O was studied by using the crossed-beam technique at several collision energies from 10.3 to 52.0 kJ mol(-1). The Y atomic beam was generated by laser vaporization and crossed with the O2 beam at a right angle. Among the energetically accessible electronic states of YO, the formation of the A2Pi and A'2Delta states was observed by their chemiluminescence at all collision energies. By analyzing the chemiluminescence spectra of YO(A2Pi(1/2,3/2)-X2Sigma+), vibrational state distributions and relative populations of spin-orbit states were determined for YO(A2Pi(1/2,3/2)). At low collision energies, the vibrational distributions agree quite well with those expected from the statistical energy partitioning, while a little deviation from the statistical expectation was observed at the highest energy, 52.0 kJ mol(-1). The populations of two spin-orbit states are in good agreement with the statistical expectations at all collision energies. The vacuum ultraviolet laser-induced fluorescence (VUV-LIF) technique was employed to determine the distributions of spin-orbit states of the product O(3P(J)) at two collision energies, 20.7 and 52.0 kJ mol(-1). The line shapes of the O atom transitions were analyzed to determine relative branching ratio of the ground state to the excited states of YO, i.e. YO(X2Sigma+)+ O(3P(J))vs. YO(A2Pi and A'2Delta)+ O(3P(J)). The results showed that the electronically excited YO was formed with comparable amount with the ground state which is statistically more favorable, and suggested the occurrence of two mechanisms taking place in the title reaction.     

Cross sections and rate constants for the C(3P)+OH(X 2Pi)-->CO(X 1Sigma+)+H(2S) reaction using a quasiclassical trajectory method

First quasiclassical trajectory calculations have been carried out for the C(3P)+OH(X 2Pi)-->CO(X 1Sigma+)+H(2S) reaction using a recent ab initio potential energy surface for the ground electronic state, X 2A', of HCO/COH. Total and state-specific integral cross sections have been determined for a wide range of collision energies (0.001-1 eV). Then, thermal and state-specific rate constants have been calculated in the 1-500 K temperature range. The thermal rate constant varies from 1.78x10(-10) cm3 s-1 at 1 K down to 5.96x10(-11) cm3 s-1 at 500 K with a maximum value of 3.39x10(-10) cm3 s-1 obtained at 7 K. Cross sections and rate constants are found to be almost independent of the rovibrational state of OH.     

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.     

Reaction pathway for the nonadiabatic reaction of Ca(4s3d 1D)+H2-->CaH(X 2Sigma+)+H

The reaction pathway and the nascent CaH product distribution in the reaction Ca(4s3d (1)D)+H(2)-->CaH(X (2)Sigma(+))+H are obtained using a pump-probe technique. The Ca atom is first prepared in the 3 (1)D state by a two-photon absorption, and then in brief time delay the laser-induced fluorescence of the reaction product CaH is monitored. The CaH(v=0,1) distributions appear to be single peaked, as characterized by Boltzmann rotational temperature of 807+/-38 K (v=0) and 684+/-77 K (v=1). The vibrational population ratio of CaH(v=0)/CaH(v=1) is determined to be 3.3+/-0.1, while the v=2 population is not detectable. The fractions of the available energy partitioning into rotation, vibration, and translation are estimated to be 0.36+/-0.05, 0.28+/-0.04, and 0.36+/-0.05, respectively. With the aid of the potential energy surfaces calculations, the current reaction should favor a near C(2v) collision configuration. The temperature dependence measurement yields a positive slope, indicative of the reaction occurrence without any potential barrier. The colliding species are anticipated to follow an attractive 1B(2) (or 2A') surface and then transit nonadiabatically to the reactive ground state surface.