Exotic decays of the 125 GeV Higgs boson at future $e^+e^-$ lepton colliders

Discovery of unexpected properties of the Higgs boson offers an intriguing opportunity of shedding light on some of the most profound puzzles in particle physics. The Beyond Standard Model (BSM) decays of the Higgs boson could reveal new physics in a direct manner. Future electron-positron lepton colliders operating as Higgs factories, including CEPC, FCC-ee and ILC, with the advantages of a clean collider environment and large statistics, could greatly enhance the sensitivity in searching for these BSM decays. In this work, we perform a general study of Higgs exotic decays at future $e^+e^-$ lepton colliders, focusing on the Higgs decays with hadronic final states and/or missing energy, which are very challenging for the High-Luminosity program of the Large Hadron Collider (HL-LHC). We show that with simple selection cuts, $O(10^{-3}\sim10^{-5})$ limits on the Higgs exotic decay branching fractions can be achieved using the leptonic decaying spectator $Z$ boson in the associated production mode $e^+e^-\rightarrow Z H$. We further discuss the interplay between the detector performance and Higgs exotic decay, and other possibilities of exotic decays. Our work is a first step in a comprehensive study of Higgs exotic decays at future lepton colliders, which is a key ingredient of Higgs physics that deserves further investigation.


Introduction
The recent discovery of the Higgs boson by the AT-LAS and CMS experiments at the Large Hadron Collider (LHC) opened a new era of particle physics [1,2]. The Standard Model (SM)-like Higgs boson has a very deep connection to many profound puzzles of fundamental physics, such as the hierarchy problem and the naturalness problem, the nature of dark matter, the origin of neutrino mass, the origin of the fermion mass hierarchy (the flavor problem), the essence of the electroweak phase transition and electroweak bayrogenesis. Essentially all beyond Standard Model (BSM) solutions to these puzzles, such as supersymmetry (SUSY) [3][4][5], composite Higgs models [6][7][8][9], and grand unified theories [10], predict modifications of the properties of the Higgs boson. Hence, the precise measurements of the properties of this Higgs boson have great potential to shed light on BSM physics.
The HL-LHC will measure many SM model decay modes of the Higgs boson to a relative precision of O(10%) [11][12][13][14][15]. The future lepton colliders operating as Higgs factories, with their clean collider environment and large statistics, would measure the Higgs boson couplings to a relative precision of O(0.1% ∼ 1%) [16][17][18]. Many discussions about the physics potential of future lepton colliders focus on the precision measurement of the Higgs properties in the effective field theory (EFT) framework [19][20][21][22][23][24]. However, new physics could manifest itself through Higgs exotic decays if some new light degrees of freedom are present, which are not described by the SM EFT. Hence, systematically searching for Higgs exotic decays would be an important physics component of the future lepton collider programs. Moreover, since many of these future facilities are currently at different stages of planning, investigation of these new physics potentials could impact their designs to achieve some more comprehensive physics goals.
In Section 2, we present an overview of exotic decay searches at lepton colliders and some general discus-sions of the Higgs exotic decays at different future lepton colliders. In Section 3, we describe our simulation framework and present our phenomenological analysis for various Higgs exotic decay modes. We summarize the physics potential from the Higgs exotic decays at the (HL-)LHC and the future lepton collider programs in Section 4. In our summary table, we include comprehensive projections and show the complementarity between future lepton collider programs and the HL-LHC. We also discuss many important future directions for the Higgs exotic decay programs.
2 Theoretical framework

