Numerical heating in particle-in-cell simulations with Monte Carlo binary collisions

Abstract

The binary Monte Carlo (MC) collision algorithm is a standard and robust method to include binary Coulomb collision effects in particle-in-cell (PIC) simulations of plasmas. Here, we show that the coupling between PIC and MC algorithms can give rise to (nonphysical) numerical heating of the system, that significantly exceeds that observed when these algorithms operate independently. We argue that this deleterious effect results from an inconsistency between the particle motion associated with MC-collisions and the work performed by the collective electromagnetic field on the PIC grid. This inconsistency manifests as the (artificial) stochastic production of electromagnetic energy, which ultimately heats the plasma particles. The MC-induced numerical heating can significantly impact the evolution of the simulated system for long simulation times ($\gtrsim 10^3$ collision periods, for typical numerical parameters). We describe the source of the MC-induced numerical heating analytically and discuss strategies to minimize it.

References in corpus (12)

• SMILEI: a collaborative, open-source, multi-purpose particle-in-cell code for plasma simulation
• Particle acceleration in relativistic collisionless shocks: Fermi process at last?
• Dense Electron-Positron Plasmas and Ultra-Intense Bursts of Gamma-Rays from Laser-Irradiated Solids
• Formation of Hard Power-laws in the Energetic Particle Spectra Resulting from Relativistic Magnetic Reconnection
• Interpenetrating plasma shells: near-equipartition magnetic field generation and non-thermal particle acceleration
• Observation of magnetic field generation via the Weibel instability in interpenetrating plasma flows
• One-to-one direct modeling of experiments and astrophysical scenarios: pushing the envelope on kinetic plasma simulations
• Weibel-instability-mediated collisionless shocks in laboratory with ultraintense lasers
• Efficient Nonthermal Particle Acceleration by the Kink Instability in Relativistic Jets
• Large-scale magnetic field generation via the kinetic Kelvin-Helmholtz instability in unmagnetized scenarios
• Statistical kinetic treatment of relativistic binary collisions
• Nonthermal ion acceleration by the kink instability in nonrelativistic jets

Cited by in corpus (3)

• Data-driven discovery of reduced plasma physics models from fully-kinetic simulations
• Exploring ultra-high-intensity wakefields in carbon nanotube arrays: an effective plasma-density approach
• Implementation of higher-order velocity mapping between marker particles and grid in the particle-in-cell code XGC