Ab initio chemical kinetics for reactions of ClO with Cl2O2 isomers

The mechanisms for the reactions of ClO with ClOClO, ClOOCl, and ClClO(2) have been investigated at the CCSD(T)/6-311+G(3df)//PW91PW91∕6-311+G(3df) level of theory. The rate constants for their low energy channels have been calculated by statistical theory. The results show that the main products for the reaction of ClO with ClOClO are ClOCl + ClOO, which can be produced readily by ClO abstracting the terminal O atom from ClOClO. This process occurs without an intrinsic barrier, with the predicted rate constant: k (ClO + ClOClO) = 7.26 × 10(-10) T(-0.15) × exp (-40/T) cm(3)molecule(-1)s(-1) for 200-1500 K. For the reactions of ClO + ClOOCl and ClClO(2), the lowest abstraction barriers are 7.2 and 7.3 kcal/mol, respectively, suggesting that these two reactions are kinetically unimportant in the Earth's stratosphere as their rate constants are less than 10(-14) cm(3)molecule(-1)s(-1) below 700 K. At T = 200-1500 K, the computed rate constants can be represented by k (ClO+ ClOOCl) = 1.11 × 10 (-14) T (0.87) exp (-3576/T) and k (ClO+ ClClO(2)) = 4.61 × 10(-14) T(0.53) exp (-3588/T) cm(3)molecule(-1)s(-1). For these systems, no experimental or theoretical kinetic data are available for comparison.     

Ab initio studies of ClO(x) reactions: prediction of the rate constants of ClO+NO2 for the forward and reverse processes.

The potential-energy surface for the reaction of ClO with NO2 has been constructed at the CCSD(T)/6-311+G(3df)//B3LYP/6-311+G(3df) level of theory. Six ClNO3 isomers are located; these are ClONO2, pc-ClOONO, pt-ClOONO, OClNO2, pt-OClONO, pc-OClONO, with predicted energies relative to the reactants of -25.6, -0.5, 1.0, 1.9, 12.2 and 13.6 kcal mol-1, and heats of formation at 0 K of 7.8, 32.9, 34.4, 35.5, 45.6 and 47.0 kcal mol-1, respectively. Isomerizations among them are also discussed. The rate constants for the low-energy pathways have been computed by statistical theory calculations. For the association reaction producing exclusively ClONO2, the predicted low- and high-pressure-limit rate constants in N2 for the temperature range of 200-600 K can be represented by: (N2)=3.19 x 10-17 T-5.54 exp(-384 K/T) cm6 molecule-2 s-1 and =3.33 x 10-7 T-1.48 exp(-18 K/T) cm3 molecule-1 s-1. The predicted low- and high-pressure-limit rate constants for the decomposition of ClONO2 in N2 at 200-600 K can be expressed, respectively, by =6.08 x 1013 T-6.54 exp(-13813 K/T) cm3 molecule-1 s-1 and =4.59 x 1023 T-2.43 exp(-13437 K/T) s-1. The predicted values compare satisfactorily with available experimental data. The reverse Cl+NO3 reaction was found to be independent of the pressure, giving exclusively ClO+NO2; the predicted rate constant can be expressed as k(Cl+NO3)=1.19 x 10-9 T-0.60 exp(58 K/T) cm3 molecule-1 s-1..     

Theoretical study of the mechanism of NO2 production from NO + ClO

The reaction of NO with ClO has been studied theoretically using density-functional and wave function methods (B3LYP and CCSD(T)). Although a barrier for cis and trans additions could be located at the RCCSD(T) and UCCSD(T) levels, no barrier exists at the B3LYP/6-311+G(d) level. Variational transition state theory on a CASPT2(12,12)/ANO-L//B3LYP/6-311+G(d) surface was used to calculate the rate constants for addition. The rate constant for cis addition was faster than that for trans addition (cis:trans 1:0.76 at 298 K). The rate constant data summed for cis and trans addition in the range 200-1000 K were fit to a temperature-dependent rate in the form kdi) = 3.30 x 10(-13)T(0.558) exp(305/T) cm3.molecule(-1).s(-1), which is in good agreement with experiment. When the data are fit to an Arrhenius plot in the range 200-400 K, an activation barrier of -0.35 kcal/mol is obtained. The formation of ClNO2 from ONOCl has a much higher activation enthalpy from the trans isomer compared to the cis isomer. In fact, the preferred decomposition pathway from trans-ONOCl to NO2 + Cl is predicted to go through the cis-ONOCl intermediate. The trans --> cis isomerization rate constant is kiso = 1.92 x 10(13) exp(-4730/T) s(-1) using transition state theory.     

