The Exclusive Vision of Rare $K$ and $B$ Decays and of the Quark Mixing in the Standard Model

The most common predictions for rare $K$ and $B$ decay branching ratios in the Standard Model are based on the CKM elements $|V_{cb}|$ and $|V_{ub}|$ resulting from global fits, that are in the ballpark of their inclusive and exclusive determinations, respectively. In the present paper we follow another route. We assume that the future true values of $|V_{cb}|$ and $|V_{ub}|$ will be both from exclusive determinations and set them equal to the most recent ones from FLAG. An unusual pattern of SM predictions results from this study with some existing tensions being dwarfed and new tensions being born. In particular using the HPQCD $B^0_{s,d}-\bar B^0_{s,d}$ hadronic matrix elements a $3.1\sigma$ tension in $\Delta M_s$ independently of $\gamma$ is found. For $60^\circ\le\gamma\le 75^\circ$ the tension in $\Delta M_d$ between $4.0\sigma$ and $1.1\sigma$ is found and in the case of $\epsilon_K$ between $5.2\sigma$ and $2.1\sigma$. Moreover, the room for new physics in $K^+\to\pi^+\nu\bar\nu$, $K_L\to\pi^0\nu\bar\nu$ and $B\to K(K^*)\nu\bar\nu$ decays is significantly increased. We compare the results in this EXCLUSIVE scenario with the HYBRID one in which $|V_{cb}|$ in the former scenario is replaced by the most recent inclusive $|V_{cb}|$ and present the dependence of all observables considered by us in both scenarios as functions of $\gamma$. We also compare the determination of $|V_{cb}|$ from $\Delta M_s$, $\Delta M_d$, $\epsilon_K$ and $S_{\psi K_S}$ using $B^0_{s,d}-\bar B^0_{s,d}$ hadronic matrix elements from LQCD with $2+1+1$ flavours, $2+1$ flavours and their average. Only for the $2+1+1$ case values for $\beta$ and $\gamma$ exist for which the same value of $|V_{cb}|$ is found: $|V_{cb}|=42.6(4)\times 10^{-3}$. This in turn implies a $2.7\sigma$ anomaly in $B_s\to\mu^+\mu^-$.


