Implications of perturbative unitarity for scalar di-boson resonance searches at LHC

We study the constraints implied by partial wave unitarity on new physics in the form of spin-zero di-boson resonances at LHC. We derive the scale where the effective description in terms of the SM supplemented by a single resonance is expected to break down depending on the resonance mass and signal cross-section. Likewise, we use unitarity arguments in order to set perturbativity bounds on renormalizable UV completions of the effective description. We finally discuss under which conditions scalar di-boson resonance signals can be accommodated within weakly-coupled models.


Introduction
Unitarity of the time evolution of an isolated quantum system and in particular of the associated S-matrix is one of the cornerstones of quantum field theory. In practical perturbative calculations however, S-matrix unitarity is always approximate and asymptotic. Nonetheless, significant violations of unitarity at low orders in perturbation theory are heralds of a stronglycoupled system and can be used to constrain the range of validity of a given (effective) quantum field theory description.
Perhaps most famously, constraints imposed by perturbative unitarity in W W scattering have been used in the past to infer an upper bound on the Higgs boson mass or, alternatively, on the scale where the standard model (SM) description of weak interactions would need to be completed in the ultraviolet (UV) in terms of some new strongly coupled dynamics [1,2]. Correspondingly it allowed to narrow down the relevant mass search window and motivate the construction of the LHC with capabilities that ensured the eventual Higgs boson discovery (cf. [3] for a review).
More generally, perturbative unitarity constraints on the validity of a certain theoretical description are applicable both in non-renormalizable as well as renormalizable models. In both cases they allow to assess the limitations of a perturbative expansion. In the non-renormalizable effective field theory (EFT) approach this amounts to a truncated power expansion in (E/Λ), where E is a typical energy in a process and Λ is the EFT cut-off scale. Violations of perturbative unitarity signal the breakdown of such an expansion, when the leading powers do not represent a good approximation to the physical result. A notable standard example is the pion scattering in chiral perturbation theory, where the loop and power expansion are adequate at low enough scattering energies but violate perturbative unitarity at higher energies and eventually need to be UV completed with the inclusion of dynamical vector resonances. On the other hand within renormalizable models, the expansion proceeds in terms of positive powers of the renormalizable couplings. Sizable violations of unitarity at leading (tree) order signal the breakdown of such an expansion and the onset of strongly coupled dynamics. Here the most renown case is that of the aforementioned W W boson scattering in presence of a heavy SM Higgs boson.
The recently rekindled interest in new physics (NP) in the form of (possibly broad) di-photon resonances [4][5][6][7][8] at the LHC prompt us to reconsider the implications of perturbative unitarity for EFT interpretations of resonances decaying to di-boson final states. In particular, focusing on promptly produced scalar SM singlets decaying to two SM gauge bosons we aim to address the following questions: at which maximal energies do we expect the effective description in terms of the SM supplemented by a single scalar to break down? What can we learn about the possible UV completions of such effective theory from unitarity arguments? In particular, whether and under which conditions can a potential di-boson signal be accommodated within weakly-coupled models?
We further motivate the endeavor with the observation that in perturbative weakly-coupled models, decays of a scalar singlet into two transverse SM gauge bosons can only arise at loop level involving massive charged and/or colored particles leading to a suppression factor of Γ V T V T /M ∝ α 2 V /16π 3 . Even in the case of QCD Γ V T V T /M 10 −4 would require large couplings and/or large multiplicies of new states contributing in the loop. Both possibilities are potentially subject to constraints coming from perturbative unitarity. In particular, we will show how they enter the amplitudes of 2 → 2 scatterings of the new degrees of freedom. Similar considerations have already triggered several studies addressing the issue of the predictivity and calculability within weakly-coupled perturbative models of di-photon resonances. 1 These include studying the renormalization group equations (RGE) of the models [10] or the actual appearance of Landau poles [11][12][13]. For marginal operators such as those corresponding to the gauge couplings, Yukawas or the scalar quartic, both effects are however only logarithmically sensitive to the UV cut-off scale of the theory. The resulting constraints can also be circumvented if the models can be UV completed into theories exhibiting an infrared (IR) fixed point behavior. In case of scalar extensions, the stability of the scalar potential has also been used [14,15]. In this case the possibility of a metastable vacuum with its intricate relations to the cosmological history of the Universe requires additional assumptions going beyond quantum field theory arguments. Some aspects of partial wave unitarity for di-photon resonances which partially overlap with our work have already been discussed in [16][17][18], however with a different focus with respect to our analysis.
The rest of the paper is structured as follows: Sect. 2 contains a brief recap of partial wave unitarity arguments, which we first apply in Sect. 3 to the EFT case where a di-boson resonance is the only new degree of freedom beyond the SM. In Sect. 4 we then consider weakly-coupled benchmark models with either new fermionic or scalar degrees of freedom coupling to a diboson resonance and inducing the EFT operators in the low-energy limit. Our main results are summarized in Sect. 5. Finally, some relevant technical details of our computations can be found in Appendix A.
2 Brief review on partial wave unitarity Let us denote by T f i ( √ s, cos θ) the matrix element of a 2 → 2 scattering amplitude in momentum space, defined via where T is the interacting part of the S-matrix, S = 1 + iT . The dependence of the scattering amplitude on cos θ is eliminated by projecting it onto partial waves of total angular momentum J (see e.g. [19][20][21]) where d J µ i µ f is the J-th Wigner d-function appearing in the Jacob-Wick expansion [22], while µ i = λ i1 − λ i2 and µ f = λ f 1 − λ f 2 are defined in terms of the helicities of the initial (λ i1 , λ i2 ) and final (λ f 1 , λ f 2 ) states. The function β(x, y, z) = x 2 + y 2 + z 2 − 2xy − 2yz − 2zx is a kinematical factor related to the momentum (to the fourth power) of a given particle in the center of mass frame. The right hand side of Eq. (2) must be further multiplied by a 1 √ 2 factor for any identical pair of particles either in the initial or final state.
When restricted to a same-helicity state (zero total spin), the Wigner d-functions reduce to the Legendre polynomials, i.e. d J 00 = P J . In practice, we will only focus on J = 0 (d 0 00 = P 0 = 1), 1 For a broad survey of such models cf. [9].
since higher partial waves typically give smaller amplitudes, unless J = 0 amplitudes are suppressed or vanish for symmetry reasons. Hence, the quantity we are interested in is In the high-energy limit, √ s → ∞, one has β where the sum over h is restricted to 2-particle states, which slightly underestimates the left hand side. For i = f Eq. (4) reduces to Hence, a J ii must lie inside the circle in the Argand plane defined by (cf. also Fig. 1) which implies Under the assumption that the tree-level amplitude is real, Eq. (7) suggests the following perturbativity criterium In fact, a Born value of Re a J ii = 1 2 and Im a J ii = 0 needs at least a correction of 40% in order to restore unitarity (cf. Fig. 1).
In reality, one expects to have issues with perturbativity even before saturating the bound in Eq. (8), which is hence understood to be a conservative one. Stronger constraints can be obtained by considering the full transition matrix connecting all the possible 2-particle states, which amount to applying Eq. (8) to the highest eigenvalue of |Re (a J if ) Born |.

