Heavy dark matter particle annihilation in dwarf spheroidal galaxies: radio signals at the SKA telescope

A weakly interacting dark matter candidate is difficult to detect at high-energy colliders like the LHC, if its mass is close to, or higher than a TeV. On the other hand, pair-annihilation of such particles may give rise to $e^+ e^-$ pairs in dwarf spheroidal galaxies (dSph), which in turn can lead to radio synchrotron signals that are detectable at the upcoming Square Kilometre Array (SKA) telescope within a moderate observation time. We investigate the circumstances under which this complementarity between collider and radio signals of dark matter can be useful in probing physics beyond the standard model of elementary particles. Both particle physics issues and the roles of diffusion and electromagnetic energy loss of the $e^\pm$ are taken into account. First, the criteria for detectability of trans-TeV dark matter are analysed independently of the particle physics model(s) involved. We thereafter use some benchmarks based on a popular scenario, namely, the minimal supersymmetric standard model. It is thus shown that the radio flux from a dSph like Draco should be observable in about 100 hours at the SKA, for dark matter masses upto 4-8 TeV. In addition, the regions in the space spanned by astrophysical parameters, for which such signals should be detectable at the SKA, are marked out.


I. INTRODUCTION
Some yet unseen weakly interacting massive particles (WIMP) are often thought of as constituents of dark matter (DM) in our universe. Many scenarios beyond the standard model (SM) of particle physics including them are regularly proposed, with various phenomenological implications. They are constrained by the data on relic density [1] of the universe as well as various direct search experiments [2,3]. It is expected that a WIMP DM candidate should also be detected at the Large Hadron Collider (LHC) in the form of missing transverse energy (MET) (see for example [4][5][6]).
Detectability at colliders and also in direct search experiments, however, depends on the mass as well as the interaction cross-section of the DM particle. In particular, detection becomes rather difficult if the WIMP mass approaches a TeV [7]. For a weakly interacting DM, production at the LHC has to depend mostly on Drell-Yan processes where the rate gets suppressed by the s, the square of the subprocess centre-of-mass energy, and also by the parton distribution function at high x. An exception is the minimal supersymmetric standard model (MSSM) [8]; there the production of colored superparticles (squarks/gluions) via strong interaction may be followed by decays in cascade leading to DM pairproduction. However, here, too, the event rates go down drastically when one looks at the current constraints on colored superparticle masses as well as the various other limits on the MSSM spectrum. On the whole, while a near-TeV DM particle still admits of some hope at the LHC in the MSSM [7,[9][10][11], the reach is considerably lower for most other scenarios where the 'dark sector' is at most weakly interacting [12,13]. It is therefore a challenge to think of additional indirect evidence, if a trans-TeV DM particle has to be explored.
While the prospect of radio fluxes unveiling DM annihilation has been explored in earlier works [29][30][31][32][34][35][36], it was pointed out in reference [37] that SKA opens up a rather striking possibility. The annihilation of trans-TeV DM pairs in dwarf spheroidal galaxies (dSph) lead to electron-positron pairs which, upon acceleration by the galactic magnetic field, produces such radio synchrotron emission. dSph's are suitable for studying DM, since star formation rates there are low [38,39], thus minimising the possibility of signal originating from astrophysical processes. Their generic faintness prompts one to concentrate on such galaxies which are satellites of the Milky Way. The SKA can ensure sufficient sensitivity required to detect the faint signal from the sources, and at the same time, will have high enough resolution to remove the foregrounds. It was shown in [37] that about 100 hours of observation at the SKA can take us above the arXiv:1905.11426v2 [hep-ph] 9 Jan 2020 detectability threshold for radio signals from the annihilation of DM particles in the 5-10 T eV range. Of course, the compatibility of such massive WIMP with the observed relic density requires a dark sector spectrum with enough scope for co-annihilation in early universe, as was demonstrated in [37] in the context of the MSSM. What one learns from such an exercise is that the probe of at least some DM scenarios should thus be possible on a time scale comparable with the running period of the LHC, and that the reach of the LHC [7,[9][10][11] for WIMP detection may be exceeded considerably through such a probe. It was also found that cases where SKA could observe radio fluxes from a dSph were consistent with limits from γ-ray observations as well as antiparticles in cosmic rays.
In order to ascertain which scenarios are more accessible in such radio probes, one needs to understand in detail the mechanisms whereby high-mass DM particles can produce higher radio fluxes, and also the effects of astrophysical processes that inevitably affect radio emission. This is the task undertaken in the present work.
The spectrum as well as the dynamics of the particle physics scenario, along with the DM profile in a dSph, is responsible for DM annihilation as well as the subsequent cascades leading to electron-positron pairs. 1 The electron(positron) energy distribution at the source function level is also the determined by the above factors [40][41][42]. However, they subsequently pass through the interstellar medium (ISM) of the galaxy, facing several additional effects. These include diffusion as well as electromagnetic energy loss in various forms, including Inverse Compton effect, Synchrotron loss, Coulomb effect and Bremsstrahlung effect [40,41,43,44]. Besides, the galactic magnetic field is operative all along. The way these affect the final e + e − energy distribution is highly inter-connected and nonlinear. For example, while the magnetic field causes electrons to lose more energy in synchrotron radiation, it puts a check on the reduction of the flux through diffusion by confining them longer within the periphery of the dSph. This check applies to electrons at lower energies at the cost of those at higher energies. Also, electromagnetic energy loss enhances the population of low-energy electrons and positrons, the enhancement being more when they have higher kinematic limits, enabled by higher mass of the DM particle. Such low-energy e + e − pairs enhance the flux in the frequency range appropriate for a radio telescope.
In the following sections, we analyse the various ingredients in the radio flux generation process, as outlined above. We first do this for fixed DM particle masses, 1 In principle, electron-positron pairs may also be directly produced in a complete model-independent scenario. We have referred to cascade since our benchmark correspond to bb, tt, W + W − and τ + τ − as dominant annihilation channels, inspired by the MSSM where direct e + e − is not found to dominate for high DM masses.
and for values of their annihilation cross-sections fixed by hand. The relative strengths of the effects of the particle theoretical scenario as well as diffusion and radiative process are thus assessed. This also serves to evolve an understanding of the dependence on the diffusion coefficient, parameters involved in electromagnetic energy loss, and, of course, the strength of the galactic magnetic field. If the nature of the DM particle(s) is known in independent channels, then the observation of the observed radio flux may be turned around to improve our understanding of these astrophysical parameters, using the results presented here.
Finally, we use some theoretical benchmark points to demonstrate the usefulness of the SKA in probing trans-TeV DM. A few sample MSSM spectra are used, largely because they offer the scope of co-annihilation that is so essential for maintaining the right relic density [37]. Using the minimum annihilation cross-section required for any DM mass for detection at SKA (in 100 hours), and the maximum value of this cross section that is compatible with limits from γ-ray and cosmic-ray data, we show that several benchmark points with DM mass in the 1-8 T eV range, which are yet to be ruled out by any observation can be investigated in 100 hours of SKA observation.
We have organised the paper as follows. In section II we have discussed in somewhat brief manner the calculations of synchrotron flux originating from DM annihilation inside a dSph. In section III we have analysed the effects of various astrophysical parameters in different steps of the production of radio flux. Section IV describes the features of heavy DM. In section V we have shown the detectability curves or threshold limits for observing radio flux in SKA, in both model independent as well as model dependent way. We have also shown the final radio fluxes for some theoretical benchmarks. Finally in section VI we conclude.

II. ESSENTIAL PROCESSES: RECAPITULATION
We start with a brief resume of sequence of processes that leads to radio flux from a dSph, as discussed, for example, in [30,31,40,41,43,44]. Dark matter pair annihilation inside of a dSph can produce SM particle pairs such as bb, τ + τ − , W + W − , tt etc. These particles then cascade and give rise to large amount of e ± flux. The resulting energy distribution of which can be obtained from the source function [41,42]: where σv and N pairs (r) = ρ 2 χ (r) 2m 2 χ are respectively the velocity averaged annihilation rate and the number density of DM pairs inside the dSph. m χ is the DM mass and ρ χ (r) is the DM density profile in the dSph as a function of radial distance r from the centre of the dSph. For our analysis we have taken Draco 2 dSph assuming an Navarro-Frenk-White (NFW) profile [48], with where ρ s = 1.4 GeV.cm −3 and r s = 1.0 kpc [41,44].
dE B f estimates the number of e ± produced with energy E per annihilation in any of the aforementioned SM channel (f ) having branching fraction B f . Figure 1 shows the electron energy distribution for four different annihilation channels (bb, τ + τ − , W + W − , tt) and for two DM masses, namely, m χ = 300 GeV (upper left panel) 5 T eV (upper right). For all cases σv has been assigned a fixed value (10 −26 cm 3 s −1 ), and the branching fraction corresponding to each channel has been set in turn at 100%. In addition, a comparison between cases with the two masses has also been shown in the lower panel for the bb and τ + τ − channel. The prediction for tt and W + W − falls in between those curves. All these energy distributions are obtained using mi-crOMEGAs [49,50]. Further discussions on these curves will be taken up in section IV.
The electrons produced in DM pair annihilation, diffuse through the interstellar medium (ISM) in the galaxy and lose energy through various electromagnetic processes. Assuming a steady state and homogeneous diffusion, the equilibrium distribution dn dE can be obtained by solving the equation [31,40,41,43,44,51] where D(E) is the diffusion parameter. The exact form of D(E) is not known for dSphs, hence for simplicity, we assume it to have a Kolmogorov spectrum D(E) = D 0 E GeV 0.3 , where D 0 is the diffusion coefficient [32,41,44]. Very little is known about the value of D 0 for dSphs due to their very low luminosity. A choice like D 0 = 3 × 10 28 cm 2 s −1 is by and large reasonable for a dSph such as Draco [44], although higher values of D 0 , too, are sometimes used as benchmark [36]. We have mostly used D 0 = 3 × 10 28 cm 2 s −1 , and demonstrated side by side the effect of other values of this largely unknown parameter. As explained in [30,31,52], assuming a proper scaling, D 0 for a dSph can also in principle have a value one order lower than that in the Milky Way [53]. However, the present study is not restricted to this value, mainly D 0 3 × 10 28 cm 2 s −1 . The detectability of the radio flux for higher D 0 , too, has been studied by us, as will be seen when we come to Figure 13.
The energy loss term b(E) takes into account all the electromagnetic energy loss processes such as Inverse Compton (IC) effect, Synchrotron (Synch) radiation, Coulomb loss (Coul), bremsstrahlung (Brem) etc. and their combined effect can be expressed as [40,41,43] b  [40,41,43]. In a dSph like Draco, the first two terms (i.e. IC and Synch) in Equation 4 dominates over the last two terms (i.e. Coul and Brem) for E > 1 GeV [31]. The energy loss term b(E) depends on the galactic magnetic field (B) through the synchrotron loss term which goes as B 2 . It is extremely difficult to get insights (say, through polarization experiments) for the magnetic field properties of dSph. The lack of any strong observational evidence suggests that the magnetic fields could be, in principle, extremely low. On the other hand, there could be various effects which can significantly contribute to the magnetic field strengths in dSph. In fact, numerous theoretical arguments have been proposed for values of B at the µG level. It has been argued, for example, that the observed fall of B from the centre of Milky Way to its peripheral region can be linearly extrapolated to nearby dSph's, leading to B 1 µG. The possibility of the own magnetic field of a dSph has also been suggested. For detailed discussions we refer the reader to reference [54]. For the dSph considered in this work, we have mostly used B = 1 µG [32,41,44], but more conservative values like B ∼ 10 −1 − 10 −3 µG have also been considered. The manner in which diffusion and electromagnetic processes affect our observables will be discussed in detail in sections III and V. The overall potential of the SKA in probing regions in the D 0 − B space in a correlated fashion will be reported in section V. Equation 3 has a solution of the form [31,40,41,44] dn with the boundary condition dn dE (r h ) = 0, where r h is the radius of the diffusion zone. For Draco, r h is taken to be 2.5 kpc [41,44]. G(r, ∆v) is the Green's function of the equation and has a form √ ∆v determines the distance traveled by an electron as it loses energy from E to E. For small galaxies like dSphs, the mean value of this length scale is expected to be larger than r h even for non-conservative choices of D 0 and B. We shall discuss this in detail in section III.
After interacting with the magnetic field B present inside the galaxy, the produced electron/positron distribution dn dE will emit synchrotron radiation (with frequency ν) at a rate governed by the synchrotron emission power P Synch (ν, E, B) [40,41,43,44,51].
In Figure 2, we have shown the dependence of the synchrotron emission power on the electron energy E for two different magnetic fields (B = 1 and 0.1 µG) and for two frequencies (ν = 10 and 10 4 M Hz). For each frequency shown, it is clear that a stronger magnetic field always intensifies the power spectrum.
The synchrotron emissivity J Synch as a function of frequency ν and radius r is obtained by folding dn Finally, the approximate radio synchrotron flux can be obtained by integrating J Synch over the diffusion size (r h ) of the dSph [40,41,43,44,51] 10 -1 10 0 10 1 10 2 10 3 10 4 where L is the luminosity distance of the dSph. For Draco, L ∼ 80 kpc [41,45]).

