Gamma rays from Dark Matter Annihilation in Three-loop Radiative Neutrino Mass Generation Models

We present the Sommerfeld enhanced Dark Matter (DM) annihilation into gamma ray for a class of three-loop radiative neutrino mass models with large electroweak multiplets where the DM mass is in O(TeV) range. We show that in this model, the DM annihilation rate becomes more prominent for larger multiplets and it is already within the reach of currently operating Imaging Atmospheric Cherenkov telescopes (IACTs), High Energy Stereoscopic System (H.E.S.S.). Furthermore, Cherenkov Telescope Array (CTA), which will begin operating in 2030, will improve this sensitivity by a factor of $\mathcal{O}{(10)}$ and may exclude a large portion of parameter space of this radiative neutrino mass model with larger electroweak multiplet. This implies that the only viable option is the model with lowest electroweak multiplets i.e. singlets of $SU(2)_{L}$ where the DM annihilation rate is not Sommerfeld enhanced and hence it is not yet constrained by the indirect detection limits from H.E.S.S. or future CTA.


I. INTRODUCTION
We are yet to identify the mass and particle nature of the Dark Matter (DM) of the universe despite having extensive astrophysical and cosmological observations supporting its existence.Recently the DM, having quantum numbers under the Standard Model (SM) gauge group SU (2) L × U (1) Y and with mass in the TeV range, has come under the focus of Imaging Atmospheric Cherenkov Telescopes (IACT) as gamma ray produced in the DM annihilation at the central region of the Milky Way galaxy is within the detection reach of IACTs [1][2][3][4].
The flux of the gamma ray photons from cosmic sources rapidly falls when their energies reach E > ∼ O(TeV).For this reason, the satellite which has the detection area of sub-m 2 , is not very sensitive for the very high energetic (VHE) gamma ray of the energy range, 100 GeV -100 TeV.On the other hand, such high energy gamma ray interacts with the upper region of the atmosphere and creates a shower of very energetic secondary charged particles which reaches at about 10 km height.These particles move faster than the speed of light in the air and therefore emit the faint blue Cherenkov light.This Cherenkov light is beamed around the direction of the incident primary photon and it illuminates the ground of about 50000 m 2 , which is often referred as Cherenkov light pool.Only 100 photons per m 2 on the ground can be seen for a primary photon of TeV energy.So if a telescope is somewhere within the light pool and has large mirror area to collect enough photons, it will observe the air shower.Therefore, the effective detection area of a Cherenkov telescope is approximately given by the area of the light pool which is much larger than that of a satellite.The prominent IACTs like MAGIC [5], VERITAS [6], CANGAROO [7] and H.E.S.S. [8] have revealed intriguing astrophysical VHE gamma ray sources of our universe.In [9,10] H.E.S.S. collaboration presented search result for gamma signal coming from DM annihilation in the inner region of Milky Way and put upper limits on the annihilation cross section σv for the DM in the TeV mass range that is not within the reach of collider searches or direct detection experiments.
In addition, the non-relativistic (NR) DM which has O(TeV) mass and electroweak charge, receives the nonperturbative Sommerfeld enhancement  which increases the DM annihilation rate into gauge bosons, i.e DM DM → W W, ZZ, γγ, γZ ,significantly.As a consequence, this class of NR DM with TeV mass range, has better prospects of being detected in currently operating H.E.S.S or future CTA [35,36] which will begin operation in about 2030.
Apart from the DM nature of the universe, the origin and smallness of the neutrino mass is yet to be concluded.The Krauss-Nasri-Trodden (KNT) model [37] ties these two issues together by radiatively generating neutrino mass 1 at three-loop with DM particle running in the loop.The Beyond Standard Model (BSM) content of the model consists of two single charged singlet scalars, S + 1 , S + 2 and three singlet RH neutrinos, N Ri , i = 1, 2, 3 under SM gauge group with masses lie in the GeV-TeV range.Here, the lightest singlet RH neutrino N R1 plays the role of DM.Subsequently, KNT model can be generalized [39] by replacing S + 2 with Φ having integer isospin and hypercharge, Y = 1 and N Ri with F i that has integer isospin and Y = 0 under SM gauge group.
In the generalized KNT model, the lightest neutral fermion component, F 0 1 is the viable DM candidate.Such replacement in KNT model with large electroweak multiplets have been studied for triplet [40], 5-plet [41] and 7plet [42] cases.In [43], we have investigated the charged lepton flavor violating processes in the generalized KNT model.In this work we focus on 5-plet and 7-plet cases because the Z 2 symmetry, {S + 2 , N Ri } → {−S + 2 , −N Ri } needed to prevent the Dirac neutrino mass term in the Lagrangian, is not required anymore for larger multiplets like in 5-plet and 7-plet cases.
The article is organized as follows.In section II, we present the generalized KNT model.Section III describes the formalism to calculate Sommerfeld enhanced DM annihilation processes in generalized KNT model.In section IV, we present the relic densities of the DM candidate, including the correction due to SE, in 5-plet and 7-plet cases via thermal freeze-out process and also via non-thermal out-of-equilibrium decay.The DM annihilation cross sections into electroweak bosons at the galactic center are presented in section V. Finally we conclude in section VI.

