Very-high energy gamma-ray astronomy: A 23-year success story in high-energy astroparticle physics

Very-high energy (VHE) gamma quanta contribute only a minuscule fraction - below one per million - to the flux of cosmic rays. Nevertheless, being neutral particles they are currently the best "messengers" of processes from the relativistic/ultra-relativistic Universe because they can be extrapolated back to their origin. The window of VHE gamma rays was opened only in 1989 by the Whipple collaboration, reporting the observation of TeV gamma rays from the Crab nebula. After a slow start, this new field of research is now rapidly expanding with the discovery of more than 150 VHE gamma-ray emitting sources. Progress is intimately related with the steady improvement of detectors and rapidly increasing computing power. We give an overview of the early attempts before and around 1989 and the progress after the pioneering work of the Whipple collaboration. The main focus of this article is on the development of experimental techniques for Earth-bound gamma-ray detectors; consequently, more emphasis is given to those experiments that made an initial breakthrough rather than to the successors which often had and have a similar (sometimes even higher) scientific output as the pioneering experiments. The considered energy threshold is about 30 GeV. At lower energies, observations can presently only be performed with balloon or satellite-borne detectors. Irrespective of the stormy experimental progress, the success story could not have been called a success story without a broad scientific output. Therefore we conclude this article with a summary of the scientific rationales and main results achieved over the last two decades.Comment: 45 pages, 38 figures, review prepared for EPJ-H special issue "Cosmic rays, gamma rays and neutrinos: A survey of 100 years of research

MAGIC and H.E.S.S. detect VHE gamma rays from the blazar OT081 for the first time: a deep multiwavelength study

https://pos.sissa.it/395/815/pdfPublished versio

A lower bound on intergalactic magnetic fields from time variability of 1ES 0229+200 from MAGIC and Fermi /LAT observations

Context. Extended and delayed emission around distant TeV sources induced by the effects of propagation of ? ray s through the intergalactic medium can be used for the measurement of the intergalactic magnetic field (IGMF). Aims. We search for delayed GeV emission from the hard-spectrum TeV ?-ray emitting blazar 1ES 0229+200, with the goal of detecting or constraining the IGMF-dependent secondary flux generated during the propagation of TeV ? rays through the intergalactic medium. Methods. We analysed the most recent MAGIC observations over a 5 year time span, and complemented them with historic data of the H.E.S.S. and VERITAS telescopes, along with a 12-year-long exposure of the Fermi/LAT telescope. We used them to trace source evolution in the GeV TeV band over a decade and a half. We used Monte Carlo simulations to predict the delayed secondary ?-ray flux, modulated by the source variability, as revealed by TeV-band observations. We then compared these predictions for various assumed IGMF strengths to all available measurements of the ?-ray flux evolution. Results. We find that the source flux in the energy range above 200 GeV experiences variations around its average on the 14-year time span of observations. No evidence for the flux variability is found in the 1100 GeV energy range accessible to Fermi/LAT. The non-detection of variability due to delayed emission from electromagnetic cascade developing in the intergalactic medium imposes a lower bound of B > 1.8 - 1017 G for the long-correlation-length IGMF and B > 1014 G for an IGMF of cosmological origin. Though weaker than the one previously derived from the analysis of Fermi/LAT data, this bound is more robust, being based on a conservative intrinsic source spectrum estimate and accounting for the details of source variability in the TeV energy band. We discuss implications of this bound for cosmological magnetic fields that might explain the baryon asymmetry of the Universe.ISSN:0004-6361ISSN:1432-074

MAGIC very large zenith angle observations of the Crab Nebula up to 100 TeV

Aims. We measure the Crab Nebula γ-ray spectral energy distribution in the ~100 TeV energy domain and test the validity of existing leptonic emission models at these high energies. Methods. We used the novel very large zenith angle observations with the MAGIC telescope system to increase the collection area above 10 TeV. We also developed an auxiliary procedure of monitoring atmospheric transmission in order to assure proper calibration of the accumulated data. This employs recording optical images of the stellar field next to the source position, which provides a better than 10% accuracy for the transmission measurements. Results. We demonstrate that MAGIC very large zenith angle observations yield a collection area larger than a square kilometer. In only ~ 56 h of observations, we detect the γ-ray emission from the Crab Nebula up to 100 TeV, thus providing the highest energy measurement of this source to date with Imaging Atmospheric Cherenkov Telescopes. Comparing accumulated and archival MAGIC and Fermi/LAT data with some of the existing emission models, we find that none of them provides an accurate description of the 1 GeV to 100 TeV γ-ray signal.ISSN:0004-6361ISSN:1432-074

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 γγ 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 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 μm range relative to the studied models, yet compatible with them within systematics.ISSN:0035-8711ISSN:1365-296

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