III. EFFECT OF VARIOUS ASTROPHYSICAL PARAMETERS
Let us now take a closer look at the astrophysical effects encapsulated in equations 3 and 4. Such effects are driven by the diffusion coefficient D 0 and the electromagnetic energy loss coefficient b(E), which determine the steady state e ± distribution dn dE (r, E) for a given dN e dE . As has been already mentioned, b(E) is implicitly dependent on the galactic magnetic field B.
We use again the two DM masses (m χ = 300 GeV and 5 T eV ) for studying the effects of D 0 and b(E). To have some idea on these effects separately, as well as their contribution in an entangled fashion, we have considered the implication of equation 3 for three different scenarios: For the NSD scenario, solution of equation 3 has a simpler form [40]: Note that, b(E) increases with E (from equation 4). Its effect on the number density in different energy regions is more prominent for this NSD scenario where we neglect the effect of distribution. The number density decreases in the high energy region due to the combined effect of the 1 b(E) suppression as well as the lower limit of the integration in equation 10. This can also be seen from Figure 3 where we have plotted the ratio of dn dE and the source function Q e (which determines the initial electron flux due to DM annihilation) against E. The term dn dE 1 Qe (in unit of s −1 ) essentially determines how the shape of the initial distribution Q e gets modified due to the effect of various astrophysical processes. As discussed in section II, the energy loss term b(E) depends on the square of the magnetic field B through synchrotron loss. Thus b(E) increases with increase of B, which in turn will reduces the dn dE at all E. This phenomenon can also be observed in Figure 3 where the red and magenta lines indicate the distributions for B = 1 and 10 µG 3 , respectively. For comparison, we have also shown the cases for B = 0.1 µG (green lines). Note that, the later case coincides with the case for B = 1 µG. This is due to the fact that the magnetic field dependence (through synchrotron loss) in the energy loss term b(E) gets suppressed by the Inverse Compton (b IC ) term for lower B (B < 1 µG), which can be clearly seen from equation 4. Also, for heavier DM masses the energetic electrons are produced in greater abundance (as seen from Figure 1) which can lead to a larger dn dE 1 Qe . The effects of two different DM masses, m χ = 300 GeV and 5 T eV , have been shown in this figure by dashed and solid lines respectively. Now, for the Nb scenario, the solution becomes with the same boundary condition as what we assumed for equation 5. Here Since the effect of b(E) is absent in this case, the distribution dn dE will not depend on the magnetic field B.  Qe for the Nb scenario will not depend on the initial spectrum Q e and its energy profile will follow the energy dependence of We should mention here that neglecting diffusion is not a bad assumption where the length scale ( √ ∆v) over which the e ± losses energy is much shorter than the typical size of the system [40,41]. However, for smaller systems like dSphs the effect of diffusion cannot be neglected. This can be justified from Figure 4 where we have plotted the SD+b case (i.e. taking diffusion into account along with energy loss effect; shown by red curves) along with the two previously mentioned scenarios. It is clear that the NSD scenario is strikingly different from the one which takes diffusion into account. Addition of diffusion in the system essentially suppresses the e ± distribution, especially in the low energy energy region. As a result, the plots for the SD+b case closely match the Nb one at lower energy, but diverges following the NSD scenario at higher energies where the effect of b(E) dominates. This phenomenon can also be explained in terms of the diffusion length scale ( √ ∆v) which is a result of combined effects of diffusion and energy loss (equation 7). It is indicative of the length covered by an electron as it loses energy from the source energy E to the in-teraction energy E which is typically larger than the size of the dSph. In order to determine the minimum energy that an electron can possess before escaping the diffusion radius, we have plotted the interaction energy E vs. √ ∆v in Figure 5 for different sets of D 0 and B. The left and right panels correspond to the source energy E = 1 T eV and 100 GeV respectively. The diffusion radius (for Draco) r h = 2.5 kpc is shown by the black solid lines for both cases. The minimum energy is thus calculated from the point where the model crosses this line. A higher D 0 will increase this minimum value of interaction energy and as a result, it reduces the e ± number density by making them leave the diffusion zone earlier. Note that, the cases with B = 0.1 µG practically coincides with that for B = 1 µG. This is due to the fact that b(E) remains unchanged for B ≤ 1 µG, as already been mentioned. On the other hand, we have explicitly checked that, a larger value of B (say B = 10 µG; not shown in the plot) will help the electron in losing energy more quickly before escaping the diffusion zone, causing a suppression in the number density of high energy electrons and an enhancement of the low energy electrons.
Similar conclusions can also be drawn from Figure  6. To illustrate the effects of D 0 and B in the SD+b scenario, we have plotted dn dE 1 Qe against E for the DM masses 5 T eV (upper panels) and 300 GeV (lower panels) at two different radii 0.1 kpc (left panels) and 2 kpc (right panels). For all cases, the dominant annihilation channel has been assumed to be bb with σv = 10 −26 cm 3 s −1 . Cases with different D 0 and B are indicated by different color coding as mentioned in the insets of corresponding figures. Similarly, Figure 7 shows the variation of the same with respect to radius r for those two DM masses (m χ = 5 T eV − left panel and m χ = 300 GeV − right panel). For each DM mass we have assumed two different combinations of energies, one being low (E = 0.1 GeV for both DM masses) and the other, high (E = 1 T eV for m χ = 5 T eV and E = 100 GeV for m χ = 300 GeV ). In both panels, the black line indicates the diffusion radius (2.5 kpc) for Draco. One can see that, large D 0 (blue curves) implies a lower density distribution compared to the corresponding cases with smaller diffusion coefficient. The difference is more prominent for low energy.
The effects of diffusion coefficient and magnetic field on the final radio flux S ν can be seen from Figure 8. The radio flux corresponding to higher diffusion coefficient (D 0 = 10 30 cm 2 s −1 ) will be suppressed compared to the lower diffusion case (D 0 = 3 × 10 28 cm 2 s −1 ) due to the escape of a large number of e ± from the stellar object, as discussed above. This suppression is slightly larger in the lower frequency region, since the number density dn dE decreases more in the low energy regime (see Figure 6). For, a constant D 0 , say, 3 × 10 28 cm 2 s −1 , a relatively lower value of B (0.1 µG) reduces the radio flux by about an order of magnitude in all frequency ranges. This is solely due to the fact that, although dn dE has similar values for different values of B (1 and 0.1 µG), the synchrotron power spectrum P Synch decreases with de-  crease in B (as evident from Figure 2). We further emphasise that, though this analysis assumes an NFW profile for Draco, the choice of other profiles such as Burkert [44,55] or Diemand et al. (2005) [56] (hereafter D05) [44] keep our predictions similar, as can be seen in Figure 9.

