Gravitational wave and collider probes of a triplet Higgs sector with a low cutoff

We study the scalar triplet extension of the standard model with a low cutoff, preventing large corrections to the quadratic masses that would otherwise worsen the hierarchy problem. We explore the reach of LISA to test the parameter space region of the scalar potential (not yet excluded by Higgs to diphoton measurements) in which the electroweak phase transition is strongly first-order and produces sizeable gravitational waves. We also demonstrate that the collider phenomenology of the model is drastically different from its renormalizable counterpart. We study the reach of the LHC in ongoing searches and project bounds for the HL-LHC. Likewise, we develop a dedicated analysis to test the key but still unexplored signature of pair-production of charged scalars decaying to third-generation quarks: pp→tb¯(t¯b),bb¯\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$pp\rightarrow t\overline{b} (\overline{t}b), b\overline{b}$$\end{document}. These results apply straightforwardly to other extensions of the Higgs sector such as the 2HDM/MSSM.


Introduction
Hyperchargeless triplet scalars arise in a variety of models of new physics. They include (i) theories of grand unification (GUT), where scalar multiplets, often transforming in the adjoint representation of the GUT group, break spontaneously the GUT symmetry. A simple example is the 24 in SU (5), which decomposes as (1, 1) 0 +(1, 3) 0 +· · · under the Standard Model (SM) gauge group SU (3) c × SU (2) L ×U (1) Y , therefore delivering a scalar triplet. Likewise, the 45 representation of SO (10) contains the 24 of SU (5) and therefore a SM triplet as well. (ii) Supersymmetric (SUSY) models. As a matter of fact, the triplet extension of the MSSM is one of the simplest options to alleviate the little hierarchy problem [1,2]. (iii) Composite Higgs models (CHM). The scalar sector of most CHMs is non minimal. It includes a hyperchargeless triplet in one of the two simplest cosets admitting an UV completion à la QCD in four dimensions, viz. SU (5)/SO (5) [3,4]. Moreover, models based on SO(7)/G 2 [5] provide exactly one triplet in addition to the Higgs boson.
Therefore, the phenomenology of such triplet is not dictated by the renormalizable Lagrangian. The latter has to be instead supplemented with effective operators encoding the effects of the heavier resonances (SUSY partners, composite states, etc.), which can modify drastically the dynamics of . To demonstrate this, we will work under the assumption that the triplet does not get a (custodial symmetry breaking) vacuum expectation value (VEV). This limit can be naturally enforced assuming the triplet is a CP-odd scalar and CP is conserved in the Higgs sector. At the renormalizable level, the Lagrangian becomes accidentally Z 2 symmetric, i.e. → − , making the neutral component of the triplet a potential dark matter candidate. The charged components are in turn long-lived. The corresponding phenomenology has been studied in Refs. [6][7][8]. However, the effective operators make all components decay promptly even if the cutoff is f ∼ several TeV at which new resonances are out of the reach of current facilities. A much larger cutoff would introduce too large corrections also to the triplet mass, worsening the hierarchy problem. 1 In this article, we study probes of current and future colliders to this more natural version of the inert triplet model (ITM).
The extended Higgs sector modifies also the electroweak (EW) phase transition (EWPT). Thus, we extend previous studies in this respect [10,11] computing the reach of future detectors for gravitational waves that originate in the production, evolution and eventual collisions of bubbles of vacuum in a first-order phase transition. The paper is organized as follows. We introduce the model in Sect. 2. We discuss the dynamics of the EWPT in Sect. 3. We explore collider signatures in Sect. 4. In Sect. 5 we propose an LHC analysis that has not been yet worked out experimentally for probing the key channel pp → tb(tb), bb. We study signal and background and provide prospects for the HL-LHC, namely the LHC running at a center of mass energy (c.m.e) √ s = 13 TeV with integrated luminosity L = 3 ab −1 . We conclude in Sect. 6.

