Broadband characterisation of the very intense TeV flares of the blazar 1ES 1959+650 in 2016

1ES 1959+650 is a bright TeV high-frequency-peaked BL Lac object exhibiting interesting features like "orphan" TeV flares and broad emission in the high-energy regime that are difficult to interpret using conventional one-zone Synchrotron Self-Compton (SSC) scenarios. We report the results from the Major Atmospheric Gamma Imaging Cherenkov (MAGIC) observations in 2016 along with the multi-wavelength data from the Fermi Large Area Telescope (LAT) and Swift instruments. MAGIC observed 1ES 1959+650 with different emission levels in the very-high-energy (VHE, E> 100 GeV) gamma -ray band during 2016. In the long-term data, the X-ray spectrum becomes harder with increasing flux and a hint of a similar trend is also visible in the VHE band. An exceptionally high VHE flux reaching similar to 3 times the Crab Nebula flux was measured by MAGIC on the 13 and 14 of June, and 1 July 2016 (the highest flux observed since 2002). During these flares, the high-energy peak of the spectral energy distribution (SED) lies in the VHE domain and extends up to several TeV. The spectrum in the gamma -ray (both Fermi-LAT and VHE bands) and the X-ray bands are quite hard. On 13 June and 1 July 2016, the source showed rapid variations in the VHE flux within timescales of less than an hour. A simple one-zone SSC model can describe the data during the flares requiring moderate to large values of the Doppler factors (delta >= 30-60). Alternatively, the high-energy peak of the SED can be explained by a purely hadronic model attributed to proton-synchrotron radiation with jet power L-jet similar to 10(46) erg s(-1) and under high values of the magnetic field strength (similar to 100 G) and maximum proton energy (similar to few EeV). Mixed lepto-hadronic models require super-Eddington values of the jet power. We conclude that it is difficult to get detectable neutrino emission from the source during the extreme VHE flaring period of 2016

Measurement of the extragalactic background light using MAGIC and Fermi-LAT gamma-ray observations of blazars up to z=1

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 gamma gamma 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 sigma 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 m range relative to the studied models, yet compatible with them within systematics

Constraints on Gamma-Ray and Neutrino Emission from NGC 1068 with the MAGIC Telescopes

Starburst galaxies and star-forming active galactic nuclei are among the candidate sources thought to contribute appreciably to the extragalactic gamma-ray and neutrino backgrounds. NGC 1068 is the brightest of the star-forming galaxies found to emit gamma-rays from 0.1 to 50 GeV. Precise measurements of the high-energy spectrum are crucial to study the particle accelerators and probe the dominant emission mechanisms. We have carried out 125 hr of observations of NGC 1068 with the MAGIC telescopes in order to search for gamma-ray emission in the very-high-energy band. We did not detect significant gamma-ray emission, and set upper limits at the 95% confidence level to the gamma-ray flux above 200 GeV f < 5.1. x. 10(-13) cm(-2) s(-1). This limit improves previous constraints by about an order of magnitude and allows us to put tight constraints on the theoretical models for the gamma-ray emission. By combining the MAGIC observations with the Fermi-LAT spectrum we limit the parameter space (spectral slope, maximum energy) of the cosmic ray protons predicted by hadronuclear models for the gamma-ray emission, while we find that a model postulating leptonic emission from a semi-relativistic jet is fully consistent with the limits. We provide predictions for IceCube detection of the neutrino signal foreseen in the hadronic scenario. We predict a maximal IceCube neutrino event rate of 0.07 yr(-1)

New Hard-TeV Extreme Blazars Detected with the MAGIC Telescopes*

Extreme high-frequency-peaked BL Lac objects (EHBLs) are blazars that exhibit extremely energetic synchrotron emission. They also feature nonthermal gamma-ray emission whose peak lies in the very high-energy (VHE, E > 100 GeV) range, and in some sources exceeds 1 TeV: this is the case for hard-TeV EHBLs such as 1ES 0229+200. With the aim of increasing the EHBL population, 10 targets were observed with the MAGIC telescopes from 2010 to 2017, for a total of 265 hr of good-quality data. The data were complemented by coordinated Swift observations. The X-ray data analysis confirms that all but two sources are EHBLs. The sources show only a modest variability and a harder-when-brighter behavior, typical for this class of objects. At VHE gamma-rays, three new sources were detected and a hint of a signal was found for another new source. In each case, the intrinsic spectrum is compatible with the hypothesis of a hard-TeV nature of these EHBLs. The broadband spectral energy distributions (SEDs) of all sources are built and modeled in the framework of a single-zone, purely leptonic model. The VHE gamma-ray-detected sources were also interpreted with a spine-layer model and a proton synchrotron model. The three models provide a good description of the SEDs. However, the resulting parameters differ substantially in the three scenarios, in particular the magnetization parameter. This work presents the first mini catalog of VHE gamma-ray and multiwavelength observations of EHBLs

