Product channels of the HCCO + NO reaction

The product branching ratio of the HCCO + NO reaction was investigated using the laser photolysis/infrared absorption technique. Ethyl ethynyl ether (C(2)H(5)OCCH) was used as the HCCO radical precursor. Transient infrared detection of CO, CO(2), and HCNO products was used to determine the following branching ratios at 296 K: phi(CO+HCNO) = 0.78 +/- 0.04 and phi(CO(2)+HCN) = 0.22 +/- 0.04. These values are in good agreement with some recent ab initio calculations.     

Direct ab initio dynamics study of the OH + HOCO reaction

The reaction between OH and HOCO has been examined using the coupled-cluster method to locate and optimize the critical points on the ground-state potential energy surface. The energetics are refined using the coupled-cluster method with basis set extrapolation to the complete basis set (CBS) limit. Results show that the OH + HOCO reaction produces H2O + CO2 as final products and the reaction passes through an HOC(O)OH intermediate. In addition, the OH + HOCO reaction has been studied using a direct dynamics method with a dual-level ab initio theory. Dynamics calculations show that hydrogen bonding plays an important role during the initial stages of the reaction. The thermal rate constant is estimated over the temperature range 250-800 K. The OH + HOCO reaction is found to be nearly temperature-independent at lower temperatures, and at 300 K, the thermal rate constant is predicted to be 1.03 x 10(-11) cm3 molecule(-1) s(-1). In addition, there may be an indication of a small peak in the rate constant at a temperature between 300 and 400 K.     

Ab initio study on the oxidation of NCN by O (3P): prediction of the total rate constant and product branching ratios

The reaction of NCN with O is relevant to the formation of prompt NO according to the new mechanism, CH+N2-->cyclic-C(H)NN- -->HNCN-->H+NCN. The reaction has been investigated by ab initio molecular orbital and transition state theory calculations. The mechanisms for formation of possible product channels involved in the singlet and triplet potential energy surfaces have been predicted at the highest level of the modified GAUSSIAN-2 (G2M) method, G2M (CC1). The barrierless association/dissociation processes on the singlet surface were also examined with the third-order Rayleigh-Schr?dinger perturbation (CASPT3) and the multireference configuration interaction methods including Davidson's correction for higher excitations (MRCI+Q) at the CASPT3(6,6)/6-311+G(3df)//UB3LYP/6-311G(d) and MRCI+Q(6,6)/6-311+G(3df)//UB3LYP/6-311G(d) levels. The rate constants for the low-energy channels producing CO+N2, CN+NO, and N(4S)+NCO have been calculated in the temperature range of 200-3000 K. The results show that the formation of CN+NO is dominant and its branching ratio is over 99% in the whole temperature range; no pressure dependence was noted at pressures below 100 atm. The total rate constant can be expressed by: kt=4.23x10(-11) T0.15 exp(17/T) cm3 molecule(-1) s(-1).     

Ab initio energies and product branching ratios for the O+C3H6 reaction

Intermediate and transition-state energies have been calculated for the O+C3H6 (propene) reaction using the compound ab initio CBS-QB3 and G3 methods in combination with density functional theory. The lowest-lying triplet and singlet potential energy surfaces of the O-C3H6 system were investigated. RRKM statistical theory was used to predict product branching fractions over the 300-3000 K temperature and 0.001-760 Torr pressure ranges. The oxygen atom adds to the C3H6 terminal olefinic carbon in the primary step to generate a nascent triplet biradical, CH3CHCH2O. On the triplet surface, unimolecular dissociation of CH3CHCH2O to yield H+CH3CHCHO is favored over the entire temperature range, although the competing H2CO+CH3CH product channel becomes significant at high temperature. Rearrangement of triplet CH3CHCH2O to CH3CH2CHO (propanal) via a 1,2 H-atom shift has a barrier of 122.3 kJ mol(-1), largely blocking this reaction channel and any subsequent dissociation products. Intersystem crossing of triplet CH3CHCH2O to the singlet surface, however, leads to facile rearrangement to singlet CH3CH2CHO, which dissociates via numerous product channels. Pressure was found to have little influence over the branching ratios under most conditions, suggesting that the vibrational self-relaxation rates for p相似文献   

