Mapping out the Space for New Physics with Leptonic and Semileptonic $B_{(c)}$ Decays

Decays of $B$ mesons with leptons in the final state offer an interesting laboratory to search for possible effects of physics from beyond the Standard Model. In view of puzzling patterns in experimental data, the violation of lepton flavour universality is an interesting option. We present a strategy, utilising ratios of leptonic and semileptonic $B$ decays, where the elements $|V_{ub}|$ and $|V_{cb}|$ of the Cabibbo--Kobayashi--Maskawa (CKM) matrix cancel, to constrain the short-distance coefficients of (pseudo)-scalar, vector and tensor operator contributions. The individual branching ratios allow us then to extract also the CKM matrix elements, even in the presence of new-physics contributions. Bounds on unmeasured leptonic and semileptonic decays offer important additional constraints. In our comprehensive analysis, we give also predictions for decays which have not yet been measured in a variety of scenarios.


Introduction
Decays of B mesons caused by semileptonic b → u ν and b → c ν quark-level transitions provide valuable insights into the quark-flavour sector of the Standard Model (SM) that is described by the Cabibbo-Kobayashi-Maskawa (CKM) matrix. These channels play a key role to determine the CKM matrix elements |V ub | and |V cb | from experimental data, where we are facing extractions from inclusive and exclusive modes that are not fully consistent with one another [1]. Decays of this kind have also received a lot of attention in view of experimental data indicating violations of lepton flavour universality, which is a key feature of the SM (see, e.g., Refs. [2][3][4] and references therein). Is it actually violated in Nature? Analyses of New Physics (NP) effects in leptonic and semileptonic modes mediated by these transitions have been performed in, e.g., Refs. [5][6][7][8][9][10][11][12][13][14][15][16][17][18][19][20].
For NP analyses to constrain the corresponding short-distance parameters, care is needed concerning the values of |V ub | and |V cb | as their determinations may be affected by contributions from physics beyond the Standard Model. In view of this feature, we have proposed a strategy for the exclusive b → u ν modes in Ref. [21], introducing suitable ratios of leptonic and semileptonic modes of this decay class, where the CKM factor |V ub | cancels. Utilising these observables, we then determined allowed regions for the short-distance coefficients of (pseudo)-scalar operators, arising from NP effects. Using the individual branching ratios, we may finally extract |V ub |, while simultaneously allowing for physics from beyond the SM.
In this paper, we build upon this strategy, addressing two new points. The first concerns the analysis of the effects of new vector and tensor operators in the exclusive b → u ν decays: the B − → ν as well as B → π ν and B → ρ ν channels, thereby complementing our previous study where we analysed only (pseudo)-scalar operators. The second -and main -point is the development of a similar strategy for exclusive B decays originating from b → c ν processes, B → D ν and B → D ( * ) ν transitions as well as the leptonic B − c → ν modes. Here we consider the NP effects of vector, (pseudo)-scalar and tensor operators, and constrain the corresponding short-distance coefficients through the currently available experimental data. Moreover, we extract the CKM matrix element |V cb | with the help of the individual branching ratios, and make predictions for decays which have not yet been measured.
The theoretical framework of our study is given by the corresponding low-energy effective Hamiltonians, allowing for NP operators. We perform also the proper renormalisation group evolution to connect the short-distance coefficients at the high-energy scale at the 1 TeV level with the low-energy scale at the µ = m b level of the b-quark mass that is relevant for the analysis of experimental data for the B-meson decays under consideration.
The outline of this paper is as follows: in Section 2, we introduce the theoretical tools, including the general expressions for the exclusive decay rates, the renormalisation group evolution and the strategy for extracting the CKM matrix elements |V ub | and |V cb |. In Section 3, we focus on the exclusive b → u ν modes, and explore the space for NP effects arising from vector and tensor operators, as well as update our previous analysis of the (pseudo)-scalar operators. In Section 4, these results are applied to determine the allowed range for |V ub | and certain branching ratios that have not been measured. The exclusive b → c ν quark-level decays are then discussed in Section 5, giving a comprehensive picture of the constraints on the short-distance coefficients of the corresponding (pseudo)-scalar, vector and tensor operators following from the current experimental data. In Section 6, we present a detailed discussion of the picture for the CKM factor |V cb |, emerging from our NP analysis and the experimental information available. There we also consider the leptonic B c decays, which have so far not been observed experimentally. Finally, we summarise our conclusions in Section 7. In two appendices, we collect details of the specific form factor parameterisations applied in our numerical analysis, and the various input parameters that are used in our calculations.

Theoretical Framework
The effective Hamiltonian relevant to our studies is where q = u, c and = e, µ, τ . The corresponding operators are where P L,R = 1 2 (1 ∓ γ 5 ). The Wilson coefficients can be decomposed in their SM and NP contributions asC q, Potential NP contributions enter through the C q, a . In view of Eq. (4), we haveC q, a = C q, a for a = {V R , S, P, T }, i.e., they are fully characterized by their NP component. We will consider CP-conserving NP effects and hence only real Wilson coefficients. As discussed in [22][23][24], within the framework of the Standard Model Effective Field Theory (SMEFT), the coefficient C q, V L depends in general on the leptonic flavour. In contrast, at the lowest order in the v 2 /Λ 2 expansion, where v is the vacuum expectation value of the Higgs boson and Λ is the NP scale, the coefficient C q, V R turns out to be lepton-flavour universal. Therefore, unless indicated otherwise, we will assume To constrain the Wilson coefficients in Eq. (1), we will consider different leptonic and semileptonic B meson decays. In the presence of pseudoscalar and vector interactions, the leptonic modes obey the following expression: where Since the tensor matrix element vanishes, 0 |qσ µν P L b | B q = 0, leptonic B q decays do not help us to constrain the coefficient C q, T . They also receive no contribution from the scalar operator.
In the case of semileptonic B decays involving a vector mesonṼ in the final state, we have the following expression [25]:  In addition, for a semileptonic B decay into a pseudoscalar meson P , we have [25] dB(B → P −ν ) In Eqs. (8) and (9), q 2 represents the four momentum squared transferred to the leptonic system, composed of the andν . In addition, the magnitude of the three momemtum of the meson in the final state has been denoted as | p P | and | pṼ |. The expression of | pṼ | in terms of q 2 is given by The equation for | p P | is analogous to Eq. (10), with the straightforward replacement of MṼ by M P . Before continuing we would like to introduce some notation. In what follows, we will use the symbol to denote either an electron, a muon or a tau lepton, on the other hand will refer only to electrons and muons. On the experimental side it is not always transparent what the branching ratios of the semileptonic modes represent. For the purpose of determining |V qb | under the assumption of the SM, it is advantageous to combine information from decays with electrons and muons in the final state. Hence, they are often presented as B(B → V (P ) −ν ), where is an electron or a muon. The exact definition of these branching ratios is often ambiguous. For example, they could be averages over both leptons, or involve just one of the two. In order to test the SM and search for NP, it would be instrumental to have measurements separately for electrons and muons. It would then also be pertinent that these analyses do not rely on theoretical assumptions such as a form factor parametrization, and are presented independently for the B − andB 0 modes, which is currently not always the case. In this work, we will employ PDG averages whenever available, averaging their values for the B − andB 0 channels, and interpret the experimental data forB →Ṽ (P ) −ν as averages over electronic and muonic modes.
Our aim is to search for potential NP effects, entering through the Wilson coefficients of the operators in Eq. (2), as well as to determine the |V qb | while allowing for new contributions. We follow the strategy introduced in Ref. [21], which is discussed in more detail in Section 2.2.
Our analysis works under the general asumption of NP in τ leptons as well as electrons and muons. Therefore, in principle, we should fit for three different coefficients C e X , C µ X and C τ X , where X = P, S, V, T . However, to simplify the analysis, we will consider scenarios where the coefficients for the electron and muon are correlated according to where f e µ is a constant. In general, there are three possibilities which can be addressed in a model-independent way: C e X < C µ X , C e X ∼ C µ X and C µ X < C e X . For the purposes of illustration, we consider three possible scenarios: To evaluate the confidence regions allowed for our NP Wilson coefficients C X at the one-σ level we consider the following equation: where O is a given physical observable with central value O Experiment , and σ Theory and σ Experiment are the theoretical and experimental uncertainties respectively. The theoretical uncertainty is estimated based on where x i are the individual theoretical inputs on which O Theory depends. Notice that the theoretical uncertainty is in general not constant and depends on the particular value assigned to C X .

