The presence of extra scalar singlets is a feature of several motivated extensions of the Standard Model, and the mixing of such a singlet with the Higgs boson is allowed to be quite large by current experiments. In this paper we perform a thorough phenomenological study of this possibility. We consider both direct and indirect searches, and we quantify the current constraints as well as the prospects for future hadron and lepton machines – from the forthcoming LHC run up to a futuristic 100 TeV protonproton collider. The direct reaches are obtained extrapolating the current limits with a technique that we discuss and check with various tests. We find a strong complementarity between direct and indirect searches, with the former dominating for lower values of the singlet mass. We also find that the trilinear Higgs coupling can have sizeable deviations from its Standard Model value, a fact for which we provide an analytical understanding. The results are first presented in a general scalar singlet extension of the Standard Model, taking advantage of the very small number of parameters relevant for the phenomenology. Finally, we specify the same analysis to a few most natural models, i.e. the NexttoMinimal Supersymmetric Standard Model, Twin Higgs and Composite Higgs.
Contents
section 1 Introduction
Is the Higgs boson [1, 2] the only light scalar particle, or is it rather the lowhanging fruit of a richer scalar sector of Nature? This question is of paramount importance for the field of particle physics. Discovering another scalar closeby would in fact have deep implications for our understanding of the concept of naturalness[3], challenging frameworks that do not aim at removing the apparent tuning of the weak scale[4].
Among the many conceivable scalar extensions of the Standard Model (SM), arguably the simplest one consists in the addition of a singlet under the SM gauge group which mixes with the already discovered Higgs. This possibility is ubiquitous in many new physics scenarios, most notably in two of the most successful proposals to address the hierarchy problem of the Fermi scale: the NexttoMinimal Supersymmetric Standard Model (NMSSM, see [5] for a review), and Twin Higgs [6]. It is also interesting that, in the context of “finite naturalness” [7], an extra singletlike scalar would also allow to dynamically generate the weak and Dark Matter scales [8]. Moreover, and independently of the hierarchy problem, the same addition to the SM is among the simples possibilities to achieve electroweak baryogengesis [9, 10, 11].
From an experimental point of view, the 8 TeV LHC run has given us a knowledge of the main Higgs properties to a 20% accuracy [12, 13], and no direct evidence for new physics whatsoever. Despite that, scenarios exist that are poorly constrained by present data, an extra scalar singlet mixed with the Higgs being a paradigmatic example of this case.
The above considerations strongly motivate a thorough experimental quest for such a particle. The goal of this paper is to provide a general and solid strategy to perform this hunt both at current and future machines.
Several studies in this direction have been performed in the literature (see e.g. the recent [14, 15, 16, 17]), however they have mostly restricted their attention to some particular direct and/or indirect search, to some particular model, or to some particular machine. In this paper our aim is to be more comprehensive, building a picture that allows to grasp the general phenomenological features, yet in an exhaustive way. In order to do so, we develop an analytical understanding of what the parameters most relevant for the phenomenology are, freeing us from the necessity of numerical scans. In this way we are able to identify the most promising ways to probe this scenario, thus providing key information for current experiments, as well as for planning future ones.
In particular, we will consider several experimental programs besides the second run of the LHC that has just started. Our study will include the high luminosity phase of the LHC itself (HLLHC [18]), a higher energy phase (HELHC), linear (ILC [19] and CLIC [20, 21]) as well as circular (FCCee [22] and CEPC [23]) electronpositron colliders, up to a futuristic 100 TeV protonproton collider (FCChh [24]).
The natural approach to study such an extended Higgs sector relies on two complementary strategies:

Indirect searches: modification of the couplings of the 125 GeV Higgs boson;

Direct searches: production and detection of the heavier scalar.
If on one hand the impact of an improvement in precision on the Higgs signal strengths is quite easy to estimate, on the other hand the study of the reach of direct searches requires generally an amount of numerical analysis. In this paper, however, we rely on a technique introduced in [25] and developed in [26] to extrapolate the reach of future colliders from the existing limits of the first LHC run. This will help to get an intuition of the capability of the new experiments, avoiding the use of detailed collider simulations, and it will allow for an easy comparison between the relevance of precision physics and direct searches.
The paper is organised as follows. In section 2 we introduce the model and present the basic phenomenological consequences, showing also what are the simplifications that can be made. In section 3 we review the extrapolation procedure presented in [25, 26] discussing its virtues and limitations, and we present our projections for the sensitivity of direct searches in the relevant channels at future machines. These results have a general validity that goes beyond the context of the singlet scenario analysed in this paper. We then apply the bounds to the case of generic singlet scalar coupled to the SM. Sections 4 and 5 are dedicated to the NMSSM, and Twin and Composite Higgs, respectively, and there we project present and future bounds into the relevant parameter space of the two models. We conclude in section 6. Technical details about the extrapolation procedure are given in appendix A.
section 2 Phenomenology of a scalar singlet
2.1 A convenient description
We add to the SM particle content a scalar field , singlet under the SM gauge group. The quantum numbers of the real component of , on which we will focus, allow it to mix with the SM Higgs boson, and we dub the mixing angle . Let and be the two neutral CPeven propagating degrees of freedom, being the one observed at the LHC, and GeV [27] and their physical masses. In terms of the gauge eigenstates
(2.1) 
the mass eigenstates read
(2.2) 
where is the vacuum expectation value (vev) of the singlet field, and GeV.
All the and couplings to SM fermions and vectors are fixed, due to the singlet nature of , only by the mixing angle . They can be written in terms of the couplings of a standard Higgs boson of the same mass as
(2.3) 
where we introduce the notation , . Notice that the branching ratios of are not modified with respect to their SM values, since all the couplings are rescaled by the same factor; the only observable deviation from the SM is therefore the reduced production crosssection. Concerning , its production crosssection is simply the one of a standard Higgs boson of mass rescaled by . Its branching ratios are again those of a SM Higgs boson of the same mass below the kinematic threshold , where the width can in principle become sizeable. All these considerations can be modified by the presence of operators of dimension greater than four in a strongly coupled theory, but in the following we will consider only weakly coupled scenarios.
The phenomenology of the two scalars can then be summarised as
(2.4)  
(2.5)  
(2.6) 
where are respectively the signal strengths of and in any SM channel (unless specified), is the signal strength of a SM Higgs boson, and is its production crosssection. The measurement of the signal strengths, and hence of its couplings, puts a model independent constraint on the mixing angle . This in turn limits the magnitude of possible direct signals of , which for large enough decays mainly into a pair of vector bosons, and , and into if the relative branching ratio is large enough. The only relevant fermionic decay mode of is a pair of top quarks, which is however subdominant with respect to the previous ones.
One reads in expressions (2.5) and (2.6) that the decays of both scalars are completely determined by just three parameters:

