Theoretical investigation on structures,stability and properties of [P,X, Y] (X=C,Si; Y = O,S) isomers

The structures, energetics and stability of the [P, X, Y] (X?=?C, Si; Y?=?O, S) radicals are explored by means of the density functional theory and ab initio levels. Seventeen [P, X, Y] isomers and 14 interconversion transition states are obtained at the B3LYP/6-311G(d) level. At the CCSD(T)/6-311?+?G(2df)//QCISD/6-311G(d)?+?ZPVE level, the lowest-lying isomers are the linear PCO 1a (0.0?kcal/mol), PCS 1b (0.0) and the three-membered ring cPSiO 1c (0.0), cPSiS 1d (0.0) on their respective potential energy surfaces. These four isomers exhibit considerably not only thermodynamic but also kinetic stabilities. Additionally, the cyclic cPCS 2b (32.8) and linear PSiS 2d (18.6) possess also high kinetic stability. All of six isomers 1a, 1b, 2b, 1c, 1d and 2d are considerably stabilized by a barrier of at least 20?kcal/mol, and may be detected in the laboratory or interstellar space. Their valence bond structures and possible formation strategies in the laboratory and space are discussed in detail. Finally, the similarities and discrepancies on structures and stabilities between [P, X, Y] (X?=?C, Si; Y?=?O, S) isomers are compared. These predicted results are highly expected to be informative for the future identification of [P, X, Y] (X?=?C, Si; Y?=?O, S) in the laboratory and space.     

Theoretical study on structures and stability of SiC3S isomers

Various levels of calculations are carried out to explore the potential energy surfaces (PES) of singlet and triplet SiC₃S, a molecule of potential interest in interstellar chemistry. At the DFT/B3LYP/6-311G(d) level, a total of 57 minimum isomers and 92 interconversion transition states are located. 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 the chain-like SiCCCS ³ 1 (0.0?kcal/mol) with a great kinetic stability of 54.1?kcal/mol. In addition, ring isomers CC-cCSSi ¹ 9 (19.8?kcal/mol), S-cCCCSi ¹ 12 (30.4?kcal/mol), S-cCCSiC ¹ 18 (9.4?kcal/mol), S-cSiCCC ¹ 21 (34.4?kcal/mol) and cage-like isomer cage-SiSCCC ¹ 23 (51.8?kcal/mol) also possess considerable kinetic stability (more than 10.0?kcal/mol). As a result, these six isomers are predicted to be possible candidates for future experimental and astrophysical detection. The bond natures and possible formation pathways in interstellar space of the SiCCCS are discussed. The predicted structure and spectroscopic properties for it are expected to be informative for the identification of SiC₃S and even larger SiC _n S species either in laboratory or in space.     

Theoretical study of B3O radical isomers and their interconversion pathways

Despite the fact that B₃O is the second simplest B _n O radical after BO, a controversy recently emerged concerning the molecular structure of its global minimum. Two recent theoretical groups predicted the linear quartet BBBO to be the ground isomer. By contrast, another recent theoretical group reported that B₃O has a doublet B₃-ring ground structure. Moreover, larger B _n O clusters usually have low-lying B₃-ring isomers. In order to determine the accurate energetic competition between linear and cyclic structures in both the doublet and quartet, and to understand the detailed isomerism between various isomers, which is vital for understanding the formation mechanism of B₃O, we report the first potential energy surface (PES) study of B₃O at various computational levels, including CCSD(T)/6-311+G(2df), CCSD(T)/aug-cc-pVTZ, CCSD(T)/aug-cc-pVQZ and G3B3 for the single-point energy, as well as B3LYP/6-311+G(d) and QCISD/6-311+G(d) for geometrical optimisation. It is shown that the isomers in the quartet state are all thermodynamically more stable than the corresponding doublet ones, and on both the quartet and doublet PESs, the linear form has the lowest energy. Therefore, our study on both linear and cyclic isomers shows that the linear quartet BBBO ⁴ 01 is definitively the ground isomer. Although being much less stable than the quartet linear BBBO global minimum by >20 kcal mol^?1, five cyclic isomers exist as local minima, with the bi-cyclic structure ⁴ 02 possessing the smallest barrier of around 15 kcal mol^?1. The dissociation energies for direct combination processes B₃ + O, B₂ + BO and B + B₂O are discussed. The present work may be helpful in obtaining a deep understanding of the doping and oxidation process of pure B _n clusters.     

