Insights into the high-energy γ-ray emission of Markarian 501 from extensive multifrequency observations in the Fermi era

We report on the γ-ray activity of the blazar Mrk 501 during the first 480 days of Fermi operation. We find that the average Large Area Telescope (LAT) γ-ray spectrum of Mrk 501 can be well described by a single power-law function with a photon index of 1.78 ± 0.03. While we observe relatively mild flux variations with the Fermi-LAT (within less than a factor of two), we detect remarkable spectral variability where the hardest observed spectral index within the LAT energy range is 1.52 ± 0.14, and the softest one is 2.51 ± 0.20. These unexpected spectral changes do not correlate with the measured flux variations above 0.3 GeV. In this paper, we also present the first results from the 4.5 month long multifrequency campaign (2009 March 15-August 1) on Mrk 501, which included the Very Long Baseline Array (VLBA), Swift, RXTE, MAGIC, and VERITAS, the F-GAMMA, GASP-WEBT, and other collaborations and instruments which provided excellent temporal and energy coverage of the source throughout the entire campaign. The extensive radio to TeV data set from this campaign provides us with the most detailed spectral energy distribution yet collected for this source during its relatively low activity. The average spectral energy distribution of Mrk 501 is well described by the standard one-zone synchrotron self-Compton (SSC) model. In the framework of this model, we find that the dominant emission region is characterized by a size ≲0.1 pc (comparable within a factor of few to the size of the partially resolved VLBA core at 15-43 GHz), and that the total jet power (≃1044 erg s-1) constitutes only a small fraction (∼10-3) of the Eddington luminosity. The energy distribution of the freshly accelerated radiating electrons required to fit the time-averaged data has a broken power-law form in the energy range 0.3 GeV-10 TeV, with spectral indices 2.2 and 2.7 below and above the break energy of 20 GeV. We argue that such a form is consistent with a scenario in which the bulk of the energy dissipation within the dominant emission zone of Mrk 501 is due to relativistic, proton-mediated shocks. We find that the ultrarelativistic electrons and mildly relativistic protons within the blazar zone, if comparable in number, are in approximate energy equipartition, with their energy dominating the jet magnetic field energy by about two orders of magnitude. © 2011. The American Astronomical Society

Observation of a sudden cessation of a very-high-energy gamma-ray flare in PKS 1510-089 with H.E.S.S. and MAGIC in May 2016

The flat spectrum radio quasar (FSRQ) PKS 1510-089 is known for its complex multiwavelength behavior, and is one of only a few FSRQs detected at very high energy (VHE, E >100 GeV) -rays. VHE -ray observations with H.E.S.S. and MAGIC during late May and early June 2016 resulted in the detection of an unprecedented flare, which reveals for the first time VHE -ray intranight variability in this source. While a common variability timescale of 1.5 hr is found, there is a significant deviation near the end of the flare with a timescale of ∼ 20 min marking the cessation of the event. The peak flux is nearly two orders of magnitude above the low-level emission. For the first time, curvature is detected in the VHE -ray spectrum of PKS 1510-089, which is fully explained through absorption by the extragalactic background light. Optical R-band observations with ATOM reveal a counterpart of the -ray flare, even though the detailed flux evolution differs from the VHE lightcurve. Interestingly, a steep flux decrease is observed at the same time as the cessation of the VHE flare. In the high energy (HE, E >100 MeV) -ray band only a moderate flux increase is observed with Fermi-LAT, while the HE -ray spectrum significantly hardens up to a photon index of 1.6. A search for broad-line region (BLR) absorption features in the -ray spectrum indicates that the emission region is located outside of the BLR. Radio VLBI observations reveal a fast moving knot interacting with a standing jet feature around the time of the flare. As the standing feature is located ∼ 50 pc from the black hole, the emission region of the flare may have been located at a significant distance from the black hole. If this correlation is indeed true, VHE rays have been produced far down the jet where turbulent plasma crosses a standing shock.Accepted manuscrip

The Proterozoic evolution of northern Siberian Craton margin: a comparison of U–Pb–Hf signatures from sedimentary units of the Taimyr orogenic belt and the Siberian platform

