Impact of cosmic-ray physics on dark matter indirect searches

The quest for the elusive dark matter (DM) that permeates the Universe (and in general the search for signatures of Physics beyond the Standard Model at astronomical scales) provides a unique opportunity and a tough challenge to the high-energy astrophysics community. In particular, the so-called DM indirect searches - mostly focused on a class of theoretically well-motivated DM candidates such as the weakly-interacting massive particles - are affected by a complex astrophysical background of cosmic radiation. The understanding and modeling of such background requires a deep comprehension of an intricate classical plasma physics problem, i.e. the interaction between high-energy charged particles, accelerated in peculiar astrophysical environments, and magneto-hydrodynamic turbulence in the interstellar medium of our Galaxy. In this review we highlight several aspects of this exciting interplay between the most recent claims of DM annihilation/decay signatures from the sky and the Galactic cosmic-ray research field. Our purpose is to further stimulate the debate about viable astrophysical explanations, discussing possible directions that would help breaking degeneracy patterns in the interpretation of current data. We eventually aim to emphasize how a deep knowledge on the physics of CR transport is therefore required to tackle the DM indirect search program at present and in the forthcoming years.

The quest for the elusive dark matter (DM) that permeates the Universe (and in general the search for signatures of Physics beyond the Standard Model at astronomical scales) provides a unique opportunity and a tough challenge to the high-energy astrophysics community.
In particular, the so-called DM indirect searches -mostly focused on a class of theoretically well-motivated DM candidates such as the weakly-interacting massive particles -are affected by a complex astrophysical background of cosmic radiation. The understanding and modeling of such background requires a deep comprehension of an intricate classical plasma physics problem, i.e. the interaction between high-energy charged particles, accelerated in peculiar astrophysical environments, and magneto-hydrodynamic turbulence in the interstellar medium of our Galaxy.
In this review we highlight several aspects of this exciting interplay between the most recent claims of DM annihilation/decay signatures from the sky and the Galactic cosmic-ray research field. Our purpose is to further stimulate the debate about viable astrophysical explanations, discussing possible directions that would help breaking degeneracy patterns in the interpretation of current data. We eventually aim to emphasize how a deep knowledge on the physics of CR transport is therefore required to tackle the DM indirect search program at present and in the forthcoming years.

EXOTICA: WHERE TO FIND THEM?
The particle dark matter (DM) [1] discovery may potentially undertake a different path than the beaten track of collider searches [2] and direct detection experiments [3]. Early Universe thermal relics, in particular, may be well-motivated DM candidates [4][5][6] expected to annihilate (or decay) even in today galactic halos, producing Standard Model (SM) particle yields. Therefore, the measurement of charged particle and gamma-ray fluxes of cosmic origin in a wide energy range -say from few MeV all the way up to the multi-TeV domain -may be a quite unique tool at our disposal in order to probe the putative particle nature of DM [7,8].
It is widely recognized that the indirect extraction of a DM signal in this context is an extremely challenging task. The DM problem in general and the indirect searches in particular have been already presented in comprehensive review papers (see e.g. [9,10]). Here, we wish to focus our attention on the important interplay between particle DM signatures and the background signals expected from astrophysics, discussing in particular the phenomenological relevance of cosmic-ray physics.
As far as charged particles are concerned, the first consideration in order is that high-energy protons, nuclei and electrons are injected in copious amounts as cosmic rays (CRs) by different classes of astrophysical sources (such as shocks associated to supernova explosions or superbubbles, or possibly accretion-powered mildly relativistic jets), which provide a huge and irreducible background. On the other hand, the paucity of anti-particles, mostly produced by secondary interactions of CRs, could in principle drastically improve the signal-to-noise ratio in favor of a putative DM detection, within a rather low expected astrophysical background; this possibility of detecting early Universe relics by studying Galactic antiparticles was first outlined in the early 1980s in several pioneering papers (e.g. [11,12]) and has been studied in much larger detail in particular during the last decade, mainly thanks to the dramatic improvement in the quality of the data provided by PAMELA [13] and AMS-02 [14] experiments. For a recent discussion on cosmic antimatter opportunities, see for example [15,16].
Despite the low background, even in the case of antimatter searches the large uncertainties involved in the modeling of both conventional astrophysical production and Galactic transport play a major role and have hindered a firm DM detection so far, although some recent tentative claims (in particular, among others, [17][18][19][20][21]) have triggered an important debate in the community. In sections 4,5,6 we will describe in detail the current status for positron, anti-proton and antinuclei indirect searches: Our purpose is to provide a case-by-case discussion mostly focused on the relevance of CR transport physics in antimatter channels for DM indirect detection.
While charged antiparticles may be promising indirect messengers of the particle DM nature, they do not retain the directionality from their emission point and, hence, cannot provide the morphological characterization of a DM signal: This possibility is accomplished instead with the analysis of gamma-ray data (for early studies, see e.g. [22][23][24]; more recently, [25,26]). If DM particles eventually decay or annihilate into gamma radiation, it is crucial to identify the most promising regions of the sky where either the expected signal is large, or the astrophysical background is low: Among the most important targets, we can certainly mention the inner Galaxy (which satisfies the first criterion), and the DM-dominated satellite galaxies orbiting around our Galaxy. The most tentative claims and interesting bounds from the study of the gamma-ray sky will be reviewed in detail in section 7.
Before going through an extensive discussion of all these channels for DM indirect searches, in the next two sections we will set up the stages of this review by briefly presenting some relevant aspects of DM models connected with the scope of the present paper, and then highlighting the several key aspects of the Galactic CR transport problem. Those concepts will be recalled all across the paper when the most relevant interpretations of CR and gamma-ray anomalies are discussed.

