No Arabic abstract
The ubiquitous relativistic jet phenomena associated with black holes play a major role in high and very-high-energy (VHE) astrophysics. In particular, observations have demonstrated that blazars show VHE emission with time-variability from days (in the GeV band) to minutes (in the TeV band), implying very compact emission regions. The real mechanism able to explain the particle acceleration process responsible for this emission is still debated, but magnetic reconnection has been lately discussed as a strong potential candidate and, in some circumstances, as the only possible one. In this work, we present the results of three-dimensional special relativistic magnetohydrodynamic simulations of the development of reconnection events driven by turbulence induced by current-driven kink instability along a relativistic jet. We have performed a systematic identification of all reconnection regions in the system, characterizing their local magnetic field topology and quantifying the reconnection rates. We obtained average rates of $0.051pm0.026$ (in units of the Alfv{e}n speed) which are comparable to the predictions of the theory of turbulence-induced fast reconnection. Detailed statistical analysis also demonstrated that the fast reconnection events follow a log-normal distribution, which is a signature of its turbulent origin. To probe the robustness of our method, we have applied our results to the blazar Mrk 421. Building a synthetic light curve from the integrated power of the magnetic reconnection events, we evaluated the time-variability from a power spectral density analysis, obtaining a good agreement with the observations in the GeV band. This suggests that turbulent fast magnetic reconnection driven by kink instability can be a possible process behind the high energy emission variability phenomena observed in blazars.
Magnetic Towers represent one of two fundamental forms of MHD outflows. Driven by magnetic pressure gradients, these flows have been less well studied than magneto-centrifugally launched jets even though magnetic towers may well be as common. Here we present new results exploring the behavior and evolution of magnetic tower outflows and demonstrate their connection with pulsed power experimental studies and purely hydrodynamic jets which might represent the asymptotic propagation regimes of magneto-centrifugally launched jets. High-resolution AMR MHD simulations (using the AstroBEAR code) provide insights into the underlying physics of magnetic towers and help us constrain models of their propagation. Our simulations have been designed to explore the effects of thermal energy losses and rotation on both tower flows and their hydro counterparts. We find these parameters have significant effects on the stability of magnetic towers, but mild effects on the stability of hydro jets. Current-driven perturbations in the Poynting Flux Dominated (PDF) towers are shown to be amplified in both the cooling and rotating cases. Our studies of the long term evolution of the towers show that the formation of weakly magnetized central jets within the tower are broken up by these instabilities becoming a series of collimated clumps which magnetization properties vary over time. In addition to discussing these results in light of laboratory experiments, we address their relevance to astrophysical observations of young star jets and outflow from highly evolved solar type stars.
Particle acceleration in magnetized relativistic jets still puzzles theorists, specially when one tries to explain the highly variable emission observed in blazar jets or gamma-ray bursts putting severe constraints on current models. In this work we investigate the acceleration of particles injected in a three-dimensional relativistic magnetohydrodynamical jet subject to current driven kink instability (CDKI), which drives turbulence and fast magnetic reconnection. Test protons injected in the nearly stationary snapshots of the jet, experience an exponential acceleration up to a maximum energy. For a background magnetic field of $B sim 0.1$ G, this saturation energy is $sim 10^{16}$ eV, while for $B sim 10$ G it is $sim 10^{18}$ eV. The simulations also reveal a clear association of the accelerated particles with the regions of fast reconnection. In the early stages of the development of the non-linear growth of CDKI in the jet, when there are still no sites of fast reconnection, injected particles are also efficiently accelerated, but by magnetic curvature drift in the wiggling jet spine. However, they have to be injected with an initial energy much larger than that required for particles to accelerate in reconnection sites. Finally, we have also obtained from the simulations an acceleration time due to reconnection with a weak dependence on the particles energy $E$, $t_A propto E^{0.1}$. The energy spectrum of the accelerated particles develops a high energy tail with a power law index $p sim$ -1.2 in the beginning of the acceleration, in agreement with earlier works. Our results provide an appropriate multi-dimensional framework for exploring this process in real systems and explain their complex emission patterns, specially in the very high energy bands and the associated neutrino emission recently detected in some blazars.
We present 3D-MHD AMR simulations of Poynting flux dominated (PFD) jets formed by injection of magnetic energy. We compare their evolution with a hydrodynamic jet which is formed by injecting kinetic energy with the same energy flux than the PFD jets. We predict characteristic emission distributions for each of these jets. Current-driven perturbations in PFD jets are amplified by both cooling and rotation for the regimes studied: Shocks and thermal pressure support are weakened by cooling, making the jets more susceptible to kinking. Rotation amplifies the toroidal magnetic field which also exacerbates the kink instability.
The ejecta composition is an open question in gamma-ray bursts (GRB) physics. Some GRBs possess a quasi-thermal spectral component in the time-resolved spectral analysis, suggesting a hot fireball origin. Others show a featureless non-thermal spectrum known as the Band function, consistent with a synchrotron radiation origin and suggesting that the jet is Poynting-flux-dominated at the central engine and likely in the emission region as well. There are also bursts showing a sub-dominant thermal component and a dominant synchrotron component suggesting a likely hybrid jet composition. Here we report an extraordinarily bright GRB 160625B, simultaneously observed in gamma-rays and optical wavelengths, whose prompt emission consists of three isolated episodes separated by long quiescent intervals, with the durations of each sub-burst being $sim$ 0.8 s, 35 s, and 212 s, respectively. Its high brightness (with isotropic peak luminosity L$_{rm p, iso}sim 4times 10^{53}$ erg/s) allows us to conduct detailed time-resolved spectral analysis in each episode, from precursor to main burst and to extended emission. The spectral properties of the first two sub-bursts are distinctly different, allowing us to observe the transition from thermal to non-thermal radiation between well-separated emission episodes within a single GRB. Such a transition is a clear indication of the change of jet composition from a fireball to a Poynting-flux-dominated jet.
We investigate the distribution of particle acceleration sites, independently of the actual acceleration mechanism, during plasmoid-dominated, relativistic collisionless magnetic reconnection by analyzing the results of a particle-in-cell numerical simulation. The simulation is initiated with Harris-type current layers in pair plasma with no guide magnetic field, negligible radiative losses, no initial perturbation, and using periodic boundary conditions. We find that the plasmoids develop a robust internal structure, with colder dense cores and hotter outer shells, that is recovered after each plasmoid merger on a dynamical time scale. We use spacetime diagrams of the reconnection layers to probe the evolution of plasmoids, and in this context we investigate the individual particle histories for a representative sample of energetic electrons. We distinguish three classes of particle acceleration sites associated with (1) magnetic X-points, (2) regions between merging plasmoids, and (3) the trailing edges of accelerating plasmoids. We evaluate the contribution of each class of acceleration sites to the final energy distribution of energetic electrons -- magnetic X-points dominate at moderate energies, and the regions between merging plasmoids dominate at higher energies. We also identify the dominant acceleration scenarios, in order of decreasing importance -- (1) single acceleration between merging plasmoids, (2) single acceleration at a magnetic X-point, and (3) acceleration at a magnetic X-point followed by acceleration in a plasmoid. Particle acceleration is absent only in the vicinity of stationary plasmoids. The effect of magnetic mirrors due to plasmoid contraction does not appear to be significant in relativistic reconnection.