Dark matter in Inert Doublet Model with one scalar singlet and $U(1)_X$ gauge symmetry

TL;DR

This paper explores a dark matter model extending the Standard Model with an extra $U(1)_X$ gauge symmetry, analyzing its parameter space, experimental constraints, and potential detection prospects at future colliders.

Contribution

It introduces a novel inert doublet model with a scalar singlet and $U(1)_X$ symmetry, providing detailed phenomenological analysis and experimental viability of the dark matter candidate.

Findings

Allowed parameter regions consistent with current experimental limits.

Dark matter mass range from light to heavy, depending on model parameters.

Future collider detection possible for dark matter masses between 10-60 GeV.

Abstract

We study Dark Matter (DM) abundance in the framework of the extension of the Standard Model (SM) with an additional $U(1)_X$ gauge symmetry. One complex singlet is included to break the $U(1)_X$ gauge symmetry, meanwhile one of the doublets is considered inert to introduce a DM candidate. The stability of the DM candidate is analyzed with a continuous $U(1)_X$ gauge symmetry as well as discrete $Z_2$ symmetry. We find allowed regions for the free model parameters which are in agreement with the most up-to-date experimental results reported by CMS and ATLAS collaborations, the upper limit on WIMP-nucleon cross section imposed by XENON1T collaboration and the upper limit on the production cross-section of a $Z^{\prime}$ gauge boson times the branching ratio of the $Z^{\prime}$ boson decaying into $\ell^-\ell^+$. We also obtain allowed regions for the DM candidate mass from the relic density reported by the PLANCK collaboration including light, intermediate and heavy masses; depending mainly on two parameters of the scalar potential, $\lambda_{2x}$ and $\lambda_{345}=\lambda_3+\lambda_4+2\lambda_5$. We find that trough $pp\rightarrow \chi\chi \gamma$ production, it may only be possible for a future hadron-hadron Circular Collider (FCC-hh) to be able to detect a DM candidate within the range of masses 10-60 GeV.

Tables

Table 1. Table 1: IDMS couplings involved in the calculations of this work. We define λ 345 = λ 3 + λ 4 + 2 ​ λ 5 subscript 𝜆 345 subscript 𝜆 3 subscript 𝜆 4 2 subscript 𝜆 5 \lambda_{345}=\lambda_{3}+\lambda_{4}+2\lambda_{5} . For Z ′ ​ f i ​ f ¯ i superscript 𝑍 ′ subscript 𝑓 𝑖 subscript ¯ 𝑓 𝑖 Z^{\prime}f_{i}\bar{f}_{i} coupling we consider the limit when the kinetic mixing term ε → 0 → 𝜀 0 \varepsilon\to 0 .

Coupling	Expression
$h{f}_i{\overline{f}}_i$	$\frac{m_{f_i}}{\upsilon }\cos {\alpha }_1$
$h{H}^{-}{H}^{+}$	$\left({\lambda }_3\cos {\alpha }_1+{\lambda }_{2x}/2\right)\upsilon$
$h{W}_{\mu }^{-}{W}_{\nu }^{+}$	$g{m}_W\cos {\alpha }_1{g}_{\mu \nu }$
$h\chi \chi$	$\left({\lambda }_{345}\cos {\alpha }_1+{\lambda }_{2x}/2\right)\upsilon$
$Z_{\mu }^{\prime }{f}_i{\overline{f}}_i$	$\frac{g_x}{2}\left(1-{\gamma }^5\right){\gamma }^{\mu }$
$Z_{\mu }^{\prime }\chi \chi$	$\frac{g_x}{2}{\left(p_{Z_{\mu }^{\prime }}-p_{\chi }\right)}^{\mu }$
$S{f}_i{\overline{f}}_i$	$\frac{m_{f_i}}{\upsilon }\sin {\alpha }_1$
$S{W}_{\mu }^{-}{W}_{\nu }^{+}$	$g{m}_W\sin {\alpha }_1{g}_{\mu \nu }$
$S\chi \chi$	${\lambda }_{2x}{\upsilon }_x\cos {\alpha }_1-{\lambda }_{345}\upsilon \sin {\alpha }_1$
$\chi \left(A\right)H^{\pm }{W}_{\mu }^{\mp }$	$i\frac{g}{\sqrt{2}}{\left(p_{H^{\pm }}-p_{\chi \left(A\right)}\right)}^{\mu }\cos {\alpha }_2$

Table 2. Table 2: Relations between fermions charges under U ​ ( 1 ) X 𝑈 subscript 1 𝑋 U(1)_{X} to guarantee the anomaly cancellation.

Field	$U{\left(1\right)}_X$
$q_L$	$x_q$
$u_R$	$x_u=4x_q$
$d_R$	$x_d=-2x_q$
$l_L$	$x_l=-3x_q$
$e_R$	$x_e=-6x_q$

Table 3. Table 3: The setting of the values for the model parameters.

Parameter	Value	Constraint
$c_{{\alpha }_1}$	$0.99$	${ℛ}_{WW}$
$g_x$	$0.4$	${ℛ}_{ZZ}$ and ${\sigma }^{SI}\left(\chi N\to \chi N\right)$
${\upsilon }_x$	$23$ TeV	${ℛ}_{ZZ}$ and ${\sigma }^{SI}\left(\chi N\to \chi N\right)$
$m_{H^{\pm }}$	$0.5$ TeV	${ℛ}_{\gamma \gamma }$
$m_{Z^{\prime }}$	$3$ TeV	$\sigma \left(pp\to Z^{\prime }\right)ℬ\left(Z^{\prime }\to \ell \ell \right)$

Equations

$${\Phi }_1$$

$${\Phi }_2$$

$${\mathcal{S}}_X$$

$${\Phi }_1$$

$${\Phi }_2$$

$${\mathcal{S}}_X$$

$$G_{\text{SM}}\otimes U(1{)}_X\stackrel{⟨{\mathcal{S}}_X⟩}{\to }{G}_{\text{SM}}\stackrel{⟨{\Phi }_1⟩}{\to }SU(3{)}_C\otimes U(1{)}_{\text{EM}},$$

$$V$$

$$M_{0}^{2}=\left(\begin{array}[]{cccc}M_{11}&M_{12}&M_{13}&0\\ M_{12}&M_{22}&M_{23}&0\\ M_{13}&M_{23}&M_{33}&M_{34}\\ 0&0&M_{34}&M_{44}\end{array}\right),$$

$$M_{11}$$

$$M_{22}$$

$$M_{33}$$

$$M_{34}$$

$$M_{44}$$

$$\left(\begin{array}[]{c}h\\ S\end{array}\right)=\left(\begin{array}[]{cc}\cos\alpha_{1}&-\sin\alpha_{1}\\ \sin\alpha_{1}&\cos\alpha_{1}\end{array}\right)\left(\begin{array}[]{c}\phi_{1}\\ s_{x}\end{array}\right)$$

