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.

Agnostic conservative down-sampling for optimizing statistical representations and PIC simulations

Efficient simulation of stochastic interactions among representative Monte Carlo particles

Efficient particle control in systems with large density gradients

Reducing Noise for PIC Simulations Using Kernel Density Estimation Algorithm

A Fine Resampling Algorithm for General Particle Filters

New Macroparticle Coalescing Models That Conserve Particle'S Phase-Space Distribution in 3-D Particle-In-Cell Simulations of Plasmas

Metropolis-adjusted interacting particle sampling

Inference and De-Noising of Non-Gaussian Particle Distribution Functions: A Generative Modeling Approach

Denoising Particle-In-Cell Data via Smoothness-Increasing Accuracy-Conserving Filters with Application to Bohm Speed Computation

Explicit energy-conserving modification of relativistic PIC method

Importance sampling method for Feldman-Cousins confidence intervals

A representative particle approach to coagulation and fragmentation of dust aggregates and fluid droplets

Efficient and Unbiased Sampling of Boltzmann Distributions via Consistency Models

Use the force! Reduced variance estimators for densities, radial distribution functions and local mobilities in molecular simulations

Particles on demand method: Theoretical analysis, simplification techniques, and model extensions

A dynamical particle merging and splitting algorithm for Particle-In-Cell simulations

High-Order Interpolation Algorithms for Charge Conservation in Particle-in-Cell Simulations

Predicting the statistical error of analog particle tracing Monte Carlo

Superposition Enhanced Nested Sampling