Higgs exotic decay modes considered in this work
The Higgs boson BSM decays have a rich variety of possibilities. To organize this study on Higgs boson BSM decays, we selectively choose a set of phenomenologically driven processes. We focus on two-body Higgs decays into BSM particles, which are allowed to subsequently decay further, up to four-body final states. We only consider the Higgs boson as an CP-even particle. CPviolation effects would affect various differential distributions, and this demands future study. These processes are well-motivated by SM+singlet extensions, two-Higgsdoublet-models, SUSY models, Higgs portals, gauge extensions of the SM, etc. These assumptions have also been emphasized in the recent overview of Higgs exotic decays [25] and the CERN yellow report [26].
We consider in general the exotic Higgs decays into BSM particles dubbed as X i , h → X 1 X 2 . The cascade decay modes are classified into four cases, schematically shown in Fig. 1. We discuss their major physics motivation and features at lepton colliders in order. h → 2: The X i s in this case are detector-stable and charge-neutral. * They could be dark matter candidates. The Higgs portal [27] to dark matter models, including various SUSY light dark matter models [28][29][30][31][32][33][34][35][36][37][38], motivates this BSM search channel. The lepton collider background for this channel are mainly from the process e + e − → ZZ → Z + νν and e + e − → W + W − → + − νν. This channel, due to its simplicity and importance, has been studied by most of the future lepton collider programs [16-18] and we will quote these results in our summary table. We include this channel here for completeness. In addition, many of the models that motivate this channel also induce other Higgs exotic decays we consider in this study.
h → 2 → 3 → 4: This is the topology in which X 1 is detector-stable and X 2 decays to two particles, with one of these decay products further decaying into two particles. A typical BSM model for such decay modes is the Higgs decaying into the lightest supersymmetric particle (LSP) plus a heavier neutralino, which subsequently decays into the LSP plus a resonant BSM particle. This resonant BSM particle could be a singletlike scalar in the Next-to-Minimal-Supersymmetric-Standard-Model (NMSSM). Many SUSY models which motivate Higgs invisible decays also induce this decay channel, e.g. [38,39]. It also commonly exists in the so-called "stealth SUSY" models [40]. This singlet-like scalar decays into SM fermion pairs, giving rise to the final state of a pair of resonant SM particles plus missing energy, dubbed h → (f f ) + / E T . † In this study, we only consider the channels which are very challenging at the LHC, h → (jj) + / E T , h → (bb) + / E T and h → (τ + τ − ) + / E T . For the hadronic channels, the major background is from the SM Higgs decay modes h → ZZ * → jj +νν and h → ZZ * → bb+νν. For h → (τ + τ − )+/ E T , in addition, the SM Higgs decay h → W W * → τ + τ − +νν also contributes to the background. h → 2 → (1 + 3): This is the topology when X 1 is detector-stable and X 2 decays into three particles. The typical BSM model is very similar to the previous decay topology. The difference comes from X 2 decaying into three particles without an intermediate on-shell resonance. This scenario takes place very naturally if the singlet-like scalar is heavy, or the particle X 2 decays through an off-shell Higgs or Z-boson. The final state of this Higgs exotic decay signature would be h → f f +/ E T . ‡ The lepton collider background sources are similar to the previous topology as well. h → 2 → 4: In this channel, the Higgs decays to a pair of BSM particles, both of which subsequently decay into two final state particles. There are a wide range of BSM models which give rise to such a decay pattern, and we selectively discuss several benchmark cases. The intermediate particle could be a pair of vectors from the "dark photon" or "dark Z-prime " models [43]. These models commonly include a new gauge boson that couples to the SM with suppressed strength. A typical example is a small kinetic mixing with the hypercharge field strengths, but more general couplings are certainly possible. They not only induce h → Z Z decays but also sizable h → ZZ decays in which the masses of the intermediate particles are uneven, if a source of mass mixing is allowed. An important feature of such models relevant for the phenomenology of Higgs exotic decay is that the Z will have sizable O(few%) decay branching fractions to SM charged leptons. Consequently, this scenario could be severely constrained by the (HL-)LHC searches, unless Z is leptophobic. This option requires more contrived model building to survives collider direct search limits [44]. The intermediate particle could also be a pair of scalars from SM+scalar, 2HDM+scalar, NMSSM models, etc [25]. In addition, a dark sector with strong dynamics, such as the Twin Higgs [45] models, could also give rise to the Higgs decays into a pair of spin-0 composite particles. § The intermediate scalar decays strongly prefer SM heavy fermions and thus bb, cc and τ + τ − decays would dominate. These decay modes are very hard to probe by (HL-)LHC searches due to the large background in the hadron collider environment. There is the interesting possibility for the intermediate scalars to decay into diphoton pairs that we consider as well. Hence, we consider many combinations of the (bb), (cc), (jj), (τ + τ − ) and (γγ) decays of the intermediate particles.
The backgrounds are again mainly from SM Higgs decays into four particles through SM gauge bosons. For final states involving photons, the SM electroweak processes at the lepton colliders dominate the background.
There are several other decay topologies in Ref.
[25] that we do not include in our current study. We comment on them here. h → 2 → 3: In this case, X 1 is detectorstable and X 2 decays promptly. For instance, in dark photon models the SM Higgs decays into h → γZ and Z could subsequently decay into SM particles via two-body decay. In certain SUSY scenarios, the Higgs could decay to the LSP and the next-to-lightest-supersymmetricparticle, which subsequently decays into a photon plus the LSP. h → 2 → 4 → 6 and h → 2 → 6: the direct decay product from the Higgs undergoes a decay chain or decays into a three-body final state. These decay topologies are common for (R-parity-violating) SUSY models. These decay modes are well-motivated and should be studied in follow-up works.

