The DFT and ab initio investigations on the mechanism of CF3O2 + H reaction

The mechanisms for the reaction of CF₃O₂ with atomic hydrogen were studied with ab initio and DFT methods. The results reveal that the reaction could take place on the singlet and triplet potential energy surfaces (PES). For the singlet PES, addition/elimination and substitution mechanisms are determined, and the former one is dominant. The most favorable channel involves the association of CF₃O₂ with H atom to form CF₃O₂H (IM1) via a barrierless process, and then the O–O bond dissociates to give out CF₃O + OH. The secondary product might be CF₃OH + O, formed from the O–O bond cleavage in the initial adduct CF₃O(H)O (IM2). Other products such as CF₃ + O₂H, HF + CF₂O₂ and O₂ + CHF₃ are of no importances because of higher barriers. On the triplet PES, only substitution mechanism is located. With higher barriers involving, the channels on the triplet PES could be negligible compared with the channels on the singlet PES.     

DFT and ab initio theoretical study for the CF3S + CO reaction

The potential energy surface for the reaction of CF₃S with CO is calculated at the G4//B3LYP/6-311++G(d,p) level of theory. The results show that F-abstraction and addition-elimination mechanisms are involved, and the latter one is dominant thermodynamically and kinetically. The dominant channel is the reactant addition to form CF₃SCO, and then decomposes to CF₃ + OCS. While the direct F-abstraction channel and CF₃SCO isomerization channel are not significant for the title reaction due to higher barriers involved. The comparisons among four reactions of CX₃Y + CO (X = H, F; and Y = O, S) are made to imply the similar and different properties and reactivities of the same family elements and the F- and S-substituted derivatives.     

Singlet-triplet energy separations in divalent five-membered cyclic conjugated C5H3X, C4H3SiX, C4H3GeX, C4H3SnX, and C4H3PbX (X = H, F, Cl, and Br)

Singlet-triplet energy gaps in cyclopenta-2,4-dienylidene, as well as its 2- or 3-halogenated derivatives, are compared and contrasted with their sila, germa, stana, and plumba analogues; at HF/6-31G* and B3LYP/ 6-311++G(3df, 2p) levels of theory. Energy gaps (ΔG_t-s), between triplet (t) and singlet (s) states, appear linearly proportional to: (i) the size of the group 14 divalent element (M = C, Si, Ge, Sn and Pb), (ii) the angle ∠C-M-C, and (iii) the ΔG_(LUMO-HOMO) of the singlet state involved. The magnitude of ΔG_t-s, for each 2- and/or 3-substituted species studied, increases with an order of: carbenes < silylenes < germylenes < stanylenes < plumbylenes. This order reverses for the barriers of the ring puckering. The puckering occurs with more ease for every singlet, compared to its corresponding triplet form.Regardless of the group 14 element (M) employed, every 3-halo-substituted species is more stable than the corresponding 2-halo-substituted isomer. For M = Pb, Sn and/or Ge; 3-halo-substituted species have higher ΔG_t-s than their corresponding 2-halo-substituted analogues. For M = Si, similar ΔG_t-s are found for 2- and 3-halogenated isomers. For M = C, 3-halo-substituted species have lower ΔG_t-s than their corresponding 2-halo-substituted analogues.Every cyclic singlet has a larger ∠C-M-C angle, than its corresponding cyclic triplet state, except for 3-halosilacyclopenta-2,4-dienylidenes where triplet has a larger ∠C-M-C angle than its corresponding singlet state.     

A quest for triplet silylenes XHSi3 at ab initio and DFT levels (X = H, F, Cl and Br)

Four ground state triplet silylenes are found among 30 possible silylenic XHSi₃ structures (X = H, F, Cl and Br), at seven ab initio and DFT levels including: B3LYP/6-311++G^∗∗, HF/6-311++G^∗∗, MP3/6-311G^∗, MP2/6-311+G^∗∗, MP4(SDTQ)/6-311++G^∗∗, QCISD(T)/6-311++G^∗∗ and CCSD(T)/6-311++G^∗∗. The latter six methods indicate that the triplet states of 3-flouro-1,2,3-trisilapropadienylidene, 1-chloro-1,2,3-trisilapropargylene and 3-chloro-1,2,3-trisilapropargylene are energy minima. These triplets appear more stable than their corresponding singlet states which cannot even exist for showing negative force constants. Also, triplet state of 1-flouro-1,2,3-trisilapropargylene is possibly accessible for being an energy minimum, since its corresponding singlet state is not a real isomer. Some discrepancies are observed between energetic and/or structural results of DFT vs. ab initio data.     

