Abstract:Numerical simulation of strong-field quantum electrodynamics (SFQED) processes is an essential step towards current and future high-intensity laser experiments. The complexity of SFQED phenomena and their stochastic nature make them extremely computationally challenging, requiring the use of supercomputers for realistic simulations. Recently, we have presented a novel approach to numerical simulation of SFQED processes based on an accurate approximation of precomputed rates, which minimizes the number of rate calculations per QED event. The current paper is focused on the high-performance implementation of this method, including vectorization of resource-intensive kernels and improvement of parallel computing efficiency. Using two codes, PICADOR and hi-$\chi$ (the latter being free and publicly available), we demonstrate significant reduction in computation time due to these improvements. We hope that the proposed approach can be applied in other codes for the numerical simulation of SFQED processes.

What problem does this paper attempt to address?

### Problems the paper attempts to solve This paper aims to solve the computational efficiency problem in the numerical simulation of strong - field quantum electrodynamics (SFQED) processes. Specifically, the paper proposes a high - performance implementation method to optimize the SFQED event generator, in order to reduce the computation time and improve the efficiency of parallel computing. The main problems include: 1. **High demand for computational resources**: SFQED processes are highly random and nonlinear, resulting in a large amount of computational resources required in the simulation. Especially in the interaction between a strong laser field and plasma, the number of particles may increase exponentially, which will significantly increase the requirements for memory and computational resources. 2. **Time - step limitation**: Traditional SFQED simulation methods usually require very small time - steps to analyze these processes, which will lead to huge computational overhead. The method proposed in the paper reduces the strict requirements for time - steps by minimizing the number of rate calculations for each QED event. 3. **Parallel computing efficiency**: In order to conduct large - scale simulations on supercomputers, an efficient parallel computing strategy is required. The paper solves this problem by vectorization and improving the efficiency of parallel computing. ### Main contributions 1. **Algorithm optimization**: The paper proposes an optimized algorithm implementation. Based on pre - calculated rate approximations, it reduces the number of rate calculations for each QED event, thereby increasing the simulation speed by 30 - 60%. 2. **Performance improvement methods**: It describes and analyzes in detail the methods for improving the performance of QED simulations, hoping that these methods can be applied to other codes. 3. **Practical applications**: These methods are implemented and tested in two codes (PICADOR and Hi - χ). Hi - χ is an open - source project and can be obtained in the source code. ### Key technologies 1. **Vectorization**: By using SIMD (Single Instruction Multiple Data) technology, the computational efficiency of resource - intensive kernels is optimized. 2. **Parallel computing**: The efficiency of parallel computing is improved, especially when dealing with newly generated particles. By saving new particles into separate arrays for each thread, data conflicts are avoided. 3. **Data structure optimization**: An appropriate memory layout (such as SoA structure) is selected to improve cache locality and the efficiency of vectorized computing. ### Experimental verification The paper verifies the effectiveness of the optimization methods through a series of experiments. The experimental results show that the optimized QED module significantly reduces the computation time while maintaining the accuracy of the simulation results. The specific experimental settings include: - **Test problems**: Simulating the QED cascade process in an extremely strong laser field, using a circularly polarized wave of a plane standing wave as the test field configuration. - **Performance analysis**: Simulations are carried out at different time - steps, and the performance differences between the baseline implementation and the optimized implementation are compared. ### Conclusion The high - performance implementation method proposed in this paper effectively solves the computational efficiency problem in the SFQED process simulation, providing strong support for the numerical simulation of future high - intensity laser experiments.

High-Performance Implementation of the Optimized Event Generator for Strong-Field QED Plasma Simulations

Optimized event generator for strong-field QED simulations within the hi-$χ$ framework

Toward a First-Principles Integrated Simulation of Tokamak Edge Plasmas

SFQEDtoolkit: a high-performance library for the accurate modeling of strong-field QED processes in PIC and Monte Carlo codes

Simulations of spin/polarization-resolved laser-plasma interactions in the nonlinear QED regime

Simulations of spin/polarization-resolved laser–plasma interactions in the nonlinear QED regime

QuickPIC: a highly efficient fully parallelized PIC code for plasma-based acceleration

Study on Extreme Plasma Dynamics by Quantum Electrodynamic Particle-in-Cell Simulations

Extended PIC schemes for physics in ultra-strong laser fields: review and developments

Implementation of a Mesh refinement algorithm into the quasi-static PIC code QuickPIC

Quantum simulations for strong-field QED

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

Mitigation of Numerical Cerenkov Radiation and Instability Using a Hybrid Finite Difference-Fft Maxwell Solver and a Local Charge Conserving Current Deposit

The Plasma Simulation Code: A modern particle-in-cell code with load-balancing and GPU support

Simulations of relativistic quantum plasmas using real-time lattice scalar QED.

Integrating a ponderomotive guiding center algorithm into a quasi-static particle-in-cell code based on azimuthal mode decomposition

Simulating Nonlinear Radiation Diffusion Through Quantum Computing

VSHPIC: A Particle-In-Cell Algorithm Based On Vector Spherical Harmonics Expansion

Particle-in-cell plasma simulation on heterogeneous cluster systems

Accelerating high-order continuum kinetic plasma simulations using multiple GPUs