Prospects for direct searches for light Higgs bosons at the ILC with 250 GeV

The particle discovered in the Higgs boson searches at the LHC with a mass of about 125 GeV is compatible within the present uncertainties with the Higgs boson predicted in the Standard Model (SM), but it could also be identified with one of the neutral Higgs bosons in a variety of beyond the SM (BSM) theories with an extended Higgs sector. The possibility that an additional Higgs boson (or even more than one) could be lighter than the state that has been detected at 125 GeV occurs generically in many BSM models and has some support from slight excesses that were observed above the background expectations in Higgs searches at LEP and at the LHC. The couplings between additional Higgs fields and the electroweak gauge bosons in BSM theories could be probed by model-independent Higgs searches at lepton colliders. We present a generator-level extrapolation of the limits obtained at LEP to the case of a future e+e-\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$e^+e^-$$\end{document} collider, both for the search where the light Higgs boson decays into a pair of bottom quarks and for the decay-mode-independent search utilising the recoil method. We find that at the ILC with a centre-of-mass energy of 250 GeV, an integrated luminosity of 500 fb-1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^{-1}$$\end{document} and polarised beams, the sensitivity to a light Higgs boson with reduced couplings to gauge bosons is improved by more than an order of magnitude compared to the LEP limits and goes much beyond the projected indirect sensitivity of the HL-LHC with 3000 fb-1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^{-1}$$\end{document} from the rate measurements of the detected state at 125 GeV.


Introduction
The properties of the Higgs boson that was discovered in 2012 are in agreement with the predictions of the Standard Former address: P. Drechsel. a e-mail: gudrid.moortgat-pick@desy.de (corresponding author) b e-mail: georg.weiglein@desy.de Model (SM) within the current experimental accuracy, but they are also compatible with a wide variety of extensions of or alternatives to the SM. Extended Higgs sectors predict the existence of additional Higgs bosons, which could be heavier but also lighter than the observed state at 125 GeV. The coupling of the SM Higgs boson to the gauge bosons W and Z is such that terms with a bad high-energy behaviour in longitudinal vector-boson scattering exactly cancel with each other. As a consequence, in an extended Higgs sector comprising the SM-like state at 125 GeV and additional Higgs bosons the couplings of the additional neutral Higgs bosons to W and Z , g φ i V V , are expected to be small. This implies that additional heavy neutral Higgs bosons may not be detectable via the search channels φ i → Z Z, W + W − and that for an additional light Higgs boson (or more than one) g φV V may be so small that such a light Higgs boson would have escaped the limits from the Higgs searches at LEP [1,2] (see also Ref. [3]) and the Tevatron [4].
The case of an extended Higgs sector containing a SMlike state that can be identified with the observed Higgs signal and further Higgs bosons of which at least one is lighter than 125 GeV can be realised, for instance, in a general two-Higgsdoublet model 1 (2HDM), and it occurs generically in extensions with a light singlet such as the next-to-minimal supersymmetric extension of the SM (NMSSM), see e.g. Refs. [7][8][9], or the N2HDM, a 2HDM with an additional real Higgs singlet, see e.g. Ref. [10]. It should be noted that extensions of the SM with a singlet-dominated state in the mass range around or just below 100 GeV are also of interest in view of the observed local excesses around 96 GeV at the 2-3 σ level in both the searches at LEP in the e + e − → Z φ, φ → bb channel [1,2] and at CMS in the φ → γ γ searches [11] (the CMS result is compatible [12] with the results of the corresponding searches at ATLAS [13]). Possible interpreta-tions have been discussed e.g. in the NMSSM [9,14,15], an inflation-inspired μNMSSM [16,17] and the N2HDM [10].
Specifically, the neutral Higgs bosons φ i of extended Higgs sectors consisting of any number of doublets and singlets fulfill the sum rule at lowest order that the squared couplings to gauge bosons of all φ i add up to the squared coupling of the SM Higgs boson to gauge bosons, and 100 GeV appears to be particularly promising, as this mass range is only mildly constrained by the existing limits from the Higgs searches at LEP and the φ → γ γ searches at the LHC (and as discussed above, in both types of these searches an interesting excess above the background expectation has been reported).
In the current paper we study the sensitivity of searches at the ILC in its first stage of √ s = 250 GeV centre-of-mass energy for light additional Higgs bosons with masses below the one of the observed signal at 125 GeV. Due to the clean environment and consequently a very favourable signal-tobackground ratio for the prospective integrated luminosity, the low beamstahlung, the precise knowledge of the beam energy and the availability of polarised beams, the ILC has a high physics potential in the direct search for such light additional Higgs bosons. We perform a generator-level extrapolation of the limits obtained at LEP to the ILC case, both for the search for the φ → bb final state and for the decaymode independent recoil technique (the latter was pioneered at LEP by the OPAL collaboration [3] and forms the basis for the total cross-section measurement for Zh(125) production at future e + e − colliders).
The two main production processes for a neutral Higgs boson at the ILC are Higgs-strahlung (e + e − → (H/φ)Z ), dominant at lower masses and lower collider energies, and W W -fusion (e + e − → (H/ )νν), dominant at higher Higgs masses and higher collider energies. In the current study, we therefore focus on the Higgs-strahlung process, for which we investigate the φ → bb channel as well as the recoil method, where only the leptonic Z -boson decays are utilised for the reconstruction of the final state (in our study we have concentrated on the Z → μ + μ − decay). We first validate our approach with the LEP results and then perform an extrapolation to the case of the ILC. The present paper builds up on a preliminary study that was carried out in Ref. [18], see also Ref. [19].
The paper is organised as follows. The adopted statistical method is described in Sect. 2. In Sect. 3 the method is validated with the results obtained at LEP, and in Sect. 4 the sensitivity of the ILC searches for light Higgs bosons is discussed. Our conclusions are given in Sect. 5.

