Simulation results: General-relativistic radiative cooling in neutron star magnetospheres

This repository contains some of the simulation data presented in the recent article titled "General-relativistic radiative cooling in neutron star magnetospheres" (https://arxiv.org/abs/2603.25949). The data available are from 3D particle pusher simulations, which investigate the effects of radiative cooling in relativistic plasmas in realistic astrophysical conditions such as curved spacetime and non-uniform electromagnetic field geometries. The simulations were performed using a parallelized particle pusher that integrates the equations of motion for a charged particle in a curved spacetime background with a 6th order Runge-Kutta method. By propagating a sufficiently large number of particles, we reconstruct the evolution of the distribution function. The datasets included here contain a random sampling of 1 out of every 100 particles in the simulation, resulting in a dataset that reflects a resolution 100 times lower than the actual simulation. The raw data includes the position in Cartesian coordinates (x,y,z), the parallel and perpendicular momentum components to the magnetic field, and the gyrophase angle (ppa, ppe, pth) for all sampled particles. The datasets included for the reduced parameter simulations extend up to three light-crossing times, while the realistic parameter simulations were run for 10 ring formation times. These datasets represent key aspects of the simulation results discussed in the paper, where the focus is on understanding the effect of realistic astrophysical conditions, such as curved spacetime and non-uniform electromagnetic field geometries, on the radiative cooling process. More details on the simulations and the analysis of these results can be found in the corresponding article.