DFT investigation of the mechanism of CH2CO + O(3P) reaction

A detailed theoretical survey of the potential energy surface (PES) for the CH₂CO + O(³P) reaction is carried out at the QCISD(T)/6‐311+G(3df,2p)//B3LYP/6‐311+G(d,p) level. The geometries, vibrational frequencies, and energies of all stationary points involved in the reaction are calculated at the B3LYP/6‐311+G(d,p) level. More accurate energy information is provided by single‐point calculations at the QCISD(T)/6‐311+G(3df,2p) level. Relationships of the reactants, transition states, intermediates, and products are confirmed by the intrinsic reaction coordinate (IRC) calculations. The results suggest that P1(CH₂+CO₂) is the most important product. This study presents highlights of the mechanism of the title reaction. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2005     

Direct dynamics study on mechanism and kinetics of the biradical self‐reaction of HOO

A theoretical study of the mechanism and the kinetics for the hydrogen abstraction reaction of the biradical hydroperoxy radical has been presented at the CCSD(T)/6‐311++G(3d,2p)//CCSD/6‐31+G(d,p) level of theory. Our theoretical calculations suppose a stepwise mechanism involving the formation of a postreactant complex in the triplet and singlet entrance channels. Four transition states of the six‐membered chain complexes (³TS1 and ¹TS1) and six‐membered ring complexes (³TS2 and ¹TS2) are located at the high dual level CCSD(T)/6‐311++G(3d,2p)//CCSD/6‐31+G(d,p) method. The rate constants of Path 1 ～ Path 4 at the CCSD(T)/6‐311++G(3d,2p)//CCSD/6‐31+G (d,p) level are calculated by means of the conventional transition state theory (TST) and canonical variational TST without and with small‐curvature tunneling (SCT) correction within the temperature range of 200–2,500 K. The calculated results show that the triplet channel is the dominating reaction channel and Path 2 is found to be the most favorable pathway. The rate constants of Path 2 are in good agreement with the experimental values at the experimentally measured temperatures. Moreover, the variational effect is not obvious in the low temperature range but is not neglectable in the high temperature range. The SCT plays an important role particularly in the low temperature range. © 2010 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

DFT studies on the multi‐channel reaction of CH3S+NO2

The mechanisms for the reaction of CH₃S with NO₂ are investigated at the QCISD(T)/6‐311++G(d,p)//B3LYP/6‐311++G(d,p) on both single and triple potential energy surfaces (PESs). The geometries, vibrational frequencies, and zero‐point energy (ZPE) correction of all stationary points involved in the title reaction are calculated at the B3LYP/6‐311++G(d,p) level. More accurate energies are obtained at the QCISD(T)/6‐311++G(d,p). The results show that 5 intermediates and 14 transition states are found. The reaction is more predominant on the single PES, while it is negligible on the triple PES. Without any barrier height for the whole process, the main channel of the reaction is to form CH₃SONO and then dissociate to CH₃SO+NO. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

Ab initio study on the reaction mechanism of ozone with bromine atom

Ab initio calculations of the potential energy surface (PES) for the Br+O₃ reaction have been performed using the MP2, CCSD(T), and QCISD(T) methods with 6‐31G(d), 6‐311G(d), and 6‐311+G(3df). The reaction begins with a transition state (TS) when the Br atom attacks a terminal oxygen of ozone, producing an intermediate, the bromine trioxide (M), which immediately dissociates to BrO+O₂. The geometry optimizations of the reactants, products, and intermediate and transition states are carried out at the MP2/6‐31G(d) level. The reaction potential barrier is 3.09 kcal/mol at the CCSD(T)/6‐311+G(3df)//MP2 level, which shows that the bromine atom trends intensively to react with the ozone. The comparison of the Br+O₃ reaction with the F+O₃ and Cl+O₃ reactions indicates that the reactions of ozone with the halogen atoms have the similar reaction mechanism. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

Model Study of the Photochemical Rearrangement Pathways of 1,2,4‐Oxadiazole