Introduction
The rare K and B decays and the quark mixing being GIM [1] suppressed in the Standard Model (SM) and simultaneously being often theoretically clean are very powerful tools for the search of New Physics (NP) [2]. Unfortunately the persistent tension between inclusive and exclusive determinations of |V cb | (see e.g. [3][4][5][6]) weakens this power significantly. As recently reemphasized by us [7] this is in particular the case of the branching ratios for rare K-meson decays and the parameter ε K that exhibit stronger |V cb | dependences than rare B decay branching ratios and the ∆M s,d mass differences. Also similar tensions in the determination of |V ub | [8] matter.
One possible solution to cope with this difficulty is to consider within the SM suitable ratios of two properly chosen observables so that the dependences on |V cb | and |V ub | are eliminated [7,9,10]. While in [9,10] B physics observables were considered, the analysis in [7] was dominated by the K system and its correlation with rare B decays and B 0 s,d −B 0 s,d mixing. In this manner we could construct 16 |V cb |-independent ratios that were either independent of the CKM parameters or only dependent on the angles β and γ, that can be determined in tree-level processes. Having one day precise experimental values for the ratios in question and also precise values on β and γ will hopefully allow one to identify particular pattern of deviations from SM expectations independently of |V cb | pointing towards a particular extension of the SM.
But these ratios, even if useful in the context of the tensions in question, are not as interesting as the observables themselves. Therefore, assuming in addition no NP in ε K , ∆M d and ∆M s and in the mixing induced CP-asymmetry S ψK S , these ratios allowed to obtain |V cb |-independent SM predictions for a number of branching ratios [7]. As these four quark mixing observables are very precisely measured and theoretically rather clean, the resulting SM predictions obtained in this manner turned out to be the most precise to date. A brief summary of the results of this analysis just appeared [11].
Another insight in this problematic has been provided recently by the authors of [12] who made a determination of |V cb | and |V ub | from loop processes, rare decays and quark mixing, by assuming no NP contributions to these observables. To this end they could use only well measured observables in the B system and ε K . This strategy has already been explored in [13] but there only ε K , ∆M d and ∆M s and S ψK S have been considered.
There is no question about that the analyses in [7,9,10] will help us to identify possible departures from SM predictions for the |V cb |-independent ratios and possible pattern of |V cb | determinations from various loop processes as analysed in [12,13], but also to some extent in [7]. See in particular Figs. 12 and 14 of the latter paper. Yet, eventually the most obvious procedure to look for NP is to determine all CKM parameters in tree-level processes under the assumption that NP contributions to these decays are negligible. This assumption is more likely to be correct than assuming no NP contributions in loop induced decays. Subsequently the resulting values of the CKM parameters inserted into SM amplitudes for loop induced processes would allow for definite predictions for GIM suppressed observables.
In this spirit in the present paper we follow a more direct but a novel route, which to our knowledge has not been explored by anybody to date, at least as far as SM predictions for theoretically clean observables like rare K and B decays, ∆M s,d and the mixing induced CP asymmetry S ψK S are concerned. Instead of using the values of |V cb | and |V ub | resulting from the global fits of the CKM matrix [14,15], we assume, in contrast to the prevailing inclusive expectations in the case of |V cb |, that the future values of both |V cb | and |V ub | will 2 be determined from exclusive tree-level decays. Therefore we use FLAG [6] averages, 1 which are based on a number of LQCD calculations that are listed after (2). With the precisely known |V us | the resulting rare decay branching ratios, ε K , ∆M d , ∆M s and S ψK S depend than only on the angles β and γ in the unitarity triangle that are moreover correlated through the CKM unitarity. An unusual pattern of SM predictions results from this study. Some present tensions are dwarfed and new tensions are born. In some cases their sizes depend sensitively on the value of γ which enhances the importance of precise measurements of this parameter, stressed in particular in [7,16].
The view that exclusive decays will eventually lead to the best determination of |V cb | is rather unusual but has been already expressed by the first author in the past [2]. The point is that precise measurements of formfactors by Lattice QCD (LQCD) accompanied by improved measurements of the relevant branching ratios should allow eventually a better control over theoretical uncertainties than it is possible in inclusive decays and consequently determinations of |V cb | and |V ub | that do not rely on quark-hadron duality. Yet, to be on the safe side, in view of the important progress in the inclusive determination of |V cb | [3][4][5][6], we compare at all stages the results in this EXCLUSIVE scenario with the HYBRID one in which |V cb | in the former scenario is replaced by the most recent inclusive |V cb | from [4].
To our knowledge in the literature only the authors of [17] performed a similar study, but only for ε K , finding, similar to us, a significant deviation of the SM prediction from the data. However, their analysis differs from ours in that for the CKM parameters they used the values obtained from global fits of the UT which can be questioned because in fact these SM global analyses used already ε K in their fits. Moreover, until now they did not incorporate the theoretical advances in ε K from [18] which have been taken by us into account in [7] and also in the present analysis.
The outline of our paper is as follows. In Section 2 we set up our strategy as far as CKM parameters are concerned. We also list the input parameters used in our numerical analysis. However, we refrain from the expressions for the observables which have been studied already by us in [7] and are collected there and in [2]. The numerical analysis is presented in Section 3 with the SM predictions for many observables resulting from the EXCLUSIVE LQCD scenario and from the HYBRID scenario. In Section 4 we calculate the impact of the hadronic matrix elements with 2 + 1 + 1 flavours from the HPQCD collaboration [19] on our results for rare B decays in [7], where the averages of HPQCD results and 2 + 1 results from Fermilab Lattice and MILC Collaborations (FNAL/MILC) [20] calculated in [10] have been used. We also illustrate how the determination of |V cb | from ∆M s , ∆M d and ε K depends on the number of flavours used in LQCD calculations of the relevant hadronic matrix elements.
We conclude in Section 5. Two short appendices list LQCD results for F Bq and F Bq B Bq for N f = 2 + 1 and N f = 2 + 1 + 1.