Introduction of New Physics Effects
In this work, we will consider NP contributions arising at a high energy scale of µ = 1 TeV. The evolution down to the bottom scale µ = m b is given by the following matrix [26]: Consequently, the short-distance coefficients at the bottom scale entering in the B-meson branching fractions are essentially functions of the coefficients at µ = 1 TeV. Notice that, as the result of the renormalization group evolution, in our basis the coefficient C q, T (1 TeV) yields contributions to C q, S and C q, P at the scale µ = m b . In the following discussion, we will determine allowed regions for the Wilson coefficients at the scale µ = 1 TeV, making the connection with the bottom scale implicitly through Eq. (15). From this point onwards and in order to simplify the notation, wherever we write C q, i (where i ∈ {V L , V R , S, P, T }), i.e., without specifying the scale µ, we are referring to the coefficients at µ = 1 TeV. There are a few exceptions, for which we will clearly indicate in the text that we are considering the scale µ = m b .
2.2 Strategy for the Determination of |V ub | and |V cb | Typically, the |V qb | are determined assuming no NP contributions (see for instance Ref. [27]). We have proposed a strategy to take into account potential NP effects arising from (pseudo)-scalar operators when determining |V ub | from exclusive b → u ν modes in Ref. [21]. In this section, we describe this strategy, generalizing it to b → q ν and |V qb | for all operators in Eq. (2).
There are analyses that have evaluated |V ub | and |V cb | considering simultaneously NP contributions in scalar, vector and tensor operators, see for instance [7,8,12,17]. These studies are based on simultaneous fits to |V qb | and the different NP Wilson coefficients, and take advantage of the correlation between these quantities established in branching fractions and differential distributions.
Our strategy follows a different approach, which allows us to extract NP effects independently of |V qb | for q = u, c. The key steps are as follows: 1. We first establish the possible regions for the different NP Wilson coefficients using ratios of branching fractions where |V qb | cancels out.
2. We then proceed with the evaluation of |V qb | by substituting the allowed NP regions established in the previous step inside specific observables sensitive to the different NP Wilson coefficients and |V qb |. We would like to stress that -since the NP contributions were already extracted in a |V qb | independent fashion -observables such as branching fractions only depend on one unknown parameter, which is precisely |V qb |. We can now finally extract this parameter, obtaining a value that takes potential NP contributions into account.
This procedure can, in certain cases, be further refined to account for correlations between the different observables and parameters considered during the analysis, as done in Ref. [21].

The b → u Transitions
The discussion in this section focuses on the determination of the b → u ν short-distance coefficients, and complements the analysis in Ref. [21]. There, the possibility of having NP effects from scalar and pseudoscalar interactions entering in different leptonic and semileptonic decays was explored. Here, we will consider in addition the vector and tensor structures, taking into account the renormalization group evolution as introduced in Section 2.1. For completeness, we will include the (pseudo)-scalar operators as well, switching on their coefficients at µ = 1 TeV, whereas the results in Ref. [21] relate to µ = m b . Our determinations will employ the following experimental measurements for B − leptonic decays [28,29]: For the corresponding process involving electrons, we will consider the upper bound obtained by the Belle collaboration in 2007 [30]: Within the SM, we have where we have used [27] |V ub | = (3.67 ± 0.15) × 10 −3 (21) and [31] f B − = 0.1900 ± 0.0013 GeV.
The value of |V ub | in Eq. (21) is determined by HFLAV from B → π ν data, where is an electron or a muon, with form factors from lattice QCD (LQCD) and light-cone sumrule (LCSR) calculations, assuming no NP in the light leptonic generations. This should be kept in mind when considering our evaluations for the branching fractions in the SM given in Eqs. (18), (19) and (20). As argued in Ref. [21], decoupling the determination of |V ub | from potential NP effects in b → u ν transitions requires a more careful approach which is developed further in this work. As a result, the values in Eqs. (18), (19) and (20) serve an illustrative purpose and should not be seen as predictions. Moreover, we would like to point out that there is a long-standing tension between exclusive, as in Eq. (21), and inclusive determinations of |V ub |. For the latter, HFLAV gives the value [27]: For the present analysis, we have also taken into account the semileptonic decay processes involving the ρ vector meson in the final state. From the results reported by the Belle collaboration [32], we obtain Here we have taken into account only the q 2 region satisfying q 2 ≤ 12 GeV 2 , as indicated by the subindex. The reason for this is that, on the theory side, we make use of form factors calculated using LCSRs, which are applicable at low values of q 2 [33]. Details on the form factors are provided in Appendix A.1. We interpret the experimental results in Eq. (24) as an average over electrons and muons, as discussed in Sec. 2. Using the isospin symmetry, we combine the results in Eq. (24) to obtain For semileptonic decays including a pseudo-scalar meson in the final state, we consider theB → π −ν processes, where = e, µ. On the experimental side, we have [29] B(B 0 → π + −ν ) [ = e,µ] = (1.50 ± 0.06) × 10 −4 , Making use of the isospin symmetry, we get In addition, we consider the following upper bound: as obtained by the Belle Collaboration [34]. In Appendix A.1 we provide details on the B → π form factors. When studying potential NP contributions to these decays, we have to take into account that they are actually used to determine |V ub | (see, e.g., Ref. [27]). Therefore, the value of |V ub | may be affected by NP, and we have to construct observables independent of |V ub | for consistency. To that end, we consider ratios of branching fractions to evaluate the possible values of the NP Wilson coefficients.
The |V ub |-independent leptonic ratios are where 1 , 2 are any of the leptons e, µ, τ . Additionally, we include the following combinations of leptonic and semileptonic branching fractions: with = e, µ, τ . Finally, we use also the ratio of semileptonic processes With the currently available branching ratios in Eqs. (16), (25) and (27), we have the following ratios at our disposal: where we use the notation defined in Eqs. (29)- (32). Given the number of independent branching fractions, at most three of these observables are independent. However, they are sensitive to different short-distance coefficients. In the following sections, we will consider each of the structures i ∈ {P, S, V L , V R , T } independently, constraining the coefficients C u,µ i and C u,τ i through three appropriate ratios of branching fractions. For the C u,e i coefficient, we will consider the three scenarios from Eq. (12). There are two exceptions to the approach discussed above: The Wilson coefficient C V R is lepton-flavour universal (see Eq. (5)), and in case of the C S short-distance coefficient, we have only sensitivity to C u,e S (1 TeV) and C u,µ S (1 TeV) through B(B → π −ν ) [ = e,µ] , so we will constrain these coefficients without assuming any relation between them.

