Possible explanation of indirect gamma ray signatures from hidden sector fermionic dark matter

We propose the existence of a hidden or dark sector besides the standard model (SM) of particle physics, whose members (both fermionic and bosonic) obey a local SU(2)H\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$_\mathrm{H}$$\end{document} gauge symmetry, while behaving like a singlet under the SM gauge group. However, the fermiomic fields of the dark sector also possess another global U(1)H\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$_\mathrm{H}$$\end{document} symmetry, which remains unbroken. The local SU(2)H\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$_\mathrm{H}$$\end{document} invariance of the dark sector is broken spontaneously when a scalar field in this sector acquires a vacuum expectation value (VEV), thereby generating masses to the dark gauge bosons and dark fermions charged under the SU(2)H\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$_\mathrm{H}$$\end{document}. The lightest fermion in this dark SU(2)H\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$_\mathrm{H}$$\end{document} sector can be a potential dark matter candidate. We first examine the viability of the model and constrain the model parameter space by theoretical constraints such as vacuum stability and by the experimental constraints such as Planck limit on relic density, LHC data, limits on spin-independent scattering cross section from dark matter direct search experiments etc. We then investigate the gamma rays from the pair annihilation of the proposed dark matter candidate at the Galactic Centre region. We also extend our calculations of gamma rays flux for the case of dwarf galaxies and compare the signatures of gamma rays obtained from these astrophysical sites.


Introduction
Experimental observations led by WMAP [1] and Planck [2] satellites reveal that only about 4 % of our Universe is made up of ordinary baryonic matter and about 26.5 % of it is constituted by the unknown nonluminous matter or dark matter. Dark matter (DM) is supposed to have a very weak interaction with the visible sector of the Universe and strong evidence of its presence come from gravitational probes only. Howa e-mail: amit.duttabanik@saha.ac.in b e-mail: debasish.majumdar@saha.ac.in c e-mail: anirban.biswas@saha.ac.in ever, the particle nature of DM and the reason for its high longevity (stability) are still unexplained. Direct searches of DM are also performed by various DM direct search experiments namely CDMS [3][4][5][6], CoGent [7], Xenon100 [8], LUX [9] etc. Although no convincing detection of DM has yet been reported by these direct search experiments, until present, LUX experiment provides most stringent bounds on DM spin-independent elastic scattering cross sections σ SI of the dark matter-nucleon with respect to its mass.
Dark matter particles can get trapped inside massive astrophysical bodies (due to their enormous gravity) like the Galactic Centre (GC), the solar core. This may happen if the velocities of DM particles fall below their escape velocities inside those massive bodies. When accumulated in considerable amount, these trapped DM particles may undergo pair annihilation and produce fermion-antifermion pairs, γ -rays etc. These γ -rays or fermions (positrons, neutrinos, antiprotons etc.), if found to be emitted in excess amount from these sites which cannot be explained by known astrophysical processes, may possibly be due to dark matter annihilation inside these sites. The detection of such excess γ -rays, positrons, neutrinos etc. thus provide valuable indirect signatures of the particle nature of dark matter. Besides the GC, the dwarf galaxies may also be rich in dark matter. The dwarf galaxies are a class of faint and small satellite galaxies of our Milky Way galaxy. The huge amount of DM content within these dwarf spherical galaxies (dsphs) is inferred from their mass to luminosity ratio ( M L ). The M L ratio for these galaxies are found to be much higher than what is expected from the estimation of their visible mass. The dark matter rich dsphs can also emit excess γ -rays due to the pair annihilation of dark matter. Nine such dwarf galaxies have recently been discovered in addition to the previously discovered 15 dwarf satellite galaxies of Milky way. Among these satellite galaxies, the γ -ray flux obtained from the Reticulum 2 (Ret2 in brief) dwarf galaxy shows an excess of γ -rays in the γ -ray energy range of 2-10 GeV [10] while the null results from other satellite galaxies provide stringent upper limits on the annihilation cross sections of DM particles for its various possible annihilation channels [11,12].
In the last few years, the analyses of Fermi-LAT publicly available data [13] by several groups [14][15][16][17][18][19][20] have confirmed the existence of a low energy (few GeV range) γray excess, which appears to be emerging from the regions close to the centre of our Milky way galaxy. The analysis of Fermi-LAT data [13] by Daylan et al. [21] shows that the γ -ray excess from the GC can be well explained by the annihilation dark matter scenario. They have also excluded all the known astrophysical processes which can act as the possible origin of this phenomenon. In Ref. [21] it is shown that the observed γ -ray spectrum from the GC can be well fitted by an annihilating dark matter particle having mass in the range ∼30−40 GeV which annihilates significantly into bb final state with an annihilation cross section σ v bb ∼ (1.4−2.0) × 10 −26 cm 3 /s with local dark matter density ρ = 0.3 GeV/cm 3 . In this work the authors have taken an angular region of 5 • around the centre of our galaxy as their region of interest (ROI) and used Navarro, Frenk and White (NFW) halo profile with γ = 1.26 for the computation of γ -ray flux. However, more recently, the authors Calore, Cholis and Weniger (CCW) of Ref. [22] have claimed to perform a detailed analysis of Fermi-LAT data along with all the possible systematic uncertainties using 60 galactic diffusion excess (GDE) models. The results obtained from the analysis of CCW provides a best fit for DM annihilation into bb final state having mass 49 +6. 4 −5.4 GeV with σ v bb = 1.76 +0. 28 −0.27 × 10 −26 cm 3 s −1 . Moreover, the analysis of CCW assumes a NFW profile with γ = 1.2 and a different region of interest (ROI) with galactic latitude |l| ≤ 20 • and longitude |b| ≤ 20 • masking out the inner region corresponding to |b| ≤ 2 • . Different particle physics models for dark matter that are simple extensions of Standard Model (SM) are proposed to account for this 1-3 GeV excess in γrays from GC . It is to be noted that apart from dark matter, non-DM sources such as millisecond pulsars may provide a feasible explanation to the excess of γ -ray observed at GC [51]. Study of unresolved point sources near GC by Lee et al. [52] suggests that point sources also contribute significantly to the gamma ray excess. However, in this work, we will consider DM as the origin of the observed excess in GC gamma ray to explore the phenomenology of dark matter. Also the study of γ -rays from previously known 15 different dwarf galaxies by Fermi-LAT [11] and eight newly discovered dwarf galaxies by Fermi-LAT with DES collaboration [12] give bound on DM mass and corresponding σ v for different annihilation channels.
In this work, we consider a "hidden sector" framework of dark matter without pretending any such discrete symmetry associated with it. We propose the existence of a hidden sector which has SU(2) H gauge structure. Dark fermions in this hidden sector are charged under this SU(2) H gauge group while all the SM particles behave like a singlet. Hence, the SM sector is decoupled from the dark sector and could interact only through the exchange of scalar bosons that exist in both sectors. Gauge bosons charged under SU(2) H are heavy and decay into dark fermions. Thus, the lightest one among dark fermions is stable and can be treated as a viable DM candidate. We show that the DM candidate in the present model that satisfy the limits from vacuum stability, LHC constraints, relic density, direct detection experiments can duely explain the Galactic centre γ -ray excess and also is in agreement with the limits on DM annihilation cross section obtained from the study of dwarf galaxies.
The paper is organised as follows: In Sect. 2 we present the hidden sector SU(2) H model. Constraints and bound on the model parameter space from vacuum stability, LHC results on SM Higgs, relic density etc. are described in Sect. 3. In Sect. 4 we show the results obtained for the available model parameter space and study of indirect searches of γ -ray is performed. Finally in Sect. 5 we summarise the work and make concluding remarks.

