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.     

Exploring the dynamics of hydrogen atom release from the radical-radical reaction of O(3P) with C3H5

The gas-phase radical-radical reaction dynamics of O(3P) + C3H5 --> H(2S) + C3H4O was studied at an average collision energy of 6.4 kcal/mol in a crossed beam configuration. The ground-state atomic oxygen [O(3P)] and allyl radicals (C3H5) were generated by the photolysis of NO2 and the supersonic flash pyrolysis of allyl iodide, respectively. Nascent hydrogen atom products were probed by the vacuum-ultraviolet-laser induced fluorescence spectroscopy in the Lyman-alpha region centered at 121.6 nm. With the aid of the CBS-QB3 level of ab initio theory, it has been found that the barrierless addition of O(3P) to C3H5 forms the energy-rich addition complexes on the lowest doublet potential energy surface, which are predicted to undergo a subsequent direct decomposition step leading to the reaction products H + C3H4O. The major counterpart C3H4O of the probed hydrogen atom is calculated to be acrolein after taking into account the factors of barrier height, reaction enthalpy, and the number of intermediates involved along the reaction pathway. The nascent H-atom Doppler profile analysis shows that the average center-of-mass translational energy of the H + C3H4O products and the fraction of the total available energy released as the translational energy were determined to be 3.83 kcal/mol and 0.054, respectively. On the basis of comparison with statistical calculations, the reaction proceeds through the formation of short-lived addition complexes rather than statistical, long-lived intermediates, and the polyatomic acrolein product is significantly internally excited at the moment of the decomposition.     

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.     

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.     

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.     

Radical-molecule reaction C3H+H2O: a mechanistic study

Despite the importance of the C(3)H radical in both combustion and interstellar space, the reactions of C(3)H toward stable molecules have never been studied. In this paper, we report our detailed mechanistic study on the radical-molecule reaction C(3)H+H(2)O at the Becke's three parameter Lee-Yang-Parr-B3LYP6-311G(d,p) and coupled cluster with single, double, and triple excitations-CCSD(T)6-311G(2d,p) (single-point) levels. It is shown that the C(3)H+H(2)O reaction initially favors formation of the carbene-insertion intermediates HCCCHOH (1a,1b) rather than the direct H- or OH-abstraction process. Subsequently, the isomers (1a,1b) can undergo a direct H- extrusion to form the well-known product propynal HCCCHO (P(5)). Highly competitively, (1a,1b) can take the successive 1,4- and 1,2-H-shift interconversion to isomer H(2)CCCHO(2a,2b) and then to isomer H(2)CCHCO(3a,3b), which can finally take a direct C-C bond cleavage to give product C(2)H(3) and CO (P(1)). The other products are kinetically much less feasible. With the overall entrance barrier 10.6 kcal/mol, the title reaction can be important in postburning processes. Particularly, our calculations suggest that the title reaction may play a role in the formation of the intriguing interstellar molecule, propynal HCCCHO. The calculated results will also be useful for the analogous C(3)H reactions such as with ammonia and alkanes.     

C2H+H2CO: a new route for formaldehyde removal

The title unknown reaction is theoretically studied at various levels to probe the interaction mechanism between the ethynyl radical (HC triple bond C) and formaldehyde (H(2)C double bond O). The most feasible pathway is a barrier-free direct H-abstraction process leading to acetylene and formyl radical (C(2)H(2)+HCO) via a weakly bound complex, and then the product can take secondary dissociation to the final product C(2)H(2)+CO+H. The C-addition channel leading to propynal plus H-atom (HCCCHO+H) has the barrier of only 3.6, 2.9, and 2.1 kcal/mol at the CCSD(T)/6-311+G(3df,2p)MP2//6-311G(d,p)+ZPVE, CCSD(T)/6-311+G(3df,2p)//QCISD/6-311G(d,p)+ZPVE, and G3//MP2 levels, respectively [CCSD(T)--coupled cluster with single, double, and triple excitations; ZPVE--zero-point vibrational energy; QCISD--quadratic configuration interaction with single and double excitations; G3//MP2-Gaussian-3 based on Moller-Plesset geometry]. The O addition also leading to propynal plus H atom needs to overcome a higher barrier of 5.3, 8.7, and 3.0 kcalmol at the three corresponding levels. The title no-barrier reaction presents a new efficient route to remove the pollutant H(2)CO, and should be included in the combustion models of hydrocarbons. It may also represent the fastest radical-H(2)CO reaction among the available theoretical data. Moreover, it could play an important role in the interstellar chemistry where the zero- or minute-barrier reactions are generally favored. Discussions are also made on the possible formation of the intriguing propynal in space via the title reaction on ice surface.     