Constraints on (Pseudo)-Scalar Wilson Coefficients
We will start our analysis of the b → u transitions by considering new contributions to the Wilson coefficients of the (pseudo)-scalar operators O u, P and O u, S introduced in Eq. (2). We have considered these operators before in Ref. [21]. The main difference with respect to our previous analysis is that here we will include RGE effects, considering the short-distance coefficients at the scale µ = 1 TeV. From Eq. (15), we can see that this involves a rescaling of the coefficients at µ = m b . In addition, we use more recent data on B(B − → µ −ν µ ) [28].

New Physics Entering through C u, P
The short-distance coefficient C u, P contributes to all branching ratios considered here with the exception of B(B → π −ν ) [ = e,µ] . Still we may consider ratios involving this decay mode, as it may be still be used to cancel |V ub |. Actually, all six observables in Eq. (33) are sensitive to the C u, P coefficient, but, as stated before, they are not all independent. We will consider the following three ratios: We have verified, in the scenario where C u,e P = C u,µ P , that a different subset of observables would not have a meaningful impact on the results. Interestingly, of the ratios in Eq. (34), the latter two do not involve a τ in the final state, so they constrain only C u,µ P , whereas R τ µ is sensitive to both C u,µ P and C u,τ P . As it turns out, there is no agreement between theory and experiment for the observable R e,µ ;ρ [q 2 ≤12] GeV 2 e,µ ;π in any of the three scenarios for C u,e P . We will address this in the next paragraph. First, we present in Fig. 1 the constraints on C u,µ P and C u,τ P from the other two observables listed in Eq. (34). We note that the results are very insensitive to the assumptions about C u,e P . In each scenario, we obtain four distinct regions where both constraints overlap, one of which includes the SM. under different assumptions about the relation between C u,e P and C u,µ P , as indicated in the plots.
Studying the ratio R e,µ ;ρ [q 2 ≤12] GeV 2 e,µ ;π further, we find that the theoretical value and experimental data only barely disagree. In fact, small variations in the uncertainty may have a significant impact in Fig. 1, yielding a contour that agrees with all four regions defined by the overlap of R τ µ and R µ e,µ ;ρ [q 2 ≤12] GeV 2 . Specifically, we find that, on the basis of this observable, we cannot exclude these solutions at the 1 σ level. µ] . Consequently, we cannot constrain C u,τ S and we will hence consider C u,e S and C u,µ S only. We will still determine the bounds from all three ratios involving the B(B → π −ν ) [ = e,µ] branching fraction; the distinction is that the cancellation of |V ub | is obtained from different modes in each ratio. Specifically, we take

New Physics Entering through
The constraints in the C u,µ S -C u,e S plane corresponding to the observables in Eq. (35) are shown in Fig. 2. Since they all overlap, we show them for clarity separately on the top row, and all together on the bottom. In case of the R µ e,µ ;π ratio, we note that there is no solution for the central value; a region only appears once the uncertainties are taken into account. Furthermore, the ratio R e,µ ;ρ [q 2 ≤12] GeV 2 e,µ ;π yields a doughnut-shaped constraint with a hole around the SM point. However, this is slightly over 1 σ, i.e., not a significant deviation.

Constraints on Vector Wilson Coefficients
We continue by considering the vector operators O u, V L and O u, V R . First we will allow for NP to enter through C u, V L , where is any of the charged leptons e, µ and τ . Next, we will consider a non-vanishing C u V R , which is lepton-flavour universal, as we discussed in Section 2.

New Physics Entering through
Allowing NP to enter only through C u, V L , the expressions in Eqs. (6), (8) and (9) take the following form: and Given that all these expressions have exactly the same dependence with respect to C u, V L , any ratio of these quantities involving the same lepton flavour will lead to the cancellation of the NP contributions. Consequently, to obtain constraints on the vector short-distance contributions, we have to include ratios of leptonic and semileptonic branching ratios with different leptonic content in numerator and denominator. Hence, from the set of observables in the list in Eq. (33) the only relevant ratios are Depending on the correlation between C u,e V L and C u,µ V L , some of the other observables in Eq. (33) may be considered as well. However, as indicated before, they are not independent. In addition, the observables in Eq. (39), with the semileptonic branching fractions in R τ e,µ ;ρ [q 2 ≤12] GeV 2 and R τ e,µ ;π replaced by B(B → ρµ −ν µ ) and B(B → πµ −ν µ ), respectively, would be the obvious candidates to constrain C u,µ V L and C u,τ V L in the future, should these measurements become available.
Let us start by making the assumption allowing us to constrain the independent coefficients C u,µ V L and C u,τ V L . The corresponding results are shown in Fig. 3, where the cross-shaped coloured areas illustrate the NP values for the Wilson coefficients that are allowed by the observables in Eq. (39). One of the four branches agrees with the SM, though the constraint from R τ e,µ ;ρ [q 2 ≤12] GeV 2 is slightly more than 1 σ away.
Unfortunately, it is not possible to use the regions in Fig. 3 to determine |V ub |, or to make predictions for branching ratios that have not yet been measured, such as B(B − → e −ν e ), since they are not further constrained. Although at some point the contours following from the three ratios may no longer overlap, other effects will start to play a role. Since |V ub | cancels in the ratios in Eq. (39), the region in Fig. 3 may correspond to values of |V ub | that are completely unrealistic when compared to other CKM constraints.
Let us illustrate this feature by considering the leptonic decays. In Fig. 4, we show again the allowed region in the C u,µ V L -C u,τ V L plane following from R τ µ . Moving along the dotted line, which indicates the central value, we obtain a relation between C u,µ V L and C u,τ V L . This information allows us to determine the corresponding correlation with |V ub | from any of the branching fractions. In this case, we employ the C τ The variation of |V ub |, moving along the central value of R τ µ , is added to Fig. 4. Here, we have restricted ourselves to a range of ten times the uncertainty around the central value in Eq. (21). We note that the value of |V ub | may differ significantly from that in Eq.   So far we have considered universality between C u,e V L and C u,µ V L . Let us continue with the other two scenarios from Eq. (12). The results can be found in Fig. 5. The regions agree with the SM, and in the case of C u,e V L = 10C u,µ V L , the regions are actually quite constrained in comparison with the other two scenarios considered here. In both plots in Fig. 5, we find agreement with the SM for one region, though at a bit more than 1 σ for R τ e,µ ;ρ [q 2 ≤12] GeV 2 , as was the case when making the assumption C u,e V L = C u,µ V L . Figure 5: Allowed regions in the C u,µ V L -C u,τ V L plane following from the observables defined in the text. The left panel corresponds to the assumption C u,e

New Physics Entering through
As indicated in Section 2, the short-distance coefficient C q V R is lepton-flavour universal. Switching on only this coefficient, we will determine the constraints from several observables. From Eq. (6), we find Considering that all the leptonic branching ratios have the same dependence on C u V R , taking ratios of purely leptonic processes cannot help us to constrain C u V R since its contribution will cancel. Nevertheless, we note that R τ µ should agree with its SM prediction of R τ µ = 1. From the experimental data, we obtain so this constraint is satisfied.
Let us now consider the other observables from the list in Eq. (33) that are also sensitive to C u V R : In Fig. 6, we plot the dependence of these ratios on C u V R , along with the region that corresponds to the experimental data. Comparing these contours in the range −15 ≤ C u V R ≤ 15, we find that each observable yields agreement with the SM value C u V R = 0 at the (1-2) σ level. In addition, each ratio allows for one or more intervals of NP values for C u V R , though not all of them overlap at 1 σ.

