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

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     

Theoretical Studies on Reaction Mechanism of CH3SO with NO

The potential energy surface （PES） of CH3SO radical with NO reaction has been studied at MP2/6-311G（2df, p） and QCISD/6-311G（2df, p） levels. Geometries of the reactants, transition states （TS） and products were optimized at B3LYP/6-311G （d,p） level. The geometries of the transition states were found for the first time. The calculated results show that the reaction can proceed via singlet-state or triplet-state PES. Because of the high energy barrier of triplet surface, the singlet surface reactions are dominant. The topological analysis of electron density shows that there are two kinds of structaral transition states （the bifurcation-type ring structure transition state and the T-shaped conflict structure transition state） in the titled reaction. The total electronic density of the reactants, TS and products and the spin electronic density on the triplet surface were also discussed in this paper.     

A theoretical study on C2HXSi silylenes (X = H,CN, NH2, and OMe)

Relative stabilities and structural characters of 30 silylenic C₂HXSi species (X = H, NH₂, CN, and OMe), with singlet (s) and/or triplet (t) states, are calculated at six levels of theory: HF/6‐311++G**, MP3/6‐31G*, B1LYP/6‐311++G**, B3LYP/6‐311++G**, MP2/6‐311++G**, and MP4(SDTQ)/6‐311++G**. The four possible isomers considered for C₂SiHX are (i) 3‐X‐1‐silacyclopropenylidene ( 1 _s‐X and 1 _t‐X ), (ii) X‐vinilydensilylene ( 2 _s‐X and 2 _t‐X ), (iii) ethynyl‐X‐silylene ( 3 _s‐X and 3 _t‐X ), and (iv) (X‐ethynyl)silylene ( 4 _s‐X and 4 _t‐X ). The GIAO–NICS calculations show that singlet cyclic structures, 1 _s‐X , are considerably more aromatic than benzene. Conversely, triplet cyclic C₂HCNSi breaks down through optimization, and transforms into a novel high‐spin acyclic carbenosilylene minimum ( 1 _t‐CN ). Singlet 3 and triplet 3 cross at a divalent angle (|XSiC) of 152°. This angle narrows to 137° for crossing of singlet 3 _s‐CN and triplet, 3 _t‐CN . The smallest |XSiC occurs at 132° for crossing of 3 _s‐H and 3 _t‐H . © 2007 Wiley Periodicals, Inc. Heteroatom Chem 18:283–293, 2007; Published online in Wiley InterScience ( www.interscience.wiley.com ). DOI 10.1002/hc.20297     

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     

Theoretical study on structures and stability of triplet SiC<Subscript>3</Subscript>O isomers

Various levels of calculations are carried out~for exploring the potential energy surface (PES) of triplet SiC₃O, a molecule of potential interest in interstellar chemistry. A total of 38 isomers are located on the PES including chain-like, cyclic and cage-like structures, which are connected by 87 interconversion transition states at the DFT/B3LYP/6-311G(d) level. The structures of the most relevant isomers and transition states are further optimized at the QCISD/6-311G(d) level followed by CCSD(T)/6-311+G(2df) single-point energy calculations. At the QCISD level, the lowest lying isomer is a linear SiCCCO 1 (0.0 kcal/mol) with the ³ ∑ electronic state, which possesses great kinetic stability of 59.5 kcal/mol and predominant resonant structure . In addition, the bent isomers CSiCCO 2 (68.3 kcal/mol) and OSiCCC 5 (60.1 kcal/mol) with considerable kinetic stability are also predicted to be candidates for future experimental and astrophysical detection. The bond natures and possible formation pathways in interstellar space of the three stable isomers are discussed. The predicted structures and spectroscopic properties for the relevant isomers are expected to be informative for the identification of SiC₃O and even larger SiC _n O species in laboratory and interstellar medium. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

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     

Isomerization study of C5H5NO molecules

The complex potential energy surface (PES) for the isomerization of C₅H₅NO species, including 18 isomers and 23 interconversion transition states, is probed theoretically at the B3LYP/6‐311++G(d,p) and MP2//B3LYP/6‐311++G(d,p) levels of theory. The geometries and relative energies for various stationary points were determined. The zero‐point vibrational energy (ZPVE) corrections have been made to calculate the reliable energy. We predicted a six‐membered ring structure as a global minima isomer I, which is 118.49 and 131.48 kcal · mol^?1 more stable than the least stable, four‐ and three‐membered ring isomer VIII at B3LYP and MP2//B3LYP levels of theory, respectively. The isomers and interconversion transition states have verified by frequency calculation. The intrinsic reaction coordinates (IRC) calculations have been performed to confirm that each transition state is linked by the desired reactants and products. The isomer stability has been studied using relative energies, chemical hardness, and chemical potential. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

