Spin-orbit effect in the energy pooling reaction O2(a 1Delta)+O2(a 1Delta)-->O2(b 1Sigma)+O2(X 3Sigma)

Five-dimensional nonadiabatic quantum dynamics studies have been carried out on two new potential energy surfaces of S(2)((1)A(')) and T(7)((3)A(")) states for the title oxygen molecules collision with coplanar configurations, along with the spin-orbit coupling between them. The ab initio calculations are based on complete active state second-order perturbation theory with the 6-31+G(d) basis set. The calculated spin-orbit induced transition probability as a function of collision energy is found to be very small for this energy pooling reaction. The rate constant obtained from a uniform J-shifting approach is compared with the existing theoretical and experimental data, and the spin-orbit effect is also discussed in this electronic energy-transfer process.     

A time-dependent wave packet study of the vibronic and spin-orbit interactions in the dynamics of Cl((2)P)+H(2)-->HCl(X (1)Sigma(g) (+))+H((2)S) reaction

We investigate the vibronic and spin-orbit (SO) coupling effects in the state-selected dynamics of the title reaction with the aid of a time-dependent wave packet approach. The ab initio potential energy surfaces of Capecchi and Werner [Science 296, 715 (2002)] have been employed for this purpose. Collinear approach of the Cl((2)P) atom to the H(2) molecule splits the degeneracy of the (2)P state and gives rise to (2)Sigma and (2)Pi electronic states. These two surfaces form a conical intersection at this geometry. These states transform as 1 (2)A('), 1 (2)A("), and 2 (2)A('), respectively, at the nonlinear configurations of the nuclei. In addition, the SO interaction due to Cl atom further splits these states into (2)Sigma(1/2), (2)Pi(3/2), and (2)Pi(1/2) components at the linear geometry. The ground-state reagent Cl((2)P(3/2))+H(2) correlates with (2)Sigma(1/2) and (2)Pi(3/2), where as the SO excited reagent Cl(*)((2)P(1/2))+H(2) correlates with (2)Pi(1/2) at the linear geometry. In order to elucidate the impact of the vibronic and SO coupling effects on the initial state-selected reactivity of these electronic states we carry out quantum scattering calculations based on a flux operator formalism and a time-dependent wave packet approach. In this work, total reaction probabilities and the time dependence of electronic population of the system by initiating the reaction on each of the above electronic states are presented. The role of conical intersection alone on the reaction dynamics is investigated with a coupled two-state model and for the total angular momentum J=0 (neglecting the electronic orbital angular momentum) both in a diabatic as well as in the adiabatic electronic representation. The SO interaction is then included and the dynamics is studied with a coupled three-state model comprising six diabatic surfaces for the total angular momentum J=0.5 neglecting the Coriolis Coupling terms of the Hamiltonian. Companion calculations are carried out for the uncoupled adiabatic and diabatic surfaces in order to explicitly reveal the impact of two different surface coupling mechanisms in the dynamics of this prototypical reaction.     

Ab initio analytical potential energy surface and quasiclassical trajectory study of the O(+)((4)S)+H(2)(X (1)Sigma(g) (+))-->OH(+)(X (3)Sigma(-))+H((2)S) reaction and isotopic variants

An analytical potential energy surface (PES) representation of the O(+)((4)S)+H(2)(X (1)Sigma(g) (+)) system was developed by fitting around 600 CCSD(T)/cc-pVQZ ab initio points. Rate constant calculations for this reaction and its isotopic variants (D(2) and HD) were performed using the quasiclassical trajectory (QCT) method, obtaining a good agreement with experimental data. Calculations conducted to determine the cross section of the title reaction, considering collision energies (E(T)) below 0.3 eV, also led to good accord with experiments. This PES appears to be suitable for kinetics and dynamics studies. Moreover, the QCT results show that, although the hypotheses of a widely used capture model are not satisfied, the resulting expression for the cross section can be applied within a suitable E(T) interval, due to errors cancellation. This could be a general situation regarding the application of this simple model to ion-molecule processes.     

Differential cross sections and product energy distributions for the C(3P)+OH(X 2Pi)-->CO(X 1Sigma+)+H(2S) reaction using a quasiclassical trajectory method

Quasiclassical trajectory calculations have been carried out for the C((3)P)+OH(X (2)Pi)-->CO(X (1)Sigma(+))+H((2)S) reaction using a recent ab initio potential energy surface for the ground electronic state X (2)A(') of COH. Differential cross sections (DCSs), and product vibrational, rotational and translational distributions have been determined for a wide range of collision energies (0.001-1 eV). The role of excitations (rotation or vibration) of the OH reactant on these quantities has been investigated. Product vibrational, rotational, and translational distributions are found to be almost independent on the rovibrational state of OH, whereas DCSs show a weak dependence on the initial rotational state of OH. We also analyze the results using a study based on the lifetime of the intermediate complex and on the kinematic constraint associated with the mass combination.     

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.     

A crossed beam study of (18)O((3)P)+NO(2) and (18)O((1)D)+NO(2): Isotope exchange and O(2)+NO formation channels

The products and dynamics of the reactions (18)O((3)P)+NO(2) and (18)O((1)D)+NO(2) have been investigated using crossed beams and provide new constraints on the structures and lifetimes of the reactive nitrogen trioxide intermediates formed in collisions of O((3)P) and O((1)D) with NO(2). For each reaction, two product channels are observed - isotope exchange and O(2)+NO formation. From the measured product signal intensities at collision energies of ～6 to 9.5 kcal∕mol, the branching ratio for O(2)+NO formation vs. isotope exchange for the O((3)P)+NO(2) reaction is 52(+6∕-2)% to 48(+2∕-6)%, while that for O((1)D)+NO(2) is 97(+2∕-12)% to 3(+12∕-2)%. The branching ratio for the O((3)P)+NO(2) reaction derived here is similar to the ratio measured in previous kinetics studies, while this is the first study in which the products of the O((1)D)+NO(2) reaction have been determined experimentally. Product energy and angular distributions are derived for the O((3)P)+NO(2) isotope exchange and the O((1)D)+NO(2)→O(2)+NO reactions. The results demonstrate that the O((3)P)+NO(2) isotope exchange reaction proceeds by an NO(3)? complex that is long-lived with respect to its rotational period and suggest that statistical incorporation of the reactant (18)O into the product NO(2) (apart from zero point energy isotope effects) likely occurs. In contrast, the (18)O((1)D)+NO(2)→O(2)+NO reaction proceeds by a direct "stripping" mechanism via a short-lived (18)O-O-NO? complex that results in the occurrence of (18)O in the product O(2) but not in the product NO. Similarly, (18)O is detected in O(2) but not NO for the O((3)P)+NO(2)→O(2)+NO reaction. Thus, even though the product energy and angular distributions for O((3)P)+NO(2)→O(2)+NO derived from the experimental data are uncertain, these results for isotope labeling under single collision conditions support previous kinetics studies that concluded that this reaction proceeds by an asymmetric (18)O-O-NO? intermediate and not by a long-lived symmetric NO(3)? complex, as earlier bulk isotope labeling experiments had concluded. Applicability of these results to atmospheric chemistry is also discussed.     

Full six-dimensional nonadiabatic quantum dynamics calculation for the energy pooling reaction O(2)(a (1)Delta)+O(2)(a (1)Delta)-->O(2)(b (1)Sigma)+O(2)(X (3)Sigma)

Six new potential energy surfaces of four singlet states and two triplet states for the title oxygen molecule reaction along with the spin-orbit coupling among them have been constructed from the complete active space second-order perturbation theory with a 6-311+G(d) basis. Accurate integral cross sections are calculated with a full six-dimensional nonadiabatic time-dependent quantum wave packet method. The thermal rate constant based on the integral cross sections agrees well with the result of the experimental measurements, and the intersystem crossing effects are also discussed in this electronic energy-transfer process.     

Reaction dynamics of Si(3PJ)+O2-->SiO(X 1 Sigma+)+O studied by a crossed-beam laser-induced fluorescence technique

Oxidation reaction of the ground state Si atom was studied by using a crossed molecular beam technique at 13.0 kJ/mol of collision energy. The Si atomic beam was generated by laser vaporization and crossed with the oxygen molecular beam at right angle. Products at the crossing region were detected by the laser-induced fluorescence (LIF). The LIF of SiO(A 1 Pi-X 1 Sigma+) was used to determine the vibrational state distribution of the electronic ground state, SiO(X 1 Sigma+). The determined distribution was inverted with the maximum population at v"=4, and in good agreement with the recent quasiclassical trajectory calculation on the singlet potential energy surface. The agreement suggested that an abstraction mechanism is dominant at the collision energy studied here. One of the counterproducts, O(3PJ), was also observed by the vacuum ultraviolet LIF and the distribution of the spin-orbit levels were determined. The formation of O(3PJ) was consistent with the significant population of v"=7 and 8 states of SiO, which could be explained by the presence of the triplet product channel with higher exothermicity.     

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

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

Chemiluminescent reaction of Ba(3P) with N2O at hyperthermal collision energies: rotational alignment of the BaO(A 1Sigma+) product

The chemiluminescent reaction Ba(6s6p (3)P)+N(2)O was studied at an average collision energy of 1.56 eV in a beam-gas arrangement. Ba((3)P) was produced by laser ablation of barium, which resulted in a broad collision energy distribution extending up to approximately 5.7 eV. A series of experiments was made to extract the Ba((3)P) contribution to chemiluminescence from that corresponding to Ba 6s(2) (1)S0 and 6s5d (3)D, which are the other two most populated states in the atomic beam. The fully dispersed polarized chemiluminescence spectra at 400-600 nm from the title reaction were recorded and assigned to a BaO molecule excited in the A (1)Sigma+ level. In addition, the average and wavelength-resolved degrees of polarization associated to the parallel BaO(A (1)Sigma+-->X (1)Sigma+) emission are reported. The analysis of the average polarization degree show that the BaO(A (1)Sigma+) product is significantly aligned, suggesting that the reaction mechanism is predominantly direct. The product rotational alignment was found to depend markedly on the emission wavelength, which revealed a negative correlation with the BaO(A (1)Sigma+) product vibrational state. On the basis of experimental and theoretical investigations on the reactions of N(2)O with both the (1)S0, (3)D, and (1)P1 states of Ba and the lighter group 2 atoms, it is suggested that the Ba((3)P) reaction involves a charge transfer at relatively short reagent separations and that restricted collision geometries at the highest velocity components of the broad distribution are necessary to rationalize the data.     

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).     

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.     

Spin-orbit coupling in O2(upsilon)+O2 collisions: I. Electronic structure calculations on dimer states involving the X 3Sigmag-, a 1Deltag, and b 1Sigmag+ states of O2

The importance of vibrational-to-electronic (V-E) energy transfer mediated by spin-orbit coupling in the collisional removal of O2(X 3Sigmag-,upsilon>or=26) by O2 has been reported in a recent communication [F. Dayou, J. Campos-Martinez, M. I. Hernandez, and R. Hernandez-Lamoneda, J. Chem. Phys. 120, 10355 (2004)]. The present work provides details on the electronic properties of the dimer (O2)2 relevant to the self-relaxation of O2(X 3Sigmag-,upsilon>0) where V-E energy transfer involving the O2(a 1Deltag) and O2(b 1Sigmag+) states is incorporated. Two-dimensional electronic structure calculations based on highly correlated ab initio methods have been carried out for the potential-energy and spin-orbit coupling surfaces associated with the ground singlet and two low-lying excited triplet states of the dimer dissociating into O2(X 3Sigmag-)+O2(X 3Sigmag-), O2(a 1Deltag)+O2(X 3Sigmag-), and O2(b 1Sigmag+)+O2(X 3Sigmag-). The resulting interaction potentials for the two excited triplet states display very similar features along the intermolecular separation, whereas differences arise with the ground singlet state for which the spin-exchange interaction produces a shorter equilibrium distance and higher binding energy. The vibrational dependence is qualitatively similar for the three studied interaction potentials. The spin-orbit coupling between the ground and second excited states is already nonzero in the O2+O2 dissociation limit and keeps its asymptotic value up to relatively short intermolecular separations, where the coupling increases for intramolecular distances close to the equilibrium of the isolated diatom. On the other hand, state mixing between the two excited triplet states leads to a noticeable collision-induced spin-orbit coupling between the ground and first excited states. The results are discussed in terms of specific features of the dimer electronic structure (including a simple four-electron model) and compared with existing theoretical and experimental data. This work gives theoretical insight into the origin of electronic energy-transfer mechanisms in O2+O2 collisions.     

Evidence for the production of N4 via the N2 A3Sigma(u)+ + N2 A3Sigma(u)+ energy pooling reaction

We report mass spectrometric evidence supporting our proposed mechanistic pathway for the production of N4 through the energy pooling reaction N2 A3Sigma(u)+ + N2 A3Sigma(u)+. N2 A3Sigma(u)+ is generated from the quenching of resonantly excited xenon in a mixture of xenon, 15N2, and 14N2 that is illuminated with xenon resonant lamps (147 nm). Mass spectra are periodically taken of the mixture. Over time, we observe significant isotopic scrambling of the 15N2 and 14N2, generating 15N14N in concentrations approaching 10% (approximately 2 Torr) of the initial 15N2 concentration. Though we do not observe the direct formation of N4, the isotopic ratios indicate that an excited complex (15N2(14)N2) exists long enough so that scrambling of the nitrogen atoms can occur, offering a possible route to the formation of tetrahedral nitrogen (1Td N4).     

Experimental and theoretical investigations of the reactions NH(X 3Sigma-) + D(2S)-->ND(X 3Sigma-) + H(2S) and NH(X 3Sigma-) + D(2S)-->N(4S) + HD(X 1Sigmag+)

The rate coefficient of the reaction NH(X (3)Sigma(-))+D((2)S)-->(k(1) )products (1) is determined in a quasistatic laser-flash photolysis, laser-induced fluorescence system at low pressures. The NH(X) radicals are produced by quenching of NH(a (1)Delta) (obtained in the photolysis of HN(3)) with Xe and the D atoms are generated in a D(2)/He microwave discharge. The NH(X) concentration profile is measured in the presence of a large excess of D atoms. The room-temperature rate coefficient is determined to be k(1)=(3.9+/-1.5) x 10(13) cm(3) mol(-1) s(-1). The rate coefficient k(1) is the sum of the two rate coefficients, k(1a) and k(1b), which correspond to the reactions NH(X (3)Sigma(-))+D((2)S)-->(k(1a) )ND(X (3)Sigma(-))+H((2)S) (1a) and NH(X (3)Sigma(-))+D((2)S)-->(k(1b) )N((4)S)+HD(X (1)Sigma(g) (+)) (1b), respectively. The first reaction proceeds via the (2)A(") ground state of NH(2) whereas the second one proceeds in the (4)A(") state. A global potential energy surface is constructed for the (2)A(") state using the internally contracted multireference configuration interaction method and the augmented correlation consistent polarized valence quadrupte zeta atomic basis. This potential energy surface is used in classical trajectory calculations to determine k(1a). Similar trajectory calculations are performed for reaction (1b) employing a previously calculated potential for the (4)A(") state. The calculated room-temperature rate coefficient is k(1)=4.1 x 10(13) cm(3) mol(-1) s(-1) with k(1a)=4.0 x 10(13) cm(3) mol(-1) s(-1) and k(1b)=9.1 x 10(11) cm(3) mol(-1) s(-1). The theoretically determined k(1) shows a very weak positive temperature dependence in the range 250< or =TK< or =1000. Despite the deep potential well, the exchange reaction on the (2)A(") ground-state potential energy surface is not statistical.     

Experimental and theoretical investigations of the inelastic and reactive scattering dynamics of O(3p) + D2

This paper presents a combined experimental and theoretical study of the dynamics of O((3)P) + D(2) collisions, with emphasis on a center-of-mass (c.m.) collision energy of 25 kcal mol(-1). The experiments were conducted with a crossed-molecular-beams apparatus, employing a laser detonation source to produce hyperthermal atomic oxygen and mass spectrometric detection to measure the product angular and time-of-flight distributions. The novel beam source, which enabled these experiments to be conducted, contributed unique challenges to the experiments and to the analysis, so the experimental methods and approach to the analysis are discussed in detail. Three different levels of theory were used: (1) quasiclassical trajectories (QCT), (2) time-independent quantum scattering calculations based on high-quality potential surfaces for the two lower-energy triplet states, and (3) trajectory-surface-hopping (TSH) studies that couple the triplet surfaces with the lowest singlet surface using a spin-orbit Hamiltonian derived from ab-initio calculations. The latter calculations explore the importance of intersystem crossing in the dynamics. Both experiment and theory show that inelastically scattered O atoms scatter almost exclusively in the forward direction, with little or no loss of translational energy. For the reaction, O((3)P) + D(2) --> OD + D, the experiment shows that, on average, approximately 50% of the available energy goes into product translation and that the OD product angular distributions are largely backward-peaked. These results may be interpreted in light of the QCT and TSH calculations, leading to the conclusion that the reaction occurs mainly on triplet potential energy surfaces with, at most, minor intersystem crossing to a singlet surface. Reaction on either of the two low-lying reactive triplet surfaces proceeds through a rebound mechanism in which the angular distributions are backward-peaked and the product OD is both vibrationally and rotationally excited. The quantum scattering results are in good agreement with QCT calculations, indicating that quantum effects are relatively small for this reaction at a collision energy of 25 kcal mol(-1).     

Repulsive double many-body expansion potential energy surface for the reactions N(4S)+H2<-->NH(X3Sigma-)+H from accurate ab initio calculations

A single-sheeted DMBE potential energy surface is reported for the reactions N(4S)+H2<-->NH(X3Sigma-)+H based on a fit to accurate multireference configuration interaction energies. These have been calculated using the aug-cc-pVQZ basis set of Dunning and the full valence complete active space wave function as reference, being semi-empirically corrected by scaling the two-body and three-body dynamical correlation energies. The topographical features of the novel global potential energy surface are examined in detail, including a conical intersection involving the two first 4A' potential energy surfaces which has been transformed into an avoided crossing in the present single-sheeted representation.     

Potential energy surface and product branching ratios for the reaction of dicarbon, C(2)(X(1)Sigma(g)(+)), with methylacetylene, CH(3)CCH(X(1)A(1)): an ab initio/RRKM study

Ab initio calculations of the potential energy surface for the C(2)(X(1)Sigma(g)(+)) + CH(3)CCH(X(1)A(1)) reaction have been carried at the G2M level of theory. The calculations show that the dicarbon molecule in the ground singlet electronic state can add to methylacetylene without a barrier producing a three-member or a four-member ring intermediate, which can rapidly rearrange to the most stable H(3)CCCCCH isomer on the C(5)H(4) singlet surface. This isomer can then lose a hydrogen atom (H) or molecular hydrogen (H(2)) from the CH(3) group with the formation of H(2)CCCCCH and HCCCCCH, respectively. Alternatively, H atom migrations and three-member-ring closure/opening rearrangements followed by H and H(2) losses can lead to other isomers of the C(5)H(3) and C(5)H(2) species. According to the calculated energetics, the C(2)(X(1)Sigma(g)(+)) + CH(3)CCH reaction is likely to be a major source of the C(5)H(3) radicals (in particular, the most stable H(2)CCCCCH and HCCCHCCH isomers, which are relevant to the formation of benzene through the reactions with CH(3)). Among heavy-fragment product channels, only C(3)H(3) + C(2)H and c-C(3)H(2) + C(2)H(2) might compete with C(5)H(3) + H and C(5)H(2) + H(2). RRKM calculations of reaction rate constants and product branching ratios depending on the reactive collision energy showed that the major reaction products are expected to be H(2)CCCCCH + H (64-66%) and HCCCHCCH + H (34-30%), with minor contributions from HCCCCCH + H(2) (1-2%), HCCCHCC + H(2) (up to 1%), C(3)H(3) + C(2)H (up to 1%), and c-C(3)H(2) + C(2)H(2) (up to 0.1%) if the energy randomization is complete. The calculations also indicate that the C(2)(X(1)Sigma(g)(+)) + CH(3)CCH(X(1)A(1)) reaction can proceed by direct H-abstraction of a methyl hydrogen to form C(3)H(3) + C(2)H almost without a barrier.     

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.