Composite Higgs models after Run2

We assess the status of models in which the Higgs is a composite pseudo-Nambu Goldstone boson, in the light of the latest 13 TeV Run 2 Higgs data. Drawing from the extensive Composite Higgs literature, we collect together predictions for the modified couplings of the Higgs, in particular examining the different predictions for $\kappa_V$ and $\kappa_F$. Despite the variety and increasing complexity of models on the market, we point out that many independent models make identical predictions for these couplings. We then look into further corrections induced by tree-level effects such as mass-mixing and singlet VEVs. We then investigate the compatibility of different models with the data, combining the Run 1 and recent Run 2 LHC data. We obtain a robust limit on the scale $f$ of 600 GeV, with stronger limits for different choices of fermion embeddings. We also discuss how a deficit in a Higgs channel could pinpoint the type of Composite Higgs model responsible for it.


INTRODUCTION
Composite Higgs models [1][2][3] offer an elegant solution to the hierarchy problem of Higgs physics. They postulate the existence of a new strongly interacting sector which confines not far above the electroweak scale. In recent years there has been significant interest in a specific class of these models -models in which the Higgs emerges as a pseudo-Nambu Goldstone boson of the strong sector. This sector is taken to be endowed with a global symmetry which is spontaneously broken in the confining phase, protecting the Higgs mass from corrections above the compositeness scale. Although the idea is reasonably straightforward, there are, as with most theories Beyond the Standard Model, many possibilities for its realisation.
Although this plethora of models offers a variety of unique and interesting predictions, those that are most immediately testable are the modifications of the Higgs couplings to the rest of the Standard Model fields. Of particular interest are the values of the coupling modifiers κ V and κ F , as defined in [4].
In this paper we summarise the predictions for these couplings in Composite Higgs (CH) models. We make the case that, despite the diversity of models in the literature, these predictions have very generic structures, and we attempt to provide some intuition for this fact.
We then investigate some simple cases in which treelevel effects can modify these generic structures. These can occur, for instance, in models with extra singlets that get vacuum expectation values (VEVs), or models with an extra SU (2) L doublet that mixes with the Higgs. We point out that to leading order the modifications to κ V and κ F are precisely as one would expect in corresponding models where all the scalars are elementary, plus the usual CH corrections.
Taking the generic structures we have identified, we then perform a χ 2 fit to the data, allowing for the possibility that different fermions couple in different ways. We place bounds on the compositeness scale f , and identify the classes of models that are most constrained.

