Quantum chemical and statistical rate study of the reaction of O(3P) with allene: O-addition and H-abstraction channels

The lowest-lying triplet and singlet potential energy surfaces for the O(3P) + CH2=C=CH2 reaction were theoretically characterized using the complete basis set model chemistry, CBS-QB3. The primary product distributions for the multistate multiwell reactions on the individual surfaces were then determined by RRKM statistical rate theory and weak-collision master equation analysis using the exact stochastic simulation method. The results predict that the electrophilic O-addition pathways on the central and terminal carbon atom are dominant up to combustion temperatures. Major predicted end-products for the addition routes include CO + C2H4, 3CH2 + H2CCO, and CH2=C*-CHO + H*, in agreement with experimental evidence. CO + C2H4 are mainly generated from the lowest-lying singlet surface after an intersystem crossing process from the initial triplet surface. Efficient H-abstraction pathways are newly identified and occur on two different electronic state surfaces, 3A' and 3A', resulting in OH + propargyl radicals; they are predicted to play an important role at higher temperatures in hydrocarbon combustion chemistry and flames, with estimated contributions of ca. 35% at 2000 K. The overall thermal rate coefficient k(O + C3H4) at 200-1000 K was computed using multistate transition state theory: k(T) = 1.60 x 10(-17) x T (2.05) x exp(-90 K/T) cm3 molecule(-1) s(-1), in good agreement with experimental data available for the 300-600 K range.     

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.     

Potential energy surfaces, product distributions and thermal rate coefficients of the reaction of O(3P) with C2H4(X1Ag): a comprehensive theoretical study

The potential energy surface for the O((3)P) + C(2)H(4) reaction, which plays an important role in C(2)H(4)/O(2) flames and in hydrocarbon combustion in general, was theoretically reinvestigated using various quantum chemical methods, including G3, CBS-QB3, G2M(CC,MP2), and MRCI. The energy surfaces of both the lowest-lying triplet and singlet electronic states were constructed. The primary product distribution for the multiwell multichannel reaction was then determined by RRKM statistical rate theory and weak-collision master equation analysis using the exact stochastic simulation method. Intersystem crossing of the "hot" CH(2)CH(2)O triplet adduct to the singlet surface, shown to account for about half of the products, was estimated to proceed at a rate of approximately 1.5 x 10(11) s(-1). In addition, the thermal rate coefficients k(O + C(2)H(4)) in the T = 200-2000 K range were computed using multistate transition state theory and fitted by a modified Arrhenius expression as k(T) = 1.69 x 10(-16) x T(1.66) x exp(-331 K/T) . Our computed rates and product distributions agree well with the available experimental results. Product yields are found to show a monotonic dependence on temperature. The major products (with predicted yields at T = 300 K/2000 K) are: CH(3) + CHO (48/37%), H + CH(2)CHO (40/19%), and CH(2)(X(3)B(1)) + H(2)CO (5/29%), whereas H + CH(3)CO, H(2) + H(2)CCO, and CH(4) + CO are all minor (< or =5%).     

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相似文献   

Radical reaction C3H+NO: a mechanistic study

Although a number of hydrocarbon radicals including the heavier C(3)-radicals C(3)H(3) and C(3)H(5) have been experimentally shown to deplete NO effectively, no theoretical or experimental attempts have been made on the reactivity of the simplest C(3)-radical towards NO. In this article, we report our detailed mechanistic study on the C(3)H+NO reaction at the Gussian-3//B3LYP/6-31G(d) level by constructing the singlet and triplet electronic state [H,C(3),N,O] potential energy surfaces (PESs). The l-C(3)H+NO reaction is shown to barrierlessly form the entrance isomer HCCCNO followed by the direct O-elimination leading to HCCCN+(3)O on triplet PES, or by successive O-transfer, N-insertion, and CN bond-rupture to generate the product (1)HCCN+CO on singlet PES. The possible singlet-triplet intersystem crossings are also discussed. Thus, the novel reaction l-C(3)H+NO can proceed effectively even at low temperatures and is expected to play an important role in both combustion and interstellar processes. For the c-C(3)H+NO reaction, the initially formed H-cCCC-NO can most favorably isomerize to HCCCNO, and further evolution follows that of the l-C(3)H+NO reaction. Quantitatively, the c-C(3)H+NO reaction can take place barrierlessly on singlet PES, yet it faces a small barrier 2.7 kcal/mol on triplet PES. The results will enrich our understanding of the chemistry of the simplest C(3)-radical in both combustion and interstellar processes, which to date have received little attention despite their importance and available abundant studies on its structural and spectroscopic properties.     