The mechanisms of photochemical isomerization reactions are investigated theoretically by using a model system of 1,2,4‐ oxadiazole with the CAS(14,9)/6‐311G(d) and MP2‐CAS‐(14,9)/ 6‐311++G(3df,3pd)//CAS(14,9)/6‐311G(d) methods. Three reaction pathways are examined, including 1) the direct mechanism, 2) the ring contraction–ring expansion mechanism, and 3) the internal cyclization–isomerization mechanism, which lead to two types of photoisomers. The theoretical findings suggest that conical intersections play a crucial role in the photorearrangement of 1,2,4‐oxadiazoles. These model investigations also indicate that the preferred reaction route for 1,2,4‐oxadiazole, which leads to phototransposition products, is as follows: reactant → Franck‐Condon region → conical intersection → photoproduct. In other words, the direct mechanism is a one‐step process that has no barrier. These theoretical results agree with the available experimental observations.     

An investigation on the reaction mechanism of the F2+Cl2→2ClF using the B3LYP method

The reaction mechanism of F₂+Cl₂→2ClF has been investigated with the density functional theory at the B3LYP/6‐311G* level. Six transition states have been found for the three possible reaction paths and verified by the normal mode vibrational and IRC analyses. Ab initio MP2/6‐311G* geometry optimizations and CCSD(T)/6‐311G(2df)//MP2/6‐311G* single‐point energy calculations have been performed for comparison. It is found that when the F₂ (or Cl₂) molecule decomposes into atoms first and then the F (or Cl) atom reacts with the molecule Cl₂ (or F₂) nearly along the molecular axis, the energy barrier is very low. The calculated energy barrier of F attacking Cl₂ is zero and that of Cl attacking F₂ is only 15.57 kJ?mol^?1 at the B3LYP level. However, the calculated dissociation energies of F₂ and Cl₂ are as high as 145.40 and 192.48 kJ?mol^?1, respectively. When the reaction proceeds through a bimolecular reaction mechanism, two four‐center transition states are obtained and the lower energy barrier is 218.69 kJ?mol^?1. Therefore, the title reaction F₂+Cl₂→2ClF is most probably initiated from the atomization of the F₂ molecule and terminated by the reaction of F attacking Cl₂ nearly along the Cl? Cl bond. MP2 calculations lead to the same conclusion, but the geometry of TS and the energy barrier are somewhat different. © 2002 John Wiley & Sons, Inc. Int J Quantum Chem, 2002     

Reaction of OH radical with mono‐, di‐, and trichloroacetaldehyde: an ab initio study

The hydrogen abstraction reactions between chlorine‐substituted acetaldehydes and OH radicals have been investigated by using ab initio molecular orbital theory. Equilibrium geometries and transition‐state structures have been optimized at the (U)MP2/6‐311G(d,p) level. Activation barriers and heats of reaction for different reaction channels have been estimated from the single‐point calculations at the (U)MP2/6‐311G(2df,2p) level. Three, two, and one hydrogen abstraction channel have been found for the mono‐, di‐, and trichloroacetaldehyde, respectively. At a higher temperature region, hydrogen abstraction from the formyl group is found to be the major reaction channel for all the three chloroacetaldehydes. The effect of halogen substitution on reactivity toward hydrogen abstraction has been discussed. © 2001 John Wiley & Sons, Inc. J Comput Chem 22: 1509–1521, 2001     

Theoretical Studies on the Aminolysis of Ester: Effects of the Catalysis and Substituent to the Reaction of Methyl Indole-3-acetate with Ammonia

A comprehensive exploration of the aminolysis mechanism for methyl indole-3-acetate with ammonia is carried out by employing the B3 LYP/6-311++G(d,p), M06-2 X/6-311++G(d,p) and MP2/6-311++G(d,p)//M06-2 X/6-311++G(d,p) levels. Two alterative reaction channels of the concerted and addition/elimination stepwise processes including the uncatalyzed, base-catalyzed reactions are taken into consideration. Subsequently, the substituent effects and solvent effects in methanol are also evaluated at the M06-2 X/6-311++G(d,p) level. The calculated results indicate that the calculated values of M06-2 X level are quite close to those of MP2, the stepwise pathway has more advantages to the concerted one for all of the reaction processes and the catalyst facilitates the proton migration and decreases the energy barriers as well. It is shown that the most preferred mechanism is the based-catalyzed stepwise process, the substituent of NH2 group slightly accelerates all the aminolysis reaction processes, and the solvent effect does not remarkably change the mechanism of the reaction.     