Higgsstrahlung process
For future lepton colliders running at the center of mass energy 240 ∼ 250 GeV, the most important Higgs production mechanism is Z-Higgs associated production through an off-shell Z boson e + e − → Z * → Zh. The Z boson with visible decays plays the role of Higgs spectator and enables Higgs tagging using the "recoil mass" technique. Given the known initial state energy ¶ , subtracting the Z-boson four-momentum enables the reconstruction of the Higgs four-momentum and thus its invariant mass. The so called recoil mass defined in this way sharply peaks at the Higgs boson mass. A selection cut around this peak would remove the majority of the SM background.
For an unpolarized electron-positron beam at the center of mass energy 240 GeV, the Higgs production rate is around 230 fb [46,47]. Both the CEPC and FCCee plan to mainly run with unpolarized beam at this energy. Here, we consider the CEPC running scenario with an effective integrated luminosity of 5 ab −1 , following its current plan of two interaction points and ten years of running with the designed beam luminosity. The CEPC will produce 1.15 million Higgs bosons in this production channel. The FCC-ee running scenario we consider here has six times more statistics than the CEPC, 30 ab −1 (equivalently 6.9 million Higgs bosons), following its current plan of four interaction points with a higher beam luminosity. For the ILC, we only consider the limits from its 250 GeV runs, in the H20 scenario ‡ To explicitly distinguish Higgs exotic decays with or without resonances in the final state particles, we put the pair of SM particles that form a resonance in parenthesis. We follow this notation in describing Higgs exotic decay final states throughout this study. § However, in many cases, this composite particle "glu-ball" would be meta-stable and requires a more sophisticated phenomenological study for displaced decays.
¶ Corrections from beamstrahlung effect [41] and ISR effect [42] need to be carefully taken into account of 2 ab −1 integrated luminosity with beam polarization of p(e − , e + ) = (−0.8, +0.3). The Higgs production rate is enhanced by a factor 1.4 due to beam polarization compared to the unpolarized beams of circular lepton colliders. In this scenario, considering the 250 GeV runs alone, the ILC will produce 0.64 million Higgs bosons in this channel. Before proceeding to the next section of detailed numerical analysis for individual channels, we comment on several instructive scenarios for the future lepton collider sensitivities here. If we consider the cleanest e + e − → Zh, Z → + − mode, the Higgs boson can be tagged with little background from the SM. Taking this leptonic-decaying spectator Z-boson alone and CEPC as an example, we will have 7.7×10 4 Higgs bosons, naively reaching a very impressive 4 × 10 −5 (2.5 × 10 −4 ) limit on the Higgs exotic branching fraction for the case of zero (one hundred) SM background. In this study, we choose to study this clean leptonic spectator Z-boson mode with various Higgs exotic decay modes. For most Higgs exotic decay modes, further inclusion of hadronic decaying Z boson and even invisible Z will definitely improve the limits significantly. The (HL-)LHC will produce more Higgs bosons, providing excellent limits on Higgs decaying into leptons, such as h → ( + − )( + − ), reaching better than O(10 −5 ) branching fractions. Hence, we do not consider these pure leptonic channels at future lepton colliders but focusing only on the channels that are hard for the LHC, involving hadronic decays and/or missing energy.