IV. HEAVY DM PARTICLES
The radio flux obtained in terms of DM annihilation from a dSph crucially depends on the source function Q e (E, r) (equation 1). Let us try to explain why one can get higher radio flux (obtained via dn/dE through integration of Q e (E, r)) for higher DM masses in some cases [37]. For this to happen, Q e corresponding to energetic e ± is intuitively expected to go up for higher m χ . One may consider contributions of several components in the expression of Q e (E, r): • σv : For a trans-TeV thermal DM χ, the factors affecting σv are its mass (m χ ) and the effective couplings to SM particle pairs. While we are concerned here with the annihilation cross sections for χ within a dSph, one also needs consistency with the observed relic density [1]. In general, the expression for the relic density (Ωh 2 ) indicates inverse proportionality to σv , and 1 m 2 χ occurs as the flux factor in the denominator of the latter [57]. Thus a factor of m 2 χ occurs in the numerator. Therefore, m χ 1 T eV may make the relic density unaccept-  possible in the MSSM through the participation of a pseudoscalar of appropriate mass [58]. This also serves to enhance the observed radio flux, as we shall see later. In addition, the annihilation of χ in the early universe may involve channels other than those in a dSph, through co-annihilation with particles spaced closely with it [58]. Once either of the above mechanisms is effective for the relevant particle spectrum, a trans-TeV DM particle remains consistent with all data including the relic density. This ensured, the pair-annihilation of a trans-TeV DM particle in dSph's enable the production of relatively energetic e ± pairs which in turn reinforce the radio signals in the desired frequency range.
The numerator is supplied by observation. For higher m χ , the denominator suppresses the number density of DM in the dSph. Thus this term tends to bring down the flux for higher DM masses.
• dN e f (E) dE B f : This has to have a compensatory effect for higher mass if the suppression caused by the previously mentioned term has to be overcome. Of course, the branching fractions B f in different channels have a role. dN e /dE is often (though not always) higher for higher DM masses. This feature is universal at higher energies. This is basically responsible for a profusion of higher energy elec-trons produced via cascades, after annihilation in any channel has been taken place. More discussion will follow on this point later. It should be noted that dN e /dE represents the probability that an electron (positron) produced from the annihilation of one pair of DM. This gets multiplied by the available number density of DM. Thus, for comparable value of B f in two benchmark scenarios, higher dN e /dE along with higher σv in the dSph (when that is indeed the case), can cause enhancement of Q e (E, r) for higher m χ .
Finally, the available electron energy distribution dn dE (which, convoluted with the synchrotron power spectrum, will yield the radio flux (i.e. equations 8 and 9)), is obtained by solving equation 3 which has a generic solution (equation 5). Therefore, higher values of the source function (Q e (E, r)) can produce higher radio flux at higher m χ .
As mentioned above, the annihilation of heavier DM particles will produce more energetic e ± , as can be seen by comparing the upper-left and upper-right panels of Figure 1 where we have shown dN e /dE in different channels (bb, τ + τ − , W + W − , tt). The lower panel shows the comparison of this e ± spectrum for two DM masses (300 GeV and 5 T eV ) arising from bb and τ + τ − annihilation channels. Note that, the values of dN e /dE for m χ = 300 GeV and 5 T eV are differently ordered for the τ + τ − and bb annihilation channels. The spectrum for bb channel is the steepest one, while τ + τ − is the flattest. In fact, dN e /dE for the bb channel is governed mostly by charged pion (π ± ) decay at various stages of cascade. For the τ + τ − channel, on the other hand, there is a relative dominance of 'prompt' electrons. There remains a difference in the degree of degradation in the two cases. Such degradation is reflected more in the low-energy part of the spectrum, and leads to different energy distributions, which in turn is dependent on m χ via the energy of the decaying b or τ . The reader is referred to see references [59][60][61] for more detailed explanations.
The presence of energetic e ± in greater abundance, which is a consequence of higher m χ , enhances the resulting radio signal obtained through equation 5 and 8. As a result of this, the final radio flux S ν (equation 9) gets positive contribution for higher DM masses compared to that for relatively lower masses for most of the annihilation channels (especially bb, W + W − , and tt). This can easily be seen if one compares the quantity S ν × 2m 2 χ (i.e. removing the effect of the multiplicative factor 1 m 2 χ present in the expression of N pairs ) for two DM masses, 300 GeV and 5 T eV , in the bb annihilation channel (left panel of Figure 10). If one chooses same annihilation rate (e.g. σv = 10 −26 cm 3 s −1 ), the higher DM mass will give higher S ν × 2m 2 χ (mainly in the high frequency region) for this channel. We have explicitly shown this for different choices of astrophysical parameters as indicated by different colors in the figure.
In scenarios where σv corresponding to a particular m χ is calculated from the dynamics of the model, it can happen that for some particular benchmark the annihilation rate ( σv ) is higher for larger m χ (as mentioned earlier and will be discussed further in the next section). In those cases larger σv can at least partially compensate the effect of 1 m 2 χ suppression (coming from N pairs ) for higher m χ . Thus one can get higher radio fluxes (S ν ) for larger DM masses compared to the smaller one.
For the τ + τ − channel the situation is somewhat different. As we have already seen from the lower panel of Figure 1, there is a large degradation of the e ± flux (dN e /dE) in this channel in the low energy region for higher m χ . This degradation in the source spectrum for higher m χ will continue to present in the equilibrium distribution dn dE . After folding this density distribution with power spectrum (see equation 8), the final frequency distribution will be suppressed for higher m χ , mainly in the low frequency range. In the higher frequency range, flux is still high for higher m χ similar to that in other annihilation channel (e.g. bb), because higher m χ corresponds to higher electron distribution in the high energy range. All these phenomena are evident from the red curves (solid and dashed) of right panel in Figure 10. Note that, these curves are for the choice of astrophysical parameters D 0 = 3 × 10 28 cm 2 s −1 , B = 1 µG. If one decreases the magnetic field B to 0.1 µG or increases the diffusion coefficient D 0 to 10 30 cm 2 s −1 , the corresponding effect can be seen from the green and blue curves respectively.
To summarise, one can expect high radio flux for trans-TeV dark matter annihilation from a dSph, based on the following considerations: • The high-mass DM candidate has to be consistent with relic density, for which σv at freeze-out is the relevant quantity.
• Such sizable σv can still be inadequate in maintaining the observed relic density. Co-annihilation channels may need to be available in these cases, though such co-annihilation does not occur in a dSph. This in turn may necessitate a somewhat compressed trans-TeV spectrum.
• The high mass DM candidate should have appropriate annihilation channels which retain a higher population of e ± . This not only off sets the suppression due to large m χ but also enhances through the energy loss term (b(E)) the e ± density at energies low enough to contribute the radio observable at the SKA [33].
• Higher magnetic fields will be more effective in producing synchrotron radiation by compensating the suppression caused by a large galactic diffusion coefficient. As we have checked through explicit calculation, this happens even after accounting for electromagnetic energy loss of the e ± through synchrotron effects.

