What is the scale of new physics behind the $B$-flavour anomalies?

Motivated by the recent hints of lepton flavour non-universality in $B$-meson semi-leptonic decays, we study the constraints of perturbative unitarity on the new physics interpretation of the anomalies in $b \to c \ell \bar \nu$ and $b \to s \ell \bar \ell$ transitions. Within an effective field theory approach we find that $2 \to 2$ fermion scattering amplitudes saturate the unitarity bound below 9 TeV and 80 TeV, respectively for $b \to c \ell \bar \nu$ and $b \to s \ell \bar \ell$ transitions. Stronger bounds, up to few TeV, are obtained when the leading effective operators are oriented in the direction of the third generation, as suggested by flavour models. We finally address unitarity constraints on simplified models explaining the anomalies and show that the new physics interpretation is ruled out in a class of perturbative realizations.


Introduction
In the recent years we have witnessed a growing pattern of experimental anomalies in flavour physics, which can be schematically summarized as follows: 1. Semi-leptonic B-decays in flavour changing neutral currents (FCNC) b → sℓℓ, suggesting a deficit of muons compared to electrons. The main observables are: i) the angular distributions of B → K * µµ [1][2][3], ii) the rate of semi-leptonic decays such as B → K * µµ [1] and B s → φµµ [4] and iii) the lepton flavour universality (LFU) violating observables R K [5] and R K * [6], which are defined by the ratios B(B → K ( * ) µµ)/B(B → K ( * ) ee). We remark that the Standard Model (SM) theoretical uncertainty for R K ( * ) is very small (few percent due to QED radiative corrections [7]). Updated fits based on effective field theory (EFT) analyses, including the most recent R K * measurement, can be found in [8][9][10][11][12].
In this case it is non-trivial that three different experiments agree well among each other. A recent EFT fit of R D ( * ) can be found for instance in [16]. Very recently, there has also been a new measurement of R D * by the LHCb collaboration [17], which is remarkably compatible with the previous ones. However, this single measurement does not affect much the global fit.
In both cases the statistical significance reaches the 4σ level, while from a theoretical standpoint it is very intriguing that both the set of anomalies can be interpreted within a coherent framework. In particular, one can envisage two orthogonal structures: i) a vertical (gauge) one: global fits seem to prefer effective operators featuring only SU(2) L doublets and ii) a horizontal (flavour) one: data hints to violation of LFU with a similar hierarchical pattern as in the SM, with new physics contributions negligible in electrons (basically no effects), sizeable in muons (observable only in b → sµµ) and large in taus (effects in b → cτ ν τ and potentially in b → sν τ ν τ ). These facts motivated the community to speculate about the simultaneous explanation of these two sets of anomalies and their connection with the origin of the SM flavour. It is then maybe not too early to dream about new physics and ask "what is the scale of new physics behind the B-flavour anomalies?" Here, we address this question by using an old tool of theoretical physics, namely perturbative unitarity. 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 SM description of weak interactions needed to be completed in the ultraviolet (UV) in terms of some new strongly coupled dynamics [18,19]. What we are going to consider here instead resembles in some sense the Fermi theory of weak interactions [20]. In fact, already in the 1930's, from the low-energy measurement of G F one could have inferred what was the scale of "new physics" behind the Fermi theory. By looking at 2 → 2 scatterings via four-fermion effective operators in the Fermi theory one finds that unitarity is violated 1 at energies of the order of Λ U = 900 GeV (see e.g. [21]). As is well known, the dy-namical degrees of freedom of the SM turned out to be weakly coupled and hence much lighter than the unitarity bound, e.g. M W ≪ Λ U .
In this paper, we do something similar to the unitarity analysis in the Fermi theory by considering the four-fermion operators of the d = 6 SM-invariant EFT (SMEFT) semi-leptonic basis, under the hypothesis of a short-distance new physics explanation of the experimental anomalies in semi-leptonic B-meson decays. 2 Note that the analysis can be independently carried out for the b → sℓℓ and b → cℓν anomalies, and we do not necessarily rely on a common explanation of the two. In short, once the Wilson coefficient of an effective operator is fixed by the fit to the anomaly, we can use it in order to extract the scale of unitarity violation without the need of passing through the ambiguous separation of mass vs. coupling. The common lore is that on-shell new degrees of freedom should appear below the scale of unitarity violation (see however [22] for exceptions), with interesting consequences for direct searches at LHC and future colliders.
A simple message that we would like to emphasize is that scattering amplitudes employing SM invariant effective operators lead to scales of unitarity violation Λ U which are typically smaller than the naive dimensional analysis (NDA) estimate of the strong coupling regime g ⋆ = 4π, i.e. M ⋆ = 4πΛ O , where Λ O denotes the scale of the SM invariant effective operator required to fit the anomaly normalized to unit Wilson coefficient. This enhancement is also in part due to the correlation of the scattering amplitudes in the gauge group space, which is important to take into account when thinking about the energy reach of LHC or future colliders. A related point is the flavour structure of the effective operators. If these are oriented along the third generation fermion families (as motivated in various flavour models), one typically predicts a strong enhancement of the unitarity bound which can even reach few TeV (cf. Table 1).
Similarly to the EFT analysis, unitarity arguments can also be used in order to set perturbativity constraints on the parameter space of simplified models explaining the flavour anomalies. Note, however, that in the latter case the scattering amplitudes do not grow with the energy but reach asymptotic values proportional to the Yukawa-like couplings of the new mediators. It is possible then to translate the unitarity bounds on the coupling into an upper bound on the mass of the new states (once the ratio coupling/mass is fixed in terms of the fit to the relevant anomaly). Remarkably, in some cases the upper bound on the new mediators' mass is so strong that the perturbative interpretation of the anomaly within a given simplified model can be ruled out, or soon tested at the LHC.
The layout of the paper is the following: in Sect. 2 we start by introducing and comparing different kind of scales in the EFT. After discussing in Sect. 3 motivated flavour structures for the effective operators, we briefly introduce the partial-wave-unitarity tool in Sect. 4. We continue in Sects. 5-6 where we derive the unitarity bounds respectively in the EFT and for simplified models addressing the B-flavour anomalies. We finally conclude in Sect. 7, where we also provide a summary of our results. In Appendix A, as a paradigmatic example, we report the details of the unitarity bound calculation in the presence of an SU(2) L triplet effective operator.

