Sampling metastable systems using collective variables and Jarzynski-Crooks paths
arXiv:2405.18160 · doi:10.1016/j.jcp.2025.113806
Abstract
We consider the problem of sampling a high dimensional multimodal target probability measure. We assume that a good proposal kernel to move only a subset of the degrees of freedoms (also known as collective variables) is known a priori. This proposal kernel can for example be built using normalizing flows. We show how to extend the move from the collective variable space to the full space and how to implement an accept-reject step in order to get a reversible chain with respect to a target probability measure. The accept-reject step does not require to know the marginal of the original measure in the collective variable (namely to know the free energy). The obtained algorithm admits several variants, some of them being very close to methods which have been proposed previously in the literature. We show how the obtained acceptance ratio can be expressed in terms of the work which appears in the Jarzynski-Crooks equality, at least for some variants. Numerical illustrations demonstrate the efficiency of the approach on various simple test cases, and allow us to compare the variants of the algorithm.
Rewritten introduction, added references, small corrections and clarifications
References in corpus (14)
- Escaping free-energy minima
- Well-Tempered Metadynamics: A Smoothly Converging and Tunable Free-Energy Method
- The Entropy Production Fluctuation Theorem and the Nonequilibrium Work Relation for Free Energy Differences
- Statistically optimal analysis of samples from multiple equilibrium states
- Equilibrium free energy differences from nonequilibrium measurements: a master equation approach
- Enhanced sampling methods for molecular dynamics simulations
- Neural Network Renormalization Group
- Nonequilibrium candidate Monte Carlo: A new tool for efficient equilibrium simulation
- Conditions for rapid mixing of parallel and simulated tempering on multimodal distributions
- Non-Reversible Parallel Tempering: a Scalable Highly Parallel MCMC Scheme
- Efficient Modelling of Trivializing Maps for Lattice Theory Using Normalizing Flows: A First Look at Scalability
- Escorted free energy simulations
- Computation of free energy differences through nonequilibrium stochastic dynamics: the reaction coordinate case
- Free energy reconstruction from steered dynamics without post-processing