Description of the applied methods
We consider the "Higgs-strahlung" process Z e e Z H/φ (2) where for our generator-level analysis we focus on the Zboson decay into muons, Z → μ + μ − , and treat the cases where either the decay of the scalar H/φ → bb is reconstructed, or only the information from the leptonic Z -boson decay is used. For the scalar either the SM Higgs boson H is considered or a light Higgs boson φ with a reduced coupling to gauge bosons compared to the case of a SM-like Higgs boson with the same mass.
Concerning the statistical treatment we follow the prescription that was outlined for S 95 in Refs. [1,20], using a simplified approach, see also [21]. We have generated event samples corresponding to the two hypotheses "all events are generated by the background only" or "all events are generated by the background plus a hypothetical signal". This has been done by comparing the events for the full e + e − → 4 fermion process for the considered mass value of the assumed signal, taking into account all interference contributions, with the events for the full e + e − → 4 fermion process for the case where the mass value of the assumed signal is beyond the kinematic reach of the collider. The quantity gives an upper limitσ on a cross section that is compatible with the "background only" hypothesis at the 95% confidence level, normalised to a reference cross section σ ref , or equivalently on the ratio of the allowed signal raten, normalised to the reference signal rate n. As reference process we use in Eq. (3) the case of the SM Higgs boson at the considered mass value. Accordingly, Eq. (4) can be interpreted as the ratio between the squared couplings of the scalars φ and H to the Z -boson, for each mass value. We have generated event samples for the process with the Monte Carlo generator Whizard-2.4.1 [22,23], which we apply for both signal processes and their reference values (with and without reconstruction of the H/φ → bb decay) as well as for the corresponding backgrounds (a simple cut has been applied to control backgrounds with photon radiation). Since we generate the full e + e − → 4 fermion process, the corresponding signal and background events are not generated separately. The obtained events are classified by either the invariant mass of the quark system or the recoil mass of the lepton system for this analysis. The signal rate in the i-th bin (we use a width of 1 GeV for the mass bins) is determined by the difference of the event rates for the background plus signal (sb i ) and the background only hypotheses (b i ), and the signal rate s i is obtained by sb i − b i . Negative rates for signal events are cut, We assume that the number of potential signal events d i in each of the i bins are distributed according to a Poisson distribution with the expected values b i for the "background only" and s i +b i for the "background plus signal" hypotheses, respectively. The test statistic orders the outcome of test experiments according to their "signal likeness". While the expectation values of the Poisson distributions have to be determined a priori, the numbers of potential signal events d i have to be determined experimentally. For this work we consider only simulated events, and thus d i = s i . The logarithm of the test-statistic yields the weights for the number of potential signal events per bin as This definition requires the presence of at least one background event per considered bin. In order to accommodate this, we choose the luminosity for the simulation large enough that at least one event can be found in each bin, i.e. b i ≥ 1, and scale the luminosity afterwards with a factor c to the desired value. The weights w i are not affected by the scaling procedure.
With the weight factors of Eq. (8), the scaling factor S 95 given in Eq. (4) can be expressed [20] as where K denotes the number of standard deviations for the required significance: for S 95 the value of K = 1.96 is used, corresponding to a confidence level (C.L.) of 5% for the signal hypothesis. The variance σ sb is given by

