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

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 the 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.

[1]  Giovanni Lapenta,et al.  Exactly energy conserving semi-implicit particle in cell formulation , 2016, J. Comput. Phys..

[2]  M. Lyutikov,et al.  Particle acceleration in explosive relativistic reconnection events and Crab Nebula gamma-ray flares , 2016, 1804.10291.

[3]  Jeremiah U. Brackbill,et al.  On energy and momentum conservation in particle-in-cell plasma simulation , 2015, J. Comput. Phys..

[4]  K. Nishikawa,et al.  COLLISIONLESS WEIBEL SHOCKS AND ELECTRON ACCELERATION IN GAMMA-RAY BURSTS , 2015, 1507.05374.

[5]  A. Spitkovsky,et al.  Simultaneous acceleration of protons and electrons at nonrelativistic quasiparallel collisionless shocks. , 2014, Physical review letters.

[6]  A. Spitkovsky,et al.  AB INITIO PULSAR MAGNETOSPHERE: THREE-DIMENSIONAL PARTICLE-IN-CELL SIMULATIONS OF OBLIQUE PULSARS , 2014, 1412.0673.

[7]  L. Hernquist,et al.  NUMERICAL CONVERGENCE IN SMOOTHED PARTICLE HYDRODYNAMICS , 2014, 1410.4222.

[8]  L. Sironi,et al.  RELATIVISTIC PAIR BEAMS FROM TeV BLAZARS: A SOURCE OF REPROCESSED GeV EMISSION RATHER THAN INTERGALACTIC HEATING , 2013, 1312.4538.

[9]  Jacob Trier Frederiksen,et al.  photon-plasma: A modern high-order particle-in-cell code , 2012, 1211.4575.

[10]  A. Broderick,et al.  THE COSMOLOGICAL IMPACT OF LUMINOUS TeV BLAZARS. I. IMPLICATIONS OF PLASMA INSTABILITIES FOR THE INTERGALACTIC MAGNETIC FIELD AND EXTRAGALACTIC GAMMA-RAY BACKGROUND , 2011, 1106.5494.

[11]  Stefano Markidis,et al.  The energy conserving particle-in-cell method , 2011, J. Comput. Phys..

[12]  Stefano Markidis,et al.  Particle acceleration and energy conservation in particle in cell simulations , 2011 .

[13]  Luis Chacón,et al.  An energy- and charge-conserving, implicit, electrostatic particle-in-cell algorithm , 2011, J. Comput. Phys..

[14]  David C. Seal,et al.  A positivity-preserving high-order semi-Lagrangian discontinuous Galerkin scheme for the Vlasov-Poisson equations , 2010, J. Comput. Phys..

[15]  Alexander S. Lipatov,et al.  The Hybrid Multiscale Simulation Technology: An Introduction with Application to Astrophysical and Laboratory Plasmas , 2010 .

[16]  A. Spitkovsky,et al.  NONLINEAR STUDY OF BELL'S COSMIC RAY CURRENT-DRIVEN INSTABILITY , 2008, 0810.4565.

[17]  Anatoly Spitkovsky,et al.  Particle Acceleration in Relativistic Collisionless Shocks: Fermi Process at Last? , 2008, 0802.3216.

[18]  A. Spitkovsky.  Simulations of relativistic collisionless shocks: shock structure and particle acceleration , 2005, astro-ph/0603211.

[19]  Hiroshi Matsumoto,et al.  A new charge conservation method in electromagnetic particle-in-cell simulations , 2003 .

[20]  Nicolas Besse,et al.  Semi-Lagrangian schemes for the Vlasov equation on an unstructured mesh of phase space , 2003 .

[21]  Alexander S. Lipatov,et al.  The Hybrid Multiscale Simulation Technology , 2002 .

[22]  T. Esirkepov,et al.  Exact charge conservation scheme for Particle-in-Cell simulation with an arbitrary form-factor , 2001 .

[23]  E. Liang,et al.  X-ray and gamma-ray Emissioons from Galactic Black Hole candidates: Observations and Theories , 1999 .

[24]  Colin J. McKinstrie,et al.  Accurate formulas for the Landau damping rates of electrostatic waves , 1999 .

[25]  Marco Brambilla,et al.  Kinetic Theory of Plasma Waves: Homogeneous Plasmas , 1998 .

[26]  John D. Villasenor,et al.  Rigorous charge conservation for local electromagnetic field solvers , 1992 .

[27]  J. W. Eastwood,et al.  The virtual particle electromagnetic particle-mesh method , 1991 .

[28]  C. Birdsall,et al.  Plasma Physics via Computer Simulation , 2018 .

[29]  C. Birdsall,et al.  Plasma Physics Via Computer , 1985 .

[30]  J. Brackbill,et al.  Multiple time scales , 1985 .

[31]  Charles K. Birdsall,et al.  Plasma self-heating and saturation due to numerical instabilities☆ , 1980 .

[32]  A. Bruce Langdon,et al.  “Energy-conserving” plasma simulation algorithms , 1973 .

[33]  H. Lewis,et al.  Energy-conserving numerical approximations for Vlasov plasmas , 1970 .

[34]  R W Hockney,et al.  Computer Simulation Using Particles , 1966 .

[35]  John M. Dawson,et al.  One‐Dimensional Plasma Model , 1962 .

[36]  O. Buneman,et al.  Dissipation of Currents in Ionized Media , 1959 .