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.     

Cyanomethylidyne: a reactive carbyne radical.

The cyanomethylidyne radical (CCN) has been a long-standing subject of extensive structural and spectroscopic studies. However, its chemical reactivity has received rather little attention. Recently, we studied the reaction of CCN with the simplest alkane, CH4, which follows a mechanism of carbyne insertion-dissociation rather than that of direct H abstraction proposed by a recent experimental study. However, we are aware that alkanes like CH4 bear no electron lone pairs and thus are not ideal diagnostic molecules for distinguishing between the carbyne-insertion and H-abstraction mechanisms. Hence, we chose a series of sigma-bonded molecules HX (X=OH, NH2, and F) which bear electron lone pairs and are better diagnostics for carbyne-insertion behavior. The new results at the CCSD(T)/6-311+G(2df,p)//B3LYP/6-311G(d,p)+ZPVE, CCSD(T)/aug-cc-pVTZ//B3LYP/6-311G(d,p)+ZPVE, G2M(CC1), and MC-QCISD//B3LYP/6-31G(d)+ZPVE levels definitively confirm the carbyne-insertion behavior of the CCN radical towards HX. In addition, we make the first attempt to understand the reactivity of the CCN radical toward pi-bonded molecules, using the CCN+C2H2 model reaction. This reaction involves carbenoid addition to the C[triple chemical bond]C bond without a potential-energy barrier to form a C3 three-membered cyclic intermediate followed by H extrusion. Therefore, the reactions of CCN with both sigma- and pi-bonded molecules conclusively show that CCN is a reactive carbyne radical and may be more reactive than the well-known CN radical. Future experimental studies, especially on product characterization, are strongly desired to test our proposed carbyne mechanism. The studied reactions of CCN with CH4, NH3, H2O, and C2H2 could be of interest to combustion science and astrophysics, and they could provide efficient routes to form novel cyano-containing molecules in interstellar space.     

Mechanism of HCS + O2 reaction: hydrogen- or oxygen-transfer?

In spite of the potential importance of the HCS radical in both combustion and interstellar processes, its chemical reactivity has not been tackled previously. In the present paper, the oxidation reaction of the HCS radical is theoretically investigated for the first time at the CCSD(T)/6-311++G(3df,2p)//BH&HLYP/6-311++G(d,p)+ZPVE and Gaussian-3//B3LYP/6-31G(d) levels. It is shown that the most feasible pathway is the O2 addition to the HCS radical forming the intermediate SC(H)OO which can undergo a subsequent O-extrusion leading to SC(H)O + 3O. This features an indirect O-transfer mechanism with the overall barrier of 4.4 and 3.5 kcal mol(-1), respectively, at the two levels. However, formation of the H-transfer product CS + HO2 is kinetically much less feasible, i.e., the direct mechanism has barriers of 14.3 and 8.7 kcal mol(-1), whereas the indirect mechanism has barriers of 12.6 and 10.7 kcal mol(-1), respectively. This result is in sharp contrast to the analogous HCO + O2 reaction, where the direct (with a barrier of 2.98 kcal mol(-1)) and indirect (2.26 kcal mol(-1)) H-transfer processes are highly competitive over the indirect O-transfer process (the least endothermicity is 19.9 kcal mol(-1)). The possible explanations and implications of the present results are provided.     

三氟化氯和环氧丙烷反应的理论研究

应用密度泛函理论对三氟化氯和环氧丙烷反应产生C3H5O和C1F2自由基的机理进行了研究。在B3PW91/6-31+G(d,p)水平优化了12个不同反应通道上各驻点(反应物、中间体、过渡态和产物) 的几何构型,并计算了它们的振动频率和零点振动能。采用CCSD(T)/6-31+ G(d,p) // B3PW91/6-31+G(d,p)单点能计算方法求得各物种的能量,并作了零点能校正。计算结果表明,三氟化氯和环氧丙烷反应可经过不同的反应路径引发C3H5O自由基和C1F2自由基,其中,三氟化氯呈对称的F原子与环氧丙烷的C(1)上与CH3在同一侧的上的H原子结合的活化能最低,仅为16.81 kJ/mol。     

Radical-molecule reactions HCO/HOC + C2H2: mechanistic study