ISSN:0035-8711ISSN:1365-296

Deep observations of the globular cluster M15 with the MAGIC telescopes

A population of globular clusters (GCs) has been recently established by the Fermi-LAT telescope as a new class of GeV γ-ray sources. Leptons accelerated to TeV energies, in the inner magnetospheres of MSPs or in their wind regions, should produce γ-rays through the inverse Compton scattering in the dense radiation field from the huge population of stars. We have conducted deep observations of the GC M15 with the MAGIC telescopes and used 165 h in order to search for γ-ray emission. A strong upper limit on the TeV γ-ray flux <3.2×10−¹³cm−²s−¹ above 300 GeV (<0.26 per cent of the Crab nebula flux) has been obtained. We interpret this limit as a constraint on the efficiency of the acceleration of leptons in the magnetospheres of the MSPs. We constrain the injection rate of relativistic leptons, ηₑ, from the MSPs magnetospheres and their surrounding. We conclude that ηₑ must be lower than expected from the modelling of high-energy processes in MSP inner magnetospheres. For leptons accelerated with the power-law spectrum in the MSP wind regions, ηₑ is constrained to be much lower than derived for the wind regions around classical pulsars. These constraints are valid for the expected range of magnetic field strengths within the GC and for the range of likely energies of leptons injected from the inner magnetospheres, provided that the leptons are not removed from the GC very efficiently due to advection process. We discuss consequences of these constraints for the models of radiation processes around millisecond pulsars.ISSN:0035-8711ISSN:1365-296

MAGIC observations of the diffuse γ-ray emission in the vicinity of the Galactic center

Aims. In the presence of a sufficient amount of target material, γ-rays can be used as a tracer in the search for sources of Galactic cosmic rays (CRs). Here we present deep observations of the Galactic center (GC) region with the MAGIC telescopes and use them to infer the underlying CR distribution and to study the alleged PeV proton accelerator at the center of our Galaxy. Methods. We used data from ≈100 h observations of the GC region conducted with the MAGIC telescopes over five years (from 2012 to 2017). Those were collected at high zenith angles (58−70 deg), leading to a larger energy threshold, but also an increased effective collection area compared to low zenith observations. Using recently developed software tools, we derived the instrument response and background models required for extracting the diffuse emission in the region. We used existing measurements of the gas distribution in the GC region to derive the underlying distribution of CRs. We present a discussion of the associated biases and limitations of such an approach. Results. We obtain a significant detection for all four model components used to fit our data (Sgr A*, “Arc”, G0.9+0.1, and an extended component for the Galactic Ridge). We observe no significant difference between the γ-ray spectra of the immediate GC surroundings, which we model as a point source (Sgr A*) and the Galactic Ridge. The latter can be described as a power-law with index 2 and an exponential cut-off at around 20 TeV with the significance of the cut-off being only 2σ. The derived cosmic-ray profile hints to a peak at the GC position and with a measured profile index of 1.2 ± 0.3 is consistent with the 1/r radial distance scaling law, which supports the hypothesis of a CR accelerator at the GC. We argue that the measurements of this profile are presently limited by our knowledge of the gas distribution in the GC vicinity.ISSN:0004-6361ISSN:1432-074

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 J073326.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 EHBL. In order to investigate the physical nature of this extreme emission, different theoretical frameworks were tested to model the broad-band 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-band SED.ISSN:0035-8711ISSN:1365-296

