MAGIC reveals a complex morphology within the unidentified gamma-ray source HESS J1857+026

Abstract not availableJ. Aleksić ... V. Stamatescu ... et.al

MAGIC observations of MWC 656, the only known Be/BH system

Context. MWC 656 has recently been established as the first observationally detected high-mass X-ray binary system containing a Be star and a black hole (BH). The system has been associated with a gamma-ray flaring event detected by the AGILE satellite in July 2010. Aims. Our aim is to evaluate whether the MWC 656 gamma-ray emission extends to very high energy (VHE > 100 GeV) gamma rays. Methods. We observed MWC 656 with the MAGIC telescopes for ~23 h during two observation periods: between May and June 2012, and in June 2013. During the last period, observations were performed contemporaneously with X-ray (XMM-Newton) and optical (STELLA) instruments. Results. We did not detect the MWC 656 binary system at TeV energies with the MAGIC telescopes in either of the two campaigns. Upper limits (ULs) to the integral flux above 300 GeV have been set, as well as differential ULs at a level of ~5% of the Crab nebula flux. The results obtained from the MAGIC observations do not support persistent emission of VHE gamma rays from this system at a level of 2.4% the Crab flux.J. Aleksić, S. Ansoldi, L. A. Antonelli, P. Antoranz, A. Babic ... Sabrina Einecke ... et al

Very-high-energy gamma-ray observations of the Type Ia Supernova SN 2014J with the MAGIC telescopes

Context. In this work we present data from observations with the MAGIC telescopes of SN 2014J detected on January 21 2014, the closest Type Ia supernova since Imaging Air Cherenkov Telescopes started to operate. Aims. We aim to probe the possibility of very-high-energy (VHE; E ≥ 100 GeV) gamma rays produced in the early stages of Type Ia supernova explosions. Methods. We performed follow-up observations after this supernova (SN) explosion for five days, between January 27 and February 2 2014. We searched for gamma-ray signals in the energy range between 100 GeV and several TeV from the location of SN 2014J using data from a total of ~5.5 h of observations. Prospects for observing gamma rays of hadronic origin from SN 2014J in the near future are also being addressed. Results. No significant excess was detected from the direction of SN 2014J. Upper limits at 95% confidence level on the integral flux, assuming a power-law spectrum, dF/dE ∝ E− Γ, with a spectral index of Γ = 2.6, for energies higher than 300 GeV and 700 GeV, are established at 1.3 × 10-12 and 4.1 × 10-13 photons cm-2 s-1, respectively. Conclusions. For the first time, upper limits on the VHE emission of a Type Ia supernova are established. The energy fraction isotropically emitted into TeV gamma rays during the first ~ 10 days after the supernova explosion for energies greater than 300 GeV is limited to 10-6 of the total available energy budget (~1051 erg). Within the assumed theoretical scenario, the MAGIC upper limits on the VHE emission suggest that SN 2014J will not be detectable in the future by any current or planned generation of Imaging Atmospheric Cherenkov Telescopes.M. L. Ahnen, S. Ansoldi, L. A. Antonelli, P. Antoranz, C. Arcaro ... Sabrina Einecke ... et al

Insights into the emission of the blazar 1ES 1011+496 through unprecedented broadband observations during 2011 and 2012

Context. 1ES 1011+496 (z = 0.212) was discovered in very high-energy (VHE, E> 100 GeV) γ rays with MAGIC in 2007. The absence of simultaneous data at lower energies led to an incomplete characterization of the broadband spectral energy distribution (SED). Aims. We study the source properties and the emission mechanisms, probing whether a simple one-zone synchrotron self-Compton (SSC) scenario is able to explain the observed broadband spectrum. Methods. We analyzed data in the range from VHE to radio data from 2011 and 2012 collected by MAGIC, Fermi-LAT, Swift, KVA, OVRO, and Metsähovi in addition to optical polarimetry data and radio maps from the Liverpool Telescope and MOJAVE. Results. The VHE spectrum was fit with a simple power law with a photon index of 3.69 ± 0.22 and a flux above 150 GeV of (1.46 ± 0.16) × 10-11 ph cm-2 s-1. The source 1ES 1011+496 was found to be in a generally quiescent state at all observed wavelengths, showing only moderate variability from radio to X-rays. A low degree of polarization of less than 10% was measured in optical, while some bright features polarized up to 60% were observed in the radio jet. A similar trend in the rotation of the electric vector position angle was found in optical and radio. The radio maps indicated a superluminal motion of 1.8 ± 0.4 c, which is the highest speed statistically significant measured so far in a high-frequency-peaked BL Lac. Conclusions. For the first time, the high-energy bump in the broadband SED of 1ES 1011+496 could be fully characterized from 0.1 GeV to 1 TeV, which permitted a more reliable interpretation within the one-zone SSC scenario. The polarimetry data suggest that at least part of the optical emission has its origin in some of the bright radio features, while the low polarization in optical might be due to the contribution of parts of the radio jet with different orientations of the magnetic field with respect to the optical emission.J. Aleksić ... S. Einecke ... et al