II. THE MODEL
Apart from the SM field content, we add the following BSM fields in the generalized KNT model which are charged under SM gauge group, SU (3) c × SU (2) L × U (1) Y as Complex scalars: S + 1 ∼ (0, 0, 1), Φ ∼ (0, j φ , 1), and Real fermions: where j φ and j F are integer isospin of SU (2) L .
In this comparative study, we focus on two set of models in this class; 5-plet model: The SM Lagrangian is augmented in the following way, where, c denotes the charge conjugation and dot sign, in shorthand, refers to appropriate SU (2) contractions.Also L α and e Rα are the LH lepton doublet and RH charged leptons respectively and Greek alphabet α stands for generation index.Moreover, [F ] αβ = f αβ and [G] iα = g iα are 3 × 3 complex antisymmetric and general complex matrices respectively.Finally, H denotes the SM Higgs doublet.The scalar potential is given by, The three-loop neutrino mass generation and the DM stability depend on the V 3 term of Eq.(3).Explicitly the 1 For a review of radiative neutrino mass generation models, please see [38].
relevant terms of V 3 for 5-plet and 7-plet models are, Here the λ term in Eq.( 4) is not invariant under Z 2 and eventually induce the decay of F 0 1 where the width is Γ DM ∼ λ 2 .But, as pointed out in [41], the bound on DM mean life-time sets λ to be very small, and in the limit when λ → 0, the Z 2 symmetry emerges.On the other hand, the λ term is absent in Eq.(5) because j φ ⊗ j φ contains symmetric (antisymmetric) irreducible representation with same isospin T φ for even (odd) integer isospin which is further contracted with Φ † to obtain a singlet.So, for two identical scalar multiplets, the antisymmetric combination is zero and hence no λ term for j φ = 3.
As pointed out in [43], the mass splittings among component fields of the scalar multiplet is controlled by λ Hφ2 (Φ † .H).(H † .Φ) ⊂ V 2 term after electroweak symmetry breaking and allowed splittings only lead to ∆m 2 ij /M 2 0 ∼ 10 −3 for invariant mass of the scalar multiplet, M 0 = 10 TeV and the ratio becomes smaller for M 0 > 10 TeV.On the other hand, the mass splittings in fermionic component fields are zero at tree-level and only receive O(100) MeV splittings due to radiative correction after electroweak symmetry breaking.Therefore such scenario is can be considered as near-degenerate case.

III. SOMMERFELD ENHANCED DARK MATTER ANNIHILATION A. DM Candidate
In the generalized KNT model, the lightest neutral component of the fermion multiplet, F 0 1 is the viable DM candidate.In comparison, the neutral component of the scalar multiplet, φ 0 = 1 √ 2 (S + iA) could have provided S to be DM but it is ruled out as it induces Z-mediated dark matter nucleon scattering of the order 10 −39 cm 2 which is much larger than the exclusion limit set by the direct detection experiments [54].One can avoid this DM-nucleon scattering channel if the splitting between S and A is large enough to make this scattering kinematically forbidden but there is no renormalizable term in the Lagrangian which can induce such splitting in a generic way.Still, higher dimensional operator can split the S and A component [44] but then it is needed to address the UV completion of the model.Therefore, we restrict ourselves only to the renormalizable Lagrangian, and therefore the DM candidate is set to F 0 1 .