Ab initio chemical kinetics for the ClOO + NO reaction: effects of temperature and pressure on product branching formation

The kinetics and mechanism for the reaction of ClOO with NO have been investigated by ab initio molecular orbital theory calculations based on the CCSD(T)/6-311+G(3df)//PW91PW91∕6-311+G(3df) method, employed to evaluate the energetics for the construction of potential energy surfaces and prediction of reaction rate constants. The results show that the reaction can produce two key low energy products ClNO + (3)O(2) via the direct triplet abstraction path and ClO + NO(2) via the association and decomposition mechanism through long-lived singlet pc-ClOONO and ClONO(2) intermediates. The yield of ClNO + O(2) ((1)△) from any of the singlet intermediates was found to be negligible because of their high barriers and tight transition states. As both key reactions initially occur barrierlessly, their rate constants were evaluated with a canonical variational approach in our transition state theory and Rice-Ramspergen-Kassel-Marcus/master equation calculations. The rate constants for ClNO + (3)O(2) and ClO + NO(2) production from ClOO + NO can be given by 2.66 × 10(-16) T(1.91) exp(341/T) (200-700 K) and 1.48 × 10(-24) T(3.99) exp(1711/T) (200-600 K), respectively, independent of pressure below atmospheric pressure. The predicted total rate constant and the yields of ClNO and NO(2) in the temperature range of 200-700 K at 10-760 Torr pressure are in close agreement with available experimental results.     

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

Protonated acetylene revisited

The potential energy surface for protonated acetylene has been re-examined with large basis sets and highly correlated methods. The energy difference of 3.6–3.8 kcal/mol between the classical structure and non-classical (bridged) structure computed with CCSD (T)/cc-pVQZ, CCSD(T)/6-311+G(3df,2pd), BD(T)/cc- pVQZ, BD(T)/6-311+G(3df,2pd) and CBS-APNO methods is in very good agreement with the best previous calculations, 3.7–4.0 kcal/mol. In contrast, BLYP, B3LYP, PW91, PBE and TPSS density functional methods do rather poorly, yielding −0.52. 0.29, 1.81, 2.16 and 0.62 kcal/mol, respectively, with the 6-311+G(3df,2pd) basis. MP2 calculations predict the classical structure to be a transition state; however, frequency calculations at the CCSD/6-311+G(3df,2pd) level of theory show that the classical structure is a local minimum. CCSD(T), BD(T) and CBS-APNO energy calculations along the MP2/6-311+G(3df,2pd) reaction path indicate that the classical structure is a shallow local minimum separated from the non-classical structure by a very small barrier of 0.11–0.13 kcal/mol. Because the barrier for proton exchange between the non-classical isomers via the classical structure is broad and nearly flat at the top, the tunneling splitting should be reduced, possibly accounting for the 15% difference between the calculated and experimental barrier heights. Contribution to the Fernando Bernardi Memorial Issue.     

Gas phase SN2 reactions of halide ions with trifluoromethyl halides: front- and back-side attack vs. complex formation

Density functional theory computations and pulsed-ionization high-pressure mass spectrometry experiments have been used to explore the potential energy surfaces for gas-phase S(N)2 reactions between halide ions and trifluoromethyl halides, X(-) + CF(3)Y --> Y(-) + CF(3)X. Structures of neutrals, ion-molecule complexes, and transition states show the possibility of two mechanisms: back- and front-side attack. From pulsed-ionization high-pressure mass spectrometry, enthalpy and entropy changes for the equilibrium clustering reactions for the formation of Cl(-)(BrCF(3)) (-16.5 +/- 0.2 kcal mol(-1) and -24.5 +/- 1 cal mol(-1) K(-1)), Cl(-)(ICF(3)) (-23.6 +/- 0.2 kcal mol(-1)), and Br(-)(BrCF(3)) (-13.9 +/- 0.2 kcal mol(-1) and -22.2 +/- 1 cal mol(-1) K(-1)) have been determined. These are in good to excellent agreement with computations at the B3LYP/6-311+G(3df)//B3LYP/6-311+G(d) level of theory. It is shown that complex formation takes place by a front-side attack complex, while the lowest energy S(N)2 reaction proceeds through a back-side attack transition state. This latter mechanism involves a potential energy profile which closely resembles a condensed phase S(N)2 reaction energy profile. It is also shown that the Cl(-) + CF(3)Br --> Br(-) + CF(3)Cl S(N)2 reaction can be interpreted using Marcus theory, in which case the reaction is described as being initiated by electron transfer. A potential energy surface at the B3LYP/6-311+G(d) level of theory confirms that the F(-) + CF(3)Br --> Br(-) + CF(4) S(N)2 reaction proceeds through a Walden inversion transition state.     

