Multicomponent Scalar Dark Matter with an Extended Gauge Sector

We consider an extension of the Standard Model of particle physics with an additional $SU(2)$ gauge sector along with an additional scalar bidoublet and a non-linear sigma field. The neutral components of the bidoublet serve as dark matter candidates by virtue of the bidoublet being odd under a $Z_2$ symmetry. Generic beyond Standard Model constraints like vacuum stability, invisible decay of higgs, Higgs alignment limit and collider constraints on heavy gauge bosons restrict the parameter space of this model. In this multicomponent dark matter scenario, we investigate the interplay between the annihilation and co-annihilation channels originating from the new gauge sector as those contribute to the relic abundance. We also inspect the direct detection constraints on scattering cross-sections of the dark matter particles with the detector nucleons and present our observations.


I. INTRODUCTION
The Standard Model (SM) of particle physics has been experimentally proven to be a very successful theory over the last few decades.However, there are a few limitations of the SM which motivate physicists to look for Beyond the Standard Model (BSM) scenarios.One such very important drawback with the SM is that it can not explain the particulate nature of the Dark Matter (DM).DM consists of more than one fourth of the energy budget of our universe [1,2].In the quest of solving the DM problem, one very popular class of models are those that incorporate a Weakly Interacting Massive Particle (WIMP) in their spectrum [3][4][5].In a typical WIMP model, there is an extended sector in addition to the SM particles which is weakly coupled to the SM and protected by a stabilising symmetry.Various scalar extensions including singlet extensions [6][7][8][9][10][11][12], two-singlet extensions [13][14][15][16][17][18][19][20][21][22], inert 2HDM [23,24], 2HDM+scalar singlet [25,26], 2HDM with gauge extension [27], and also various gauge extensions mostly of different kinds of additional U (1) [28][29][30][31][32][33][34] with an extended scalar and/or fermion sector have been studied in this context.
In this work, we extend the gauge sector of SM with an additional SU (2) and also the scalar sector with an additional scalar bidoublet and a non-linear sigma field.The SM-like scalar doublet and the sigma field contribute to the spontaneous symmetry breaking whereas the scalar bidoublet provides the DM candidates protected by an imposed Z 2 symmetry.An extended gauge sector introduces co-annihilation channels for the relic calculation as well as modifies the direct detection constraints -the main goal of the present work is to study in this specific model context how these phenomenological considerations affect the parameter space of the model.
The paper is organised as follows: in Sec.II we introduce the model.In Sec.III we discuss the model constraints arising from vacuum stability, invisible decay of the Higgs particle, Higgs alignment, and collider searches for heavy gauge bosons.In Sec.IV, we investigate the DM phenomenology motivating the DM candidates and calculating the relic abundance and the direct detection constraints before finally presenting our conclusions in Sec.V.All the relevant couplings are organised in Appendix B after a few comments on perturbativity in Appendix A. A brief note on the loop level calculations of Higgs decay to two photons is provided in Appendix C with appropriate loops relevant to the model under study.The Boltzmann equations and Feynman diagrams relevant for the relic are given in Appendix D.

