Theoretical investigation on stability and isomerizations of CH3SO isomers

The stability and isomerizations of CH3SO isomers have been investigated at B3LYP/6-311G(d,p), MP2/6-311G(d,p), QCISD/6-311G(d,p), and CCSD(T)/6-311G(d,p) levels. Geometries of isomers and transition states (TS) have been optimized at the B3LYP/6-311G(d, p) level. Vibration analysis and the intrinsic reaction coordinate (IRC) calculated at the same level have been applied to validate the connection of the stationary points. The four different methods give similar results: 11 isomers and 9 isomerization channels were found. CH3SO and CH2(S)OH are the most stable species among the 11 isomers. Furthermore, the breakage and formation of the chemical bonds in isomerization reactions have been discussed by the topological analysis method of electronic density. The "energy transition state (ETS)" and the "structure transition state (STS)" of all the isomerizations have been found. The topological analysis shows that the relative positions of ETS and STS are determined by reaction energy. The nonplanar four-member ring structure transition state (STS), which was first found in this paper, extended the concept of ring STS.     

Topological studies on IRC paths of the isomerization reaction of silicon methyl-nitrene

On the basis of kinetic study of isomerization reaction of H_3Si-N, ab initio (RHF, UHF/6-31G) calculations on some points of the singlet and triplet reaction paths were carried out. The breakage and formation of chemical bond in the reaction are discussed. The calculated results show that there is a transitional structure of three-membered ring on each of reaction paths. A 'structural transition region' and a 'structural transition state' in both of studied reaction are found. Our previous conclusion that the structure transition state (STS) always appears before the energy transition state (ETS) in endothermic reaction and after ETS in exothermic reaction is further confirmed. The relationship between the change of spin density distribution and the structural transition state are investigated.     

CH2XH→CH3X (X＝O, S, Se)异构化的量子拓扑研究

利用从头算和量子拓扑方法讨论了CH₂XH→CH₃X (X＝O, S, Se)异构化过程的反应机理. 着重从电子密度拓扑分析计算了反应进程中的各点, 讨论了反应进程中键的断裂和生成, 上述反应都经历了三元环过渡结构, 找到了这类反应的"能量过渡态"和"结构过渡态", 且结构过渡态均在能量过渡态之后出现. 三元结构过渡态结构出现的范围与反应热成正比.     

IClO2异构化反应机理及电子密度拓扑研究

利用密度泛函理论方法研究了IClO2异构化反应机理. 优化得到了七种异构体, 其中OIClO和IClOO还未见报道, 对各异构体的热力学稳定性进行了比较. 找到了异构化过程的过渡态, 并通过内禀反应坐标(IRC)计算确认了各个异构体之间的相互转化关系. 从量子拓扑学的角度, 对典型异构化反应通道IRC途径上的各点进行了电子密度拓扑分析, 讨论了反应过程中化学键的断裂、生成以及化学键的变化规律, 找到了反应途径上的能量过渡态(ETS)和结构过渡态(STS).     

HCHO+X(X=F、Cl、Br)的反应机理

采用CCSD/6-311++G(d,p)//B3LYP/6-311++G(d,p)方法研究了HCHO与卤素原子X(X=F、Cl、Br)的反应机理. 计算结果表明, 卤素原子X(X=F、Cl、Br)主要通过直接提取HCHO中的H原子生成HCO+HX(X=F、Cl、Br). 另外还可以生成稳定的中间体, 中间体再通过卤原子夺氢和氢原子直接解离两个反应通道分别生成HCO+HX(X=F、Cl、Br)和H+XCHO(X=F、Cl、Br). 其中卤原子夺氢通道为主反应通道, HCO和HX(X=F、Cl、Br)为主要的反应产物; 且三个反应的活化能均较低, 说明此类反应很容易进行, 计算结果与实验结果符合很好. 电子密度拓扑分析显示, 在HCHO+X反应通道(b)中出现了T型结构过渡态, 结构过渡态(STS)位于能量过渡态(ETS)之后. 并且按F、Cl、Br的顺序, 结构过渡态出现得越来越晚.     

