Cherenkov Telescope Array Research Articles

Prospects for annihilating dark matter from M31 and M33 observations with the Cherenkov Telescope Array

M31 and M33 are the closest spiral galaxies and the largest members (together with the Milky Way) of the Local group, which makes them interesting targets for indirect dark matter searches. In this paper we present studies of the expected sensitivity of the Cherenkov Telescope Array (CTA) to an annihilation signal from weakly interacting massive particles from M31 and M33. We show that a 100 h long observation campaign will allow CTA to probe annihilation cross-sections up to 〈συ〉 ≈ 5·10-25 cm3 s-1 for the τ + τ - annihilation channel (for M31, at a DM mass of 0.3 TeV), improving the current limits derived by HAWC by up to an order of magnitude. We present an estimate of the expected CTA sensitivity, by also taking into account the contributions of the astrophysical background and other possible sources of systematic uncertainty. We also show that CTA might be able to detect the extended emission from the bulge of M31, detected at lower energies by the Fermi/LAT.

The effect of night sky background on source reconstruction in atmospheric Cherenkov telescope arrays and implications for image cleaning

Abstract Gamma-ray astronomy through imaging atmospheric Cherenkov telescope (IACTs) has enabled the study of numerous astrophysical sources in the GeV-TeV domain. The presence of night sky background (NSB) in IACT images affects the reconstruction of the gamma-ray being observed and is usually dealt through image cleaning. In this work, we propose that the image cleaning approach might be adapted to the reconstruction method and NSB levels. In particular, we study the effect of NSB on source reconstruction for gamma-rays observed through IACT arrays through Monte Carlo simulations. We find that the presence of NSB can have a smoothing effect on the shower images enabling better source reconstruction under certain circumstances. We also find that there is an optimum cleaning threshold depending on the energy of the initial gamma-ray. This threshold shows very little dependence on the NSB level or core position. Based on these results, we propose and test a hybrid image cleaning method for source reconstruction. We find that this hybrid image cleaning method results in better source reconstruction when observations are made in high NSB conditions. Whereas at low NSB levels, the application of a single image cleaning threshold yields better results. We finally conclude that this implies that different image cleaning methods might be needed to best carry out source reconstructions at different NSB levels.

Gamma-ray emission from spectrally resolved cosmic rays in galaxies

ABSTRACT Cosmic rays (CRs) are ubiquitous in the interstellar medium (ISM) of nearby galaxies, but many of their properties are not well-constrained. Gamma-ray observations provide a powerful tool in this respect, allowing us to constrain both the interaction of CR protons with the ISM and their transport properties. To help better understand the link between observational signatures and CR physics, we use a series of magneto-hydrodynamical arepo simulations of isolated galaxies performed using spectrally resolved CR transport in every computational cell, with subsequent gamma-ray emission calculated using the Crayon + (Cosmic RAY emissiON) code. In each of our simulated haloes, modelling the energy-dependent spatial diffusion of CRs leads to a more extended distribution of high-energy (~100 GeV) gamma rays compared to that predicted by a ‘grey’ steady-state model, which is especially visible in the corresponding emission maps and radial profiles. Despite this, the total gamma-ray spectra can often be well approximated by the steady-state model, although recovering the same spectral index typically requires a minor variation of the energy dependence of the diffusion coefficient. Our simulations reproduce the observed shapes of the gamma-ray spectra of nearby star-forming galaxies and also match recent observations of the far infrared–gamma-ray relation. We find, however, that the spectrally resolved model yields marginally smaller luminosities for lower star formation rates compared to grey simulations of CRs. Our work highlights the importance of modelling spectrally resolved CR transport for an accurate prediction of spatially resolved high-energy gamma-ray emission, as will be probed by the upcoming Cherenkov Telescope Array observatory.

Development of a SiPM Pixel prototype for the Large-Sized Telescope of the Cherenkov Telescope Array

