Top-antitop production from $W_L^+ W_L^-$ and $Z_L Z_L$ scattering under a strongly-interacting symmetry-breaking sector

By using a Non-linear Electroweak Chiral Lagrangian, including the Higgs, coupled to heavy quarks, and the Equivalence Theorem, we compute (in the regime $M_t^2/v^2\ll\sqrt{s}M_t/v^2\ll s/v^2$) the one-loop amplitudes $W^+W^-\to t\bar t$, $ZZ\to t\bar t$ and $hh\to t\bar t$ (to NLO in the effective theory). We calculate the scalar partial-wave helicity amplitudes which allow us to check unitarity at the perturbative level in both $M_t/v$ and $s/v$. As with growing energy perturbative unitarity deteriorates, we also introduce a new unitarization method with the right analytical behavior on the complex $s$-plane and that can support poles on the second Riemann sheet to describe resonances in terms of the Lagrangian couplings. Thus we have achieved a consistent phenomenological description of any resonant $t\bar t$ production that may be enhanced by a possible strongly interacting Electroweak Symmetry Breaking Sector.

√ sM t /v 2 s/v 2 and to NLO in the effective theory). We calculate the scalar partial-wave helicity amplitudes which allow us to check unitarity at the perturbative level in both M t /v and s/v. As with growing energy perturbative unitarity deteriorates, we also introduce a new unitarization method with the right analytical behavior on the complex s-plane and that can support poles on the second Riemann sheet to describe resonances in terms of the Lagrangian couplings. Thus we have achieved a consistent phenomenological description of any resonant tt production that may be enhanced by a possible strongly interacting Electroweak Symmetry Breaking Sector.

Introduction
The Higgs-like particle with a mass of 125 GeV found at the Large Hadron Collider (LHC) [1,2] completes a possible framework of the fundamental interactions, as this new boson has quantum numbers and couplings compatible with those expected for the Higgs of the Standard Model (SM) in its minimal version. In addition, new scalar-resonances associated to new physics effects have been constrained roughly up to 600-700 GeV [3]. For new vector bosons, the lowest energy for a possible resonance to lie at is even higher [4]. The discrepancy among the Higgs mass scale and that of any new physics appearance is suggestive of a Goldstone boson (GB) interpretation of the Higgs that (together with the Goldstone bosons associated with the W ± L and Z L components of vector bosons), may be related to some global spontaneous symmetry breaking that in turn prompts a breaking of the electroweak gauge symmetry SU(2) L ×U(1) Y → U(1) Q .
To describe such pseudo-Goldstone behavior of the Higgs boson, some effective description of the Electroweak Symmetry Breaking Sector (EWSBS) of the SM must be taken into account [5][6][7][8][9][10][11][12][13]. These Effective Field Theory (EFT) descriptions are useful even when the Higgs is not a GB. In consequence, EFTs are a convenient way of parametrizing the EWSBS.
The energy gap may also favor a non-linear Lagrangian description of the symmetry breaking, which is a very general approach to the EWSBS in the EFT. The old Electroweak Chiral Lagrangian (ECL) technique [14], built up on standard Chiral Perturbation Theory for hadron physics [15], can be extended to include the scalar Higgslike particle h transforming as a singlet of custodial SU(2) C to give the so-called Higgs Effective Field Theory (HEFT). Meanwhile, the longitudinal gauge bosons transform as a triplet. This pattern is analogous to low energy hadron physics, where pions fall in a triplet and the η meson is embedded in a singlet representation of the strong SU(2) V isospin group. The global symmetry breaking scheme, SU(2) L × SU(2) R → SU(2) C , is common to both effective field theories of the strong and electroweak interactions.
As the HEFT theories are derivative expansions, for most of parameter space (saliently excluding that of the Standard Model and perhaps other very carefully tuned sets), the interactions will generically become strong at sufficiently high energy, and we have argued that a second, very broad scalar pole is expected [16,17]. This motivates theoretical studies of new resonances with energies 700 GeV < E < 4πv ∼ 3 TeV that require methods extending perturbation theory in the HEFT Lagrangianthat we explote to Next to Leading Order, (NLO)-. One strategy is extending the low-energy amplitudes through dispersion relations (DR) compatible with analyticity and unitarity. Resonances can then be found as poles in the second Riemann sheet due to the proper analytical behavior of the amplitudes. Such unitarization methods introduce some level of arbitrariness, as unitarity, analyticity, and the low-energy behavior are not sufficient to determine a scattering amplitude with arbitrary accuracy. Nevertheless, in [18] we showed that the analytical and unitary description of higher energy dynamics provided by DRs extending the one-loop results, is essentially unique qualitatively; at least so up to the first resonance in each spin-isospin channel. Other groups have recently pursued related unitarization methods in the context of the EWSBS [19].
The top quark is quite strongly coupled to the EWSBS and offers an opportunity for numerous analysis [20][21][22][23]. Current experimental efforts have studied in detail processes where heavy quarks are produced as intermediate (subsequently decaying into jets) or final states [24][25][26][27].
It is then reasonable to introduce fermions in the theory for energy scales compatible with those where resonances may appear at the LHC and may be described within the EFT framework. Our work analyzes the coupling of the pure Goldstone sector to top quarks. Newphysics fermionic couplings in the HEFT entitle us to flexibly describe the amplitudes W + L W − L → tt, Z L Z L → tt and hh → tt in the regime In the high energy limit s M 2 Z,W ∼ M 2 h and by means of the Equivalence Theorem (ET) [28], we can compute all amplitudes V L V L → tt substituting the V L longitudinal vector bosons by GBs (denoted ω in what follows). Therefore, in this article we will take 0 = M 2 W = M 2 Z = M 2 h consistently. The last equality, for M h , can also be a consequence of a new symmetry breaking pattern (such as in Composite Higgs Models) and is within the philosophy of EFT, but holds (approximately) anyway because of the experimental Higgs mass value, which is close to that of the electroweak gauge bosons, below the TeV scale that we explore.
Watson's final state interaction theorem, implemented in our unitarization method, guarantees that amplitudes with final tt pairs feature poles in the second Riemann sheet in the same position as the elastic GB amplitudes. Dynamical resonances are thus linked to the parameter space of chiral couplings in the DR-unitarized HEFT.
We have organized the presentation as follows: section 2 discusses the introduction of a heavy fermion in the Effective Lagrangian for Electroweak-chiral interactions. The amplitudes for ω i ω j → tt and hh → tt processes at tree and one-loop levels are computed and summarized in sections 3 and 4, respectively. We dedicate sections 5 and 6 to the helicity amplitudes and to the scalar partial-wave computation. Section 7 is dedicated to the study of unitarity and the Inverse Amplitude Method (IAM [29][30][31]) implementation, both for single and coupled channels (while a derivation of the later is deferred to an appendix). Section 8 offers our final remarks and discussion.

