Benchmarking the Inert Doublet Model for e+e− colliders

We present benchmarks for the Inert Doublet Model, a Two Higgs Doublet Model with a dark matter candidate. They are consistent with current constraints on direct detection, including the most recent bounds from the XENON1T experiment and relic density of dark matter, as well as with known collider and low-energy limits. We focus on parameter choices that promise detectable signals at lepton colliders via pair-production of H+H− and HA. For these we choose a large variety of benchmark points with different kinematic features, leading to distinctly different final states in order to cover the large variety of collider signatures that can result from the model.


Introduction
The results from the ATLAS and CMS collaborations from Run I and ongoing Run II are in good agreement with the predictions of the Standard Model (SM) [1,2]. Although the discovered Higgs particle appears to be consistent with the expectations for a SM Higgs boson, both the experimental uncertainties and theoretical speculations still leave room for new physics. In particular the scalar sector can provide intriguing scenarios in this respect which should be further scrutinized.
Although experimental collider data is in good agreement with predictions of the Standard Model alone, a number of non-collider observations can only be described in models containing additional (new physics) constituents. A prime example for this is dark matter (DM). Within the standard model of cosmology, the Planck mission data [3] implies that nearly 85% of the total matter content in the universe is dark. However, so far only the gravitational interactions of these hypothetical particles have been detected, and a fundamental nature of DM remains largely unknown. Since the Standard Model of elementary particles does not contain a viable DM candidate, any evidence of DM in the direct detection or indirect detection experiments or production at colliders would be a signal of new physics, the discovery of which is arguably one of the most important goals in the field.
An intriguing extension of the SM scalar sector is the Inert Doublet Model (IDM) which features a dark matter candidate [4][5][6]. In this two Higgs doublet model a discrete Z 2 JHEP12(2018)081 symmetry (called D-symmetry) is imposed, with the following transformation properties: (1.1) where the φ S doublet plays the same role as the corresponding doublet in the SM, providing the SM-like Higgs particle. This doublet is even under the D-symmetry, while the second doublet, the inert (or dark) φ D , is D-odd and contains four scalars, two charged and two neutral ones, labelled H ± and H, A, respectively. The above symmetry renders the additional SU(2) L doublet φ D inert, i.e. prevents its couplings to the SM matter sector, thereby providing a dark matter candidate. In the rest of this work, we consider cases where H serves as the dark matter candidate of the model. The IDM was first discussed in [4] and later in [5,6]. The model was further studied in [7][8][9][10][11], followed with further analyses of the IDM at colliders .
We here present a set of the IDM benchmarks proposed for detailed studies at the future e + e − colliders ILC and CLIC. They have been selected from updates of the scan presented in [23,30] and represent distinct features for two prominent production processes at linear colliders, e + e − → H + H − and e + e − → AH. Our benchmarks are designed to cover all interesting parameter space, featuring different mass splittings between H and other dark particles, leading to distinct collider signatures.
The outline of this paper is as follows. We start with a short description of the IDM in section 2, followed with a discussion of theoretical constraints and current experimental limits in section 3. In section 4 we define benchmark points that pass all constraints and discuss their applications to further studies at the first stage of linear colliders running at 250, 380 and 500 GeV center-of-mass energies. Benchmarks suitable for testing at high energy stages of 1, 1.5 and 3 TeV linear colliders are presented in section 5. We also comment on the impact of future XENON-nT measurements and prospects of testing the IDM at the LHC in section 6.

The IDM
The scalar sector of the IDM consists of two SU(2) L doublets of complex scalar fields, φ S and φ D , with the D-symmetric potential: (2.1) Exact D-symmetry implies that only φ S can acquire a nonzero vacuum expectation value (v). As a result the scalar fields in φ D do not mix with the SM-like field from φ S , and the lightest particle of the dark sector is stable. The dark sector contains four new particles: H, A and H ± . We here choose H to denote the dark matter candidate. A priori, any of the new scalars can function as a dark matter candidate. However, we neglect the choice of a charged dark matter candidate, as these are strongly constrained [32]. As inert scalars do not couple to fermions, it is not possible to assign a definite CP-property to them.

JHEP12(2018)081
Therefore, choosing A instead of H to be the lightest particle changes only the meaning of λ 5 , with rephasing of λ 5 → −λ 5 , but not the overall phenomenology of the model, cf. appendix A in [23]. After electroweak symmetry breaking, the model contains seven free parameters. Agreement with the Higgs boson discovery and electroweak precision observables fixes the SM-like Higgs mass M h and v, and we are left with five free parameters, which we take as where the λ's refer to couplings within the dark sector and to the SM-like Higgs respectively. In the following, we will use the abbreviation λ 345 = λ 3 + λ 4 + λ 5 .

