SiCS分子结构及其稳定性的理论研究

采用DFT,QCISD及CCSD(T)方法对单重态、三重态SiCS的分子体系势能面进行理论计算,在QCISD/6-311G(d)水平上得到由3个过渡态连接的5个稳定构型.经动力学及热力学分析均是稳定的三重态线型分子SiCS(³1)、单重态线型分子SiCS(¹1)以及单重态的环状分子cSiCS(¹2).     

预测[C,O,S]体系的稳定异构体

采用 DFT, QCISD及CCSD(T)方法对单重态和三重态的COS分子体系势能面进行了理论计算, 在 QCISD/6-311G(d) 水平上得到4个过渡态连接的6个稳定构型. 经动力学及热力学分析发现只有一个单重态线性的分子11(O—C—S) 能够稳定存在.     

NC2S+离子的结构和稳定性的理论研究

在密度泛函和从头算理论水平下计算了单重态的NC2S+离子的结构、能量、光谱以及稳定性. 在B3LYP/6-311G(d)水平下, 得到8个异构体, 它们由15个过渡态相连接. 在CCSD(T)/6-311+G(2df)//QCISD/6-311G(d)+ZPVE水平下, 得到能量最低的异构体是直线型的具有1Σ电子态的NCCS+(1)(0.0 kJ/mol), 其次是直线型的异构体CNCS+(2)(54.8 kJ/mol). 两个低能量的异构体1和2及另外一个高能量的直线型异构体CCNS+(3)(323.8 kJ/mol)都具有相当大的动力学稳定性, 这三个异构体在具备一定条件的实验室和星际条件下是可以进行观测的. 分析了这3个异构体的成键性质.     

三重态SiCP2异构体的结构和稳定性的理论研究

采用DFT,QCISD和CCSD(T)等理论计算方法对三重态SiCP2异构体的结构和稳定性进行了理论研究.在B3LYP/6-311G(d)水平下,共计算得到由17个过渡态相连接的15个异构体.在CCSD(T)/6-311 +G(2df)//QCISD/6-311G(d)水平下,考虑重点振动能相对能量最低的三元环状异构体P-cCSiP 8(0.0 kJ/mol)及四元环状结构的cPCSiP 4具有相当大的动力学稳定性,在一定的实验室和星际条件下可能被检测到.另外,对它们的成键性质也进行了分析.     

密度泛函理论研究SiC与C2H4反应的产物通道

在B3LYP/6-311++G(d,p)水平上研究了SiC与乙烯的单重态和三重态反应机理,优化得到了反应物、过渡态、中间体和产物的几何构型;通过振动分析对过渡态和中间体构型进行了确认。在CCSD(T)/cc-pVTZ水平上对计算得到的构型进行了能量校正。计算结果表明,SiC+C2H4反应在单重态和三重态条件下均可发生,其中单重态反应为主反应通道,1P5为主产物。     

HPS2异构体结构与相对稳定性

在MP2/6-311++G(d,p)和QCISD(t)/6-311++G(3df,2p)(单点)水平下计算得到了包括9个异构体和15个过渡态的HPS2势能面.其中,异构体trans-HSPS(E1)的能量最低, 其次是cis-HSPS(E2)和HPS(S)(C2_V,E3),能量分别比trans-HSPS(E1)高3.43和14.17 kJ/mol.计算结果表明,异构体E1,E2,E3和立体的三元环结构的异构体HP(S)S(C_S,E4), 及与E4共存的trans-HPSS(E5)和cis-HPSS(E6)具有一定的稳定性.在QCISD(t)/6-311++G(3df,2p)//MP2/6-311++G(d,p)并包含零点能水平下,PH和S₂反应生成的E6和E5分别越过65.75和71.73 kJ/mol的势垒就可以异构化为具有较高动力学稳定性的产物E4, 计算结果对实验认定的产物是cis-HPSS(E6)的结论进行了修正.     

SiC(X~3Ⅱ)的结构与势能函数的从头计算

应用密度泛函理论B3LYP和B3P86以及组态相互作用方法 CCSD,CCSD(T),QCISD和QCISD(T),采用6-311g,6-311G(df),6-311+G(d,p),6-311++G(3df,3pd),aug-cc-pvdz和D95(d)多种基组,优化计算了SiC分子的平衡结构和能量.通过优化计算结果和实验数据R=0.171 82nm进行对比,选择B3LYP/6-311G(df),CCSD/6-311G(df)和QCISD/6-311G(df)方法对SiC(X3Ⅱ)分子进行单点能扫描,同时计算其光谱参数(Be,αe,ωe,ωexe)和力常数(f2,f3,f4),这些计算结果与实验数值相吻合,为研究SiC/SiC复合材料提供了理论数据参考.     

乙硼烷离子和自由基异构体重排与体系势能面的能量子化学计算研究

在B3LYP/6-311G(d,p)和QCISD(T)/6-311++G(3df,2p)(单点)水平下, 对B2H+5阳离子和B2H*5自由基全优化得到9个几何异构体: B2H+5单态体系(D3h, C1), B2H+5三重态(Cs, Cs, C1), B2H*5自由基(C2v, Cs, Cs, Cs). 得到势能面上与体系异构化过程相联系的5种过渡态.     

