Gamma-ray astronomy has become one of the main experimental ways to test the modified dispersion relations (MDRs) of photons in vacuum, obtained in some attempts to formulate a theory of Quantum ...Gravity. The MDRs in use imply time delays which depend on the energy, and which increase with distance following some function of redshift. The use of transient, or variable, distant and highly energetic sources, already allows us to set stringent limits on the energy scale related to this phenomenon, usually thought to be of the order of the Planck energy, but robust conclusions on the existence of MDR-related propagation effects still require the analysis of a large population of sources. In order to gather the biggest sample of sources possible for MDR searches at teraelectronvolt energies, the H.E.S.S., MAGIC and VERITAS collaborations enacted a joint task force to combine all their relevant data to constrain the Quantum Gravity energy scale. In the present article, the likelihood method used, to combine the data and provide a common limit, is described in detail and tested through simulations of recorded data sets for a gamma-ray burst, three flaring active galactic nuclei and two pulsars. Statistical and systematic errors are assessed and included in the likelihood as nuisance parameters. In addition, a comparison of two different formalisms for distance dependence of the time lags is performed for the first time. In a second article, to appear later, the method will be applied on all relevant data from the three experiments.
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.
The Lorentz Invariance Violation (LIV), a proposed consequence of certain quantum gravity (QG) scenarios, could instigate an energy-dependent group velocity for ultra-relativistic particles. This ...energy dependence, although suppressed by the massive QG energy scale \(E_\mathrm{QG}\), expected to be on the level of the Planck energy \(1.22 \times 10^{19}\) GeV, is potentially detectable in astrophysical observations. In this scenario, the cosmological distances traversed by photons act as an amplifier for this effect. By leveraging the observation of a remarkable flare from the blazar Mrk\,421, recorded at energies above 100 GeV by the MAGIC telescopes on the night of April 25 to 26, 2014, we look for time delays scaling linearly and quadratically with the photon energies. Using for the first time in LIV studies a binned-likelihood approach we set constraints on the QG energy scale. For the linear scenario, we set \(95\%\) lower limits \(E_\mathrm{QG}>2.7\times10^{17}\) GeV for the subluminal case and \(E_\mathrm{QG}> 3.6 \times10^{17}\) GeV for the superluminal case. For the quadratic scenario, the \(95\%\) lower limits for the subluminal and superluminal cases are \(E_\mathrm{QG}>2.6 \times10^{10}\) GeV and \(E_\mathrm{QG}>2.5\times10^{10}\) GeV, respectively.
Clusters of galaxies are the largest known gravitationally bound structures in the Universe, with masses around 1015 M⊙, most of it in the form of dark matter. The ground-based Imaging Atmospheric ...Cherenkov Telescope MAGIC made a deep survey of the Perseus cluster of galaxies using almost 400 h of data recorded between 2009 and 2017. This is the deepest observational campaign so far on a cluster of galaxies in the very high energy range. We search for gamma-ray signals from dark matter particles in the mass range between 200 GeV and 200 TeV decaying into standard model pairs. We apply an analysis optimized for the spectral and morphological features expected from dark matter decays and find no evidence of decaying dark matter. From this, we conclude that dark matter particles have a decay lifetime longer than ∼1026 s in all considered channels. Our results improve previous lower limits found by MAGIC and represent the strongest limits on decaying dark matter particles from ground-based gamma-ray instruments.
The MAGIC stereoscopic system collected 69 hours of Crab Nebula data between October 2009 and April 2011. Analysis of this data sample using the latest improvements in the MAGIC stereoscopic software ...provided an unprecedented precision of spectral and night-by-night light curve determination at gamma rays. We derived a differential spectrum with a single instrument from 50 GeV up to almost 30 TeV with 5 bins per energy decade. At low energies, MAGIC results, combined with Fermi-LAT data, show a flat and broad Inverse Compton peak. The overall fit to the data between 1 GeV and 30 TeV is not well described by a log-parabola function. We find that a modified log-parabola function with an exponent of 2.5 instead of 2 provides a good description of the data (χred2=35/26). Using systematic uncertainties of the MAGIC and Fermi-LAT measurements we determine the position of the Inverse Compton peak to be at (53±3stat+31syst−13syst) GeV, which is the most precise estimation up to date and is dominated by the systematic effects. There is no hint of the integral flux variability on daily scales at energies above 300 GeV when systematic uncertainties are included in the flux measurement. We consider three state-of-the-art theoretical models to describe the overall spectral energy distribution of the Crab Nebula. The constant B-field model cannot satisfactorily reproduce the VHE spectral measurements presented in this work, having particular difficulty reproducing the broadness of the observed IC peak. Most probably this implies that the assumption of the homogeneity of the magnetic field inside the nebula is incorrect. On the other hand, the time-dependent 1D spectral model provides a good fit of the new VHE results when considering a 80 μG magnetic field. However, it fails to match the data when including the morphology of the nebula at lower wavelengths.
In recent years, a new generation of optical intensity interferometers has emerged, leveraging the existing infrastructure of Imaging Atmospheric Cherenkov Telescopes (IACTs). The MAGIC telescopes ...host the MAGIC-SII system (Stellar Intensity Interferometer), implemented to investigate the feasibility and potential of this technique on IACTs. After the first successful measurements in 2019, the system was upgraded and now features a real-time, dead-time-free, 4-channel, GPU-based correlator. These hardware modifications allow seamless transitions between MAGIC's standard very-high-energy gamma-ray observations and optical interferometry measurements within seconds. We establish the feasibility and potential of employing IACTs as competitive optical Intensity Interferometers with minimal hardware adjustments. The measurement of a total of 22 stellar diameters are reported, 9 corresponding to reference stars with previous comparable measurements, and 13 with no prior measurements. A prospective implementation involving telescopes from the forthcoming Cherenkov Telescope Array Observatory's northern hemisphere array, such as the first prototype of its Large-Sized Telescopes, LST-1, is technically viable. This integration would significantly enhance the sensitivity of the current system and broaden the UV-plane coverage. This advancement would enable the system to achieve competitive sensitivity with the current generation of long-baseline optical interferometers over blue wavelengths.
We present the first multi-wavelength study of Mrk 501 including
very-high-energy (VHE) gamma-ray observations simultaneous to X-ray
polarization measurements from the Imaging X-ray Polarimetry ...Explorer (IXPE).
We use radio-to-VHE data from a multi-wavelength campaign organized between
2022-03-01 and 2022-07-19. The observations were performed by MAGIC, Fermi-LAT,
NuSTAR, Swift (XRT and UVOT), and several instruments covering the optical and
radio bands. During the IXPE pointings, the VHE state is close to the average
behavior with a 0.2-1 TeV flux of 20%-50% the emission of the Crab Nebula.
Despite the average VHE activity, an extreme X-ray behavior is measured for the
first two IXPE pointings in March 2022 with a synchrotron peak frequency >1
keV. For the third IXPE pointing in July 2022, the synchrotron peak shifts
towards lower energies and the optical/X-ray polarization degrees drop. The
X-ray polarization is systematically higher than at lower energies, suggesting
an energy-stratification of the jet. While during the IXPE epochs the
polarization angle in the X-ray, optical and radio bands align well, we find a
clear discrepancy in the optical and radio polarization angles in the middle of
the campaign. We model the broad-band spectra simultaneous to the IXPE
pointings assuming a compact zone dominating in the X-rays and VHE, and an
extended zone stretching further downstream the jet dominating the emission at
lower energies. NuSTAR data allow us to precisely constrain the synchrotron
peak and therefore the underlying electron distribution. The change between the
different states observed in the three IXPE pointings can be explained by a
change of magnetization and/or emission region size, which directly connects
the shift of the synchrotron peak to lower energies with the drop in
polarization degree.
Axion-like particles (ALPs) are pseudo-Nambu-Goldstone bosons that emerge in various theories beyond the standard model. These particles can interact with high-energy photons in external magnetic ...fields, influencing the observed gamma-ray spectrum. This study analyzes 41.3 hrs of observational data from the Perseus Galaxy Cluster collected with the MAGIC telescopes. We focused on the spectra the radio galaxy in the center of the cluster: NGC 1275. By modeling the magnetic field surrounding this target, we searched for spectral indications of ALP presence. Despite finding no statistical evidence of ALP signatures, we were able to exclude ALP models in the sub-micro electronvolt range. Our analysis improved upon previous work by calculating the full likelihood and statistical coverage for all considered models across the parameter space. Consequently, we achieved the most stringent limits to date for ALP masses around 50 neV, with cross sections down to \(g_{a\gamma} = 3 \times 10^{-12}\) GeV\(^{-1}\).
Monthly Notices of the Royal Astronomical Society, Volume 527,
Issue 3, January 2024, Pages 5856-5867 Gamma-ray bursts (GRBs) are explosive transient events occurring at
cosmological distances, ...releasing a large amount of energy as electromagnetic
radiation over several energy bands. We report the detection of the long
GRB~201216C by the MAGIC telescopes. The source is located at $z=1.1$ and thus
it is the farthest one detected at very high energies. The emission above
\SI{70}{\GeV} of GRB~201216C is modelled together with multi-wavelength data
within a synchrotron and synchrotron-self Compton (SSC) scenario. We find that
SSC can explain the broadband data well from the optical to the
very-high-energy band. For the late-time radio data, a different component is
needed to account for the observed emission. Differently from previous GRBs
detected in the very-high-energy range, the model for GRB~201216C strongly
favors a wind-like medium. The model parameters have values similar to those
found in past studies of the afterglows of GRBs detected up to GeV energies.