Experimental and theoretical constraints
In this work, we make use of the tool chain used in [23], and explicitly follow the scan procedure described therein; however, we update several constraints as discussed below. Explicit benchmark points (BPs) for the IDM, taking all constraints viable at that time into account, have been presented in [23,33], which focus on processes at the LHC; some of these were already investigated in a linear collider context in [24]. Ref. [30] shows how the available parameter space for certain scenarios is further limited by more recent constraints.
As the experimental constraints have evolved significantly since that time, we decided to define the new set of benchmark points, fulfilling the updated constraints and focused on the detailed analysis of e + e − collider sensitivity. Unless stated otherwise, all considered BPs fulfil the latest experimental limits; following [23] (see also the discussion in [34]), we do not require the IDM to provide 100% of the dark matter relic density. Below we briefly summarize the imposed constraints, emphasizing on updates with respect to [23], and describe the set-up to find good BPs and discuss the obtained limits.

Theoretical and experimental constraints
Positivity constraints: we require that the potential is bounded from below, therefore no field configuration leads to V → −∞, resulting in tree-level relations [35] (3.1) These relations hold on tree level and in this work we do not consider higher-order contributions, which in principle could lead to change in stability of the electroweak vacuum [10,36].
Perturbative unitarity: we require the scalar 2 → 2 scattering matrix to be unitary, i.e. all eigenvalues of scattering matrices for scalars with specific hypercharge and isospin should satisfy |L i | ≤ 16 π [37,38]. Furthermore, we require all quartic scalar couplings to be perturbative, i.e. to take absolute values ≤ 4 π.

JHEP12(2018)081
Global minimum: in the IDM two neutral minima can coexist even at tree level. Unless the following relation is satisfied the inert minimum is only a local one, with the global vacuum corresponding to the case of massless fermions [39]. We impose the above relation in our scan.
Higgs mass and signal strengths: the mass of the SM-like Higgs boson h is set to in agreement with limits from ATLAS and CMS experiments [40], while the total width of the SM-like Higgs boson obeys an upper limit of [41] Γ tot ≤ 9 MeV. In the IDM, the total width of the SM-like state can obtain modifications from the following two contributions. For dark matter masses M H ≤ M h /2, invisible decays of the 125 GeV resonance can lead to large additional contributions. Therefore, in these scenarios the above bound poses one of the most dominant constraints, especially affecting λ 345 [23,30]. Furthermore, the partial decay width of h to diphoton final states can be altered significantly [11,19], as the new physics corrections are formally of the same order as the SM process. This leads to a clear distinction between allowed and forbidden regions in the (λ 345 , M H ± ) plane [23,30]. The Run I combined ATLAS and CMS limit for h → γγ signal strength is given by µ γγ = 1.14 +0. 38 −0.36 [42]. In our analysis we use both the upper limit (3.3) and require agreement within 2 σ for the prediction of h → γ γ. We furthermore check agreement with all other branching ratios of the 125 Higgs on the 2 σ level using the publicly available tool HiggsSignals-2.2.1beta [43], and require ∆χ 2 ≤ 11.31, corresponding to a 95% confidence level.
Gauge bosons width: introduction of light new particles could in principle significantly change the total width of electroweak gauge bosons (cf. e.g. [44]). To ensure that W ± → HH ± and Z → HA, H + H − decay channels are kinematically forbidden we set: Electroweak precision tests (EWPT): we call for a 2 σ (i.e. 95% C.L.) agreement with electroweak precision observables, parametrized through the electroweak oblique parameters S, T, U [45][46][47][48]. In our work, calculations were done through the routine implemented in the Two Higgs Doublet Model Calculator (2HDMC) tool [49], which checks whenever model predictions fall within the observed parameter range [50].
Charged scalar mass and lifetime: we take a conservative lower estimate on the mass of M H ± following analysis in [51] to be