Experimental and theoretical studies of the C2F4 + O reaction: nonadiabatic reaction mechanism

In this work, the C(2)F(4)(X(1)A(g)) + O((3)P) reaction was investigated experimentally using molecular beam-threshold ionization mass spectrometry (MB-TIMS). The major primary products were observed to be CF(2)O (+ CF(2)) and CF(3) (+ CFO), with measured approximate yields of % versus %, respectively, neglecting minor products. Furthermore, the lowest-lying triplet and singlet potential energy surfaces for this reaction were constructed theoretically using B3LYP, G2M(UCC, MP2), CBS-QB3, and G3 methods in combination with various basis sets such as 6-31G(d), 6-311+G(3df), and cc-pVDZ. The primary product distribution for the multiwell multichannel reaction was then determined by RRKM statistical rate theory and weak-collision master equation analysis. It was found that the observed production of CF(3) (+ CFO) can only occur on the singlet surface, in parallel with formation of ca. 5 times more CF(2)O(X) + CF(2)(X(1)A(1)). This requires fast intersystem crossing (ISC) from the triplet to the singlet surface at a rate of ca. 4 x 10(12) s(-1). The theoretical calculations combined with the experimental results thus indicate that the yield of triplet CF(2)(?(3)B(1)) + CF(2)O formed on the triplet surface prior to ISC is < or =35%, whereas singlet CF(2)(X(1)A(1)) + CF(2)O is produced with yield > or =60%, after ISC. In addition, the thermal rate coefficients k(O + C(2)F(4)) in the T = 150-1500 K range were computed using multistate transition state theory and can be expressed as k(T) = 1.67 x 10(-16) x T(1.48) cm(3) molecule(-1) s(-1); they are in agreement with the available experimental results in the T = 298-500 K range.     

New insight into the gas-phase bimolecular self-reaction of the HOO radical

The singlet and triplet potential energy surfaces (PESs) for the gas-phase bimolecular self-reaction of HOO*, a key reaction in atmospheric environments, have been investigated by means of quantum-mechanical electronic structure methods (CASSCF and CASPT2). All the reaction pathways on both PESs consist of a first step involving the barrierless formation of a prereactive doubly hydrogen-bonded complex, which is a diradical species lying about 8 kcal/mol below the energy of the reactants at 0 K. The lowest energy reaction pathway on both PESs is the degenerate double hydrogen exchange between the HOO* moieties of the prereactive complex via a double proton transfer mechanism involving an energy barrier of only 1.1 kcal/mol for the singlet and 3.3 kcal/mol for the triplet at 0 K. The single H-atom transfer between the two HOO* moieties of the prereactive complex (yielding HOOH + O2) through a pathway keeping a planar arrangement of the six atoms involves a conical intersection between either two singlet or two triplet states of A' and A" symmetries. Thus, the lowest energy reaction pathway occurs via a nonplanar cisoid transition structure with an energy barrier of 5.8 kcal/mol for the triplet and 17.5 kcal/mol for the singlet at 0 K. The simple addition between the terminal oxygen atoms of the two HOO* moieties of the prereactive complex, leading to the straight chain H2O4 intermediate on the singlet PES, involves an energy barrier of 7.3 kcal/mol at 0 K. Because the decomposition of such an intermediate into HOOH + O2 entails an energy barrier of 45.2 kcal/mol at 0 K, it is concluded that the single H-atom transfer on the triplet PES is the dominant pathway leading to HOOH + O2. Finally, the strong negative temperature dependence of the rate constant observed for this reaction is attributed to the reversible formation of the prereactive complex in the entrance channel rather than to a short-lived tetraoxide intermediate.     