B. SE Annihilation Cross-sections
When the DM is non-relativistic, v DM c and m W,Z m DM , the exchange of massive W and Z gauge bosons between DM components will induce Yukawa potential and γ exchange will induce Coulomb potential which in turn significantly modifies the wavefunction of the incoming DM states and enhances the annihilation cross-sections.This phenomenon is known as Sommerfeld Enhancement (SE).The calculation of Sommerfeld enhanced DM annihilation cross section is well studied subject so here we follow the prescriptions given in [34,44].In the following we briefly review them to set up our notation.
As the Sommerfeld enhancement is considered for 2 → 2 processes, we first define 2-particle states which consist of incoming component fields of fermion multiplet.Sommerfeld enhancement takes place in DM (co)annihilation processes with final states W ± , Z and γ bosons so we only consider the 2-particle states which are CP-even and have total charges Q = 0, ±1, ±2.In the case of DM annihilation in the galaxy halo at present times, only 2-particle states with Q = 0 are relevant.Moreover, we have considered M F1 < M F2,3 , therefore only component fields of F 1 multiplet enter into 2-particle states.We define 2-particle state vector corresponding to F 1 as The ordering of the component 2-particle states, as above, is arbitrary.Therefore, one can choose any other ordering of states.The modification of the wavefunction is determined by solving the radial Schrodinger equation with effective potential, where r is the magnitude of the relative distance between two component fields in their center-of-mass frame, the kinetic energy of the incoming DM states, i.e.
The wavefunction Ψ jj ,ii gives the transition amplitude from |ii states to |jj states in the presence of effective potential, V .The double indices ii , jj and kk run over the states of the 2-particle state vector defined in Eq.( 6).
We primarily focus on the S-wave annihilation so we set l = 0 and have Here, is the momentum associated with the 2-particle state, |jj and d jj = m j + m j − 2m S denotes the mass differences between DM and other states of the multiplet.Q kk is the electric charge associated with state |kk .Also, α W = α and n W = 1 for W boson exchange and α Z = α/ cos 2 θ W and n Z = 1/ cos θ W for Z boson exchange.Finally, f jj ,kk is the group theoretical factor associated with SU (2).Now by using dimensionless variables defined as x = αm F1 r, φ = (m W /m F1 )/α, v = (v/c)/α and d ii = d ii /m F1 /α, we re-write the coupled radial Schrodinger equations as where the dimensionless momentum, k2 jj = 2 v − 2 d jj .At large x, Ψ jj ,ii behaves as Ψ jj ,ii ∼ T jj ,ii e i kjj x where T jj ,ii is the transition amplitude provided the effective potential is dominated by Yukawa potential.Now if the annihilation matrix for final state f is given by Γ (f ) jj ,ii , the annihilation cross section is, where c = 2 for |F 0 1 F 0 1 state as it consists of identical fields.

IV. DARK MATTER RELIC DENSITY
The relic density of DM in the universe is measured by Planck Collaboration as Ω DMh 2 = 0.1199 ± 0.0022 (68% C. L.) [45].The fermionic DM in the generalized KNT model can achieve this relic density either by the thermal freeze-out process or non-thermal process as we will describe below.

