A Sinking Approach to Explore Arbitrary Areas in Free Energy Landscapes

Abstract

To address the time-scale limitations in molecular dynamics (MD) simulations, numerous enhanced sampling methods have been developed to expedite the exploration of complex free energy landscapes. A commonly employed approach accelerates the sampling of degrees of freedom associated with pre-defined collective variables (CVs), which typically tends to traverse the entire CV range. However, in many scenarios, the focus of interest is on specific regions within the CV space. This paper introduces a novel "sinking" approach that enables enhanced sampling of arbitrary areas within the CV space. We begin by proposing a gridded convolutional approximation that productively replicates the effects of metadynamics, a powerful CV-based enhanced sampling technique. Building on this, we present the SinkMeta method, which "sinks" the interior bias potential to create restraining potential "cliffs" at the grid edges. This technique can confine the exploration of CVs in MD simulations to a preset area. Our experimental results demonstrate that SinkMeta requires minimal sampling steps to estimate the free energy landscape for CV subspaces of various shapes and dimensions, including irregular two-dimensional regions and one-dimensional pathways between metastable states. We believe that SinkMeta will pioneer a new paradigm for sampling partial phase spaces, especially offering an efficient and flexible solution for sampling minimum free energy paths in high-dimensional spaces.

arXiv comment29 pages, 6 figures

References in corpus (21)

• Escaping free-energy minima
• Well-Tempered Metadynamics: A Smoothly Converging and Tunable Free-Energy Method
• PLUMED: a portable plugin for free-energy calculations with molecular dynamics
• Simulated Tempering: A New Monte Carlo Scheme
• String Method for the Study of Rare Events
• Identification of slow molecular order parameters for Markov model construction
• Rethinking Metadynamics: from bias potentials to probability distributions
• Enhanced sampling methods for molecular dynamics simulations
• A Variational Approach to Enhanced Sampling and Free Energy Calculations
• Caliber based spectral gap optimization of order parameters (SGOOP) for sampling complex molecular systems
• Data-Driven Collective Variables for Enhanced Sampling
• A self-learning algorithm for biased molecular dynamics
• A variational conformational dynamics approach to the selection of collective variables in metadynamics
• Homogeneous nucleation of ice
• A Perspective on Deep Learning for Molecular Modeling and Simulations
• Gaussian Mixture Based Enhanced Sampling
• Making the best of a bad situation: a multiscale approach to free energy calculation
• Combining Metadynamics and Integrated Tempering Sampling
• Enhanced sampling of transition states
• Atomistic View of Homogeneous Nucleation of Water into Polymorphic Ices
• A Generalized Nucleation Theory for Ice Crystallization