gen.parRep: a first implementation of the Generalized Parallel Replica dynamics for the long time simulation of metastable biochemical systems
arXiv:1807.02431 · doi:10.1016/j.cpc.2019.01.005
Abstract
Metastability is one of the major encountered obstacle when performing long molecular dynamics simulations, and many methods were developed to address this challenge. The "Parallel Replica" (ParRep) dynamics is known for allowing to simulate very long trajectories of metastable Langevin dynamics in the materials science community, but it relies on assumptions that can hardly be transposed to the world of biochemical simulations. The later developed "Generalized ParRep" variant solves those issues, but it was not applied to significant systems of interest so far. In this article, we present the program gen.parRep, the first publicly available implementation of the Generalized Parallel Replica method (BSD 3-Clause license), targeting frequently encountered metastable biochemical systems, such as conformational equilibria or dissociation of protein-ligand complexes. It will be shown that the resulting C++ implementation exhibits a strong linear scalability, providing up to 70 % of the maximum possible speedup on several hundreds of CPUs.
arXiv commentThird version of the article. 37 pages (including references), 13 Figures, 3 Tables
References in corpus (9)
- PLUMED 2: New feathers for an old bird
- String Method for the Study of Rare Events
- Elaborating Transition Interface Sampling Methods
- Improved neighbor list algorithm in molecular simulations using cell decomposition and data sorting method
- A multiple replica approach to simulate reactive trajectories
- A Generalized Parallel Replica Dynamics
- Maximum Flux Transition Paths of Conformational Change
- Accelerated dynamics: Mathematical foundations and algorithmic improvements
- The parallel replica method for simulating long trajectories of Markov chains