Cherenkov Telescope Array Observatory Research Articles

The population of Galactic supernova remnants in the TeV range

Supernova remnants (SNRs) are likely to be significant sources of cosmic rays up to the knee of the local cosmic-ray (CR) spectrum. They produce gamma rays in the very-high-energy (VHE) ($E>0.1$ TeV) range mainly via two mechanisms: hadronic interactions of accelerated protons with the interstellar medium and leptonic interactions of accelerated electrons with soft photons. Observations with current instruments have lead to the detection of about a dozen SNRs emitting VHE gamma rays and future instruments should significantly increase this number. Yet, the details of particle acceleration at SNRs and of the mechanisms producing VHE gamma-rays at SNRs remain poorly understood. We aim to study the population of SNRs detected in the TeV range and its properties and confront it to simulated samples in order to address fundamental questions concerning particle acceleration at SNR shocks. Such questions concern the spectrum of accelerated particles, the efficiency of particle acceleration, and the gamma-ray emission being dominated by hadronic or leptonic interactions. By means of Monte Carlo methods, we simulated the population of SNRs in the gamma--ray domain and confronted our simulations to the catalogue of sources from the High Energy Stereoscopic System (H.E.S.S.) Galactic Plane Survey (HGPS). We systematically explored the parameter space defined in our model, including for example, the slope of accelerated particles alpha , the electron-to-proton ratio $K_ ep $, and the efficiency of particle acceleration $ $. In particular, we found possible sets of parameters for which $ 90$<!PCT!> of Monte Carlo realisations are found to be in agreement with the HGPS. These parameters are typically found at $ 4.2 ep $, and $0.03 0.1 $ . We are able to strongly argue against some regions of the parameter space describing the population of Galactic SNRs in the TeV range, such as $ 4.35$, or ep Our model is so far able to explain the SNR population of the HGPS. Our approach, when confronted with the results of future systematic surveys, such as the Cherenkov Telescope Array Observatory, will help remove degeneracy from the solutions and to better understand particle acceleration at SNR shocks in the Galaxy.

Characterization and performance of an upgraded front-end-board for NectarCAM

This paper presents an analysis of the updated version of the Front-End Board (FEB) for the NectarCAM camera, developed for the Cherenkov Telescope Array Observatory (CTAO). The FEB is a critical component responsible for reading and converting signals from the camera’s photo-multiplier tubes into digital data and generating module-level trigger signals. This study provides an overview of the design and performance of the new FEB version, including the use of an improved NECTAr3 chip with advanced features. The NECTAr3 chip contains a switched capacitor array for sampling signals at 1 GHz and a 12-bit analog-to-digital converter (ADC) for digitization upon receiving a trigger signal. The integration of the new NECTAr3 chip results in a significant reduction of NectarCAM’s deadtime by an order of magnitude compared to the previous version. The paper also presents the results of laboratory testing, including measurements of timing performance, linearity, dynamic range, and deadtime, to characterize the new FEB’s performance.

The ASTRI Mini-Array: A New Pathfinder for Imaging Cherenkov Telescope Arrays

The ASTRI Mini-Array is an Istituto Nazionale di Astrofisica (INAF) project to build and operate an array of nine Imaging Atmospheric Cherenkov Telescopes (IACTs) at the Teide Astronomical Observatory of the Instituto de Astrofisica de Canarias in Tenerife (Spain) based on a host agreement with INAF and, as such, it will be the largest IACT array until the Cherenkov Telescope Array Observatory starts operations. Implementing the ASTRI Mini-Array poses several challenges from technical, logistic, and management points of view. Starting from the description of the innovative technologies adopted to build the telescopes, we will discuss the solutions adopted to overcome these challenges, making the ASTRI Mini-Array a great instrument to perform deep observations of the galactic and extra-galactic sky at very high energies.

Optical spectroscopy of blazars for the Cherenkov Telescope Array – III