JHEP12(2018)081
We also set an upper limit on the charged scalar lifetime of τ ≤ 10 −7 s, (3.6) in order to evade bounds from quasi-stable charged particle searches. This translates to a lower bound on the total decay width of the charged scalar H ± of Γ tot ≥ 6.58 × 10 −18 GeV. In [52], the authors show that limits from a recast of heavy charged particle searches at the LHC [53] can set bounds on the lifetimes of H ± between 10 −9 and 10 −5 s, with the lower value prevailing for charged masses M H ± 500 GeV. For most of the parameter space, the limit (3.6) therefore corresponds to a conservative upper limit.
Collider searches for new physics: we require agreement with the null-searches from the LEP, Tevatron, and LHC experiments. We use the publicly available tool HiggsBounds-5.2.0beta [54][55][56][57]. In addition the reinterpreted LEP II searches for supersymmetric particles analysis exclude the region of masses in the IDM where simultaneously [12] M as it would lead to a visible di-jet or di-lepton signal. After taking into account all the above limits we are outside of the region excluded due to the reinterpretation of the supersymmetry analysis from LHC Run I [21].
Dark matter phenomenology: we apply dark matter relic density limits obtained by the Planck experiment [3]: For a DM candidate that provides 100% of observed DM in the Universe we require the above bound to be fulfilled within the 2σ limit. However, we also allow for the case where H is only a subdominant DM candidate, with (see [23,30,34]). In such a scenario, additional dark matter candidates would be needed in order to account for the missing relic density. In the results presented here, we apply XENON1T limits [58]. 1 These supersede previous bounds applied e.g. in [30] in relevant regions of parameter space, and are therefore crucial for a correct determination of the available parameter space, especially for low dark masses. Results from indirect detection experiments, e.g. Fermi-LAT [60], give less stringent constraints than collider and direct detection experiments discussed above [10,34]. A number of planned DM indirect detection experiments, mainly the Cherenkov Telescope Array, will be able to probe the heavy mass region, with DM particle heavier than 500 GeV [61]. Furthermore, if the reported gamma-ray excess from the Galactic center is of DM origin [62], it can be explained by the IDM with dark matter masses near the Higgs resonance or around 72 GeV [63]. All dark matter variables were calculated with the use of micrOmegas version 4.3.5 [64].

Scan setup and limits
After fixing the value of M h , and hence both λ 1 and m 2 11 , we are left with five independent input parameters for the scan: three masses of dark scalars M A,H,H ± and two couplings, λ 345 , λ 2 . In the initial setup of our scan, masses take values between 0 and 1 TeV, with M H always being the lightest and M H ± ≥ 70 GeV. Unless stated otherwise, scalar couplings fall in the range of λ 2 ∈ [0; 4.5], λ 345 ∈ [−1.5; 4π].
In order to get interesting benchmark points we follow the procedure described in [23]. All constraints described in the previous section are checked in steps, with the aid of publicly available tools. In the process we track the impact of each exclusion criterion.
The first step contains a check through the 2HDMC in order to establish agreement with theoretical constraints (positivity, stability, perturbative unitarity). The same code is used to check the SM-like Higgs and electroweak gauge bosons widths, the decay rates of h → invisible, and h → γγ, properties of charged scalar (lifetime, mass) as well as the EWPT observables. Points that have passed the first step are then checked against limits from collider searches, in particular the Higgs signal strength limits, with the use of HiggsBounds and HiggsSignals. Points which passed the collider test are then confronted with limits from DM phenomenology, i.e. the relic density constraints with upper Planck bound (points that correspond to 100% of measured Ω c h 2 are selected at later stage) and direct detection limits from XENON1T, with the use of micrOmegas.
Although the IDM is one of the simplest extensions of the SM and has a limited number of parameters, it is still difficult to establish which constraint has the greatest impact on excluding a given point in parameter space. In the following, we point to some generic features that however can be used for a clear distinction in the model parameters space (see also [23]).
It is important to emphasize that the couplings that govern the production and decay processes at e + e − colliders are mainly determined by the electroweak parameters of the SM; the additional parameters in the potential do not play any significant role for the allowed scenarios. On the other side, the parameter λ 345 is especially sensitive to constraints from dark matter observations. This nicely demonstrates the important complementarity of collider and astrophysical measurements in constraining the IDM parameter space.
We divide the discussion in two different subsections (i) constraints on the masses of the dark scalars, with the constraint given in eq. (3.9) for the dark matter relic density; (ii) additional limits we obtain when requiring exact relic density.
We only briefly discuss the second point here and refer to the literature (see e.g. [34]) for further details.