THE NON-LINEAR COMPOSITE HIGGS
In Composite Higgs models, the Higgs is realised as a pseudo-Nambu Goldstone boson (pNGB) of a broken global symmetry. This symmetry is a symmetry of a new strongly interacting sector, out of which the Higgs emerges as a composite.
Let the global symmetry be denoted G and the subgroup to which it spontaneously breaks be denoted H. Then the Higgs and the other pNGBs (denoted collectively by φ a , one for each broken generator X a ), are parametrised via where f is an energy scale associated with the spontaneous symmetry breaking. U transforms non-linearly under the global symmetry G: where g ∈ G and h ∈ H. By non-linear we mean that the transformation h is field-dependent: h = h(g, φ a ). In cases where the coset G/H is symmetric 1 we are allowed to construct an object (which we will label Σ) whose transformation under G is linear. In all the models considered here [3,[5][6][7][8][9][10][11][12][13][14][15][16][17][18][19][20][21][22], and in the vast majority of models in the literature, G/H will be symmetric. This reduces the task of writing down a low-energy effective theory for the pNGBs to a relatively trivial search for invariant combinations of Σ and the other relevant fields. 1 If T a and X a are the unbroken and broken generators respectively, then the Lie algebra of a symmetric coset obeys the schematic relations

arXiv:1703.10190v2 [hep-ph] 2 Feb 2018
We will assume that the Higgs boson is a doublet under SU (2) L , which, along with U (1) Y , must be embedded as an unbroken subgroup of G. Although data strongly supports the doublet scenario (e.g. see LHC constraints on the ratio of couplings to W and Z bosons [4]), nonlinear models have been studied in which the four scalar fields are actually a singlet and a triplet under SU (2) L [23-26] 2 .

Gauge couplings
The couplings of the Higgs to the gauge bosons come from the kinetic term for Σ, which in the CCWZ prescription [29] is: where D µ = ∂ µ − igA µ , with A µ = A a µ T a for each gauged generator T a . We assume that the Higgs is embedded in a bidoublet (2, 2) of a custodial SO(4) SU (2) L × SU (2) R ∈ H -this is necessary in order to protect the ρ parameter from unwanted corrections [30]. Note that this imposes the non-trivial requirement that H must contain an unbroken factor of SO(4).
Since we are interested in the couplings of the physical Higgs boson to SM fields, we will expand Σ along the direction in which the Higgs will get a VEV, and set all other pNGB fields to zero. The term in (3) will generically 3 lead to a Higgs-gauge coupling of the form: which is valid as a series expansion around H/f . Expanding around the Higgs VEV H → H +h (where h is the physical exictation of the Higgs field) we find the gauge boson masses and couplings: we find Thus Since κ V is defined as g W W h /g SM W W h , we find Since the structure of (3) is generic, so too is this result, at leading order, across all Composite Higgs models.

Fermion couplings
In Composite Higgs models the SM fermions usually couple to the strong sector via the partial compositeness mechanism [8,31,32]. As far as this mechanism pertains to the construction of the low energy effective theory, it involves embedding the SM fermions in representations of the global symmetry G, and then constructing G invariant operators out of these multiplets and Σ. Such an embedding is sometimes called a spurion -the term spurion refers to the 'missing' elements of the multiplet, since after all, the SM particles do not come in full multiplets of the new symmetry G. The incompleteness of these spurious multiplets contributes to the explicit breaking of G and allows the Higgs to acquire a potential via loops of SM fermions.
The appropriate representation in which to embed the SM particles would, in principle, depend on the UV completion of the model. Some attempts towards UV completions of Composite Higgs models have been made (see, for example [5,7,9]), however for the purposes of most model building the choice of representation is a 'free parameter' of the model. There is, however, good cause to restrict the choice of representation into which the SU (2) L quark doublet is embedded. As shown in [33], embedding q L into a bidoublet (2, 2) of the custodial SO(4) SU (2) L × SU (2) R can prevent anomalous contributions to the Z → bb coupling. This restriction forces one to choose representations that contain a bidoublet in their decomposition under the custodial SO(4) subgroup of G.
To treat the EFT in full generality, one should embed q L , t R and b R into different multiplets Ψ q , Ψ t and Ψ b . The kind of representation that the three quarks are embedded into need not be the same. Thus, even for each coset G/H, there are a bewildering number of possibilities. However, for the vast majority of models the form of κ F is actually quite restricted. We tabulate a few examples in Table I.
It might seem strange that so many distinct models lead to so few possibilities for κ F . In fact, when one examines the structure of the allowed terms in the effective Lagrangian, a general pattern emerges: the lowest order coupling of the Higgs to fermions will generally contain either one or two factors of Σ. For example, in the Minimal Composite Higgs Model (MCHM), the coset group is SO(5)/SO(4), and one can define a linearly transforming Σ in the 5 of SO (5), which, expanded along the H direction can be expressed as With q L and t L embedded in the 5, Yukawa couplings come from the SO(5) invariant effective operator leading to a term proportional to sin(H/f ) cos(H/f ). Alternatively one could embed q L into a 10, the t R into a 5 -in this case the Yukawa term originates from an operator like and the interaction is proportional to sin(H/f ) 5 .
In general the structure must be such that the leading term in the trigonometric expansion is H/f . In almost all cases the relevant term will be proportional to either sin(H/f ) or sin(H/f ) cos(H/f ). This argument is certainly not intended to be rigorous -we merely hope to provide some intuition for the fact the non-linear nature of a pNGB Higgs boson leads to repeated structures even across different models and choices of representations 6 .
Following the same procedure as in equation (5), we can expand around the Higgs VEV to find the expression for κ F , defined by yv/m F . A coupling of the form ψψ sin(H/f ) leads to while a coupling of the form ψψ sin(H/f ) cos(H/f ) leads to As we stated above, the representation into which we embed t R and b R might not be the same -in this case it is quite possible (depending on the details of the model) that the top and bottom couplings to the Higgs have different structures. For instance, in the second example above, although the t R is embedded into a 5, the b R might be embedded into a 10. As a result the top coupling would scale with 1 − 1 2 ξ while the bottom coupling would scale with 1 − 3 2 ξ. There are (as always) some interesting exceptions. For example, in [19], with q L in a 5 and t R in a 14, one can derive κ F ≈ 1−3ξ, see also Ref. [36]. In some models (for some examples, see [11,19]) more than one operator can be constructed which contributes to the same Yukawa coupling. The degree to which each operator contributes will, in such cases, be a free parameter and will lead to more complex expressions for κ F . Such models are interesting insofar as they are exceptions -however more minimal scenarios will follow the structure we have outlined above.
No mention has been made so far of the leptonic sector. In theory the lepton Yukawas can also be generated via the partial compositeness mechanism (see for instance [12]). This means that κ τ (for instance) would also receive corrections, and in minimal scenarios would depend on ξ like κ A F or κ B F , as defined in Table I.

TREE-LEVEL EFFECTS
In this section we will briefly look at two interesting scenarios that can lead to tree-level corrections to κ V and κ F from the integrating-out of heavier states. We will describe these corrections as leading to a new effective ξ eff to be compared with the vanilla prediction for ξ.
The first possibility is that in models with an extra singlet pNGB (such as the SU (4)/Sp(4) and SU (5)/SO(5) cosets), the pNGB potential could induce a VEV for the singlet. This can modify κ F and κ V in two ways -firstly a VEV for the singlet η will induce singlet-doublet mixing between η and H. Singlet-doublet mixing (in the elementary case) and its effect on Higgs couplings was studied in detail in [37]. The fact that the H mixes with another scalar means that the couplings will be modified by a factor of cos θ, where θ is the mixing angle between H and η. For small mixing angles: In this and in the following we are assuming that the singlet is heavier than the Higgs and that it makes sense to integrate it out. Generally, in the absense of further tuning, one expects the extra pNGBs to be heavier than the Higgs by a factor of ξ = v 2 /f 2 , since this is the amount by which the mass of the Higgs has to be tuned to satisfy electroweak precision test [38]. Thus, in models with around 10% tuning, values for the extra pNGB masses of around 300 − 500 GeV are not unreasonable.
There could also be effects similar to those studied above, arising from higher-dimensional terms in the nonlinear effective theory. As an example we will look at the SU (4)/Sp(4) model. The gauge boson coupling to the Higgs and η (the equivalent of equation (4)) will be (neglecting hypercharge) As expected, there is no dimension-4 coupling of η to the SU (2) L gauge bosons, but there are higher order terms involving η which could modify the hW W coupling if η gets a VEV. However one should also note that the kinetic term in (3) corrects the Higgs kinetic term: After canonically normalising the Higgs field and expanding around small values of ξ η = v 2 η /f 2 we find that the O(ξ η ) correction to κ V actually cancels. To leading order in ξ, ξ η and θ we have: The correction due to the singlet VEV thus neatly "factorises" into the mass-mixing correction O(θ 2 ) plus the usual compositeness correction O(ξ). We can thus define a ξ eff = ξ + θ 2 , such that κ V ≈ 1 − ξ eff /2. One finds a similar result for κ F . The singlet VEV modifies κ F from ≈ 1 − 3 2 ξ to and in this case our effective ξ eff = ξ + 1 3 θ 2 . In the regime where m η and v are both v, the mixing will be small and will scale approximately as where we have related m η to f via some coupling: m η = g η f . The amount of tuning present in such a model was analysed in [39]. This coset was also investigated in a cosmological setting in [34,40], where the singlet η plays the role of the inflaton. In such a scenario the size of the singlet VEV has important implications for the scale of inflation, and the mass-mixing of the inflaton would be important also for the process of reheating. Moreover, the singlet η and a non-zero value of ξ η could be a key component of a solution to the matter-antimatter asymmetry in the Universe [41].
If the value of ξ eff were the same for all couplings (i.e. the modifications to κ V and κ Fi were the same), then the theory would resemble a CH model without any mixing, only with an apparent rescaling of f . However it is interesting to note that in the above case the inferred values of ξ eff from the measurements of κ V and κ F are different, which would in principle allow us to experimentally distinguish between these two scenarios.
Another possibility is that the spontaneous breaking leads to another pNGB doublet of SU (2) L (a composite two Higgs doublet model). In principle, explicit breaking effects could lead to a mixing between the two doublets. This possibility is discussed in [14,15], and in a different context in [18], in which the two doublets appear from two different spontaneous breakings at different scales.
In this case we will obtain similar results to our expressions above for ξ eff , with a correction from the massmixing at O(θ 2 ) that will be present in the elementary case, and the usual correction at O(ξ) coming from higher dimensional operators (see [42] for a review of the elementary two Higgs doublet model, and [37] for an analysis of the Higgs EFT in such a scenario).
Since we have looked at tree-level corrections to κ V and κ F coming from new states in the composite sector, one should also talk about loop level modifications. In principle loops of scalar, fermionic and vector resonances of the strong sector can modify the Higgs couplings. These will arise from higher dimensional (d ≥ 6) operators in the effective theory, suppressed by factors of f 4−d .

STATUS AFTER RUN 2
In this section we study the impact of Run 1 LHC data on Composite Higgs models, as well as the improvement which results when adding the 13 TeV results recently released by the collaborations. In Table II we summarize the channels considered in the combination of Run 1 and 2 data from ATLAS and CMS, as well as indicate the coupling modifiers that one would obtain in Composite Higgs models, as discussed previously.
The couplings of the Composite Higgs to gluons and photons, κ g and κ γ , are functions of the modifications of the couplings to fermions and gauge bosons, which  II: List of 13 TeV channels considered in the fit, with the corresponding κ modifyiers. Note that the 7+8 TeV Run 1 data was included using the results of the combination of ATLAS and CMS data in Ref. [4].
We then perform a χ 2 fit to the ATLAS and CMS data 7 , with the restriction ξ > 0. 7 When two measurements of the same channel were available, we discarded the worse measurement, or kept both if they were of similar significance. Results from [54,55] were considered but not included in the fit.
The dependence of the χ 2 function with the scale of new physics f is shown in Fig. 1. The green and yellow bands correspond to the one-and two-sigma regions of the fit, and the left and right panels correspond to Run 1 and the combination of Run 1 and Run 2, respectively. Different choices of fermion representations κ A,B F (as shown in Table I) lead to different χ 2 dependences.
The model-independent limit on f improves from 450 GeV (Run 1) to 600 GeV (Run 1+2) at 95% CL, and we see that the most constrained screnario is κ t = κ A F , κ b = κ τ = κ B F . Moreover, one can see that the spread of limits on the scale f due to these fermion choices increases with the addition of more data. This is a signal that the data is increasingly sensitive to these choices, due to better determination of the Higgs couplings to the heavy fermions. To illustrate this point, assume that at some point in the future a deficit in one channel is observed, whereas other channels remain consistent with the SM. For example, assume that the signal strength of the ttH processes was found to be a third of the SM rate, whereas other processes involving the coupling of the Higgs to vector bosons remained consistent with the SM. In this case, certain representations for fermion embeddings of the top and bottom quarks would be preferred by data, see Fig. 2.
These limits on f should be compared with the limits of direct searches for new resonances. One would typically expect a set of new resonances, e.g. new massive W and Z , to appear at some scale related to f , m W = g ρ f , with g ρ O(4π). The value of g ρ is an input to the effective theory, but can be obtained by performing a lattice simulation of the theory and investigating the spectrum of resonances. Its value depends on the specific pattern of breaking as well as the possible electroweak effects. As an indicator of the value of g ρ in these kind of models, we draw attention to the work done in the coset for (κt, κ b ). In this case, the choice κt = κ b = κ B F would be preferred by data. We assume a 20% uncertainty in these channels, except in gg → H → γγ where a 10% accuracy is assumed.
SO(6)/SO(5) [56], and in others scenarios [57], where g ρ was found to be O (10). In this case, a limit on f ∼ 600 GeV, would correspond to a Z and W in the multi-TeV scale, certainly competitive with direct searches for these resonances.
Besides vector resonances, one would expect a tower of fermion resonances, or techni-baryons. Typically, these techni-baryons are heavier than the vector bound states by a factor of N T C , with N T C the number of colours in the new strongly coupled sector [58,59]. Hence, naively one would expect fermion resonances again in the multi-TeV scale. Yet, in most Composite Higgs models the mechanism of electroweak symmetry breaking depends on the existence of light techni-baryons (top partners) with masses of the order of f , contrary to the large-N expectation. This mechanism is being tested by direct searches of heavy partners of the top, with recent Run 2 results already sensitive to the 1.2 TeV region [60], clearly more competitive than the indirect searches in Higgs data if one believed this is the correct mechanism in place. Note, though, that the mass of the top partner is also linked to the amount of fine-tuning in these models. From this point of view the strong limits in top-partners may lead one to consider alternative constructions, such as Composite Twin Higgs models [21,22,61,62], or models involving the see-saw mechanism devloped in [18]. In such models the top partners can be significantly heavier without introducing more fine-tuning.

CONCLUSIONS
In this paper we have summarised the structure of the Higgs couplings (parameterised by κ V and κ F ) in Composite Higgs models. Although different CH models have very different predictions for the UV theory and the spectrum of higher mass resonances, we have identified generic forms for κ V and κ F which hold for many different choices of the coset group and fermion representations.
We have also looked into tree level effects on these couplings coming from extra states. In particular we studied the interesting possibility that an extra singlet pNGB may acquire a VEV. The modifications to κ V and κ F are to leading order just a sum of the corrections in elementary singlet + doublet models, and the usual correction expected in composite models. The same can be said for the case in which the Higgs mixes with an extra doublet.
We combined the Run 1 and recent Run 2 LHC data to set limits on CH models, finding that different choices for fermion representations lead to a spread of limits but a lower bound on the scale f can be set to 600 GeV. We also discussed how an observed deficit in a Higgs channel such as ttH could pinpoint the type of CH model responsible for it.