Abstract Simulations of relativistic plasmas traditionally focus on the dynamics of two-species mixtures of charged particles under the influence of external magnetic fields and those generated by particle currents. However, the extreme conditions of astrophysical plasmas near compact objects, such as black holes and neutron stars, are often characterized by mixtures of electrons, protons, and positrons, whose dynamics can differ significantly, because of the considerable mass contrast. We present the first two-dimensional particle-in-cell simulations of relativistic turbulence and magnetic reconnection in a three-species plasma, varying the relative abundances of electrons, protons, and positrons, while employing realistic mass ratios to achieve unprecedented accuracy. We find that turbulence leads to the formation of magnetic islands, current sheets, and plasmoids. Reconnection occurs between these structures, with plasma composition playing a key role in determining the number of reconnection sites and their energy conversion efficiency. In particular, as the proton fraction increases, very small-scale features of the turbulence are washed out, while global dissipative effects are amplified. Finally, using a novel generalization of Ohm’s law for a relativistic multispecies plasma, we find that the reconnection rate is primarily governed by the electric fields associated with the divergence of the positron and electron pressure tensors. These results provide new insights into dissipation and particle acceleration in turbulent relativistic plasmas, such as those near black holes and neutron stars, and can be used to interpret their high-energy emission and phenomenology.
Imbrogno et al. (Mon,) studied this question.