the mass of the singletlike scalar ,

the mixing angle ,

the branching ratio BR.
We now discuss them in order to understand better and simplify the description of the phenomenology, starting with BR.
For , one knows that decays of into vector bosons are determined by their Goldstone nature, and that they dominate over all the other SM decay modes. This implies the asymptotic relation
(2.7) 
For large enough the number of free parameters drops then from three to two, further simplifying the description of the model. It is natural to choose as one of them.
Concerning the other free parameter, notice that for any given scalar potential, the mixing angle between the two CPeven scalars can be expressed in the form
(2.8) 
where is the (1,1) element of the mass matrix of the scalar system before diagonalisation, in the basis. Notice that one has three possibilities to obtain a relation between , the physical masses and the parameters of the Lagrangian: one possibility for each independent entry of the mass matrix . We find it convenient to choose , because it is always of the order of the electroweak scale , despite its explicit form as a function of the Lagrangian parameters depends on the model under consideration. While this statement will be made more explicit in the discussion of the potential in section 2.2, it can be understood qualitatively by noting that, in general, it applies to all cases where the only field responsible for EW symmetry breaking (EWSB) is the one associated with that single entry of .^{1}^{1}1The only dimensionful parameter of the SM potential is a trivial example of the above property: its size is related to by an order one factor, the Higgs quartic coupling. This property holds in particular for the doubletsinglet system under consideration, since the scalar field does not play any role in EWSB. We conclude that the use of as a free parameter to describe the phenomenology, instead of , is a convenient choice.
2.2 Potential and couplings among the Higgses
The discussion performed so far allows to treat the decay only in an approximate manner. Moreover, it has nothing to say about possible deviations of the trilinear Higgs coupling from its SM value. It is known that such deviations are allowed to be quite large, at least in some specific models containing singlets[28]. It would then be interesting to compare them with the other indirect signals, possibly in a most general way.
With this in mind, we write down the most general form of a renormalisable scalar potential involving a doublet and a singlet as
(2.9) 
By deriving the potential with respect to , one sees that the condition for the singlet vev requires . This further implies , as one can infer from (2.9). This is not the case we are interested in, and hence we assume in the rest of our discussion. Since spontaneously breaks the symmetry , we assume in general that also and are different from zero. A non zero positive is then a necessary condition for the potential to be bounded from below (neglecting runnings at higher energies that could change this limitation, see e.g. [8]).
We can then substitute a subset of the seven parameters of the potential in (2.9), with the vevs and and with the physical masses and . We start by removing the mass parameters and , in favour of the vevs. In this way one realises that the following relation holds
(2.10) 
substantiating the discussion of the previous section. We then proceed to exchange the last two dimensionful parameters of (2.9), and , in favour of the physical masses and . The most general addition of a scalar singlet to the SM can then be described in terms of five free parameters:
(2.11) 
The mixing angle , which determines most of the phenomenology of the model, is a function of only two of them, see (2.8).
The couplings among the scalars, on the other hand, explicitly depend also on the other three parameters. They read
(2.12)  
(2.13) 
where is as in (2.8) and is the SM value of the triple Higgs coupling. One would like to know how , and impact the phenomenology, compared to and : if all the parameters of (2.11) were equally important, as one could expect in principle from (2.12) and (2.13), the description of this part of the phenomenology would be somehow obscured. We show here that this is not the case, and that the impact of actually dominates over the one of and in most of the relevant parameter space. We do so by performing a large expansion of the interesting phenomenological quantities, BR and , which reads
(2.14)  
(2.15)  
Here we are assuming that and are of order of the weak scale, and we collectively denote by orders of , , and . We do not keep track of the explicit dependence on the couplings in the higher order terms, since we are considering only weakly coupled scenarios where these couplings are never much greater than one. A few remarks are in order:

both quantities have the correct asymptotic behaviour if the singlet is decoupled;

