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

在B3LYP／6－311G（d,p)和QCISD（T）／6－311＋＋G（3df,2p)（单点）水平下，对B2H6^ 阳离子和B2H5自由基全优化得到9个几何异构体，B2H5^＋单态体系（D3h,C1),B2H5^ 三重态（Cs,Cs,C1),B2H5自由基（C2v,Cs,Cs,Cs)，得到势能面上与体系异构化过程相联系的5种过渡态。     

乙烯基自由基(C2H·3)与羟基自由基(OH·)反应势能面的理论研究

在QCISD（T）/6-311G（2df,p)//B3LYP/6-311G（d,p)水平上对自由基反应C2H3^. OH^.进行了计算，结果表明，经过缔合、多步H转移、CH3转移和离解等复杂过程，最终要得到8种产物（P1-P8），茯中产物P2（H2CCO H2）和P6（CH3CO^. H^.）是主要产物。本文得到的CH2CHOH（1或1‘），CH3CHO（2）和CH3COH（3）之间的过渡态TS1/2，TS1‘/3和TS2/3的能量顺序与Wesdemiotis等的实验推测相反，而与Smith等的计算结果一致。     

乙烯酮自由基(HCCO·)与氧气(O2)反应势能面的理论研究

在G2(B3LYP／MP2／CC)水平上对反应HCCO O2进行了计算,得到了反应势能面,提出了3种可能的反应机理：(1)四元环反应机理得到产物P1(HCO CO2);(2)三元环反应机理得到产物P2(CO HCO2);(3)O-O键断裂反应机理得到产物P3(O OCC(O)H)和P4(O CO HCO)．由反应势能面推测产物P1(HCO CO2)为主要产物,产物P2(CO HCO2),P3(O OCC(O)H)P4(O CO HCO)为次要产物．     

HCO+NO2反应势能面的理论研究

在B3LYP/6-311G(d,p)和CCSD(T)/6-311G(d,p)水平上给出了HCO+NO₂反应详细的势能面信息.计算结果表明,该反应采用两种无垒进攻方式,分别得到两种加合物H(O)CNO₂和H(O)CONO.找到7种能量低于反应物且合理的产物及相应的反应路径.通过对热力学和动力学的分析,产物HONO+CO(P2,P3),HNO+CO₂(P1)和H+CO₂+NO(P6)的形成更为有利.计算结果同实验相符,且有助于深入了解HCO自由基的化学行为.     

亚甲基氰自由基(HCCN)与一氧化氮(NO)反应势能面的理论研究

在CCSD(T)/6-311G(d,p)//B3LYP/6-311G(d,p)+ZPE水平上对反应HCCN+NO的二重态反应势能面进行了计算,得到了4种产物:P1(HCN+NCO),P2(OH+NCCN),P₃[HCN+(CNO)]和P4(HCN+CNO).其中产物P1为主要产物,P2为次要产物,P₃和P4很难得到.在G2(B3LYP/MP2/CC)水平,对产物P1和P2的反应通道的单点能量进行了校正.     

乙烯酮自由基(HCCO·)与氧气(O2)反应势能面的理论研究

在G2(B3LYP/MP2/CC)水平上对反应HCCO+O2进行了计算,得到了反应势能面,提出了3种可能的反应机理:(1)四元环反应机理得到产物P1(HCO+CO2);(2)三元环反应机理得到产物P2(CO+HCO2);(3)O—O键断裂反应机理得到产物P3(O+OCC(O)H)和P4(O+CO+HCO).由反应势能面推测产物P1(HCO+CO2)为主要产物,产物P2(CO+HCO2),P3(O+OCC(O)H)和P4(O+CO+HCO)为次要产物.     

O+HCNO反应势能面的理论研究

采用CCSD(T)/6-311G(d,p)//B3LYP/6-311G(d,p)+ZPVE方法对反应O＋HCNO进行了研究. 通过反应势能面揭示了该反应的机理, 通过H或O迁移等多步反应路径得到3种产物, 其中, P1(HCO+NO)为主要产物, P2(HNO+CO)和P3(NCO+OH)为次要产物. 为进一步实验研究提供了参考.     