Strategy
The CKM parameters entering our analysis will be with γ one of the angles in the UT, shown in Fig. 1. It is equal, within an excellent accuracy, to the single phase in the standard parametrization of the CKM matrix [21,22].
As the input parameters we will use λ = 0.225 and the FLAG values for |V cb | and |V ub | extracted from the exclusive tree-level decays. Now these values, as given in the latest FLAG's report, read [6] |V cb | = 39.48 (68) These results are based on a number of different LQCD analyses as summarized in Fig. 38 of [6]. These are from FNAL/MILC [23][24][25][26], HPQCD [27][28][29] and RBC/UKQCD [30] with further details given in the original papers and [6]. However, the value for |V cb | in (2) does not include the most recent one from Fermilab/MILC [31] that is significantly lower 38.40(74) × 10 −3 . Fortunately, we were able to obtain from FLAG a preliminary result for |V cb | 2 that includes the latter result. Our basic values for |V cb | and |V ub | obtained from the overall 2022 FLAG's (|V cb |,|V ub |) fit will be then as follows 3 with |V ub | practically unchanged. Larger values for |V cb | from exclusive decays using LQCD, in the ballpark of 41.0 × 10 −3 , have been reported in [32][33][34][35] and we are looking forward to the 2023 FLAG report incorporating these results. We will also compare our EXCLUSIVE scenario with the HYBRID one in which the value for |V cb | is the inclusive one from [4] and the exclusive one for |V ub | as above: For γ we will use a broad range 60 • ≤ γ ≤ 75 • . Using CKM unitarity the angle β in the UT can then be determined through the correlation of β and γ On the other hand the sides of the UT, R t and R b , can be solely expressed in terms of the angles β and γ, as follows [36] We observe that R t depends dominantly on γ, while R b on β. These approximations follow from the experimental fact that β + γ ≈ 90 • and it is an excellent approximation to set sin(β + γ) = 1 in the formulae below although we will not do it in the numerical evaluations. The values of |V cb | and |V ub | in (3) imply then and using [37] that both differ mildly from the measured values [22] On the other hand in the HYBRID scenario we find and β = (21.74 ± 0.84) • , S ψK S = sin(2β) = 0.688 ± 0.022, (HYBRID) (12) in perfect agreement with the experimental measurements in (10). This brief exercise is an overture to the new tensions emerging from the exclusive strategy. In Fig. 2 we show β as a function of γ for the values of |V ub |/|V cb | in two scenarios in question and compare them to the one-sigma range for β in (10). We observe that in the case of the EXCLUSIVE strategy there is indeed a mild tension. To remove this tension a negative NP phase has to be added to β in the formula for S ψK S in (9). See (23). This negative phase originates in a NP phase in B 0 d −B 0 d mixing. Significantly larger tensions will be found in most observables analyzed by us. Figure 3: |V td | and Imλ t as functions of γ in the EXCLUSIVE and HYBRID scenarios.
In this context we would like to mention the analysis in [38] in which the ratio |V ub |/|V cb | was proposed as a useful test of the SM because of reduced hadronic uncertainties combined with the fact that this ratio is almost the same for the exclusive and inclusive determinations of |V cb | and |V ub |.
Finally, useful are also the following expressions (λ t = V td V * ts ) where the values in (3) have been used. For the HYBRID scenario we have Moreover As |V td | and Imλ t play an important role in rare K and B decays we show in Fig. 3 their dependence on γ in both scenarios.
A recent review of tree-level determinations of β and γ can be found in Chapter 8 of [2]. See also [39,40]. But here we will use γ as a free parameter and β as an output by means of (5).  [22] m K 0 = 497.61(1)MeV [22] S ψK S = 0.699 (17) [22] Table 1: Values of the experimental and theoretical quantities used as input parameters. For future updates see FLAG [6], PDG [22] and HFLAV [42].