A tale of scales
In what follows we will focus for simplicity on purely left-handed operators, since they provide the best fit for both the anomalies in b → sµµ and b → cτ ν transitions. The analysis can be easily generalized to scenarios including more operators by using the results given in Sect. 5. In order to start the discussion it is useful to identify and compare four (conceptually different) scales in the EFT: 3 1. Λ A : the "Fermi constant" of the process. This is the scale required to explain the anomaly, to be evaluated at the typical energy of the process which is fixed by the B-meson mass. The low-energy EFT description is based on SU(3) C × U(1) EM invariant operators. The index A on Λ A runs over the anomalies, schematically A = {R D ( * ) , R K ( * ) }, and the EFT Lagrangian featuring purely left-handed operators reads where we assumed alignment with the phases of the CKM elements that appear in the corresponding SM operators. Note that the fit of the R D ( * ) and R K ( * ) anomalies requires an opposite sign interference with the SM contribution. We also included an extra factor of 2 in the definition of the charged-current operator, so that the latter has the same normalization of the neutral-current operator when considering a SMEFT. The best fit values of the R D ( * ) [23] and R K ( * ) [10] anomalies yield respectively where the errors are at 1σ. In the following we will only consider central values.
2. Λ O : the scale of the SMEFT operator. This is the scale required to explain the anomaly using an EFT at higher energies 4 (SU(3) C × SU(2) L × U(1) Y invariant), with Wilson coefficient normalized to one. The index O on Λ O is associated with an operator of the SMEFT semi-leptonic basis and runs over all the possible Lorentz and flavour structures. For definiteness we will consider here an SU(2) L triplet operator (Q and L denoting SU(2) L doublets) and two reference flavour structures such that the operator is aligned in the direction of the flavour eigenstates responsible for the anomalies, namely O = Q 23 L 33 (for b → cτ ν transitions) and O = Q 23 L 22 (for b → sµµ transitions). The matching with Eq. (1) yields 3 Some of the results presented here will be derived in the following sections. 4 QCD running effects on the Wilson coefficients are of the order of 1 + αs 4π × log ΛO m b . For Λ O = 1 TeV, this corresponds to an O(5%) correction that will be neglected in the following.
As we will discuss in detail in Sect. 3, depending on the specific flavour ansatz, the scale Λ O can be effectively reduced with respect to the "Fermi constant" of the process. For example, the transition b → cτ ν could originate from the operator O = Q 33 L 33 , where the 3 → 2 transition in the up sector is due to a CKM mixing (in the basis where 3. Λ U : the scale of unitarity violation. This is the scale where the EFT description breaks down. The important point is that it can be expressed in terms of the scale Λ O , without passing through the ambiguous separation between coupling and mass. Using the results of Sect. 5 (which are based on a non-trivial calculation of the scattering amplitude, including gauge group multiplicity factors) we obtain which yields These are the most conservative bounds on the scale of new physics responsible for the anomalies in b → cτ ν and b → sµµ.
4. M ⋆ : the NDA mass scale in the strongly coupled regime. This is the mass scale associated with the effective operator when saturating perturbativity. After reintroducing in the NDA (see e.g. [24][25][26]), one can formally distinguish among scales (Λ), masses (M) and couplings (g), and set M = gΛ. By naively saturating perturbativity at |g ⋆ | = 4π, we can write which leads to Note that M ⋆ is a factor 5 larger than the scale of unitarity violation in Eqs. (8)-(9).
Our results for the EFT analysis are summarized in Table 1 (cf. also Sects. 3-5 for more details on the flavour structure of the effective operators and the unitarity bounds), where we report the values of the four different scales discussed above for the anomalies in either b → cτ ν or b → sµµ transitions, and depending on the flavour structure of the operator O. The two main points to be observed are the following: i) Λ U is sizably smaller than M ⋆ and ii) depending on the flavour structure of the operator O, the scale Λ U approaches the energy reach of LHC. This motivates an interesting interplay of the flavour anomalies with direct searches, which is further explored in Sect. 6 by employing simplified models.

