3-甲基环状乙撑磷酸二酯醇解的反应途径研究

采用密度泛函理论和MP2方法研究了3-甲基环状乙撑磷酸二酯(MEP)与甲醇的反应途径: (I) CH3O－＋MEP; (II) CH3OH＋MEP; (III) CH3O－＋HMEP (MEP的质子化形式); (IV) CH3OH＋HMEP. 在B3LYP/6-31＋＋G(d,p)水平上优化了四条反应途径的反应物、中间体、过渡态及产物的几何构型, 并在同水平上进行了自然电荷分析, 然后在MP2/6-311＋＋G(3df,2p)水平上计算了各驻点的单点能. 采用极化连续介质模型(PCM)研究了各途径在苯、甲醇和水溶液中的溶剂化效应. 计算结果表明, 溶剂效应使途径(I)的自由能垒降低, 而使途径(II)和(IV)的决速步骤的自由能垒升高. 在气相和苯溶剂中途径(IV)是反应的优势途径, 在甲醇和水溶剂中途径(I)则成为最优. 研究结果进一步表明实验条件下途径(II)与(IV)对总醇解反应的贡献相当.     

大气中OH引发的有机磷农药MTMPDA降解机理的研究

采用量子化学方法在ωB97XD/6-311++G(3df,2pd)//ωB97XD/6-31+G(d,p)和M062X/6-311++G(3df,2pd)//M062X/6-31+G(d,p)水平下对大气中MTMPDA在OH自由基引发下的反应机理及其产物的后继反应进行理论研究.结果表明MTMPDA+OH反应主要通过OH抽提与N原子相连的—CH3上的H原子、以及抽提与O原子相连的—CH3上的H原子的方式进行.此外,抽提与N原子相连的—CH3上的H原子生成的产物与大气中的O2和NO发生进一步反应并没有导致MTMPDA毒性明显地降低;而抽提与O原子相连的—CH3上的H原子的产物发生一系列后继反应可以降低MTMPDA的毒性.     

S→N烷基重排反应的机理和溶剂化效应

采用密度泛函理论的B3LYP方法、从头算的MP2方法和自洽反应场极化连续模型(PCM),在6-311++G(2d,2p)基组水平上研究了N,N’-二甲基-S-异苯并呋喃在气相和溶液中发生S→N烷基重排反应的机理、溶剂效应和取代基效应.结果表明:该反应通过四元环机理和双位迁移机理生成产物,在气相和溶剂水中,双位迁移途径的能垒均比四元环途径低,反应主要通过双位迁移途径生成产物.在气相,苯环上发生-Cl,-NO2和-OCH3取代时,双位迁移途径的能垒在MP2/6-311++G(2d,2p)水平上比没有取代时分别低4.18,7.61,4.96kJ/mol,反应的取代基效应不明显.而在溶剂水中,苯环上发生-Cl,-NO2和-OCH3取代时,双位迁移途径的能垒在PCM-MP2/6-311++G(2d,2p)水平上比气相时分别低37.73,39.96和37.17kJ/mol,反应的溶剂化效应非常明显.理论研究结果与实验观察结果一致.     

硝基甲烷异构化反应势能面的ab initio研究

在B3LYP/6-311+ +G(2d,2p)水平上,优化得到硝基甲烷CH3NO2的10种异构体和23个异构化反应过渡态,并用G2MP2方法进行了单点能计算.根据计算得到的G2MP2相对能量,探讨了CH3NO2势能面上异构化反应的微观机理.研究表明,反应初始阶段的CH3NO2异构化过程具有较高的能垒,其中CH3NO2的两个主要异构化反应通道,即CH3NO2→CH3ONO和CH3NO2→CH2N(O)OH的活化能分别为270.3和267.8 kJ/mol,均高于CH3NO2的C-N键离解能.因而,从动力学角度考虑, CH3NO2的异构化反应较为不利.     

CH3+HNCO反应机理的理论研究

