The particle discovered in the Higgs boson searches at the LHC with a mass of about 125 GeV can be identified with one neutral Higgs boson in a variety of Beyond the Standard Model (BSM) theories with an extended Higgs sector. Limits on the couplings between additional Higgs fields to the electroweak gauge-bosons in such theories can be obtained by model-independent Higgs searches at lepton colliders. We present an extrapolation of the limits obtained at LEP for a future lepton collider and can show that the ILC with polarized beams and a total luminosity of fb is very sensitive to such reduced Higgs-gauge-boson couplings up to about 4%. We apply the extrapolated limits on BSM models with an extended Higgs-sector.
Sensitivity of the ILC to light Higgs masses
DESY, Notkestraße 85, D–22607 Hamburg, Germany
II. Institut für Theoretische Physik, Universität Hamburg,
Luruper Chaussee 149, 22761 Hamburg, Germany
Talk presented at the International Workshop on Future Linear Colliders (LCWS2017), Strasbourg, France, 23-27 October 2017. C17-10-23.2.
With the discovery of a Higgs boson in the year 2012, the Standard Model (SM) seems to be completed. Within the current experimental accuracy, the observed properties of the Higgs boson coincide very well with the predictions for a Higgs boson in the SM. Additional Higgs bosons appear in a large class of new physics models as, for instance, in a singlet extended SM, two-Higgs-doublet model (2HDM) and in supersymmetric models (MSSM, NMSSM, etc.). These new Higgs fields can mix, but their couplings to gauge bosons have to fulfill a sum rule:
In physics models Beyond the Standard Model (BSM) the Higgs boson couplings get further contributions, however, only within the per cent level. Therefore, measurements with high precision will be mandatory to be sensitive to deviations from the SM. The measurements of all Higgs bosons to gauge bosons will therefore be an important consistency test for the closure of the SM or for opening the window to new physics.
As already mentioned, in many BSM models additional Higgs bosons appear. In the Minimal Supersymmetric extension of the Standard Model (MSSM), for instance, a Higgs boson at 125 GeV would uniquely point to a SM-like Higgs boson at this mass without any further lighter Higgs masses. However, already in the Next-to-Minimal Supersymmetric Model (NMSSM) the additional Higgs singlet relieves the bounds and offers light Higgs bosons below 125 GeV. Since such bosons necessarily have a reduced couplings to gauge bosons, very light Higgs bosons —that escaped LEP and Tevatron bounds— are still not yet excluded and provide a very interesting framework for new physics.
In the current paper we study the sensitivity of the ILC in its first stage of GeV cms energy to such light Higgs bosons far below the LEP limit. Due to the clean environment, the low beamstahlung, the precise knowledge on the beam energy and the availability of polarized beams, the ILC is perfectly suited to be sensitive to such light additional Higgs bosons.
The two main production processes for Higgs bosons are Higgs-strahlung (), dominant at lower masses and energies, and via -fusion () dominant at higher Higgs masses and cms energies. In the current study, we concentrate therefore on the Higgs-strahlung process and use two approaches: we either study the Higgs decay process and use the -decay into for validation (perfomed at LEP) or use the recoil method and expoit only the -decays, here , (perfomed at the OPAL experience).
As criteria whether a light Higgs boson is accessible, the statistical method has been applied for the two approaches, the details are explained in Sect. 2. In Sect. 3 the method is validated with the results obtained at LEP and in Sect. 4 the discovery potential of the ILC for such light Higgs bosons is discussed. Conclusions are given in Sec. 5.
2 Details on the applied methods
We consider the ’Higgs-strahlung’ process
with the -boson and the scalar decaying into fermions,
In order to study the sensitivity to lighter Higgs masses, we evaluate a hypothetical signal process, where the SM-Higgs in Eq. 3 is replaced by a new scalar . The new scalar is assumed to have the same couplings as the SM-Higgs scalar, but to have a different mass .
We use the statistical method called , described in refs. [1, 2], here used in a simplified approach, see also , where event samples for the two hypotheses ’all events are generated by the background only’ or ’all events are generated by the background plus a hypothetical signal’. The quantity
gives an upper limit on a cross-section, that is compatible with the ’background only’ hypothesis at the 95% confidence level, normalized on a reference cross section . As reference process we regard in Eq. 3 the scalar to be the Higgs boson of the SM. This quantity is equal to the maximal allowed signal rate normalized on the total signal rate . Since it is assumed that has SM-coupling but only a different mass, Eq. 4 can be interpreted as the ratio between the squared couplings of the scalars () with the -boson, (), respectively:
The event samples for
are generated with the Monte Carlo generator Whizard-2.4.1 and contain the signal processes of Eq. 3.
For the given process gauge invariance forbids to generate signal and background events separately. The obtained events are ordered by either the invariant or recoil mass of the quark- or lepton-system for this analysis. The obtained signal rate in the -th bin is obtained by the difference of the event rates for the background plus signal () and the background only hypotheses (), and the signal rate is obtained by . Negative rates for signal events are cut, .
The obtained events are identified by either the invariant or recoil mass of the quark- or lepton-system for this analysis.
We assume that the number of potential signal events in each of the bins are distributed according to a Poisson distribution with the expected values for the ’background only’ and 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 have to be determined experimentally. For this work we consider only simulated events, and thus . The logarithm of the test-statistic reveals 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. , and scale the luminosity afterwards with a factor to the desired value. The weights are not affected by the scaling procedure.
3 Validation of the methods with results from LEP I and LEP II
In order to test the method we apply it in a first step on the obtained data from LEP and try to recover the published results and the interpretation for new physics [1, 4]. Furthermore, in a second step, we reproduce the results of the OPAL collaboration that performed also Higgs analysis with the recoil method .
We simulated the process with the SM implementation of Whizard 2.4.1 [6, 7]. This process contains the decays of the scalar and -boson in the Higgs-strahlung process, where we consider only the dominant process with the Higgs decaying into -quarks and the -boson decaying into the -leptons. Our ’background only’ and ’background plus signal’ hypotheses are generated by adapting the mass of the implemented SM-Higgs field, correspondingly, i.e. concerning the ’background only’ hypothesis we chose the mass of the implemented scalar beyond the kinematic accessibility of the LEP experiments.
For LEP I we consider the integrated luminosity at the center-of-mass energy of 91.2 GeV. For LEP II we simulate the experiment with the integrated luminosity per year with the highest achieved energy in each period, with exception of the last year, where we assume the lowest energy of 206 GeV . The considered luminosity per experiment is given in tab. 1.
For the simulation, however we generated a luminosity that was 400 times the LEP-luminosity in order to obtain finite weight factors (Eq. 8) and to reduce the statistical error.
We do not consider hadronization effects of the -quarks or detector effects. This simplification leads to an over-optimistic estimate for the signal efficiency and thus for the observed signal rates. This effect is countered, however, by the restraint to only one decay mode for the scalar and -bosons.
Since the signal rates are not explicitly given for the LEP experiments , we scale the event rates of our results by a factor , which yields an effective luminosity of
where is the actual luminosity derived from the values given in tab. 1. The scaling factors depend on the parameters of the experiment and the analysis. We determined the factors via comparison with the official numbers from LEP for . Such a scaling factor can only approximate the effects of an exact analysis. In case that the decay modes and are considered —we call this the ’traditional’ method at LEP— it is expected to obtain the luminosity (for an ideal detector). Assuming a median signal efficiency of 50% for the signal rates, the coefficients should therefore be not larger than . Applying the recoil method, used by the OPAL experiment, only the decay mode has been used and explicit efficiency have been published, therefore the corresponding coefficients are smaller. The coefficients are given in tab. 2.
3.1 Comparison with LEP data using the / decay mode – ’LEP Traditional’
In order to compare with LEP data, event samples for the process are generated. However, this process depends on the width via . Changing the mass would therefore also induce to adjust the width. However, we obtain the best consistency results when assuming a constant width of the scalar
The events for ’background only’, , and ’background plus signal’, , are taken from the decay channel. The -decay channel is used only for validation that one observes both a signal from a scalar and a -boson. A challenging region for this analysis is if the scalar is close to , since a large number of background events close to the -boson mass , weakens the result for ; 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 the and fit these points with a third order polynomial to obtain the number of events in the bins around as the value of the fit function at the central mass of each bin.
Our result concerning the reconstruction of the LEP expectation for are given in Fig.1 including scaling factors , corresponding to Table 2, that accounts for detector effects and the restriction to the -channel only. As can be seen from Fig.1, the expected limits are very well reconstructed with our method.
3.2 Comparison with OPAL data using the recoil method – ’LEP Recoil’
The OPAL experiment also used the recoil method for analyzing the data, i.e. exploiting the recoil of on the -boson and analyzing the Z-decay only. This method has the great advantage of being completely independent of the decay channels.
The OPAL analysis used the decay modes , . We restrict our analysis to the decay mode only, but extrapolated in a second step the luminosity taken by OPAL to the full LEP luminosity, see Table 1. We include only bins close to (interval [84 GeV, 98 GeV]), so that the weighted mean of the central masses of the bins is in a small interval around
with the number of events in the -th of bins with the central mass .
The events from the -lepton pairs are ordered by their total energy into bins and for each bin we calculated the respective recoil mass ,
to obtain the event rates for background , , and signal,
A mass-dependent width is taken into account,
Since the OPAL experiment published the used signal efficiency in ref. , we estimated a mean signal efficiency of for our approach and obtained the effective luminosity as
where the obtained values for the scaling factor are given in Tab. 2. Practically , since the used efficiency is given. It turned out that neglecting hadronization as well as detector effects practically compensated our disregard of the -channel.
The reconstructed results for for OPAL and our comparison are shown in Fig.2, as well as the distribution to the full LEP luminosities. A you can see, also the reconstruction of the recoil results seem to be very promising.
4 Discovery potential at the ILC for light Higgs
After having validated the method with the LEP data, we turn our analysis now into expectation for the ILC at GeV. We assume a beam polarization of for the electron beam and for the positron beam, corresponding to the baseline design . We assume a total luminosity of fb expected to be collected within 15 years .
We compare the projections for the ILC both for the case, where the -decays have been taken into account, corresponding to the ’LEP traditional’ method, and the recoil method, where only has been used. Of course, a mass-dependent width for has been calculated and taken into account. We compare these new results with the derived curve for LEP. As one can see in Fig. 3, we gain about a factor 10 compared to LEP with the ILC, already with only a cms of GeV.
From Fig. 3 we derive the following limits for :
With respect to Eq. 5, these limits predict a sensitivity of the ILC to deviations from Standard-Model gauge boson-Higgs couplings up to
These limits for a reduced coupling in the light scalar sector can be easily accommodated within a large class of new physics models, 2HDM, MSSM, NMSSM. Taking the NMSSM, for instance, scanning through the plane and performing a global fit and comparing with the Standard Model as done in , see Fig. 4, allows perfectly well these range of limits. As can be seen from Fig. 4, where the fit value of the reduced squared coupling of a singlet-like Higgs state to the Z-boson is given in dependence of the light singlet Higgs mass, the ILC at GeV would be sensitive to the light singlet mass region down to about 60 GeV and a reduced Higgs-gauge-boson coupling of up to only a few percent.
Such a sensitivity greatly enhances the physics potential of the ILC at GeV. If, for instance, a signal around GeV, see Fig. 5 (left) (corresponding to Fig. 7 in ), would become reality, then one could provide strict bounds for the new Scalar-gauge-boson couplings in different BSM models. One should remember that also at LEP, there was a slight excess for a light scalar GeV, see Fig. 5 (right) (corresponding to Fig. 10 in ).
At the ILC, one therefore would expect at this mass a limit of , cf. 3, which would allow that a gauge-boson-scalar-coupling of such a new scalar at 95 GeV is only about 3.9% of the SM-coupling. Such reduced couplings in the light scalar sector can be accommodated with a large class of new physics models, as, for instance, in the NMSSM, see discussion above.
In the current study we could show that an ILC at GeV has a large, still unexplored, physics potential for the discovery of light, non-SM-like Higgs masses. The sum rule for Higgs-gauge-boson-couplings is very restrictive w.r.t. new physics models and allows only small admixtures. We use the method to estimate the sensitivity to non SM-like couplings and cross sections. We validate the method with LEP data both the traditional method, i.e. using Higgs decays and Z-decays, as well as the recoil method where only the decays of the Z-boson are exploited within the Higgs-strahlungs process. With an ILC at GeV and an integrated luminosity of about fb and with a polarization configuration of , one has an impressive higher sensitivity up to one order of magnitude in this channel, comparing LEP/LHC data for masses in the range GeV. Furthermore one can probe the hypothesis of an current excess at GeV with regard to the new scalar-gauge-boson and we show that the ILC could be sensitive to a respective couplings of only . That’s an impressive improvement compared to LHC results. The ILC already at its first energy stage of GeV will therefore provide essential steps towards new physics, is highly sensitive to such non-SM-like Higgs couplings opening a window to a large class of new phyics models as 2HDM, MSSM, NMSSM.
The authors acknowledge support by the DFG through the SFB 676 “Particles, Strings and the Early Universe”. This research was supported in part by the European Commission through the “HiggsTools” Initial Training Network PITN-GA-2012-316704.
- email: email@example.com
- email: firstname.lastname@example.org
- email: email@example.com
- OPAL, DELPHI, LEP Working Group for Higgs boson searches, ALEPH, L3, R. Barate et al., Phys. Lett. B565 (2003) 61, hep-ex/0306033.
- P. Bock, JHEP 01 (2007) 080, hep-ex/0405072.
- P. Bechtle, SUSY Higgs Boson Searches at LEP and SUSY Parameter Measurements at TESLA, PhD thesis, Hamburg U., 2004.
- DELPHI, OPAL, ALEPH, LEP Working Group for Higgs Boson Searches, L3, S. Schael et al., Eur. Phys. J. C47 (2006) 547, hep-ex/0602042.
- OPAL, G. Abbiendi et al., Eur. Phys. J. C27 (2003) 311, hep-ex/0206022.
- W. Kilian, T. Ohl and J. Reuter, Eur. Phys. J. C71 (2011) 1742, 0708.4233.
- M. Moretti, T. Ohl and J. Reuter, (2001) 1981, hep-ph/0102195.
- R. Assmann, M. Lamont and S. Myers, Nucl. Phys. Proc. Suppl. 109B (2002) 17, [,17(2002)].
- C. Adolphsen et al., (2013), 1306.6328.
- K. Fujii et al., (2017), 1710.07621.
- M. Klute et al., EPL 101 (2013) 51001, 1301.1322.
- D. Lopez-Val, T. Plehn and M. Rauch, JHEP 10 (2013) 134, 1308.1979.
- F. Domingo and G. Weiglein, JHEP 04 (2016) 095, 1509.07283.
- CMS, C. Collaboration, (2017).