Computational mechanistic study of methanol and molecular oxygen reaction on the triplet and singlet potential energy surfaces

The reaction of CH₃OH with the O₂ on the triplet and singlet potential energy surfaces (PES) was carried out using the B3LYP, MP2, and CCSD(T)//B3LYP theoretical approaches in connection with the 6-311++G(3df–3pd) basis set. Three pre-reactive complexes, ¹C1, ¹C2, and ³C1, on the singlet and triplet PES were formed between methanol and molecular oxygen. From a variety of the complexes, seven types of products are obtained, of which four types are found to be thermodynamically stable. Results reveal that there exists one intersystem crossing between triplet and singlet PES. For P4 adduct that is the main and kinetically the most favorable product, the rate constants are calculated in the temperature range of 200–1,000 K in the reliable pathway.     

DFT investigation of the mechanism of CH2CO + O(3P) reaction

A detailed theoretical survey of the potential energy surface (PES) for the CH₂CO + O(³P) reaction is carried out at the QCISD(T)/6‐311+G(3df,2p)//B3LYP/6‐311+G(d,p) level. The geometries, vibrational frequencies, and energies of all stationary points involved in the reaction are calculated at the B3LYP/6‐311+G(d,p) level. More accurate energy information is provided by single‐point calculations at the QCISD(T)/6‐311+G(3df,2p) level. Relationships of the reactants, transition states, intermediates, and products are confirmed by the intrinsic reaction coordinate (IRC) calculations. The results suggest that P1(CH₂+CO₂) is the most important product. This study presents highlights of the mechanism of the title reaction. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2005     

Stabilization Studies on Divalent Five‐Membered Rings C4H4C,C4H4Si,C4H4Ge,C4H4Sn,and C4H4Pb

Energy differences, ΔX_S‐t (X = E, H and G) (ΔX_S‐t = X_(singlet)‐X_(triplet)) between singlet (s) and triplet (t) states are calculated at B3LYP/6‐311++G (3df,2p). The DFT calculations show that the triplet state of C₄H₄C is a ground state with planar conformer respect to its corresponding nonplanar singlet state. Both singlet and triplet states of C₄H₄M (M = Si, Ge, Sn and Pb) have a planar conformer with the singlet ground state. Four isodesmic reactions are presented for determining the stability energies, SE. NICS calculations are carried out for C₄H₄M to determine the aromatic character.     

Ab initio study on the reaction mechanism of ozone with bromine atom

Ab initio calculations of the potential energy surface (PES) for the Br+O₃ reaction have been performed using the MP2, CCSD(T), and QCISD(T) methods with 6‐31G(d), 6‐311G(d), and 6‐311+G(3df). The reaction begins with a transition state (TS) when the Br atom attacks a terminal oxygen of ozone, producing an intermediate, the bromine trioxide (M), which immediately dissociates to BrO+O₂. The geometry optimizations of the reactants, products, and intermediate and transition states are carried out at the MP2/6‐31G(d) level. The reaction potential barrier is 3.09 kcal/mol at the CCSD(T)/6‐311+G(3df)//MP2 level, which shows that the bromine atom trends intensively to react with the ozone. The comparison of the Br+O₃ reaction with the F+O₃ and Cl+O₃ reactions indicates that the reactions of ozone with the halogen atoms have the similar reaction mechanism. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

Ti, Na与O2反应机理的理论研究

基于MP2/6-311+G(d)水平, 分别对过渡金属Ti和碱金属Na与O2的反应机理进行了研究. 比较了Ti和Na分别以垂直O—O键和沿着O—O键的方向逼近O2, 以及中性体系Ti/Na+O2和带1个负电荷的体系(Ti/Na+O2)-的情况. 详细分析了不同反应路径的结合能和电荷变化的曲线, 预测了最佳反应方式. 结果表明, 垂直接近方式要比水平接近方式更具有优势; 体系带一个负电荷(Ti/Na+O2)-有利于金属与O2的结合. 同时, 计算结果表明在Ti+O2和(Ti+O2)-体系中Ti容易与单态的O2结合; 在中性体系中Na也容易与单态O2结合, 而在(Na+O2)-体系中Na与三态O2的结合更稳定. 在CCSD(T)/6-311++G(3df)//MP2/6-311+G(d)水平下, 计算了Ti+O2和(Ti+O2)-的反应势能面.     

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     

