A quasi-static particle-in-cell algorithm based on an azimuthal Fourier decomposition for highly efficient simulations of plasma-based acceleration: QPAD

The 3D quasi-static particle-in-cell (PIC) algorithm is a very efficient method for modeling short-pulse laser or relativistic charged particle beam-plasma interactions. In this algorithm, the plasma response to a non-evolving laser or particle beam is calculated using Maxwell's equations based on the quasi-static approximate equations that exclude radiation. The plasma fields are then used to advance the laser or beam forward using a large time step. The algorithm is many orders of magnitude faster than a 3D fully explicit relativistic electromagnetic PIC algorithm. It has been shown to be capable to accurately model the evolution of lasers and particle beams in a variety of scenarios. At the same time, an algorithm in which the fields, currents and Maxwell equations are decomposed into azimuthal harmonics has been shown to reduce the complexity of a 3D explicit PIC algorithm to that of a 2D algorithm when the expansion is truncated while maintaining accuracy for problems with near azimuthal symmetry. This hybrid algorithm uses a PIC description in r-z and a gridless description in $\phi$. We describe a novel method that combines the quasi-static and hybrid PIC methods. This algorithm expands the fields, charge and current density into azimuthal harmonics. A set of the quasi-static field equations are derived for each harmonic. The complex amplitudes of the fields are then solved using the finite difference method. The beam and plasma particles are advanced in Cartesian coordinates using the total fields. Details on how this algorithm was implemented using a similar workflow to an existing quasi-static code, QuickPIC, are presented. The new code is called QPAD for QuickPIC with Azimuthal Decomposition. Benchmarks and comparisons between a fully 3D explicit PIC code, a full 3D quasi-static code, and the new quasi-static PIC code with azimuthal decomposition are also presented.

[1]  Joyce,et al.  Electron-hose instability in the ion-focused regime. , 1991, Physical review letters.

[2]  Jean-Luc Vay,et al.  A spectral, quasi-cylindrical and dispersion-free Particle-In-Cell algorithm , 2015, Comput. Phys. Commun..

[3]  Jean-Luc Vay,et al.  Numerical stability analysis of the pseudo-spectral analytical time-domain PIC algorithm , 2014, J. Comput. Phys..

[4]  A. E. Dangor,et al.  Monoenergetic beams of relativistic electrons from intense laser–plasma interactions , 2004, Nature.

[5]  Viktor K. Decyk,et al.  Enhancing parallel quasi-static particle-in-cell simulations with a pipelining algorithm , 2009, J. Comput. Phys..

[6]  W. Mori,et al.  Hosing instability in the blow-out regime for plasma-wakefield acceleration. , 2007, Physical review letters.

[7]  A. Jullien,et al.  Relativistic electron beams driven by kHz single-cycle light pulses , 2016, Nature Photonics.

[8]  L. Soby,et al.  Acceleration of electrons in the plasma wakefield of a proton bunch , 2018, Nature.

[9]  Wei Lu,et al.  Modeling of laser wakefield acceleration in Lorentz boosted frame using EM-PIC code with spectral solver , 2014, J. Comput. Phys..

[10]  Eric Esarey,et al.  Tunable laser plasma accelerator based on longitudinal density tailoring , 2011 .

[11]  Su,et al.  Acceleration and focusing of electrons in two-dimensional nonlinear plasma wake fields. , 1991, Physical review. A, Atomic, molecular, and optical physics.

[12]  P. P. Rajeev,et al.  Gamma-rays from harmonically resonant betatron oscillations in a plasma wake , 2011 .

[13]  W. Mori,et al.  Relativistic single-cycle tunable infrared pulses generated from a tailored plasma density structure , 2017, Relativistic Plasma Waves and Particle Beams as Coherent and Incoherent Radiation Sources III.

[14]  J. Vay,et al.  Noninvariance of space- and time-scale ranges under a Lorentz Transformation and the implications for the study of relativistic interactions. , 2007, Physical review letters.

[15]  Thomas M. Antonsen,et al.  An improved iteration loop for the three dimensional quasi-static particle-in-cell algorithm: QuickPIC , 2013, J. Comput. Phys..

[16]  C. Nielson,et al.  Particle-code models in the nonradiative limit , 1976 .

[17]  W. Mori,et al.  Nonlinear theory for relativistic plasma wakefields in the blowout regime. , 2006, Physical review letters.

[18]  B. Mahieu,et al.  High-Brilliance Betatron γ-Ray Source Powered by Laser-Accelerated Electrons. , 2017, Physical review letters.

[19]  Wei Lu,et al.  Energy doubling of 42 GeV electrons in a metre-scale plasma wakefield accelerator , 2007, Nature.

[20]  Robert D. Falgout,et al.  hypre: A Library of High Performance Preconditioners , 2002, International Conference on Computational Science.