Investigating the Blazar TXS 0506+056 through Sharp Multiwavelength Eyes During 2017–2019

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted reuse, distribution, and reproduction in any medium, provided the original work is properly cited.--Full list of authors: Acciari, V. A.; Aniello, T.; Ansoldi, S.; Antonelli, L. A.; Arbet Engels, A.; Artero, M.; Asano, K.; Baack, D.; Babić, A.; Baquero, A.; Barres de Almeida, U.; Barrio, J. A.; Batković, I.; Becerra González, J.; Bednarek, W.; Bernardini, E.; Bernardos, M.; Berti, A.; Besenrieder, J.; Bhattacharyya, W.; Bigongiari, C.; Biland, A.; Blanch, O.; Bökenkamp, H.; Bonnoli, G.; Bošnjak, Ž.; Busetto, G.; Carosi, R.; Ceribella, G.; Cerruti, M.; Chai, Y.; Chilingarian, A.; Cikota, S.; Colombo, E.; Contreras, J. L.; Cortina, J.; Covino, S.; D'Amico, G.; D'Elia, V.; Vela, P. Da; Dazzi, F.; De Angelis, A.; De Lotto, B.; Del Popolo, A.; Delfino, M.; Delgado, J.; Mendez, C. Delgado; Depaoli, D.; Di Pierro, F.; Di Venere, L.; Do Souto Espiñeira, E.; Dominis Prester, D.; Donini, A.; Dorner, D.; Doro, M.; Elsaesser, D.; Fallah Ramazani, V.; Fariña, L.; Fattorini, A.; Font, L.; Fruck, C.; Fukami, S.; Fukazawa, Y.; García López, R. J.; Garczarczyk, M.; Gasparyan, S.; Gaug, M.; Giglietto, N.; Giordano, F.; Gliwny, P.; Godinović, N.; Green, J. G.; Green, D.; Hadasch, D.; Hahn, A.; Hassan, T.; Heckmann, L.; Herrera, J.; Hoang, J.; Hrupec, D.; Hütten, M.; Inada, T.; Iotov, R.; Ishio, K.; Iwamura, Y.; Jiménez Martínez, I.; Jormanainen, J.; Jouvin, L.; Kerszberg, D.; Kobayashi, Y.; Kubo, H.; Kushida, J.; Lamastra, A.; Lelas, D.; Leone, F.; Lindfors, E.; Linhoff, L.; Lombardi, S.; Longo, F.; López-Coto, R.; López-Moya, M.; López-Oramas, A.; Loporchio, S.; Machado de Oliveira Fraga, B.; Maggio, C.; Majumdar, P.; Makariev, M.; Mallamaci, M.; Maneva, G.; Manganaro, M.; Mannheim, K.; Mariotti, M.; Martínez, M.; Mas Aguilar, A.; Mazin, D.; Menchiari, S.; Mender, S.; Mićanović, S.; Miceli, D.; Miener, T.; Miranda, J. M.; Mirzoyan, R.; Molina, E.; Moralejo, A.; Morcuende, D.; Moreno, V.; Moretti, E.; Nakamori, T.; Nava, L.; Neustroev, V.; Nievas Rosillo, M.; Nigro, C.; Nilsson, K.; Nishijima, K.; Noda, K.; Nozaki, S.; Ohtani, Y.; Oka, T.; Otero-Santos, J.; Paiano, S.; Palatiello, M.; Paneque, D.; Paoletti, R.; Paredes, J. M.; Pavletić, L.; Peñil, P.; Persic, M.; Pihet, M.; Prada Moroni, P. G.; Prandini, E.; Priyadarshi, C.; Puljak, I.; Rhode, W.; Ribó, M.; Rico, J.; Righi, C.; Rugliancich, A.; Sahakyan, N.; Saito, T.; Sakurai, S.; Satalecka, K.; Saturni, F. G.; Schleicher, B.; Schmidt, K.; Schmuckermaier, F.; Schweizer, T.; Sitarek, J.; Šnidarić, I.; Sobczynska, D.; Spolon, A.; Stamerra, A.; Strišković, J.; Strom, D.; Strzys, M.; Suda, Y.; Surić, T.; Takahashi, M.; Takeishi, R.; Tavecchio, F.; Temnikov, P.; Terzić, T.; Teshima, M.; Tosti, L.; Truzzi, S.; Tutone, A.; Ubach, S.; van Scherpenberg, J.; Vanzo, G.; Vazquez Acosta, M.; Ventura, S.; Verguilov, V.; Viale, I.; Vigorito, C. F.; Vitale, V.; Vovk, I.; Will, M.; Wunderlich, C.; Yamamoto, T.; Zarić, D.; Hodges, M.; Hovatta, T.; Kiehlmann, S.; Liodakis, I.; Max-Moerbeck, W.; Pearson, T. J.; Readhead, A. C. S.; Reeves, R. A.; Lähteenmäki, A.; Tornikoski, M.; Tammi, J.; D'Ammando, F.; Marchini, A.; MAGIC Collaboration; OVRO Collaboration; Metsahovi Collaboration.The blazar TXS 0506+056 got into the spotlight of the astrophysical community in 2017 September, when a high-energy neutrino detected by IceCube (IceCube-170922A) was associated at the 3σ level with a γ-ray flare from this source. This multi-messenger photon-neutrino association remains, as per today, the most significant association ever observed. TXS 0506+056 was a poorly studied object before the IceCube-170922A event. To better characterize its broadband emission, we organized a multiwavelength campaign lasting 16 months (2017 November to 2019 February), covering the radio band (Metsähovi, OVRO), the optical/UV (ASAS-SN, KVA, REM, Swift/UVOT), the X-rays (Swift/XRT, NuSTAR), the high-energy γ rays (Fermi/LAT), and the very high-energy (VHE) γ rays (MAGIC). In γ rays, the behavior of the source was significantly different from the behavior in 2017: MAGIC observations show the presence of flaring activity during 2018 December, while the source only shows an excess at the 4σ level during the rest of the campaign (74 hr of accumulated exposure); Fermi/LAT observations show several short (on a timescale of days to a week) flares, different from the long-term brightening of 2017. No significant flares are detected at lower energies. The radio light curve shows an increasing flux trend that is not seen in other wavelengths. We model the multiwavelength spectral energy distributions in a lepto-hadronic scenario, in which the hadronic emission emerges as Bethe-Heitler and pion-decay cascade in the X-rays and VHE γ rays. According to the model presented here, the 2018 December γ-ray flare was connected to a neutrino emission that was too brief and not bright enough to be detected by current neutrino instruments. © 2022. The Author(s). Published by the American Astronomical Society.The financial support of the German BMBF, MPG, and HGF; the Italian INFN and INAF; the Swiss National Fund SNF; the ERDF under the Spanish Ministerio de Ciencia e Innovación (MICINN) (PID2019-104114RB-C31, PID2019-104114RB-C32, PID2019-104114RB-C33, PID2019-105510GB-C31,PID2019-107847RB-C41, PID2019-107847RB-C42, PID2019-107847RB-C44, PID2019-107988GB-C22); the Indian Department of Atomic Energy; the Japanese ICRR, the University of Tokyo, JSPS, and MEXT; the Bulgarian Ministry of Education and Science, National RI Roadmap Project DO1-400/18.12.2020 and the Academy of Finland grant nr. 320045 is gratefully acknowledged. This work was also supported by the Spanish Centro de Excelencia "Severo Ochoa" (SEV-2016-0588, SEV-2017-0709, CEX2019-000920-S), the Unidad de Excelencia "María de Maeztu" (CEX2019-000918-M, MDM-2015-0509-18-2) and by the CERCA program of the Generalitat de Catalunya; by the Croatian Science Foundation (HrZZ) Project IP-2016-06-9782 and the University of Rijeka Project 13.12.1.3.02; by the DFG Collaborative Research Centers SFB823/C4 and SFB876/C3; the Polish National Research Centre grant UMO-2016/22/M/ST9/00382; and by the Brazilian MCTIC, CNPq and FAPERJ. E.P. acknowledges funding from the Italian Ministry of Education, University and Research (MIUR) through the "Dipartimenti di eccellenza" project Science of the universe. M.C. has received financial support through the Postdoctoral Junior Leader Fellowship Programme from la Caixa Banking Foundation, grant No. LCF/BQ/LI18/11630012. This research has made use of data from the OVRO 40 m monitoring program, which was supported in part by NASA grants NNX08AW31G, NNX11A043G, and NNX14AQ89G, and NSF grants AST-0808050 and AST-1109911, and private funding from Caltech and the MPIfR. This publication makes use of data obtained at Metsähovi Radio Observatory, operated by Aalto University in Finland.Peer reviewe