Previous detections of individual astrophysical sources of neutrinos are limited to the Sun and the supernova 1987A, whereas the origins of the diffuse flux of high-energy cosmic neutrinos remain unidentified. On 22 September 2017, we detected a high-energy neutrino, IceCube-170922A, with an energy of ~290 tera-electron volts. Its arrival direction was consistent with the location of a known γ-ray blazar, TXS 0506+056, observed to be in a flaring state. An extensive multiwavelength campaign followed, ranging from radio frequencies to γ-rays. These observations characterize the variability and energetics of the blazar and include the detection of TXS 0506+056 in very-high-energy γ-rays. This observation of a neutrino in spatial coincidence with a γ-ray-emitting blazar during an active phase suggests that blazars may be a source of high-energy neutrinos.

A neutrino with energy ∼290 TeV, IceCube-170922A, was detected in coincidence with the BL Lac object TXS0506+056 during enhanced gamma-ray activity, with chance coincidence being rejected at ∼3σ level. We monitored the object in the very-high-energy (VHE) band with the Major Atmospheric Gamma-ray Imaging Cherenkov (MAGIC) telescopes for ∼41 hr from 1.3 to 40.4 days after the neutrino detection. Day-timescale variability is clearly resolved. We interpret the quasi-simultaneous neutrino and broadband electromagnetic observations with a novel one-zone lepto-hadronic model, based on interactions of electrons and protons coaccelerated in the jet with external photons originating from a slow-moving plasma sheath surrounding the faster jet spine. We can reproduce the multiwavelength spectra of TXS 0506+056 with neutrino rate and energy compatible with IceCube-170922A, and with plausible values for the jet power of 10 4 10 erg s 45 46 1 -´-. The steep spectrum observed by MAGIC is concordant with internal γγ absorption above ∼100 GeV entailed by photohadronic production of a ∼290 TeV neutrino, corroborating a genuine connection between the multi-messenger signals. In contrast to previous predictions of predominantly hadronic emission from neutrino sources, the gamma-rays can be mostly ascribed to inverse Compton upscattering of external photons by accelerated electrons. The X-ray and VHE bands provide crucial constraints on the emission from both accelerated electrons and protons. We infer that the maximum energy of protons in the jet comoving frame can be in the range ∼10 14 -10 18 eV.

MAGIC, a system of two imaging atmospheric Cherenkov telescopes, achieves its best performance under dark conditions, i.e. in absence of moonlight or twilight. Since operating the telescopes only during dark time would severely limit the duty cycle, observations are also performed when the Moon is present in the sky. Here we develop a dedicated Moon-adapted analysis to characterize the performance of MAGIC under moonlight. We evaluate energy threshold, angular resolution and sensitivity of MAGIC under different background light levels, based on Crab Nebula observations and tuned Monte Carlo simulations. This study includes observations taken under non-standard hardware configurations, such as reducing the camera photomultiplier tubes gain by a factor â\u88¼ 1.7 (reduced HV settings) with respect to standard settings (nominal HV) or using UV-pass filters to strongly reduce the amount of moonlight reaching the cameras of the telescopes. The Crab Nebula spectrum is correctly reconstructed in all the studied illumination levels, that reach up to 30 times brighter than under dark conditions. The main effect of moonlight is an increase in the analysis energy threshold and in the systematic uncertainties on the flux normalization. The sensitivity degradation is constrained to be below 10%, within 15â\u80\u9330% and between 60 and 80% for nominal HV, reduced HV and UV-pass filter observations, respectively. No worsening of the angular resolution was found. Thanks to observations during moonlight, the maximal duty cycle of MAGIC can be increased from â\u88¼ 18%, under dark nights only, to up to â\u88¼ 40% in total with only moderate performance degradation