II. THE MODEL
We consider an SU (2) extension of the gauge sector of the SM.Furthermore, the scalar sector is extended with an additional scalar bi-doublet odd under a Z 2 symmetry and a non-linear sigma field (for details see Tab.I).The gauge sector can now be written as SU (2) 0 × SU (2) 1 × U (1) 2 (where the subscripts help us to keep track of different gauge
charges).We will see that three out of the four degrees of freedom of the scalar doublet and three degrees of freedom of the sigma field are eaten up by the six gauge bosons of the gauge sector, thus giving rise to massive gauge bosons, while the photon remains massless due to the residual U (1) EM symmetry.The electrically neutral components of the scalar bidoublet protected by a Z 2 symmetry furnish the dark matter candidates of the universe.The various scalar fields can be conveniently expressed in the following form After a few redefinitions, the latter two can be recast as where and similarly for the π's.The scalars transform under these gauge symmetries in the following manner: Thus the Φ 2 and Σ transform as bidoublets under the SU (2) 0 × SU (2) 1 whereas the Φ 1 is a doublet under SU (2) 1 .
The U i signify the unitary transformation matrices corresponding to the gauge symmetry indexed i and the covariant derivatives accordingly take the form With these definitions in place, the kinetic energy term of the scalar sector can be written as The SU (2) 0 × SU (2) 1 × U (1) 2 symmetry is broken down to SU (2) × U (1) 2 as Σ develops a vacuum expectation value F .We identify this gauge group with SU (2) L × U (1) Y of the SM.Then this symmetry is further broken down to U (1) EM as the Φ 1 develops a vacuum expectation value v (= 246 GeV).After symmetry breaking, collecting the relevant terms, the neutral gauge boson mass matrix can be written as where g1 g0 = x ≪ 1, g2 g1 = t, and v F = r.This matrix can be diagonlized perturbatively in the small parameter x.The eigenvalues corresponding to the eigenstates of the photon, Z and Z ′ bosons are given by and the mass eigenstates are (writing t = s c ) ) (2.9) The charged gauge boson mass matrix can similarly be computed and given by . (2.10) The mass eigenvalues are The mass eigenstates of the charged gauge bosons are where the one with the smaller mass is identified as the W boson and the heavier one as the W ′ boson.It is to be noticed that the massive SM gauge bosons get their contribution predominantly from SU (2) 1 .However, for the BSM gauge bosons, SU (2) 0 is the major contributor.In the limit g 0 → ∞ (or equivalently x → 0), the m W ′ ,Z ′ decouple and the gauge spectrum at low energies thus reduces to that of the SM as expected.For all practical purposes x can be parametrized as the ratio of the light and heavy gauge boson masses (to leading order) as which can be readily seen from Eqns.2.11 and 2.6.After taking into account the gauge charges of the scalars and the imposed Z 2 symmetry, the most general scalar potential can be written as (2.15) After spontaneous symmetry breaking the higgs mass (m 2 h ) is given by the expression 2λ 1 v 2 and the charged scalar mass (m 2 ϕ ± ) is given by m 2 2 −λ 23 F 2 .The two neutral components of Φ 2 have masses given by the following expressions: It is important to note here that the Z 2 symmetry of the Φ 2 bi-doublet works as the stabilising symmetry required for the DM under the WIMP scenario.However, this discrete symmetry can not forbid the number changing processes among the components of this bidoublet.Thus for a positive (λ 23 + λ23 ), the ϕ 3 component will be the natural DM candidate and for negative values of this quantity, the ϕ 0 component serves the same purpose from the stability argument.However, if the couplings are such that these two are degenerate, then we have an interesting possibility of a multicomponent DM scenario where both neutral components are cosmologically stable DM candidates -we describe how this will play out in more detail in Sec.IV.The Yukawa sector is assumed to mimic that of the SM as Φ 1 plays a similar role of the SM Higgs doublet.However, the fermions can still couple to the heavy gauge bosons due to gauge mixing.All the relevant couplings can be straightforwardly computed (again in a perturbative expansion in x) and are listed in Appendix B. Specifically, the purely scalar sector couplings are given in Tables III and IV.The scalar couplings with the gauge bosons are listed in Tables V,VI,VII,VIII,IX,X and XI.The relevant couplings of the gauge bosons with the fermions can be found in Tables XII and XIII and that of the gauge bosons among themselves in Table XIV.

III. MODEL CONSTRAINTS
In this section, we impose all the theoretical and experimental constraints on the model.We will not only identify all the relevant parameters for carrying out further analysis, but also assign values to some and constrain the rest with various experimental observations.Among the many model parameters, the following are the most relevant for our analysis: v, m h , m ϕ3 , λ 12 , λ 23 , λ23 , m W ′ (or m Z ′ ) and F (or equivalently r).The other important parameters can be represented in terms of these above mentioned ones.For example, x can be represented in terms of m Z ′ and F following Eqn.2.6 or more conveniently through the mass ratio in Eqn.2.11 (for a given v, F and r can be used interchangeably).We also set m h = 125 GeV and v = 246 GeV in this paper.We now look at four key constraints on the model parameters coming from vacuum stability, the branching ratio of invisible decay of the Higgs, the Higgs alignment limit and collider constraints on the heavy gauge bosons.

