Abstract

We study the Inert Higgs Doublet Model and its inert scalar Higgs as the only source for dark matter. It is found that three mass regions of the inert scalar Higgs can give the correct dark matter relic density. The low mass region (between 3 and 50 GeV) is ruled out. New direct dark matter detection experiments will probe the intermediate (between 60 and 100 GeV) and high (heavier than 550 GeV) mass regions. Collider experiments are advised to search for decay in the two jets plus missing energy channel.

1. Introduction

Astrophysical observations provide strong evidence for the existence of dark matter (DM) [1] and its abundance in the current phase of the Universe [2]. According to the newest results from Planck collaboration, there is 74% approximately of matter which is not directly visible but is observed due to its gravitational effects on visible matter. Additional evidence for the existence of DM comes from study of the rotation curves of spiral galaxies [3], the analysis of the bullet cluster [4], and the study of baryon acoustic oscillations [5]. One of the most common hypotheses used to explain these phenomena is to postulate the existence of weakly interacting massive particles (WIMPs) [6].

In order to explain DM, we study a simple extension of the Standard Model (SM), called the Inert Higgs Doublet Model (IHDM). This model, which was originally proposed for studies on electroweak (EW) symmetry breaking [7], introduces an additional doublet and a discrete symmetry. These two characteristics of the model modify the SM phenomenology, but there are some regions of the parameter space which predict only small deviations from the SM. Nevertheless, one of the most attractive characteristics of the IHDM is the presence of a stable neutral particle which can be a DM candidate.

In this work the IHDM is revisited, considering various restrictions, focusing our analysis on the zones of the parameter space which reproduce the correct DM relic density according to the newest measurements [8, 9]. In this context, three not connected mass regimes for the lightest inert particle are found. These regimes are also analyzed using LHC observables like branching ratios to invisible particles and a specific SM-like Higgs boson decay mode. Additionally, we study some inert decays modes. Finally, we use a direct detection approach to rule out one of the mass regimes. It is further shown that this regime can also be ruled out by constraints from collider physics [10].

This paper is organized as follows. After a short introduction in Section 1 the model is introduced in Section 2 by formulating the associated potential and constraints of the model, and by exploring the parameter space and its characteristics in Section 3. In Section 4 the behavior of the model is presented from a collider physics perspective, studying the modifications of the SM and its implications for the IHDM. In Section 5 the results of this study are complemented by an analysis from the dark matter perspective. Finally, in Section 6 we remark on the most important conclusions of our work.

2. The Inert Higgs Doublet Model

Consider an extension of Standard Model (SM), which contains two Higgs doublets and a discrete symmetry [7] ( and ). All fields of the SM are invariants under the discrete symmetry, and is completely analogous to the SM Higgs doublet.

The most general renormalizable invariant Higgs potential that also preserves the discrete symmetry iswhere , , are given byThe parameters and are intrinsically real, and will be assumed to be real [11]. After Spontaneous Symmetry Breaking, the vacuum expectation values of the Higgs doublets arewhere is forced by the discrete symmetry and  GeV. Expanding the fields around those vacua, we definewhere and are the neutral and charged Goldstone bosons and is the SM-like Higgs boson. The fields in the second doublet belong to the so-called dark or inert sector. They are the scalar and pseudoscalar , both neutral, and the charged scalar . As in the SM, the parameter is related to by the tree-level tadpole condition . The masses of physical states [12, 13] arewithAs independent and free parameters we take the masses , , and at tree-level and the couplings and . The SM-like Higgs boson mass is fixed now thanks to the measurement  GeV [14, 15].

The constraints we initially impose include vacuum stability at tree-level, where constraints on the couplings and mass terms appear [1618]; perturbation () [19, 20] and unitarity [21], where we impose that the scalar potential is unitary and that several scattering processes between scalar and gauge bosons are bounded; electroweak precision tests through the , , and parameters [22] applied to the IHDM [21, 23], with values given bywith a correlation coefficient of [24]; and collider constraints [16, 2527], where we satisfy lower bonds on the Higgs boson masses.