Limits on masses
The collider phenomenology of the IDM is mainly determined by the dark scalar masses, as all relevant production and decay channels are governed by electroweak couplings and the JHEP12(2018)081 corresponding SM parameters. As dominant constraints for the regions for M H ≤ M h /2 and M H ≥ M h /2 originate from different sources, we will discuss these separately.
• In general, constraints arising from EWPT are of a great importance to our studies. As found in [23,34,65], mass splittings between inert particles are heavily limited by the S,T,U parameters. First, only moderate mass splittings are allowed by EWPT data, with the preferred value of M H ± −M A below 100 GeV. Also, there is a hierarchy between masses, with the charged scalar being the heaviest particle. The reverse relation is not excluded, however it will lead to a larger tension with combined S and T limits. In general, points with moderate mass splittings are preferred, especially for dark matter candidates with masses ≥ 300 GeV, where splittings between M A and M H are typically of order 10 % or lower. For masses M H ≤ 100 GeV, we found that relatively large mass differences are allowed between the two neutral scalars.
The mass hierarchy and constraints on mass splittings can influence cascade decays of inert particles.
• In addition to eq. (3.5), the measurement of h → γγ puts a lower limit on the charged Higgs mass as a function of λ 345 , cf. e.g. [30].
• Masses of DM particles below 45 GeV are excluded [23,34]. For M H ≤ M h /2 limits from the invisible branching ratio of the SM Higgs lead to relatively low values of λ 345 , which in turn result in the relic density exceeding the measured value by orders of magnitude or the mass splittings changing electroweak gauge bosons widths 2 significantly.
• For dark scalars with masses 100 GeV, additional specific constraints are given by eq. (3.7).
The parameters λ 2 and λ 345 are also constrained from combined positivity, unitarity and global minimum conditions [66]. The exact value of λ 2 parameter would matter in studies that consider interaction between inert particles, for example for the astrophysical implications of self-interacting dark matter, or if the loop processes like HH → γγ through the H + H − loop are considered, this is however beyond the scope of this work.
Dark matter constraints give major bounds on allowed values of λ 345 . In general, |λ 345 | 1.5 for dark masses up to 1 TeV. For M H M h /2, this limit decreases to |λ 345 | 0.006, due to the inclusion of results presented in [58]. Limits from Higgs signal strength measurements also limit this parameter, albeit less strongly than dark matter constraints. If we allow for a subdominant dark matter candidate, the parameter space largely opens up for regions where the dark matter relic density is much lower than the Planck value. Direct detection limits are then rescaled and therefore considerably relaxed (see also [23]). A prominent example for this is the region where M H ∼ M h /2, leading to large annihilation cross sections. This broadens especially the allowed range for λ 345 .
In summary, M H ± and M A are relatively degenerate throughout the allowed parameter space. The splittings between these and the dark matter candidate depend on the DM mass JHEP12(2018)081 and can become quite large, especially for low M H . The couplings λ 2 and λ 345 do not play a significant role in collider phenomenology.

Requiring exact relic density
Requiring the model to render exact relic density, as specified by eq. (3.8), puts additional constraints on the model. These are mainly on the coupling combination λ 345 , but also affect the values of possible dark scalar masses in such a scenario.
• If H is the only source of DM in the Universe then it is possible to find good points only for masses 55 GeV M H 75 GeV, where the allowed range of dark scalar masses is determined by electroweak constraints, without additional fine-tuning.

• In the region 75 GeV
M H 500 GeV DM direct detection searches exclude all points that have exact relic density. Points that have passed all other constraints provide only a subdominant dark matter candidate.
• For M H larger than 500 GeV extreme fine-tuning of the dark masses is required in order to obtain exact relic density; cf. the discussion in [34]. The authors of that reference find that O (GeV) mass splittings between the dark scalars would be required in order to obtain the correct relic density; in turn, this can lead to long-lived particles at the LHC, where minimal allowed mass splittings are O (0.2 GeV) [52]. In our work, we do not study this particularly fine-tuned region in more detail.

Benchmark points
As the aim of this note is to present benchmarks useful for further studies at linear colliders, we are particularly interested in points that provide an observable signal. For the possible signal we take the pair-production processes e + e − → H + H − and e + e − → AH (4.1) for charged and neutral scalar production. The s-channel production e + e − → h → AA is also possible, however it is suppressed by a small electron Yukawa coupling. The dark scalars H ± and A will then decay into a virtual or on-shell electroweak gauge boson and the dark matter candidate. Here, we first concentrate on center-of-mass energies of √ s ∈ {250; 380; 500} GeV, therefore constraining ourselves to scenarios where either M H + ≤ 250 GeV or M H + M A ≤ 500 GeV (or both). In section 5, we extend the analysis to contain higher dark scalar masses, in order to investigate the collider reach of the high energy ILC at 1 TeV and higher energy stages of CLIC at 1.5 TeV and 3 TeV [67]. In this work, we consider benchmarks with both on-shell and off-shell intermediate gauge bosons, as these differ in the collider phenomenology, and different cut strategies have to be applied.
Our benchmark selection was done in the following steps: 1. Benchmark candidates were generated by employing the scan presented in [23], with updated experimental constraints as discussed above. We selected about 7500 bench-