Long-term multi-wavelength variability and correlation study of Markarian 421 from 2007 to 2009

Aims: We study the multi-band variability and correlations of the TeV blazar Mrk 421 on year timescales, which can bring additional insight on the processes responsible for its broadband emission. Methods: We observed Mrk 421 in the very high energy (VHE) γ-ray range with the Cherenkov telescope MAGIC-I from March 2007 to June 2009 for a total of 96 h of effective time after quality cuts. The VHE flux variability is quantified using several methods, including the Bayesian Block algorithm, which is applied to data from Cherenkov telescopes here for the first time. The 2.3 yr long MAGIC light curve is complemented with data from the Swift/BAT and RXTE/ASM satellites and the KVA, GASP-WEBT, OVRO, and Metsähovi telescopes from February 2007 to July 2009, allowing for an excellent characterisation of the multi-band variability and correlations over year timescales. Results: Mrk 421 was found in different γ-ray emission states during the 2.3 yr long observation period: The flux above 400 GeV spans from the minimum nightly value of (1:3 ± 0:4) × 10⁻¹¹ cm⁻² s⁻¹ to the maximum flux, that is about 24 times higher, at (3:1 ± 0:1) × 10⁻¹⁰ cm⁻² s⁻¹. Flares and different levels of variability in the γ-ray light curve could be identified with the Bayesian Block algorithm. The same behaviour of a quiet and active emission was found in the X-ray light curves measured by Swift/BAT and the RXTE/ASM, with a direct correlation in time. The behaviour of the optical light curve of GASP-WEBT and the radio light curves by OVRO and Metsähovi are different as they show no coincident features with the higher energetic light curves and a less variable emission. Overall, the fractional variability increases with energy. The comparable variability in the X-ray and VHE bands and their direct correlation during both high- and low-activity periods spanning many months show that the electron populations radiating the X-ray and γ-ray photons are either the same, as expected in the synchrotron-self-Compton mechanism, or at least strongly correlated, as expected in electromagnetic cascades.M. L. Ahnen ... S. Einecke ... et al

Very high-energy γ -ray observations of novae and dwarf novae with the MAGIC telescopes

Context. In the last five years the Fermi Large Area Telescope (LAT) instrument detected GeV γ-ray emission from five novae. The GeV emission can be interpreted in terms of an inverse Compton process of electrons accelerated in a shock. In this case it is expected that protons in the same conditions can be accelerated to much higher energies. Consequently they may produce a second component in the γ-ray spectrum at TeV energies. Aims. We aim to explore the very high-energy domain to search for γ-ray emission above 50 GeV and to shed light on the acceleration process of leptons and hadrons in nova explosions. Methods. We have performed observations, with the MAGIC telescopes of the classical nova V339 Del shortly after the 2013 outburst; optical and subsequent GeV γ-ray detections triggered the MAGIC observations. We also briefly report on VHE observations of the symbiotic nova YY Her and the dwarf nova ASASSN-13ax. We complement the TeV MAGIC observations with the analysis of contemporaneous Fermi-LAT data of the sources. The TeV and GeV observations are compared in order to evaluate the acceleration parameters for leptons and hadrons. Results. No significant TeV emission was found from the studied sources. We computed upper limits on the spectrum and night-by-night flux. The combined GeV and TeV observations of V339 Del limit the ratio of proton to electron luminosities to Lp ≲ 0.15 Le.M. L. Ahnen, S. Ansoldi, L. A. Antonelli, P. Antoranz, A. Babic ... Sabrina Einecke ... et al

Unprecedented study of the broadband emission of Mrk 421 during flaring activity in March 2010

