Sensitivity to invisible scalar decays at CLIC

Abstract

We studied the possibility of constraining production of new scalar particles at CLIC running at 380 GeV and 1.5 TeV, assuming the associated production of Higgs-like neutral scalar with $\mathrm{Z}{}{}$ boson and its invisible decays. The analysis is based on the Whizard event generation and fast simulation of the CLIC detector response with Delphes. We considered ${\mathrm{e}{}{}}^{+} {\mathrm{e}{}{}}^{-}$ background processes but also relevant $\upgamma {}{} \upgamma {}{}$ and $\upgamma {}{} \mathrm{e}{}{} ^{\pm }$ interactions. The approach consisting of a two-step analysis was used to optimise separation between signal and background processes. First, a set of preselection cuts was applied; then, multivariate analysis methods were employed to optimise the significance of observations. We first estimated the expected limits on the invisible decays of the 125 GeV Higgs boson, which were then extended to the cross section limits for production of an additional neutral scalar, assuming its invisible decays, as a function of its mass. Extracted model-independent branching ratio and cross section limits were then interpreted in the framework of the Higgs-portal models to set limits on the mixing angle between the SM-like Higgs boson and the new scalar of the “dark sector”.

1 Introduction

All available experimental results seem to confirm that the new particle discovered in 2012 by the ATLAS and CMS experiments at LHC [1, 2] is the last missing constituent of the Standard Model (SM), the Higgs boson. The Standard Model predicts that the Higgs boson with a mass of about 125 GeV should decay predominantly to $\mathrm{b}{}{} \overline{\hbox {b}}{}{} {}{}$ (about 58% of all decays), but also to $\mathrm{W}{}{} \mathrm{W}{}{} ^*$ (21%), $\tau ^{+} \tau ^{-}$ (6,3%) or $\mathrm{Z}{}{} \mathrm{Z}{}{} ^*$ (2,6%) [3], the latter channel resulting also in about 0.1% of ‘invisible’ decays (with both $\mathrm{Z}{}{}$ bosons decaying to neutrinos). Some extensions of the Standard Model predict additional channels for invisible Higgs boson decays—into new, unobservable particles. These particles could contribute to the dark matter (DM) density of the Universe. As of today, the best direct limits on invisible Higgs boson decays come from ATLAS and CMS experiments at the LHC—at 95% C.L. the branching fraction is less than 13% [4] and 19% [5], respectively.

Experimental constraints on the invisible Higgs boson decays can be set either directly, by searching for such decays in channels where Higgs boson production can be tagged independently from the decay mode (e.g. via vector boson fusion in $\mathrm{p}{}{} \mathrm{p}{}{}$ collisions or production together with a $\mathrm{Z}{}{}$ boson in ${\mathrm{e}{}{}}^{+} {\mathrm{e}{}{}}^{-}$ collisions), or indirectly, based on the global fit to all production and decay channel measurements, assuming that the total width of the Higgs boson can also be directly measured. Prospects for constraining properties of the 125 GeV Higgs boson at high luminosity LHC (HL-LHC) and different future collider projects were recently summarised in [6]. For all considered electron–positron Higgs factories, the Future Circular Collider (FCC-ee) [7], the Circular Electron Positron Collider (CEPC) [8], the International Linear Collider (ILC) [9] and the Compact Linear Collider (CLIC) [10], expected indirect constraints on the “untagged” Higgs boson decays are on the level of 1–2%. Stronger constraints are expected from direct searches at ${\mathrm{e}{}{}}^{+} {\mathrm{e}{}{}}^{-}$ Higgs factories, where the Higgs-strahlung process, ${\mathrm{e}{}{}}^{+} {\mathrm{e}{}{}}^{-} \rightarrow \mathrm{Z}{}{} \mathrm{H}{}{}$, allows for decay-mode independent tagging of Higgs boson production events. Both leptonic and hadronic decay modes of the $\mathrm{Z}{}{}$ boson can be considered. While the leptonic final state allows for more precise reconstruction of the event and better background suppression, sensitivity to invisible Higgs boson decays is dominated by hadronic $\mathrm{Z}{}{}$ boson decay channel providing an order of magnitude higher statistics [11].

Invisible decays of the SM-like Higgs boson are considered as one of the possible signatures for existence of DM particles and physics Beyond the Standard Model (BSM) [12,13,14]. In the so-called Higgs-portal models, existence of additional fundamental scalars of the “hidden sector” is assumed [15,16,17]. These new particles could mix with the SM Higgs, and thus, open new decay channels of the SM-like 125 GeV state into DM particles. The new scalars, if they are relatively light, could also be produced in ${\mathrm{e}{}{}}^{+} {\mathrm{e}{}{}}^{-}$ collisions in the process corresponding to the Higgs-strahlung process in the SM: ${\mathrm{e}{}{}}^{+} {\mathrm{e}{}{}}^{-} \rightarrow \mathrm{Z}{}{} \mathrm{H}{}{} '$. We extend our search for invisible decays of the SM-like Higgs boson to the search for production and invisible decays of a scalar particle, $\mathrm{H}{}{} '$, with arbitrary mass. We then interpret our results in terms of the limits on the scalar sector mixing angle in the Higgs-portal scenarios. For illustration, we employed the vector-fermion dark matter model (VFDM) [18, 19], a simple extension of the Standard Model with one extra scalar, two Majorana fermions and one gauge boson. The choice was primarily based on our previous experience with this model [20, 21].