On the flavour structure of the effective operators
The R D ( * ) and R K ( * ) anomalies can be interpreted via new physics contributions in quark flavour transitions involving the third and second generation, respectively b → c for FCCC and b → s for FCNC. In models with motivated flavour structures, it is natural to expect sizable effects in channels not directly related to the flavour anomalies. In particular, it may happen that operators involving fermions of the third family are enhanced compared to flavour violating ones. This implies that a stronger unitary bound can be derived from 2 → 2 scatterings of fermions of the third generation. For example, when considering the channel related to the anomaly in b → cτ ν τ we always get a unitarity bound from the scattering bc → τ ν τ , but we can reasonably expect that scatterings of the form bb → τ τ give stronger unitarity constraints. In order to create a link between the different channels, a flavour structure has to be assumed. In the following, we review some well-known frameworks:

Minimal Flavour Violation (MFV)
The MFV hypothesis [27] states that the strength of new physics effects are linked to the SM Yukawa couplings, which act as sources of breaking of the enlarged symmetry of the gauge-kinetic terms for fermions, SU(3) 3 for quarks. In particular, for quark doublets we get that flavour violating interactions are generated at the leading order (in powers of where we omitted SU(2) L and Lorentz indices. Here, a and b are coefficients of similar size. This implies a suppression of flavour violating quark currents compared to flavour conserving ones 2. SU(2) Q flavour symmetry In the limit of vanishing SM Yukawas for the first two quark generations, an SU(2) 3 global symmetry is restored. This approximate symmetry (or a subgroup of it) might be promoted to be a fundamental symmetry in the UV. In particular, there might be an SU(2) Q symmetry that distinguishes the quark doublets of the first two generations from the third one, and which has to be eventually broken in order to reproduce the observed pattern of SM masses and mixings. If the breaking is achieved via a spurion field X that transforms as the fundamental representation of SU(2) Q , we get (see e.g. [28]) that the the typical size of | X| is of O(λ 2 ), where λ ∼ 0.2 is the Cabibbo angle. In this case, we also expect that new physics effects in FCCC and FCNC scale like

Partial compositeness (PC)
A dynamical explanation of the flavour structure of the SM is provided by the paradigm of PC [29] in the context of composite Higgs models. In this framework the SM fields are linear combinations of elementary and composite states. The admixture elementarycomposite of every SM state is regulated by a parameter ǫ A i , where A runs over the SM fermion fields (A = Q, L, u, d, e) and i is a family index. In terms of the mixing parameters, the Yukawas of the SM are given by ( It is possible to show (see e.g. [30]) that the ǫ Q i are linked to the size of the CKM matrix elements, i.e. ǫ 2 /ǫ 3 ∼ λ 2 and ǫ 1 /ǫ 3 ∼ λ 3 . New physics effects are hence related to the size of the ǫ A i coefficients, and for quark left-handed currents one expects a similar scaling for FCCC and FCNC as in Eq. (15).
We conclude that for all the three frameworks above the transition between the third and second generation is suppressed by a factor O(λ 2 ) compared to the diagonal case involving only the third family. This implies that stronger unitarity bound can be derived from 2 → 2 scattering of the third family. For the presentation of our results in Table 1 we fix the numerical values to the MFV case, leading to a |V cb | suppression in FCCC and a |V ts | one in FCNC.
On the other hand, the situation in the lepton sector crucially depends on the unknown origin of neutrino masses. Note that the new physics effects required by the B-flavour anomalies do not violate the accidental U(1) e, µ, τ symmetry of the SM which arises in the m ν → 0 limit (or, equivalently, in the decoupling limit of lepton-number-violating effective operators). It is hence reasonable to assume that the source of LFU breaking required by the B-flavour anomalies is connected to the charged lepton masses. Two structures can be easily motivated: where the first option corresponds to PC with ǫ L i ∼ ǫ e i (implying (Y E ) ij ∼ ǫ L i ǫ L j ) and the second one to MFV in the charged lepton sector. In Table 1 we use these two benchmarks, though different patterns can be of course envisaged.

Partial wave unitarity
Here we briefly recap the partial-wave unitarity formalism. More details can be found e.g. in [31,32]. Let us denote by M f i ( √ s, cos θ) the matrix element of a 2 → 2 scattering amplitude in momentum space, where √ s is the center of mass energy and θ is the azimuthal angle of the scattering. The dependence from cos θ can be eliminated by projecting the amplitude onto partial waves of total angular momentum J. In our case it suffices to consider the lowest partial wave, defined by This expression is only valid in the high-energy limit, since we neglected kinematical factors ensuring that the partial wave is zero at threshold (see e.g. [33]). The right hand side of Eq. (17) must be further multiplied by a 1 √ 2 factor for any identical pair of particles either in the initial or final state. The unitarity of the S-matrix implies where the inequality originates from the fact that we restricted the sum over h to 2-particle states. For at the level of the Born amplitude. Although the choice in Eq. (19) is somewhat arbitrary, it yields a reasonable indication of the range of validity of the perturbative expansion. In fact, a Born value of Re a 0 ii = 1 2 and Im a 0 ii = 0 needs at least a higher-order correction of 40% in order to restore unitarity (see e.g. [33]), thus signalling the breakdown of the expansion itself.
It is also useful to note that in order to optimize the unitarity bound one can look for correlations in the partial-wave matrix (e.g. in the gauge group or flavour space). This corresponds to diagonalizing the partial-wave matrix and setting the bound on the largest eigenvalue, Reã 0 ii < 1/2, with the forward scattering i = f understood to correspond to a superposition of states which is an eigenvector ofã 0 ii . Note that neglecting a scattering channel for the partialwave matrix corresponds to removing the associated row/column. Thanks to the Cauchy interlacing theorem, we also know that the largest eigenvalue of the reduced matrix is always ≤ than the largest eigenvalue of the full matrix. Hence, by neglecting a scattering channel the unitarity bound still (conservatively) applies.

Unitarity bounds in the EFT
where flavour indices have been suppressed. Here, Q L and L L denote SU(2) L doublets, while u R , d R and e R are SU(2) L singlets. An important aspect to be taken into account for the determination of the unitarity bound is the correlation of the scattering amplitude in the SU(3) C × SU(2) L space. Let us consider, for instance, the scattering ( Similarly, in SU(2) L space we can represent it via a 4×4 matrix in the basis {ψ 1 ψ 1 , ψ 1 ψ 2 , ψ 2 ψ 1 , ψ 2 ψ 2 }, where ψ a (a = 1, 2 being an SU(2) L index) denotes either (Q L ) α or L L . A stronger unitarity bound can be hence obtained by preparing the initial and final states of the scattering in the eigenstate corresponding to the highest eigenvalue of a 0 both in SU(3) C and SU(2) L space (cf. also the discussion at the end of Sect. 4). By looking at the different scattering channels displayed in the first column of Table 2 Table 2 we also show the enhancement of the a 0 eigenvalue due to the SU(3) C × SU(2) L group structure of the partial wave. The full calculation of the unitarity bound for the triplet operator O = QL (3) (including a detailed discussion of the gauge group enhancement) is exemplified in Appendix A, while the bounds for the other cases are obtained in a similar way. We finally observe that since the tensor operator does not contribute to the J = 0 partial wave, in order to apply the unitarity bound from Λ QuLe ′ one would need to inspect higher partial waves.

Unitarity bounds in simplified models
We continue by applying unitarity constraints on the parameter space of simplified models for the explanation of the R D ( * ) and R K ( * ) anomalies. Note that this case is slightly different from the unitarity bounds in the EFT, since the scattering amplitudes do not grow with the energy. Still, one can examine the 2 → 2 scatterings of SM fermions in order to set perturbativity limits on the renormalizable couplings of the new mediators and, in turn, translate them into an upper bound on the mass of the new states (once the ratio coupling/mass is fixed in terms of the fit to the relevant observable). As two representative classes of simplified models, we consider colorless spin-1 mediators and scalar/vector leptoquarks. However, some comments are in order about the phenomenological viability of the simplified models. The criterium that we are going to follow in order to select the suitable representations  for the new mediators is that after integrating them out they are able to generate triplet and singlet left-handed operator, namely those associated with the coefficients Λ QL (3) and Λ QL (1) in Eq. (20). In all the cases that we are going to consider the phenomenologically disfavoured right-handed and scalar/tensor operator of Eq. (20) can be set to zero by a proper choice of the mediator's coupling. Given these conditions, the full set of simplified models is displayed in Table 3. Table 3: Overview of simplified models which can possibly contribute to R D ( * ) or R K ( * ) via a singlet/triplet left-handed operator. Only for specific values of the ratio of the Wilson coefficients c 1 /c 3 (obtained by integrating out a given mediator) the dangerous d i → d j νν operators are not generated (U 1 case).
From the SU(2) L decomposition (neglecting flavour indices and reinserting the Wilson co-efficients explicitly) it is evident that for c 1 /c 3 = −1 there are no b → sµµ transitions. Similarly, for c 1 /c 3 = 1 processes of the type d i → d j νν are absent. The latter are particularly dangerous, since decays like B → K ( * ) νν or K → πνν are very constraining [34,35]. From this point of view U 1 is phenomenologically favoured, since it automatically ensures the absence of d i → d j νν operators at the scale of the threshold. 5 For an incomplete list of references addressing both R D ( * ) and R K ( * ) with this leptoquark see [37][38][39][40]. Other phenomenological issues that have to be taken into account when considering a simplified model are electroweak precision tests and the radiative generation of LFU breaking effects in Z and τ decays [41,42]. In order to avoid those bounds one has to assume either a certain level of tuning within the couplings of the simplified model or rely on some non-generic features of the UV completion of the simplified model. For an example of a leptoquark model where all these bounds have been consistently addressed see e.g. [43]. Finally, one has to consider direct searches that we briefly address in Sect. 6.3. Our results on the unitarity bounds for colorless vectors and leptoquarks, which are summarized in Tables 4-5, provide an extra constraint which has to be satisfied within perturbative models.

Colorless vectors
Let us first consider the case of a real electroweak vector, V ′ µ ∼ (1, 3, 0), which couples to the SM fermions via At energies √ s ≫ M V ′ the partial-wave scattering matrix in the (Q j Q i , L l L k ) basis is given by 6 where we also took into account the SU(3) C × SU(2) L multiplicity factors. The formalism for extracting the correlation in the gauge group space follows very closely the sample calculation of the scattering with the effective triplet operator, which is detailed in Appendix A. The largest eigenvalue of Eq. (23) is and the associated unitarity bound reads Note that this is stronger than the perturbativity bound sometimes quoted in the literature, e.g. |λ Q,L ij | < √ 4π [44]. In the following, we exemplify the unitarity bounds in the case where the couplings of V ′ are aligned with the operators responsible for R D ( * ) and R K ( * ) , respectively λ Q 23 λ L 33 and λ Q 23 λ L 22 . This actually yields the most conservative bounds without flavour enhancements. The generalization to non-aligned cases is straightforward and it is reported in Table 4 for some representative cases. In this respect, we note that the multiple coupling configuration with λ Q 33 ∼ λ Q 23 might help in relaxing the bounds from Refs. [41,42]. Integrating out the V ′ and matching with Eq. (1), we obtain It is convenient to define the auxiliary functions so that the bound in Eq. (25) can be recast as (using also Eq. (26)) where A = {R D ( * ) , R K ( * ) }. The most conservative bound is obtained by maximizing the function f (r) at r = √ 3, which yields for the case of R D ( * ) (R K ( * ) ). The analysis for the Z ′ is basically identical to that of the V ′ as far as concerns neutral currents. So we do not repeat it here. The unitarity bounds for both the cases are collected in Table 4.

Leptoquarks
Let us start by first discussing the flavour structure of the leptoquark Lagrangian. Neglecting Lorentz and gauge indices, we have L LQ ⊃ y ij QL Q i L j Φ + h.c., where Φ denotes one of the four leptoquarks in Table 3. The simplest way to generate a contribution for either R D ( * ) or R K ( * ) is to switch on a single coupling, e.g. y 3j QL , with the lepton index j aligned either along the third or second generation. The 3 → 2 transition in the quark sector can be then obtained either via a V cb or V ts suppression. However, within such an approach the sign of the Wilson coefficient, which goes either like y 3j QL 2 V cb or y 3j QL 2 V ts (recall that V cb > 0 and V ts < 0 in the standard parametrization), is fixed and does not always correspond to the one necessary Table 4: Summary of unitarity bounds for colorless spin-1 mediators. FS Q and FS L denote the flavour suppression factors in the quark and lepton sectors (same as in Table 1).
to reproduce the anomaly. 7 Hence, in the following we define our simplified models based on the two leptoquark couplings y 3j QL and y 2j QL , so that the sign of the contribution can be always matched. We further assume the scaling y 2j QL ∼ y 3j QL λ 2 , as suggested by motivated flavour structures. Given the hierarchy y 3j QL ≫ y 2j QL , the scattering amplitudes are dominated by y 3j QL and the correlation of the partial wave in flavour space can be safely neglected. Indeed, for a leptoquark-mediated processes in the t-channel one should make the following replacement in the bound: |y 3j QL | 2 → |y 3j QL | 4 + |y 2j QL | 4 ∼ |y 3j QL | 2 √ 1 + λ 4 , while no such flavour enhancement is even present for an s-channel scattering. Given these considerations, for each leptoquark of Table 2 we compute the unitarity constraints on its couplings and the matching condition with the effective operators in Eq. (1). Following the conventions of Ref. [45] we have: The strongest unitarity bound comes from the t-channel mediated Q c 3 Q c 3 → L j L j scattering, which in the limit √ s ≫ M S 1 gives where we included a √ 3 enhancement factor due to the correlation of the partial wave in color space. Integrating out S 1 and matching with the operators in Eq. (1) we obtain . (31) Analogously to the previous case we consider the t-channel mediated Q c 3 Q c 3 → L j L j scattering, from which we get the unitarity bound 7 With the single coupling y 3j QL we find that S 3 cannot explain neither of the anomalies, while U 3 cannot explain R D ( * ) .
where we included a √ 3 × 3 enhancement factor due to the correlation of the partial wave in the SU(3) C × SU(2) L space. Integrating out S 3 and matching with the operators in Eq. (1) we obtain . (33) where we included a factor 2 enhancement from the correlation of the partial wave in the SU(2) L space (while there is no SU(3) C enhancement since the color flows through the diagram). Integrating out U 1 and matching with the operators in Eq. (1) we obtain . (35) Analogously to the previous case, from the s-channel process Q 3 L j → Q 3 L j we obtain Integrating out U 3 and matching with the operators in Eq. (1) we obtain After saturating the matching condition required to reproduce the anomalies, we can translate the unitarity bounds on the leptoquark couplings into an upper bound on the leptoquark masses. As a reference value we fix λ 2 = |V cb | (|V ts |) for R D ( * ) (R K ( * ) ). The results are displayed in Table 5, depending on the flavour structure of the leptoquark couplings.

Direct searches at the LHC
We will now briefly discuss the bounds from direct searches for the simplified models of Table 3 and compare them with the unitarity bounds on the new mediators' masses from Tables 4-5. We will focus in particular on decay channels involving the third family, since these are theoretically motivated by flavour models and because it is precisely in those cases that the upper bounds on the mass of the new states are more stringent. Table 5: Summary of unitarity bounds for leptoquarks. FS Q and FS L indicate the flavour suppression factors in the quark and lepton sectors (same as in Table 1). Red (light-red) boxes denote the cases excluded (disfavoured) by direct searches (see Sect. 6.3).
Let us discuss in turn the various cases. Ref. [46] considered vector triplet V ′ exclusion limits by recasting pp(bb) → τ τ searches. The conclusion is that for relatively heavy vectors M V ′ 500 GeV the resolution of the R D ( * ) anomaly with dominant third generation couplings requires a very large Z ′ decay width (where here Z ′ denotes here the neutral component of V ′ ), which is beyond the perturbative regime. This is somehow compatible with our unitarity bound M V ′ < 1.3 TeV in Table 4. Note, however, that a large Z ′ width also implies extra modeldependent decay channels which would in principle contribute to our scattering amplitudes, and would yield in turn a bound stronger than 1.3 TeV. On the other hand, for light masses M Z ′ 400 GeV a perturbative window with a relatively small Z ′ width is not yet excluded by τ τ searches. However, this requires a suppression of electroweak precision observables which are generically quite constraining [47].
If leptoquarks are light enough they can be pair-produced at LHC with sizable cross-section via QCD interactions. As already stated we assume that decay channels are dominated by third generation SM fermions. S 1 has the same quantum numbers of a sbottom and decays into either S 1 → bν τ or S 1 → tτ (both with B = 50%). Using the results of [48] we obtain M S 1 > 570 GeV, which is still compatible with the unitarity bound in Table 5. On the other hand, S 3 comprises three charge eigenstates, respectively with charges 4/3, 1/3 and −2/3. The predominant decays are S 4/3 3 → bτ , S 1/3 3 → bν τ or S 1/3 3 → tτ (both with B = 50%) and S −2/3 3 → tν τ . There will be electroweak mass splittings between the three leptoquark states, allowing the heavier ones to decay to the lighter ones, but these decays will be subdominant to those through the leptoquark couplings, if the mass splittings are small. In fact, by using the results of Ref. [49] we find that electroweak precision data exclude mass splittings within  there is a dedicated leptoquark search for third generation final states [51], which yields S 4/3 3 850 GeV. The bound on S 1/3 3 basically corresponds to the previous one for S 1 . All in all, when comparing the limits from direct searches with the unitarity bounds in Table 5, we conclude that a leptoquark S 3 with couplings dominantly aligned along the third generation cannot explain within a perturbative framework the R D ( * ) anomaly (and R K ( * ) as well, under the hypothesis of MFV in the lepton sector).
We finally discuss vector leptoquarks. Under the assumption of leading third generation couplings one can look at pp(bb) → τ τ searches, which however are not yet sensitive enough to rule out the explanation of R D ( * ) via U 1 [46]. On the other hand, vector leptoquarks can also be efficiently pair-produced at LHC via their coupling to gluons. This interaction depends however on the UV completion of the vector. The most general CP-conserving Lagrangian describing the interaction of the vector U µ with gluons (including operators up to d = 4) is given by [52] where D µ = ∂ µ −ig s t a G a µ is the QCD covariant derivative and G a µν = ∂ µ G a ν −∂ ν G a µ +g s f abc G b µ G c ν is the usual QCD field strength. As two benchmark scenarios we consider the minimal coupling (MC) and the Yang Mills (YM) type of coupling of Ref. [52]. The former case (κ G = 1), refers to the interaction stemming purely from the QCD covariant derivative of the vector, while the latter (κ G = 0) includes non-minimal interactions between the vector and the gluons which arise when the vector has a gauge origin. 9 We remark, however, that κ G is an unknown parameter.
In the exact U(2) flavour limit U 1 decays in either U 1 → tν τ or U 1 → bτ (both with B = 50%). By revisiting a √ s = 8 TeV ATLAS search [54] for QCD pair-produced third generation scalar leptoquark in the ttνν channel, Ref. [38] Table 5, we remark that a dedicated experimental search in this case would be very helpful. Finally, since U 3 contains an isospin component with the same charge of U 1 we expect similar bounds, though optimized searches for the other charge eigenstates might yield better constraints.

Conclusions
In this work we have investigated the constraints of partial-wave unitarity for the new physics interpretation of the recent hints of LFU violation in B-meson semi-leptonic decays, both within an EFT approach and by employing simplified models. In order to simplify the discussion we focussed on a single SU(2) L triplet operator (cf. Eq. (4)) which can contribute to both R D ( * ) or R K ( * ) , but without necessarily relying on the common explanation of both the anomalies. This can be straightforwardly extended to the more general situation involving multiple operators, by employing the results of Sect. 5 in which we derived the connection between the scale of unitarity violation Λ U and the coefficients Λ O of the semi-leptonic SMEFT operator basis.
The results of the EFT analysis are summarized in Table 1. In particular, we find that the most conservative bound on the scale of unitarity violation is Λ U = 9.2 TeV and 84 TeV, respectively for R D ( * ) or R K ( * ) . This corresponds to the case when the effective operators are aligned in the direction of the flavour eigenstates responsible for the anomalies. On the other hand, motivated frameworks like e.g. MFV, U(2) flavour models and PC suggest an alignment of the effective operators along the third generation, thus implying that stronger unitarity bounds can be actually extracted by considering third generation fermions' scatterings. For instance, in the case of third generation alignment in the quark sector the previous bounds become Λ U = 1.9 TeV and 17 TeV, with the latter reaching even few TeV in the case of hierarchical flavour structures also in the lepton sector.
In a similar way one can use the tool of perturbative unitarity to set constraints on the parameter space of simplified models explaining the B-flavour anomalies. As a representative class of models we considered colorless vectors and scalar/vector leptoquarks (cf. Table 3). In all those cases, it was possible to use partial-wave unitarity in order to set bounds on the renormalizable couplings of the new mediators with the SM fermions. By fixing the ratio coupling/mass in order to fit the anomaly, the unitarity bound was hence translated into an upper bound on the mass of the simplified model's mediator. The results are collected in Tables 4-5 for some reference flavour structures.
While for the anomalies in b → sℓℓ transitions it is much easier to accomodate direct searches, that is not the case for R D ( * ) . Simplified models for explaining the latter are problematic for various reasons: d i → d j νν transitions, electroweak precision observables, radiative generation of LFU breaking effects in Z and τ decays, etc. On top of that, one should take into account unitarity constraints within perturbatively calculable models. The vector leptoquark U 1 seems phenomenologically in a better shape for explaining R D ( * ) , since it is automatically free from issues like d i → d j νν transitions and also because, being an SU(2) L singlet, bounds from electroweak precision data are more easily evaded. For this specific case we provided a new bound by rescaling recent searches at LHC Run-2 with full dataset, finding that M U 1 1 ÷ 1.3 TeV (depending on the UV completion of the vector). In those cases where the leptoquark couplings are dominantly aligned along the third generation, the open window between direct searches and perturbativity is quite reduced and might be eventually closed in the near future.
Since in the massless limit the fermions in Eq. (39) are helicity eigenstates, at energies √ s ≫ v only the + − −+ polarization survives, yielding 10 The J = 0 partial-wave scattering matrix is obtained via In order to maximize the unitarity bound one can prepare the scattering eigenstates in such a way that they correspond to the highest eigenvalues of a 0 in the gauge group space. Let us discuss in turn the SU(3) C and SU(2) L structures. In the former case the partial wave can be rep- whose eigenvalues are ( √ 3, − √ 3, 0, 0). Thus, by preparing the initial and final states of the scattering in the eigenstate 1 √ 6 (1, 1, 1, √ 3), the color enhancement factor corresponds to √ 3. On the other hand, the partial wave in the SU(2) L space has the matrix form defined on the basis {ψ 1 ψ 1 , ψ 1 ψ 2 , ψ 2 ψ 1 , ψ 2 ψ 2 }, where ψ a (a = 1, 2 being an SU(2) L index) denotes either (Q L ) α or L L . In order to derive Eq. (45) it is convenient to use the Fierz identity (σ A ) ab (σ A ) cd = 2δ ad δ cb −δ ab δ cd . Since the eigenvalues of a 0 SU (2) L are (2, 2, −2, 0), by preparing the initial and final states of the scattering in the eigenstate 1 √ 2 (0, 1, 1, 0), the SU(2) L enhancement factor is 2. Summarizing, the gauge group enhancement leads to an extra √ 3 × 2 factor in the partial-wave eigenvalue, and including the latter we obtain From the condition in Eq. (19) it finally follows the unitarity bound √ s < Λ U , where As a final remark, we briefly mention an alternative way to work out the gauge group enhancement which employs irreducible representations for the scattering amplitude [33]. Denoting by ψ i (ψ j ) the fundamental (anti-fundamental) representation of an SU(N) group, a general two-particle state |ψ i ψ j can be decomposed into a singlet and an adjoint channel where T A , with A = 1, . . . , N 2 −1, are SU(N) generators (in the normalization Tr T A T B = δ AB ) and we properly normalized the states to unitary norm. The scattering amplitude in Eq. (42) has both SU(3) C and SU(2) L components. In the former case the S-matrix elements in the (color) singlet and adjoint channels are where M SU (3) C denotes the matrix element in Eq. (42) stripped from the color structure. For the SU(2) L case istead let us collectively denote the doublets (either Q or L) by ψ a , with a = 1, 2 being an SU(2) L index. Then the singlet and adjoint scattering channels are where M SU (2) C denotes the matrix element in Eq. (42) stripped from the SU(2) L structure. Hence, by considering the singlet channel in color space and the adjoint channel in SU(2) L space, we gain respectively a factor √ 3 and 2 in the partial wave.