A. Thermal Freeze-out of DM
The thermal freeze-out of fermionic DM of the 5-plet and 7-plet, both proceed mainly through • gauge interactions in dominant S-wave and sub-dominant P wave channels as the DM is non-relativistic and they are controlled by gauge coupling g and receive non-negligible Sommerfeld enhancement in mainly S-wave annihilation cross-sections.
• yukawa interactions in sub-dominant P -wave channels which are controlled by g iα couplings and are less significant because of large gauge annihilation as pointed out in subsequent discussion.
For DM in TeV mass range, the thermal freeze-out may take place either in the broken or symmetric phase of the SM depending on its mass.The critical temperature where the cross-over between high-temperature symmetric phase and low-temperature broken phase of the SM takes place is T c = 159 ± 1 GeV [47].But the gauge singlet with mass at the electroweak scale ( > ∼ M W )and coupled with Higgs may change the dynamics of the phase transition.In that case, the phase transition can be first order and critical temperature can be T c ∼ 100 GeV [46].So considering the freeze-out temperature of DM as T F > ∼ T c ∼ 100 GeV, the freeze-out condition, z F = M DM /T F ∼ 20 − 30 implies that when M DM > ∼ 2 − 3 TeV, the DM freezes out in the symmetric phase of the SM.
When the DM freezes out in the broken phase, the calculation of the Sommerfeld enhanced cross-sections of (co)-annihilation processes which enter into the Boltzmann equation, requires solving the Schrodinger matrix equation with effective potential given in Eq. (9).But this calculation is greatly simplified if the freeze-out takes place in the symmetric phase because then SU (2) L is a good symmetry and therefore we can express the twoparticle states, which are direct product states, in terms of definite states of SU (2) irreducible representations.Moreover, as SU (2) gauge bosons are massless in this phase, the effective potential is the Coulomb potential and therefore the Sommerfeld enhanced cross-section is where σ 0 is the perturbative cross-section and S is given by, Here, β = v DM /c and α is the corresponding coupling of the SU (2) group.Moreover, Y = 0 for F 1 so no need to consider U (1) Y contribution.Following [50][51][52], the effective potential is, potential for total iso-spin is given as, The DM abundance is calculated by solving the following Boltzmann equation [48,49], Here, z = M DM /T .H and s are the Hubble rate in the radiation dominated era and the entropy density respectively.They are given by, where M pl = (8πG N ) −1/2 = 2.44 × 10 18 GeV is the reduced Planck Mass.Also g * and g * s are the total and effective relativistic degrees of freedom respectively.In addition, Y = n/s where n is the number density and Y eq is value of Y in thermal equilibrium.And hence The thermal rate of 2 → 2 scattering that involves (co)annihilating component fields of F 1 into SM fields, ij → a b, at temperature T , is denoted by γ and given as, where K 1 (x) is the modified Bessel function of first kind and the reduced cross section, summed over all (co)annihilation channels, is σ(s) = 2sλ(1, M 2 DM /s, M 2 DM /s)σ(s).Also, λ(a, b, c) = (a − b − c) 2 − 4bc is the Kallen function.
At the freeze-out, the DM and other component fields are non-relativistic, therefore we can decompose the reduced cross section up to P wave contribution as follows, where β is the DM velocity (here, c = 1) in the center of mass frame given by β = 1 − 4M 2 DM /s.Now as V is in isospin-1 representation, the 2-particle state V V will be in 1 ⊗ 1 = 0 ⊕ 1 ⊕ 2 representations with definite isospin.Therefore we only have to consider 0, 1 and 2 definite representations coming from 2 ⊗ 2 and 3 ⊗ 3 in the case of 5-plet and 7-plet respectively.Therefore the potentials for definite representations in the case of 5-plet are, whereas, for 7-plet, they are, V Finally, we have the Sommerfeld enhanced S-wave coefficients for 5-plet and 7-plet, c s and c s respectively as follows, where S(x) T symbolically denotes the result of integration on S(x) in Eq.( 19) because it is an implicit function of s.Also, the first, second and third terms of above equations denote the 0, 2 and 1 definite representations respectively.On the other hand, the gauge contribution to P-wave coefficient, c p is given as Also there will be P wave contribution from yukawa terms, where the sum is taken over the charged lepton flavors.As the components of the fermion multiplets are exactly degenerate or almost-degenerate in the symmetric and broken phases respectively, all the (co)annihilation channels contribute equally.The set of parameters of generalized KNT model which is relevant for DM relic density calculation via thermal freeze-out is, {M F1 , M F2 , M F3 , M φ , g 1α } apart from the SM gauge couplings of scalar and fermion multiplets.In this analysis we scan over M F1 ∈ (1, 50) TeV, M F2,3 ∈ M F1 + (1, 10) TeV, M φ ∈ M F1 + (10, 100) TeV.Moreover, g 1α couplings are chosen so that they satisfy the neutrino constraints as described in [43].
From Fig. 1, we can see that the inclusion of Sommerfeld enhanced S-wave contribution significantly changes the mass of the DM for which the correct relic density can be obtained.Moreover, the dominant contribution in thermal freeze-out comes from the gauge contribution as it involves large S-wave and P-wave coefficients in comparison to the P-wave contribution coming from the Yukawa sector of generalized of KNT model, as shown in Fig. 2 for two temperature values.Apart from thermal freeze-out of DM which is mainly controlled by the gauge interactions as seen from Fig. 2, it is possible to set the DM relic density non-thermally by the out-of-equilibrium decay of φ scalar via φ + → F 0 1 e + R in generalized KNT model.But as both Φ and F 1 are charged under the gauge group, the processes φ i φ j ↔ V V and F 1i F 1j ↔ V V will keep them in the thermal equilibrium.Therefore, one important condition is that the temperature where the decay takes place must be smaller than the temperatures where the gauge reactions of Φ and F 1 decouple.
The decay temperature is given as, where the decay width Γ φ for φ + → F 0 1 e + α process is Moreover, the thermal rate of the decay at temperature T is given as where, n eq is the number density of φ + in equilibrium.Finally, the gauge reaction density, ij ↔ ab at temperature T is given by, where s min = Max{(M i + M j ) 2 , (m a + m b ) 2 } and σA is the reduced 2 → 2 scattering cross section of component fields i, j of the scalar and fermion multiplets, Φ and F 1 to SM fields a, b via gauge interactions.The decoupling conditions for both scattering and decay rates are set to be γ neqH ≤ 1.From Fig. 3, we can see that the gauge reaction densities of Φ component fields decouple when temperature becomes small.One the other hand, the decoupling of the inverse decay process F 0 1 e + R → φ + which would deplete the amount of F 0 1 , is necessary.This condition sets the corresponding decay width to be very small, at the order of ∼ 10 −18 GeV so that the inverse process remains decoupled throughout the whole thermal history of the universe as shown in the Fig. 4 (left).Besides, Fig. 4 (right) represents the δ − |g 1α | plane bounded by the constraint Γ φ < ∼ 10 −18 GeV so that the inverse decay process remains out of equilibrium during the thermal evolution of the universe.We can see that, such small decay width of φ + implies that the mass difference between φ + and F 1 0 needs to be of the order O(1 − 10) MeV and |g 1α | ∼ 10 −4 for M φ ∼ 10 TeV.Therefore from this estimates, we can infer that the out-ofequilibrium decay of φ + to generate DM content of the universe only holds for a fine-tuned parameters of the model.Nevertheless, one can extend the generalized KNT model with another sector which can non-thermally produce the DM without any fine-tuning.The detailed construction of such extended model and its connection to the thermal history of the universe are beyond the scope of this paper and left for future investigation.

