Linelike features in TeV $\gamma$ rays constitute a “smoking gun” for TeV-scale particle dark matter and new physics. Probing the Galactic Center region with ground-based Cherenkov telescopes enables the search for TeV spectral features in immediate association with a dense dark matter reservoir at a sensitivity out of reach for satellite $\gamma$-ray detectors, and direct detection and collider experiments. We report on 223 hours of observations of the Galactic Center region with the MAGIC stereoscopic telescope system reaching $\gamma$-ray energies up to 100 TeV. We improved the sensitivity to spectral lines at high energies using large-zenith-angle observations and a novel background modeling method within a maximum-likelihood analysis in the energy domain. No linelike spectral feature is found in our analysis. Therefore, we constrain the cross section for dark matter annihilation into two photons to $⟨\sigma v⟩\lesssim 5\times {10}^{-28}{\mathrm{cm}}^3{\mathrm{s}}^{-1}$ at 1 TeV and $⟨\sigma v⟩\lesssim 1\times {10}^{-25}{\mathrm{cm}}^3{\mathrm{s}}^{-1}$ at 100 TeV, achieving the best limits to date for a dark matter mass above 20 TeV and a cuspy dark matter profile at the Galactic Center. Finally, we use the derived limits for both cuspy and cored dark matter profiles to constrain supersymmetric wino models.

Monochromatic gamma-ray signals constitute a potential smoking gun signature for annihilating or decaying dark matter particles that could relatively easily be distinguished from astrophysical or instrumental backgrounds. We provide an updated assessment of the sensitivity of the Cherenkov Telescope Array (CTA) to such signals, based on observations of the Galactic centre region as well as of selected dwarf spheroidal galaxies. We find that current limits and detection prospects for dark matter masses above 300 GeV will be significantly improved, by up to an order of magnitude in the multi-TeV range. This demonstrates that CTA will set a new standard for gamma-ray astronomy also in this respect, as the world's largest and most sensitive high-energy gamma-ray observatory, in particular due to its exquisite energy resolution at TeV energies and the adopted observational strategy focussing on regions with large dark matter densities. Throughout our analysis, we use up-to-date instrument response functions, and we thoroughly model the effect of instrumental systematic uncertainties in our statistical treatment. We further present results for other potential signatures with sharp spectral features, e.g. box-shaped spectra, that would likewise very clearly point to a particle dark matter origin.

We study the broadband emission of Mrk 501 using multiwavelength observations from 2017 to 2020 performed with a multitude of instruments, involving, among others, MAGIC, Fermi's Large Area Telescope (LAT), NuSTAR, Swift, GASP-WEBT, and the Owens Valley Radio Observatory. Mrk 501 showed an extremely low broadband activity, which may help to unravel its baseline emission. Nonetheless, significant flux variations are detected at all wave bands, with the highest occurring at X-rays and very-high-energy (VHE) γ-rays. A significant correlation (>3σ) between X-rays and VHE γ-rays is measured, supporting leptonic scenarios to explain the variable parts of the emission, also during low activity. This is further supported when we extend our data from 2008 to 2020, and identify, for the first time, significant correlations between the Swift X-Ray Telescope and Fermi-LAT. We additionally find correlations between high-energy γ-rays and radio, with the radio lagging by more than 100 days, placing the γ-ray emission zone upstream of the radio-bright regions in the jet. Furthermore, Mrk 501 showed a historically low activity in X-rays and VHE γ-rays from mid-2017 to mid-2019 with a stable VHE flux (>0.2 TeV) of 5% the emission of the Crab Nebula. The broadband spectral energy distribution (SED) of this 2 yr long low state, the potential baseline emission of Mrk 501, can be characterized with one-zone leptonic models, and with (lepto)-hadronic models fulfilling neutrino flux constraints from IceCube. We explore the time evolution of the SED toward the low state, revealing that the stable baseline emission may be ascribed to a standing shock, and the variable emission to an additional expanding or traveling shock.

Context.Certain types of supernova remnants (SNRs) in our Galaxy are assumed to be PeVatrons, capable of accelerating cosmic rays (CRs) to ~ PeV energies. However, conclusive observational evidence for this has not yet been found. The SNR G106.3+2.7, detected at 1–100 TeV energies by different γ-ray facilities, is one of the most promising PeVatron candidates. This SNR has a cometary shape, which can be divided into a head and a tail region with different physical conditions. However, in which region the 100 TeV emission is produced has not yet been identified because of the limited position accuracy and/or angular resolution of existing observational data. Additionally, it remains unclear as to whether the origin of the γ-ray emission is leptonic or hadronic.Aims.With the better angular resolution provided by new MAGIC data compared to earlierγ-ray datasets, we aim to reveal the acceleration site of PeV particles and the emission mechanism by resolving the SNR G106.3+2.7 with 0.1° resolution at TeV energies.Methods.We observed the SNR G106.3+2.7 using the MAGIC telescopes for 121.7 h in total – after quality cuts – between May 2017 and August 2019. The analysis energy threshold is ~0.2 TeV, and the angular resolution is 0.07−0.1°. We examined theγ-ray spectra of different parts of the emission, whilst benefitting from the unprecedented statistics and angular resolution at these energies provided by our new data. We also used measurements at other wavelengths such as radio, X-rays, GeVγ-rays, and 10 TeVγ-rays to model the emission mechanism precisely.Results.We detect extended γ-ray emission spatially coincident with the radio continuum emission at the head and tail of SNR G106.3+2.7. The fact that we detect a significantγ-ray emission with energies above 6.0 TeV from only the tail region suggests that the emissions above 10 TeV detected with air shower experiments (Milagro, HAWC, Tibet ASγand LHAASO) are emitted only from the SNR tail. Under this assumption, the multi-wavelength spectrum of the head region can be explained with either hadronic or leptonic models, while the leptonic model for the tail region is in contradiction with the emission above 10 TeV and X-rays. In contrast, the hadronic model could reproduce the observed spectrum at the tail by assuming a proton spectrum with a cutoff energy of ~1 PeV for that region. Such high-energy emission in this middle-aged SNR (4−10 kyr) can be explained by considering a scenario where protons escaping from the SNR in the past interact with surrounding dense gases at present.Conclusions.Theγ-ray emission region detected with the MAGIC telescopes in the SNR G106.3+2.7 is extended and spatially coincident with the radio continuum morphology. The multi-wavelength spectrum of the emission from the tail region suggests proton acceleration up to ~PeV, while the emission mechanism of the head region could either be hadronic or leptonic.