Considered in this paper is sensitivity of CLIC running at 380 GeV and 1.5 TeV to invisible decays of the scalar particle produced in association with the $\mathrm{Z}{}{}$ boson. Only hadronic $\mathrm{Z}{}{}$ boson decays are considered, following the approach described in [22, 23]. The expected sensitivity of CLIC running at 380 GeV to invisible decays of the 125 GeV Higgs boson is weaker than that expected for other Higgs factories running at energies of 240–250 GeV, close to the maximum of the Higgs-strahlung cross section [6]. However, larger collision energy allows to search for production of new scalars in a wider mass range. For CLIC running at 1.5 TeV, scalar masses of up to about 1 TeV can be probed. Production of new scalars in the Higgs-strahlung process has not been considered before for CLIC. Similar study, based on the leptonic $\mathrm{Z}{}{}$ boson decays, has been recently presented for ILC [24].

2 Event generation and fast simulation of the detector response

Results presented in this paper are based on a realistic simulation of events, including fast simulation of the CLICdet [25] detector response. Signal and background event samples were generated using Whizard 2.7.0 [26, 27], taking into account the beam energy profile expected for CLIC running at 380 GeV and 1.5 TeV parametrised with the CIRCE2 package [28]. As the signal, we considered the Higgs-strahlung process: Higgs boson production (with decay into an invisible final state) together with a $\mathrm{Z}{}{}$ boson decaying into a quark–antiquark pair. Samples of events with production of the new, hypothetical scalar particle $\mathrm{H}{}{} '$ with the emission of two quarks (${\mathrm{e}{}{}}^{+} {\mathrm{e}{}{}}^{-} \rightarrow \mathrm{H}{}{} '\mathrm{Z}{}{} \rightarrow inv+\mathrm{q}{}{} \mathrm{q}{}{}$) were also generated using the Standard Model implementation in Whizard (SM_CKM model), but with the modified Higgs boson mass and width (according to SM predictions for given mass).¹ Masses of the new scalar in the range 120–280 GeV (for the first stage of CLIC) and 150–1200 GeV (for the second stage) were considered.

For the background, we studied processes both with and without Higgs boson production. We also took into account the possible background contribution from $\upgamma {}{} \upgamma {}{}$ and $e^{\pm }\gamma$ interactions, where both beamstrahlung photons ($\upgamma {}{} ^{_\text {BS}}$) and photon radiation by the incoming electrons, as described by the effective photon approximation² ($\upgamma {}{} ^{_\text {EPA}}$), were taken into account. The baseline design of CLIC [10] includes polarisation of the electron beam only. For the first running stage, at 380 GeV, the same integrated luminosity of data is expected to be collected with negative and positive electron beam polarisation. We could therefore analyse this data as a single, unpolarised, dataset with an integrated luminosity³ of 1000 fb$^{-1}$ [29]. We also consider the alternative running scenario presented recently [30], where up to 4000 fb$^{-1}$ of data can be collected at the initial CLIC stage. At 1.5 TeV we consider the two electron beam polarisations separately, assuming 2000 fb$^{-1}$ to be collected with -80% polarisation and 500 fb$^{-1}$ with +80% polarisation [29]. Cross sections for processes taken into account in the presented study⁴ calculated by Whizard and numbers of generated events are shown in Tables 1 (for 380 GeV running) and 2 (for 1.5 TeV). The cross sections include standard cuts at generator level, as applied in [23]. Only for the $\upgamma {}{} ^{_\text {BS}}\upgamma {}{} ^{_\text {BS}}\rightarrow \mathrm{q}{}{} \mathrm{q}{}{}$ sample, the cut on the invariant mass of the produced quark pair was increased to 50 GeV, to avoid large contribution from soft events. Similarly, four-momentum transfer value between the incoming and outgoing electron (or positron) was required to be greater than 100 GeV for ${\mathrm{e}{}{}}^{+} {\mathrm{e}{}{}}^{-} \rightarrow {\mathrm{e}{}{}}^{+} {\mathrm{e}{}{}}^{-}$ events. Numbers of Monte Carlo events generated for some of the background channels are below the event statistics expected in the actual experiment. However, it was verified that the resulting statistical uncertainties of the presented results are on percent level. The statistical cross section uncertainties resulting from the Whizard integration are of the order of one permille or below and have been neglected in the analysis.

