Steady-state simulations using weighted ensemble path sampling

Abstract

We extend the weighted ensemble (WE) path sampling method to perform rigorous statistical sampling for systems at steady state. The straightforward steady-state implementation of WE is directly practical for simple landscapes, but not when significant metastable intermediates states are present. We therefore develop an enhanced WE scheme, building on existing ideas, which accelerates attainment of steady state in complex systems. We apply both WE approaches to several model systems confirming their correctness and efficiency by comparison with brute-force results. The enhanced version is significantly faster than the brute force and straightforward WE for systems with WE bins that accurately reflect the reaction coordinate(s). The new WE methods can also be applied to equilibrium sampling, since equilibrium is a steady state.

References in corpus (4)

• Forward Flux Sampling-type schemes for simulating rare events: Efficiency analysis
• Peptide Folding Kinetics from Replica Exchange Molecular Dynamics
• Computing stationary distributions in equilibrium and non-equilibrium systems with Forward Flux Sampling
• Efficient path sampling on multiple reaction channels

Cited by in corpus (19)

• Enhanced sampling methods for molecular dynamics simulations
• Simulating rare events using a Weighted Ensemble-based string method
• Efficient Stochastic Simulation of Chemical Kinetics Networks using a Weighted Ensemble of Trajectories
• Weighted Ensemble Milestoning (WEM): A Combined Approach for Rare Event Simulations
• Non-Stationary Forward Flux Sampling
• Steered Transition Path Sampling
• Optimizing weighted ensemble sampling of steady states
• Analysis and optimization of weighted ensemble sampling
• Weighted ensemble: Recent mathematical developments
• Constant pH Simulation with FMM Electrostatics in GROMACS. (A) Design and Applications
• Investigating the Role of Non-Covalent Interactions in Conformation and Assembly of Triazine-Based Sequence-Defined Polymers
• Computing transition path theory quantities with trajectory stratification
• Weighted-ensemble Brownian dynamics simulation: Sampling of rare events in non-equilibrium systems
• Efficiently Sampling Conformations and Pathways Using the Concurrent Adaptive Sampling (CAS) Algorithm
• An ergodic theorem for the weighted ensemble method
• WeTICA: A directed search weighted ensemble based enhanced sampling method to estimate rare event kinetics in a reduced dimensional space
• Transient probability currents provide upper and lower bounds on non-equilibrium steady-state currents in the Smoluchowski picture
• A Study on the Performance of Distributed Training of Data-driven CFD Simulations
• Sampling Rare Event Energy Landscapes via Birth-Death Augmented Dynamics