at first order in , only determines BR and (besides and ), and at second order only enters the expressions. We have thus obtained an analytical understanding of which extra parameters will impact more on the phenomenology;
In the following, we will always use the complete expressions (2.12), (2.13) in our phenomenological studies. We notice that, choosing a very small value for , the width of the singletlike scalar becomes larger than its mass, and the triple Higgs coupling grows as , thus making the model effectively strongly coupled. Furthermore, also grows for small , introducing a hierarchy among the dimensionful parameters in the potential.
We conclude this section observing that the minimum is always the lowest one, for a given choice of the parameters. We verified this in a range that widely covers our phenomenological study. For a more detailed study of the theoretical constraints on a singletHiggs potential, see [29].
2.3 Indirect signals: triple Higgs coupling vs. signal strengths
We have now developed enough tools to study the impact of indirect searches in the scenario under consideration. Before exposing quantitative predictions, it is important to know the present and future experimental situation. We thus begin this section with an overview of the current exclusions, as well as the expectations for future sensitivities in precision Higgs measurements.
Current and future sensitivities to Higgs couplings.
A global fit to all the Higgs signal strengths as measured at the 8 TeV LHC constrains at 95% C.L. [30, 31]. Concerning future projections, we rely on the comprehensive summary performed at Snowmass in 2013 [32], and we refer the interested reader to it for more details.^{2}^{2}2More recent results, when available (see e.g. [13, 33]), substantially agree with the numbers given here.
1 reach in  