Table 1 Cross sections, $\sigma$, and numbers of generated events, N$_{^\text {GEN}}$, for each final state considered at 380 GeV; the reference cross section for $\mathrm{H}{}{} _\text {inv}+ \mathrm{q}{}{} \mathrm{q}{}{}$ channel is the SM Higgs boson production cross section for $\mathrm{H}{}{} _\text {SM}+ \mathrm{q}{}{} \mathrm{q}{}{}$ final state

Table 2 Cross sections, $\sigma ^{neg}$ and $\sigma ^{pos}$, for $-80\%$ and +80% electron beam polarisation, respectively, at 1.5 TeV CLIC, and numbers of generated events, N$_{^\text {GEN}}$, for each considered final state; the reference cross sections for $\mathrm{H}{}{} _\text {inv}+ \mathrm{q}{}{} \mathrm{q}{}{}$ channel are the SM Higgs boson production cross sections for $\mathrm{H}{}{} _\text {SM}+ \mathrm{q}{}{} \mathrm{q}{}{}$ final state

To simulate detector response, the fast simulation framework Delphes was used [31]. Control cards prepared for the new detector model CLICdet [32] were modified to make Higgs particles ‘invisible’ in the simulation (ignored when generating detector response), so that the invisible Higgs boson decays can be modelled by defining the Higgs boson as stable in Whizard and Pythia. An event reconstruction begins with searching for isolated electrons, muons and photons (assuming reconstruction efficiency resulting from the full detector simulation). For CLICdet, Delphes identifies isolated electrons and photons with energy of at least 2 GeV and muons of 3 GeV. Jet clustering was carried out with the VLC algorithm [33] run in exclusive mode for reconstruction of two hadronic jets with minimal transverse jet momentum of 20 GeV. The algorithm was run with parameters $\beta$ = 1 and $\gamma$ = 1, and the optimal value for the parameter describing the jet cone size was found to be $R=1.5$ for 380 GeV and $R=0.7$ for 1.5 TeV CLIC running. While the algorithm is configured to reconstruct two hadronic jets, VLC distance values at which an event transitions from four to three jets and from three to two jets are stored as parameters y$_{34}$ and y$_{23}$, respectively. We further require that each jet includes at least two charged particles.

Because of the 0.5 ns bunch spacing in the CLIC beams, the pile-up of beam-induced backgrounds can affect the event reconstruction. Realistic levels of pile-up from the most important beam-induced background, the $\upgamma {}{} \upgamma {}{} \rightarrow$hadrons process, were included in full simulation studies to verify the impact on event reconstruction results [23]. A selection based on the time stamp and transverse momentum of the reconstructed particles can be used to reduce this background significantly. Further suppression results from use of the VLC algorithm which was designed to be more resilient to the impact of backgrounds [33]. Effects of the beam-induced backgrounds turned out to be negligible for CLIC running at 380 GeV. To take into account their contribution at 1.5 TeV, additional energy smearing is applied for jets reconstructed at 1.5 TeV: for central jets ($|\cos \theta |<0.64$), the pile-up contribution is expected to result in an additional 1% jet energy smearing, while for more forward jets ($|\cos \theta |\ge 0.64$) 5% smearing is assumed [32]. Jet momentum is scaled by the same factor as the jet energy⁵ and possible impact of pile-up on the reconstruction of the jet direction is neglected. Validity of the proposed approach for the reconstruction of hadronic vector boson decays has been recently demonstrated in a similar study [34]. Impact of the jet energy smearing on the expected cross section limits presented in this work for CLIC running at 1.5 TeV was found to be small, at the level of 3–5%.

3 Signal event selection

3.1 Preselection of events

A main purpose of the preselection is to remove all background events which are not consistent with the expected signature of the signal process. For the process ${\mathrm{e}{}{}}^{+} {\mathrm{e}{}{}}^{-} \rightarrow \mathrm{H}{}{} \mathrm{Z}{}{} \rightarrow inv+\mathrm{q}{}{} \mathrm{q}{}{}$, the two reconstructed jets are expected to have an invariant mass consistent with the mass of the Z boson. In the initial preselection step, all events which were not consistent with this signature were rejected. In particular, events with isolated leptons (electrons or muons) or isolated energetic photons (with energy greater than 5 GeV) were excluded. For a significant fraction of events, the difference between the energy sum of the reconstructed jets and the energy sum of all identified particles in the event was sizable, indicating an incomplete event reconstruction (e.g. additional jet with transverse momentum below the required threshold of 20 GeV or deposits removed by the VLC algorithm run in exclusive mode due to the small beam distance). To avoid such events, we required this difference to be less than 10 GeV.