Ab initio/density functional theory and multichannel RRKM study for the ClO + CH2O reaction

The potential energy surface for the CH(2)O + ClO reaction was calculated at the QCISD(T)/6-311G(2d,2p)//B3LYP/6-311G(d,p) level of theory. The rate constants for the lower barrier reaction channels producing HOCl + HCO, H atom, OCH(2)OCl, cis-HC(O)OCl and trans-HC(O)OCl have been calculated by TST and multichannel RRKM theory. Over the temperature range of 200-2000 K, the overall rate constants were k(200-2000K) = 1.19 x 10(-13)T(0.79) exp(-3000.00/T). At 250 K, the calculated overall rate constant was 5.80 x 10(-17) cm(3) molecule(-1) s(-1), which was in good agreement with the experimental upper limit data. The calculated results demonstrated that the formation of HOCl + HCO was the dominant reaction channel and was exothermic by 9.7 kcal/mol with a barrier of 5.0 kcal/mol. When it retrograded to the reactants CH(2)O + ClO, an energy barrier of 14.7 kcal/mol is required. Furthermore, when HOCl decomposed into H + ClO, the energy required was 93.3 kcal/mol. These results suggest that the decomposition in both the forward and backward directions for HOCl would be difficult in the ground electronic state.     

Theoretical study on the hydrogen abstraction reaction of CH3CH2F (HFC-161) with Cl atom

A direct dynamics method is employed to study the hydrogen abstraction reaction of CH₃CH₂F+Cl. Three distinct transition states are located, one for -H abstraction and two for β-H abstraction. The potential energy surface (PES) information is obtained at the QCISD(T)/6-311+G(3df,2p)//MP2/6-311G(d,p), CCSD(T)/6-311+G(3df,2p)//MP2/6-311G(d,p) and G2//MP2/6-311G(d,p) level. Based on the QCISD(T)/6-311+G(3df,2p)//MP2/6-311G(d,p) results, the rate constants of the three reaction channels are evaluated by using the canonical variational transition state theory (CVT) with small-curvature tunneling (SCT) contributions over the temperature range of 220–2800 K. The calculated results indicate that -H abstraction dominates the total reaction almost over the whole temperature range.     

A 193 nm laser photofragmentation time-of-flight mass spectrometric study of chloroiodomethane

The photodissociation dynamics of chloroiodomethane (CH2ICl) at 193 nm has been investigated by employing the photofragment time-of-flight (TOF) mass spectrometric method. Using tunable vacuum ultraviolet undulator synchrotron radiation for photoionization sampling of nascent photofragments, we have identified four primary dissociation product channels: CH2Cl + I(2P(1/2))/I(2P(3/2)), CH2I + Cl(2P(1/2))/Cl(2P(3/2)), CHI + HCl, and CH2 + ICl. The state-selective detection of I(2P(3/2)) and I(2P(1/2)) has allowed the estimation of the branching ratio for I(2P(1/2)):I(2P(3/2)) to be 0.73:0.27. Theoretical calculations based on the time-dependent density-functional theory have been also made to investigate excited electronic potential-energy surfaces, plausible intermediates, and transition structures involved in these photodissociation reactions. The translation energy distributions derived from the TOF measurements suggest that at least two dissociation mechanisms are operative for these photodissociation processes. One involves the direct dissociation from the 2 1A' state initially formed by 193 nm excitation, leading to significant kinetic-energy releases. For the I-atom and Cl-atom elimination channels, the fragment kinetic-energy releases observed via this direct dissociation mechanism are consistent with those predicted by the impulsive dissociation models. Other mechanisms are likely predissociative or statistical in nature from the lower 1 1A' and 1 1A' states and/or the ground X 1A' state populated by internal conversion from the 2 1A' state. On the basis of the maximum kinetic-energy release for the formation of CH2Cl + I(2P(1/2)), we have obtained a value of 53+/-2 kcal/mol for the 0 K bond dissociation energy of I-CH2Cl. The intermediates and transition structures for the CHI + HCl and CH2 + ICl product channels have been also investigated by ab initio quantum calculations at the MP2(full)/6-311G(d) and B3LYP(full)/6-11G(d) levels of theory. The maximum kinetic-energy releases observed for the CHI + HCl and CH2 + ICl channels are consistent with the interpretation that the formation of CHI and CH2 in their ground triplet states is not favored.     