TARGETING DM INDIRECT SEARCHES ON HIGH ENERGY PHYSICS
The quest for the fundamental origin of DM may require to consider a priori an impressive range of energy scales (see [27,28] for interesting historical retrospectives). For instance, sitting on the extremes of the viable mass window for DM searches, today we may be looking for imprints on the cosmological matter power spectrum of super-light candidates such as axion-like particles [29][30][31] from the string landscape [32,33], as well as aiming to detect the gravitational echoes of massive black hole merging [34,35], possibly originated from primordial density fluctuations in the early Universe [36,37].
However, the phenomenology of DM candidates may be intriguingly correlated to the long-standing puzzles pertaining to the realm of the SM (see, e.g., [10,38] for broad reviews on the topic). Of particular significance, one of the main driving forces of research on High Energy Physics has been the quadratic UV sensitivity of the Higgs boson mass to any New Physics energy scale above the electroweak one [39]. On general grounds, the attention of this review is mostly paid to DM candidates motivated by New Physics at the electroweak scale. Many extensions of the SM theory addressing the electroweak hierarchy problem can indeed accommodate such DM candidate in their spectrum; see, e.g., [5] for a review on the widely mentioned case of supersymmetry; [40,41] for higher dimensional theories aimed at stabilizing the electroweak scale; [42][43][44][45] for examples of strongly coupled UV completions. Moreover, DM candidates sitting around the electroweak scale may also be well-motivated in the context of non-natural theories addressing other possible SM issues, such as proton stability and/or gauge coupling unification; see e.g. [46][47][48] for the case of supersymmetry, [49][50][51][52] for interesting universal/warped extra-dimensional proposals.
In all these scenarios, the DM particle is usually stable due to its charge under a (discrete) symmetry of the new theory, while a benchmark range of masses and couplings can be eventually individuated on the basis of the UV guiding principles. Most importantly, the emerging DM phenomenology from these studies typically falls in the experimental window of sensitivities for antimatter and gamma-ray searches discussed in the next sections.
Marginalizing over the specific details of UV models, the DM reference framework we mainly refer to, in this review, corresponds to early Universe cold thermal relics. Assumed to be in thermal equilibrium with the primordial plasma at the very early stages, DM decoupling as non-relativistic species eventually leads to [53] (for more details, see also [54]): Ω is the DM cosmological relative abundance observed today [55], σv is the DM thermally-averaged particle annihilation cross section and in the last step we make use of the approximation σv σ c/3. By means of dimensional analysis, we can naively estimate the DM annihilation cross section in terms of its mass m and dimensionless coupling constant g, obtaining: Eq. (1) in conjunction with the estimate in Eq. (2) characterizes the so-called Weakly Interacting Massive Particle (WIMP) miracle: A cold thermal relic charged under weak interactions, and with mass close by the electroweak scale, naturally accounts for the present DM abundance. The WIMP mass range may be bracketed from below according to the seminal works in [56,57] as m 10 GeV, while from the above unitarity arguments on the DM cross section pinpoint to m 10 2 TeV [58] (but caveats exist [59]). Typical expectations from the WIMP paradigm may be notably disregarded by the presence of resonant regimes, near mass thresholds, co-annihilation with other particles in the thermal bath [60]. Many other examples of variants to the standard WIMP freezeout scenario have been investigated in literature [61][62][63][64][65][66][67][68][69]. Non-perturbative effects such as Sommerfeld enhancement [70][71][72][73] and bound-state formation [74][75][76][77][78] have also more recently acknowledged to be of possible dramatic impact in the broad context of WIMP phenomenology. Importantly, Eq. (1) shows that the DM relic abundance is mostly sensitive to the annihilation cross section: It follows that a WIMP-less miracle is perfectly conceivable [79], opening on general grounds a broader range of viable mass scales and couplings for the phenomenology of DM thermal relics. Therefore, indirect signals from WIMP-like scenarios -whose signatures have been comprehensively inspected in [80,81] and are of particular importance for this review -are after all: 1) intimately connected to the possibility of being effectively visible today (in particular, non exhibiting important velocity suppression); 2) possibly connected to the tantalizing discovery of New Physics near the electroweak scale [82].

Origin of CRs
A clear identification of the classes of astrophysical sources able to accelerate particles from GeV all the way up to PeV energies (and, for extragalactic accelerators, up to ∼ 10 20 eV) is crucial for DM indirect detection.
We will not digress here on the long debate about the origin of CRs and their acceleration mechanisms (see e.g. [83] for an excellent review). For the purpose of this review, let us start by mentioning the supernova paradigm as the main guideline. Supernova remnants (SNRs) were proposed as potential sources out of energy budget arguments in [84,85]; the picture was better defined later in terms of SNRs located in our own Galaxy [86,87]; however, a physical process capable of such a powerful CR acceleration had not been proposed yet at that time. Later, the theory of diffusive shock acceleration was presented in four famous seminal papers [88][89][90][91], and is currently considered the main reference framework in the field.
However, it is important to keep in mind that other classes of sources (e.g. pulsar wind nebulae for leptonic CRs [92], OB associations [93], X-ray binaries [94]), and several other acceleration mechanisms, have been proposed as well. As we will see below, many of these potential accelerators are expected to play a role in the explanation of some tentative claims of DM detection.

CR transport: Preliminary considerations
Let us now turn our attention to a crucial aspect of CR physics, which has an extremely relevant impact on DM indirect searches, i.e. CR propagation in our Galaxy.
The usual starting point is a collection of several key observations that characterize the cosmic-ray flux.
• The isotropy of the arrival direction (at the level of ∼ 10 −3 in the TeV -PeV range, recently measured with high precision by many experiments [95][96][97][98][99][100]); • The much larger abundance of Lithium, Beryllium and Boron compared to the solar system abundances, which is naturally interpreted as the signature of the interaction of primary species such as protons and heavy nuclei with a column density of interstellar gas as large as few g/cm 2 : Such a grammage implies that the primary species have crossed the Galactic disk many times; • The presence of a diffuse gamma-ray emission across the whole Galactic disk, already predicted in the early 1960s [101], and first measured by pioneering satellite experiments such as OSO-3 (1967) and SAS-2 (1972); nowadays, Fermi-LAT has provided a state-of-the art description of this emission in the 300 MeV -300 GeV range, as detailed below.
These pieces of information, combined together, suggest a "conventional scenario" for CR transport that was shaped by the pioneering work of Ginzburg and colleagues (see [102] and references therein, and [103]), based on a random walk through the Galaxy governed by the quasi-linear theory of pitch-angle scattering on Alfvénic turbulence (QLT), first presented in the 1960s [104,105].