CVD Reactions of TiCl4 with Ammonia： a Quantum Chemical Study

Based on density functional calculations, the mechanism and the energetic course of the chemical vapor deposition （CVD） reaction of TiCl4 with NH3 were studied at the level of B3LYP with 6-311g（d） basis set. Furthermore, the polymerization processes of dimerization, trimerization and tetramerization were investigated. The calculation results indicate that the formation of polymers is favored at the elimination reaction. On the basis of the calculated energetics, a possible mechanism of the reduction reaction has been proposed.     

二氟硅烯与甲醛环加成反应机理的理论研究

The mechanism of the cycloadditohn reaction of singlet difluorosilylene with formaldehyde have been studied by RHF/6-311G* gradient method. The electron correlation energy corrections of energies for all the structures were computed using second-order Moller-Plesset perturbation theory(MP2). The results show that this reaction proceeds via two steps:1)Difluorosilylene and formaldehyde form an intermediate complex, it is an exothermal reaction with no barrier.2) The intermediate complex isomerizes to form the product, after being corrected by zero-point energies, the barrier is 127.28 kJ•mol-1 (MP2/6-311G* 6-311G*).     

A theoretical study on decomposition and rearrangement reaction mechanism of trichloroacetyl chloride (CCl3COCl)

In the present study, the possible decomposition and rearrangement reaction profile of trichloroacetyl chloride have been studied using UMP2/6‐311++G (2d, 2p) level of ab initio and UB3LYP/6‐311++G (2d, 2p) level of density functional theory methods. The harmonic vibrational frequencies were calculated at the same level of theory used for the characterization of stationary points and zero‐point vibrational energy corrections. The potential energy barrier and activation energy between each step of the reaction have been calculated for the seven possible reaction pathways (Ia–c, IIa–b, IIIa–b). The trichloroacetyl chloride is an asymmetric ketone where the two α bonds of acetyl chloride, the C? C and C? Cl bonds are strong with dissociation energy of 72 kcal/mol. The phosgene (COCl₂), dichloroketene (CCl₂CO), carbon dichloride (CCl₂), carbon tetrachloride (CCl₄), and carbon monoxide (CO) are the major dominant products on the decomposition of the trichloroacetyl chloride. These resultant products are more hazards than the parent trichloroacetyl chloride molecules. The positive value of the reaction energy indicate that the overall reaction profile is found to be endothermic at the UMP2 and UBLYP/6‐311++G(2d, 2p) levels of theory, respectively, at UMP2/6‐311G** optimized geometry. © 2010 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

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.     

Transition states for deactivation reactions in the modeled 2,2,6,6‐tetramethyl‐1‐piperidinyloxy‐mediated free‐radical polymerization of acrylonitrile

The geometries and energetics of transition states (TS) for radical deactivation reactions, including competitive combination and disproportionation reactions, have been studied for the modeled 2,2,6,6‐tetramethyl‐1‐piperidinyloxy (TEMPO)‐mediated free‐radical polymerization of acrylonitrile with quantum mechanical calculations at the DFT/UB3‐LYP/6‐311+G(3df,2p)//(U)AM1 level of theory (where DFT is density functional theory, AM1 is Austin model 1, and UAM1 is unrestricted Austin model 1). A method providing reasonable starting geometries for an effective search for TS between the TEMPO radical and 1‐cyanopropyl radical mimicking the growing polyacrylonitrile macroradical is shown. For the hydrogen atom abstraction reaction by the TEMPO radical from the 1‐cyanopropyl radical, practically one TS has been found, whereas for the combination reaction of the radicals, several TS have been found, mainly differing in out‐of‐plane angle α of the N? O bond in the TEMPO structure. α in the TS is correlated with the activation energy, ΔE, determined from the single‐point calculation at the DFT UB3‐LYP/6‐311+G(3df, 2p)//UAM1 level for the combination reaction of CH₃AN· with the TEMPO radical. The theoretical activation energy for the coupling reaction from DFT UB3‐LYP/6‐311+G(3df, 2p)//UAM1 calculations has been estimated to be 11.6 kcal mol^?1, that is, only about 4.5 times smaller than ΔE for the disproportionation reaction obtained with the DFT UB3‐LYP/6‐311+G(3df, 2p)//(U)AM1 approach. © 2005 Wiley Periodicals, Inc. J Polym Sci Part A: Polym Chem 44: 914–927, 2006     