Context. Blazars, which include BL Lacs and flat-spectrum radio quasars, represent the brightest persistent extragalactic sources in the high-energy (HE; 10 MeV–100 GeV) and very-high-energy (VHE; E > 100 GeV) γ-ray sky. Due to their almost featureless optical/UV spectra, it is challenging to measure the redshifts of BL Lacs. As a result, about 50% of γ-ray BL Lacs lack a firm measurement of this property, which is fundamental for population studies, indirect estimates of the extragalactic background light, and fundamental physics probes (e.g., searches for Lorentz-invariance violation or axion-like particles). Aims. This paper is the third in a series of papers aimed at determining the redshift of a sample of blazars selected as prime targets for future observations with the next generation, ground-based VHE γ-ray astronomy observatory, Cherenkov Telescope Array Observatory (CTAO). The accurate determination of the redshift of these objects is an important aid in source selection and planning of future CTAO observations. Methods. Promising targets were selected following a sample selection obtained with Monte Carlo simulations of CTAO observations. The selected targets were expected to be detectable with CTAO in observations of 30 h or less. We performed deep spectroscopic observations of 41 of these blazars using the Keck II, Lick, SALT, GTC, and ESO/VLT telescopes. We carefully searched for spectral lines in the spectra and whenever features of the host galaxy were detected, we attempted to model the properties of the host galaxy. The magnitudes of the targets at the time of the observations were also compared to their long-term light curves. Results. Spectra from 24 objects display spectral features or a high signal-to-noise ratio (S/N). From these, 12 spectroscopic redshifts were determined, ranging from 0.2223 to 0.7018. Furthermore, 1 tentative redshift (0.6622) and 2 redshift lower limits at z > 0.6185 and z > 0.6347 were obtained. The other 9 BL Lacs showed featureless spectra, despite the high S/N (≥100) observations. Our comparisons with long-term optical light curves tentatively suggest that redshift measurements are more straightforward during an optical low state of the active galactic nucleus. Overall, we have determined 37 redshifts and 6 spectroscopic lower limits as part of our programme thus far.

The Cherenkov Telescope Array Observatory workflow management system

The Cherenkov Telescope Array Observatory (CTAO) is the next generation ground-based observatory for gamma-ray astronomy at very high energies. It is expected to produce about 2 PB of raw data each year and to manage a global data volume which will grow through the years to reach about 100 PB in 2030. In addition, CTAO will require a high computing capacity for data processing and Monte Carlo simulations, of the order of hundreds of millions of CPU HS06 hours per year. To meet these requirements, CTAO will adopt a distributed computing model using 4 academic data centers, and will use the DIRAC framework as its workload management system. In the past ten years, to optimize the instrument design and study its performances, CTAO has used the European Grid Infrastructure (EGI) to run massive Monte Carlo campaigns. In order to handle these campaigns and to automatize simulation and data processing workflows, we have developed a production system prototype based on DIRAC. Recently, we have also developed a user interface allowing for the configuration and submission of complex workflows. In this contribution we present the production system prototype, its user interface for workflow management as well as its application to CTAO workflows.

Silicon photomultipliers for the SST camera of the Cherenkov Telescope Array

The Cherenkov Telescope Array Observatory (CTAO) will be the major global observatory for gamma-ray astronomy over the next decade and beyond. It will consist of two arrays of telescopes of different sizes, one for each hemisphere, and will be sensitive to gamma rays in the energy range from a few tens of GeV to hundreds of TeV. The Small-Sized Telescopes (SSTs) are a crucial component of the southern array, as they will extend the sensitivity of the observatory to the highest energies. Their focal plane will be equipped with 2048 Silicon Photomultiplier (SiPM) pixels, each one read independently by a state-of-the-art full waveform sampling readout. These solid-state sensors offer advantages over the traditional photomultiplier tubes, such as lower operating voltage, higher photon detection efficiency, and tolerance to bright illumination. In particular, they are the best choice for a small and compact camera such as the SST one. After a detailed comparative study, LVR3-type SiPMs from Hamamatsu Photonics were selected, with an active area of 6 mm × 6 mm, a microcell of 50 μm and without a protective coating, for optimum performance. The sensors demonstrated to have a higher photon detection efficiency and a lower crosstalk compared to the competitors, alongside a low dark count rate. In this contribution, we present the selection process and the latest measurements performed on the SiPMs mounted on the SST Camera module.

Performance of the joint LST-1 and MAGIC observations evaluated with Crab Nebula data

Aims. Large-Sized Telescope 1 (LST-1), the prototype for the Large-Sized Telescope at the upcoming Cherenkov Telescope Array Observatory, is concluding its commissioning phase at the Observatorio del Roque de los Muchachos on the island of La Palma. The proximity of LST-1 to the two MAGIC (Major Atmospheric Gamma Imaging Cherenkov) telescopes makes it possible to carry out observations of the same gamma-ray events with both systems. Methods. We describe the joint LST-1+MAGIC analysis pipeline and used simultaneous Crab Nebula observations and Monte Carlo simulations to assess the performance of the three-telescope system. The addition of the LST-1 telescope allows for the recovery of events in which one of the MAGIC images is too dim to survive analysis quality cuts. Results. Thanks to the resulting increase in the collection area and stronger background rejection, we found a significant improvement in sensitivity, allowing for the detection of 30% weaker fluxes in the energy range between 200 GeV and 3 TeV. The spectrum of the Crab Nebula, reconstructed in the energy range between ~60 GeV and ~10 TeV, is in agreement with previous measurements.

