Monte Carlo simulation of particle interactions at high dynamic range: Advancing beyond the Googol

Abstract

We present a method which extends Monte Carlo studies to situations that require a large dynamic range in particle number. The underlying idea is that, in order to calculate the collisional evolution of a system, some particle interactions are more important than others and require more resolution, while the behavior of the less important, usually of smaller mass, particles can be considered collectively. In this approximation groups of identical particles, sharing the same mass and structural parameters, operate as one unit. The amount of grouping is determined by the zoom factor -- a free parameter that determines on which particles the computational effort is focused. Two methods for choosing the zoom factors are discussed: the `equal mass method,' in which the groups trace the mass density of the distribution, and the `distribution method,' which additionally follows fluctuations in the distribution. Both methods achieve excellent correspondence with analytic solutions to the Smoluchowski coagulation equation. The grouping method is furthermore applied to simulations involving runaway kernels, where the particle interaction rate is a strong function of particle mass, and to situations that include catastrophic fragmentation. For the runaway simulations previous predictions for the decrease of the runaway timescale with the initial number of particles $\mathcal{N}$ are reconfirmed, extending $\mathcal{N}$ to $10^{160}$. Astrophysical applications include modeling of dust coagulation, planetesimal accretion, and the dynamical evolution of stars in large globular clusters. The proposed method is a powerful tool to compute the evolution of any system where the particles interact through discrete events, with the particle properties characterized by structural parameters.

18 pages, 10 figures. Re-submitted to ApJ with comments of the referee included

References in corpus (14)

• Coagulation, fragmentation and radial motion of solid particles in protoplanetary disks
• Dust coagulation in protoplanetary disks: a rapid depletion of small grains
• Closed-form expressions for particle relative velocities induced by turbulence
• Scaling Theory and Exactly Solved Models In the Kinetics of Irreversible Aggregation
• Planet Formation by Coagulation: A Focus on Uranus and Neptune
• Dust coagulation in protoplanetary disks: porosity matters
• Monte-Carlo Simulations of Globular Cluster Evolution - I. Method and Test Calculations
• Resolving the Formation of Protogalaxies. I. Virialization
• A Hybrid N-body--Coagulation Code for Planet Formation
• Runaway collisions in young star clusters. I. Methods and tests
• A Survey of Numerical Solutions to the Coagulation Equation
• Stochastic Coagulation and the Timescale for Runaway Growth
• Solving the Coagulation Equation by the Moments Method
• A runaway collision in a young star cluster as the origin of the brightest supernova

Cited by in corpus (31)

• The outcome of protoplanetary dust growth: pebbles, boulders, or planetesimals? II. Introducing the bouncing barrier
• Dust Evolution and the Formation of Planetesimals
• Dust coagulation and fragmentation in molecular clouds. I. How collisions between dust aggregates alter the dust size distribution
• Planetesimal formation by sweep-up: How the bouncing barrier can be beneficial to growth
• On the water delivery to terrestrial embryos by ice pebble accretion
• Adding particle collisions to the formation of asteroids and Kuiper belt objects via streaming instabilities
• Numerical Modeling of the Coagulation and Porosity Evolution of Dust Aggregates
• Dust growth and evolution in protoplanetary disks
• Erosion and the limits to planetesimal growth
• Electrostatic Barrier against Dust Growth in Protoplanetary Disks. I. Classifying the Evolution of Size Distribution
• Residence Times of Particles in Diffusive Protoplanetary Disk Environments I. Vertical Motions
• A new condition for the transition from runaway to oligarchic growth
• Accretion among preplanetary bodies: the many faces of runaway growth
• Effect of nucleation on icy pebble growth in protoplanetary discs
• Modeling dust growth in protoplanetary disks: The breakthrough case
• A method for coupling dynamical and collisional evolution of dust in circumstellar disks: the effect of a dead zone
• Dust diffusion and settling in the presence of collisions: Trapping (sub)micron grains in the midplane
• Influence of grain growth on the thermal structure of protoplanetary discs
• Self-sustained Levitation of Dust Aggregate Ensembles by Temperature Gradient Induced Overpressures
• Dust growth, fragmentation and self-induced dust traps in PHANTOM
• Coagulation Instability in Protoplanetary Disks: A Novel Mechanism Connecting Collisional Growth and Hydrodynamical Clumping of Dust Particles
• Haze Formation on Triton
• Breakthrough revisited: investigating the requirements for growth of dust beyond the bouncing barrier
• Dust Accumulation near the Magnetospheric Truncation of Protoplanetary Discs around T Tauri Stars
• Dust Growth and Dynamics in Protoplanetary Nebulae: Implications for Opacity, Thermal Profile and Gravitational Instability
• Massive Protostellar Disks as a Hot Laboratory of Silicate Grain Evolution
• Nonlinear Outcome of Coagulation Instability in Protoplanetary Disks II: Dust Ring Formation Mediated by Backreaction and Fragmentation
• Nonlinear Outcome of Coagulation Instability in Protoplanetary Disks I: First Numerical Study of Accelerated Dust Growth and Dust Concentration at Outer Radii
• Conglomeration of kilometre-sized planetesimals
• SHAMPOO: A stochastic model for tracking dust particles under the influence of non-local disk processes
• Evading the dust fragmentation barrier with the streaming instability in protoplanetary disks