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 ∼ 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–30% 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 ∼ 18%, under dark nights only, to up to ∼ 40% in total with only moderate performance degradation.
Clusters of galaxies are expected to be reservoirs of cosmic rays (CRs) that should produce diffuse γ-ray emission due to their hadronic interactions with the intra-cluster medium. The nearby Perseus ...cool-core cluster, identified as the most promising target to search for such an emission, has been observed with the MAGIC telescopes at very-high energies (VHE, E ≳ 100 GeV) for a total of 253 h from 2009 to 2014. The active nuclei of NGC 1275, the central dominant galaxy of the cluster, and IC 310, lying at about 0.6° from the centre, have been detected as point-like VHE γ-ray emitters during the first phase of this campaign. We report an updated measurement of the NGC 1275 spectrum, which is described well by a power law with a photon index Γ = 3.6 ± 0.2stat ± 0.2syst between 90 GeV and 1200 GeV. We do not detect any diffuse γ-ray emission from the cluster and so set stringent constraints on its CR population. To bracket the uncertainties over the CR spatial and spectral distributions, we adopt different spatial templates and power-law spectral indexes α. For α = 2.2, the CR-to-thermal pressure within the cluster virial radius is constrained to be ≲ 1−2%, except if CRs can propagate out of the cluster core, generating a flatter radial distribution and releasing the CR-to-thermal pressure constraint to ≲ 20%. Assuming that the observed radio mini-halo of Perseus is generated by secondary electrons from CR hadronic interactions, we can derive lower limits on the central magnetic field, B0, that depend on the CR distribution. For α = 2.2, B0 ≳ 5−8 μG, which is below the ~25 μG inferred from Faraday rotation measurements, whereas for α ≲ 2.1, the hadronic interpretation of the diffuse radio emission contrasts with our γ-ray flux upper limits independently of the magnetic field strength.
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-11 cm-2 s-1 to the maximum flux, that is about 24 times higher, at (3.1 ± 0.1)×10-10 cm-2 s-1. 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.
Context. Blazars are variable sources on various timescales over a broad energy range spanning from radio to very high energy (>100 GeV, hereafter VHE). Mrk 501 is one of the brightest blazars at TeV ...energies and has been extensively studied since its first VHE detection in 1996. However, most of the γ-ray studies performed on Mrk 501 during the past years relate to flaring activity, when the source detection and characterization with the available γ-ray instrumentation was easier toperform. Aims. Our goal is to characterize the source γ-ray emission in detail, together with the radio-to-X-ray emission, during the non-flaring (low) activity, which is less often studied than the occasional flaring (high) activity. Methods. We organized a multiwavelength (MW) campaign on Mrk 501 between March and May 2008. This multi-instrument effort included the most sensitive VHE γ-ray instruments in the northern hemisphere, namely the imaging atmospheric Cherenkov telescopes MAGIC and VERITAS, as well as Swift, RXTE, the F-GAMMA, GASP-WEBT, and other collaborations and instruments. This provided extensive energy and temporal coverage of Mrk 501 throughout the entire campaign. Results. Mrk 501 was found to be in a low state of activity during the campaign, with a VHE flux in the range of 10%–20% of the Crab nebula flux. Nevertheless, significant flux variations were detected with various instruments, with a trend of increasing variability with energy and a tentative correlation between the X-ray and VHE fluxes. The broadband spectral energy distribution during the two different emission states of the campaign can be adequately described within the homogeneous one-zone synchrotron self-Compton model, with the (slightly) higher state described by an increase in the electron number density. Conclusions. The one-zone SSC model can adequately describe the broadband spectral energy distribution of the source during the two months covered by the MW campaign. This agrees with previous studies of the broadband emission of this source during flaring and non-flaring states. We report for the first time a tentative X-ray-to-VHE correlation during such a low VHE activity. Although marginally significant, this positive correlation between X-ray and VHE, which has been reported many times during flaring activity, suggests that the mechanisms that dominate the X-ray/VHE emission during non-flaring-activity are not substantially different from those that are responsible for the emission during flaring activity.
Aims. The radio galaxy NGC 1275, recently identified as a very high energy (VHE, >100 GeV) γ-ray emitter by MAGIC, is one of the few non-blazar active galactic nuclei detected in the VHE regime. The ...purpose of this work is to better understand the origin of the γ-ray emission and locate it within the galaxy. Methods. We studied contemporaneous multifrequency observations of NGC 1275 and modeled the overall spectral energy distribution. We analyzed unpublished MAGIC observations carried out between October 2009 and February 2010, and the previously published observations taken between August 2010 and February 2011. We studied the multiband variability and correlations by analyzing data of Fermi-LAT in the 100 MeV–100 GeV energy band, as well as Chandra (X-ray), KVA (optical), and MOJAVE (radio) data taken during the same period. Results. Using customized Monte Carlo simulations corresponding to early MAGIC stereoscopic data, we detect NGC 1275 also in the earlier MAGIC campaign. The flux level and energy spectra are similar to the results of the second campaign. The monthly light curve above 100 GeV shows a hint of variability at the 3.6σ level. In the Fermi-LAT band, both flux and spectral shape variabilities are reported. The optical light curve is also variable and shows a clear correlation with the γ-ray flux above 100 MeV. In radio, three compact components are resolved in the innermost part of the jet. One of these components shows a similar trend as the Fermi-LAT and KVA light curves. The γ-ray spectra measured simultaneously with MAGIC and Fermi-LAT from 100 MeV to 650 GeV can be well fitted either by a log-parabola or by a power-law with a subexponential cutoff for the two observation campaigns. A single-zone synchrotron-self-Compton model, with an electron spectrum following a power-law with an exponential cutoff, can explain the broadband spectral energy distribution and the multifrequency behavior of the source. However, this model suggests an untypical low bulk-Lorentz factor or a velocity alignment closer to the line of sight than the parsec-scale radio jet.
Context. We present the results of a multi-year monitoring campaign of the Galactic center (GC) with the MAGIC telescopes. These observations were primarily motivated by reports that a putative gas ...cloud (G2) would be passing in close proximity to the super-massive black hole (SMBH), associated with Sagittarius A*, located at the center of our galaxy. This event was expected to give astronomers a unique chance to study the effect of in-falling matter on the broad-band emission of a SMBH. Aims. We search for potential flaring emission of very-high-energy (VHE; ≥100 GeV) gamma rays from the direction of the SMBH at the GC due to the passage of the G2 object. Using these data we also study the morphology of this complex region. Methods. We observed the GC region with the MAGIC Imaging Atmospheric Cherenkov Telescopes during the period 2012–2015, collecting 67 h of good-quality data. In addition to a search for variability in the flux and spectral shape of the GC gamma-ray source, we use a point-source subtraction technique to remove the known gamma-ray emitters located around the GC in order to reveal the TeV morphology of the extended emission inside that region. Results. No effect of the G2 object on the VHE gamma-ray emission from the GC was detected during the 4 yr observation campaign. We confirm previous measurements of the VHE spectrum of Sagittarius A*, and do not detect any significant variability of the emission from the source. Furthermore, the known VHE gamma-ray emitter at the location of the supernova remnant G0.9+0.1 was detected, as well as the recently discovered VHE source close to the GG radio arc.
The Geminga pulsar, one of the brighest gamma-ray sources, is a promising candidate for emission of very-high-energy (VHE > 100 GeV) pulsed gamma rays. Also, detection of a large nebula has been ...claimed by water Cherenkov instruments. We performed deep observations of Geminga with the MAGIC telescopes, yielding 63 h of good-quality data, and searched for emission from the pulsar and pulsar wind nebula. We did not find any significant detection, and derived 95% confidence level upper limits. The resulting upper limits of 5.3 × 10-13 TeV cm-2 s-1 for the Geminga pulsar and 3.5 × 10-12 TeV cm-2 s-1 for the surrounding nebula at 50 GeV are the mostconstraining ones obtained so far at VHE. To complement the VHE observations, we also analyzed 5 yr of Fermi-LAT data from Geminga, finding that the sub-exponential cut-off is preferred over the exponential cut-off that has been typically used in the literature. We also find that, above 10 GeV, the gamma-ray spectra from Geminga can be described with a power law with index softer than 5. The extrapolation of the power-law Fermi-LAT pulsed spectra to VHE goes well below the MAGIC upper limits, indicating that the detection of pulsed emission from Geminga with the current generation of Cherenkov telescopes is very difficult.
A search for tau neutrino induced showers with the MAGIC telescopes is presented. The MAGIC telescopes located at an altitude of 2200 m a.s.l. in the Canary Island of La Palma, can point towards the ...horizon or a few degrees below across an azimuthal range of about 80°. This provides a possibility to search for air showers induced by tau leptons arising from interactions of tau neutrinos in the Earth crust or the surrounding ocean. In this paper we show how such air showers can be discriminated from the background of very inclined hadronic showers by using Monte Carlo simulations. Taking into account the orography of the site, the point source acceptance and the event rates expected have been calculated for a sample of generic neutrino fluxes from photo-hadronic interactions in AGNs. The analysis of about 30 h of data taken towards the sea leads to a 90% C.L. point source limit for tau neutrinos in the energy range from 1.0 × 1015 eV to 3.0 × 1018 eV of about Eντ2×ϕ(Eντ)<2.0×10−4 GeV cm−2 s−1 for an assumed power-law neutrino spectrum with spectral index γ=−2. However, with 300 h and in case of an optimistic neutrino flare model, limits of the level down to Eντ2×ϕ(Eντ)<8.4×10−6 GeV cm−2 s−1 can be expected.