[C, N, N, P]体系异构体结构与稳定性

在B3LYP/6-311G(d)和QCISD(t)/6-311+G(2df)水平下计算得到了[C, N, N, P]体系的9个异构体、20个过渡态及一些相关的解离碎片, 并研究了异构体之间的异构化过程, 讨论了该体系异构体的结构与稳定性. 结果表明在得到的9个异构体中, 具有2A′电子态的、非线型链状NCNP, NCPN, CNPN和CNNP 4个异构体是较稳定的. 其余异构体由于易于重排或是解离, 稳定性较低. 研究同时给出了这4个异构体的振动频率、转动常数、偶极矩、第一绝热电离能和绝热电子亲合能等相关数据, 为可能的实验和星际光谱观测提供相关的理论数据. 得到的计算结果与等价电子体系SiNNP进行了比较.     

三重态类硅烯HB＝SiLiF的结构及其与R—H(R=F,OH,NH_2)的插入反应

用密度泛函理论(DFT)和二次组态相互作用(QCISD)方法研究了三重态类硅烯HB=SiLiF的结构及其与RH(R=F,OH,NH2)的插入反应.计算结果表明,类硅烯HB=SiLiF有三种平衡构型,其中四元环构型能量最低,是其存在的主要构型.HB=SiLiF与HF,H2O和NH3发生插入反应的机理相同.QCISD/6-311++G(d,p)//B3LYP/6-311+G(d,p)计算的三个反应的势垒分别为124.85,140.67和148.16kJ·mol-1,反应热分别为-2.22,20.08和23.22kJ·mol-1.相同条件下发生插入反应时,反应活性都是H—FH—OHH—NH2.     

Fourier transform infrared isotopic study of ZnC3: identification of the ν1(a') mode

An isomer of ZnC(3) with bent structure in its (1)A(') electronic state has been detected in the products from the dual laser ablation of carbon and zinc rods that were trapped in solid Ar at ~12 K. Measurements of (13)C isotopic shifts have enabled the identification of the ν(1)(a(')) asymmetric carbon stretching fundamental at 1858.9 cm(-1). The experimental results are in good agreement with the predictions of DFT-B3LYP/6-311G(d) calculations that indicate a singlet bent isomer ground state structure with triplet linear and singlet cyclic isomers lying slightly higher in energy. This is the first optical detection of any isomer of ZnC(3).     

Structure and Stability of Interstellar Molecule C_3S

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

Lowest singlet and triplet potential energy surfaces of S2N2

Forty four stationary points have been located on the lowest singlet and triplet potential energy surfaces of S(2)N(2). Ten minima and ten saddle points on the lowest singlet surface and eleven minima and thirteen saddle points on the lowest triplet surface were found. All saddle points were connected to minima or lower-order saddle points by following the intrinsic reaction coordinate. Renner-Teller effects in the linear isomers were studied by examining their bending curves. The S(2)N(2) polymerization mechanism was investigated by first locating the transition state corresponding to ring opening and then considering all species connected to it that are close in energy. The commonly accepted mechanism is problematic due to the number of species that would lead to dissociation to SN + SN. Other possible isomers that are consistent with the experimental evidence but do not connect to SN radicals in the dissociation limit were examined. A mechanism of polymerization to (SN)(x)() is proposed that involves excitation of the square planar singlet molecule to the triplet surface. The triplet species then undergoes a puckering, and polymerization occurs in a direction approximately perpendicular to the S(2)N(2) plane. Consideration of the predicted vibrational frequencies suggests the structure of the second isomer of S(2)N(2). This isomer has a trans-NSSN structure with a long SS bond. The energetics of trans-NSSN are consistent with the observed temperature effects in the dimerization of SN. Analysis of the bending curves of linear NSSN and NSNS indicates that trans-NSSN is the only isomer which has a small yet significant barrier to that dimerization.     

Singlet and triplet potential energy surfaces of C3H2

The detailed singlet and triplet potential energy surfaces of C₃H₂ involving nine isomers and 13 transition structures are studied at the G3 level of theory. The rearrangement mechanisms and the electronic properties of various isomers in a broad energy range have been studied in both singlet and triplet states. Cyclopropenylidene and propargylene are found to be the most stable isomers in the singlet and triplet states, respectively. The singlet isomers are found to be more kinetically stable species as a result of high conversion barriers through which they pass. The calculations indicate that cyclopropyne in its triplet state is the least kinetically stable isomer. It is realized that the G3 method comprises both computational cost and accuracy and thus can be applied to investigation of potential energy surface of small molecules.     

A coupled-cluster study of linear and rhombic boron nitride dimers: a revisit

