The limits on the strong Higgs sector parameters in the presence of new vector resonances

In this paper, we investigate how the LHC data limit the Higgs-related couplings in the effective description of a strongly interacting extension of the Standard model. The Higgs boson is introduced as a scalar composite state and it is followed in the mass hierarchy by an $SU(2)$ triplet of vector composites. The limits are calculated from the constraints obtained in the recent ATLAS+CMS combined analysis of the data from 2011 and 2012. We find that the data prefer the scenario where the Higgs couplings to the electroweak gauge bosons differ from its couplings to the vector triplet. We also investigate the unitarity limits of the studied effective model for the experimentally preferred values of the Higgs couplings. We find from the $\pi\pi\rightarrow\pi\pi$ scattering amplitudes that for the vector resonance masses between one and two TeV significant portions of the experimentally allowed regions are well below the unitarity limit. We also evaluate how the existing ATLAS and CMS Run-2 data restrict our model with the upper bounds on the resonance production cross section times its branching ratio for various decay channels. The masses in the range $1\;\mbox{TeV}\leq M_\rho\leq 2\;\mbox{TeV}$ are not excluded in parts or even full parameter space of our theory.


Introduction
Even though the LHC experiments ATLAS and CMS achieved a spectacular success by discovering the 125 GeV Higgs boson [1,2] it was more the beginning rather than the end of the struggle to uncover the character of physics beyond the Standard model (SM). To this moment, it has not even been settled down whether new a e-mail: gintner@fyzika.uniza.sk b e-mail: josef.juran@utef.cvut.cz physics takes the form of weakly coupled supersymmetry or strongly coupled composites.
If the Higgs is generated as a composite state by new strong interactions the extension of the SM can be effectively described by higher dimensional operators that do not decouple in the low-energy limit. Presumably, they would modify the SM couplings of the Higgs boson with the heavy SM fields, such as the electroweak (EW) gauge bosons and/or the third quark generation. However, while the light SM Higgs boson can guarantee unitarity of the SM to virtually arbitrary high energies, this is not true anymore if the Higgs couplings become anomalous [27,28]. Nevertheless, the least one could require from the successful effective description of the composite state phenomenology is that it will not break down at energy below the compositeness scale. Meeting this expectation might be assisted with by the presence of additional new composite states which naturally occur in strongly interacting theories, anyway. Any further progress in dealing with these questions largely depends on an experimental input. Therefore, it is interesting to find what the most recent measurements of the Higgs boson couplings imply for various effective descriptions of strongly interacting extensions of the SM. 2 In this paper, we calculate the LHC limits on the effective theory describing possible early signs of strongly interacting physics beyond the SM. The effective description we work with is a rather simplistic view of what might be observed at the LHC beyond the 125 GeV Higgs boson. We work with the vision where the Higgs boson is a scalar composite state followed in the mass hierarchy by a vector composite SU (2) triplet state. Our approach is closely related to the formalism used in [29].
In particular, the Higgs sector of the effective Lagrangian under consideration is based on the non-linear sigma model with the 125-GeV SU (2) L+R scalar singlet complementing the non-linear triplet of the Nambu-Goldstone bosons. The new vector resonances are explicitly present in the form of an SU (2) L+R triplet. The vector triplet is introduced as a gauge field via the hidden local symmetry approach [30]. Consequently, it mixes with the EW gauge bosons. It results in the appearance of the mixing-generated (indirect) couplings of the vector triplet with all SM fermions. This setup fits the situation when the global SU (2) L ×SU (2) R symmetry is broken down to SU (2) L+R . The gauge sector of this effective description is equivalent to the gauge sector of highly deconstructed Higgsless model with only three sites [31].
The above-mentioned effective scenario has also been a basis for the effective Lagrangian we introduced and investigated in [32,33] 1 . Therein, the vector triplet couples directly to the third quark doublet only and to none of the other SM fermions. In addition, even the interactions to the right quarks are not necessarily universal. Similar interaction patterns can be found in various recent extensions of the SM, including the partial compositeness and extra-dimensional scenarios. As was shown in [32,33] this arrangement helps relax the tight restrictions placed by the electroweak precision data on the vector triplet coupling to the light SM fermions, the bottom quark included. While the study in this paper follows our tBESS model (and, in the relevant parts, matches the model considered in [29]), its conclusions will also be partly, or fully, applicable to a wider class of effective models, e.g., with different fermionic sectors.
The question we address in this paper concerns a possible structure of the interactions between the new vector triplet and the Higgs boson. Under the considered symmetry, the Higgs boson coupling to the new vector triplet can differ from the Higgs couplings to the EW gauge bosons. This splitting might appear as an unwelcome complication. Nevertheless, as we will 1 We call it the tBESS model to stress its connection to the BESS model [34]. The name modifying "t" suggests a special standing of the top-quark related doublet in the model. demonstrate in the paper, the ATLAS and CMS data support it. In particular, the H → γγ constraint is the key component of the data that makes it difficult for the no-splitting scenario to satisfy the LHC measurements.
Using the results of [29,32,33] we also analyze the tree-level unitarity limits of our model resulting from the scattering of the longitudinal EW gauge bosons for the allowed values of the Higgs couplings obtained in this work. When authors of [29] addressed unitarity questions the discovery of the Higgs boson was not confirmed yet. Thus, while they used the correct mass of the Higgs boson in their analysis, they were lacking any of the experimental input on the Higgs interactions available to us today. This was compensated for by the usage of the sum rules. However, one of the used sum rules imposes the no-spitting condition. The results of our analysis compel us to abandon this assumption and calculate the unitarity limits under different conditions. This paper is organized as follows. In Section 2, we introduce the necessary components of the effective Lagrangian and work out the consequences of the Higgsto-gauge-boson coupling splitting. Section 3 is devoted to the calculations of the constraints for the Higgs-togauge boson couplings. Particularly, in Subsection 3.1, we set up the framework for the constraint calculations. In Subsection 3.2, we demonstrate the tension between the H → γγ and other LHC measurement constraints for the model parameters when the universality of the Higgs-to-gauge-boson couplings is assumed. In Subsection 3.3, we calculate the best fits and constraints for the Higgs-to-gauge-boson parameters when the Higgs coupling to the new vector triplet can differ from the Higgs coupling to the EW gauge bosons. In Subsection 3.4 we investigate the unitarity limits for our model. Finally, in Section 4, we investigate the applicability of the existing LHC limits on the masses of new vector resonances to our model. We compare the predictions of our model for the production cross section of the vector resonance times its branching ratio for various decay channels with the existing experimental upper bounds obtained by the ATLAS and CMS Collaborations. Section 5 presents the conclusions of the paper.

