Enhanced Stopping of Macro-Particles in Particle-in-Cell Simulations

Abstract

We derive an equation for energy transfer from relativistic charged particles to a cold background plasma appropriate for finite-size particles that are used in particle-in-cell simulation codes. Expressions for one-, two-, and three-dimensional particles are presented, with special attention given to the two-dimensional case. This energy transfer is due to the electric field of the wake set up in the background plasma by the relativistic particle. The enhanced stopping is dependent on the $q^2/m$, where $q$ is the charge and $m$ is the mass of the relativistic particle, and therefore simulation macro-particles with large charge but identical $q/m$ will stop more rapidly. The stopping power also depends on the effective particle shape of the macro-particle. These conclusions are verified in particle-in-cell simulations. We present 2D simulations of test particles, relaxation of high-energy tails, and integrated fast ignition simulations showing that the enhanced drag on macro-particles may adversely affect the results of these simulations in a wide range of high-energy density plasma scenarios. We also describe a particle splitting algorithm which can potentially overcome this problem and show its effect in controlling the stopping of macro-particles.

29 pages, 7 figures; v2 - added reference; v3 - improved fig. 1&2

References in corpus (6)

• On the structure of relativistic collisionless shocks in electron-ion plasmas
• One-to-one direct modeling of experiments and astrophysical scenarios: pushing the envelope on kinetic plasma simulations
• Exploiting multi-scale parallelism for large scale numerical modelling of laser wakefield accelerators
• Laser-plasma interactions for fast ignition
• Statistical kinetic treatment of relativistic binary collisions
• Energy Loss of High-Energy Particles in Particle-in-Cell Simulation

Cited by in corpus (7)

• Persistence of precursor waves in two-dimensional relativistic shocks
• The energetics of relativistic magnetic reconnection: ion-electron repartition and particle distribution hardness
• Mitigation of numerical Cerenkov radiation and instability using a hybrid finite difference-FFT Maxwell solver and a local charge conserving current deposit
• Non-thermal particle acceleration at highly oblique non-relativistic shocks
• Two-stage electron acceleration by 3D collisionless guide field magnetic reconnection
• Extended particle absorber for efficient modeling of intense laser-solid interactions
• Stopping Power Enhancement From Discrete Particle-Wake Correlations in High Energy Density Plasmas