| Autor: |
Novikov IS; Skolkovo Institute of Science and Technology, Skolkovo Innovation Center, Nobel St. 3, Moscow 143026, Russia., Shapeev AV; Skolkovo Institute of Science and Technology, Skolkovo Innovation Center, Nobel St. 3, Moscow 143026, Russia., Suleimanov YV; Computation-Based Science and Technology Research Center, Cyprus Institute, 20 Kavafi Street, Nicosia 2121, Cyprus. |
| Jazyk: |
angličtina |
| Zdroj: |
The Journal of chemical physics [J Chem Phys] 2019 Dec 14; Vol. 151 (22), pp. 224105. |
| DOI: |
10.1063/1.5127561 |
| Abstrakt: |
Ring polymer molecular dynamics (RPMD) has proven to be an accurate approach for calculating thermal rate coefficients of various chemical reactions. For wider application of this methodology, efficient ways to generate the underlying full-dimensional potential energy surfaces (PESs) and the corresponding energy gradients are required. Recently, we have proposed a fully automated procedure based on combining the original RPMDrate code with active learning for PES on-the-fly using moment tensor potential and successfully applied it to two representative thermally activated chemical reactions [I. S. Novikov et al., Phys. Chem. Chem. Phys. 20, 29503-29512 (2018)]. In this work, using a prototype insertion chemical reaction S + H 2 , we show that this procedure works equally well for another class of chemical reactions. We find that the corresponding PES can be generated by fitting to less than 1500 automatically generated structures, while the RPMD rate coefficients show deviation from the reference values within the typical convergence error of the RPMDrate. We note that more structures are accumulated during the real-time propagation of the dynamic factor (the recrossing factor) as opposed to the previous study. We also observe that a relatively flat free energy profile along the reaction coordinate before entering the complex-formation well can cause issues with locating the maximum of the free energy surface for less converged PESs. However, the final RPMD rate coefficient is independent of the position of the dividing surface that makes it invulnerable to this problem, keeping the total number of necessary structures within a few thousand. Our work concludes that, in the future, the proposed methodology can be applied to realistic complex chemical reactions with various energy profiles. |
| Databáze: |
MEDLINE |
| Externí odkaz: |
|
