Nonequilibrium candidate Monte Carlo: A new tool for efficient equilibrium simulation

Abstract

Metropolis Monte Carlo simulation is a powerful tool for studying the equilibrium properties of matter. In complex condensed-phase systems, however, it is difficult to design Monte Carlo moves with high acceptance probabilities that also rapidly sample uncorrelated configurations. Here, we introduce a new class of moves based on nonequilibrium dynamics: candidate configurations are generated through a finite-time process in which a system is actively driven out of equilibrium, and accepted with criteria that preserve the equilibrium distribution. The acceptance rule is similar to the Metropolis acceptance probability, but related to the nonequilibrium work rather than the instantaneous energy difference. Our method is applicable to sampling from both a single thermodynamic state or a mixture of thermodynamic states, and allows both coordinates and thermodynamic parameters to be driven in nonequilibrium proposals. While generating finite-time switching trajectories incurs an additional cost, driving some degrees of freedom while allowing others to evolve naturally can lead to large enhancements in acceptance probabilities, greatly reducing structural correlation times. Using nonequilibrium driven processes vastly expands the repertoire of useful Monte Carlo proposals in simulations of dense solvated systems.

References in corpus (18)

• A nonequilibrium equality for free energy differences
• The Entropy Production Fluctuation Theorem and the Nonequilibrium Work Relation for Free Energy Differences
• Statistically optimal analysis of samples from multiple equilibrium states
• Simulated Tempering: A New Monte Carlo Scheme
• Path ensembles averages in systems driven far-from-equilibrium
• Measuring thermodynamic length
• Optimal finite-time processes in stochastic thermodynamics
• Optimal protocols for minimal work processes in underdamped stochastic thermodynamics
• Optimized free energies from bidirectional single-molecule force spectroscopy
• Computing the optimal protocol for finite-time processes in stochastic thermodynamics
• Equilibrium free energies from fast-switching trajectories with large time steps
• Far-from-Equilibrium Measurements of Thermodynamic Length
• Taking Bigger Metropolis Steps by Dragging Fast Variables
• Optimal estimators and asymptotic variances for nonequilibrium path-ensemble averages
• Free energy reconstruction from steered dynamics without post-processing
• Hamiltonian derivation of a detailed fluctuation theorem
• Estimating equilibrium ensemble averages using multiple time slices from driven nonequilibrium processes: theory and application to free energies, moments, and thermodynamic length in single-molecule pulling experiments
• Langevin dynamics with constraints and computation of free energy differences

Cited by in corpus (22)

• The geometry of thermodynamic control
• Using nonequilibrium fluctuation theorems to understand and correct errors in equilibrium and nonequilibrium discrete Langevin dynamics simulations
• Boundary layers in stochastic thermodynamics
• Annealing glasses by cyclic shear deformation
• Skipping the Replica Exchange Ladder with Normalizing Flows
• Stochastic Normalizing Flows for Inverse Problems: a Markov Chains Viewpoint
• Preserving Correlations Between Trajectories for Efficient Path Sampling
• Accelerating Langevin Sampling with Birth-death
• Stochastic Normalizing Flows
• Computer Simulations of a Heterogeneous Membrane with Enhanced Sampling Techniques
• Dynamical computation of the density of states and Bayes factors using nonequilibrium importance sampling
• Reducing the effect of Metropolization on mixing times in molecular dynamics simulations
• Policy-guided Monte Carlo on general state spaces: Application to glass-forming mixtures
• Grand-canonical molecular dynamics simulations powered by a hybrid 4D nonequilibrium MD/MC method: Implementation in LAMMPS and applications to electrolyte solutions
• A rigorous sequential update strategy for parallel kinetic Monte Carlo simulation
• Donnan equilibrium in charged slit-pores from a hybrid nonequilibrium Molecular Dynamics / Monte Carlo method with ions and solvent exchange
• Transient probability currents provide upper and lower bounds on non-equilibrium steady-state currents in the Smoluchowski picture
• Scaling of Stochastic Normalizing Flows in $\mathrm{SU}(3)$ lattice gauge theory
• Contact geometry and quantum thermodynamics of nanoscale steady states
• Learned Mappings for Targeted Free Energy Perturbation between Peptide Conformations
• Sampling metastable systems using collective variables and Jarzynski-Crooks paths
• Neural Mode Jump Monte Carlo