Phenomenological analysis
Following the discussion in the previous section, we perform a numerical study of the future lepton collider reach for selected Higgs exotic decay modes. Though operating at slightly different center of mass energies in the range of 240 ∼ 250 GeV, the proposed future lepton colliders in general have similar detector performance. For simplicity, we choose the CEPC as the benchmark accelerator and detector model in this section for the analysis. We will extrapolate the sensitivity for FCC and ILC in the next section.
For numerical analyses, we generate both the signal and the background events for an 240 GeV electronpositron collider with MadGraph5 at parton level [48] and impose the detector acceptance, energy and momentum smearing, and separation cuts with our own analysis code.
We describe here our parameter choices for the detector effects, and our pre-selection cuts that are universal for the analyses for all Higgs exotic decay mode. To reach a high particle identification efficiency, all of the visible particles in the final state are required to have | cos θ| < 0.98 (following Ref. [16]), or equivalently |η| < 2.3. The final state particles are required to be well separated with We only study the case where the Z boson decays into + − final state where ± = e ± , µ ± , and leave the study of other decay modes of the Z boson for future works. The signal events are required to contain at least a pair of opposite-sign same-flavor charged leptons with an opening angle greater than 80 • , and satisfy and |m − m Z | < 10 GeV.
Furthermore, the recoil mass is defined as To suppress the ISR contribution to the backgrounds, for Higgs exotic decay modes without missing energy, we require the events to have the total visible energy E vis > 225 GeV. (6) In this work, we mimic the detector resolution effect by adding Gaussian smearing effects on the fourmomentum of the particles, following the performance described in Ref. [16]. For photons in the final state, the energy resolution is determined by the electromagnetic calorimeter, which performs approximately as The energy resolution of jets is affected by the hadron calorimeter, and performs approximately as For electrons and muons in the final state, we include the momentum resolution effect of the track system with the approximate performance of Next, we discuss the phenomenology of individual Higgs exotic decay channels.
This final state appears if the SM-like Higgs boson decays into X 2 X 1 with X 2 → X 1 s and s → jj. In the NMSSM, for example, the particles X 1 , X 2 and s could be identified as the light neutralinosχ 0 1 ,χ 0 2 and the light singlet-like (pseudo-)scalar h 1 (a 1 ), respectively. We generate the irreducible SM background e + e − → + − ν ν jj with MadGraph5. Beyond the pre-selection cut and the recoil mass cut, we require that there are two additional jets which satisfy E j > 10 GeV and | cos θ j | < 0.98.
After these cuts, the invariant mass distribution of the dilepton system is shown in Fig. 2. The dominant background after the recoil mass cut will clearly be the Higgsstrahlung process with h → ZZ * → qqνν. After the recoil mass cut, the SM background cross section is 0.063 fb. The dijet invariant mass (m jj ) distribution and the two-dimensional differential distribution of m jj versus / E T of the SM background after the recoil mass cut are shown in Fig. 3. There is a clear valley in the distribution between 35 to 75 GeV, in which none of the Z bosons from the SM-like Higgs boson decay are on-shell and thus the h → qqν ν is doubly suppressed. This property could be used to optimize the cut and increase the sensitivity to the signal if the invariant mass of the light (pseudo)scalar falls in this range. We use the likelihood function of the m jj -/ E T distribution to give the exclusion limit. We show the 95% C.L. exclusion limit in Fig. 4 in the plane of X 1 , mass m 1 , and the mass splitting between X 2 and X 1 , m 2 −m 1 , for two benchmark intermediate scalar masses of 10 GeV (a) and 40 GeV (b). In most of the parameter space, the Higgs branching fraction to this exotic decay channel h → (jj) + / E T can be excluded to the level of 2-6×10 −4 . We discuss several kinematical features here that impact the exclusion limit. When m s = 10 GeV, the larger the m 2 − m 1 , the better the reach. This is due to the fact that the signal events populate the lower-left corner in the m jj -/ E T plane with low SM background for large mass splitting when the MET is small. Consequently, the highest sensitivity is reached when m 1 = 10 GeV, m 2 = 100 GeV. When m 2 − m 1 is small, the signal events will tend to have a large / E T and look like the SM background events e + e − → Z(h → ZZ * ) where the on-shell Z from the Higgs boson decay decays to νν. Furthermore, for a light intermediate scalar mass, small mass splitting also results in soft jets in the final states which are more likely to fail the pre-selection cuts of two jets. When m h 1 = 40 GeV (so m 2 − m 1 40 GeV), the m jj falls in the valley and the sensitivity is high (around 2 × 10 −4 ) as we expected. When the intermediate scalar mass m h 1 is close to the Z boson mass, the sensitivity is low, due to the relatively large SM background e + e − → Z(h → ZZ * ) where the on-shell Z from the Higgs boson decay decays to dijets. Fig. 4. The 95% C.L. upper limit on the Higgs exotic decay branching fractions into (jj) + / ET for various lightest detector-stable particle X1 with mass m1 and mass splittings m2 −m1. The results for the benchmark cases of the dijet mother particle mass of 10 GeV and 40 GeV are shown in (a) and (b), respectively.

