Decomposition reaction rate of BCl3-C3H6(propene)-H2 in the gas phase

The decomposition reaction rate in the BCl(3)-C(3)H(6)-H(2) gas phase reaction system in preparing boron carbides was investigated based on the most favorable reaction pathways proposed by Jiang et al. [Theor. Chem. Accs. 2010, 127, 519] and Yang et al. [J. Theor. Comput. Chem. 2012, 11, 53]. The rate constants of all the elementary reactions were evaluated with the variational transition state theory. The vibrational frequencies for the stationary points as well as the selected points along the minimum energy paths (MEPs) were calculated with density functional theory at the B3PW91/6-311G(d,p) level and the energies were refined with the accurate model chemistry method G3(MP2). For the elementary reaction associated with a transition state, the MEP was obtained with the intrinsic reaction coordinates, while for the elementary reaction without transition state, the relaxed potential energy surface scan was employed to obtain the MEP. The rate constants were calculated for temperatures within 200-2000 K and fitted into three-parameter Arrhenius expressions. The reaction rates were investigated by using the COMSOL software to solve numerically the coupled differential rate equations. The results show that the reactions are, consistent with the experiments, appropriate at 1100-1500 K with the reaction time of 30 s for 1100 K, 1.5 s for 1200 K, 0.12 s for 1300 K, 0.011 s for 1400 K, or 0.001 s for 1500 K, for propene being almost completely consumed. The completely dissociated species, boron carbides C(3)B, C(2)B, and CB, have very low concentrations, and C(3)B is the main product at higher temperatures, while C(2)B is the main product at lower temperatures. For the reaction time 1 s, all these concentrations approach into a nearly constant. The maximum value (in mol/m(3)) is for the highest temperature 1500 K with the orders of -13, -17, and -23 for C(3)B, C(2)B, and CB, respectively. It was also found that the logarithm of the overall reaction rate and reciprocal temperature have an excellent linear relationship within 700-2000 K with a correlation coefficient of 0.99996. This corresponds to an apparent activation energy 337.0 kJ/mol, which is comparable with the energy barrier 362.6 kJ/mol of the rate control reaction at 0 K but is higher than either of the experiments 208.7 kJ/mol or the Gibbs free energy barrier 226.2 kJ/mol at 1200 K.     

H3PO→H2POH异构化反应的直接动力学研究

在QCISD(T)/6-311C++G(2df,2pd)//QCISD/6-311C++G(d,p)+ZPE水平上,对H3PO的异构化反应H3PO→(1)H2POH(trans)→(2)H2POH(cis)进行了计算研究.结果表明,H原子由P原子向O原子迁移反应(1)的能垒为250.0kJ/mol,是反应速率控制步骤,而O_H键绕P_O键旋转的构型转化反应(2)的能垒只为12.3kJ/mol.利用经典过渡态理论(TST)与变分过渡态理论(CVT)分别计算了反应(1)在200～2000K温度区间内的速率常数kTST和kCVT,获得了经小曲率隧道效应(SCT)及Eckart模型校正后的速率常数kTST/Eckart和kCVT/SCT.对只涉及H原子迁移的反应(1),量子力学隧道效应的影响在低温段非常明显,而变分效应对反应速率常数的影响很小.     

A computational study of conformational interconversions in 1,4-dithiacyclohexane (1,4-dithiane)