III.1. Vacuum Stability
It is the quartic part of the scalar potential that contributes to the vacuum stability conditions [35].After spontaneous symmetry breaking when Σ develops a vacuum expectation value, the quartic part of the scalar potential in our model can be written in the unitary gauge ⟨Σ⟩ = 1 as, The vacuum stability conditions can now be found through the standard procedure described in [35].Let's assume that Φ • Let us consider first the b = 0 direction (a → ∞).We demand V 4 ≥ 0, which implies λ 1 a 2 ≥ 0. The condition is then λ 1 ≥ 0.
• Next we consider the a = 0 direction (b → ∞).Similarly by demanding V 4 to be positive, one can derive the condition λ 2 ≥ 0.
• Lastly, the direction a = λ2 λ1 b (where a, b → ∞) can be considered.Again demanding V 4 ≥ 0 one can find the condition Now collecting all the equations together the vacuum stability conditions for this model are given by Another relevant theoretical constraint, the perturbativity bound, is discussed in Appendix A.

III.2. Higgs Invisible Decay
The Higgs invisible decay branching ratio is given by Among ATLAS [36] and CMS [37] measurements of this quantity, the more stringent bound comes from the ATLAS experiment and is 11%.In Fig. 1, assuming ϕ 3 to be the dark matter candidate, we present this bound in the parameter space m ϕ3 − λ 12 for dark matter masses smaller than half of m h (in order, of course, that h can decay into a pair of them).There is an upper bound on λ 12 as a function of m ϕ3 as indicated by the shaded region in the figure.
For higher DM masses the restriction will be drastically lenient due to off-shell suppression of the decay width.

III.3. Higgs Alignment Limit
Given the extended gauge and scalar sector in our model, the standard Higgs boson couplings to various constituents of the SM are modified due to gauge mixing.For example, the coupling between Higgs and the pair of W boson (Table V ) will take the form One can notice that the deviation of this coupling from the SM counterpart is proportional to x 2 .A similar deviation can be found in other SM-like gauge couplings in Tables V and VI.On the other hand the couplings between the SM Higgs and different fermion pairs in this model maintain the same structure as SM due to the Yukawa sector of both these models being the same.To analyse how new physics can modify Higgs properties we look into the branching ratio of h → γγ.Apart from the usual SM fermion and W boson loops, the heavy W ′ boson as well as the charged scalar ϕ ± will contribute to this loop mediated process.After incorporating these additional contributions we have calculated the branching ratio BR (h → γγ) (for a detailed calculation, see Appendix C) in this model for a range of charged Higgs and W ′ masses.In Fig. 2, we present the allowed region (shaded region is disallowed) of the parameter space after recasting ATLAS [38] and CMS [39] measurements -we have set the ratio between two vev's r at 0.25 and presented our results for two different λ 12 benchmark values given by 0.1 and 0.01.It is found that for both the scenarios, the region m W ′ > 1 TeV and m ϕ ± > 200 GeV is allowed.
The dark cyan shaded region in the m ϕ ± vs. m W ′ plane is disallowed from the branching ratio measurements of h → γγ channel.

III.4. Collider Constraints on Heavy Gauge Bosons
Due to its extended gauge sector, the model admits heavy charged and neutral gauge bosons.The Z ′ couples to different SM fields primarily through gauge mixing.As a result, this Z ′ , if produced in colliders, can be discovered from its decays to different SM final states.Specifically, the Z ′ can decay into ℓ L lL (ℓ represents all possible charged leptons), ν L νL and q L qL (q represents all possible SM coloured particles) pairs where the the corresponding coupling is suppressed by the parameter x w.r.t the SM counterpart.One can notice in Tables XII and XIII that the Z ′ specifically couples to the left-handed chiral fermions since it is induced purely from a gauge mixing in the neutral gauge sector and does not receive any contribution from the U (1) 2 (see Eqn. 2.9).In addition to the fermion pairs, the Z ′ can decay into two other SM final states Zh and W W . Apart from these SM final states the Z ′ boson does couple to the dark sector particles ϕ + ϕ − and ϕ 0 ϕ 3 .The three point vertices corresponding to these two decay modes are proportional to 1 x at leading order (see Table XI).As a consequence the decay width of these two modes are significantly larger than SM final states.In Table II we present the analytic formulas for these decay widths.