The DM particle must be neutral. In our analysis we assume it is the boson; thus and , which due to (5) translates to

We do not consider as the DM candidate because it is analogous to consider as the DM candidate defining instead of .

3. IHDM Parameter Space

We randomly scan the parameter space of the IHDM, taking into account all the constraints mentioned in the previous section. Additionally, we compute some astrophysical properties of the model using the micrOMEGAs software [28]. We consider masses satisfying 1 GeV 1 TeV, where . In addition, we consider cosmological measurements: the DM relic density is a property related to its abundance in the current phase of the Universe. This quantity is well measured by WMAP [29] and Planck [30] experiments. Following [31] to combine both measurements we obtain

In Figure 1 the coupling is shown as a function of the Higgs boson mass varying (, , , , and ). We work with the hypothesis that the light inert Higgs boson is providing the complete DM density given in (9). The color code is as follows: red points (dark gray) produce a relic density above the limit given in (9); blue points (black) produce a relic density within the region; green points (light gray) produce a relic density below the limit. Regarding the points that satisfy the relic density we see three clear regions [12, 32], one for low ( GeV), another one for medium ( GeV), and finally one for high values of ( GeV). The explanation for the gap is related to annihilation processes and it will be given later. At  GeV the IHDM can no longer be compatible with vacuum existence and stability [33, 34].

In Figure 2 we have for the same scan and color code the mass of the heavy pseudoscalar (a) and the mass of the charged Higgs (b) as a function of the mass of the DM candidate . Due to dedicated pre-LHC collider searches, a bound that captures most of the features is  GeV and  GeV. This is so with the exception of a small strip for  GeV seen in Figure 2(a), where due to the proximity of the and masses the search loses sensitivity. In this figure the gap in values of when the relic density is imposed is apparent. Notice that the density of solutions is larger when the masses for , , and are close to each other and that this feature is more pronounced when the masses of these particles are near the TeV scale (due to the logarithmic scale). Finally we notice that is more sensitive to the parameters and (Figure 1) than the masses of the other inert particles (Figure 2).

4. Collider Physics

As we mentioned before, the Higgs boson discovered at CERN in 2012 is the SM-like Higgs boson of our model from the non-Inert Higgs Doublet field . This particle couples to the charged Higgs pair (), which contributes to the diphoton decay width (in this minimal scenario, the IHDM cannot account for the reported excess of diphoton events by ATLAS [35] and CMS [36] collaborations in their Run-II 13 TeV analyses, because the symmetry prevents the extra Higgs bosons of the model from decaying into just two photons) [37]. For the same reason also contribute to .

It is convenient to work with the parameter [21, 33, 38]The value we use for the SM is MeV [39]. ATLAS [40] and CMS [41] collaborations have studied this decay mode, and if we combined both results [31] we obtain .

In Figure 3 we have the parameter as a function of the DM candidate mass (a) and as a function of the coupling (b). The points in parameter space that produce a correct relic density can be divided into three groups. In the case of very light masses for the DM candidate ( GeV approximately) the decay mode is open, and is close to zero ruling those masses out [42]. In the intermediate mass case, approximately between 60 and 100 GeV, there is a region with acceptable solutions characterized by . This region is characterized by increasingly heavier values for and . In the large mass region ( GeV approximately), the charged Higgs gives a negligible contribution to the decay such that is close to unity. Interestingly, Figure 3(b) shows that perturbative values () are preferred.

In addition, if the inert particles are light enough, there are two other two-body decays which areThere is no phase space for a two-body decay . It is possible to define the parameter , in analogy to defined in (10). It is interesting that even though the decay is not well measured, it still can give additional insight to the model [43, 44]. As Figure 4 shows, there appears a very narrow correlation between and , which is a common feature for versus plots [17, 18, 45]. The bisector branch only contains points that satisfy (inert invisible decay channel open). Points of parameter space which satisfy the relic density and have low DM candidate mass are ruled out, because they produce a very small value for . By analyzing the characteristics of those data points one finds that the larger branch includes only points with , and the ones that also satisfy relic density are close to , as was mentioned before. The two branches seen in the versus relation also appear in the nonnormalized versus relation (not shown). But it is reduced only to the long (green) branch in the versus relation. Further one sees that when the decay channel is closed, the loop transforms into the long branch, the otherwise SM dot (the intersection point between the two branches). If the decay channel is open, the second branch appears because the channel tends to dominate [46].

