A fast, very-high-energy γ -ray flare from BL Lacertae during a period of multi-wavelength activity in June 2015

The mechanisms producing fast variability of the γ-ray emission in active galactic nuclei (AGNs) are under debate. The MAGIC telescopes detected a fast, very-high-energy (VHE, E > 100 GeV) γ-ray flare from BL Lacertae on 2015 June 15. The flare had a maximum flux of (1.5 ± 0.3) × 10-10 photons cm-2 s-1 and halving time of 26 ± 8 min. The MAGIC observations were triggered by a high state in the optical and high-energy (HE, E > 100 MeV) γ-ray bands. In this paper we present the MAGIC VHE γ-ray data together with multi-wavelength data from radio, optical, X-rays, and HE γ rays from 2015 May 1 to July 31. Well-sampled multi-wavelength data allow us to study the variability in detail and compare it to the other epochs when fast, VHE γ-ray flares have been detected from this source. Interestingly, we find that the behaviour in radio, optical, X-rays, and HE γ-rays is very similar to two other observed VHE γ-ray flares. In particular, also during this flare there was an indication of rotation of the optical polarization angle and of activity at the 43 GHz core. These repeating patterns indicate a connection between the three events. We also test modelling of the spectral energy distribution based on constraints from the light curves and VLBA observations, with two different geometrical setups of two-zone inverse Compton models. In addition we model the γ-ray data with the star-jet interaction model. We find that all of the tested emission models are compatible with the fast VHE γ-ray flare, but all have some tension with the multi-wavelength observations

Constraining Lorentz Invariance Violation Using the Crab Pulsar Emission Observed up to TeV Energies by MAGIC

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 E-QG1 > 5.5 x 10(17) GeV (4.5 x 10(17) GeV) for a linear, and E-QG2 > 5.9 x 10(10) GeV (5.3 x 10(10) GeV) 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

MAGIC observations of the microquasar V404 Cygni during the 2015 outburst

The microquasar V404 Cygni underwent a series of outbursts in 2015, June 15-31, during which its flux in hard X-rays (20-40 keV) reached about 40 times the Crab nebula flux. Because of the exceptional interest of the flaring activity from this source, observations at several wavelengths were conducted. The MAGIC telescopes, triggered by the INTEGRAL alerts, followed-up the flaring source for several nights during the period June 18-27, for more than 10 h. One hour of observation was simultaneously conducted on a giant 22 GHz radio flare and a hint of signal at GeV energies seen by Fermi-LAT. The MAGIC observations did not show significant emission in any of the analysed time intervals. The derived flux upper limit, in the energy range 200-1250 GeV, is 4.8 x 10(-12) photons cm(-2) s(-1). We estimate the gamma-ray opacity during the flaring period, which along with our non-detection points to an inefficient acceleration in the V404 Cyg jets if a very high energy emitter is located further than 1 x 10(10) cm from the compact object

Search for very high-energy gamma-ray emission from the microquasar Cygnus X-1 with the MAGIC telescopes

The microquasar Cygnus X-1 displays the two typical soft and hard X-ray states of a black hole transient. During the latter, Cygnus X-1 shows a one-sided relativistic radio-jet. Recent detection of the system in the high energy (HE; E greater than or similar to 60 MeV) gamma-ray range with FermiLAT associates this emission with the outflow. Former MAGIC observations revealed a hint of flaring activity in the very high-energy (VHE; E greater than or similar to 100 GeV) regime during this X-ray state. We analyse similar to 97 h of Cygnus X-1 data taken with the MAGIC telescopes between July 2007 and October 2014. To shed light on the correlation between hard X-ray and VHE gamma rays as previously suggested, we study each main X-ray state separately. We perform an orbital phase-folded analysis to look for variability in the VHE band. Additionally, to place this variability behaviour in a multiwavelength context, we compare our results with Fermi-LAT, AGILE, Swift-BAT, MAXI, RXTE-ASM, AMI and RATAN-600 data. We do not detect Cygnus X-1 in the VHE regime. We establish upper limits for each X-ray state, assuming a power-law distribution with photon index Gamma = 3.2. For steady emission in the hard and soft X-ray states, we set integral upper limits at 95 per cent confidence level for energies above 200 GeV at 2.6 x 10(-12) photons cm(-2) s(-1) and 1.0 x 10(-11) photons cm(-2) s(-1), respectively. We rule out steady VHE gamma-ray emission above this energy range, at the level of theMAGIC sensitivity, originating in the interaction between the relativistic jet and the surrounding medium, while the emission above this flux level produced inside the binary still remains a valid possibility

Observation of the black widow B1957+20 millisecond pulsar binary system with the MAGIC telescopes

B1957+20 is a millisecond pulsar located in a black-widow-type compact binary system with a low-mass stellar companion. The interaction of the pulsar wind with the companion star wind and/or the interstellar plasma is expected to create plausible conditions for acceleration of electrons to TeV energies and subsequent production of very high-energy gamma-rays in the inverse Compton process. We performed extensive observations with the Major Atmospheric Gamma Imaging Cherenkov Telescopes (MAGIC) telescopes of B1957+20. We interpret results in the framework of a few different models, namely emission from the vicinity of the millisecond pulsar, the interaction of the pulsar and stellar companion wind region or bow shock nebula. No significant steady very high-energy gamma-ray emission was found. We derived a 95 per cent confidence level upper limit of 3.0 x 10(-12) cm(-2) s(-1) on the average gamma-ray emission from the binary system above 200 GeV. The upper limits obtained with the MAGIC constrain, for the first time, different models of the high-energy emission in B1957+20. In particular, in the inner mixed wind nebula model with mono-energetic injection of electrons, the acceleration efficiency of electrons is constrained to be below similar to 2-10 per cent of the pulsar spin-down power. For the pulsar emission, the obtained upper limits for each emission peak are well above the exponential cut-off fits to the Fermi-LAT data, extrapolated to energies above 50 GeV. The MAGIC upper limits can rule out a simple power-law tail extension through the sub-TeV energy range for the main peak seen at radio frequencies