Structures,vibrational spectra,and relative energetics of FC(O)ONO and FC(O)NO2 isomers at DFT and ab initio levels

The structural and vibrational parameters of FC(O)ONO and FC(O)NO₂ isomers were examined theoretically using the B3LYP/6-311+G(3df) and CCSD(T)/6-311G(d) methods. Four conformers of FC(O)ONO isomer and one FC(O)NO₂ isomer are found here. Among them, the trans–cis and cis–cis FC(O)ONO configuration are new conformers. The energetics were refined with G3//B3LYP and CBS-QB3 calculations. The trans–trans conformer of the FC(O)ONO isomer is found to be the lowest energy structure, with an estimated heat of formation of ?104.9 kcal mol^?1 at 0 K as determined from CBS-QB3 theory. The next lowest structure is the cis–trans FC(O)ONO lying 1.7 kcal mol^?1 above the trans–trans structural form. The highest energy structure is the FC(O)NO₂ isomer with a predicted heat of formation of ?84.8 kcal mol^?1. A comparison of the relative stability of the FCNO₃ isomers with the isomers of ClCNO₃ shows that the Cl analogues follow the same pattern of stability, as do the F isomers. However, the chlorine isomers are unstable relative to their fluorine analogues.     

Isomers of OCS and their reaction with H2O on singlet potential energy surface

The isomers of the carbonyl sulfide (OCS) molecule are investigated in detail at CCSD(T)/cc-pVTZ//MP2/6-311++G(2d,2p) level of theory. One cyclic isomer was identified along with three different linear minima of the OCS molecule. Three interconversion transition states were also located between cyclic and linear forms of OCS. Among these four isomers, the singlet potential energy surface (PES) for the molecule–molecule reaction between the three most energetically favoured isomers of OCS and H₂O has been explored theoretically at the CCSD(T)/cc-pVTZ//MP2/6-311++G(2d,2p) level. This singlet PES comprises of three paths. Path 1 is the reaction of linear OCS molecule with water producing the major product P1 (CO₂?+?H₂S), minor product P2 (S?+?HCOOH) and two isomers via 14 minima and 15 transition states. The Path 2 is an isomerization process in which cyclic isomer of OCS reacts with water molecule via another initial barrierless aduct producing five isomers of the OCS–H₂O system through five interconversion transition states. The reaction of linear COS isomer with water is shown in Path 3. This path produces the radicals SH and COOH from another COS–H₂O complex via a transition state. Among these three products, the product P1 is energetically most favoured. The overall exothermicity of the product channels for the formation of major product P1 on PES is calculated to be about 10.60?kcal/mol possessing initial high entrance barriers of 45.48 and 55.47?kcal/mol in two possible pathways. As the process is favoured thermodynamically but not kinetically, the reaction is expected to be very slow.     

甲基乙烯醚和OH反应机理的理论研究

利用量子化学从头计算的方法对甲基乙烯醚的两个异构体之间的转化,羟基与顺式-甲基乙烯醚和反式-甲基乙烯醚的加成反应,以及羟基提取甲基上的氢原子的反应机理进行了研究.研究结果表明:顺式-甲基乙烯醚比反式-甲基乙烯醚更加稳定,在QCISD/6-31G(d,P)//BHandHLYP/6.311 G(d,P)理论水平下,OH加到顺式-甲基乙烯醚1号住的碳原子上需要跨越的能垒比其它反应通道需要跨越的能垒少7.5～34 KJ/mol,因此是主要的反应通道,而OH加在反式.甲基乙烯醚2号位的碳原子上所需要跨越的能垒比其它反应路径所需要跨越的能垒少8.3～26.7 kJ/mol,因此是主要的反应路径.利用经典过渡态理论计算了总的速率常数     

Theoretical Study on Reaction Mechanism of Aluminum-Water System

在B3LYP/6-311+G(d,p)和QCISD(T)/6-311++G(d,p) 水平上,对气相中铝与水的反应进行了计算. 理论计算结果表明,铝水反应经过了三条可能的反应通道并且涉及了四个异构体、七个过渡态和两种产物. 其中两条反应通道所涉及的中间物和产物已被实验测得. 同时,在298 K和2000 K计算了铝水反应的焓变和吉布斯自由能,一些计算结果与以前的计算和实验结果吻合较好.     

过氧硝酸分解势能面的理论研究