Theoretical study on the reaction mechanism of CH4 with CaO

The reaction pathways and energetics for the reaction of methane with CaO are discussed on the singlet spin state potential energy surface at the B3LYP/6-311+G(2df,2p) and QCISD/6-311++G(3df,3pd)//B3LYP/6-311+G(2df,2p) levels of theory. The reaction of methane with CaO is proposed to proceed in the following reaction pathways: CaO + CH₄ → CaOCH₄ → [TS] → CaOH + CH₃, CaO + CH₄ → OCaCH₄ → [TS] → HOCaCH₃ → CaOH + CH₃ or [TS] → CaCH₃OH → Ca + CH₃OH, and OCaCH₄ → [TS] → HCaOCH₃ → CaOCH₃ + H or [TS] → CaCH₃OH → Ca + CH₃OH. The gas-phase methane–methanol conversion by CaO is suggested to proceed via two kinds of important reaction intermediates, HOCaCH₃ and HCaOCH₃, and the reaction pathway via the hydroxy intermediate (HOCaCH₃) is energetically more favorable than the other one via the methoxy intermediate (HCaOCH₃). The hydroxy intermediate HOCaCH₃ is predicted to be the energetically most preferred configuration in the reaction of CaO + CH₄. Meanwhile, these three product channels (CaOH + CH₃, CaOCH₃ + H and Ca + CH₃OH) are expected to compete with each other, and the formation of methyl radical is the most preferable pathway energetically. On the other hand, the intermediates HCaOCH₃ and HOCaCH₃ are predicted to be the energetically preferred configuration in the reaction of Ca + CH₃OH, which is precisely the reverse reaction of methane hydroxylation.     

Ab initio quantum chemical studies of the reactions of CF3CFHO2 with HO2

The potential energy surface (PES) for the CF₃CFHO₂+HO₂ reaction has been theoretically investigated using the DFT [B3LYP/6‐311G(d,p)] and B3LYP/6‐311++G(3df,3pd)//B3LYP/6‐311G(d,p) levels of theory. Both singlet and triplet PESs are investigated. The reaction mechanism on the triplet surface is simple. It is revealed that the formation of CF₃CFHOOH+³O₂ is the dominant channel on the triplet surface. On the basis of the ab initio data, the total rate constants for the reaction CF₃CFHO₂+HO₂ in the T = 210–500 K range have been computed using conventional transition state theory with Wigner's tunneling correction and have been fitted by a rate constant expression as k = 1.04 ×10^?12(cm³ molecule^?1 s^?1) exp (700.33/T). Calculated transition state rate constants with Wigner's tunneling correction for the reaction CF₃CFHO₂+HO₂ are in good agreement with the available experimental values. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

Direct dynamics studies for the reactions of CF3CHFCF3 and CF3CF2CHF2 with H atoms

The rate constants of the hydrogen abstraction reactions of CF₃CHFCF₃ + H (R1) and CF₃CF₂CHF₂ + H (R2) have been calculated by means of the dual-level direct dynamics method. Optimized geometries and frequencies of stationary points and extra points along the minimum-energy path (MEP) are obtained at the MPW1K/6-311+G(d,p) level, and the classical energetic information is further corrected with the interpolated single-point energy (ISPE) approach by the G3(MP2) level of theory. Using the canonical variational transition state theory (CVT) with small-curvature tunneling corrections (SCT), the rate constants are evaluated over a wide temperature range of 200-2000 K. The calculated CVT/SCT rate constants are in good agreement with available experimental values. It is found that the variational effect is very small and almost negligible over the whole temperature region. However, the small-curvature tunneling correction plays an important role in the lower temperature range. Furthermore, the heats of formation of species CF₃CF₂CHF₂ (SC₁ or SC₂) and CF₃CF₂CF₂ are studied using isodesmic reactions to further elucidate the thermodynamic properties.     

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.     

