No Arabic abstract
Particle-in-cell (PIC) simulations have shown that relativistic collisionless magnetic reconnection drives nonthermal particle acceleration (NTPA), potentially explaining high-energy (X-ray/$gamma$-ray) synchrotron and/or inverse Compton (IC) radiation observed from various astrophysical sources. The radiation back-reaction force on radiating particles has been neglected in most of these simulations, even though radiative cooling considerably alters particle dynamics in many astrophysical environments where reconnection may be important. We present a radiative PIC study examining the effects of external IC cooling on the basic dynamics, NTPA, and radiative signatures of relativistic reconnection in pair plasmas. We find that, while the reconnection rate and overall dynamics are basically unchanged, IC cooling significantly influences NTPA: the particle spectra still show a hard power law (index $geq -2$) as in nonradiative reconnection, but transition to a steeper power law that extends to a cooling-dependent cutoff. The steep power-law index fluctuates in time between roughly $-$3 and $-$5. The time-integrated photon spectra display corresponding power laws with indices $approx -0.5$ and $approx -1.1$, similar to those observed in hard X-ray spectra of accreting black holes.
First results are presented from kinetic numerical simulations of relativistic collisionless magnetic reconnection in pair plasma that include radiation reaction from both synchrotron and inverse Compton (IC) processes, motivated by non-thermal high-energy astrophysical sources, including in particular blazars. These simulations are initiated from a configuration known as ABC fields that evolves due to coalescence instability and generates thin current layers in its linear phase. Global radiative efficiencies, instability growth rates, time-dependent radiation spectra, lightcurves, variability statistics and the structure of current layers are investigated for a broad range of initial parameters. We find that the IC radiative signatures are generally similar to the synchrotron signatures. The luminosity ratio of IC to synchrotron spectral components, the Compton dominance, can be modified by more than one order of magnitude with respect to its nominal value. For very short cooling lengths, we find evidence for modification of the temperature profile across the current layers, no systematic compression of plasma density, and very consistent profiles of E.B. We decompose the profiles of E.B with the use of the Vlasov momentum equation, demonstrating a contribution from radiation reaction at the thickness scale consistent with the temperature profile.
Magnetic reconnection, especially in the relativistic regime, provides an efficient mechanism for accelerating relativistic particles and thus offers an attractive physical explanation for nonthermal high-energy emission from various astrophysical sources. I present a simple analytical model that elucidates key physical processes responsible for reconnection-driven relativistic nonthermal particle acceleration (NTPA) in the large-system, plasmoid-dominated regime in two dimensions. The model aims to explain the numerically-observed dependencies of the power-law index $p$ and high-energy cutoff $gamma_c$ of the resulting nonthermal particle energy spectrum $f(gamma)$ on the ambient plasma magnetization $sigma$, and (for $gamma_c$) on the system size $L$. In this self-similar model, energetic particles are continuously accelerated by the out-of-plane reconnection electric field $E_{rm rec}$ until they become magnetized by the reconnected magnetic field and eventually trapped in plasmoids large enough to confine them. The model also includes diffusive Fermi acceleration by particle bouncing off rapidly moving plasmoids. I argue that the balance between electric acceleration and magnetization controls the power-law index, while trapping in plasmoids governs the cutoff, thus tying the particle energy spectrum to the plasmoid distribution.
We present large scale 3D particle-in-cell (PIC) simulations to examine particle energization in magnetic reconnection of relativistic electron-positron (pair) plasmas. The initial configuration is set up as a relativistic Harris equilibrium without a guide field. These simulations are large enough to accommodate a sufficient number of tearing and kink modes. Contrary to the non-relativistic limit, the linear tearing instability is faster than the linear kink instability, at least in our specific parameters. We find that the magnetic energy dissipation is first facilitated by the tearing instability and followed by the secondary kink instability. Particles are mostly energized inside the magnetic islands during the tearing stage due to the spatially varying electric fields produced by the outflows from reconnection. Secondary kink instability leads to additional particle acceleration. Accelerated particles are, however, observed to be thermalized quickly. The large amplitude of the vertical magnetic field resulting from the tearing modes by the secondary kink modes further help thermalizing the non-thermal particles generated from the secondary kink instability. Implications of these results for astrophysics are briefly discussed.
Particle energization in shear flows is invoked to explain non-thermal emission from the boundaries of relativistic astrophysical jets. Yet, the physics of particle injection, i.e., the mechanism that allows thermal particles to participate in shear-driven acceleration, remains unknown. With particle-in-cell simulations, we study the development of Kelvin-Helmholtz (KH) instabilities seeded by the velocity shear between a relativistic magnetically-dominated electron-positron jet and a weakly magnetized electron-ion ambient plasma. We show that, in their nonlinear stages, KH vortices generate kinetic-scale reconnection layers, which efficiently energize the jet particles, thus providing a first-principles mechanism for particle injection into shear-driven acceleration. Our work lends support to spine-sheath models of jet emission - with a fast core/spine surrounded by a slower sheath - and can explain the origin of radio-emitting electrons at the boundaries of relativistic jets.
Using fully kinetic simulations, we study the scaling of the inflow speed of collisionless magnetic reconnection from the non-relativistic to ultra-relativistic limit. In the anti-parallel configuration, the inflow speed increases with the upstream magnetization parameter $sigma$ and approaches the light speed when $sigma > O(100)$, leading to an enhanced reconnection rate. In all regimes, the divergence of pressure tensor is the dominant term responsible for breaking the frozen-in condition at the x-line. The observed scaling agrees well with a simple model that accounts for the Lorentz contraction of the plasma passing through the diffusion region. The results demonstrate that the aspect ratio of the diffusion region remains $sim 0.1$ in both the non-relativistic and relativistic limits.