Constraints on Tensor Wilson Coefficients
As a final step in our b → u analysis, we consider the tensor operator. From Eq. (15) we can see that taking a non-zero vale for C T at the scale µ = 1 TeV yields contributions at the low energy scale µ = m b not only in C T , but also in the coefficients C P and C S . In consequence, even though the branching ratios of the leptonic modes do not depend on C u, T (m b ), they are sensitive to C u, T (1 TeV) through C u, P (m b ). Hence, all six ratios in Eq. (33) can be employed to constrain the C u, T coefficients at the high-energy scale. Since they are not all independent, we will consider again three ratios of branching fractions, to constrain the coefficients C u,µ T and C u,τ T at µ = 1 TeV. We verified, in the scenario where C u,e T = C u,µ T , that a different choice of observables would not meaningfully affect the results.
We find, similar to the pseudoscalar case, that there is no agreement between theory and experiment for R e,µ ;ρ [q 2 ≤12] GeV 2 e,µ ;π . The contours from the remaining two ratios are shown in Fig. 7. Each plot corresponds to a scenario from Eq. (12). Note that when we assume C u,e T = C u,µ T and C u,e T = 0.1C u,µ T , we obtain four regions where all contours overlap, one of which includes the SM. On the other hand, in the scenario where C u,e T = 10C u,µ T , the ratio R µ e,µ ;ρ [q 2 ≤12] GeV 2 does not yield a solution for the central value. Still, a region appears once the uncertainties are taken into account, yielding two solutions from the combination of both observables, one of which agrees with the SM. plane for the observables listed in the text. The three plots correspond to different assumptions about the relation between C u,e T and C u,µ T , as indicated.
As for the ratio R e,µ ;ρ [q 2 ≤12] GeV 2 e,µ ;π , the constraint is slightly stronger than in the case of C u, P . One the one hand, in the scenarios where C u,e T = C u,µ T and C u,e T = 0.1C u,µ T , we find that the constraint from this observable is quite unstable, as was the case when we considered the C u, P coefficients. In particular, at the 1 σ level, it is actually still in agreement with all four regions following from the overlap of the contours from R τ µ and R µ e,µ ;ρ [q 2 ≤12] GeV 2 . On the other hand, while the R e,µ ;ρ [q 2 ≤12] GeV 2 e,µ ;π ratio agrees with the SM at the 1 σ level in the scenario where C u,e T = 10C u,µ T , it does constrain the NP parameter space in this scenario.

Determination of |V ub | and Predictions
Now that we have determined allowed regions for the different short-distance coefficients in several scenarios, we may use these results to determine |V ub | and make predictions for branching ratios that have not yet been measured. We can consider any of the measured branching ratios to determine |V ub | while accounting simultaneously for NP. Then, we may use the results for |V ub | and the Wilson coefficients together to predict the branching fractions . Since experimental limits exist for two of these branching ratios, B(B − → e −ν e ) and B(B 0 → π + τ −ν τ ), we may find that (part of) our predicted range for these modes is already excluded. In that case, the bounds will be used to constrain the allowed parameter space even further.
The purpose of this section is to demonstrate our strategy for the determination of |V ub |, as well as the constraining power and potential enhancements of branching ratios such as B(B − → e −ν e ), which have not yet been measured. Therefore, we have chosen to limit ourselves to a selection of scenarios, which are: • NP through the pseudoscalar coefficients with the assumption C u,e P = C u,µ P , • NP through the left-handed vector coefficients with the assumption C u,e V L = 10C u,µ V L , • NP through the tensor coefficients with the assumption C u,e T = 0.1C u,µ T .

The Pseudoscalar Coefficients
The first step in our strategy is to obtain bounds on the short-distance coefficients C u,µ P and C u,τ P from the overlap of the contours from R τ µ and R µ e,µ ;ρ [q 2 ≤12] GeV 2 in the first plot in Fig. 1. As indicated in Sec. 3.1.1, the ratio R e,µ ;ρ [q 2 ≤12] GeV 2 e,µ ;π just about excludes the four solutions from the former two observables. However, as was discussed there, this constraint is quite unstable, so we will still apply our strategy using R τ µ and R µ e,µ ;ρ [q 2 ≤12] GeV 2 . The resulting allowed ranges are given in Table 1, where the regions are numbered 1 to 4 clockwise, starting at the SM-like solution in the top right corner of the first plot in Fig. 1. Table 2: Predictions for branching ratios that have not yet been measured corresponding to the allowed ranges of |V ub | and C u,µ P or C u,τ P , in the scenario where C u,e P = C u,µ P .
The next step is then to determine Here we scan the ranges for C u,µ P in Table 1 together with values for |V ub | to find agreement with the measurement in Eq. (25). The resulting allowed ranges for |V ub | are added in the fourth column of Table 1. Since the branching ratio B(B → ρ −ν ) [ = e,µ], q 2 ≤12 GeV 2 is very stable on such a small scale for C u,µ P , we find the same result from all four regions. The range contains the HFLAV result for |V ub | from exclusive decays in Eq. (21) at the upper end, but not |V ub | incl. as given in Eq. (23). Now that the short-distance coefficients and |V ub | have been determined, we can apply them to predict branching ratios that have not yet been measured. Of the three modes that were introduced for this purpose in the introduction of this section, two are sensitive to C u, It turns out that the branching ratio prediction of B(B − → e −ν e ) is larger than the experimental bound in Eq. (17) for both regions 3 and 4, thereby excluding these solutions. The allowed ranges for regions 1 and 2 are given in Table 2. As we found in Ref. [21], as well as in Ref. [35] for neutral leptonic B-meson decays, the C P coefficient potentially enhances the branching ratio of the B − → e −ν e decay by several orders of magnitude because it lifts the helicity suppression.

The Left-Handed Vector Coefficients
Let us now apply the same strategy to the Wilson coefficient C u, V L , specifically the scenario where we assume C u,e V L = 10C u,µ V L . The relevant contours are given in the left panel of Fig. 5. There are a few key differences with respect to the C u, P scenario: • The regions are now defined by the overlap of three different contours.
• Within the uncertainties, two solutions merge into one region, thereby giving two instead of four regions in total.
In this instance, we find that neither the upper bound on B(B − → e −ν e ) nor the one on B(B 0 → π + τ −ν τ ) gives any further constraints. The results for the short-distance coefficients and |V ub | are given in Table 3. Region 1 is the SM-like solution on the top in the left panel of Fig. 5, and region 2 is the other solution on the bottom. We find quite a large range for |V ub |, which includes the HFLAV value from exclusive modes in Eq. (21), as well as the inclusive result in Eq. (23). The large size of the range is dominated by the NP effects, i.e., the variation of |V ub | with respect to C u,µ V L . This is in contrast to the C u P scenario, where |V ub | was very stable over the range of the short-distance coefficient. On the other hand, this means that improving the determination of C u,µ V L has the potential to drastically reduce the allowed range for |V ub | as well.
The predictions for the branching ratios are given in Table 4. In contrast to the results for C P , the helicity supression of the leptonic decays is not lifted. Consequently, the potential enhancement of B(B − → e −ν e ) is not as dramatic as when NP would enter through the pseudoscalar coefficient. In Table 5, we give the allowed ranges of the short-distance coefficients and |V ub | that are in agreement with the upper bound on B(B 0 → π + τ −ν τ ). The corresponding predictions for the branching ratios are given in Table 6. We have numbered the regions again in a clockwise fashion, starting with 1 for the SM-like solution in the top right corner of the relevant plot in Fig. 7.

