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.     

Quantum mechanical investigation of the atmospheric reaction CH3O2 + NO

The important stationary points on the potential energy surface of the reaction CH(3)O(2) + NO have been investigated using ab initio and density functional theory techniques. The optimizations were carried out at the B3LYP/6-311++G(d,p) and MP2/6-311++G(d,p) levels of theory while the energetics have been refined using the G2MP2, G3//B3LYP, and CCSD(T) methodologies. The calculations allow the proper characterization of the transition state barriers that determine the fate of the nascent association conformeric minima of methyl peroxynitrite. The main products, CH(3)O + NO(2), are formed through either rearrangement of the trans-conformer to methyl nitrate and its subsequent dissociation or via the breaking of the peroxy bond of the cis-conformer to CH(3)O + NO(2) radical pair. The important consequences of the proposed mechanism are (a) the allowance on energetic grounds for nitrate formation parallel to radical propagation under favorable external conditions and (b) the confirmation of the conformational preference of the homolytic cleavage of the peroxy bond, discussed in previous literature.     

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.     

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.     

Quantum chemical study of the reaction of trichloroethylene with O(3P)

The adiabatic mechanism of the reaction of trichloroethylene with O(³P), exploring the various O-atom addition and H-atom abstraction channels, is theoretically studied at the MP2/6-311++G(2d, 2p), MP2/aug-cc-pVTZ, CCSD/6-31G(d), G3, and CBS-QB3 levels of theory. From a kinetic point of view, the addition to the less substituted carbon atom of the double bond is more favorable than the addition to the more substituted carbon. Such O-atom addition reactions are favored over the one possible hydrogen-abstraction reaction. Calculations of the present study showed that five products are obtained: HCCl + C(O)Cl₂ (P1), Cl + ClC(O)CHCl (P2), H + ClC(O)CCl₂ (P3), Cl + HC(O)CCl₂ (P4), and CH(O)Cl + CCl₂ (P5). The products P2 and P4 are found to be the most favored ones. The kinetic calculations of rate constant in the range of 285–395 K are performed at the CBS-QB3 level of theory and are in conformity with the experimental outcomes.     

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

The reaction for SiH3+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 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 pathway is the SiH3+O(3P)→IM1→TS3→IM2→TS8→HOSi+H2. The other minor products include the HSiOH+H, H2SiO+H and HSiO+H2. Furthermore, the products HOSi, HSiO and HSiOH(cis) can undergo dissociation into the product SiO. In addition, the calculations provide a possible interpretation for disagreement about the mechanism of the reaction SiH4+O(3P). It suggests that the products HSiOH, H2SiO and SiO observed by Withnall and Andrews are produced from the secondary reaction SiH3+O(3P) and not from the reaction SiH4+O(3P).     

Transient infrared spectra of CH3SOO and CH3SO observed with a step-scan Fourier-transform spectrometer

A step-scan Fourier-transform spectrometer coupled with a multipass absorption cell was employed to monitor time-resolved infrared absorption of transient species produced upon irradiation at 248 nm of a flowing mixture of CH(3)SSCH(3) and O(2) at 260 K. Two transient bands observed with origins at 1397±1 and 1110±3 cm(-1) are tentatively assigned to the antisymmetric CH(3)-deformation and O-O stretching modes of syn-CH(3)SOO, respectively; the observed band contour indicates that the less stable anti-CH(3)SOO conformer likely contributes to these absorption bands. A band with an origin at 1071±1 cm(-1), observed at a slightly later period, is assigned to the S=O stretching mode of CH(3)SO, likely produced via secondary reactions of CH(3)SOO. These bands fit satisfactorily with vibrational wavenumbers and rotational contours simulated based on rotational parameters of syn-CH(3)SOO, anti-CH(3)SOO, and CH(3)SO predicted with density-functional theories B3LYP/aug-cc-pVTZ and B3P86/aug-cc-pVTZ. Two additional bands near 1170 and 1120 cm(-1) observed at a later period are tentatively assigned to CH(3)S(O)OSCH(3) and CH(3)S(O)S(O)CH(3), respectively; both species are likely produced from self-reaction of CH(3)SOO. The production of SO(2) via secondary reactions was also observed and possible reaction mechanism is discussed.     

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.     

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

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

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.     