使用密度泛函理论B3LYP/6-311+ G(2d,2p)研究了过氧硝酸的最低能量结构．采用耦合簇方法CCSD(T)/aug-cc-pVDZ首次分别扫描了过氧硝酸沿氧-氮和氧-氧键的分解势能面．计算结果表明在氧-氮势能面上,当O3—N4键长是2.82 ?时,对应的疏松过渡态的能垒是25.6 kcal/mol;在氧$-$氧键的势能面上,当O2—O3键长是2.35 ?时,对应的疏松过渡态的能垒是37.4 kcal/mol．这表明过氧硝酸更容易分解为HO₂和NO₂．     

Li原子修饰笼型Si6团簇的结构和储氢性能

利用杂化密度泛函B3LYP方法, 在6-311+G(d, p)基组水平上对Si6和Li修饰的Si6团簇的几何结构和电子性质及储氢性能进行模拟计算和理论研究. 结果表明, Si6团簇最低能量构型为笼型结构, 纯Si6团簇不能有效吸附氢分子. Li原子的引入显著改善了Si6团簇的储氢能力. 以两个Li原子端位修饰Si6团簇为载体, 其氢分子的平均吸附能为1.692~2.755 kcal/mol, 每个Li原子周围可以有效吸附五个氢分子, 储氢密度可达9.952wt%. 合适的吸附能和较高储氢密度表明Li修饰Si6团簇有望成为理想的储氢材料.     

Li原子修饰笼型Si_6团簇的结构和储氢性能

利用杂化密度泛函B3LYP方法,在6-311+G(d,p)基组水平上对Si_6和Li修饰的Si_6团簇的几何结构和电子性质及储氢性能进行模拟计算和理论研究.结果表明,Si_6团簇最低能量构型为笼型结构,纯Si_6团簇不能有效吸附氢分子.Li原子的引入显著改善了Si_6团簇的储氢能力.以两个Li原子端位修饰Si_6团簇为载体,其氢分子的平均吸附能为1.692~2.755 kcal/mol,每个Li原子周围可以有效吸附五个氢分子,储氢密度可达9.952 wt%.合适的吸附能和较高储氢密度表明Li修饰Si_6团簇有望成为理想的储氢材料.     

Conformational analysis of 4,4‐dimethyl‐4‐silathiane and its S‐oxides

4,4‐Dimethyl‐4‐silathiane and its S‐oxides [n = 0 ( 1 ), 1 ( 2 ), 2 ( 3 )] were studied experimentally by variable temperature dynamic NMR spectroscopy down to 103 K and the frozen ring inversion was revealed for all three compounds. The barriers for the degenerate ring inversion in 1 and 3 were measured to be 4.8 and 5.0 kcal/mol at the coalescence temperatures of 111 and 116 K, respectively, and practically coincide with the calculated barriers of 4.60 kcal/mol in 1 and 4.46 kcal/mol in 3 . The frozen equilibrium mixture 2‐ax/2‐eq contains 37% of the 2‐ax and 63% of the 2‐eq conformer. The ring inversion barrier proved to be ca. 4.8 kcal/mol. Calculations at the B3LYP/6‐311+G(d,p) level of theory showed the 2‐ax conformer to be 0.90 kcal/mol more stable than the 2‐eq conformer in the gas phase whereas in solution the relative stability of the conformers calculated using the PCM model at the same level of theory is inverted to become 0.19 (in CHCl₃) or 0.36 kcal/mol (in DMSO) in favor of the 2‐eq conformer. The chair–chair interconversion mechanism of sulfoxide 2 includes two intermediate energetically equivalent 1,4‐twist forms and the 2,5‐boat transition state: 2‐ax (chair) ? 2 (1,4‐twist) ? [ 2 (2,5‐boat)]^≠ ? 2 (1,4‐twist) ? 2‐eq (chair). The calculated ring inversion barriers are 5.1 ( 2‐ax → 2‐eq ) and 4.2 kcal/mol ( 2‐eq → 2‐ax ) in the gas phase, and 4.03 and 4.22 kcal/mol, respectively, in chloroform. Copyright © 2011 John Wiley & Sons, Ltd.     

Si3Xn (X=C,O,N;n=1,2)团簇的密度泛函研究