The Tensor Coefficients
Although we reach the upper bound of B(B 0 → π + τ −ν τ ) in all four regions, it mostly    Table 6: Allowed ranges for branching ratios that have not yet been measured corresponding to the allowed ranges for |V ub | and C u,µ T or C u,τ T , in the scenario where C u,e T = 0.1C u,µ T .
restricts solution 2 and 3. For these regions, the ranges of |V ub | are also quite small, and do not contain the HFLAV value from exclusive decays in Eq. (21), whereas the ranges found in regions 1 and 4 do. None of them reach the inclusive result in Eq. (23). In all four regions, the variation of |V ub | with respect to the short-distance coefficient has only a minor impact on the range that we obtain. As for the predictions of branching ratios that have not yet been measured, we find again potentially spectacular enhancements for B(B − → e −ν e ), though the maximum values are somewhat smaller than in the case of C u, P . It is not surprising to find such values, because the C u, T (1 TeV) yields a contribution to C u, P (m b ) as stated before, thereby also lifting the helicity supression of the leptonic modes.
Due to the potential enhancement of B(B − → e −ν e ) by several orders of magnitude with respect to the SM value, we have decided to summarize these results in Fig. 8. There we show the range of values larger than the SM that may be obtained for B(B − → e −ν e ) in the two allowed regions for C u, P and the four regions for C u, T . We also indicate the SM value from Eq. (20) and the experimental limit from Eq. (17). The results from the left-handed vector coefficient are not included because on the scale of this plot they would be essentially at the SM level.

The b → c Transitions
To continue our discussion, we focus on semileptonic processes induced by the b → c quark transition. More specifically, we will concentrate on the semileptonic decays of B-mesons B → D ν and B → D * ν with = e, µ, τ . Since the D and D * are pesudoscalar and vector meson states, respectively, this will allow us to probe different Lorenz structures involving b → c processes.
For the semileptonic decays involving the B → D transition, the following experimental averages calculated from independent measurements reported by the BaBar, Belle, CLEO and ALEPH collaborations [29] are available: For channels with a τ lepton in the final state, we have considered the processes where the first and the second ratios were obtained by the BaBar [36] and the Belle [37] experiments, respectively. On the other hand, for the B → D * ν decays we have the following averages which combine partial measurements by Belle, BaBar, DELPHI, CLEO, OPAL, ARGUS and ALEPH [29]: For the analogous processes involving a τ in the final state, we have the following experimental results available [29]: which are based on measurements performed by the Belle and BaBar Collaborations. We can average both determinations in Eq. (48), obtaining To exemplify the application of our strategy for the extraction of the CKM element |V cb | in the presence of possible NP effects, we need observables where |V cb | itself is not present. One possibility is to construct new quantities based on ratios of the branching fractions in Eqs. (45), (46), (47) and (48). Fortunately, ratios of this nature are already available directly from experimental determinations.
Indeed, the ratios satisfy the requirement of |V cb | independence. Furthermore, they posses the feature of showing a combined tension with the SM determination at the 3 σ level, as reported in Ref. [27]. The most recent averages include the latest measurement reported by Belle [38,39], and reads as follows [27]: R(D) = 0.340 ± 0.027 ± 0.013 = 0.340 ± 0.030, R(D * ) = 0.295 ± 0.011 ± 0.008 = 0.295 ± 0.014, being in tension with the SM at 3.08 σ, as reported in [27]. We would like to emphasize that R(D) and R(D * ) are defined in terms of physical quantities which depend on the CKM factor |V cb |. However, this factor actually cancels out in the ratios which define R(D) and R(D * ) themselves. This is precisely the desirable feature we need to decouple the determination of NP effects in the different Wilson coefficients and the independent extraction of |V cb |.
Our SM evaluations of the observables R(D) and R(D * ) are where we have performed the corresponding theoretical evaluations using the form factors for the B → D and B → D * transitions obtained from QCD sum-rule (QCDSR) and LQCD calculations. The technical details of the corresponding paramterizations are summarized in Appendix A.2.
In addition, different |V cb |-independent polarization observables can, in principle, serve our purposes as well [14]. Nevertheless, for most of them, their implementation suffers from large experimental uncertainties or from the lack of experimental information at all.
In this analysis, we consider the following polarization ratio which allows us to obtain useful constraints on NP coefficients associated with semileptonic processes involving a D * meson in the final state: Within the SM, we find which is compatible with the current experimental result [40] F L (D * ) = 0.60 ± 0.08 ± 0.04, at the 1.6 σ level, in agreement with other studies reported in the literature [41]. Finally, the Belle Collaboration [42] has performed an experimental test which quantifies the flavour universality between electrons and muons through the semileptonic ratio The corresponding experimental result reads [42] R e µ (D * ) = 1.01 ± 0.01 ± 0.03, which can be compared with our SM evaluation R e µ (D * )| SM = = 1.0045(1).
A striking feature of R e µ (D * ) is its small uncertainty on both the theoretical and the experimental side. This quantity will play an important role in restricting NP in the light generations of leptons.
For completeness, we present a set of semi-numerical formulae which display the dependence of the different observables in our study on the presence of the different NP contributions. To simplify the presentation, we show two cases separately depending on whether NP enters either through the τ leptons or through the light generations, i.e., muons and electrons. In addition, we don't give the uncertainty of the numerical coefficients, which arises from, e.g., the form factors. However, it is important to mention that we do take these uncertainties into account in our analysis. Finally, we would like to stress that, in order to compare with other sources available in the literature, we present results where the Wilson coefficients are evaluated at the scale µ = m b . The connection with the 1 TeV value can be established straightforwardly using the renormalization group evolution given in Eq. (15).
In the case of NP entering through the τ leptons, we have On the other hand, should NP affect the light generations, the effects on R D * are described as with For R(D), the corresponding contributions to the light generations are given by with For the light semileptonic ratio R e µ (D * ) introduced in Eq. (56), we have with In order to avoid any potential confusion, let us emphasize again that the short-distance coefficients in Eqs.