Quantum study on the electrocyclic reactions of CH2CHCHCHX (XH,F, Cl)

The microscopic mechanisms of the electrocyclic reactions for cis‐1,3‐butadiene and its monofluoro‐, monochloroderivatives have been studied by density functional theory (DFT), using the B3LYP method and 6‐311++G** basis sets. We optimized the geometric configurations of reactants, transition states, and products; verified all the probable transition states through vibrational analysis; and calculated the relative single‐point energies at the QCISD(T)/6‐311++G**//B3LYP/6‐311++G**. The results show that the monofluoro‐, monochloroderivatives of cis‐1,3‐butadiene both have two conformers; the reactant favors the electrocyclic reaction when one outboard hydrogen atom of the CH₂ groups is substituted by the fluorine or chlorine atom. © 2004 Wiley Periodicals, Inc. Int J Quantum Chem, 2004     

Deactivation reactions in the modeled 2,2,6,6‐tetramethyl‐1‐piperidinyloxy‐mediated free‐radical polymerization of styrene: A comparative study with the 2,2,6,6‐tetramethyl‐1‐piperidinyloxy/acrylonitrile system

The competitiveness of the combination and disproportionation reactions between a 1‐phenylpropyl radical, standing for a growing polystyryl macroradical, and a 2,2,6,6‐tetramethyl‐1‐piperidinyloxy (TEMPO) radical in the nitroxide‐mediated free‐radical polymerization of styrene was quantitatively evaluated by the study of the transition geometry and the potential energy profiles for the competing reactions with the use of quantum‐mechanical calculations at the density functional theory (DFT) UB3‐LYP/6‐311+G(3df, 2p)//(unrestricted) Austin Model 1 level of theory. The search for transition geometries resulted in six and two transition structures for the radical combination and disproportionation reactions, respectively. The former transition structures, mainly differing in the out‐of‐plane angle of the N? O bond in the transition structure TEMPO molecule, were correlated with the activation energy, which was determined to be in the range of 8.4–19.4 kcal mol^?1 from a single‐point calculation at the DFT UB3‐LYP/6‐311+G(3df, 2p)//unrestricted Austin Model 1 level. The calculated activation energy for the disproportionation reaction was less favorable by a value of more than 30 kcal mol^?1 in comparison with that for the combination reaction. The approximate barrier difference for the TEMPO addition and disproportionation reaction was slightly smaller for the styrene polymerization system than for the acrylonitrile polymerization system, thus indicating that a β‐proton abstraction through a TEMPO radical from the polymer backbone could diminish control over the radical polymerization of styrene with the nitroxide even more than in the latter system. © 2006 Wiley Periodicals, Inc. J Polym Sci Part A: Polym Chem 45: 232–241, 2007     

Theoretical Study on Gas—phase Pyrolytic Reactions of N—Ethyl,Isopropyl and N—t—Butyl Substituted 2—Aminopyrazine

Density functional theory (DFT) and ab initio methods were used to study gas‐phase pyrolytic reaction mechanisms of iV‐ethyl, N‐isopropyl and N‐t‐butyl substituted 2‐aminopyrazine at B3LYP/6–31G* and MP2/6–31G*, respectively. Single‐point energies of all optimized molecular geometries were calculated at B3LYP/6–311 + G(2d,p) level. Results show that the pyrolytic reactions were carried out through a unimolecular first‐order mechanism which were caused by the migration of atom H(17) via a six‐member ring transition state. The activation energies which were verified by vibrational analysis and correlated with zero‐point energies along the reaction channel at B3LYP/6–311 + G(2d,p) level were 252.02 kJ. mo^?1 (N‐ethyl substituted), 235.92 kJ‐mol^?1 (N‐t‐isopropyl substituted) and 234.27 kJ‐mol^?1 (N‐t‐butyl substituted), respectively. The results were in good agreement with available experimental data.     

D-苯丙氨酸与Cu＋(1S0, 3d10)气相反应理论研究