The effective Lagrangian
The effective Lagrangian is built to respect the global SU (2) L × SU (2) R × U (1) B−L × SU (2) HLS symmetry of which the SU (2) L × U (1) Y × SU (2) HLS subgroup is also a local symmetry. The SU (2) HLS symmetry is an auxiliary gauge symmetry invoked to accommodate the SU (2) triplet of vector resonances. Beside the scalar singlet h(x) and the vector triplet the effective Lagrangian is built out of the SM fields only. The Lagrangian can be split in three terms 2 where L GB describes the gauge-boson sector including the SU (2) HLS triplet, L ESB is the scalar sector responsible for spontaneous breaking of the electroweak and hidden local symmetries, and L ferm is the fermion Lagrangian of the model. Let us express L ESB as a sum of two terms, contains the kinetic term and the self-interactions of the Higgs boson with the mass M h = 125 GeV and free parameters c h and c ′ h . Furhter, is responsible for the masses of all gauge bosons including the new vector triplet, and describes their interactions with the Higgs boson. The interactions are parameterized by the free parameters a V , a ρ , a ′ V , a ′ ρ , . . .. Below, only the interaction terms of (4) that are at most linear in h will be considered 3 . Terms with higher powers of h, which are not important for higgs phenomenology at the LHC, can be restored at any time without affecting our conclusions.
The fermion sector L ferm in its minimal formulation contains the fermionic kinetic terms and the fermion interactions with the EW gauge bosons as well as terms responsible for the couplings of the SM fermions to the Higgs boson. While the EW part is kept identical to the SM one, the interactions of the Higgs boson with the fermions can assume non-SM values. Their parameterization is based on the following interaction Lagrangian: 2 While the full formulation of the model can be found in [33], the definitions of basic quantities used in Eqs. (2) through (5) are, for the reader's convenience, summarized in Appendix A. 3 Although the quadratic terms would be needed if one wished to maintain the possibility to eliminate the linear growth in s from the scattering amplitude ππ → hh and thus improve the unitarity limit for the model.
where M k f is a 2 × 2 diagonal matrix with the masses of the upper and bottom k th fermion doublet components on its diagonal, and U = ξ(π) · ξ(π) = exp(2iπτ /v). Note that when c k = 1, ∀k, and the rest of c's are zeros the scalar resonance interactions with fermions imitate the corresponding interactions of the SM Higgs boson. Again, only the interaction terms that are at most linear in h will be considered below.
Possible direct interactions of the fermions with the vector triplet do not play a role in the calculation of the limits investigated in this paper. Thus, this part of the effective Lagrangian is left unspecified. Nevertheless, we would like to mention the example of the setup of this sector where only the third quark generation couples directly to the vector triplet and the interaction of the right top quark is disentangled from the interaction of the right bottom quark. We suggested and analyzed the effective model with this kind of the fermion sector in [32,33]. For the sake of completeness, we would like to point out that in this model the vector resonances do interact with the SM fermions, including the light ones, even if there are no direct interactions introduced. This is due to the mixing between the EW gauge bosons and the vector resonance triplet. Of course, the mixinginduced couplings are suppressed by the transformation matrix elements; they are proportional to 1/g ′′ .
The mixing of the gauge fields occurs in the process of diagonalization of the gauge-boson mass matrix. After gauging out all six Goldstone bosons the Lagrangian L hV reads where r = a ρ /a V and where is the gauge-boson mass term. Further, X ± = (W ± f , V ± ) T , X 0 = (W 3 , B, V 3 ) T , and M 2 C (α) and M 2 N (α) are the squared-mass matrices of the charged and neutral gauge bosons, respectively, The diagonalization process results in the transformation of the gauge-boson basis, from the mass one to the flavor one, {Y } → {X}: where Y ± = (W ± m , ρ ± ) T , Y 0 = (A, Z, ρ 0 ) T . Note that we use V a µ to denote the vector resonance components in the flavor basis and ρ a µ for the vector resonance components in the mass basis 4 .
In the limit M W ± , M Z ≪ M ρ 0 , equivalent to the condition g, G ≪ √ αg ′′ , the mixing matrices read In the same limit, the next-to-leading order approximations of the gauge-boson masses read and The leading order approximation for the partial width of the ρ decay to the EW gauge bosons is The interactions of the gauge bosons with the Higgs can be read off from L hV in the mass basis with the anomalousness factors 5 c Y Y ′ . If the flavor basis splitting factor r equals to 1 (a ρ = a V ) then the mass basis couplings c Y Y ′ follow a simple pattern Fig. 1 The mass basis splitting factors ζ XY as functions of the splitting r = a ρ /a V . The graphs are plotted for M ρ 0 = 1.5 TeV and g ′′ = 15.
For a more general situation, a ρ = a V , the relations of c Y Y ′ 's to a V and a ρ become more intricate. Let us introduce the mass basis splitting factors ζ Y Y ′ such that Then The ζ Y Y ′ factors for the individual gauge bosons are summarized in Table 1. There, we have introduced auxiliary variables and Note that k N, Fig. 1 helps to understand how the flavor basis splitting, r = a ρ /a V , translates into the anomalous factors c Y Y ′ . The essential role in this issue is played by the ζ factors. Therefore, the graphs in the figure depict the dependences of ζ's on r. There, the plots of ζ W,Z are almost perfect horizontal lines at 1 which complies with ζ W,Z = 1 − (1 − r)O(x 4 ). On the other hand, ζ ρ ± ,ρ 0 = r + (1 − r)O(x 4 ) suggests the straight line of the 45 degree slope for ζ ρ ± ,ρ 0 (r). Finally, the ζ factors of the mixed interaction terms are negligible when compared to the other ζ's, which is in agreement with the finding that ζ W ρ ± ,Zρ 0 = (1 − r)O(x 2 ).
While all ζ's depend on g ′′ and M ρ these dependences are very weak. When changing (g ′′ , M ρ ) from (10, 1 TeV) to (25, 2 TeV), ζ W and ζ Z vary no more Table 1 The mass basis splitting factors ζ Y Y ′ (α, r) for the individual gauge bosons Y , Y ′ ∈ {W ± , ρ ± , A, Z, ρ 0 }. than by about 10 −4 for −3 ≤ r ≤ 3. The same conclusion applies to ζ ρ ± and ζ ρ 0 . The actual size of ζ W ρ ± and ζ Zρ 0 as well as their dependence on g ′′ and M ρ can be seen in Fig. 2. Consequently, for many phenomenological considerations the relations represent satisfactory approximations over quite a large region of r values and for all relevant values of g ′′ and M ρ . Now, let us turn our attention to the Higgs-to-gaugeboson couplings. Beside being proportional to the split-ting factors the couplings are amplified by the (squares of) gauge-boson masses. Namely, Here, the whole dependence on new physics is contained in the O(x 4 ) term. Thus, the splitting can modify the SM expectation for the ratio (30) only very slightly. On the other hand, owing to the new vector triplet's large mass new physics becomes manifest the most in the Higgs interaction with the new vector triplet as can be seen in Note that while the vector mass affects the ratios significantly their dependences on g ′′ are completely ignorable. In Fig. 3, the ratios g hρ + ρ − /g hW W and g hρ 0 ρ 0 /g hZZ as functions of r and for various vector resonance masses are depicted. If the splitting takes place new physics also manifests via the emergence of two new vertices, hW ± ρ ± and hZρ 0 , not present either in the SM or in the Lagrangian L hV when a ρ = a V . Even though the new coupling strengths lag far behind the strengths of the hρ + ρ − and hρ 0 ρ 0 couplings their presence would introduce new phenomena. The relative strengths of the hW ± ρ ± and hZρ 0 couplings with respect to the hW W and hZZ couplings, respectively, are given by Note that while these ratios are affected significantly by g ′′ their dependences on the vector mass are completely ignorable. In Fig. 4, the ratios g hW ρ ± /g hW W and g hZρ 0 /g hZZ as functions of r and for various g ′′ values are depicted. The splitting of a ρ from a V will also cause c W = c Z . Thus, in principle, r affects the model prediction for the custodial symmetry protected rho parameter. Fortunately, the effect is negligible for a wide range of r values. To demonstrate it we inspect the prediction for the ratio Numerically, our model predicts |λ W Z − 1| < 6 × 10 −5 at tree level when −3 ≤ r ≤ 3, 10 ≤ g ′′ ≤ 25, and M ρ ∈ (1, 2) TeV. The current experimental limit is λ W Z = 0.88 +0.10 −0.09 [35].
3 Limits on the Higgs-related parameters