Theoretical investigation on the GaCl3-catalyzed ring-closing metathesis reaction of N-2,3-butadienyl-2-propynyl-1-amine: three-membered ring versus four-membered ring mechanism

The gallium chloride (GaCl(3))-catalyzed ring-closing metathesis reaction mechanism of N-2,3-butadienyl-2-propynyl-1-amine has been studied at the Becke three-parameter hybrid functional combined with Lee-Yang-Parr correlation functional (B3LYP)/6-31G(d), B3LYP/6-31+G(d,p), B3LYP/6-311++G(d,p)//B3LYP/ 6-31G(d) and the second-order M?ller-Plesset perturbation (MP2)/6-311++G(d,p)//B3LYP/6-31+G(d,p) levels. It was found that the final metathesis product can be yielded via a three-membered or four-membered ring mechanism. The three-membered ring pathway is favorable due to its low energy barrier at the rate determining step. The whole reaction is stepwise and strongly exothermic.     

CH3CH2S自由基H迁移异构化及裂解反应的理论研究

采用密度泛函方法(B3LYP)在6-311+G(d,p)基组水平上研究了CH3CH2S自由基H迁移异构化以及裂解反应的微观动力学机理. 在QCISD(T)/6-311++G(d,p)//B3LYP/6-311+G(d,p)+ZPE水平上进行了单点能校正. 利用经典过渡态理论(TST)与变分过渡态理论(CVT)分别计算了在200~2000 K温度区间内的速率常数kTST和kCVT, 同时获得了经小曲率隧道效应模型(SCT)校正后的速率常数kCVT/SCT. 研究结果表明, CH3CH2S自由基1,2-H迁移、1,3-H迁移、C—C键断裂和β-C—H键断裂反应的势垒ΔE≠分别为149.74, 144.34, 168.79和198.29 kJ/mol. 当温度低于800 K时, 主要发生1,2-H迁移反应, 高于1800 K时, 主要表现为C—C键断裂反应, 在1300—1800 K范围内, 1,3-H迁移反应是优势通道, 在计算的整个温度段内, β-C—H键断裂反应可以忽略.     

The shape of gaseous n-butylbenzene: assessment of computational methods and comparison with experiments

Conformational landscape of neutral and ionized n-butylbenzene has been examined. Geometries have been optimized at the B3LYP/6-31G(d), B3LYP/6-31+G(d,p), B3LYP-D/6-31+G(d,p), B2PLYP/6-31+G(d,p), B2PLYP-D/6-31+G(d,p), B97-D/6-31+G(d,p), and M06-2X/6-31+G(d,p) levels. This study is complemented by energy computations using 6-311++G(3df,2p) basis set and CBS-QB3 and G3MP2B3 composite methods to obtain accurate relative enthalpies. Five distinguishable conformers have been identified for both the neutral and ionized systems. Comparison with experimentally determined rotational constants shows that the best geometrical parameters are provided by B3LYP-D and M06-2X functionals, which include an explicit treatment of dispersion effects. Composite G3MP2B3 and CBS-QB3 methods, and B2PLYP-D, B3LYP-D, B97-D, and M06-2X functionals, provide comparable relative energies for the two sets of neutral and ionized conformers of butyl benzene. An exception is noted however for conformer V(+) the stability of which being overestimated by the B3LYP-D and B97-D functionals. The better stability of neutral conformers I, III, and IV, and of cation I(+) , demonstrated by our computations, is in perfect agreement with conclusions based on micro wave, fluorescence, and multiphoton ionization experiments.     

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~＇)自由 基与臭氧之间的反应活性很强。     

