Radiative turbulent flares in magnetically-dominated plasmas

TL;DR

This study uses kinetic simulations to explore how radiative cooling influences particle acceleration and emission spectra in turbulent, magnetically-dominated plasmas relevant to neutron stars and black holes.

Contribution

It provides new insights into the effects of radiative cooling on turbulence-driven particle acceleration and emission in magnetized plasmas, including the transition between different cooling regimes.

Findings

Weak cooling leads to broad, anisotropic spectra dominated by stochastic acceleration.

Strong cooling suppresses stochastic acceleration, favoring thermal particle emission.

Particle distribution is highly intermittent and cannot be modeled by a universal diffusion coefficient.

Abstract

We perform 2D and 3D kinetic simulations of reconnection-mediated turbulent flares in a magnetized electron-positron plasma, with weak and strong radiative cooling. Such flares can be generated around neutron stars and accreting black holes. We focus on the magnetically-dominated regime where tension of the background magnetic field lines exceeds the plasma rest mass by a factor $\sigma_0 > 1$. In the simulations, turbulence is excited on a macroscopic scale $l_0$, and we observe that it develops by forming thin, dynamic current sheets on various scales. The deposited macroscopic energy dissipates by energizing thermal and nonthermal particles. The particle energy distribution is shaped by impulsive acceleration in reconnecting current sheets, gradual stochastic acceleration, and radiative losses. We parameterize radiative cooling by the ratio $A$ of light-crossing time $l_0/c$ to a cooling timescale, and study the effect of increasing $A$ on the flare. When radiative losses are sufficiently weak, $A<\sigma_0^{-1}$, the produced emission is dominated by stochastically accelerated particles, and the radiative power depends logarithmically on $A$. The resulting flare radiation spectrum is broad and anisotropic. In the strong-cooling regime, $A>\sigma_0^{-1}$, stochastic acceleration is suppressed while impulsive acceleration in the current sheets continues to operate. As $A$ increases further, the emission becomes dominated by thermal particles. Our simulations offer a new tool to study particle acceleration by turbulence, especially at high energies, where cooling competes with acceleration. We find that the particle distribution is influenced by strong intermittency of dissipation, and stochastic acceleration cannot be described by a universal diffusion coefficient.