$$\left(\begin{array}[]{c}\chi\\ A\end{array}\right)=\left(\begin{array}[]{cc}\cos\alpha_{2}&-\sin\alpha_{2}\\ \sin\alpha_{2}&\cos\alpha_{2}\end{array}\right)\left(\begin{array}[]{c}\phi_{2}\\ \eta_{2}\end{array}\right),$$

$$m_{S,h}^2={\lambda }_1{\upsilon }^2+{\lambda }_x{\upsilon }_x^2\pm ({\lambda }_1{\upsilon }^2+{\lambda }_x{\upsilon }_x^2)\sqrt{1+r_1^2},$$

$$m_{H^{\pm }}^2$$

$$m_{A,\chi }^2$$

$$m_h^2$$

$$m_S^2$$

$${\mathcal{L}}_{Kin}=-\frac{1}{4}{\hat{B}}_{\mu \nu }{\hat{B}}^{\mu \nu }+\frac{1}{2}\frac{\epsilon }{\cos {\theta }_W}{\hat{B}}^{\mu \nu }{\hat{Z}}_{0\mu \nu }^{\prime }-\frac{1}{4}{\hat{Z}}_{0\mu \nu }^{\prime }{\hat{Z}}^{\prime 0\mu \nu },$$

$$\displaystyle\left(\begin{array}[]{c}Z^{\prime}_{0\mu}\\ B_{\mu}\end{array}\right)=\left(\begin{array}[]{cc}\sqrt{1-\varepsilon^{2}/\cos^{2}\theta_{W}}&0\\ -\varepsilon/\cos^{2}\theta_{W}&1\end{array}\right)\left(\begin{array}[]{c}\hat{Z^{\prime}}_{0\mu}\\ \hat{B}_{\mu}\end{array}\right),$$

$$A_{\mu }={\hat{A}}_{\mu }-\epsilon {\hat{Z}}_{0\mu }^{\prime },$$

$$Z_{0\mu }={\hat{Z}}_{0\mu }+\epsilon \tan {\theta }_W{\hat{Z}}_{0\mu }^{\prime },$$

$$Z_{0\mu }^{\prime }={\hat{Z}}_{0\mu }^{\prime }.$$

$${\mathcal{L}}_{\text{scalar}}$$

$$D_{\mu }=\left({\partial }_{\mu }+i{g}^{\prime }Y{\hat{B}}_{\mu }+ig{T}_3{\hat{W}}_{3\mu }+i{g}_x{Q}_i^{\prime }{\hat{Z}}_{0\mu }^{\prime }\right),$$

$${\mathcal{L}}_{\text{scalar}}=\frac{1}{2}{m}_{Z^{\prime }}^2{Z}^{\prime 0}{Z}^{\prime 0}+\frac{1}{2}{m}_Z^2{Z}^0{Z}^0-{\Delta }^2{Z}^0{Z}^{\prime 0}+...,$$

$$m_{Z^{\prime }}^2$$

$${\Delta }^2=\frac{1}{2}{g}_Z\left(\frac{g^{\prime }\epsilon }{2\cos {\theta }_W}+g_{x}x\right){\upsilon }^2,$$

$$m_Z^2=g^2\frac{{\upsilon }^2}{4{\cos }^2{\theta }_W}.$$

$$\displaystyle\left(\begin{array}[]{c}Z\\ Z^{\prime}\end{array}\right)=\left(\begin{array}[]{cc}\cos\xi&-\sin\xi\\ \sin\xi&\cos\xi\end{array}\right)\left(\begin{array}[]{c}Z^{0}\\ Z^{\prime 0}\end{array}\right),$$

$${\mathcal{L}}_{Yukawa}=$$

$$\mp x_1-x_q+x_{u,d}=0$$

Full text

††journal: Eur. Phys. J. C

1 Introduction

Cosmological observations have shown anomalies that establish the existence of non-luminous matter as a possible solution. This non-luminous matter was called Dark Matter (DM) by F. Zwicky zwicky . F. Zwicky applied the virial theorem to the Coma Cluster and concluded that a large amount of non-luminous matter must be considered to keep the system bound together. Forty years later, V. Rubin and N. Thonnard found gravitational evidence through the rotation curve of spiral galaxies rubin ; rubin1 ; rubin2 ; rubin3 . Several proposals arose to explain the observations, namely, modified gravity milgrom , a dark component of matter bradac ; bradac1 ; bradac2 ; bradac3 , non-baryonic DM bergstrom . Nowadays observations suggest the existence of non-baryonic DM as the most viable solution. The PLANCK collaboration reveals that cold non-baryonic content of the matter density is $\Omega h^{2}=0.120\pm 0.001$ PLANCK , constituting about $25\%$ of the energy content of the universe.

As it is well known, the Standard Model (SM) of particle physics SM ; SM1 ; SM2 ; Cruz:2019vuo does not provide answer to fundamental issues; in particular, we can highlight the absence of a DM candidate, which motivates to extend the SM, opening the door to possible new physics beyond SM (BSM). BSM can include one or more scalar fields to introduce DM candidate, which corresponds to the simplest type of DM known as Weak Interacting Massive Particle (WIMP) IDM1 ; IDM2 ; IDM3 ; IDM4 ; IDM5 ; SUSYDM ; Bertone1 ; Bertone2 . The scalar particle as DM candidate must satisfy experimental and theoretical constraints 10testpoint ; for instance: it must have the right relic density, neutral particle and it must be consistent with direct DM searches. On the experimental side, searches for WIMPs based on different methodologies are realized by collaborations such as given in Refs. DirectDM ; DirectDM1 ; DirectDM2 , the CDMS cdms , CoGeNT cogent , Xenon xenon and LUX lux . The second one is through indirect searches IndirectDM by PAMELA pamela , ATIC atic , and Fermi LAT fermilat experiments for particles resulting from WIMP annihilation, for example, positron-electron pairs. Finally, DM search at colliders DMcollidersATLAS , DMcollidersCMS , such as the LHC, the WIMPs can be produced in pairs in association with other particles. A process to study DM at colliders is $pp\to\chi\chi+P$, where $\chi$ is a DM candidate and $P=g,\,\gamma,\,W,\,Z,\,H$.

The simplest proposal for a DM candidate is to extend the SM by introducing a singlet scalar field McDonald:1993ex . An interesting and simple model with a scalar field as DM candidate is the Inert Higgs Doublet Model (IDM) Deshpande:1977rw which contains a neutral scalar particle to play the role of WIMP Dolle:2009fn . The IDM shows an important dependence on the mass splitting parameter defined as the masses difference between pseudo-scalar and scalar coming from second doublet, the inert doublet. The heavy DM mass region for small values of the mass splitting parameter is obtained for masses from 500 GeV to 1000 GeV, meanwhile, the light DM mass region for the mass splitting parameter of the order of 50 GeV - 90 GeV is obtained for masses from 30 GeV to 80 GeV.