Identifying the cratonic affinity of Neoproterozoic crust that surrounds the northern margin of the Siberian Craton (SC) is critical for determining its tectonic evolution and placing the Craton in Neoproterozoic supercontinental reconstructions. Integration of new U–Pb–Hf detrital zircon data with regional geological constraints indicates that distinct Neoproterozoic arc-related magmatic belts can be identified within the Taimyr orogen. Sedimentary rocks derived from 970 to 800 Ma arc-related suites reveal abundant Archean and Paleoproterozoic detritus, characteristic of the SC. The 720–600 Ma arc-related zircon population from the younger Cambrian sedimentary rocks is also complemented by an exotic juvenile Mesoproterozoic zircon population and erosional products of older arc-related suites. Nonetheless, numerous evidences imply that both arcs broadly reworked Siberian basement components. We suggest that the early Neoproterozoic (ca. 970–800 Ma) arc system of the Taimyr orogen evolved on the active margin of the SC and probably extended along the periphery of Rodinia into Valhalla orogen of NE Laurentia. We also suggest the late Neoproterozoic (750–550 Ma) arc system could have been part of the Timanian orogen, which linked Siberia and Baltica at the Precambrian/Phanerozoic transition

MAGIC very large zenith angle observations of the Crab Nebula up to 100 TeV

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

Observation of inverse Compton emission from a long γ-ray burst

Long-duration γ-ray bursts (GRBs) originate from ultra-relativistic jets launched from the collapsing cores of dying massive stars. They are characterized by an initial phase of bright and highly variable radiation in the kiloelectronvolt-to-megaelectronvolt band, which is probably produced within the jet and lasts from milliseconds to minutes, known as the prompt emission1,2. Subsequently, the interaction of the jet with the surrounding medium generates shock waves that are responsible for the afterglow emission, which lasts from days to months and occurs over a broad energy range from the radio to the gigaelectronvolt bands1–6. The afterglow emission is generally well explained as synchrotron radiation emitted by electrons accelerated by the external shock7–9. Recently, intense long-lasting emission between 0.2 and 1 teraelectronvolts was observed from GRB 190114C10,11. Here we report multi-frequency observations of GRB 190114C, and study the evolution in time of the GRB emission across 17 orders of magnitude in energy, from 5 × 10−6 to 1012 electronvolts. We find that the broadband spectral energy distribution is double-peaked, with the teraelectronvolt emission constituting a distinct spectral component with power comparable to the synchrotron component. This component is associated with the afterglow and is satisfactorily explained by inverse Compton up-scattering of synchrotron photons by high-energy electrons. We find that the conditions required to account for the observed teraelectronvolt component are typical for GRBs, supporting the possibility that inverse Compton emission is commonly produced in GRBs

Searches for Dark Matter annihilation signatures in the Segue 1 satellite galaxy with the MAGIC-I telescope

We report the results of the observation of the nearby satellite galaxy Segue 1 performed by the MAGIC-I ground-based gamma-ray telescope between November 2008 and March 2009 for a total of 43.2 hours. No significant gamma-ray emission was found above the background. Differential upper limits on the gamma-ray flux are derived assuming various power-law slopes for the possible emission spectrum. Integral upper limits are also calculated for several power-law spectra and for different energy thresholds. The values are of the order of 10^{-11} ph cm^{-2}$ s^{-1} above 100 GeV and 10^{-12} ph cm^{-2} s^{-1} above 200 GeV. Segue 1 is currently considered one of the most interesting targets for indirect dark matter searches. In these terms, the upper limits have been also interpreted in the context of annihilating dark matter particles. For such purpose, we performed a grid scan over a reasonable portion of the parameter space for the minimal SuperGravity model and computed the flux upper limit for each point separately, taking fully into account the peculiar spectral features of each model. We found that in order to match the experimental upper limits with the model predictions, a minimum flux boost of 10^{3} is required, and that the upper limits are quite dependent on the shape of the gamma-ray energy spectrum predicted by each specific model. Finally we compared the upper limits with the predictions of some dark matter models able to explain the PAMELA rise in the positron ratio, finding that Segue 1 data are in tension with the dark matter explanation of the PAMELA spectrum in the case of a dark matter candidate annihilating into tau+tau-. A complete exclusion however is not possible due to the uncertainties in the Segue 1 astrophysical factor.Comment: 26 pages, 10 figures. Matched to published versio

Unravelling the complex behavior of Mrk 421 with simultaneous X-ray and VHE observations during an extreme flaring activity in April 2013