Computational studies on the kinetics and mechanisms for NH3 reactions with ClOx (x = 0-4) radicals

Kinetics and mechanisms for NH3 reactions with ClOx (x = 0-4) radicals have been investigated at the G2M level of theory in conjunction with statistical theory calculations. The geometric parameters of the species and stationary points involved in the reactions have been optimized at the B3LYP/6-311+G(3df,2p) level of theory. Their energetics have been further refined with the G2M method. The results show that the H-abstraction process is the most favorable channel in each reaction and the barriers predicted in decreasing order are OClO > ClO > Cl > ClO3 > ClO4. All reactions were found to occur by hydrogen-bonding complexes; the rate constants for these complex metathetical processes have been calculated in the temperature range 200-2000 K by the microcanonical VTST and/or RRKM theory (for ClO4 + NH3) with Eckart tunneling and multiple reflection corrections. The predicted rate constants are in good agreement with the available experimental data.     

Transition state structure and energetics of the N(2)O + X (X = Cl,Br) reactions

The structural and vibrational properties of the transition state of the N(2)O + X (X = Cl,Br) reactions have been characterized by ab initio methods using density functional theory. We have employed Becke's hybrid functional (B3LYP), and transition state optimizations were performed with 6-31G(d), 6-311G(2d,2p), 6-311+G(3d,2p), and 6-311+G(3df,2p) basis sets. For the chlorine atom reaction the coupled-cluster method (CCSD(T)) with 6-31G(d) basis set was also used. All calculations resulted in transition state structures with a planar cis arrangement of atoms for both reactions. The geometrical parameters of transition states at B3LYP are very similar, and the reaction coordinates involve mainly the breaking of the N-O bond. At CCSD(T)/6-31G(d) level a contribution of the O-Cl forming bond is also observed in the reaction coordinate. In addition, several highly accurate ab initio composite methods of Gaussian-n (G1, G2, G3), their variations (G2(MP2), G3//B3LYP), and complete basis set (CBS-Q, CBS-Q//B3LYP) series of models were applied to compute reaction energetics. All model chemistries predict exothermic reactions. The G3 and G2 methods result in the smallest deviations from experiment, 1.8 and 0 kcal mol(-1), for the enthalpies of reaction for N(2)O reaction with chlorine and bromine, respectively. The G3//B3LYP and G1 methods perform best among the composite methods in predicting energies of the transition state, with a deviation of 1.9 and 3.0 kcal mol(-1), respectively, in the activation energies for the above processes. However, the B3LYP/6-311+G(3df,2p) method gives smaller deviations of 0.4 and -1.0 kcal mol(-1), respectively. The performance of the methodologies applied in predicting transition state energies was analyzed.     

A theoretical study of the H-abstraction reactions from HOI by moist air radiolytic products (H, OH, and O (3P)) and iodine atoms (2P(3/2))

The rate constants of the reactions of HOI molecules with H, OH, O ((3)P), and I ((2)P(3/2)) atoms have been estimated over the temperature range 300-2500 K using four different levels of theory. Geometry optimizations and vibrational frequency calculations are performed using MP2 methods combined with two basis sets (cc-pVTZ and 6-311G(d,p)). Single-point energy calculations are performed with the highly correlated ab initio coupled cluster method in the space of single, double, and triple (pertubatively) electron excitations CCSD(T) using the cc-pVTZ, cc-pVQZ, 6-311+G(3df,2p), and 6-311++G(3df,3pd) basis sets. Reaction enthalpies at 0 K were calculated at the CCSD(T)/cc-pVnZ//MP2/cc-pVTZ (n = T and Q), CCSD(T)/6-311+G(3df,2p)//MP2/6-311G(d,p), and CCSD(T)/6-311++G(3df,3pd)//MP2/6-311G(d,p) levels of theory and compared to the experimental values taken from the literature. Canonical transition-state theory with an Eckart tunneling correction is used to predict the rate constants as a function of temperature. The computational procedure has been used to predict rate constants for H-abstraction elementary reactions because there are actually no literature data to which the calculated rate constants can be directly compared. The final objective is to implement kinetics of gaseous reactions in the ASTEC (accident source term evaluation code) program to improve speciation of fission products, which can be transported along the reactor coolant system (RCS) of a pressurized water reactor (PWR) in the case of a severe accident.     