Other possibilities are Supersymmetry, which provides a WIMP candidate through the lightest neutralino susy ; susy1 , or universal extra dimension models with the lightest Kaluza-Klein partner as DM candidate ED ; ED1 ; ED2 ; ED3 ; ED4 .

In this work, we consider a model with an additional $U(1)_{X}$ gauge symmetry which includes two doublets and one complex singlet of scalar fields. One doublet is inert of which we identify a degree of freedom as a DM candidate, meanwhile, the other doublet is the usual SM doublet. The stability of the DM candidate is ensured by imposing a discrete $Z_{2}$ symmetry or by $U(1)_{X}$ gauge symmetry. Models with extra $U(1)_{X}$ gauge symmetries as extensions of the SM has many motivations. For example, Grand Unified and superstring theories contain additional $U(1)_{X}$ factors in the effective low energy limit. Supersymmetric extensions include theoretical and phenomenological aspects such as flavor physics, neutrino physics and DM Hewett:1988xc ; Suematsu:1994qm ; Demir:2005ti ; Langacker:2008yv . Extended models with a $U(1)_{X}$ gauge symmetry also have phenomenological importance because they predict a heavy vector gauge boson $Z^{\prime}$ derived from the spontaneous symmetry breaking (SSB) Martinez:2014rea , Camargo:2018klg . Besides, the $U(1)_{X}$ gauge symmetries can be incorporated in extended models that are free from triangle anomalies adding new fermions.

Previously, one of our authors published a paper with a similar approach, it can be found in the Ref. Martinez:2014ova , in which DM candidate mass of the order of $1.3$ GeV to $70$ GeV are allowed, depending on the assignment of the free parameters associated. The experimental data from LEP and relic density observation are considered to find an allowed mass of DM candidate of the order of 70 GeV in a scenario that assigns the parameters of the model as Higgs-phobic type, in which the $Z^{\prime}$ boson provides the channel of annihilation for DM suppressing the participation of the Higgs channel. The decay signal of Higgs diphoton also imposes strong restrictions through recent data from the CERN-LHC collider Martinez:2015kmn . When it is combined with the observed value of DM relic density, the allowed mass region is obtained such that $5$ GeV $\leq m_{\chi}\leq 62$ GeV for values of the order of $0.02\sim 0.08$ of the quartic coupling between doublets and singlet scalar in a model with $U(1)$ gauge symmetry Martinez:2014rea . In reference DMseveral ; DMseveral1 ; DMseveral2 ; DMseveral3 ; DMseveral4 ; DMseveral6 ; DMseveral7 ; DMseveral8 ; DMseveral9 ; DMseveral10 we include extensive literature to be consulted.

The organization of our research is as follows. In Sec. 2 we give a general view of the model. Sec. 3 is focused to constrain the free model parameters. In Sec. 4, we show the branching ratios for the $Z^{\prime}$ and the neutral scalar associated with the singlet field. Sec. 5 is devoted to the analysis of the relic density, we present our results and an analysis of them. The DM production at future colliders through $pp\to\chi\chi\gamma$ is presented in Sec. 6. Finally, in Sec. 7 conclusions are presented.

2 Inert Doublet Model plus a complex Singlet scalar (IDMS)

The IDMS incorporates a local $U(1)_{X}$ gauge symmetry and a $SU(2)$ scalar doublet to the SM gauge symmetry $G_{\text{SM}}=SU(3)_{C}\otimes SU(2)_{L}\otimes U(1)_{Y}$. The $Z^{\prime}$ gauge boson associated with $U(1)_{X}$ will provide an additional channel to the production and annihilation in scattering processes. On the other hand, the singlet scalar field is included to break down the $U(1)_{X}$ symmetry to $G_{\text{SM}}$. The DM candidate arises from the second doublet scalar field, which has a Vacuum Expectation Value (VEV) equal to zero to guarantee the stability of the DM candidate. But not only with a null value of VEV can achieve the stability of DM candidate, but it is also necessary a mechanism to control the couplings responsible for the DM candidate decays. Two possible options to control the stability of the DM are considered: a discrete $Z_{2}$ symmetry or the $U(1)_{X}$ gauge symmetry Haber:2018iwr ; Bonilla:2014xba ; Krawczyk:2015xhl .

2.1 Scalar fields

The scalar fields and their assignments under the $G_{\text{SM}}\otimes U(1)_{X}$ group are given by:

[TABLE]

where two first entries denote the representation under $SU(3)_{C}$ and $SU(2)_{L}$, respectively, meanwhile the hypercharge and charge under $U(1)_{X}$ are written in the last two entries. The scalar fields are written as follows:

[TABLE]

The spontaneous symmetry breaking (SSB) is achieved as

[TABLE]

where $\langle\mathcal{S}_{X}\rangle=\upsilon_{x}/\sqrt{2}$ and $\langle\Phi_{1}\rangle^{T}=(0,\upsilon/\sqrt{2})$ with $\upsilon=246$ GeV. Note that $\Phi_{2}$ must have VEV equal to zero to guarantee the stability of the DM candidate. The most general, renormalizable and gauge invariant potential is

[TABLE]

where $\mu_{1,\,2}^{2}$, $\lambda_{1,\,2,\,3,\,4,\,1x,\,2x}$ are real parameters and $\mu_{12}^{2}$, $\lambda_{5,\,6,\,7,\,12x}$ can be complex parameters. Note that $\mu_{2}^{2}>0$ because the DM candidate arises from $\Phi_{2}$, which has $<\Phi_{2}>=0$. The terms in the scalar potential that are proportional to $\Phi_{1}^{{\dagger}}\Phi_{2}\mathcal{S}_{X}$ or $\Phi_{2}^{{\dagger}}\Phi_{1}\mathcal{S}_{X}$ can generate a decay of DM candidate into two neutral scalars. We assume that $x_{2}-x_{1}\pm x\neq 0$ to leave these terms non-invariant under gauge symmetry. Thus, the parameters that accompany these terms must be zero to recover the gauge invariance and at the same time eliminate the couplings that are responsible for a decay of DM candidate at two neutral scalars.

After SSB the mass matrix for scalars in the $\{\phi_{1},\,s_{x},\,\phi_{2},\,\eta_{2}\}$ basis is

[TABLE]

where

[TABLE]