Model
The Lagrangian of the CP-odd scalar triplet when it is assumed embedded in an UV theory, takes the form: where We will assume flavour diagonal couplings: c with c a constant and y Yukawa. We will comment on departures from this assumption in the conclusions. The relevant parameter is therefore the ratio c/ f . The productH (H ) stands for the doublet in the SU (2) L ×U (1) Y decomposition 2 −(+)1/2 × 3 0 = 2 −(+)1/2 + 4 −(+)1/2 . Explicitly: 1 (The fine-tuning scales roughly as ∼ m 2 h / f 2 , with m h the Higgs mass [9]. f = 1 TeV gives already ∼ 1 %, and this falls below the permile level for f > 4 TeV.) and analogously for H upon the replacement h * 0 → h + and h − → −h 0 . Therefore, in the unitary gauge after EWSB, we obtain: with v ∼ 246 GeV and the sum extends to the first and second families of quarks, that we will denote collectively by q. Couplings to the leptons could be also present. We neglect them in this analysis. φ 0(±) can decay into SM quarks. Likewise, for m t > m , the top quark can decay into the triplet and a bottom quark. (m t stands for the top quark mass, whereas m is the physical mass of ; at tree level m 2 = μ 2 +λ H v 2 /2.) The following relations hold: In the second equation we are assuming m q m q , which is normally the case if the former is an up quark and the later a down quark. Note that decays into the light quarks are dominant if channels involving the top quark are kinematically closed. Partial widths as a function of m are depicted in Fig. 1.
Note also that, had we assumed a CP violating trilinear term in the potential, ∼ κ H 2 , this term would induce a VEV for the triplet, v ∼ κv 2 /m 2 . The triplet could also  decay into the Higgs degrees of freedom, the corresponding width scaling as ∼ (v /v 2 ) 2 m 3 . Given that v modifies the ρ parameter, it is bounded to be v GeV [12]. Therefore, the corresponding decay would still be subdominant with respect to those suppressed by v/ f even for f ∼ 10 TeV.
This setup can be easily accommodated in a SUSY framework. The minimal model consists of the MSSM extended with a supermultiplet with quantum numbers (1, 3) 0 ; see Ref. [13]. The most general and renormalizable superpotential is the MSSM superpotential extended by with H 1 and H 2 the two doublet superfields. Likewise, the extra soft-breaking Lagrangian reads In the limit λ → 0, λμ → finite, the fermionic partner in decouples and the model is approximately inert with In the previous expression, κ stands for the trilinear coupling in κ H H 2 . (With a slight abuse of notation, we are denoting here by H 2 the scalar component of this superfield.) After integrating H 2 out, we get with m 2 the mass of H 2 and y u(d) 2 its Yukawa couplings. The Lagrangian above can also arise naturally in CHMs, in which both H and are pNGBs originated in the (approximate) symmetry breaking pattern G → H at a scale f at which a new strong sector confines. The smallest realization of this setup relies on SO(7) → G 2 . The global symmetry is only approximate because it is explicitly broken by loops of SM gauge bosons, as well as by linear mixings between the left-and right-handed top quark fields and composite operators O L ,R . The latter transform in representations of SO (7). If O L ∼ 35 and O R ∼ 1, one obtains the Lagrangian [14] L ∼ y u q LH and similarly for other quarks.  The elementary ITM can also get 1/ f corrections provided it is extended with new vector-like quarks with quantum numbers ; and/or with scalar doublets with quantum numbers (1, 2) 1/2 . The effective operator in Eq. 1 is then generated after integrating the heavy modes at the mass scale M ∼ f ; see Fig. 2. This list exhausts the possible tree-level weakly-coupled UV completions of the ITM that fit into our phenomenological framework.
LEP operated at √ s = 209 GeV, excluding masses below ∼ 100 GeV. (This limit is however slightly model dependent; other sources of new physics could weaken it to even ∼ 75 GeV [15].) We will therefore restrict our analysis to the mass range 100 < m < 500 GeV.

