Electronic quenching of OH A 2Sigma+ radicals in single collision events with H2 and D2: a comprehensive quantum state distribution of the OH X 2Pi products

A pump-probe laser-induced fluorescence technique has been used to examine the nascent OH X (2)Pi product state distribution arising from non-reactive quenching of electronically excited OH A (2)Sigma(+) by molecular hydrogen and deuterium under single-collision conditions. The OH X (2)Pi products were detected in v'=0, 1 and 2; the distribution peaks in v'=0 and decreases monotonically with increasing vibrational excitation. In all vibrational levels probed, the OH X (2)Pi products are found to be highly rotationally excited, the distribution peaking at N'=15 when H(2) was used as the collision partner and N'=17 for D(2). A marked propensity for production of Pi(A') Lambda-doublet levels was observed, while both OH X (2)Pi spin-orbit manifolds were equally populated. These observations are interpreted as dynamical signatures of the nonadiabatic passage of the OH + H(2)/D(2) system through the seams of conical intersection that couple the excited state (2 (2)A') and ground state (1 (2)A') surfaces.     

Cross sections and thermal rate constants for the isotope exchange reaction: D(2S)+OH(2Pi)-->OD(2Pi)+H(2S)

We report state-to-state and overall thermal rate constants for the isotope exchange reaction D((2)S)+OH((2)Pi)-->OD((2)Pi)+H((2)S) for 0 K相似文献   

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.     

A new mechanism for the production of highly vibrationally excited OH in the mesosphere: an ab initio study of the reactions of O2(A 3Sigmau+ and A' 3Deltau)+H

In an attempt to explain the observed nightglow emission from OH(v=10) in the mesosphere that has the energy greater than the exothermicity of the H+O(3) reaction, potential energy surfaces were calculated for reactions of high lying electronic states of O(2)(A (3)Sigma(u) (+) and A' (3)Delta(u)) with atomic hydrogen H((2)S) to produce the ground state products OH((2)Pi)+O((3)P). From collinear two-dimensional scans, several adiabatic and nonadiabatic pathways have been identified. Multiconfigurational single and double excitation configuration interaction calculations show that the adiabatic pathways on a (4)Delta potential surface from O(2)(A' (3)Delta)+H and a (4)Sigma(+) potential surface from O(2)(A (3)Sigma(u) (+))+H are the most favorable, with the zero-point corrected barrier heights of as low as 0.191 and 0.182 eV, respectively, and the reactions are fast. The transition states for these pathways are collinear and early, and the reaction coordinate suggests that the potential energy release of ca. 3.8 eV (larger than the energy required to excite OH to v=10) is likely to favor high vibrational excitation.     

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.     

A combined crossed beam and theoretical investigation of O(3P) + C3H3 --> C3H2 + OH

The radical-radical reaction dynamics of ground-state atomic oxygen [O(3P)] with propargyl radicals (C3H3) has first been investigated in a crossed beam configuration. The radical reactants O(3P) and C3H3 were produced by the photodissociation of NO2 and the supersonic flash pyrolysis of precursor propargyl bromide, respectively. A new exothermic channel of O(3P) + C3H3 --> C3H2 + OH was identified and the nascent distributions of the product OH in the ground vibrational state (X 2Pi:nu" = 0) showed bimodal rotational excitations composed of the low- and high-N" components without spin-orbit propensities. The averaged ratios of Pi(A')/Pi(A") were determined to be 0.60 +/- 0.28. With the aid of ab initio theory it is predicted that on the lowest doublet potential energy surface, the reaction proceeds via the addition complexes formed through the barrierless addition of O(3P) to C3H3. The common direct abstraction pathway through a collinear geometry does not occur due to the high entrance barrier in our low collision energy regime. In addition, the major reaction channel is calculated to be the formation of propynal (CHCCHO) + H, and the counterpart C3H2 of the probed OH product in the title reaction is cyclopropenylidene (1c-C3H2) after considering the factors of barrier height, reaction enthalpy and structural features of the intermediates formed along the reaction coordinate. On the basis of the statistical prior and rotational surprisal analyses, the ratio of population partitioning for the low- and high-N" is found to be about 1:2, and the reaction is described in terms of two competing addition-complex mechanisms: a major short-lived dynamic complex and a minor long-lived statistical complex. The observed unusual reaction mechanism stands in sharp contrast with the reaction of O(3P) with allyl radical (C3H5), a second significant conjugated hydrocarbon radical, which shows totally dynamic processes [J. Chem. Phys. 117, 2017 (2002)], and should be understood based upon the characteristic electronic structures and reactivity of the intermediates on the potential energy surface.     

Evidence for excited spin-orbit state reaction dynamics in F+H2: theory and experiment

We describe fully quantum, time-independent scattering calculations of the F+H2-->HF+H reaction, concentrating on the HF product rotational distributions in v'=3. The calculations involved two new sets of ab initio potential energy surfaces, based on large basis set, multireference configuration-interaction calculations, which are further scaled to reproduce the experimental exoergicity of the reaction. In addition, the spin-orbit, Coriolis, and electrostatic couplings between the three quasidiabatic F+H2 electronic states are included. The calculated integral cross sections are compared with the results of molecular beam experiments. At low collision energies, a significant fraction of the reaction is due to Born-Oppenheimer forbidden, but energetically allowed reaction of F in its excited (2P 1/2) spin-orbit state. As the collision energy increases, the Born-Oppenheimer allowed reaction of F in its ground (2P 3/2) spin-orbit state rapidly dominates. Overall, the calculations agree reasonably well with the experiment, although there remains some disagreement with respect to the degree of rotational excitation of the HF(v'=3) products as well as with the energy dependence of the reactive cross sections at the lowest collision energies.     

Isotope branching and tunneling in O(3P)+HD-->OH+D; OD+H reactions

The O((3)P)+HD and O((3)P)+D(2) reactions are studied using quantum scattering calculations and chemically accurate potential energy surfaces developed for the O((3)P)+H(2) system by Rogers et al. [J. Phys. Chem. A 104, 2308 (2000)]. Cross sections and rate coefficients for OH and OD products are calculated using accurate quantum methods as well as the J-shifting approximation. The J-shifting approach is found to work remarkably well for both O+HD and O+D(2) collisions. The reactions are dominated by tunneling at low temperatures and for the O+HD reaction the hydrogen atom transfer leading to the OH product dominates at low temperatures. Our result for the OH/OD branching ratio is in close agreement with previous calculations over a wide range of temperatures. The computed OH/OD branching ratios are also in close agreement with experimental results of Robie et al. [Chem. Phys. Lett. 134, 579 (1987)] at temperatures above 400 K but the theoretical results do not reproduce the rapid rise in the experimental values of the branching ratio for temperatures lower than 350 K. We believe that new measurements could resolve the long-standing discrepancy between experiment and theory for this benchmark reaction.     

Quasi-classical trajectory study of the S + OH → SO + H reaction: from reaction probability to thermal rate constant

First quasi-classical trajectory calculations have been carried out for the S((3)P) + OH(X?(2)Π) → SO(X?(3)Σ(-)) + H((2)S) reaction on an ab initio global potential energy surface for the ground electronic state, X?(2)A', of HSO. Cross sections, computed for collision energies up to 1 eV, show no energy threshold and decrease with the increasing collision energy. Rate constants have been calculated in the 5-500 K temperature range. The thermal rate constant is in good agreement with approximate quantum results, while a disagreement is found at 298 K with the experimental data. Product energy distributions have also been reported at four collision energies from 0.001 to 0.5 eV. The shapes of the rovibrational and angular distributions suggest the formation of an intermediate complex that is more and more long-lived as the collision energy increases.     

Electronic quenching of OH A 2Sigma+ radicals in single collision events with molecular hydrogen: quantum state distribution of the OH X 2Pi products

We report a combined experimental and theoretical investigation of the nonreactive quenching channel resulting from electronic quenching of OH A 2Sigma+ by molecular hydrogen. The experiments utilize a pump-probe scheme to determine the OH X 2Pi population distribution following collisional quenching in a pulsed supersonic expansion. The pump laser excites OH A 2Sigma+ (nu'=0, N'=0), which has a significantly reduced fluorescence lifetime due to quenching by H2. The probe laser monitors the OH X 2Pi (nu", N") population via laser-induced fluorescence on various A-X transitions under single collision conditions. The experiments reveal a high degree of rotational excitation (N") of the quenched OH X 2Pi products observed in nu"=1 and 2 as well as a pronounced propensity for quenching into the Pi(A') Lambda-doublet level. These experiments have been supplemented by extensive multireference, configuration-interaction calculations aimed at exploring the topology of the relevant potential energy surfaces. Electronic quenching of OH A 2Sigma+ by H2 proceeds through conical intersections between two potentials of A' reflection symmetry (in planar geometry) that correlate with the electronically excited A 2Sigma+ and ground X 2Pi states of OH. The conical intersections occur in high-symmetry geometries, in which the O side of OH points toward H2. Corroborating and extending earlier work of Hoffman and Yarkony [J. Chem. Phys. 113, 10091 (2000)], these calculations reveal a steep gradient away from the OH-H2 conical intersection as a function of both the OH orientation and interfragment distance. The former will give rise to a high degree of OH rotational excitation, as observed for the quenched OH X 2Pi products.     

Two photon dissociation dynamics of NO2 and NO2 + H2O

To explore the dynamics of OH formation from two photon absorbed NO(2) with H(2)O, a high-level multiconfigurational perturbation theory was used to map the potential energy profiles of NO(2) dissociation to O ((1)D) + NO (X(2)Π), and subsequent hydrogen abstraction producing 2OH (X(2)Π) + NO (X(2)Π) in the highly excited S(PP) (?(2)A', (2)ππ*) state. The ground state NO(2) is promoted to populate in the S(NP1) (?(2)A", (2)nπ*) intermediate state by one photon absorption at ~440 nm, one thousandth of which is further excited to S(PP) (?(2)A', (2)ππ*) state and undergoes a medium-sized barrier (~11.0 kcal/mol) to give rise to OH radicals. In comparison with the hydrogen abstraction reaction in highly vibrationally excited NO(2) ground state, two photon absorption facilitates NO(2) dissociation to O ((1)D) and O ((1)D) + H(2)O → 2OH (X(2)Π) but results in low quantum yield of NO(2)** since there is a weak absorption upon the second beam light at ~440 nm. It can be concluded that the reaction of two photon absorbed NO(2) with H(2)O makes negligible contributions to the formation of OH radicals. In contrast, single photon absorption at <554 nm is a possible process on the basis of the present and previous computations.     

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.     

Comparison of experimental time-of-flight spectra of the HF products from the F+H2 reaction with exact quantum mechanical calculations

High resolution HF product time-of-flight spectra measured for the reactive scattering of F atoms from n-H2(p-H2) molecules at collision energies between 69 and 81 meV are compared with exact coupled-channel quantum mechanical calculations based on the Stark-Werner ab initio ground state potential energy surface. Excellent agreement between the experimental and computed rotational distributions is found for the HF product vibrational states v'=1 and v'=2. For the v'=3 vibrational state the agreement, however, is less satisfactory, especially for the reaction with p-H2. The results for v'=1 and v'=2 confirm that the reaction dynamics for these product states is accurately described by the ground electronic state 1 (2)A' potential energy surface. The deviations for HF(v'=3, j' > or =2) are attributed to an enhancement of the reaction resulting from the 25% fraction of excited ((2)P(12)) fluorine atoms in the reactant beam.     

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.     

Quasi-classical trajectory study of Si+O2-->SiO+O reaction

Quasi-classical trajectory calculations for the Si(3P)+O2(X 3Sigmag-)-->SiO(X 1Sigma+)+O(1D) reaction have been carried out using the analytical ground 1A' potential energy surface (PES) recently reported by Dayou and Spielfiedel [J. Chem. Phys. 119, 4237 (2003)]. The reaction has been studied for a wide range of collision energies (0.005-0.6 eV) with O2 in its ground rovibrational state. The barrierless PES leads to a decrease of the total reaction cross section with increasing collision energy. It has been brought to evidence that the reaction proceeds through different reaction mechanisms whose contributions to reactivity are highly dependent on the collision energy range. At low collision energy an abstraction mechanism occurs involving the collinear SiOO potential well. The associated short-lived intermediate complex leads to an inverted vibrational distribution peaked at v'=3 and low rotational excitation of SiO(v',j') with a preferentially backward scattering. At higher energies the reaction proceeds mainly through an insertion mechanism involving the bent and linear OSiO deep potential wells and associated long-lived intermediate complexes, giving rise to nearly statistical energy disposals into the product modes and a forward-backward symmetry of the differential cross section.     

Direct evidence for nonadiabatic dynamics in atom+polyatom reactions: crossed-jet laser studies of F+D2O-->DF+OD

Quantum-state-resolved reactive-scattering dynamics of F+D(2)O-->DF+OD have been studied at E(c.m.)=5(1) kcal/mol in low-density crossed supersonic jets, exploiting pulsed discharge sources of F atom and laser-induced fluorescence to detect the nascent OD product under single-collision conditions. The product OD is formed exclusively in the v(OD)=0 state with only modest rotational excitation ( =0.50(1) kcal/mol), consistent with the relatively weak coupling of the 18.1(1) kcal/mol reaction exothermicity into "spectator" bond degrees of freedom. The majority of OD products [68(1)%] are found in the ground ((2)Pi(32) (+/-)) spin-orbit state, which adiabatically correlates with reaction over the lowest and only energetically accessible barrier (DeltaE( not equal) approximately 4 kcal/mol). However, 32(1)% of molecules are produced in the excited spin-orbit state ((2)Pi(12) (+/-)), although from a purely adiabatic perspective, this requires passage over a DeltaE( not equal) approximately 25 kcal/mol barrier energetically inaccessible at these collision energies. This provides unambiguous evidence for nonadiabatic surface hopping in F+D(2)O atom abstraction reactions, indicating that reactive-scattering dynamics even in simple atom+polyatom systems is not always isolated on the ground electronic surface. Additionally, the nascent OD rotational states are well fitted by a two-temperature Boltzmann distribution, suggesting correlated branching of the reaction products into the DF(v=2,3) vibrational manifold.     

Global potential energy surfaces for O((3)P) + H2O((1)A1) collisions

Global analytic potential energy surfaces for O((3)P) + H(2)O((1)A(1)) collisions, including the OH + OH hydrogen abstraction and H + OOH hydrogen elimination channels, are presented. Ab initio electronic structure calculations were performed at the CASSCF + MP2 level with an O(4s3p2d1f)/H(3s2p) one electron basis set. Approximately 10(5) geometries were used to fit the three lowest triplet adiabatic states corresponding to the triply degenerate O((3)P) + H(2)O((1)A(1)) reactants. Transition state theory rate constant and total cross section calculations using classical trajectories to collision energies up to 120?kcal mol(-1) (～11?km s(-1) collision velocity) were performed and show good agreement with experimental data. Flux-velocity contour maps are presented at selected energies for H(2)O collisional excitation, OH + OH, and H + OOH channels to further investigate the dynamics, especially the competition and distinct dynamics of the two reactive channels. There are large differences in the contributions of each of the triplet surfaces to the reactive channels, especially at higher energies. The present surfaces should support quantitative modeling of O((3)P) + H(2)O((1)A(1)) collision processes up to ～150?kcal mol(-1).     

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.