In the next steps, quantities describing event topology were considered. First, we analysed the distributions of parameters y$_{23}$ and y$_{34}$ describing the results of jet clustering with VLC algorithm. While the algorithm was forced to reconstruct two jets in each event, these distributions allowed us to distinguish actual two-jet events from events with a larger number of underlying jets in the final state. The distributions of -log$_{10}$y$_{23}$ and -log$_{10}$y$_{34}$ are shown in Fig. 1a and b, respectively. The double-peak structure, clearly visible in the background sample of SM Higgs boson decays ($\mathrm{H}{}{} _\text {SM}$), corresponds to pure two- and four-jet events. The distribution for the other background channels is more uniform. For the signal sample, with two hadronic jets, we expect values of y$_{23}$ and y$_{34}$ to be relatively small. We selected events for which y$_{23}<0.01$ (-log$_{10}$y$_{23}>$ 2.0) and y$_{34}<0.001$ (-log$_{10}$y$_{34}>$3.0).

Fig. 1

Distributions of variables used for event preselection: jet clustering parameters, y$_{23}$ (a) and y$_{34}$ (b), dijet invariant mass, m$_{jj}$ (c) and the cosine of the dijet emission angle, $\cos (\theta )$ (d). Compared are distributions for different event samples: total background (thin solid black line), ${\mathrm{e}{}{}}^{+} {\mathrm{e}{}{}}^{-}$ background without Higgs boson production (long-dashed pink line), background of SM Higgs boson production and decays (short-dashed red line), photon interactions (dash-dotted blue line) and signal (thick solid green line), assuming integrated ${\mathrm{e}{}{}}^{+} {\mathrm{e}{}{}}^{-}$ luminosity of 1000 fb$^{-1}$ collected at 380 GeV. For easier shape comparison, the signal sample is normalised to BR$(\mathrm{H}{}{} \rightarrow inv)=100\%$. Black vertical lines indicate preselection cuts, see text for details

The next quantity considered for the preselection of signal events was the invariant mass of the two-jet final state—m$_{jj}$. It should correspond to the mass of the $\mathrm{Z}{}{}$ boson, so only events for which this value was in the range of 80–100 GeV were selected for further analysis. The distribution of the invariant masses for different event samples is shown in Fig. 1c. For ${\mathrm{e}{}{}}^{+} {\mathrm{e}{}{}}^{-}$ SM background events (SM bg), dijet mass distribution is also peaked at the mass of the $\mathrm{Z}{}{}$ boson, so only the tails of the distribution can be rejected. However, photon-induced backgrounds ($\upgamma {}{} \upgamma {}{}$ and $\mathrm{e}^{\pm }\gamma$) are dominated by hadronic decays of the $\mathrm{W}{}{}$ boson and significant suppression is possible with the dijet mass cut. The peak visible in the channel of SM Higgs boson decays ($\mathrm{H}{}{} _\text {SM}$) around 120 GeV corresponds to the process ${\mathrm{e}{}{}}^{+} {\mathrm{e}{}{}}^{-} \rightarrow \mathrm{H}{}{} \mathrm{Z}{}{} \rightarrow \mathrm{q}{}{} \mathrm{q}{}{} \upnu {}{} \upnu {}{}$ which, considering only its topology (two jets and missing energy), is indistinguishable from the signal.⁶ We also studied the distribution of the dijet emission angle, $\theta$, defined as the angle between the beam axis and a sum of the jet four-momenta (the emission angle of the Z boson for the signal events). For the majority of background events, small emission angles are reconstructed, while for signal events the distribution is almost flat (angles close to 90$^\circ$ are slightly preferred). Therefore, events for which $|\cos (\theta )|$ was greater than 0.8 were excluded from the analysis. The distribution of cosine of the angle $\theta$ is shown in Fig. 1d.

The results of the preselection are presented in Tables 3 and 4. Shown in Fig. 2 is the distribution expected for CLIC running at 380 GeV of the so-called recoil mass, the invariant mass of the Higgs boson produced together with the Z boson, after preselection cuts, reconstructed from the energy–momentum conservation (missing mass). For the background sample, the distribution has two maxima: at around 300 GeV, which is the maximum recoil mass allowed (as we require two jets to have an invariant mass of at least 80 GeV), and at around 90 GeV, which is mainly due to invisible Z boson decays. For signal events, normalised in Fig. 2 to BR$(\mathrm{H}{}{} \rightarrow inv)=1\%$, the expected recoil mass distribution is consistent with the SM Higgs boson mass of 125 GeV. The slight shift of the maxima and longer tail in the reconstructed recoil mass distribution towards higher mass values is most likely due to the influence of the beam energy spectra, which is not accounted for in the recoil mass reconstruction.

Table 3 Efficiency and expected event numbers after preselection, N$_\text {pre}$, corresponding to the integrated ${\mathrm{e}{}{}}^{+} {\mathrm{e}{}{}}^{-}$ luminosity of 1000 fb$^{-1}$ collected at 380 GeV, for each considered sample or events. Results for the signal sample were calculated for SM Higgs boson production cross section and BR$(\mathrm{H}{}{} \rightarrow inv)=100\%$