After the $M_{0}^{2}$ matrix is diagonalized and neutral scalars are rotated to physical states, the $M_{13}$ and $M_{23}$ matrix elements allow the mixing between neutral scalars and DM candidate, as shown in equation (4). This means that terms proportional to $\Phi_{1}^{{\dagger}}\Phi_{2}$ in the scalar potential must be eliminated, otherwise DM candidate will be unstable. For $x_{1}=x_{2}$ the terms proportional to $\Phi_{1}^{{\dagger}}\Phi_{2}$ in the scalar potential are gauge invariant. Then, it is required to introduce an additional discrete $Z_{2}$ symmetry for the doublets to eliminate these terms in the potential. Moreover, for $x_{1}\neq x_{2}$ the gauge invariance of the $U(1)_{X}$ symmetry guarantees the stability for the DM candidate. In either case, we will assume that $\lambda_{6}=\lambda_{7}=\lambda_{12x}=0$ in order to maintain the invariance under $Z_{2}$ or $U(1)_{X}$ symmetries.

2.2 $Z_{2}$ symmetry and $x_{1}=x_{2}$ case

The terms proportional to $\Phi_{1}^{{\dagger}}\Phi_{2}$ in the potential are invariant under $U(1)_{X}$; then it is necessary to introduce a $Z_{2}$ discrete symmetry to eliminate them. The proper assignment is $\Phi_{1}\rightarrow\Phi_{1}$ and $\Phi_{2}\rightarrow-\Phi_{2}$. Under the last assignment for the doublet, the $M_{13}=M_{23}=0$ and the mass matrix for the neutral scalar, Eq. (8), can be diagonalized by

[TABLE]

and

[TABLE]

where $\tan\alpha_{1,2}=\frac{r_{1,2}}{1+\sqrt{1+r_{1,2}^{2}}}$ with $r_{1}=\frac{\lambda_{1x}\upsilon\upsilon_{x}}{\lambda_{1}\upsilon^{2}-\lambda_{x}\upsilon_{x}^{2}}$ and $r_{2}=\frac{-\textrm{Im}[\lambda_{5}]}{\textrm{Re}[\lambda_{5}]}$ Cabral-Rosetti:2017mai . Therefore the masses for the scalars are

[TABLE]

while the $H^{\pm}$ charged scalar, $A$ pseudoscalar and $\chi$ masses are given, respectively, by

[TABLE]

We assume that $\chi$ plays the role of DM.

Two interesting limits can arise when approximations are realized about the $\lambda_{5}$ quartic couplings involved in $r_{2}$. The LHC results imply that $\upsilon\ll\upsilon_{x}$ and $\lambda_{1x}\sim 1$, then $r_{1}\approx-\frac{\upsilon}{\upsilon_{x}}$ and $\tan\alpha_{1}\approx-\frac{\upsilon}{2\upsilon_{x}}$. By considering $\textrm{Im}[\lambda_{5}]\sim\textrm{Re}[\lambda_{5}]$, then $r_{2}\approx-1$ and $\tan\alpha_{2}\approx-\frac{1}{1+\sqrt{2}}$. Moreover, if $\textrm{Im}[\lambda_{5}]=0$, which is the CP conservation case, then $r_{2}=0$ and $\tan\alpha_{2}=0$. By considering the previous approximation on $r_{1}$ and $\tan\alpha_{1}$, we can write Eq. (12) as

[TABLE]

An important fact is that the model allows the $\chi\to H^{\pm}W^{\mp}$ decay, whose $\chi H^{\pm}W^{\mp}$ coupling is shown in table 1. To avoid the instability of the DM candidate, we demand that the masses must satisfy $m_{H^{\pm}}^{2}\,(m_{A}^{2})>m_{\chi}^{2}$. To achieve this, from eqs. (13) and (14), we impose the following constraint: $\lambda_{4}>2|\lambda_{5}|$.

2.3 $U(1)_{X}$ gauge symmetry, $x_{1}\neq x_{2}$ case

The DM candidate can also be stable when $x_{1}\neq x_{2}$. In this case, the same parameters in the scalar potential, as the previous case, must be eliminated and $\lambda_{5}$ must be also zero. In addition, $\phi_{2}$ and $\eta_{2}$ are not mixing since $M_{34}=0$.

2.4 Gauge bosons interactions

The kinetic terms for the $U(1)_{Y}$ and $U(1)_{X}$ gauge symmetries are given by:

[TABLE]

where, $\hat{B}^{\mu\nu}$ and $\hat{Z^{\prime}_{0}}^{\mu\nu}$ are the field strength tensors defined by $\hat{F}_{\mu\nu}=\partial_{\mu}\hat{F}_{\nu}-\partial_{\nu}\hat{F}_{\mu}$ for $\hat{F}_{\nu}=\hat{B}_{\nu},\,\hat{Z}_{0\nu}$ Lee:2013fda ; Davoudiasl:2012ag . The mixing term between $\hat{B}_{\mu\nu}$ and $\hat{Z}_{0\mu\nu}^{\prime}$ is allowed by the gauge invariance. However, this mixing term can be eliminated by the field redefinition

[TABLE]

where the fields with hat notation contain the kinetic mixing term and $\varepsilon$ must be small to agree with the experiment. After SSB, the gauge bosons in the mass basis are

[TABLE]

[TABLE]

[TABLE]

The parameter $\varepsilon$ is assumed to be $\varepsilon\ll\cos\theta_{W}$ in order to ignore terms higher or equal to $\mathcal{O}(\varepsilon^{2})$. The term $\varepsilon$ is constrained experimentally with values smaller than $10^{-3}$ Abel:2008ai .

The interaction between gauge and scalar fields is

[TABLE]

where the covariant derivative $D_{\mu}$ for neutral gauge bosons is defined as

[TABLE]

where $g_{x}$ and $Q^{\prime}_{i}$ are the coupling constants and the charge for $U(1)_{X}$, respectively. When the SSB is achieved not only the mass terms are generated but also mixing terms are obtained:

[TABLE]

where

[TABLE]

and

[TABLE]

meanwhile, the $Z$ gauge boson mass retains the same value set by the SM,

[TABLE]

In order to cancel the mixing term, the following rotation is required

[TABLE]

where the mixing angle $\xi$ satisfies the expression $\tan 2\xi=\frac{2\Delta^{2}}{m^{2}_{Z^{0}}-m^{2}_{Z^{{}^{\prime}0}}}$, and has been constrained to the interval $|\xi|<10^{-3}$ Bouchiat:2004sp .

2.5 Fermion interactions

The most general Yukawa Lagrangian is

[TABLE]

where $Y_{a}^{0f}$ are the $3\times 3$ Yukawa matrices, for $f=u,d,l$. $q_{L}$ and $l_{L}$ denote the left-handed fermion doublets under $SU(2)_{L}$, while $u_{R}$, $d_{R}$, $l_{R}$ correspond to the right-handed singlets. The zero superscript in fermion fields stands for the interaction basis. The DM stability is lost if the couplings $Y_{2ij}^{0f}$ appear in the Eq. (40). These Yukawa couplings can be eliminated by the correct assignment of values for charges under the $Z_{2}$ and $U(1)_{X}$ symmetries, as previously done.