Constraints Pseudo-scalar Wilson Coefficients
In order to constrain the possible values of the pseudo-scalar Wilson coefficients for τ and light leptons e, µ, we make a combined scan which includes the observables R(D * ) R e µ (D * ) and and  Table 7: Bounds for C c,µ P and C c,τ P in different scenarios correlating C c,e P and C c,µ P .
and F L (D * ) introduced in Eq. (50) and Eq. (53), respectively. The corresponding results are presented in Fig. 9 from which we can read off the intervals allowed for C c,µ P and C c,τ P at 1 TeV. We summarize these results in the 2nd and the 4th columns of Table 7. Here we see how the regions for C c,τ P fall into any of the two intervals (−3.78, −2.32) or (0.48, 1.94) regardless of the assumed NP scenario correlating C c,e P with C c,µ P .  Figure 9: Constraints in the C c,τ P -C c,µ P plane considering different correlations between C c,e P and C c,µ P .
In the case of the light generations, the NP intervals depend on the scenario under consideration, see the 2nd column in Table 7, and can be relatively sizeable, reaching values as large as |C c,µ P | ∼ 3. We can, however, refine these bounds by including the light universality ratio R e µ (D * ). Unlike R(D * ) and F L (D * ), the observable R e µ (D * ) is only sensitive to potential NP in light generations and has the potential of further constraining C c,µ P by one order of magnitude in most of the cases. We show the effect of R e µ (D * ) in Figure 10, and present our numerical results in the 3rd column of Table 7. Interestingly, since C c,µ P and C c,τ P are correlated by R(D * ) and F L (D * ), we can return to these ratios with our new findings for C c,µ P to improve the available regions for C c,τ P . We present the corresponding results in the 5th column of Table 7. We have seen how the observable R e µ (D * ), introduced in Eq. (56), plays an important role in constraining the NP contributions to the light generations of leptons. This results from its high experimental precision, as given in Eq. (57). An interesting question is whether this observable can impose any restriction on the correlation factor f e µ which relates the short distance contributions for e and µ presented in Eq. (11). To answer this, we study the relationship between f e µ and |C c,µ P (m b )|, finding where the Wilson coefficients are evaluated at µ = m b . A careful analysis of the expression inside the square root shows that as long as the quantity R e µ (D * )/R e,SM µ (D * ) is compatible with one, there is no upper bound on |f e µ |. If we take R e µ (D * ) equal to the experimental value given in Eq. (57) we see that, within the current theoretical and experimental uncertainties, we have R e µ (D * ) = R e,SM µ (D * ). Hence, we conclude that, in spite of the precision in theory and experiment available for R e µ (D * ), it is not possible to restrict the values for f e µ . For completeness, we present graphically the relationship between f e µ and C c,µ P in Fig. 11.

Constraints on Scalar Wilson Coefficients
The scalar Wilson coefficients for semileptonic b → c transitions enter in our analysis only through the observable R(D). Even though R(D) correlates C c,τ S with C c,µ S , we find that this observable does not lead to strong bounds on these coefficients. We illustrate the interplay between C c,τ S and C c,µ S as imposed by R(D) in Fig. 12 for the different assumed scenarios between C c,e S and C c,µ S .

Constraints on Vector Wilson Coefficients
At first sight, the constraints on the left-handed Wilson coefficients can be derived using R(D), R(D * ), F L (D * ) and R e µ (D * ). The observables R(D) and R(D * ) are sensitive to the effects of τ as well as e and µ. On the other hand, F L (D * ) is only sensitive to NP in τ leptons in general. Interestingly, if the vector-left coefficient is the only effect entering in F L (D * ), then they cancel exactly since they have equal contributions in numerator and denominator, thereby rendering F L (D * ) insensitive to C c,τ V L . This can be verified using Eq. (59). Therefore, the limits on C c,τ V L can only be computed through R(D) and R(D * ), which as shown in Figure 13 cannot yield bounded regions in the C c,τ V L -C c,µ V L plane. Nevertheless, the inclusion of R e µ (D * ), which is sensitive to C c,e V L and C c,µ V L , improves this situation. Indeed, by providing bounds on C c,e V L and C c,µ V L as shown in Fig. 14, the observable R e µ (D * ) is also able to constrain C c,τ V L because the NP Wilson coefficients are not independent from each other but are correlated through R(D) and R(D * ), as can be seen in Fig. 13. Our findings are summarized in Table 8. As seen in Fig. 14 Table 8: Bounds for C c,µ V L and C c,τ V L in different scenarios correlating C c,e V L and C c,µ V L . For the scenario C c,e V L = C c,µ V L the allowed NP regions are not bounded and thus not presented. In order to study the behaviour of C V R , we consider the universality behaviour with respect to the different lepton flavours as established in Eq. (5). However, as seen in Fig. 15, the combination of observables included in this study cannot explain the current experimental data if C V R is the only NP contribution.

Constraints on Tensor Wilson Coefficients
The tensor Wilson coefficient enters in R(D), R(D * ), F L (D * ) and R e µ (D * ). In the case of F L (D * ), the tensor Wilson coefficient C c,τ T does not reproduce the experimental data within 1 σ. As for R(D) and R(D * ), the resulting regions are not bounded in the C c,τ T -C c,µ T plane. This can be seen in Figure 16. For completeness, we perform the analysis for R e µ (D * ) finding bounded regions for C c,µ T in Figure 17, we also present the corresponding numerical intervals in Table 9. However, our results for R(D) and R(D * ) already exclude the possibility of explaining data with a tensor interaction at the 1 σ level.  Table 9: Bounds for C c,µ T in different NP scenarios correlating it with C c,e T . The observable considered is R e µ (D * ). 6 Extraction of |V cb | and Implications for B(B c → −ν ) In this section, we proceed with the extraction of the values of |V cb | which are compatible with the NP contributions determined in the previous subsections. In analogy with the determination of |V ub | presented in Section 4, we use the bounds for the NP Wilson coefficients derived from the |V cb |-independent observables R(D), R(D * ), F L (D * ) and R e µ (D * ) to extract |V cb | through the interplay of an extra observable which is |V cb |dependent. To accomplish this goal and for the purposes of illustration, we will focus on the averaged branching fraction B(B → D * τ −ν τ ) which is sensitive to both the NP short-distance contributions and |V cb | itself. Moreover, this particular average is not affected by the ambiguity concerning the leptonic combination in the final state which we encountered while interpreting the experimental measurement of B(B → D * −ν ) , where the contributions from electrons and muons were not transparent.
The possible regions for the short-distance coefficients and |V cb | can be further refined by imposing upper bounds on the leptonic B c branching ratios as extra constraints, especially in scenarios with non-vanishing pseudoscalar Wilson coefficients at the bottom scale, as they lift the helicity suppression of the leptonic modes. Unfortunately, no direct constraints are available from the LHC. In view of this, an estimate of a bound on B(B c → τ −ν τ ) has been derived from LEP data at the Z peak in Ref. [43]. In Refs. [14,15], it was argued that this bound is too strict and values up to 60% cannot be excluded. Thus we will take Information on the lifetime τ Bc of the B c mesons can also be converted into bounds for the leptonic branching ratios [22,44], giving results in the same ballpark [14]. However, as this observable depends on any possible B c decay channel, the interpretation is more complex since other NP contributions may enter. In order to incorporate Eq. (67) in our studies, we use the theoretical expression withB and the Wilson coefficients are evaluated at µ = m b . For the purposes of numerical comparison, we make a small digression from our strategy and evaluate the leptonic branching fractions for the B c meson using the value of |V cb | obtained in [27] from exclusiveB → D −ν decays: Assuming this value of |V cb |, we may calculate the following "SM" leptonic branching fractions: We would like to stress that since the results in Eq. (71) to Eq. (73) use an external determination for |V cb |, they were not derived from our strategy and presented only as reference values for our future discussion. For completeness and for the purposes of comparison we also provide the inclusive value of |V cb | as presented in [45] |V cb | = 0.04162 +0.00026 −0.00080 .
We are now ready to discuss the extraction of |V cb | in the presence of NP contributions. In view of the results obtained in Section 5, we can only fulfill this task for the Wilson coefficients C c,τ P and C c,τ V L where the resulting NP regions are bounded. We describe both cases in the following subsections.