Theoretical study and rate constant computation on the reaction HFCO + OH --> CFO + H2O

The potential energy surface, including the geometries and frequencies of the stationary points, of the reaction HFCO + OH is calculated using the MP2 method with 6-31+G(d,p) basis set, which shows that the direct hydrogen abstraction route is the most dominating channel with respect to addition and substitution channels. For the hydrogen abstraction reaction, the single-point energies are refined at the QCISD(T) method with 6-311++G(2df,2pd) basis set. The calculated standard reaction enthalpy and barrier height are -17.1 and 4.9 kcal mol(-1), respectively, at the QCISD(T)/6-311++G(2df,2pd)//MP2/6-31+G(d,p) level of theory. The reaction rate constants within 250-2500 K are calculated by the improved canonical variational transition state theory (ICVT) with small-curvature tunneling (SCT) correction at the QCISD(T)/6-311++G(2df,2pd)//MP2/6-31+G(d,p) level of theory. The fitted three-parameter formula is k = 2.875 x 10(-13) (T/1000)1.85 exp(-325.0/T) cm(3) molecule(-1) s(-1). The results indicate that the calculated ICVT/SCT rate constant is in agreement with the experimental data, and the tunneling effect in the lower temperature range plays an important role in computing the reaction rate constants.     

On the chaperon mechanism: application to ClO + ClO (+N2) --> ClOOCl (+N2)

The dynamics of the ClO + ClO (+N(2)) radical complex (or chaperon) mechanism is studied by electronic structure methods and quasi-classical trajectory calculations. The geometries and frequencies of the stationary points on the potential energy surface (PES) are optimized at the B3LYP/6-311+G(3df) level of theory, and the energies are refined at the CCSD(T)/6-311+G(3df) (single-point) level of theory. Basis set superposition error (BSSE) corrections are applied to obtain 1.5 kcal mol(-1) for the binding energy of the ClO.N(2) van der Waals (VDW) complex. A model PES is developed and used in quasi-classical trajectory calculations to obtain the capture rate constant and nascent energy distributions of ClOOCl* produced via the chaperon mechanism. A range of VDW binding energies from 1.5 to 9.0 kcal mol(-1) are investigated. The anisotropic PES for the ClO.N(2) complex and a separable anharmonic oscillator approximation are used to estimate the equilibrium constant for formation of the VDW complex. Rate constants, branching ratios to produce ClOOCl, and nascent energy distributions of excited ClOOCl* are discussed with respect to the VDW binding energy and temperature. Interestingly, even for weak VDW binding energies, the N(2) usually carries away enough energy to stabilize the nascent ClOOCl*, although the VDW equilibrium constant is small. For stronger binding energies, the stabilization efficiency is reduced, but the capture rate constant is increased commensurately. The resulting rate constants for forming ClOOCl* from the title reaction are only weakly dependent on the VDW binding energy and temperature. As a result, the relative importance of the chaperon mechanism is mostly dependent on the VDW equilibrium constant. For the calculated ClO.N(2) binding energy of 1.5 kcal mol(-1), the VDW equilibrium constant is small, and the chaperon mechanism is only important at very high pressures.     

Ab initio studies of CIO, reactions: prediction of the rate constants of CIO + NO for the forward and reverse processes.