A detailed computational study is performed on the unknown radical-molecule reactions between HCO/HOC and acetylene (C2H2) at the CCSD(T)/6-311G(2d,p)//B3LYP/6-311G(d,p)+ZPVE, Gaussian-3//B3LYP/6-31G(d), and Gaussian-3//MP2(full)/6-31G(d) levels. For the HCO + C2H2 reaction, the most favorable pathway is direct C-addition forming the intermediate HC=CHCH=O followed by a 1,3-H-shift leading to H2C=CHC=O, which finally dissociates to the product C2H3 + CO. The overall reaction barrier is 13.8, 10.5, and 11.3 kcal/mol, respectively, at the three levels. The quasi-direct H-donation process to produce C2H3 + CO with barriers of 14.0, 14.1, and 14.1 kcal/mol is less competitive. Thus only at higher temperatures could the HCO + C2H2 reaction play a role. In contrast, the HOC + C2H2 reaction can barrierlessly generate C2H3 + CO via the quasi-direct H-donation mechanism proceeding via a prereactive complex with OH...C2 hydrogen bonding. This is suggestive of the potential importance of the HOC + C2H2 reaction in both combustion and interstellar processes. However, the direct C-addition channel is much less competitive. For both reactions, the possible formation of the intriguing interstellar molecules propadiene and propynal is also discussed. The present theoretical study represents the first attempt to probe the reaction mechanism between HOC and pi-systems. Future laboratory investigations on both reactions (particularly HOC + C2H2) are recommended.     

Kinetics of hydrogen abstraction reactions of butene isomers by OH radical

The rate coefficients of H-abstraction reactions of butene isomers by the OH radical were determined by both canonical variational transition-state theory and transition-state theory, with potential energy surfaces calculated at the CCSD(T)/6-311++G(d,p)//BH&HLYP/6-311G(d,p) level and CCSD(T)/6-311++G(d,p)//BH&HLYP/cc-pVTZ level and quantum mechanical tunneling effect corrected by either the small-curvature tunneling method or the Eckart method. While 1-butene contains allylic, vinylic, and alkyl hydrogens that can be abstracted to form different butene radicals, results reveal that s-allylic H-abstraction channels have low and broad energy barriers, and they are the most dominant channels which can occur via direct and indirect H-abstraction channels. For the indirect H-abstraction s-allylic channel, the reaction can proceed via forming two van der Waals prereactive complexes with energies that are 2.7-2.8 kcal mol(-1) lower than that of the entrance channel at 0 K. Assuming that neither mixing nor crossover occurs between different reaction pathways, the overall rate coefficient was calculated by summing the rate coefficients of the s-allyic, methyl, and vinyl H-abstraction paths and found to agree well with the experimentally measured OH disappearance rate. Furthermore, the rate coefficients of p-allylic H abstraction of cis-2-butene, trans-2-butene, and isobutene by the OH radical were also determined at 300-1500 K, with results analyzed and compared with available experimental data.     

Theoretical study on structures and stability of C4P isomers

The structures, energetics, spectroscopies, and stabilities of doublet C(4)P isomeric species are explored at the DFT/B3LYP, QCISD, and CCSD(T) (singlet-point) levels. A total of 12 minimum isomers and 27 interconversion transition states are located. At the CCSD(T)/6-311G(2df)//QCISD/6-311G(d)+ZPVE level, the lowest-lying isomer is a floppy CCCCP 1 (0.0 kcal/mol) mainly featuring a cumulenic structure |C=C=C=C=P*|, which differs much from the analogous C4N radical (|*C-C[triple bond]C-C[triple bond]N|). The quasi-linearity and the low bending mode of 1 are in contrast to the previous prediction. The second energetically followed isomer PC-cCCC 3 (14.9 kcal/mol) possesses a CCC ring-bonded to CP. The two low-lying isomers are separated by a high-energy ring-closure/open transition state (26.5 kcal/mol) and thus are very promising candidates for future laboratory and astrophysical detection. Furthermore, four high-energy isomers, that is, two bent isomers CCPCC 2 (68.4 kcal/mol) and CCPCC 2' (68.5 kcal/mol) and two cagelike species 10 (56.0 kcal/mol) and 11 (67.9 kcal/mol), are also stabilized by considerable barriers. The present work is the first detailed potential energy survey of CnP clusters and can provide useful information for the investigation of larger CnP radicals and for understanding the isomerism of P-doped C vaporization processes.     

Theoretical study on structures and stability of Si2CP isomers