|V cb | Compatible with NP Pseudoscalar Coefficients
Using the values for C c,τ P given in Table 7 and the branching fraction B(B → D * τ −ν τ ) , we obtain the bounds for |V cb | shown in the second column of Table 10. Notice that we present our results for each of the possible correlations between C c,e P and C c,µ P . We observe that values for |V cb | in the range 0.032 ≤ |V cb | ≤ 0.042 are compatible with the experimental data considered so far. This seemingly large interval can be further reduced once the constraint in Eq. (67) is included. Indeed, this extra restriction has two effects. Firstly, it restricts the values of |V cb | to the subinterval 0.039 ≤ |V cb | ≤ 0.042. Secondly, it leads to the following region for the Wilson coefficient C c,τ P : 0.49 ≤ C c,τ P ≤ 0.60. The interplay between C c,τ P , |V cb | and B(B → D * τ −ν τ ) is shown in Fig. 18. We would like to stress that our values for |V cb | are compatible with the exclusive determination shown in Eq. (70) as well as with the inclusive value shown in Eq. (74). Interestingly, since the Wilson coefficients for tau leptons and muons are correlated through the observable R(D * ), the constraint imposed by the bound in Eq. (67) over C c,τ P has effects on the light generations as well. In the 5th column of Table 10 we show the new intervals for C c,τ P after the inclusion of B(B c → τ −ν τ ) in our analysis. In addition we recalculated the new intervals for C c,µ P correspondingly and present them in the 4th-column of the same table.
We would like to highlight that, as can be seen in Table 10, the interval for C c,τ P does not include the SM value C c,τ P = 0, which indicates that the experimental measurements involving τ leptons require NP effects to be addressed theoretically. In contrast, for muons C c,µ P = 0 is compatible with the current experimental results which in our case also implies C c,e P = 0 in view of the different scenarios that we are studying. Another interesting feature that distinguishes τ leptons from its lighter counterparts is the fact that our final intervals for C c,µ P and C c,e P translate into values for B(B c → µν µ ) and B(B c → eν e ) that are at most about 30%. To the best of our knowledge, this does not break any constraint analogous to the one in Eq. (67). Therefore, we can consider the values for B(B c → µν µ ) and B(B c → eν e ) obtained from our regions for C c,µ P and C c,e P as predictions, which we present in Table 11. Notice that the NP contributions can potentially lead possible cancellations which imply branching fractions below the Scenario |V cb | C c,µ P (1 TeV) C c,τ P (1 TeV)  Table 11: Predictions of B(B c → −ν ) in the presence of a NP Wilson coefficient C c, P for different scenarios correlating C c,e P and C c,µ P . The zero minima predicted for the different leptonic branching fractions is the consequence of a perfect cancellation between the NP contributions and the purely SM ones in Eq. (6). This implies that the effect of a pseudoscalar NP contributions can lead to values of the branching fractions of leptonic decays below their corresponding SM values.
corresponding SM values. We find that the masses of the light leptons can enhance the branching fractions B(B c → µν µ ) and B(B c → eν e ) by several orders of magnitude through NP effects entering in pseudoscalar couplings. This effect is most dramatic for electrons, where the enhancement can reach up to eight orders of magnitude. In the case of muons, the largest enhancement is about four orders of magnitude. This situation is analogous to the one explored in Ref. [35] for B leptonic decays. A previous model-dependent study, which accounts for enhancements in the branching fractions of the processes B c → ν , can be found in Ref. [5]. The dependence of B(B c → eν e ) and B(B c → µν µ ) on C c,e P and C c,µ P , respectively, is presented in Fig. 19 showing the maximal effects achievable. Moreover, in Fig. 20 we illustrate how the enhancements for B(B c → ν ), with = e, µ, compare to each other for each one of our NP scenarios.  Table 12: Values of |V cb | in different NP scenarios correlating C c,e V L and C c,µ V L . We do not present the scenario C c,e V L = C c,µ V L since in this case there are not bounded regions.  Figure 19: Dependence of B(B c → ν ) normalized to its SM value on C c,e P and C c,µ P for = e, µ.

|V cb | Compatible with NP Left-Handed Vector Coefficients
As discussed in Sec. 5.3, should NP enter through a left-handed vector interaction at 1 TeV, the available constraints given by the observables R(D), R(D * ) and R e µ (D * ) lead to four independent bounded regions for C c,e V L , C c,µ V L and C c,τ V L only when either C c,e V L = 10C c,µ V L or C c,e V L = 0.1C c,µ V L . However, once we include the bound for B(B c → τ −ν τ ) given in Eq. (67), we find that there are actually only two regions in the C c,τ V L -C c,µ V L plane compatible with data, which can be read in the 3rd and 4th columns of Table 12.
From our maximally constrained regions for C c,µ V L and C c,τ V L we proceed with the extraction of |V cb | using the CKM-dependent observable B(B → D * τ −ν τ ) . Just as for the pseudoscalar NP interaction, we do this by assessing the combinations of values for  Table 13: Predictions for B(B c → −ν ) in the presence of a NP Wilson coefficient C V L for different scenarios correlating C c,e V L and C c,µ V L . The scenario C c,e V L = C c,µ V L does not yield bounded NP regions and it is thus not presented.
C c,τ V L and |V cb | compatible with the measurement in Eq. (49). Our final result for |V cb | is then 0.038 ≤ |V cb | ≤ 0.043, whose correlation with C c,τ V L is shown in Fig. 21. Once more, we find that our determination for |V cb | is compatible with the exclusive and inclusive values presented in Eq. (70) and Eq. (74) respectively. To conclude our study, we analyze the interplay between C V L and B(B c → −ν ). As we have discussed at the beginning of this section, B(B c → τ −ν τ ) reduces the possible values that C c,µ V L and C c,τ V L can take by eliminating two out of the four regions allowed by the CKM-independent experimental ratios R(D), R(D * ) and R e µ (D * ). However, it is interesting to notice that the remaining solutions have C c,τ V L components which imply values for B(B c → τ −ν τ ) that are below the experimental bound in Eq. (67). Hence, we can interpret these effects as predictions for B(B c → τ −ν τ ). Something similar occurs with the projections on C c,µ V L with respect to B(B c → µ −ν µ ), and with the corresponding transformation to C c,e V L with respect to B(B c → e −ν e ). In contrast to the pseudoscalar case discussed in the previous section, the helicity enhancement does not occur with the vector interaction. In other words, unlike NP entering through pseudo-scalar interactions, the contribution of NP through vector interactions are suppressed by the mass of the lepton in the final state in the decay processes B(B c → −ν ) for = e, µ, τ . Consequently, our NP scenarios lead to values for B(B c → −ν ) that are close to their SM values, as shown in Table 13.