h → (bb) + / E T
The background and the benchmark model for this mode are the same as the h → (jj) + / E T case discussed in the previous section. The signal event is required to contain two b-tagged jets. The b-tagging efficiency is conservatively chosen to be 80%, and the charm mis-tagging rate and the light flavor mis-tagging rate are set to be 9% and 1%, respectively. Similarly, we use the likelihood function of the m bb -/ E T distribution to derive the exclusive limit.
We show the 95% C.L. exclusion limit in Fig. 5 in the plane of X 1 , mass m 1 , and the mass splitting between X 2 and X 1 , m 2 −m 1 , for two benchmark intermediate scalar masses of 10 GeV (a) and 40 GeV (b). In most of the parameter space, the Higgs branching fraction to this exotic decay channel h → (bb) + / E T can be excluded to the level of 5 × 10 −5 ∼ 1.5 × 10 −4 . The features for various kinematical regions are also similar to the analysis in the previous section. The limits are roughly a factor of 4 better than h → (jj)+/ E T across the whole parameter region, due to the b-tagging reducing the flavor universal quark jets background. In the highest sensitivity benchmark points, the 95% C.L. exclusion bound can reach 6×10 −5 (e.g., m 1 = 10 GeV, m 2 = 60 ∼ 100 GeV, m s = 40 GeV). This impressive result nearly reaches the statistical limit of CEPC, after folding in a factor 0.64 on signal strength from the requirement of double b-tagging. ET for various lightest detector-stable particle with mass m1 and mass splittings m2 −m1. The results for the benchmark cases of a dijet mother particle mass of 10 GeV and 40 GeV are shown in (a) and (b), respectively.

h → jj + / E T and h → bb + / E T
Although the final states of these channels are the same as the h → (jj) + / E T and h → (bb) + / E T cases discussed in the previous sections, the distributions of the kinetic variables are quite different since there is no dijet resonance in this case. We assume the intermediate particle (e.g., scalar s) to be very heavy so that the decay of the X 2 can be fully described by a four-fermion contact operator. Similarly, we use the likelihood function of the m bb -/ E T distribution to derive the exclusive limit. The results are shown in Fig. 6 in the plane of X 1 , mass m 1 , and the mass splitting between X 2 and X 1 , m 2 −m 1 , for h → jj+/ E T (a) and h → bb+/ E T (b). Comparing with the Higgs exotic decays with intermediate resonance in the previous section, the exclusion limits on the branching fraction are only slightly worse in the bulk region of the parameter space, reaching generally 3 × 10 −4 ∼ 8 × 10 −4 and 2 × 10 −4 ∼ 4 × 10 −4 for h → jj + / E T and h → bb + / E T , respectively. Fig. 6. The 95% C.L. upper limit on the Higgs exotic decay branching fractions into jj +/ ET (a) and bb+/ ET (b) for various lightest detector-stable particle mass m1 and mass splittings m2 − m1.