The electroweak phase transition
The new scalar potential modifies the EWPT, which in the SM is a cross over. In the region of the parameter space where it is first order and strong, gravitational waves can be produced via nucleation and eventual collision of bubbles of symmetrybreaking vacuum. In order to explore this phenomenology, we study the evolution of the one-loop effective potential at finite temperature: C is a constant fixed so that V vanishes at the origin of the field space. V tree stands for the tree-level potential. V CW is the one-loop correction at zero temperature in MS and Landau gauge, namely where i runs over all bosons (+) and fermions (−). The factor n i denotes de number of degrees of freedom of the field i, while c i is 5/6 for gauge bosons and 3/2 otherwise. The field-dependent masses squared of the spectator fields are In addition, we have two more field dependent masses squared, m 2 1 and m 2 2 , given by the eigenvalues of the mixing matrix Finally, the finite temperature corrections read As input parameters, we take m , λ H and λ . The remaining three parameters in the tree level potential are numerically obtained after requiring V tree + V CW to have a extreme at h = v, φ 0 = 0, at which the physical Higgs and masses are m h ∼ 125 GeV and m , respectively. In other words: At tree level, (v, 0) is guaranteed to be an extreme provided λ , λ H > 0 and μ 2 > −1/2v 2 λ H . A comparison between tree and loop level values of μ H , μ and λ H in a set of benchmark inputs can be seen in Table 1. At high temperatures, the EW symmetry is restored. For certain values of the model parameters, the transition between h = 0 → h = v n is not smooth as in the SM, but rather first order. One example is given in Fig. 3. In this case, the EWPT proceeds in two steps. An example of a first-order EWPT in one step is shown in Fig. 4. We will denote by T n the nucleation temperature, namely the temperature at which the Higgs first order phase transition takes place. This is determined by the condition S 3 /T n ∼ 100, where S 3 stands for the action of the thermal transition between vacua [16,17].  Fig. 3 Evolution of the VEV with the temperature for the parameter space point m = 120 GeV, λ h = 1.1, λ = 1.1. At high temperatures, the EW symmetry is restored. It is spontaneously broken to ( h , φ 0 ) ∼ (0, 10) GeV at T ∼ 160 GeV, evolving until T n ∼ 85 GeV at which the step (0, 120) → (220, 0) GeV takes place. This latter transition is clearly strong; h /T n > 1. The shape of the potential at T n is also shown Note that at T > 0, the triplet squared term reads ∼ μ 2 + T 2 , which can not be negative for any value of μ 2 > 0. Therefore, the 2-step EWPT can only occur if the triplet minimum is present at T = 0. Moreover, for a fixed m , there is a minimum λ H below which μ 2 is not negative. Likewise, there is a maximum value of the coupling above which the potential at the triplet minimum, V (0, φ 0 ) ∼ −|μ | 4 /λ , is deeper than the Higgs one, the theory being therefore unstable. Altogether, they explain the bounded shape of the figure above.
It is also well known that strong first order phase transition, resulting from non-standard Higgs sectors, produce gravitational waves . They are roughly characterized by the normalized latent heat of the phase transition (with (T n ) the latent heat at T n ), and by the inverse duration time of the phase transition, We computed these quantities using CosmoTransitions [48]. (S 3 /T ) is estimated finding the two values of T for which S 3 /T = 100, 200 GeV, respectively. (Therefore, (S/T ) = 100 GeV.) We warn that, due to the rapid growth of S 3 /T with T , the linear estimation of the derivative can be sensibly overestimated. Given that small values of β/H give rise to stronger gravitational waves, our results are conservative. The region of the parameter space that we estimate it can be tested by the future gravitational wave observatory LISA lies above the dotted green line in Fig. 5. The points in this are lead to α, β within the region "C1" of Ref. [21] for T n = 100 GeV. (The bubble velocity is close to unity in good approximation. Also, we have neglected the effect that sounds waves might be not "long-lasting", what could weaken the gravitational wave signal [49].) Large values of λ H can be also probed in the h → γ γ channel. Indeed, the width of the later in this case reads with α the electromagnetic constant and τ i = 4m 2 i /m 2 h . The last ATLAS+CMS combined measurement of the Higgs decay into photons was provided in Ref. [50], (h → γ γ )/ (h → γ γ ) SM = 1.14 +0. 19 −0.18 . The region in the plane (m , λ H ) that is consequently excluded at the 95 % CL is enclosed by the solid orange line in Fig. 5. The expectation at the HL-LHC is that ratios outside the range 1.0 ± 0.1 will be excluded [51]. The corresponding region is enclosed by the dashed orange line. It is clear that, if departures from the SM prediction on the Higgs to diphoton rate are not observed, only one-step EWPT (i.e. single peak signatures) could be detected by LISA.
Finally, let us very briefly comment on the possibility of EW baryogenesis [52][53][54][55]. In our scenario, CP is violated spontaneously during the second transition in the two-step case, when both h and φ 0 change VEV and therefore the top mass acquires a CP violating phase. In related models [56] (see also Refs. [32,57,58]), EW baryogenesis has been shown successful provided c v/ f 0.1, with v the change in VEV during the EWPT. In our case, v can be easily 100 GeV (see Fig. 3) and therefore c v/ f 0.1 for c/ f ∼ 1 TeV −1 .
A small explicit CP violating potential V /T 4 n H/T n ∼ 10 −16 , with the Hubble scale H, is only needed to avoid domain wall problems [56]. In our setup, this can be triggered by a small CP-violating term in the potential, ∼ κ H 2 . At leading order in κ, it reflects in the (finitetemperature) potential as V ∼ κ T 3 /(4π) 2 . Avoiding domain walls then implies κ 10 −12 GeV. Let us show that this amount of CP violation evades easily neutron and electron dipole moment (EDM) constraints. Indeed, the neutron EDM arises mainly from the running of the diagram on the right panel of Fig. 6. It gives [59] d n e ∼ 20 MeV For the values of c/ f and κ stated previously we obtain |d n | ∼ 10 −38 e cm, much smaller than the current 90 % CL bound 2.9× ∼ 10 −26 e cm [60].
An electron EDM will be generated mainly via two-loop diagrams as that depicted in the left panel of Fig. 6. Using the expressions of Ref. [61], we find that the value of the electron EDM in our case reads with We obtain |d e | ∼ 10 −42 e cm, much smaller than the latest measurement by ACME [62], d < 1.1 × 10 −29 e cm.
Regarding c/ f , more stringent bounds could be set at colliders. We dedicate next section to this point.