几个P-ylide反应机理的量子拓扑研究

采用MP2(FC)/6-311++G(d,p)对磷叶立德和类磷叶立德自由基反应机理进行了探讨.优化了中间体、过渡态和产物的几何构型,并采用内禀反应坐标法进行追踪.侧重从量子拓扑学的角度,对反应过程中各点进行电子密度拓扑分析,讨论了反应过程中化学键的断裂、生成和化学键的变化规律.上述两个反应都经历三员环过渡态,找到了这类反应的能量过渡态和结构过渡态.     

Topological characteristics of electron density distribution in SSXY --> XSSY (X or Y = F, Cl, Br, I) isomerization reactions

The reactions of SSXY to XSSY (X or Y = F, Cl, Br, I) have been studied at B3LYP/6-311++G(2df) and MP2/6-311++G(2df) levels based on B3LYP/6-311++G(2df) optimized geometries. There are two pathways (X or Y atom transferring) for each reaction of SSXY to XSSY(X or Y = F, Cl, Br, I). The "Atoms in Molecules" theory has been applied to analyze the topological characteristics of electron density distribution along the reaction path. (a) The S-S-X or S-S-Y ring structure transition region (STR) is contributed to S(1)-X or S(1)-Y bond formation and S(2)-X or S(2)-Y bond annihilation. The STR and structure transition state (STS) are defined. (b) Comparing the two pathways of each SSXY --> XSSY (X, Y = F, Cl, Br, I) reaction, X or Y atom transferring, the broader the STR is, the later the STS appears, and the pathway is easier. (c) When X linked to the same S site of a three-member ring (S-S-F, S-S-Cl, S-S-Br, or S-S-I), the ring STR becomes broader and the STS appears later according to the sequence of X = F, Cl, Br, and I. And in these cases, for exothermic reactions, the higher the DeltaE(P-->TS1) is, the broader the STR is and the later the STS appears. (d) When the same Y atom linked to different three-member-ring S-S-X, the STR becomes broader and the STS appears later according to the sequence of X = F, Cl, Br, and I in the ring. And in these cases, for exothermic reactions, the lower the DeltaE(P-->TS2) is, the broader the STR is and the later the STS appears.     

ClO与ClO自由基反应机理及电子密度拓扑分析