The structures, energetics, spectroscopies, and isomerization of various doublet Si2CP species are explored theoretically. In contrast to the previously studied SiC2N and SiC2P radicals that have linear SiCCN and SiCCP ground states, the title Si2CP radical has a four-membered-ring form cSiSiPC 1 (0.0 kcal/mol) with Si-C cross-bonding as the ground-state isomer at the CCSD(T)/6-311G(2df)//B3LYP/6-311G(d)+ZPVE level, similar to the Si2CN radical. The second low-lying isomer 2 at 11.6 kcal/mol has a SiCSiP four-membered ring with C-P cross-bonding, yet it is kinetically quite unstable toward conversion to 1 with a barrier of 3.5 kcal/mol. In addition, three cyclic species with divalent carbene character, i.e., cSiSiCP 7, 7' with C-P cross-bonding and cSiCSiP 8 with Si-Si cross-bonding, are found to possess considerable kinetic stability, although they are energetically high lying at 44.4, 46.5, and 41.4 kcal/mol, respectively. Moreover, a linear isomer SiCSiP 5 at 44.3 kcal/mol also has considerable kinetic stability and predominantly features the interesting cumulenic /Si=C=Si=P/* form with a slight contribution from the silicon-phosphorus triply bonded form /Si=C*-Si[triple bond]P/. The silicon-carbon triply bonded form *Si[triple bond]C-Si[triple bond]P/ has negligible contribution. All five isomers are expected to be observable in low-temperature environments. Their bonding nature and possible formation strategies are discussed. For relevant species, the QCISD/6-311G(d) and CCSD(T)/6-311+G(2df) (single-point) calculations are performed to provide more reliable results. The calculated results are compared to those of the analogous C3N, C3P, SiC2N, and Si2CN radicals with 17 valence electrons. Implications in interstellar space and P-doped SiC vaporization processes are also discussed.     

