Resetting Metadynamics

Ofir Blumer,Shlomi Reuveni,Barak Hirshberg
2023-07-12
Abstract:Metadynamics is a powerful method to accelerate molecular dynamics simulations, but its efficiency critically depends on the identification of collective variables that capture the slow modes of the process. Unfortunately, collective variables are usually not known a priori, and finding them can be very challenging. We recently presented a collective variables-free approach to enhanced sampling using stochastic resetting. Here, we combine the two methods for the first time, showing that it can lead to greater acceleration than either of them separately. We also demonstrate that resetting Metadynamics simulations performed with suboptimal collective variables can lead to speedups comparable with those obtained with optimal collective variables. Therefore, the application of stochastic resetting can be an alternative to the challenging task of improving suboptimal collective variables, at almost no additional computational cost. Finally, we propose a method to extract unbiased mean first-passage times from Metadynamics simulations with resetting, resulting in an improved tradeoff between speedup and accuracy. This work opens the way for combining stochastic resetting with other enhanced sampling methods to accelerate a broad range of molecular simulations.
Chemical Physics,Statistical Mechanics,Computational Physics
What problem does this paper attempt to address?
The paper primarily explores how to enhance the sampling efficiency in molecular dynamics simulations by combining two methods—Metadynamics (MetaD) and Stochastic Resetting (SR). Specifically, the paper attempts to address the following core issues: 1. **Overcoming the timescale problem of slow processes**: Molecular dynamics simulations are a powerful tool at atomic-level time and spatial resolution, but their application is limited to simulating microsecond-scale processes. Phenomena such as protein folding, which occur over longer timescales, are difficult to simulate. 2. **Optimizing the selection of Collective Variables (CVs)**: The Metadynamics method relies on the selection of effective collective variables that can capture slow-changing modes. These variables should distinguish between the metastable states of interest and the transition dynamics between them. However, finding suitable collective variables in practice is highly challenging. 3. **Improving sampling efficiency**: Although Metadynamics is an effective enhanced sampling technique, its efficiency heavily depends on the quality of the collective variables. When the collective variables are suboptimal, it may lead to poor unbiased free energy surfaces or kinetic inference results. The main contribution of the paper is the proposal of a new strategy that combines stochastic resetting with the Metadynamics method to achieve higher acceleration. Specifically: - **Combining SR with MetaD**: The paper demonstrates that even with the optimal collective variables for Metadynamics simulations, adding stochastic resetting can further accelerate the simulation process. More importantly, for cases using suboptimal collective variables, combining stochastic resetting can achieve acceleration effects similar to those using optimal collective variables. This provides an alternative solution to the challenge of selecting collective variables in complex systems. - **Evaluating the effect of stochastic resetting**: The researchers propose a simple method to evaluate whether stochastic resetting can accelerate the simulation and what the optimal resetting rate is. Additionally, they developed a method to extract unbiased mean first passage times from Metadynamics simulations that include stochastic resetting, thereby achieving a better balance between speed and accuracy. In summary, this paper aims to significantly enhance the efficiency of molecular dynamics simulations by combining stochastic resetting and Metadynamics, especially when dealing with complex chemical and physical systems.