Collider signatures
The scalar triplet can be produced at pp colliders in a variety of ways; see Fig. 7. The corresponding cross sections at √ s = 8, 13 TeV are given in Fig. 8. For completeness, we also provide numbers for 27 and 100 TeV center-of-mass energy.
The triplet can be singly produced in qq initiated processes. The Yukawa suppression, together with the 1/ f factor, makes the production cross section in this channel very small, though. Still, φ 0 can be singly produced in gluon fusion. In the regime m < 2m t , the most constraining searches are those looking for single production of bb resonances. The most up-to-date such analysis was recently released by CMS; see Ref. [66]. It is based on 35.9 fb −1 of integrated luminosity collected at √ s = 13 TeV. The region of the plane (m , c/ f ) that is excluded by this analysis is enclosed by the solid blue line in Fig. 9. It is expected that they become a factor of √ 3 fb/35.9 ab ∼ 9 stronger at the HL-LHC. The projected bound on the plane is enclosed by the dashed blue line in the same figure. For m > 2m t , φ 0 decays mostly into tt. There are however no resonant searches for invariant masses below 500 GeV, neither at √ s = 8 TeV nor √ s = 13 TeV. Moreover, the scalar triplet can be produced in association with top and bottom quarks, namely pp → φ + tb. For m > m t , the most updated and constraining search is the ATLAS study of Ref. [67], which uses 36.1 fb −1 of LHC data collected at 13 TeV. It combines both the semi-and dileptonic channels. The limits on σ ( pp → tbφ ± ) × B(φ ± → tb) translate into the bounded region delimited by the green solid line in Fig. 9. A naive rescaling with the luminosity enhancement suggests that cross sections a factor of ∼ 0.1 smaller can be tested at the HL-LHC. Translated to the plane (m , c/ f ), the corresponding bound is given by the region enclosed by the dashed line of the same colour.
In addition, for m t > m , the triplet can be also produced in the decay of the top quark. Current searches for tt production with t → φ ± b, φ ± → j j have been carried out in CMS at √ s = 8 TeV with an integrated luminosity of L = 19.7 fb −1 [65]. This latter reference sets an upper bound on this rare top decay of B(t → φ ± b, φ ± → j j) < 1-2 % for m ∼ 100-160 GeV. Using Eq. 5, this constraint translates into the region enclosed by the solid red line in Fig. 9. The projected bound at the HL-LHC is depicted, too.
Finally, irrespectively of the value of c/ f , the scalar triplets can be always pair-produced via EW charged currents (CC), pp → W ±( * ) → φ ± φ 0 , as well as via neutral currents (NC), pp → Z /γ → φ + φ − . (Note that φ 0 does not interact with the Z boson and therefore it can not be pair-produced via NCs.) For m < m t , the new charged and neutral scalars decay mainly into q q and bb, respectively. Searches for pairproduced dijet resonances might therefore be sensitive to this regime. The most constraining such search is the CMS analysis presented in Ref. [68]. At this mass scale, each pair of quarks is very collimated and manifests as a single jet. The experimental analysis uses boosted techniques, including jet grooming to remove QCD radiation. The analysis considers L = 35.9 fb −1 at 13 TeV of c.m.e. The current limits on the total cross section range from ∼ 170 pb (100 GeV) to ∼ 20 pb (170 GeV). Therefore, the parameter space of our model is not constrained. Furthermore, a naive rescaling with the larger luminosity shows that this analysis will not be even constraining at the HL-LHC.
For m t < m < 2m t , the NC process gives rise to the final state tb, tb. The latest analysis exploring this channel for masses below 500 GeV was performed by CMS at √ s = 8 we have computed the corresponding ratio in MadGraph, by using the Higgs EFT of the ggHFullLoop model. This ratio turns out to be a good approximation of the ratio in the full theory; see Ref. [64] TeV; see Ref.
[69]. Unfortunately, the corresponding limits range from ∼ 2.5 pb (250 GeV) to ∼ 0.5 pb (500 GeV). No region in our parameter space can be even constrained at the HL-LHC. Likewise, the CC gives tb(tb), bb. To the best of our knowledge, there is however no dedicated search for pair produced resonances decaying to these final states.
Being this channel c/ f independent, we perform a signal and background simulation of this process in Sect. 5. For m > 2m t , the NCs still give resonant tb, bt. The CC channel instead results in tt, tb(bt). Once more, no dedicated analysis exists for this final state. (The lack of analyses sen-sitivity to similar final states has been also recently pointed out in Ref. [70] in the context of composite dark sectors.) However, in comparison to this one, the tb(bt), bb analysis is much cleaner. Furthermore, it probes the mass range where the 2-step EWPT, and therefore EW baryogenesis, can occur.