The Electroweak Chiral Lagrangian with massive fermions
There are several equivalent forms of the universal Electroweak Chiral Lagrangian employing only the experi-mentally known particles. At leading order we adopt the gauged SU(2) L × SU(2) R /SU(2) C = SU(2) S 3 HEFT. Adding the chiral fermionic interactions, the Lagrangian reads where the U(x) ∈ SU(2) can be parametrized in terms of the would-be Goldstone fields as The Higgs potential can be expanded as We recover the SM with In most models of interest that the low energy theory formulated as HEFT is supposed to describe, the coefficients of the Higgs self-potential scale in the same way, as powers of the Higgs mass. It is a reasonable hypothesis to maintain this scaling as the constraints on these couplings have so far been found to be close to their SM values. As discussed in the introduction, we are neglecting M h [32] because we work in the M 2 W ∼ M 2 h s limit, therefore the potential V (h) is negligible and we further set it to zero.
In the Yukawa sector (last line) of the Lagrangian in Eq. (1), the quark doublets are where the two Q entries are made of the different up and down quark sectors and the Yukawa-coupling matrix in Eq. (1) has the form This matrix can be diagonalized by transforming independently the right-and left-handed up and down quarks as where V U,D L,R are four 3 × 3 unitary matrices. Thus the Yukawa part of the Lagrangian can be written as where is the Cabibbo-Kobayashi-Maskawa matrix and the new quark fields are mass eigenstates with M U and M D being the corresponding diagonal and real mass matrices.
In keeping with the m t /v √ s/v philosophy, lighter quark masses are completely irrelevant so we focus only on the heaviest quark generation, for which where the matrix element V tb has now been omitted since it is very close to unity. As can be seen, this part of the Lagrangian explicitly breaks custodial symmetry because of the (very) different values of the t and b quark masses. The F and G functions appearing in the Lagrangian are arbitrary analytical functions on the Higgs field h, which are usually parametrized as and In this work, these functions are only needed up to the quadratic terms. Also we will consider the limit of vanishing mass for the bottom quark (M b = 0). Then, the expanded Yukawa Lagrangian is where we have kept only O(ω 2 ) terms. Finally, the relevant HEFT Lagrangian that couples the EWSBS to the 3rd fermion generation, so as to describe the ωω → tt and hh → tt processes in the regime As we will see later, in order to properly unitarize the ωω → tt and hh → tt amplitudes, one has to consider also the amplitudes ωω → ωω, ωω → hh and hh → hh. The one-loop divergences appearing in all of them can be absorbed in the couplings corresponding to the Lagrangian (as we will explicitly show for the top amplitudes later) 3 Tree level and one loop contributions for ω a ω b → tt In this section we address the process V L V L → tt (where V = W, Z) at energies that are high when compared with M Z , M W and M h ; then we can use the ET and concentrate only in the GB ω i , h and the b and t quarks. More specifically, we will consider the regime In earlier work we have set all masses to zero from the start, M h = M Z = M W = M t = 0 since we were interested in the high energy regime, appropriate for LHC resonance searches. However in this work we deal with tt production and in that strict limit the amplitude vanishes and the minimal non-vanishing contribution must be at least linear in M t . More precisely, the lowest order (tree level) V L V L → tt is of the order of √ sM t /v 2 . At the one-loop level one should in principle include diagrams with ω, h and t loops. However, diagrams with t loops are higher order in diagrams with t loops can safely be ignored, even if M t /v is not a very small parameter. On the other hand, one-loop diagrams with ω and h loops are order √ sM t /v 2 , as the tree level ones, and must consistently be taken into account. Consequently, we will ignore diagrams as those in Fig. 3. This will not only make the computation manageable (because of the significantly smaller number of Feynman diagrams to be taken into account) but also it will make the renormalization of the amplitudes much simpler, so that only two new counter-terms must be introduced and the corresponding two couplings renormalized. We consider this a very sensible approximation to the V L V L → tt reaction in the M 2 t /v 2 √ sM t /v 2 s/v 2 regime and, in any case, a necessary first step to a more complete future computation that should be performed if more accuracy was ever needed.
At tree level (see Fig. 1), the scattering amplitude is given by where i, j are the custodial isospin indices of the incoming GB; p 1 , p 2 and λ 1 , λ 2 are top, antitop momenta and he-licities respectively. The √ 3 factor is a color factor since the tt pair is produced in a color singlet state. Next we consider the one-loop terms. The Feynman diagrams contributing to this order can be seen in Fig. 2. By using dimensional regularization with dimension D = 4 − ε, the result is where as usual and µ is an arbitrary renormalization scale. Hence, the sum of the two contributions is and The divergence in Eq. (21) can be absorbed by renormalizing the g t coupling. Using the MS renormalization scheme we define and consequently the next to leading order contribution (NLO) is given by Notice that we do not have any wavefunction or mass renormalization. To see it, observe the typical diagrams contributing in Fig. 3. All diagrams in the figure are clearly proportional to M t or higher orders. Since all our amplitudes have been computed to linear order in , attaching them to any of our external legs or propagators would increase the order in the chiral M t √ s/v 2 counting of Fig. 6. In the absence of wave or mass renormalization, squared amplitudes must be observable and hence µindependent. Then, we require the total derivatives of the NLO amplitude with respect to log µ 2 to vanish, Hence, the renormalization equation for the coupling g r which can be integrated to give On the other hand, the different spinor combinations of helicities appearing in the above amplitudes are, to the Leading Order (LO) in M t / √ s expansion, where the helicity indices + and − refer to λ = +1/2 and λ = −1/2, respectively. Therefore, the tree level amplitude in Eq.   Figure 3: The contribution of these diagrams to ωω → tt scattering (through wavefunction or mass renormalization) does not need to be considered, as all our amplitudes are already linear in M t . Attaching any of these corrections would increase the order in the M t /v expansion by at least one unit, as exposed in Fig. 6.
Thus, the amplitude ω i ω j → tt is given by 4 hh → tt process In a similar way to the triplet-states annihilation, we may consider hh → tt annihilation. The contributing diagram (direct vertex) to the LO hh → tt amplitude is depicted in Fig. 4 and the tree level amplitude is easily obtained, Combining the tree plus the one-loop amplitude yields where the renormalized coupling g r t is obviously defined as and Again, the lack of wave function renormalization at this level requires this amplitude to be scale independent. Thus the coupling dependence on µ is given by By integrating Eq. (38), the renormalized coupling evolves with the scale as In a similar way than in the would-be GB case we get, for the amplitude in Eq. (35), where