异氰酸(HNCO)分解引发的一系列自由基反应是氮氧化物快速消除机理[1,2](RAPRENOX)所研究的领域,该反应涉及到燃烧化学中氮氧化物NOX的消除,所以获得这些反应准确的位垒就成为实验化学和理论化学所要解决的问题。本文中我们重点研究CH3+HNCO反应机理,探讨CH3自由基是否也能象氮氢自由基一样,在异氰酸(HNCO)分解反应中起作用。1　计算方法用量子化学MP2方法,在6 311++G 水平上计算了CH3自由基与HNCO反应的反应物、产物、中间体和过渡态的几何构型,用QCISD(T)方法在6 311++G 水平上计算了它们的能量。通过振动分析确定…     

氨和甲醇团簇中的质子转移反应

应用激光多光子电离质谱和分子束技术研究了氨和甲醇二元团簇,实验观测到两个系列质子化的团簇离子: (CH3OH)nH+和(CH3OH)nNH4+(1≤n≤14 ),其产生是经过二元团簇内的质子转移反应。同时也研究了氘代甲醇CH3OD和氨混合团簇,结果表明OD原子团中的D转移概率比CH3原子团中的质子转移概率大几倍。在HF/STO-3G和MP2/6-31G*　*水平上对氨和甲醇二元团簇进行了计算,结果表明与CH3相比OH中的质子转移更加容易,因为CH3中的质子转移过程要克服高度约120 kJ/mol的能垒。     

CH3NO2+X(X=H,OH, CH3, CH2[^3B1], O[^3P])→CH2NO2+XH吸氢反应过渡态 结 构和反应位垒的DFT计算研究

采用DFT(B3LYP)方法,分别在6-311g(d,p),6-311++g(d,p)和自洽相关基组cc-pVIZ水平上优化了基态硝基甲烷和自由基H,OH,CH3,CH2[^3B1]以及O[^3P]等发生吸氢反应时的过渡态结构,并计算了反应的位垒。研究表明,对同一反应,不同基组下优化得到的过渡态几何结构基本一致;反应位垒数值的大小也基本接近,经校正,硝基甲烷同自由基反应位垒的理论计算值同实验结果基本吻合。     

氢甲酰化反应溶剂效应的量子化学研究

在B3LYP/6-31G(d, p)(Rh和P采用LANL2DZ + Polar)水平下,利用自洽反应场（SCRF）的Onsager模型,系统研究了C6H12、C6H6、THF、CH2Cl2、CH3OH和H2O（相对介电常数ε分别为2.02、2.25、7.58、8.93、32.63和78.39）六种溶剂对乙烯氢甲酰化反应中各驻点的结构和反应机理的影响.分别在上述溶剂中优化了反应物、产物、中间体和过渡态的结构,并计算了自由能和活化自由能.计算结果表明,不同溶剂对反应势能面上各驻点的结构和能量均有一定的影响,而且随着相对介电常数ε的增加,反应的活化自由能下降.水是以上六种溶剂中最佳的溶剂,这与大量的实验研究结果一致.     

甲烷-空气-H2O模拟光化学反应体系中CH3OO

在1.013×105 Pa,(298±2)K及O2-N2气氛下,研究了羟基自由基*OH引发的甲烷光化学反应体系中过氧甲基自由基CH3OO*自身复合反应.反应物和产物采用长光路Fourier红外光谱(LP-FTIR)和高效液相色谱(HPLC)测定.证实产物中有甲基过氧化氢(CH3OOH,MHP)和过氧甲醚(CH3OOCH3,DMP)存在并首次在该体系中发现了羟甲基过氧化氢(HOCH2OOH,HMHP).HMHP的检出表明,CH3OO*自身复合的可能途径之一生成了Criegee中间体过氧次甲基双自由基*CH2OO*.采用G2,G2(MP2)和G2(ful)方法对一些反应的标准焓变和标准Gibbs自由能变化进行了理论计算.结果表明CH3OO*自身复合生成*CH2OO*及*CH2OO*与H2O反应生成HMHP的途径在热力学上是可能的.     

F+CH3OH碰撞反应机理和反应势能面

