Abstract:Context: Interaction processes between discrete particles are often modelled with stochastic methods such as the Representative Particle Monte Carlo (RPMC) method which simulate mutual interactions (e.g.\ chemical reactions, collisions, gravitational stirring) only for a representative subset of $n$ particles instead of all $N$ particles in the system. However, in the traditionally employed computational scheme the memory requirements and the simulation runtime scale quadratically with the number of representative particles. Aims: We want to develop a computational scheme that has significantly lower memory requirements and computational costs than the traditional scheme, so that highly resolved simulations with stochastic processes such as the RPMC method become feasible. Results: In this paper we propose the bucketing scheme, a hybrid sampling scheme that groups similar particles together and combines rejection sampling with a coarsened variant of the traditional discrete inverse transform sampling. For a nu -partite bucket grouping, the storage requirements scale with $n$ and $ and the computational cost per fixed time increment scales with $n both thus being much less sensitive to the number of representative particles $n$. Extensive performance testing demonstrates the higher efficiency and the favourable scaling characteristics of the bucketing scheme compared to the traditional approach, while being statistically equivalent and not introducing any new requirements or approximations. With this improvement, the RPMC method can be efficiently applied not only with very high resolution but also in scenarios where the number of representative particles increases over time, and the simulation of high-frequency interactions (such as gravitational stirring) as a Monte Carlo process becomes viable.

Reweighting in nonequilibrium simulations

Reducing finite-size effects with reweighted renormalization group transformations

Reducing the time-step errors in diffusion Monte Carlo

RiteWeight: Randomized Iterative Trajectory Reweighting for Steady-State Distributions Without Discretization Error

Girsanov reweighting for simulations of underdamped Langevin dynamics. Theory

Optimized Monte Carlo Methods

On reweighting for twisted boundary conditions

A Comment on the Reweighting Method for Accelerated Molecular Dynamics Simulations

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

Reweighted Jarzynski sampling: Acceleration of rare events and free energy calculation with a bias potential learned from nonequilibrium work

Reweighted ensemble dynamics simulations: theory, improvement, and application

Improved Reweighting of Accelerated Molecular Dynamics Simulations for Free Energy Calculation

Multipoint reweighting method and its applications to lattice QCD

Efficient simulation of stochastic interactions among representative Monte Carlo particles

Equilibrium sampling by re-weighting non-equilibrium simulation trajectories

One flavor mass reweighting: foundations

Efficient particle control in systems with large density gradients

Implementation of Girsanov Reweighting in OpenMM and Deeptime