V. RESULT AND DISCUSSION
In section 1, we have identified the masses of F 1 0 for which standard thermal freeze-out process gives the current relic density in both 5-plet and 7-plet cases.Moreover, if we relax the thermal freeze-out scenario and consider the non-thermal production of DM in the generalized KNT model, one can achieve the current DM density for a wide range of masses, 5−50 TeV but to achieve that, fine-tunings in mass difference |M φ −M F1 | and |g 1α are required.Still one can have extended KNT model with a dark sector which can assist the non-thermal DM production but will not be subject to any fine-tuning.Therefore, in the subsequent analysis of Sommerfeld enhanced DM annihilation, we have considered the range 1 − 50 TeV because this range has better sensitivity in IACTs.

A. DM Direct Detection
The DM candidate F 0 1 does not couple to quarks at tree-level because of its vanishing hypercharge.However at one loop level, due to exchange of W boson, it has effective coupling with the quarks which leads to both spin-dependent and spin-independent contribution in DM-nucleon scattering.The spin-dependent cross-section is suppressed by the mass of the DM which is at O(TeV).On the other hand, the spin-independent cross-section for fermionic multiplet with integer isospin j, which doesn't depend on the DM mass, is given by [49], where, M Nucl is the mass of the target nucleus, f parametrizes nucleon matrix element as n| q m q qq|n = f m n nn and from lattice result, f = 0.347131 [53].FIG.5: Spin-independent cross section of DM-nucleon interaction.The shaded region is excluded by XENON1T (2017) data [54].The exclusion limit on DM mass in [54] is given up to 10 TeV.Here we have extrapolated this exclusion limit up to 50 TeV.As we can see, the thermal DM scenario with MF 1 = 9.9 TeV for 5-plet and MF 1 = 22.85 TeV for 7-plet, are almost at the verge of exclusion by the XENON1T (2017).

