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.     

Mechanism and rate constants for ammonia photochemical oxidation with O2 on the singlet and triplet potential energy surfaces

The mechanism and thermodynamic of NH₃ + O₂ reaction on the singlet and triplet potential energy surfaces (PES), were carried out using the RMP2 and CCSD (T)//RMP2 theoretical approaches in connection with the 6-311++G(d, p) basis set. Three pre-reactive complexes, ¹C1, ¹C2, and ³C1 on the singlet and triplet PES were formed between ammonia and molecular oxygen. With variety of pre-reactive complexes, six types of products are obtained, of which two types are found to be thermodynamically stable. The mechanistic properties of all products channels are discussed. Results show that production of HONO + H₂ and HN(OH)₂ are the main reaction channels in thermodynamic viewpoint with the Gibbs free energy of ?G° = ?34.681 and ?27.153 kcal/mol, respectively. Rate constants of the title reaction over the temperature range of (200–1000 K) show kinetic products are different from thermodynamic products.     

Computational study on the reaction mechanism of the gas-phase atom-negative ion of S?+?NO2 ?: comparative study of mechanism with S?+?O3 reaction as isoelectronic and isostructure systems

The reaction mechanism of sulfur vapor (S) with nitrite ion (NO₂ ⁻) has been investigated theoretically on the triplet and singlet potential energy surfaces (PESs). All stationary points for the title reaction have been optimized at the B3LYP/6-311+G(3df) level. The energetic data have been obtained at the CCSD(T)//B3LYP level employing the 6-311+G(3df) basis set. Five stable collision complexes, ³IN1 (S–ONO⁻), ³IN2 (cyclic SONO⁻), ¹IN1 (cis S–ONO⁻), ¹IN2 (S–NO₂ ⁻), and ¹IN3 (trans S–ONO⁻), have been considered on the triplet and singlet PESs through barrier-less and exothermic processes. By starting from these complexes, a simple mechanism has been obtained on the triplet PES while a complex mechanism has been considered on the singlet PES. The calculated results show that there are no favorable paths for the reaction of S with NO₂ ⁻ on the singlet PES. Therefore, the S + NO₂ ⁻ reaction proceeds only on the triplet PES to produce ³SO + ³NO⁻ as main products. The results from the comparative study of S + NO₂ ⁻ reaction mechanism with S + O₃ (as isoelectronic and isostructure reactions) on the singlet PES show similarities in the overall trend of reaction mechanism and atom connectivity and differences in the stability of intermediates and the energy barriers of transition states.     

A barrier-free molecular radical-molecule reaction: {^{3}C_{2} (a^{3}\Pi) {+} O_{2} (X^{3} \Sigma)}

The reaction of ³C₂ (a³Π) radical with O₂ (X³Σ) molecule has been studied theoretically using ab initio Quantum Chemistry method. Both singlet and triplet potential energy surfaces (PES) are calculated at the CCSD(T)/aug-cc-pVDZ//B3LYP/6-311+G(d) + ZPE and G3B3 levels of theory. On the singlet PES of the title reaction, it is shown that the most feasible pathway should be the O-atom of O₂ attacking the C-atom of the  ³C₂ molecule first to form the adduct 1 CCOO, followed by the O-shift to give intermediate 2 CC(OO), and then to the major products P1 (2CO). Alternatively, 1 can be directly dissociated to P1 via transition state TS1-P1. The other reaction pathways are less competitive due to thermodynamical or kinetic factors. On the other hand, the pathways on the triplet PES are less competitive than those on the singlet PES in low temperature range, whereas it is not the case in high temperature ranges. On the basis of the analysis of the kinetics of all pathways through which the reactions proceed, we expect that the competitive power of reaction pathways may vary with experimental conditions for the title reaction. The reaction heats of formation calculated are in good agreement with that obtained experimentally.     

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.     

Theoretical study on the reaction mechanism of CH2SH + NO2

The mechanisms of CH₂SH with NO₂ reaction were investigated on the singlet and triplet potential energy surfaces (PES) at the BMC-CCSD//B3LYP/6-311 + G(d,p) level. The result shows that the title reaction is more favourable on the singlet PES thermodynamically, and it is less competitive on the triplet PES. On the singlet PES, the initial addition of CH₂SH with NO₂ leads to HSCH₂NO₂ (IM2) without any transition state, followed by a concerted step involving C–N fission and shift of H atom from S to O giving out CH₂S + trans-HONO, which is the major products of the title reaction. With higher barrier height, the minor products are CH₂S + HNO₂, formed by a similar concerted step from the initial adduct HSCH₂ONO (IM1). The direct abstraction route of H atom in SH group abstracted by O atom might be of some importance. It starts from the addition of the reactants to form a weak interaction molecular complex (MC3), subsequently, surmounts a low barrier height leading to another complex (MC2), which gives out CH₂S + trans-HONO finally. Other direct hydrogen abstraction channels could be negligible with higher barrier heights and less stable products.     