In the case of discrete $Z_{2}$ symmetry with $x_{1}=x_{2}$, the couplings $Y_{2ij}^{0f}$ must be equal to zero in order to respect the discrete $Z_{2}$ symmetry. The couplings $Y_{1ij}^{0f}$ are allowed if the assignment of the $U(1)_{X}$ charges for the fermions satisfy

[TABLE]

and

[TABLE]

where $x_{q,l}$ are the $U(1)_{X}$ charges of left-handed doublet fermions, meanwhile, $x_{u,d,e}$ are the $U(1)_{X}$ charges of right-handed fermions.

In the case of $x_{1}\neq x_{2}$ we set the $U(1)_{X}$ charges such that $\mp x_{2}-x_{q}+x_{u,d}\neq 0$ and $x_{2}-x_{l}+x_{e}\neq 0$ in order to eliminate the couplings $Y_{2ij}^{0f}$ in Eq. (40). Obviously, $\Phi_{1}$ also satisfies Eqs. (41) and (42) to provide the masses of the fermions as in SM. Feynman rules of IDMS are shown in table 1.

It is important to mention that the fermion charges under $U(1)_{X}$ must satisfy the triangle anomaly equations, which can be reviewed in Mantilla:2016lui , in order to garantize an anomaly free model. The anomaly cancellation requirements for fermion charges $x_{f}$, for $f=q,u,d,l,e$, are shown in table 2 as a function of $x_{q}$.

3 Constraints on free model parameters

In this section, we obtain the experimentally allowed regions for the free model parameters involved in our analysis by considering the most up-to-date experimental collider results reported by CMS Sirunyan:2018koj and ATLAS ATLAS:2018doi collaborations, namely, signal strengths, denoted by $\mathcal{R}_{x\bar{x}}$. In this work we consider the production of $H_{i}$ via gluon fusion and we use the narrow width approximation. Then, $\mathcal{R}_{x\bar{x}}$ can be written as follows:

[TABLE]

where $\Gamma(H_{i}\to gg)$ is the decay width of $H_{i}$ into gluon pair, with $H_{i}=h^{\text{IDMS}}$ and $h^{\text{SM}}$. Here $h^{\text{IDMS}}$ is the SM-like Higgs boson coming from IDMS and $h^{\text{SM}}$ is the SM Higgs boson; $\mathcal{B}(H_{i}\to x\bar{x})$ is the branching ratio of $H_{i}$ decaying into a $x\bar{x}$, where $x\bar{x}=b\bar{b},\;\tau^{-}\tau^{+},\;\mu^{-}\mu^{+},\;WW^{*},\;ZZ^{*},\;\gamma\gamma$. Besides to measurements of colliders, we use the most-up-date upper limit on WIMP-nucleon cross-section, for the spin independent case, reported by XENON1T collaboration Aprile:2018dbl and whose value for a DM candidate mass of $30$ GeV is given by:

[TABLE]

On the other side, the free parameters of the IDMS involved in our analysis are the following:

• •

Mixing angle $\alpha_{1}$.

• •

Vacuum Expectation Value of the scalar singlet, $\upsilon_{x}$.

• •

$U(1)_{X}$ coupling constant, $g_{x}$.

• •

$Z^{\prime}$ gauge boson mass, $m_{Z^{\prime}}$.

• •

Scalar mass, $m_{S}$.

• •

Charged scalar boson mass, $m_{H^{\pm}}$.

• •

Dark matter boson mass, $m_{\chi}$.

In order to constrain the ${Z^{\prime}}$ gauge boson mass, $m_{Z^{\prime}}$, the upper limit on the production cross-section of a $Z^{\prime}$ gauge boson times the branching ratio of the $Z^{\prime}$ decaying into $\ell^{-}\ell^{+}$ Aaboud:2017buh , with $\ell=e,\,\mu$, was considered.

3.1 Constraint on mixing angle $\alpha_{1}$

Due to the coupling $g_{hPP}^{\text{IDMS}}=\cos{\alpha_{1}}\cdot g_{hPP}^{\text{SM}}$, with $P=f_{i},W$, allowed regions for $\cos\alpha_{1}=c_{\alpha_{1}}$ can be extracted experimentally from $\mathcal{R}_{x\bar{x}}$. We find that $\mathcal{R}_{WW^{*}}$ is the most stringent way of limiting $c_{\alpha_{1}}$. In the fig. 1 we show the $c_{\alpha_{1}}-\mathcal{R}_{WW^{*}}$ plane, where the dark area (orange online) is the allowed region for $\mathcal{R}_{WW^{*}}$ at $2\sigma$. The graph was generated via SpaceMath SpaceMath .

We note that the allowed interval for $c_{\alpha_{1}}$ is between $\sim 0.99-1$. This is to be expected since $c_{\alpha_{1}}$ must be closed to the unit to have small deviations of the SM couplings. In particular, when $c_{\alpha_{1}}=1$ the SM is recovered. From now on we will consider $c_{\alpha_{1}}=0.99$.

3.2 Constraint on the $Z^{\prime}$ gauge boson mass $m_{Z^{\prime}}$

In order to constrain the $Z^{\prime}$ gauge boson mass, we now turn to analyze the $Z^{\prime}$ production cross-section times the branching ratio of $Z^{\prime}$ decaying into $\ell^{-}\ell^{+}$ ($\sigma_{Z^{\prime}}\mathcal{B}_{Z^{\prime}}$), with $\ell=e,\,\mu$. The ATLAS and CMS collaborations Aaboud:2017buh , Sirunyan:2018exx searched for a new resonant and non-resonant high-mass phenomena in dilepton final states at $\sqrt{s}=13$ TeV with an integrated luminosity of 36.1 fb*-1* and 36 fb*-1*, respectively. Nevertheless no significant deviation from the SM prediction was observed. Lower limits excluded on the resonant mass was reported depending on specific models.

Figure 2 shows $\sigma_{Z^{\prime}}\mathcal{B}_{Z^{\prime}}$ as a function of the $Z^{\prime}$ gauge boson mass for $g_{x}=0.4,\,0.5$ and $2m_{Z}/\upsilon.$ The last value is related to the coupling of $Z$ gauge boson to fermions. We present two regions, the largest (magenta online) represents the results reported by ATLAS collaboration for a center-of-mass energy of $\sqrt{s}=13$ TeV and 36.1 fb*-1* as mentioned above, while the smallest (yellow online) area corresponds to $\sqrt{s}=14$ TeV and 3000 fb*-1*, which is the goal of the High Luminosity Large Hadron Collider ATL-PHYS-PUB-2018-044 ; these analyses are based on generator-level information with parameterized estimates applied to the final state particles to simulate the response of the upgraded ATLAS detector and pile-up collisions.