利用密度泛函理论对ClO与ClO自由基反应机理进行了深入理论探讨,在B3LYP/6-311＋＋G（3df）水平上对该反应体系的反应物、中间体、过渡态及产物进行了几何构型优化,对反应通道进行了IRC(内禀反应坐标）路径解析,计算了沿各反应通道的能垒和离解能,并进行了零点能校正.从量子拓扑学的角度,对反应通道IRC途径上一些重要点进行了电子密度拓扑分析,讨论了反应过程中化学键的断裂、生成以及键的变化规律,找到了反应途径的能量过渡态和结构过渡态.     

几个硫-叶立德反应机理的量子拓扑研究

采用MP2(FC)/6-311++G(d,p)对硫叶立德和类硫叶立德自由基反应机理进行了探讨。优化了中间体、过滤态和产物的几何构型。本文侧重从量子拓扑学的角度,对IRC(内禀反应坐标)反应进程中各点进行电子密度拓扑分析,讨论了反应过程中化学键的断裂、生成和化学键的变化规律。找到了这类反应的能量过渡态和结构过渡态,上述两个反应都是先经历一个没有形成三元环拓扑结构的能量过渡态,再经历一个形成了三元环拓扑结构的结构过渡态,最后到达产物。     

Studies on reactions INCX --> IXCN (X = O, S, and Se)

The reactions of INCO to IOCN, INCS to ISCN, and INCSe to ISeCN have been studied at the MP2/6-311++G(2df)//B3LYP/6-311++G(2df) level. Geometries of reactants, transition states, and products have been optimized, and geometries of transition states are reported for the first time. The reasons that INCO, ISCN, and ISeCN are easily detected instead of IOCN, INCS, and INCSe have been explained successfully. The breakage and formation of chemical bonds along the reaction paths have been discussed by the topological analysis method of electronic density. The calculated results show that there are two kinds of structure transition states (STS) in the reactions studied.     

BrO与CH3SH反应机理的量子化学及拓扑研究

利用密度泛函和电子密度拓扑分析方法对BrO与CH₃SH反应的微观机理进行了理论研究. 在B3LYP/6-311G (d, p)水平上对反应势能面上的各驻点进行几何构型的全优化; 振动分析和IRC计算证实了中间体和过渡态的真实性和相互连接关系; 计算得到了各反应通道的活化能, 并进行了零点能校正. 计算结果表明: 该反应存在7个反应通道, 其中生成CH₃S＋HOBr和CH₃SO＋HBr的通道为主要反应通道. 通过对反应过程中部分驻点的电子密度拓扑分析, 首次发现了接近平面的四元环状过渡态, 从而拓展了原来对环状结构过渡态定义的适用范围.     

RN (R＝CH3, CH3CH2)异构化及脱氢反应的量子拓扑研究

利用量子化学从头算CASSCF方法在6-311＋G (d, p)基组水平上对单线态和三线态RN (R＝CH₃, CH₃CH₂)异构化反应及RN脱氢反应的微观机理进行了理论研究. 在MP2/6-311＋G (d, p)和CCSD/6-311＋G (d, p)水平上进行了单点能校正. 单态和三态势能面的交叉点(ISC)的存在清楚地说明了基态反应物³RN异构化为基态产物¹R'NH (R'＝CH₂, CH₃CH)的过程. 电子密度拓扑分析显示在整个异构化过程中有两种类型的结构过渡态: 单态反应通道为T型过渡态, 三态反应通道为环状过渡态. 单线态RN脱氢反应通道中“原子-分子键”的存在说明两个H原子是以H₂的形式从RN中脱去的.     

CH3S自由基H迁移异构化及脱H2反应的直接动力学研究

采用密度泛函方法(MPW1PW91)在6-311G(d,p)基组水平上研究了CH₃S自由基H迁移反应CH₃S→CH₂SH (R1), 脱H₂反应CH₃S→HCS＋H₂ (R2)以及脱H₂产物HCS异构化反应HCS→CSH (R3)的微观动力学机理. 在QCISD(t)/6- 311＋＋G(d,p)//MPW1PW91/6-311G(d,p)＋ZPE水平上进行了单点能校正. 利用经典过渡态理论(TST)与变分过渡态理论(CVT)分别计算了各反应在200～2000 K温度区间内的速率常数k^TST和k^CVT, 同时获得了经小曲率隧道效应模型(SCT)校正后的速率常数k^CVT/SCT. 结果表明, 反应 R1, R2 和R3的势垒△E^≠分别为160.69, 266.61和241.63 kJ/mol, R1为反应的主通道. 低温下CH₃S比CH₂SH稳定, 高温时CH₂SH比CH₃S更稳定. 另外, 速率常数计算结果显示, 量子力学隧道效应在低温段对速率常数的计算有显著影响, 而变分效应在计算温度段内对速率常数的影响可以忽略.     

HNCX→HXCN(X=O, S, Se)反应的量子拓扑研究

The reactions of HNCO to HOCN, HNCS to HSCN and HNCSe to HSeCN have been studied at MP2/6-311 ＋ ＋ G（2df, pd）//B3LYP/6-311 ＋ ＋G（2df, pd） level. Geometries of the reactants, transition states and products have been optimized and geometries of the transition states are reported for the first time. The reasons why HNCO and HNCS instead of HOCN and HSCN were easily detected have been explained. It was predicted that HNCSe will be more easily detected than HSeCN. The breakage and formation of the chemical bonds in the reactions have been discussed by the topological analysis method of electronic density. The calculated results show that there are two kinds of structure transition states （STS） in reactions studied.     

HCN→HNC异构化反应过渡态附近的拓扑特性