H+HCNO反应势能面的理论研究

在CCSD(T)/6-311G(d,p)//B3LYP/6-311G(d,p)+ZPVE水平下, 对反应H＋HCNO进行了研究. 建立了反应势能面, 揭示了该反应的反应机理, 通过H迁移、N—O键或C—N键断裂等多步反应, 得到4种产物, 其中最主要产物为P1(HCN+OH).     

乙烯酮自由基(HCCO·)与二氧化氮(NO2)反应势能面的理论研究

在CCSD(T)/6-311G(d,p)//MP2/6-311G(d,p)+ZPE水平上对反应HCCO+NO2进行了计算, 建立了反应势能面. 此反应由反应物通过三步反应到达产物. 首先, NO2的O原子进攻HCCO自由基中与H相邻的C原子, 形成异构体1[ONOC(H)CO]或2[H(CONOC)O]. 然后, 异构体1和2通过N-O键的断裂形成产物NO和OC(H)CO. 最后, 产物中的OC(H)CO可以通过C-C键的断裂进一步分解为HCO和CO. 由HCCO+NO2反应得到产物NO+HCO+CO.     

氧气和CS自由基反应势能面的密度泛函理论研究

应用量子化学从头计算和密度泛函理论(DFT)对O_2和CS自由基的反应进行了研 究。在B3LYP/6-311G~(**)水平上计算出了各物种的优化构型、振动频率和零点振 动能(ZPVE)。各物种的总能量由CCSD（T）/6-311G~(**)//B3LYP/6-311G~(**)给出 ，并对总能量进行了零点能校正。计算结果表明：CS自由基中的C端沿着O_2的双键 中线方向进攻，进行加成反应，反应的第一步放出大量的热量(450 kJ/mol)，推动 反应继续进行，从稳定的中间体4(Cs)出发，反应主要通过O的相邻位置的迁移生成 P_1(CO+SO)和P_3(COS+O)；通过S的相邻位置的迁移生成了重要的反应复合物 (complex 1)，进一步离解为产物P_2(CO_2+S)。计算结果可以很好地解释反应机理 。     

Theoretical study of the Si2NO potential energy surface

The structures, spectroscopies, and stabilities of the doublet Si₂NO radical are explored at the density functional theory (DFT) and ab initio levels. Seventeen isomers are located, connected by 26 interconversion transition states. At the CCSD(T)/6‐311+G(2df)//QCISD/6‐311G(d)+ZPVE level, three low‐lying isomers are predicted, that is, one bent species SiNSiO 3 (5.1 kcal/mol) containing the important Si?N triple bonding and two four‐membered ring isomers including cyclic cSiNSiO 1 (0.0) with Si? Si cross‐bonding with C_2v symmetry and puckered cSiNSiO 1′ (11.9) with divalent carbene character. Three low‐lying isomers 1, 1′, and 3 have reasonable kinetic stabilities and might be observable either experimentally or astrophysically. The possible formation strategies of 1, 1′, and 3 in laboratory and in space are discussed in detail. The calculated vibrational frequencies and possible formation processes of 3 are consistent with recent experimental observations. In light of the fact that no cyclic nitrogen‐containing species have been detected in space, two cyclic isomers 1 and 1′ could be promising candidates. Furthermore, the bonding nature of three isomers 1, 1′, and 3 is analyzed. The calculated results are also compared with those of the analogue C₂NO radical. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

Theoretical study on the potential energy surface and vibrational energy levels of HXeI

The potential energy surface for the electronic ground state of the HXeI molecule is constructed by using the internally contracted multi-reference configuration interaction with the Davidson correction（icMRCI＋Q）method and large basis sets. The stabilities and dissociation barriers are identified from the potential energy surfaces.The three-body dissociation channel is found to be the dominate dissociation channel for HXeI.Based on the obtained potentials,vibrational energy levels of HXeI are calculated using the Lanczos algorithm.Our theoretical results are in excellent agreement with the available observed values.     