JHEP12(2018)081
mark candidates which fulfill all constraints (we allow for under-abundant dark matter, therefore relaxing bounds from direct detection considerably in certain regions of parameter space).
2. We then calculated the production cross sections at 250, 380 and 500 GeV center-ofmass energies for processes (4.1); we required a minimal production cross section of 10 fb to classify a point as "accessible" for a certain process/energy stage; 3. if possible, we required a high-enough relic density, providing at least 50% of the value observed by the Planck collaboration; 4. finally we selected benchmark points corresponding to different accessibility at the subsequent energy stages and different kinematical configurations, namely on-shell vs off-shell intermediate gauge bosons.
The calculation of cross sections was performed using input files generated by SARAH 4.13.0 [68][69][70] and SPheno 4.0.3 [71,72], which were passed to WHizard 2.2.8 [73,74]. Initial state radiation was taken into account but not the beam luminosity spectra. Point 2 above allows us to categorize different classes of benchmarks, which we label "XYZ" in the following, with X=1(2) corresponding to 1 (2) production processes accessible at 250 GeV, while Y and Z are defined accordingly for 380 GeV and 500 GeV. Within a certain category, we then considered different mass splitting configurations, in order to cover all possible typical parameter configurations leading to distinct collider signatures that can be generated by the IDM. The final selection of benchmark points, focusing on kinematic properties and selection criteria, is given in table 1. The table also contains the complete set of independent parameters for each point, as well as the relic density.
In our selection, BP1 is an example for a relatively small mass splitting, that forces the intermediate gauge boson to be off-shell. Other benchmark points, as e.g. BP9, allow for on-shell decays in both channels. However, these force the corresponding dark scalar masses to be O (150 GeV), therefore leading to smaller production cross sections.  figure 1 (red points) cover mass gaps up to about 250 GeV only, due to the required minimal cross section (see point 2 above). Notice that for most selected benchmark points the DM candidate is relatively light, with a mass below 80 GeV. Tables 2 and 3 show the production cross sections at various center-of-mass energies for all benchmark scenarios. These indicate promising prospects of detection at future linear colliders. Because A is always a lighter particle, neutral channels are usually accessible at lower energies than charged ones. If both channels are accessible, the cross section for charged scalar pair-production is usually larger than for the neutral ones.

Extending to higher mass scales
So far, we have discussed IDM scenarios that are accessible at center-of-mass energies up to 500 GeV, with a mass range of M 250 GeV for all dark scalars. However, the parameter space of the IDM largely opens up for higher dark scalar masses, especially as direct detection constraints are less strict and direct collider searches pose less stringent constraints. 3 Therefore, in this section we consider high-mass benchmark points (HP) that can be explored at higher center-or-mass energies, e.g. at 1 TeV ILC or at the high-energy stages of the CLIC collider, at 1.5 TeV or 3 TeV.
As before, all benchmark points have passed theoretical and experimental constraints discussed in section 3. We now consider the nominal collider energies of 1 TeV, 1.5 TeV and 3 TeV. From about 6000 parameter points, we selected 20 benchmark points. The selection was done analogously to the low-energy case, aiming in addition to cover the whole mass range of the dark scalars up to 1 TeV. Benchmark points are summarized in table 4 and the production cross sections for the considered IDM scalar pair-production processes in table 5. The maximum mass splitting between the dark scalars is about 140 GeV, as the relic density decreases rapidly with the increasing mass difference (see section 3.2.2).
M  Table 2. Production cross sections in fb for on-shell charged scalar pair-production, e + e − → H + H − , for the center-of-mass energies considered in this work. We only list benchmark points with at least one non-zero production cross section. We also display the branching ratio H + → W + H (the other possible channel, H + → W + A, is suppressed for most BPs). The displayed branching ratios were calculated using 2HDMC.