Table 4 Efficiency and expected event numbers after preselection, N$_\text {pre}$, corresponding to the integrated ${\mathrm{e}{}{}}^{+} {\mathrm{e}{}{}}^{-}$ luminosity of 2000 fb$^{-1}$ (negative polarisation) and 500 fb$^{-1}$ (positive polarisation) collected at 1.5 GeV, for each considered sample of events. Results for the signal were calculated for SM Higgs boson production cross section and BR$(\mathrm{H}{}{} \rightarrow inv)=100\%$

Fig. 2

Distributions of the reconstructed invariant mass of the invisible Higgs boson decay products expected for different event samples after preselection cuts: total background (thin solid black line), ${\mathrm{e}{}{}}^{+} {\mathrm{e}{}{}}^{-}$ background without Higgs boson production (long-dashed pink line), background of SM Higgs boson production and decays (short-dashed red line), photon interactions (dash-dotted blue line) and signal (thick solid green line), assuming integrated ${\mathrm{e}{}{}}^{+} {\mathrm{e}{}{}}^{-}$ luminosity of 1000 fb$^{-1}$ collected at 380 GeV. The signal sample is normalised to BR$(\mathrm{H}{}{} \rightarrow inv)=1\%$

3.2 Final selection

The second stage of the analysis was based on multivariate analysis and machine learning. The boosted decision tree (BDT)[35] algorithm, as implemented in TMVA framework[36], was used, with 1000 trees and 5 input variables. The following parameters were selected as the BDT input variables:

• 1. E$_{j\!j}$—dijet energy,

• 2. m$_{j\!j}$—dijet invariant mass,

• 3. m$^{\text {miss}}$—reconstructed recoil mass,

• 4. p$^{\text {miss}}_{t}$—missing transverse momentum,

• 5. $\alpha _{j\!j}$—angle between the two reconstructed jets in the LAB frame.

This choice of parameters was selected as optimal, resulting in the most efficient event classification, from a large number of different parameter sets considered. Distributions of the selected variables for the signal and the background samples are shown in Fig. 3.

Fig. 3

Distributions of the input variables selected for the BDT algorithm: dijet energy, E$_{j\!j}$ (a), dijet invariant mass, m$_{j\!j}$ (b), reconstructed recoil mass, m$^{\text {miss}}$ (c), missing transverse momentum, p$^{\text {miss}}_{t}$ (d) and angle between the two reconstructed jets, $\alpha _{jj}$ (e). Expected signal (thick blue line) and background (thin red line) event distributions for CLIC running at 380 GeV are normalised to 1000 fb$^{-1}$ and BR$(\mathrm{H}{}{} \rightarrow inv)=1\%$ for signal events

Fig. 4

Distribution of BDT response for the signal (thick blue line) and the background (thin red line) events (left) and the expected significance of signal observation as a function of the BDT response cut (right), for 1000 fb$^{-1}$ collected at 380 GeV and BR$(\mathrm{H}{}{} \rightarrow inv)=1\%$ for signal events

Distributions of the BDT algorithm response for considered signal and background event samples (after preselection cuts) are shown in Fig. 4a. Most of the background events can be easily distinguished from the signal, but there is also a significant contribution of background events for which BDT response values are positive, consistent with the response expected for signal events. This indicates that it is not possible to achieve full separation between the signal and the background processes. One should note that about 0.1% of SM Higgs boson decays result in fact in the invisible final state ($\mathrm{H}{}{} \rightarrow \mathrm{Z}{}{} \mathrm{Z}{}{} ^* \rightarrow \upnu {}{} \upnu {}{} \upnu {}{} \upnu {}{}$), which is included in the background simulation.

In the final step of the analysis, we select the cut on the BDT response which gives the highest significance for the expected signal. The dependence of the signal significance at 380 GeV on the BDT algorithm response cut is shown in Fig. 4b, for the signal sample normalised to BR$(\mathrm{H}{}{} \rightarrow inv)=1\%$. The highest significance for invisible Higgs decays at 380 GeV CLIC is obtained for a BDT response cut of about 0.14, corresponding to a BDT selection efficiency for signal events of about 50% and background rejection efficiency of about 95%. The background remaining after the BDT cut is dominated by contributions from ${\mathrm{e}{}{}}^{+} {\mathrm{e}{}{}}^{-} \rightarrow \mathrm{q}{}{} \mathrm{q}{}{} \upnu {}{} \upnu {}{}$ (68% of expected background events), ${\mathrm{e}{}{}}^{+} {\mathrm{e}{}{}}^{-} \rightarrow \mathrm{q}{}{} \mathrm{q}{}{} \mathrm{l}{}{} \upnu {}{}$ (14%) and $\gamma \mathrm{e}^{\pm } \rightarrow \mathrm{q}{}{} \mathrm{q}{}{} \upnu {}{}$ (17%). While we do not discuss systematic uncertainties of the measurement here, one has to note that they can be significantly constrained based on measurements of ${\mathrm{e}{}{}}^{+} {\mathrm{e}{}{}}^{-} \rightarrow \mathrm{Z}{}{} \mathrm{Z}{}{}$ and ${\mathrm{e}{}{}}^{+} {\mathrm{e}{}{}}^{-} \rightarrow \mathrm{W}{}{} \mathrm{W}{}{}$ processes in other decay channels (data-driven approach).