电荷密度拓扑分析方法(topological analysis of charge distribution)已经被广泛应用于研究化合物的各种静态性质及反应特性.对IRC反应途径上过渡态附近构型的拓扑性质研究报导较少.本工作是用电荷密度分布的拓扑分析方法对HCN→HNC异构化反应IRC过渡态附近的性质进行了研究.对电荷密度分布拓扑分析和IRC途径确定选用相同的基组,得到了一些新的结论。     

Approach to the atmospheric chemistry of methyl nitrate and methylperoxy nitrite. Chemical mechanisms of their formation and decomposition reactions in the gas phase

Potential energy surfaces, minimum energy reaction paths, minima, transition states, reaction barriers, and conical intersections for the most important atmospheric reactions of methyl nitrate (CH(3)ONO(2)) and methylperoxy nitrite (C(3)HOONO) on the electronic ground state have been studied (i) with the second-order multiconfigurational perturbation theory (CASPT2) by computation of numerical energy gradients for stationary points and (ii) with the density functional theory (DFT). The proposed mechanism explains the conversion of unreactive alkyl peroxy radicals into alkoxy radicals: CH(3)O(2) + NO <=> CH(3)OONO <=> CH(3)O + NO(2) left arrow over right arrow CH(3)ONO(2). Additionally, several discrepancies found in the comparison of the results obtained from the two employed approaches are analyzed. CASPT2 predicts that all dissociation reactions into radicals occur without an extra exit energy barrier. In contrast, DFT finds transition states for the dissociations of cis- and trans-methylperoxy nitrite into CH(3)O + NO(2). Furthermore, multiconfigurational methods [CASPT2 and complete active space SCF (CAS-SCF)] predict the isomerization of CH3ONO2 to CH3OONO to occur in a two-step mechanism: (i) CH(3)ONO(2) --> CH(3)O + NO(2); and (ii) CH(3)O + NO(2) --> CH(3)OONO. The reason for this has to do with the coupling of the ground electronic state with the first excited state. Therefore, it is demonstrated that DFT methods based on single determinantal wave functions give an incorrect picture of the aforementioned reaction mechanisms.     

卤代硅烯异构化为卤代甲硅烷基硅烯的热力学与动力学分析

利用密度泛函理论(DFT) 中的B3LYP方法,全参数优化了卤代硅烯HXSi=SiXH（X=F、Cl、Br、I）异构化反应的反应物、产物及过渡态的几何构型,计算出了它们的振动频率、零点振动能(ZPVE)和总能量,并对它们进行了振动分析,以确定过渡态的真实性。又计算了反应的热力学函数变化,平衡常数及速率常数,比较了不同卤素对反应的影响。热力学与动力学计算结果表明,该异构化反应过程是一个放热的、且在低温下可自发进行的反应,但对于溴代硅烯和碘代硅烯而言,当温度达到1000 K时,反应开始转化为非自发反应。     

Modeling the structure,absorption spectra,and cis-trans isomerization of thiacarbocyanine dyes

The relative stability of the trans-and cis-isomers of 3,3′-diethylthiacarbocyanine (Dye1) and 3,3′-diethyl-9-methylthiacarbocyanine (Dye2)¹, as well as sections of the potential energy surfaces along the internal coordinate of the isomerization reaction, were studied using the density functional theory. Calculation of the minimum energy pathway for the isomerization reaction showed that the barrier for rotation about the C₈–C₉ bond is higher for Dye1 than for Dye2. Local minimums were found for the singlet excited state of the 8,9-cis-and trans-isomers of the dyes. In the case of the trans-isomers, substantial changes in the dye structure do not occur and the local minimum of the excited state corresponds to the geometry of the starting trans-isomers, which favors efficient fluorescence. A search for the nearest local minimum of the singlet excited state of the 8,9-cis-isomers leads to structures, which differ significantly from the starting structures, and the intensity of the S₁ → S₀ transition in those structures appears to be practically zero. The results are in agreement with experimental data on the absorption, fluorescence, and fluorescence excitation spectra of the dyes.     

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.