Decay Modes
Decay rate expressions The relevant Z ′ decay rates for the purposes of this section -the expressions for the couplings given herein can be found in Appendix B.
Here N c stands for the colour factor (3 and 1 for quarks and leptons respectively) and the λ's denote the relevant cubic vertices.We have neglected the mass of the fermions in comparison with the Z ′ mass for the fermionic decay modes.In principle, there exists a cubic vertex of the form λ Z ′ W W ′ .However, in the present model m Z ′ and m W ′ are degenerate and the Z ′ decay via W W ′ channel is kinematically forbidden.In Fig. 3 (left) we present the corresponding BRs of these channels for m Z ′ ranging from 1 TeV to 5 TeV.As expected, the BR of the ϕ + ϕ − and ϕ 0 ϕ 3 modes dominate for the entire m Z ′ range.Among the SM final states the W W has the maximum and Zh the minimum branching ratio.For the present calculation we have set r = 0.25 and m ϕ0 = m ϕ3 = m ϕ ± = 400 GeV.In Fig. 3 (right) we present the There is a visible growth in this ratio at the mass range m Z ′ > 2.5 TeV.However, the overall value remains below ≲ 30% for the entire range.This is to emphasize that for the current benchmark choice the perturbative treatment is applicable.For different parameter choices, one should explicitly calculate the Γ Z ′ m Z ′ and ensure that the corresponding value is less than unity.From the above discussion, it is evident that the Z ′ in this model primarily decays via exotic scalars coming from the scalar bidoublet and that the branching ratios in the SM final states are significantly low.Nevertheless, it is useful to check the direct detection limits to understand the allowed parameter space.Before the LHC era, the LEP and Tevatron colliders have extensively searched for the heavy neutral gauge bosons.The LEP searches [40] can be divided into two categories: i) precision measurement around the Z pole which can effectively constraint the mixing angle between the Z − Z ′ and ii) direct searches for the heavy neutral boson in the e + e − → f f (here f is all possible SM fermions) at the mass range above Z-pole.On other hand, at Tevatron the Z ′ boson was primarily looked for in the pp → Z ′ X → ℓ + ℓ − X channel and using the data, both the CDF [41,42] and D0 [43,44] collaborations have put bounds on m Z ′ for different U (1) extensions of the SM.The commonality between the LEP and Tevatron is the corresponding bounds are strictly dependent on the Z ′ couplings to SM fermions.For a comprehensive analysis on how these results can be translated in a specific model, the interested reader can consult Ref. [45,46].At LHC, both the Drell-Yan (DY) and Vector Boson Fusion (VBF) production mechanisms are considered for the Z ′ searches.
After production, the Z ′ can decay to dileptons, Zh, and W W final states.For the present study, we will consider these three final states and will recast the corresponding ATLAS and CMS bounds in our model.To calculate the signal cross section in our model we have used the prescription illustrated in Ref. [47,48].We begin with the process pp → Z ′ Zh, where the Z ′ is produced via DY and decays into Zh.The Z further decays leptonically and the Higgs boson decays to bb.Both the ATLAS [49] and CMS [50] have searched the Z ′ boson in this mode.The model dependent upper bounds on the σ × BR is displayed in Fig. 4 (upper panel left) by red dashed (ATLAS) and blue dotted curve (CMS) curves.The black solid line represents the corresponding number for the present model.
In the case of VBF mechanism the Z ′ boson can decay via W W mode.After that, one of the W bosons decays hadronically and the other one decays leptonically.In Fig. 4 (upper panel right), we display the corresponding ATLAS, CMS and model specific σ × BR values via red dashed, blue dashed and black solid line respectively.We have also translated the LHC limits in dilepton final state where the Z ′ is produced via DY mechanism.The corresponding limits [53,54] are displayed in Fig. 4 (lower panel).Considering these search channels one can see that the most stringent bound can be obtained from the Zh mode and to satisfy current collider limits one have to set m Z ′ ≳ 2.5 TeV for the specific benchmark points that we have chosen.Apart from the heavy neutral gauge bosons the LHC also looked for the additional charged gauge bosons.However the number of searches is comparatively lower as well as less stringent in the latter case.

