The physics of gamma-ray bursts & relativistic jets

We provide a comprehensive review of major developments in our understanding of gamma-ray bursts, with particular focus on the discoveries made within the last fifteen years when their true nature was uncovered. We describe the observational properties of photons from the radio to multi-GeV bands, both in the prompt emission and the afterglow phases. Mechanisms for the generation of these photons in GRBs are discussed and confronted with observations to shed light on the physical properties of these explosions, their progenitor stars and the surrounding medium. After presenting observational evidence that a powerful, collimated, jet moving at close to the speed of light is produced in these explosions, we describe our current understanding regarding the generation, acceleration, and dissipation of the jet and compare these properties with jets associated with AGNs and pulsars. We discuss mounting observational evidence that long duration GRBs are produced when massive stars die, and that at least some short duration bursts are associated with old, roughly solar mass, compact stars. The question of whether a black-hole or a strongly magnetized, rapidly rotating neutron star is produced in these explosions is also discussed. We provide a brief summary of what we have learned about relativistic collisionless shocks and particle acceleration from GRB afterglow studies, and discuss the current understanding of radiation mechanism during the prompt emission phase. We discuss theoretical predictions of possible high-energy neutrino emission from GRBs and the current observational constraints. Finally, we discuss how these explosions may be used to study cosmology, e.g. star formation, metal enrichment, reionization history, as well as the formation of first stars and galaxies in the universe.

/pdf/the-physics-of-gamma-ray-bursts-and-relativistic-jets-571z4a98ye.pdf

The Physics of Gamma-Ray Bursts and Relativistic Jets

A high-energy neutrino event detected by IceCube on 22 September 2017 was coincident in direction and time with a gamma-ray flare from the blazar TXS 0506+056. Prompted by this association, we investigated 9.5 years of IceCube neutrino observations to search for excess emission at the position of the blazar. We found an excess of high-energy neutrino events, with respect to atmospheric backgrounds, at that position between September 2014 and March 2015. Allowing for time-variable flux, this constitutes 3.5σ evidence for neutrino emission from the direction of TXS 0506+056, independent of and prior to the 2017 flaring episode. This suggests that blazars are identifiable sources of the high-energy astrophysical neutrino flux.

/pdf/neutrino-emission-from-the-direction-of-the-blazar-txs-0506-15ewkjze3i.pdf

Neutrino emission from the direction of the blazar TXS 0506+056 prior to the IceCube-170922A alert

We describe new implementations of leptonic and hadronic models for the broadband emission from relativistic jets in active galactic nuclei in a temporary steady state. For the leptonic model, a temporary equilibrium between particle injection/acceleration, radiative cooling, and escape from a spherical emission region is evaluated, and the self-consistent radiative output is calculated. For the hadronic model, a temporary equilibrium between particle injection/acceleration, radiative and adiabatic cooling, and escape is evaluated for both primary electrons and protons. A new, semianalytical method to evaluate the radiative output from cascades initiated by internal γγ pair production is presented. We use our codes to fit snapshot spectral energy distributions (SEDs) of a representative set of Fermi-LAT-detected blazars. We find that the leptonic model provides acceptable fits to the SEDs of almost all blazars with parameters close to equipartition between the magnetic field and the relativistic electron population. However, the hard γ-ray spectrum of AO 0235+164, in contrast to the very steep IR-optical-UV continuum, poses a severe problem for the leptonic model. If charge neutrality in leptonic models is provided by cold protons, the kinetic energy carried by the jet should be dominated by protons. We find satisfactory representations of the snapshot SEDs of most blazars in our sample with the hadronic model presented here. However, in the case of two quasars the characteristic break at a few GeV energies cannot be well modeled. All of our hadronic model fits require powers in relativistic protons in the range Lp ~ 1047-1049 erg s–1.

https://iopscience.iop.org/article/10.1088/0004-637X/768/1/54/pdf

Leptonic and Hadronic Modeling of Fermi-detected Blazars

The blazar Markarian 501 (Mrk 501) was observed at energies above 100 GeV with the MAGIC telescope from May through July 2005. The high sensitivity of the instrument enabled the determination of the flux and spectrum of the source on a night-by-night basis. Throughout our observational campaign, the flux from Mrk 501 was found to vary by an order of magnitude, and to be correlated with spectral changes. Intra-night flux variability with flux-doubling times down to 2 minutes was also observed. The strength of variability increased with the energy of the {gamma}-ray photons. The energy spectra were found to harden significantly with increasing flux, and a spectral peak clearly showed up during very active states. The position of the spectral peak seems to be correlated with the source luminosity.

/pdf/variable-very-high-energy-g-ray-emission-from-markarian-501-2h7i02gzhh.pdf

Variable very high energy γ-ray emission from markarian 501

