Ring Polymer Molecular Dynamics and Active Learning of Moment Tensor Potential for Gas-Phase Barrierless Reactions: Application to S + H2

Abstract

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, Y. V. Suleimanov, A. V. Shapeev, Phys. Chem. Chem. Phys., 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 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 relatively flat free energy profile of the 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 future, the proposed methodology can be applied to realistic complex chemical reactions with various energy profiles.

References in corpus (14)

• Moment Tensor Potentials: a class of systematically improvable interatomic potentials
• Less is more: sampling chemical space with active learning
• Active learning of linearly parametrized interatomic potentials
• Phase transitions of hybrid perovskites simulated by machine-learning force fields trained on-the-fly with Bayesian inference
• A KInetic Database for Astrochemistry (KIDA)
• Active Learning of Uniformly Accurate Inter-atomic Potentials for Materials Simulation
• Accelerating crystal structure prediction by machine-learning interatomic potentials with active learning
• Impact of lattice relaxations on phase transitions in a high-entropy alloy studied by machine-learning potentials
• Machine learning of molecular properties: locality and active learning
• Chemical Reaction Rates from Ring Polymer Molecular Dynamics: Theory and Practical Applications
• Automated Calculation of Thermal Rate Coefficients using Ring Polymer Molecular Dynamics and Machine-Learning Interatomic Potentials with Active Learning
• An experimental and theoretical investigation of the C(1D) + D2 reaction
• Thermal Rate Coefficients for the Astrochemical Process C + CH${}^{+}$ $\to$ C${}_2^{+}$ + H by Ring Polymer Molecular Dynamics
• The low temperature D${}^{+}$ + H${}_2$ $\to$ HD + H${}^{+}$ reaction rate coefficient: a ring polymer molecular dynamics and quasi-classical trajectory study

Cited by in corpus (3)

• A Deep Dive into Machine Learning Density Functional Theory for Materials Science and Chemistry
• Dynamical Strengthening of Covalent and Non-Covalent Molecular Interactions by Nuclear Quantum Effects at Finite Temperature
• Couplings for Andersen Dynamics