Quantum calculations of the O(3P)+H2-->OH+H reaction

Quantum scattering calculations are reported for the O(3P)+H2(v=0,1) reaction using chemically accurate potential energy surfaces of 3A' and 3A" symmetry. We present state-to-state reaction cross sections and rate coefficients as well as thermal rate coefficients for the title reaction using accurate quantum calculations. Our calculations yield reaction cross sections that are in quantitative accord with results of recent crossed molecular beam experiments. Comparisons with results obtained using the J-shifting calculations show that the J-shifting approximation is quite reliable for this system. Thermal rate coefficients from the exact calculations and the J-shifting approximation agree remarkably well with experimental results. Our calculations also reproduce the markedly different OH(v'=0)/OH(v'=1) branching in O(3P)+H2(v=1) reaction, observed in experiments that use different O(3P) atom sources. In particular, we show that the branching ratio is a strong function of the kinetic energy of the O(3P) atom.     

Quantum mechanical investigations of the N(4S)+O2(X 3Sigmag -)-->NO(X 2Pi)+O(3P) reaction

The reaction between energetic nitrogen atoms and oxygen molecules has received important attention in connection with nitric oxide chemistry in the lower thermosphere. We report time-independent quantum mechanical calculations of the N(4S)+O2-->NO+O reaction employing the X 2A' and a 4A' electronic potential energy surfaces of Sayos et al. [J. Chem. Phys. 117, 670 (2002)]. We confirm the production of highly vibrationally excited NO molecules, consistent with previous semiclassical and more recent time-dependent quantum wave packet studies. Calculations are carried out for total angular momentum quantum number J=0 and cross sections and rate coefficients are extracted using the J-shifting approximation. The results are in good agreement with available experimental and theoretical data.     

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.     

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

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

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.     

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.     

Vibrational state-dependent predissociation dynamics of ClO (A (2)Pi(3/2)): Insight from correlated fine structure branching ratios

We have studied the v'-dependent predissociation dynamics of the ClO A (2)Pi(3/2) state using velocity-map ion-imaging. Experimental final correlated state branching ratios, i.e. Cl((2)P(J=3/2,1/2)) + O((3)P(J=2,1,0)) channels, have been measured for v' = 6-11. We find that the branching ratios are highly variable and depend strongly on v', providing a window into the v'-dependent predissociation mechanism. A comparison of the experimental results with the recent model of Lane et al. (I. C. Lane, W. H. Howie and A. J. Orr-Ewing, Phys. Chem. Chem. Phys., 1999, 1, 3087) in both the diabatic and adiabatic limits suggests that the dynamics are closer to the diabatic limit. The overall Cl((2)P(J)) branching ratios are in good agreement with the diabatic model results. There are significant differences, however, between theory and experiment at the correlated state level, demonstrating the sensitivity of correlated measurements to the role of the exit channel coupling in the predissociation dynamics. The results highlight the need for more sophisticated quantum dynamical calculations to describe the correlated fine structure branching ratios in this system.     

Reactivity enhancement of ultracold O(3P)+H2 collisions by van der Waals interactions

The role of van der Waals forces in O((3)P)+H(2)(upsilon=1,j=0) collisions is investigated theoretically at low and ultralow temperatures. Quantum scattering calculations have been performed for zero total angular momentum using the lowest London-Eyring-Polanyi-Sato double-polynomial (3)A(") potential-energy surface reported by [Rogers et al., J. Phys. Chem. A 104, 2308 (2000)] and its recent BMS1 and BMS2 extensions developed by [Brandao et al., J. Chem. Phys. 121, 8861 (2004)] which provide a more accurate treatment of the van der Waals interaction. Our calculations show that van der Waals forces strongly influence chemical reactivity at ultracold translational energies. The presence of a zero-energy resonance for the BMS1 surface is found to enhance reactivity in the ultracold regime and shift the Wigner threshold to lower temperatures.     

Rate constant for OH(2 Pi)+O(3P)-->H(2S)+O2(3 Sigma g-) reaction on an improved ab initio potential energy surface and implications for the interstellar oxygen problem

The authors report a global potential energy surface for the ground electronic state of HO(2)(X (2)A(")), which improves upon the XXZLG potential [Xu and et al., J. Chem. Phys. 122, 244305 (2005)] with additional high-level ab initio points for the long-range interaction potential in the O+OH channel. Exact J=0 quantum mechanical reaction probabilities were calculated on the new potential and the rate constant for the title reaction was obtained using a J-shifting method. The calculated rate constant is in good agreement with available experimental values and our results predict a significantly lower rate at temperature range below 30 K, offering a possible explanation for the "interstellar oxygen problem."     

Influence of rotation and isotope effects on the dynamics of the N((2)D)+H(2) reactive system and of its deuterated variants

Integral cross sections and thermal rate constants have been calculated for the N((2)D)+H(2) reaction and its isotopic variants N((2)D)+D(2) and the two-channel N((2)D)+HD by means of quasiclassical trajectory and statistical quantum-mechanical model methods on the latest ab initio potential-energy surface [T.-S. Ho et al., J. Chem. Phys. 119, 3063 (2003)]. The effect of rotational excitation of the diatom on the dynamics of these reactions has been investigated and interesting discrepancies between the classical and statistical model calculations have been found. Whereas a net effect of reagent rotation on reactivity is always observed in the classical calculations, only a very slight effect is observed in the case of the asymmetric N((2)D)+HD reaction for the statistical quantum-mechanical method. The thermal rate constants calculated on this Potential-Energy Surface using quasiclassical trajectory and statistical model methods are in good agreement with the experimental determinations, although the latter are somewhat larger. A reevaluation of the collinear barrier of the potential surface used in the present study seems timely. Further theoretical and experimental studies are needed for a full understanding of the dynamics of the title reaction.     

Quantum reactive scattering of H + hydrocarbon reactions

A practical quantum-dynamical method is described for predicting accurate rate constants for general chemical reactions. The ab initio potential energy surfaces for these reactions can be built from a minimal number of grid points (average of 50 points) and expressed in terms of analytical functionals. All the degrees of freedom except the breaking and forming bonds are optimised using the MP2 method with a cc-pVTZ basis set. Single point energies are calculated on the optimised geometries at the CCSD(T) level of theory with the same basis set. The dynamics of these reactions occur on effective reduced dimensionality hyper-surfaces accounting for the zero-point energy of the optimised degrees of freedom. Bonds being broken and formed are treated with explicit hyperspherical time independent quantum dynamics. Application of the method to the H + CH(4)--> H(2)+ CH(3), H + C(2)H(6)--> H(2)+ C(2)H(5), H + C(3)H(8)--> H(2)+n-C(3)H(7)/H(2)+i-C(3)H(7) and H + CH(3)OH --> H(2)+ CH(3)O/H(2)+ CH(2)OH reactions illustrate the potential of the approach in predicting rate constants, kinetic isotope effects and branching ratios. All studied reactions exhibit large quantum tunneling in the rate constants at lower temperatures. These quantum calculations compare well with the experimental results.     

A global 12-dimensional ab initio potential energy surface and dynamical studies for the SiH4+H-->SiH3+H2 reaction

A global 12-dimensional ab initio interpolated potential energy surface (PES) for the SiH(4)+H-->SiH(3)+H(2) reaction is presented. The ab initio calculations are based on the unrestricted quadratic configuration interaction treatment with all single and double excitations together with the cc-pVTZ basis set, and the modified Shepard interpolation method of Collins and co-workers [K. C. Thompson et al., J. Chem. Phys. 108, 8302 (1998); M. A. Collins, Theor. Chem. Acc. 108, 313 (2002); R. P. A. Bettens and M. A. Collins, J. Chem. Phys. 111, 816 (1999)] is applied. Using this PES, classical trajectory and variational transition state theory calculations have been carried out, and the computed rate constants are in good agreement with the available experimental data.     

Reaction rates for O + HD → OH + D and O + HD → OD + H

We have calculated reaction rates for the reactions O + HD → OH + D and O + DH → OD + H using improved canonical variational transition state theory and least-action ground-state transmission coefficients with an ab initio potential energy surface. The kinetic isotope effects are in good agreement with experiment. The optimized tunneling paths and properties of the variational transition states and the rate enhancement for vibrationally excited reactants are also presented and compared with those for the isotopically unsubstituted reaction O + H₂ → OH + H. The thermal reactions at low and room temperature are predicted to occur by tunneling at extended configurations, i.e., to initiate early on the reaction path and to avoid the saddle point regions. Tunneling also dominates the low and room temperature reactions for excited vibrational states, but in these cases the results are not as sensitive to the nature of the tunneling path. Overbarrier mechanisms dominate for both thermal and excited-vibrational state reactions for T > 600 K. For the excited-state reaction (with initial vibrational quantum number n > 0) a transition state switch occurs for T > 1000 K for the O + HD(n = 1) → OD + H case and for T > 1500 K for the O + DH(n = 1) → OD + H reaction, and this may be a general phenomenon for excited-state reactions at higher temperature. In the present case the switch occurs from an early variational transition state where the vibrationally adiabatic approximation is expected to be valid to a tighter variational transition state where nonadiabatic effects are probably important and should be included.     

A reinvestigation of the kinetics and the mechanism of the CH3C(O)O2 + HO2 reaction using both experimental and theoretical approaches

The kinetics and the mechanism of the reaction CH(3)C(O)O(2)+ HO(2) were reinvestigated at room temperature using two complementary approaches: one experimental, using flash photolysis/UV absorption technique and one theoretical, with quantum chemistry calculations performed using the density functional theory (DFT) method with the three-parameter hybrid functional B3LYP associated with the 6-31G(d,p) basis set. According to a recent paper reported by Hasson et al., [J. Phys. Chem., 2004, 108, 5979-5989] this reaction may proceed by three different channels: CH(3)C(O)O(2)+ HO(2)--> CH(3)C(O)OOH + O(2) (1a); CH(3)C(O)O(2)+ HO(2)--> CH(3)C(O)OH + O(3) (1b); CH(3)C(O)O(2)+ HO(2)--> CH(3)C(O)O + OH + O(2) (1c). In experiments, CH(3)C(O)O(2) and HO(2) radicals were generated using Cl-initiated oxidation of acetaldehyde and methanol, respectively, in the presence of oxygen. The addition of amounts of benzene in the system, forming hydroxycyclohexadienyl radicals in the presence of OH, allowed us to answer that channel (1c) is <10%. The rate constant k(1) of reaction (1) has been finally measured at (1.50 +/- 0.08) x 10(-11) cm(3) molecule(-1) s(-1) at 298 K, after having considered the combination of all the possible values for the branching ratios k(1a)/k(1,)k(1b)/k(1,)k(1c)/k(1) and has been compared to previous measurements. The branching ratio k(1b)/k(1), determined by measuring ozone in situ, was found to be equal to (20 +/- 1)%, a value consistent with the previous values reported in the literature. DFT calculations show that channel (1c) is also of minor importance: it was deduced unambiguously that the formation of CH(3)C(O)OOH + O(2) (X (3)Sigma(-)(g)) is the dominant product channel, followed by the second channel (1b) leading to CH(3)C(O)OH and singlet O(3) and, much less importantly, channel (1c) which corresponds to OH formation. These conclusions give a reliable explanation of the experimental observations of this work. In conclusion, the present study demonstrates that the CH(3)C(O)O(2)+ HO(2) is still predominantly a radical chain termination reaction in the tropospheric ozone chain formation processes.     

Angular distributions and angular momentum alignment of O((3)P(J)) atoms formed in the photolysis of O(2)via the Herzberg continuum

Sliced velocity-map imaging has been used to measure photofragment scattering distributions for the O((3)P(2)) and O((3)P(1)) products of O(2) photolysis following laser excitation into the Herzberg continuum between 205 and 241 nm. The images were analysed to extract the photofragment spatial anisotropy parameter, β, together with the alignment parameters a(∥), a(⊥), a(⊥), and Re[a(∥, ⊥)]. Our alignment measurements bridge the gap between the recent 193 nm measurement of Brouard et al., Phys. Chem. Chem. Phys., 2006, 8, 5549 and those of Alexander et al., J. Chem. Phys., 2003, 118, 10566 at 222 and 237 nm, and extend out to the threshold at 241 nm. Our measured parameters show no strong dependence on photolysis wavelength. Near the threshold we were able to separate the contributions from the O((3)P(2)) + O((3)P(2)) and O((3)P(2)) + O((3)P(1)) channels, and found significantly different photofragment alignments for the two cases.     

Recombination of ozone via the chaperon mechanism

The recombination of ozone via the chaperon mechanism, i.e., ArO+O2 --> Ar+O3 and ArO2+O --> Ar+O3, is studied by means of classical trajectories and a pairwise additive Ar-O3 potential energy surface. The recombination rate coefficient has a strong temperature dependence, which approximately can be described by T(-n) with n approximately 3. It is negligible for temperatures above 700 K or so, but it becomes important for low temperatures. The calculations unambiguously affirm the conclusions of Hippler et al. [J. Chem. Phys. 93, 6560 (1990)] and Luther et al. [Phys. Chem. Chem. Phys. 7, 2764 (2005)] that the chaperon mechanism makes a sizable contribution to the recombination of O3 at room temperature and below. The dependence of the chaperon recombination rate coefficient on the isotopomer, studied for two different isotope combinations, is only in rough qualitative agreement with the experimental data. The oxygen atom isotope exchange reaction involving ArO and ArO2 van der Waals complexes is also investigated; the weak binding of O or O2 to Ar has only a small effect.     

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

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

The relaxation of OH (v=1) and OD (v=1) by H2O and D2O at temperatures from 251 to 390 K

We report rate coefficients for the relaxation of OH(v=1) and OD(v=1) by H2O and D2O as a function of temperature between 251 and 390 K. All four rate coefficients exhibit a negative dependence on temperature. In Arrhenius form, the rate coefficients for relaxation (in units of 10(-12) cm3 molecule-1 s-1) can be expressed as: for OH(v=1)+H2O between 263 and 390 K: k=(2.4+/-0.9) exp((460+/-115)/T); for OH(v=1)+D2O between 256 and 371 K: k=(0.49+/-0.16) exp((610+/-90)/T); for OD(v=1)+H2O between 251 and 371 K: k=(0.92+/-0.16) exp((485+/-48)/T); for OD(v=1)+D2O between 253 and 366 K: k=(2.57+/-0.09) exp((342+/-10)/T). Rate coefficients at (297+/-1 K) are also reported for the relaxation of OH(v=2) by D2O and the relaxation of OD(v=2) by H2O and D2O. The results are discussed in terms of a mechanism involving the formation of hydrogen-bonded complexes in which intramolecular vibrational energy redistribution can occur at rates competitive with re-dissociation to the initial collision partners in their original vibrational states. New ab initio calculations on the H2O-HO system have been performed which, inter alia, yield vibrational frequencies for all four complexes: H2O-HO, D2O-HO, H2O-DO and D2O-DO. These data are then employed, adapting a formalism due to Troe (J. Troe, J. Chem. Phys., 1977, 66, 4758), in order to estimate the rates of intramolecular energy transfer from the OH (OD) vibration to other modes in the complexes in order to explain the measured relaxation rates-assuming that relaxation proceeds via the hydrogen-bonded complexes.