Aims. We aim to characterize the multiwavelength emission from Markarian 501 (Mrk 501), quantify the energy-dependent variability, study the potential multiband correlations, and describe the temporal evolution of the broadband emission within leptonic theoretical scenarios. Methods. We organized a multiwavelength campaign to take place between March and July of 2012. Excellent temporal coverage was obtained with more than 25 instruments, including the MAGIC, FACT and VERITAS Cherenkov telescopes, the instruments on board the Swift and Fermi spacecraft, and the telescopes operated by the GASP-WEBT collaboration. Results. Mrk 501 showed a very high energy (VHE) gamma-ray flux above 0.2 TeV of ∼0.5 times the Crab Nebula flux (CU) for most of the campaign. The highest activity occurred on 2012 June 9, when the VHE flux was ∼3 CU, and the peak of the high-energy spectral component was found to be at ∼2 TeV. Both the X-ray and VHE gamma-ray spectral slopes were measured to be extremely hard, with spectral indices < 2 during most of the observing campaign, regardless of the X-ray and VHE flux. This study reports the hardest Mrk 501 VHE spectra measured to date. The fractional variability was found to increase with energy, with the highest variability occurring at VHE. Using the complete data set, we found correlation between the X-ray and VHE bands; however, if the June 9 flare is excluded, the correlation disappears (significance < 3σ) despite the existence of substantial variability in the X-ray and VHE bands throughout the campaign. Conclusions. The unprecedentedly hard X-ray and VHE spectra measured imply that their low- and high-energy components peaked above 5 keV and 0.5 TeV, respectively, during a large fraction of the observing campaign, and hence that Mrk 501 behaved like an extreme high-frequency-peaked blazar (EHBL) throughout the 2012 observing season. This suggests that being an EHBL may not be a permanent characteristic of a blazar, but rather a state which may change over time. The data set acquired shows that the broadband spectral energy distribution (SED) of Mrk 501, and its transient evolution, is very complex, requiring, within the framework of synchrotron self-Compton (SSC) models, various emission regions for a satisfactory description. Nevertheless the one-zone SSC scenario can successfully describe the segments of the SED where most energy is emitted, with a significant correlation between the electron energy density and the VHE gamma-ray activity, suggesting that most of the variability may be explained by the injection of high-energy electrons. The one-zone SSC scenario used reproduces the behavior seen between the measured X-ray and VHE gamma-ray fluxes, and predicts that the correlation becomes stronger with increasing energy of the X-rays.

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 which 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 γγ 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σ upper bounds, including systematic uncertainties, are between 13% and 23% 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 µm range relative to the studied models, yet compatible with them within systematics.

We analyze Chandra X-ray observatory data for a sample of 63 luminous infrared galaxies (LIRGs), sampling the lower-infrared luminosity range of the Great Observatories All-Sky LIRG survey (GOALS), which includes the most luminous infrared selected galaxies in the local Universe. X-rays are detected for 84 individual galaxies within the 63 systems, for which arcsecond resolution X-ray images, fluxes, infrared and X-ray luminosities, spectra and radial profiles are presented. Using X-ray and mid-infrared (MIR) selection criteria, we find AGN in (31 ± 5)% of the galaxy sample, compared to the (38 ± 6)% previously found for GOALS galaxies with higher infrared luminosities (C-GOALS I). Using MIR data, we find that (59 ± 9)% of the X-ray selected AGN in the full C-GOALS sample do not contribute significantly to the bolometric luminosity of the host galaxy. Dual AGN are detected in two systems, implying a dual AGN fraction in systems that contain at least one AGN of (29 ± 14)%, compared to the (11 ± 10)% found for the C-GOALS I sample. Through analysis of radial profiles, we derive that most sources, and almost all AGN, in the sample are compact, with half of the soft X-ray emission generated within the inner ∼1 kpc. For most galaxies, the soft X-ray sizes of the sources are comparable to those of the MIR emission. We also find that the hard X-ray faintness previously reported for the bright C-GOALS I sources is also observed in the brightest LIRGs within the sample, with LFIR > 8 × 1010 L⊙.

The dwarf spheroidal galaxy Ursa Major II (UMaII) is believed to be one of the most dark-matter dominated systems among the Milky Way satellites and represents a suitable target for indirect dark matter (DM) searches. The MAGIC telescopes carried out a deep observation campaign on UMaII between 2014 and 2016, collecting almost one hundred hours of good-quality data. This campaign enlarges the pool of DM targets observed at very high energy (E 50 GeV) in search for signatures of DM annihilation in the wide mass range between ∼100 GeV and ∼100 TeV. To this end, the data are analyzed with the full likelihood analysis, a method based on the exploitation of the spectral information of the recorded events for an optimal sensitivity to the explored DM models. We obtain constraints on the annihilation cross-section for different channels that are among the most robust and stringent achieved so far at the TeV mass scale from observations of dwarf satellite galaxies.

It is widely believed that the bulk of the Galactic cosmic rays is accelerated in supernova remnants (SNRs). However, no observational evidence of the presence of particles of PeV energies in SNRs has yet been found. The young historical SNR Cassiopeia A (Cas A) appears as one of the best candidates to study acceleration processes.