Relevant measurements
In this paper, we would like to identify the restrictions that the current LHC measurements provide for the 10 15 25 Fig. 4 The relative strengths of the hW ρ ± and hZρ 0 couplings with respect to the hW + W − and hZZ couplings, respectively, as functions of the splitting r = a ρ /a V . The solid lines correspond to g ′′ = 10, the dashed lines to g ′′ = 15, and the dotted lines to g ′′ = 25.
free parameters of the effective Lagrangian under consideration. While the LHC cannot compete yet with the low-energy data from the LEP, SLC, and Tevatron in setting a limit on the value of g ′′ , it certainly plays the key role in restricting the Higgs-related couplings of the EW gauge bosons, fermions, and the vector resonance.
In [33], we calculated the indirect limits on g ′′ and other free parameters of the tBESS phenomenological Lagrangian, a special case of the Lagrangian considered here, fitting the observables ǫ 1 , ǫ 2 , ǫ 3 , Γ b (Z → bb + X), and BR(B → X s γ) under the assumption a V = a ρ = 1 and M ρ = 1 TeV and 2 TeV. In addition to the setup investigated here, the tBESS model contains three independent direct interactions of the vector triplet with fermions: one with the left top-bottom quark doublet, another with the right top quark, and yet another with the right bottom quark. The analysis of the limits lead to the conclusion about the preference of higher g ′′ values, namely g ′′ > 12 at 95% CL when combined with preferred values of other free parameters. Further, the analysis seemed to support the assumption of some models of partial fermion compositeness that the new strong physics resonances couple stronger to the right top quark than to the right bottom quark. 7 The requirement that our Lagrangian be treatable perturbatively bounds the values of g ′′ from above by the naive perturbativity limit, g ′′ /2 < ∼ 4π, implying g ′′ < ∼ 25. If we took this value as the final say in this issue it would not be reasonable to use g ′′ higher than about 20 in our calculations. However, one can imagine that a more rigorous analysis of the perturbativity limit could somehow modify its value one way or the other.
The limits on the Higgs-related couplings of the studied Lagrangian can be obtained from the existing measurements of the LHC experiments. In this paper, we use constraints from the ATLAS+CMS Collaborations analysis [35] based on the combination of six decay channels, namely H → γγ, ZZ * , W W * , bb, ττ and µμ, and of five production processes, namely gluon and vector boson fusions, and associated productions with W , Z, or a pair of top quarks. In our analysis, we assume that there is only one SM-like Higgs boson state at about 125 GeV of a negligible decay width.