Considering the results reported by LHC (HL-LHC),

$\sigma_{Z^{\prime}}\mathcal{B}_{Z^{\prime}}>10^{-4}$ pb ($\sim 10^{-6}$ pb) excludes $m_{Z^{\prime}}\lesssim 3$ TeV ($m_{Z^{\prime}}\lesssim 5$ TeV) for $g_{x}=0.4$, while $m_{Z^{\prime}}\lesssim 3.4$ TeV ($m_{Z^{\prime}}\lesssim 5.4$ TeV) for $g_{x}=0.5$ are excluded. Finally, we explored the case in which $g_{x}=g_{Z}=2m_{Z}/\upsilon$ and we observe a similar behavior as reported in the refs. Aaboud:2017buh -ATL-PHYS-PUB-2018-044 , excluding $m_{Z^{\prime}}\lesssim 4.5$ TeV ($m_{Z^{\prime}}\lesssim 6.5$ TeV) .

3.3 Constraint on $\upsilon_{x}$, $g_{x}$

In the fig. 3 we show the $g_{x}-\upsilon_{x}$ plane, in which allowed regions for $\mathcal{R}_{ZZ^{*}}$ and the upper limit on WIMP-nucleon cross section, $\sigma^{SI}(\chi N\to\chi N)$, are displayed. We generate the Feynman rules of the IDMS via LanHEP lanhep and we evaluate $\sigma^{SI}(\chi N\to\chi N)$ through CalcHep calchep . We observe that the consistent zone with both $\mathcal{R}_{ZZ^{*}}$ and $\sigma^{SI}(\chi N\to\chi N)$ allows values for $\upsilon_{x}$ in the interval from $\sim 11$ to $\sim 16$ TeV for $g_{x}=0.4$, while the white area represents the excluded region.

3.4 Constraint on the charged scalar boson mass $m_{H^{\pm}}$

We use $\mathcal{R}_{\gamma\gamma}$ in order to constrain the charged scalar boson mass $m_{H^{\pm}}$. In addition to the SM contributions, the $h\to\gamma\gamma$ decay receives contributions at one-loop level of charged scalar bosons predicted by the IDMS. The $b\to s\gamma$ decay is another process that can also impose strong restrictions on the charged scalar boson mass. However, since this particle arises from the inert doublet, its Yukawa couplings with fermions are absent, so the $b\to s\gamma$ decay is not a way to restrict the charged scalar boson mass.

Figure 4 shows the $m_{H^{\pm}}-\mathcal{R}_{\gamma\gamma}$ plane. We note that $\mathcal{R}_{\gamma\gamma}$ imposes a lower bound on the charged scalar boson mass as $330$ GeV $\lesssim m_{H^{\pm}}$ at 1$\sigma$ with $\lambda_{2x}=0$. Figure 4 also shows bounds for $\lambda_{2x}=0.005$, which are less restrictive than the previous case. We find that $170$ GeV $\lesssim m_{H^{\pm}}$ ($75$ GeV $\lesssim m_{H^{\pm}}$) at 1$\sigma$ (2$\sigma$), respectively. Finally, the total decay width of the Higgs boson Tanabashi:2018oca excludes $60$ GeV $\lesssim m_{H^{\pm}}$. The values for $\lambda_{2x}$ parameter are select such that they are compatibles with viable values for relic density within the framework of the IDMS. We generated random values for $\lambda_{3}$ between $0.01-0.0105$ ($0.0297-0.03$) for $\lambda_{2x}=0.005$ ($\lambda_{2x}=0$), respectively, for the same reason.

In table 3 we present a summary of the values for the model parameters used in our following analysis.

4 Phenomenology for $S$ and $Z^{\prime}$

We now analyze the behavior of the branching ratio of the dominant decay channels of the particles coming from complex singlet, namely, scalar $S$ and the $Z^{\prime}$ gauge boson. Analytical formulas of the partial decay widths are presented in A. Figure 5 shows the branching ratios for relevant decays of the $S$ neutral scalar at tree and one-loop level. In Figure 6 the relevant decay channels are also presented but for the $Z^{\prime}$ gauge boson predicted by the IDMS.

We observe that the dominant $S$ decay modes are $S\to VV$, with $V=W,\,Z$, and $S\to hh$. These processes are of the order of $10^{-1}$. Once the $t\bar{t}$ channel became open for $m_{S}\geq 2m_{t}$, its branching ratio is of the order of $S\to VV$ channel up to a $S$ mass of about $800$ GeV. Later when $m_{S}$ increases, the value of the $BR(S\to t\bar{t})$ decreases such that $BR(S\to t\bar{t})\sim\mathcal{O}(10^{-2})$ for $m_{S}=2$ TeV. Another relevant decay mode is $S\to b\bar{b}$ whose branching ratio range decreases from $10^{-3}$ to $10^{-5}$. At one-loop level, the dominant channel is $S\to gg$ with a branching ratio up to $\mathcal{O}(10^{-4})$ for $m_{S}=200$ GeV and $\mathcal{O}(10^{-5})$ for $m_{S}=2$ TeV.

As far as the $Z^{\prime}$ gauge boson is concerned, their dominant decay modes are into type-up quarks, whose sum is about $3.5\times 10^{-1}$, followed by type-down quarks and finally by charged leptons with a branching ratio of the order $10^{-1}$.

5 Relic density

Once the model parameters are bounded by experimental and theoretical constraints, we now turn to analyze if the model can help us to understand the relic density, which is the current experimental quantity of DM particle that remains after of freeze-out process. The observed value for non-baryonic matter reported by PLANCK collaboration PLANCK is

[TABLE]

where $h$ is the Hubble constant in units of $100\frac{km}{s.Mpc}$. Our analysis is based on selecting $\chi$ as DM candidate. The relic density is obtained by solving the Boltzmann equation for the number density rate which is given by

[TABLE]

where $n$ is the DM number density and $a$ is a scale factor. All information about the model is contained in the thermally averaged cross section $\langle\sigma v\rangle$. The relic density, $\Omega h^{2}$, is obtained by using the micrOmegas package micromegas ; micromegas1 , which requires all information about the IDMS for which we implement the model via the LanHep package lanhep .

In Fig. 7 we present a scattering plot of the relic density as a function of the DM candidate mass, $m_{\chi}$. We show three scenarios to note the sensitivity of the relic density on $\lambda_{2x}$. These scenarios are classified by their random value intervals:

• •

$R_{1}:\lambda_{2x}\sim\mathcal{O}(10^{-3}-10^{-2})$,

• •

$R_{2}:\lambda_{2x}\sim\mathcal{O}(10^{-4}-10^{-3})$,