V. DETECTABILITY CURVES AND FINAL RADIO FLUX
The upcoming Square Kilometre Array (SKA) radio telescope will play an important role in detecting DM induced radio signal from dwarf spheroidal galaxies. Due to its large effective area and better baseline coverage, the SKA has a significantly higher surface brightness sensitivity compared to existing radio telescopes. Since the focus of this work is on the diffuse synchrotron signal from dSphs, we estimate the instrument sensitivity corresponding to the surface brightness and compare with the predicted signal. To calculate the sensitivity for a given source dSph, we need to know the baseline distribution of the telescope. The sensitivity can also be affected by the properties of the sky around the source, e.g., whether there exists any other bright source in the field of view.
Since the detailed design of the SKA is yet to be finalised, it is difficult to estimate the noise in the direction of a given source. For our purpose, however, it is sufficient to estimate the approximate values of the sensitivity using the presently accepted baseline design. In order to do so, we make use of the documents provided in the SKA website [33]. The calculations presented in the document allow us to compute the expected sensitivity near zenith in a direction well away from the Galactic plane for all the frequencies relevant to the SKA, i.e., for both SKA-MID and SKA-LOW. These include contribution from the antenna receiver, the spill-over and the sky background (consisting of the CMB, the galactic contribution and also the atmospheric contribution). Note that we have assumed that all other observational systematics, such as the presence of other point sources in the field, the effect of the primary beam, pointing directiondependence of the sky temperature etc are already cor-rected for and hence are not included in the calculation of the thermal noise.
From the above calculation, we find that typical values of the SKA surface brightness sensitivity in the frequency range 50 MHz -50 GHz for 100 hours of observation time is 10 −6 -10 −7 Jy with a bandwidth of 300 MHz [33,34,37]. This may allow one to observe very low intensity radio signal coming from ultra-faint dSph's. While comparing the predicted signal with the telescope sensitivity, we have assumed that the SKA field of view is larger than the galaxy sizes considered here and hence all the flux from the galaxy will contribute to the detected signal. This assumption need not be true for the SKA precursors like the Murchison Widefield Array (MWA) where the effect of the primary beam needs to be accounted for while computing the expected signal [47].
Till date, people do not have clear understanding about either the DM particle physics model which governs the production of initial stage e ± spectrum, or the astrophysical parameters like galactic diffusion coefficient (D 0 ), magnetic field (B) etc. which are responsible for the creation of radio synchrotron flux. Thus in this analysis, we have constrained the DM parameter space, responsible for detecting the radio signal at SKA, for particular choices of astrophysical parameters which are on the conservative side for a typical dSph like Draco [44]. On the other hand, assuming some simplified DM model scenarios with trans-TeV DM masses, we have estimated the limits on the B − D 0 plane, required for the radio signal from Draco to be observed at SKA. Figure 11 shows the minimum σv required for any m χ in four different annihilation channels (bb, τ + τ − , W + W − , tt) for detection of DM annihilation induced radio signal from Draco with 100 hours of observation at SKA assuming a typical band width 300 M Hz. The radio signal has been taken to be detectable when the observed flux in any frequency bin rises three times above the noise so as to ensure that the detection is statistically significant and is not affected by spurious noise features. We vary the DM mass over a wide range of 10 GeV to 50 T eV , assuming 100% branching fraction (B f ) in one annihilation channel at a time. The predictions here are for a conservative choice of the diffusion coefficient (D 0 = 3 × 10 28 cm 2 s −1 ) [44]. The bands are due to the variation of the galactic magnetic field B from 1 µG (lower part of the band) to a more conservative value 0.1 µG (upper part of the band). As expected, the minimum σv required will be larger for lower magnetic fields as lower B reduces the the radio synchrotron frequency distribution (S ν ). These limits are the detection thresholds for SKA to observe radio signal from Draco. For m χ ∼ 1 T eV this limit could be as low as σv ∼ 3×10 −29 cm 3 s −1 . For lower values of the diffusion coefficient such as D 0 = 3 × 10 26 cm 2 s −1 [30,31,52], the detectability threshold band comes down even further. Moreover, we have not considered any halo substructure contributions which are expected to enhance the radio flux [43] or lower the threshold limits even more. Along with these lower limits, we have also shown the modelindependent upper limits (in 95% C.L.) on σv in various channels from cosmic-ray (CR) antiproton observation (dashed curves) [28] and from 6 years of Fermi LAT (FL) γ-ray data (dotted curves) [27]. Note that the upper bounds from cosmic-ray antiproton observations are the strongest ones. For each annihilation channel, the area bounded by the upper and lower limits represents the region in the σv −m χ plane which can be probed or constrained by SKA with 100 hours of observation. It is clear from the figure that even with conservative choices of astrophysical parameters for m χ ≈ 50 T eV , there are significantly large regions of the parameter space which can be probed in SKA. In such extreme cases, it is of course necessary to have a dark sector that allows high co-annihilation rates, so that the observed relic density bound is not excluded.
In the context of model-dependent analysis, three benchmark points, named as Model A1a, B2a and E, from our earlier work [37] have been considered here. These benchmarks correspond to the minimal supersymmetric standard model (MSSM) scenario where the lightest neutralino (χ 0 1 ) is the DM candidate (χ). Table I contains the possible annihilation channels with branching fractions, DM masses (m χ 0 1 ) and annihilation rates ( σv ) calculated in those benchmark points. All these quantities have been calculated using publicly available package micrOMEGAs [49,62]. Masses of neutralino and all other supersymmetric particles in these three cases are in the trans-TeV range. All of these benchmarks produce relic densities within the expected observational limits [1,63,64] and satisfy constraints coming from direct DM searches [2,3], collider study [65], lightest neutral Higgs mass measurements [66] and other experiments [67,68]. These benchmarks have been discussed in further de-tailed in reference [37]. Figure 12 shows the detectability of these benchmarks for Draco dSph in SKA after 100 hours of observations. The predicted flux in each case falls within the area which is still allowed by the data and is at least 3 times above the observation threshold of SKA. The upper limits here are taken from cosmic-ray antiproton observation at 95% C.L [28] and the lower limits have been calculated for D 0 = 3 × 10 28 cm 2 s −1 and B = 1 µG. Along with these we have indicated the total σv (listed in third column of table I) by different points for each benchmark. It should be noted that each such benchmark allows more than one annihilation channel. Consistently, the cross-sections are obtained by taking a sum over different channels appropriately weighted by the branching ratios: It is clear that all these high mass cases, which are allowed by cosmic-ray antiproton data, can easily be probed in SKA with 100 hours of observations. This is due to the following reason: • σv is more effective in offsetting 1 m 2 χ suppression (present in the expression for number density of DM pair; equation 1) partially, though not fully.
• A greater abundance of high-energy e ± is created by high m χ . This, via electromagnetic energy loss driven by the term proportional to b(E) in equation 3, generates a bigger flux of radio synchrotron emission, as evinced in the expression for J Synch (ν, r) (equation 8).
• The cases where one predicts more intense radio flux for higher m χ have DM annihilation mostly in the bb channels, as against the τ + τ − channel.
The corresponding cascade branching ratios as well as the three-body decay matrix elements and their energy integration limits are responsible for bigger radio flux. 4 Figure 11 shows how the detectability threshold at the SKA varies with the galactic magnetic field B, for a fixed value of D 0 . However, the values of these two parameters for the dSphs are quite uncertain. In view of this, it is also important to find out which regions in the astrophysical parameter space are within the scope of the SKA with 100 hours of observation, when both B and D 0 vary over substantial ranges. The detectability of the signal in the B − D 0 space is shown in Figure 13 for some illustrative DM scenarios. These correspond to DM masses of 1 (red curve) and 5 (blue curve) T eV .  ) and its pair annihilation rate ( σv ) inside a dSph along with branching fractions (B f in percentage) in different annihilation channels for the selected MSSM benchmark points from [37].
Cases where the dominant annihilation channel is bb are shown in the left panel, while the curves on the right panel capture the corresponding situations with τ + τ − as the main channel. The solid line of each color corresponds to the maximum value of σv for the chosen m χ , con-sistent with the cosmic-ray upper limit. All points above and on the left of the curve correspond to the combinations of B and D 0 that make the radio signals detectable over 100 hours of observations with the SKA. Points for higher B in this region correspond to models that are   [37] in the σv − m χ plane. The upper bars represent the 95% C.L. upper limits on σv corresponding to these benchmark points from cosmic-ray (CR) antiproton observation [28]. The lower bars show the minimum σv required for those benchmark points for the observation of radio flux from Draco dSph at SKA with 100 hrs of observations. The diffusion coefficient and the magnetic field have been assumed as, D 0 = 3 × 10 28 cm 2 s −1 and B = 1 µG respectively. detectable even with lower values of σv . Similarly, the detectable models having progressively lower σv are arrived at, as one moves to lower D 0 for a fixed value of B. It is evident from both the bb and τ + τ − channels (chosen for illustrations in Figure 13) that consistent values of σv can lead to detectability with 100 hours at SKA, for B 5 × 10 −4 µG for D 0 ∼ 10 27 cm 2 s −1 . On the other hand, for larger values of D 0 ∼ 10 30 cm 2 s −1 , we need B 10 −2 µG for the signal to be detectable.
Finally, we show in Figure 14 the final radio fluxes in two of the MSSM benchmarks listed in table I. For illustration we have taken the benchmark points A1a (m χ ∼ 1 T eV ) and E (m χ ∼ 8.5 T eV ). The annihilation is bbdominated for both cases. The yellow band here represents the SKA sensitivity [33], with a bandwidth of 300 M Hz. The band is due to the variation of the observation time from 10 hours (upper part of the band) to 100 hours (lower part of the band) [33,37]. The choice of the diffusion coefficient and the magnetic field is on the conservative side (D 0 = 3 × 10 28 cm 2 s −1 and B = 1 µG). We have already shown the detectability of these benchmarks in Figure 12. From this figure we can see that, even for 10 hours of observation, these high-m χ benchmarks can easily be observed in SKA over most of its frequency range.
One important feature of Figure 14 is that, in a wide range of frequency (300 M Hz -50 GHz), suitable for SKA, S ν for very high mass case (Model E) is higher than that for low mass case (Model A1a). Though the DM mass in Model E is larger, annihilation rate in this case is higher (∼ 9 × 10 −26 cm 3 s −1 ) than in Model A1a (∼ 3 × 10 −27 cm 3 s −1 ). In general, σv should have a 1 m 2 χ suppression, because of the energies of the colliding DM particles [57]. In spite of that, Model E (higher m χ ) has greater σv 5 than Model A1a (lower m χ ) mainly due to the closer proximity to a s-channel resonance mediated by CP-odd pseudoscalar in the annihilation process like χ 0 1 χ 0 1 → bb. For detailed information reader is referred to see reference [37] As discussed earlier, The higher value of σv partially offsets the effect of larger DM mass and consequently we get σv m 2 χ (which appears in the source function; equation 1) for Model A1a and Model E as 2.7 and 1.26 (in units of 10 −33 GeV −2 cm 3 s −1 ), respectively, though the latter model has larger m χ (∼ 8.5 T eV ) compared to the former (∼ 1 T eV ). Thus a 72-fold suppression due to m 2 χ results in a suppression just by a factor of ≈ 2 at the level of σv m 2 χ . Also, note that the bb annihilation channel dominates for both models. These observations, together with the discussion in section IV and the contents of Figure 10, explains a higher mass DM particle generating higher radio flux for scenarios where the bb annihilation channel dominates over τ + τ − .