IV. DARK MATTER PHENOMENOLOGY
In this section we consider our model described in Sec.II as a plausible solution to the dark matter problem.In what follows, we discuss the choice of the DM candidate, calculate the relic abundance for a multicomponent scenario, and address the direct detection constraints.

IV.1. Choice of Dark Matter Candidate
For generic choices of parameters in this model, the couplings mentioned in Table XI allow the decays ϕ 0 → ϕ ± W ′∓ , ϕ 3 Z ′ for a positive λ 23 + λ23 , and ϕ 3 → ϕ 0 Z ′ for negative values of λ 23 + λ23 .Such decays, despite being suppressed by the heavy gauge boson mass, would still make one of the scalars unstable.In other words, constraining such decays with the age of universe would push the heavy gauge boson masses to very high values essentially decoupling them from SM gauge sector for all practical purposes.While this is certainly a reasonable path to take, in this work our aim is to treat the multi-component DM scenario while maintaining a meaningful way of discussing the gauge sector in the context of present experiments and constraints.Thus, we work under the mass degeneracy limit of the Z 2 odd sector that makes the neutral components cosmologically stable DM candidates.Referring to Eqn. 2.16, we see that this corresponds to the particular parameter choice λ23 = −λ 23 .The charged components will completely annihilate via photon couplings.For the mass degenerate scenario, all components of the Z 2 odd sector are stable at tree level in the context of collider production as well.However, in the absence of higher dimensional terms in our scalar Lagrangian, loop level effects can give rise to a sub-GeV mass non-degeneracy between the charged and neutral components [55] leading to charged scalar decays via MeV scale charged pions as in reference [24], thus avoiding any constraints pertaining to the production of an otherwise stable charged particle in colliders which could potentially disallow charged scalar masses (and thus the DM masses) upto a few TeV.Given that the present model admits an extra SU (2) gauge group in addition to the SM one, we expect such loop effects to be even more pronounced in our case in providing a mass splitting between the charged and the neutral scalars.

