Maria Petropoulou

Relativistic Jets Shine through Shocks or Magnetic Reconnection

Observations of gamma-ray-bursts and jets from active galactic nuclei reveal that the jet flow is characterized by a high radiative e ffi ciency and that the dissipative mechanism must be a powerful accelerator of non-thermal particles. Shocks and magnetic reconnection have long been considered as possible candidates for powering the jet emission. Recent progress via fully-kinetic particle-in-cell simulations allows us to revisit this issue on firm physical grounds. We show that shock models are unlikely to account for the jet emission. In fact, when shocks are effi cient at dissipating energy, they typically do not accelera te particles far beyond the thermal energy, and vice versa. In contrast, we show that magnetic reconnection can deposit more than 50% of the dissipated energy into non-thermal leptons as long as the energy density of the magnetic field in the bulk flow is larger tha n the rest mass energy density. The emitting region, i.e., the reconnection downstream, is characterized by a rough energy equipartition between magnetic fields and radiating partic les, which naturally accounts for a commonly observed property of blazar jets.

Photohadronic origin of

The recent IceCube discovery of 0.1-1 PeV neutrinos of astrophysical origin opens up a new era for high-energy astrophysics. Although there are various astrophysical candidate sources, a firm association of the detected neutrinos with one (or more) of them is still lacking. A recent analysis of plausible astrophysical counterparts within the error circles of IceCube events showed that likely counterparts for nine of the IceCube neutrinos include mostly BL Lacs, among which Mrk 421. Motivated by this result and a previous independent analysis on the neutrino emission from Mrk 421, we test the BL Lac-neutrino connection in the context of a specific theoretical model for BL Lac emission. We model the spectral energy distribution (SED) of the BL Lacs selected as counterparts of the IceCube neutrinos using a one-zone leptohadronic model and mostly nearly simultaneous data. The neutrino flux for each BL Lac is self-consistently calculated, using photon and proton distributions specifically derived for every individual source. We find that the SEDs of the sample, although different in shape and flux, are all well fitted by the model using reasonable parameter values. Moreover, the model-predicted neutrino flux and energy for these sources are of the same order of magnitude as those of the IceCube neutrinos. In two cases, namely Mrk 421 and H 1914-194, we find a suggestively good agreement between the model prediction and the detected neutrino flux. Our predictions for all the BL Lacs of the sample are in the range to be confirmed or disputed by IceCube in the next few years of data sampling.

\boldsymbol {\gamma }

Blobs, or quasi-spherical emission regions containing relativistic particles and magnetic fields, are often assumed ad hoc in emission models of relativistic astrophysical jets, yet their physical origin is still not well understood. Here, we employ a suite of large-scale two-dimensional particle-in-cell simulations in electron-positron plasmas to demonstrate that relativistic magnetic reconnection can naturally account for the formation of quasi-spherical plasmoids filled with high-energy particles and magnetic fields. Our simulations extend to unprecedentedly long temporal and spatial scales, so we can capture the asymptotic physics independently of the initial setup. We characterize the properties of the plasmoids that are continuously generated as a self-consistent by-product of the reconnection process: they are in rough energy equipartition between particles and magnetic fields; the upper energy cutoff of the plasmoid particle spectrum is proportional to the plasmoid width w, corresponding to a Larmor radius ~0.2 w; the plasmoids grow in size at ~0.1 of the speed of light, with most of the growth happening while they are still non-relativistic (first they grow); their growth is suppressed once they get accelerated to relativistic speeds by the field line tension, up to the Alfven speed (then they go). The largest plasmoids, whose typical recurrence interval is ~2.5 L/c, reach a characteristic size w ~ 0.2 L independently of the system length L, they have nearly isotropic particle distributions and they contain the highest energy particles, whose Larmor radius is ~0.03 L. The latter can be regarded as the Hillas criterion for relativistic reconnection. We briefly discuss the implications of our results for the high-energy emission from relativistic jets and pulsar winds.

-ray BL Lac emission: implications for IceCube neutrinos

We investigate the origin of high-energy emission in blazars within the context of the leptohadronic one-zone model. We find that γ -ray emission can be attributed to synchrotron radiation either from protons or from secondary leptons produced via photohadronic processes. These possibilities imply differences not only in the spectral energy distribution (SED) but also in the variability signatures, especially in the X- and γ -ray regime. Thus, the temporal behaviour of each leptohadronic scenario can be used to probe the particle population responsible for the high-energy emission as it can give extra information not available by spectral fits. In this work, we apply these ideas to the non-thermal emission of Mrk 421, which is one of the best monitored TeV blazars. We focus on the observations of 2001 March, since during that period Mrk 421 showed multiple flares that have been observed in detail both in X-rays and γ -rays. First, we obtain pre-flaring fits to the SED using the different types of leptohadronic scenarios. Then, we introduce random-walk-type, small-amplitude variations on the injection compactness or on the maximum energy of radiating particles and follow the subsequent response of the radiated photon spectrum. For each leptohadronic scenario, we calculate the X-ray and γ -ray fluxes and investigate their possible correlation. Whenever the ‘input’ variations lead, apart from flux variability, also to spectral variability, we present the resulting relations between the spectral index and the flux, both in X-rays and γ -rays. We find that proton synchrotron models are favoured energetically but require fine tuning between electron and proton parameters to reproduce the observed quadratic behaviour between X-rays and TeV γ -rays. On the other hand, models based on pion decay can reproduce this behaviour in a much more natural way.

