CH_3O与CIO双自由基反应机理的量子化学研究

在QCISD(T)/6-311+G(d,p)//B3LYP/6-311+G(3df,3pd)水平上,对CH3O与ClO双自由基反应进行了理论研究.结果表明,该反应共有三个反应通道,产物分别为HOCI+CH2O,CH2O2+HCl和CH3Cl+O2(1△).不论从动力学角度,还是从热力学角度看,形成产物HOCl+CH2O的通道均是最有利的,因此为主要反应通道,这与实验观察到的结果是一致的.     

CH_2CH(~2A')自由基与臭氧反应机理的理论研究

用量子化学MP2（full）方法，在6－311+ +G~(**)基组水平上研究了CH_2CH (~2A~＇)自由基与臭氧反应的机理，全参数优化了反应过程中反应物、中间体、过 渡态和产物的几何构型，在QCISD（T,full）/6-311+ +G~(**)水平上计算了它们的 能量，并对它们进行了振动分析，以确定中间体和过渡态的真实性，研究结果表明 ：CH_2CH(~2A~＇)自由基与臭氧反应有两条可行的反应通道，分别为：CH_2CH (~2A~＇)+O_3→TS1→M1→TS2→O_2+OCH_2CH→TS4+O_2→O_2(~3∑_g)+CH_2CHO (~2A~＂)和CH_2CH(~2A~＇)+O_3→M2→TS3→O_2(~3∑_g)+CHO(~2A~＂),后一个反 应通道较容易发生，而且反应活化能小（2.97kJ/mol），说明CH_2CH(~2A~＇)自由 基与臭氧之间的反应活性很强。     

氧气和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)。计算结果可以很好地解释反应机理 。     

CH3O与ClO双自由基反应机理的量子化学研究

在QCISD(T)/6-311＋G(d,p)//B3LYP/6-311＋G(3df,3pd)水平上, 对CH₃O与ClO双自由基反应进行了理论研究. 结果表明, 该反应共有三个反应通道, 产物分别为HOCl＋CH₂O, CH₂O₂＋HCl和CH₃Cl＋O₂(¹Δ). 不论从动力学角度, 还是从热力学角度看, 形成产物HOCl＋CH₂O的通道均是最有利的, 因此为主要反应通道, 这与实验观察到的结果是一致的.     