Multireference configuration interaction calculations for the F(2P)+HCl-->HF+Cl(2P) reaction: a correlation scaled ground state (1 2A') potential energy surface

This paper presents a new ground state (1 (2)A(')) electronic potential energy surface for the F((2)P)+HCl-->HF+Cl((2)P) reaction. The ab initio calculations are done at the multireference configuration interaction+Davidson correction (MRCI+Q) level of theory by complete basis set extrapolation of the aug-cc-pVnZ (n=2,3,4) energies. Due to low-lying charge transfer states in the transition state region, the molecular orbitals are obtained by six-state dynamically weighted multichannel self-consistent field methods. Additional perturbative refinement of the energies is achieved by implementing simple one-parameter correlation energy scaling to reproduce the experimental exothermicity (DeltaE=-33.06 kcalmol) for the reaction. Ab initio points are fitted to an analytical function based on sum of two- and three-body contributions, yielding a rms deviation of <0.3 kcalmol for all geometries below 10 kcalmol above the barrier. Of particular relevance to nonadiabatic dynamics, the calculations show significant multireference character in the transition state region, which is located 3.8 kcalmol with respect to F+HCl reactants and features a strongly bent F-H-Cl transition state geometry (theta approximately 123.5 degrees ). Finally, the surface also exhibits two conical intersection seams that are energetically accessible at low collision energies. These seams arise naturally from allowed crossings in the C(infinityv) linear configuration that become avoided in C(s) bent configurations of both the reactant and product, and should be a hallmark of all X-H-Y atom transfer reaction dynamics between ((2)P) halogen atoms.     

Ab initio chemical kinetics for the OH+HNCN reaction

The kinetics and mechanism of the reaction of the cyanomidyl radical (HNCN) with the hydroxyl radical (OH) have been investigated by ab initio calculations with rate constants prediction. The single and triplet potential energy surfaces of this reaction have been calculated by single-point calculations at the CCSD(T)/6-311+G(3df,2p) level based on geometries optimized at the B3LYP/6-311+G(3df,2p) and CCSD/6-311++G(d,p) levels. The rate constants for various product channels in the temperature range of 300-3000 K are predicted by variational transition-state and Rice-Ramsperger-Kassel-Marcus (RRKM) theories. The predicted total rate constants can be represented by the expressions ktotal=2.66 x 10(+2)xT-4.50 exp(-239/T) in which T=300-1000 K and 1.38x10(-20)xT2.78 exp(1578/T) cm3 molecule(-1) s(-1) where T=1000-3000 K. The branching ratios of primary channels are predicted: k1 for forming singlet HON(H)CN accounts for 0.32-0.28, and k4 for forming singlet HONCNH accounts for 0.68-0.17 in the temperature range of 300-800 K. k2+k7 for producing H2O+NCN accounts for 0.55-0.99 in the high-temperature range of 800-3000 K. The branching ratios of k3 for producing HCN+HNO, k6 for producing H2N+NCO, k8 for forming 3HN(OH)CN, k9 for producing CNOH+3NH, and k5+k10 for producing NH2+NCO are negligible. The rate constants for key individual product channels are provided in a table for different temperature and pressure conditions.     

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.     

Dynamics on the HOCO potential energy surface studied by dissociative photodetachment of HOCO- and DOCO-

An experimental study of the dissociative photodetachment (DPD) dynamics of HOCO(-) and DOCO(-) at a photon energy of 3.21 eV has been carried out to probe the potential energy surface of the HOCO free radical and the dynamics of the OH+CO-->H+CO(2) reaction. These photoelectron-photofragment coincidence experiments allow the identification of photodetachment processes leading to the production of stable HOCO free radicals and both the H+CO(2) and OH+CO dissociation channels on the neutral surface. Isotopic substitution by deuterium in the parent ion is observed to reduce the product branching ratio for the D+CO(2) channel, consistent with tunneling playing a role in this dissociation pathway. Other isotope effects on the detailed partitioning of kinetic energy between photoelectrons and photofragments are also discussed. The results are compared to recent theoretical predictions of this DPD process, and evidence for the involvement of vibrationally excited HOCO(-) anions is discussed.     

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.     

The reaction of phenyl radical with molecular oxygen: a G2M study of the potential energy surface

Ab initio G2M calculations have been performed to investigate the potential energy surface for the reaction of C6H5 with O2. The reaction is shown to start with an exothermic barrierless addition of O2 to the radical site of C6H5 to produce phenylperoxy (1) and, possibly, 1,2-dioxaspiro[2.5]octadienyl (dioxiranyl, 8) radicals. Next, 1 loses the terminal oxygen atom to yield the phenoxy + O products (3) or rearranges to 8. The dioxiranyl can further isomerize to a seven-member ring 2-oxepinyloxy radical (10), which can give rise to various products including C5H5 + CO2, pyranyl + CO, o-benzoquinone + H, and 2-oxo-2,3-dihydrofuran-4-yl + C2H2. Once 10 is produced, it is unlikely to go back to 8 and 1, because the barriers separating 10 from the products are much lower than the reverse barrier from 10 to 8. Thus, the branching ratio of C6H5O + O against the other products is mostly controlled by the critical transition states between 1 and 3, 1 and 8, and 8 and 10. According to the calculated barriers, the most favorable product channel for the decomposition of 10 is C5H5 + CO2, followed by pyranyl + CO and o-benzoquinone + H. Since C6H5O + O and C5H5 + CO2 are expected to be the major primary products of the C6H5 + O2 reaction and thermal decomposition of C6H5O leads to C5H5 + CO, cyclopentadienyl radicals are likely to be the major product of phenyl radical oxidation, and so it results in degradation of the six-member aromatic ring to the five-member cyclopentadienyl ring. Future multichannel RRKM calculations of reaction rate constants are required to support these conclusions and to quantify the product branching ratios at various combustion conditions.     

Direct ab initio molecular dynamics study on a microsolvated SN2 reaction of OH-(H2O) with CH3Cl

Reaction dynamics for a microsolvated SN2 reaction OH-(H2O)+CH3Cl have been investigated by means of the direct ab initio molecular dynamics method. The relative center-of-mass collision energies were chosen as 10, 15, and 25 kcal/mol. Three reaction channels were found as products. These are (1) a channel leading to complete dissociation (the products are CH3OH+Cl- +H2O: denoted by channel I), (2) a solvation channel (the products are Cl-(H2O)+CH3OH: channel II), and (3) a complex formation channel (the products are CH3OH...H2O+Cl-: channel III). The branching ratios for the three channels were drastically changed as a function of center-of-mass collision energy. The ratio of complete dissociation channel (channel I) increased with increasing collision energy, whereas that of channel III decreased. The solvation channel (channel II) was minor at all collision energies. The selectivity of the reaction channels and the mechanism are discussed on the basis of the theoretical results.     

Theory of multichannel thermal unimolecular reactions. 2. Application to the thermal dissociation of formaldehyde

The thermal dissociation of formaldehyde proceeds on three channels, the molecular-elimination channel H2CO --> H2 + CO (1), the radical-forming bond-fission channel H2CO --> H + HCO (2), and the bond-fission-initiated, intramolecular-hydrogen-abstraction channel H2CO --> H...HCO --> H2 + CO (3) which also forms molecular products. The kinetics of this system in the low-pressure range of the unimolecular reaction is shown to be governed by a subtle superposition of collisional channel coupling to be treated by solving a master equation, of rotational channel switching accessible through ab initio calculations of the potential as well as spectroscopic and photophysical determinations of the threshold energies and channel branching above the threshold energy for radical formation which can be characterized through formaldehyde photolysis quantum yields as well as classical trajectory calculations. On the basis of the available information, the rate coefficients for the formation of molecular and radical fragments are analyzed and extrapolated over wide ranges of conditions. The modeled rate coefficients in the low-pressure range of the reaction (neglecting tunneling) over the range 1400-3200 K in the bath-gas Ar in this way are represented by k0,Mol/[Ar] approximately 9.4 x 10(-9) exp(-33,140 K/T) cm3 molecule(-1) s(-1) and k0,Rad/[Ar] approximately 6.2 x 10(-9) exp(-36,980 K/T) cm3 molecule(-1) s(-1). The corresponding values for the bath-gas Kr, on which the analysis relies in particular, are k0,Mol/[Kr] approximately 7.7 x 10(-9) exp(-33,110 K/T) and k0,Rad/[Kr] approximately 4.1 x 10(-9) exp(-36 910 K/T) cm3 molecule(-1) s(-1). While the threshold energy E0,2 for channels 2 and 3 is taken from spectroscopic measurements, the threshold energy E0,1 for channel 1 is fitted on the basis of experimental ratios k0,Rad/k0,Mol in combination with photolysis quantum yields. The derived value of E0,1(1) = 81.2 (+/-0.9) kcal mol(-1) is in good agreement with results from recent ab initio calculations, 81.9 (+/-0.3) kcal mol(-1), but is higher than earlier results derived from photophysical experiments, 79.2 (+/-0.8) kcal mol(-1). Rate coefficients for the high-pressure limit of the reaction are also modeled. The results of the present work markedly depend on the branching ratio between channels 2 and 3. Expressions of this branching ratio from classical trajectory calculations and from photolysis quantum yield measurements were tested. At the same time, a modeling of the photolysis quantum yields was performed. The formaldehyde system so far presents the best characterized multichannel dissociation reaction. It may serve as a prototype for other multichannel dissociation reactions.     

Quantum dynamics study of the dissociative photodetachment of HOCO-

Six-dimensional wave packet calculations are carried out to study the behavior of HOCO subsequent to the photodetachment of an electron from the negative anion, HOCO-. It is possible to form stable and/or long-lived HOCO complexes, as well as the dissociative products OH+CO and H+CO2. A variety of observables are determined: the electron kinetic energy (eKE) distributions associated with the OH+CO and H+CO2 channels, the correlated eKE and product translational energy distribution for the OH+CO channel, and product branching ratios. Most of our results are in good accord with the experimental results of Clements, Continetti, and Francisco [J. Chem. Phys. 117, 6478 (2002)], except that the calculated eKE distribution for the H+CO2 channel is noticeably colder than experiment. Reasons for this discrepancy are suggested.     

Quantum chemical and theoretical kinetics study of the O(3P) + C2H2 reaction: a multistate process

The potential energy surfaces of the two lowest-lying triplet electronic surfaces 3A' and 3A' for the O(3P) + C2H2 reaction were theoretically reinvestigated, using various quantum chemical methods including CCSD(T), QCISD, CBS-QCI/APNO, CBS-QB3, G2M(CC,MP2), DFT-B3LYP and CASSCF. An efficient reaction pathway on the electronically excited 3A' surface resulting in H(2S) + HCCO(A2A') was newly identified and is predicted to play an important role at higher temperatures. The primary product distribution for the multistate multiwell reaction was then determined by RRKM statistical rate theory and weak-collision master equation analysis using the exact stochastic simulation method. Allowing for nonstatistical behavior of the internal rotation mode of the initial 3A' adducts, our computed primary-product distributions agree well with the available experimental results, i.e., ca. 80% H(2S) + HCCO(X2A' + A2A') and 20% CH2(X3B1) + CO(X1sigma+) independent of temperature and pressure over the wide 300-2000 K and 0-10 atm ranges. The thermal rate coefficient k(O + C2H2) at 200-2000 K was computed using multistate transition state theory: k(T) = 6.14 x 10(-15)T (1.28) exp(-1244 K/T) cm3 molecule(-1) s(-1); this expression, obtained after reducing the CBS-QCI/APNO ab initio entrance barriers by 0.5 kcal/mol, quasi-perfectly matches the experimental k(T) data over the entire 200-2000 K range, spanning 3 orders of magnitude.     

Quantum force molecular dynamics study of the reaction of O atoms with HOCO

The reaction of HOCO with O atoms has been studied using a direct ab initio dynamics approach based on the scaling all correlation UCCD/D95(d,p) method. Ab initio calculations point to two possible reaction mechanisms for the O+HOCO-->OH+CO2 reaction. They are a direct hydrogen abstraction and an oxygen addition reaction through a short-lived HOC(O)O intermediate. The dynamics results show that only the addition mechanism is important under the conditions considered here. The lifetime of the HOC(O)O complex is predicted to be 172+/-15 fs. This is typical of a direct and fast radical-radical reaction. At room temperature, the calculated thermal rate coefficient is 1.44 x 10(-11) cm(3) mol(-1) s(-1) and its temperature dependence is rather weak. The two kinds of reactive trajectories are illustrated in detail.     

Experimental and theoretical studies of rate coefficients for the reaction O(3P)+C2H5OH at high temperatures

Rate coefficients of the reaction O(3P)+C2H5OH in the temperature range 782-1410 K were determined using a diaphragmless shock tube. O atoms were generated by photolysis of SO2 at 193 nm with an ArF excimer laser; their concentrations were monitored via atomic resonance absorption. Our data in the range 886-1410 K are new. Combined with previous measurements at low temperature, rate coefficients determined for the temperature range 297-1410 K are represented by the following equation: k(T)=(2.89+/-0.09)x10(-16)T1.62 exp[-(1210+/-90)/T] cm3 molecule(-1) s(-1); listed errors represent one standard deviation in fitting. Theoretical calculations at the CCSD(T)/6-311+G(3df, 2p)//B3LYP/6-311+G(3df) level predict potential energies of various reaction paths. Rate coefficients are predicted with the canonical variational transition state (CVT) theory with the small curvature tunneling correction (SCT) method. Reaction paths associated with trans and gauche conformations are both identified. Predicted total rate coefficients, 1.60 x 10(-22)T3.50 exp(16/T) cm3 molecule(-1) s(-1) for the range 300-3000 K, agree satisfactorily with experimental observations. The branching ratios of three accessible reaction channels forming CH3CHOH+OH (1a), CH2CH2OH+OH (1b), and CH3CH2O+OH (1c) are predicted to vary distinctively with temperature. Below 500 K, reaction 1a is the predominant path; the branching ratios of reactions 1b,c become approximately 40% and approximately 11%, respectively, at 2000 K.     

Structural isomers of C2N+: a selected-ion flow tube study

Reactivities of the structural isomers CCN+ and CNC+ were examined in a selected-ion flow tube at 300 +/- 5 K. The less reactive CNC+ isomer was identified as the product of the reactions of C(+) + HCN and C(+) + C2N2; in these reactions only CNC+ can be produced because of energy constraints. Rate coefficients and branching ratios are reported for the reactions of each isomer with H2, CH4, NH3, H2O, C2H2, HCN, N2, O2, N2O, and CO2. Ab initio calculations are presented for CCN+ and CNC+; a saddle point for the reaction CCN+ --> CNC+ is calculated to be 195 kJ mol-1 above the CNC+. The results provide evidence that the more reactive CCN+ isomer is unlikely to be present in measurable densities in interstellar clouds.     

Determination of the CH + O2 product channels

The multichannel CH + O2 reaction was studied at room temperature, in a low-pressure fast-flow reactor. CH radical was obtained from the reaction of CHBr3 with potassium atoms. The overall rate constant was determined from the decay of CH with distance, O2 being introduced in excess. The result, after corrections for axial and radial diffusion, is k = (3.6 +/- 0.5) x 10(-11) cm3 molecule-1 s-1. The OH(A2 sigma +) chemiluminescence was observed, confirming the existence of the OH + CO channel. The vibrational population distribution of OH(A2 sigma +) is 32% in the v' = 1 level and 68% in the v' = 0 level (+/- 5%). The relative atomic concentrations were determined by resonance fluorescence in the vacuum ultraviolet. A ratio of 1.4 +/- 0.2 was found between the H atom density (H atoms being produced from the H + CO2 channel and from the HCO dissociation) and the O atom density (O + HCO). Ab initio calculations of the transition structures have been performed, associated with statistical estimations. The estimated branching ratios are: O + HCO, 20%; O + H + CO, 30%; H + CO2, 30%; and CO + OH, 20%.     

Quasi-classical trajectory study of the H + CO2 → HO + CO reaction on a new ab initio based potential energy surface

A detailed quasi-classical trajectory study of the H + CO(2) → HO + CO reaction is reported on an accurate potential energy surface based on ab initio data. The influence of the vibrational and rotational excitations of CO(2) was investigated up to the collision energy of 2.35 eV. It was found that the total reaction integral cross section increases monotonically with the collision energy, consistent with experimental results. The excitation of the CO(2) bending vibration enhances the reaction, while the excitation in its asymmetric stretching vibration inhibits the reaction. The calculated thermal rate constants are in excellent agreement with experiment. At the state-to-state level, the rotational state distributions of the HO product are in good agreement with experimental results, while those for the CO product are much hotter than measurements. The calculated differential cross sections are dominated by forward scattering, suggesting that the lifetime of the HOCO intermediate may not be sufficiently long to render the reaction completely statistical.