Multiwavelength observations of the lensed quasar PKS 1830-211 during the 2019 γ-ray flare

ABSTRACT PKS 1830 -211 is a γ-ray emitting, high-redshift (z =2.507 ± 0.002), lensed flat-spectrum radio quasar. During the period 2019 mid-February to mid-April, this source underwent a series of strong γ-ray flares that were detected by both AGILE-GRID (Gamma-Ray Imaging Detector) and Fermi Large Area Telescope (Fermi-LAT), reaching a maximum γ-ray flux of $F_{\rm E\gt 100\, MeV}\approx 2.3\times 10^{-5}$ photons cm−2 s−1. Here, we report on a coordinated campaign from both on-ground [Medicina, Owens Valley Radio Observatory (OVRO), Rapid Eye Mount (REM), and Sardinia Radio Telescope (SRT)] and orbiting facilities (AGILE, Fermi, INTEGRAL, NuSTAR, Swift, and Chandra), with the aim of investigating the multiwavelength properties of PKS 1830-211 through nearly simultaneous observations presented here for the first time. We find a possible break in the radio spectra in different epochs above 15 GHz, and a clear maximum of the 15 GHz data approximately 110 d after the γ-ray main activity periods. The spectral energy distribution shows a very pronounced Compton dominance (> 200) which challenges the canonical one-component emission model. Therefore, we propose that the cooled electrons of the first component are re-accelerated to a second component by, for example, kink or tearing instability during the γ-ray flaring periods. We also note that PKS 1830-211 could be a promising candidate for future observations with both Compton satellites [e.g. enhanced ASTROGAM (e-ASTROGAM)] and Cherenkov arrays [Cherenkov Telescope Array Observatory (CTAO)] which will help, thanks to their improved sensitivity, in extending the data availability in energy bands currently uncovered.

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.

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.

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.

Sensitivity of the Cherenkov Telescope Array to spectral signatures of hadronic PeVatrons with application to Galactic Supernova Remnants

The local Cosmic Ray (CR) energy spectrum exhibits a spectral softening at energies around 3 PeV. Sources which are capable of accelerating hadrons to such energies are called hadronic PeVatrons. However, hadronic PeVatrons have not yet been firmly identified within the Galaxy. Several source classes, including Galactic Supernova Remnants (SNRs), have been proposed as PeVatron candidates. The potential to search for hadronic PeVatrons with the Cherenkov Telescope Array (CTA) is assessed. The focus is on the usage of very high energy γ-ray spectral signatures for the identification of PeVatrons. Assuming that SNRs can accelerate CRs up to knee energies, the number of Galactic SNRs which can be identified as PeVatrons with CTA is estimated within a model for the evolution of SNRs. Additionally, the potential of a follow-up observation strategy under moonlight conditions for PeVatron searches is investigated. Statistical methods for the identification of PeVatrons are introduced, and realistic Monte-Carlo simulations of the response of the CTA observatory to the emission spectra from hadronic PeVatrons are performed. Based on simulations of a simplified model for the evolution for SNRs, the detection of a γ-ray signal from in average 9 Galactic PeVatron SNRs is expected to result from the scan of the Galactic plane with CTA after 10 h of exposure. CTA is also shown to have excellent potential to confirm these sources as PeVatrons in deep observations with O(100) hours of exposure per source.

Observation of the Cumbre Vieja volcano plume above the Observatorio del Roque de los Muchachos with the Barcelona Raman LIDAR

The Cherenkov Telescope Array Observatory (CTAO), currently under construction, is the next-generation very-high-energy gamma-ray observatory, providing the coverage for photons in the energy range 20GeV to 300TeV. CTAO will increase detection sensitivity in the 100 GeV to 10TeV range by a factor of 5 — 10 with respect to present experiments. CTAO retrieves the properties of very-high-energy gamma-rays by measuring Cherenkov light emitted by atmospheric showers of secondary particles that incident gamma rays produce in upper layers of the atmosphere. The key for reaching the required energy measurement accuracy is a precise knowledge of the atmospheric transmittance for Cherenkov light, which can be obtained using a dedicated Raman LIDAR. The device should operate at 355nm (near the maximum of Cherenkov light spectrum) and have the capability of taking data at specific azimuth and zenith angles up to distances of 30 km, so that atmospheric transmission along all possible air-shower directions can be determined. The Barcelona Raman LIDAR (BRL) is the official CTAO Pathfinder prototype, developed for atmospheric characterization of the Northern CTAO Site at the Observatorio del Roque de los Muchachos (ORM) on the Canary island of La Palma. BRL was deployed at ORM for extensive on-field tests between February 2021 and May 2022. We report on the commissioning results, including the remote operation capabilities of the system and its contribution to the understanding of atmospheric phenomena during its deployment period. In particular, we report on the properties of the volcanic plume from the eruption of the Cumbre Vieja volcano on 22 September 2021.