The interim framework
In [35], the combined ATLAS and CMS measurements were utilized to calculate fits and limits on the free kappa parameters of the interim framework introduced in [36]. In this framework, kappas parameterize possible deviations of the Higgs interactions from their SM expectations. The kappas are introduced so that the squared free parameters κ i scale the SM Higgs production cross sections σ SM i and/or the partial Higgs decay widths Γ SM i associated with the ith SM particle. When κ i = 1 then the best available SM predictions for σ i ·BR i are recovered so that for the true SM Higgs boson no artificial deviations caused by ignored higher-order corrections are present.
The particular interim framework scenario that suits our effective Lagrangian is that of Section 6.1 of [35]. In this scenario, it is assumed that there are no non-SM decays of the Higgs and that the branching ratio of invisible and/or undetected decay products is zero. New particles in loops are allowed. Assuming that κ c = κ t , κ s = κ b , and κ µ = κ τ there are generally seven free parameters in this scenario: five tree-level kappas (κ W , κ Z , κ t , κ b , κ τ ) and two loop-level kappas (κ g and κ γ ). The loop-level kappas are associated with the higher-order effective vertices Hgg and Hγγ.
We will derive the limits on a V , a ρ , and c t fitting the experimental values of kappas from Table 17 of [35] (the B BSM = 0 part). The kappa values were obtained assuming κ t > 0 and using the following input parameters for the calculation of the SM quanti-ties 6 [37]: M W = 80.398 GeV, M Z = 91.1876 GeV, G F = 1.16637 × 10 −5 GeV −2 , M t = 172.5 GeV, M h = 125.5 GeV. There are five kappas relevant for calculations of our limits: κ γ , κ W , κ Z , κ g , and κ t . The remaining kappas, κ b and κ τ do not depend on a V , a ρ , or c t , at leading order.
The experimental uncertainties on kappas quoted in [35] are asymmetric. In our analysis, we will simplify the situation by equalizing both sides of the uncertainty to the one that is bigger. In addition, the 1-sigma interval of κ Z consists of two disconnected regions. In this case, we take into account only the region where the bestfit value of κ Z is placed. Thus modified experimental limits on κ's that will be used in our calculations are shown in Table 2. Table 2 The experimental limits for κ's that will be fitted to find the restrictions and the best values of a V , a ρ , and c t . Note that the kappa limits were obtained under the assumption that κ t > 0.

Fitting kappas with the effective Lagrangian parameters
To restrict the free Higgs parameters via the fit to kappas the relations between the related kappas and c i 's have to be established. For that purpose we use some of the kappa defining quantities utilized in fitting the experiment and equate them to the predictions of our effective theory. With the input parameters listed above we obtain where σ ttH is the cross section of associated production of the Higgs boson with a pair of top quarks, σ ggF is the gluon-fusion Higgs-production cross section, and Γ jj 's are the partial Higgs decay widths to the dibosons, jj = W W * , ZZ * . The index "SM" denotes the SM values.
As is well known, the H → γγ decay in the SM occurs at the loop level only. In the SM, two dominant contributions originate from the Feynman diagrams with the W boson and top quark loops. Beyond the SM the anomalous couplings of the Higgs boson to W and top quark are parameterized by the factors c W and c t , respectively. In addition, the H → γγ decay can be modified by the extra diagram with the ρ ± resonance in the loop. Thus, our effective theory predicts where c SM γ is the SM coupling of the Hγγ effective vertex and c γ (c W , c t , c ρ ± ) is its anomalous analog. Further, α EM is the electromagnetic coupling constant which, in the case of the effective Lagrangian and for the given set of the input parameters [37], depends on new physics. The leading order approximation that dominates the ratio of the squared α EM 's when M ρ ≫ M Z reads (38) Numerically, when g ′′ varies between 10 and 25 then (α EM /α SM EM ) 2 changes from 0.992 to 0.999. Thus, this effect will be ignored in our further analysis. Also, note that in the leading order the ratio is not affected by the mass of the vector resonance.
The contribution of the ρ ± resonance to h → γγ mimics the contribution of W ± ; the only difference comes from different masses and couplings of the vector particles. Thus, considering the principal contributions only -from the top, W , and ρ loops -c γ reads where where For the given input values we get F ferm (x t ) = −1. Note that, for the displayed decimal places, the numerical coefficients in this formula are not sensitive to varying the ρ mass between 1 and 2 TeV. The SM value of c γ is obtained when c t = c W = 1 and c ρ ± = 0: c SM γ = 0.81. The effective Lagrangian predictions for observables do not change under the simultaneous sign change of all c i parameters. Therefore, a sign of one of c i 's can be fixed without losing physically distinguishable configurations of theory. Owing to that we choose c t > 0 throughout the paper.
3.2 Failure of the a V = a ρ scenario Before performing the full fit on κ γ , κ W , κ Z , κ g , and κ t with three free parameters a V , a ρ , c t let us provide a simple demonstration that the no-splitting scenario, a V = a ρ , has a hard time to satisfy the experimental restrictions on the kappas under consideration. Applying (37) to this situation we obtain Using the experimental restriction on κ γ (see Table 2) the Eq. (44) results in the allowed region comprised of two parallel stripes crossing the a V − c t plane as shown in Fig. 5. Let us recall that we work under the assumption c t > 0 which reduces the a V − c t plane to a half-plane.
In addition, the graphs in Fig. 5 display rectangular intersections of the experimentally allowed regions for a V and c t . The limit on a V is derived by fitting κ 2 W and κ 2 Z while taking into account the no-splitting relations c W = c Z = a V . The relevant χ 2 -function consists of the last two terms of (46). The obtained fit reads |a V | = 0.93 ± 0.08 where the absolute value originate in the fitting of squares of variables. Analogically, the c t constraint is obtained by fitting κ 2 t and κ 2 g . The corresponding χ 2 -function contains the second and third terms of (46). The obtained best value of |c t | is c t = 0.85 +0.11 −0.12 . The graphs in Fig. 5 indicate a tension between the h → γγ limits and the combined κ W − κ Z and κ t − κ g restrictions in the no-splitting version of our effective Lagrangian. In Section 3.3, we will assess if and to what extent this tension can be removed by allowing a V and a ρ to become independent.