We calculate the spectral energy distribution (SED) of electromagnetic radiation and the spectrum of high energy neutrinos from BL Lac objects in the context of the Synchrotron Proton Blazar Model. In this model, the high energy hump of the SED is due to accelerated protons, while most of the low energy hump is due to synchrotron radiation by co-accelerated electrons. To accelerate protons to sufficiently high energies to produce the high energy hump, rather high magnetic fields are required. Assuming reasonable emission region volumes and Doppler factors, we then find that in low-frequency peaked BL Lacs (LBLs), which have higher luminosities than high-frequency peaked BL Lacs (HBLs), there is a significant contribution to the high frequency hump of the SED from pion photoproduction and subsequent cascading, including synchrotron radiation by muons. In contrast, in HBLs we find that the high frequency hump of the SED is dominated by proton synchrotron radiation. We are able to model the SED of typical LBLs and HBLs, and to model the famous 1997 flare of Markarian 501. We also calculate the expected neutrino output of typical BL Lac objects, and estimate the diffuse neutrino intensity due to all BL Lacs. Because pion photoproduction is inefficient in HBLs, as protons lose energy predominantly by synchrotron radiation, the contribution of LBLs dominates the diffuse neutrino intensity. We suggest that nearby LBLs may well be observable with future high-sensitivity TeV gamma-ray telescopes.

BL Lac Objects in the Synchrotron Proton Blazar Model

(abr.) The spectral energy distribution (SED) of blazars typically has a double-humped appearance usually interpreted in terms of synchrotron self-Compton models. In proton blazar models, the SED is instead explained in terms of acceleration of protons and subsequent cascading. We discuss a variation of the Synchrotron Proton Blazar model, first proposed by Mucke & Protheroe (1999), in which the low energy part of the SED is mainly synchrotron radiation by electrons co-accelerated with protons which produce the high energy part of the SED mainly asproton synchrotron radiation. Using a Monte Carlo/numerical technique to simulate the interactions and subsequent cascading of the accelerated protons, we are able to fit the observed SED of Markarian 501 during the April 1997 flare. We find that the emerging cascade spectra initiated by gamma-rays from $\pi^0$ decay and by $e^\pm$ from $\mu^\pm$ decay turn out to be relatively featureless. Synchrotron radiation produced by $\mu^\pm$ from $\pi^\pm$ decay, and even more importantly by protons, and subsequent synchrotron-pair cascading, is able to reproduce well the high energy part of the SED. For this fit we find that synchrotron radiation by protons dominates the TeV emission, pion photoproduction being less important with the consequence that we predict a lower neutrino flux than in other proton blazar models.

A Proton Synchrotron Blazar Model for Flaring in Markarian~501

A new Monte Carlo program for photohadronic interactions of relativistic nucleons with an ambient photon radiation field is presented. The event generator is designed to fulfil typical astrophysical requirements, but can also be used for radiation and background studies at high energy colliders such as LEP2 and HERA, as well as for simulations of photon induced air showers. 
We consider the full photopion production cross section from the pion production threshold up to high energies. 
It includes resonance excitation and decay, direct single pion production and diffractive and non-diffractive multiparticle production. 
The cross section of each individual process is calculated by fitting experimental data, while the kinematics is determined by the underlying particle production process. We demonstrate that our model is capable of reproducing known accelerator data over a wide energy range.

https://arxiv.org/pdf/astro-ph/9903478.pdf

Monte-Carlo simulations of photohadronic processes in astrophysics

The April 1997 giant flare of Mkn 501 is modelled by a stationary Synchrotron-Proton-Blazar model. Our derived model parameters are consistent with X-ray-to-TeV-data in the flare state and diffusive shock acceleration of electron and proton in a Kolmogorov/Kraichnan turbulence sprectrum. While the emerging pair-synchrotron cascade spectra initiated by photons from neutral pi-decay and electrons from pi+- --> mu+- --> e+- decay turn out to be relatively featureless, muon and proton synchrotron radiation and their cascade radiation produce a double-humped spectral energy distribution. For the present model we find proton synchrotron radiation to dominate the TeV emission, while the contribution from the synchrotron radiation of the pairs, produced by the high energy hump, is only minor.

Modeling the April 1997 flare of Mkn 501

The Synchrotron Proton Blazar model is a promising model to explain high energy emission from gamma-ray loud BL Lac objects like Mkn 421. In contrast to leptonic models, the hadronic explanation of gamma-ray emission predicts ultrahigh energy neutrinos. 
The predicted neutrino spectra from a typical High-energy cutoff BL Lac Object (HBL) and a Low-energy cutoff BL Lac Object (LBL) are presented. We find that cooling due to muon synchrotron radiation causes a cutoff of the neutrino spectrum at $\sim 10^{18}$ eV, with the exception of $\nu_\mu$ from kaon decay which may extend to higher energies if meson production takes place in the secondary resonance region of the cross section. 
The impact of the neutrino output from both source populations to the diffuse neutrino background is discussed.

A. Muecke

Papers

BL Lac Objects in the Synchrotron Proton Blazar Model

A Proton Synchrotron Blazar Model for Flaring in Markarian~501

Monte-Carlo simulations of photohadronic processes in astrophysics

Modeling the April 1997 flare of Mkn 501

Neutrino Emission from HBLs and LBLs