Assembly and performance of SiPM arrays for the prototype SCT proposed for CTA

The Near Ultraviolet High Density (NUV-HD) SiPMs produced by Fondazione Bruno Kessler have been employed to develop 16-pixel optical units to equip the focal plane of the prototype Schwarzschild–Couder Telescope (pSCT) proposed as a possible design for the Medium-Sized Telescope of the Cherenkov Telescope Array Observatory. After the assembly procedure, the optical units were tested and characterized to study their performance and homogeneity in terms of gain and dark count rate. In this work, we report on the assembly procedure and on the laboratory tests performed on different production of NUV-HD and the selection we made for the best quality sensors to be used in the installation on the telescope camera. Currently 36 NUV-HD3 optical units have been successfully integrated on the pSCT camera, together with 64 HAMAMATSU MPPCs. An upgrade of the pSCT camera is foreseen over the next years when the full focal plane is expected to be equipped entirely with FBK NUV-HD3 SiPMs, for a total of 11328 pixels.

Very-high-energy flat spectral radio quasar candidates

ABSTRACT The attenuation of very-high-energy (VHE) photons by the extragalactic background light (EBL) prevents the observation of high-redshift flat spectrum radio quasars (FSRQs). However, the correlation of the VHE spectral index with source redshift suggests that the EBL intensity may be less than what is predicted. This deviation can draw new constraints on the opacity of the Universe to VHE gamma-rays. Therefore, more FSRQs may fall above the sensitivity of the forthcoming VHE telescopes than the ones predicted by the existing EBL models. In order to account for the lower EBL intensity predicted by the index–redshift correlation, we introduce a redshift-dependent correction factor to the opacity, estimated from a commonly used cosmological EBL model. Considering this modified opacity, we identify the plausible VHE FSRQ candidates by linearly extrapolating the Fermi gamma-ray spectrum at 10 GeV to the VHE regime. Our study suggests that among 744 FSRQs reported in the Fermi Fourth Catalogue Data Release 2, 32 FSRQs will be detectable by the Cherenkov Telescope Array Observatory (CTAO). Because FSRQs are proven to be highly variable, we assume a scenario where the average Fermi gamma-ray flux increases by a factor of 10, and this predicts an additional 90 FSRQs that can be detected by the CTAO.

The Detection of GRBs at VHE: A Challenge Lasting for More than Two Decades, What Is Next?

Unveiling the mystery of gamma-ray bursts (GRBs) has been the target of many multi-waveband observational and theoretical efforts during the last decades. The results collected by current and past space-based instruments have provided important insights into the mechanisms at the origin of their prompt and afterglow phases. On the other hand, many questions, such as the the origin of the multi-GeV signal observed in a large number of events, remained unanswered. Within this framework, the first firm detections of a very-high-energy (VHE, E≳100 GeV) emission component by MAGIC and H.E.S.S. collaborations represented an important, long-awaited result for the VHE astrophysics community. However, while such discoveries opened a new era in the study of GRBs, they also provided an unexpected complexity due to the differences between the phenomenology of the observed events. This revealed that we still have an incomplete comprehension of GRB physics. In the nearby future, observations by the Cherenkov Telescope Array Observatory (CTAO), with unprecedented sensitivity in the VHE band, will have a key role in the study of these enigmatic objects and their interactions with the surrounding environment. In this review we will cover the recent GRB history, highlighting the efforts of follow-up campaigns by the VHE community that led to the first VHE GRB detection, and outlining what we can expect from future facilities in the next decades.

Gamma-Ray Bursts Afterglow Physics and the VHE Domain