A new scheme for determining the intramolecular seven-membered ring N-H...O=C hydrogen-bonding energies of glycine and alanine peptides

In this paper a new scheme was proposed to calculate the intramolecular hydrogen-bonding energies in peptides and was applied to calculate the intramolecular seven-membered ring N-H...O=C hydrogen-bonding energies of the glycine and alanine peptides. The density-functional theory B3LYP6-31G(d) and B3LYP6-311G(d,p) methods and the second-order Moller-Plesset perturbation theory MP26-31G(d) method were used to calculate the optimal geometries and frequencies of glycine and alanine peptides and related structures. MP26-311++G(d,p), MP26-311++G(3df,2p), and MP2/aug-cc-pVTZ methods were then used to evaluate the single-point energies. It was found that the B3LYP6-31G(d), MP26-31G(d), and B3LYP6-311G(d,p) methods yield almost similar structural parameters for the conformers of the glycine and alanine dipeptides. MP2/aug-cc-pVTZ predicts that the intramolecular seven-membered ring N-H...O=C hydrogen-bonding strength has a value of 5.54 kcal/mol in glycine dipeptide and 5.73 and 5.19 kcal/mol in alanine dipeptides, while the steric repulsive interactions of the seven-membered ring conformers are 4.13 kcal/mol in glycine dipeptide and 6.62 and 3.71 kcal/mol in alanine dipeptides. It was also found that MP26-311++G(3df,2p) gives as accurate intramolecular N-H...O=C hydrogen-bonding energies and steric repulsive interactions as the much more costly MP2/aug-cc-pVTZ does.     

Vibrational frequencies and infrared intensities of the hydrogen-bonded complexes of nitrous acid with ethers: ab initio and DFT studies

The vibrational spectra of the binary complexes formed by HONO-trans and HONO-cis with dimethyl and diethyl ethers 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 optimisation was made for the complexes studied. 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 ab initio and DFT calculations show that the changes in the vibrational characteristics (vibrational frequencies and infrared intensities) upon hydrogen bonding for the hydrogen-bonded complex (CH3)2O...HONO-trans are larger than for the complex (CH3)2O...HONO-cis.     

Energetics of the O-H bond and of intramolecular hydrogen bonding in HOC6H4C(O)Y (Y = H, CH3, CH2CH=CH2, C[triple bond]CH, CH2F, NH2, NHCH3, NO2, OH, OCH3, OCN, CN, F, Cl, SH, and SCH3) compounds

The energetics of the phenolic O-H bond in a series of 2- and 4-HOC 6H 4C(O)Y (Y = H, CH3, CH 2CH=CH2, C[triple bond]CH, CH2F, NH2, NHCH 3, NO2, OH, OCH3, OCN, CN, F, Cl, SH, and SCH3) compounds and of the intramolecular O...H hydrogen bond in 2-HOC 6H 4C(O)Y, was investigated by using a combination of experimental and theoretical methods. The standard molar enthalpies of formation of 2-hydroxybenzaldehyde (2HBA), 4-hydroxybenzaldehyde (4HBA), 2'-hydroxyacetophenone (2HAP), 2-hydroxybenzamide (2HBM), and 4-hydroxybenzamide (4HBM), at 298.15 K, were determined by micro- or macrocombustion calorimetry. The corresponding enthalpies of vaporization or sublimation were also measured by Calvet drop-calorimetry and Knudsen effusion measurements. The combination of the obtained experimental data led to Delta f H m (o)(2HBA, g) = -238.3 +/- 2.5 kJ.mol (-1), DeltafHm(o)(4HBA, g) = -220.3 +/- 2.0 kJ.mol(-1), Delta f H m (o)(2HAP, g) = -291.8 +/- 2.1 kJ.mol(-1), DeltafHm(o)(2HBM, g) = -304.8 +/- 1.5 kJ.mol (-1), and DeltafHm(o) (4HBM, g) = -278.4 +/- 2.4 kJ.mol (-1). These values, were used to assess the predictions of the B3LYP/6-31G(d,p), B3LYP/6-311+G(d,p), B3LYP/aug-cc-pVDZ, B3P86/6-31G(d,p), B3P86/6-311+G(d,p), B3P86/aug-cc-pVDZ, and CBS-QB3 methods, for the enthalpies of a series of isodesmic gas phase reactions. In general, the CBS-QB3 method was able to reproduce the experimental enthalpies of reaction within their uncertainties. The B3LYP/6-311+G(d,p) method, with a slightly poorer accuracy than the CBS-QB3 approach, achieved the best performance of the tested DFT models. It was further used to analyze the trends of the intramolecular O...H hydrogen bond in 2-HOC 6H 4C(O)Y evaluated by the ortho-para method and to compare the energetics of the phenolic O-H bond in 2- and 4-HOC 6H 4C(O)Y compounds. It was concluded that the O-H bond "strength" is systematically larger for 2-hydroxybenzoyl than for the corresponding 4-hydroxybenzoyl isomers mainly due to the presence of the intramolecular O...H hydrogen bond in the 2-isomers. The observed differences are, however, significantly dependent on the nature of the substituent Y, in particular, when an intramolecular H-bond can be present in the radical obtained upon cleavage of the O-H bond.     