Context. Because of its proximity, Mrk 421 is one of the best sources on which to study the nature of BL Lac objects. Its proximity allows us to characterize its broadband spectral energy distribution (SED). Aims. The goal is to better understand the mechanisms responsible for the broadband emission and the temporal evolution of Mrk 421. These mechanisms may also apply to more distant blazars that cannot be studied with the same level of detail. Methods. A flare occurring in March 2010 was observed for 13 consecutive days (from MJD 55 265 to MJD 55 277) with unprecedented wavelength coverage from radio to very high energy (VHE; E> 100 GeV) γ-rays with MAGIC, VERITAS, Whipple, Fermi-LAT, MAXI, RXTE, Swift, GASP-WEBT, and several optical and radio telescopes. We modeled the day-scale SEDs with one-zone and two-zone synchrotron self-Compton (SSC) models, investigated the physical parameters, and evaluated whether the observed broadband SED variability can be associated with variations in the relativistic particle population. Results. The activity of Mrk 421 initially was high and then slowly decreased during the 13-day period. The flux variability was remarkable at the X-ray and VHE bands, but it was minor or not significant at the other bands. The variability in optical polarization was also minor. These observations revealed an almost linear correlation between the X-ray flux at the 2–10 keV band and the VHE γ-ray flux above 200 GeV, consistent with the γ-rays being produced by inverse-Compton scattering in the Klein-Nishina regime in the framework of SSC models. The one-zone SSC model can describe the SED of each day for the 13 consecutive days reasonably well, which once more shows the success of this standard theoretical scenario to describe the SEDs of VHE BL Lacs such as Mrk 421. This flaring activity is also very well described by a two-zone SSC model, where one zone is responsible for the quiescent emission, while the other smaller zone, which is spatially separated from the first, contributes to the daily variable emission occurring at X-rays and VHE γ-rays. The second blob is assumed to have a smaller volume and a narrow electron energy distribution with 3 × 104<γ< 6 × 105, where γ is the Lorentz factor of the electrons. Such a two-zone scenario would naturally lead to the correlated variability at the X-ray and VHE bands without variability at the optical/UV band, as well as to shorter timescales for the variability at the X-ray and VHE bands with respect to the variability at the other bands. Conclusions. Both the one-zone and the two-zone SSC models can describe the daily SEDs via the variation of only four or five model parameters, under the hypothesis that the variability is associated mostly with the underlying particle population. This shows that the particle acceleration and cooling mechanism that produces the radiating particles might be the main mechanism responsible for the broadband SED variations during the flaring episodes in blazars. The two-zone SSC model provides a better agreement with the observed SED at the narrow peaks of the low- and high-energy bumps during the highest activity, although the reported one-zone SSC model could be further improved by varying the parameters related to the emitting region itself (δ, B and R), in addition to the parameters related to the particle population.Fil: Aleksic, J.. IFAE; EspañaFil: Ansoldi, S.. Università di Udine; ItaliaFil: Antonelli, L. A.. INAF National Institute for Astrophysics; ItaliaFil: Antoranz, P.. Università di Siena; ItaliaFil: Babic, A.. University of Rijeka; CroaciaFil: Bangale, P.. Max-Planck-Institut für Physik; AlemaniaFil: Barres de Almeida, U.. Max-Planck-Institut für Physik; AlemaniaFil: Barrio, J. A.. Universidad Complutense de Madrid; EspañaFil: Becerra Gonzalez, J.. Inst. de Astrofísica de Canarias; EspañaFil: Bednarek, W.. University of Łódź,; PoloniaFil: Bernardini, E.. Deutsches Elektronen-Synchrotron (DESY); AlemaniaFil: Biasuzzi, B.. Università di Udine; ItaliaFil: Biland, A.. ETH Zurich; SuizaFil: Blanch, O.. IFAE; EspañaFil: Boller, A.. ETH Zurich; SuizaFil: Bonnefoy, S.. Universidad Complutense de Madrid; EspañaFil: Bonnoli, G.. INAF National Institute for Astrophysics ; ItaliaFil: Borracci, F.. Max-Planck-Institut für Physik; AlemaniaFil: Bretz, T.. Universität Würzburg ; AlemaniaFil: Carmona, E.. Centro de Investigaciones Energéticas, Medioambientales y Tecnológicas ; EspañaFil: Carosi, A.. INAF National Institute for Astrophysics; ItaliaFil: Colin, P.. Max-Planck-Institut für Physik; AlemaniaFil: Colombo, E.. Inst. de Astrofísica de Canarias; EspañaFil: Contreras, J. L.. Universidad Complutense; EspañaFil: Cortina, J.. IFAE; EspañaFil: Covino, S.. INAF National Institute for Astrophysics; ItaliaFil: Da Vela, P.. Università di Siena; ItaliaFil: Dazzi, F.. Max-Planck-Institut für Physik; AlemaniaFil: Pichel, Ana Carolina. Consejo Nacional de Investigaciónes Científicas y Técnicas. Oficina de Coordinación Administrativa Ciudad Universitaria. Instituto de Astronomía y Física del Espacio. - Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Instituto de Astronomía y Física del Espacio; ArgentinaFil: Rovero, Adrian Carlos. Consejo Nacional de Investigaciónes Científicas y Técnicas. Oficina de Coordinación Administrativa Ciudad Universitaria. Instituto de Astronomía y Física del Espacio. - Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Instituto de Astronomía y Física del Espacio; ArgentinaFil: The Veritas Collaboration.Fil: The MAGIC Collaboration