Switching of global minima of novel germylenic reactive intermediates via halogens (X): C2GeH2 vs. C2GeHX at ab initio and DFT levels

For 30 C₂GeHX germylenic isomers, one cyclic structure, X-germacyclopropenylidene, and three acyclics are considered, which include: ethynyl-X-germylene, X-vinylidenegermylene, and (X-ethynyl)germylene (X = H, F, Cl, and Br). The global minimum among six isomeric C₂GeH₂ (where X = H), is found to be cyclic, aromatic, singlet germacyclopropenylidene. In contrast, among the 24 corresponding halogermylenes, C₂GeHX (where X = F, Cl, and Br), the global minima switch to acyclic, singlet ethynylhalogermylenes, at eight reasonably high ab initio and DFT levels. The direct resonance interaction between X and the divalent center Ge in the singlet acyclic ethynylhalogermylene structures, is claimed to justify switching of the calculated global minima in the halo derivatives. GIAO-NICS calculations indicate that the X-germacyclopropenylidene isomer is more aromatic for X = H than X = F, Cl, or Br. The angle ∠XGeC bending potential energy curves show the singlet and triplet ethynylgermylene crossing at ≈146°, for X = H.     

Theoretical study on the reaction mechanisms of CH2SH + NO reaction

The mechanisms for the CH₂SH + NO reaction were investigated on both of the singlet and triplet PES at the BMC-CCSD//B3LYP/6-311+G(d,p) level. The results indicate that the singlet PES is much lower than the triplet PES energetically; therefore, the reaction occurs on the singlet PES dominantly. The most favorable channel on the singlet PES takes place by a barrierless addition of N atom to CH₂SH radical to form HSCH₂NO. Subsequently, the rearrangement of the initial adduct HSCH₂NO (IM1) to form another intermediate IM3 via a four-center transition state, followed by the C–O bond fission in IM3 leading to the major product CH₂S + HNO. Due to high barriers, other product including HC(N)SH + HO, HON + CH₂S, and HNO + CHSH could be negligible. The direct abstraction channel was also determined to yield CH₂S + HON. With high barrier (33.3 kcal/mol), it is not competitive with the addition channel, in which all stationary points are lower than reactant energetically. While on the triplet PES, with the lowest barrier height (18.8 kcal/mol), the direct N-abstracted channel to form CH₂S + HNO is dominant. However, it is not competitive with the channels on the singlet PES. Our results are in good accordance with experimental conclusions that the reaction proceeds via addition mechanism.     

Theoretical study on the potential energy surface of NC<Subscript>3</Subscript>P isomers

Density functional theory (DFT) calculations have been used to study the isomerization process in the NC₃P system. At the DFT/B3LYP/6-311G(d) level, 28 triplet and 28 singlet minima were obtained on their respective potential energy surfaces. The linear triplet ³NCCCP is the lowest-energy structure among the isomers. On the triplet PES, only linear isomers ³NCCCP, ³CNCCP, ³CCCNP, and ³CCNCP possess great kinetic and thermodynamic stabilities to exist under low-temperature conditions (such as in the dense interstellar clouds). At the same time, one chain-like and four three-membered-ring isomers on the singlet PES have been located with high kinetic and thermodynamic stabilities. Further CCSD(T)/6-311G(2df)//QCISD/6-311G(d), CCSD(T)/cc-pVTZ//DFT/B3LYP/cc-pVTZ, and CASPT2(14,12)/cc-pVQZ//CASSCF(14,12)/cc-p VQZ calculations are performed on the structures, frequencies, and energies of the relevant species. The bonding natures were analyzed and the results were compared with the analogous NC₃N and NC₂P molecules so as to aid their future experimental or astrophysical detection.     

Ab initio and density functional theory study of the enthalpies of formation of F2SOx and FClSOx (x=1, 2)