The CR transport equation
Magneto-hydrodynamic turbulence -which is ubiquitous in the interstellar medium (ISM) and covers a very wide range of scales from astronomical units (AUs) to kpc [106] -is widely considered as the main responsible for this diffusive regime. In more detail, the relativistic motion of charged particles in our Galactic environment is affected by the presence of both a coherent large-scale magnetic field component, B [107,108], on top of which magnetic inhomogeneities, δ B are propagating. These fluctuations in the magnetized interstellar medium (ISM) are associated to a turbulent cascade that is believed to be either initiated at large scales ∼ 10 2 pc (by supernova explosions, differential rotation of the Galactic disk, or other mechanisms [109]), or (especially at small scales) triggered by CR themselves via streaming instability [110,111]. This cascade has been usually considered, in the basic scenario of QLT, as isotropic and mainly composed of Alfvén waves, i.e. transverse magnetic perturbations moving at the Alfvén speed 1 : Following in part the approach of [113], we recap here the main features of QLT (see also [102,103]).
The rationale of QLT is to consider the interaction of a charged particle of momentum p = m v with magnetic inhomogeneities δ B that are sufficiently small (with respect to the regular field B) at the scale of interest. The process is well described by a stochastic equation for the pitch angle, defined as µ = cos(p ∧ B). On average the variance of the pitch angle can be shown to feature a resonance condition [103]. According to it, the particle only interacts with the inhomogeneities of wavelength ∼ 2π/k matching the particle Larmor gyroradius r L 2 : Let us now consider an ensemble of particles described by a phase-space distribution f ( x, p, t), with probability density Ψ( p, ∆ p ) for transitions p → p + ∆ p in momentum space, due to interactions with stochastic fluctuations in the magnetized environment. We can state that, after a lapse ∆t, in the Alfvén wave rest frame (primed), the evolved phase-space distribution must be equal to: being v the CR particle velocity in the wave frame. We assume that detailed balance holds, i.e. a transition p → p − ∆ p described by a probability Ψ( p , −∆ p ), is equivalent to the one described by Ψ( p − ∆ p , ∆ p ). Applying this principle in the limit of ∆p /p 1 (originating from |δ B |/| B | 1) that characterizes QLT, we can write: where . . . ∆ p ≡ d∆ p Ψ( p , ∆ p ), and a Taylor expansion of Ψ( p , ∆ p ) has been performed.
Assuming the static limit, i.e. ∆t/t 1, starting from a Taylor expansion of Eq. (5), with the help of Eq. (6) we finally get, in the same wave frame: This is a Boltzmann equation where in the right-hand side the "collision operator" qualifies Brownian motion in momentum space with diffusion coefficient describing indeed the momentum isotropization due to CR stochastic scattering with Alfvén waves. 3 Let us now perform a transformation to the Galactic rest frame. In this frame of reference, Eq. (7) features a spatial diffusion operator along the direction of the regular field as well: This is the most important term that governs CR transport in the Galaxy. The spatial diffusion coefficient D zz is related to D pp by [114]: for B = Bẑ . Inspired and guided by the results of QLT, a general transport equation is usually considered, mainly based on the aforementioned process of diffusion in both position and momentum space, but featuring a wider set of terms associated to other physical phenomena. The full equation reads: with N ( x, p, t) and Q( x, p, t) being respectively the CR density species and CR injecting density source per unit of momentum. In the left-hand side, the diffusion term is usually isotropic and described by a scalar, positionindependent coefficient, despite the fact that QLT predicts a highly anisotropic transport along the regular field direction (see e.g. the discussion in [115]). The scalar spatial diffusion coefficient is generally taken as where F(k) is defined as the (normalized) power associated to the turbulent modes with wave number k ∝ 1/p resonating with the particles carrying momentum p. Since the turbulent power scales as a power law, the rigidity dependence of the diffusion coefficient is usually parametrized as: with D 0 and δ as free parameters to be fixed by comparison with CR data. The spatial dependence of such normalization, that stems from the spatial variations of the turbulent power, is usually neglected with some relevant exceptions [115][116][117][118][119][120].
The momentum diffusion is also called stochastic reacceleration, and the relation 10 is assumed to hold.
The left-hand side also involves an advection term originally present in Eq. (7) as well, now characterized by the bulk velocity of the plasma in the lab frame: Galactic winds affecting CR motion may be described by such a term, together with adiabatic energy losses, involving velocity gradients even of O(10 2 ) km s −1 kpc −1 perpendicularly to the Galactic disk.
The physics of advective-diffusive transport is enriched by two more relevant phenomena included in Eq. (11): net energy losses and spallation. In fact, we need to consider CR energy loss processes, characterized by the continuous loss rate dp/dt, particularly important e.g. at high energies for light charged species such as leptons (see e.g. [121] and [122]).
Eventually, on the right-hand side of Eq. (11): • The primary source term Q 0 captures the primary accelerators of CRs: As mentioned above, while the supernova paradigm is still the most accredited one, other classes of sources can certainly be at work; • The secondary source term Q sec describes the production of a given species from spallation of the heavier ones onto interstellar gas; • A loss term due to inelastic collisions characterized by an interaction time τ N is also introduced.
An important remark is needed at this point. Our picture of MHD turbulence has dramatically improved during the latest decades: According to the current scenarios [123,124], MHD turbulence is composed of an anisotropic cascade of both Alfvén waves, and isotropic fast magnetosonic modes, as theoretically demonstrated and numerically confirmed by several simulations. As a consequence of the anisotropy of the Alfvénic cascade, the scattering efficiency on Alfvén waves turns out to be very low [125]. posing a tough challenge to the whole scenario discussed above. Among others, a possible solution [126,127] is that magnetosonic modes dominate gyroresonance interaction for most of the pitch-angle range.
However, although the actual microphysics underlying the CR random walk is still far from being exhaustively addressed, the QLT can still be considered a useful guideline to be taken as a reference, and Eq. (11) should be understood as a phenomenological tool to tame the complexity of the plasma physics problem, allowing us to make predictions against a plethora of data.
While a detailed and realistic study of Galactic CR propagation requires the extensive use of those numerical or complex semi-analytical methods, we can extract some physical insight useful for the next sections looking at a simplified version of Eq. (11).
Indeed, at the basis of CR transport may be conceived an important hierarchic game of scales: convective and re-acceleration effects are typically related to low-energy regimes, while energy-loss rates are negligible for highenergy hadrons such as protons or heavier nuclei. So, for energies E few GeV, the Galactic motion of heavy species can be approximately described to be in a purely diffusive regime: Then, we may trade the spatial diffusion operator for an effective time of confinement, τ D , i.e.
and treat the Galaxy as a box where CRs perform a random walk up to the box boundaries, beyond which they leak out. In the steady-state limit of this simplified framework, CR secondary species, produced by the interaction of source-injected (namely, primary) CR particles, can be estimated as: Hence, the timescale for confinement of CRs in the Galaxy is intimately linked to relative abundances of secondary and primary species. As anticipated in the above, anti-matter species are typically produced as secondaries, and therefore Eq. (15) is of direct relevance for DM indirect searches in anti-matter channels.
Note that by dimensional analysis, τ D ∼ H 2 /D, where H captures the typical size of the box. Therefore, secondary fluxes are sensitive to both the spatial diffusion coefficient and the height of the CR propagation halo. As we will discuss below, our current poor knowledge of the value of H generally translates into an important source of uncertainty in DM indirect detection studies.
As stated above, the spatial diffusion coefficient scales as a power law with rigidity. Supplementing Eq.(15) with such expectation, we have a theoretical prediction that nicely fits the trend of available experimental data [138]: at energies above few GeV, the measurements of local CR secondary-over-primary observables like B/C are indeed compatible with a power-law behavior (with index δ) of the diffusion coefficient. Typical estimates of this parameter are in the range 0.3 < δ < 0.6 [139][140][141][142][143], with a normalization at GeV corresponding to D 0 10 28 cm 2 s −1 , but pertain only to a local measurement (see [117] for a possible indirect inference of δ across the Galaxy).