• •

$R_{3}:\lambda_{2x}=0\cup\mathcal{O}(10^{-7}-10^{-3})$.

In all cases, we use random values for $\lambda_{345}\sim\mathcal{O}(10^{-2}-10^{-1})$ and values for $m_{\chi}$ from $1$ to $3000$ GeV. The most favored scenario is $R_{3}$, which contains the special case $\lambda_{2x}=0$. When this occurs, the IDM $h\chi\chi$ coupling is recovered, as the Table 1 shows. Nevertheless, in the IDMS two new portals contribute to relic density, namely, $Z^{\prime}$ gauge boson and the neutral scalar $S$; both arise from the complex singlet $\mathcal{S}_{X}$. We observe that, depending on $\lambda_{345}$ and $\lambda_{2x}$, masses from a few GeV to about 2 TeV are in agreement with the results for the relic density of the PLANCK collaboration PLANCK . It is worth mentioning that we include a pair of photons in the final state in the process of annihilation of the DM particles, i.e., $\chi\chi\to h\to\gamma\gamma$.

Figure 8 shows the representatives values of $\lambda_{345}$ in the interval $0.0297-0.03$ ($0.01-0.0105$) and $\lambda_{2x}=0$ ($\lambda_{2x}=0.005$), respectively. The scattering process $\sigma^{SI}(\chi N\to\chi N)$ excludes an important region of allowed values of $g_{x}$ and $\upsilon_{x}$ for $\lambda_{345}>0.0105$ ($\lambda_{345}>0.03$) by assuming $\lambda_{2x}=0$ ($\lambda_{2x}=0.005$). Under these considerations we find intervals for the masses of DM candidates:

For $\lambda_{2x}=0:$

• •

Light masses: $1\lesssim m_{\chi}\lesssim 80$ GeV.

• •

Heavy masses: $1100\lesssim m_{\chi}\lesssim 1600$ GeV. 2. 2.

For $\lambda_{2x}=0.005:$

• •

Light and intermediate masses: $1\lesssim m_{\chi}\lesssim 700$ GeV.

• •

Heavy masses: $1800\lesssim m_{\chi}\lesssim 2100$ GeV.

6 Dark Matter production at hadron colliders

The ATLAS and CMS collaborations Aaboud:2017dor searched for the reaction $pp\to\chi\chi\gamma$ through events that contain an energetic photon and large missing transverse momentum, corresponding to an integrated luminosity of 36.1 fb*-1* at centre-of-mass energy of 13 TeV. However, only the exclusion limits were reported. In order to motivate a potential and sophisticated study of the production of DM particles at hadron colliders, we evaluate the $pp\to\chi\chi\gamma$ production cross section and the main SM background processes via MadGraph5 Alwall:2011uj . Our study is focused on future hadron colliders, namely:

• •

High-Luminosity Large Hadron Collider Apollinari:2017cqg (HL-LHC). The HL-LHC is a new stage of the LHC starting about 2026 with a center-of-mass energy of 14 TeV. The upgrade aims at increasing the integrated luminosity by a factor of ten ($\sim$3000 fb*-1*) with respect to the final stage of the LHC ($300$ fb*-1*).

• •

High-Energy Large Hadron Collider Benedikt:2018ofy (HE-LHC). The HE-LHC is a possible future project at CERN. The HE-LHC will be a 27 TeV $pp$ collider being developed for the 100 TeV Future Circular Collider. This project is designed to reach up to 12000 fb*-1* which opens a large window for new physics research.

• •

Future Circular hadron-hadron Collider Arkani-Hamed:2015vfh (FCC-hh). The FCC-hh is a future 100 TeV $pp$ hadron collider which will be able to discover rare processes, new interactions up to masses of around 30 TeV and search for a possible substructure of the quarks. The FCC-hh will reach up to an integrated luminosity of 30000 fb*-1* in its final stage.

6.1 Signal and background events

The main SM background to the $\gamma+E_{T}^{miss}$ final state are events containing either a true photon or an object misidentified as a photon. The dominant background processes are the electroweak production of $Z(\to\nu\nu)\gamma$, $W(\to\ell\nu)\gamma$ and $Z(\to\ell\ell)\gamma$ with unidentified charged leptons, $e,\,\mu$, or with $\tau\to$hadrons+$\nu_{\tau}$.

As far as our computation scheme is concerned, we first use the LanHEP lanhep routines to obtain the IDMS Feynman rules for MadGraph5 Alwall:2011uj . Secondly, we evaluated the production cross section of the signal and background processes (PCSS and PCSB) and we generated $10^{5}$ events for both reactions.

In fig. 9, we present the PCSS and PCSB (axis left) and number of events (axis right) for the different future hadron colliders, i.e., 9(a) HL-LHC, 9(b) HE-LHC and finally 9(c) for the FCC-hh, with integrated luminosities 3000 fb*-1*, 12000 fb*-1*, 30000 fb*-1*, respectively. In all graphics, horizontal lines represent the potential SM background processes. We observe that light masses for the DM candidate are favored producing up to about $10^{5}$ ($10^{6}$, $10^{7}$) events at the HL-LHC (HE-LHC, FCC-hh) by considering a DM mass of 10 GeV. However, the intermediate regimen of masses ($\sim$ 500 GeV) is disadvantaged by this channel, even at the FCC-hh only one event will be produced. Therefore, we analyze the range of masses 10-100 GeV for event reconstruction.

6.2 Event reconstruction

We closely follow the strategy by ATLAS Collaboration Aaboud:2017dor in which the photon identification is based on energy deposited at the electromagnetic calorimeter. Candidate photons are required to have $E_{T}^{\gamma}>150$ GeV, to be within $|\eta|<1.37$ and be isolated by demanding energy in the calorimeter in a cone of size $\Delta R=\sqrt{(\Delta\eta)^{2}+(\Delta\phi)^{2}}=0.4$. Due to the elusive nature of the DM candidates, these particles are characterized by missed energy transverse and therefore we demand for $E_{T}^{\text{miss}}$> 150 GeV. It is also required that the photon and the $E_{T}^{\text{miss}}$ do not overlap in the azimuthal plane, then is required the condition $\Delta\phi(\gamma,\,E_{T}^{\text{miss}})>0.4$. In our analysis, the above requirements work well for intermediate masses, 100 GeV $\lesssim m_{\chi}$. In the analysis performed in Section 6.1, masses in the interval of 10-100 GeV (10-300 GeV), for HL-LHC and HE-LHC (for FCC-hh), respectively, are favored. Therefore, for light DM masses we apply slightly different cuts, namely, 10<$E_{T}^{\text{miss}}$<150 GeV and 10<$E_{T}^{\gamma}$<150 GeV. In Fig. 10 we present the $E_{T}^{\text{miss}}$ distribution to both signal and background processes, while in figure 11 we show the photon transverse energy.