063102-7
From the exclusion limits shown in Fig. 6, we find that when the mass splitting m 2 −m 1 is around 80 GeV, the future lepton colliders have the strongest sensitivities on these Higgs exotic channels, reaching around 3.1×10 −4 and 1.6×10 −4 for h → jj+/ E T and and h → bb+/ E T , respectively. When X 1 is light and m 2 −m 1 is large, the energy is shared by the two jets and the X 1 . Consequently, when the mass splitting m 2 − m 1 is around 80 GeV, the dijet invariant mass will be around 40∼60 GeV, falling in the "valley" of low SM background as shown in Fig. 3. For heavier X 1 , the MET will be lower due to less momentum available for the LSP. The optimal limits will be reached for an even smaller mass splitting. The bb + / E T case has a higher sensitivity again by roughly a factor of two since the b-tagging suppresses the SM background.

h → (jj)(jj), h → (cc)(cc) and h → (bb)(bb)
For this class of Higgs exotic decays, we consider the scalar mediator (s), the pseudoscalar (a) and the vector (Z µ ) mediator. We assume the effective interactions between the SM-like Higgs boson and the mediators are hss, haa and hZ µ Z µ , respectively. The (pseudo)scalar mediator can decay into dijet final states via sf f (af γ 5 f ) or sG µν G µν (aG µνG µν ) interactions. For the vector mediator case, we consider both the vector-like and righthanded interaction with the SM fermions. In addition to the pre-selection cuts and the recoil mass cut, we require that there are at least four jets that satisfy E j > 5 GeV.
The most important background from the SM is where the four jets could be either from the hadronic decay of the SM vector bosons in h → V V * , or from h → jj with jet-splitting. We show the recoil mass distribution of the SM background with all but the recoil mass cut applied in Fig. 7. We also include in the red curves the background distribution with the inclusion of ISR effect. Its effect is negligibly small after the relatively large window of the recoil mass cut. We hence neglect the ISR effect in this analysis. Another kinematical variable which is useful for separating the signal and background is With the combination which gives δm, we calculate the likelihood function of the m j 1 j 2 +m j 3 j 4 versus δm distribution and get the 95% C.L. exclusive bound shown in

063102-8
If the mediator decays to the light flavors (u, d, s, c, g) but not the bottom quark, a b-veto could be used to suppress the SM background. In this case, we require that there is no b-tagged jet in the final state. The b-veto does not increase the sensitivity significantly, as shown in Fig. 8(b), because the SM background is dominated by the Z(h → V V * → jjjj), which has a similar composition of quark flavors.
The future lepton collider with 5 ab −1 integrated luminosity could exclude 10 −3 branching fractions of h → (jj)(jj) in a wide range of the mediator mass. We do not consider the mediator mass below 10 GeV as a different analysis strategy should then be taken because the jets are more likely to fail the separation cuts, manifested in the left-hand corner of Fig. 7.   The sensitivity in this channel is worse than in h → (jj) + / E T and h → jj + / E T . This is because we miss the information of the correct combination of jets in the final state. The wrong combination of jets in the final state from the SM background will mimic the signal.
For the case of di-bottom pair resonances h → (bb)(bb) and di-charm pair resonances h → (cc)(cc), the simulation is nearly identical to the h → (jj)(jj) case. We require at least three b-tagged jets and c-tagged jets in the final state, respectively. For the charm-tagging, we assume a tagging efficiency of 60% and mis-tag rate from b-jets 15% and light jets 10%. The relatively larger fake rate assumed here leads to slightly worse sensitivity for h → (cc)(cc) when comparing with h → (bb)(bb). The results are shown in Fig. 9. The future lepton collider with 5 ab −1 integrated luminosity could exclude branching fractions of h → (bb)(bb) and h → (cc)(cc) down to 3 × 10 −4 ∼ 4 × 10 −4 and 7 × 10 −4 ∼ 9 × 10 −4 , respectively, in a wide range of the mediator mass.

