Long live the Higgs portal!

In Higgs portal models of fermion dark matter, scalar couplings are unavoidably suppressed by strong bounds from direct detection experiments. As a consequence, thermal dark matter relics must coexist with mediators in a compressed spectrum of dark particles. Small couplings and small mass splittings lead to slow mediator decays, leaving signatures with displaced vertices or disappearing tracks at colliders. We perform a comprehensive analysis of long-lived mediators at the LHC in the context of a minimal dark matter model with a naturally small Higgs portal, also known as the wino-bino scenario in supersymmetry. Existing searches for disappearing charged tracks and displaced hard leptons already exclude tiny portal couplings that cannot be probed by current direct and indirect detection experiments. For larger portal couplings, we predict new signatures with displaced soft leptons, which are accessible with run-II data. Searches for displaced particles are sensitive to weakly coupling mediators with masses up to the TeV scale, well beyond the reach of prompt signals.


Introduction
The hypothesis of thermal Higgs portal dark matter had almost been declared a dead end. If dark matter interacts with the standard model through the Higgs boson, the portal strength is indeed strongly constrained by null results from direct detection experiments [1][2][3][4][5][6][7]. A dark matter relic from thermal freeze-out would thus be overabundant today, unless its annihilation rate in the early universe was enhanced, for instance, by annihilation through a resonance or co-annihilation with mediators, or altered due to a non-standard cosmological history. In Higgs-portal scenarios with fermion dark matter, mediators are a necessary prediction of a UV-complete theory. In this work, we assume that dark matter is a fermion with no charges under the standard model gauge group, interacting only through electroweak mixing. Fermionic mediators with renormalizable Higgs couplings induce generally strong dark matter-nucleon scattering through electroweak mixing [8][9][10][11][12]. Mediator triplets or higher multiplets of the weak gauge group lead to non-renormalizable Higgs JHEP02(2019)140 couplings, which are naturally suppressed by a high cutoff scale [13]. The small portal coupling implies slow decays of the heavy dark states. This can result in a decay length of several centimeters or more, before the mediator decays into dark matter and leptons or hadrons. At the LHC, long-lived mediators thus leave characteristic signatures of displaced vertices or disappearing charged tracks in the detector.
Higgs-portal scenarios with small dark matter couplings are thus a perfectly viable option for thermal dark matter that can be tested at colliders. In this work, we explore the prospects to discover thermal Higgs-portal dark matter through signatures with long-lived mediators at the LHC. We focus on a non-renormalizable Higgs portal with a Majorana fermion singlet and a fermion triplet in the adjoint representation of the weak gauge group. The portal coupling is naturally small and thus evades direct detection bounds. Our model can be UV-completed for instance by a fermion doublet or a scalar triplet. The former construction is known as the wino-bino scenario in the context of supersymmetry [14], the latter is similar to type-II seesaw models for neutrino masses [15,16]. Due to its rich dark matter phenomenology and promising collider signatures, the small coupling regime of this and similar models has recently received increased interest [17][18][19][20][21][22]. We perform a comprehensive analysis of two scenarios: a scalar and a pseudo-scalar Higgs portal. While the scalar scenario has been the focus of most previous collider studies, the phenomenology of long-lived mediators in the pseudo-scalar scenario is still largely unexplored. We point out the characteristic differences between both scenarios, regarding their dark matter and collider phenomenology. We propose new signatures with displaced particles and show that they can probe mediator masses well beyond the reach of prompt signatures at the LHC.
This article is organized as follows. In section 2, we introduce our model and discuss the relations between the scalar and pseudo-scalar Higgs portals. Section 3 is a brief interlude on the interpretation of our model in the minimal supersymmetric standard model (MSSM). In section 4, we derive constraints on the portal strength from direct detection experiments. Section 5 is devoted to a detailed discussion of the relic dark matter abundance in different regimes of a small portal coupling. In section 6, we eventually explore the collider phenomenology of the two scenarios under the assumption of a thermal dark matter candidate. We conclude in section 7.

The singlet-triplet Higgs portal
We extend the standard model by two self-adjoint fermion fields with vector-like weak interactions, Here χ S is a standard-model singlet Majorana fermion and χ T transforms under the weak gauge group as a triplet with zero hypercharge. We assume a discrete Z 2 symmetry, under which χ S and χ T are odd and all standard-model particles are even. The lightest fermion state is stable and a dark matter candidate. In this scenario, there are no renormalizable