Unprecedented study of the broadband emission of Mrk 421 during flaring activity in March 2010

Context. Because of its proximity, Mrk 421 is one of the best sources on which to study the nature of BL Lac objects. Its proximity allows us to characterize its broadband spectral energy distribution (SED). Aims. The goal is to better understand the mechanisms responsible for the broadband emission and the temporal evolution of Mrk 421. These mechanisms may also apply to more distant blazars that cannot be studied with the same level of detail. Methods. A flare occurring in March 2010 was observed for 13 consecutive days (from MJD 55 265 to MJD 55 277) with unprecedented wavelength coverage from radio to very high energy (VHE; E > 100 GeV) γ-rays with MAGIC, VERITAS, Whipple, Fermi-LAT, MAXI, RXTE, Swift, GASP-WEBT, and several optical and radio telescopes. We modeled the day-scale SEDs with one-zone and two-zone synchrotron self-Compton (SSC) models, investigated the physical parameters, and evaluated whether the observed broadband SED variability can be associated with variations in the relativistic particle population. Results. The activity of Mrk 421 initially was high and then slowly decreased during the 13-day period. The flux variability was remarkable at the X-ray and VHE bands, but it was minor or not significant at the other bands. The variability in optical polarization was also minor. These observations revealed an almost linear correlation between the X-ray flux at the 2–10 keV band and the VHE γ-ray flux above 200 GeV, consistent with the γ-rays being produced by inverse-Compton scattering in the Klein-Nishina regime in the framework of SSC models. The one-zone SSC model can describe the SED of each day for the 13 consecutive days reasonably well, which once more shows the success of this standard theoretical scenario to describe the SEDs of VHE BL Lacs such as Mrk 421. This flaring activity is also very well described by a two-zone SSC model, where one zone is responsible for the quiescent emission, while the other smaller zone, which is spatially separated from the first, contributes to the daily variable emission occurring at X-rays and VHE γ-rays. The second blob is assumed to have a smaller volume and a narrow electron energy distribution with 3 × 10⁴ < γ < 6 × 10⁵, where γ is the Lorentz factor of the electrons. Such a two-zone scenario would naturally lead to the correlated variability at the X-ray and VHE bands without variability at the optical/UV band, as well as to shorter timescales for the variability at the X-ray and VHE bands with respect to the variability at the other bands. Conclusions. Both the one-zone and the two-zone SSC models can describe the daily SEDs via the variation of only four or five model parameters, under the hypothesis that the variability is associated mostly with the underlying particle population. This shows that the particle acceleration and cooling mechanism that produces the radiating particles might be the main mechanism responsible for the broadband SED variations during the flaring episodes in blazars. The two-zone SSC model provides a better agreement with the observed SED at the narrow peaks of the low- and high-energy bumps during the highest activity, although the reported one-zone SSC model could be further improved by varying the parameters related to the emitting region itself (δ, B and R), in addition to the parameters related to the particle population.J. Aleksić ... S. Einecke ... V. Stamatescu ... et al

Multiband variability studies and novel broadband SED modeling of Mrk 501 in 2009