Full fit analysis
Assuming that a V and a ρ are free independent parameters of the effective Lagrangian under consideration, we calculate limits on a V , a ρ , and c t and search for their best values minimizing the χ 2 -function built out of the relevant measured observables and their predictions. We use the experimental values of κ W , κ Z , κ t , κ γ , and κ g shown in Table 2. Theoretical predictions of the kappas in terms of c Z , c W , c ρ ± , and c t of the Lagrangian interaction term (18) and (5), respectively, are given in (36) and (37). The dependence of c γ on c W , c ρ ± , and c t in (37) is given by (39). The c Z , c W , and c ρ ± are, in turn, related to a V and a ρ via (20).
In the interim framework the squares of kappas rather than the kappas themselves are defined and related to measured observables. Therefore, our theoretical predictions are also related to squared kappas and the best fit should be sought for minimizing the following χ 2function: whereσ i 's are the experimental errors for κ 2 i 's. The fitting of squares introduces degeneracy of solutions caused by the insensitivity of χ 2 to the relative signs between the kappas and theory parameters. Recall that c t > 0 by assumption.
If we approximate c W = c Z = a V and c ± ρ = a ρ the Eq. (45) can be simplified without any significant impact on the best-fit values as can be inferred from our conclusions obtained in Section 2 (viz., the Eq. (26)). Then, substituting (43) into (45) we obtaiñ It is obvious thatχ 2 has degenerate minima. If thẽ χ 2 -function assumes its minimum valueχ 2 min for some triplet a V , a ρ , c t then alsoχ 2 (−a V , −a ρ , −c t ) =χ 2 min . Nevertheless, by fixing c t > 0 we eliminate a half of the degenerate minima.
A simple inspection of the Eq. (46) implies thatχ 2 can be minimized when, at the same time, c t assumes the value of "weighted average" of κ t and κ g , |a V | assumes a "weighted average" of κ W and κ Z , and a ρ sets to zero the first term with the before-obtained values of c t and a V substituted in. Zeroing the first term of (46) amounts to solving a quadratic equation in a ρ . Thus, there are generally two solutions for (a V , c t ) and two solutions for (−a V , c t ), i.e. four solutions in total. Since the first term of Eq. (46) has zero contribution toχ 2 min and because the following two terms depend only on c t and two last terms only on |a V |, all four solutions result in the same value ofχ 2 min . In particular, by fitting the kappa values of Table 2 we get four minimizing triplets of {a V , a ρ , c t } with the same minimal values ofχ 2 ,χ 2 min = 4.17. Having d.o.f. = 5 − 3 = 2, the value corresponds to the hypothesis backing of 12%. The values of the degenerate best-fit triplets along with the corresponding parameter constraints at 20, 68, and 95% CL are shown in Table 3.
The graphic representation of the best-fit values is depicted in Fig. 6. There, the two-dimensional cut of the c t − a ρ − a V allowed regions by the c t = 0.85 plane is shown. Note that 0.85 is the best-fit value of c t . The contours depicted in the graph correspond to the 68 and 95% CL regions in the c t − a ρ − a V space. The splitting factors r for the best-fit points A, B, C, and D have the values −0.09, −1.81, −0.56, and −2.27, respectively. The straight line indicates the points of the no-splitting scenario, a ρ = a V . The full 95% CL region in the c t − a ρ − a V space around the best-fit point A is shown in Fig. 7. The allowed regions around the best-fit points Table 3 The best-fit values of c t , a V , and a ρ corresponding to the four minima of theχ 2 -function (46) with the constraints at 20% CL (1-sigma), 68% CL, and 95% CL. All minima (labeled as A, B, C, and D) have the same backing of 12%. There is a good reason why a V a ρ < 0 for all four best-fit points. It is because the combined contributions of a V and c t to |c γ | -the values of a V and c t being determined by the other terms of the χ 2 -function (46)overshoot the optimal value of |c γ |. Since the a V contribution dominates the c t contribution, the a ρ has to have a sign opposite to the sign of a V in order to counterbalance the excess. In addition, since we optimize |c γ |, rather than c γ itself, a ρ resulting in the optimal c γ plays as well as a ρ resulting in −c γ . Thus, we end up with two a ρ 's for each of the two best-fit values of a V .
It is not a difficult exercise to impose the no-splitting condition, a ρ = a V , onχ 2 of (46) in order to obtain a more rigorous quantification and justification of the conclusion we have reached in Section 3.2. The fitting of χ 2 results in two minima of unequal depths. This asymmetry results from the sensitivity ofχ 2 on the relative sign of a V and c t in the κ γ term. The global minimum assumes the valueχ 2 min = 21.3 which for d.o.f. = 5−2 = 3 corresponds to the hypothesis backing of 0.9%. The minimum resides at a V = 0.55±0.05 and c t = 0.89 +0. 10 −0.11 . The local minimum at (a V , c t ) = (−0.33, 0.81) amounts to the value higher by ∆χ 2 = 10.7 above the global minimum.
The fit under the no-splitting condition and the related CL contours are depicted in Fig. 8. Besides, the graph contains reminiscence of Fig. 5 where the limits obtained by the cruder approach of Section 3.2 were shown. From a naive visual inspection we would say that the results of the cruder analysis do not contradict the more sophisticated fit performed here. Not only the position of the minimum coincides with the guessed expectation based on Fig. 5, but the obtained hypothesis backing also confirms our qualitative conclusion of Section 3.2 about a low support of the data for the no-splitting scenario.
Having the rho-to-Higgs coupling disconnected from the Higgs interactions with the EW gauge bosons and acknowledging that a V = c t = 1 is compatible with the experiment, one might wonder how successful the parameter a ρ alone would be in accommodating the existing data if a V and c t were kept on their SM values. This simple exercise amounts to minimizing thẽ χ 2 -function (46) when a V = c t = 1. There are two Fig. 7 The 95% CL allowed region of the parameters a V , a ρ , and c t of theχ 2 minima A and the region projections to two-parameter planes. Dots indicate the best-fit value and its projections to the planes. minima, at a ρ = −0.12 and a ρ = −1.72, both with χ 2 min = 7. Since d.o.f. = 5 − 1 = 4 our effective Lagrangian with the SM couplings of the Higgs boson to the EW gauge bosons and top quark has 13.8% backing by the data. Raising χ 2 above its minimum by 1 the 1sigma (68% CL) constraints read +0. 12 −0.14 and +0.14 −0.12 for the former and latter a ρ best values, respectively. Hence, if we had a theoretical reason to demand a V = c t = 1 the experimental support for the vector triplet of particular a ρ values would be as good as in the model with loose a V and c t .
There are two kappas in the combined ATLAS+CMS measurements [35] not utilized in our analysis: In the approximation used in this paper, κ τ is solely related to c τ , κ 2 τ = c 2 τ . Thus, it has no impact on the fits of other parameters. Regarding κ b , while it shares the dependence on c b with κ g , its effect on the Higgs production via gluon-gluon fusion is negligible. Consequently, we had dropped the κ b term from the χ 2 -function.
On the other hand, the effective coupling κ Zγ ≡ Γ Zγ /Γ SM Zγ has a similar structure as κ γ : it depends on a V , a ρ , and c t parameters at leading order. Thus, it has a potential to affect our fits significantly. Unfortunately, the existing measurements restrict κ Zγ very poorly [38]. Because of that, neither κ Zγ was included in the χ 2function (45).