The same analysis procedure was applied for signal and background samples generated for CLIC running at 1.5 TeV, separately for two considered electron beam polarisation settings. Distributions of the BDT algorithm response for considered signal and background event samples (after preselection cuts) are shown in Fig. 5.

Fig. 5

Distribution of BDT response for the signal (thick blue line) and the background (thin red line) for CLIC running at 1.5 TeV, with electron beam polarisation of $-80\%$ (left) and $+80\%$ (right)

A similar level of signal–background separation is obtained for each polarisation. The background remaining after the BDT cut is again dominated by $\mathrm{q}{}{} \mathrm{q}{}{} \upnu {}{} \upnu {}{}$ contribution (about 74%) and $\gamma \mathrm{e}^{\pm }$ processes (17%), while contribution from $\mathrm{q}{}{} \mathrm{q}{}{} \mathrm{l}{}{} \upnu {}{}$ channel decreases and is at the level of about 7%.

The analysis procedure developed to discriminate between the background of different SM processes and the signal of invisible scalar decays, described above for the case of the 125 GeV SM-like Higgs boson, was used to estimate the expected sensitivity of CLIC experiment to production and invisible decays of a new scalar state. Same preselection cuts and same set of input BDT variables were used for each scenario, and the BDT algorithm was trained separately for each considered scalar mass (each generated signal sample). With the BDT response cut optimised for signal significance, one can extract an expected limit on the cross section for the production of the new scalar $\mathrm{H}{}{} '$ as a function of its mass, assuming its invisible decays, BR$(\mathrm{H}{}{} '\rightarrow inv) = 100\%$.

For the scalar mass of 125 GeV, expected numbers of background events and efficiency of signal event selection can also be translated into a constraint on the invisible branching ratio of the SM-like Higgs boson. For the first stage of the CLIC accelerator, assuming that the measured event distributions are consistent with the predictions of the Standard Model, the expected 95% C.L. limit⁷ is:

$$\begin{aligned} \hbox {BR}(\mathrm{H}{}{} \rightarrow inv)<1.0\% \;\;\; (0.5\%) \end{aligned}$$

for the integrated luminosity of 1000 fb$^{-1}$ (4000 fb$^{-1}$).⁸ A significance above $5 \sigma$, necessary to confirm the discovery of a new decay channel (and therefore also existence of new, invisible particles), is expected for an invisible Higgs boson branching ratio above 3.0% (1.5%). Presented results seem to be consistent with previous estimates presented by CLICdp Collaboration, BR$(\mathrm{H}{}{} \rightarrow inv)<0.69\% \; (0.34\%)$ at 90% C.L. [29, 30]. However, direct comparison is not possible, as these estimates are based on the study assuming CLIC running at 350 GeV [22, 23], with higher expected Higgs production cross section and lower cross sections for main background channels. Beamstrahlung and EPA photon interactions were also not taken into account in [22, 23]. Presented results on searches for invisible Higgs boson decays at 380 GeV CLIC show that even for 4000 fb$^{-1}$ of integrated luminosity the expected sensitivity is weaker than at other Higgs factories, running at energies of 240–250 GeV, where the exclusion limits of the order of 0.22–0.27% are expected (for first stages of FCC-ee, ILC and CEPC) [6]. That is why we focus on the search for production of new scalars in the following.

4 Results

In Fig. 6a and b presented are 95% C.L. limits on the cross section for the production of the new scalar $\mathrm{H}{}{} '$ in association with the $\mathrm{Z}{}{}$ boson, relative to the expected cross section for the production of the SM-Higgs boson (for given mass), as a function of the assumed scalar mass. The results shown in Fig. 6a and b are obtained for CLIC running at 380 GeV and 1.5 TeV, respectively.

Fig. 6

Expected limits on the production cross section of the new scalar $\mathrm{H}{}{} '$, relative to the expected SM Higgs production cross section, as a function of its mass, for CLIC running at 380 GeV (left) and 1.5 TeV (right). New scalar is assumed to have only invisible decay channels, BR$(\mathrm{H}{}{} '\rightarrow inv) = 100\%$