Enthalpies of formation of F₂SO, F₂SO₂, FClSO and FClSO₂ molecules have been determined using ab initio molecular orbital theory and density functional theory (DFT) calculations. Different DFT approaches and levels of the Gaussian-3 and the complete basis set (CBS) ab initio model chemistries have been employed to calculate enthalpies of formation from both total atomization energies and isodesmic reaction schemes. The best values at 298 K for F₂SO, F₂SO₂, FClSO and FClSO₂ as derived from an average of G3, G3B3, CBS-Q and CBS-QB3 isodesmic energies are −140.6, −181.1, −92.6 and −132.3 kcal mol⁻¹, respectively. The results obtained suggest that the accumulated small component errors found in the DFT-based methods are significantly reduced at the ab initio levels employed. Structural properties, harmonic vibrational frequencies, mode assignations and infrared intensities derived from B3LYP and mPW1PW91 functional with the 6-311+G(3df) basis set are presented.     

Ab initio/DFT theory and multichannel RRKM study on the mechanisms and kinetics for the CH3S + CO reaction

The potential energy surface for the reaction of CH₃S with CO was calculated at the G3MP2//B3LYP/6-311++G(d,p) level. The rate constants for feasible channels leading to several products were calculated by TST and multichannel-RRKM theory. The results show that addition–elimination mechanism is dominant, while hydrogen abstraction mechanism is uncompetitive. The major channel is the addition of CO to CH₃S leading to an intermediate CH₃SCO which then decomposes to CH₃ + OCS. In the temperature range of 200–3000 K, the overall rate constants are positive temperature dependence and pressure independence, and it can be described by the expression as k = 1.10 × 10⁻¹⁶T^1.57exp(−3359/T) cm³ molecule⁻¹ s⁻¹. At temperature between 208 and 295 K, the calculated rate constants are in good agreement with the experimental upper limit data. At T = 1000 and 2000 K, the major product is CH₃ + OCS at lower pressure; while at higher pressure, the stabilization of IM1 is dominant channel.     

Theoretical investigation on the reaction of HS+ with CH3NH2

The singlet and triplet potential energy surfaces for the reaction of HS⁺ with the simplest primary amine, CH₃NH₂, were determined at the CCSD(T)/6-311+G(d,p) level using the B3LYP/6-311G(d,p) and QCISD/6-311G(d,p) geometries. All possible reaction channels were explored. The results show that three paths on the singlet potential energy surface and one path on the triplet potential energy surface are competitive. These four feasible paths provide products which are presented in the paper and they are consistent with previous experimental results. On the other hand, the stationary points involved in the most favourable path all lie below those of the reactant and thus the title reaction is expected to be rapid, which is also consistent with the experiment.     

Theoretical study on the (O2-HF) complex

The potential energy surface of the (O₂-HF)⁺ complex has been investigated theoretically. Equilibrium structure has been calculated at the UCCSD/6-311++G(2d,2p) and UCCSD/aug-cc-pVTZ levels of theory. Two stable complexes, T-shape and collinear, have been found on the surface whose interaction energies are calculated to be −12.36 and −10.19 kcal mol⁻¹, respectively, at the UCCSD(T)/aug-cc-pVTZ level after correcting for the basis set superposition error (BSSE). Natural bond orbital (NBO) analysis revealed intermolecular charge transfers occur followed by intramolecular charge rearrangement. A large contribution from the n_F to π* _OO is found.     

全氟烷基磺酸酯C—O键断裂的同面SN2反应

用密度泛函理论研究了CF₃SO₃CF₂CF₃＋F^－的碳氧键断裂反应的机理. 首先, 用DFT方法优化了反应物、中间体、过渡态、产物的平衡构型, 分析了碳氧键断裂反应的势能面变化. 发现在S_N2反应机理中, 除了S—O断裂S_N2反应外, 引起C—O键断裂的同面进攻也是一个可能的反应途径. 理论计算表明, 最终反应的产物是受热力学控制的, S—O键的断裂绝对地优于C—O的断裂. 因此, C—O断裂的同面机理虽然是可能的, 但却难以被实验观察到. 本文还讨论了端基 —F₃在同面S_N2反应中的邻位效应, 以及基组对这个效应的影响.     

Ab initio direct dynamics studies on the hydrogen abstraction reactions of CF2H2 and CF3H with F atom