Ab initio molecular orbital theory with the 6-31G(d), 6-31G(d,p), 6-31+G(d), 6-31+G(d,p), 6-31+G(2d,p), 6-311G(d), 6-311G(d,p), and 6-311+G(2d,p) basis sets and density functional theory (BLYP, B3LYP, B3P86, B3PW91) have been used to locate transition states involved in the conformational interconversions of 1,4-dithiacyclohexane (1,4-dithiane) and to calculate the geometry optimized structures, relative energies, enthalpies, entropies, and free energies of the chair and twist conformers. In the chair and 1,4-twist conformers the C-Hax and C-Heq bond lengths are equal at each carbon, which suggest an absence of stereoelectronic hyperconjugative interactions involving carbon-hydrogen bonds. The 1,4-boat transition state structure was 9.53 to 10.5 kcal/mol higher in energy than the chair conformer and 4.75 to 5.82 kcal/mol higher in energy than the 1,4-twist conformer. Intrinsic reaction coordinate (IRC) calculations showed that the 1,4-boat transition state structure was the energy maximum in the interconversion of the enantiomers of the 1,4-twist conformer. The energy difference between the chair conformer and the 1,4-twist conformer was 4.85 kcal/mol and the chair-1,4-twist free energy difference (deltaG degrees (c-t)) was 4.93 kcal/mol at 298.15 K. Intrinsic reaction coordinate (IRC) calculations connected the transition state between the chair conformer and the 1,4-twist conformer. This transition state is 11.7 kcal/mol higher in energy than the chair conformer. The effects of basis sets on the 1,4-dithiane calculations and the relative energies of saturated and unsaturated six-membered dithianes and dioxanes are also discussed.     

HNCO+HCO→NCO+CH2O氢转移反应的从头算及动力学研究

在UMP2(Full)/6-311G(d,p)计算水平上,优化了标题反应的反应物、过渡态、产物的几何结构,沿最小能量途径讨论了异氰酸(HNCO)和甲酰自由基(HCO)发生氢转移反应位能面上驻点的结构以及相互作用分子结构变化.指出该反应是一个N-H键断裂和C-H键生成的协同反应.进一步采用UQCISD(T,Full)方法对反应途径上的驻点进行了单点能量校正,得出该反应的计算位垒是91.47 kJ/mol,与实验值108.92 kJ/mol接近在500～2500K实验温度范围内,运用变分过渡态理论(CVT)计算得到的速率常数与实验观测值进行了比较     

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更稳定. 另外, 速率常数计算结果显示, 量子力学隧道效应在低温段对速率常数的计算有显著影响, 而变分效应在计算温度段内对速率常数的影响可以忽略.     

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键断裂反应可以忽略.     

Reaction pathways of propene pyrolysis

The gas‐phase reaction pathways in preparing pyrolytic carbon with propene pyrolysis have been investigated in detail with a total number of 110 transition states and 50 intermediates. The structure of the species was determined with density functional theory at B3PW91/6‐311G(d,p) level. The transition states and their linked intermediates were confirmed with frequency and the intrinsic reaction coordinates analyses. The elementary reactions were explored in the pathways of both direct and the radical attacking decompositions. The energy barriers and the reaction energies were determined with accurate model chemistry method at G3(MP2) level after an examination of the nondynamic electronic correlations. The heat capacities and entropies were obtained with statistical thermodynamics. The Gibbs free energies at 298.15 K for all the reaction steps were reported. Those at any temperature can be developed with classical thermodynamics by using the fitted (as a function of temperature) heat capacities. It was found that the most favorable paths are mainly in the radical attacking chain reactions. The chain was proposed with 26 reaction steps including two steps of the initialization of the chain to produce H and CH₃ radicals. For a typical temperature (1200 K) adopted in the experiments, the highest energy barriers were found in the production of C₃ to be 203.4 and 193.7 kJ/mol. The highest energy barriers for the production of C₂ and C were found 174.1 and 181.4 kJ/mol, respectively. These results are comparable with the most recent experimental observation of the apparent activation energy 201.9 ± 0.6 or 137 ± 25 kJ/mol. © 2010 Wiley Periodicals, Inc. J Comput Chem, 2010     

Ab initio study of the reaction path of the reaction HNCO+CN

HNCO is a convenient photolytic source of NCO and NH radicals for laboratory kinetics studies of elementary reaction[1] and plays an important role in the combustion and atmosphere chemistry. It can re- move deleterious compounds rapidly from exhausted ga…     

An investigation of the lowest reaction pathway of propene + BCl3 decomposition in chemical vapor deposition process