Helicity amplitudes
In order to study the unitarity of the strongly interacting ωω, hh and tt processes it is quite convenient to consider partial waves of the corresponding helicity amplitudes, as the unitarity relations do not couple different J nor custodial isospin I. For example, for elastic Goldstone-boson scattering ωω → ωω there are three custodial isospin A I amplitudes (I = 0, 1, 2), analogous to those in pion-pion scattering in hadron physics, which are defined in terms of the amplitude A(s,t, u)δ i j δ kl + A(t, s, u)δ ik δ jl + A(u,t, s)δ il δ jk .
These amplitudes can be expanded as loop + . . . (45) The projection over definite orbital angular momentum (the GBs carry zero spin) is then These partial waves also accept a chiral expansion which takes the general form The constants K, D and E and the function B(µ) depend on the different channels IJ = 00; 11; 20; 02; 22, as is shown in [16,18]. We will use the notation of that paper for the inelastic and pure-h scattering reactions too. As A IJ (s) must be scale independent we have This B(µ) function depends on the NLO chiral constants (a, b, a 4 , a 5 , etc.) and from now on we omit the superindices r on the renormalized coupling constants for simplicity.
Since the Higgs boson is assigned zero custodial isospin, ωω → hh and hh → hh occur only in the isospin zero channel I = 0. The corresponding partial waves can also be expanded as and respectively. Both ωω and hh may couple to the tt state. The hh pair is always produced in an I = 0 state as h is a custodial symmetry singlet. On the other hand, as t is a member of a custodial isospin doublet (t, b) T , a tt doublet can be projected to both I = 0 and I = 1, with and The ωω state with I = 0 is defined as As the t and b quarks interact proportionally to their masses which are so different, the reactions involving fermions considered here are not custodial invariant (and in fact we are neglecting all the time bb pair production since M t M b 0). The initial I = 0 ωω or hh states couple to |tt in a superposition of |I = 0, I z = 0 and |I = 1, I z = 0 .
If we concentrate in the J = 0 case, the parity of the I = J = 0 ωω or hh pairs is positive; so must be that of the tt state. The reason is that, though Eq. (14) contains parity-violating terms, the Feynman diagrams in Fig. 2 and 5 only employ the parity-conserving pieces of that equation.
As the product of intrinsic fermion and antifermion parities is P = −1, their relative orbital angular momentum must be L = 1 (to obtain P = +1) and J = 0 then needs spin S = 1. Therefore we denote the tt helicity states by |λ 1 , λ 2 and build up the S = 1, S z = 0 state The orthogonal spin state |S = 0, S z = 0 , having negative parity in an s-wave, does not couple to the I = 0 ωω or hh states. Putting it in a p-wave would entail one more order in the chiral (M t , √ s) counting (depicted in Fig. 6).