[21]  G. White,et al.  High-efficiency acceleration of an electron beam in a plasma wakefield accelerator , 2014, Nature.

[22]  T. M. Antonsen,et al.  QUICKPIC: A highly efficient particle-in-cell code for modeling wakefield acceleration in plasmas , 2006, J. Comput. Phys..

[23]  David H. Whittum,et al.  Transverse two-stream instability of a beam with a Bennett profile , 1997 .

[24]  J. Cary,et al.  High-quality electron beams from a laser wakefield accelerator using plasma-channel guiding , 2004, Nature.

[25]  Carl Schroeder,et al.  HiPACE: a quasi-static particle-in-cell code , 2014 .

[26]  Hans Werner Meuer,et al.  Top500 Supercomputer Sites , 1997 .

[27]  Y. Glinec,et al.  A laser–plasma accelerator producing monoenergetic electron beams , 2004, Nature.

[28]  W. Press,et al.  Numerical Recipes: The Art of Scientific Computing , 1987 .

[29]  K. A. Marsh,et al.  Multi-gigaelectronvolt acceleration of positrons in a self-loaded plasma wakefield , 2015, Nature.

[30]  William H. Press,et al.  Numerical recipes: the art of scientific computing, 3rd Edition , 2007 .

[31]  G. G. Stokes.  "J." , 1890, The New Yale Book of Quotations.

[32]  A. J. Gonsalves,et al.  Multistage coupling of independent laser-plasma accelerators , 2016, Nature.

[33]  Esarey,et al.  Nonlinear interaction of intense laser pulses in plasmas. , 1989, Physical review. A, Atomic, molecular, and optical physics.

[34]  Sanjiva K. Lele,et al.  A highly accurate technique for the treatment of flow equations at the polar axis in cylindrical coordinates using series expansions , 2002 .

[35]  Jean-Luc Vay,et al.  PPPS-2013: Topic 1.2: A domain decomposition method for pseudo-spectral electromagnetic simulations of plasmas , 2013, 2013 Abstracts IEEE International Conference on Plasma Science (ICOPS).

[36]  Wei Lu,et al.  Elimination of the numerical Cerenkov instability for spectral EM-PIC codes , 2015, Comput. Phys. Commun..

[37]  Warren B. Mori,et al.  On numerical errors to the fields surrounding a relativistically moving particle in PIC codes , 2020, J. Comput. Phys..

[38]  Mori,et al.  Group velocity of large-amplitude electromagnetic waves in a plasma. , 1994, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[39]  Jian Liu,et al.  Why is Boris algorithm so good , 2013 .

[40]  W. Mori,et al.  Preservation of beam emittance in the presence of ion motion in future high-energy plasma-wakefield-based colliders. , 2010, Physical review letters.

[41]  Wei Lu,et al.  OSIRIS: A Three-Dimensional, Fully Relativistic Particle in Cell Code for Modeling Plasma Based Accelerators , 2002, International Conference on Computational Science.

[42]  Wei Lu,et al.  Numerical instability due to relativistic plasma drift in EM-PIC simulations , 2012, Comput. Phys. Commun..

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

[44]  J. Rosenzweig,et al.  Effects of ion motion in intense beam-driven plasma wakefield accelerators. , 2005, Physical review letters.

[45]  J. Dawson.  Particle simulation of plasmas , 1983 .

[46]  Thomas M. Antonsen,et al.  Kinetic modeling of intense, short laser pulses propagating in tenuous plasmas , 1997 .

[47]  Zulfikar Najmudin,et al.  Bright spatially coherent synchrotron X-rays from a table-top source , 2010 .

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

[49]  Erik Lefebvre,et al.  Particle-in-Cell modelling of laser-plasma interaction using Fourier decomposition , 2009, J. Comput. Phys..

[50]  K. Lotov,et al.  Fine wakefield structure in the blowout regime of plasma wakefield accelerators , 2003 .

[51]  Wei Lu,et al.  Enabling Lorentz boosted frame particle-in-cell simulations of laser wakefield acceleration in quasi-3D geometry , 2016, J. Comput. Phys..

[52]  R. A. Fonseca,et al.  Implementation of a hybrid particle code with a PIC description in r-z and a gridless description in ϕ into OSIRIS , 2015, J. Comput. Phys..

[53]  P. Sprangle,et al.  Nonlinear theory of intense laser-plasma interactions. , 1990, Physical review letters.

[54]  W. Mori,et al.  Ion Motion Induced Emittance Growth of Matched Electron Beams in Plasma Wakefields. , 2017, Physical review letters.

[55]  T. Ditmire,et al.  Quasi-monoenergetic laser-plasma acceleration of electrons to 2 GeV , 2013, Nature Communications.