The Cherenkov Telescope Array (CTA) will be the next generation ground-based gamma-ray observatory. CTA consists of different telescope types of which the largest ones (Large-Sized Telescopes, LSTs) cover the lower energy range, between 20GeV and 200GeV. The first LST is currently being commissioned at the Roque de los Muchachos Observatory, La Palma, Canary Islands. Its camera has 1855 photomultipliers (PMTs) with 1.5 inch cathodes. Silicon Photomultipliers (SiPMs) are increasingly becoming valid alternatives to PMTs also in gamma-ray astronomy. In the context of the LST project, there is an effort to study a novel Advanced Camera, equipped with SiPMs and a completely redesigned electronics based on a fully digital approach. To study and develop solutions on the sensors of these camera, we built a prototype camera module with a fully re-designed pre-amplifying stage and sensor bias control while re-using the digitizing and triggering stages of the existing LST camera module. We report on the design choices made to achieve the highest performance in terms of timing and charge resolution and the laboratory measurements validating those choices.

Probing the IC/CMB interpretation for the X-ray knots of AGNs through VHE observations

ABSTRACT The detection of hard X-ray spectra (spectral index <2) from the kiloparsec-scale jet of active galactic nuclei cannot be accounted for by the synchrotron emission mechanism from the electron distribution responsible for the radio/optical emission. Alternate explanations are the inverse Compton scattering of cosmic microwave background photons (IC/CMB) or synchrotron emission from a second electron population. When the X-ray emission is interpreted as an IC/CMB process, the Compton spectrum often peaks at GeV energy and many sources are predicted to be Fermi candidate sources. The absence of significant gamma-ray flux from some of these galaxies by Fermi disfavours the IC/CMB interpretation of the high-energy emission. We extend this study to predict the very-high-energy (VHE) gamma-ray emission due to the IC/CMB model, which can be investigated by the Cherenkov Telescope Array Observatory (CTAO). The model parameters deciding the broad-band spectral energy distribution are estimated using an analytical approximation of the emissivity functions. The emission model is extrapolated to VHE and then compared with the CTAO sensitivity. For this particular study, we have selected 18 knots with harder X-ray spectra and for which the IC/CMB model for X-ray emission has been suggested.

Photon-ALP oscillations at CTA energies

Axions and axion-like-particles (ALPs) are well motivated beyond the standard model particles that can explain a variety of unsolved problems in physics, such as the strong CP problem and the nature of dark matter. These particles are characterised by their two-photon coupling, which leads to so-called photon-ALP oscillation as photons propagate through an external magnetic field. Such oscillations lead to characteristic signatures in the energy spectrum of high-energy photons from astrophysical sources, allowing us to probe the existence of ALPs and possibly dark matter. We review the signatures from ALPs in photon spectra and discuss a new method that can be used to directly search for the energy dependence of the oscillations. The focus is on photons at TeV-energies relevant for the upcoming Cherenkov Telescope Array (CTA).

Modeling the propagation of very-high-energy γ-rays with the CRbeam code: Comparison with CRPropa and ELMAG codes

Context. Very-high-energy γ-rays produce electron positron pairs in interactions with low-energy photons of extragalactic background light during propagation through the intergalactic medium. The electron-positron pairs generate secondary γ-rays detectable by γ-ray telescopes. This secondary emission can be used to detect intergalactic magnetic fields (IGMF) in the voids of large-scale structure. Aims. A new γ-ray observatory, namely, Cherenkov Telescope Array (CTA), will provide an increase in sensitivity for detections of these secondary γ-ray emission and enable the measurement of its properties for sources at cosmological distances. The interpretation of the CTA data, including the detection of IGMF and study of its properties and origins, will require precision modeling of the primary and secondary γ-ray fluxes. Methods. We assess the precision of the modeling of the secondary γ-ray emission using model calculations with publicly available Monte Carlo codes CRPropa and ELMAG and compare their predictions with theoretical expectations and with model calculations of a newly developed CRbeam code. Results. We find that model predictions of different codes differ by up to 50% for low-redshift sources, with discrepancies increasing up to order-of-magnitude level with the increasing source redshifts. We identify the origin of these discrepancies and demonstrate that after eliminating the inaccuracies found, the discrepancies between the three codes are reduced to 10% when modeling nearby sources with z ~ 0.1. We argue that the new CRbeam code provides reliable predictions for the spectral, timing, and imaging properties of the secondary γ-ray signal for both nearby and distant sources with z ~ 1. Thus, it can be used to study gamma-ray sources and IGMF with a level of precision that is appropriate for the prospective CTA study of the effects of γ-ray propagation through the intergalactic medium.