The lowest reaction pathway or one of the favored possible paths in the CVD process of preparing boron carbides with BCl₃-C₃H₆(propene)-H₂ precursors was searched theoretically, which involves 95 transition states and 103 intermediates. The geometries of the species were optimized by employing the B3PW91/6-311G(d,p) method. The transition states as well as their linked intermediates were confirmed with frequency and the intrinsic reaction coordinates analyses. The energy barriers and the reaction energies were evaluated with the accurate model chemistry method at G3(MP2) level after a non-dynamical electronic correlation detection. The heat capacities and entropies were obtained with statistical thermodynamics, and the heat capacities were fitted into analytical equations. The Gibbs free energies at 298.15 K for all of the reaction steps were reported. The energies at any temperature could be derived classically by using the analytical heat capacities. All the possible elementary reactions, including both direct decomposition and the radical attacking dissociations, for each reaction step were examined, and the one with the lowest energy or energy barrier was further studied in the next step. It was found that there are 19 reaction steps in the lowest path to produce the final BC₃ cluster including two steps of initializing the reaction chain of producing H and Cl radicals. The highest energy in the lowest reaction pathway is 215.1 kJ/mol at 298.15 K and that for 1,200 K is 275.1 kJ/mol. The results are comparable with the most recent experimental observation of the apparent activation energy 208.7 kJ/mol.     

3,6-二羟基哒嗪质子转移过程的理论研究

为了探索3,6-二羟基哒嗪分子醇式和酮式结构互变异构化的反应机理,本研究组采用DFT B3LYP/6-311+G(d)方法对标题化合物异构化反应的势能面进行了研究,在探讨各种可能的反应途径中,发现至少有34种异构体和43种过渡态.结果表明,6-羟基-3(2H)-哒嗪酮不论是单体,与水形成配合物,还是二聚体,比其相对应的异构体能量低,表明在通常情况下是以6-羟基-3(2H)-哒嗪酮形式稳定存在的,这与前人通过实验数据对互变异构体的比率进行预测的结果是一致的;在考察的可能反应途径中,直接进行的分子内质子转移过程需要的活化能为142.2 kJ·mol^-1,水助催化时,反应活化能为41.3 kJ·mol^-1,考虑溶剂效应后,其活化能为59.2 kJ·mol^-1,二聚体双质子转移的活化能为16.8 kJ·mol^-1,二聚体双质子转移所需活化能最低,在室温下就可以进行.由此可见氢键在降低反应活化能方面起着重要的作用.     

黄嘌呤异构体质子转移异构化反应机理的理论研究

采用量子化学密度泛函B3LYP/6-31G(d,p)和M06-2X/6-311++G(d,p)方法对黄嘌呤两种酮式异构体X(1,3,7)与X(1,3,9)间质子转移引起的互变异构反应机理进行了计算研究,获得了异构化反应过程的反应焓﹑活化吉布斯自由能和质子转移反应的速率常数等参数。水相计算采用极化连续(PCM)模型。结果表明,由于可能的氢迁移顺序差异,分子内由X(1,3,7)向X(1,3,9)异构化可能共有16条反应通道,涉及11个中间体和20个过渡态,其主反应通道速控步骤的活化吉布斯自由能为183.10k J/mol,速率常数为5.17×10-20s-1,其余各通道速控步骤活化吉布斯自由能均较高,而且整体水溶剂效应不利于质子转移的发生。     

Conformers and intramolecular hydrogen bonding of the oxalic acid monomer and its anions

Density functional theory, B3LYP/6‐31G** and B3LYP/6‐311+G(2d,p), and ab initio MP2/6‐31G** calculations have been carried out to investigate the conformers, transition states, and energy barriers of the conformational processes of oxalic acid and its anions. QCISD/6‐31G** geometrical optimization is also performed in the stable forms. Its calculated energy differences between the two most stable conformers are very near to the related observed value at 7.0 kJ/mol. It is found that the structures and relative energies of oxalic acid conformers predicted by these methods show similar results, and that the conformer L1 (C_2h) with the double‐interfunctional‐groups hydrogen bonds is the most stable conformer. The magnitude of hydrogen bond energies depends on the energy differences of various optimized structures. The hydrogen bond energies will be about 32 kJ/mol for interfunctional groups, 17 kJ/mol for weak interfunctional groups, 24 kJ/mol for intra‐COOH in (COOH)₂, and 60 kJ/mol for interfunctional groups in (COOH)COO⁻¹ ion if calculated using the B3LYP/6‐311+G(2d,p) method. © 2000 John Wiley & Sons, Inc. Int J Quant Chem 76: 541–551, 2000     