在MP2(full)/6-311++g(d,p)水平上详细研究了氟原子与甲醇抽氢反应的多通道反应机理,得到了各条通道中涉及的驻点的构型和振动频率及其能量,给出了两张完整的反应势能面.结果表明,氟原子从C原子上抽氢时有一条明显的最低能量通道,而从氧原子上抽氢时要涉及多条分支通道和多个驻点构型,给出了各分支通道的势能面示意图,结果表明以形成五元环状过渡态通道为优势通道.计算得到经途径1生成CH2OH时反应放热170.62kJ/mol,经分支途径6生成CH3O自由基时反应放热119.41 kJ/mol,此结果与实验值一致.     

Reaction of CH3 with NO2

The elementary reaction of the CH3 radical with NO2 was investigated by time-resolved FTIR spectroscopy and quantum chemical calculations. The CH3 radical was produced by laser photolysis of CH3Br or CH3I at 248 nm. Vibrationally excited products OH, HNO and CO2 were observed by the time-resolved spectroscopy for the first time. The formation of another product NO was also verified. According to these observations, the product channels leading to CH3O+NO, CH2NO+OH and HNO+H2CO were identified. The channel of CH3O+NO was the major one. The reaction mechanisms of the above channels were studied by quantum chemical calculations at CCSD(T)/6-311++G(df,p)//MP2/6-311G(d,p) level. The calculated results fit with the experimental observations well.     

CH3CH2+O(3P)反应机理的理论研究

The reaction for CH3CH2+O(3P) was studied by ab initio method. The geometries of the reactants, intermediates, transition states and products were optimized at MP2/6-311+G(d,p) level. The corresponding vibration frequencies were calculated at the same level. The single-point calculations for all the stationary points were carried out at the QCISD(T)/6-311+G(d,p) level using the MP2/6-311+G(d,p) optimized geometries. The results of the theoretical study indicate that the major products are the CH2O＋CH3, CH3CHO+H and CH2CH2+OH in the reaction. For the products CH2O＋CH3 and CH3CHO+H, the major production channels are A1: (R)→IM1→TS3→(A) and B1: (R)→IM1→TS4→(B), respectively. The majority of the products CH2CH2+OH are formed via the direct abstraction channels C1 and C2: (R)→TS1(TS2)→(C). In addition, the results suggest that the barrier heights to form the CO reaction channels are very high, so the CO is not a major product in the reaction.     

Electronic and molecular structure of aminimides (1-acyl-2,2,2-trimethyldiazan-2-ium-1-ide). 1. Formaminimide (HCON- N+ Me3)

The electronic structure and geometries of (Z)- and (E)-H-CON- N+(CH3)3 have been examined at two levels of theory: B3LYP (basis sets 6-311+G(d,p), 6-311++G(d,p), and 6-311G(3df,3pd)) and MP2(full)/6-311++G(d,p). The (Z) conformation about the C(O)-N(-) bond is thermodynamically preferred over the (E) configuration. Natural bond orbital calculation locates one lone pair of the N- in the HOMO, which is the p(z) natural hybrid orbital (perpendicular to the O=CN- N+ plane). The second lone pair (of lower energy) of N- occupies the HOMO-3, which is the natural hybrid orbital sp(1.12) (sp(1.01) for the (E) conformation, sp(1.74) in the rotational transition state). The carbonyl pi bond is the HOMO-2. The charge-transfer ability of the negative nitrogen in H-CON- N+ (CH3)3 is more powerful than that of the neutral amidic nitrogen in dimethylformamide. The following facts convincingly sustain this view: (1) the higher rotational barrier (stronger C-N(-) bond) in the case of H-CON- N+ (CH3)3, (2) natural resonance theory analysis predicts almost equal weights for the (Z)-H-C(=O)N- N+ (CH3)3 and the (Z)-H-C(O-)=NN+ (CH3)3 canonical resonance structures whereas the weight of the HCON(CH3)2 structure is almost twice as large as that of HC(O-)=N+ (CH3)2, and (3) the second-order perturbation stabilization, as a result of the donor (N-)/acceptor (carbonyl) interaction, is 101.3 kcal/mol for H-CON- N+ (CH3)3 and only 64.4 kcal/mol for dimethylformamide.     

CH3 CH2 +N(4S)反应机理的理论研究