Context. We present an extensive study of the BL Lac object Mrk 501 based on a data set collected during the multi-instrument campaign spanning from 2009 March 15 to 2009 August 1, which includes, among other instruments, MAGIC, VERITAS, Whipple 10 m, and Fermi-LAT to cover the γ-ray range from 0.1 GeV to 20 TeV; RXTE and Swift to cover wavelengths from UV tohard X-rays; and GASP-WEBT, which provides coverage of radio and optical wavelengths. Optical polarization measurements were provided for a fraction of the campaign by the Steward and St. Petersburg observatories. We evaluate the variability of the source and interband correlations, the γ-ray flaring activity occurring in May 2009, and interpret the results within two synchrotron self-Compton (SSC) scenarios. Aims. The multiband variability observed during the full campaign is addressed in terms of the fractional variability, and the possible correlations are studied by calculating the discrete correlation function for each pair of energy bands where the significance was evaluated with dedicated Monte Carlo simulations. The space of SSC model parameters is probed following a dedicated grid-scan strategy, allowing for a wide range of models to be tested and offering a study of the degeneracy of model-to-data agreement in the individual model parameters, hence providing a less biased interpretation than the “single-curve SSC model adjustment” typically reported in the literature. Methods. We find an increase in the fractional variability with energy, while no significant interband correlations of flux changes are found on the basis of the acquired data set. The SSC model grid-scan shows that the flaring activity around May 22 cannot be modeled adequately with a one-zone SSC scenario (using an electron energy distribution with two breaks), while it can be suitably described within a two (independent) zone SSC scenario. Here, one zone is responsible for the quiescent emission from the averaged 4.5-month observing period, while the other one, which is spatially separated from the first, dominates the flaring emission occurring at X-rays and very-high-energy (>100 GeV, VHE) γ rays. The flaring activity from May 1, which coincides with a rotation of the electric vector polarization angle (EVPA), cannot be satisfactorily reproduced by either a one-zone or a two-independent-zone SSC model, yet this is partially affected by the lack of strictly simultaneous observations and the presence of large flux changes on sub-hour timescales (detected at VHE γ rays). Results. The higher variability in the VHE emission and lack of correlation with the X-ray emission indicate that, at least during the 4.5-month observing campaign in 2009, the highest energy (and most variable) electrons that are responsible for the VHE γ rays do not make a dominant contribution to the ~1 keV emission. Alternatively, there could be a very variable component contributing to the VHE γ-ray emission in addition to that coming from the SSC scenario. The studies with our dedicated SSC grid-scan show that there is some degeneracy in both the one-zone and the two-zone SSC scenarios probed, with several combinations of model parameters yielding a similar model-to-data agreement, and some parameters better constrained than others. The observed γ-ray flaring activity, with the EVPA rotation coincident with the first γ-ray flare, resembles those reported previously for low frequency peaked blazars, hence suggesting that there are many similarities in the flaring mechanisms of blazars with different jet properties.M. L. Ahnen ... S. Einecke ... et al

The 2009 multiwavelength campaign on Mrk 421: variability and correlation studies

Aims: We perform an extensive characterization of the broadband emission of Mrk 421, as well as its temporal evolution, during the non-flaring (low) state. The high brightness and nearby location (z = 0.031) of Mrk 421 make it an excellent laboratory to study blazar emission. The goal is to learn about the physical processes responsible for the typical emission of Mrk 421, which might also be extended to other blazars that are located farther away and hence are more difficult to study. Methods: We performed a 4.5-month multi-instrument campaign on Mrk 421 between January 2009 and June 2009, which included VLBA, F-GAMMA, GASP-WEBT, Swift, RXTE, Fermi-LAT, MAGIC, and Whipple, among other instruments and collaborations. This extensive radio to very-high-energy (VHE; E > 100 GeV) γ-ray dataset provides excellent temporal and energy coverage, which allows detailed studies of the evolution of the broadband spectral energy distribution. Results: Mrk421 was found in its typical (non-flaring) activity state, with a VHE flux of about half that of the Crab Nebula, yet the light curves show significant variability at all wavelengths, the highest variability being in the X-rays. We determined the power spectral densities (PSD) at most wavelengths and found that all PSDs can be described by power-laws without a break, and with indices consistent with pink/red-noise behavior. We observed a harder-when-brighter behavior in the X-ray spectra and measured a positive correlation between VHE and X-ray fluxes with zero time lag. Such characteristics have been reported many times during flaring activity, but here they are reported for the first time in the non-flaring state. We also observed an overall anti-correlation between optical/UV and X-rays extending over the duration of the campaign. Conclusions: The harder-when-brighter behavior in the X-ray spectra and the measured positive X-ray/VHE correlation during the 2009 multiwavelength campaign suggests that the physical processes dominating the emission during non-flaring states have similarities with those occurring during flaring activity. In particular, this observation supports leptonic scenarios as being responsible for the emission of Mrk 421 during nonflaring activity. Such a temporally extended X-ray/VHE correlation is not driven by any single flaring event, and hence is difficult to explain within the standard hadronic scenarios. The highest variability is observed in the X-ray band, which, within the one-zone synchrotron self-Compton scenario, indicates that the electron energy distribution is most variable at the highest energies.J. Aleksić ... S. Einecke... et al