Current status and operation of the H.E.S.S. array of imaging atmospheric Cherenkov telescopes

The High Energy Stereoscopic System (H.E.S.S.) is an array of five imaging atmospheric Cherenkov telescopes (IACTs) to study gamma-ray emission from astrophysical objects in the Southern hemisphere. It is the only hybrid array of IACTs, composed of telescopes with different collection area and footprint, individually optimised for a specific energy range. Collectively, the array is most sensitive to gamma rays in the range of 100GeV to 100TeV. The array has been in operation since 2002 and has been upgraded with new telescopes and cameras multiple times. Recent hardware upgrades and changes in the operational procedures increased the amount of observing time, which is of key importance for time-domain science. H.E.S.S. operations saw record data taking in 2020 and 2021 and we describe the current operations with specific emphasis on system performance, operational processes and workflows, quality control, and (near) real-time extraction of science results. In light of this, we will briefly discuss the early detection of gamma-ray emission from the recurrent nova RS Oph and alert distribution to the astrophysics community.

The prototype Schwarzschild–Couder Telescope: A medium-sized telescope for the Cherenkov Telescope Array

The Schwarzschild–Couder Telescope (SCT) is a dual mirror medium-sized telescope proposed for the Cherenkov Telescope Array (CTA), the next-generation very-high energy (from about 20 GeV to 300 TeV) gamma-ray observatory. The SCT design consists of a dual-mirror optics and a high resolution camera with a field of view (FoV) of 8 degrees squared, which will allow exceptional performance in terms of angular resolution and background rejection. A prototype telescope (named pSCT) has been installed at the Fred Lawrence Whipple Observatory in Arizona, USA. Its camera is partially equipped and covers a FoV of 2.7°. The pSCT has recently successfully detected the Crab Nebula with a statistical significance of 8.6 standard deviations. The upgrade of the pSCT focal plane is now ongoing, aimed to equip the full camera with upgraded sensors and electronics, enhancing the telescope field of view from the current 2.7°to the final 8°. In this presentation, an overview of the pSCT project and obtained results will be given, together with the camera upgrade status and expected performance.

The silicon photomultiplier-based camera for the Cherenkov Telescope Array small-sized telescopes

The Cherenkov Telescope Array (CTA) is the next generation ground-based gamma-ray astronomy observatory, planned to comprise two arrays of imaging air Cherenkov telescopes (IACTs) located in the northern and southern hemispheres. Three telescope sizes are required to cover the CTA gamma-ray energy range from 20 GeV to 300 TeV. An array of several tens of Small-Sized Telescopes (SSTs) at the southern site situated in the Andes at Paranal in Chile, will provide unprecedented sensitivity above 1 TeV and up to 300 TeV, and offer the highest angular resolution of any instrument at these energies. Following a down selection from three prototype telescopes, the design finally selected for SST comprises a dual mirror Schwarzschild–Couder optic with a 4.3 m diameter primary mirror and a 1.8 m secondary mirror imaged by a silicon photomultiplier (SiPM)-based camera with a ∼9°field of view (FoV). The dual mirror optics produce a smaller plate-scale aplanatic focal plane allowing a small, low-cost camera to be employed, compared to that required for the conventional single mirror Davies-Cotton IACT design. The camera comprises an array of 2048 SiPM pixels, configured as 32 sensor and electronics modules each with an 8 × 8 pixel2 tile populated with 6 × 6 mm2 SiPM pixels. Full waveform capture on every channel is provided by the TARGET ASIC which performs the dual function of event triggering and waveform digitization of the full camera at 1 GSample/s. We describe the finalized SST camera design including its optimization for the production phase of the project anticipated to begin in 2023.

The NectarCAM timing system

NectarCAM is a Cherenkov camera which is going to equip the Medium-Sized Telescopes (MST) of the northern site of the Cherenkov Telescope Array Observatory (CTAO). NectarCAM is equipped with 265 modules, each consisting of 7 photo-multiplier tubes (PMTs), a Front-End Board and a local camera trigger system used for data acquisition. This paper addresses the timing performance of NectarCAM which are crucial to reduce the noise in shower images and improve image cleaning as well as to discriminate between gamma-ray photons and cosmic-ray background and finally to allow coincidence identification with neighboring telescopes for stereoscopic operations. Verification tests of the system have been performed in a dark room using various light sources to illuminate the first NectarCAM unit. The resulting timing precision and accuracy of the trigger arrival relative to a laser source, of individual and multiple pixel signals have been studied and are shown to comply to CTAO requirements.