Nucleophilic substitution at the imidoyl carbon atom: intermediate mechanistic and reactivity behavior between carbonyl and vinyl carbon substitution.

Gas-phase nucleophilic substitution reactions at the imidoyl carbon have been investigated using chloride exchanges, Cl- + RY=CHCl right harpoon over left harpoon RY=CHCl + Cl- with Y = N and R = F, H or CH3, at the MP2, B3LYP and G2(+) levels using the MP2/6-311+G geometries. The results are compared with those for the vinyl (Y = CH) and carbonyl (Y = O) carbon substitution. The mechanism and reactivity of substitution at the imidoyl carbon are intermediate between those of carbonyl (SNpi) and vinyl carbon (SNsigma) substitution, which is directly related to the electronegativity of Y, CH < N < O. The prediction of competitive SNsigma with SNpi path for the imidoyl chloride is consistent with the S(N)1-like mechanism proposed for reactions in solution. The important factors in favor of an in-plane concerted SN2 (SNsigma) over an out-of-plane pi-attack (SNpi) path are (i) lower proximate sigma-sigma* charge-transfer energies (DeltaECT), (ii) stronger electrostatic stabilization (DeltaENCT), and (iii) larger lobe size on C(alpha) for the sigma*- than pi*-LUMO despite the higher sigma* than pi* level. The electron correlation energy effects at the MP2 level are overestimated for the relatively delocalized structure (S(N)pi TS) but are underestimated for the localized structure (SNsigma TS) so that the MP2 energies lead to a wrong prediction of preferred reaction path for the vinyl chloride. The DFT at the B3LYP level predicts correct reaction pathways but overestimates the electron correlation effects.     

Theoretical Investigation of CH3CF2O2+HOO Reaction

The reactions of CH3CF2O2 with HOO are important chemical cyclic processes of photochemical contamination. In this paper, the reaction pathways and reaction mechanism of CH3CF2O2+HOO are investigated extensively with the Gaussian 98 package at the B3LYP/6-311++G** basis sets. The use of vibrational mode analysis and electron population analysis to reveal the reaction mechanism is firstly reported. The study shows that CH3CF2CO2+HOO→IM1→TS1→CH3CF2O2H+O2 channel is the energetically most favorable, CH3CF2CO2H and O2 are the principal products, and the formation of CH3OH and CF2O is also possible.     

FC(O)O自由基与NO2反应的量子化学研究

用量子化学DFT, MP2, G3和G3MP2方法对FC(O)O自由基与NO₂的反应机理进行了理论研究. 优化了反应势能面上各驻点的几何结构, 通过内禀反应坐标(IRC)计算和振动分析, 确认了反应中的过渡态, 并用过渡态理论(TST)计算了相关反应的速率常数.