LHC sensitivity
EW pair production of tb(tb), bb resonances occurs naturally in broadly-studied models, such as the 2HDM. Moreover, the  Fig. 9 Current (solid) and future (dashed) bounds on the parameter space (m , c/ f ). The red region to the left comes from LHC searches for rare top decays, t → bφ ± , φ ± → jets; see Ref. [65]. The upper blue region comes from LHC searches for single production of φ 0 → bb via gluon fusion; see Ref. [66]. The green region to the right comes from LHC searches for single production of φ ± in association with a top and a bottom quark; see Ref. [67]. The middle vertical orange band is the region that could be tested using our dedicated analysis cross section is independent of scalar to fermion couplings, provided the former decay promptly. It is therefore surprising that no experimental search has explored this channel at the LHC yet. A plausible explanation is that the majority of these analyses are based on the 2HDM of the MSSM. In that case, one Higgs doublet gives mass to the up fermions, while the second gives mass to the down fermions. The physical charged components then couple to the top and bottom quarks with effective c/ f ∼ (y t cot β + y b tan β)/v, with tan β the ratio of the two doublet VEVs. Therefore, c/ f 2 √ y b y t /v ∼ 1 TeV −1 , the inequality being saturated at tan β = √ y t /y b ∼ 1.3. As it can be seen in Fig. 9, this value is at the reach of bb resonant searches. Therefore, there is a priori no need for further analyses. The situation in our model and in other versions of the 2HDM is very different, though, because the coupling of the new scalars to the SM fermions can be (and typically is) very small. In the composite Higgs scenario, the effective scalarfermion coupling can also be small, i.e. c/ f 1 TeV −1 , according to Eq. 13 for natural values of 0 < γ ≤ 1.
With the aim of filling this gap, we perform a dedicated analysis for the process pp → φ ± φ 0 → tb(tb), bb. We generate the hard processes for signal and background using MadGraph v5 [71] and shower them to a fully hadronised final state using Pythia v8 [72]. No parton-level cuts are applied. The main backgrounds are tt + jets, ttbb, t (t) + 3b and W + 4b.
At the reconstruction level, a lepton is considered isolated if the hadronic energy deposit within a cone of size R = 0.3 is smaller than 10 % of the lepton candidate's p T . Jets are The longitudinal component of the missing neutrino is reconstructed using the requirement m 2 W = ( p l + p ν ) 2 , with m W the mass of the W boson. The neutrino and lepton fourmomenta, p ν and p l , are then added to the jet wich gives a total invariant mass closest to the top quark mass. After this,