h → (γγ)(γγ)
For this scenario, we consider the (pseudo)scalar mediator. In addition to the preliminary cuts, we require there are at least four hard photons that satisfy The most important background from the SM is We show the SM background distribution in the recoil mass after imposing all the other kinematical cuts in Fig. 10. The contribution from e + e − → Z + h + 2γ → + − + 4γ is highly suppressed by the h → γγ decay branching ratio. The SM Z boson decaying into a four photon final state has a tiny branching ratio and thus it The cross section of e + e − → Z + h → + − + 2γ is less than 40 ab and is severely suppressed by two more hard photons in the final state, so it is neglected here.
does not contribute. Therefore, there is no sizable resonance structure in the recoil mass distribution in the SM background of this channel, as shown in this figure.
Beyond the pre-selection cuts, the four photons should form two pairs of diphoton resonances with similar masses and thus we further require for a signal event. The SM background is 0.14×10 −3 fb. Around three signal events are needed to be excluded at 95% C.L. for 5 ab −1 integrated luminosity. We tabulate the signal efficiency for various intermediate masses and show the derived 95% C.L. limits on this Higgs exotic branching fraction Br(h → (γγ)(γγ)) in Table 1. The signal acceptance for this analysis increases as the mediator mass increases, as a result of the requirement on the photon energy in Eq. (14). The limits can reach up to 4.7 × 10 −4 on the exotic branching fractions of this channel for a mediator mass of 50 GeV.

h → (jj)(γγ)
For this scenario, we consider the (pseudo)scalar mediator. In additional to the preliminary cuts, we require there are at least two hard photons and two hard jets that satisfy E γ,j > 10 GeV.
The most important background from the SM is The contribution from e + e − → Z + h → + − + bb → + − + bb + 2γ is suppressed by the fine structure constant and the separation cut on y in Eq. (1). We show the SM background distribution in the recoil mass with the all but the recoil mass cut in Fig. 11. The background is mostly from e + e − → ZZ + 2γ where the photons are from the ISR. There is also a sizable contribution from e + e − → ZZ where one of the Z boson decays into a qq final state and the photons are from the final state radiation from the jets. We can find a small peak at the Z-pole in the recoil mass distribution in this figure. We also require the dijet invariant mass and diphoton invariant mass to be close in value, to further suppress the background. The SM background is 0.31 × 10 −3 fb. Around four signal events are needed to be excluded at 95% C.L. for 5 ab −1 integrated luminosity. We tabulate the signal efficiency for various intermediate masses and show the derived 95% C.L. limits on this Higgs exotic branching fraction Br(h → (jj)(γγ)) in Table 2. We can see the signal selection efficiency for this analysis increases as the mediator mass increases, as a result of the requirement on final state particle energies in Eq. (17). The limits can reach up to 5.6×10 −4 on the exotic branching fractions of this channel for a mediator mass of 50 GeV.