We observe that $E_{T}^{\gamma}$ and $E_{T}^{\text{miss}}$ grow as $m_{\chi}$ increase. For a better illustration, Fig. 12 shows the normalized $E_{T}^{\gamma}$ and $E_{T}^{\text{miss}}$ distributions for $m_{\chi}=$10, 100 and 500 GeV.

6.3 Signal significance

We compute the signal significance defined as $\textbf{S}=\frac{N_{S}}{\sqrt{N_{S}+N_{B}}}$, where $N_{S}$ ($N_{B}$) are the number of signal (background) events after the kinematic cuts were applied. For colliders considered (HL-LHC, HE-LHC and FCC-hh), we find that through $pp\to\chi\chi\gamma$ production only at the FCC-hh will be possible to claim detection of the DM candidate in the range of masses 10-60 GeV once a center-of-mass energy of 100 TeV and an integrated luminosity about 22000 fb*-1* are reached. This is illustrated in Fig. 13 which shows the signal significance as a function of the DM candidate mass.

7 Conclusions

In this work, we study an extension of the SM with $U(1)_{X}$ gauge symmetry that includes two doublets and one complex singlet scalar field in order to introduce a WIMP as DM candidate. The proposed candidate as DM in this extension arises from one inert doublet, whose VEV is equal to zero. In order to ensure the stability of the DM candidate, we consider two scenarios to control the scalar couplings: a discrete $Z_{2}$ symmetry and the $U(1)_{X}$ symmetry.

In the constrained IDMS Bonilla:2014xba , the parameters associated with the singlet cubic terms in the scalar potential are responsible for the source of CP violation, since three neutral scalars are mixed to generate the physical states. In the IDMS with local gauge $U(1)_{X}$ symmetry this type of mixture, which produces an undefined state of CP for neutral scalars, can also occur when $\lambda_{6}\neq 0$ and $\lambda_{12x}\neq 0$. However, this analysis is out of the objective at the moment. Then, the study of explicit CP violation can be considered with this model.

We explore the allowed regions for free model parameters of the IDMS taking into account the most up-to-date experimental collider and astrophysical results. We find that the signal strength $\mathcal{R}_{WW^{*}}$ is the most stringent, allowing an interval for the neutral scalar mixing angle such that $0.99\lesssim\cos{\alpha_{1}}\lesssim 1$. The analysis of the $\sigma(pp\to Z^{\prime})$ production cross-section times $\mathcal{B}(Z^{\prime}\to\ell^{-}\ell^{+})$, with $\ell=e,\,\mu$, excludes regions for $m_{Z^{\prime}}\lesssim 3$ TeV with $g_{x}=0.4$.

Regions for the masses of DM candidate in the order of light ($\mathcal{O}(10)$ GeV), intermediate ($\mathcal{O}(100)$ GeV), and heavy ($\mathcal{O}(2)$ TeV) are in agreement with the upper limit on $\sigma^{SI}(\chi N\to\chi N)$ and relic density reported by XENON1T and PLANCK collaborations, depending mainly on $\lambda_{345}$ and $\lambda_{2x}$. We find that the allowed interval for the DM candidate mass is highly sensitive to $\lambda_{345}$ and $\lambda_{2x}$. For instance, for the values of $\lambda_{2x}=0$ and $\lambda_{345}=0.03$, the allowed values for DM mass are obtained such that $1.1$ TeV$\lesssim m_{\chi}\lesssim 1.6$ TeV meanwhile for $\lambda_{2x}=0.005$ and $\lambda_{345}=0.01$ the result is $m_{\chi}\lesssim 0.7$ TeV. Additionally, IDMS presents an improvement in the mass region due to the portals associated with the $Z^{\prime}$ gauge boson and a scalar boson $S$, both portals are predicted by the IDMS, which are absent in models as IDM.

We conclude that the IDMS is a viable model for the study of DM which provides an improvement in the allowed regions of the DM candidate mass. The IDMS has a rich phenomenology through processes involving $Z^{\prime}$, $S$ and $H^{\pm}$ bosons that could be tested at hadron colliders. Also, the IDMS predict DM particle masses in the interval 10-60 GeV that could be detectable at the FCC-hh through the $pp\to\chi\chi\gamma$ process. However, other processes could also be analyzed to complement the search for DM particles. On the other hand, we also find restrictions for $\lambda_{4,5}$ parameters of the model that prohibit the decay $\chi\to W^{\pm}H^{\mp}$.

Acknowledgements.

M. A. Arroyo-Ureña especially thanks to PROGRAMA DE BECAS POSDOCTORALES DGAPA-UNAM for postdoctoral funding and thankfully acknowledge computer resources, technical advise and support provided by Laboratorio Nacional de Supercómputo del Sureste de México. This work was supported by projects Programa de Apoyo a Proyectos de Investigación e Innovación Tecnológica (PAPIIT) with registration codes IA107118 and IN115319 in Dirección General de Asuntos de Personal Académico de Universidad Nacional Autónoma de México (DGAPA-UNAM), and Programa Interno de Apoyo para Proyectos de Investigación (PIAPI) with registration code PIAPIVC07 in FES-Cuautitlán UNAM and Sistema Nacional de Investigadores (SNI) of the Consejo Nacional de Ciencia y Tecnología (CONACYT) in México.

Appendix A Decay widths of scalar $S$ and pseudoscalar $A$ bosons

A.1 Scalar boson decays

The most relevant decays of both $CP$-even and $CP$-odd scalar bosons have been long studied in the literature. We will present the decay width formulas for the sake of completeness. The tree-level two-body widths are given as follows:

[TABLE]

with $\tau_{i}=4m_{i}^{2}/m_{S}^{2}$ and $N_{c}$ is the color number. From here we easily obtain the flavor conserving decay width. The $CP$-even scalar boson decays into pairs of real electroweak gauge bosons can also be kinematically allowed. The corresponding decay width is

[TABLE]

with $n_{V}=1\;(2)$ for $V=W\;(Z)$ and $g_{S\bar{f}_{i}f_{j}}$-$g_{SVV}$ given in the table 1.

Additionally to the tree level decays, other relevant channels arise at one-loop, such as $S\to\gamma\gamma$ and $S\to gg$, whose decay widths are given by:

[TABLE]

with the subscript $s$ standing for the spin of the charged particle circulating into the loop. The $A_{s}^{S\gamma\gamma}$ function is given by

[TABLE]

where

[TABLE]

The two-gluon decay can only receive contributions from quarks and its decay width can be obtained from (49) by only summing over quarks and making the replacements $\alpha^{2}\to 2\alpha^{2}_{S}$, $N_{c}Q_{f}^{2}\to 1$.