Afterglow radiation in gamma-ray bursts (GRB), extending from the radio band to GeV energies, is produced as a result of the interaction between the relativistic jet and the ambient medium. Although in general the origin of the emission is robustly identified as synchrotron radiation from the shock-accelerated electrons, many aspects remain poorly constrained, such as the role of inverse Compton emission, the particle acceleration mechanism, the properties of the environment and of the GRB jet itself. The extension of the afterglow emission into the TeV band has been discussed and theorized for years, but has eluded for a long time the observations. Recently, the Cherenkov telescopes, MAGIC and H.E.S.S., have unequivocally proven that afterglow radiation is also produced above 100 GeV, up to at least a few TeV. The accessibility of the TeV spectral window will largely improve with the upcoming facility CTA (the Cherenkov Telescope Array). In this review article, we first revise the current model for afterglow emission in GRBs, its limitations and open issues. Then, we describe the recent detections of very high energy emission from GRBs and the origin of this radiation. Implications on the understanding of afterglow radiation and constraints on the physics of the involved processes will be deeply investigated, demonstrating how future observations, especially by the CTA Observatory, are expected to give a key contribution in improving our comprehension of such elusive sources.

Performance of a small array of Imaging Air Cherenkov Telescopes sited in Australia

AbstractAs TeV gamma-ray astronomy progresses into the era of the Cherenkov Telescope Array (CTA), there is a desire for the capacity to instantaneously follow up on transient phenomena and continuously monitor gamma-ray flux at energies above$10^{12}\,\mathrm{eV}$. To this end, a worldwide network of Imaging Air Cherenkov Telescopes (IACTs) is required to provide triggers for CTA observations and complementary continuous monitoring. An IACT array sited in Australia would contribute significant coverage of the Southern Hemisphere sky. Here, we investigate the suitability of a small IACT array and how different design factors influence its performance. Monte Carlo simulations were produced based on the Small-Sized Telescope (SST) and Medium-Sized Telescope (MST) designs from CTA. Angular resolution improved with larger baseline distances up to 277 m between telescopes, and energy thresholds were lower at 1 000 m altitude than at 0 m. The${\sim} 300\,\mathrm{GeV}$energy threshold of MSTs proved more suitable for observing transients than the${\sim}1.2\,\mathrm{TeV}$threshold of SSTs. An array of four MSTs at 1 000 m was estimated to give a5.7$\sigma$detection of an RS Ophiuchi-like nova eruption from a 4-h observation. We conclude that an array of four MST-class IACTs at an Australian site would ideally complement the capabilities of CTA.

Predictions of TeV emission for a set of hard BL Lac objects

ABSTRACT We focus our analysis on 55 BL Lac objects with a hard Fermi gamma-ray spectrum, and for which a redshift or a lower limit to it has been determined by a previous study of ours. We extrapolate the spectral fits given by the 4FGL catalogue to the VHE band (>0.1 TeV), which can be explored by imaging atmospheric Cherenkov telescopes. Furthermore, we take into account the absorption due to the extragalactic background light, strongly depending on the redshift. Finally, we compare our results with publicly available sensitivity curves for a selection of imaging atmospheric Cherenkov telescope arrays currently operating or under construction. From our extrapolations and simulations we find a large number of promising candidates for observation with the forthcoming Cherenkov Telescope Array observatory.

Sensitivity of the Cherenkov Telescope Array to dark subhalos

In this work, we study the potential of the Cherenkov Telescope Array (CTA) for the detection of Galactic dark matter (DM) subhalos. We focus on low-mass subhalos that do not host any baryonic content and therefore lack any multiwavelength counterpart. If the DM is made of weakly interacting massive particles (WIMPs), these dark subhalos may thus appear in the gamma-ray sky as unidentified sources. A detailed characterization of the instrumental response of CTA to dark subhalos is performed, for which we use the {\it ctools} analysis software and simulate CTA observations under different array configurations and pointing strategies, such as the scheduled extragalactic survey. This, together with information on the subhalo population as inferred from N-body cosmological simulations, allows us to predict the CTA detectability of dark subhalos, i.e., the expected number of subhalos in each of the considered observational scenarios. In the absence of detection, for each observation strategy we set competitive limits to the annihilation cross section as a function of the DM particle mass, that are at the level of $\langle\sigma v\rangle\sim4\times10^{-24}$ ($7\times10^{-25}$) $\mathrm{cm^3s^{-1}}$ for the $b\bar{b}$ ($\tau^+\tau^-$) annihilation channel in the best case scenario. Interestingly, we find the latter to be reached with no dedicated observations, as we obtain the best limits by just accumulating exposure time from all scheduled CTA programs and pointings over the first 10 years of operation. This way CTA will offer the most constraining limits from subhalo searches in the intermediate range between $\sim 1-3$ TeV, complementing previous results with \textit{Fermi}-LAT and HAWC at lower and higher energies, respectively.