Quantum chemical study of the mechanism of reaction between NH (X 3sigma-) and H2, H2O, and CO2 under combustion conditions

Reactions of ground-state NH (3sigma-) radicals with H2, H2O, and CO2 have been investigated quantum chemically, whereby the stationary points of the appropriate reaction potential energy surfaces, that is, reactants, products, intermediates, and transition states, have been identified at the G3//B3LYP level of theory. Reaction between NH and H2 takes place via a simple abstraction transition state, and the rate coefficient for this reaction as derived from the quantum chemical calculations, k(NH + H2) = (1.1 x 10(14)) exp(-20.9 kcal mol(-1)/RT) cm3 mol(-1) s(-1) between 1000 and 2000 K, is found to be in good agreement with experiment. For reaction between triplet NH and H2O, no stable intermediates were located on the triplet reaction surface although several stable species were found on the singlet surface. No intersystem crossing seam between triplet NH + H2O and singlet HNO + H2 (the products of lowest energy) was found; hence there is no evidence to support the existence of a low-energy pathway to these products. A rate coefficient of k(NH + H2O) = (6.1 x 10(13)) exp(-32.8 kcal mol(-1)/RT) cm3 mol(-1) s(-1) between 1000 and 2000 K for the reaction NH (3sigma-) + H2O --> NH2 (2B) + OH (2pi) was derived from the quantum chemical results. The reverse rate coefficient, calculated via the equilibrium constant, is in agreement with values used in modeling the thermal de-NO(x) process. For the reaction between triplet NH and CO2, several stable intermediates on both triplet and singlet reaction surfaces were located. Although a pathway from triplet NH + CO2 to singlet HNO + CO involving intersystem crossing in an HN-CO2 adduct was discovered, no pathway of sufficiently low activation energy was discovered to compare with that found in an earlier experiment [Rohrig, M.; Wagner, H. G. Proc. Combust. Inst. 1994, 25, 993.].     

A crossed molecular beams and ab initio study on the formation of C6H3 radicals. an interface between resonantly stabilized and aromatic radicals

The crossed molecular beams reaction of dicarbon molecules, C(2)(X(1)Σ(g)(+)/a(3)Π(u)) with vinylacetylene was studied under single collision conditions at a collision energy of 31.0 kJ mol(-1) and combined with electronic structure calculations on the singlet and triplet C(6)H(4) potential energy surfaces. The investigations indicate that both reactions on the triplet and singlet surfaces are dictated by a barrierless addition of the dicarbon unit to the vinylacetylene molecule and hence indirect scattering dynamics via long-lived C(6)H(4) complexes. On the singlet surface, ethynylbutatriene and vinyldiacetylene were found to decompose via atomic hydrogen loss involving loose exit transition states to form exclusively the resonantly stabilized 1-hexene-3,4-diynyl-2 radical (C(6)H(3); H(2)CCCCCCH; C(2v)). On the triplet surface, ethynylbutatriene emitted a hydrogen atom through a tight exit transition state located about 20 kJ mol(-1) above the separated stabilized 1-hexene-3,4-diynyl-2 radical plus atomic hydrogen product; to a minor amount (<5%) theory predicts that the aromatic 1,2,3-tridehydrobenzene molecule is formed. Compared to previous crossed beams and theoretical investigations on the formation of aromatic C(6)H(x) (x = 6, 5, 4) molecules benzene, phenyl, and o-benzyne, the decreasing energy difference from benzene via phenyl and o-benzyne between the aromatic and acyclic reaction products, i.e., 253, 218, and 58 kJ mol(-1), is narrowed down to only ～7 kJ mol(-1) for the C(6)H(3) system (aromatic 1,2,3-tridehydrobenzene versus the resonantly stabilized free radical 1-hexene-3,4-diynyl-2). Therefore, the C(6)H(3) system can be seen as a "transition" stage among the C(6)H(x) (x = 6-1) systems, in which the energy gap between the aromatic isomer (x = 6, 5, 4) is reduced compared to the acyclic isomer as the carbon-to-hydrogen ratio increases and the acyclic isomer becomes more stable (x = 1, 2).     