Partial waves in perturbation theory
For the I = 0 projection of the processes computed in this work we introduce the partial wave helicity amplitudes where we only need the case J = 0 and λ = λ 1 −λ 2 = 0 (at an energy lower than intrinsic new physics scales in the 1 TeV region, only one or at most a few partial waves suffice to accurately represent the whole amplitude). Then, the rotation is simply represented by the identity matrix. These partial waves are trivially obtained from our amplitudes in Eqs. (29)(30)(31) and (40-42). Taking into account Eq. (30), Then, the partial wave corresponding to the |S = 1, S z = 0 tt state in Eq. (56) is given by This partial wave can be expanded as where the first two contributions to Q have the form respectively. Projecting Eq. (24), the coefficients are given by The real part of this scalar partial-wave is shown in Fig. 7. We can see in the figure the effect of the parameters g t and c 1 . At present, g t is not very constrained, due to the low cross section of hh production. But, according to [33], c 1 ∈ (1.0, 1.7) at 2-σ confidence level. For comparison, we also give the line corresponding to c 1 = 1, g t = 0 (at µ = 3TeV) and with a = 1, so that the coupling to the ωω sector is as in the Standard Model. As visible, the amplitude may grow with √ s and may eventually violate perturbative unitarity (see section 7 below for an extensive discussion). The parameters from the top sector can easily enhance this behavior: for example, a value g t = 0.03 − 0.05 will already cause trouble with unitarity below 3 TeV, as the amplitude is seen to approach 1 rapidly. Likewise, example imaginary parts of this partial wave are shown in Fig. 8 (as this is not an elastic amplitude, the imaginary part can actually be negative depending on the parameter set).
In a similar way it is possible to obtain the J = 0 partial wave for the hh → tt reaction, also accepting the chiral expansion where the first two terms share the general form of Eq. (62-63), and the constants are given by With this partial wave it is possible to describe J = 0 scattering including the ωω, hh and tt states. We collect them all in a partial-wave amplitude-matrix, in the order just quoted, where S is the appropriate tt → tt partial wave. This partial wave is of order M 2 t /v 2 , one order higher than we have retained, so we may consistently set it to zero against others that are M t √ s/(v 2 ). As the interactions considered here are T -reversal invariant, this matrix is symmetric. Each of the elements has a right unitary cut starting at s = 0 associated with the threshold for producing ωω, hh and tt (which are all considered massless here in accordance with the use of the Equivalence Theorem in the mid-to high-energy region). The physical partial waves have support on this cut along s = E 2 + iε, where E is the reaction's center of mass energy. For these physical values, the unitarity condition for the F matrix reads This matrix equation can be expanded to The ellipsis stand for third terms that are higher order in M t /v. Thus, this system is only approximately equivalent to exact unitarity. As shown in the appendix (Fig. 12) our unitarization method in Eq. (98), as implemented on a computer, satisfies this system with very good accuracy. However, the NLO computations shown in Refs. [16,17] for A(s), M(s) and T (s); and those of this work for Q(s) and N(s), are unitary only in the perturbative sense. We expand the matrix F, using an obvious notation, as Substitution in Eq. (74) yields the perturbative unitarity relations Im The set of Eq. (75a-75f) become (noting that the second, third and fourth equations simplify, as the last term drops because T (0) 0 = 0 for the scalar hh → hh channel), The first three equations were obtained in [16,17]. The last two are equivalent to (with K 0 the LO constant in the M amplitude). These two equations can be explicitly checked and are a very non trivial consistency test of the results found in this work. However, as Fig. 9 shows, perturbative unitarity in a chiral expansion separates from exact unitarity as energy increases, invalidating perturbation theory.