Effective field theory of a scalar resonance
We consider the EFT of a gauge singlet spin-0 resonance, S with mass M S , coupled to the SM fields. Assuming CP invariance, we choose S to transform as a scalar. 2 The only renormalizable terms couple S to the Higgs in the scalar potential where µ S s α m 2 S /v m 2 S /600 GeV. In the inequality we have introduced v 246 GeV and s α 0.4 [23,24] as the sine of the mixing angle between S and the Higgs, h (in the . While for a CP-even S the µ S term can contribute to the S → hh, W L W L , Z L Z L widths, it is not relevant for unitarity bounds in the high-energy limit. The d = 5 Lagrangian instead reads where we have suppressed flavor indices. This parametrization makes it clear that apart from the µ S term in Eq. (9), the interactions of a scalar singlet with the SM fields, directly relevant for di-boson resonances at the LHC, are all due to non-renormalizable d = 5 operators. Their effects are thus expected to be enhanced at high energies eventually leading to the breakdown of perturbative unitarity. In order to quantify this simple observation in the following subsections we evaluate the relevant scattering amplitudes involving SM gauge bosons, Higgs and quarks at the respective leading orders in perturbation theory. Moreover, since we are interested in studying 2 → 2 scattering processes at energies √ s M S v, we can safely set all the massive parameters (including M S ) to zero and work within the unbroken SM theory. This also implies that we can neglect any h-S mixing effects and set the masses of the final state SM particles to zero. We distinguish between two classes of tree-level processes characterized by a different energy scaling of the amplitude: scalar mediated scatterings and d = 5 contact interactions.

Scalar mediated boson scattering
Let us start, as an example, by considering the γγ → γγ scattering amplitude due to the effective operator whose matching with the operators in Eq. (10) is given by The calculation is detailed in Appendix A.1. In the (++, −−) helicity basis we find where in the last step we took the high-energy limit. Note that only the s-channel survives at high energies. The projection on the J = 0 partial waves is obtained by applying Eq. (3) and by multiplying by a 1/2 factor which takes into account the presence of identical particles both in the initial and final states. In the high-energy limit we get which, confronted with Eq. (8), leads to the tree-level unitarity bound As a matter of fact, the bound above can be made stronger if one considers the full V V → V V scattering matrix, where V and V are any of the 8 + 3 + 1 (transversely polarized) SM gauge bosons of the effective Lagrangian in Eq. (10). In such a case, the previous calculation is readily generalized in the high-energy limit where only the s-channel survives. To this end, we note that a scattering amplitude in the s-channel can be written as where a i and a j are obtained by cutting any i → j diagram in two parts along the s-channel propagator. The matrix in Eq. (16) has rank 1 and its non-zero eigenvalue is given by the trace. Hence, denoting byã 0 the eigenvalue of the V V → V V scattering matrix, in the high-energy limit we getã Correspondingly, the tree-level unitarity bound is given by We remark that in deriving these bounds we consider only the transverse polarizations of the W and Z gauge bosons. Generally, scattering amplitudes involving longitudinally polarized massive vector bosons can grow as positive powers of E/m W,Z implying apparently stronger dependence on s. However, as it can be easily verified (through an explicit calculation of the processes at hand or more generally via a clever gauge choice [25]), the scattering amplitudes involving longitudinally polarized states sourced by the gauge field strengths in Eq. (10) are suppressed by powers of m W,Z /E and thus do not lead to relevant unitarity constraints at high s. In the v → 0 limit there is just one additional tree-level s-channel contribution leading to 2 → 2 scatterings of SM particles from Eq. (10), that is due to the operator where we have neglected vertices with 4 or 5 particles and Imposing the unitarity bound on the highest eigenvalue we get Note that in the EW broken vacuum the constraint corresponds to scattering of both the physical Higgs bosons as well as the longitudinally polarized massive EW gauge bosons. Considering thus also |H † i H i as possible initial and final states,

Fermion-scalar contact interactions
Next we consider the contact interaction where we have explicitly factored out the color and SU (2) L group structure. In this case the leading scattering process is Qd → SH. By explicitly writing the polarization and gauge indices in the amplitude, one finds Only the ++ and −− polarizations survive. By explicit evaluation (cf. Appendix A.2 for the expression of the spinor polarizations) we get At high energies the J = 0 partial wave is obtained by considering the color singlet channel for a state in the linear combination 1 √ 2 |Qd + |SH , which gives Correspondingly, the tree-level unitarity bound reads √ s 8πΛ d .
Similarly, from the other two contact interactions in the last row of Eq. (10) we get √ s 8πΛ u and √ s 8πΛ e .

Unitarity bounds
As an exemplification we consider a scalar resonance S with mass M S and total width Γ S appearing in a di-photon final state at the LHC. 3 Expanding the effective Lagrangian in Eq. (10) around the broken electroweak (EW) vacuum, the part relevant for S production at the LHC is whose operators give rise to the decay widths The matching between the operators in Eq. (29) and Eq. (10) then yields In the narrow width approximation the prompt S production at the LHC can also be fully parametrized in terms of the relevant decay widths where √ s is the LHC pp collision energy and C PP parametrize the relevant parton luminosities. For illustration purposes in the following we consider in turn either gg and γγ induced processes or alternatively bb and γγ rates at a benchmark mass of M S = 750 GeV. The remaining possibilities lie in between these two limiting cases considering the values of relevant parton luminosities (their values at √ s = 8 TeV and 13 TeV LHC can be found e.g. in [12]). In the former case given a 13 TeV cross-section σ γγ ≡ σ(pp → S)B γγ one obtains the relation while for the latter we obtain These relations define the phenomenological benchmarks for the resonance partial widths into gauge boson and quark final states, to be subjected to constraints from perturbative unitarity.
To make contact with the EFT unitarity discussion of the preceding subsections we use Eqs. (30)-(31) and trade Λ g , Λ W and Λ B for Γ gg , Γ γγ and the ratio r ≡ Λ B /Λ W . In particular, we get with Barring the fine-tuned region around r = −1 (corresponding to 1/Λ γ = 0), the function f (r) has the global minimum 1.6 for r = 0.030 and reaches asymptotically the maximum 57 for r → ±∞. Hence, we can set the following unitarity bounds where the values 13 and 79 in the last equation correspond respectively to the boundary values r → ±∞ and 0.030. Generally, these bounds can be interpreted as the indication of the mass scale of new degrees of freedom UV completing the effective low-energy description and regularizing (unitarizing) the amplitude growth. If S is a member of a new strongly coupled sector (i.e. a composite state) [12,[26][27][28][29][30][31][32][33][34][35][36], the above results imply upper bounds on its compositeness scale. 4 Unfortunately, in this context unless a prospective O(TeV) mass di-boson resonance would have a very large di-boson decay width, the bounds do not appear strong enough to guarantee observable effects at LHC energies and a prospective future 50-100 TeV hadron-hadron collider [37,38] would be called for. On the other hand, in perturbative weakly-coupled realizations discussed in the next section, where S remains an elementary particle in the UV, its couplings to SM gauge field strengths cannot be generated at the tree level. Thus one expects new dynamics to appear much below the above conservative unitarity estimates.
In the case of quark scattering, we use Eq. (32) and Eq. (33). Thus the bound in Eq.
where on the r.h.s. we have normalized the partial width in qq to a broad resonance scenario. Contrary to S couplings to SM gauge field strengths, its couplings to SM fermions can be easily realized in weakly-coupled renormalizable models already at the tree level. In particular, this requires (a) S mixing with the SM Higgs doublet, (b) embedding S into an EW doublet with the quantum numbers of the SM Higgs, or (c) the introduction of new massive fermions mixing with the SM quarks and/or leptons. Case (a) is constrained by Higgs coupling measurements [23,24].
In both remaining cases, the above result can be interpreted as an upper bound on the mass scale of the extra EW (and color) charged states present in the UV completions. Unfortunately, unless S decay channels to SM quarks induce a sizable width, LHC energies will not necessarily be sufficient to exhaust these possibilities directly within the EFT. One should thus consider explicit UV realizations. In the case (b) which goes beyond the scope of this paper, precision Higgs boson and EW measurements can be used to provide additional handles [39][40][41][42][43][44]. Case (c) on the other hand, is covered in the next section.
In Fig. 2 we display the scale of unitarity violation Λ U [TeV] in the B γγ vs. σ γγ plane, for either gg or bb production and assuming either a broad or narrow resonance. In particular, for gg production we have while for bb production As reference values we take M S = 750 GeV and f (r) = 30. The horizontal lines from top to bottom indicate a cross-section signal of 6, 0.6 and 0.2 fb, assuming the same significance of the signal over the three integrated luminosities L = 3.2, 300 and 3000 fb −1 . The red curve denotes instead the reference value Λ = 20 TeV, corresponding to the typical squark-gluino reach of a futuristic 100 TeV collider [45], which applies in the case of coloured new physics generating the effective operators. Hence, if a signal is observed above the red curve it basically means that a 100 TeV collider could potentially probe the physics responsible for the restoration of unitarity. We observe that such low-scale violation of unitarity are more readily obtained in the large width scenario and that for any given σ γγ and B γγ , unitarity violation sets in earlier for bb induced production, compared to gluon fusion processes, due to much smaller PDFs.

Weakly-coupled models
In this section we consider explicit UV completions of the effective operators of Sect. 2, capturing the main features of several proposed NP models, which have recently appeared in the litera-ture. In particular, we will assume either fermion or scalar mediators 5 and CP-even couplings (the CP-odd case leads to similar conclusions as far as concerns unitarity bounds). Moreover, we restrict ourselves to the cases of bb, gg and/or γγ decays and postulate different sets of fields which separately contribute to the relevant partial widths. Note that as far as perturbativity is concerned, this latter hypothesis leads to conservative bounds. Colored mediators are experimentally much more constrained, and their masses generally need to lie close to or above the TeV scale. On the other hand, much lighter uncolored mediators are still allowed, potentially leading to resonantly enhanced one-loop contributions to radiative S decays [47,48]. The first model comprises new fermionic mediators (see e.g. [10]), all singlets under SU (2) L . To this end, we introduce N Q copies of electromagnetic (EM) neutral vector-like QCD triplets Q A ∼ (3, 1, 0) (with A = 1, . . . , N Q ) as well as N E copies of colorless vector-like fermions E B (with B = 1, . . . , N E ), with (hyper)charge Y (E B ∼ (1, 1, Y )). We assume the theory to be invariant under a U (N Q ) ⊗ U (N E ) global symmetry and the di-boson resonance is represented by a real scalar field S. The Lagrangian featuring the new fermions reads where the details of the scalar potential are not needed for our discussion.
The second model we are going to consider involves instead new scalar mediators. In analogy to the previous case, we introduce NQ copies of EM neutral QCD scalar tripletsQ A ∼ (3, 1, 0) and NẼ copies of colorless charged scalarsẼ B ∼ (1, 1, Y ), again all singlets under SU (2) L . We also assume the theory to be invariant under a U (NQ) ⊗ U (NẼ) global symmetry and the di-boson resonance is represented by a real scalar field S. The Lagrangian featuring the new scalars reads where the ellipses stand for additional terms in the scalar potential which are irrelevant for our discussion. Focusing on the CP-even couplings, the contributions to Γ γγ and Γ gg can now be written as [12] where τ i = 4m 2 i /M 2 S (for i = E,Ẽ, Q,Q), I Q = IQ = 1/2 is the index of the QCD representation, while Q E(Ẽ) is the EM charge of E(Ẽ). The loop functions read In particular, in the limit of heavy particles (τ → ∞), they decouple as S(τ ) 2/(3τ ) and F(τ ) 1/(3τ ). As a reference value we fix M S = 750 GeV, α s (M S /2) = 0.1, α EM = 1/137 and set the masses of the mediators close to the current experimental bounds from direct searches, 6 m E,Ẽ = 400 GeV and m Q,Q = 1 TeV, thus getting where we have separately considered the cases of new fermions and scalars. For heavier mediator masses the rates decouple as powers of 1/τ i = M 2 S /(4m 2 i ) and thus even larger couplings are required. For this reason, perturbativity bounds extracted using Eqs. (53)-(54) are understood to be conservative.
Finally, we also consider a special case of the fermionic model, where at least one colored fermionic mediator has the SM gauge quantum numbers of the down-like right-handed SM quarks B ∼ (3, 1, −1/3) and mixes with the b-quark, in turn inducing Sbb interactions. 7 The relevant b − B mixing Lagrangian is where Q 3 = (t L , b L ), we have used reparametrization invariance to rotate away a possible Bb R mass-mixing term, and have also neglected small CKM induced mixing terms with the first two SM generations. In the following we assume all couplings to be real in accordance with the CP-even nature of S. After EW symmetry breaking, the physical eigenstates B and b are then given in terms of the above weak eigenstates as where 6 Stable charged leptons must be heavier than about 400 GeV in order to avoid excessive Drell-Yan production [49,50], while the bounds on long-lived colored particles are more model dependent due to non-perturbative QCD uncertainties and typically range from few hundreds of GeV to 1 TeV [51,52]. 7 Analogous cases for vector-like fermions mixing with other quark flavors can easily be derived using the results of [53]. and the masses are related via In this basis, the S interactions with b and B are The θ L Bb mixing angle is constrained by EW precision measurements to sin θ L Bb = 0.05(4) [53], while θ R Bb is parametrically further suppressed as θ R Bb ∼ (m b /m B )θ L Bb . The S → bb decay width can thus be written compactly as up to terms suppressed as m 2 b / {M 2 S , m 2 B } . Note that contrary to the loop induced decay modes, Γ bb does not explicitly depend on the mediator mass. On the other hand, its implicit dependence through θ L Bb ∼ v/m B is well constrained experimentally. The resulting unitarity constraints based on Eq. (61) and saturating the upper bound on θ L Bb can thus again be considered as conservative.

Single fermion case
Let us first consider a simplified model featuring a real scalar singlet S and a non-colored Dirac fermion ψ, with the interaction Lagrangian We denote the masses of S and ψ, respectively as M S and m ψ . Focusing on the J = 0 sector, the most relevant scattering amplitude is given by ψψ → ψψ (cf. Appendix A.2). In particular, the matrix of scattering amplitudes in the (++, −−) helicity basis 8 is found to be where in the last step we took the high-energy limit. The projection on the J = 0 partial waves is readily obtained by applying Eq. (3). We report here the expression in the high-energy limit (for the full expression see Eqs. The tree-level unitarity bound (cf. Eq. (8)) relative to the highest eigenvalue of the partial wave matrix yields The behaviour of |Re a 0 ++++ | (|Re a 0 ++−− |) with the full kinematical dependence is displayed in Fig. 3, for the reference values M S = 750 GeV, m ψ = 400 GeV and y = √ 8π (y = √ 4π). In Fig. 4 we show the tree-level unitarity bound, restricted for simplicity to the ++++ helicity channel and for the three reference values m ψ = 250, 400 and 1000 GeV. The above discussion prompts us to investigate resonance width effects, which can also become important very close to the scattering poles and effectively regulate the formally diverging tree-level amplitudes. Since such effects necessarily go beyond the tree-level approximation (they can be viewed as the absorptive part of the resummed self-energy contributions of S), we do not attempt to include them explicitly. 9 Instead we superimpose contours of constant s (in shades of grey) where the (on-shell) width effects parametrized as 10 are expected to become important. Unitarity constraints derived in such regions cannot be considered meaningful. The parameter α in Eq. (66) can be viewed as a measure of the relative error ∆ introduced by using the tree-level propagator in the squared amplitude instead of one corrected in a Breit-Wigner approximation. In particular, we have α = 1/∆ − 1. So, for example, α = 3 corresponds to ∆ = 10%. For concreteness we fix Γ S /M S = 0.10. Note that due to the scaling of Eq. (66), smaller S decay widths can only lead to more stringent constraints (derived closer to the resonance poles). The bounds derived in this way can thus be considered conservative.
For m ψ = 250 GeV, S can directly decay into ψψ, thus giving the following contribution to the total decay rate In fact the requirement Γ S /M S < 10% is always more constraining than the tree-level unitarity bound whenever the s-pole resonance is above threshold, M S > 2m ψ (cf. shaded light-green region in the first plot of Fig. 4). On the other hand, for cases where the s-pole resonance is below threshold, tree-level unitarity is violated well above the region where resonance width effects are relevant. It is interesting to compare the tree-level unitarity bounds in Fig. 4 with those obtained via the RGE criterium [10] β y y = 5y 2 16π 2 < 1 .
The latter agrees up to an O(1) factor with the bound based on tree-level unitarity in the asymptotic high-energy regime y < √ 8π. Finally we note that in addition to ψψ scattering, in bounding tree-level unitarity within the fermionic mediator model one can also consider other elastic channels, such as ψS or ψψ. It turns out however, that the corresponding J = 0 partial wave amplitudes vanish in the √ s → ∞ limit and also do not receive possible enhancements due to nearby s-channel resonance poles, thus leading to no additional constraints.

Single scalar case
Let us next consider the scalar resonance S interacting with a complex scalar field φ via where A is a massive coupling and the masses of S and φ are denoted as M S and m φ , respectively. The amplitude for the φφ * → φφ * scattering reads (70) Correspondingly, the J = 0 partial wave is found to be  Note that, differently from the fermion mediators' case, the unitarity bound is never relevant in the high-energy regime √ s M S , m φ . Such situation is expected since the scalar interaction in Eq. (69) is in the form of a relevant operator, whose tree-level contribution to a 0 vanishes as 1/s in the s → ∞ limit. Thus tree-level unitarity in this case cannot be used to bound the validity of the leading order perturbative description at high energies. It can nonetheless identify problematic kinematical regions in vicinity of scattering poles. Fig. 6 shows the unitarity bound for the three reference values m φ = 250, 400 and 1000 GeV. For m φ = 250 GeV, the S → φφ * decay channel contributes to the total width of S via Analogously to the fermionic case, whenever the s-pole resonance is above threshold, M S > 2m φ , the requirement Γ S /M S < 10% is always more constraining than the tree-level unitarity bound (cf. light-green shaded area in the first plot of Fig. 6). Below threshold, the issue of the s-pole resonance width is treated in a similar way as for the fermionic case, by identifying and avoiding kinematical regions in √ s via Eq. (66) where finite width effects can become important. For m φ = 400 (1000) GeV, tree-level unitarity is then violated for values of A/M S 6.6 (11), at scales of √ s 920 GeV (2.2 TeV). Comparing the above tree-level unitarity bound with a complementary perturbativity criterium, we notice that in this case the RGEs cannot be used since, A being associated to a relevant operator, by dimensional reasons it cannot enter its beta function alone. However, A does give a finite perturbative correction to the trilinear scalar vertex Sφφ * . By evaluating the one-loop correction at zero external momentum we find In the m φ M S limit we have We can hence define a perturbativity criterium via the relation ∆A/A < 1. In any of the two limits above, the bound ∆A/A < 1 is approximately given by 11 which agrees within an O(1) factor with the bound based on tree-level unitarity (cf. also Fig. 6).
We also note that a conceptually different bound could be inferred by requiring that A does not destabilize too much the d = 2 operators. 12 For instance, by inspecting the beta function of M 2 S (see e.g. [59]) we might require β M 2 S /M 2 S = A 2 /8π 2 < 1, which yields a bound very similar to that in Eq. (76). On the other hand, an interesting feature of the mass-hierarchy bound is that, unlike the one obtained via the finite vertex correction, it gets enhanced by the multiplicity N of fields φ coupling to S, via the replacement A 2 → N A 2 .
Finally, in addition to the φφ * channel, one could also consider the φS or φφ scatterings. However, for reasons similar to the fermionic case, these processes do not lead to additional constraints and we do not discuss them any longer.

Generalization in flavor space
The results of the previous two subsections can be readily generalized to the case of N copies of the mediators. The same conclusions apply for fermion and scalar mediators, but for definiteness we are going to explicitly discuss them for fermions only. To this end, let us consider N copies of fermion fields, ψ i (i = 1, . . . , N ), interacting via the Lagrangian term where y ij is understood in the mass basis. Let us assume then some flavor structures for y ij and study the corresponding form of the unitarity bound: 1. y ij = y (∀ i and j) In such a case the amplitude matrix in Eq. (63) gets generalized into where ⊗ denotes Kronecker product and J N is the N -dimensional matrix made all by 1's. Since the only non-zero eigenvalue of J N is equal to N (recall that J N is a rank-1 matrix), all the results of the previous section are readily generalized by the replacement y → √ N y.

y ij = yδ ij
This case corresponds to the weakly-coupled models discussed at the beginning of Sect. 4. The Lagrangian features an extra U (N ) global symmetry which can be conveniently used to label the irreducible sectors of the ψψ → ψψ scattering amplitudes. Since N ⊗ N = 1 ⊕ Adj N , a general two-particle state |ψ i ψ j can be decomposed into a singlet channel and an adjoint one where T A , with A = 1, . . . , N 2 − 1, are SU (N ) generators in the fundamental representation (in the normalization Tr T A T B = δ AB ) and we properly took into account the normalization of the states.
Due to the specific flavor structure, y ij = yδ ij , one has where T s and T t denote respectively the s-and t-channel contribution to the scattering amplitudes in Eq. (63).
Let us hence discuss in turn the non-zero scattering amplitudes. For the singlet-singlet channel one finds In the asymptotic limit, √ s M S , m ψ , the t-channel decouples and one recovers the same multiplicity suppression in the unitarity bound, as in case 1. The results in the lowenergy region are instead displayed in Fig. 7, which shows the tree-level unitarity bound in the ( √ s, √ N y) plane, for different values of N . Notice that, in this normalization, the s-channel contribution is not affected by N , while the t-channel contribution is suppressed like 1/N (cf. Eq. (83)). Hence, for large enough N the unitarity bound coincides with the s-channel one and becomes relevant only in the asymptotic region √ s M S , m ψ .
The other non-zero scattering amplitude is the adjoint-adjoint one, which is found to be Hence, we conclude that the adjoint-adjoint scattering is phenomenologically less relevant: only the subleading t-channel contributes, without the high-multiplicity enhancement.
This is the most general case relevant for a di-boson resonance, for which the mediators' couplings enter the partial width Γ γγ as | i y i | 2 . On the other hand, the unitarity bound on the 2 → 2 scatterings applies to the combination i |y i | 2 . Hence, at fixed value of | i y i | 2 , the sum that enters in the amplitude for the 2 → 2 scattering is minimized when y i = y (∀ i). In this way the bound from unitarity is minimized too.
Finally, we briefly discuss the case where the mediators carry extra gauge quantum numbers, as e.g. color. This exactly matches the identity-y scenario and thus all the previous results carry over. In particular, given an N R -dimensional irreducible representation of the gauge group, the state corresponding to the gauge singlet combination always features an N R enhancement in the s-channel.

Application to mediator models
We are now ready to discuss the implication of the unitarity bounds on the required partial widths needed to reproduce any given γγ signal at the LHC. In particular, in the case of gg-initiated production processes (at M S = 750 GeV) the constraints to be fulfilled are the following: • Fermion mediators (model in Eq. (47), cf. also Eq. (65)): • Scalar mediators (model in Eq. (48)): The values in Eqs. (88)-(89) refer to the s-channel bounds of Fig. 6, for which the flavor and color enhancement apply, while in deriving Eq. (90) we have used Eq. (35) and Eq. (54). On the other hand, the following constraints (obtained by looking at the full partial wave amplitude in Fig. 6) hold irrespectively of the flavor and color copies. Note that the bounds on A Q are weaker then on A E because the partial wave amplitudes are decreasing fast for heavy mediators (away from the poles). Thus, contrary to the fermionic case, unitarity bounds on these scalar couplings crucially depend on the assumed mediator masses. Nevertheless, the bounds cannot be circumvented by decoupling the mediator masses (for fixed M S ) since the decoupling of the partial rates in Eqs. (49)-(50) is faster than that of the partial wave amplitude (cf. Eq. (71)).
In the case of fermion mediators we have 5 parameters (y E , y Q , N E , N Q and Q E ) entering the expression in Eq. (87) corresponding to a particular di-photon signal strength. Hence, a possible way to display the tree-level unitarity bounds in Eqs. (85)-(86) is to choose a value of Q E and fix y Q = y E . Fig. 8 (upper side plots) displays iso-curves reproducing the benchmark signal of σ γγ = 1 fb and Γ S /M S = 0.1 in the N Q vs. N E plane and the associated perturbativity bounds for different values of Q E . A very similar discussion applies to the case of scalar mediators (cf. lower side plots). As it emerges from Fig. 8, the only possibilities to accommodate the benchmark di-photon signal within weakly-coupled models are either via exotically-large EM charges 13 and/or a very large number of mediators' copies. These two latter options are also bounded by usual RGE arguments, which however are not sufficient to exclude such possibilities (see e.g. [10]).
We finally discuss the case of the model in Eq. (55) where the production of S is due to bb-initiated processes. Using Eq.  to be confronted with the tree-level unitarity bound where we also took into account the color enhancement of the s-channel. In this case, the perturbative unitarity constraint is very severe (see Fig. 9). In particular for our benchmark it excludes the possibility for S → bb decays to saturate a large decay width.

Conclusions
Perturbative unitarity is a powerful theoretical tool for inferring the range of validity of a given EFT, with notable examples of applications both in the physics of strong and electroweak interactions. The continued interest in di-boson resonances at the LHC motivated us to investigate the implications of partial wave unitarity for the theoretical description of such possible signals both in the minimal EFT extension of the SM as well as in its renormalizable UV completions. In the case of a TeV-scale scalar di-boson resonance observable at the LHC we have, under some very basic and natural assumptions on the structure of the EFT (mainly that S is a spin-0 SM gauge singlet and that the dim = 5 operators in Eq. (10) are the most relevant ones for the decay of S), demonstrated a potential violation of tree-level unitarity in the scattering of SM fields at energy scales of few tens of TeV. One should stress, however, that in many models (both weakly and strongly coupled) predicting observable di-boson resonances, new states are typically predicted to lie much below our energy estimates.
In a similar way one can use perturbative unitarity in order to estimate the range of validity of perturbation theory in explicit renormalizable UV completions of the low-energy EFT and accordingly set perturbativity bounds on the relevant model couplings. Especially in the case of a large total S width, the inferred bounds are typically very constraining, and are in particular endangering the calculability of many weakly-coupled models present in the literature.
Interestingly, tree-level unitarity bounds are important not only at high energies but also close to thresholds of new physics. This is especially crucial for scalars interacting via relevant operators, since the corresponding unitarity bounds are always saturated at finite scattering energies relatively close to threshold. Other perturbativity criteria such as those based on Landau poles are only logarithmically sensitive to the energy scale and typically need a few decades of running before hitting the singularity of the Landau pole.
Finally, we find that our perturbative bounds are sensitive not only to the strengths of the couplings (y) of the mediators to a di-boson resonance but also to the multiplicity N of the mediator states. For example, for fermions the bounds scale as N y 2 , exhibiting a similar 't Hooft scaling as the perturbative bounds obtained by analyzing the RGE flow of the couplings [10].
We conclude that in the event of an experimental observation of a scalar di-boson resonance at the LHC, while our estimates cannot provide a guarantee to see on-shell effects of additional new degrees of freedom at the LHC, they would immediately imply the existence of additional phenomena within the energy reach of the next generation 50-100 TeV hadron colliders, thus making a strong physics case for their construction.

Note added
While completing this paper we came across Ref. [18]. Though part of our work overlaps with it, we reach different conclusions.
spinors. The latter are chosen according to the following convention on the definition of the spinors' helicities Σ ·p u r (p) = ru r (p) , Σ ·p v s (p) = −sv s (p) , where Σ = diag( σ, σ) denotes the spin operator. Note that for anti-particles the helicity is defined with the opposite sign. A standard basis, for the two component spinors, which satisfies Eq. (109) is provided by In order to evaluate Eqs. (106)-(107), we need the rotated spinors is the rotation matrix of a bi-spinor in the 1-3 plane by an angle θ with respect to the 2nd axis. R(π) and R(θ + π) are instead the relevant rotation matrices for v(k) and v(k ), respectively. The helicity amplitudes for the ψψ → ψψ scattering (Eqs. (106)-(107)) are displayed in Table 1. Finally, we report for completeness the analytical expression of the J = 0 partial wave for the + + ++ and + + −− helicity-state scattering amplitudes, which read respectively