Efficient simulation of stochastic interactions among representative Monte Carlo particles

• Published in: Astronomy and astrophysics (Berlin) 2024-04, Vol.684, p.A2
• Main Authors: Beutel, M., Dullemond, C. P., Strzodka, R.
• Format: Article
• Language: English
• Subjects: Chemical reactions; Computing costs; Run time (computers); Sampling; Simulation; Stirring; Stochastic processes
• ISSN: 0004-6361; 1432-0746
• DOI: 10.1051/0004-6361/202347304

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 v-partite bucket grouping, the storage requirements scale with n and v 2 , and the computational cost per fixed time increment scales with n ⋅ v , 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.