Partial-wave unitarization
In this section we will show how it is possible to promote the perturbative unitarity obtained from the one-loop computations to the exact (in a sense that will be clarified below) unitarity. To show how that can be done we will start first with the particular case a 2 = b. From the results in [16,17] it is very easy to realize that the ωω → hh amplitude vanishes, and hh decouples from ωω elastic scattering. Thus, the reaction matrix can be written in 2 × 2 form as where A 00 , Q 0 and S 0 are the elastic ωω → ωω, crosschannel ωω → tt, and elastic tt → tt I = J = 0 partial waves respectively. As we saw in the previous section the A, Q and S amplitudes can be expanded as On the (RC) right cut (the physical region) the unitarity relation Im F = FF † applies, entailing These equations, if expanded in perturbation theory, return Eqs. (78a), (78d) and (78f) (setting M = 0 there). In order to fulfill these relations we proceed as follows. First, we consider the elastic scattering ωω amplitude A.
As shown in [16,17], Eq. (78d) can be satisfied by using the Inverse Amplitude Method (IAM) which introduces the unitarized amplitudẽ This ensures elastic unitarity, ImÃ =ÃÃ * , provided we have perturbative unitarity, i.e. Im A (1) = (A (0) ) 2 , as is the case here. Now, in order to unitarize Q we introducẽ Again, it is very easy to show that this partial wave fulfills the unitarity relation ImQ = AQ * by using the perturbative result Im Q (1) = Q (0) A (0) as follows, Thus, we have that Hence, we recover Eq. (82b) as announced. Therefore we are left with two unitarized amplitudesÃ andQ for the processes ωω → ωω and ωω → tt, respectively. These amplitudes also respect the perturbative expansioñ and, in addition they feature the proper analytical structure on the whole complex plane. In particular they have a right (unitarity) cut and also the expected left cut. Moreover, they can be analytically extended to the second Riemann sheet beyond the unitarity cut and they can have poles there that can be understood as resonances (whether "dynamic" or "intrinsic") developing in some regions of the chiral coupling space. Those resonances are typical of strongly interacting scenarios for the symmetry breaking sector of the SM and are under active research at the LHC. The unitarity condition for the Q amplitude linking ωω and tt introduced in Eq. (84) can now be checked numerically, and we have do so (not shown). Our numeric precision is, for the entire energy interval of interest up to 3 TeV, of order 10 −5 without any particular effort (and this small error probably stems from our setting b not quite equal to a 2 to avoid numerical problems elsewhere, so that a tiny leak to the hh channel may be there), so that Watson's final state theorem is well satisfied and the phase of the Q amplitude is correctly set to that of the strongly interacting A(ωω → ωω).
We now exemplify the power of a method by generating a resonance in the elastic A(ωω → ωω) amplitude and feeding it to the tt channel. So that a comparison with many other theory works can be made, that were inspired by a presumed narrow LHC excess at around 0.75 TeV, we choose µ M = 750GeV, a = 0.81, a 5 = 0.0023, and all other parameters from the ωω sector as in the SM (particularly a 4 = 0 and b a 2 ). This generates a relatively narrow resonance with mass around 750 GeV, Γ/M 0.06 similar to what the community was considering before it was clear that it had been a statistical fluctuation. The resonance is shown in the top plot of Fig. 10. The lower plot shows its effect on the real part of the Q amplitude, where we have set differently from zero the parameter c 1 (to ±1). As can be seen in the figure, the resonance, a textbook Breit-Wigner in the A elastic channel (EWSBS) appears as a dip due to its interference with the background in the Figure 10: Top plot: a relatively narrow, scalar EWSBS resonance in the elastic A (ωω → ωω amplitude). Lower plot: its appearance, for the indicated c 1 , in the tt channel. Other parameters are indicated in the text. Note that this resonance interferes destructively with the background Q amplitude.
Q(ωω → tt) amplitude. Of course, such dips will appear broadened and lessened after convolution with the parton distribution functions producing the top-antitop system, the hard kernels, and the reconstruction efficiency (and detector acceptance) of the final product decays. Though a full simulation is beyond the scope of this work, it is possible that they may be observable, providing a signal that is not so often expected (as practitionners often seek excess cross-sections). A similar phenomenon has been observed by [21] in interference between perturbative SM production WW → tt (the equivalent of our Q amplitude) and the s-channel production of tt via a new resonance of O(TeV) mass. The coincidence suggests that this may be a robust result. These interference phenomena of backgrounds and narrow resonances are well known in hadron physics [34,35] and it would be interesting to discover them in the EWSBS.
The entire discussion can now be extended to the more general case a 2 = b where the cross-channels ωω → hh and hh → tt are active again. Then the reaction matrix is 3 × 3, Here, A 00 , Q 0 and S 0 are the scalar partial waves as in Eq. (81a) and following. And now, since the crossed channel amplitude ωω → hh is present, we also include its J = 0 partial wave, and those for hh → hh and hh → tt.
All of them accept a chiral expansion as Once again, the partial waves Q, N and S are suppressed by M t /v factors. In particular, where X = A, M or T . On the RC, the unitarity relation Im F = FF † applies, which leads to the set of Eqs. (75), where we omit terms suppressed by higher powers of M t /v in such a way that all equations are correct up to . This is essential to be able to decouple the unitarization of the WBGBs sector (the A, M and T partial waves) from the tt amplitudes. For the unitarization of the WBGBs sector we can use again the (coupled) IAM method. For this, we first define the 2 × 2 matrix which as usual admits a chiral expansion K = K (0) + K (1) + . . . with Im K (1) = K (0) K (0) on the RC (perturbative unitarity). The corresponding unitarized matrixK is provided by the IAM method (basically, a dispersive analysis for this matrix that employs the chiral expansion on the LC and everywhere for small s, and exact two-channel unitarity on the RC), that generalizes Eq. (83) By construction, the unitarity relation ImK =KK † holds. Now, the remaining unitarity conditions on the RC, Eqs. (75d-75f) can be written in a condensed way as This can also be expanded in perturbation theory, A solution to Eq. (96) can be written down in generalizing the simpler a 2 = b discussion. The unitarized amplitudes are then (a simple demonstration is relegated to the appendix) Notice that we are using the notation K 0 ≡ K (0) , and that Eq. (98) is a generalization of Eq. (84). In the particular case a 2 = b we have M (0) =M = T (0) =T = 0 and we then recover the previous definitions of the unitarizedÃ andQ. In the general case, the amplitudes obtained from Eq. (98) feature all the good properties mentioned above as analyticity in the whole complex plane, left and right cuts, the possibility for developing poles in the second Riemann sheet, etc.