Structure and Stability of Isomers of the Promising Interstellar Molecule PC3O

DFT/B3LYP/6-311G(d) and CCSD(T)/6-311G(2d) single-point calculations are carried out for exploring the doublet potential energy surface (PES) of PC₃O, a molecule of potential interest in interstellar chemistry. A total of 29 minima connected by 65 interconversion transition states are located. The structures of the most relevant isomers and transition states are further optimized at the QCISD level followed by CCSD(T) single-point energy calculations. At the CCSD(T)/6-311G(2df)//QCISD/6-311G(d)+ZPVE level, the global minimum is the quasi-linear structure PCCCO 1 (0.0 kcal/mol) with a great kinetic stability of 47.9 kcal/mol, and the cumulenic form features largely in its resonance structures. Moreover, the chainlike isomer OPCCC 3 (64.5) and five-membered-ring species cPCCCO 19 (77.8) possess considerable kinetic stability of about 18.0 kcal/mol. All these three isomers are very promising candidates for future experimental and astrophysical detection. Additionally, a three-membered-ring isomer CC-cCOP 10 (69.6) has slightly lower kinetic stability of around 15 kcal/mol and may also be experimentally observable. Possible formation mechanisms of the four stable isomers in interstellar space are discussed. The present research is the first attempt to study the isomerization and dissociation mechanisms of PC _n O series. The predicted spectroscopic properties, including harmonic vibrational frequencies, dipole moments and rotational constants for the relevant isomers, are expected to be informative for the identification of PC₃O in laboratory and interstellar medium.     

Structure and Stability of Interstellar Molecule C_3S

ＩｎｔｒｏｄｕｃｔｉｏｎＳｍａｌｌｃｌｕｓｔｅｒｓｃｏｎｔａｉｎｉｎｇｃａｒｂｏｎａｎｄｓｕｌｆｕｒ,ｓｕｃｈａｓＣＳ ,Ｃ2 ＳａｎｄＣ3 Ｓ ,ｗｈｉｃｈｐｏｓｓｅｓｓｌａｒｇｅｐｅｒｍａｎｅｎｔｄｉｐｏｌｅｍｏｍｅｎｔｓａｎｄｈａｖｅｂｅｅｎｉｄｅｎｔｉｆｉｅｄｉｎｔｈｅｃａｒｂｏｎｓｔａｒＩＲＣ+ 10°2 16ａｎｄｉｎｔｈｅＴａｕｒｕｓｃｏｌｄｍｏｌｅｃｕｌａｒｄｅｎｓｅｃｌｏｕｄＴＭＣ 1,1 7ｈａｖｅａｔｔｒａｃｔｅｄｍｕｃｈａｔｔｅｎｔｉｏｎｂｅｃａｕｓｅｏｆｔｈｅｉｒｉｍｐｏｒｔａｎｔｒｏｌｅｓ…     

A density functional theory study on the atmospheric reaction of CH3O2 with HS: Mechanism and kinetics

The reaction mechanism of CH₃O₂ and HS was systematically investigated by density functional theory (DFT). Six singlet pathways and seven triplet ones are located on the potential surface (PES). The result indicates that the main products are CH₃O and HSO both on the singlet and triplet PES, different from the CH₃O₂ + OH reaction. Moreover, deformation density (ρ_def) and atoms in molecules (AIM) analyses were carried out to further uncover the nature of chemical bonding evolution in the primary pathways. Furthermore, reaction rate constants were calculated in the temperature range from 200 to 1000 K using the transition state theory with the Wigner and Eckart tunneling corrections. Our results can shed light on the title reaction and offer instructions for analogous atmospheric reactions, as well as experimental research in the future.     

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.     

Structures and stabilities of HPO<Subscript>2</Subscript> isomers

The potential energy surface of HPO₂ system including eight isomers and twelve transition states is predicated at MP2/6-311++G(d, p) and QCISD(t)/6-311++G(3df,2p)(single-point) levels of theory. On the potential energy surface, cis-HOPO(E1) is found to be thermodynamically and kinetically most stable isomer followed by trans-HOPO(E2) and HPO(O)(C_2v, E3) at 10.99 and 48.36 kJ/mol higher, respectively. Based on the potential energy surface, only E1 and E3 are thermodynamically stable isomers, and should be experimentally observable. The products cis-HPOO(E5) and trans-HPOO(E6) in the first-step reaction of HP with O₂ can isomerize into isomer E1 that has higher stability. The reaction of OH with PO will directly lead to the formation of isomer E1. The computed results are well consistent with the previous experimental studies.