JHEP02(2019)140
Higgs couplings to dark fermions. At energies below a cutoff scale Λ, the scalar sector is described by the effective Lagrangian 1 where D µ = ∂ µ − ig τ a W a µ is the covariant derivative, τ a are the generators of the weak SU(2) gauge group, and H = (h + , is the Higgs doublet of the standard model. Gauge-invariant contractions of the fields are implicitly assumed. Higgs couplings to dark fermions are of mass dimension five and thus naturally small at energies well below Λ. Since we will focus on the parameter region with a high cutoff scale Λ M W , the impact of operators with higher mass dimensions can be neglected. We assume all parameters to be real in order to preserve CP invariance. The dark matter and collider phenomenology of the Higgs-portal couplings κ S and κ T has been extensively studied [23][24][25][26]. In general, a complete theory that induces a singlet-triplet Higgs portal can also generate singlet-singlet and/or triplet-triplet Higgs couplings. However, since κ S and κ T , κ T will play essentially no role in our phenomenological analysis, we neglect them by setting κ S = κ T = κ T = 0.
After electroweak symmetry breaking, χ 0 S and χ 0 T mix through the Higgs portal. The scalar Lagrangian for the neutral fermions reads We can choose m T positive without losing generality. The sign of µ is not observable in tree-level processes, so we assume it to be positive as well. Furthermore, we always request m T > −m S , so that the singlet fermion is lighter than the triplet. Due to electroweak mixing, the gauge eigenstates are not mass eigenstates of our model. We define a mass matrix M by By diagonalizing M , the physical eigenstates χ and χ h are readily obtained as The corresponding mass eigenvalues for the neutral states χ , χ h and the charged states

JHEP02(2019)140
For m S > µ 2 /m T , the mass m of the lightest state is positive. The physical Lagrangian then reads (neglecting interactions with two Higgs bosons) where θ w denotes the weak mixing angle. Singlet-triplet mixing induces scalar Higgs couplings and charged vector currents of the lightest state χ . Neutral currents of χ and χ h are absent due to the Majorana nature of the neutral fermions. In the limit of small mixing, the lightest state is mostly a weak singlet with suppressed couplings to the standard model. A pseudo-scalar singlet-triplet Higgs portal can be obtained through a chiral rotation of the singlet fermion, (2.8) By applying this transformation to the Lagrangian in eq. (2.3), we obtain The chiral rotation turns the scalar portal into a pseudo-scalar portal and flips the sign of m S . In the mass eigenbasis, we obtain The physical mass terms and interactions are given by Now the Higgs interactions with one heavy and one light neutral fermion are pseudo-scalar, while heavy-heavy and light-light couplings remain scalar. The gauge couplings of χ have an axial-vector structure. The chiral rotation also flips the sign of the lightest mass eigenvalue in the spectrum, see eq. (2.10). The parameter space m S < µ 2 /m T with negative mass m < 0 in the scalar Lagrangian L S thus corresponds to positive mass m > 0 in the pseudo-scalar Lagrangian L P . We identify two physical scenarios, scalar scenario: L S with m S > µ 2 /m T ↔ m > 0, (2.12) pseudo-scalar scenario: L P with m S < µ 2 /m T ↔ m > 0. It is instructive to study the features of our model in the limit of small singlet-triplet mixing. In this limit, the couplings of the lightest state χ are approximated by

JHEP02(2019)140
(2.13) For fixed values of m T and µ, the mixing θ in the pseudo-scalar scenario is smaller than in the scalar scenario, see eq. (2.12). Gauge and diagonal Higgs couplings of the lightest state are thus weaker in the pseudo-scalar scenario. The masses of the neutral fermions are approximated as (2.14) The lightest state χ is thus mostly a gauge singlet, while the heavier states χ + and χ h approximately correspond to the charged and neutral components of a weak triplet. In both scenarios, the mass splitting between the heavier states is given by The first contribution is due to singlet-triplet mixing and the second contribution due to electroweak loop corrections [27]. The mass difference between the charged state and the lightest state is different in the two scenarios. Neglecting electroweak corrections and assuming small fermion mixing, one finds In figure 1, we show two typical scenarios of dark fermion spectra. For very small singlettriplet mixing, electroweak corrections dominate the mass splitting ∆m hc and the charged fermion is the heaviest state (right panel). For larger mixing, χ + and χ h become degenerate in mass and eventually flip their positions (left panel). The spectra in the scalar and pseudoscalar scenarios look very similar, unless the mixing is large. The mass hierarchy plays a crucial role for the collider phenomenology of the dark fermions (see section 6).

JHEP02(2019)140
3 In case of supersymmetry Our model can be interpreted as the so-called wino-bino scenario in the MSSM with conserved R parity, where the superpartners of the gauge bosons reside around the electroweak scale and the higgsinos, as well as all scalars but the standard-model-like Higgs boson, are much heavier. This scenario is commonly referred to as split supersymmetry [28] with heavy higgsinos. Here we briefly discuss this supersymmetric realization of our Higgs-portal scenario. We adopt the notation for the wino-bino scenario from ref. [21] and relate it to ours. The phenomenology at energies around the weak scale is described by an effective Lagrangian, where W A and B are the wino and bino fields, C 1 and C 2 are Wilson coefficients of the dimension-five operators. 2 The standard-model-like Higgs field H is a linear combination of the two Higgs doublets H u and H d [28]. This Lagrangian is to be compared with the Higgs-portal Lagrangian in eq. (2.2). The mass parameters m S and m T in our model correspond with the gaugino masses M 1 and M 2 . The Higgs-portal coupling µ/v is related to the higgsino mass parameter µ through where tan β = v u /v d is the ratio of vacuum expectation values of the two Higgs fields H u and H d , and M W is the mass of the W boson. The cutoff Λ can thus be interpreted as the higgsino mass µ, provided that κ ST and sin(2β) are not suppressed. For tan β = 1, a portal coupling of µ/v = 0.01 corresponds to a higgsino mass of µ = 1.4 TeV. As we will see in section 6, collider searches with long-lived mediators probe parameter regions with µ/v 0.01 and thus higgsino masses well beyond resonant production at the LHC. The scalar and pseudo-scalar scenarios correspond to different choices of the phases in the diagonalization of the bino-wino mass matrix (see eq. (17) in [21]), pseudo-scalar scenario: α = β = 0, φ = π/2.
Using these relations, one can directly interpret the results of our work in the wino-bino scenario. The pseudo-scalar scenario is a special case of the wino-bino scenario the complex MSSM with flipped bino mass M 1 → −M 1 . Signatures of long-lived winos scenario have been discussed in the scalar scenario [20,29,30]. To the best of our knowledge, a comprehensive analysis of LHC signatures with displaced particles has not been performed yet. The phenomenology of mediator decays in the pseudo-scalar scenario has been much less explored [31]. Our work can serve as a framework to systematically search for supersymmetric gauginos in the limit of higgsino decoupling. Notice that in the MSSM the dark JHEP02(2019)140 matter phenomenology and the resulting collider signatures can be altered by the presence of other light superpartners or additional scalars. The interpretation of our results in a specific scenario should thus be done with care.

Dark matter scattering off atomic nuclei
From now on, we will interpret the lightest neutral fermion χ as a dark matter candidate. An important bound on our model is derived from searches for dark matter scattering off atomic nuclei in direct detection experiments. Since χ does not couple to the Z boson, spin-independent scattering is mediated by Higgs boson exchange via a scalar current. 3 Thanks to the small momentum transfer, the interaction of dark matter with the quarks inside the nucleon can be described by an effective Lagrangian, where G F is the Fermi constant and M h is the Higgs mass. Notice that the Higgs coupling to two dark matter states is not affected by a chiral rotation and is scalar, regardless of whether the Higgs portal has a scalar or pseudo-scalar structure. In terms of the effective interaction, the coupling of dark matter to protons can be expressed as and analogously for neutrons with p → n. The contribution of a particular quark flavor to the proton mass is m p f (p) T q = p|m qq q|p , determined experimentally. The cross section for spin-independent dark matter scattering off a nucleus with Z protons and A nucleons in total is given by is the reduced dark matter-nucleon mass and σ n is the cross section for spin-independent dark matter-nucleon scattering. The currently strongest upper bound on nucleon scattering has been obtained by the Xenon1T collaboration [32], In our model, this translates to a strong bound on the Higgs coupling to dark matter, which we show in figure 2 as a function of m . For instance, for a dark matter mass of 100 GeV we obtain For sizeable mass splittings m T − m S µ, fermion mixing is small. In this regime, the Xenon1T results can be interpreted as an upper bound on the mass splitting between the mediators, The mediator states are thus nearly degenerate in mass. For smaller mass splittings m T − m S ≈ 15-30 GeV, as favored by co-annihilation, µ/v must be suppressed to evade direct detection. 4 In this regime, fermion mixing can still be close to maximal. As µ/v is lowered below the Xenon1T bounds, fermion mixing decreases and dark matter decouples from the standard model. Future direct detection experiments are expected to probe even smaller dark matter-nucleon scattering cross sections. If they were able to test rates comparable to coherent neutrino scattering (see figure 2), this would probe Higgs couplings at the permil level.

Dark matter annihilation and relic abundance
The interpretation of the lightest state as dark matter strongly depends on the Higgsportal strength µ/v. We consider couplings that are large enough for dark matter to be in thermal equilibrium with the primordial plasma before freeze-out. The relic abundance is then determined by the freeze-out of processes that change the dark matter number density.
In table 1, we show the annihilation and scattering processes relevant around the freeze-out temperature and their dependence on the model parameters µ/v and θ. Throughout our analysis, we consider dark matter masses below 1 TeV, corresponding to the region that JHEP02(2019)140 process scaling process scaling can be probed at the LHC. In this mass range, the effect of Sommerfeld enhancement on pair annihilation is mild and will be neglected in our analysis [35]. Scattering off standardmodel fermions keeps the dark fermions in thermal equilibrium. Co-scattering, mediator decays and inverse decays ensure chemical equilibrium among the dark fermions.

Pair annihilation
Dark matter pairs can annihilate through the Higgs boson or through weak charged currents. Since direct detection results set stringent bounds on the Higgs coupling and also on the mixing of dark fermions below the TeV scale, pair annihilation is suppressed. In figure 3, we show the parameter regions of our model that are excluded by Xenon1T in the ∆m c − m plane for three values of µ/v. The grey areas are excluded by searches for pair-produced charged fermions at LEP [36]. For µ/v = 0.2, Xenon1T excludes large parts of the parameter space, while for µ/v = 0.02 Xenon1T sets weaker bounds than LEP. Notice that in the pseudo-scalar scenario the bounds from direct detection are much weaker, since the mixing θ for fixed µ/v, m , and ∆m c is smaller than in the scalar scenario. The colored curves show the observed relic abundance Ω χ h 2 = 0.1199 [37] for fixed µ/v. To obtain these predictions, we have implemented our model in micrOMEGAs [34]. The vertical lines around m = 63 GeV correspond to the freeze-out of pair annihilations through the Higgs resonance. As can be seen in the figure, resonant annihilation is compatible with the observed relic abundance even for small µ/v. For larger dark matter masses, the relic abundance is determined by pair annihilation through gauge interactions, provided that fermion mixing is sizeable. In the scalar scenario, this is the case for µ/v = 0.2, while in the pseudo-scalar scenario pair annihilation is too small to provide the correct relic abundance.

Co-annihilation
For smaller µ/v, pair annihilation is suppressed by θ 4 and becomes irrelevant around the freeze-out temperature. The relic abundance is now set by co-annihilation processes like χ χ + → ff [38], which scale as θ 2 . Since the thermally averaged co-annihilation cross section is Boltzmann-suppressed, a moderate mass difference of ∆m c /m ≈ 10 % is required JHEP02(2019)140 to prevent overabundance. In the scalar model, for large dark matter masses co-annihilation becomes relevant at larger µ/v, since the pair annihilation rate is smaller than for lighter dark matter. Similarly, in the pseudo-scalar model co-annihilation starts at larger µ/v than in the scalar model, because the mixing suppression of pair annihilation is stronger.

Mediator annihilation
For even smaller µ/v, the co-annihilation rates become inefficient and dark matter decouples earlier than the heavier dark fermions. The relic density in this case is determined by mediator annihilation χ + χ − → ff (and similar processes involving χ h ), which can still change the dark matter number density through mediator decays χ + → χ ff . Mediator annihilation is not suppressed by fermion mixing (see table 1), but the thermally averaged rate is Boltzmann-suppressed by the mediator mass. Therefore χ + and χ h should not be much heavier than χ . Mediator annihilation thus predicts a compressed spectrum, similar to co-annihilation but with smaller dark matter couplings.
In summary, direct detection sets an upper bound of µ/v 0.2 (0.6) on viable thermal relics below 1 TeV in the scalar (pseudo-scalar) scenario. In the pseudo-scalar scenario, all processes decouple for comparably larger Higgs-portal couplings, due to the smaller fermion mixing in this scenario. In particular, for all three displayed values of µ/v the relic abundance away from the Higgs resonance is set by mediator annihilation.

Co-scattering
In the region of mediator annihilation, dark matter is still in chemical equilibrium with the heavier dark fermions through co-scattering and decays. While decays, co-scattering and co-annihilation all scale as θ 2 or (µ/v) 2 , the latter is relatively suppressed by the number density of the non-relativistic mediators and thus decouples earlier. In this regime, determining the relic abundance requires solving the coupled system of Boltzmann equations for the number density evolution of all dark fermions, taking into account co-scattering and decay processes. This approach, however, is not included in automated programs such as micrOMEGAs, DarkSUSY and MadDM [34,39,40]. While a detailed numerical analysis of the JHEP02(2019)140 non-equilibrium processes is beyond the scope of our work, we offer a qualitative discussion of the dark matter phenomenology in this regime.
Mediator annihilation determines the relic abundance, as long as mediator decays are prompt around the freeze-out temperature. For a fixed value of µ/v, decays are still rapid in our model when co-annihilation processes have already decoupled. Once the mediator decays drop below the Hubble rate, the dark matter number density can only decrease through co-scattering processes χ SM ↔ χ + SM and χ SM ↔ χ h SM, followed by mediator annihilation. This happens only for very small µ/v, where mediator decays decouple before co-scattering processes. Eventually, the latter decouple as well and dark matter departs from chemical equilibrium, while the mediators remain in equilibrium. The relic abundance is now driven by the freeze-out of co-scattering processes and thus very sensitive to the strength of the Higgs portal. Similar scenarios, dubbed conversion-driven freeze-out, have been identified in refs. [41,42]. In non-standard cosmological scenarios with early matter domination, the observed relic abundance can also be obtained out of thermal equilibrium [43].
In figure 4, we illustrate the different phases of dark matter freeze-out for small Higgsportal couplings in the scalar scenario (blue) and the pseudo-scalar scenario (green). For concreteness, we choose a benchmark point with m = 500 GeV and a mass splitting of ∆m c = 30 GeV, as it is typical for efficient co-annihilation and mediator annihilation. For other dark matter masses, the main features of the phase diagram are very similar. In the pseudo-scalar scenario, each phase is reached at a larger coupling µ/v. This is due to the fact that the freeze-out of the relevant processes is very sensitive to the fermion mixing, which is smaller than in the scalar scenario [44]. The different Lorentz structure of Higgs and gauge couplings in both scenarios (cf. eq. (2.13)) has only a subleading effect on the annihilation rates. In particular, all processes that dominate the dark matter abundance can proceed in an s-wave in either scenario (cf. ref. [45]).

JHEP02(2019)140
For large values of µ/v, the relic abundance is set mostly by pair annihilation and co-annihilation. Part of this region is excluded by Xenon1T (shaded in grey). For µ/v 0.02 (0.2) in the scalar (pseudo-scalar) scenario, the relic abundance is determined by mediator annihilation, mostly by χ h χ + → SM SM. Part of this region might be probed by future direct detection experiments, but the region below the neutrino floor (indicated by a dashed line) is not accessible with current methods. At µ/v 3 × 10 −4 (10 −3 ), the thermally averaged co-scattering rate becomes smaller than the (thermally averaged) mediator annihilation rate (see the dashed line). Mediator decays χ + → χ ff remain fast, so that co-scattering does not affect the relic abundance yet. At µ/v 2 × 10 −7 (5 × 10 −6 ), charged mediator decays drop below the Hubble rate and the relic abundance is determined by co-scattering. 5 Thermal equilibrium is preserved for couplings well below the range considered in this work. The regimes of mediator annihilation and co-scattering thus provide us with a thermal dark matter candidate that cannot be tested by direct detection experiments. Therefore colliders play an important role in probing Higgs portal dark matter with tiny portal couplings.

Long-lived mediators at the LHC
The hypothesis of a Higgs-portal dark matter relic with small couplings is directly testable at colliders. In this section, we investigate the LHC phenomenology of our model, constraining the parameter space with existing searches and predicting new observables that test regions that have not been explored yet. At the LHC, the mediators are pair-produced through Drell-Yan-like processes and subsequently decay into the lightest dark fermion. Two examples of such processes are shown in figure 5. The production rate is set by the invariant mass of the mediator pair and the weak gauge coupling. 6 LEP has set a lower bound on the mass of the charged fermion, m c 100 GeV. Direct detection results imply that the Higgs-portal coupling must be small, µ/v 0.2 for mediators below the TeV scale. Furthermore, direct detection sets an upper bound on the mass splitting between mediators, ∆m hc . Viable scenarios of dark matter freeze-out favor a moderate mass splitting between dark matter and the mediators, ∆m c . These requirements determine the parameter region of interest as m c 100 GeV, ∆m hc few GeV, In this parameter range, fermion mixing is small. At the LHC, we search for a compressed spectrum of dark fermions, featuring potentially long-lived mediators with soft decay products and missing energy. In what follows, we first determine the lifetimes of the mediators and then discuss the resulting signatures and how to test them. 5 The heavy neutral fermions become long-lived already at larger values of µ/v ∼ 10 −4 (0.005), due to the smaller decay rate for χ h (see section 6). However, since processes involving χ h contribute less than about 10 % to mediator annihilation in this region, we expect that the relic abundance can still be reliably obtained without including co-scattering in the Boltzmann equations. 6 For large fermion mixing, the production rate of pp → χ + χ h is suppressed by a factor (g cos θ) 2 .

Mediator decays
The lifetime of the mediators depends sensitively on the Higgs-portal coupling, µ/v, and on the mass difference between the dark fermions in the initial and final states, ∆m. Neglecting a potential phase-space suppression, the mediator decay width scales as where x and y depend on the decay process. For decays into the lightest state, χ + → χ and χ h → χ , the mass splitting is sizeable, ∆m c ≈ ∆m h ≈ 15 − 30 GeV. The mediator is thus long-lived on collider scales only if its decay is suppressed by a small portal coupling µ/v. In the decay χ + → χ h , in turn, ∆m hc can be arbitrarily small, as we discussed in section 2. In this case, the mediator decay is suppressed by the mass splitting. Let us first focus on the charged state χ + . Depending on the Higgs-portal strength, χ + can decay either via the two-body decay χ + → χ h π + or the three-body decays χ + → χ + ν and χ + → χ + hadrons. The two-body decay is kinematically allowed for mass splittings larger than the pion mass. In our model this is fulfilled for very small µ/v, where the mixing-induced mass difference is negligible and ∆m hc ≈ −160 MeV, see eq. (2.15). For m π < |∆m hc | m c , the decay width is given by where f π 130 MeV is the pion decay constant and V ud is a CKM matrix element. The decay rate is strongly suppressed by the small mass difference ∆m hc , as well as by the limited kinematic phase space. Since fermion mixing is tiny in this regime, cos θ ≈ 1 and the decay rates in the scalar and pseudo-scalar scenarios are the same. The nominal decay length of the charged fermion for |∆m hc | ≈ 160 GeV is given by where τ χ is the proper lifetime of χ + . Charged particles with a decay length in the centimeter range leave tracks in the inner layers of the ATLAS and CMS detectors, which we will discuss in more detail below. As µ/v increases, the mass splitting ∆m hc drops below the pion mass and the two-body decay is forbidden. In this regime the three-body decays χ + → χ + ν and χ + → χ + hadrons dominate the decay rate. In the scalar and pseudoscalar scenarios, three-body decays proceed through vector currents (2.7) and axial-vector currents (2.11), respectively. For ∆m c M W , m c , the partial width of the leptonic decay χ + → χ + ν is given by [46] In our analysis, we estimate hadronic three-body decays using the branching ratios from ref. [47].

JHEP02(2019)140
in the scalar and pseudo-scalar scenarios. The mass splitting is very similar in both scenarios, ∆m c ≈ ∆m c . Contrary to pion decay, three-body decays are not suppressed by the mass splitting, but by the small fermion mixing θ. Notice that in the pseudo-scalar scenario the three-body decays are smaller than in the scalar scenario, leading to a longer mediator lifetime. This is due to the smaller mixing angle, which reduces the weak coupling to dark matter, see eq. (2.13). In the region where three-body decays dominate, the nominal decay length reaches up to cτ χ + ≈ 1.5 cm (4 cm) in the scalar (pseudo-scalar) scenario. For larger portal couplings, the decay length decreases. Leptonic three-body decays lead to displaced soft leptons, which can in principle be observable at the LHC for displacements larger than 200 µm. Hadronic decays lead to soft jets, which are difficult to detect.
The lifetime of the heavy neutral state χ h depends on the mass hierarchy of the dark fermions. For sizeable µ/v, corresponding to a normal mass hierarchy ( figure 1, left), χ h decays dominantly via χ h → χ + π − , provided that the channel is kinematically allowed. In this case, the decay rate of χ h is the same as for χ + in the limit µ/v → 0, given by eq. (6.3). For smaller µ/v, the mass hierarchy is inverted ( figure 1, right), and χ h decays dominantly through χ h → χ bb via an off-shell Higgs boson. For ∆m h m h , the partial decay width in the respective scenarios is given by The decay rate is not suppressed by mixing, but by the small Higgs-portal and bottom Yukawa couplings. 8 In the pseudo-scalar scenario, the decay proceeds through a pseudoscalar current. Parity conservation requires that χ and the bb pair are emitted in a relative p-wave near the kinematic endpoint, resulting in an additional suppression of (∆m h /m h ) 2 . Therefore the heavy neutral fermion lives longer in the pseudo-scalar scenario than in the scalar scenario. Compared with χ + , the strong suppression of χ h decays results in a longer lifetime of the neutral mediator. For µ/v → 0, χ h can be arbitrarily long-lived. The typical signature of a long-lived neutral fermion is a pair of b-jets with displaced vertices, which is expected to be detectable by ATLAS and CMS in the range of 1 cm < cτ χ h < 1 m [21].

LHC signatures with displaced particles
We are now prepared to investigate the predicted LHC signatures in detail. In figure 6, we present the parameter region of our model that can be tested with prompt and long-lived mediators in the m c versus µ/v plane. For each parameter point, the mass difference ∆m c is determined by requiring the observed dark matter abundance of Ω χ h 2 = 0.1199 ± 0.0022. In most of the parameter space, co-annihilation or mediator annihilation set the relic abundance (cf. figure 4). At the lower end of the plots, the relic abundance starts to be determined by co-scattering, as indicated by a dashed line.

JHEP02(2019)140
Null results of Xenon1T exclude the upper left corner in the scalar scenario. 9 In the pseudo-scalar scenario, direct detection bounds are weaker and do not appear in the figure. In either case, current direct detection experiments do not probe parameter regions with long-lived mediators yet. If future experiments became sensitive to scattering rates comparable to coherent neutrino scattering, they would test the region with µ/v 0.01. Indirect detection cannot probe small portal couplings, as dark matter pair annihilation is strongly suppressed. Colliders are thus the only terrestrial instruments to date that can test the hypothesis of dark matter from a tiny Higgs portal.
The various signatures are classified as follows. Green regions are already excluded by existing searches; orange regions correspond to predictions for the LHC with full run-II data (displaced b-jet pairs) or for the HL-LHC (disappearing charged tracks); red regions have not been explored yet, but can be probed with new signatures we predict in this work.
In what follows, we discuss the signatures one by one, starting from tiny portal couplings and moving upwards in the parameter space.

Disappearing charged tracks
When the portal coupling is tiny, the mass difference between χ + and χ h is induced radiatively by electroweak corrections and decreases with increasing µ/v. In this parameter region, the spectrum is inverted and χ + is the heaviest particle in the dark sector (see figure 1, right). If the splitting is saturated by electroweak corrections, the decay χ + → χ h π + is kinematically allowed and dominant. The charged mediator decays with a nominal length up to 7 cm (see eq. (6.4)). It leaves a track in the innermost layers of the detector and decays before reaching the outer tracking layers [30,47,48]. The ATLAS collaboration has performed a dedicated search for supersymmetric winos with similar decay length [49], which is directly applicable to our case. Notice that the analysis assumes that the byproduct of chargino production, the neutral wino, is stable at collider scales. In our model, χ h decays via χ h → χ bb, which is strongly suppressed by the tiny portal coupling in this regime. This assumption is thus fulfilled by our heavy neutral fermion.
Interpreting the ATLAS results in our scenarios, we find that it excludes mediator masses up to m c ≈ 460 (480) GeV and portal couplings of µ/v 10 −6 (10 −5 ) in the scalar (pseudo-scalar) scenario. Notice that parts of the parameter space with disappearing charged tracks correspond to co-scattering. In this region, the mediator's decay length can deviate from that predicted from co-annihilation. Since co-scattering saturates the relic abundance at smaller couplings µ/v, the mass splitting tends to be larger, but without exceeding ∆m hc ≈ 160 MeV. The shown results are thus not expected to change much in the co-scattering phase. In particular, we expect the area below the plot to be excluded by disappearing charged track searches as well. At the upper edge of the excluded region, the rates for two-and three-body decays are equal, Γ(χ + → χ h π + ) = Γ(χ + → χ + ν)+Γ(χ + → χ +hadrons). In the pseudo-scalar scenario, this condition is met at larger µ/v, due to the smaller three-body decay rate (see eq. (6.5)). 9 The bound is conservative in the sense that the region is excluded for any dark matter mass below the TeV scale. Above this boundary, three-body decays dominate. In the semi-transparent regions, the branching ratio for χ + → χ h π + is less than 90 %, which may weaken the bounds derived assuming 100 % decay into pions.

JHEP02(2019)140
In ref. [47], a dedicated study of the HL-LHC prospects for disappearing charged tracks has been performed for supersymmetric higgsinos. Rescaling the production rate in these predictions for our wino-bino scenario, we estimate that with a luminosity of 3 ab −1 the HL-LHC can extend the search to mediator masses up to 1 TeV.

Displaced soft leptons
In the region where three-body decays dominate the lifetime of χ + , a soft charged lepton or soft jets are a typical signature of χ + → χ ff decays. Here we focus on signatures with displaced soft leptons. The dominant process with a final state of two soft leptons and missing energy is shown in the left panel of figure 5. Due to the small mass splitting 15 ∆m c 30 GeV, the transverse momenta of the leptons typically range around p T, ≈ 5 − 40 GeV. For a sufficiently small portal coupling, the charged mediator can be long-lived and the soft leptons originate from a displaced vertex. The longest decay length of χ + is obtained for small portal couplings, cτ χ + = 1.5 cm (4 cm) in the scalar (pseudoscalar) scenario. The smallest decay length that can be observed, cτ χ + ≈ 200 µm, is limited by the vertex resolution of the detector [50]. The corresponding parameter region is shown in figure 6.
Searches for displaced leptons have been performed by CMS at 8 TeV with a lepton momentum cut of p T, > 25 GeV [51] and at 13 TeV for p T, > 40 GeV [50]. While the momentum cut of the 13-TeV analysis is too strong for our model, a fraction of events with displaced soft leptons falls into the signal region of the 8-TeV analysis. In ref. [17], this analysis has been recasted for quintuplet dark matter, which leads to the same final state of displaced soft lepton pairs through a decay chain of doubly-charged fermions. We rescale the quintuplet event rates by a factor of 1/4 and derive the bounds for our model, requiring that the decay length of χ + is equal to the decay length of the doubly-charged quintuplet fermion. The resulting 95 % CL exclusion bound is shown as a green area in JHEP02(2019)140 figure 6. Displaced lepton searches exclude mediator masses up to m c ≈ 200 GeV for portal couplings µ/v 3 × 10 −5 (3 × 10 −4 ) in the scalar (pseudo-scalar) model. The upper boundary on µ/v is determined by the minimal decay length that can be probed with the CMS search. The sensitivity of the displaced lepton search extends up to decay lengths of cτ χ + 7 cm (and thus to smaller µ/v), where disappearing charged tracks can be observed. In the light green region, the mass splitting ∆m c is up to 25 % smaller than in the quintuplet model, and only very few soft leptons still pass the p T, cut. The sensitivity to our model is thus strongly limited in this region.
We therefore suggest to extend searches for displaced lepton pairs to lower transverse momenta. As the lepton momenta spread over a certain range, it could be experimentally beneficial to first lower the threshold for one displaced lepton, i.e., to search for signals with one soft and one harder displaced lepton. Soft displaced leptons can potentially probe much larger mediator masses, as indicated by the red area. The highest accessible masses are experimentally limited by the signal-background discrimination efficiency for small event rates.

Displaced b-jet pairs
In the region where χ + decays are prompt, parameter regions with heavy mediators can still be probed by displaced signatures if χ h is long-lived. This is indeed the case, because the decay χ h → χ bb through an off-shell Higgs boson is small (see eq. (6.6)). Produced via pp → χ + χ h , the slow χ h decay leaves a signature of a pair of displaced b-jets. Due to the small mass splitting ∆m h , the b-jets are rather soft. In ref. [21], such a signature has been analyzed in the context of a supersymmetric wino-bino scenario. Projections for the LHC running at 14 TeV are made, assuming an integrated luminosity of 300 fb −1 . The detection criteria for displaced vertices are adapted from an ATLAS study based on 8 TeV data. Accordingly, a good detection efficiency can be achieved for a decay length in the range 1 cm cτ χ h 1 m, reaching its maximum around cτ χ h ≈ 10 cm. Reinterpreting these predictions for our scenarios, we derive the parameter region that can be probed with displaced b-jet pairs, as indicated by the orange area in figure 6. The lower edge of this area corresponds to the largest detectable decay length, cτ χ h = 1 m. The upper edge is set either by the sensitivity limit cτ χ h = 1 cm or, in the pseudo-scalar scenario for small m c , by m h = m c + m π , where the spectrum is inverted and χ h decays dominantly via χ h → χ + π − . The hatched regions indicate the (large) uncertainties on the displaced vertex reconstruction [21]. Overall, it is expected that mediator masses up to m h ≈ 800 GeV can be probed for appropriate couplings µ/v. Neutral mediators with longer lifetimes, corresponding to a smaller portal coupling µ/v, escape the LHC detector before decaying. Proposed surface detectors might be able to probe this scenario [19] and thus cover parts of the area below the yellow band in figure 6.

Displaced and prompt soft leptons
For m h > m c , the mass splitting due to fermion mixing dominates over electroweak corrections. In this region, the portal coupling is too large to cause displaced χ + decays. However, χ h can be long-lived for m h m c + m π , due to a strong phase-space suppression JHEP02(2019)140 by the small mass difference ∆m hc . From pair-produced mediators via pp → χ h χ + , one expects one prompt decay, χ + → χ + ν, and one slow decay, χ h → π − /ff (χ + → χ + ν), which produces a displaced soft lepton ( figure 5, right). We thus predict a signature with large missing energy and two soft leptons, one displaced and one prompt. The leptons can be either of opposite or of same electric charge, appearing at equal rates. The pion from the decay χ h → π − χ + is very soft and not detected.
In figure 6, we show the parameter space for displaced and prompt soft lepton signals in red. The upper edge of the area is determined by requiring a minimal nominal decay length of cτ χ h 200 µm. Above this line, the mass splitting ∆m hc is larger and decays become rapid. The lower edge is set by the kinematic threshold for a two-body decay, m h = m c + m π . Below this threshold, the heavy neutral fermion decays mostly via χ h → χ bb, and no soft leptons are produced.

Prompt soft leptons
Mediators with a decay length cτ χ + < 200 µm leave signatures with prompt soft leptons. Both χ + and χ h decays can produce prompt soft leptons via χ + → χ + ν or χ h → π − (χ + → χ + ν), respectively. Signatures with two prompt leptons and missing energy are thus expected from pp → χ + χ − and pp → χ + χ h production. Which process dominates depends on the χ h decay branching ratio into leptons. If χ h is long-lived, pairs of prompt soft leptons can only be produced via pp → χ + χ − with subsequent prompt decays. Searches for pairs of prompt soft leptons have been performed by ATLAS [52] and CMS [53] in the context of pure wino production via pp → χ + χ h . The analyses assume a decay via χ + χ h → χ W * χ Z * → χ + ν χ + − to 100 %. However, in our model χ h → χ Z * → χ + − is loop-and mixing-suppressed and thus small compared to χ h → χ h * → χ bb. In ref. [17], the CMS analysis [53] has been recasted under the assumption of pure χ + χ − production for a triplet dark matter model, which corresponds to our scalar scenario. We reinterpret their results and derive a lower bound of m c 130 GeV for 5 × 10 −5 < µ/v < 0.02 in the JHEP02(2019)140 The mass difference ∆m c is varied within the range that is compatible with the observed dark matter abundance. scalar scenario and 7 × 10 −4 < µ/v < 0.07 in the pseudo-scalar scenario. The excluded parameter region is shown in figure 6 as a green band.
In case of a normal mass hierarchy, χ h decays are prompt for µ/v 0.02 (0.07) in the scalar (pseudo-scalar) scenario. Soft di-leptons are now produced from both pp → χ + χ − and pp → χ + χ h . Adding the contribution from the latter process enhances our signal rate. However, our pp → χ + χ h contribution favors a different kinematic regime and produces only two charged leptons in the final state, compared to three leptons in the experimental analyses. We therefore do not expect a significantly stronger bound than m c 130 GeV when including χ + χ h contributions. The sensitivity could be enhanced by optimizing the signal region for signatures with exactly two soft leptons in the final state. All in all, however, searches for prompt soft leptons are very limited in their mass reach, due to the small production rates and large backgrounds. Searches for displaced and prompt leptons or displaced b-jet pairs are expected to be sensitive to much higher mediator masses.
In figure 7, we summarize the various displaced signatures that are predicted from long-lived mediator decays at the LHC. They are classified according to the nominal decay length of the mediator χ + (salmon) or χ h (purple) in both scenarios. Colored areas correspond to regions of experimental sensitivity with the ATLAS and CMS detectors. Hatched regions could be probed with an extended sensitivity. All other edges are theory bounds, which were explained above for each individual signature. Disappearing charged track searches have already excluded the blue and green parameter regions in the scalar and pseudo-scalar scenario, respectively. It is apparent that our model predicts signatures with different final states in basically all accessible layers of the LHC detectors. The search for dark matter from a small Higgs portal is thus most efficiently done by gathering all these signatures in a combined interpretation.

JHEP02(2019)140 7 Conclusions
In this work, we have investigated singlet-triplet fermion dark matter interacting with the standard model through a scalar or pseudo-scalar Higgs portal. The nature of the Higgs portal also has implications on the weak charged currents, which are vector-like in the scalar scenario and axial-vector-like in the pseudo-scalar scenario. Both scenarios are related through a chiral rotation of the fermion singlet. Besides the different Lorentz structure of the Higgs and W -boson couplings, the fermion mixing in the pseudo-scalar scenario is generically smaller than in the scalar scenario. This leads to different lifetimes of the mediators, with observable consequences for the dark matter and collider phenomenology.
Due to the strong bounds on the portal coupling from direct detection experiments, the thermal relic dark matter abundance relies on co-annihilation with mediators. In the regime of very small portal couplings, co-scattering and mediator decays have a crucial impact on the dark matter number density during freeze-out. For a reliable prediction of the relic abundance, co-scattering and decay processes thus need to be taken into account when solving the coupled system of Boltzmann equations. Since the phase of co-scattering is a common prediction in models with small portal couplings, we suggest to include coscattering and mediator decays in existing automated tools for relic density calculations.
Thermal dark matter with tiny Higgs couplings and a compressed dark sector implies long-lived mediators. At the LHC, this leads to a plethora of signatures with both prompt and displaced vertices, as well as disappearing tracks. Due to the small electroweak production rates and the softness of the visible decay products, searches for prompt signatures are limited in their mass reach. Current searches for prompt soft leptons cannot probe mediator masses above about 150 GeV. Displaced signatures, in turn, leave clear signatures that can be distinguished from the background with a handful of events. Disappearing charged track searches already exclude mediator masses up to m c ≈ 460 (480) GeV for tiny portal couplings of µ/v 10 −6 (10 −5 ) in the scalar (pseudo-scalar) scenario. Parameter regions with larger couplings can be probed with signatures of displaced b-jet pairs, displaced soft lepton pairs, or one displaced and one prompt lepton. Existing searches for displaced leptons probe portal couplings µ/v 3×10 −5 (3×10 −4 ). They exclude mediator masses up to m c ≈ 200 GeV, but are very limited in their sensitivity, due to the cut on the lepton transverse momentum. Lowering the momentum threshold will strongly enhance the sensitivity of displaced lepton searches to soft decay products, which are typical for compressed dark sectors. Notice that the parameter region with mediator masses not much larger than m c ≈ 150 GeV and larger portal couplings can be probed with prompt leptons, as well as displaced and prompt leptons and displaced b-jet pairs. A potential discovery of one of these signatures in this region can thus be confirmed by a complementary search for the other signature.
Our analysis shows that mediator masses of several hundred GeV can be accessible at the ATLAS and CMS detectors using the run-II data set. With the data set expected at the HL-LHC, the reach can be extended to probe mediator masses up to the TeV scale and even beyond with future colliders [29,47,[54][55][56]. While neither indirect nor direct detection will be able to test dark matter scenarios with tiny Higgs-portal couplings in the foreseeable JHEP02(2019)140 future, collider signatures are perfectly characteristic probes. Displaced soft objects in association with missing energy are sensitive to the mass splitting and the coupling of the dark sector at the same time. We enthusiastically encourage the ATLAS and CMS collaborations to exploit the lifetime frontier in the search for Higgs-portal dark matter.