Plasmoids in relativistic reconnection, from birth to adulthood: first they grow, then they go

Powerful flares from blazars with short (

Mrk 421 as a case study for TeV and X-ray variability in leptohadronic models

\sim

Blazar flares powered by plasmoids in relativistic reconnection

min) variability timescales are challenging for current models of blazar emission. Here, we present a physically motivated ab initio model for blazar flares based on the results of recent particle-in-cell (PIC) simulations of relativistic magnetic reconnection. PIC simulations demonstrate that quasi-spherical plasmoids filled with high-energy particles and magnetic fields are a self-consistent by-product of the reconnection process. By coupling our PIC-based results (i.e., plasmoid growth, acceleration profile, particle and magnetic content) with a kinetic equation for the evolution of the electron distribution function we demonstrate that relativistic reconnection in blazar jets can produce powerful flares whose temporal and spectral properties are consistent with the observations. In particular, our model predicts correlated synchrotron and synchrotron self-Compton flares of duration of several hours--days powered by the largest and slowest moving plasmoids that form in the reconnection layer. Smaller and faster plasmoids produce flares of sub-hour duration with higher peak luminosities than those powered by the largest plasmoids. Yet, the observed fluence in both types of flares is similar. Multiple flares with a range of flux-doubling timescales (minutes to several hours) observed over a longer period of flaring activity (days or longer) may be used as a probe of the reconnection layers orientation and the jets magnetization. Our model shows that blazar flares are naturally expected as a result of magnetic reconnection in a magnetically-dominated jet.

A combined radio and GeV γ-ray view of the 2012 and 2013 flares of Mrk 421

In 2012 Markarian 421 underwent the largest flare ever observed in this blazar at radio frequencies. In the present study, we start exploring this unique event and compare it to a less extreme event in 2013. We use 15 GHz radio data obtained with the Owens Valley Radio Observatory 40-m telescope, 95 GHz millimetre data from the Combined Array for Research in Millimeter-Wave Astronomy, and GeV γ-ray data from the Fermi Gamma-ray Space Telescope. The radio light curves during the flaring periods in 2012 and 2013 have very different appearances, in both shape and peak flux density. Assuming that the radio and γ-ray flares are physically connected, we attempt to model the most prominent sub-flares of the 2012 and 2013 activity periods by using the simplest possible theoretical framework. We first fit a one-zone synchrotron self-Compton (SSC) model to the less extreme 2013 flare and estimate parameters describing the emission region. We then model the major γ-ray and radio flares of 2012 using the same framework. The 2012 γ-ray flare shows two distinct spikes of similar amplitude, so we examine scenarios associating the radio flare with each spike in turn. In the first scenario, we cannot explain the sharp radio flare with a simple SSC model, but we can accommodate this by adding plausible time variations to the Doppler beaming factor. In the second scenario, a varying Doppler factor is not needed, but the SSC model parameters require fine-tuning. Both alternatives indicate that the sharp radio flare, if physically connected to the preceding γ-ray flares, can be reproduced only for a very specific choice of parameters.

Temporal signatures of leptohadronic feedback mechanisms in compact sources

The hadronic model of active galactic nuclei and other compact high-energy astrophysical sources assumes that ultra-relativistic protons, electron–positron pairs and photons interact via various hadronic and electromagnetic processes inside a magnetized volume, producing the multiwavelength spectra observed from these sources. A less studied property of such systems is that they can exhibit a variety of temporal behaviours due to the operation of different feedback mechanisms. We investigate the effects of one possible feedback loop, where γ-rays produced by photopion processes are being quenched whenever their compactness increases above a critical level. This causes a spontaneous creation of soft photons in the system that result in further proton cooling and more production of γ-rays, thus making the loop operate. We perform an analytical study of a simplified set of equations describing the system, in order to investigate the connection of its temporal behaviour with key physical parameters. We also perform numerical integration of the full set of kinetic equations verifying not only our analytical results but also those of previous numerical studies. We find that once the system becomes ‘supercritical’, it can exhibit either a periodic behaviour or a damped oscillatory one leading to a steady state. We briefly point out possible implications of such a supercriticality on the parameter values used in active galactic nuclei spectral modelling, through an indicative fitting of the VHE emission of blazar 3C 279.

Implications of a PeV neutrino spectral cut-off in gamma-ray burst models

The recent discovery of extragalactic PeV neutrinos opens a new window to the exploration of cosmic-ray accelerators. The observed PeV neutrino flux is close to the Waxman-Bahcall upper bound implying that gamma-ray bursts (GRBs) may be the source of ultra-high energy cosmic rays (UHECRs). Starting with the assumption of the GRB-UHECR connection, we show using both analytical estimates and numerical simulations that the observed neutrinos can originate at the jet as a result of photopion interactions with the following implications: the neutrino spectra are predicted to have a cutoff at energy