Ab initio quantum chemical studies on the reactions of CF3O2 with OH

The potential energy surface for the CF₃O₂ + OH reaction has been theoretically investigated using the DFT (B3LYP/6-311G(d,p)) level of theory. Both singlet and triplet potential energy surfaces are investigated. The reaction mechanism on the triplet surface is simple. However, the reaction mechanism on the singlet surface is more complicated. It is revealed that the formation of CF₃O + HO₂ is the dominant channel on the triplet surface. The potential energy surface (PES) for this reaction has been given according to the relative energies calculated at the DFT/B3LYP/6-311G(d,p) level. Because this reaction involves both triplet and singlet states, triplet–singlet intersystem crossing (ISC) crossing also have been investigated in this paper.     

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     

The DFT study on mechanisms for NCO with C2H5 reaction

A detailed investigation has been performed at the QCISD(T)/6‐311++G(d,p)//B3LYP/6‐311+G(d,p) level for the reaction of NCO with C₂H₅ by constructing singlet and triplet potential energy surfaces (PES). The results show that the title reaction is more favorable on the singlet PES than on the triplet PES. On the singlet PES, the initial addition processes are barrierless and release lots of energy. The dominant channel occurs via the fragmentations of the initial adduct C₂H₅NCO and C₂H₅OCN to form C₂H₄ + HNCO and HOCN, respectively. With higher barrier heights, other products such as CH₄ + HNC + CO, CH₃CHNH + CO, CH₃CH + HNCO, and CH₃CN + H₂ + CO are less competitive. On the triplet PES, the entrance reactions surpass significant barriers; therefore, it could be negligible at the normal atmospheric condition. However, the most feasible channel on the triplet PES is the direct hydrogen abstraction channel to form CH₂CH₂ + HNCO. © 2010 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

Ab Initio Study on the Potential Energy Surfaces of B4 Cluster

The singlet and triplet potential energy surfaces (PES) for the isomerization and dissociation reactions of B₄ isomers have been investigated using ab initio methods. Ten B4 isomers have been identified and of these 10 species, 4 have not been reported previously. The singlet rhombic structure ¹1 is found to be the most stable on the B₄ surface, in agreement with the results of previous reports. Several isomerization and dissociation pathways have been found. On the singlet PES, the linear ¹3b can rearrange to rhombus ¹1 directly, while ¹3c rearranges to ¹1 through two‐step reactions involving a cyclic intermediate. On the triplet PES, the capped triangle structure ³2 undergoes ring opening to the linear isomer ³3b with a barrier of 34.8 kcal/mol and 44.9 kcal/mol, and the latter undergoes ring closure to the square structure ³1 with a barrier of 30.4 kcal/mol and 33.0 kcal/mol at the MP4/6–311+G(3df)//MP2/6–311G(d) and CCSD/aug‐cc‐pVTZ//MP2/6–311G(d) levels of theory, respectively. The direct decomposition of singlet B₄ yielding to B₃+B is shown to have a large endothermicity of 87.3 kcal/mol (CCSD), and that producing 2B₂ to have activation energy of 133.4 kcal/mol (CCSD).     

Theoretical study on the radical–radical reaction mechanism of CH2I + NO2

The mechanisms of CH₂I with NO₂ reaction were investigated on the singlet and triplet potential energy surfaces (PESs) by the UB3LYP method. The energetic information is further refined at the UCCSD(T) and UQCISD(T) levels of theory. Our results indicated that the title reaction is more favorable on the singlet PES thermodynamically, and less competitive on the triplet one. On the singlet PES, the title reaction is most likely to be initiated by the carbon-to-oxygen approach forming the adduct IM1 (H₂ICONO-trans) without any transition state, which can isomerizes to IM2 (H₂ICNO₂) and IM3 (H₂ICONO-cis), respectively. The most feasible pathway is the 1, 3-I shift with C–I and O–N bonds cleavage along with the N–I bond formation of IM1 lead to the product P1 (CH₂O + INO), which can further dissociate to give P3 (CH₂O + I + NO). The competitive pathway is 1, 3-H shift associated with O–N bond rupture of IM1 to form P2 (CHIO + HNO). The theoretically obtained major product CH₂O and adducts IM1 and IM2 are in good agreement with the kinetic detection in experiment. The similarities and discrepancies between CH₂I + NO₂ and CH₂Br + NO₂ reactions are discussed in terms of the electronegativity of halogen atom and the barrier height of the rate-determining process. The present study may be helpful for further experimental investigation of the title reaction.     

Theoretical study on the atmospheric reaction of CH3O2 with OH