Theoretical study on structures and stabilities of [H,Ge,C,N

Theoretical investigations are performed for the first time on the simplest hydrogenated germanium cyanide [H,Ge,C,N], whose analogs [H,C(2),N] and [H,Si,C,N] have been detected in space and laboratory, respectively. The detailed potential energy surfaces in both singlet and triplet states are constructed at the CCSD(T)/6-311+G(3df,2p)//B3LYP/6-31G(d)+ZPVE level, including 18 minimum isomers and 26 interconversion transition states. The former three low-lying and kinetically stabilized isomers are HGeCN (1)1 (0.0 kcal/mol), HGeNC (1)2 (5.1 kcal/mol), and cyclic cCHNGe(1)7 (11.1 kcal/mol). In addition, five isomers HCNGe (1)3 (33.8), HNCGe (1)5 (29.8), cNHCGe (1)8 (37.9), HGeCN (3)1 (30.1), and HNCGe (3)5 (26.5) each have considerable barriers, despite their high energies. Future laboratory characterization and astrophysical detection of the eight [H,Ge,C,N] isomers, especially the former three low-lying species (1)1, (1)2, and (1)7, are highly recommended. The accurate spectroscopic data at the QCISD/6-311G(d,p) level are provided. For some species, the CBS-QB3 calculations are also performed. Wherever possible, comparisons with the analogous [H,C(2),N] and [H,Si,C,N] are made on the structural, energetic, and bonding properties.     

NC2S+离子的结构和稳定性的理论研究

在密度泛函和从头算理论水平下计算了单重态的NC2S+离子的结构、能量、光谱以及稳定性. 在B3LYP/6-311G(d)水平下, 得到8个异构体, 它们由15个过渡态相连接. 在CCSD(T)/6-311+G(2df)//QCISD/6-311G(d)+ZPVE水平下, 得到能量最低的异构体是直线型的具有1Σ电子态的NCCS+(1)(0.0 kJ/mol), 其次是直线型的异构体CNCS+(2)(54.8 kJ/mol). 两个低能量的异构体1和2及另外一个高能量的直线型异构体CCNS+(3)(323.8 kJ/mol)都具有相当大的动力学稳定性, 这三个异构体在具备一定条件的实验室和星际条件下是可以进行观测的. 分析了这3个异构体的成键性质.     

Theoretical study on stability and properties of NC2O isomers

The structures, energetics, spectroscopies, and stabilities of the doublet NC(2)O radical are explored at density functional theory and ab initio levels. Nine minimum isomers are located connected by 22 interconversion transition states. At the CCSD(T)/6-311+G(2df)//QCISD/6-311G(d)+ZPVE level, the lowest-lying isomer is bent NCCO 1 (0.0 kcal/mol) with (2)A' state followed by bent isomer CNCO 2 (16.7). Two isomers (1 and 2) and another high-lying species CCNO 4 (99.4) with bent structure are considerably stabilized by a barrier of at least 20 kcal/mol. All of the three isomers should be experimentally or astrophysically observable. This result is consistent with their indication of neutralization-reionization mass spectrometry experiments. Also, the calculated spectroscopic properties and bond distances of known NCCO 1 are consistent with recent experimental observations and theoretical studies. The bonding natures of the isomers 1, 2, and 4 are analyzed. Their molecular properties including the heats of formation, adiabatic ionization potentials, and adiabatic electronic affinities are calculated at the higher levels G3//B3LYP, G3(MP2)//B3LYP, QCISD, and CCSD(T) (single-point). Possible formation strategies of the isomers 1, 2, and 4 in laboratory and space are also discussed in detail.     

乙烯酮自由基(HCCO·)与二氧化氮(NO2)反应势能面的理论研究

在CCSD(T)/6-311G(d,p)//MP2/6-311G(d,p)+ZPE水平上对反应HCCO+NO2进行了计算, 建立了反应势能面. 此反应由反应物通过三步反应到达产物. 首先, NO2的O原子进攻HCCO自由基中与H相邻的C原子, 形成异构体1[ONOC(H)CO]或2[H(CONOC)O]. 然后, 异构体1和2通过N-O键的断裂形成产物NO和OC(H)CO. 最后, 产物中的OC(H)CO可以通过C-C键的断裂进一步分解为HCO和CO. 由HCCO+NO2反应得到产物NO+HCO+CO.     

Isomerization pathways from the norbornadiene to the cycloheptatriene radical cation by opening a bridgehead-methylene bond: a theoretical investigation

Three skeletal rearrangement channels for the norbornadiene (N*+) to the 1,3,5-cycloheptatriene (CHT*+) radical cation conversion, initialized by opening a bridgehead-methylene bond in N*+, are investigated using the quantum chemical B3LYP, MP2 and CCSD(T) methods in conjunction with the 6-311 +G(d,p) basis set. Two of the isomerizations proceed through the norcaradiene radical cation (NCD*+), either through a concerted path (N*+ - NCD*+), or by a stepwise mechanism via a stable intermediate (N*+ - I1 - NCD*+). At the CCSD(T)/6-311 +G(d,p)//B3LYP/6-311 +G(d,p) level, the lowest activation energy, 28.9 kcal mol(-1), is found for the concerted path whereas the stepwise path is found to be 2.3 kcal mol(-1) higher. On both pathways, NCD*+ rearranges further to CHT*+ with significantly less activation energy. The third channel proceeds from N*+ through I1 and then directly to CHT*+, with an activation energy of 37.1 kcal mol(-1). The multi-step channel reported earlier by our group, which proceeds from N*+ to CHT*+ via the quadricyclane and the bicyclo[2.2.1]hepta-2-ene-5-yl-7-ylium radical cations, is 4.6 kcal mol(-1) lower than the most favorable path of the present study. If the methylene group is substituted with C(CH3)2, however, the concerted path is estimated to be 5.6 kcal mol(-1) lower than the corresponding substituted multi-step path at the B3LYP/6-311+(d,p) level. This shows that substitution of particular positions can have dramatic effects on altering reaction barriers in the studied rearrangements. We also note that identical energies are computed for CHT*+ and NCD*+ whereas, in earlier theoretical investigations, the former was reported to be 6-17 kcal mol(-1) more stable than the latter. Finally, a bent geometry is obtained for CHT*+ with MP2/6-311 +G(d,p) in contradiction with the planar conformation reported for this cation in earlier computational studies.     

Reaction mechanism between carbonyl oxide and hydroxyl radical: a theoretical study

The reaction mechanism of carbonyl oxide with hydroxyl radical was investigated by using CASSCF, B3LYP, QCISD, CASPT2, and CCSD(T) theoretical approaches with the 6-311+G(d,p), 6-311+G(2df, 2p), and aug-cc-pVTZ basis sets. This reaction involves the formation of H2CO + HO2 radical in a process that is computed to be exothermic by 57 kcal/mol. However, the reaction mechanism is very complex and begins with the formation of a pre-reactive hydrogen-bonded complex and follows by the addition of HO radical to the carbon atom of H2COO, forming the intermediate peroxy-radical H2C(OO)OH before producing formaldehyde and hydroperoxy radical. Our calculations predict that both the pre-reactive hydrogen-bonded complex and the transition state of the addition process lie energetically below the enthalpy of the separate reactants (DeltaH(298K) = -6.1 and -2.5 kcal/mol, respectively) and the formation of the H2C(OO)OH adduct is exothermic by about 74 kcal/mol. Beyond this addition process, further reaction mechanisms have also been investigated, which involve the abstraction of a hydrogen of carbonyl oxide by HO radical, but the computed activation barriers suggest that they will not contribute to the gas-phase reaction of H2COO + HO.     

Theoretical study of the addition and abstraction reactions of hydroxyl radical with uracil

The addition as well as abstraction reactions of hydroxyl radical (OH) with the nucleic acid base, uracil (U), in the gas phase has been explored at the B3LYP/6-31+G(d,p) level of density functional theory (DFT). The energy barrier of the OH addition to both the C5 and C6 positions of the uracil is less than 1 kcal/mol while the hydrogen abstractions (H-abstractions) from either the N1 or the N3 positions are ∼9.5 kcal/mol. Further the energetics of these reactions are assessed by applying the effect of aqueous medium through the polarizable continuum model (PCM). Both the gas and the solution phase data established that the thermodynamic and kinetic factors are more favorable for the OH addition to either C5 or C6 positions of the uracil than the H-abstraction reactions. Moreover, calculations at the MPW1K/6-31+G(d,p), CCSD(T)/6-31+G(d,p)//B3LYP/6-31+G(d,p) and CCSD(T)/6-31+G(d,p)//MPW1K/6-31+G(d,p) levels of theoretical methods qualitatively supported the B3LYP/6-31+G(d,p) results.     

H-Atom abstraction from CH(3)NHNH(2) by NO(2): CCSD(T)/6-311++G(3df,2p)//MPWB1K/6-31+G(d,p) and CCSD(T)/6-311+G(2df,p)//CCSD/6-31+G(d,p) calculations

Stationary points of paths for H atom abstraction from CH(3)NHNH(2) (monomethylhydrazine) by NO(2) were characterized via CCSD(T)/6-311++G(3df,2p)//MPWB1K/6-31+G(d,p) and CCSD(T)/6-311+G(2df,p)//CCSD/6-31+G(d,p) calculations. Five transition states connecting CH(3)NHNH(2)-NO(2) complexes to a manifold that includes CH(3)NHNH-HONO, CH(3)NNH(2)-HONO, CH(3)NNH(2)-HNO(2), and CH(3)NHNH-HNO(2) complexes were identified. Transition states that connect CH(3)NHNH-HONO, CH(3)NNH(2)-HONO, CH(3)NNH(2)-HNO(2), and CH(3)NHNH-HNO(2) complexes to each other via H atom exchange and/or hindered internal rotation were also identified. The high point in the minimum energy path from the CH(3)NHNH(2) + NO(2) reactant asymptote to the manifold of HONO-containing product states is a transition state 8.6 kcal/mol above the reactant asymptote. From a kinetics standpoint, this value is considerably higher than the 5.9 kcal/mol value that was estimated for it based on theoretical results for H atom abstraction from NH(3) by NO(2).     

Thermodynamic and ab initio analysis of the controversial enthalpy of formation of formaldehyde.

There are two values, -26.0 and -27.7 kcal mol(-1), that are routinely reported in literature evaluations for the standard enthalpy of formation, Delta(f) H(o)(298), of formaldehyde (CH(2)=O), where error limits are less than the difference in values. In this study, we summarize the reported literature for formaldehyde enthalpy values based on evaluated measurements and on computational studies. Using experimental reaction enthalpies for a series of reactions involving formaldehyde, in conjunction with known enthalpies of formation, its enthalpy is determined to be -26.05+/-0.42 kcal mol(-1), which we believe is the most accurate enthalpy currently available. For the same reaction series, the reaction enthalpies are evaluated using six computational methods: CBS-Q, CBS-Q//B3, CBS-APNO, G2, G3, and G3B3 yield Delta(f) H(o)(298)=-25.90+/-1.17 kcal mol(-1), which is in good agreement to our experimentally derived result. Furthermore, the computational chemistry methods G3, G3MP2B3, CCSD/6-311+G(2df,p)//B3LYP/6-31G(d), CCSD(T)/6-311+G(2df,p)//B3LYP/6-31G(d), and CBS-APNO in conjunction with isodesmic and homodesmic reactions are used to determine Delta(f) H(o)(298). Results from a series of five work reactions at the higher levels of calculation are -26.30+/-0.39 kcal mol(-1) with G3, -26.45+/-0.38 kcal mol(-1) with G3MP2B3, -26.09+/-0.37 kcal mol(-1) with CBS-APNO, -26.19+/-0.48 kcal mol(-1) with CCSD, and -26.16+/-0.58 kcal mol(-1) with CCSD(T). Results from heat of atomization calculations using seven accurate ab initio methods yields an enthalpy value of -26.82+/-0.99 kcal mol(-1). The results using isodesmic reactions are found to give enthalpies more accurate than both other computational approaches and are of similar accuracy to atomization enthalpy calculations derived from computationally intensive W1 and CBS-APNO methods. Overall, our most accurate calculations provide an enthalpy of formation in the range of -26.2 to -26.7 kcal mol(-1), which is within computational error of the suggested experimental value. The relative merits of each of the three computational methods are discussed and depend upon the accuracy of experimental enthalpies of formation required in the calculations and the importance of systematic computational errors in the work reaction. Our results also calculate Delta(f) H(o)(298) for the formyl anion (HCO(-)) as 1.28+/-0.43 kcal mol(-1).     

Theoretical study on the [Si, C, N, O] potential energy surface

The structures, energetics, spectroscopies, and stabilities of the doublet [Si, C, N, O] radical are explored at the density functional theory and ab initio levels. Sixteen isomers are located, connected by 29 interconversion transition states. At the CCSD(T)/6-311+G(2df)//QCISD/6-311G(d)+ZPVE level, the lowest lying isomer is a linear SiNCO 1 (0.0 kcal/mol) mainly featuring a cumulene | . Si = N = C = O. The second and third low-lying isomers are bent OSiCN 2 (8.8) and bent OSiNC 3 (11.1), respectively. All the three low-lying isomers 1, 2, 3, and another high-lying species 5 (75.4) with a linear SiCNO structure are shown to have considerable kinetic stability and may be experimentally observable. The predicted results of isomers 1 and 2 are consistent with the previous mass spectrometry experiments. Moreover, the fourth low-lying species SiOCN 4 (23.9) with bent structure is expected to be observable in low-temperature environments. The bonding nature of the five isomers 1, 2, 3, 4, and 5 is analyzed. The calculated results are compared with those of the analogous molecules C(2)NO and Si(2)NO. Implications in interstellar space and N,O-doped SiC vaporization processes are also discussed.     

Theoretical study on reactions of cyclic-N3 with NO, NO2 and Cl2

In this article, we report our detailed mechanistic study on the reactions of cyclic-N₃ with NO, NO₂ at the G3B3//B3LYP/6-311+G(d) and CCSD(T)/aug-cc-pVTZ//QCISD/6-311+G(d)+ZPVE levels; the reactions of cyclic-N₃ with Cl₂ was studied at the G3B3//B3LYP/6-311+G(d) and CCSD(T)/aug-cc-pVTZ//QCISD/6-31+G(d)+ZPVE levels. Both of the singlet and triplet potential-energy surfaces (PESs) of cyclic-N₃ + NO, cyclic-N₃ + NO₂ and the PES of cyclic-N₃ + Cl₂ have been depicted. The results indicate that on singlet PESs cyclic-N₃ can undergo the barrierless addition–elimination mechanism with NO and NO₂ forming the respective dominant products N₂ + ¹cyclic-NON and ¹NNO(O) + N₂. Yet the two reactions on triplet PESs are much less likely to take place under room temperature due to the high barriers. For the cyclic-N₃ + Cl₂ reaction, a Cl-abstraction mechanism was revealed that results in the product cyclic-N₃Cl + Cl with an overall barrier as high as 14.7 kcal/mol at CCSD(T)/aug-cc-pVTZ//QCISD/6-31+G(d)+ZPVE level. So the cyclic-N₃ radical could be stable against Cl₂ at low temperatures in gas phase. The present results can be useful for future experimental investigation on the title reactions.     

O+HCNO反应势能面的理论研究

采用CCSD(T)/6-311G(d,p)//B3LYP/6-311G(d,p)+ZPVE方法对反应O＋HCNO进行了研究. 通过反应势能面揭示了该反应的机理, 通过H或O迁移等多步反应路径得到3种产物, 其中, P1(HCO+NO)为主要产物, P2(HNO+CO)和P3(NCO+OH)为次要产物. 为进一步实验研究提供了参考.