THE ANTIPROTON CHANNEL
In the conventional scenario sketched above, antiprotons are produced in the Galactic environment by spallation of heavy nuclei and protons onto interstellar gas.
The early measurements of the antiproton flux date back to the 1970s and early 1980s; a first tentative claim of anomaly with respect to the conventional expectations based on the picture of nuclear spallation goes back to [145]. However, DM connections were not outlined at that time: Cosmic antiprotons were considered a promising channel for DM searches only some years later, in several seminal papers (see, for instance, [11,146,147] and, more recently, [148]).
A dramatic improvement in the accuracy of the data was provided by the PAMELA collaboration in 2009 [149] (with further refinement in [150]): The measurement showed a reasonably good agreement with conventional models based on purely secondary origin, as confirmed by [151] and [152]. A note here is in order about the meaning of conventional model: In what follows, we will use this expression for a model based on the simplest version of Eq. (11), taking a single class of sources (SNRs) at work, with antiparticles produced only as secondary products from primary spallation, and featuring constant and homogeneous diffusion, tuned on local CR data.
Given the absence of significant unexpected spectral features (such as bumps), the potential constraining power of PAMELA dataset for DM searches was soon demonstrated in a series of papers (from the early ones as [153,154] to the more recent [155]), which provided a comprehensive discussion on the upper bound on the WIMP annihilation cross section from the detection of cosmic antiprotons. The most relevant point made in those papers is the crucial role of CR transport. A major source of uncertainty, in particular, is the size of the diffusive halo H introduced in the previous section, i.e. the volume where Galactic CRs are effectively confined by the presence of a turbulent magnetic field. Models of CR transport based on larger diffusion halos usually feature larger average values of the diffusion coefficient in order to correctly reproduced the secondary/primary ratio data: Therefore, in these scenarios, the antiproton flux probe a larger region of the Galaxy, and exhibits a larger constraining power. On the other hand, cases where one assumes a very thin halo (smaller than 2 kpc) turn out to be much less restricting on the particle DM properties indirectly probed.
Another source of uncertainty certainly lies on the properties of CR transport in a much smaller environment: the Heliosphere. We refer to [155] and references therein for a comprehensive discussion on this aspect.
As a consequence of this complicated puzzle, it was not possible to firmly exclude some relevant tentative DM claims made in other channels (see the gamma-ray section for more details), and the most severe limitation came from the poor constraints we actually have on the size of the Galactic CR diffusion halo. A possible improvement in this direction may come from more accurate measurements of the Beryllium isotopes, and from a careful analysis of the current and forthcoming data on the vertical profile (with respect to the Galactic plane) of the synchrotron emission from the Galaxy (following [156,157], we point out that the current data seem to favor large values for the size of the diffusion halo).
In 2015 a much more accurate dataset was published by the AMS-02 collaboration [19]. The debate on the antiproton channel has then included tentative claims of anomalies with respect to the conventional scenario, possibly explained in terms of DM indirect detection.
First of all, the AMS-02 collaboration itself initially claimed the presence of an excess at high energies over 100 GeV. Right after, the significance of this anomaly was better characterized [144,158], pointing only to a mild overshooting of the expected background. The relevance of this putative discrepancy an the estimated impact of the different sources of uncertainty on the model predictions is well depicted in Fig. 1 (taken from [144]). 4 Taking at face-value the original tentative claim from AMS, the interpretation of the excess requires scenarios beyond the conventional on of CR transport pictured in the mini-guide of section 3. For instance, a mechanism that may be at work and explain the discrepancy is the secondary production at the accelerator: The idea, proposed before AMS data in [166] as a possible explanation to the positron ratio anomaly (see section 6), is that secondary products of hadronic interactions inside the sources can participate in the acceleration process and subsequently escape into the interstellar medium as an extra component featuring a very flat spectrum. In [167] this scenario was applied to (pre-AMS) antiproton data as well, and the authors demonstrated that the boronover-carbon ratio has much constraining power for this interpretation.
More recently, refined Bayesian analyses that include the official AMS data on this observable as well (see e.g. [142]) confirm the presence of a mild discrepancy between the regions of the parameter space pointed by AMS antiproton and B/C data.
However, spatial-dependent diffusion setups (e.g. the phenomenological two-zone models as those considered in [168] and [169], designed to capture both CR transport in pre-existing SNR-driven Galactic turbulence, and confinement by CR-driven turbulence via streaming instability) seem to solve the discrepancy, as well as the latest scenarios that include secondary production at the accelerator [170].
DM interpretations for the high-energy discrepancy are also still viable for quite large values of the DM particle mass (as shown, e.g., in [171][172][173]), in particular for light mediator scenarios. At the same time, the constraining power of AMS data in the energy range where no relevant feature or anomaly is present has been most recently explored in [174,175].
On the other hand, on the low-energy side, a possible indication of a DM signal for DM masses near 80 GeV has been found [175][176][177], with a hadronic annihilation cross section close to the thermal value: Interestingly, this tentative claim is compatible with the DM interpretation of the Galactic center gamma-ray excess (see section 7). Again, more investigation of the transport uncertainties (both in standard and beyond-standard scenarios), and more detailed combined studies of this signal together with constraints from other probes (e.g., the observation of dwarf spheroidals in the gamma-ray band) will be crucial in order to confirm the existence of this anomaly.

