Inert Doublet Model ============================================================================ *Jon Butterworth, Louie Corpe, Tania Robens, Gustavs Zilgalvis* The following is a study of the Inert Doublet Model :cite:`PhysRevD.18.2574,Cao:2007rm,Barbieri:2006dq`, a simple two-Higgs doublet model with the scalar potential .. math:: \begin{array}{c} V=-\frac{1}{2}\left[m_{11}^{2}\left(\phi_{S}^{\dagger} \phi_{S}\right)+m_{22}^{2}\left(\phi_{D}^{\dagger} \phi_{D}\right)\right]+\frac{\lambda_{1}}{2}\left(\phi_{S}^{\dagger} \phi_{S}\right)^{2}+\frac{\lambda_{2}}{2}\left(\phi_{D}^{\dagger} \phi_{D}\right)^{2}+\lambda_{3}\left(\phi_{S}^{\dagger} \phi_{S}\right)\left(\phi_{D}^{\dagger} \phi_{D}\right) \\ +\lambda_{4}\left(\phi_{S}^{\dagger} \phi_{D}\right)\left(\phi_{D}^{\dagger} \phi_{S}\right)+\frac{\lambda_{5}}{2}\left[\left(\phi_{S}^{\dagger} \phi_{D}\right)^{2}+\left(\phi_{D}^{\dagger} \phi_{S}\right)^{2}\right] \end{array} that obeys a discrete :math:`Z_2` symmetry and provides a dark matter (DM) candidate :cite:`Kalinowski:2018ylg,Ilnicka:2015jba`. The UFO files come from the Feynrules model library https://feynrules.irmp.ucl.ac.be/wiki/InertDoublet and can also be found here: https://gitlab.com/hepcedar/contur/-/tree/main/data/Models/2HDM/InertDoublet Events are generated with Herwig7 :cite:`Bellm:2015jjp`, including all 2-to-2 processes involving a beyond-the-SM (BSM) particle in the matrix-element or an outgoing leg. In this model, :math:`H^{0}` is considered to be the DM candidate, and only regions where it is the lightest new particle are considered. Taking into account that larger :math:`\lambda` values and larger mass splittings between :math:`A^{0}` and the charged Higgs are excluded by precision measurements of electroweak observables and theory considerations, we use parameter values :math:`\lambda_{2} = 2`, :math:`\lambda_{L} = \frac{1}{2}(\lambda_{3}+\lambda_{4}+\lambda_{5}) = 0.005`, :math:`m_{H^{\pm}}=m_{A^0}+20` GeV, and scan :math:`m_{H^{0}}` and :math:`m_{A^{0}}` over the range [50, 500] GeV in which we have the most sensitivity. A 10x10 grid is used, with 30,000 simulated events for each grid point. The regions where :math:`H^{0}` is not the lightest particle are excluded. .. image:: images/p_p__to__H0_H0.png :scale: 40% .. image:: images/p_p__to__H0_Hpm.png :scale: 40% .. image:: images/cbar.png :scale: 50% .. image:: images/p_p__to__A0_H0.png :scale: 40% .. image:: images/p_p__to__A0_A0.png :scale: 40% .. image:: images/cbar.png :scale: 50% The figures above show some of the dominant cross sections at each point according to the Herwig event generator. The lower half of the plane, coloured in orange, is not considered. White sections of the plots indicate regions where the cross-section was not sampled because it was negligible compared to other processes. The largest cross-section is for low-mass :math:`m_{H^{0}}` pair-production, shown in the upper-left figure. In this case the DM candidate is produced with no other SM objects, and would therefore be a challenging missing-transverse-energy (:math:`E_T^{\rm miss}`) -only signature. In other parts of the plane, :math:`H^{0}` may also be produced with other new particles such as :math:`A^{0}`, or :math:`H^{\pm}`, or pair-production of these other BSM particles may occur. These BSM particles then decay back to a :math:`H^{0}` in association with a SM vector boson, which follow their normal decays to quarks or leptons. For instance, pair-produced :math:`A^{0}` particles may decay to pairs of :math:`H^{0}` and two SM Z bosons. The figures below show the results of propagating these production and decays to final states, which could be observed at the LHC. The final states with the highest cross-sections are shown. Although the :math:`E_T^{\rm miss}`-only signature would have the highest cross-sections, it is experimentally very challenging to analyse this final state at the LHC. It is clear that the signatures more easily observable at the LHC would come from :math:`H^{0}` production associated with jets, leading to :math:`E_T^{\rm miss}` +jets final states. Only few such measurements (as opposed to searches) are available at time of writing. Signatures with one or more leptons with :math:`E_T^{\rm miss}` and jets may also play a role. .. image:: images/nLeptons0_MET.png :scale: 50% .. image:: images/nLeptons0_nJets2_MET.png :scale: 50% .. image:: images/nLeptons1_MET.png :scale: 50% .. image:: images/cbar.png :scale: 50% The heatmaps below show the results of applying the CONTUR method to this model, using for all currently-available LHC data preserved in Rivet (7, 8 and 13 TeV runs in Rivet as of 07/07/2020)]. .. image:: images/dominantPools0.png :scale: 60% .. image:: images/combinedOverlay.png :scale: 50% .. image:: images/combinedMeshcbar.png :scale: 50% As expected from the cross section plots, there is very little sensitivity, with only a small amount of exclusion at about the 60% level showing at low :math:`m_{H^0}`. As can be seen in the left-hand plot, the :math:`E_T^{\rm miss}` +jets analyses are the most sensitive across the plane. However, other signatures play a role in particular at low :math:`m_{H^{0}}`, where the overall sensitivity is greatest. These additional signatures include the EEJET and WW pools, which collect lepton+jet and lepton+ :math:`E_T^{\rm miss}` +jets signatures, and the four-lepton invariant mass measurement, which is sensitive in the rare cases when DM candidates are produced with four leptons, such as :math:`{A^{0}}{A^{0}} \rightarrow {H^{0}}{H^{0}} ZZ \rightarrow {H^{0}}{H^{0}} \ell \ell \ell \ell`. The :math:`E_T^{\rm miss}` +jets measurements and searches are available in Contur only for 3.2/fb, so the 139/fb full-run-2 updates are expected to boost sensitivity. The CLs exclusion is low in all regions of the heatmap, showing that the Inert Doublet Model is very challenging to exclude, as it appears to evade many of the measurements made at the LHC thus far. It is therefore a good candidate for dedicated searches with higher integrated luminosity.