CH_2CH(~2A')自由基与臭氧反应机理的理论研究

用量子化学MP2（full）方法,在6－311+ +G~(**)基组水平上研究了CH_2CH (~2A~＇)自由基与臭氧反应的机理,全参数优化了反应过程中反应物、中间体、过 渡态和产物的几何构型,在QCISD（T,full）/6-311+ +G~(**)水平上计算了它们的 能量,并对它们进行了振动分析,以确定中间体和过渡态的真实性,研究结果表明 ：CH_2CH(~2A~＇)自由基与臭氧反应有两条可行的反应通道,分别为：CH_2CH (~2A~＇)+O_3→TS1→M1→TS2→O_2+OCH_2CH→TS4+O_2→O_2(~3∑_g)+CH_2CHO (~2A~＂)和CH_2CH(~2A~＇)+O_3→M2→TS3→O_2(~3∑_g)+CHO(~2A~＂),后一个反 应通道较容易发生,而且反应活化能小（2.97kJ/mol）,说明CH_2CH(~2A~＇)自由 基与臭氧之间的反应活性很强。     

(H_2O)_n(n=0～2)催化HS和HOCl反应机理的理论研究

在CCSD(T)/6-311+G(3df,2p)//M06-2X/6-311+G(3df,2p)水平上研究了(H_2O)n(n=0~2)催化HS和HOCl的反应机理.结果表明,HS与HOCl反应中HS夺取HOCl上的H原子形成产物H_2S和ClO.在无水催化时,该反应存在2种不同的路径(分别经过过渡态TS1和TS2,二者互为顺反结构),对应的能垒分别为100.28和100.91kJ/mol,到达产物(H_2S+ClO)需吸收18.99kJ/mol能量,反应不易发生;在单个水分子参与时,水分子可通过形成弱相互作用或者作为H原子转移桥梁影响反应机理,获得了4种水催化路径,能垒(间于53.97~92.39kJ/mol之间)均低于无水催化过程.同时发现,在反应到达产物前,水分子可以与产物形成中间体IM,IM相对能仅为0.46kJ/mol,有利于产物形成;有2个水分子参与反应时,找到了3条催化路径,最优反应路径过渡态TS7的能垒为45.05kJ/mol,低于无水催化过程,相比单个水分子最优路径能垒(53.97kJ/mol)并无显著降低.     

卟吩内氢迁移反应的理论研究

用B3LYP/6-31G~(* *)方法计算卟吩内氢迁移反应中的反应物、产物、中间体 、过渡态及其二级鞍点的结构与能量、结果表明，分步反应历程的几率较大，与实 验的预测一致，同时对几种理论研究的结果进行了分析比较。     

ClO/ClO~-耦合体系结构性质及反应通道研究

利用从头算法和密度泛函理论对ClO/ClO~-体系进行了研究。首先应用密度泛 函理论的六种方法(B3LYP,BLYP,B3P86,BP86,BHLYP,LSDA)和从头算理论的CCSD方法 在6-3+G~*,6-311+G~*及6-311+G(3df)基组水平上对单体进行了优化,选出最适合 该体系的方法和基组B3LYP/6-311+G(3df);然后在B3LYP/6-311+G(3df)水平上计算 了沿各种反应通道的离解能,并且进行了校正。找出了存在的中间体及中间体异构 化的过渡态,并进行了IRC路径解析。对各稳定体进行了频率分析和成键分析。结 果表明,单体ClO和单体ClO~-结合为稳定的中间体后,其离解方式主要是向着 ClOO+Cl~-和ClOO~-+Cl两种方式进行,两种离解方式的离解能分别为-33.39和82. 88kJ/mol,并且前者是一个离解性电子转移过程,经历一个电子转移过渡态。     

O（~1D）与CF_2HCl反应的理论研究

用量子化学密度泛函理论（DFT）和G3B3方法，对为（~1D）与CF_2HCl的反应 进行了研究，在B3LYP/6-311+G(d)，B3LYP/6-311+G(2df,2pd)和G3B3计算水平上， 优化了反应热能面上各驻点的几何结构，通过内禀反应坐标（IRC）计算和振动分 析，对反应过渡态进行了确认，并确定了反应机理。     

甲苯定向硝化的理论研究

运用密度泛函理论(DFT)B3LYP方法，在6-31G~(* *)基组水平上，全优化计算 了硝酰阳离子NO_2~+对苯和从0-，m-,p-位进攻甲苯亲电取代硝化后应，求得4条反 应途径上包括反应物、过渡态和Wheland中间体共(4×3)12个反应驻点σ络合物的 分子几何、电子结构、能量和IR光谱等性质，阐明了反应中无同位素效应的实验事 实，求得各反应途径的活化能排序：p-＞o-＞PhH＞m-和σ络合物(R,TSak或INT)的 相对稳定化能排序：p-ArCH_3-NO_2~+＞o-ArCH_3-NO_2~+＞m-ArCH_3-NO_2~+＞ PhH-NO_2~+，从而阐明了甲基对苯环致活(或致钝)以及增加甲苯硝化络合物稳定性 的双重功能，对甲苯定向硝化的理论预示与实验结果相吻合。     

Ab initio study of the ClO + NH2 reaction: prediction of the total rate constant and product branching ratios

The mechanism for ClO + NH2 has been investigated by ab initio molecular orbital and transition-state theory calculations. The species involved have been optimized at the B3LYP/6-311+G(3df,2p) level and their energies have been refined by single-point calculations with the modified Gaussian-2 method, G2M(CC2). Ten stable isomers have been located and a detailed potential energy diagram is provided. The rate constants and branching ratios for the low-lying energy channel products including HCl + HNO, Cl + NH2O, and HOCl + 3NH (X(3)Sigma(-)) are calculated. The result shows that formation of HCl + HNO is dominant below 1000 K; over 1000 K, Cl + NH2O products become dominant. However, the formation of HOCl + 3NH (X(3)Sigma(-)) is unimportant below 1500 K. The pressure-independent individual and total rate constants can be expressed as k1(HCl + HNO) = 4.7 x 10(-8)(T(-1.08)) exp(-129/T), k(2)(Cl + NH2O) = 1.7 x 10(-9)(T(-0.62)) exp(-24/T), k3(HOCl + NH) = 4.8 x 10(-29)(T5.11) exp(-1035/T), and k(total) = 5.0 x 10(-9)(T(-0.67)) exp(-1.2/T), respectively, with units of cm(3) molecule(-1) s(-1), in the temperature range of 200-2500 K.     

OXO(X=Cl,Br)与H(2S)反应机理的量子化学研究

卤素二氧化物自由基(OXO)是重要的同温层物种，在极地的夜间含量尤其丰富，由于其光解产生的卤素原子会造成臭氧的严重损耗，因此愈来愈受到人们的广泛关注嘲．OXO自由基可以通过自由基之间的反应或微波放电等方式获得，有关它与其它大气物种的反应已经进行了许多研究，Bemand等在实验上研究了Cl，H，O与OClO的反应动力学，对于OBrO与H之间的反应目前尚未见实验研究的报道，只有Guhua等在理论上对部分XBrO2(X=H，Cl，Br)同分异构体进行了研究．     

Ab initio/density functional theory and multichannel RRKM study for the ClO + CH2O reaction

The potential energy surface for the CH(2)O + ClO reaction was calculated at the QCISD(T)/6-311G(2d,2p)//B3LYP/6-311G(d,p) level of theory. The rate constants for the lower barrier reaction channels producing HOCl + HCO, H atom, OCH(2)OCl, cis-HC(O)OCl and trans-HC(O)OCl have been calculated by TST and multichannel RRKM theory. Over the temperature range of 200-2000 K, the overall rate constants were k(200-2000K) = 1.19 x 10(-13)T(0.79) exp(-3000.00/T). At 250 K, the calculated overall rate constant was 5.80 x 10(-17) cm(3) molecule(-1) s(-1), which was in good agreement with the experimental upper limit data. The calculated results demonstrated that the formation of HOCl + HCO was the dominant reaction channel and was exothermic by 9.7 kcal/mol with a barrier of 5.0 kcal/mol. When it retrograded to the reactants CH(2)O + ClO, an energy barrier of 14.7 kcal/mol is required. Furthermore, when HOCl decomposed into H + ClO, the energy required was 93.3 kcal/mol. These results suggest that the decomposition in both the forward and backward directions for HOCl would be difficult in the ground electronic state.     

Structures and Electronic Properties of HOCl···HCOCl Complexes

B3LYP/6-311++G** and MP2/6-311++G** calculations were used to analyze the interaction between hypochlorous acid (HOCl) and formyl chloride (HCOCl). The results showed that there were four equilibrium geometries (S1, S2, S3, and S4) optimized at B3LYP/6-311++G** level, and all the equilibrium geometries were confirmed to be in stable states by analytical frequency calculations. Complexes S1 and S3 use the 5H atom of HOCl as proton donor and the terminal 1O atom of HCOCl as acceptor to form red shift hydrogen bond systems. However, the blue-shifted hydrogen bond (2C-3H···6O) coexists with 4Cl···5O interaction in structures S2. As for S4, it uses the 7Cl atom of HOCl as proton donor and the terminal 1O atom of HCOCl as acceptor to form red shift halogen bond system. Interaction energies between monomers in the four complexes corrected with basis set superposition error (BSSE) and zero-point vibrational energy (ZPVE) lie in the range from −5.05 to −14.76 kJ·mol⁻¹ at MP2/6-311++G** level. The natural bond orbital (NBO) and atoms in molecules (AIM) theories have also been applied to explain the structures and the properties of the complexes.     

HOCl…HCOCl复合物的结构和电子性质

在DFT-B3LYP/6-311++G**水平上求得HOCl+HCOCl复合物势能面上的四种稳定构型(S1, S2, S3和S4). 其中, 在复合物S1和S3中, HOCl单体的5H原子作为质子供体, 与HCOCl单体中作为质子受体的1O原子相互作用, 形成红移氢键复合物; 在复合物S4中, HOCl单体的7Cl原子作为质子供体, 与HCOCl单体中作为质子受体的1O原子相互作用, 形成红移卤键复合物; 而在复合物S2中, 同时存在2C—3H…6O蓝移氢键和4Cl…5O相互作用. 在MP2/6-311++G**水平上计算的单体间的相互作用能考虑了基组重叠误差(BSSE)和零点振动能(ZPVE)校正, 其值在-5.05与-14.76 kJ·mol-1之间. 采用自然键轨道理论(NBO)对两种单体间相互作用的本质进行了考查, 并通过分子中原子理论(AIM)分析了复合物中氢键和卤键键鞍点处的电子密度拓扑性质.     

OXO (X=Cl,Br)与X (2P3/2)反应机理的理论研究

用密度泛函B3LYP/6-311+G**和高级电子相关的组态相互作用QCISD(T)/6-311+G**方法研究了OXO与X (2P3/2)双自由基反应的微观机理.研究结果表明该反应存在两个反应通道,产物分别为XO和X2+O2.由于形成产物XO的活化势垒较低,因而是主要反应通道,这与实验观察到的结果是一致的.而形成X2+O2的通道从动力学上看是不利的.     

On the aromatic character of the heterocyclic bases of DNA and RNA

Studies based on ab initio optimized geometries (at B3LYP/6-311+G** and MP2/6-311+G** levels) and on experimental structures retrieved from the Cambridge Structural Database (CSD) reveal that the nucleobases constituting DNA and RNA differ significantly in their aromatic character, as shown by the geometry-based index of aromaticity HOMA that ranges from 0.466 for thymine to 0.917 for adenine, based on B3LYP/6-311+G** calculations, and 0.495-0.926, respectively, if based on the MP2/6-311+G** level. Aromaticity of the bases decreases markedly with an increase of the number of double-bond C=X (X = N, O) substituents at the rings. H-bonds involving C=O groups in Watson-Crick pairs cause an increase of the aromatic character of the rings.     

Theoretical studies on the reaction of pentafulvenone with water in gas phase and aqueous solvent

In gas phase, the hydrations of pentafulvenone to generate three types of cyclopentadienyl carboxylic acids are studied theoretically at the MP2/6-311+G**//B3LYP/6-311+G** level. A water molecule attacking the C=O double bond of pentafulvenone can yield cyclopentadienyl carboxylic acids via the formation of fulvenediols, and attacking the C=C double bond of pentafulvenone can directly yield cyclopentadienyl carboxylic acid. The barriers of rate-determining transition states are 42.2 and 30.4 kcal mol⁻¹, respectively. The barriers of rate-determining transition states for two water molecules system are 20.2 and 19.6 kcal mol⁻¹, respectively. The products can isomerize to each other. In aqueous solvent, the hydrations of pentafulvenone are investigated using PCM-UAHF model at the MP2 (PCM)/6-311+G**// B3LYP (PCM)/6-311+G** and MP2 (PCM)/6-311+G**// B3LYP/6-311+G** levels. The barriers of all rate-determining transition states are decreased. The added water molecule acts as catalyst in both gas phase and aqueous solvent. Electronic Supplementary Material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Theoretical Study of Reaction Mechanism of 1-Propenyl Radical with NO

The reaction system of 1-propenyl radical with NO is an ideal model for studying the intermolecular and intramolecular reactions of complex organic free radicals containing C=C double bonds. On the basis of the full optimization of all species with the Gaussian 98 package at the B3LYP/6-311++G** level, the reaction mechanism was elucidated extensively using the vibrational mode analysis. There are seven reaction pathways and five sets of small molecule end products: CH2O+CH3CN, CH2CHCN+H2O, CH3CHO+HCN, CH3CHO+HNC, and CH3CCH+HNO. The channel of C3H5￠+NO→ IM1→TS1→IM2→TS2→IM3→TS3→CH3CHO+HCN is thermodynamically most favorable.