Summary and outlook
We summarize the set of Higgs exotic decays in Table 3, including current collections and projections of LHC constraints, and limits from our study for the CEPC, ILC and FCC-ee. For the LHC constraints, we tabulate both the current limits and projected limits on these exotic decay channels from various references. Due to the subtleties in the projections at the LHC from 063102-10 possible systematics, we only quote the projected limits from studies with highest integrated luminosity and label them accordingly, instead of naively scaling these results to the projected luminosity of (HL-)LHC. The current LHC limits are collected in the second column, and the square brackets represent 7 and 8 TeV LHC results alone. The projections for HL-LHC are collected in the third column, where the limits for 100 fb −1 and 300 fb −1 are shown in parentheses and square brackets respectively. The reference for the current limits and projected limits are also included. For projections on the reaches of HL-LHC, it should be understood that they are from Ref. [25] unless specified otherwise. The limits for the case of the ILC and FCC-ee are extrapolated from our CEPC analyses in Section 3 following the running scenario discussion in Section 2, assuming the reach is dominated by statistics. The decay channels involving τ -leptons are extrapolated and validated by our preliminary simulation, assuming a 40% τ -tagging efficiency with more background considered. We mark these extrapolated limits with an asterisk. Table 3. The current and projected limits on selected Higgs exotic decay modes for the (HL-)LHC, CEPC, ILC, and FCC-ee. Throughout our analysis, we assume SM production rates for the Higgs boson. The current LHC limits are collected in the second column, with the square brackets indicating 7 and 8 TeV LHC results alone. The projections for the HL-LHC are collected in the third column, where the limits for 100 fb −1 and 300 fb −1 alone are shown in parentheses and square brackets respectively. The references for the current limits and projected limits are included; for projections on the HL-LHC limits we omit the references if they are from Ref. [25]. We also omit the references for projections on the future lepton collider programs if the results are from this study. Limits extracted from related searches with reasonable assumptions are marked with asterisks following the results with details explained in the text. For more details about the benchmark parameter choices for some of the decay modes, see discussion in the main text.
Similar to the lepton collider, the LHC reach can depend on model parameters. While a comparison of the reach throughout the full parameter space is possible, it is tedious and not illuminating. We choose to focus on comparison for particular benchmark points, which is good enough to demonstrate the qualitative difference between the LHC and future lepton colliders. These choices are to make the limits more or less representative of the intermediate ("average") limits across wide parameter regions without being in the extremely good or bad kinematical points.
Photons are one of the "clean" objects in the LHC collider environment and the large statistics at the HL-LHC can constrain decays with photons to good precision. However, these is some discrepancy between the theoretical projections in Ref. [25] of 3×10 −5 at 300 fb −1 , and the extrapolation from the current limit at the 8 TeV LHC [54]. We choose to use the more conservative extrapolation by us and thus put an asterisk after the projected limit in this summary table.
For the Higgs invisible decays at lepton colliders, we quote the limits from current studies [16][17][18]. These limits do not depend on the invisible particle mass using the recoil mass technique at lepton colliders.  Table 3. We put several vertical lines in this figure to divide different types of Higgs exotic decays. Table 3 and the corresponding Fig. 12, we can clearly see the improvement in exotic decays from the lepton collider Higgs factories. These exotic Higgs decay channels are selected such that they are hard to be constrained at the LHC but important for probing BSM decays of the Higgs boson. The improvements on the limits of the Higgs exotic decay branching fractions vary from one to four orders of magnitude for these channels. The lepton colliders can improve the limits on the Higgs invisible decays beyond the HL-LHC projection by one order of magnitude, reaching the SM invisible decay branching fraction of 0.12% from h → ZZ * → νννν [56]. For the Higgs exotic decays into hadronic particle plus missing energy, (bb) + / E T , (jj)+/ E T and (τ + τ − )+/ E T , the future lepton colliders improve on the HL-LHC sensitivity for these channels by roughly four orders of magnitude. This great advantage benefits a lot from low QCD background and the Higgs tagging from recoil mass technique at future lepton colliders. As for the Higgs exotic decays without missing energy, the improvement varies between two to three orders of magnitude, except for the one order of magnitude improvement for the (γγ)(γγ) channel. Being able to reconstruct the Higgs mass from the final state particles at the LHC does provide additional signal-background discrimination power and hence the future lepton colliders improvement on Higgs exotic decays without miss-ing energy is less impressive than for those with missing energy. Furthermore, as discussed earlier, leptons and photons are relatively clean objects at the LHC and the sensitivity at the LHC on these channels will be very good. Future lepton colliders complement the HL-LHC for hadronic channels and channels with missing energies.

From this summary in
There are many more investigations to be carried out under the theme of Higgs exotic decays. For our study, we take the cleanest channel of e + e − → ZH with Z → + − and h →exotics up to four-body final state, but further inclusion of the hadronic decaying spectator Z-boson and even invisible decays of the Z-boson would definitely improve the statistics and consequently result in better limits. As a first attempt to evaluate the Higgs exotic decay program at future lepton colliders, we do not include the case of very light intermediate particles whose decay products will be collimated, but postpone this for future study when the detector performance is more clearly defined. There are many more exotic Higgs decay modes to consider, such as Higgs decaying to a pair of intermediate particles with un-even masses [25], Higgs CP property measurements from its decay differential distributions [57][58][59][60], flavor violating decays, decays to light quarks [61], decays into meta-stable particles, and complementary Higgs exotic productions [62]. Our work is a first systematic study evaluating the physics potential of future lepton colliders to probe Higgs exotic decays. More work will be needed to obtain a complete picture.
LTW would like to thank Matt Strassler for useful discussions, and the Institute of High Energy Physics in Beijing for hospitality. ZL and LTW would like to thank the Kavli Institute for Theoretical Physics for hospitality.