Optimized routines for event generators in QED-PIC codes

Abstract

In recent years, the prospects of performing fundamental and applied studies at the next-generation high-intensity laser facilities have greatly stimulated the interest in performing large-scale simulations of laser interaction with matter with the account for quantum electrodynamics (QED) processes such as emission of high energy photons and decay of such photons into electron-positron pairs. These processes can be modeled via probabilistic routines that include frequent computation of synchrotron functions and can constitute significant computational demands within accordingly extended Particle-in-Cell (QED-PIC) algorithms. In this regard, the optimization of these routines is of great interest. In this paper, we propose and describe two modifications. First, we derive a more accurate upper-bound estimate for the rate of QED events and use it to arrange local sub-stepping of the global time step in a significantly more efficient way than done previously. Second, we present a new high-performance implementation of synchrotron functions. Our optimizations made it possible to speed up the computations by a factor of up to 13.7 depending on the problem. Our implementation is integrated into the PICADOR and Hi-Chi codes, the latter of which is distributed publicly (https://github.com/hi-chi/pyHiChi).

References in corpus (13)

• SMILEI: a collaborative, open-source, multi-purpose particle-in-cell code for plasma simulation
• Modelling Gamma Ray Emission and Pair Production in High-Intensity Laser-Matter Interactions
• QED cascades induced by circularly polarized laser fields
• Laser absorption via QED cascades in counter propagating laser pulses
• Warp-X: a new exascale computing platform for beam-plasma simulations
• Exploiting multi-scale parallelism for large scale numerical modelling of laser wakefield accelerators
• Electron dynamics, gamma and electron-positron production by colliding laser pulses
• Quantum quenching of radiation losses in short laser pulses
• Particle-in-Cell Laser-Plasma Simulation on Xeon Phi Coprocessors
• Laser-driven plasma pinching in $e^{-}{e}^{+}$ cascade
• Narrowing of the emission angle in high-intensity Compton scattering
• Analytical Fits to the Synchrotron Functions
• Particle dynamics and spatial $e^{-}{e}^{+}$ density structures at QED cascading in circularly polarized standing waves

Cited by in corpus (3)

• Charged particle motion and radiation in strong electromagnetic fields
• Optimized event generator for strong-field QED simulations within the hi-$\chi$ framework
• High-Performance Implementation of the Optimized Event Generator for Strong-Field QED Plasma Simulations