THE AVENUE FOR ANTINUCLEI
A milestone campaign for imprints of particle DM on the observable CR radiation may correspond to the discovery of Galactic light antinuclei such as antideuteron ( 2 H) [178][179][180][181] and antihelium-3 ( 3 He) [182,183]. At present, no compelling evidence for a detection of antimatter with mass number A ≥ 2 has been experimentally corroborated in the measurement of Galactic CR fluxes. A notable upper limit on the flux of cosmic antideuteron has been set by the BESS facility -reporting at the 95% of confidence level Φ 2 H 2 × 10 −4 (m 2 s sr GeV/n) −1 for kinetic energy per nucleon 0.17 ≤ T ≤ 1.15 GeV/n [184] -while the BESS-Polar collaboration currently constrains antihelium-to-helium flux to be smaller than 10 −7 in the interval of probed rigidities, 1.6 -14 GV [185].
Today, the AMS-02 mission is operating in the direction to perform the first historical observation of cosmic light antinuclei, with promising projected sensitivities [186]. A tentative claim of few antihelium eventspossibly measured by the AMS-02 collaboration -might be already at hand, waiting for a firmer experimental response in the upcoming years; see, for instance, ref. [187].
The phenomenological relevance of the avenue for antinuclei in relation to DM indirect searches stems from the kinematics of spallation processes producing CR secondaries: the energy threshold associated to the production of one antideuteron is roughly 2.5 greater than the one required to produce a secondary antiproton. Moreover, the same energy threshold is a monotonic increasing function of the mass number A. Therefore, the low energy flux of cosmic antinuclei is expected to be small [188], opening a low-energy window -related to kinetic energies per nucleon between ∼0.1 and few GeV/n -for large signal-to-background ratios from exotica.
The physics of CR accelerators and of propagation impacting the prediction for antinuclei fluxes at Earth sits on the very similar grounds of the one discussed for the antiproton channel. On the one hand, the unknowns stemming from Galactic CR propagation affecting the predicted Galactic fluxes for antideuteron and antihelium should be correlated to the Galactic antiproton spectrum and boron-to-carbon ratio data. On the other one, the injection of antinuclei in the Galactic interstellar medium is intimately connected to the antiproton production mechanism. As discussed in the previous section, the latter may be sourced from in situ acceleration in the downstream region of supernova remnant shockwaves. Such mechanism generally yields a source term distinguished into a A-component, referred to CR secondaries accelerated by the shock, and a B-component, related to standard spallation processes sourced by CR primaries. Most recent analyses on the topic have included also these extra contributions in their predictions [189][190][191]. Fig. 2, from [189], reports the maximum contribution to the antideuteron flux with respect to the typical unknowns from CR physics, properly calibrated on antiproton AMS-02 data. The components from production in supernovae (SNR-A and SNR-B bands) are found to be subdominant, giving (at most) a 10% effect with respect to the main standard component, obtained by considering the interactions of CR primaries with the interstel- Maximum flux for the expected astrophysical background is also shown: in black color the contribution from CR spallation processes with the interstellar medium, and the subleading SNR-A and SNR-B ones. Bands bracket uncertainties from force-field approximated solar modulation and an event-byevent coalescence model. Credit to [189], Figure 3. lar medium (black line). Fig. 2 also shows predictions for some benchmark exotica such as annihilating DM particles with masses of O(100) GeV and 100% branching ratios into bb or W + W − final states, together with the expectations from exotica that gained recent interest [34,35] such as primordial black holes [192,193]. The plot clearly highlights the importance in the aforementioned low-energy window. However, contrary to common wisdom [186], Fig. 2 also underlines the unlucky possibility that within the forthcoming years current and future experimental facilities may not be sensitive to antideuteron yields from DM/exotica production.
At this point, it is important to stress that unknowns stemming from the physics of CR accelarators and Galactic transport should not be retained to be the major source of uncertainty in the prediction of cosmic antinuclei fluxes. The bands reported in Fig. 2, while including the effects of solar modulation -typically studied within the force field approximation [194] and better investigated in [195] by means of numerical tools [196,197] -are most importantly related to the coalescence model adopted to establish antideuteron formation. At present, the prediction of antinuclei fluxes seems to strongly depend on the assumptions made to describe antiproton-antineutron fusion [186]. Nowadays, state-ofthe-art analyses can avoid to rely on simplistic analytical modeling [198], making instead use of Monte Carlo event generators and available data from colliders. Large systematics on the formation of antideuteron and antihelium can be understood on the basis of the sensitivity of these studies on the hadronization model implemented and the experimental dataset considered, see e.g. [186,199]. In-terestingly, a recent work focused on the description of nuclear coalescence via a physical modeling for the fusion of nucleons into composite nuclei [200], and exploiting two-particle correlation measurements [201], has pointed out the possibility that the production cross section in pp collisions for antihelium-3 may have been underestimated by up to two orders of magnitude. Consequently, a putative detection of cosmic 3 He events related to kinetic energies greater than 1 GeV/n may actually be within the reach of AMS-02 in few years [16]. f Therefore, we may conclude that searches for antinuclei, while being potentially exposed to CR propagation details, are mainly plagued by the assumptions and the systematics involved in the estimate of the poorly known production cross section. This can have a dramatic impact on our present ability to make projections for a concrete signal detection of antideuterons and antihelium-3. At the same time, these uncertainties leave us hope to foresee a spectacular discovery of exotica such as DM in the peculiar window of low-energy antinuclei events.