Performance and systematic uncertainties of CTA-North in conditions of reduced atmospheric transmission

The Cherenkov Telescope Array (CTA) is the next-generation stereoscopic system of Imaging Atmospheric Cherenkov Telescopes (IACTs). In IACTs, the atmosphere is used as a calorimeter to measure the energy of extensive air showers induced by cosmic gamma rays, which brings along a series of constraints on the precision to which energy can be reconstructed. The presence of clouds during observations can severely affect Cherenkov light yield, contributing to the systematic uncertainty in energy scale calibration. To minimize these systematic uncertainties, a calibration of telescopes is of great importance. For this purpose, the influence of cloud transmission and altitude on CTA-N performance degradation was investigated using detailed Monte Carlo simulations for the case where no action is taken to correct for the effects of clouds. Variations of instrument response functions in the presence of clouds are presented. In the presence of clouds with low transmission (≤ 80%) the energy resolution is aggravated by 30% at energies below 1 TeV, and by 10% at higher energies. For higher transmissions, the energy resolution is worse by less than 10% in the whole energy range. The angular resolution varies up to 10% depending both on the transmission and altitude of the cloud. The sensitivity of the array is most severely reduced at lower energies, even by 60% at 40 GeV, depending on the clouds' properties. A simple semi-analytical model of sensitivity degradation has been introduced to summarize the influence of clouds on sensitivity and provide useful scaling relations.

Design and development of laser temperature control system of LHAASO

The Large High Altitude Air Shower Observatory (LHAASO) is located in Haizi Mountain, Daocheng County, Sichuan Province, with an average altitude of 4410 m. The Wide Field of view Cherenkov Telescope Array (WFCTA) is one of its three observation arrays. Laser calibration system for WFCTA consists of 2 nitrogen lasers and 1 Nd:YAG laser system. During the operation of WFCTA, the outdoor temperature may dropped from 17 °C to about -30 °C, which will bring serious challenges to the operation of laser calibration system operating in the field. In order to ensure the stable working environment and remote operation of the laser, we have developed the LHAASO-WFCTA laser temperature control, remote control and data monitoring system. The fluctuation of temperature in the box is within 1 °C around the target. The laser temperature control system was officially put into operation in Sept. 2020. In this paper, the design scheme, fabrication scheme, comparison of heat dissipation coefficient and test of temperature control system are introduced. The working data shows that both nitrogen molecule lasers and Nd:YAG laser temperature control system can achieve stable and long-term operation in high-altitude extreme field.

Mesotuned WIMPs and avoided deconfinement