The model
We consider the existence of a "dark sector" that governs the particle candidate of dark matter. Just as the "visible sector" related to the known fundamental particles successfully explained by the Standard Model, we propose the existence of a hidden "dark sector" that relates the dark matter particles. We also presume that the Lagrangian of this hidden sector remains invariant under the transformations of a local SU (2) Therefore, the most general Lagrangian of the present proposed model contains the following gauge invariant terms: 1 In order to cancel the Witten anomaly [106] we need at least two (even numbers) of left handed fermionic SU(2) H doublets in our model. 2 Although, in order to keep similarity with the expression of the Standard Model Higgs doublet H , we have introduced the notation of three scalar fields, in the expression of , as G + 2 , φ • and G • 2 , however, the symbols + and 0 appearing in the superscript of dark sector scalar fields do not represent the electric charge of the corresponding scalar field as electric charge itself is not defined in the dark sector, which is invariant only under SU(2) H . with the covariant derivatives of the SU(2) L × U(1) Y doublet H and the SU(2) H doublets , χ iL , respectively, whilẽ = iσ 2 with σ 2 is the Pauli spin matrix. Moreover, g, g and g H are the respective gauge couplings corresponding to the gauge groups SU(2) L , U(1) Y and SU(2) H . In the above equation (Eq. (2)) F μν is the field strength tensor for the gauge fields A i μ (i =1-3) of the SU(2) H gauge group while H is the usual SM Higgs doublet. The global U(1) H invariance of the dark sector Lagrangian forbids the presence of any Majorana type mass terms of the fermionic fields ( f i , i = 1, 4) in Eq. (2). We have assumed at the beginning that the dark sector fermions are charged under a global U(1) H symmetry. Therefore invariance of the dark sector Lagrangian (Eq. (2)) under this U(1) H symmetry requires equal and opposite U(1) H charges between each fermion and its antiparticle. Thus we can say that there is some conserved quantum number in the theory which can differentiate between a fermion and its antiparticle. In other words it can be stated that the dark sector fermions in the present theory are Dirac type fermions. We have also assumed that the dark sector fermions ( f i , i = 1, 4) are in "mass basis" or "physical basis" so that the Lagrangian (Eq. (2)) does not contain any mixing term between these fermionic states. 3 The dark sector fermions can interact among themselves by exchanging dark gauge bosons A i μ and due to the presence of these interaction modes all the heavier fermions such as f i (i = 2-4) decay into the lightest one ( f 1 ). Consequently, the lightest fermion f 1 is stable and can be a viable dark matter candidate. Like the hidden sector gauge fields A iμ , the dark matter candidate f 1 also gets mass when the postulated SU(2) H symmetry of the hidden sector breaks spontaneously by the VEV of . Thus, the expression of mass of the fermionic dark matter candidate can easily be obtained using Eq. (2); it is We have already mentioned before that, due to the presence of the gauge invariant term λ 3 H † H † , the neutral components of both scalar doublets, namely h 0 and φ 0 , possess mass mixing between themselves. The mass squared mixing matrix between these two real scalar fields are given by After diagonalising the mass squared matrix M 2 scalar , we obtain two physical eigenstates h 1 and h 2 , which are related to the old basis sates h 0 and φ 0 by an orthogonal transformation matrix O(α) where α is the mixing angle between the resulting physical scalars. The relation between the physical scalars h 1 and h 2 with the scalar fields h 0 and φ 0 are given as The expressions of the mixing angle α and the masses of the physical real scalars h 1 and h 2 are given by We assume the physical scalar h 1 is the SM-like Higgs boson which has been observed by the ATLAS and the CMS detector [107,108]. Therefore we have adopted the mass (m 1 ) of h 1 and VEV v of h 0 to be ∼125.5 and 246 GeV, respectively. Thus, we have three unknown model parameters which control the interactions of the dark matter candidate f 1 in the early Universe, namely the mixing angle α, the mass (m 2 ) of the extra physical scalar boson h 2 and more importantly, the mass m f 1 of the dark matter particle f 1 . In the rest of our work we have computed the allowed ranges of these model parameters using various theoretical, experimental as well as observational results. Throughout the work, for simplicity we take mass of fermionic DM candidate ( f 1 ) to be m.

