Non-empirical weighted Langevin mechanics for the potential escape problem: parallel algorithm and application to the Argon clusters

Abstract

Recently a non-empirical stochastic walker algorithm has been developed to search for the minimum-energy escape paths (MEP) from the minima of the potential surface [J. Phys. Soc. Jpn. 87, 063801 (2018)]. This algorithm is novel in that it tracks the MEP monotonically and does not use the whole Hessian matrix but only gradient and Laplacian of the potential. In this work, we implement an MPI-parallelized version of this algorithm in a simple way. We also explore efficient ways to reduce the number of walkers required for the accurate tracking of the MEP and generate initial positions automatically. We apply the whole scheme to the Lennard-Jones argon cluster with 7-38 atoms to demonstrate the successful tracking of the reaction paths. This achievement paves the path to non-empirical simulation of rare reactions without coarse-graining or artificial potential.

arXiv comment30 pages (single column), 13 figures

References in corpus (15)

• Escaping free-energy minima
• Deep Potential Molecular Dynamics: a scalable model with the accuracy of quantum mechanics
• Sampling rare switching events in biochemical networks
• Direct evaluation of large-deviation functions
• Nonequilibrium Markov processes conditioned on large deviations
• The double-funnel energy landscape of the 38-atom Lennard-Jones cluster
• Evolution of the Potential Energy Surface with Size for Lennard-Jones Clusters
• Simulating rare events in dynamical processes
• Phase changes in 38 atom Lennard-Jones clusters. II: A parallel tempering study of equilibrium and dynamic properties in the molecular dynamics and microcanonical
• Variational and optimal control representations of conditioned and driven processes
• Probing rare physical trajectories with Lyapunov weighted dynamics
• The Gentlest Ascent Dynamics
• Effective interactions and large deviations in stochastic processes
• Efficient characterisation of large deviations using population dynamics
• Stochastic formalism for thermally driven distribution frontier: A nonempirical approach to the potential escape problem

Cited by in corpus (1)

• AtomREM: Non-empirical seeker of the minimum energy escape paths on many-dimensional potential landscapes without coarse graining