LHC8  0.2  – 
LHC14  0.080.12  – 
HLLHC  48  0.5 
HELHC  –  0.2 
FCChh  –  0.08 
ILC  2  0.210.83 
ILCup  4  0.130.46 
CLIC  23  0.10.21 
CEPC  2  – 
FCCee  1  – 
We report in table 1 the expected precisions in coming from the measurement of a single coupling; for comparison we provide also the current best precision in a single channel from the LHC8 data. In the case of the mixing with a singlet, these numbers provide at least a good indication of the possible outcome of a global fit. Notice that this is not the case for more complex Higgs sectors, like when the Higgs mixes with another doublet (e.g. MSSM), where a global fit is usually much more constraining (see e.g. [28]). The hadronic collider configurations that we consider are LHC at 14 TeV with 300 fb (LHC14) and 3000 fb (HLLHC) of integrated luminosity. Assuming one will measure SM central values of the Higgs signal strengths, the expected precision will be in the ballpark of 10% and 5% respectively. The ranges reported in table 1 come from different assumptions on the future theory errors [33], and from different systematics scenarios [34]. Future projections for higher energy hadronic colliders are currently unavailable, needless to say they would constitute an essential ingredient in the present analysis. For leptonic machines, we report the expected precisions for the ILC and its luminosity upgrade (ILCup), for CLIC and for circular colliders such as FCCee or CEPC. For CLIC we report two numbers, the first assuming a data taking up to 1.4 TeV with 1.5 fb, the second adding 2 fb at 3 TeV. Assuming SM central values, could be constrained from the percent down to the permille level, FCCee being the most sensitive collider in this respect. Concerning the trilinear Higgs coupling, for definiteness we stick again to the review [32], and refer the interested reader to it for details. The expected precisions achievable at future machines are also reported in table 1. The measurement of the trilinear Higgs coupling benefits from higher energy collisions, thus the projection for FCCee is not reported. The benchmarks we consider for hadronic () colliders are a 33 TeV one (HELHC) and a 100 TeV one (FCChh), both with 3 ab of integrated luminosity. The two numbers given for ILC and ILCup correspond to the ILC500 and ILC1000 cases, those given for CLIC have the same meaning as for the signal strengths, and assume electron beam polarisation (for unpolarised beams, the worst expected precision would increase to 0.28). Again we warn the reader that such numbers are to be taken with a grain of salt, possibly even more than in the case of signal strengths [32]. In any case it looks very challenging to push the sensitivity in below the 10% level.
Results and interplay of different measurements.
The predictions for and are shown in figures 1 and 2. The shaded pink region is excluded by a combined fit to Higgs coupling measurements [30, 31], while continuous pink lines display some values of . Given that the estimates for future sensitivities are only indicative, and that those for some machines (like FCChh) are not yet available, we prefer not to present lines corresponding to specific future colliders. This will make it easier to visualise the reach due both to new precision assessments, and to new data. With this in mind and to avoid clutter, in the lefthand plots we focus on a region of parameters interesting for LHC and its upgrades, and in the righthand ones we zoom into an area more interesting for future leptonic machines. One can see that reaching a permille precision on the Higgs signal strengths would constitute an impressive probe of this scenario. This will become even clearer when compared with direct searches in section 3.3.
We show the predictions for as dashed black lines. As discussed in the previous section, this quantity is not univocally fixed in terms of and , so we show it for the benchmark values GeV and GeV. One can see that the allowed parameter space can accommodate very large deviations from one in , with both signs and up to an order of magnitude. From (2.15) one infers that the large size of such deviations is quite general, in other words it is not typical of the particular values chosen for , at least until does not become too small. We explicitly verified this by exploring other values of not too far from . We further check numerically that variations of both and have a minor impact in the size of , confirming the insight from our analytical formula (2.15) (the plots being for ). Finally, another check we perform is that, in the whole interval of parameters that we consider, one never has . This would be the case e.g. if we have chosen a value of closer to zero.
An important message of figures 1 and 2 and table 1 is that, at each stage of the future experimental program (except for FCCee and CEPC), the trilinear Higgs coupling and the Higgs signal strengths will be both crucially important as probes of the model. Moreover, under some circumstances, could even be the only observable where deviations will show up.^{3}^{3}3This can be expected in a bottomup perspective from the EFT analysis of [35]. Perhaps more important than that, in the currently allowed parameter space one could have deviations which are so large, to induce an increase of an order of magnitude (or more) of double Higgs production (see e.g. [36]), likely already in the reach of current machines.
The presence of so large deviations in the trilinear Higgs coupling is a rather peculiar property: the opposite is in fact true in many other models (see e.g. [37]), for example when the Higgs mixes with another doublet. Our study shows an explicit and rather general model where such large deviations are actually typical. It then motivates a further effort in the community to measure the Higgs selfcoupling, already with the available data, and to assess the relative sensitivity of future machines.
Electroweak precision tests.
In the model under consideration the only contribution to the electroweak precision tests (EWPT) is due to the IRlogarithms in the and parameters, that can be estimated as [38]
(2.16) 
where is the finestructure constant and and are the sine and cosine of the Weinberg angle. and are proportional to the reduction of the couplings to the weak bosons, i.e. to the mixing angle. They also vanish when , being the mixing unphysical in this limit. EWPT do not have a strong impact on the picture under consideration, unless does not go to zero when is decoupled. For this to be possible, one needs to be proportional to (see (2.8)): while this is not the case in Supersymmetry, Twin and Composite Higgs models realise such a behaviour. We will then return on the EWPT in section 5.2.
section 3 Direct searches
Searches for scalar resonances decaying into a pair of SM particles have been performed by the ATLAS and CMS experiments in a large variety of channels. Being based on the full dataset of the first LHC run, many of these searches are becoming sensitive to crosssections that are significantly smaller than the ones of a standard Higgs boson. Therefore they already set relevant constraints for scalar singlets. With the increasing centerofmass energy and luminosity, all these bounds will strengthen further. In this section we discuss the searches that are relevant for our purposes, and present extrapolations of their limits to future colliders.
3.1 Extrapolation of the bounds
In order to get a projection of the mass and crosssection reach of future collider experiments, we extrapolate the expected limits from the 8 TeV LHC relying on a parton luminosity rescaling, as introduced in [25] and developed in [26]. This method, which we are going to recall here, is subject to a number of assumptions that are potentially not satisfied in a generic search. It nevertheless constitutes a quick and efficient tool to get a reasonable estimate of the masses and crosssections that could be probed, avoiding to perform a full MonteCarlo simulation of the experimental analysis.
The main underlying idea is that the number of signal events that can be excluded with a certain confidence level at a given collider, for a given mass of the resonance, depends exclusively on the corresponding number of background events . In other words, starting from the LHC bounds, it is possible to obtain a projection for the exclusions at a different collider, as a function of the resonance mass, studying the scaling of the background with the collider energy and luminosity.
Extrapolation procedure.
Concretely, one proceeds as follows [26]. Take a value of the mass, , and a corresponding value of the crosssection times branching ratio, , that can be excluded at 95% C.L. at a collider with centerofmass energy and integrated luminosity . It is possible to find a new value of the mass, , for which the same number of background events is produced at a different collider with energy and luminosity , by solving
(3.1) 
The number of signal events that can be excluded for the new mass with the new collider will then simply be equal to the one of the original search, .
If the signal acceptances and efficiencies do not change significantly between the two experiments – an assumption which we shall discuss further below – one can finally get the new value of the excluded crosssection dividing the number of signal events by the luminosity , and therefore, imposing (3.1),
(3.2) 
The extrapolation of the backgrounds to higher energies and luminosities, which is needed to determine the mass from (3.1), is performed in the following way. The number of background events of a certain process can be written as
(3.3) 
where is the parton luminosity for a channel involving two initial partons with a centerofmass energy , and is the partonic crosssection for that process. The sum is performed over all the partonic reactions that contribute to the background. An explicit expression for the parton luminosities in terms of the parton distribution functions is given in (A.1) in the appendix.
Two simplifications can now be made. We assume the resonance to be narrow, hence the integral in (3.3) is performed over a small energy range close to the resonance peak, with . The integrand can then be taken approximately constant over the integration range. If furthermore one performs the searches in a window of constant relative width , the equality (3.1) simplifies to
(3.4) 
At energies much higher than the SM thresholds, the crosssections decrease as , and one can finally write
(3.5) 
where the coefficients determine the composition of the background. While the behaviour of the crosssections is clearly a good approximation for highmass resonances, as was the case in [26], one may wonder whether it can be used also in our context, where one is interested in scalar particles that are not too much heavier than the Higgs boson. We verified that (3.5) is actually valid down to a few hundreds of GeV, depending on the process. In any case, we do not use the experimental constraints below the region of validity for our extrapolations, and we do not provide extrapolated exclusions for masses below that limit. More details on this issue are given in appendix A.
Since the extrapolated values of obtained from (3.5) are higher than the original for a collider with increased energy or luminosity, the range of masses where one gets a projection of the bounds will be shifted to higher values. Often, however, the physically more interesting region – at or below the TeV scale– lies outside that range, especially in the case of the colliders with highest luminosity or energy. In order to get an extrapolated exclusion in the full region of interest, in [26] the following method was used. One notices that, for a given energy of the new collider, the values of the masses where an exclusion is provided grows with the luminosity . This follows immediately from (3.5) and from the fact that the parton luminosities are a decreasing function of . Then, one can easily extend the extrapolation of the bounds to lower masses, decreasing the luminosity of the collider from to ; lowering the luminosity enough, one will eventually be able to get an exclusion for arbitrary masses down to the initial values. However, since the excluded crosssection increases with the inverse of , the bounds for low masses will be much less powerful, and can be interpreted only as upper bounds on the actual reach of a collider. This is what has been done in [26].
In order to get more realistic projections for low masses, in our analysis we rescale the crosssection limits, for each value of the mass, by the square root of the ratio between the nominal of the new collider and the luminosity that was actually used for the extrapolation, . This is a quite accurate approximation of the scaling of the excluded crosssection with the luminosity when the systematic errors are small (see (3.7)).
To summarise, starting from an LHC exclusion , the corresponding projection for a different collider, for a given mass , reads
(3.6) 
where is determined from (3.5), substituting with .
Discussion of the main assumptions.
Let us finally make a few comments on our prime assumption that the exclusion is driven by the background alone, which is not necessarily very accurate. First of all, the composition of the background must not change dramatically between the two experiments: the background that is dominant at the LHC must still be the leading one at higher energies. Next, the assumption is justified only if the analysis is done as a counting experiment of signal vs. background events, assuming the resonance to be narrow. If effects due to the shape of the resonance peak are taken into account, the results depend on signal and background kinematical distributions in a nontrivial way.^{4}^{4}4This is however what is done by the experimental collaborations at the LHC; we assume that the results do not deviate much from those obtained with a cutandcount analysis. We always work in the narrow width approximation. The significance of a signal over a given background then is
(3.7) 
where and are the systematical errors on the background and signal, respectively. In the hypothesis that the background dominates over the signal, which is always verified in our case, this last error can safely be neglected. However, in order to have the same significance at different collider experiments when is fixed, one must assume that the systematic error , when it is not negligible, remains approximately constant.
A further issue is represented by the fact that the values of and in (3.7) in general depend on the acceptances and efficiencies, due to the experimental analysis and detector effects. They could vary for different experiments, especially when extrapolating over a wide range of energies – e.g. from 8 TeV to 100 TeV. We will nevertheless assume, following [26], all the acceptances and efficiencies to be independent of the collider energy and the resonance mass. This assumption is another possible limitation of the method, but it is quite reasonable for the second stage of the LHC. For higher energy colliders, our results should be viewed as an estimate of what limits could be achieved assuming an experimental setup similar to the present one, and at the same time as a goal for the future detectors.
3.2 Current and future constraints
Given that we are interested in models with a scalar singlet, the direct searches relevant to us are those for a heavy scalar decaying into , and final states. The summary of all the searches that we have considered is provided in table 2. These are only a partial list of all the searches for extra scalars carried out by LHC experiments, but they are the most sensitive ones for scenarios where the extra scalar is a singlet.
Please note that all the extrapolations that we present here are valid for a generic spin0 resonance decaying into , and , and have therefore a general validity beyond the singlet scenario that we consider in this paper.
The Zz and Ww channels.
Many different searches for scalar resonances decaying into and in fully leptonic and semileptonic final states have been presented both by ATLAS and CMS and are summarised in table 2. The strongest individual constraint comes always from the and channels, which dominate respectively below and above about 500 GeV. The and the other channels give a bound on the crosssection which is always weaker by a factor of a few. However, a recent analysis by CMS [42] shows that the inclusion of all these final states is able to significantly improve the limit over all the mass range.^{5}^{5}5In the same paper, the results are explicitly interpreted in the case of a scalar singlet extension of the SM. We will use only this global combination for our analysis in the divector channel.
The background in the channel, which dominates the exclusion in the region of interest for our extrapolation, comes to a large extent from nonresonant production, and is mainly produced through the partonic reaction [40]. We therefore extrapolate this background using the parton luminosities, using everywhere the MSTW2008 parton distribution functions [51].
The experimental exclusion is provided in the mass range from 145 GeV up to 1 TeV, but we consider only masses above 200 GeV because below that value the SM thresholds influence the shape of the background and the extrapolation procedure breaks down. We explicitly checked that above 200 GeV the excluded crosssection is driven by the parton luminosity to a very good approximation (see appendix A).
In figure 3 we show our results for the projected limits on . The lefthand panel shows the projections for the LHC at 13 TeV with an integrated luminosity of 100 fb, and at 14 TeV with a luminosity of 300 fb and 3000 fb, together with the current constraint from CMS. The righthand plot shows the extrapolations for the luminosity and energy upgrades of the LHC, respectively with a centerofmass energy of 14 TeV and 33 TeV and both with a luminosity of 3 ab, and for the 100 TeV FCChh collider, also considering a luminosity of 3 ab. The reach of FCChh extends up to masses of about 11 TeV.
We have checked that our limits are in very reasonable agreement with those obtained in the Snowmass study [52], where a full collider simulation of the signal and the background has been carried out, and with the ones for the HLLHC presented by CMS in [53].
One sees that the next phase of the LHC will improve the present bounds by more than a factor of 2 – with some improvement to be expected already in the first 13 TeV run – while another factor of 3 will be achieved by the high luminosity phase. In terms of the absolute limit on the crosssection, the highenergy colliders HELHC and FCChh will give weaker constraints, but one has to keep in mind that the crosssections for a given mass will be much larger at such colliders.
The hh channels.
Searches for resonant production have been performed by both the ATLAS and CMS collaborations. The two most sensitive channels at present are [47, 48], which gives the strongest exclusions for masses GeV, and [49, 50], which is more constraining at low masses. Other decay channels, e.g. involving leptons, could also become relevant for the next LHC run and future experiments [54, 55, 56].
The background in the channel is dominated by and multijet events [48]. In terms of partonic processes, both and contribute, but the gluonfusion background is expected to dominate at high energies, hence we will use the parton luminosities for the extrapolation. As in the previous case, we verified that the inclusion of the other parton luminosities would change the result only marginally. We have checked explicitly that the extrapolation procedure is robust above about GeV, where the threshold starts to influence the background (see appendix A), and therefore we restrict the exclusion to masses above this value.
In figure 4 we show our results for the projected limits on from the channel, extrapolating the CMS exclusions [48]. The energy and luminosity benchmarks and the notation are the same as in figure 3. The improvements from LHC8 to the future machines is roughly the same as for the channel.
Concerning the channel, its extrapolation is trustworthy only above 400 GeV, and in this region the exclusion is always more powerful. For this reason we do not consider this channel for our future projections. We however include in our analysis the present exclusion which is relevant for low masses. We use also here the CMS data [50], since ATLAS sees a slight excess in this channel [49].
3.3 Results and comparison of bounds
We now interpret the bounds presented in the previous section in the context of a generic scalar singlet extension of the SM.
The total production crosssection times branching ratio for is given by (2.5) and (2.6). When dealing with current bounds we always use the SM Higgs total crosssection and branching ratios given by the Higgs Cross Section Working Group [57], which are provided for masses up to 1 TeV. On the other hand, since the extrapolated limits extend up to much higher energies, where we cannot rely on the results of [57], we use the LO quantities in our projections. In particular, we compute the gluonfusion crosssection with the program higlu [58, 59], and we use madgraph [60] to calculate the vector boson fusion and VH associated production crosssections. The sum of these three components constitutes a good approximation of the total Higgs production crosssection. Concerning the branching ratios, we use the results of [57] for masses below 250 GeV, and we match them to the LO (plus NLO QCD) results above that scale, where the two quantities are practically equal.
In figure 5 and figure 6 we show the present and future exclusions in the same plane of figure 1, for two different values of the singlet vev . We fixed , however, as discussed above, the impact of their precise value on the results is very weak. The lefthanded panels of both figures show the projections for the present and future stages of the LHC, up to the luminosity upgrade, while in the righthanded panels we show the range relevant for future highenergy colliders – HELHC and the 100 TeV FCChh. The constraints from the channel are drawn in red, and the ones from the two and channels in purple and magenta, respectively. As discussed above, we do not show future projections for the channel, since they are never relevant in the region where we perform the extrapolations.
One sees immediately that, depending on the value of , the relative importance of the and searches can be very different. In particular, we show the cases of GeV and GeV as benchmarks of a large enhancement and suppression of , respectively. In the first case the searches in the channel dominate both at the 8 TeV and 14 TeV LHC, but for the larger masses that can be reached at future colliders their reach becomes comparable with the one from the channel. In the second case of figure 6, on the other hand, is so strongly suppressed that searches in will barely be sensitive to this scenario even at FCChh, while always gives a much stronger exclusion.
The message that emerges from these figures is therefore pretty clear: at a highenergy collider, be it HELHC or FCChh, the searches will always put the strongest bound on the production of a singletlike scalar. On the contrary, at the various phases of the LHC the process can be the most powerful probe, given some particular choice of the parameters. This is easily understood in terms of the fact that the enhancement (or suppression) of the branching ratio is particularly effective at low masses, where the deviations from the limit (2.7) can be larger (see (2.14)), but in the region of parameter space relevant for future hadron colliders (high masses and small mixing angle) it is not possible to modify more than a certain amount. Despite this strong dependence on the choice of parameters of the bounds in the and channels, the combined reach of the two exclusions is very similar in all cases: a suppression in will enhance the branching ratio into vectors, and vice versa, and the constraints on the crosssection in the two channels are of comparable order of magnitude.
We further point out that the enhancement of is correlated with the modification of the triple Higgs coupling, as one can see from (2.14) and (2.15), and comparing figure 2 with figure 5. This implies in particular that, if a resonance in the channel will be discovered in the future, in the singlet scenario under consideration one expects to see an enhancement also in nonresonant double Higgs production.
In figure 7 we show a comparison of all the direct searches with the indirect constraints from the measurement of Higgs signal strengths. Here we fix to the asymptotic value of 0.25, and the direct exclusion is always dominated by . As in figures 1 and 2, we draw contours of constant instead of explicit indirect exclusions for each future collider, since the projections for the precisions that can be attained on the Higgs couplings are still quite uncertain. One can however easily read the numbers off table 1. Notice that in the righthanded plot the range of the vertical axis is changed in order to better show the region of small mixing angle most relevant for future colliders.
The direct searches are always more powerful than the Higgs signal strength measurements for low values of the resonance mass, while the opposite is true for high masses. While this is a general feature of direct and indirect searches, such an interplay was a priori not guaranteed in the physically interesting range of masses, and for the actual sensitivities of the future experiments. At the LHC, the point where the two searches become comparable varies from GeV for LHC8 to slightly more than 1 TeV at HLLHC. Also a linear collider, with a sensitivity to deviations in the Higgs couplings in the 1% ballpark, would become more sensitive than the HELHC upgrade above about 1 TeV (500 GeV for HLLHC). The same mass of a TeV sets also the limit between FCCee and FCChh.
section 4 Non minimal supersymmetry
We can now apply the previous general discussions to some concrete models that are particularly motivated scenarios for physics beyond the SM. In this section we consider the NexttoMinimal Supersymmetric Standard Model (NMSSM), which is the minimal deformation of the MSSM that includes a scalar singlet in the spectrum. A chiral superfield , singlet under the SM gauge groups, is coupled to the two Higgs doublets and of the MSSM through the superpotential
(4.1) 
where is a generic polynomial up to third order in .
This deformation has important consequences for the naturalness of the theory, given the measured value of the Higgs mass. Indeed, in the MSSM a 125 GeV Higgs requires the largest contribution to to originate from SUSYbreaking effects, , where is the radiative correction mainly due to the topstop sector (see e.g. [61]). This in turn requires either top squarks in the multiTeV range, or large trilinear terms, both cases implying a large tuning of the electroweak scale. The ultimate reason for this is that the supersymmetric couplings in the MSSM are controlled by the weak gauge couplings. The situation is different in the NMSSM, where a supersymmetric contribution to the Higgs mass originates at treelevel from the Yukawalike term ,
(4.2) 
where is the ratio of the vev’s of the two Higgs doublets . As manifest from (4.2), the most natural scenario occurs when the coupling is of order one, so that the largest contribution to has a supersymmetric origin. We here briefly sketch the major features of such a picture:

It gives a large supersymmetric contribution to the Higgs mass (with a different dependence on than the Dterms) that is sufficient to achieve 125 GeV at treelevel.

It modifies the dependence of the electroweak vacuum expectation value on the Lagrangian parameters, allowing for parametrically larger soft masses of stops and gluinos – by a factor of – than in the MSSM, for a given amount of tuning [62, 63, 64]. While a large softens the tuning in , it increases the tuning in the Higgs mass because it overshoots its measured value.

The coupling is not asymptotically free and if the theory undergoes a strong coupling regime at a scale lower than the GUT scale [65]. We will consider both regimes in the following, as they both have been deeply studied in the literature.
The NMSSM with a sizeable – which is the leading parameter also for phenomenology – opens then the interesting possibility that the lightest new particles be the extra scalar bosons of its extended Higgs sector (possibly with the exception of the LSP). These new physical degrees of freedom consist in three neutral scalars, two pseudoscalars, and one charged Higgs boson.
The phenomenology of the NMSSM is controlled by several parameters, so that its study is often performed via scans on some motivated ranges (see e.g. [66]), which need the specification of a particular version of the NMSSM potential, i.e. the function and all the soft terms. Here we adopt a more general approach to have an analytic understanding of the phenomenological properties of the extra Higgs bosons, focussing on the CPeven ones, in the spirit of [28, 67].
The SMlike Higgs field of (2.1) reads . As the present data suggest (see [28, 67, 68]), a particularly interesting limit is the one where the orthogonal Higgs doublet, which does not participate in EWSB, is decoupled.^{6}^{6}6The case of a nondecoupled second doublet is mostly probed via Higgs coupling measurements, and its phenomenology is similar to that of an MSSM Higgs[28]. Moreover, such possibility is further constrained by the presence of a charged Higgs boson. Under the assumption of CP conservation, we can single out the CPeven component from the complex field . In this case the system is equal to the twostate scenario described in section 2, and can be studied along the same lines.
The main parameters are again and , with the former depending on the parameters of the model like
(4.3) 
Similarly to section 2, the mass eigenstates are defined as in (2.2), while the mixing angle between the Higgs and the singlet is now given by
(4.4) 
The phenomenology is mainly set by and but, as discussed in section 2.1, quantities like the trilinear couplings depend on other parameters. For simplicity, in the following we will fix BR to its asymptotic value 1/4, so that there is no need to enter the details of the scalar potential.^{7}^{7}7The coefficients of the potential in (2.9) can be computed in terms of the parameters of the model. This amounts to have an explicit expression for in (4.1) and for the soft term lagrangian . As an example in the scale invariant NMSSM (see [5] for a review), The matching for the coefficients, upon the decoupling of one of the two Higgs doublets, is given by Notice however that all the analysis of section 2.2 applies, in general, also to the present case. In particular, from the point of view of phenomenology is again the leading extra parameter.
4.1 Results and discussion
The quantity now depends on parameters with a defined physical interpretation. Given the expression (4.4) we have chosen to fix and let free to vary. Notice that the value of is almost irrelevant for sizeable . We consider the following two cases:

Perturbative coupling, . In this case we use a value for consistent with perturbative unification at the GUT scale, and we take GeV in order to more easily reproduce the Higgs mass.
Figures 8 and 9 show the present and future constraints on this scenario for both indirect and direct searches. We restrict ourselves to low values of since this is the region most relevant for naturalness [62, 63]. Contrary to the generic Singlet case, the direct searches always dominate at the LHC and its upgrades. In the context of figure 7, this can be interpreted with the observation that, for the physically motivated ranges of NMSSM parameters, cannot be too large, see (4.3). The two benchmarks with different values of have significantly different behaviours: the larger is the smaller is the allowed parameter space. Moreover, while for a larger leptonic colliders will start to be competitive with hadronic ones, for a smaller only a permille level of precision in the Higgs couplings could be competitive with direct searches. This means that the observation of a deviation in the Higgs signal strengths, without any direct indication of an extra scalar, will strongly disfavour this scenario. It is interesting that, for supersymmetry at the LHC, our findings indicate that the strategies to best look for a singletlike Higgs and a doubletlike one look quite different: direct searches in the first case, and mainly indirect ones in the second case (barring the large case [72, 73]).
From the point of view of indirect searches, it is notable that also here large deviations from the SM value are possible in the triple Higgs coupling, as shown in [28] and confirmed in a particular model by [74], and contrary for example to the MSSM case. In general there are no particular differences with respect to the analysis of section 2, given that , and are free parameters also in the NMSSM. However, as we have already discussed, the NMSSM prefers low values of . Accordingly, the deviations from one in are expected to be milder than in the generic case (see figures 1 and 2), but still allowing for observable effects at the LHC.
section 5 Twin & Composite Higgs
The presence of an additional scalar singlet is also typical of another class of models, namely those where the Higgs is a pseudoGoldstone boson. In this type of models the Higgs boson emerges as an approximate Goldstone boson associated to the breaking of some global symmetry . In the case where there is only one physical pseudo Goldstone Higgs, the only other additional scalar singlet is the “radial mode”, i.e. the degree of freedom associated with vacuum expectation value , with , that breaks the global symmetry. Throughout this section we dub the extra physical scalar instead of .
If the above description derives from a weakly coupled model at higher energies, then the radial mode could be relatively light, significantly below . This is for example the case of the Twin Higgs (TH) idea [6], with a weakly coupled description [75, 76, 77, 78, 79], and we now turn to this possibility.
Twin Higgs.
Twin Higgs models solve the little hierarchy problem up to some scale , without requiring the presence of coloured degrees of freedom below such scale.^{8}^{8}8Above , a theory solving the big hierarchy between and the Planck (or GUT) scale is needed. This is contrary to the usual realisations of SUSY and composite Higgs models (CHM), where new coloured bosons or fermions are expected to lie within the LHC reach. In TH models the low energy degrees of freedom that cancel the “quadratic divergences” are total singlets of the SM and hence very difficult to detect at the LHC, while coloured particles are expected to be naturally heavier. In TH models even a null result from the LHC will leave the tuning of the weak scale at the level of (for recent studies see [79, 80, 81, 82, 83]).
In this paper we would like to emphasise that, if the TH has a weakly coupled description, an interesting possibility consists in trying to detect the “radial mode”, the mirror partner of the Higgs. In fact, differently from other mirror particles, it can be singly produced through a mixing with the Higgs. This is complementary to the idea that the only modelindependent probe of such scenarios consists of the Higgs signal strengths. In order to discuss the phenomenology of this picture we here briefly review the basics of the model.