Validation of the methods with results from LEP I and LEP II
In order to validate the described method we apply it in a first step to the search in the φ → bb channel and compare with the results obtained by the LEP combination [1,2]. Furthermore, in a second step, we validate our approach for the case of the search via the recoil method by comparing with the results obtained by the OPAL collaboration in their Higgs analysis employing the recoil method [3]. We simulated the process e + +e − → b+b+μ + +μ − with the SM implementation of Whizard 2.4.1 [22,23]. This process contains for the signal events the Higgs-strahlung process where the produced Higgs boson decays into bb and the Z boson decays into μ + μ − . As described above, for simplicity we restrict to the bbμ + μ − final state, i.e. we do not perform a separate analysis for the bbe + e − final state and for our analysis via the recoil method we also do not simulate additional decay modes of the Higgs boson. As explained above, our "background only" and "background plus signal" hypotheses are generated by appropriately adapting the mass of the implemented SM-like Higgs boson, i.e. for the "background only" hypothesis the mass m φ of the scalar is chosen beyond the kinematic accessibility of the LEP experiments.
For LEP 1 we consider the integrated luminosity of the four experiments at the centre-of-mass energy of 91.2 GeV. For LEP 2 we use for the comparison the integrated luminosity that the four LEP experiments recorded in a particular year and assign it to the highest energy that was achieved during that year, with the exception of the last year of LEP running. For the latter we use the energy of 206 GeV [24], where the bulk of the luminosity was recorded during that year. The luminosity per experiment that we have used in our comparison is given for the different LEP energy stages in Table 1. The number of events that we actually generated for our simulation is in fact 400 times higher than the LEP luminosity in order to facilitate the determinaton of the weight factors w i , see Eq. (8), and to reduce the statistical error. The results were then scaled down to the appropriate luminosities as described above.
In our generator-level analysis we do not take into account hadronisation effects of the b-quarks and we also do not simulate detector effects. This simplification leads to an overoptimistic estimate for the signal efficiencies and thus for the observed signal rates. Since for the searches making use of the φ → bb final state the signal efficiencies have not explicitly been given for the LEP experiments [21], we introduce a scale factor multiplying the luminosity that we use for the comparison with the LEP 1 results and we also allow for such a factor for the comparison of the LEP 2 results. While our over-optimistic treatment of the signal efficiencies implies the need to scale down our effective luminosities, we also need to appropriately scale them up to account for the fact that we have simulated only events for the Z → μ + μ − final state, whereas the LEP analyses also incorporate the other decay modes of the Z boson. Both effects can be combined into a single scale factor. We introduce such a factor, c bb , for the φ → bb analyses both at LEP 1 and LEP 2. The corresponding factors for the recoil method analyses at LEP 1 and LEP 2 are denoted as c recoil . Specifically, we obtain the effective luminosities L eff for our comparison as where L is the luminosity derived from the values given in Table 1. The scaling factors c bb , c recoil depend on the parameters of the experiment and the analysis. We have determined those factors from a comparison with the LEP-combined results for S 95 in the case of the φ → bb analyses and with the OPAL results for the case of the recoil method analyses. The four approximate scaling factors, rounded to integer values, that we have obtained in this way are listed in Table 2.
Obviously such a simple scaling factor can only roughly approximate an actual experimental analysis. For the case of the φ → bb analysis the scaling factor for incorporating besides the Z → μ + μ − final state also the events for the other Z decay modes would roughly correspond to an increase of the effective luminosity by a factor of 30 for an ideal detector. Assuming for simplicity a median signal efficiency of 50% for the signal rates leads to the conclusion that the coefficients c bb in Table 2 should not be larger than ≈ 15. The comparison with the listed results for c bb shows that the value that we obtained for the LEP 1 case (91.2 GeV) is rather close to the expectation for this idealised case. For the LEP 2 case (> 91.2 GeV), on the other hand, the detector effects play a larger role, and our generator-level extrapolation corresponds to a more optimistic estimate in comparison with the analysis incorporating a realistic treatment of the backgrounds and the experimental efficiencies. Table 2 shows that for the recoil method the required correction factors c recoil are smaller. This is expected from the facts that the OPAL analysis only made use of the decay modes Z → μ + μ − , e + e − , that for this analysis explicit signal efficiencies have been published, which we have taken into account in Eq. (10), and that generally the detector effects are expected to be less important for the recoil analysis as it only relies on a leptonic final state.