Multiwavelength observations of a VHE gamma-ray flare from PKS 1510-089 in 2015

Context. PKS 1510 089 is one of only a few flat spectrum radio quasars detected in the very-high-energy (VHE, > 100 GeV) gamma-ray band.Aims. We study the broadband spectral and temporal properties of the PKS 1510 089 emission during a high gamma-ray state.Methods. We performed VHE gamma-ray observations of PKS 1510 089 with the Major Atmospheric Gamma Imaging Cherenkov (MAGIC) telescopes during a long, high gamma-ray state in May 2015. In order to perform broadband modeling of the source, we have also gathered contemporaneous multiwavelength data in radio, IR, optical photometry and polarization, UV, X-ray, and GeV gamma-ray ranges. We construct a broadband spectral energy distribution (SED) in two periods, selected according to VHE gamma-ray state.Results. PKS 1510 089 was detected by MAGIC during a few day-long observations performed in the middle of a long, high optical and gamma-ray state, showing for the first time a significant VHE gamma-ray variability. Similarly to the optical and gamma-ray high state of the source detected in 2012, it was accompanied by a rotation of the optical polarization angle and the emission of a new jet component observed in radio. However, owing to large uncertainty on the knot separation time, the association with the VHE gamma-ray emission cannot be firmly established. The spectral shape in the VHE band during the flare is similar to those obtained during previous measurements of the source. The observed flux variability sets constraints for the first time on the size of the region from which VHE gamma rays are emitted. We model the broadband SED in the framework of the external Compton scenario and discuss the possible emission site in view of multiwavelength data and alternative emission models

Multi-wavelength characterization of the blazar S5 0716+714 during an unprecedented outburst phase

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, Metsahovi, 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 y-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.02(stats) +/- 0.05(sys), 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

Detection of persistent VHE gamma-ray emission from PKS 1510-089 by the MAGIC telescopes during low states between 2012 and 2017

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 gammaray 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 emissio-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 x10(-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.50-sigma ) VHE gamma-ray emission at the level of (4.27 +/- 0.61(stat)) x 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

Periastron Observations of TeV Gamma-Ray Emission from a Binary System with a 50-year Period

We report on observations of the pulsar/Be star binary system PSR J2032+4127/MT91 213 in the energy range between 100 GeV and 20 TeV with the Very Energetic Radiation Imaging Telescope Array and Major Atmospheric Gamma Imaging Cherenkov telescope arrays. The binary orbit has a period of approximately 50 years, with the most recent periastron occurring on 2017 November 13. Our observations span from 18 months prior to periastron to one month after. A new point-like gamma-ray source is detected, coincident with the location of PSR J2032+4127/MT91 213. The gamma-ray light curve and spectrum are well characterized over the periastron passage. The flux is variable over at least an order of magnitude, peaking at periastron, thus providing a firm association of the TeV source with the pulsar/Be star system. Observations prior to periastron show a cutoff in the spectrum at an energy around 0.5 TeV. This result adds a new member to the small population of known TeV binaries, and it identifies only the second source of this class in which the nature and properties of the compact object are firmly established. We compare the gamma-ray results with the light curve measured with the X-ray Telescope on board the Neil Gehrels Swift Observatory and with the predictions of recent theoretical models of the system. We conclude that significant revision of the models is required to explain the details of the emission that we have observed, and we discuss the relationship between the binary system and the overlapping steady extended source, TeV J2032+4130

MAGIC and Fermi-LAT gamma-ray results on unassociated HAWC sources

The HAWC Collaboration released the 2HWC catalogue of TeV sources, in which 19 show no association with any known high-energy (HE; E greater than or similar to 10 GeV) or very-high-energy (VHE; E greater than or similar to 300 GeV) sources. This catalogue motivated follow-up studies by both the Major Atmospheric Gamma-ray Imaging Cherenkov (MAGIC) and Fermi-LAT (Large Area Telescope) observatories with the aim of investigating gamma-ray emission over a broad energy band. In this paper, we report the results from the first joint work between High Altitude Water Cherenkov (HAWC), MAGIC, and Fermi-LAT on three unassociated HAWC sources: 2HWC J2006+341, 2HWC J1907+084*, and 2HWC J1852+013*. Although no significant detection was found in the HE and VHE regimes, this investigation shows that a minimum 1 degrees extension (at 95 per cent confidence level) and harder spectrum in the GeV than the one extrapolated from HAWC results are required in the case of 2HWC J1852+013*, whilst a simply minimum extension of 0.16 degrees (at 95 per cent confidence level) can already explain the scenario proposed by HAWC for the remaining sources. Moreover, the hypothesis that these sources are pulsar wind nebulae is also investigated in detail