THE POSITRON CHANNEL
The positron channel has been under the spotlight for a long time in the DM indirect detection community. As mentioned in section 1, the paucity of those particles makes them an ideal target for DM searches, and the presence of a significant anomaly (with respect to the conventional expectations) has further increased the interest around this observable during the latest decade.
Let us start by clarifying, also in this case, what we mean by "conventional" predictions. In the context of a simplified treatment of acceleration and transportas described in the mini-guide of CR physics -a slightly different discussion is actually required for leptons in general. In fact, high-energy leptons feature relevant energy losses especially at high energies, which implies a new timescale competing with the diffusion one. The simplified transport equation governing e ± propagation is therefore of this kind: where inverse Compton scattering on the photon background and synchrotron radiation define the typical timescale for energy losses, τ E loss (p) ∼ p/|dp/dt| ∝ p −1 . The solution of Eq. (16) is given by a Green function that boils down to: For primary electrons this result implies a scaling N ∝ p αe− δ  FIG. 3. The most recent leptonic data from AMS and DAMPE, interpreted within a pulsar scenario. Credit to [202], Figure 2.
In the standard scenarios, positrons in the Galaxy are believed to originate, like antiprotons, as an entirely secondary component arising from the collisions of relativistic protons with the ISM gas according to a chain like p + H → · · · → π ± → µ ± + · · · → e ± + . . . . The source function of positrons is then expected to scale as: where we have implemented the rather simplistic approximation of an energy independent cross section σ p→e + , and used the fact that the relevant timescale for propagation of high-energy protons is the diffusion timescale, and defined α p the spectral index for the proton injection source function. Plugging this result in Eq. (18), we can predict the scaling of the propagated positron flux over the electron one: with α e the spectral index at injection for the electron source distribution.
In the framework of diffusive shock acceleration, the injected spectral index should not differ much among different species. Consequently, the ratio of secondary positrons over primary electrons is predicted to decrease with increasing energy, unless a (very unlikely) large difference between the source spectral indexes for protons and electrons is assumed ad hoc.
The rise at high energy in the positron fraction originally discovered by PAMELA in 2009 [17], and subsequently confirmed by Fermi-LAT and AMS-02 [18] collaborations, constitutes then a substantial deviation from the standard prediction of Eq. (19) and appears robust with respect to uncertainties in CR transport models, implemented in a more realistic way (see, however, [203]). The release of the data on the absolute positron spectrum [204] confirmed and strengthened this conclusion.
The detection immediately triggered a debate in the community (see e.g. [205,206] and references therein). A natural explanation in terms of nearby astrophysical accelerators of primary e + +e − pairs, e.g. pulsar wind nebulae (already invoked in [92] as potential contributors to the leptonic flux), was soon considered as a very promising one; see Fig. 3 for a recent realization of this scenario, compared to up-to-date experimental data. Other astrophysical interpretations were proposed (see, e.g., [207]), including the already mentioned secondary production at accelerators [166].
On the other hand, many DM scenarios were invoked as well: The tough challenges for model building are in this case the large annihilation cross section required to sustain the measured positron flux at high energy, and the strong constraints originating from other channels (including gamma rays, CMB, and antiprotons); we refer [208] for an early review on the topic.
More interestingly, nowadays it is possible to challenge the widely debated pulsar hypothesis in several ways, and the uncertainties in CR transport play a major role.
First of all, it is possible to look for an anisotropy in the arrival direction of high-energy leptons; moreover, the gamma-ray observatories may now allow to identify the emission from the leptons leaving nearby known pulsars. Along this track, a detection of a TeV halo around Geminga has recently been reported in [209]: In that paper a naive estimate of the diffusion coefficient in the vicinity of Geminga is presented, which turns out to be much smaller than the average Galactic one inferred by secondary-to-primary ratios, posing a challenge both to CR transport models and to the pulsar interpretation of the positron anomaly as well; see also the follow-up detailed discussion in [210].
Very recently, the antiproton and positron channels were critically re-examined in [211]. In that paper it is noticed that the ratio between the positron (or antiproton) flux to the proton one is consistent with the secondary production rates in the conventional picture.
Based on these considerations, the author suggests that Galactic positrons and antiprotons may have a common origin as secondaries in hadronic interactions, probably produced in the local interstellar environment, so that diffusion and energy loss do not act for enough time to leave an observable imprint on the spectrum. If confirmed, this would imply a completely different propagation scenario characterized by a much lower residence time (∼ 1 Myr) compared to current benchmark values: Such scenario would also accommodate the spectral break in the electron spectrum reported by H.E.S.S. [212], but not the electron slope, which is actually steeper than the antiproton and positron one: In such alternative framework, the e − /p discrepancy could in principle be generated at source, by not-yet identified mechanisms.
The take-home message of this discussion is that the positron channel is far from being understood, and the nature of the emission above 30 GeV remains mysterious. However, DM interpretations of this anomaly seems disfavored with respect to several alternative astrophysical scenarios, in particular the pulsar hypothesis. Possible avenues towards a clearer understanding of these issues are: 1) More detailed studies of the leptonic CR anisotropy (that can in principle provide a smoking gun of the pulsar scenario, or in alternative strongly constrain the scenario itself); 2) More data beyond TeV energies: experiments such as DAMPE 5 and CALET 6 are already operating, and the first results from DAMPE [213] already showed some interesting features to be confirmed and further investigated [202] (as shown in Fig. 3); 3) More investigations of the interplay with the high-energy gamma-ray observations, such as the TeV halo around Geminga, aimed at characterizing the diffusion properties of leptons in the vicinity of the sources; 4) Also in this case, a better characterization of the diffusion halo size, for instance by means of analyses focused on the radio emission, following e.g. the approach of [156,157].

