Abstract
The coincident detection of GW170817 in gravitational waves and electromagnetic radiation spanning the radio to MeV gamma-ray bands provided the first direct evidence that short gamma-ray ...bursts (GRBs) can originate from binary neutron star (BNS) mergers. On the other hand, the properties of short GRBs in high-energy gamma-rays are still poorly constrained, with only ∼20 events detected in the GeV band, and none in the TeV band. GRB 160821B is one of the nearest short GRBs known at
z
= 0.162. Recent analyses of the multiwavelength observational data of its afterglow emission revealed an optical-infrared kilonova component, characteristic of heavy-element nucleosynthesis in a BNS merger. Aiming to better clarify the nature of short GRBs, this burst was automatically followed up with the MAGIC telescopes, starting from 24 s after the burst trigger. Evidence of a gamma-ray signal is found above ∼0.5 TeV at a significance of ∼ 3
σ
during observations that lasted until 4 hr after the burst. Assuming that the observed excess events correspond to gamma-ray emission from GRB 160821B, in conjunction with data at other wavelengths, we investigate its origin in the framework of GRB afterglow models. The simplest interpretation with one-zone models of synchrotron-self-Compton emission from the external forward shock has difficulty accounting for the putative TeV flux. Alternative scenarios are discussed where the TeV emission can be relatively enhanced. The role of future GeV–TeV observations of short GRBs in advancing our understanding of BNS mergers and related topics is briefly addressed.
Context.
Markarian 501 (Mrk 501) is a very high-energy (VHE) gamma-ray blazar located at
z
= 0.034, which is regularly monitored by a wide range of multi-wavelength instruments, from radio to VHE ...gamma rays. During a period of almost two weeks in July 2014, the highest X-ray activity of Mrk 501 was observed in ∼14 years of operation of the
Neil Gehrels Swift
Gamma-ray Burst Observatory.
Aims.
We characterize the broadband variability of Mrk 501 from radio to VHE gamma rays during the most extreme X-ray activity measured in the last 14 years, and evaluate whether it can be interpreted within theoretical scenarios widely used to explain the broadband emission from blazars.
Methods.
The emission of Mrk 501 was measured at radio with Metsähovi, at optical–UV with KVA and
Swift
/UVOT, at X-ray with
Swift
/XRT and
Swift
/BAT, at gamma ray with
Fermi
-LAT, and at VHE gamma rays with the FACT and MAGIC telescopes. The multi-band variability and correlations were quantified, and the broadband spectral energy distributions (SEDs) were compared with predictions from theoretical models.
Results.
The VHE emission of Mrk 501 was found to be elevated during the X-ray outburst, with a gamma-ray flux above 0.15 TeV varying from ∼0.5 to ∼2 times the Crab nebula flux. The X-ray and VHE emission both varied on timescales of 1 day and were found to be correlated. We measured a general increase in the fractional variability with energy, with the VHE variability being twice as large as the X-ray variability. The temporal evolution of the most prominent and variable segments of the SED, characterized on a day-by-day basis from 2014 July 16 to 2014 July 31, is described with a one-zone synchrotron self-Compton model with variations in the break energy of the electron energy distribution (EED), and with some adjustments in the magnetic field strength and spectral shape of the EED. These results suggest that the main flux variations during this extreme X-ray outburst are produced by the acceleration and the cooling of the high-energy electrons. A narrow feature at ∼3 TeV was observed in the VHE spectrum measured on 2014 July 19 (MJD 56857.98), which is the day with the highest X-ray flux (>0.3 keV) measured during the entire
Swift
mission. This feature is inconsistent with the classical analytic functions to describe the measured VHE spectra (power law, log-parabola, and log-parabola with exponential cutoff) at more than 3
σ
. A fit with a log-parabola plus a narrow component is preferred over the fit with a single log-parabola at more than 4
σ
, and a dedicated Monte Carlo simulation estimated the significance of this extra component to be larger than 3
σ
. Under the assumption that this VHE spectral feature is real, we show that it can be reproduced with three distinct theoretical scenarios: (a) a pileup in the EED due to stochastic acceleration; (b) a structured jet with two-SSC emitting regions, with one region dominated by an extremely narrow EED; and (c) an emission from an IC pair cascade induced by electrons accelerated in a magnetospheric vacuum gap, in addition to the SSC emission from a more conventional region along the jet of Mrk 501.
Context. PKS 1510–089 is a flat spectrum radio quasar strongly variable in the optical and GeV range. To date, very high-energy (VHE, > 100 GeV) emission has been observed from this source either ...during long high states of optical and GeV activity or during short flares. Aims. We search for low-state VHE gamma-ray emission from PKS 1510–089. We characterize and model the source in a broadband context, which would provide a baseline over which high states and flares could be better understood. Methods. PKS 1510–089 has been monitored by the MAGIC telescopes since 2012. We use daily binned Fermi-LAT flux measurements of PKS 1510–089 to characterize the GeV emission and select the observation periods of MAGIC during low state of activity. For the selected times we compute the average radio, IR, optical, UV, X-ray, and gamma-ray emission to construct a low-state spectral energy distribution of the source. The broadband emission is modeled within an external Compton scenario with a stationary emission region through which plasma and magnetic fields are flowing. We also perform the emission-model-independent calculations of the maximum absorption in the broad line region (BLR) using two different models. Results. The MAGIC telescopes collected 75 hr of data during times when the Fermi-LAT flux measured above 1 GeV was below 3 × 10−8 cm−2 s−1, which is the threshold adopted for the definition of a low gamma-ray activity state. The data show a strongly significant (9.5σ) VHE gamma-ray emission at the level of (4.27 ± 0.61stat) × 10−12 cm−2 s−1 above 150 GeV, a factor of 80 lower than the highest flare observed so far from this object. Despite the lower flux, the spectral shape is consistent with earlier detections in the VHE band. The broadband emission is compatible with the external Compton scenario assuming a large emission region located beyond the BLR. For the first time the gamma-ray data allow us to place a limit on the location of the emission region during a low gamma-ray state of a FSRQ. For the used model of the BLR, the 95% confidence level on the location of the emission region allows us to place it at a distance > 74% of the outer radius of the BLR.
Aims.
We measure the Crab Nebula
γ
-ray spectral energy distribution in the ~100 TeV energy domain and test the validity of existing leptonic emission models at these high energies.
Methods.
We used ...the novel very large zenith angle observations with the MAGIC telescope system to increase the collection area above 10 TeV. We also developed an auxiliary procedure of monitoring atmospheric transmission in order to assure proper calibration of the accumulated data. This employs recording optical images of the stellar field next to the source position, which provides a better than 10% accuracy for the transmission measurements.
Results.
We demonstrate that MAGIC very large zenith angle observations yield a collection area larger than a square kilometer. In only ~ 56 h of observations, we detect the
γ
-ray emission from the Crab Nebula up to 100 TeV, thus providing the highest energy measurement of this source to date with Imaging Atmospheric Cherenkov Telescopes. Comparing accumulated and archival MAGIC and
Fermi
/LAT data with some of the existing emission models, we find that none of them provides an accurate description of the 1 GeV to 100 TeV
γ
-ray signal.
Abstract
The blazar
TXS 0506+056
got into the spotlight of the astrophysical community in 2017 September, when a high-energy neutrino detected by IceCube (IceCube-170922A) was associated at the 3
σ
...level with a
γ
-ray flare from this source. This multi-messenger photon-neutrino association remains, as per today, the most significant association ever observed. TXS 0506+056 was a poorly studied object before the IceCube-170922A event. To better characterize its broadband emission, we organized a multiwavelength campaign lasting 16 months (2017 November to 2019 February), covering the radio band (Metsähovi, OVRO), the optical/UV (ASAS-SN, KVA, REM, Swift/UVOT), the X-rays (Swift/XRT, NuSTAR), the high-energy
γ
rays (Fermi/LAT), and the very high-energy (VHE)
γ
rays (MAGIC). In
γ
rays, the behavior of the source was significantly different from the behavior in 2017: MAGIC observations show the presence of flaring activity during 2018 December, while the source only shows an excess at the 4
σ
level during the rest of the campaign (74 hr of accumulated exposure); Fermi/LAT observations show several short (on a timescale of days to a week) flares, different from the long-term brightening of 2017. No significant flares are detected at lower energies. The radio light curve shows an increasing flux trend that is not seen in other wavelengths. We model the multiwavelength spectral energy distributions in a lepto-hadronic scenario, in which the hadronic emission emerges as Bethe-Heitler and pion-decay cascade in the X-rays and VHE
γ
rays. According to the model presented here, the 2018 December
γ
-ray flare was connected to a neutrino emission that was too brief and not bright enough to be detected by current neutrino instruments.
Context.
QSO B1420+326 is a blazar classified as a flat-spectrum radio quasar (FSRQ). At the beginning of the year 2020, it was found to be in an enhanced flux state and an extensive multiwavelength ...campaign allowed us to trace the evolution of the flare.
Aims.
We search for very high-energy (VHE) gamma-ray emission from QSO B1420+326 during this flaring state. We aim to characterize and model the broadband emission of the source over different phases of the flare.
Methods.
The source was observed with a number of instruments in radio, near-infrared, optical (including polarimetry and spectroscopy), ultraviolet, X-ray, and gamma-ray bands. We use dedicated optical spectroscopy results to estimate the accretion disk and the dust torus luminosity. We performed spectral energy distribution modeling in the framework of combined synchrotron-self-Compton and external Compton scenario in which the electron energy distribution is partially determined from acceleration and cooling processes.
Results.
During the enhanced state, the flux of both SED components of QSO B1420+326 drastically increased and the peaks were shifted to higher energies. Follow-up observations with the MAGIC telescopes led to the detection of VHE gamma-ray emission from this source, making it one of only a handful of FSRQs known in this energy range. Modeling allows us to constrain the evolution of the magnetic field and electron energy distribution in the emission region. The gamma-ray flare was accompanied by a rotation of the optical polarization vector during a low -polarization state. Also, a new superluminal radio knot contemporaneously appeared in the radio image of the jet. The optical spectroscopy shows a prominent FeII bump with flux evolving together with the continuum emission and a MgII line with varying equivalent width.
Context.
It has become evident that one-zone synchrotron self-Compton models are not always adequate for very high-energy (VHE) gamma-ray-emitting blazars. While two-component models perform better, ...they are difficult to constrain due to the large number of free parameters.
Aims.
In this work, we make a first attempt at taking into account the observational constraints from very long baseline interferometry (VLBI) data, long-term light curves (radio, optical, and X-rays), and optical polarisation to limit the parameter space for a two-component model and test whether or not it can still reproduce the observed spectral energy distribution (SED) of the blazars.
Methods.
We selected five TeV BL Lac objects based on the availability of VHE gamma-ray and optical polarisation data. We collected constraints for the jet parameters from VLBI observations. We evaluated the contributions of the two components to the optical flux by means of decomposition of long-term radio and optical light curves as well as modelling of the optical polarisation variability of the objects. We selected eight epochs for these five objects based on the variability observed at VHE gamma rays, for which we constructed the SEDs that we then modelled with a two-component model.
Results.
We found parameter sets which can reproduce the broadband SED of the sources in the framework of two-component models considering all available observational constraints from VLBI observations. Moreover, the constraints obtained from the long-term behaviour of the sources in the lower energy bands could be used to determine the region where the emission in each band originates. Finally, we attempt to use optical polarisation data to shed new light on the behaviour of the two components in the optical band. Our observationally constrained two-component model allows explanation of the entire SED from radio to VHE with two co-located emission regions.
ABSTRACT
We present a measurement of the extragalactic background light (EBL) based on a joint likelihood analysis of 32 gamma-ray spectra for 12 blazars in the redshift range z = 0.03–0.944, ...obtained by the MAGIC telescopes and Fermi-LAT. The EBL is the part of the diffuse extragalactic radiation spanning the ultraviolet, visible, and infrared bands. Major contributors to the EBL are the light emitted by stars through the history of the Universe, and the fraction of it that was absorbed by dust in galaxies and re-emitted at longer wavelengths.
The EBL can be studied indirectly through its effect on very high energy photons that are emitted by cosmic sources and absorbed via γγ interactions during their propagation across cosmological distances. We obtain estimates of the EBL density in good agreement with state-of-the-art models of the EBL production and evolution. The 1σ upper bounds, including systematic uncertainties, are between 13 per cent and 23 per cent above the nominal EBL density in the models. No anomaly in the expected transparency of the Universe to gamma-rays is observed in any range of optical depth. We also perform a wavelength-resolved EBL determination, which results in a hint of an excess of EBL in the 0.18–0.62 $\mu\mathrm{ m}$ range relative to the studied models, yet compatible with them within systematics.
Although Galactic cosmic rays (protons and nuclei) are widely believed to be mainly accelerated by the winds and supernovae of massive stars, definitive evidence of this origin remains elusive nearly ...a century after their discovery. The active regions of starburst galaxies have exceptionally high rates of star formation, and their large size—more than 50 times the diameter of similar Galactic regions—uniquely enables reliable calorimetric measurements of their potentially high cosmic-ray density. The cosmic rays produced in the formation, life and death of massive stars in these regions are expected to produce diffuse -ray emission through interactions with interstellar gas and radiation. M82, the prototype small starburst galaxy, is predicted to be the brightest starburst galaxy in terms of -ray emission. Here we report the detection of >700-GeV -rays from M82. From these data we determine a cosmic-ray density of 250 eV cm-3 in the starburst core, which is about 500 times the average Galactic density. This links cosmic-ray acceleration to star formation activity, and suggests that supernovae and massive-star winds are the dominant accelerators.