The dynamic properties of the hydrogen abstraction reactions of CF₂H₂ and CF₃H with F atom are investigated in the temperature range of 182–2000 K. The minimum-energy path (MEP) is optimized at MP2/6-311 G(d, p) level, then the energy profiles are refined at the CCSD(T)/6-311++G(3df, 2pd) level (single-point). The theoretical rate constants, which are calculated by the variational transition state theory (VTST) including the small curvature tunneling (SCT) correction, are in good agreement with the experimental ones. It is found that the rate constant of the CF₂H₂ + F reaction are larger than that of the CF₃H + F reaction and the activation energies exhibit in the just opposite order. This phenomenon can be rationalized by the hardness η of the halomethane molecules. The comparison of the two reactions with the CFH₃ + F reaction is made. It is found that the rate constants decrease in the order of CFH₃ + F > CF₂H₂ + F > CF₃H + F. The effect of fluorine substitution leads to a dramatic increase in the activation energy and a decrease in the preexponential factor. We hope that present theoretical studies for these compounds can give further information concerning how fluorine substitution affects the rate constants of hydrogen abstraction reactions.     

Gas-phase thermolysis reaction of formaldehyde diperoxide. Kinetic study and theoretical mechanisms

Gas-phase thermolysis reaction of formaldehyde diperoxide (1,2,4,5-tetroxane) was performed in an injection chamber of a gas chromatograph at a range of 463-503 K. The average Arrhenius activation energy and pre-exponential factor were 29.3 ± 0.8 kcal/mol and 5.2 × 10¹³ s⁻¹, respectively. Critical points and reaction paths of the ground singlet and first triplet potential energy surfaces (PES) were calculated, using DFT method at BHANDHLYP/6-311+G^∗∗ level of the theory. Also, G3 calculations were performed on the reactant and products. Reaction by the ground-singlet and first-triplet states turned out to be endothermic and exothermic, respectively. The mechanism in three steps seemed to be the most probable one. An electronically non-adiabatic process appeared, in which a crossing, at an open diradical structure, from the singlet to the triplet state PES occurred, due to a spin-orbit coupling, yielding an exothermic reaction. Theoretical kinetic constant coming from the non- adiabatic transition from the singlet to the triplet state agrees with the experimental values.     

Computational mechanistic investigation of the gas phase C2H4 + CO reaction on the singlet and triplet potential energy surfaces

Reaction pathways of ethylene and carbon monoxide on the singlet and triplet potential energy surfaces (PESs) have been calculated at B3LYP/6-311++G (3df, 3dp), G3B3 and CCSD(T)//B3LYP levels. Reaction mechanisms have been investigated by analysis of various structures. Suggested reaction mechanisms reveal that ³P3(CH₂CHCHO) and ³P4(CH₃CCHO) are thermodynamically stable adducts with the negative value in Gibbs free energies on the triplet PES. In addition, results show that one intersystem crossing exists between triplet and singlet PESs, which are obtained by scanning of the C–C bond length in ¹IN3 and ³IN7 species.     

Effects of acid treatment on activated carbon used as a support for Rb and K catalyst for C2F5I synthesis and its mechanism

In the present work, the activated carbon (AC) support was treated with HCl, HNO₃ and HF solution. The order of catalyst dispersion was as follows: Rb-K/AC-HNO₃ > Rb-K/AC-HF > Rb-K/AC-HCl > Rb-K/AC. The same sequence was also observed for the amount of the acid surface oxygen groups on AC, but not for the basicity of the catalyst. The key role of acid treatment on AC surface chemistry and the basic sites, which are closely related to catalyst dispersion and basicity, is examined to rationalize these findings. On the other hand, a consideration of the reaction mechanism suggests that the reaction proceeds via CF₂ carbenes formed on the catalyst surface as intermediates, followed by carbine disproportionation to CF₃ radicals and CF₃CF₂ radicals, followed by reaction with I₂ to produce CF₃CF₂I, and it was also found that the Rb-K/AC-HCl catalyst with a high dispersion and moderate basicity was helpful for the enhancement of catalytic activity for C₂F₅I synthesis.