A quantum chemical investigation on the reaction mechanism of CH₃O₂ with OH has been performed. Based on B3LYP and QCISD(T) calculations, seven possible singlet pathways and seven possible triplet pathways have been found. On the singlet potential energy surface (PES), the most favorable channel starts with a barrierless addition of O atom to CH₃O₂ leading to CH₃OOOH and then the O? O bond dissociates to give out CH₃O + HO₂. On the triplet PES, the calculations indicate that the dominant products should be ³CH₂O₂ + H₂O with an energy barrier of 29.95 kJ/mol. The results obtained in this work enrich the theoretical information of the title reaction and provide guidance for analogous atmospheric chemistry reactions. © 2015 Wiley Periodicals, Inc.     

Theoretical study of [Si,O,C,O] species: Prediction of new species on triplet potential energy surface

The intermediates [Si,O,C,O] of the Si + CO₂ reaction have been studied in detail using high level ab iniitio methods. Both singlet and triplet [Si,O,C,O] species are characterized structurally and energetically. On the singlet potential energy surface (PES), the vdw‐OSi–CO isomer and in the triplet PES, the bent‐SiOCO isomer is found to be thermodynamically as well as kinetically most stable species. All possible isomerization transition states (TS) are located on both singlet and triplet potential surfaces. On the triplet surface, the stability of the trans‐OSiCO isomer is comparable with that of the bent‐SiOCO isomer. A non‐planar cis‐SiOCO isomer is located on the triplet PES, which is predicted for the first time. Heats of formation at 0 K (Δ_fH°, 0 K) for all singlet and triplet species are computed using G3B3, G3MP2, and CBS‐Q theories. The discrepancy between G3B3 and the other two methods for the heat of formation value for triplet trans‐OSiCO is discussed. The PESs for singlet as well as triplet species with their dissociation asymptotes are explored at the CCSD(T)/6‐311G(d,p)//MP2/6‐311G(d,p) level of theory. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

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.     

Copper complexes with N-aminotriazolethione azomethines: Structures and magnetochemical properties

Azomethine derivatives of 4-amino-1,2,4-triazole-3-thiones (H₂L) and their metal complexes were obtained. The stabilities of the ligand conformers were calculated using quantum-chemical techniques. Ab initio (B3LYP/LANL2DZ) calculations of the complexes in the lower singlet and triplet states were performed with full geometry optimization. The structures and magnetochemical properties of the chelate complexes obtained were examined over a wide temperature range.     

Ab initio MO study of potential energy surface of NH2 with CN reaction

An extensive quantum chemical study of the potential energy surfaces (PES) for the association reaction of NH₂ with CN and the subsequent isomerization and dissociation reactions has been carried out using density functional theory (DFT)/B3LYP/6‐311++G(3df,2p) level of theory on both singlet and triplet states. The reaction mechanism on the triplet surface is more complicated than that on the singlet surface. A total of 19 isomers and 46 transition states have been identified and characterized on the triplet PES. Among them, IM2 (IM2a), IM3 (IM3a, IM3b), and IM10 are the lowest‐lying isomers with thermodynamic stability. Twenty available dissociation channels, depending on the different initial isomers, have been identified. On the singlet surface, only 12 isomers and 16 transition states have been found, and among them IM1(S) and IM2(S) are the lowest‐lying isomers. The higher isomerization and dissociation barriers on the singlet surface indicate that the addition and the subsequent reactions of NH₂+CN are most likely to occur on the triplet PES because of the lower barriers. A prediction can be made for the possible mechanism explaining the production of H+HNCN. Besides HNCN, other major products are NH+HCN and NH+HNC, which are produced by direct dissociation reactions from triplet IM2 and IM3, respectively. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2006     

Bonding and correlation analysis of various Si2CO isomers on the potential energy surface

At various levels of theory, singlet and triplet potential energy surfaces (PESs) of Si₂CO, which has been studied using matrix isolation infrared spectroscopy, are investigated in detail. A total of 30 isomers and 38 interconversion transition states are obtained at the B3LYP/6‐311G(d) level. At the higher CCSD(T)/6‐311+G(2d)//QCISD/6‐311G(2d)+ZPVE level, the global minimum ¹1 (0.0 kcal/mol) corresponds to a three‐membered ring singlet O‐cCSiSi (¹A′). On the singlet PES, the species ¹2 (0.2 kcal/mol) is a bent SiCSiO structure with a ¹A′ electronic state, followed by a three‐membered ring isomer Si‐cCSiO (¹A′) ¹3 (23.1 kcal/mol) and a linear SiCOSi ¹4 (¹Σ⁺) (38.6 kcal/mol). The isomers ¹1, ¹2, ¹3 , and ¹4 possess not only high thermodynamic stabilities, but also high kinetic stabilities. On the triplet PES, two isomers ³1 (³B₂) (18.8 kcal/mol) and ³7 (³A″) (23.3 kcal/mol) also have high thermodynamic and kinetic stabilities. The bonding natures of the relevant species are analyzed. The similarities and differences between C₃O, C₃S, SiC₂O, and SiC₂S are discussed. The present results are also expected to be useful for understanding the initial growing step of the CO‐doped Si vaporization processes. © 2008 Wiley Periodicals, Inc. Int J Quantum Chem, 2009     

