SHARP: A Spatially Higher-order, Relativistic Particle-in-Cell Code

Abstract

Numerical heating in particle-in-cell (PIC) codes currently precludes the accurate simulation of cold, relativistic plasma over long periods, severely limiting their applications in astrophysical environments. We present a spatially higher-order accurate relativistic PIC algorithm in one spatial dimension, which conserves charge and momentum exactly. We utilize the smoothness implied by the usage of higher-order interpolation functions to achieve a spatially higher-order accurate algorithm (up to fifth order). We validate our algorithm against several test problems -- thermal stability of stationary plasma, stability of linear plasma waves, and two-stream instability in the relativistic and non-relativistic regimes. Comparing our simulations to exact solutions of the dispersion relations, we demonstrate that SHARP can quantitatively reproduce important kinetic features of the linear regime. Our simulations have a superior ability to control energy non-conservation and avoid numerical heating in comparison to common second-order schemes. We provide a natural definition for convergence of a general PIC algorithm: the complement of physical modes captured by the simulation, i.e., those that lie above the Poisson noise, must grow commensurately with the resolution. This implies that it is necessary to simultaneously increase the number of particles per cell and decrease the cell size. We demonstrate that traditional ways for testing for convergence fail, leading to plateauing of the energy error. This new PIC code enables us to faithfully study the long-term evolution of plasma problems that require absolute control of the energy and momentum conservation.

26 pages, 19 figures, discussion about performance is added, published in ApJ

References in corpus (4)

• Particle acceleration in relativistic collisionless shocks: Fermi process at last?
• Simultaneous Acceleration of Protons and Electrons at Nonrelativistic Quasiparallel Collisionless Shocks
• Ab-initio pulsar magnetosphere: three-dimensional particle-in-cell simulations of oblique pulsars
• Numerical Convergence in Smoothed Particle Hydrodynamics

Cited by in corpus (20)

• Electron Heating in Low Mach Number Perpendicular shocks. I. Heating Mechanism
• A new cosmic ray-driven instability
• Importance of resolving the spectral support of beam-plasma instabilities in simulations
• Growth of beam-plasma instabilities in the presence of background inhomogeneity
• The growth of the longitudinal beam-plasma instability in the presence of an inhomogeneous background
• Lepton-driven Non-resonant Streaming Instability
• Deciphering the physical basis of the intermediate-scale instability
• A Particle-In-Cell Code Comparison for Ion Acceleration: EPOCH, LSP, and WarpX
• Refining pulsar radio emission due to streaming instabilities: Linear theory and PIC simulations in a wide parameter range
• Radio Emission by Soliton Formation in Relativistically Hot Streaming Pulsar Pair Plasmas
• Parallel transport modeling of linear divertor simulators with fundamental ion cyclotron heating
• A guiding center implementation for relativistic particle dynamics in the PLUTO code
• Energy Dissipation in Strong Collisionless Shocks: The Crucial Role of Ion-to-Electron Scale Separation in Particle-in-Cell Simulations
• AHKASH: a new Hybrid particle-in-cell code for simulations of astrophysical collisionless plasma
• Constraining blazar heating with the 2<z<3 Lyman-$\alpha$ forest
• Parallel Collisionless Shocks in strongly Magnetized Electron-Ion Plasma. I. Temperature anisotropies
• Thermal fluid closures and pressure anisotropies in numerical simulations of plasma wakefield acceleration
• Bunch Expansion as a cause for Pulsar Radio Emissions
• The mechanism of efficient electron acceleration at parallel non-relativistic shocks
• Kinetic Simulations of Strongly-Magnetized Parallel Shocks: Deviations from MHD Jump Conditions