Theoretical mechanistic study on the reaction of CN radical with HNCN

The mechanism for the reaction of the cyanogen radical (CN) with the cyanomidyl radical (HNCN) has been investigated theoretically. The electronic structure information of the singlet and triplet potential energy surfaces (PESs) is obtained at the B3LYP/6-311+G(3df,2p) level, and the single-point energies are refined at the CCSD(T)/6-311+G(3df,2p) level as well as by multilevel MCG3-MPWB method. The calculations show that the C atom of CN additions to middle- and end-N atoms of HNCN are two barrierless association processes leading to the energy-rich intermediates IM1 HN(CN)CN and IM2 HNCNCN, respectively, on the singlet PES. The higher barriers of the subsequent isomerization and dissociation channels from IM1 and IM2 indicate that these two intermediates, which have considerably thermodynamic and kinetic stability, are the dominant product at high pressure. While at low pressure, the most favorable product is P(2) H + NCNCN, which will be formed from both IM1 and IM2 via direct dissociation processes by the H-N bond rupture, and the secondary feasible product is P(4) HCN + (1) NCN, while P(5) HCCN + N(2) and P(6) HCNC + N(2) are the least competitive products. On the triplet PES, P(14) NCNC + HN may be a comparable competitive product at high temperature. In addition, the comparison between the mechanisms of the CN + HNCN and OH + HNCN reactions is made. The present results will enrich our understanding of the chemistry of the HNCN radical in combustion processes and interstellar space.     

自由基HO2与C2H2反应势能面的理论研究

应用量子化学从头计算和密度泛函理论(DFT)对HO2+C2H2反应体系的反应机理进行了研究.在B3LYP/6-311G**和CCSD(T)/6-311G**水平上计算了HO2+ C2H2反应的二重态反应势能面.计算结果表明,主要反应方式为自由基HO2的H原子和C2H2分子中的C原子结合,经过一系列异构化,最后分解得到主要产物P1 (CH2O+ HCO).此反应是放热反应,化学反应热为-321.99 kJ·mol-1.次要产物为P2 (CO2 +CH3),也是放热反应.     

Theoretical study on the [Si, C, N, O] potential energy surface

The structures, energetics, spectroscopies, and stabilities of the doublet [Si, C, N, O] radical are explored at the density functional theory and ab initio levels. Sixteen isomers are located, connected by 29 interconversion transition states. At the CCSD(T)/6-311+G(2df)//QCISD/6-311G(d)+ZPVE level, the lowest lying isomer is a linear SiNCO 1 (0.0 kcal/mol) mainly featuring a cumulene | . Si = N = C = O. The second and third low-lying isomers are bent OSiCN 2 (8.8) and bent OSiNC 3 (11.1), respectively. All the three low-lying isomers 1, 2, 3, and another high-lying species 5 (75.4) with a linear SiCNO structure are shown to have considerable kinetic stability and may be experimentally observable. The predicted results of isomers 1 and 2 are consistent with the previous mass spectrometry experiments. Moreover, the fourth low-lying species SiOCN 4 (23.9) with bent structure is expected to be observable in low-temperature environments. The bonding nature of the five isomers 1, 2, 3, 4, and 5 is analyzed. The calculated results are compared with those of the analogous molecules C(2)NO and Si(2)NO. Implications in interstellar space and N,O-doped SiC vaporization processes are also discussed.     

Theoretical studies on the potential energy surfaces and vibrational energy levels of HXeF and HXeCl

The potential energy surfaces for the electronic ground state of the HXeCl and HXeF molecules areconstructed by using the internally contracted multi-reference configuration interaction with theDavidson correction(icMRCI Q)method and large basis sets.The stabilities and dissociation barriersare identified from the potential energy surfaces.The three-body dissociation channel is found to bethe dominate dissociation channel for HXeCl,while two dissociation channels are possible and com-petitive for HXeF.Based on the obtained potentials,vibrational energy levels of HXeCl and HXeF arecalculated using the Lanczos algorithm.Our theoretical results are in good agreement with the avail-able observed values.Particularly,the calculated fundamental frequency of the H—Xe stretching vi-bration including the Xe matrix effect of HXeCl is found to be 1666.6 cm-1,which is only 17.6 cm-1higher than the recently observed value of 1649 cm-1.     