B. Sommerfeld Enhanced Cross-sections
In this section we present the Sommerfeld enhanced cross-sections of F 0 1 F 0 1 → W + W − and F 0 1 F 0 1 → γγ which are sensitive to IACTs.Moreover, in the analysis we set, v DM = 10 −3 which is the scale of DM average velocity in the galactic halo.In Fig. 6, we can see the resonance and dips occurring for both 5-plet (red) and 7-plet (black) at particular mass values of DM due to SE in the presence of Yukawa potential induced by the exchange of massive W and Z bosons in the limit of non-relativistic velocity.Apart from the dips at 3.8 TeV (5-plet) and 4 TeV (7-plet), σv ww is larger than its tree-level value (blue and brown lines for 5-plet and 7-plet respectively) for almost all of the DM mass range, 1-50 TeV.In fact for this mass range, it is large enough to be almost excluded by the H.E.S.S. limit (orange dashed line) provided that the F 0 1 is the dominant DM of the universe and follows the Einasto density profile [9].In addition we can see from Fig. 6 that the future CTA will improve the exclusion limit by a factor of O(10) (purple dashed line) [36].
In addition, Fig. 7 represents the Sommerfeld enhanced cross section, σv γγ for the process, F 0 1 F 0 1 → γγ.At tree-level, this process does not take place because the DM is charge neutral but due to the multiple exchange of gauge bosons i.e W ± , Z and γ and charged states, F (±Q) 1 in the ladder diagrams, the effective coupling with the photons is possible when the DM is non-relativistic.If φ > ∼ 1 and/or v > ∼ 1, the Sommerfeld enhancement will be suppressed and in that case, F 0 1 F 0 1 → γγ proceeds through one-loop process that gives, for 1-50 TeV mass range, σv γγ of the order 10 −28 − 10 −31 cm 3 s −1 .Again, we can see from Fig. 7 that apart from some dips, for almost all of 1-20 TeV mass range, H.E.S.S. (orange dashed line) can exclude the DM in case of 5-plet (red line) and 7-plet (black line) using gamma-line searches.In passing, we would like to point out that, unlike the case of scalar DM with larger electroweak multiplet which is focused in [44], the mass splittings among the fermionic component fields do not suppress the Sommerfeld enhancement as they are nearly degenerate for O(TeV) mass range.

VI. CONCLUSION AND OUTLOOK
In this work, we have presented the Sommerfeld enhanced DM annihilation cross sections to gamma rays for a class of three loop radiative neutrino mass generation models with large electroweak fermion and scalar multiplets i.e 5-plet and 7-plet under SU (2) L .We have shown that larger multiplet leads to more enhanced annihilation cross section because more charged component fields of the fermion multiplet contribute to the ladder diagrams.As a consequence, the potential matrix element, V ii ,jj induced by W exchange depends roughly on the group theoretical factor, j(j + 1) − m(m + 1) where j is the total isospin and m is the T 3 value of the component field in the multiplet.Moreover, the annihilation matrix element into final states W W , Γ ii ,jj also depends on factor, j(j + 1) − m 2 .Therefore, larger multiplets will naturally give rise to larger cross sections that we have already observed by comparing σv ww and σv γγ for 5-plet and 7-plet in Fig. 6 and Fig. 7 respectively.
Finally, putting astrophysical uncertainties aside, because of the larger Sommerfeld enhancement for the higher electroweak multiplets, we can see from Fig. 6 and Fig. 7 that for almost all of mass range, 1-50 TeV, the constraints from H.E.S.S. can exclude the F 0 1 being the DM, provided it is the dominant DM component and follows the Einasto density profile in both cases of 5-plet and 7-plet.That leaves only the singlet and triplet case as a viable DM candidate in the generalized KNT model.The singlet fermion, N R1 of the KNT model which is electroweak neutral, does not receive any Sommerfeld enhancement.On the other hand, the triplet case, where the DM candidate is the neutral component of the fermion multiplet with isospin, j = 1, will have enhanced annihilation processes due to exchange of electroweak bosons but being smaller representation than the 5-plet or 7-plet, it may have larger potion of parameter space yet to be excluded by the H.E.S.S. limit but will be within the reach of future CTA sensitivity limits.The detailed analysis for gamma rays coming from DM annihilation in different astrophysical environments (galactic center and dwarf spheroidal galaxies) for the triplet case is left for future investigation [55].

FIG. 3 :
FIG.3: Decoupling of the gauge reaction densities of Φ with temperature.
Sommerfeld enhanced cross-section σv WW at the galactic halo of the Milky Way for 5-plet (red) and 7-plet (black).Moreover, σvww without SE is given for 5-plet (blue) and 7-plet (brown).Here orange and purple dashed lines are H.E.S.S observed limit and CTA sensitivity limit on σvww respectively