Numerical Analysis
Our numerical analysis uses the formulae for various branching ratios that we have collected in [7]. The parameters, other than the CKM ones, entering the formulae in [2,7] are collected in the Table 1. Except for the values of F Bs B s and F B d B d that are taken this time from the HPQCD collaboration [19] 4 , other parameters are unchanged. These two inputs from N f = 2 + 1 + 1 LQCD calculations are only slightly lower than the ones used in [2,7] but are significantly lower than the ones from N f = 2 + 1 LQCD average given by FLAG in [6]. These differences are summarized in Appendix B. Their impact on the determination of |V cb | from ∆M s and ∆M d will be analysed in Section 4.
Here we only recall that the dependence of the observables considered by us on |V us | is negligible. As far as β and γ are concerned, the angle β is already known from the mixing induced CP-asymmetry S ψK S with respectable precision as given in (10) and there is a significant progress by the LHCb collaboration on the determination of γ from treelevel strategies [37] so that we have presently from tree-level decays the value given in (8). Moreover, in the coming years the determination of γ by the LHCb and Belle II collaborations should be significantly improved so that precision tests of the SM using the strategy in [7] and the one presented here will be possible.
However we emphasize that we do not use the value of γ above as an input parameter. Our strategy will be to treat γ as a free parameter in the rather broad range 60 • ≤ γ ≤ 75 • and in view of the future measurements of γ by LHCb and Belle II to exhibit the γ dependence of the observables considered by us. We recall that the angle β is the output by means of the unitarity relation (5) as given in Fig. 2.
In Table 2 we show SM predictions for a number of rare K and B branching ratios and ∆F = 2 observables resulting from the EXCLUSIVE input in (3) setting γ = 65.4 • , the central LHCb value. The uncertainties appearing therein, thus, do not include any error on the γ determination. We also show our results in the HYBRID scenario defined in (4). The latter are not far from the ones obtained in [7] where the absence of NP contributions to   (4) and in the fourth the experimental data. ∆F = 2 observables was assumed. We do not consider the decays like B → K(K * ) + − that have larger theoretical uncertainties than the observables considered by us. Their |V cb | dependence has been investigated recently in [12]. The decay B s → X s γ was not considered in [7]. The result for B s → X s γ in both scenarios is obtained here from [48] that effectively corresponds to the inclusive |V cb | = 42.0 × 10 −3 . We just rescaled it using the exclusive and inclusive values of |V cb | for EXCLUSIVE and HYBRID scenarios, respectively.
In Figs.4, 5 and 6 we show the γ dependence of the following observables in both scenarios for |V cb | and |V ub |. B(K S → µ + µ − ) has the same γ dependence as B(K L → π 0 νν) and the remaining observables do not depend or only weakly depend on γ.
Concentrating first on the EXCLUSIVE scenario we observe: • As seen in Table 2 for γ = 65.4 • , the largest tensions are found for ε K (4.1σ), ∆M s (3.1σ) and ∆M d (2.6σ).
• As ∆M s is practically independent of γ this tension remains for other values of γ.   Figure 6: ε K and the CP asymmetry S ψK S as functions of γ, in the EXCLUSIVE and HYBRID scenarios. The bands represent the uncertainties related to |V cb |, |V ub |, |V us | and to the non-CKM parameters. The red band in the upper panel represents the experimental value for ε K , with its 1σ uncertainty. The same for S ψK S .
• As seen in Fig. 4, the branching ratios for K + → π + νν and K L → π 0 νν are significantly suppressed below the values found in the literature that are in the ballpark of 8.5×10 −11 and 3.0 × 10 −11 , respectively [7]. But this suppression decreases with increasing γ.
• As seen in Fig. 5, the tension in ∆M d decreases to 0.6 σ for γ = 75 • but is as large as 4σ for γ = 60 • . The branching ratio for B d → µ + µ − shows a similar behaviour because its ratio to ∆M d is CKM parameters independent. The uncertainty in ∆M d is a bit larger because of the additional hadronic uncertainty in the parameterB d . The red band in the upper panel represents the experimental value for ∆M d , with its 1σ uncertainty.
• As seen in Fig. 6, the tension for ε K (with the experimental measurement shown in red) is practically linear in γ and in the range of γ considered varies from 2.0σ for γ = 75 • to 5.2σ for γ = 60 • . While significant tension in ε K in the EXCLUSIVE scenario has been already identified in [17], our analysis differs in several respects from that paper as we already stated at the beginning of this writing. On the other hand the tension for S ψK S is practically independent of γ and in the ballpark of 1.0σ so that in this case one really cannot talk about an anomaly.
On the other hand in the HYBRID scenario all these tensions disappear but the one in B s → µ + µ − is independently of γ in the ballpark of 2.1σ [10].  [7] for different values of γ within the SM. The SM area corresponds to the HP QCD 2 + 1 + 1 input (left panel) and to the 2 + 1 + 1 and 2 + 1 average (right panel) used in [7]. The green area represents the EXCLUSIVE scenario and the gray area the present experimental situation.

