A new radiation reaction approximation for particle dynamics in the strong field regime

Jérôme Pétri

DOI: https://doi.org/10.1051/0004-6361/202346515

2023-08-02

Abstract:(Abridged) Following particle trajectories in the intense electromagnetic field of a neutron star is prohibited by the large ratio between the cyclotron frequency $\omega_{\rm B}$ and the stellar rotation frequency $\Omega$. No fully kinetic simulations on a macroscopic scale and with realistic field strengths have been performed so far due to the huge computational cost implied by this enormous scale of separation. In this paper, we derive new expressions for the particle velocity subject to strong radiation reaction that are intended to be more accurate than the current state-of-the-art expression in the radiation reaction limit regime, the so-called Aristotelian regime. We shortened the timescale hierarchy by solving the particle equation of motion in the radiation reaction regime, where the Lorentz force is always and immediately balanced by the radiative drag, and including a friction not necessarily opposite to the velocity vector, as derived in the Landau-Lifshitz approximation. Starting from the reduced Landau-Lifshitz equation (i.e. neglecting the field time derivatives), we found expressions for the velocity depending only on the local electromagnetic field configuration and on a new parameter related to the field strength that controls the strength of the radiative damping. The new velocity expressions produce results similar in accuracy to the radiation reaction limit approximation. We therefore do not expect this new method to improve the accuracy of neutron star magnetosphere simulations. The radiation reaction limit is a simple but accurate, robust, and efficient way to follow ultra-relativistic particles in a strong electromagnetic field.

High Energy Astrophysical Phenomena,High Energy Physics - Theory

What problem does this paper attempt to address?

This paper attempts to solve the problem that in a strong electromagnetic field environment (such as near a neutron star), due to the influence of the radiation reaction force on the motion trajectory of charged particles, existing simulation methods cannot accurately track the particle motion. Specifically, the paper focuses on how to overcome the computational challenges brought about by the huge difference in time scales between the gyro frequency and the rotation frequency of the neutron star. This difference makes it extremely difficult to perform full - dynamic simulations on a macroscopic scale, because the range of time scales to be processed is very wide, resulting in extremely high computational costs.

To meet this challenge, the author proposes a new radiation - reaction approximation method, aiming to more accurately describe the particle velocity under strong radiation - reaction conditions. This method is based on the Landau - Lifshitz approximation and takes into account that the radiation resistance is not only in the opposite direction of the particle velocity, but may also contain components along the electric field, magnetic field, and electric drift directions. In this way, the author hopes to provide a more accurate method to simulate the behavior of ultra - relativistic particles in the neutron star magnetosphere, especially under strong electromagnetic field conditions.

The main contributions of the paper are as follows:
1. **Proposing new particle velocity expressions**: These expressions only depend on the local electromagnetic field configuration and a parameter related to the field strength, which controls the intensity of radiation damping.
2. **Verifying the effectiveness of the new method**: By comparing with the existing radiation - reaction limit approximation (RRL), it is found that the results of the new method are basically consistent with RRL at ultra - relativistic speeds (γ ≳ 10), indicating that the new method does not significantly improve in accuracy, but provides another theoretical possibility.
3. **Exploring the applicability and limitations of different approximation methods**: By comparing different approximation methods (such as IRR, VRR, and LLR), the author finds that these methods can give similar results in most cases, but differences may occur under certain specific conditions, such as the particle trajectories near the neutron star surface.

Overall, this paper provides a new theoretical framework for more accurately simulating the motion of charged particles in a strong electromagnetic field environment and lays the foundation for further research.