Constraints
In this section we will discuss various constraints and bounds on model parameters that arise from both theoretical aspects and the experimental observations.
• Vacuum stability To ensure the stability of the vacuum, the scalar potential for the model must remain bounded from below. The quartic terms of the scalar potential is given as where H is the SM Higgs doublet and is the hidden sector Higgs doublet. Conditions for the vacuum stability in this framework is given as • LHC Phenomenology In the present model of hidden sector (SU(2) H ) fermionic dark matter discussed earlier in Sect. 2, an extra Higgs doublet is added to the SM. This dark SU(2) H Higgs doublet provides an additional Higgs-like scalar that mixes up with the SM Higgs. Large Hadron Collider (LHC) performing the search of Higgs particle (ATLAS and CMS collaboration) have already discovered a Higgs-like particle having mass about 125 GeV. The excess in γ γ channel reported independently by ATLAS [107] and CMS [108] confirmed the existence of Higgs-like bosons. In the case of Hidden sector SU(2) H model, the mixing between SM Higgs with Dark Higgs results in two Higgs-like scalars. In the present scenario we take one of the scalar (h 1 ) as the SM Higgs with mass m 1 = 125 GeV. We further assume that the signal strength of scalar h 1 also satisfies the limits on the same obtained for the newly discovered boson. Thus, h 1 in the present framework is identical with the SM-like Higgs as reported by LHC Higgs search experiments (ATLAS and CMS). The signal strength of Higgs boson (h), decaying into a particular final state (x x, x is any SM particle), is defined as where σ ( pp → h) and Br(h → xx) are the Higgs production cross section and its branching ratio of any particular decay mode (x = quark, lepton or gauge boson), obtained from LHC experiments. The corresponding quantities computed using Standard Model of electroweak interaction are denoted by σ SM ( pp → h) and Br SM (h → x x), respectively. For the present model, the signal strength of the SM-like scalar h 1 is then defined by where the quantities in the numerator of Eq. (11) are the production cross section and the branching ratio of SM-like Higgs boson h 1 , which are computed using the present formalism. Now due to the mixing of scalar bosons, the coupling of SM-like Higgs boson to the SM fermions and gauge bosons are modified with respect to the SM Higgs boson (h) by the cosine of mixing angle α, whereas the couplings of non-SM scalar boson h 2 to SM particles are multiplied by a factor sin α. Hence the ratio σ (pp→h 1 ) σ SM ( pp→h) = cos 2 α and from a similar argument one can see that σ (pp→h 2 ) σ SM ( pp→h) = sin 2 α. The SM branching ratio can be expressed as Br SM is the decay width of SM Higgs boson h into any final state particles and SM is the total SM Higgs decay width having mass m 1 = 125 GeV. Similarly one can derive the expression for branching ratio of h 1 into any specific decay channel in the present model Br( is the decay width of h 1 decaying into x x final state while 1 is the total decay width of h 1 in the present model. Hence, the signal strength of h 1 in Eq. (11) can be written in the form where we have denoted cos α as c α . It is to be noted that apart from the decay into SM particles the SM-like scalar h 1 can also have invisible decay mode into dark matter particles. Therefore the total decay width of h 1 , in the present model, can be written In Eq. (13), inv 1 is the invisible decay width h 1 for the channel h 1 → f 1f1 . For m 1 > 2 m the expression of invisible decay width of h 1 is given by since the coupling between h 1 and the dark matter candidate is proportional to m v s s α . In the above, v s is the VEV of SU(2) H Higgs doublet and s α = sin α. Similarly for the other scalar involved in our model, the signal strength R 2 is expressed as where SM (m h = m 2 ) is the total decay width of SM Higgs boson if it has mass m h = m 2 while 2 is the total decay width for the non-SM scalar boson h 2 The coupling between dark matter and h 2 depends on the factor m v s cos α . Hence, the invisible decay width of h 2 ( inv 2 ) for m 2 > 2m is given as As stated earlier, we consider h 1 with mass m 1 = 125 GeV to be the Higgs-like scalar and infer R 1 > 0. 8 [109] and invisible decay branching ratio Br 1 inv ≤ 0.2 [110] where Br 1 inv = 1 inv / 1 is defined as the ratio of invisible decay width to the total decay width.
• Dark matter relic density The DM relic density as measured by Planck satellite experiment is given as [2] DM h 2 = 0.1199 ± 0.0027.
In Eq. (19), h is the Hubble parameter measured in units of 100 km s −1 Mpc −1 . We calculate the relic density for the fermionic (SU(2) H ) dark matter candidate in the assumed dark sector in our model by solving the Boltzmann equation. The relic density of the DM candidate is obtained by solving the Boltzmann equation [111] dn dt where n is the number density of DM particle and n eq is the same in equilibrium. In Eq. (20), σ v is the thermal where g * is the effective number of d.o.f. (degrees of freedom), M Pl is the Planck mass (∼1.22 × 10 19 GeV) and x f = m/T f with T f being the freeze out temperature of the DM species, respectively. We compute the freeze out temperature T f by solving iteratively the following equation: In order to obtain the freeze out temperature of DM and hence its relic density using Eqs. (21) and (22) we need to calculate the thermal average of the product between total DM annihilation cross section (σ ) and the relative velocity (v) of two annihilating DM particles. The expression for the thermally averaged DM annihilation cross section into all possible final states is given as where the factors K i , (i = 1, 2) are the modified Bessel functions and √ s is the centre of mass energy. In the present formalism dark matter candidate f 1 can annihilate into the SM particles through s-channel processes mediated by the scalar bosons h 1 and h 2 . In Eq. (23), σ (s) denotes the total annihilation cross section of dark matter into all possible final states which are allowed by the Lagrangian given in Eq. (2). Feynman diagrams for different annihilation channels of f 1 are shown in Fig. 1. The expressions of σ v for different final state annihilation of dark matter into SM particles are derived from the Feynman diagrams shown in Fig. 1. The value of σ v obtained for DM annihilation into SM fermion and antifermion pairs ( ff ) in the final state is of the form where m is the DM mass and m f is the mass of specific fermion ( f = quark or lepton). In Eq. (24) v and v s are the vacuum expectation values of SM Higgs doublet and dark Higgs doublet, N c is the colour quantum number (3 for quarks and 1 for leptons). Moreover, 1 , 2 are the total decay widths of the scalar bosons h 1 , h 2 and the expressions of 1 and 2 are given in Eqs. (13) and (17). We also calculate σ v for W + W − and Z Z channels which proceed through the s-channel exchange of scalar bosons h 1 , h 2 (see Fig. 1). The expressions of σ v W + W − and σ v Z Z are furnished here: In the above, m W and m Z denotes the respective masses of W and Z bosons. Annihilations of DM particles into scalar bosons h 1 and h 2 are also taken into account. The process of DM annihilation into scalars h 1 or h 2 is also scalar mediated; it depends on scalar couplings between h 1 and h 2 . The s-channel annihilation cross sections of f 1 annihilating into the pairs of h 1 and h 2 , calculated using the f 1f1 → h i h i , i = 1, 2 annihilation diagram, take the following form: where λ i jk is the coupling for the vertex involving three scalar fields h i h j h k . The expressions for the scalar couplings λ 111 , λ 211 , λ 122 and λ 222 are given in Appendix A. We calculate the thermally averaged annihilation cross section of the present DM candidate using Eqs. (23)- (28). We then compute the freeze out temperature T f by solving Eq. (22) and finally obtain the relic density of f 1 at the present epoch from Eq. (21).
• DM direct detection Direct detection of DM particle is based on the scattering of the DM particle with the target nucleus of the detector material. Fermionic dark matter in the present model can undergo elastic scattering with the detector nucleus. This elastic scattering of the DM and the nucleus will transfer a recoil energy to the target nucleus which is then calibrated. From the nonobservance of such elastic scattering events the direct detection experiments give the upper bound of elastic scattering cross sections for different possible masses of dark matter. The scattering cross section is expressed as cross section per nucleon for enabling direct comparison of the results from different experiments. In the present model DM fermion of mass m can interact with the tar-get nucleus through t-channel Higgs mediated processes through both h 1 and h 2 . The spin-independent (SI) elastic scattering cross section off the detector material normalised to cross section per nucleon can be written [66] σ SI = sin 2 2α 4π where m r = mm p m+m p is the reduced mass for the DMnucleon system and λ p [66] is given in terms of the form factors f q and the proton mass m p as Using Eqs. (29)- (30), we calculate the spin-independent elastic scattering cross section of the DM fermion off the nucleon and compare it with the experimental bounds from LUX [9]. Note that both the DM annihilation cross section and the DM-nucleon scattering cross section depend on an effective coupling g eff = | m v s s α c α | (Eqs. (24)- (29). This effective coupling is a useful parameter to explain the dark matter phenomenology in the present framework. Further discussions on the effective coupling are given in Sect. 4.
• DM indirect detection The existence of DM has now been well established from gravitational evidence in astrophysical scale. Indirect search of DM focuses on the non-gravitational search of DM candidate and explores the particle physics nature of DM. The astrophysical sites such as the Galactic Centre (GC), dwarf galaxies etc. are of great interest, since dark matter can be trapped and accumulate at GC due to the enormous gravity in the region of GC and the mass to luminosity ratio of dwarf galaxies indicates the presence of dark matter in large magnitude. These sites are suitable for indirect search of DM as DM particles trapped in these regions can undergo annihilation into various SM particles which can further produce gamma rays, neutrinos etc. Thus any observed excess in the fluxes of γ -rays, positrons, anti-protons from such sites can indicate DM annihilation processes in those sites if other astrophysical phenomena cannot explain the observed excess. Fermi-LAT [112] searches for the excess emission of γ -rays originating from GC and dwarf galaxies. Observation of the excess in e + /e − and p/p flux is performed by the AMS-02 [113] experiment. In this section we will study Fermi-LAT observed gamma ray flux results from the centre of the Milky Way and surrounding dwarf spheroidal galaxies (dSphs).
The expression for the differential γ -ray flux obtained from a region of interest (ROI) subtends a solid angle d centered at GC is given as where σ v f is the average thermal annihilation cross section of DM particles annihilating into final state particle f and dN f dE γ is the photon energy spectrum of DM annihilation into the same. The factor J appearing in Eq. (31) is related to the quantity of dark matter present at the astrophysical site considered and is expressed in terms of dark matter density as In Eq. (32) the line of sight (los) integral is performed over an angle θ , is the angular aperture between the line connecting GC to the Earth and the direction of line of sight. In the above Eq. (32), r = r 2 + s 2 − 2r s cos θ where r = 8.5 kpc, is the distance to the Sun from GC. It is clear from the expression of Eq. (32) that the value of the J factor is dependent on the nature of the chosen ρ(r ) factor, i.e., the DM halo density profile ρ(r ). In the present work, we consider the Navarro-Frenk-White (NFW) [114] halo profile. The DM density distribution for the NFW halo profile is given as where r s = 20 kpc is the characteristic distance and ρ 0 is normalised to the local DM density, i.e., ρ = 0.4 GeV cm −3 at a distance r from GC. The analysis by Daylan et al. [21] of Fermi-LAT data suggests an excess in γ -ray in the γ energy range of 2-3 GeV at GC. The same analysis demonstrates that this excess can be explained by the annihilation of 31-40 GeV DM into bb with σ v bb = 1.4-2.0 × 10 −26 cm 3 s −1 . In this work [21], the inner galaxy gamma ray flux (5 0 from GC) is calibrated using a NFW halo profile with γ = 1.26 and a local DM density ρ = 0.3 GeV cm −3 . In a recent work by Calore, Cholis and Weniger (CCW) [22] a detailed analysis is performed for the GC γ -rays along with the systematic uncertainties using 60 galactic diffusion excess (GDE) models. Results from CCW analysis provides a best fit for DM annihilation into bb having mass 49 +6.4 −5.4 GeV with σ v bb = 1.76 +0.28 −0.27 × 10 −26 cm 3 s −1 . However, CCW analysis of Galactic Centre excess (GCE) for gamma ray have also considered generalised NFW profile (γ = 1.2, ρ = 0.4 GeV cm −3 ) for a different region of interest (ROI) with galactic latitude |l| ≤ 20 • and longitude |b| ≤ 20 • masking out inner |b| ≤ 2 • . In another work Agrawal et al. [115] reported that annihilation of heavier dark matter (up to 165 GeV for bb channel) can also explain the observed GCE in γ -ray when the uncertainties in DM halo profile (NFW) and the J -factor are taken into account. However, in the present work, we do not consider any such uncertainties in halo profiles or J values and we use the canonical NFW halo profile used in the CCW analysis. Using Eqs. (31)-(33), we calculate the γ -ray flux (in GeV cm −2 s sr −1 ) for the ROI described in CCW analysis for Fermi-LAT data. As mentioned earlier we consider for our calculations the NFW profile with γ = 1.2 and ρ = 0.4 GeV cm −3 . Apart from the GC region, dwarf galaxies of the Milky Way are also of great significance for indirect search of DM as these galaxies are supposed to be rich in dark matter. Recent analyses of γ -ray fluxes from 15 Milky-Way dSphs reported by Fermi-LAT [11] provide a limit on the DM mass and corresponding thermally averaged annihilation cross section σ v f into different channels f (τ and b). Fermi-LAT have used their 6 year data collected by the Fermi Large area Telescope and performed an analysis for 15 dSphs using "pass-8 event level analysis" (see [11] and references therein). In another work [12] Fermi-LAT in collaboration along with Dark Energy Survey (DES) collaboration also provide similar bound on σ v f where they include data for 8 new dSphs. For both the analyses presented in [11,12] a canonical NFW halo profile (γ = 1) is considered, and the astrophysical J factors are measured over a solid angle = 2.4 × 10 −3 sr with angular radius 0.5 • . Independent searches carried out by Fermi-LAT [11] and DES-Fermi-LAT collaboration on 15 previously discovered and eight recently discovered different dSphs reported no significant excess in observed γ -ray. Results from the DES dSphs [12] also predict an upper bound to the observed γ -ray energy flux with 95 % confidence limit (C.L.) for eight newly found dSphs. Gamma ray flux for dwarf galaxies when integrated for an energy range extending over a region of solid angle is expressed as where dN dE γ is the γ -ray. The expression of the flux presented in Eq.  r (s, θ))ds, (35) calculated over a solid angle = 2.4 × 10 −3 sr subtended by the ROI (0.5 • angular radius) for the NFW halo profile (γ = 1). The density distribution function for the NFW profile with γ = 1 is then where r s is the NFW scale radius and ρ 0 represents the characteristic density for the dSphs. In the case of a Fermi-LAT analysis, the J factors for different dSphs are adopted from Ref. [11]. We use the values of the J factor from [12] for computing the gamma ray flux for eight DES dSphs for the dark matter candidates in our model. However, it is to be noted that J factors for DES dSphs candidates are obtained assuming a point-like dSphs instead of having a spatial extension (as in the case of [11]) to avoid the uncertainties in halo profile arising from spatial extension. Calculation of gamma ray flux is also based on the assumption that the spectrum dN dE γ follows the conventional power law dN dE γ ∼ 1 E 2 . As mentioned earlier, study of 15 dSphs by Fermi-LAT and eight other dSphs by DES-Fermi-LAT collaboration found no significant excess in γ -ray from these dwarf galaxies. However, a recent search on a newly discovered dwarf galaxy Reticulum 2 (Ret2) in a work by Geringer-Sameth et al. [10] has reported an excess in observed γ -ray signal. In the present work, we calculate the γ -ray flux for annihilation of hidden SU(2) H fermionic dark matter into γ -ray through different SM final states and explore whether the model can account for GCE in γ -ray and also satisfies the bounds on gamma ray flux from dwarf satellite galaxies. As mentioned earlier, in the present model dark matter candidate ( f 1 ) is fermionic in nature and it interacts with the visible world (SM particles) through the exchange of two real scalar bosons h 1 and h 2 . As a result the annihilation cross sections of the DM dark candidate f 1 into the final states that composed of SM particles (mainly light quarks and leptons) are proportional to the square of relative velocity (v 2 ) between the annihilating dark matter particles (p wave process). Now the averaged DM relative velocity is proportional to ∼ √ 3/x [116,117] with x = m T is a dimensionless quantity and T being the temperature of the Universe. Hence, in our model, the thermally averaged annihilation cross section used for computing DM relic density, at x ∼ 20 − −30, is different from the annihilation cross section (for x ∼ 3 × 10 6 [116,117]) needed to calculate the γ -ray flux at the Galactic Centre and dwarf galaxies. The latter quantity is velocity suppressed as the average DM relative velocity is ∼10 −3 when the annihilation of DM occurs at the GC. Among all the annihilation channels of f 1 , the annihilation mode f 1f1 → bb plays a significant role for the γ -ray excess observed from GC and dwarfs satellite galaxies as it is the most dominant annihilation channel for the considered mass range of DM. In order to explain the GC gamma-excess by DM annihilation to bb, the annihilation cross section should be ∼1.76 +0. 28 −0.27 × 10 −26 cm 3 /s [22]. Although in the present case, the thermally averaged annihilation cross section for the bb annihilation is quite small, however, the quantity σ v bb can be significantly enhanced using the Breit-Wigner resonant enhancement mechanism [116,117]. Breit-Wigner enhancement occurs only when the mass of the dark matter (m) is nearly equal to half of the mediator mass (in our case it is the mass of h 2 ). Therefore, we have defined the mass of the hidden sector scalar boson (h 2 ) and the centre of mass energy √ s in the following way: where δ < 0 represents the physical pole and z is the measure of excess centre of momentum energy scaled by 4m 2 . In terms of z, Eq. (23) for the bb annihilation channel, can now be written with the expression of σ (z) bb is given by 4 and where γ 2 = 2 m 2 , 2 being the total decay width of h 2 of mass m 2 . It is to be noted that the upper limit 4 Since the Breit-Wigner enhancement occurs when m ∼ m 2 /2, as a result only the term proportional to of the above integration should be ∞ (see Eq. (23)), however, the integrand becomes negligibly small when z approaches z e f f ∼ max [4/x, 2|δ|] for δ < 0 [37,117].
Using the above prescription, we calculate the thermally averaged annihilation cross section σ v bb of the dark matter candidate f 1 for GC and dwarf spheroidal galaxies. The actual values of σ v bb , γ 2 and δ for the two chosen bench mark points (BP1, BP2) are given in Table 1 of Sect. 4. We have found that for |δ| ∼ 10 −3 the annihilation cross section σ v bb ∼ 1.9 × 10 −26 cm 3 /s which can explain the excess of gamma ray flux in GC. 5

Calculational procedures and results
In this section we present the computation of dark matter annihilation cross sections as also the DM-nucleon elastic scattering cross sections. They are required for the calculation of relic densities and the comparison of the latest DM scattering cross section bound given by the LUX direct detection experiment. The invisible decay widths and signal strengths for the SM-like scalar is also calculated in order to constrain the model parameter space. The gamma ray flux is then computed within the framework of SU(2) H fermionic dark matter for galactic centre as also for dwarf galaxies and the results are compared with the experimental analysis.

Constraining the model parameter space
The fermionic dark matter in the present model can annihilate through scalar mediated (h 1 and h 2 ) s-channel processes. As mentioned in Sect. 3, the model parameter space is first constrained by the vacuum stability conditions given in Eq. (9). The signal strengths R 1 and R 2 for the Higgs doublets h 1 (SM) and h 2 (dark sector) are then computed using Eqs. (19)- (23). With the chosen constraints on R 1 (R 1 ≥ 0.8, Ref. [109]) the invisible decay branching ratio of SM-like Higgs Br 1 inv is calculated and the parameter space is further constrained by the LHC experiment limit of Br 1 inv (Br 1 inv ≤ 0.2 [70]). The parameter space thus constrained is then used to compute the thermal averaged annihilation cross section σ v of the present fermionic dark matter candidate and the relic density is obtained by solving the Boltzmann equation (using Eqs. (20)-(23)). The annihilation cross sections are computed with the calculated analytical formulae given in Eqs. (24)-(28) with two choices of VEV for (dark Higgs doublet) namely v s = 246 and 500 GeV. In our calculation we consider the mass m 1 of the SM-like Higgs boson h 1 to be 125 GeV. The calculation is performed for two values of the dark sector scalar h 2 masses and they are m 2 = 100 and 110 GeV. These relic densities are compared with the dark matter relic density given by Planck [2]. Thus the Planck result further constrains the parameter space of our model. With this available parameter space we evaluate the dark matter-nucleon spin-independent scattering cross section (σ SI ) for the purpose of comparing our results with those given by the dark matter direct detection experiments such as LUX, XENON100 etc. In this way we restrict our model parameter space by different experimental results.
In Fig. 2a and b we show the calculated values of σ SI with different DM mass in the present model where the conditions from vacuum stability, bound on SM Higgs signal strength and DM relic density results from Planck have been imposed. We first choose certain values of m 1 and m 2 and vary the couplings λ i , i = 1-3 (satisfying vacuum stability conditions given in Eq. (9) for two different values of v s , which also constrain the mixing angle α through the Eq. (7).
Here we want to mention that we have varied λ 1 and λ 2 in the range 0-0.2 with the values of both λ's being evenly spread within the considered range. Consequently the value of the parameter λ 3 becomes fixed by the vacuum stability criteria given in Eq. (9), which is also varied with an equal interval in the range |λ 3 | < 2 √ λ 1 λ 2 . The model parameter space thus obtained is then further constrained by imposing the conditions R 1 > 0.8 and Br 1 inv < 0.2 from the LHC results. Using this restricted model parameter space satisfying both vacuum stability and LHC bounds, we therefore calculate the relic density of the dark matter candidate f 1 by solving the Boltzmann equation (Eq. (20)) for different values of DM mass. Finally, we consider specific range of model parameter space which is in agreement with DM relic density reported by Planck experiment and for these parameter space we compute the spin-independent direct detection cross section using Eqs. (29)- (30). In this way the viable model parameter space for the dark matter candidate f 1 is obtained. Figure 2a is for the case m 2 = 100 GeV, while Fig. 2b is for the case m 2 = 110 GeV. The upper limits on σ SI for different values of DM mass, obtained from the LUX DM direct search experiment, are also shown in Fig. 2a, b by the blue line for comparison. The red and green scattered regions as shown in Fig. 2a, b correspond to the two choices of v s = 246 and 500 GeV, respectively. From Fig. 2a it can be observed that only the region near the resonances of scalar bosons h 1 and h 2 is in agreement with the upper limit on σ SI predicted by LUX. It is also seen from Fig. 2a that the choice of v s does not alter the allowed range of parameter space. Observation of Fig. 2b shows that, apart from the SM Higgs resonance region (m ∼ m 1 /2) there exists another allowed range of the m-σ SI parameter space in the vicinity of the non-SM scalar resonance (m ∼ m 2 /2). Note that variation of m with σ SI depicted in Fig. 2b depends only on the masses of scalar bosons and does not suffer any significant change due to change in v s . The non-SM Higgs signal strength R 2 (cal- Higgs and relic abundance of DM obtained from Planck [2]. Limits on DM-nucleon cross section from LUX [9] is also plotted (blue line) for comparison culated using Eq. (16)) for the valid m-σ SI parameter space shown in Fig. 2a, b is very small and R 2 < 0.2.
In this work, we assumed two values for the VEV v s (246 and 500 GeV) for the hidden sector Higgs doublet HS . From Eq. (6), we observe that the mixing between the scalars h 1 and h 2 depends on the VEV of HS and H . Hence, the choice of v s may change the range of available model parameter space. In Fig. 3a, b, we plot the variation of the Higgs mixing angle α between h 1 and h 2 with λ 3 for m 2 = 100 and 110 GeV with m 1 = 125 GeV (mass of SM-like Higgs). Needless to mention that the regions of the α-λ 3 space shown in Fig. 3a, b are consistent with the bounds from vacuum stability, SM Higgs signal strength from LHC, relic abundance of DM from Planck and limits on DM-nucleon scattering cross section from the LUX direct DM search experiment. Plots in Fig. 3 are produced using a similar method; we have applied previously to obtain viable model parameter space for Fig. 2. However, the α-λ 3 plane in Fig. 3 is further constrained by imposing LUX DM direct detection bound. The plots in Fig. 3a are for the case when m 1 = 125 GeV and m 2 = 100 GeV while plots in Fig. 3b represent the allowed α − λ 3 parameter space when m 2 = 110 GeV for the fixed value of m 1 = 125 GeV. The green and blue regions in Fig. 3a, b correspond to two different values for the VEV of the dark Higgs doublet, v s = 246 and v s = 500 GeV, respectively. From Fig. 3a (m 2 = 100 GeV case) one observes that for both the considered values of VEV v s , the mixing parameter λ 3 remains small and is confined within the region |λ 3 | < 0.01. For the case when v s = 246 GeV (the red region of Fig. 3a), the limit of mixing angle α ranges between −0.1 and 0.1. However, these range (of mixing angle) varies within the limit |α| ≤ 0.2 when v s = 500 GeV is chosen (green region shown in Fig. 3a). Study of the λ 3 -α plots in Fig. 3b (plotted for m 2 = 110 GeV) reveals that for both values of v s considered in Fig. 3, the mixing parameter is small (|λ 3 | < 0.01). The mixing angle α is bounded in the range |α| < 0.15 and |α| ≤ 0.30 for v s = 246 and 500 GeV, respectively. 4.2 Calculation of gamma ray signals from galactic centre and dwarf galaxies In this section, we calculate the γ -ray flux from the galactic centre and dwarf galaxies for the fermionic dark matter in the framework of the present model and compare our results with the experimental observations. For these calculations we consider two benchmark points (BPs) from the restricted parameter space that satisfy both theoretical and the experimental bounds (mainly vacuum stability, LHC constraints on the SM Higgs signal, Planck results for relic abundance and a direct detection limit on m-σ SI from LUX) for two choices of the h 2 mass, mainly, m 2 = 100 and 200 GeV. In Table 1 we tabulate the chosen BPs along with the model parameters. There are two chosen sets of benchmark points in Table 1 and we denote them BP1 and BP2. The GC gamma ray flux is calculated using Eqs. (31)- (33) for the BPs tabulated in Table 1. The annihilation cross section σ v bb for the dark matter particle is calculated using the Breit-Wigner enhancement technique using Eqs. (37)-(40) discussed in   . 4 Comparison of the GC γ -ray flux data from [22] with those calculated for benchmark points in Table 1 Sect. 3. The gamma ray spectrum dN dE in Eq. (31) is obtained from Ref. [118] for the annihilation of DM into any specific channel. The gamma ray spectra for BP1 and BP2 are then calculated for the specified region of interest adopted from Ref. [22] (|l| ≤ 20 • , 2 • ≤ |b| ≤ 20 • ) using the NFW halo profile (with γ = 1.2, ρ = 0.4 GeV cm −3 ). In Fig. 4, we show the calculated GC gamma ray flux (in GeV cm −2 sr −1 ) for our proposed DM candidate with BP1 and BP2. We also show in Fig. 4 the CCW data for comparison. Green and blue lines in Fig. 4 represent the calculated γ -ray spectra for BP1 and BP2, respectively. Both benchmarks points are in agreement with the findings from GC gamma ray study presented in CCW [22]. From Fig. 4 it can be observed that flux calculated using the set BP1 (m = 50 GeV) is in better agreement with the findings from CCW analysis.
We now further investigate how well the DM candidate in our model can explain the observed extragalactic γ -ray signatures from various dwarf galaxies. From their six years observations on 15 dwarf galaxies, the Fermi-LAT experiment did not obtain any significant excess of γ -rays. The Fermi-LAT collaboration [11], however, in a recent work provides combined bound on DM mass and thermally averaged DM annihilation cross section into SM particles for these 15 dSphs. A similar bound in the m-σ v f plane is also presented recently in another work [12] for eight new dSphs jointly by the Fermi-LAT and the DES collaborations. In this work we calculate the thermally averaged annihilation cross section of DM annihilating into SM sector in our model and compare them with experimental results given by [11,12]. The allowed range of m-σ v space (for annihilation into bb) along with the bounds on σ v (into bb channel only) obtained from GC γ ray search results CCW [22] and dwarf galaxies [11,12] compared with the same obtained from benchmark points in Table 1 In Fig. 5, we plot the bounds on DM annihilation cross section σ v bb (for the annihilation channel DMDM→ bb) with dark matter mass m obtained from galactic [22] and extragalactic [11,12] γ -ray search experiments. We calculate the variations of the same plotted in Fig. 5 for the benchmark points BP1 (for m 2 = 100 GeV) and BP2 (for m 2 = 110 GeV) considered in our model.
Black contours shown in Fig. 5 are the 1σ , 2σ and 3σ contours given by the CCW [22] analysis of GC gamma ray excess observations. The blue line in Fig. 5 describes the bounds in the m-σ v bb plane given by the analysis of gamma rays from previously discovered 15 dSphs and they are adopted from [11]. Also shown in Fig. 5, the yellow band which is the 95 % confidence limit (C.L.) region adopted from the analysis in Ref. [11] for DM annihilation into bb. The combined bounds on σ v bb for different DM mass m from a recent study of the newly discovered eight DES dwarf galaxies [12] are given by the pink coloured line in Fig. 5. The green horizontal line in Fig. 5 shows the annihilation cross section for thermal dark matter that may yield the right DM relic abundance obtained from the Planck experiment.
From Fig. 5 one readily observes that the calculated values of σ v bb for the benchmark points BP1 and BP2 in our model broadly agrees with the 1σ , 2σ and 3σ allowed regions in m-σ v bb plane obtained from the experimental results. This can also be noted from Fig. 5 that these benchmark points are consistent the combined limit from DES dwarf satellite data and falls within the 95 % C.L. limit predicted by Fermi-LAT for 15 dSphs. Also the calculated values of σ v bb for the benchmark points considered in our work lie below the upper bound on thermal DM annihilation cross section. Hence, DM fermion in the present model can account for the galactic centre excess in γ -ray and is also consistent with the bounds on gamma ray flux from Milky-Way dwarf satellite galaxies.
We now calculate the gamma ray flux for eight new dwarf satellite galaxies discovered by the DES experiment for the hidden sector fermionic dark matter candidate proposed in this work. These calculations are performed with each of the benchmark parameter sets BP1 and BP2 given in Table 1. The Gamma ray flux for each of these eight dSphs in the work [12] is computed using Eq. (34) and the values of the J factors (Eq. (35)) for each of the eight dSphs adopted from Ref. [12]. In Ref. [12] these J factors are estimated by integrating the dark matter density (adopting the NFW halo profile for the DM density distribution) along the line of sight over a solid angle = 2.4 × 10 −4 sr −1 . As previously mentioned the gamma ray spectrum dN dE is also obtained from Ref. [118] for this calculation. The calculated flux for each of the eight dSphs are shown in eight plots (a-h) of Fig. 6. Also shown in each of the eight plots of Fig. 6, the respective upper bounds of the flux given by the experimental observations of gamma rays from each of the eight dSphs. These are shown as red coloured points, while the computed flux in this work for the respective dSphs are given by continuous lines in Fig. 6. The green and blue continuous lines in each of the plots (a-h) of Fig. 6 correspond to the calculated flux using the benchmark points BP1 and BP2, respectively. It is clear from Fig. 6 that the fluxes calculated, assuming the annihilation of the DM candidate in our proposed model, for all the eight dSphs do not exceed the upper limit of γ flux set by the experimental observations of DES collaboration.
Besides the 15 dwarf galaxies investigated earlier and the eight other recently explored dwarf galaxies, one more dwarf galaxy, namely Reticulum 2 (Ret2), has been probed very recently. Geringer-Sameth et al. [10], after an analysis of observed gamma rays from Ret2 dwarf galaxy reported an excess of gamma ray emission from Ret2. From their analysis of the Ret2 data Geringer-Sameth et al. provide different C.L. allowed contours in the m-J 19 σ v −26 plane where m is the mass of the dark matter and J 19 σ v −26 is the product of the J factor in the units of 10 19 GeV 2 cm −5 and thermal averaged product σ v of annihilation cross section and relative velocity in the units of 10 −26 cm 3 s −1 for various final state SM channels. As mentioned earlier in this work the DM candidate primarily annihilates into bb; only the contours for the DM pair annihilation into bb channel are adopted. For the present dark matter model with the constrained parameter space discussed earlier we compute the quantity J 19 σ v −26 for different dark matter mass m annihilating into bb channel. However, the value of the J factor for Ret2 has been adopted from [10]. In their work Geringer-Sameth et al. [10] estimated the J values by performing line of sight integral over a circular region with angular radius 0.5 • surrounding the dwarf and  Table 1 over a solid angle  [10] annihilating into bb pair) are presented by green coloured patches in both plots of Fig. 7. From Fig. 7, it can easily be observed that J 19 σ v −26 in the present model calculated for DM annihilating into bb channel (for benchmark points with m 2 = 100 and 110 GeV) is within the 3σ C.L. limit. Hence fermionic DM candidate in the present framework can also explain the observed excess in γ -ray from Ret2.

Discussions and conclusions
In this work, we have proposed the existence of a hidden sector which obeys a local SU(2) H and a global U(1) H gauge symmetries. In order to introduce fermions which are charged under this SU(2) H gauge group one should have at least two fermion doublets in order to avoid the "Witten anomaly". The particle and the antiparticle of these dark fermions are different as they possess equal and opposite U(1) H charges. Similar to the usual Higgs doublet in the visible sector, this hidden sector also has an SU(2) H scalar doublet , which, however, does not have any U(1) H charge. The SU(2) H gauge symmetry breaks spontaneously when the neutral component of the scalar doublet gets a VEV and thereby generates masses of all the dark gauge bosons ( A μ ) and dark fermions ( f i ). Since the dark sector fermions interact among themselves through the dark gauge bosons, therefore all the heavier fermions as well as the dark gauge bosons can decay into the lightest fermion and hence the lightest fermion in this dark sector can be treated as a particle for the viable dark matter candidate. In fact, in this model this lightest fermion is the only dark matter candidate. The dark fermions and dark gauge bosons do not mix with the SM fermions and gauge bosons due to the non-abelian nature of two SU(2) groups. However, the dark sector scalar field can interact with the SM Higgs-like scalar in the visible sector and only through this interaction two sectors are mutually connected.
We therefore test the viability of the present model by using theoretical and experimental constraints on the relevant model parameters, such as vacuum stability conditions, bounds on relic abundance of DM from the Planck experiment, direct detection limits on DM-nucleon scattering cross section from LUX experiment. LHC bounds on signal strength and invisible decay width of the SM Higgs, are also used to constrain the parameter space. From such analyses we find that only a small region of the parameter space near the scalar resonances (when m m 1 2 and m 2 2 ), is consistent with the current experimental bounds. Study of the model parameters, thus constrained, shows that the mixing between the two scalars (h 1 , h 2 ) of the model is very small (mixing angle α ≤ 0.3 deg) and depends on the VEV of the dark scalar doublet ( ). With the allowed regions of parameter space thus obtained, for the present DM candidate (dark fermion f 1 ) we compute the gamma ray flux from the GC region. While calculating the gamma ray flux from GC we have used the Breit-Wigner enhancement mechanism for the computation of DM annihilation cross section into the bb final state ( f 1f1 → bb). These computational results are then compared with the experimental analyses of the Fermi-LAT GC gamma ray flux data considering that the dark matter at the GC primarily pair annihilates into the bb channel. Our proposed DM candidate can indeed explain the results from these experimental analyses.
In search of indirect evidence of dark matter from astrophysical sources, the gamma rays from various dwarf satellite galaxies are also explored for a possible signature of excess gamma rays from these sites. To this end 15 such dwarf galaxies have earlier been investigated and more recently the gamma ray observation has also been reported from eight more newly discovered dSphs. From the analyses of these observational results different C.L. bounds have been given in the parameter space of the σ v bb -m plane. We compare our computational results with these experimental bounds and found that the γ -rays that the DM candidate in our model produce on pair annihilation can simultaneously satisfy the observational results from GC and dwarf galaxies. We also demonstrate that the calculated fluxes in our model for each of the recently discovered eight dwarf galaxies lie below the corresponding upper limits of the fluxes obtained from the observational results of these dwarf galaxies. We further demonstrate that our calculations are also in good agreement with the analysis of Ret2 dwarf galaxy observations.
Our work clearly demonstrates that the dark matter candidate proposed in this work is a viable one to explain the γ -rays from both the GC region and dwarf galaxies simultaneously. However, the dark matter can also pair annihilate into fermion-antifermion pairs and there are experiments such as AMS-02 that look for the excess of e + /e − or pp in cosmos. In a recent work, the AMS-02 collaboration have reported their first measurement of p/p flux [119]. A modelindependent analysis of this AMS-02 p/p data is performed by Jin et al. [120]. In this work [120], the upper limits in σ v value for DM annihilation into SM particles (quarks and gauge bosons) for different considered DM halo profiles (NFW, Isothermal, Moore) are obtained. The analysis presented in the work [120] also considered four different propagation models, namely the conventional, MED, MIN and MAX models. 6 We have also checked that the DM in our model satisfies upper bound on σ v bb given in Ref. [120] when NFW profile is considered. This is found to be true for both the cases of dark sector scalar mass m 2 = 100 GeV and m 2 = 110 GeV. Hence, fermionic dark matter explored in the present model can serve as a potential candidate for dark matter. Upcoming results from LHC and also DM direct and indirect search experiments may provide stringent limits on the available model parameter space.