Prediction of product branching ratios in the C(3P)+C2H2-->l-C3H+H/c-C3H+H/C3+H2 reaction using ab initio coupled clusters calculations extrapolated to the complete basis set combined with Rice-Ramsperger-Kassel-Marcus and radiationless transition theories

Ab initio CCSD(T) calculations of intermediates and transition states on the singlet and triplet C3H2 potential energy surfaces extrapolated to the complete basis set limit are combined with statistical computations of energy-dependent rate constants of the C(3P)+C2H2 reaction under crossed molecular beam conditions. Rice-Ramsperger-Kassel-Marcus theory is applied for isomerization and dissociation steps within the same multiplicity and radiationless transition and nonadiabatic transition state theories are used for singlet-triplet intersystem crossing rates. The calculated rate constants are utilized to predict product branching ratios. The results demonstrate that, in qualitative agreement with available experimental data, c-C3H+H and C3+H2 are the most probable products at low collision energies, whereas l-C3H+H becomes dominant at higher Ec above approximately 25 kJ/mol.     

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.     

A study of the radical-radical reaction dynamics of O(3P)+t-C4H9 --> OH+iso-C4H8

The radical-radical reaction dynamics of ground-state atomic oxygen [O(3P)] with t-butyl radicals (t-C4H9) in the gas phase were investigated using high-resolution laser spectroscopy in a crossed-beam configuration, together with ab initio theoretical calculations. The radical reactants, O(3P) and t-C4H9, were produced by the photodissociation of NO2 and the supersonic flash pyrolysis of the precursor, azo-t-butane, respectively. A new exothermic channel, O(3P)+t-C4H9 --> OH+iso-C4H8, was identified and the nascent rovibrational distributions of the OH (X 2Pi: upsilon" = 0,1,2) products were examined. The population analyses for the two spin-orbit states of F1(2Pi32) and F2(2Pi12) showed that the upsilon" = 0 level is described by a bimodal feature composed of low- and high-N" rotational components, whereas the upsilon" = 1 and 2 levels exhibit unimodal distributions. No noticeable spin-orbit or Lambda-doublet propensities were observed in any vibrational state. The partitioning ratio of the vibrational populations (Pupsilon") with respect to the low-N" components of the upsilon" = 0 level was estimated to be P0:P1:P2 = 1:1.17+/-0.24:1.40+/-0.11, indicating that the nascent internal distributions are highly excited. On the basis of the comparison of the experimental results with the statistical theory, the reaction mechanism at the molecular level can be described in terms of two competing dynamic pathways: the major, direct abstraction process leading to the inversion of the vibrational populations, and the minor, short-lived addition-complex process responsible for the hot rotational distributions. After considering the reaction exothermicity, the barrier height, and the number of intermediates along the addition reaction pathways on the lowest doublet potential energy surface, the formation of CH3COCH3(acetone)+CH3 was predicted to be dominant in the addition mechanism.     

Singlet and triplet potential surfaces for the O2+C2H4 reaction

Electronic structure calculations at the CASSCF and UB3LYP levels of theory with the aug-cc-pVDZ basis set were used to characterize structures, vibrational frequencies, and energies for stationary points on the ground state triplet and singlet O(2)+C(2)H(4) potential energy surfaces (PESs). Spin-orbit couplings between the PESs were calculated using state averaged CASSCF wave functions. More accurate energies were obtained for the CASSCF structures with the MRMP2/aug-cc-pVDZ method. An important and necessary aspect of the calculations was the need to use different CASSCF active spaces for the different reaction paths on the investigated PESs. The CASSCF calculations focused on O(2)+C(2)H(4) addition to form the C(2)H(4)O(2) biradical on the triplet and singlet surfaces, and isomerization reaction paths ensuing from this biradical. The triplet and singlet C(2)H(4)O(2) biradicals are very similar in structure, primarily differing in their C-C-O-O dihedral angles. The MRMP2 values for the O(2)+C(2)H(4)→C(2)H(4)O(2) barrier to form the biradical are 33.8 and 6.1 kcal/mol, respectively, for the triplet and singlet surfaces. On the singlet surface, C(2)H(4)O(2) isomerizes to dioxetane and ethane-peroxide with MRMP2 barriers of 7.8 and 21.3 kcal/mol. A more exhaustive search of reaction paths was made for the singlet surface using the UB3LYP/aug-cc-pVDZ theory. The triplet and singlet surfaces cross between the structures for the O(2)+C(2)H(4) addition transition states and the biradical intermediates. Trapping in the triplet biradical intermediate, following (3)O(2)+C(2)H(4) addition, is expected to enhance triplet→singlet intersystem crossing.     

Crossed beams and theoretical studies of the O(3P)+CH4-->H+OCH3 reaction excitation function

The excitation function for the reaction, O(3P)+CH4-->H+OCH3, has been measured in a crossed molecular beams experiment and determined with direct dynamics calculations that use the quasiclassical trajectory method in conjunction with a recently developed semiempirical Hamiltonian. Good agreement is found between experiment and theory, enabling us to address two fundamental issues for the O(3P)+CH4 reaction that arise for all O(3P)+saturated hydrocarbon reactions: (1) the importance of triplet excited states that correlate adiabatically to ground-state reactants and products and (2) the importance of intersystem crossing processes involving the lowest singlet surface [corresponding to reaction with O(1D)]. Our results indicate that the first excited triplet surface contributes substantially to the cross section when the collision energy exceeds the reaction barrier (approximately 2 eV) by more than 0.5 eV. Although triplet-singlet crossings may occur at all energies, we have found that their effect on the excitation function is negligible for the collision energies studied-up to 1.5 eV above threshold.     

Reaction dynamics of OH+(3Sigma-)+C2H2 studied with crossed beams and density functional theory calculations

The reactions between OH+(3Sigma-) and C2H2 have been studied using crossed ion and molecular beams and density functional theory calculations. Both charge transfer and proton transfer channels are observed. Products formed by carbon-carbon bond cleavage analogous to those formed in the isoelectronic O(3P)+C2H2 reaction, e.g., 3CH2 + HCO+, are not observed. The center of mass flux distributions of both product ions at three different energies are highly asymmetric, with maxima close to the velocity and direction of the precursor acetylene beam, characteristic of direct reactions. The internal energy distributions of the charge transfer products are independent of collision energy and are peaked at the reaction exothermicity, inconsistent with either the existence of favorable Franck-Condon factors or energy resonance. In proton transfer, almost the entire reaction exothermicity is transformed into product internal excitation, consistent with mixed energy release in which the proton is transferred with both the breaking and forming bonds extended. Most of the incremental translational energy in the two higher-energy experiments appears in product translational energy, providing an example of induced repulsive energy release.     

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.     

Quasiclassical trajectory calculations of the reaction C+C2H2-->l-C3H, c-C3H+H, C3+H2 using full-dimensional triplet and singlet potential energy surfaces

Full-dimensional, density functional theory (B3LYP/6-311g(d,p))-based potential energy surfaces (PESs) are reported and used in quasi-classical calculations of the reaction of C with C(2)H(2). For the triplet case, the PES spans the region of the reactants, the complex region (with numerous minima and saddle points) and the products, linear(l)-C(3)H+H, cyclic(c)-C(3)H+H and c-(3)C(3)+H(2). For the singlet case, the PES describes the complex region and products l-C(3)H+H, c-C(3)H+H and l-(1)C(3)+H(2). The PESs are invariant under permutation of like nuclei and are fit to tens of thousands of electronic energies. Energies and harmonic frequencies of the PESs agree well the DFT ones for all stationary points and for the reactant and the products. Dynamics calculations on the triplet PES find both l-C(3)H and c-C(3)H products, with l-C(3)H being dominant at the energies considered. Limited unimolecular reaction dynamics on the singlet PES find both products in comparable amounts as well as the C(3)+H(2) product.     

B(2)A(')-X(2)A(') detection of vibrationally excited HCO produced by the O((3)P)+C(2)H(4) reaction

The distribution of rotational and vibrational energy in HCO produced by the O((3)P)+C(2)H(4) reaction has been measured using laser-induced fluorescence detection via the B(2)A(')-X(2)A(') transition. Over a detection wavelength range of 248-290 nm, our experiments have shown that HCO is formed in both the ground state and in at least six vibrationally excited states with up to two quanta of energy in the C-O stretch and the bending mode. Dispersed fluorescence experiments were conducted to positively assign all of the HCO vibrational bands. The experiments confirmed that many bands, including the B(000)-X(000) band, are affected by overlap with other HCO bands. Spectral modeling was used to separate the contributions of overlapping HCO B-X bands and to determine a nascent HCO rotational temperature of approximately 600 K, corresponding to approximately 6% of the total energy from the O((3)P)+C(2)H(4) reaction. HCO vibrational distributions were determined for two different average collision energies and were fit with vibrational temperatures of 1850+/-80 K and 2000+/-100 K, corresponding to approximately 15% of the total energy. The observed Boltzmann distribution of vibrational energy in HCO indicates that HCO and CH(3) are formed by the dissociation of an energized intermediate complex.     

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.     

Products of the gas-phase reactions of OH radicals with (C2H5O)2P(S)CH3 and (C2H5O)3PS

Products of the gas-phase reactions of OH radicals with O,O-diethyl methylphosphonothioate [(C2H5O)2P(S)CH3, DEMPT] and O,O,O-triethyl phosphorothioate [(C2H5O)3PS, TEPT] have been investigated at room temperature and atmospheric pressure of air using in situ atmospheric pressure ionization mass spectrometry (API-MS) and, for the TEPT reaction, gas chromatography and in situ Fourier transform infrared (FT-IR) spectroscopy. Combined with products quantified previously by gas chromatography, the products observed were: from the DEMPT reaction, (C2H5O)2P(O)CH3 (21+/-4% yield) and C2H5OP(S)(CH3)OH or C2H5OP(O)(CH3)SH (presumed to be C2H5OP(O)(CH3)SH by analogy with the TEPT reaction); and from the TEPT reaction, (C2H5O)3PO (54-62% yield), SO2 (67+/-10% yield), CH3CHO (22-40% yield) and, tentatively, (C2H5O)2P(O)SH. The FT-IR analyses showed that the formation yields of HCHO, CO, CO2, peroxyacetyl nitrate [CH3C(O)OONO2], organic nitrates, and acetates from the TEPT reaction were <5%, 3+/-1%, <7%, <2%, 5+/-3%, and 3+/-2%, respectively. Possible reaction mechanisms are discussed.     

Ab initio Study of the Potential Energy Surface and Product Branching Ratios for Reaction of O（^1D） with C2H5CI

The potential energy surface of O（^1D）＋C2H5Cl reaction was studied using QCISD（T）/6- 311＋＋G（d,p）//MP2/6-31G（d,p） method. The calculations reveal an insertion-elimination mechanism. The insertion reaction of O（^1D） and C2H5Cl produces two energy-rich intermediates, IM1 and IM2, which subsequently decompose into various products. The calculations of the branching ratios of various products formed through the two intermediates were carried out using RRKM （Rice-Ramsperger-Kassel-Marcus） theory at the collision energies of 0, 20.9, 41.8, 62.7, 83.6, 104.5, and 125.4 kJ/mol. HCl is the main decomposition product for IM1; CH2OH is the main decomposition product for IM2. Since IM1 is more stable than IM2, HCl is probably the main product of the O（^1D）＋C2H5Cl reaction.     

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

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