An implementation of the variational quantum RRKM program is presented to utilize the direct ab initio dynamics approach for calculating k(tE, J), k(E) and k(T) within the framework of the microcanonical transition st...An implementation of the variational quantum RRKM program is presented to utilize the direct ab initio dynamics approach for calculating k(tE, J), k(E) and k(T) within the framework of the microcanonical transition state (μTST) and microcanonical variational TST (μVT) theories. An algorithm including tunneling contributions in Beyer-Swinehart method for calculating microcanonical rate constants is also proposed. An efficient piece-wise interpolation method is developed to evaluate the Boltzmann integral in calculation of thermal rate constants. Calculations on several test reactions, namely the H(D)2CO→ H(D)2 + CO, CH2CO→CH2 + CO and CH4 + H→CH3 + H2 reactions, show that the results are in good agreement with the previous rate constants calculations. This approach would require much less computational resource.展开更多
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 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.展开更多
基金supported by the Pre-research Special Project(Grant No.2001CCD03500)of the 973 Fundamental Research Project.
文摘An implementation of the variational quantum RRKM program is presented to utilize the direct ab initio dynamics approach for calculating k(tE, J), k(E) and k(T) within the framework of the microcanonical transition state (μTST) and microcanonical variational TST (μVT) theories. An algorithm including tunneling contributions in Beyer-Swinehart method for calculating microcanonical rate constants is also proposed. An efficient piece-wise interpolation method is developed to evaluate the Boltzmann integral in calculation of thermal rate constants. Calculations on several test reactions, namely the H(D)2CO→ H(D)2 + CO, CH2CO→CH2 + CO and CH4 + H→CH3 + H2 reactions, show that the results are in good agreement with the previous rate constants calculations. This approach would require much less computational resource.
基金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 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.