The potential energy surfaces of both singlet and triplet B₂N₂ have been investigated computationally at the coupled-cluster level with a polarized triple zeta basis set augmented with diffuse functions. Calculated vibrational frequencies and intensities are also reported. The triplet species are consistently more stable than their singlet analogs and the stabilities of the linear B₂N₂ isomers increase with increasing number of B–N bonds. The most stable isomer is the linear triplet BNBN isomer with a rhombic form with a short diagonal BB distance close in energy. Our results are consistent with the results of the matrix IR studies of Andrews et al. nucleus-independent chemical shift (NICS) values were calculated for the singlet D_2h rhombic form and its C_2v dication, and these were compared to those of the D_2h cyclobutadiene and its D_2d dication, respectively. Electron density plots for the linear and rhombic B₂N₂ minima showed similar distributions for the singlet and triplet states. These plots confirmed weak BB bonding interactions in both rhombic forms but larger BN bond orders.     

Two-state reactivity, electromerism, tautomerism, and "surprise" isomers in the formation of compound II of the enzyme horseradish peroxidase from the principal species, compound I

QM and QM/MM calculations on Compound II, the enigmatic species in the catalytic cycle of the horseradish peroxidase enzyme, reveal six low-lying isomers. The principal isomer is the triplet oxo-ferryl form (PorFe(IV)=O) that yields the hydroxo-ferryl isomer (PorFe(IV)-OH+). These are the only forms observed in experimental studies. Theory shows, however, that these are the least stable isomers of Compound II. The two most stable forms are the singlet and triplet states of the Por+*Fe(III)-OH electromer. In addition, theory reveals species never considered in heme enzymes: the singlet and triplet states of the Por+*Fe(III)-OH2 electromer. The computational results reproduce the experimental features of the known isomers and enable us to draw relationships and make predictions regarding the missing ones. For example, while the "surprise" species, singlet and triplet Por+*Fe(III)-OH2, have never been considered in heme chemistry, the calculated Fe-O bond lengths indicate that these isomers may have, in fact, been observed in one of the two opposing EXAFS studies reported previously. Furthermore, these ferric-aqua complexes could be responsible for the reported 18O exchange with bulk water. It is clear, therefore, that the role of Compound II in the HRP cycle is considerably more multi-faceted than has been revealed so far. Our suggested multi-state reactivity scheme provides a paradigm for Compound II species. The calculated M?ssbauer parameters may be helpful toward eventual characterization of these missing isomers of Compound II.     

Accurate ab initio quartic force fields of cyclic and bent HC2N isomers

Highly correlated ab initio quartic force fields (QFFs) are used to calculate the equilibrium structures and predict the spectroscopic parameters of three HC(2)N isomers. Specifically, the ground state quasilinear triplet and the lowest cyclic and bent singlet isomers are included in the present study. Extensive treatment of correlation effects were included using the singles and doubles coupled-cluster method that includes a perturbational estimate of the effects of connected triple excitations, denoted as CCSD(T). Dunning's correlation-consistent basis sets cc-pVXZ, X = 3,4,5, were used, and a three-point formula for extrapolation to the one-particle basis set limit was used. Core-correlation and scalar relativistic corrections were also included to yield highly accurate QFFs. The QFFs were used together with second-order perturbation theory (PT) (with proper treatment of Fermi resonances) and variational methods to solve the nuclear Schro?dinger equation. The quasilinear nature of the triplet isomer is problematic, and it is concluded that a QFF is not adequate to describe properly all of the fundamental vibrational frequencies and spectroscopic constants (though some constants not dependent on the bending motion are well reproduced by PT). On the other hand, this procedure (a QFF together with either PT or variational methods) leads to highly accurate fundamental vibrational frequencies and spectroscopic constants for the cyclic and bent singlet isomers of HC(2)N. All three isomers possess significant dipole moments, 3.05 D, 3.06 D, and 1.71 D, for the quasilinear triplet, the cyclic singlet, and the bent singlet isomers, respectively. It is concluded that the spectroscopic constants determined for the cyclic and bent singlet isomers are the most accurate available, and it is hoped that these will be useful in the interpretation of high-resolution astronomical observations or laboratory experiments.     

Photochemical E (trans)-->Z (cis) isomerization in substituted 1-naphthylacrylates.

Substituted naphthylacrylates, 1-3, not showing rotamerism have been synthesized with a view to study photochemical E (trans)-->Z (cis) isomerization. Photostationary state composition of the isomers upon direct excitation, triplet sensitized isomerization, quantum yield of isomerization, and steady state and time-resolved fluorescence behavior have been studied for these naphthylacrylates. The direct excitations of the compounds yield high Z (approximately 80%) isomer composition, whereas the triplet sensitization results in less Z (approximately 20%) isomer composition. This indicates that the singlet pathway is very efficient in converting the E isomer to the Z isomer. The naphthylacrylates 1 and 2 exhibit structured fluorescence at room temperature in hexane and upon changing the solvent to CH3CN; the structure of the fluorescence is lost, indicating that the singlet excited-state develops a polar character in a polar environment. The polar nature of the singlet excited state becomes more clear in the case of 3 from its fluorescence solvatochromism. The naphthylacrylates did not exhibit excitation wavelength-dependent fluorescence at room temperature suggesting that the ground state conformers (rotamers) are not involved. Fluorescence lifetimes measured for these compounds displayed biexponential behavior, which is explained using a two-state model.