Journal
JOURNAL OF CHEMICAL PHYSICS
Volume 151, Issue 22, Pages -Publisher
AMER INST PHYSICS
DOI: 10.1063/1.5127561
Keywords
-
Funding
- Russian Science Foundation [18-13-00479]
- European Regional Development Fund
- Republic of Cyprus through the Research Promotion Foundation [INFRASTRUCTURE/1216/0070, Cy-Tera NEA YPiODeltaOMH/SigmaTPATH/0308/31]
- Russian Science Foundation [18-13-00479] Funding Source: Russian Science Foundation
Ask authors/readers for more resources
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.
Authors
I am an author on this paper
Click your name to claim this paper and add it to your profile.
Reviews
Recommended
No Data Available