Theoretical study of low-lying triplet states of aniline

Multireference complete active space self-consistent-field CASSCF(10,12)/ANO and second-order perturbation theory MS-CASPT2 calculations were performed to determine the vertical low-lying singlet and triplet states of aniline. The sequence of the seven lower lying triplet states is T1(1(3)A'), T2(1(3)A' '), T3(2(3)A'), T4(3(3)A'), T5(2(3)A' '), T6(4(3)A'), and T7(3(3)A' '). The 3(3)A', 4(3)A', and 3(3)A' ' states are assigned as 3s, 3py, and 3pz Rydberg states, respectively, while other states correspond to pi <-- pi excitations. Both the T1 and T2 states are found to be below at the lowest-lying singlet S1 (1(1)A' ') state. Geometry, vibrational modes, and electron distribution of the lowest lying T1 state were determined using UB3LYP calculations. The vertical and adiabatic singlet-triplet energy gaps DeltaE(S0-T1) amount to 3.7 and 3.5 +/- 0.2 eV, respectively. In clear contrast with the S0 state, the triplet aniline is no longer aromatic, and its protonation occurs preferentially at the ring meta-carbon site, with a proton affinity PA = 243 +/- 3 kcal/mol.     

Gas phase hydration and deprotonation of the cyclic C3H3+ cation. Solvation by acetonitrile, and comparison with the benzene radical cation

The binding energies of the first 5 H2O molecules to c-C3H3+ were determined by equilibrium measurements. The measured binding energies of the hydrated clusters of 9-12 kcal/mol are typical of carbon-based CH+...X hydrogen bonds. The ion solvation with the more polar CH3CN molecules results in stronger bonds consistent with the increased ion-dipole interaction. Ab initio calculations show that the lowest energy isomer of the c-C3H3+(H2O)4 cluster consists of a cyclic water tetramer interacting with the c-C3H3+ ion, which suggests the presence of orientational restraint of the water molecules consistent with the observed large entropy loss. The c-C3H3+ ion is deprotonated by 3 or more H2O molecules, driven energetically by the association of the solvent molecules to form strongly hydrogen bonded (H2O)nH+ clusters. The kinetics of the associative proton transfer (APT) reaction C3H3+ + nH2O --> (H2O)nH+ + C3H2* exhibits an unusually steep negative temperature coefficient of k = cT(-63+/-4) (or activation energy of -37 +/- 1 kcal mol(-1)). The behavior of the C3H3+/water system is exactly analogous to the benzene+*/water system, suggesting that the mechanism, kinetics and large negative temperature coefficients may be general to multibody APT reactions. These reactions can become fast at low temperatures, allowing ionized polycyclic aromatics to initiate ice formation in cold astrochemical environments.     

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.     

Formation and decomposition of chemically activated cyclopentoxy radicals from the c-C5H9 + O reaction