使用密度泛函理论(DFT)的杂化密度泛函B3LYP方法在6-31G*基组水平上对Si3Xn(X=C,O,N;n=1,2)团簇各种可能的构型进行几何结构优化,预测了各团簇的最稳定结构.并对最稳定结构的电子结构、振动特性、成键特性和电荷特性等进行了理论研究.结果表明团簇的几何结构都是平面结构,通常Si3X2出现是Si-X键,较少出现X-X键;而Si3X中出现Si-X键和Si-Si键共存,Si3Xn(X=C,O,N;n=1,2)团簇的电荷布局分布表明这种电荷转移的作用使得团簇中所有X原子呈负电性,Si原子显正电性.处于不同位置的Si原子呈不同大小布局数,而且由于Si3X2的对称性,2个X负电性相同.     

DFT study on the mechanism of trimolecular radical reactions: isomerisation in small clusters

ABSTRACT

Structural and thermodynamic properties of 48 trimolecular clusters containing one radicl and two protic molecules (H₂O, NH₃, H₂O₂, CH₃OH, HOCl) were studied at B3LYP/6-311++G(3df,3pd) level of theory. These radical-clusters have non-cyclic structures and are stabilised via two inter-molecular hydrogen bonding interactions. The calculated enthalpies of formation of the radical-clusters were generally in the range of ?30 to ?50 kJ/mol. The calculated activation energies (E_a) of the intra-cluster hydrogen transfers were smaller than 70 kJ/mol. Also, structures and thermodynamics of 15 cyclic molecular clusters as well as multi-hydrogen transfers in them were investigated. The results showed that the stability of the cyclic clusters and activation energies of the multi-hydrogen transfers depend on the cluster size.     

Structure and stability of neutral and cationic AlnO clusters

In this paper various structural possibilities for AlnO neutral and cationic isomers were investigated by using the B3LYP/6-311G（3df） method. Calculations of this paper predicted the existence of a number of previously unknown isomers. The stabilities of the AlnO （n = 2 - 7） clusters with even n are greater than those with odd n, however the stabilities Of cationic ions have the opposite odd-even alternation. The mass spectra observations of Al17O^＋ and Al19O^＋ ions support our theoretical predictions on their stabilities.     

Isomerization of C3H3NO isomers: ab initio study

The structures and isomerization process of C₃H₃NO species have been explored at the MP2/6–311++G(d,p) level of theory of the ab initio method. Eleven minima and four interconversion transition states have been identified. The zero-point vibrational energy corrections were made to predict reliable energies. We predict a five-membered ring-like structure to be the lowest energy isomer, which is 177.73?kcal?mol^?1 more stable than the least stable isomer X found on the potential energy surface. The transition states and minima isomers were verified by frequency calculation. Intrinsic reaction coordinate (IRC) calculations have been performed to confirm that each transition state is linked by the desired reactants and products. The isomer stabilities have been studied using the relative energies, chemical hardness and chemical potential. The MHP principle could not predict the order of stability for C₃H₃NO isomers as arrived at with the relative energies. The role of intramolecular hydrogen bonds on the equilibrium structure has been discussed. The energy barrier and reaction enthalpy have been calculated during isomerization.     

Computational study on structures,thermochemical properties,and bond energies of disulfide oxygen (S–S–O)‐bridged CH3SSOH and CH3SS(=O)H and radicals

Cleavage of disulfide bonds is a common method used in linking peptides to proteins in biochemical reactions. The structures, internal rotor potentials, bond energies, and thermochemical properties (Δ_fH°, S°, and Cp(T)) of the S–S bridge molecules CH₃SSOH and CH₃SS(=O)H and the radicals CH₃SS?=O and C?H₂SSOH that correspond to H‐atom loss are determined by computational chemistry. Structure and thermochemical parameters (S° and Cp(T)) are determined using density functional Becke, three‐parameter, Lee–Yang–Parr (B3LYP)/6‐31++G (d, p), B3LYP/6‐311++G (3df, 2p). The enthalpies of formation for stable species are calculated using the total energies at B3LYP/6‐31++G (d, p), B3LYP/6‐311++G (3df, 2p), and the higher level composite CBS–QB3 levels with work reactions that are close to isodesmic in most cases. The enthalpies of formation for CH₃SSOH, CH₃SS(=O)H are ?38.3 and ?16.6 kcal mol^?1, respectively, where the difference is in enthalpy RSO–H versus RS(=O)–H bonding. The C–H bond energy of CH₃SSOH is 99.2 kcal mol^?1, and the O–H bond energy is weaker at 76.9 kcal mol^?1. Cleavage of the weak O–H bond in CH₃SSOH results in an electron rearrangement upon loss of the CH₃SSO–H hydrogen atom; the radical rearranges to form the more stable CH₃SS· = O radical structure. Cleavage of the C–H bond in CH₃SS(=O)H results in an unstable [CH₂SS(=O)H]* intermediate, which decomposes exothermically to lower energy CH₂ = S + HSO. The CH₃SS(=O)–H bond energy is quite weak at 54.8 kcal mol^?1 with the H–C bond estimated at between 91 and 98 kcal mol^?1. Disulfide bond energies for CH₃S–SOH and CH3S–S(=O)H are low: 67.1 and 39.2 kcal mol^?1. Copyright © 2011 John Wiley & Sons, Ltd.     