Exploring the Reaction Paths on the Potential Energy Surfaces of the S1 and T1 States in Methylenecyclopropane

The reaction paths of methylenecyclopropane 1 on the potential energy surfaces (PESs) of the lowest triplet (T₁) state and the lowest excited singlet (S₁) state, as well as that of the ground state (S₀), were explored by using the nudged elastic band method at the MRMP2//MCSCF/6‐31++G(d,p) and DFT(B3LYP)/6‐31++G(d,p) levels of theory. After vertical excitation of 1, three transition states on the PES of the lowest triplet state and one transition state on the S₁ PES were found along the reaction path to produce a carbene, cyclobutylidene 2. All of these transition states are lower in energy than the S₁ state produced by vertical excitation at the S₀ energy minimum in 1. Fast transition is predicted to occur from the T₁ state or from the S₁ state to the S₀ state due to strong spin‐orbit coupling or nonadiabatic coupling in the geometrical vicinity of 2. On the MRMP2 S₀ PES, the energy barriers of 5.0, 10.3 and 13.5 kcal mol^?1 were obtained for C migration reaction (backward reaction), 1,2‐H migration reaction to cyclobutene 3, and 1,3‐H migration reaction to bicyclopropane 4, respectively, started at 2. The introduction of phenyl groups makes the energy barriers smaller due to the π conjugation between the carbene center and phenyl groups.     

Theoretical study on the gas‐phase reaction mechanism between nickel monoxide and methane for syngas production

The comprehensive mechanism survey on the gas‐phase reaction between nickel monoxide and methane for the formation of syngas, formaldehyde, methanol, water, and methyl radical has been investigated on the triplet and singlet state potential energy surfaces at the B3LYP/6‐311++G(3df, 3pd)//B3LYP/6‐311+G(2d, 2p) levels. The computation reveals that the singlet intermediate HNiOCH₃ is crucial for the syngas formation, whereas two kinds of important reaction intermediates, CH₃NiOH and HNiOCH₃, locate on the deep well, while CH₃NiOH is more energetically favorable than HNiOCH₃ on both the triplet and singlet states. The main products shall be syngas once HNiOCH₃ is created on the singlet state, whereas the main products shall be methyl radical if CH₃NiOH is formed on both singlet and triplet states. For the formation of syngas, the minimal energy reaction pathway (MERP) is more energetically preferable to start on the lowest excited singlet state other than on the ground triplet state. Among the MERP for the formation of syngas, the rate‐determining step (RDS) is the reaction step for the singlet intermediate HNiOCH₃ formation involving an oxidative addition of NiO molecule into the C? H bond of methane, with an energy barrier of 120.3 kJ mol^?1. The syngas formation would be more effective under higher temperature and photolysis reaction condition. © 2009 Wiley Periodicals, Inc. J Comput Chem, 2009     

Halogen switching of azacarbenes C2NH ground states at ab initio and DFT levels

Relative stabilities and singlet–triplet energy differences are calculated for 24 C₂NX azacarbenes (where X is H, F, Cl, and Br). Three skeletal arrangements are employed including azacyclopropenylidene, [(imino)methylene]carbene, and cyanocarbene. Halogens appear to alternate the electronic ground states of C₂NH azacarbenes, from triplet to singlet states, at MP3/6‐311++G**, B1LYP/6‐311++G**, B3LYP/6‐311++G**, MP2/6‐311++G**, MP4(SDTQ)/6‐311++G**, QCISD(T)/6‐311++G**, CCSD(T)/6‐311++G**, CCSD(T)/cc‐pVTZ, G1, and G2 levels of theory. The aromatic characters of singlet cyclic azacyclopropenylidenes are measured using GIAO–NICS calculations. Linear correlations are found between the B3LYP/6‐311++G** calculated LUMO–HOMO energy gaps (ΔE_{HOMO ‐ LUMO}) of the singlet carbenes versus their corresponding singlet–triplet energy separations (ΔE). Electrophilic characters are found for all singlet azacarbenes in their addition reactions to alkenes with the highest electrophilicity being exhibited for X = F. © 2008 Wiley Periodicals, Inc. Heteroatom Chem 19:377–388, 2008; Published online in Wiley InterScience ( www.interscience.wiley.com ). DOI 10.1002/hc.20442