Conclusions
We have mapped out the space for NP effects in exclusive B decays originating from b → u ν and b → c ν quark-level transitions, developing and implementing strategies which utilise ratios of various decay observables to avoid input on the CKM matrix elements |V ub | and |V cb |. In the first case, considering B − → ν as well as B → π ν and B → ρ ν modes, we updated a previous analysis of NP effects arising from (pseudo)-scalar operators and complemented it with a study of vector and tensor operator contributions. The corresponding constraints for the Wilson coefficients are consistent with the SM at the (1-2) σ level. Since our approach attempts to investigate the effect of NP in all the leptonic generations, we have assumed three possible scenarios which correlate the short-distance contributions between electrons and muons. In particular, we explored specific cases which fall within the following categories: C e X C µ X , C µ X C e X and C µ X ∼ C e X . We have concluded our studies of the b → u ν modes by determining |V ub | from the branching ratio of the B → ρ ν channel, where = e, µ, while allowing simultaneously for NP in a variety of scenarios. When we considered NP entering through the pseudoscalar or tensor coefficients, we found quite similar results, obtaining 3.0 × 10 −3 ≤ |V ub | ≤ 3.7 × 10 −3 and 2.9 × 10 −3 ≤ |V ub | ≤ 3.7 × 10 −3 , respectively. Both these ranges include the exclusive determination from HFLAV, which assumes the SM, of |V ub | excl. = (3.67 ± 0.15) × 10 −3 but do not reach the inclusive value of |V ub | incl. = (4.32 ± 0.12 +0. 12 −0.13 ) × 10 −3 [27]. Furthermore, we found in both cases that the value of |V ub | only shows a minor variation with respect to the short-distance coefficient on the relevant scales. On the other hand, we found that the variation with respect to the Wilson coefficient was a lot stronger in the case of new contributions to the left-handed vector coefficients, yielding a larger range of 2.8×10 −3 ≤ |V ub | ≤ 5.8×10 −3 , which includes both the exclusive and inclusive determination from HFLAV. We then applied our results for the short-distance coefficients and |V ub | to make predictions for the branching ratios of the B − → eν e , B → πτν τ and B → ρτν τ decay modes, which have not yet been measured. For the first two, experimental limits are available, which excluded (part of) the space for the Wilson coefficients and |V ub |, which was already taken into account in the results for |V ub | presented above. The largest effects were found in the tensor and pseudoscalar NP scenarios. For the former, the branching ratio of B → πτν τ severely constrains the allowed regions. As a result, the current data allow for values of this observable up to the experimental limit B(B 0 → π + τ −ν τ ) < 2.5 × 10 −4 [34]. Additionally, the branching ratio of B − → eν e may reach values of up to 2.3 × 10 −8 , being enhanced with respect to its SM value by three orders of magnitude. NP entering through the pseudoscalar operator yields even larger effects for the leptonic mode, with a potential enhancement by four orders of magnitude with respect to the SM, corresponding to B(B − → e −ν e ) = 1.4 × 10 −7 . We have extended this framework to b → c ν modes, considering B → D ν and B → D ( * ) ν transitions. Moreover, we had a closer look at the leptonic B − c → ν channels which provide very interesting probes for NP effects. We have provided a stateof-the-art discussion of the SM predictions of the key observables, including the R D ( * ) ratios, and the constraints on contributions from (pseudo)-scalar, vector and tensor operators. Useful constraints arise also from the polarisation observable F L (D * ) and the ratio R e µ (D * ), putting an interesting experimental limit on the flavour universality between electrons and muons. Just as for the b → u transitions, we have performed our studies within the context of three scenarios correlating the Wilson coefficients for electrons and muons. We found that the experimental results for our different observables can be explained through pseudoscalar and left-handed vector NP effects. Following a procedure analogous to the one applied in the case of the b → u transitions, we have extracted |V cb | using the CKM-dependent branching fraction B(B → D * τ −ν τ ) , finding values within the range 0.038 ≤ |V cb | ≤ 0.043 which are compatible with the exclusive [27] |V cb | = 0.03958±0.00117 and the inclusive [45] |V cb | = 0.04162 +0.00026 −0.00080 determinations. In addition, we have determined potential NP enhancements of the branching fractions for the leptonic decays B − c → µ −ν µ and B − c → e −ν e . Interestingly, if NP is present in pseudo-scalar Wilson coefficients, the branching fraction B(B − c → e −ν e ) can be enhanced by up to eight orders of magnitude in our scenarios reaching B(B − c → e −ν e ) = 0.31. In the case of B(B − c → µ −ν µ ), we may obtain an increase as large as four orders of magnitude leading to B(B − c → µ −ν µ ) = 0.39. In the future, it will be especially useful to have more measurements of the semileptonic decay modes with light leptons in the final state available, separately for electrons and muons, as well as (further) searches for decays that have so far not been measured. We look forward to seeing how the picture for the NP constraints following from leptonic and semileptonic B (c) decays and the methods proposed in this paper will evolve, exploiting the wealth of experimental data to become available. and for the B → P transitions. Note that the | pṼ | and | p P | are functions of q 2 . The seven form factors V , A 0 , A 1 , A 2 , T 1 , T 2 and T 3 for B →Ṽ and the three form factors F 0 , F 1 and F T for B → P have to be determined through non-perturbative techniques. In Appendices A.1 and A.2, we discuss these form factors for B → {ρ, π} and B → D ( * ) modes, respectively.

A.1 B → {ρ, π} Form Factors
For the calculation of the branching ratio of theB → ρ −ν decay, we require the form factors V , A 0 , A 1 , A 2 , T 1 , T 2 and T 3 . As we did in Ref. [21], we make use of the LCSR calculation from Ref. [33]. The form factors are given by where with The factor k (ρ 0 ,u) = f takes into account that we have a b → u instead of a b → d transition.
The form factors can now be determined using the coefficients α i k in Table 14 and the mass terms m R,i in Table 15. We do not consider the correlations between the α i k coefficients in our evaluations. The form factors are valid in the range 0 ≤ q 2 ≤ 14 GeV 2 . We note that Ref. [33] considers the form factors A 1 , A 12 instead of A 1 , A 2 , and T 2 , T 23 instead of T 2 , T 3 . We can convert the form factors using the following expressions [33]:  Table 14: Coefficients for the B → ρ form factors [33].
The z π and t π 0 are defined in analogy to z ρ and t ρ 0 with m ρ replaced by m π . The required parameters b i n are given in Table 16. We do not consider correlations between these parameters. These form factors are valid over the full kinematic range.

A.2 B → D ( * ) Form Factors
In order to determine the different hadronic form factors required for the calculation of the B(B → D ν) and B(B → D * ν) branching ratios, we use the parameterization provided in Ref. [48] that includes O(Λ QCD /m c,b ) and O(α s ) contributions.
In the case of the B → D transition, the connection between the current form factors  Table 16: Coefficients for the B → π form factors [46,47]. F 1 , F 0 , F T and the HQET hadronic form factors h + , h − and h T is given by where w depends on the four momentum transfer q 2 to the system composed by the andν , so that it can recoil against the D mesons according to Of particular interest is the zero-recoil point, corresponding to the maximum fourmomentum transferred to the lepton pair, for which Eq. (96) trivially reduces to w 0 ≡ w(q 2 0 ) = 1.
Notice that in Eqs. (101) and (102), the power corrections in m b and m c are included through the terms multiplying The required inputs in Eq.
and as indicated in [48], the prescription has to be used everywhere except in those terms where theΛ/m b,c factors do not multiply subleading Isgur-Wise functions in Eq. (101) and Eq. (102) to ensure the cancellation of leading renormalons. It can be seen how at leading order in α s and Λ QCD /m b,c the different HQET form factors either reduce to a common Isgur-Wise function ξ or vanish.
Since we want to explore the possibility of having NP effects in light leptons, we consider the set of input parameters shown in Table 17 and present the corresponding correlation matrix in Table 18. They were obtained from a fit to QCDSR and LQCD calculations as presented in [49]. This set of inputs avoids the usage of the differential distributions for the processes B → D ( * ) ν , for = e, µ, to fit the data.

B Input Parameters
For convenience, we summarize the numerical values and sources of the input parameters used in this paper in Table 19.