A simple model based on the statistics of individual atoms [Europhys. Lett. 94 40002 (2011)] or molecules [Chin. Phys. Lett. 29 080504 (2012)] was used to predict chemical reaction rates without empirical paramete...A simple model based on the statistics of individual atoms [Europhys. Lett. 94 40002 (2011)] or molecules [Chin. Phys. Lett. 29 080504 (2012)] was used to predict chemical reaction rates without empirical parameters, and its physical basis was further investigated both theoretically and via MD simulations. The model was successfully applied to some reactions of extensive experimental data, showing that the model is significantly better than the conventional transition state theory. It is worth noting that the prediction of the model on ab initio level is much easier than the transition state theory or unimolecular RRKM theory.展开更多
The reaction of disilane with atomic hydrogen has been studied. This reaction involves both substitution and abstraction. Calculations show that the hydrogen abstraction is the strongest competing channel. The canonic...The reaction of disilane with atomic hydrogen has been studied. This reaction involves both substitution and abstraction. Calculations show that the hydrogen abstraction is the strongest competing channel. The canonical variational transition state theory with a small curvature tunneling correction (SCT) has been used for the kinetic calculation. The theoretical results are in good agreement with the available experimental data. Comparing the reactions of atomic hydrogen with disilane and silane, it can be seen that the reactivity of the Si-H bond is higher in Si2H6than that in SiH4.展开更多
The conventional symmetry numbers σ_≠~' of activated complexes may lead to error in the rate constant expression of transition state theory, whereas the statistical factor ι~≠ or ι may violate the principle o...The conventional symmetry numbers σ_≠~' of activated complexes may lead to error in the rate constant expression of transition state theory, whereas the statistical factor ι~≠ or ι may violate the principle of detailed balance. A mathematically precise definition of the symmetry number σ_≠ of activated complex is given, i.e. σ_≠=_iN_4(?)/m, m is the number of physically distinct configurations of labelled transition state and N_i is the identical atoms in the activated complex. The identical atoms must belong to the same molecule of reactants and products. The present symmetry numbers σ_≠ of activated complexes assure not only obtaining correct rate constant expressions but also obeying the principle of detailed balance. It can be used with the statistical factor ι to construct the structures of transition states for unimolecular and bimolecular exchange reactions.展开更多
The reaction of C3H8+O(^3p)→C3HT+OH is investigated using ab initio calculation and dynamical methods. Electronic structure calculations for all stationary points are obtained using a dual-level strategy. The geo...The reaction of C3H8+O(^3p)→C3HT+OH is investigated using ab initio calculation and dynamical methods. Electronic structure calculations for all stationary points are obtained using a dual-level strategy. The geometry optimization is performed using the unrestricted second-order Moller-Plesset perturbation method and the single-point energy is computed us- ing the coupled-cluster singles and doubles augmented by a perturbative treatment of triple excitations method. Results indicate that the main reaction channel is C3Hs+O(^3p)→i- C3HT+OH. Based upon the ab initio data, thermal rate constants are calculated using the variational transition state theory method with the temperature ranging from 298 K to 1000 K. These calculated rate constants are in better agreement with experiments than those reported in previous theoretical studies, and the branching ratios of the reaction are also calculated in the present work. Furthermore, the isotope effects of the title reaction are calculated and discussed. The present work reveals the reaction mechanism of hydrogenabstraction from propane involving reaction channel competitions is helpful for the understanding of propane combustion.展开更多
The reaction mechanism of SiCl4 with H2 has been studied theoretically using Gaussian 98 program at B3LYP/6-311G^* level. Three different reaction paths (a, b, c) in the gas phase were obtained. The geometries, vib...The reaction mechanism of SiCl4 with H2 has been studied theoretically using Gaussian 98 program at B3LYP/6-311G^* level. Three different reaction paths (a, b, c) in the gas phase were obtained. The geometries, vibrational frequencies and energies of every stagnation point in the reaction channel were calculated and the mechanisms have been confirmed. The results show that path a has an activation energy of 79.12 kcal/mol, which was considered as the main reaction path. Comparably, paths b and c have the energy barriers of 125.07 and 136.25 kcal/mol, respectively. The reaction rate constant was calculated by TST method over a wide temperature range of 900~1600 K, which further confirmed that path a was the main reaction channel展开更多
The interconversion between the two distinct isomers of methyl vinyl ether (MVE), the formation of the primary ozonides from O3-initated reactions of MVE, the transformation between the primary ozonides, and the sub...The interconversion between the two distinct isomers of methyl vinyl ether (MVE), the formation of the primary ozonides from O3-initated reactions of MVE, the transformation between the primary ozonides, and the subsequent fragmentation were studied using quantum chemical methods at the BHandHLYP/6311++G(d,p) level of theory for optimized geometries and frequency calculations and at the QCISD/631G(d,p) level for the single point energy calculations. The rate coefficients were calculated for the temperature range 280-440 K by using the canonical transition state theory (TST). For ozone addition to MVE, there are two different possibilities discussed on the basis of two different possible orientations for ozone attack. The results of the theoretical study indicate that although the synperiplanar-MVE is 7.11 kJ/mol more stable than the antiperiplanar-MVE, the antiperiplanar-MVE plays a more important role in formation of the primary ozonides because the primary ozonides formed from the ozone addition antiperiplanar-MVE are more stable and the energy barriers corresponding to transition states are lower. The intereonversion between the primary ozonides formed from the ozone addition to antiperiplanar-MVE is the most accessible compared with the transformations between other primary ozonides. The cleavage of the primary ozonides mainly leads to the formation of the CH2OO, which is in agreement with the experimental estimates. The calculated overall rate constant for the ozone-initiated reactions is 4.8× 10^-17 cm^3/(molecule.s) at 298.15 K, which agrees with the experimental value for ethyl vinyl ether.展开更多
The abstraction reaction of H with (CH_3)_3SiH was investigated at the high levels of ab initio molecule orbital theory. The geometries were optimized at the MP2 level with 6-31G( d ) basis set, and G2MP2 level was us...The abstraction reaction of H with (CH_3)_3SiH was investigated at the high levels of ab initio molecule orbital theory. The geometries were optimized at the MP2 level with 6-31G( d ) basis set, and G2MP2 level was used for the final energy calculations. The theoretical analysis provides the conclusive evidence that the main process is the hydrogen abstraction from the Si-H bond, leading to the formation of H_2 and silyl radicals; the hydrogen abstraction from the C-H bond has a higher barrier and is difficult to react. The kinetics was calculated with canonical variational transition-state theory (CVT) over the temperature range 200-1 000 K, and the theoretical rate constants match well with the later experimental values.展开更多
Hydrogen abstraction reaction, H+C2H4 --H2+C2H2 was studied by using A initio SCF method. Ge-ometries were fully optimized at SCF level and energies were computed at STO-3G basis set for reactants and transition state...Hydrogen abstraction reaction, H+C2H4 --H2+C2H2 was studied by using A initio SCF method. Ge-ometries were fully optimized at SCF level and energies were computed at STO-3G basis set for reactants and transition state. Vibrational analysis was performed thereupon. Finally, the rate constant calculations were carried out at different temperatures for all range of reaction temperature according to Eyring's sbwlute reaction rate theory. The calculated activation energy is 12. 68 kcal/mol, lower than observed value (H. S kcal/mol) by 1. 82 kcal/mol only. The agreement of the calculated rate constants with the experiments is satisfactory.展开更多
UMP2 method was performed to study the reaction mechanism of CH2 with HNCO.The geometrical configurations of reactants,intermediates,transition states and products were optimized by UMP2 method at the 6-311++G(d,p)l...UMP2 method was performed to study the reaction mechanism of CH2 with HNCO.The geometrical configurations of reactants,intermediates,transition states and products were optimized by UMP2 method at the 6-311++G(d,p)level.The energies of stationary points along the pathways were calculated at QCISD(T)/6-311++G(d,p).Intermediates and transition states were confirmed by the results of vibration analysis and IRC calculation.From the results of the reaction mechanism of CH2 with HNCO,it is shown that in this reaction there exist three pathways,which are CH2(S)+HNCO→IM1→TS1→CH3+NCO;CH2(S)+HNCO→IM2→TS2→CH2CONH,CH2(T)+HNCO→IM3→TS3→CH2NH+CO,and CH2(T)+HNCO→IM3→TS3→CH2NH+CO.展开更多
The potential energy surfaces (PES) of unimolecular dissociation reactions for di-ethyl beryl- lium and di-t-butyl beryllium are investigated by B3LYP, CCSD(T), and G3B3 approaches. Possible reaction pathways thro...The potential energy surfaces (PES) of unimolecular dissociation reactions for di-ethyl beryl- lium and di-t-butyl beryllium are investigated by B3LYP, CCSD(T), and G3B3 approaches. Possible reaction pathways through either the radical or transition state (TS) of the molecules are considered. The geometries, vibrational frequencies and relative energies for various sta- tionary points are determined. From the study of energetics, the TS pathways arising from concerted molecular eliminations are indicated to be the main dissociation pathways for both molecules. The PES differences of the dissociation reactions are investigated. The activation energies and rate constants will be helpful for investigating the predictive ability of the reaction in further theoretical and experimental research.展开更多
基金Project supported by the National Natural Science Foundation of China(Grant No.11274073)the Leading Academic Discipline Project of Shanghai,China(Grant No.B107)
文摘A simple model based on the statistics of individual atoms [Europhys. Lett. 94 40002 (2011)] or molecules [Chin. Phys. Lett. 29 080504 (2012)] was used to predict chemical reaction rates without empirical parameters, and its physical basis was further investigated both theoretically and via MD simulations. The model was successfully applied to some reactions of extensive experimental data, showing that the model is significantly better than the conventional transition state theory. It is worth noting that the prediction of the model on ab initio level is much easier than the transition state theory or unimolecular RRKM theory.
基金the Research Fund for the Doctoral Program of High Education of China (Grant No. 1999042201).
文摘The reaction of disilane with atomic hydrogen has been studied. This reaction involves both substitution and abstraction. Calculations show that the hydrogen abstraction is the strongest competing channel. The canonical variational transition state theory with a small curvature tunneling correction (SCT) has been used for the kinetic calculation. The theoretical results are in good agreement with the available experimental data. Comparing the reactions of atomic hydrogen with disilane and silane, it can be seen that the reactivity of the Si-H bond is higher in Si2H6than that in SiH4.
文摘The conventional symmetry numbers σ_≠~' of activated complexes may lead to error in the rate constant expression of transition state theory, whereas the statistical factor ι~≠ or ι may violate the principle of detailed balance. A mathematically precise definition of the symmetry number σ_≠ of activated complex is given, i.e. σ_≠=_iN_4(?)/m, m is the number of physically distinct configurations of labelled transition state and N_i is the identical atoms in the activated complex. The identical atoms must belong to the same molecule of reactants and products. The present symmetry numbers σ_≠ of activated complexes assure not only obtaining correct rate constant expressions but also obeying the principle of detailed balance. It can be used with the statistical factor ι to construct the structures of transition states for unimolecular and bimolecular exchange reactions.
文摘The reaction of C3H8+O(^3p)→C3HT+OH is investigated using ab initio calculation and dynamical methods. Electronic structure calculations for all stationary points are obtained using a dual-level strategy. The geometry optimization is performed using the unrestricted second-order Moller-Plesset perturbation method and the single-point energy is computed us- ing the coupled-cluster singles and doubles augmented by a perturbative treatment of triple excitations method. Results indicate that the main reaction channel is C3Hs+O(^3p)→i- C3HT+OH. Based upon the ab initio data, thermal rate constants are calculated using the variational transition state theory method with the temperature ranging from 298 K to 1000 K. These calculated rate constants are in better agreement with experiments than those reported in previous theoretical studies, and the branching ratios of the reaction are also calculated in the present work. Furthermore, the isotope effects of the title reaction are calculated and discussed. The present work reveals the reaction mechanism of hydrogenabstraction from propane involving reaction channel competitions is helpful for the understanding of propane combustion.
基金This work was supported by the Foundation of Education Committee of Liaoning Province (No.990321076)
文摘The reaction mechanism of SiCl4 with H2 has been studied theoretically using Gaussian 98 program at B3LYP/6-311G^* level. Three different reaction paths (a, b, c) in the gas phase were obtained. The geometries, vibrational frequencies and energies of every stagnation point in the reaction channel were calculated and the mechanisms have been confirmed. The results show that path a has an activation energy of 79.12 kcal/mol, which was considered as the main reaction path. Comparably, paths b and c have the energy barriers of 125.07 and 136.25 kcal/mol, respectively. The reaction rate constant was calculated by TST method over a wide temperature range of 900~1600 K, which further confirmed that path a was the main reaction channel
文摘The interconversion between the two distinct isomers of methyl vinyl ether (MVE), the formation of the primary ozonides from O3-initated reactions of MVE, the transformation between the primary ozonides, and the subsequent fragmentation were studied using quantum chemical methods at the BHandHLYP/6311++G(d,p) level of theory for optimized geometries and frequency calculations and at the QCISD/631G(d,p) level for the single point energy calculations. The rate coefficients were calculated for the temperature range 280-440 K by using the canonical transition state theory (TST). For ozone addition to MVE, there are two different possibilities discussed on the basis of two different possible orientations for ozone attack. The results of the theoretical study indicate that although the synperiplanar-MVE is 7.11 kJ/mol more stable than the antiperiplanar-MVE, the antiperiplanar-MVE plays a more important role in formation of the primary ozonides because the primary ozonides formed from the ozone addition antiperiplanar-MVE are more stable and the energy barriers corresponding to transition states are lower. The intereonversion between the primary ozonides formed from the ozone addition to antiperiplanar-MVE is the most accessible compared with the transformations between other primary ozonides. The cleavage of the primary ozonides mainly leads to the formation of the CH2OO, which is in agreement with the experimental estimates. The calculated overall rate constant for the ozone-initiated reactions is 4.8× 10^-17 cm^3/(molecule.s) at 298.15 K, which agrees with the experimental value for ethyl vinyl ether.
基金Supported by the Research Foundation for the Doctoral Program of Higher Education of China
文摘The abstraction reaction of H with (CH_3)_3SiH was investigated at the high levels of ab initio molecule orbital theory. The geometries were optimized at the MP2 level with 6-31G( d ) basis set, and G2MP2 level was used for the final energy calculations. The theoretical analysis provides the conclusive evidence that the main process is the hydrogen abstraction from the Si-H bond, leading to the formation of H_2 and silyl radicals; the hydrogen abstraction from the C-H bond has a higher barrier and is difficult to react. The kinetics was calculated with canonical variational transition-state theory (CVT) over the temperature range 200-1 000 K, and the theoretical rate constants match well with the later experimental values.
文摘Hydrogen abstraction reaction, H+C2H4 --H2+C2H2 was studied by using A initio SCF method. Ge-ometries were fully optimized at SCF level and energies were computed at STO-3G basis set for reactants and transition state. Vibrational analysis was performed thereupon. Finally, the rate constant calculations were carried out at different temperatures for all range of reaction temperature according to Eyring's sbwlute reaction rate theory. The calculated activation energy is 12. 68 kcal/mol, lower than observed value (H. S kcal/mol) by 1. 82 kcal/mol only. The agreement of the calculated rate constants with the experiments is satisfactory.
文摘UMP2 method was performed to study the reaction mechanism of CH2 with HNCO.The geometrical configurations of reactants,intermediates,transition states and products were optimized by UMP2 method at the 6-311++G(d,p)level.The energies of stationary points along the pathways were calculated at QCISD(T)/6-311++G(d,p).Intermediates and transition states were confirmed by the results of vibration analysis and IRC calculation.From the results of the reaction mechanism of CH2 with HNCO,it is shown that in this reaction there exist three pathways,which are CH2(S)+HNCO→IM1→TS1→CH3+NCO;CH2(S)+HNCO→IM2→TS2→CH2CONH,CH2(T)+HNCO→IM3→TS3→CH2NH+CO,and CH2(T)+HNCO→IM3→TS3→CH2NH+CO.
基金ACKNOWLEDGMENTS This work was supported by the NationM Nature Science Foundation of China (No.11104256) and the Open Project of State Key Laboratory Cultivation base for Nonmetal Composites and Functional Mate- rials (No.11zxfk19). We thank Dr. Shuang-lin Hu from the Chemistry Department of Uppsala University in Sweden for helpful suggestion. We would also thank the Hefei National Laboratory for Physical Sciences at the Microscale in University of Science and Technology of China for the computational facilities (Gaussian 09).
文摘The potential energy surfaces (PES) of unimolecular dissociation reactions for di-ethyl beryl- lium and di-t-butyl beryllium are investigated by B3LYP, CCSD(T), and G3B3 approaches. Possible reaction pathways through either the radical or transition state (TS) of the molecules are considered. The geometries, vibrational frequencies and relative energies for various sta- tionary points are determined. From the study of energetics, the TS pathways arising from concerted molecular eliminations are indicated to be the main dissociation pathways for both molecules. The PES differences of the dissociation reactions are investigated. The activation energies and rate constants will be helpful for investigating the predictive ability of the reaction in further theoretical and experimental research.