Can cytosine, thymine and uracil be formed in interstellar regions? A theoretical study

This theoretical study investigates possible synthetic routes to cytosine, uracil and thymine in the gas phase from precursor molecules that have been detected in interstellar media. Studies at the CCSD(T)/6-311++G(d,p)//B3LYP/6-311++G(d,p) level of theory suggest that: The reactions between :CCCNH and :CCCO with monosolvated urea may constitute viable interstellar syntheses of cytosine and uracil. No low energy equilibration between cytosine and uracil has been demonstrated. The interaction of :CH(2) with the 5 C-H bond of uracil may form thymine in an energetically favourable reaction, but competing reactions where :CH(2) reacts with double bonds and other CH and NH bonds of uracil, reduce the effectiveness of this synthesis. The reaction between the hydrated propional enolate anion and isocyanic acid may produce thymine, in a reaction sequence where ΔG(reaction)(298 K) is -22 kJ mol(-1) and the maximum energy requirement (barrier to the first transition state) is only 47 kJ mol(-1).     

Direct dynamic study on the hydrogen abstraction reaction CH3CN + OH-->CH2CN + H2O

The reaction of acetonitrile with hydroxyl has been studied using the direct ab initio dynamics methods. The geometries, vibrational frequencies of the stationary points, as well as the minimum energy paths were computed at the BHandHLYP and MP2 levels of theory with the 6-311G(d, p) basis set. The energies were further refined at the PMP4/6-311+G(2df, 2pd) and QCISD(T)/6-311+G(2df, 2pd) levels of theory based on the structures optimized at BHandHLYP/6-311G(d, p) and MP2/6-311G(d, p) levels of theory. The Polyrate 8.2 program was employed to predict the thermal rate constants using the canonical variational transition state theory incorporating a small-curvature tunneling correction. The computed rate constants are in good agreement with the available experimental data.     

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.     

A DFT study on the double bond migration of butene catalyzed by ionic pair of 1-ethyl-3-methyl-imidazolium fluoride

The double bond migration of butene catalyzed by 1-ethyl-3-methyl-imidazolium fluoride (EmimF) has been studied using quantum chemical method. The geometries of reactant, transition state and product for the isomerization have been optimized by density functional theory (DFT) at the B3PW91/6-31G(d,p), 6-311++G(d,p) and aug-cc-PVDZ levels. The computed results show that the 4-H atom on imidazole ring of EmimF has a good catalytic activity to the double bond migration of butene and the catalytic reaction of 1-butene to 2-butene is a synergetic and elementary process. The apparent activation energy of isomerization is about 197 kJ/mol.     

Transition state structure and energetics of the N(2)O + X (X = Cl,Br) reactions