Prospects for future dark matter and LHC experiments
Experiments at future e + e − colliders will be able to probe the described benchmark scenarios on the time scales of ten to twenty years. By that time, the sensitivity of direct DM search experiments will have improved significantly, and much larger samples of data will also be collected by the LHC experiments. It is therefore an important question whether the benchmarks presented here are already accessible at the LHC or via dark matter direct detection. As already described above, significant constraints of the IDM scenarios are set by the recent XENON1T measurement [58]. The experiment continues to collect data and it is expected to improve its sensitivity by a factor of about 3 [75]. In this range, most of the points providing dominant contribution to the DM density in the Universe (see tables 1 and 4) can be probed. Four of these scenarios (BP6, BP16, BP21 and HP10) will remain JHEP12(2018)081  Table 3. Production cross sections in fb for on-shell neutral scalar pair-production, e + e − → HA, for the center-of-mass energies considered in this work. We only list benchmark points with at least one non-zero production cross section. The branching ratio BR(A → HZ ( ) ) ≈ 100%.
inaccessible at XENON1T, but can be probed at XENONnT, when the sensitivity will be enhanced by another order of magnitude. Only one scenario, BP2, will not be probed after the future XENON detector upgrade. Scenarios resulting in subdominant contribution to the relic density are less constrained by direct search experiments. Still, half of the presented benchmark points can be probed with the full XENON1T data and only four scenarios (BP18, BP19, BP20 and HP3) remain inaccessible at XENONnT. This shows that direct DM search experiments will continue to set important constraints on the proposed class of models and can result in the need to redefine the set of benchmark points in the future. On the other hand, if an excess in direct detection is observed in the XENON detector, the IDM scenarios will provide a perfect test framework for interpretation of these results. Direct searches will only indicate the possible mass and coupling range for the lightest IDM scalar. Other parameters of the model need to be constrained in collider experiments. Searches for pair production of IDM scalars have also been considered for current and future runs of the LHC. Many different final states can be considered, including mono-jet, mono-Z, mono-Higgs and Vector-Boson-Fusion + missing transverse energy signatures [25,26,29,34,76]. Especially multi-lepton final states [26] are promising. 4 In addition, also multi-jet final states [25,76] and combinations [29] have been considered. For example, in [26,29,76] the authors present phenomenological studies that render scenarios with dark masses 300 GeV accessible at the HL-LHC. Regarding other channels, the analysis of mono-jet signature considered in [31,34] indicates that LHC has limited sensitivity to probe the IDM. Additional assumptions are needed to set the scalar mass limits and they range only up to about 200 GeV for HL-LHC. Also for the di-jet plus missing transverse energy signature, as studied e.g. in [25], LHC sensitivity is significantly affected by the large background, which can not be sufficiently  Table 5. Production cross sections in fb for high-mass benchmark points at 1 TeV, 1.5 TeV and 3 TeV, for the production processes considered here.

JHEP12(2018)081
suppressed even with strong kinematic selection. For the benchmark scenarios considered, the maximum signal significance at an integrated luminosity of 3000 fb −1 was about 2σ. However, we want to state that these analyses have not been on the same level as the current study, as they did not make use of more advanced analysis techniques. In general, production cross sections can be in the O (500 − 700fb) range for a center-of-mass energy of 13 TeV [77]. Without detailed analyses, projections of reachability are however difficult to make. We therefore strongly encourage the experimental collaborations to investigate the benchmark points presented here at current and future LHC runs.

Conclusions
In this paper we have revisited and updated the available parameter space of the Inert Doublet Model. The model features an exact Z 2 symmetry which results in dark scalars that do not interact with SM fermions, and the lightest neutral scalar can serve as a JHEP12(2018)081 promising dark matter particle. We took into account most recent experimental constraints from relic density and direct dark matter searches, including the latest XENON1T 2018 results, as well as collider bounds and theoretical constraints. Based on these updated results, we have provided benchmark scenarios accessible at the initial stages of future linear e + e − colliders (250 and 500 GeV ILC and 380 GeV CLIC) as well as benchmarks that can be tested at high-energy stages (1 TeV ILC and 1.5 and 3 TeV CLIC). In doing so we pursued the philosophy of covering the widest range of parameters and experimental signatures. We provide predictions of production cross sections at these energies, and supplement these with information about the branching fractions of the relevant decay modes. We encourage the LC groups to make use of these benchmark scenarios. Although the benchmarks have been defined with e + e − physics in mind, we strongly encourage our LHC experimental colleagues to consider these scenarios in the analysis of the current and upcoming LHC data.