VI. CONCLUSION
The aim of this paper has been to show that radio signal arising from a high mass (trans-TeV) WIMP DM can be detectable as radio synchrotron flux from a dSph, to be recorded by the upcoming SKA Telescope. We have analysed not only the particle physics aspects of DM annihilation and subsequent cascades leading to e ± pairs, but have also included the relevant astrophysical processes the electrons/positrons pass through before emitting radio waves, upon acceleration by the galactic magnetic field. We have set out to identify the mechanism whereby a trans-TeV DM candidate can thus be visible in radio search. We found that in the SKA frequency range, enhancement of the radio flux for this case is possible mainly due to the following reasons: • Larger cross section or annihilation rate (required to maintain relic density under the observed limit for a trans-TeV DM) facilitated by the dynamics of the particle physics model. This helps in compensating the 1 m 2 χ suppression due to large DM mass.
• Presence of energetic e ± in the DM annihilation spectrum in greater abundance. This partially re- suppression effect and on the other hand enhances, through the energy loss term, the electron/positron density at low energies which helps to produce large radio flux.
• Dominance of the annihilation channel bb which yields a comparatively larger abundance of e ± in all of the energy range of the spectrum produced by DM annihilation.
Simultaneously, effect of various astrophysical parameters (e.g. D 0 , b(E), B) on the radio synchrotron flux produced from the annihilation of a trans-TeV DM particle has been studied in detail.
Using SKA sensitivity we have drawn the limits in the σv − m χ plane to observe radio flux from Draco with 100 hours of observation. We found that, these limits are much more stronger than the previously obtained bounds on σv from Fermi-LAT γ-ray [27] or AMS-02 cosmicray antiproton observation [28]. Even for a conservative choice of astrophysical parameters (D 0 = 3 × 10 28 cm 2 s −1 and B = 1 µG), we found that these limits can go as low as σv ∼ 3 × 10 −29 cm 3 s −1 for a m χ ∼ 1 T eV . This indicates towards a large region of WIMP parameter space which can be probed through the upcoming SKA. Along with these we have also shown the limits in the B − D 0 plane. We found that, for a DM mass in the trans-TeV range, magnetic field as low as B ∼ 10 −3 µG and diffusion coefficient as high as D 0 ∼ 10 30 cm 2 s −1 are well enough to produce radio flux above SKA sensitivity in 100 hours of observation time.
Taking the minimal supersymmetric standard model (MSSM) as an illustration, we have shown that benchmark points with lightest neutralino masses (m χ 0 1 ) ∼ 1 -8 T eV , which satisfy all the constraints from observed relic density, direct DM searches and collider searches, can lead to detectable signals at the SKA with 100 hours observation, even with conservative choice of B and D 0 . We have illustrated how the effects mentioned earlier can lead to a larger radio flux for a high mass DM benchmark point compared to a low mass case, thus establishing the credibility of the search for heavy DM through radio observation in SKA 6 .