Unitarity limits for the preferred values of the Higgs couplings
In this Section, we would like to determine how the usability of our phenomenological Lagrangian is restricted by the unitarity limits when the data preferred Higgs couplings obtained in Subsection 3.3 are considered.
Opposite to the SM, our effective Lagrangian is not renormalizable and its applicability is limited to a finite range of energies, not exceeding the point where scattering amplitudes violate unitarity. Considering the low-energy phenomenology of underlying new strong physics it is natural to demand that its successful effective description does not break down below the scale significant to these new interactions. Thus, the results of the investigation of the unitarity limits could be suggestive either of the new physics scale or, in a less fortunate case, of the defects in our Lagrangian. The least we can, and need, to deduce from such an analysis is the range of usability of the effective Lagrangian we play with.
The unitarity of our Lagrangian was investigated in detail in our previous publications [32,33]. Our analysis was based on the scattering amplitudes of the longitudinal EW gauge bosons to the EW gauge bosons. We studied the unitarity of the amplitudes using the Equivalence theorem approximation where the concerned amplitudes were replaced by the pionic scattering amplitudes of the non-linear sigma model. While the approximation corresponds to the limit g, g ′ → 0 (no pion to EW gauge-boson vertices), the exchange of the Higgs and vector resonances was included.
Even more thorough investigation of the unitarity limits of the same (in relevant sectors) effective Lagrangian was performed in [29]. Beside the elastic ππ → ππ amplitudes, authors of [29] analyzed the unitarity limits implied by the non-elastic ππ → hh, ρ L ρ L , hρ L processes. Note that wherever applicable the conclusions of [29] coincide with our conclusions [32,33]. In the following, we utilize the findings of these papers for calculation of the tree-level unitarity limits for our Lagrangian. Skipping details of the very standard calculations, below we summarize and discuss the obtained results.
In Fig. 9, we show the ππ → ππ tree-level unitarity limits 7 . In our case, the ππ → ππ scattering amplitudes depend only on one of the three Higgs couplings investigated in this paper, namely |a V |. Besides, the amplitudes also depend on M ρ and g ′′ . There are three graphs in Fig. 9 corresponding to three different masses of the vector resonance: M ρ = 1, 1.5, and 2 TeV. In each graph, there are regions shown where unitarity holds up to Λ = 3 TeV, 4 TeV, and 5 TeV. The regions are superimposed by the experimentally allowed 2-sigma interval for a V obtained in Subsection 3.3: 0.77 ≤ a V ≤ 1.09, the best fit being a V = 0.93.
We can see that for the 1 TeV vector resonance unitarity holds up to at least Λ = 3 TeV when a V and g ′′ assumes allowed values 8 , excluding a small region where g ′′ → 12 and a V → 1.1. When raising M ρ , the region where unitarity holds shrinks toward higher g ′′ . At the same time, it slightly shifts toward smaller a V .
If we dropped the vector resonance from our effective Lagrangian the Higgs resonance alone could unitarize the elastic ππ amplitudes up to some Λ h that depends on a V ; for example, if a V assumes its best value of 0.93 the tree-level unitarity holds up to Λ h = 4.6 TeV. Perhaps, it might also be worth mentioning that Λ h → ∞ when a V = 1. One wonders how adding the vector triplet to the Higgs-only setup changes the Λ h limit. To assess it we plot the graph in Fig. 10. There, for given a V and M ρ , the values of g ′′ are divided into two intervals: the one where Λ < Λ h (a V ), and the other, where Λ > Λ h (a V ). We can see in the graph that if 7 All our calculations use the Γ ρ = 0 approximation. This is justifiable as long as M ρ is far below the unitarity limit in terms of the ρ width. 8 Recall that 12 ≤ g ′′ ≤ 8π. a V ≥ 1 adding the vector resonance will always lower the unitarity limit. On the other hand, if a V < 1 there is always g ′′ 0 such that for all g ′′ ≥ g ′′ 0 the unitarity limit gets bigger.
In general, the ππ scattering amplitudes are plagued with the linear growths in s. Any added ingredient or assumption that removes the linear growths has a good chance to improve unitarity limits. It can be shown [29] that assuming special relations among the parameters of our Lagrangian (the sum rules) or adding new pseudoscalar and/or axial-vector fields to our Lagrangian can eliminate the linear terms from the scattering amplitudes. For example, the relation removes the linear dependence on s from the ππ → ππ amplitude. The sum rule (47) links a V with g ′′ and M ρ via α = [2M ρ /(vg ′′ )] 2 . The resulting sets of (a V , g ′′ ) points for given M ρ 's are also shown in Fig. 9. Obviously, the vector resonance satisfying the existing experimental limits is able to follow the sum rule if its mass stays close or below 1.5 TeV. It is not possible to satisfy all sum rules resulting from the elastic and inelastic ππ scatterings at the same time [29]. In addition, the sum rule motivated by the ππ → hρ L amplitude reads Unfortunately, this is the no-splitting condition not preferred by the data. However, the unitarity-based conclusions reached in [29] rely on the assumption (48). Thus, the results obtained in our analysis call for the investigation of the inelastic unitarity limits when a V = a ρ . However, this is beyond the scope of this paper and currently work in progress. 4 The vector resonance mass limits and the upper bounds on σ(pp → ρ + X) × BR Undeniably, the search for new vector (and other) resonances has its rightful and important place in the AT-LAS and CMS Collaboration's activities. While no discovery has been made, the direct exclusion limits constantly improve. Unfortunately, the obtained limits are strongly model and parameter dependent. No wonder that the mass exclusion limits found in the literature cover only some of the interesting cases. To the best of our knowledge, there are no exclusion mass limits applicable to the vector resonance of the model studied in this paper.
One of the crucial factors on which the exclusion mass limit depends is the value of the vector resonance gauge coupling g ′′ . As we saw in Fig. 9, in our model, the superposition of the unitarity limits over the experimentally preferred region 0.77 ≤ a V ≤ 1.09 suggests that considering g ′′ below about 10 is not well justified. We cannot rely on the predictions of our Lagrangian with experimentally allowed values of a V when g ′′ = O(1). This is not of great concern to us if the motivation for our model stems from strongly interacting physics. Independently of this restriction, if we narrowed our considerations to the tBESS interaction pattern of the vector resonance to fermions the low-energy limit from the EW precision measurements reads g ′′ ≥ 12 at 95% CL [33]. This is an additional motivation for considering g ′′ values above 10.
We can evaluate how the existing ATLAS and CMS data restrict our model when we compare the predictions of our model with the upper bounds on the resonance production cross section times its branching ratio for various decay channels. The bounds are rather model independent once spin of the resonance under consideration is specified. Of course, one should keep in mind that the calculations involved proceed under the assumption of a narrow-width resonance.
In Table 4, we present the cross section times branching ratios for various decay channels of the considered model at the LHC collision energy of 13 TeV. The predictions are given for three different values of the resonance masses, namely 1, 1.5, and 2 TeV, and three values of g ′′ , namely 10, 15, and 20. The g ′′ values were chosen to span the region allowed by the combination of Table 4 The production cross section times the branching ratio for different decay channels of the 1, 1.5, and 2 TeV vector resonances of our model considering three different values of g ′′ . The predictions for the ZH and W H channels are given for a V = 1 and a ρ = 0. No direct interactions of the vector resonance with fermions are assumed. The cross sections in the table are calculated for the 13 TeV pp collisions. the limits considered and obtained in the previous sections. Note that for g ′′ > 20 our theory not only runs into its perturbativity limit, but heavier vector resonances depart from the narrow width requirement as can be read off of the upper x-axes of Fig 9. All model predictions quoted in Table 4 correspond to the scenario when the decay of the vector resonance to fermions is negligible. In the tBESS-like fermion sector, this would correspond to turning off the direct coupling of the vector resonance to the third quark generation 9 , i.e. b L = b R = 0. Consequently, in this scenario, all fermions can couple to the vector resonance through the mixing-induced interactions only and the decay widths of the neutral/charged vector resonances are dominated by their decays to the EW gauge bosons; BR(W W/W Z) > 99%.
The experimental upper bounds on the cross section times branching ratio in W W channel read 0.033, 9 In the tBESS model [33], there is no direct interaction of the vector resonance with the light fermions. The direct coupling of the vector resonance to the third quark generation is parameterized by b L and b R and by the parameter p which enables the splitting of the direct interactions of the right top and the right bottom quarks. Additional fermion sector parameters λ L,R that have been introduced in the tBESS model are, for simplicity, kept at zero values throughout this paper. 0.012, and 0.005 pb for M ρ = 1, 1.5, and 2 TeV, respectively [39]. This excludes the 1 TeV resonance for g ′′ 16 and 1.5 TeV resonance for g ′′ 12. The 2 TeV resonance is unrestricted for g ′′ ≥ 10. Nevertheless, setting b L,R to their maximally low-energy precision data allowed values -b L,R ≈ 0.1, as found in [33] 10 -can lower BR(W W ) of the 1 TeV resonance down to about 70% for g ′′ = 10, to 30% for g ′′ = 15, and to 12% for g ′′ = 20. In the 2 TeV resonance case, BR(W W ) would be lowered to about 97, 87, and 67%, respectively. Thus, we can see that the fermionic interactions of the vector resonance can noticeably decrease the predictions (and, thus, release the experimental restrictions) of the model in this channel. The same effect can be expected to occur in the ZW channel.
The exclusion limits from the W W /W Z channels mentioned in the previous two paragraphs were obtained from 13. At tree level, the ZH and W H decays of the vector resonances occur only when a V = a ρ . In Table 4, we present the predictions of our model for a V = 1 and a ρ = 0. They can be used to derive the predictions for preferred values of a V and a ρ that have been found in Section 3.3 and quoted in Table 3. Following the formalism introduced in Section 2 we find that Then, taking into account the negligibility of the contributions of Γ ZH and Γ W H to the total decay widths of the neutral and charged resonances, the (a ρ − a V ) 2 scaling applies to the values of the production cross section times the branching ratio for these channels. Thus, when (a V , a ρ ) = (0.93, −0.08), the numbers in the ZH/W H sector of Table 4 Table 4 correspond to the scenario with no direct fermion interactions with the vector resonance. In the tBESS-like fermion sector the tb channel production would generally be higher.
Finally, there are also the upper bounds for the tt channel based on 3.2 fb −1 of data [50]: 1.20, 0.33, and 0.17 pb for M ρ = 1, 1.5, and 2 TeV, respectively. When we compare these bounds with the predictions in Table 4 we can see that for g ′′ ∈ (10, 20) they do not restrict our model. It applies even in the case of the tBESS-like fermion sector with b L,R = 0.1 and p = 1. Then the predictions in the tt channel will rise to about 10 −2 , 10 −3 , and 10 −4 pb for M ρ = 1, 1.5, and 2 TeV, respectively.