气相中自由基OH在不同催化剂下对甲酰氟FCHO进行提氢反应的理论性研究

分别在水、甲酸和硫酸存在的情况下,通过CCSD(T)//M06-2X/6-311++G(3df,3pd)的理论方法,对大气中自由基OH提取甲酰氟FCHO上的氢进行了反应机理和动力学的研究.计算结果表明相对于反应物,加入催化剂的过渡态的能垒从3.64 kcal/mol分别下降到-2.89、-6.25和-7.76 kcal/mol,表明水、甲酸和硫酸在甲酰氟FCHO和自由基OH提氢反应中起了重要作用.通过运用具有Eckart隧道校正的传统过渡态理论计算出的动力学数据表明通道X...FCHO+OH(X=H₂O, HCOOH,或者H₂SO₄)要比通道X...OH+FCHO更有利于反应的发生.催化剂水、甲酸和硫酸的加入使甲酰氟FCHO 和自由基OH提氢反应的速率常数要比不加催化剂时小,说明了催化剂的加入不能促进大气中甲酰氟FCHO和自由基OH的反应.     

Reaction of C2HCl2+O2: Combined TR-FTIR Spectroscopy and Electronic Structure

用时间分辨傅立叶变换红外发射光谱(TR-FTIR)和G3MP2//B3LYP/6-311G(d,p)水平的电子结构计算研究了环境化学中重要的二氯代乙烯自由基C₂HCl₂和O₂分子的基元反应通道和机理. 通过0.5 cm^-1高分辨的TR-FTIR发射光谱观察到三种振动激发态产物CO₂、CO和HCl,由光谱拟合得到CO和HCl的振动态分布,结合电子结构计算的反应势能曲线,提出反应机理和能量上最可能的反     

小型二氧化钛团簇(TiO2)n (n=2-8)的低能异构体研究

钛和氧之间存在多种成键方式,但迄今为止,二氧化钛团簇均只有少数几种异构体被报道. 与广泛使用的全局优化方法不同,本工作通过优化大量的随机初始结构,获得(TiO₂)_n (n=2-8)团簇稳定异构体. 首先利用PM6半经验方法对高达一万个以上的初始结构进行初步的优化筛选,并对筛选出的结构进行进一步的DFT计算以获得二氧化钛团簇的稳定异构体. 利用这种策略,发现了大量未经报道的稳定异构体,并提出了(TiO₂)₅和(TiO₂)₈新的最稳定异构体. 这些结构中包括含3个末端氧原子的异构体、含5配位氧原子和6配位钛原子的异构体等未经报道的新结构类型. 与丰富成键特征相对应,发现异构体数目随团簇尺寸的增大而急剧增加,对于(TiO₂)₇和(TiO₂)₈,能量在30 kcal/mol以内的异构体都在50个以上. 该工作发现了大量的二氧化钛小型团簇异构体,并凸显了其多样的结构特征,增进了对二氧化钛纳米团簇的结构、成键的理解,并为进一步的理论模拟、力场优化等提供了理论基础.     

GaH(D,T)分子基态结构与势能函数

采用量子力学从头计算方法,运用单双取代二次组态相互作用和单双（三）取代二次组态相互作用并考虑基组6-311++G（3df,3pd）计算优化了GaH（D,T）分子基态X¹Σ⁺的结构和离解能.并采用最小二乘法拟合改进的Murrell-Sorbie函数得到了相应的势能函数.计算得到的光谱常数与实验光谱数据符合得很好. 关键词： GaH（D T）分子基态 势能函数 Murrell-Sorbie函数