The structural and vibrational properties of the transition state of the N(2)O + X (X = Cl,Br) reactions have been characterized by ab initio methods using density functional theory. We have employed Becke's hybrid functional (B3LYP), and transition state optimizations were performed with 6-31G(d), 6-311G(2d,2p), 6-311+G(3d,2p), and 6-311+G(3df,2p) basis sets. For the chlorine atom reaction the coupled-cluster method (CCSD(T)) with 6-31G(d) basis set was also used. All calculations resulted in transition state structures with a planar cis arrangement of atoms for both reactions. The geometrical parameters of transition states at B3LYP are very similar, and the reaction coordinates involve mainly the breaking of the N-O bond. At CCSD(T)/6-31G(d) level a contribution of the O-Cl forming bond is also observed in the reaction coordinate. In addition, several highly accurate ab initio composite methods of Gaussian-n (G1, G2, G3), their variations (G2(MP2), G3//B3LYP), and complete basis set (CBS-Q, CBS-Q//B3LYP) series of models were applied to compute reaction energetics. All model chemistries predict exothermic reactions. The G3 and G2 methods result in the smallest deviations from experiment, 1.8 and 0 kcal mol(-1), for the enthalpies of reaction for N(2)O reaction with chlorine and bromine, respectively. The G3//B3LYP and G1 methods perform best among the composite methods in predicting energies of the transition state, with a deviation of 1.9 and 3.0 kcal mol(-1), respectively, in the activation energies for the above processes. However, the B3LYP/6-311+G(3df,2p) method gives smaller deviations of 0.4 and -1.0 kcal mol(-1), respectively. The performance of the methodologies applied in predicting transition state energies was analyzed.     

Theoretical study and rate constant calculation for reaction of CF(3)CH(2)OH with OH

The reaction mechanism of CF(3)CH(2)OH with OH is investigated theoretically and the rate constants are calculated by direct dynamics method. The potential energy surface (PES) information, which is necessary for dynamics calculation, is obtained at the B3LYP/6-311G (d, p) level. The single-point energy calculations are performed at the MC-QCISD level using the B3LYP geometries. Complexes, with the energies being less than corresponding reactants and products, are found at the entrance and exit channels for methylene-H-abstraction channel, while for the hydroxyl-H-abstraction channel only entrance complex is located. By means of isodesmic reactions, the enthalpies of the formation for the species CF(3)CH(2)OH, CF(3)CHOH, and CF(3)CH(2)O are estimated at the MC-QCISD//B3LYP/6-311G (d, p) level of theory. The rate constants for two kinds of H-abstraction channels are evaluated by canonical variational transition state theory with the small-curvature tunneling correction (CVT/SCT) over a wide range of temperature 200-2000 K. The calculated results are in good agreement with the experimental values in the temperature region 250-430 K. The present results indicate that the two channels are competitive. Below 289 K, hydroxyl-H-abstraction channel has more contribution to the total rate constants than methylene-H-abstraction channel, while above 289 K, methylene-H-abstraction channel becomes more important and then becomes the major reaction channel.     

TDDFT investigation on nucleic acid bases: comparison with experiments and standard approach

A comprehensive theoretical study of electronic transitions of canonical nucleic acid bases, namely guanine, adenine, cytosine, uracil, and thymine, was performed. Ground state geometries were optimized at the MP2/6-311G(d,p) level. The nature of respective potential energy surfaces was determined using the harmonic vibrational frequency analysis. The MP2 optimized geometries were used to compute electronic vertical singlet transition energies at the time-dependent density functional theory (TDDFT) level using the B3LYP functional. The 6-311++G(d,p), 6-311(2+,2+)G(d,p), 6-311(3+,3+)G(df,pd), and 6-311(5+,5+)G(df,pd) basis sets were used for the transition energy calculations. Computed transition energies were found in good agreement with the corresponding experimental data. However, in higher transitions, the Rydberg contaminations were also obtained. The existence of pisigma* type Rydberg transition was found near the lowest singlet pipi* state of all bases, which may be responsible for the ultrafast deactivation process in nucleic acid bases.     

污泥热解中HCN与CaO的反应机理:密度泛函理论研究

采用密度泛函理论对污泥热解中CaO与HCN在低温段的反应进行了研究。在B3LYP/6-311++(3df,2p)水平上计算得到了反应路径上各驻点的几何构型与频率,并在此构型上使用CCSD(T)/cc-pVQZ进行单点能计算。结果表明,两个HCN分子吸附于CaO后,质子发生转移时出现反应路径中最大能垒(310.33 kJ/mol)。使用经典过渡态理论拟合了反应中各步骤的阿累尼乌斯公式,计算了三种典型温度下各步骤的反应速率,发现质子转移为该反应的决速步骤,且温度越高CaO对HCN的作用效果越好。