We are also interested in the invisible decay of the SM-like Higgs boson. If the DM candidate mass satisfies , the two-body decay channel is open, which is invisible for the LHC detectors and shows only as missing momentum. There are measurements for the invisible decay of the SM-like Higgs from the LHC experiments. Taking a simple average of the upper bounds to the invisible decay rate from ATLAS [47] and CMS [48] gives for the SM-like Higgs boson. In Figure 5 we show the branching ratio for the invisible decay of the SM-like Higgs boson , as a function of the mass of the DM candidate (a) and as a function of the parameter (b). In (a) we also have a horizontal line that shows the upper bound for mentioned above. The threshold appears clearly in (a). Most of the points with correct relic density satisfying  GeV are ruled out because they produce a very large invisible branching ratio for . On the contrary, most of the points with higher DM candidate mass are fine because they produce an invisible branching ratio equal to zero. In (b), where we have as dashed lines the bounds from LHC experiments, we see a very strong relation between and . The points that satisfy relic density with a low mass for the DM candidate are simultaneously excluded from and from . The rest of the points satisfy the bounds.

To finalize this section we discuss the branching ratios for the two observable inert Higgs bosons and . In Figure 6 we have the branching ratios for the pseudoscalar Higgs boson (a) as a function of its mass and the branching ratios for the charged Higgs boson (b) as a function of its mass. In (a) we show the decays of the inert pseudoscalar Higgs , which are (solid line) and (dashed line). As can be seen from the Feynman rules, the only unknown parameters the branching ratios depend on are the masses , , and . Therefore, four scenarios are considered: (i) , (ii) , (iii) , and (iv)  GeV. The gauge boson can be off shell, although we consider the inert Higgs bosons always on shell. The oscillation near the threshold is due to different increasing rates for the decay rates when the gauge boson is off shell. The branching ratio is always large (because ) while can be low near thresholds. There is a crossing point where . In (b) we show the decays of the charged Higgs . Analogous scenarios are considered, but replacing by . In solid line we have the branching ratio for the decay and in dash we have . In the case of , there is no crossing point; thus is always larger than . We remind the reader that the presence of a Higgs boson in a final state is seen as missing momentum at the LHC.

5. Cosmology and Dark Matter

The existence of dark matter seems to be well established now [49]. There are several candidates for DM; among them are the previously mentioned WIMPs. A good particle candidate for DM must be neutral and stable (or quasi-stable). The discrete symmetry in the model studied in this paper ensures that the lightest of the inert Higgs bosons is stable. Observation implies it is either or (or in a fine tuned scenario both). In this paper we study the former case. An important restriction this candidate must satisfy is that its mass density must agree with experimental observations. We calculate the relic density of our DM candidate using micrOMEGAs software [28]. To better understand the results on the relic density, we calculate also the thermal averaged annihilation cross section times the relative velocity , or annihilation cross section for short.

In Figure 7 we plot the relic density as a function of the DM candidate mass (a), and the annihilation cross section also as a function of the DM candidate mass (b), calculated with [28]. For the relic density case, we also show, as an horizontal dashed line, the experimentally measured value for as given in (9). The scan shows a large distribution with differences that can reach more than 10 orders of magnitude. For this reason most of the points in the scan are ruled out if one demands that is actually the only WIMP responsible for the observed DM signatures. There are two mass gaps that divide the mass region in three: low mass ( GeV approximately), medium mass ( GeV approximately), and high mass ( approximately). The origin of these mass gaps is better understood with the aid of the right frame. In the right frame of Figure 7 we show the annihilation cross section as a function of the DM candidate mass , with vertical lines denoting different thresholds. The first gap is near the threshold where the annihilation channel becomes very efficient due to the fact that the SM-like Higgs is on-shell. The second gap starts at the thresholds and , where become available, and continues later with the threshold where the channel opens up. The annihilation channel also helps. All these new annihilation channels make the DM annihilation very efficient, and it is not possible to obtain a relic density according to observations. On the other hand, for larger , it is possible to get a correct relic density if the difference between the three inert scalar masses is not so large and remains small enough [50] (see Figure 1).