Testing emission models on the extreme blazar 2WHSP J073326.7+515354 detected at very high energies with the MAGIC telescopes

Extreme high-energy-peaked BL Lac objects (EHBLs) are an emerging class of blazars. Their typical two-hump-structured spectral energy distribution (SED) peaks at higher energies with respect to conventional blazars. Multiwavelength (MWL) observations constrain their synchrotron peak in the medium to hard X-ray band. Their gamma-ray SED peaks above the GeV band, and in some objects it extends up to several TeV. Up to now, only a few EHBLs have been detected in the TeV gamma-ray range. In this paper, we report the detection of the EHBL 2WHSP J073320,7+515354, observed and detected during 2018 in TeV gamma rays with the MAGIC telescopes. The broad-band SED is studied within an MWL context, including an analysis of the Fermi-LAT data over 10 yr of observation and with simultaneous Swift-XRT, Swift-UVOT, and KVA data. Our analysis results in a set of spectral parameters that confirms the classification of the source as an EIME. In order to investigate the physical nature of this extreme emission, different theoretical frameworks were tested to model the broadband SED. The hard TeV spectrum of 2WHSP J073326.7+515354 sets the SED far from the energy equipartition regime in the standard one-zone leptonic scenario of blazar emission. Conversely, more complex models of the jet, represented by either a two-zone spine-layer model or a hadronic emission model, better represent the broad-hand SED

Monitoring of the radio galaxy M 87 during a low -emission state from m 2012 to 2015 with MAGIC

M 87 is one of the closest (z = 0.004 36) extragalactic sources emitting at very high energies (VHF, E > 100 GeV). The aim of this work is to locale the region of the VHF gamma-ray emission and to describe the observed broad-band spectral energy distribution (SED) during the low VHF gamma-ray state. The data from M 87 collected between 2012 and 2015 as part of a MAGIC monitoring programme are analysed and combined with multiwavelength data from Fermi-LAT, Chandra, HST, FVN, VLBA, and the Liverpool Telescope. The averaged VHE gamma-ray spectrum can be fitted from 100 GeV to 10 TeV with a simple power law with a photon index of (-2.41 0.07), while the integral flux above 300 GeV is (1.44 0.13) x 10-12 cm 2 s I. During the campaign between 2012 and 2015, M87 is generally found in a low-emission state at all observed wavelengths. The VIIE gamma-ray flux from the present 2012-2015M 87 campaign is consistent with a constant flux with some hint of variability ( 3 a) on a daily time-scale in 2013. The low-state gamma-ray emission likely originates from the same region as the flare-state emission. Given the broad-band SED, both a leptonic synchrotron self-Compton and a hybrid photohadronic model reproduce the available data well, even if the latter is preferred. We note, however, that the energy stored in the magnetic field in the leptonic scenario is very low, suggesting a matter-dominated emission region

Search for Very High-energy Emission from the Millisecond Pulsar PSR J0218+4232

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) gamma-ray emission. Using 11.5 yr of Fermi Large Area Telescope (LAT) data between 100 MeV and 870 GeV, and similar to 90 hr of Major Atmospheric Gamma Imaging Cherenkov (MAGIC) observations in the 20 GeV to 20 TeV range, we searched for the highest energy gamma-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 gamma-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

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

Unraveling the Complex Behavior of Mrk 421 with Simultaneous X-Ray and VHE Observations during an Extreme Flaring Activity in 2013 April*

We report on a multiband variability and correlation study of the TeV blazar Mrk 421 during an exceptional flaring activity observed from 2013 April 11 to 19. The study uses, among others, data from GLAST-AGILE Support Program (GASP) of the Whole Earth Blazar Telescope (WEBT), Swift, Nuclear Spectroscopic Telescope Array (NuSTAR), Fermi Large Area Telescope, Very Energetic Radiation Imaging Telescope Array System (VERITAS), and Major Atmospheric Gamma Imaging Cherenkov (MAGIC). The large blazar activity and the 43 hr of simultaneous NuSTAR and MAGIC/VERITAS observations permitted variability studies on 15 minute time bins over three X-ray bands (3-7 keV, 7-30 keV, and 30-80 keV) and three very-high-energy (VHE; >0.1 TeV) 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 of 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-rays 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 characteristics of the correlation change substantially over time and across energy bands. Our findings favor 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

Investigation of the correlation patterns and the Compton dominance variability of Mrk 421 in 2017

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 Metsahovi. 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 sigma. 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 gamma(min)' = 2 x 10(4) to gamma(max)' = 6 x 10(5).</p