The formation and the decomposition of chemically activated cyclopentoxy radicals from the c-C5H9 + O reaction have been studied in the gas phase at room temperature. Two different experimental arrangements have been used. Arrangement A consisted of a laser-flash photolysis set up combined with quantitative Fourier transform infrared spectroscopy and allowed the determination of the stable products at 4 mbar. The c-C5H9 radicals were produced via the reaction c-C5H10 + Cl with chlorine atoms from the photolysis of CFCl3; the O atoms were generated by photolysis of SO2. Arrangement B, a conventional discharge flow-reactor with molecular beam sampling, was used to determine the rate coefficient. Here, the hydrocarbon radicals (c-C5H9, C2H5, CH2OCH3) were produced via the reaction of atomic fluorine with c-C5H10, C2H6, and CH3OCH3, respectively, and detected by mass spectrometry after laser photoionization. For the c-C5H9 + O reaction, the relative contributions of intermediate formation (c-C5H9O) and direct abstraction (c-C5H8 + OH) were found to be 68 +/- 5 and 32 +/- 4%, respectively. The decomposition products of the chemically activated intermediate could be identified, and the following relative branching fractions were obtained: c-C5H8O + H (31 +/- 2%), CH2CH(CH2)2CHO + H (40 +/- 5%), 2 C2H4 + H + CO (17 +/- 5%), and C3H4O + C2H4 + H (12 +/- 5%). Additionally, the product formation of the c-C5H8 + O reaction was studied, and the following relative yields were obtained (mol %): C2H4, 24%; C3H4O, 18%; c-C5H8O, 30%; c-C5H8O, 23%; 4-pentenal, 5%. The rate coefficient of the c-C5H9 + O reaction was determined relative to the reactions C2H5 + O and CH3OCH2 + O leading to k = (1.73 +/- 0.05) x 10(14) cm3 mol(-1) s(-1). The experimental branching fractions are analyzed in terms of statistical rate theory with molecular and transition-state data from quantum chemical calculations, and high-pressure limiting Arrhenius parameters for the unimolecular decomposition reactions of C5H9O species are derived.     

Theoretical study of the reaction of ethane with oxygen molecules in the ground triplet and singlet delta States

Quantum chemical calculations are carried out to study the reaction of ethane with molecular oxygen in the ground triplet and singlet delta states. Transition states, intermediates, and possible products of the reaction on the triplet and singlet potential energy surfaces are identified on the basis of the coupled-cluster method. The basis set dependence of coupled-cluster energy values is estimated by the second-order perturbation theory. The values of energy barriers are also refined by using the compound CBS-Q and G3 techniques. It was found that the C(2)H(6) + O(2)(X(3)Σ(g)(-)) reaction leads to the formation of C(2)H(5) and HO(2) products, whereas the C(2)H(6) + O(2)(a(1)Δ(g)) process produces C(2)H(4) and H(2)O(2) molecules. The appropriate rate constants of these reaction paths are estimated on the basis of variational and nonvariational transition-state theories assuming tunneling and possible nonadiabatic transitions in the temperature range 500-4000 K. The calculations showed that the rate constant of the C(2)H(6) + O(2)(a(1)Δ(g)) reaction path is much greater than that of the C(2)H(6) + O(2)(X(3)Σ(g)(-)) one. At the same time, the singlet and triplet potential surface intersection is detected that leads to the appearance of the nonadiabatic quenching channel O(2)(a(1)Δ(g)) + C(2)H(6) → O(2)(X (3)Σ(g)(-)) + C(2)H(6). The rate constant of this process is estimated with the use of the Landau-Zener model. It is demonstrated that, in the case of the existence of thermal equilibrium in the distribution of molecules over the electronic states, at low temperatures (T < 1200 K) the main products of the reaction of C(2)H(6) with O(2) are C(2)H(4) and H(2)O(2), rather than C(2)H(5) and HO(2). At higher temperature (T > 1200 K) the situation is inverted.     

Reaction mechanism of the CCN radical with nitric oxide