The Impact of the HPQCD Results
It is of interest to see how the use of 2 + 1 + 1 hadronic matrix elements from the HPQCD collaboration [19] used in the present paper, instead of the ones used in [7] (the average of 2+1 and 2+1+1 results), would modify our results for rare B decays of the latter paper in which no NP in ∆M s and ∆M d has been assumed. We make this comparison in Table 3. For completeness we list there also results for rare K decays which remain unchanged. This also allows the comparison with the results obtained in EXCLUSIVE and HYBRID scenarios in Table 2.
We observe that all B decays branching ratios in Table 3 are larger than our results obtained in [7] that were rather close to the ones in the HYBRID scenario. In view of still sizable experimental errors this impact of the HPQCD results cannot be fully appreciated with the exception of B s → µ + µ − . With the branching ratio for this decay in Table 3 and the experimental data in Table 2, assuming no NP in ∆M s , the anomaly in B s → µ + µ − of 2.1σ found in [10] is raised to 2.7σ.
It should be emphasized at this point that all the correlations found in [7] that do not involve ∆M s and ∆M d remain unchanged but the predictions for B decay branching ratios change as we have just seen. This then implies a different SM region in the correlation between K + → π + νν and B s → µ + µ − that we illustrate in Fig. 7. There the result using HPQCD 2 + 1 + 1 input (left panel) is compared with the one of [7] (right panel) where the average of 2 + 1 + 1 and 2 + 1 matrix elements from [10] has been used. This difference shows the importance of charm contribution in LQCD calculations. Note that the result for K + → π + νν did not change relative to [7].
In fact our results for B s,d → µ + µ − are rather consistent with the ones obtained by the HPQCD collaboration [19] in 2019. But since then the experimental accuracy of B s,d → µ + µ − significantly increased [52][53][54] allowing a better estimate of the anomaly in question.
Finally, we would like to address another important issue. In [7], using the averages of B 0 s,d −B 0 s,d hadronic matrix elements from LQCD calculations with 2 + 1 and 2 + 1 + 1 flavours, as given in (32), we found that there are no values of β and γ for which the same value for |V cb | can be obtained from ε K , ∆M d and ∆M s when imposing the experimental constraint from S ψK S . It is then of interest to investigate what happens when this analysis is repeated separately for the 2 + 1 and 2 + 1 + 1 matrix elements as given in (26) and (29), respectively.
The result of this exercise is shown in Fig. 8. We observe that only in the case of 2 + 1 + 1 flavours consistent result for |V cb | from all observables considered by us is obtained which in turn provides unique values of |V cb | and γ. The determination of γ and |V cb | can be further improved by considering first the |V cb |-independent ratio ∆M d /∆M s from which one derives an accurate formula for sin γ with the value for ξ from HPQCD [19] and where (17) has been used. The advantage of using this ratio for the determination of γ over studying ∆M s and ∆M d separately is its |V cb |-independence and the reduced error on ξ from LQCD relative to the individual errors of hadronic parameters in ∆M s and ∆M d . See Appendix B.
Combining then ε K , ∆M d , ∆M s , S ψK S and using (19) we obtain finally the following values of the CKM parameters The value of |V cb | is somewhat larger than in the HYBRID scenario in (4) but consistent with it. It should be noted that the determination of γ in this manner is more accurate than its present determination from tree-level decays in (8). The corresponding value of |V ub | is which is slightly larger than the FLAG determination but consistent with it. We observe that ε K dominates this determination of |V cb |. This can be traced back to its larger senitivity to |V cb | than it is the case for ∆M s,d . While ∆M s,d are proportional to |V cb | 2 , |ε K | exhibits approximately |V cb | 3.4 dependence [7]. Reducing the error on β, represented by the green band, and decreasing the error on γ from its tree-level measurements will provide a determination of |V cb | with an error below 1%. The 2 + 1 case demonstrates significant inconsistencies between |V cb | values from ∆M d,s and ε K . The average in (32) considered by us in [7] is in a better shape but also various tensions are identified that we discussed in detail in the latter paper. The message from this exercise is clear. The inclusion of charm in the evaluation of B 0 s,d −B 0 s,d hadronic matrix elements by LQCD is mandatory and it is important that in addition to HPQCD [19] a second LQCD collaboration includes charm in the evaluation of these matrix elements.
Assuming that the HPQCD values will be confirmed by another LQCD group, the SM predictions in the left panel of Fig. 7 will be favoured implying But it should be remembered that in contrast to the EXCLUSIVE scenario discussed in previous sections these results assume that the SM predictions for ∆M s and ε K agree with the data. If the EXCLUSIVE scenario will turn out to be true, the predictions above will be invalid and will be replaced by the ones in Table 2 and the green area in Fig. 7.

