Search for: All records

We present the results of 3D particle-in-cell simulations that explore relativistic magnetic reconnection in pair plasma with strong synchrotron cooling and a small mass fraction of nonradiating ions. Our results demonstrate that the structure of the current sheet is highly sensitive to the dynamic efficiency of radiative cooling. Specifically, stronger cooling leads to more significant compression of the plasma and magnetic field within the plasmoids. We demonstrate that ions can be efficiently accelerated to energies exceeding the plasma magnetization parameter, ≫σ, and form a hard power-law energy distribution,f_i∝γ⁻¹. This conclusion implies a highly efficient proton acceleration in the magnetospheres of young pulsars. Conversely, the energies of pairs are limited to eitherσin the strong cooling regime or the radiation burnoff limit,γ_syn, when cooling is weak. We find that the high-energy radiation from pairs above the synchrotron burnoff limit,ε_c≈ 16 MeV, is only efficiently produced in the strong cooling regime,γ_syn<σ. In this regime, we find that the spectral cutoff scales asε_cut≈ε_c(σ/γ_syn) and the highest energy photons are beamed along the direction of the upstream magnetic field, consistent with the phenomenological models of gamma-ray emission from young pulsars. Furthermore, our results place constraints on the reconnection-driven models of gamma-ray flares in the Crab Nebula.

We present high-resolution 2D and 3D simulations of magnetized decaying turbulence in relativistic, resistive magnetohydrodynamics. The simulations show dynamic formation of large-scale intermittent long-lived current sheets being disrupted into plasmoid chains by the tearing instability. These current sheets are locations of enhanced magnetic-field dissipation and heating of the plasma. We find magnetic energy spectra ∝k^−3/2, together with strongly pronounced dynamic alignment of Elsässer fields and of velocity and magnetic fields, for strong guide-field turbulence, whereas we retrieve spectra ∝k^−5/3for the case of a weak guide-field.