Comparison with LEP data using the H /φ decay mode
In order to compare our analysis for the search using the H /φ decay mode with the LEP data, we use the event samples for the process e + + e − → b +b + μ + + μ − and employ the information from the reconstruction of the bb system for the identification of the samples for "background only", b i , and "background plus signal", b i + s i . The information from the μ + μ − system is only used for validating that the signal events are compatible with the production of a Higgs boson together with a Z boson. A challenging region for this analysis is where the scalar mass m φ is close to M Z , since the presence of a large number of background events close to the Z -boson mass M Z weakens the limit on S 95 ; in this region we fit the expected background events to emulate a more detailed analysis: we extrapolate the expected background events in an interval around M Z and fit these points with a third-order polynomial to obtain the number of events in the bins around M Z as the value of the fit function at the central mass of each bin. The result using our method for approximately reconstructing the expected LEP limit is shown in Fig. 1 in com-Fig. 1 The expected limit from LEP for the analyses using the φ → bb decay information is compared with our method for approximately reconstructing the expected limit for S 95 from the process e + e − → Z (H/φ) → μ + μ − bb. Our approximate result makes use of the two scaling factors c bb given in Table 2 parison with the expected limit that was published by the LEP collaborations. As explained above, our method makes use of the two scaling factors c bb given in Table 2. The comparison in Fig. 1 shows that with those two factors as input our method reproduces very well not only the normalisation of the expected LEP limit but also its shape as a function of m φ .

Comparison with OPAL data using the recoil method
The OPAL collaboration also used the recoil method for analysing the data, i.e. exploiting the recoil of H, φ from the Z -boson and analysing the Z decay only. This method has the great advantage of being completely independent of the H, φ decay modes.
The OPAL analysis used the decay modes Z → μ + μ − , e + e − [3]. We restrict our analysis to the decay mode Z → μ + μ − only, but extrapolate in a second step the luminosity taken by OPAL to the full LEP luminosity, see Table 1 The events from the μ-lepton pairs are ordered by their total energy E i into bins, and for each bin we calculated the respective recoil mass m Fig. 2 The expected limit from OPAL based on the recoil method is compared with our method for approximately reconstructing the expected limit for S 95 , which makes use of the two scaling factors c recoil given in Table 2. Also shown is our result where the luminosity has been extrapolated to the full LEP luminosity to obtain the event rates for background, ). Since the OPAL collaboration published the signal efficiency of their analysis in Ref. [3], we could make use of this information to estimate a mean signal efficiency of ε = 30% for our approach and obtained in this way the effective luminosity L eff according to Eq. (10). The obtained values for the two scaling factors c recoil ( √ s) are given in Table 2. Our approximate reconstruction of the expected limit for S 95 from the recoil method analysis is compared with the expected limit published by the OPAL collaboration in Fig.  2, showing overall a good agreement. We furthermore display our result where the luminosity has been extrapolated to the full LEP luminosity.