GAMMA-RAY OPPORTUNITIES
The ubiquitous flux of high-energy CR nuclei and leptons may be able to transform the Galaxy into a huge pion factory and an efficient machine to up-scatter diffuse photons emitted by stars and reprocessed by dust grains. These processes -namely π 0 decay and Inverse Compton scattering, with the addition of the (usually sub-dominant) bremsstrahlung emission -yield a diffuse flux of high-energy photons from the MeV to the multi-TeV energy domain, reaching the current sensitivity of space missions such as Fermi-LAT and AGILE, and ground-based facilities such as H.E.S.S., MAGIC, VER-ITAS and HAWC.
Therefore, the observable gamma-ray sky can give us today a quite unique diagnostics of CR transport far from the solar system environment. The Galactic diffuse gamma-ray emission stemming from CR interactions with the ambient gas and radiation field constitute indeed the bulk of photons measured along the Galactic plane region [214]. It depends on our observational knowledge of emitting targets [215], namely the indirect tracing of gas column densities [216,217] through, e.g., observed CO emissivities, or the characterization of the low-energy photon background, the so-called interstellar radiation field [218]. Moreover, it crucially relies on the details about CR propagation across the Galaxy [219].
In the last few years important progress has been pursued in the development of phenomenological viable models for Galactic CR propagation, able to match the observed GeV -TeV photon data from the Galactic plane region, while reproducing local CR measurements [116,117]. In the next future, a more systematic study of gamma-ray data in symbiosis with the analysis of local CR observables may offer us the most important chance to pin down the exact features underlying Eq. (11) (see, for example, [115,220]) in a data-driven fashion [169,[221][222][223].
Looping over uncertainties both on the side of emitting targets and possibly also on the underlying CR transport properties is a challenging task [214,224]. However, such an attempt is particularly welcome in order to constrain particle DM properties [219,225,226]. As previously mentioned, anomalies in the gamma-ray sky may be extremely compelling for indirect DM searches. State-of-the-art N-body simulations (see, e.g., [227,228]) predict an extended DM halo embedding and surrounding the Milky Way, with a central density peaking in correspondence to the Galactic center (GC). Then, within a scenario where DM particles pair-annihilate (or decay) eventually to gamma-ray photon yields, the GC region is likely the brightest possible target for DM indirect searches, being relatively close to us (∼ 8 kpc).
The kinematics of DM thermal relics annihilating today in the halo may actually give rise to peculiar photon energy spectra. In particular, gamma-ray photon lines would not have any well-known astrophysical counterpart. Therefore, mono-chromatic lines in the GeV -TeV range produced by DM particle pairs annihilating into two-body final state channels with one or two photons, give rise to a potential smoking-gun signature. However, in a standard scenario of electrically neutral particles [229], DM rate to mono-energetic photons will exhibit loop-suppression (see, e.g., [230,231] for the benchmark of neutralino DM). Interestingly, a hint in favor of such a spectacular signature has been found in 2012 from dedicated analyses of Fermi-LAT data in an extended region of the Galactic Center [232,233], showing a peak in the photon spectrum at an energy ∼130 GeV with significance of the excess at ∼ 3σ level. After an optimal observational strategy has been carried out for the purpose [234,235], the updated analysis from the Fermi-LAT collaboration does not support any longer the original evidence for such a spectral feature in the dataset [236], suggesting previous claims to be related only to a statistical fluke. While current lack of detection of gamma-ray lines place important upper-limits on today's DM annihilation cross-section/decay rate into mono-chromatic photons in the GeV -TeV energy window [236,237], the search for pronounced spectral features in the gammaray sky remains one of the most tantalizing observational programs within the WIMP mass reach [26].
Interestingly enough, within almost an entire decade, an increasing number of studies focused on the GC region, repeatedly showing the existence of a statistically significant signal in Fermi-LAT data possibly correlated with spectral and morphological features of prompt gamma-ray emission from DM pair annihilation [238,239,[241][242][243][244][245][246][247][248][249][250][251][252][253][254][255][256][257]. Originally, the analysis of Fermi-LAT data in 2009 concerning the innermost two degrees around the GC lead the authors of [242,243] to discover the existence of a bump-like feature in the photon spectrum exceeding the estimated astrophysical background with high statistical significance. Peaking around ∼ 3 GeV and fitting an approximately spherical morphology, the signal found could be immediately associated to annihilating DM particles with mass range and cross-section remarkably within the WIMP ballpark. Strengthening such compelling interpretation, [247], and successively [249], focused on a region of interest (ROI) extended up to ten degrees in latitude, in correspondence to the low-latitude part of the Fermi bubbles [258], finding new evidence for the gamma-ray excess even at few kpc of distance from the GC. Eventually, in [252] an optimized analysis implementing specific cuts to Fermi-LAT events -improving the resolution of the gamma-ray mapscould reach a statistical preference at the ∼ 30σ level for the inclusion of a WIMP-like template in the fit to the dataset. These claims triggered several phenomenological studies on the New Physics direction to undertake in order to explain this anomaly (see, e.g., [259][260][261][262][263][264][265][266]).
It is important to note that a close investigation about the impact of the Galactic diffuse emission -related to CR physics -on the robustness of the aforementioned evidence in favor of a DM indirect detection came only four years after [242], with the studies in [253,254]. In particular, [254] analyzed the uncertainties related to the galactic diffuse modeling explicitly inspecting the systematics of 60 models with different characterization of CR transport physics and emitting target properties. Most importantly, the authors in [254] derived in a more data-driven way the overall systematics in the ROI of 20 • × 20 • centered at the GC by looking at at the residuals obtained from a large number of control ROI along the whole Galactic plane, using their 60 galactic diffuse emission models for a principal component analysis. The systematics found in [254] associated to the GC excess signal enlarged the set of viable DM bestfit scenarios [263,264], and opened a new window for the interpretation of the anomaly in terms of a population of point sources such as millisecond pulsars (MPs) [267], a scenario originally proposed in [268] and successively supported in [246,248,251,269], while criticized in [270,271]. Alternative astrophysical explanations in terms of outburst events of hadronic [272] or leptonic [273] origin turned out to be not favored by data [274].
Notably, the different CR injection terms (Q 0 appearing in Eq. (11) adopted in all the aforementioned analyses, correspond to the radial distributions shown by the thin and dashed colored lines in the lower-left panel of Fig. 4. Their behavior in the GC proximity does not correlate well with the expected high-star formation rate present in the inmost few hundred parsecs around the GC [275], the so-called Central Molecular Zone (CMZ). Multi-wavelength observations of the CMZ [276] point indeed to an environment with large amount of molecular gas [277], and hosting high-mass OB stars [278], potential progenitors for standard CR acceleration sites such as supernova remnants [279]. These pieces of informations motivated [238] to re-analyze the GC excess implementing a novel steady-state CR source term capturing the CMZ star-forming activity. The three upper panels in Fig. 4 from [238] show the residual count maps of photons falling in the energy range of 1 -10 GeV, within a ROI of 10 • × 10 • centered at the GC. From left to right, the result within a benchmark scenario originally identified in [254] to optimize the study of the GC anomaly at the GeV, and the remarkable improvement in the outcome of the fit to Fermi-LAT data when adding a DM component (central panel), or when re-visiting the whole Galactic diffuse emission background on the basis of the novel steady source term peaked at the GC (right panel).
Notably, this result has been successively corroborated and refined in [239,240], with the implementation of the high-resolution Galactic gas distribution obtained in [216], able to resolve spiral arms and the Galactic bar, and most importantly providing kinematic resolution towards the GC, with the help of hydrodynamic simulations. Then, these HII density maps were correlated to the CR source injection term by means of a simple model of star formation [280]. The thick light-blue line in the left-bottom panel of Fig. 4, taken from [239], shows a benchmark realization of the novel CR source radial distribution, that correctly does not fall off any longer at small Galactocentric radii. The red dashed (blue dotdashed) line in the right-bottom panel, from Fig. 4 of [240], identifies the improvement, ∆χ 2 < 0, of the description of gamma-ray data in the inner Galaxy region as a function of the fraction f H2 of CRs injected with a spatial distribution tracing the HII density maps of [216], FIG. 4. Importance of the cosmic-ray contribution to the Galactic Center excess at few GeV from the analysis of Fermi-LAT data. The top three panels refer to the residual photon count maps shown in Fig. 3 of [238]. In the left-bottom panel, different source distributions for primary CR injection (credit to [239], Fig. 1). The right-bottom panel shows how the quality of the fit to gamma-ray data can improve according to a better modeling of the CR primary source distribution (credit to [240], Fig. 4).
in the presence (absence) of the DM component, against the best-fit scenario of a Galactic diffuse emission model built with a standard CR injection term plus the DM component. Inset numbers report the statistical preference of the red line over the blue one.
The relevance of the Galactic diffuse emission modeling in assessing the significance of the GeV excess at the GC has been further acknowledged by the Fermi-LAT collaboration [281]. In [281], injected CR electrons as in [239], together with the uncertainties on the interstellar radiation field, have been shown to play a major role in the characterization of the gamma-ray anomaly. As also marked by the right-bottom panel of Fig. 4, noisy residual photon counts around the GC in the analyses of [238,239,281] still leave room for an extra-component in the description of gamma-ray data. Notably, two independent studies implementing different statistical techniques for clustering patterns in the observed photon count maps have followed, showing very high evidence in favor of a hitherto undetected population of point sources, able to fully account for the GeV gamma-ray excess [282,283]. The results of these two works reached remarkably similar conclusions, giving substantial credit to the MSP-like interpretation of the gamma-ray signal at the GC (see [284,285] for discussions about the luminosity function of the putative MSPs at the GC) and triggering relevant dedicated searches [286]. Non-negligible mismodeling of astrophysical backgrounds and foregrounds may affect the details on the prediction for such unre-solved population of point sources within the Galactic bulge, see e.g. [255]. A recent novel tool, SkyFACT, developed in [287], based on image reconstruction and adaptive spatio-spectral template regression, has allowed for dramatic improvements in the quality of the fits to gamma-ray data through fine-grained variations of galactic diffuse emission modeling. Exploiting this powerful package, a novel investigation on the morphology of the excess in connection to the stellar distribution of the bar/boxy bulge in the inner Galaxy provides strong support to the MSP hypothesis [288].
In conclusion, the GC GeV excess remains a widely studied signal in the astro-particle community. The updated comprehensive analysis carried out by the Fermi-LAT collaboration [256] fairly summarizes the most relevant factors affecting the characterization of the signal: 1) the details on CR production and propagation, especially in the GC proximity; 2) the templates for the interstellar gas and radiation field in the inner Galaxy; 3) the emissivity and morphology of the Fermi bubbles towards low latitudes; 4) the list of point sources near the GC identified within a given background model.
The take-home message of this long debate is likely twofold: The existence of an extended emission from the inner Galaxy peaked at few GeV is well established; however, the characterization of this emission, i.e. the morphology and its intensity, strongly depends on the assumptions of the CR source distribution. Interestingly, the interpretation in terms of unresolved point sources, possibly associated to a population of millisecond pulsars -currently supported by wavelet and photon-count statistics analyses -is testable in the future with more sensitive radio facilities [289]. As a final outlook, let us mention that a recent spatially extended gamma-ray signal from the center of M31 galaxy [290] has further renewed the interest on the GC anomaly and, for instance, may possibly shed new light on the existence of galactic bulge MSP populations [291].
Let us now move away from the GC region, and consider another very relevant potential discovery window. Contrary to the complex astrophysical environment characterizing the CMZ and the GC, dwarf spheroidal galaxies of the Milky Way (dSphs) stand out as very promising targets in the gamma-ray band due to the corresponding low astrophysical background and foreground [292,293]. Being relatively close to us and associated to fairly large DM densities, the gamma-ray campaign on dSphs has been soon realized to be one of the most potentially sensitive probes to particle DM properties [294][295][296]. At present, dSph gamma-ray upper-bounds are remarkably probing the benchmark thermal relic scenario within the WIMP mass window [297][298][299][300]. Such upper-limits may be at odd with naive DM interpretations of the GC excess [301,302], while depend crucially on the estimated DM content in these galaxies, potentially affected by several systematics [303][304][305][306].
Interestingly, hints for a gamma-ray signal possibly compatible with the one observed at the GC have been found in the analysis of some of the most recently discovered Milky Way ultra-faint satellites, see, e.g., the case of Reticulum II [307][308][309]. According to the latest joint analysis of Fermi-LAT and Dark Energy Survey collaborations, the significance of these excesses remain at present well below the 3σ level [299]. Moreover, dedicated searches in the radio-band have not found any counterpart of the putative gamma-ray DM signal [310]. Note that -from the perspective of a signal detection -a broad multi-wavelength program for indirect DM searches in dSphs would be indeed promising [311,312]. However, in opposition to the case of gamma-rays, constraints on DM annihilation/decay derived from the observation of dSphs in the radio and/or X band turn out to be subject to a larger set of astrophysical uncertainties, including CR transport physics [313]. The future of indirect searches is particularly bright.
In the gamma-ray band, two new regions of the spectrum will be explored. On the low-energy side, the MeV-GeV domain can be probed by planned experiments such as e-ASTROGAM [314] 7 , and AMEGO 8 (other previously proposed missions include, e.g., COMPAIR [315] and ADEPT [316]), which could be realized in the midand long-term future in the late 2020s). All those experiments will feature a 2-3 order-of-magnitude increase in sensitivity [317,318] with respect to previous instruments operating in this window like COMPTEL and EGRET [319,320], and a remarkable energy resolution especially below 10 GeV, where the detection principle is based on Compton scattering instead of pair production. This will guarantee a high constraining power as far as sharp spectral features (e.g. lines originating from DM annihilation) are concerned [321]. On theoretical grounds, vanilla WIMP scenario in this domain of energy scales may not be adequate due to the Lee-Weinberg bound [56,57]. Interesting models designed to offer an explanation for a 511 keV emission line detected in the inner Galaxy may be instead well probed [322,323].
On the other end, in the TeV domain, while many Imaging Atmospheric Cherenkov Telescopes (MAGIC [324], H.E.S.S. 9 and VERITAS 10 ) and air-shower arrays (such as Milagro 11 and HAWC 12 [325]) have already been providing relevant results, and have provided stringent constraints on DM annihilation [326], CTA will provide a further, unprecedented increase in sensitivity above the TeV. We want to emphasize that, as mentioned in Section 2, the most naive version of the WIMP miracle naturally yields a predicted mass for the DM candidate in that ballpark. CTA can probe the DM sector in different ways [327][328][329][330].
In connection with the main topic of this review, we stress once again that the GC region appears promising also in this context, which implies that a more careful modeling of the γ-ray diffuse emission from the inner Galaxy will be required. The diffusion properties, again, play then a central role: In particular, the CR spectrum in the inner Galaxy is currently a matter or debate, and Fermi-LAT data seem to hint towards a hardening in the inner Galaxy [117,222,223]. DRAGON-based models [115,117] featuring inhomogeneous or anisotropic diffusion, which are key features of this specific numerical package, allow to reproduce Fermi-LAT data and provide TeV predictions [118,221] which seem in accord with a collection of data from different experiments, including the bright, diffuse Galactic Ridge emission measured by H.E.S.S. over the latest decade [331][332][333]. These kind of studies will be even more important in the forthcoming years, due to the necessity, that we have stressed several times along this paper, of a good characterization of the astrophysical backgrounds. The study of dwarf spheroidal galaxies will also benefit from the expected increase in sensitivity in the TeV domain, and the expected performance of CTA in this channel is remarkable [334][335][336].
As far as charged particle channels are concerned, the antiparticle/antinuclear avenue is particularly promising despite the relevant uncertainties. In fact, new highquality data are expected by many experiments: besides AMS-02, the already operating DAMPE [337] space experiment and CALET [338], on board the ISS, have the opportunity to probe in particular the leptonic channel all the way up to the multi-TeV range with unprecedented energy resolution and sensitivity; planned experiments such as HERD 13 will provide a further, very relevant extension in the covered energy range and sensitivity, and -on the antinuclei side -we are looking forward to experiments such as GAPS [339], a balloon experiment expected to operate in 2020 and optimized specifically for low-energy antinuclei signatures, thanks to a novel detection technique based on exotic atom capture and decay.
To conclude this grand tour, let us emphasize once again the opportunity offered by the tremendous increase in experimental accuracy we are already witnessing in these days, and we will be further exploiting with the upcoming years. In our opinion, in order to take full advantage of these developments, interdisciplinarity is the main avenue. The guaranteed outcome of the research program aimed at indirectly detecting particle dark matter by means of astronomical data certainly involves a more profound understanding of the CR physics outlined across this review. It is therefore crucial to deepen and broaden the connections between the research fields we have described above: Both the high-energy astrophysics and the particle physics community would unquestionably benefit from a cutting-edge increasing crossover.