To investigate the possibility of the carbyne radical CCN in removal of nitric oxide, a detailed computational study is performed at the Gaussian-3//B3LYP/6-31G(d) level on the CCN + NO reaction by constructing the singlet and triplet electronic state [C(2)N(2)O] potential energy surfaces (PESs). The barrierless formation of the chain-like isomers NCCNO (singlet at -106.5, triplet cis at -48.2 and triplet trans at -47.6 kcal/mol) is the most favorable entrance attack on both singlet and triplet PESs. Subsequently, the singlet NCCNO takes an O-transfer to form the branched intermediate singlet NCC(O)N (-85.6), which can lead to the fragments CN + NCO (-51.2) via the intermediate singlet NCOCN (-120.3). The simpler evolution of the triplet NCCNO is the direct N-O rupture to form the weakly bound complex triplet NCCN...O (-56.2) before the final fragmentation to NCCN + (3)O (-53.5). However, the lower lying products (3)NCN + CO (-105.6) and (3)CNN + CO (-74.6) are kinetically much less competitive. All the involved transition states for generation of CN + NCO and NCCN + (3)O lie much lower than the reactants. Thus, the novel reaction CCN + NO can proceed effectively even at low temperatures and is expected to play a role in both combustion and interstellar processes. Significant differences are found on the singlet PES between the CCN + NO and CH + NO reaction mechanisms.     

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.     

A theoretical study of the mechanism and kinetics of F+N3 reactions.

Both the singlet(1A') and triplet(3A') potential energy surfaces (PESs) of F+N(3) reactions are investigated using the complete-active-space self-consistent field (CASSCF) and the multireference configuration interaction (MRCI) methods with a proper active space. The minimum energy crossing point (MECP) at the intersection seam between the 1A' and 3A' PESs is located and used to clarify the reaction mechanisms. Two triplet transition states are found, with one in the cis form and the other one in the trans form. Further kinetic calculations are performed with the canonical unified statistical (CUS) theory on the singlet PES and the improved canonical variational transition-state (ICVT) method on the triplet PES. The rate constants are also reported. At 298 K, the calculated rate constant is in reasonably good agreement with experimental values, and spin-orbit coupling effects lower it by 28 %. The spectroscopic constants derived from the fitted potential-energy curves for the singlet and triplet states of NF are in very good agreement with experimental values. Our calculations indicate that the adiabatic reaction on the singlet PES leading to NF(a(1)Delta)+N(2) is the major channel, whereas the nonadiabatic reaction through the MECP, which leads to NF(X(3)Sigma(-))+N(2), is a minor channel.     

Ab initio kinetics of the reaction of HCO with NO: abstraction versus association/elimination mechanism

The kinetics and mechanism for the reaction of HCO with NO occurring by both singlet and triplet electronic state potential-energy surfaces (PESs) have been studied at the modified Gaussian-2 level of theory based on the geometric parameters optimized by the Becke-3 Lee-Yang-Parr/6-311G(d,p) method. There are two major reaction channels on both singlet and triplet PESs studied: one is direct H abstraction producing CO+HNO and the other is association forming a stable HC(O)NO (nitrosoformaldehyde) molecule. The dominant reaction is predicted to be the direct H abstraction occurring primarily by the lowest-energy path via a loose hydrogen-bonding singlet molecular complex, ON...HCO, with a 2.9-kcal/mol binding energy and a small decomposition barrier (1.9 kcal/mol). The commonly assumed HC(O)NO intermediate, predicted to lie below the reactants by 27.7 kcal/mol, has a high HNO-elimination barrier (34.5 kcal/mol). Bimolecular rate constants for the formation of the singlet products and their branching ratios have been calculated in the temperature range of 200-3000 K. The rate constant for the disproportionation process producing HNO+CO, found to be affected strongly by multiple reflections above the well of the complex at low temperature, is predicted to be k(HNO)=3.08 x 10(-12) T(0.10) exp(242T) for 200-500 K, and 1.72 x 10(-16) T(1.47) exp(888T) for 500-3000 K in units of cm(3) molecule(-1) s(-1). The high- and low-pressure rate constants for the association process forming HC(O)NO can be represented by k(infinity)=4.42 x 10(-11) T(0.25) exp(-28T) cm(3) molecule(-1) s(-1) (200-3000 K) and k(0)=7.30x10(-16) T(-5.75) exp(-719T) (200-1000 K) and 1.82 x 10(2) T(-11.92) exp(1846T) (1000-3000 K) cm(6) molecule(-2) s(-1) for N(2)-buffer gas. The absolute values of total rate constant, predicted to be weakly dependent negatively on temperature but positively on pressure, are in close agreement with most experimental data within their reported errors.     