IV.2. Relic Abundance
In the early universe all the particles were in equilibrium in the high temperature cosmic soup.Over time as the universe expanded and the temperature of the universe started to drop, the interaction rate of WIMPs with the other particles of the cosmic plasma became comparable to the Hubble parameter.It became difficult for the WIMPs to find another WIMP to annihilate with and thus they were frozen out remaining in the universe as a relic which one can compare with the observed DM relic abundance.The value of relic abundance as measured by Planck satellite is given by Ωh 2 = 0.120 ± 0.001 [56].
There are two DM candidates in this model as argued in Sec.IV.1.The dark matter candidates can annihilate mainly via higgs mediated channels to SM states, for example b b, τ τ , W + W − , ZZ, hh and t t, as allowed by the kinematics in various DM mass range.The Feynman diagrams are given in Appendix D. Among the individual contributions of the DM candidates to the relic, the contribution from ϕ 3 is much less than that of ϕ 0 for the DM mass range starting from m W upto around 300 GeV.This is due to the fact that in this model there exists a ϕ 3 ϕ 3 W + W − coupling (see Table VII) and also ϕ 3 ϕ ± W ∓ couplings (see Table XI) which is absent for ϕ 0 , and thus ϕ 3 annihilates more through those channels leaving less relic in the universe.These couplings come into effect around a DM mass of m W and finally go away as the co-annihilation channels start to dominate at a DM mass of around 300 GeV. 1 As we move into the higher DM mass range, co-annihilation among ϕ 0 and ϕ 3 via Z ′ mediated diagrams into W W , Zh and quark and lepton pairs start to dictate the relic abundance.All these properties are evident in Fig. 5 where the relic abundance contribution of each dark matter candidate as well as the total contribution have been displayed as a function of dark matter mass for a mass degenerate DM scenario.The Z ′ mass is fixed at 5 TeV as a benchmark choice and the relic varies with scalar sector parameter λ 12 in the annihilation dominated low DM mass region whereas at high DM mass region the λ 12 dependence is absent since the Z ′ mediated co-annihilation starts dictating the relic.There are two peaks in the velocity averaged cross section due to the Breit-Wigner resonance which are related to the mass of the two mediators h and Z ′ .DM relic is dictated by various annihilation and co-annihilation channels in various DM mass ranges: • The DM annihilation to bb is the dominant channel dictating the relic in the DM mass range upto around the W boson mass i.e. m W ∼ 80.4 GeV.The τ τ is the second most dominant process in this DM mass region.The amplitude for the process to bb final state is given by • One can observe the sudden decline in the relic curve for both the individual contribution and the total relic around the DM mass of m W as the Higgs mediated annihilation channel to W W final state becomes kinematically accessible.The amplitude for the same is given by for ϕ 0 and • There is a further sudden decline in the relic is evident in the ϕ 0 contribution and the total relic around the Higgs mass (i.e.m h ∼ 125 GeV) as the annihilation to hh final state becomes dominant.The case for ϕ 3 contribution to relic is slightly different as mentioned above.The amplitude goes as (4.4) The exact DM mass value can be different than 300 GeV depending upon relevant couplings and the heavy gauge boson mass.
• Finally, at DM masses higher than a few hundred GeV the Z ′ mediated co-annihilation cross-section takes over and the Z ′ mediated channel to W W final state emerges as the dominant contribution.The corresponding amplitude is proportional to where ψ is the angle between the incident and emergent particle momenta to be summed over.The region above the black line is disallowed by overabundance as per Planck data.The relevant Boltzmann equations are given in Appendix D. The relative contribution of the individual DM components to the total relic abundance can be defined as a relic fraction given by for the ϕ i component.The relic fraction of each of the components is shown in Fig. 6.One can observe that the contribution to the relic from individual components are of the same order for DM masses < 80 GeV due to similar diagrams contributing to both of them.For DM masses > 80 GeV, the diagrams with W W as final states become kinematically accessible.The diagrams shown in Fig. 14 only exist for the ϕ 3 component and are absent for the ϕ 0 component.Hence, the ϕ 3 component annihilates dominantly through these channels, especially through the ϕ 3 ϕ 3 W W vertex, leading to a very small contribution to the relic.The ϕ 0 on the other hand furnishes the major contribution to the relic.As one explores higher DM masses, the Z ′ mediated diagrams start to become relevant.The contributions are again comparable at the Breit-Wigner resonance region corresponding to the Z ′ mediator (i.e.around m ϕ3 ∼ m Z ′ 2 ∼ 2.5 TeV).Next, we turn to the total relic abundance as a function of DM mass in Fig. 7 but this time for various values of Z ′ mass and two benchmark values of λ 12 to complement our results in Fig. 5.For the higher DM mass range not much difference can be noticed in this logarithmic scale plot for various m Z ′ values apart from the shift in the Breit-Wigner resonance.However, in the lower DM mass range there is strong competition from the Higgs mediated channels.A heavier Z ′ mediated process generates less cross section compared to the Higgs mediated processes in low DM mass region.However, a lighter Z ′ mediated process can produce comparable cross section to the higgs mediated ones and thus contribute to the relic significantly even in the low DM mass range as evident in Fig. 7 (right).