We report on a multi-band variability and correlation study of the TeV blazar Mrk 421 during an exceptional flaring activity observed from 2013 April 11 to 2013 April 19. The study uses, among others, data from GASP-WEBT, Swift, NuSTAR, Fermi-LAT, VERITAS, and MAGIC. The large blazar activity, and the 43 hours of simultaneous NuSTAR and MAGIC/VERITAS observations, permitted variability studies on 15 minute time bins, and over three X-ray bands (3-7 keV, 7-30 keV and 30-80 keV) and three very-high-energy (>0.1 TeV, hereafter VHE) gamma-ray bands (0.2-0.4 TeV, 0.4-0.8 TeV and >0.8 TeV). We detected substantial flux variations on multi-hour and sub-hour timescales in all the X-ray and VHE gamma-ray bands. The characteristics of the sub-hour flux variations are essentially energy-independent, while the multi-hour flux variations can have a strong dependence on the energy of the X-ray and the VHE gamma rays. The three VHE bands and the three X-ray bands are positively correlated with no time-lag, but the strength and the characteristics of the correlation changes substantially over time and across energy bands. Our findings favour multi-zone scenarios for explaining the achromatic/chromatic variability of the fast/slow components of the light curves, as well as the changes in the flux-flux correlation on day-long timescales. We interpret these results within a magnetic reconnection scenario, where the multi-hour flux variations are dominated by the combined emission from various plasmoids of different sizes and velocities, while the sub-hour flux variations are dominated by the emission from a single small plasmoid moving across the magnetic reconnection layer

Multi-messenger characterization of Mrk 501 during historically low X-ray and γ\gamma-ray activity

We study the broadband emission of the TeV blazar Mrk501 using multi-wavelength (MWL) observations from 2017 to 2020 performed with a multitude of instruments, involving, among others, MAGIC, Fermi-LAT, NuSTAR, Swift, GASP-WEBT, and OVRO. During this period, Mrk501 showed an extremely low broadband activity, which may help to unravel its baseline emission. Despite the low activity, significant flux variations are detected at all wavebands, with the highest variations occurring at X-rays and VHE $\gamma$-rays. A significant correlation (>3$\sigma$) between X-rays and VHE $\gamma$-rays is measured, supporting leptonic scenarios to explain the variable parts of the spectral energy distribution (SED), also during low activity states. Extending our data set to 12-years (from 2008 to 2020), we find significant correlations between X-rays and HE $\gamma$-rays, indicating, for the first time, a common physical origin driving the variability between these two bands. We additionally find a correlation between HE $\gamma$-rays and radio, with the radio emission lagging the HE $\gamma$-ray emission by more than 100 days. This is consistent with the $\gamma$-ray emission zone being located upstream of the radio-bright regions of the Mrk501 jet. Furthermore, Mrk501 showed a historically low activity in both X-rays and VHE $\gamma$-rays from mid-2017 to mid-2019 with a stable VHE flux (>2TeV) of 5% the emission of the Crab Nebula. The broadband SED of this 2-year long low-state, the potential baseline emission of Mrk501, can be adequately characterized with a one-zone leptonic model, and with (lepto)-hadronic models that fulfill the neutrino flux constraints from IceCube. We explore the time evolution of the SED towards the historically low-state, revealing that the stable baseline emission may be ascribed to a standing shock, and the variable emission to an additional expanding or traveling shock

Multi-messenger characterization of Mrk 501 during historically low X-ray and γ\gamma-ray activity

International audienceWe study the broadband emission of the TeV blazar Mrk501 using multi-wavelength (MWL) observations from 2017 to 2020 performed with a multitude of instruments, involving, among others, MAGIC, Fermi-LAT, NuSTAR, Swift, GASP-WEBT, and OVRO. During this period, Mrk501 showed an extremely low broadband activity, which may help to unravel its baseline emission. Despite the low activity, significant flux variations are detected at all wavebands, with the highest variations occurring at X-rays and VHE $\gamma$-rays. A significant correlation (>3$\sigma$) between X-rays and VHE $\gamma$-rays is measured, supporting leptonic scenarios to explain the variable parts of the spectral energy distribution (SED), also during low activity states. Extending our data set to 12-years (from 2008 to 2020), we find significant correlations between X-rays and HE $\gamma$-rays, indicating, for the first time, a common physical origin driving the variability between these two bands. We additionally find a correlation between HE $\gamma$-rays and radio, with the radio emission lagging the HE $\gamma$-ray emission by more than 100 days. This is consistent with the $\gamma$-ray emission zone being located upstream of the radio-bright regions of the Mrk501 jet. Furthermore, Mrk501 showed a historically low activity in both X-rays and VHE $\gamma$-rays from mid-2017 to mid-2019 with a stable VHE flux (>2TeV) of 5% the emission of the Crab Nebula. The broadband SED of this 2-year long low-state, the potential baseline emission of Mrk501, can be adequately characterized with a one-zone leptonic model, and with (lepto)-hadronic models that fulfill the neutrino flux constraints from IceCube. We explore the time evolution of the SED towards the historically low-state, revealing that the stable baseline emission may be ascribed to a standing shock, and the variable emission to an additional expanding or traveling shock