Conclusions
We have studied the experimental and unitarity limits on the parameters of the strong Higgs sector of the phenomenological Lagrangian where beside the composite 125 GeV Higgs boson the SU (2) L+R triplet of composite vector resonances is explicitly present. The ESB sector of our effective Lagrangian has been based on the SU (2) L ×SU (2) R → SU (2) L+R non-linear sigma model while the scalar resonance has been introduced as the SU (2) L+R singlet. The vector resonance has been built in employing the hidden local symmetry approach.
For the interactions of the Higgs boson with the EW gauge fields, the vector triplet, and the top quark the assumed symmetries allow one to introduce three free independent modification factors, a V , a ρ , and c t . After the transformation from the flavor to mass gauge-boson eigenstate basis, the first two factors combine into modifiers c Z , c W , c ρ 0 , c ρ ± , c Zρ 0 , and c W ρ ± of the vertices hZZ, hW + W − , hρ 0 ρ 0 , hρ ± ρ ± , hZρ 0 , and hW ± ρ ∓ , respectively. If r = a ρ /a V = 1, then c Z differs from c W . Nevertheless, for quite a large interval of r's around one, e.g., |r| ≤ 3, the effect is very small. The corrections to the custodial symmetry protected rho parameter induced by the differing a V and a ρ are negligible and well within the experimental limits. For many phenomenological considerations, the approximations c Z = c W = a V , c ρ 0 = c ρ ± = a ρ , and c Zρ 0 = c W ρ ± = 0 are satisfactory over quite a large region of r values and for all relevant values of g ′′ and M ρ .
The limits on the free Higgs coupling factors a V , a ρ , and c t have been calculated using constraints on the kappa parameters of the interim framework. The constraints were obtained in the recent ATLAS+CMS Collaborations analysis of various Higgs-related processes based on data from 2011 and 2012. We have used the fitting scenario where no non-SM decays of the Higgs is assumed and where the branching ratio of invisible and/or undetected decay products is zero. In addition, new particles in loops are allowed. Out of seven free parameters in this scenario -loop-level κ g , κ γ , and tree-level κ W , κ Z , κ t , κ b , κ τ -we have used all but the last two to find the restrictions on a V , a ρ , and c t ; κ b and κ τ have been ignored because their impact on the fitting parameters was negligible.
By fitting the kappas we have established that the simple case of a ρ = a V is strongly disfavored by the data. There are four triplets of the best-fit values of a V , a ρ , and c t that can satisfy the fitted data. Namely, (a V , a ρ , c t ) = (0.93, −0.08, 0.85), (0.93, −1.68, 0.85), (−0.93, 0.52, 0.85), and (−0.93, 2.11, 0.85). The 1-sigma deviation (i.e., χ 2 min + 1) for a V at all four best values is the same: ±0.08. It also applies to c t . Its 1-sigma errors read +0.11 −0.12 for each best value of c t . As far as a ρ is concerned its 1-sigma deviations differ slightly at each of the best-fit values. However, overall they do not exceed ±0.17. The hypothesis backing for all four fits is the same and amounts to 12%. The tied score might be tilted in favor of one of the fits once κ Zγ gets measured more precisely.
Using scattering amplitudes of the longitudinal EW gauge bosons in the Equivalence theorem approximation we have studied the restrictions of the usability of our phenomenological Lagrangian imposed by the unitarity limits when the data preferred Higgs couplings obtained in this paper are considered. We have found from the ππ → ππ scattering that for M ρ = 1 TeV unitarity holds up to at least Λ = 3 TeV when 0.77 ≤ a V ≤ 1.09 and 12 ≤ g ′′ ≤ 25. As M ρ grows the region where unitarity holds shrinks toward higher g ′′ and lower a V . Even if M ρ = 2 TeV, the considered model is well below the unitarity limit at significant portions of the experimentally allowed region of a V . Nevertheless, these conclusions should be complemented by a similar analysis of the ππ scattering with h and ρ L in the final state where a V and a ρ would be independent parameters, thus allowing also for a V = a ρ preferred by the data.
Our calculations show that the masses in the range 1 TeV ≤ M ρ ≤ 2 TeV are not excluded in parts or even full parameter space of our theory. When the model's predictions face the upper bounds on the production cross section times branching ratio in different decay channels the 1 TeV resonance gets excluded in the W W channel when g ′′ 16 and in the W Z channel when g ′′ 17. Both, the W W and W Z channel measurements exclude the 1.5 TeV resonance when g ′′ 14. These restriction can get weakened once the direct interactions of the vector resonance with the fermion sector are introduced. None of the other reviewed decay channels excludes our model, at least when 10 ≤ g ′′ ≤ 20.
In the view of the results obtained in this paper we would conclude that even such a simplistic effective description of possible early phenomenology of strong BSM physics as the one studied here is capable to accommodate the existing data.