A kinetic study of Ca-containing ions reacting with O, O2, CO2 and H2O: implications for calcium ion chemistry in the upper atmosphere

A series of gas-phase reactions involving molecular Ca-containing ions was studied by the pulsed laser ablation of a calcite target to produce Ca(+) in a fast flow of He, followed by the addition of reagents downstream and detection of ions by quadrupole mass spectrometry. Most of the reactions that were studied are important for describing the chemistry of meteor-ablated calcium in the earth's upper atmosphere. The following rate coefficients were measured: k(CaO(+) + O --> Ca(+) + O(2)) = (4.2 +/- 2.8) x 10(-11) at 197 K and (6.3 +/- 3.0) x 10(-11) at 294 K; k(CaO(+) + CO --> Ca(+) + CO(2), 294 K) = (2.8 +/- 1.5) x 10(-10); k(Ca(+).CO(2) + O(2) --> CaO(2)(+) + CO(2), 294 K) = (1.2 +/- 0.5) x10(-10); k(Ca(+).CO(2) + H(2)O --> Ca(+).H(2)O + CO(2)) = (13.0 +/- 4.0) x 10(-10); and k(Ca(+).H(2)O + O(2) --> CaO(2)(+) + H(2)O, 294 K) = (4.0 +/- 2.5) x 10(-10) cm(3) molecule(-1) s(-1). The quoted uncertainties are a combination of the 1sigma standard errors in the kinetic data and the systematic errors in the models used to extract the rate coefficients. Rate coefficients were also obtained for the following recombination (also termed association) reactions in He bath gas: k(Ca(+).CO(2) + CO(2) --> Ca(+).(CO(2))(2), 294 K) = (2.6 +/- 1.0) x 10(-29); k(Ca(+).H(2)O + H(2)O --> Ca(+).(H(2)O)(2)) = (1.6 +/- 1.1) x 10(-27); and k(CaO(2)(+) + O(2) --> CaO(2)(+).O(2)) < 1 x 10(-31) cm(6) molecule(-2) s(-1). These recombination rate coefficients, as well as those for the ligand-switching reactions listed above, were then interpreted using a combination of high level quantum chemistry calculations and RRKM theory using an inverse Laplace transform solution of the master equation. The surprisingly slow reaction between CaO(+) and O was explained using quantum chemistry calculations on the lowest (2)A', (2)A' and (4)A' potential energy surfaces. These calculations indicate that reaction mostly occurs on the (2)A' surface, leading to production of Ca(+)((2)S) + O(2)((1)Delta(g)). The importance of this reaction for controlling the lifetime of Ca(+) in the upper mesosphere and lower thermosphere is then discussed.     

Computational study of carbon atom (3P and 1D) reaction with CH2O. Theoretical evaluation of 1B1 methylene production by C (1D)

Singlet and triplet free energy surfaces for the reactions of C atom ((3)P and (1)D) with CH(2)O are studied computationally to evaluate the excited singlet ((1)B(1)) methylene formation from deoxygenation of CH(2)O by C ((1)D) atom as suggested by Shevlin et al. Carbon atoms can react by addition to the oxygen lone pair or to the C=O double bond on both the triplet and singlet surfaces. Triplet C ((3)P) atoms will deoxygenate to give CO plus CH(2) ((3)B(1)) as the major products, while singlet C ((1)D) reactions will form ketene and CO plus CH(2) ((1)A(1)). No definitive evidence of the formation of excited singlet ((1)B(1)) methylene was found on the singlet free energy surface. A conical intersection between the (1)A' and (1)A' ' surfaces located near an exit channel may play a role in product formation. The suggested (1)B(1) state of methylene may form via the (1)A' ' surface only if dynamic effects are important. In an effort to interpret experimental observation of products trapped by (Z)-2-butene, formation of cis- and trans-1,2-dimethylcyclopropane is studied computationally. The results suggests that "hot" ketene may react with (Z)-2-butene nonstereospecifically.