Conclusions
The EXCLUSIVE vision of rare decays and quark mixing is still not excluded and could become reality in the coming years. The present paper shows, similarly to our analysis in [7], how important is the determination of |V cb | for rare decays, in particular for rare Kaon decays. A precise determination of the γ in tree-level decays in the coming years will shed additional light on the tensions identified by us. As we have seen, an unusual pattern of SM predictions results from this study with some existing tensions disappearing or being dwarfed and new ones being born. In particular the B s → µ + µ − tension disappears and instead the anomalies at the level of (2 − 5)σ are present in ∆M s , ∆M d and in particular in ε K . While the 3.1σ tension in ∆M s is practically independent of γ, the one in ∆M d increases from 0.6 σ to 4σ when γ is decreased from 75 • to 60 • . In the case of ε K the corresponding variation is from 2.0σ to 5.2σ.
Moreover, the room left for NP in K + → π + νν, K L → π 0 νν and B → K(K * )νν is significantly increased but as seen in Figs.4 and 5 it depends sensitively on γ. The tension in B → X s γ is also interesting.
It should be recalled that in 2018 with the values of γ ≈ 74 • from the LHCb, with the inclusive |V cb | and the N f = 2 + 1 hadronic B 0 s,d −B 0 s,d matrix elements, ∆M d in the SM was found significantly above the data with a smaller enhancement in ∆M s [16]. In 2022 with lower values for γ ≈ 65 • from the LHCb [37] and N f = 2 + 1 + 1 hadronic B 0 s,d −B 0 s,d matrix elements from HPQCD [19] the inclusive values of |V cb | imply good agreement of the SM with the data on ∆M s,d . But for the exclusive values of |V cb | used in the present paper ∆M s is significantly below the experimental data. This also applies to ∆M d unless γ is chosen above 70 • that is not yet excluded by experiments.
In this context we should emphasize that the R(K) and R(K * ) anomalies being independent of |V cb | remain. On the other hand, as analysed in [12], lowering the value of |V cb | decreases the anomalies in B → Kµ + µ − , B → K * µ + µ − and B s → φµ + µ − decay branching ratios but to remove them completely values of |V cb | significantly lower than the exclusive ones are required.
It is premature to make a detailed analysis of possible BSM scenarios that could remove the anomalies in the EXCLUSIVE scenario considered by us. Despite of this let us close our paper with a few observations.
As in the EXCLUSIVE scenario NP is required to enhance ∆M s , ∆M d and ε K , a natural scenario would be at first sight the constrained Minimal Flavour Violation scenario [57] because, as pointed out in [58], in this scenario the ∆F = 2 observables can only be enhanced. However, the fact that a new phase ϕ new ≈ −1.3 • is required to fit the data for S ψK S , a more apprioprate here would be the U(2) 3 scenario [59][60][61]. As pointed out in [62], in the U(2) 3 scenario the CP-asymmetry S ψK S is anti-correlated with the CP-asymmetry S ψφ so that with |β s | ≈ 1 • an enhancement of the latter asymmetry from the SM prediction 0.0363±0.0013 to 0.080±0.020 would follow. Somewhat above the present data 0.054±0.020 [22] but consistent with it.
As a byproduct we have investigated in Section 4 the impact of the hadronic matrix elements from the HPQCD collaboration [19] on our results for rare B decays in [7]. The most interesting result is the increase of the B s → µ + µ − anomaly from 2.1σ to 2.7σ. Moreover we compared the determination of |V cb | from ∆M s , ∆M d , ε K and S ψK S using B 0 s,d −B 0 s,d hadronic matrix elements from LQCD with 2+1+1 flavours, 2+1 flavours and their average. As seen in Fig. 8 only for the 2 + 1 + 1 case values for β and γ can be found for which the same value of |V cb | is found. The resulting |V cb |, γ and β are given in (20) and |V ub | in (21). In any case the coming years will hopefully reveal for us which scenario for |V cb | and |V ub | has been chosen by nature. The measurement of γ combined with the 16 |V cb |-independent ratios constructed in [7] and with γ-dependence of various observables presented here will also play an important role in the search for NP. The importance of rare K decays in the search for NP has been recently summarized in [63] and the prospects for reducing hadronic uncertainties in K decays through intensive LQCD computations in the coming years are very good [64].
The EXCLUSIVE scenario appears to us to be more interesting than the HYBRID one because it implies more tensions between the SM predictions and the data. On the other the proponents of the inclusive determinations of |V cb | could consider the tensions found by us as an argument against exclusive determinations of |V cb |.
For N f = 2 + 1 + 1 one finds [19] In this case there are significant differences between the N f = 2 + 1 and N f = 2 + 1 + 1 results. Moreover, the latter ones are more accurate and we use them in the present paper. In [7] we have used the averages of both results given by [10] These are consistent with the ones from [19] alone but higher. However FLAG-2021 advices not to make such averages so that this time we use the N f = 2 + 1 + 1 values.