Discussion
The possibility of coupling of the top-antitop quark pair to the longitudinal gauge bosons has long been considered [36]. We have here carried out a study that, while keeping perturbation theory relatively simple (by neglecting masses and transverse gauge couplings as well as lighter quarks, that is, concentrating on the electroweak symmetry breaking sector where new strong interactions appear), proceeds beyond it by implementing unitarity in the spirit of the final state interaction theorem.
In this article we have adopted an Effective Field Theory approach extending the ECL to incorporate a light Higgs, namely the HEFT, and coupled the resulting system to the top-antitop sector at NLO in a double expansion in M t /v and √ s/v. It is clear that for large enough values of the c 1 , g t or other parameters of the top sector, the coupling can be so intense than the approximation that all the strong interactions are contained in the EWSBS may fail. In that case, expressions such as Eq. (75d) are no longer a reliable guide, and a full coupled channel unitarization must be attempted, which we refrain from at the present time. Fig. 11 gives a feeling as to when this is expected to happen. We see in the figure that for values of c 1 further than one and half units from its SM value (c 1 = 1) or values of g t of order 0.01, the interactions coupling ωω and tt become about 1/2 of the elastic ωω ones, and more care is required in studying the amplitudes.
As long as the top-sector parameters are smaller, we may use the natural counting represented above in Fig. 6. Then, we have obtained the NLO amplitudes and studied their perturbative unitarity in Eqs. (78a-78f) up to M 2 t /v 2 terms. The satisfaction of these unitarity relation deteriorates with increasing energy.
We have also shown the effect of employing those NLO scattering amplitudes from perturbation theory as the lowenergy information for a dispersive analysis that can reach the resonance region (E ∼ 0.5 − 3 TeV) as encoded in the Inverse Amplitude Method. This may prove useful if the LHC finds new resonances in the TeV region that it is exploring. It would be natural to take as starting point that any such resonances are related to electroweak symmetry breaking (otherwise why would they lie in this energy range), and their coupling to tt would be a promising alley of experimental investigation. Figure 11: For high enough values of c 1 or g t (or other parameters not depicted), the M t counting can be overruled. Solid line: ImÃ IAM (ωω → ωω) elastic amplitude. Others: Im Q. From top to bottom: dashed lines correspond to c 1 = 1, 0.5, 0, −0.5, 1 respectively, while the dotted lines have been generated with c 1 = 1 but g t = 10 −3 , 5 × 10 −3 , and 2.5 × 10 −2 respectively. In all cases the nonvanishing EWSBS parameters are a = 0.81, b = a 2 as in the SM and a 4 (3 TeV) = 4 × 10 −4 . When Im Q ∼ ImÃ, the M t /v expansion is no more a reliable guide and a full coupledchannel unitarization is necessary.
The IAM can reproduce broad, σ -like resonances as that depicted in Fig. 11, driven by the LO parameter a, or narrow resonances such as that in Fig 10. We look forward to good-statistics LHC data to guide theory in the choice of HEFT parameters.

Appendix: unitarity in coupled channels
Let us now prove that the amplitudes defined in Eq. (98) indeed satisfy the unitarity relations in Eq. (96) by taking the imaginary part of the former equation, Because of Eq. (97), Eq. (99) turns into so that we recover Eq. (96), as was to be demonstrated. Notice also that the unitarization of Eq. (98) reproduces the correct low-energy behavior given by the chiral expansion, where the unitarized partial wavesÃ,M andT , which correspond to the WBGBs sector, are computed with the IAM procedure.
In Fig. 12 we show a check of our approximate coupled-channel unitarity relations from Eq. (75) above. The imaginary parts of Q and N linking the ωω and hh channels to tt, respectively, are shown against the two separate terms on the RHS of those equations and against their sum, which perfectly match the imaginary part.