IV.3. Direct Detection Constraints
Direct detection experiments like Xenon1T [57,58], PandaX-4T [59] and LUX-ZEPLIN (LZ) [60] constrain the parameter space by non-observation as they provide upper limits on the scattering cross-section of the DM candidates with the detector nuclei in the non-relativistic limit.The expression2 for spin independent scattering cross section of ϕ 0 with the nucleons of the detector material is given by [22] where, µ 0 = m ϕ 0 m N m ϕ 0 +m N .The matrix elements (f N T i ) are given in reference [61].In addition to the Higgs mediated ϕ 0 N → ϕ 0 N scattering process mentioned above, this model also exhibits a Z ′ mediated ϕ 0 N → ϕ 3 N process.In the mass degeneracy limit of the DM candidates, both the contributions are to be added together to arrive at the combined cross-section.The relevant couplings λ ϕ0ϕ3Z ′ = −λ ϕ3ϕ0Z ′ and λ Z ′ ūL u L = −λ Z ′ dL d L are given in Tables XI and XII.In the multicomponent dark matter scenario, the quantity of relevance to realise the constraints from direct detection experiments by non-observation is given by the following effective cross-section3 [63][64][65], where the DM-nucleon scattering cross-section is modified with each DM candidate's respective contribution to the relic abundance and σ SI N,i is the summation of the contributions from both diagrams ϕ i N → ϕ i N and ϕ i N → ϕ j N .Figure 8 describes the experimental constraints coming from Xenon1T, PandaX-4T and LZ experiments on the effective cross-section of the DM candidate ϕ 0 as a function of DM mass where LZ is clearly the most restricting of them all.A similar set of constraints can be obtained for the ϕ 3 as well following the same procedure.It is evident from Fig. 8 that for a Z ′ mass of 2.5 TeV (motivated by the lower bound from Fig. 4), much of the DM mass parameter space is already excluded upto around 2 TeV for both the benchmark values of coupling parameter λ 12 = 10 −2 and 10 −3 .Keeping every other benchmark value same and modifying only the Z ′ mass to 5 TeV brings the DM mass exclusion region down to ∼ 100 − 200 GeV range and very interestingly opens up some parameter space in the lower DM mass region of around 10 − 15 GeV for λ 12 = 10 −3 .A more careful look at the figures makes it apparent that for λ 12 = 10 −2 , the Higgs mediated diagram plays an important role upto the DM masses around the weak scale and beyond that the heavy gauge boson mediated diagram takes over.However, the contribution of the Higgs mediated diagram remains subdominant to the one mediated by the heavy gauge boson throughout the DM mass range for λ 12 = 10 −3 .As expected, the cross section is more suppressed by a higher value of the Z ′ mass compared to a lower value.The slightly raised plateau like shape of the curves in the mid-range of DM masses depicted here are the result of the fractional contributions of the individual DM components to the relic.One also has to consider the neutrino floor [66] and the neutrino fog [67] in this regard as for cross-sections too small the DM scattering will be difficult to tell apart due to the presence of neutrino scattering.Similar results are expected for ϕ 3 upto the relic fraction as is realized in Fig. 9.The rise (dip) in the curves are indicative of the substantial (negligible) contribution of ϕ 3 to the relic in that DM mass region.Even though in this work we have worked under the mass degeneracy limit of the DM candidates (i.e.m ϕ0 = m ϕ3 ), there is also in principle a scope for inelastic scattering between the two DM components by softly breaking their mass degeneracy.However, that would lead to some problems with the stability of the DM components and this multicomponent DM scenario will effectively reduce to a single component scenario as discussed in Sec.IV.1.

V. CONCLUSION
In this work we have extended the SM with an SU(2) gauge sector and the scalar sector with a bidoublet and a non-linear sigma field.The neutral components of the Z 2 protected bidoublet are suitable DM candidates.In the DM mass range starting from a few GeV upto around half of the Higgs mass, the constraint from invisible decay of higgs becomes important and restricts the most relevant scalar sector coupling parameter λ 12 with an upper bound of O(10 −3 − 10 −2 ).However, the sensitivity of direct detection experiments provide a much more stringent bound for a typical scalar extension under the WIMP scenario.In many scalar extensions, one has to resort to the Breit-Wigner resonance to account for the DM relic abundance while simultaneously satisfying the direct detection constraints.In our multicomponent DM model the extended gauge sector mitigates those strict constraints to a great extent.The gauge sector couplings provide additional channels for the DM annihilation and co-annihilation thus increasing the velocity averaged cross-section and reducing the relic to avoid overabundance.In the similar coupling range this model can also be constrained with the direct detection effective cross-section values from the experimental observations.Furthermore, we found that a Z ′ mediated scattering diagram can dominate a scalar mediated diagram in a sizeable