A copy of the SM gauge, matter and scalar content is added, related to the SM by a symmetry. A scalar portal coupling is allowed between the two sectors.

An important assumption is that the scalar potential is accidentally SO(8) invariant,^{9}^{9}9For our discussion also SU(4)/SU(3) is fine. in such a way we can have a spontaneous symmetry breaking SO(8)/SO(7) parametrized by
(5.1) where is the SM Higgs doublet, and is its symmetric partner, singlet under the SM gauge group. The global symmetry SO(8) is manifest when the above expression is written in terms of the real components. From (5.1), and from the fact that the two copies of the SM are gauged, out of the 7 Goldstones only one remains as a physical fluctuation, the other six being “eaten” by the vector fields. In addition there is also a radial mode of mass . At this level one would like to explicitly break both SO(8), to give a mass also to the other degree of freedom, and , in order to get a hierarchy between the two vev’s in the two sectors.

This can be achieved by adding to the above potential other gaugeinvariant terms, so that the full treelevel potential reads
(5.2) where the first term breaks SO(8) but preserves , and the second one breaks both groups. A nice feature of TH models is that the radiative corrections do not introduce any sensitivity to high scales, and the finite contributions are suppressed by the fourth power of SM couplings (as an example [6]).^{10}^{10}10A term can be reabsorbed in a redefinition of and of the coefficients in (5.1). The same is true for quadratic terms containing a different combination of and .
These three ingredients are sufficient to describe the phenomenology. The parameters in (5.2) are related to the vev’s and via
(5.3) 
so that is the symmetry breaking scale of a stronglycoupled scenario, which indeed is obtained in the limit of large , i.e. when the mass of the radial mode is close to the cutoff and it can be consistently integrated out.
From the diagonalisation of the mass matrix we can then fix and trade the residual free parameter for , the mass eigenvalue of the singletlike state. The mass eigenstates and are related to the gauge ones via a rotation by an angle as in (2.2), with
(5.4) 
In the notation of section 2, the parameter is now given by
(5.5) 
A crucial observation is that the only two free parameters are (, ), which completely determine the phenomenology. Like in the generic singlet case of section 2, the mixing angle is also the value of the normalised couplings of and to SM particles. Other important informations are provided by the trilinear couplings of the mass eigenstates
(5.6) 
Composite Higgs models.
The above model only describes two real scalars, so that it closely resembles the case of a linearised minimal composite Higgs model with an SO(5)/SO(4) symmetry breaking [84]. In this context the global symmetry is spontaneously broken by some strong interaction, and the Higgs is again a pseudoGoldsone boson of this breaking. The main phenomenological difference with the previous case is the presence of vector and fermion resonances related with the strong sector, with a mass of a few times the scale .^{11}^{11}11See e.g. [85] and references therein.
There is one single radial mode associated with this symmetry breaking pattern, which forms a quintuplet of SO(5) together with the Higgs doublet ; at energies higher than its mass one can describe this scenario through a linear model. Although this might be an improper description of a stronglyinteracting scenario, where the radial mode is generically expected to have a mass close to the cutoff of the theory , it provides a useful and simple parametrisation [38]. Using the same notation as above, and trading with the pure real singlet , the potential can be written as [86, 68]