The reaction for CH3CH2+N(4S) was studied by ab initio method. The geometries of the reactants, intermediates, transition states and products were optimized at MP2/6-311+G(d,p) level. The corresponding vibration frequencies were calculated at the same level. The single point calculations for all the stationary points were carried out at the QCISD(T)/ 6-311+G(d,p) level using the MP2/6-311+G(d,p) optimized geometries. The results of the theoretical study indicate that the major products are the CH2CH2+3NH and H2CN＋CH3, and the minor products are the CH3CHN+H in the reaction. The majority of the products CH2CH2+3NH are formed via a direct hydrogen abstraction channel. The products H2CN＋CH3 are produced via an addition/dissociation channel. The products CH3CHN+H are produced via an addition/dissociation channel.     

N(4S)+CH3X(X=Cl、Br)反应机理的理论研究

The reaction of N(4S)+CH3X(X=Cl、Br) was studied by the ab initio method. The geometries of the reactants, transition states and products were optimized at the MP2/6-311+G(d,p) level. The corresponding vibration frequencies were calculated at the same level. The single-point calculations for all the stationary points were carried out at the MP2/6-311++G(3df,2p) and the QCISD(T)/6-311+G(d,p) levels using the MP2/6-311+G(d,p) optimized geometries. The energies of all the stationary points were calculated by the G2MP2 method. The results of this theoretical study indicate that the reaction has three reaction channels: H abstraction reaction channel a, Cl or Br abstraction reaction channel b and substitution reaction channel c. For the N(4S)+CH3Cl reaction, reaction channel a is the main reaction channel. Reaction channels b and c may have a slight contribution in the reaction. For the N(4S)+CH3Br reaction, reaction channel a is the main reaction channel. Reaction channels b and c may have some contribution in the reaction.     

Theoretical study and rate constant computation on the reaction HFCO + OH --> CFO + H2O

The potential energy surface, including the geometries and frequencies of the stationary points, of the reaction HFCO + OH is calculated using the MP2 method with 6-31+G(d,p) basis set, which shows that the direct hydrogen abstraction route is the most dominating channel with respect to addition and substitution channels. For the hydrogen abstraction reaction, the single-point energies are refined at the QCISD(T) method with 6-311++G(2df,2pd) basis set. The calculated standard reaction enthalpy and barrier height are -17.1 and 4.9 kcal mol(-1), respectively, at the QCISD(T)/6-311++G(2df,2pd)//MP2/6-31+G(d,p) level of theory. The reaction rate constants within 250-2500 K are calculated by the improved canonical variational transition state theory (ICVT) with small-curvature tunneling (SCT) correction at the QCISD(T)/6-311++G(2df,2pd)//MP2/6-31+G(d,p) level of theory. The fitted three-parameter formula is k = 2.875 x 10(-13) (T/1000)1.85 exp(-325.0/T) cm(3) molecule(-1) s(-1). The results indicate that the calculated ICVT/SCT rate constant is in agreement with the experimental data, and the tunneling effect in the lower temperature range plays an important role in computing the reaction rate constants.     

Mechanistic and kinetic study of the CH3CO + O2 reaction