We finally study the direct detection prospects of our DM candidate. We do that through the tree-level spin-independent DM-nucleon interaction cross section [51], which applied to our case,

Here is the mass of the SM-like Higgs boson, is the mass of the DM candidate, is the nucleon mass, taken here to be  GeV as the average of the proton and neutron masses, is the combined coupling defined in (6), and is a form factor that depends on hadronic matrix elements [52, 53].

In Figure 8 we show the DM-nucleon cross section as a function of the DM candidate mass for a correct value of . We consider three values for : a central value (0.326) from a lattice calculation [54], and extreme values (0.260 and 0.629) from the MILC collaboration [55]. Lower bounds for past experiments and prospects of measurements for future experiments are also shown [5659]. Notice that the dispersion of points for high can be understood from analogous dispersion seen in Figure 1, and the same situation occurs with the line-like distribution for light . We also show the coherent neutrino scattering upper limit [60]. This curve represents the threshold below which the detector sensitivity is such that not only can the possible DM scattering effects be observed, but also the indistinguishable scattering effects are associated with neutrinos. Thus, this indicates a region where the neutrino background becomes dominant and little information can be obtained on DM effects. Current direct detection of DM excludes all the low DM mass points, and most of the medium DM mass points. Allowed are a narrow region near 60 GeV and all the high mass region. Note that the absence of points in the range of ~100–550 GeV in this plot is due to the fact that the plotted points are only those that give the right dark matter density (“blue points”). Future experiments will be able to test large parts of these two regions but will not be able to rule them out entirely if there is no signal.

6. Conclusions

In this paper the Inert Higgs Doublet Model is studied, with the inert Higgs boson as a DM candidate, using the latest results for DM relic density, annihilation cross section, and collider searches. As a summary we highlight the following:(i)The branching ratios for the charged Higgs and for the pseudoscalar Higgs are studied and shown in Figure 6. The symmetry strongly reduces the number of different decay channels. Considering the Higgs boson on shell and allowing the gauge boson to be off shell (a different choice would produce different decay channels, but with smaller branching ratios), we find that as opposed to the decays, where there is a crossing point not far from the threshold. For this reason, in collider searches we recommend to look for a signal for a : two jets (consistent with a ) and missing energy (from the DM candidate ).(ii)Three distinct mass regions are found that produced the correct relic density (i.e., is the only source for DM) as can be seen in Figure 7. The low mass region (between 3 and 50 GeV approximately) is already ruled out because it produces a very small value for (Figure 3), because it produces a very high value for (Figure 5) and because of direct DM searches (Figure 8). The intermediate mass region (between 60 and 100 GeV approximately) and the high mass region (heavier than 550 GeV approximately) are allowed.(iii)In Figure 8 we study the DM candidate direct detection. The low mass region is also ruled out by present experiments. In addition, future experiments will probe intermediate and high mass regions. Nevertheless, in absence of signals, it will not be possible to rule out these two regions. Notice the proximity of the region where the coherent neutrino scattering is an irreducible background.

With this one sees that the Inert Higgs Doublet Model gives a still viable DM candidate, which will most likely be tested by direct DM detection experiments.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

The authors are grateful for conversations with Drs. Germán Gómez, Nicolás Viaux, and Edson Carquín. Marco Aurelio Díaz was partly supported by Fondecyt 1141190. Benjamin Koch was partly supported by Fondecyt 1120360. Sebastián Urrutia-Quiroga was partly supported by postgraduate CONICYT grant. The work of all of the authors was also partly financed by CONICYT Grant ANILLO ACT-1102.