Decays of the new scalar are assumed to be dominated by invisible channels, BR$(\mathrm{H}{}{} '\rightarrow inv) = 100\%$. The results indicate that the experiment at CLIC will be able to exclude new scalar production with rate of about 1% of the SM production cross section for masses up to about 200 GeV, assuming 4000 fb$^{-1}$ of data collected at 380 GeV. For higher masses the experimental sensitivity decreases mainly due to the decreasing production cross section.

For the second CLIC stage, sensitivity to production and invisible decays of the light Higgs-like scalars is smaller than at 380 GeV, mainly due to the decreasing signal cross section and higher background levels. The expected limit on the invisible decays of SM Higgs boson is about 3%. Assuming the production cross section given by the SM predictions, the second stage of CLIC will be sensitive to the new ‘invisible’ scalars up to about 1 TeV.

Fig. 7

Expected sensitivity of CLIC running at 380 GeV and 1.5 TeV compared to the existing limit from LEP [37] and the expected sensitivity of ILC running at 250 GeV and 500 GeV [24]. Limits on the production cross section of the new scalar $\mathrm{H}{}{} '$, relative to the expected SM Higgs production cross section, are shown as a function of its mass. For CLIC limits, new scalar is assumed to have invisible decay channel only, BR$(\mathrm{H}{}{} '\rightarrow inv) = 100\%$, while LEP and ILC results are decay-mode independent

In Fig. 7, the expected sensitivity of CLIC running at 380 GeV and 1.5 TeV is compared to the existing limit from LEP [37] and the expected sensitivity of ILC for 2000 fb$^{-1}$ collected at 250 GeV and 4000 fb$^{-1}$ collected at 500 GeV [24]. LEP and ILC limits were evaluated in a decay-mode independent approach, based on the reconstruction of leptonic $\mathrm{Z}{}{}$ boson decays ($\mathrm{Z}{}{} \rightarrow \mathrm{e}{}{} \mathrm{e}{}{}$ and $\mathrm{Z}{}{} \rightarrow \upmu {}{} \upmu {}{}$). The two approaches are complementary. Stronger limits can be obtained when considering hadronic $\mathrm{Z}{}{}$ boson decays, if we can assume that invisible decay channels dominate for $\mathrm{H}{}{} '$. Production of new scalars with masses below 125 GeV was not considered in the presented study, as signal–background separation becomes more difficult in this range and fast detector simulation with Delphes could be not detailed enough to model this correctly.

5 Interpretation

The expected limits on invisible decays of the 125 GeV Higgs boson and limits on the production of new ‘invisible’ scalars, which were obtained in a model-independent approach, can also be used to constrain different BSM scenarios. We demonstrate the possibility of constraining parameters of the Higgs-portal models taking the VFDM model [18, 19] as an example. The Standard Model (SM) is extended by the spontaneously broken extra $U(1)_X$ gauge symmetry and a Dirac fermion. To generate mass for the dark vector $X_\mu$, the Higgs mechanism with a complex singlet S is employed in the dark sector. Dark matter candidates are the massive vector boson $X_\mu$ and two Majorana fermions $\psi _\pm$. The spontaneous symmetry breaking in the dark sector results in an additional scalar state $\phi$. This state can mix with the SM Higgs field h implying existence of two mass eigenstates:

$$\begin{aligned} \left( \begin{array}{c} \mathrm{H}{}{} \\ \mathrm{H}{}{} ' \end{array} \right) = \left( \begin{array}{cc} \cos \alpha &{} \sin \alpha \\ -\sin \alpha &{}\cos \alpha \end{array} \right) \left( \begin{array}{c} h\\ \phi \end{array} \right) \; , \end{aligned}$$

where we assume that $\mathrm{H}{}{}$ is the observed 125 GeV state. If $\alpha \ll 1$, it is SM-like, but it can also decay invisibly (to dark sector particles) via the $\phi$ component (BR$(\mathrm{H}{}{} \rightarrow inv)\sim \sin ^2\alpha$). If $\mathrm{H}{}{} '$ is also light, it can be produced in e$^+$e$^-$ collisions in the same way as the SM-like Higgs boson. We assume in the following that invisible decays to dark matter sector particles dominate for $\mathrm{H}{}{} '$ (BR$(\mathrm{H}{}{} '\rightarrow inv) \approx 100\%$). If this is the case,⁹ the cross section for new scalar production corresponding to the limits presented in the previous section can be written as:

$$\begin{aligned} \sigma _{\text {H}'} = \sigma ^{m_{\text {H}_{\text {SM}}} = \; m_{\text {H}'}}_{\text {SM}} \cdot \sin ^{2}(\alpha ) \; , \end{aligned}$$

where $\sigma _{\text {H}'}$ is the cross section for the production of a new scalar of mass $m_{\text {H}'}$, and $\sigma ^{m_{\text {H}_{\text {SM}}} = m_{\text {H}'}}_{\text {SM}}$ is the cross section for the production of the Higgs boson in the Standard Model with the same mass. Limits on the sine of the mixing angle, $\sin \alpha$, resulting from the cross section limits presented in Fig. 6, are shown in Fig. 8.

Fig. 8

Expected limits on the sine of the scalar sector mixing angle, $\sin \alpha$, as a function of the $\mathrm{H}{}{} '$ mass, for CLIC running at 380 GeV (left and right plots) and 1.5 TeV (right plot)

The mixing angle in the VFDM model can also be constrained by analysing the limit on the invisible branching ratio for the SM-like Higgs boson BR$(\mathrm{H}{}{} \rightarrow inv)$. When the contribution of the $\mathrm{H}{}{} '\mathrm{H}{}{} '$ decay channel can be neglected, the invisible partial width of the Higgs boson, $\Gamma _{\text {inv}}$, is proportional to $\sin ^{2}(\alpha )$, but depends also on other model parameters, in particular on the dark sector coupling constant, $g_{X}$, the mass of the vector dark matter, $m_X$, and the masses of the fermionic dark matter particles, $m_{\Psi _{-}}$ and $m_{\Psi _{+}}$.¹⁰ Constraints on the scalar sector mixing angle, resulting from the limits on the invisible decays of the SM-like Higgs boson expected at 380 GeV CLIC, are shown in Fig. 9. Expected limits on $\sin \alpha$, plotted as a function of the $\Psi _{-}$ particle mass, are based on the invisible decay widths¹¹ calculated with Whizard, for $g_X$ = 1 and $m_X$ = $m_{\Psi _{+}}$ = 200 GeV.¹² Also indicated in Fig. 9 are the indirect limits on the mixing angle, which can be set at CLIC from the analysis of the Higgs coupling measurements. Due to the mixing with the $\phi$ state, all couplings of the SM-like Higgs boson to SM particles are scaled by $\cos (\alpha )$. In particular, the coupling of the Higgs boson to the $\mathrm{Z}{}{}$ bosons, $g_{_{\text {HZZ}}}$, is given by:

$$\begin{aligned} g_{_{\text {HZZ}}} = g_{_{\text {HZZ}}}^{_{(\text {SM})}} \cdot \cos (\alpha ). \end{aligned}$$

It is expected that the experiment at 380 GeV CLIC will be able to measure $g_{_{\text {HZZ}}}$ in a model-independent approach with an accuracy of 0.6% [29]. If no deviations from SM are observed, the corresponding 95% C.L. limit on the mixing angle in the VFDM model is:

$$\begin{aligned} |\sin (\alpha )| < 0.14. \end{aligned}$$

Fig. 9

Expected upper limits on the sine of the scalar sector mixing angle, $\sin \alpha$, as a function of the $\Psi _{-}$ particle mass of the VFDM model, for 1 ab$^{-1}$ (magenta) and 4 ab$^{-1}$ (green) of integrated luminosity collected at 380 GeV CLIC. Indicated by blue and red curves are limits corresponding to the direct limits from ATLAS [4] and CMS [5] experiments, respectively. The horizontal lines indicate the indirect limits expected from the measurement of the $g_{_{HZZ}}$ and $\kappa$ couplings [29, 38]

If the Higgs coupling fit is performed with the assumption that the Higgs boson couplings to all SM particles scale by the same factor, $\kappa$, much stronger constraints can be set [38]. After three CLIC running stages, the overall scaling of the Higgs boson couplings should be known to

$$\begin{aligned} \Delta \kappa = 0.06\,\% \, . \end{aligned}$$

This corresponds to 95% C.L. limit on the mixing angle in the VFDM model of:

$$\begin{aligned} |\sin (\alpha )| < 0.044. \end{aligned}$$

Also indicated in Fig. 9 are constraints resulting from the ATLAS and CMS limits on invisible Higgs boson decays, BR$(\mathrm{H}{}{} \rightarrow inv) < 13\%$ [4] and $< 19\%$ [5], respectively. For masses of dark matter particles up to about 60 GeV, direct search for invisible Higgs boson decays will allow to set much better constraints on the mixing angle in the scalar sector than it will be possible with indirect methods.

6 Summary

We studied the sensitivity to invisible Higgs boson decays and the possibility of constraining production of new scalar particles at CLIC running at 380 GeV and 1.5 TeV. We assumed associated production of Higgs-like neutral scalar with $\mathrm{Z}{}{}$ boson and invisible scalar decays. The analysis was based on the Whizard event generation and fast simulation of the CLIC detector response with Delphes, taking into account beam energy profile as well as background contributions from photon–photon and electron–photon interactions. An approach consisting of a two-step analysis, with multivariate analysis methods employed at the second step, was used to optimise separation between signal and background processes. Expected limits on the production cross section of the new scalar $\mathrm{H}{}{} '$ were presented as a function of its mass, for CLIC running at 380 GeV and 1.5 TeV. For 1000 fb$^{-1}$ of data collected at initial CLIC stage invisible Higgs boson decays at the level of 1.0% can be excluded at 95% C.L. The limits obtained in the model-independent approach can also be used to set limits on the different extensions of the Standard Model. Constraints at the percent level can be set on the scalar sector mixing angle in Higgs-portal models.