Potential-energy surface of the CH3CO + O2 reaction has been calculated by ab initio quantum chemistry methods. The geometries were optimized using the second-order Moller-Plesset theory (MP2) with the 6-311G(d,p) basis set and the coupled-cluster theory with single and double excitations (CCSD) with the correlation consistent polarized valence double zeta (cc-pVDZ) basis set. The relative energies were calculated using the Gaussian-3 second-order Moller-Plesset theory with the CCSD/cc-pVDZ geometries. Multireference self-consistent-field and MP2 methods were also employed using the 6-311G(d,p) and 6-311++G(3df,2p) basis sets. Both addition/elimination and direct abstraction mechanisms have been investigated. It was revealed that acetylperoxy radical [CH3C(O)OO] is the initial adduct and the formation of OH and alpha-lactone [CH2CO2(1A')] is the only energetically accessible decomposition channel. The other channels, e.g., abstraction, HO2 + CH2CO, O + CH3CO2, CO + CH3O2, and CO2 + CH3O, are negligible. Multichannel Rice-Ramsperger-Kassel-Marcus theory and transition state theory (E-resolved) were employed to calculate the overall and individual rate coefficients and the temperature and pressure dependences. Fairly good agreement between theory and experiments has been obtained without any adjustable parameters. It was concluded that at pressures below 3 Torr, OH and CH2CO2(1A') are the major nascent products of the oxidation of acetyl radicals, although CH2CO2(1A') might either undergo unimolecular decomposition to form the final products of CH2O + CO or react with OH and Cl to generate H2O and HCl. The acetylperoxy radicals formed by collisional stabilization are the major products at the elevated pressures. In atmosphere, the yield of acetylperoxy is nearly unity and the contribution of OH is only marginal.     

Kinetics of the multichannel reaction of methanethiyl radical (CH3S*) with 3O2

The CH3S* + O2 reaction system is considered an important process in atmospheric chemistry and in combustion as a pathway for the exothermic conversion of methane-thiyl radical, CH3S*. Several density functional and ab initio computational methods are used in this study to determine thermochemical parameters, reaction paths, and kinetic barriers in the CH3S* + O2 reaction system. The data are also used to evaluate feasibility of the DFT methods for higher molecular weight oxy-sulfur hydrocarbons, where sulfur presents added complexity from its many valence states. The methods include: B3LYP/6-311++G(d,p), B3LYP/6-311++G(3df,2p), CCSD(T)/6-311G(d,p)//MP2/6-31G(d,p), B3P86/6-311G(2d,2p)//B3P86/6-31G(d), B3PW91/6-311++G(3df,2p), G3MP2, and CBS-QB3. The well depth for the CH3S* + 3O2 reaction to the syn-CH3SOO* adduct is found to be 9.7 kcal/mol. Low barrier exit channels from the syn-CH3SOO* adduct include: CH2S + HO2, (TS6, E(a) is 12.5 kcal/mol), CH3 + SO2 via CH3SO2 (TS2', E(a) is 17.8) and CH3SO + O (TS17, E(a) is 24.7) where the activation energy is relative to the syn-CH3SOO* stabilized adduct. The transition state (TS5) for formation of the CH3SOO adduct from CH3S* + O2 and the reverse dissociation of CH3SOO to CH3S* + O2 is relatively tight compared to typical association and simple bond dissociation reactions; this is a result of the very weak interaction. Reverse reaction is the dominant dissociation path due to enthalpy and entropy considerations. The rate constants from the chemical activation reaction and from the stabilized adduct to these products are estimated as functions of temperature and pressure. Our forward rate constant and CH3S loss profile are in agreement with the experiments under similar conditions. Of the methods above, the G3MP2 and CBS-QB3 composite methods are recommended for thermochemical determinations on these carbon-sulfur-oxygen systems, when they are feasible.     

Theoretical study of the structures, stability and vibrational spectra of the nitrous acid complexes with CH4

The structures, stability and vibrational spectra of the binary complexes CH4...HONO-trans and CH4...HONO-cis have been investigated using ab initio calculations at the SCF and MP2 levels with 6-311++G(d,p) basis set and B3LYP calculations with 6-31G(d,p) and 6-31+G(d,p) basis sets. Full geometry optimization was made for the complexes studied. It was established that the complex CH4...HONO-trans is more stable by 0.41 kcal mol(-1) than the complex CH4...HONO-cis. The accuracy of the ab initio calculations have been estimated by comparison between the predicted values of the vibrational characteristics (vibrational frequencies and infrared intensities) and the available experimental data. It was established, that the methods, used in this study are well adapted to the problem under examination. The predicted values with the B3LYP calculations are very near to the results, obtained with 6-311++G(d,p)/MP2. The changes in the vibrational characteristics of methane and trans-, cis-nitrous acid upon formation of the hydrogen bond show that the complexes CH4...HONO-trans and CH4...HONO-cis have geometry in which the OH group interacts with a methane molecule forming a single hydrogen bond. This fact is confirmed by relatively strong perturbation of the OH stretching vibration to lower frequencies and an increase of the infrared intensity of this vibration up to three times upon hydrogen bonding.