Discovery potential at the ILC for a light Higgs boson
After having validated our method with the existing results from the analyses at LEP data, we are now in a position to apply this method in order to derive the expected limits for the ILC at √ s = 250 GeV for the two types of analyses making use of the φ → bb decay information and employing the recoil method. For our ILC analysis we use the two scaling factors c bb and c recoil that we determined for LEP 2 as given in Table 2. Concerning the ILC we assume a beam polarisation of P e − = −80% for the electron beam and P e + = +30% for the positron beam, corresponding to the baseline design [25]. At the ILC with √ s = 250 GeV a total luminosity of L = 2000 fb −1 is expected to be collected within 15 years [26]. For our study, however, we exploit only the polarisation configuration (−80%, +30%) and assume the rather modest luminosity of 500 fb −1 . Fig. 3 Projections for the expected limits on S 95 at the ILC [19] with √ s = 250 GeV, L = 500 fb −1 and polarised beams (P e − = −80% and P e + = +30%) for the analysis using the decay φ → bb and for the recoil method analysis, based only on the recoil from the Z boson that is reconstructed from the Z → μ + μ − decay. The impact of the Higgs signal at 125 GeV is not shown. The ILC projections are compared with the observed limit at LEP from the φ → bb searches [1]. The solid horizontal lines denote the indirect sensitivities at the HL-LHC and the ILC with √ s = 250 GeV from the projection for the accuracies on the couplings of the observed Higgs boson at 125 GeV [27] (see text) In Fig. 3 we show our results for the S 95 projections at the ILC for the analysis making use of the decay φ → bb decay and and for the recoil method, where only the information from Z → μ + μ − has been used [19] (the expected limits do not take into account the impact of the Higgs signal at 125 GeV). These results are compared with the observed limit from the φ → bb search at LEP [1]. As one can see from up to higher-order corrections to the sum rule. The indirect sensitivity achievable at the HL-LHC that is indicated in Fig. 3 is based on the projection for the 2σ accuracy of the g h(125)Z Z coupling from Table 4 of Ref. [27], where an integrated luminosity of 3000 fb −1 was assumed. While the indirect sensitivity at the HL-LHC exceeds the observed limit from LEP for m φ > ∼ 85 GeV, Fig. 3 clearly demonstrates that the sensitivity of the direct search at the ILC will much surpass the indirect sensitivity at the HL-LHC even for the ultimate accuracy reachable at the HL-LHC.
The indirect sensitivity of the ILC with √ s = 250 GeV from the rate measurements of the state at 125 GeV, where as for the HL-LHC we have used the projections listed in Table 4 of Ref. [27], is indicated in Fig. 3 as the area above the solid horizontal black line that is labelled as "ILC". Fig. 3 shows that this indirect sensitivity at the ILC would be similar to the direct reach via the recoil method, which could provide important complementary information for determining the nature of a possible excess in the direct searches for additional light Higgs bosons at the ILC. It should be noted that the sensitivity of the direct ILC searches for additional Higgs bosons making use of the φ → bb decay mode will significantly improve even on the indirect sensitivity of the ILC based on its high-precision measurement of the coupling of h(125) to Z bosons.

Conclusions
In this paper we have pointed out that the ILC at √ s = 250 GeV has a large physics potential in the direct search for additional light Higgs bosons. Such a light Higgs boson would be expected to have a heavily suppressed coupling to gauge bosons as compared to a SM-like Higgs boson at the same mass, which could therefore be below the existing limits from LEP. Via those direct searches the ILC would probe favoured parameter regions of various extensions of the SM, which have recently received considerable attention also in view of the excess over the background expectation that has been reported for the φ → γ γ searches at CMS in the vicinity of a long-standing excess in the LEP Higgs searches.
We have performed a generator-level study for searches via the Higgs-strahlung process at a future e + e − collider for both the analysis type using the φ → bb decay and for the decay-mode independent search via the recoil method. In a first step we have validated our approach with the existing results for the search in the φ → bb channel, from the LEP combination, and for the search utilising the recoil method, from the OPAL collaboration. We determined normalisation factors for the effective luminosities that we employ to approximately account for signal efficiencies and detector effects. We demonstrated that with this input our method reproduces the expected limits from LEP and OPAL for the two types of analyses very well.
After this validation we applied our method for deriving the expected limits for the ILC at √ s = 250 GeV for both types of analyses. We used the normalisation factors c bb and c recoil that we determined for LEP 2 but assumed as ILC conditions a centre-of-mass energy of 250 GeV, beam polarisation of P e − = −80% for the electron beam and P e + = +30% for the positron beam and used as a very conservative approach the rather modest luminosity of 500 fb −1 . Our results show that the ILC at √ s = 250 GeV will improve the LEP limits in the sensitivity to a light Higgs boson with reduced couplings to gauge bosons for the most interesting mass range between about 60 GeV and 100 GeV by more than an order of magnitude. This sensitivity of the direct search for additional light Higgs bosons at the ILC will go much beyond the indirect sensitivity of the HL-LHC from the rate measurements of the detected state at 125 GeV even for the projected ultimate accuracy reachable at the HL-LHC with 3000 fb −1 . It is interesting to note in this context that the indirect sensitivity of the ILC with √ s = 250 GeV from the rate measurements of the state at 125 GeV is similar to the direct search reach of the ILC via the recoil method, while the direct searches for additional light Higgs bosons at the ILC utilising the φ → bb decay will significantly improve even on the ultimate indirect ILC sensitivity.
The physics potential for the direct searches at the ILC with √ s = 250 GeV discussed in this paper complements and significantly enhances the ILC physics progmamme for precision measurements in the Higgs and the electroweak sector.