Ring polymer molecular dynamics and active learning of moment tensor potential for gas-phase barrierless reactions: Application to S + H2.

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 + H2, 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.