We report on the detection of flaring activity from the Fanaroff-Riley I radio galaxy NGC 1275 in very-high-energy (VHE, E > 100 GeV) gamma rays with the Major Atmospheric Gamma Imaging Cherenkov ...(MAGIC) telescopes. The observations were performed between 2016 September and 2017 February, as part of a monitoring programme. The brightest outburst, with ∼1.5 times the Crab Nebula flux above 100 GeV (C.U.), was observed during the night between 2016 December 31 and 2017 January 1. The flux is fifty times higher than the mean flux previously measured in two observational campaigns between 2009 October and 2010 February and between 2010 August and 2011 February. Significant variability of the day-by-day light curve was measured. The shortest flux-doubling timescale was found to be of (611 ± 101) min. The spectra calculated for this period are harder and show a significant curvature with respect to the ones obtained in the previous campaigns. The combined spectrum of the MAGIC data during the strongest flare state and simultaneous data from the Fermi-LAT around 2017 January 1 follows a power law with an exponential cutoff at the energy (492 ± 35) GeV. We further present simultaneous optical flux density measurements in the R-band obtained with the Kungliga Vetenskaps Akademien (KVA) telescope and investigate the correlation between the optical and gamma-ray emission. Due to possible internal pair-production, the fast flux variability constrains the Doppler factor to values that are inconsistent with a large viewing angle as observed in the radio band. We investigate different scenarios for the explanation of fast gamma-ray variability, namely emission from magnetospheric gaps, relativistic blobs propagating in the jet (mini-jets), or an external cloud (or star) entering the jet. We find that the only plausible model to account for the luminosities here observed would be the production of gamma rays in a magnetospheric gap around the central black hole, only in the eventuality of an enhancement of the magnetic field threading the hole from its equipartition value with the gas pressure in the accretion flow. The observed gamma-ray flare therefore challenges all the discussed models for fast variability of VHE gamma-ray emission in active galactic nuclei.
Aims.
We present a detailed characterisation and theoretical interpretation of the broadband emission of the paradigmatic TeV blazar Mrk 421, with a special focus on the multi-band flux correlations.
...Methods.
The dataset has been collected through an extensive multi-wavelength campaign organised between 2016 December and 2017 June. The instruments involved are MAGIC, FACT,
Fermi
-LAT,
Swift
, GASP-WEBT, OVRO, Medicina, and Metsähovi. Additionally, four deep exposures (several hours long) with simultaneous MAGIC and
NuSTAR
observations allowed a precise measurement of the falling segments of the two spectral components.
Results.
The very-high-energy (VHE;
E
> 100 GeV) gamma rays and X-rays are positively correlated at zero time lag, but the strength and characteristics of the correlation change substantially across the various energy bands probed. The VHE versus X-ray fluxes follow different patterns, partly due to substantial changes in the Compton dominance for a few days without a simultaneous increase in the X-ray flux (i.e., orphan gamma-ray activity). Studying the broadband spectral energy distribution (SED) during the days including
NuSTAR
observations, we show that these changes can be explained within a one-zone leptonic model with a blob that increases its size over time. The peak frequency of the synchrotron bump varies by two orders of magnitude throughout the campaign. Our multi-band correlation study also hints at an anti-correlation between UV-optical and X-ray at a significance higher than 3
σ
. A VHE flare observed on MJD 57788 (2017 February 4) shows gamma-ray variability on multi-hour timescales, with a factor ten increase in the TeV flux but only a moderate increase in the keV flux. The related broadband SED is better described by a two-zone leptonic scenario rather than by a one-zone scenario. We find that the flare can be produced by the appearance of a compact second blob populated by high energetic electrons spanning a narrow range of Lorentz factors, from
γ
′
min
=2×10
4
to
γ
′
max
=6×10
5
.
We report on variability and correlation studies using multiwavelength observations of the blazar Mrk 421 during the month of 2010 February, when an extraordinary flare reaching a level of ∼27 Crab ...Units above 1 TeV was measured in very high energy (VHE) γ-rays with the Very Energetic Radiation Imaging Telescope Array System (VERITAS) observatory. This is the highest flux state for Mrk 421 ever observed in VHE γ-rays. Data are analyzed from a coordinated campaign across multiple instruments, including VHE γ-ray (VERITAS, Major Atmospheric Gamma-ray Imaging Cherenkov), high-energy γ-ray (Fermi-LAT), X-ray (Swift, Rossi X-ray Timing Experiment, MAXI), optical (including the GASP-WEBT collaboration and polarization data), and radio (Metsähovi, Owens Valley Radio Observatory, University of Michigan Radio Astronomy Observatory). Light curves are produced spanning multiple days before and after the peak of the VHE flare, including over several flare "decline" epochs. The main flare statistics allow 2 minute time bins to be constructed in both the VHE and optical bands enabling a cross-correlation analysis that shows evidence for an optical lag of ∼25-55 minutes, the first time-lagged correlation between these bands reported on such short timescales. Limits on the Doppler factor (δ 33) and the size of the emission region ( ) are obtained from the fast variability observed by VERITAS during the main flare. Analysis of 10 minute binned VHE and X-ray data over the decline epochs shows an extraordinary range of behavior in the flux-flux relationship, from linear to quadratic to lack of correlation to anticorrelation. Taken together, these detailed observations of an unprecedented flare seen in Mrk 421 are difficult to explain with the classic single-zone synchrotron self-Compton model.
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.
Spontaneous breaking of Lorentz symmetry at energies on the order of the Planck energy or lower is predicted by many quantum gravity theories, implying non-trivial dispersion relations for the photon ...in vacuum. Consequently, gamma-rays of different energies, emitted simultaneously from astrophysical sources, could accumulate measurable differences in their time of flight until they reach the Earth. Such tests have been carried out in the past using fast variations of gamma-ray flux from pulsars, and more recently from active galactic nuclei and gamma-ray bursts. We present new constraints studying the gamma-ray emission of the galactic Crab Pulsar, recently observed up to TeV energies by the Major Atmospheric Gamma-ray Imaging Cherenkov (MAGIC) collaboration. A profile likelihood analysis of pulsar events reconstructed for energies above 400 GeV finds no significant variation in arrival time as their energy increases. Ninety-five percent CL limits are obtained on the effective Lorentz invariance violating energy scale at the level of ( ) for a linear, and ( ) for a quadratic scenario, for the subluminal and the superluminal cases, respectively. A substantial part of this study is dedicated to calibration of the test statistic, with respect to bias and coverage properties. Moreover, the limits take into account systematic uncertainties, which are found to worsen the statistical limits by about 36%-42%. Our constraints would have been much more stringent if the intrinsic pulse shape of the pulsar between 200 GeV and 400 GeV was understood in sufficient detail and allowed inclusion of events well below 400 GeV.
Context. The BL Lac object S5 0716+714, a highly variable blazar, underwent an impressive outburst in January 2015 (Phase A), followed by minor activity in February (Phase B). The MAGIC observations ...were triggered by the optical flux observed in Phase A, corresponding to the brightest ever reported state of the source in the R-band. Aims.The comprehensive dataset collected is investigated in order to shed light on the mechanism of the broadband emission. Methods. Multi-wavelength light curves have been studied together with the broadband spectral energy distributions (SEDs). The sample includes data from Effelsberg, OVRO, Metsähovi, VLBI, CARMA, IRAM, SMA, Swift-UVOT, KVA, Tuorla, Steward, RINGO3, KANATA, AZT-8+ST7, Perkins, LX-200, Swift-XRT, NuSTAR, Fermi-LAT and MAGIC. Results. The flaring state of Phase A was detected in all the energy bands, providing for the first time a multi-wavelength sample of simultaneous data from the radio band to the very-high-energy (VHE, E > 100 GeV). In the constructed SED, the Swift-XRT+NuSTAR data constrain the transition between the synchrotron and inverse Compton components very accurately, while the second peak is constrained from 0.1 GeV to 600 GeV by Fermi+MAGIC data. The broadband SED cannot be described with a one-zone synchrotron self-Compton model as it severely underestimates the optical flux in order to reproduce the X-ray to γ-ray data. Instead we use a two-zone model. The electric vector position angle (EVPA) shows an unprecedented fast rotation. An estimation of the redshift of the source by combined high-energy (HE, 0.1 GeV < E < 100 GeV) and VHE data provides a value of z = 0.31 ± 0.02stats ± 0.05sys, confirming the literature value. Conclusions. The data show the VHE emission originating in the entrance and exit of a superluminal knot in and out of a recollimation shock in the inner jet. A shock–shock interaction in the jet seems responsible for the observed flares and EVPA swing. This scenario is also consistent with the SED modeling.
Abstract
We report on a long-lasting, elevated gamma-ray flux state from VER J0521+211 observed by VERITAS, MAGIC, and Fermi-LAT in 2013 and 2014. The peak integral flux above 200 GeV measured with ...the nightly binned light curve is (8.8 ± 0.4) × 10
−7
photons m
−2
s
−1
, or ∼37% of the Crab Nebula flux. Multiwavelength observations from X-ray, UV, and optical instruments are also presented. A moderate correlation between the X-ray and TeV gamma-ray fluxes was observed, and the X-ray spectrum appeared harder when the flux was higher. Using the gamma-ray spectrum and four models of the extragalactic background light (EBL), a conservative 95% confidence upper limit on the redshift of the source was found to be
z
≤ 0.31. Unlike the gamma-ray and X-ray bands, the optical flux did not increase significantly during the studied period compared to the archival low-state flux. The spectral variability from optical to X-ray bands suggests that the synchrotron peak of the spectral energy distribution (SED) may become broader during flaring states, which can be adequately described with a one-zone synchrotron self-Compton model varying the high-energy end of the underlying particle spectrum. The synchrotron peak frequency of the SED and the radio morphology of the jet from the MOJAVE program are consistent with the source being an intermediate-frequency-peaked BL Lac object.
Context.
The BL Lac object 1ES 0647+250 is one of the few distant
γ
-ray emitting blazars detected at very high energies (VHEs; ≳100 GeV) during a non-flaring state. It was detected with the MAGIC ...telescopes during a period of low activity in the years 2009−2011 as well as during three flaring activities in the years 2014, 2019, and 2020, with the highest VHE flux in the last epoch. An extensive multi-instrument data set was collected as part of several coordinated observing campaigns over these years.
Aims.
We aim to characterise the long-term multi-band flux variability of 1ES 0647+250, as well as its broadband spectral energy distribution (SED) during four distinct activity states selected in four different epochs, in order to constrain the physical parameters of the blazar emission region under certain assumptions.
Methods.
We evaluated the variability and correlation of the emission in the different energy bands with the fractional variability and the Z-transformed discrete correlation function, as well as its spectral evolution in X-rays and
γ
rays. Owing to the controversy in the redshift measurements of 1ES 0647+250 reported in the literature, we also estimated its distance in an indirect manner through a comparison of the GeV and TeV spectra from simultaneous observations with
Fermi
-LAT and MAGIC during the strongest flaring activity detected to date. Moreover, we interpret the SEDs from the four distinct activity states within the framework of one-component and two-component leptonic models, proposing specific scenarios that are able to reproduce the available multi-instrument data.
Results.
We find significant long-term variability, especially in X-rays and VHE
γ
rays. Furthermore, significant (3−4
σ
) correlations were found between the radio, optical, and high-energy (HE)
γ
-ray fluxes, with the radio emission delayed by about ∼400 days with respect to the optical and
γ
-ray bands. The spectral analysis reveals a harder-when-brighter trend during the non-flaring state in the X-ray domain. However, no clear patterns were observed for either the enhanced states or the HE (30 MeV <
E
< 100 GeV) and VHE
γ
-ray emission of the source. The indirect estimation of the redshift yielded a value of
z
= 0.45 ± 0.05, which is compatible with some of the values reported in the literature. The SEDs related to the low-activity state and the three flaring states of 1ES 0647+250 can be described reasonably well with the both one-component and two-component leptonic scenarios. However, the long-term correlations indicate the need for an additional radio-producing region located about 3.6 pc downstream from the gamma-ray producing region.
Abstract
PSR J0218+4232 is one of the most energetic millisecond pulsars known and has long been considered as one of the best candidates for very high-energy (VHE; >100 GeV)
γ
-ray emission. Using ...11.5 yr of Fermi Large Area Telescope (LAT) data between 100 MeV and 870 GeV, and ∼90 hr of Major Atmospheric Gamma Imaging Cherenkov (MAGIC) observations in the 20 GeV to 20 TeV range, we searched for the highest energy
γ
-ray emission from PSR J0218+4232. Based on the analysis of the LAT data, we find evidence for pulsed emission above 25 GeV, but see no evidence for emission above 100 GeV (VHE) with MAGIC. We present the results of searches for
γ
-ray emission, along with theoretical modeling, to interpret the lack of VHE emission. We conclude that, based on the experimental observations and theoretical modeling, it will remain extremely challenging to detect VHE emission from PSR J0218+4232 with the current generation of Imaging Atmospheric Cherenkov Telescopes, and maybe even with future ones, such as the Cherenkov Telescope Array.