The mechanisms for ClO+NO and its reverse reactions were investigated by means of ab initio molecular orbital and statistical theory calculations. The species involved were optimized at the B3LYP/6-311 +G(3df) level, and their energies were refined at the CCSD(T)/6-311+ G(3df)//B3LYP/6-311 + G(3df) level. Five isomers and the transition states among them were located. The relative stability of these isomers is ClNO2 > cis-ClONO > trans-ClONO > cis-OClNO>trans-OClNO. The heats of formation of the three most-stable isomers were predicted using isodesmic reactions by different methods. The predicted bimolecular reaction rate constant shows that, below 100 atm, the formation of Cl+NO2 is dominant and pressure-independent. The total rate constant can be expressed as: k(ClO+NO)= 1.43 x 10(-9)T(-083)exp(92/ T) cm3 molecule(-1)s(-1) in the temperature range of 200-1000 K, in close agreement with experimental data. For the reverse reaction, Cl+NO2-->ClNO2 and ClONO (cis and trans isomers), the sum of the predicted rate constants for the formation of the three isomers and their relative yields also reproduce the experimental data well. The predicted total third-order rate constants in the temperature range of 200-1000 K can be represented by: k0(He) = 4.89 x 10(-6)T(-5.85) exp(-796/T) cm6 molecule(-1)s(-1) and k0(N2) =5.72 x 10(-15)T(-5.80) exp(-814/T) cm6 molecule(-1)s(-1). The predicted high- and low-pressure limit decomposition rates of CINO2 in Ar in the temperature range 400-1500 K can be expressed, respectively, by: k-(ClNO2) = 7.25 x 10(19)T(-1.89) exp(-16875/T) s(-1) and kd(ClNO2) = 2.51 x 10(38)T(-6.8) exp(-18409/T) cm3 molecule(-1) s(-1). The value of k0(ClNO2) is also in reasonable agreement with available experimental data.     

B2F2分子异构体结构的量子化学计算研究

采用abinitio方法对B₂F₂分子异构体的结构进行了计算研究,并与Al₂F₂分子进行了比较.结果表明,B₂F₂具有D_∞h对称性,³Σ_g电子态的直线型结构FBBF是B₂F₂分子的最稳定异构体,对文献的结果进行了修正.在UCCSD(T,full)/6－311+G(2df)水平下,F-B和B-B键长分别为0.12942和0.14820nm,振动频率出现在1860.00和1320.62cm^-1处.在UQCISD(T,full)/6－311+G(2df)//UMP2(full)/6－311+G(d)+ZPVE水平下,³Σ_g态的线性FBBF分子的垂直电离势为848.58kJ/mol,而由³Σ_g电子态的BF二聚为³Σ_g态的线性FBBF分子的焓变为59.86kJ/mol,此二聚化反应是放热反应,说明二聚化过程在能量上是有利的.     

HNCS与Cl原子的反应机理及电子密度拓扑分析

在QCISD(T)/6-311++G(d,p)和B3LYP/6-311++G(d,p)级别上研究了HNCS与Cl原子的反应机理. 并应用经典过渡态理论和正则变分过渡态理论结合小曲率隧道效应, 计算了200-2500 K温度范围内各反应通道的速率常数. 结果表明, HNCS与Cl原子反应存在3个反应通道. 当温度低于294 K时, 生成HCl+NCS的夺氢反应(a)是优势通道, 温度高于294 K时, 生成HNC(Cl)S的加成反应(c)为主反应通道, Cl进攻N的反应通道(b)因能垒较高而难以进行.     

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.     

Theoretical Study on Gas Phase Reactions of OH Hydrogen-Abstraction from Formyl Fluoride with Di erent Catalysts

The mechanisms and kinetics of the gas phase reactions that the hydrogen atom in formyl uoride (FCHO) abstracted by OH in the presence of water, formic acid (FA), or sulfuric acid (SA) are theoretically investigated at the CCSD(T)/6-311++G(3df, 3pd)//M06-2X/6-311++G(3df, 3pd) level of theory. The calculated results show that the barriers of the transition states involving catalysts are lowered to -2.89, -6.25, and -7.76 kcal/mol from 3.64 kcal/mol with respect to the separate reactants, respectively, which re ects that those catalysts play an important role in reducing the barrier of the hydrogen abstraction reaction of FCHO with OH. Additionally, using conventional transition state theory with Eckart tun-neling correction, the kinetic data demonstrate that the entrance channel X FCHO+OH (X=H₂O, FA, or SA) is signi cantly more favorable than the pathway X OH+FCHO. More-over, the rate constants of the reactions of FCHO with OH radical with H₂O, FA, or SA introduced are computed to be smaller than that of the naked OH+FCHO reaction because the concentration of the formed X FCHO or X OH complex is quite low in the atmosphere.