Appendix B: Couplings
The couplings are arranged in a tabular form as the following (where θ is the Weinberg angle): TABLE VIII.Quartic couplings for charged scalar decay.The CP-even Higgs boson can decay into di-photon channel through loop mediated processes.In case of SM, the W boson loop and different fermion loops participate in this process.For the present model, in addition to these loop diagrams the heavy charged gauge boson W ′ and the charged scalar ϕ ± can also engage in the loop induced processes.In Fig. 12 we show the all the Feynman diagrams that can contribute to these process.Ref. [68] have calculated this where I i h represent the individual loop contribution from different particles that flowing in the loop.In Eq.C2, we present the explicit expression for I i h .
Here N c and Q f represents the color factor and electromagnetic charges for each SM fermions.The R h i s are the relative coupling strength of the three point vertices λ hii (where i can vary based on the particle flowing in the loop) w.r.t corresponding SM coupling.For the present model R h f = 1 as the fermion sector is identical to the SM.In case of gauge bosons and the charged scalar the relative coupling strength would be R h W = 1 + x 2 sin 2 θ W , R h W ′ = x 2 and R h ϕ ± = 2λ ab g 2 respectively.To get the sense of these values one should the check the corresponding vertices that are present in appendix B. The F j (τ i ) are the form factors corresponding to the gauge, fermion and scalar loops which can take following form. where and depending on the value of τ the f (τ ) can take following form

12 FIG. 1 .
FIG.1.The shaded region in the λ12 − m ϕ 3 space is disallowed by constraints from Higgs invisible branching ratio.

FIG. 3 .
FIG. 3. Left: The branching ratio of Z ′ for all possible two body final states.Here we have set r = 0.25 and m ϕ 0 = m ϕ 3 = m ϕ ± = 400 GeV.The branching ratio for Z ′ → qq mode (solid orange line) is summed over all six SM quark final states.Similarly, the Z ′ → ℓ + ℓ − (blue solid) and Z ′ → νν (green dashed) modes are calculated while considering all three lepton flavours.

FIG. 4 .
FIG. 4. The direct search limits in the channels DY → Z ′ → Zh (upper panel left), VBF → Z ′ jj → W W jj (upper panel right) and DY → ℓ + ℓ − (lower panel).The red-dashed, blue-dotted and black solid line denote the ATLAS, CMS and the model specific curves respectively.The arXiv IDs which are written beside the ATLAS and CMS highlight the corresponding articles [49-54] which are used in this calculation.

FIG. 5 .
FIG.5.Relic abundance as a function of dark matter mass for various couplings (lagrangian parameters) in the multi-component dark matter scenario with annihilation via h portal and co-annihilation mostly via Z ′ portal for m Z ′ = 5 TeV.

2 For λ 12
FIG.7.Relic abundance as a function of dark matter mass for various Z ′ masses in the multi-component dark matter scenario.

FIG. 10 .m
FIG.10.The variation of x in M Z ′ versus r plane.For details please see the text.

FIG. 12 .
FIG.12.All possible loop mediated processes that can contribute to h → γγ decay mode.

TABLE III .
Cubic scalar couplings.

TABLE VII .
Quartic gauge couplings for DM annihilation.

TABLE IX .
Quartic couplings for charged scalar annihilation.

TABLE X .
A few more couplings between the scalar and gauge sectors.

TABLE XII .
Cubic Couplings of gauge fields with quarks.

TABLE XIII .
Cubic Couplings of gauge fields with leptons.

TABLE XIV .
Relevant couplings of gauge bosons among themselves.
Appendix C: h → γγ Decay Width