用量子化学密度泛函(DFT)方法研究D-苯丙氨酸与一价基态金属阳离子Cu^＋在气相中反应的机理. 在B3LYP/6-31G*水平上, 优化了反应包含的4个反应通道的反应物、中间体、过渡态和产物的几何构型, 并采用B3LYP/DZVP, B3LYP/[6-311＋G**(C,H,O)＋Lanl2dz(Cu)], B3LYP/6-311＋G**, MP2/6-311＋G** 等方法对各驻点进行了单点能计算. 通过对计算结果的分析, 获得了其单重态反应势能面的一般轮廓、各驻点几何构型优化参数, 明确了其反应机理.     

马来酰亚胺阴离子聚合链引发反应机理研究

以氢氧根离子催化马来酰亚胺聚合反应为例,通过理论计算探讨了在DMF溶剂(ε=36.7)中阴离子聚合链引发阶段的微观机理,在DFT(B3LYP/6-311G*)水平搜索得到反应路径,相同水平上对反应物中间体、过渡态和产物进行几何优化和振动分析,并经过IRC验证.发现在反应过程中氢氧根先进攻酰亚胺的碳,经过一过渡态最终生成氢氧根连接在碳碳双键上的阴离子,证明反应机理为阴离子聚合.由于活化能较小,所以聚合反应在较温和的条件下即可进行.     

Direct dynamic study on the hydrogen abstraction reaction of H2CO with NCO

A direct ab initio dynamics method is used to investigate the hydrogen‐abstraction reaction of H₂CO with NCO. The potential energy surface information is obtained at the MP2/6‐311G(d,p) level. More accurate single‐point energy is refined at the G3(MP2)//MP2/6‐311G(d,p) level. Furthermore, the rate constants of reaction H₂CO + NCO are evaluated by using the canonical variational transition state theory with small‐curvature tunneling contributions over a wide temperature range of 200–2000 K. The calculated reaction enthalpy and rate constants are in good agreement with the available experimental values. © 2009 Wiley Periodicals, Inc. Int J Chem Kinet 41: 394–400, 2009     

A Theoretical Study of the Kinetics of the Hydrogen Atom Abstraction Reactions from Cyclopropane by H,O (3P), and Cl (2P3/2) Atoms and OH Radicals

The rate constants of the H‐abstraction reactions from cyclopropane by H, O (³P), Cl (²P_3/2), and OH radicals have been calculated over the temperature range of 250?2500 K using two different levels of theory. Calculations of optimized geometrical parameters and vibrational frequencies are performed using the MP2 method combined with the cc‐pVTZ basis set and the 6–311++G(d,p) basis set. Single‐point energy calculations have been carried out with the highly correlated ab initio coupled cluster method in the space of single, double, and triple (perturbatively) electron excitations CCSD(T) using either the cc‐pVTZ, aug‐cc‐pVTZ, and aug‐cc‐pVQZ basis sets or the 6–311++G(3df,3pd) basis set. The CCSD(T) calculated potential energies have been extrapolated to the complete basis limit (CBS) limit. The Full Configuration Interaction (FCI) energies have been also estimated using the continued‐fraction approximation as proposed by Goodson (J. Chem. Phys., 2002, 116, 6948–6956). Canonical transition‐state theory combined with an Eckart tunneling correction has been used to predict the rate constants as a function of temperature using two kinetic models (direct abstraction or complex mechanism) at two levels of theory (CCSD(T)‐cf/CBS//MP2/cc‐pVTZ and CCSD(T)‐cf/6–311++G(3df,3pd)//MP2/6–311++G(d,p)). The calculated kinetic parameters are in reasonable agreement with their literature counterparts for all reactions. In the light of these trends, the use of the Pople‐style basis sets for studying the reactivity of other systems such as larger cycloalkanes or halogenated cycloalkanes is recommended because the 6–311++G(3df,3pd) basis set is less time consuming than the aug‐cc‐pVQZ basis set. Based on our calculations performed at the CCSD(T)‐cf/CBS//MP2/cc‐pVTZ level of theory, the standard enthalpy of formation at 298 K for the cyclopropyl radical has been reassessed and its value is (290.5 ± 1.6) kJ mol^?1.