Replica exchange and expanded ensemble simulations as Gibbs sampling: Simple improvements for enhanced mixing

Abstract

The widespread popularity of replica exchange and expanded ensemble algorithms for simulating complex molecular systems in chemistry and biophysics has generated much interest in enhancing phase space mixing of these protocols, thus improving their efficiency. Here, we demonstrate how both of these classes of algorithms can be considered a form of Gibbs sampling within a Markov chain Monte Carlo (MCMC) framework. While the update of the conformational degrees of freedom by Metropolis Monte Carlo or molecular dynamics unavoidably generates correlated samples, we show how judicious updating of the thermodynamic state indices---corresponding to thermodynamic parameters such as temperature or alchemical coupling variables---associated with these configurations can substantially increase mixing while still sampling from the desired distributions. We show how state update methods in common use lead to suboptimal mixing, and present some simple, inexpensive alternatives that can increase mixing of the overall Markov chain, reducing simulation times necessary to obtain estimates of the desired precision. These improved schemes are demonstrated for several common applications, including an alchemical expanded ensemble simulation, parallel tempering, and multidimensional replica exchange umbrella sampling.

17 pages, 2 figures

References in corpus (7)

• Equilibrium free energies from non-equilibrium metadynamics
• Feedback-optimized parallel tempering Monte Carlo
• Optimized parallel tempering simulations of proteins
• Generalized Ensemble and Tempering Simulations: A Unified View
• Conditions for rapid mixing of parallel and simulated tempering on multimodal distributions
• Phase transition in the three dimensional Heisenberg spin glass: Finite-size scaling analysis
• Optimizing Replica Exchange Moves For Molecular Dynamics

Cited by in corpus (24)

• Accelerated weight histogram method for exploring free energy landscapes
• End-to-End Differentiable Molecular Mechanics Force Field Construction
• Non-Reversible Parallel Tempering: a Scalable Highly Parallel MCMC Scheme
• From data to noise to data: mixing physics across temperatures with generative artificial intelligence
• Machine-learned molecular mechanics force field for the simulation of protein-ligand systems and beyond
• Perturbation Potentials to Overcome Order/Disorder Transitions in Alchemical Binding Free Energy Calculations
• Redox potential replica exchange molecular dynamics at constant pH in AMBER: Implementation and validation
• Riemann metric approach to optimal sampling of multidimensional free-energy landscapes
• Statistically optimal continuous free energy surfaces from biased simulations and multistate reweighting
• Skewed Thermodynamic Geometry and Optimal Free Energy Estimation
• Alchemical Grid Dock (AlGDock): Binding Free Energy Calculations between Flexible Ligands and Rigid Receptors
• Ab initio RNA folding
• Estimating thermodynamic expectations and free energies in expanded ensemble simulations: systematic variance reduction through conditioning
• Simulated Tempering and Magnetizing: An Application of Two-Dimensional Simulated Tempering to Two-Dimensional Ising Model and Its Crossover
• Towards an optimal flow: Density-of-states-informed replica-exchange simulations
• Mass-scaling replica-exchange molecular dynamics optimizes computational resources with simpler algorithm
• Designed-walk replica-exchange method for simulations of complex systems
• Multiple Topology Replica Exchange of Expanded Ensembles (MT-REXEE) for Multidimensional Alchemical Calculations
• How cooperatively folding are homopolymer molecular knots?
• Application of Simulated Tempering and Magnetizing to a Two-Dimensional Potts Model
• Deterministic replica-exchange method without pseudo random numbers for simulations of complex systems
• Simulated tempering with irreversible Gibbs sampling techniques
• Multidimensional Replica Exchange simulations for Efficient constant pH and Redox Potential Molecular Dynamics
• High-throughput free energies and water maps for drug discovery by molecular density functional theory