Theoretical study on reaction mechanism of the vinyl radical with nitrogen atom

The complex triplet potential energy surface of the C₂H₃N system is investigated at the UB3LYP and CCSD(T) (single-point) levels in order to explore the possible reaction mechanism of C₂H₃ radical with N(⁴S). Eleven minimum isomers and 18 transition states are located. Possible energetically allowed reaction pathways leading to various low-lying dissociation products are obtained. Starting from the energy-rich reactant C₂H₃+N(⁴S), the first step is the attack of the N atom on the C atom having one H atom attached in C₂H₃ radical and form the intermediate C₂H₃N(1). The associated intermediate 1 can lead to product P₁ CH₂CN+H and P₂ ³CH₂+³HCN by the cleavage of C–H bond and C–C bond, respectively. The most favorable pathway for the C₂H₃+N(⁴S) reaction is the channel leading to P₁, which is preferred to that of P₂ due to the comparative lower energy barrier. The formation of P₃ ³C₂H₂+³NH through hydrogen-abstraction mechanism is also feasible, especially at high temperature. The other pathways are less competitive comparatively.     

Exploring the potential energy surface of retinal, a comparison of the performance of different methods

The ground state structure of retinal has been investigated. We found that DFT and CASSCF produce different results for the bond length alternation in a model system of retinal. Quantum mechanics/molecular mechanics calculations including the closest surrounding amino acids have been performed, using DFT and CASSCF to calculate the structure of retinal in the protein cavity. The planarity of the retinal molecule is affected by the surrounding protein. DFT and CASSCF produce different twist angles. The difference between CASSCF and DFT appears to be related to the positively charged nitrogen of the Schiff base, which leads to different pi-bond orders produced by the two methods.     

Theoretical study on reaction mechanism of the CF radical with nitrogen dioxide

The singlet potential energy surface of the [CFNO₂] system is investigated at the B3LYP and CCSD(T) (single‐point) levels to explore the possible reaction mechanism of CF radical with NO₂. The top attack of C‐atom of CF radical at the N‐atom of NO₂ molecule first forms the adduct isomer FCNO₂ 1 followed by oxygen‐shift to give trans‐OC(F)NO 2 and then to cis‐OC(F)NO 3 . Subsequently, the most favorable channel is a direct dissociation of 2 and 3 to product P₁ FCO+NO. The second and third less favorable channels are direct dissociation of 3 to product P₂ FNO+CO and isomerization of 3 to a complex NOF?CO 4 , which can easily dissociate to product P₃ FON+CO, respectively. The large exothermicity released in these processes further drives most of the three products P₁ , P₂ , and P₃ to take secondary dissociation to the final product P₁₂ F+CO+NO. Another energetically allowed channel is formation of product P₄ ¹NF+CO₂, yet it is much less competitive than P₁ , P₂ , P₃ , and P₁₂ . The present calculations can well interpret one recent experimental fact that the title reaction is quite fast yet still much slower than the analogous reaction CH+NO₂. Also, the results presented in this article may be useful for future product distribution analysis of the title reaction as well as for the analogous CCl and CBr reactions. © 2001 John Wiley & Sons, Inc. J Comput Chem 22: 1907–1919, 2001     

Theoretical study of nonlinear triatomic molecular potential energy surfaces: Lie algebraic method

Triatomic molecular potential energy surfaces (PES) are obtained by using coherent state to take the classical limits of algebraic Hamiltonian. The algebraic Hamiltonian for bent tria-tomic molecules can be obtained using Lie algebraic method (the expansion coefficients are obtained by fitting spectroscopic data). This PES is applied to H2O molecule, and good results are obtained.