We discuss the production of dark matter with mass of order $10$-$100$ TeV and $\mathcal{O}(1)$ couplings to the standard model, a scenario we refer to as the `mesotuned WIMP'. Given the lack of new physics observed at the LHC, indicating we live in a world with some degree of fine-tuning, this scale is a natural scale to expect new particles to appear. However, dark matter with mass in this range and large couplings is difficult to produce with the observed relic abundance via the usual mechanisms. Here we discuss the nonthermal production of dark matter and show that this mechanism can produce the observed dark matter in this mass range, but requires a specific relationship between the dark matter mass and the reheat temperature. It is shown that the avoided deconfinement model removes the dependence on the reheat temperature and predicts a dark matter candidate that can be tested at the Cherenkov Telescope Array and direct detection experiments, offering a natural setting for mesotuned WIMP dark matter.

Sensitivity of the Cherenkov Telescope Array to TeV photon emission from the Large Magellanic Cloud

ABSTRACT A deep survey of the Large Magellanic Cloud at ∼0.1–100 TeV photon energies with the Cherenkov Telescope Array is planned. We assess the detection prospects based on a model for the emission of the galaxy, comprising the four known TeV emitters, mock populations of sources, and interstellar emission on galactic scales. We also assess the detectability of 30 Doradus and SN 1987A, and the constraints that can be derived on the nature of dark matter. The survey will allow for fine spectral studies of N 157B, N 132D, LMC P3, and 30 Doradus C, and half a dozen other sources should be revealed, mainly pulsar-powered objects. The remnant from SN 1987A could be detected if it produces cosmic-ray nuclei with a flat power-law spectrum at high energies, or with a steeper index 2.3–2.4 pending a flux increase by a factor of >3–4 over ∼2015–2035. Large-scale interstellar emission remains mostly out of reach of the survey if its >10 GeV spectrum has a soft photon index ∼2.7, but degree-scale 0.1–10 TeV pion-decay emission could be detected if the cosmic-ray spectrum hardens above >100 GeV. The 30 Doradus star-forming region is detectable if acceleration efficiency is on the order of 1−10 per cent of the mechanical luminosity and diffusion is suppressed by two orders of magnitude within <100 pc. Finally, the survey could probe the canonical velocity-averaged cross-section for self-annihilation of weakly interacting massive particles for cuspy Navarro–Frenk–White profiles.

Probing gamma-ray burst afterglows with the Cherenkov Telescope Array

ABSTRACT Detection of delayed sub-TeV photons from gamma-ray bursts (GRBs) by MAGIC and HESS has proven the promising future of GRB afterglow studies with the Cherenkov Telescope Array (CTA), the next-generation gamma-ray observatory. With the unprecedented sensitivity of CTA, afterglow detection rates are expected to increase dramatically. In this paper, we explore the multidimensional afterglow parameter space to see the detectability of sub-TeV photons by CTA. We use a one-zone electron synchrotron and synchrotron self-Compton model to obtain the spectral energy distribution. We consider bursts going off in a medium of homogenous density. The blast wave is assumed to be radiatively inefficient and evolving adiabatically. Considering that the electron acceleration is not efficient if the acceleration time-scale exceeds the radiative cooling time-scale, we find that the sub-TeV emission is always due to the self-Compton process. We find that jets with high kinetic energy or large bulk Lorentz factor decelerating into a dense ambient medium offer better detection prospects for CTA. For relatively lower values of the downstream magnetic field, electrons are slow-cooling, and the emitted radiation is positively correlated with the magnetic field. For larger magnetic fields, the electron population enters the fast-cooling phase where the radiated flux is inversely proportional to the magnetic field. We apply our results in the context of bright TeV afterglows detected in recent years. Our results indicate that cosmological short GRBs have only moderate prospects of detection by CTA while local neutron star merger counterparts can be detected if the jet is launched towards the observer.

Probing gamma-ray bursts observed at very high energies through their afterglow

ABSTRACT A growing number of gamma-ray burst (GRB) afterglows is observed at very-high energies (VHE, ≳ 100 GeV). Yet, our understanding of the mechanism powering the VHE emission remains baffling. We make use of multiwavelength observations of the afterglow of GRB 180720B, GRB 190114C, and GRB 221009A to investigate whether the bursts exhibiting VHE emission share common features. We assume the standard afterglow model and microphysical parameters consistent with a synchrotron self-Compton (SSC) scenario for the VHE radiation. By requiring that the blastwave should be transparent to γ–γ pair production at the time of observation of the VHE photons and relying on typical prompt emission efficiencies and data in the radio, optical, and X-ray bands, we infer for those bursts that the initial energy of the blastwave is $\tilde{E}_{k, \rm {iso}} \gtrsim \mathcal {O}(10^{54})$ erg and the circumburst density is $n_0 \lesssim \mathcal {O}(10^{-1})$ cm−3 for a constant circumburst profile [or $A_\star \lesssim \mathcal {O}(10^{-1})$ cm−1 for a wind scenario]. Our findings thus suggest that these VHE bursts might be hosted in low-density environments, if the SSC radiation is responsible for the VHE emission. While these trends are based on a small number of bursts, the Cherenkov Telescope Array has the potential to provide crucial insight in this context by detecting a larger sample of VHE GRBs. In addition, due to the very poor statistics, the non-observation of high-energy neutrinos cannot constrain the properties of these bursts efficiently, unless additional VHE GRBs should be detected at distances closer than 15 Mpc when IceCube-Gen2 radio will be operational.

The diffuse gamma-ray flux from clusters of galaxies

The origin of the diffuse gamma-ray background (DGRB), the one that remains after subtracting all individual sources from observed gamma-ray sky, is unknown. The DGRB possibly encompasses contributions from different source populations such as star-forming galaxies, starburst galaxies, active galactic nuclei, gamma-ray bursts, or galaxy clusters. Here, we combine cosmological magnetohydrodynamical simulations of clusters of galaxies with the propagation of cosmic rays (CRs) using Monte Carlo simulations, in the redshift range z ≤ 5.0, and show that the integrated gamma-ray flux from clusters can contribute up to 100% of the DGRB flux observed by Fermi-LAT above 100 GeV, for CRs spectral indices α = 1.5 − 2.5 and energy cutoffs {E}_{max }=1{0}^{16}-1{0}^{17} eV. The flux is dominated by clusters with masses 1013 ≲ M/M⊙ ≲ 1015 and redshift z ≲ 0.3. Our results also predict the potential observation of high-energy gamma rays from clusters by experiments like the High Altitude Water Cherenkov (HAWC), the Large High Altitude Air Shower Observatory (LHAASO), and potentially the upcoming Cherenkov Telescope Array (CTA).

New determination of the production cross section for γ rays in the Galaxy

The flux of $\gamma$ rays is measured with unprecedented accuracy by the $\textit{Fermi}$ Large Area Telescope from 100 MeV to almost 1 TeV. In the future, the Cherenkov Telescope Array will have the capability to measure photons up to 100 TeV. To accurately interpret this data, precise predictions of the production processes, specifically the cross section for the production of photons from the interaction of cosmic-ray protons and helium with atoms of the ISM, are necessary. In this study, we determine new analytical functions describing the Lorentz-invariant cross section for $\gamma$-ray production in hadronic collisions. We utilize the limited total cross section data for $\pi^0$ production channels and supplement this information by drawing on our previous analyses of charged pion production to infer missing details. In this context, we highlight the need for new data on $\pi^0$ production. Our predictions include the cross sections for all production channels that contribute down to the 0.5% level of the final cross section, namely $\eta$, $K^+$, $K^-$, $K^0_S$, and $K^0_L$ mesons as well as $\Lambda$, $\Sigma$, and $\Xi$ baryons. We determine the total differential cross section $d\sigma(p+p\rightarrow \gamma+X)/dE_{\gamma}$ from 10 MeV to 100 TeV with an uncertainty of 10% below 10 GeV of $\gamma$-ray energies, increasing to 20% at the TeV energies. We provide numerical tables and a script for the community to access our energy-differential cross sections, which are provided for incident proton (nuclei) energies from 0.1 to $10^7$ GeV (GeV/n).

Deep learning techniques for imaging air Cherenkov telescopes

Very-high-energy (VHE) gamma rays and charged cosmic rays (CCRs) provide an observational window into the acceleration mechanisms of extreme astrophysical environments. One of the major challenges at imaging air Cherenkov telescopes (IACTs) designed to look for VHE gamma rays, is the separation of air showers initiated by CCRs which form a background to gamma-ray searches. Two other less well-studied problems at IACTs are (a) the classification of different primary nuclei among the CCR events, and (b) identification of anomalous events initiated by beyond-the-Standard-Model (BSM) particles that could give rise to shower signatures which differ from the standard images of either gamma rays or CCR showers. The problems of categorizing the primary particle that initiates a shower image, or the problem of tagging anomalous shower events in a model-independent way, are problems that are well suited to a machine learning approach. Traditional studies that have explored gamma-ray/CCR separation have used a multivariate analysis based on derived shower properties, which contains significantly reduced information about the shower. In our work, we address the problems outlined above by using machine learning architectures trained on full simulated shower images, as opposed to training on just a few derived shower properties. We illustrate the techniques of binary and multicategory classification using convolutional neural networks, and we also pioneer the use of autoencoders for anomaly detection at VHE gamma-ray experiments. The latter technique has been studied previously in the context of collider physics, to tag anomalous BSM candidates in a model-independent way. In this study, for the first time, we demonstrate the efficacy of these techniques in the domain of VHE gamma-ray experiments. As a case study, we apply our techniques to the High Energy Stereoscopic System experiment. However, the real strength of the techniques that we broach here in the context of VHE gamma-ray observatories, is that these methods can be applied broadly to any other IACT---such as the upcoming Cherenkov Telescope Array---or can even be suitably adapted to CCR experiments.