The invariant mass of this top is required to be within 50
GeV of the top mass; 4. Three b-tagged jets are to be found among the jets not coming from the leptonic top. 5. We require the reconstructed masses of φ 0 and φ ± to be similar, i.e. |m φ 0 ,rec − m φ ± ,rec | ≤ 50 GeV.   To decide which b-tagged jet is assigned to φ ± and which two are assigned to φ 0 , we compute all possible combinations and choose the one resulting in the minimum difference between m φ 0 ,rec and m φ ± ,rec . The normalized distribution of this former variable in the main background (tt+jets) and in the signal for m = 185 GeV and m = 310 GeV is depicted in Fig. 10. In Fig. 11, we also show the normalized distribution of the p T of the reconstructed φ 0 . However, cutting on this variable is costly in cross section, but would allow to improve the ratio of signal over background further. Finally, 6. We reconstruct the resonances φ 0 and φ ± in the mass window of ±30 and ±40 GeV, respectively. As shown in Fig. 10 the experimental width of the resonances depends on their masses, thus, the central value of the mass window in the (m φ 0 ,rec , m φ ± ,rec ) plane has to be optimised for each m separately.
The cut-flow for the signal and the relevant backgrounds is given in Table 2. The sensitivity estimates are conservative, as the reconstruction relies on fairly inclusive cuts using large mass windows. Further, it is likely that some of the background was counted twice, as b-quarks from the parton shower in tt and from the matrix element in ttbb are both contributing to the total background. Thus, we expect that the sensitivity can be improved further using a combination of multi-variate techniques [73][74][75] and highp T final states [76].
We estimate the sensitivity at the HL-LHC as S = s/ √ s + b, with s and b the number of signal and background events after all cuts, respectively. It ranges from 2.7 to 7.3 for m between 185 and 340 GeV. Thus, at the LHC ( √ s = 13 TeV) we can probe the entire mass interval with 3000 fb −1 . The corresponding region in the plane (m , c/ f ) is therefore a vertical band, the one enclosed by the dashed orange line in Fig. 9.

Conclusions
Natural scalar extensions of the SM must have a low cutoff f preventing large corrections to the scalar masses. However, such models are usually studied neglecting 1/ f terms. Basing on the real triplet extension of the SM, we have highlighted that, if these terms are taken into account, the phenomenology can be drastically different. In particular, the only renor-malizable interaction allowing the new scalars to decay is so suppressed by the measurement of the ρ parameter, that decays mediated by effective operators dominate.
We have studied the reach of current LHC analyses. We have found that, despite being f independent, searches for EW pair-produced charged scalars decaying to third generation quarks are absent. This is particularly surprising given that such signals appear in a plethora of new physics models, including the 2HDM/MSSM. Therefore, we have developed a dedicated analysis to probe the cleanest of these channels: pp → φ ± φ 0 → tb(tb), bb. We have shown that the whole range of masses ∼ 185-340 GeV can be tested at the HL-LHC.
For this analysis, we have neglected new scalar couplings to the leptons. Under the sort of Minimal Flavour Violation [77] assumed after Eq. 1, explicitly reproduced in concrete models as shown in Eq. 13, the only relevant lepton would be the tau. Still, the decay of φ 0 into τ + τ − would involve only a ∼ m 2 τ /(m 2 b N c ) ∼ 5 % of its width; our results being effectively unaffected. If couplings to the leptons are accidentally larger, the scalar could be better seen elsewhere; see Ref. [78].
We also stress that, had he assume flavour-violating couplings c u(d) i j , they would give rise to a plethora of signals in meson decays [79]. They could be also seen in top decays as in the singlet extension of the Higgs sector [80].
On another front, we have studied the reach of the future gravitational wave observatory LISA to the gravitational waves produced in the EWPT for certain region of the parameter space of the model. In particular, we have demonstrated that regions not yet excluded by Higgs to diphoton measurements will be testable. In this region, the EWPT proceeds mainly in one step, and therefore only one signal peak might be expected.
Finally, it is worth mentioning that in concrete models c and λ H related; normally λ H ∝ c. This is evident for example in CHMs, in which the former (latter) is induced by integrating out heavier resonances at tree level (one loop). To exhaust this point, we plot in Fig. 12 the current and future bounds on the plane (m , c) considering all collider searches and gravitational wave signatures for different simple assumptions on the relation λ H = λ H (c).