Emerging patterns of New Physics with and without Lepton Flavour Universal contributions

We perform a model-independent global fit to b→sℓ+ℓ-\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$b\rightarrow s\ell ^+\ell ^-$$\end{document} observables to confirm existing New Physics (NP) patterns (or scenarios) and to identify new ones emerging from the inclusion of the updated LHCb and Belle measurements of RK\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$R_K$$\end{document} and RK∗\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$R_{K^*}$$\end{document}, respectively. Our analysis, updating Refs. Capdevila et al. (J Virto JHEP 1801:093, 2018) and Algueró et al. (J Matias Phys Rev D 99(7):075017, 2019) and including these new data, suggests the presence of right-handed couplings encoded in the Wilson coefficients C9′μ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{{\mathcal {C}}}}_{9'\mu }$$\end{document} and C10′μ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{{\mathcal {C}}}}_{10'\mu }$$\end{document}. It also strengthens our earlier observation that a lepton flavour universality violating (LFUV) left-handed lepton coupling (C9μV=-C10μV\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{{\mathcal {C}}}}_{9\mu }^{\mathrm{V}}=-\,{{{\mathcal {C}}}}_{10\mu }^{\mathrm{V}}$$\end{document}), often preferred from the model building point of view, accommodates the data better if lepton-flavour universal (LFU) NP is allowed, in particular in C9U\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{{\mathcal {C}}}}_{9}^{\mathrm{U}}$$\end{document}. Furthermore, this scenario with LFU NP provides a simple and model-independent connection to the b→cτν\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$b\rightarrow c\tau \nu $$\end{document} anomalies, showing a preference of ≈7σ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\approx 7\,\sigma $$\end{document} with respect to the SM. It may also explain why fits to the whole set of b→sℓ+ℓ-\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$b\rightarrow s\ell ^+\ell ^-$$\end{document} data or to the subset of LFUV data exhibit stronger preferences for different NP scenarios. Finally, motivated by Z′\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$Z^\prime $$\end{document} models with vector-like quarks, we propose four new scenarios with LFU and LFUV NP contributions that give a very good fit to data.


Introduction
The flavour anomalies in b → s + − processes are at present among the most promising signals of New Physics a e-mail: bcapdevila@ifae.es (NP). Their analyses can be efficiently and consistently performed in a model-independent effective field theory (EFT) framework (see, for instance, [1][2][3]), where all short-distance physics (including NP) is encoded in Wilson coefficients, i.e. the coefficients of higher-dimension operators. A central open question is then which pattern(s) in the space of the Wilson coefficients is (are) preferred by b → s + − observables. More precise measurements, in particular for the observables showing deviations from the Standard Model (SM) expectations (P 5 [4], R K ,K * ,φ , Q 5 [5]…), help us to improve the results of this EFT analysis, which can then be used as a guideline for the construction of phenomenologically accurate NP models.
In this context we present here an update and extension of our recent works in Refs. [1,2], in the light of new measurements of key observables involved in b → s + − anomalies. We update the experimental value of the ratio probing lepton flavour universality (LFU) defined as R [1.1,6] K LHCb = 0.846 +0.060 +0.016 −0.054 −0.014 , R [1,6] K Belle = 0.98 +0.27 −0.23 ± 0.06 ,  [6], corresponding to the average of Run-1 and part of Run-2 (2015-2016) measurements, and the Belle collaboration [7], combining the data from charged and neutral modes. The correlations with the (finely binned) measurements of B(B → K μ + μ − ) [8] are tiny and therefore neglected here. In addition the Belle collaboration has also presented new results for R K * , the equivalent LFU-violating (LFUV) ratio for B → K * , in three bins [9], again considering both charged and neutral channels:  [15,19] K * = 1.18 +0.52 −0.32 ± 0.10.
Our treatment for the Belle observables within the global fit follows the same strategy as described in Ref. [1] for Q 4,5 where we introduced a nuisance parameter accounting for the relative weight of each isospin component. We have also updated our average for B(B s → μ + μ − ) including the latest measurement from the ATLAS collaboration [10] and taking into account the most recent lattice update of f B s for N f = 2 + 1 + 1 simulations collected in Ref. [11].
A relatively small numerical impact of such updates has been found. As in Ref. [1], our analysis also includes the lat-est update of P 4,5 from the Belle collaboration [12] where the muon and electron modes are considered separately (averaging charged and neutral modes), superseding the previous measurement in Ref. [13] where an average over both leptonic modes is presented. This allows us to include an additional measurement P 5μ (exhibiting a 2.6 σ discrepancy with respect to the SM) as well as the LFUV observable Q 5 in our analysis (see Ref. [14] for another recent analysis including this update).
In addition to updating the experimental inputs, our analysis explores new emerging directions in the parameter space spanned by the effective operators driven by data within two different frameworks. First, following Ref. [1] we assume in Sect. 2 that NP affects only muons and is thus purely Lepton-Flavour Universality Violating (LFUV). In Sect. 3 we follow the complementary approach discussed in Ref. [2], where we consider the consequences of removing the frequently made hypothesis that NP is purely LFUV. We then explore the implications of allowing both LFU and LFUV NP contributions to the Wilson coefficients C 9 ( ) and C 10 ( ) . The p-value of the SM hypothesis is 11.0% for the fit "All" and 8.0% for the fit LFUV Table 2 Most prominent 2D patterns of NP in b → sμ + μ − . The last five rows correspond to Hypothesis 1: (C NP 9μ = − C 9 μ , C NP 10μ = C 10 μ ), 2: Motivated by the new emerging directions in the LFUV case we also extend our analysis of NP scenarios to allow for the presence of LFU NP right handed-currents (RHC). In Sect. 4, we focus on a particular scenario (scenario 8) which can, within an EFT framework, link the flavour anomalies in b → s + − and b → c ν processes. Furthermore, we consider new patterns, motivated by Z models with vector-like quarks, which naturally predict LFU effects in C 10 ( ) complemented by LFUV ones. Finally, we sum up our results in Sect. 5. An appendix is devoted to the description of the correlations obtained for the various Wilson coefficients in the most relevant scenarios considered in this article.

Global fits in presence of LFUV NP
We start by considering the fits for NP scenarios which affect muon modes only. Tables 1, 2 and 3 and Fig. 1 update the corresponding tables and figures of Ref. [1] based on fits to the full set of data ("All") or restricted to quantities assessing LFUV. While we do not observe any significant difference in the 1D scenarios with "All" data compared to Ref. [1], some of the Pulls (with respect to the SM) for the LFUV 1D fits get reduced by half a standard deviation. A few other comments are in order: Table 3 1 and 2 σ confidence intervals for the NP contributions to Wilson coefficients in the 6D hypothesis allowing for NP in b → sμ + μ − operators dominant in the SM and their chirally-flipped counterparts, for the fit "All" The Pull SM is 5.1 σ and the p-value is 81.6%

Fig. 1
From left to right: allowed regions in the (C NP 9μ , C NP 10μ ), (C NP 9μ , C 9 μ ) and (C NP 9μ , C NP 9e ) planes for the corresponding 2D hypotheses, using all available data (fit "All") upper row or LFUV fit lower row Table 4 Coefficients for the polynomial parameterisation of the numerator and denominator of R [1.1,6] K in the vicinity of the SM point   [2,15], has an increased significance in the "All" fit compared to our earlier analysis. 2. The scenario C NP 9μ has the largest p-value in the "All" fit while C NP 9μ = − C NP 10μ has the largest p-value in the LFUV fit, a difference which can be solved through the introduction of LFU NP (see Ref. [2] and next section). 3. The best-fit point for the scenario C NP 9μ coincides now in the "All" and LFUV fits. 4. The scenario with only C NP 10μ has a significance in the "All" fit of only 4.0σ level and 3.9σ for the LFUV fit, which explains its absence from Table 1 as happens in Ref. [1].
Concerning the 2D scenarios collected in Table 2, the same picture arises as in Ref. [1], except that C NP 9e is now basically zero and small contributions to RHC seem slightly favoured (C 9 μ > 0, C 10 μ < 0). 1 with the coefficients provided in Table 4 (for linearised expressions, see Refs. [2,16]). We introduce a new Hyp. 5 in Table 2. The comparison between Hyps. 4 and 5 shows that the scenario C 9 μ = − C 10 μ (left-handed lepton coupling for right-handed quarks) prefers to be associated with C NP 9μ (vector lepton coupling for left-handed quarks) rather than C NP 9μ = − C NP 10μ (left-handed lepton coupling for left-handed quarks). Finally, no significant changes are observed in the 6D fit, except for the slight increase in the Pull SM , see Table 3. 1 Interestingly, these small contributions also reduce slightly the mild tension between P 5 at large and low recoils pointed out in Ref. [15] compared to the scenario with only C NP 9μ .
With the updated data, little change is observed among the preferred 2D NP models. Nevertheless, with an R [1.1,6] K value closer to one, scenarios with right-handed currents (RHC), namely (C NP 9μ , C 9 μ ) and (C NP 9μ , C 10 μ ), seem to emerge. The first scenario is naturally generated in a Z model with opposite couplings to right-handed and left-handed quarks and was proposed in Ref. [17] within the context of a gauged L μ − L τ symmetry with vector-like quarks. The latter (of masses m D and m Q ) are charged under L μ − L τ and have the same SM quantum numbers as right-handed down quarks and left-handed quark doublets, respectively. The vector-like quarks couple to the SM ones and to a scalar φ which breaks the L μ − L τ symmetry with couplings Y D,Q . We show the update of Fig. 2 of Ref. [17] assuming Y D,Q = 1 in Fig. 2. Since the current fit allows for C 9 μ = 0 at the two sigma level, the SU (2) singlet vector-like quark can still be decoupled [18].

Global fits in presence of LFUV and LFU NP
We turn to scenarios that allow also for the presence of LFU NP [2,15] (in addition to LFUV contributions to muons only), leading to the value of the Wilson coefficients (with i = 9, 10) for b → se + e − and b → sμ + μ − transitions respectively. We update some of the scenarios considered in Ref. [2] in Table 5. Concerning new directions in parameter space we allow for RHC, motivated by the results of the previous section, and focus on scenarios that could be fairly easily obtained in simple NP models.
With the updated experimental inputs, we confirm our earlier result [2] that a LFUV left-handed lepton coupling structure (corresponding to C V 9 = − C V 10 and preferred from a model-building point of view) yields a better description of data with the addition of LFU-NP in the coefficients C 9,10 , as shown by the scenarios 6, 8 in Table 5 with p-values larger than 70%.
We observe a very slight decrease in significance for the scenarios 5-7, with the exception of scenario 8 which exhibits one of the most significant pulls with respect to the SM. Scenario 8 of Ref. [2] can actually be realized via off-shell photon penguins [19] in a leptoquark model explaining also b → cτ ν data (we will return to this point in the following section).
Updated plots of the 2D LFU-LFUV scenarios discussed in Ref. [2] are shown in Fig. 3.
The new scenarios 9-13 are characterized by a C U 10( ) contribution. This arises naturally in models with modified Z couplings (to a good approximation C U 9( ) can be neglected). The pattern of scenario 9 occurs in Two-Higgs-Doublet models where this flavour universal effect can be supplemented by a C V 9 = − C V 10 effect [20]. In case of scenarios 11-13, one can invoke models with vector-like quarks where modified Z couplings are even Scenarios 5-8 were introduced in Ref. [2]. Scenarios 9 (motivated by 2HDMs [20]) and 10-13 (motivated by Z models with vector-like quarks [21]) are new induced at tree level. The LFU effect in C U 10( ) can be accompanied by a C V 9,10( ) effect from Z exchange [21]. Vectorlike quarks with the quantum numbers of right-handed down quarks (left-handed quarks doublets) generate effect in C U 10 and C V 9 (C U 10( ) and C V 9 ) for a Z boson with vector couplings to muons [21].
The comparison of scenarios 10 and 12 illustrates that C V 9μ plays an important role in LFU NP scenarios and can-not be swapped for its chirally-flipped counterpart without consequences. Finally, the allowed regions for the new LFU scenarios are displayed in Fig. 4.

Model-independent connection to b → c ν
In complement with the above EFT analysis, we focus now on the NP interpretation of scenario 8. Indeed, this sce- Fig. 4 Updated plots of Ref. [2] corresponding to the new scenarios 10, 11, 12 The red contour lines show the corresponding regions once R D ( * ) is included in the fit (for = 2 TeV). The horizontal blue (vertical yellow) band is consistent with R D ( * ) (R K ) at the 2 σ level and the contour lines show the predicted values for these ratios. Right: Impact of favoured NP scenarios on the observable P 5 . Only central values for the NP scenarios are displayed. The most interesting scenarios cluster together while traditional scenarios like C 9μ = − C 10μ or the scenario C 10μ considered in Ref. [28] fail to explain this anomaly nario allows for a model-independent connection between the anomalies in b → s + − and those in b → cτ ν, which are now at the 3.1σ level [22].
Such a correlation arises in the SMEFT scenario where C (1) = C (3) expressed in terms of gauge-invariant dimension-6 operators [23,24]. This scenario stems naturally from models with an SU (2) singlet vector leptoquark [25][26][27]. The operator involving-third generation leptons explains R D ( * ) and the one involving the second generation gives a LFUV effect in b → sμ + μ − processes. The constraint from b → cτ ν and SU (2) L invariance leads generally to large contributions to the operatorsγ μ P L bτ γ μ P L τ , which enhances b → sτ + τ − processes [24], but also mixes into O 9 and generates C U 9 at μ = m b [19]. Note that not all models addressing the charged and neutral current anomalies simultaneously have an anarchic flavour structure. In fact, in the case of alignment in the down-sector [29,30] one does not find large effects in b → sτ + τ − or C U 9 . Therefore, scenario 8 is reproduced in this setup with an additional correlation between C U 9 and R D ( * ) . Assuming a generic flavour structure so that small CKM elements can be neglected [19,24], we get Realizations of this scenario in specific NP models yield also an effect in C 7 generally [19]. However, since this effect is model dependent (and in fact small in some UV complete models [31,32]), we neglect it here, leading to the left plot in Fig. 5, where we include the recent update of Ref. [33] to draw the band for R D( * ) . Note that this scenario has a pull of 7.0 σ due to the inclusion of R D ( * ) , which increases χ 2 by ∼ 20.

Conclusions
In summary, including recent updates (R K , R K * and B(B s → μ + μ − )) our global model-independent analysis yields a very similar picture to the one previously found in Refs. [1,2] for the various NP scenarios of interest with some important peculiarities. In presence of LFUV NP contributions only, the 1D fits to "All" observables remain basically unchanged showing the preference for C NP 9μ scenario over C NP 9μ = − C NP 10μ . If only LFUV observables are considered the situation is reversed, as already found in Ref. [1], but now with an increased gap between the significances. This difference between the preferred hypotheses, depending on the data set used, can be solved introducing LFU NP contributions [2].
The main differences arise for the 2D scenarios: the cases including RHC, (C NP 9μ , C 10 μ ), (C NP 9μ , C 9 μ ) or (C NP 9μ , C 9 μ = − C 10 μ ), can accommodate better the recent updates, which enhances the significance of these scenarios compared to Ref. [1], pointing to new patterns including RHC. A more precise experimental measurement of the observable P 1 [34,35] would be very useful to confirm or not the presence of RHC NP encoded in C 9 μ and C 10 μ .
We also observe interesting changes in the 2D fits in the presence of LFU NP, where new scenarios (not considered in Ref. [2]) give a good fit to data with C U 10 ( ) and additional LFUV contributions. For example scenario 11 (C V 9μ , C 10 μ ) can accommodate b → s + − data very well, at the same level as scenario 8. Scenarios including LFU NP in left-handed currents (discussed in Ref. [2]) stay practically unchanged but with some preference for scenarios 6 and 8, which have a (V − A) structure for the LFUV-NP and a V or (V + A) structure for the LFU-NP. Furthermore, we have included additional scenarios 9 and 10 that exhibit a significance of 5.0σ and 5.5σ respectively.
We note that the amount of LFU NP is sensitive to the structure of the LFUV component. For instance, in scenario 7 (C V 9μ and C U 9 ) the LFU component is negligible at its best fit point. On the contrary, if the LFUV-NP has a (V − A) structure, the LFU-NP component (C U 9 ) is large, as illustrated by scenarios 6, 8 and 9. Scenarios with NP in RHC (either LFU or LFUV) prefer such contributions at the 2σ level (see scenarios 11 and 13) with the exception of scenario 12 with negligible C V 9 μ . The new values of R K and R K * seem thus to open a window for RHC contributions while the new B(B s → μμ) update (theory and experiment) helps only marginally scenarios with C NP 10μ . Finally, we showed that scenario 8, which allows for a model-independent connection between the b → cτ ν anomalies and the ones in b → s + − , can explain all data consistently and is preferred over the SM by 7 σ . Figure 5 illustrates the impact on the largest anomaly (P 5 ) of some of the most significant scenarios. Interestingly, several of the scenarios currently favoured cluster around the same values for the bins showing deviations with respect to the SM.
We have thus identified a number of NP scenarios with similarly good p-values and pulls with respect to the SM, which are able to reproduce the b → s + − data very well. Hierarchies among these scenarios can be identified, but additional data and reduced uncertainties are required to come to a final conclusion. The full exploitation of LHC run-2 data by the LHCb experiment (as well as by ATLAS and CMS) and the forthcoming results from the Belle and Belle II collaborations are expected to improve the situation very significantly in the forthcoming years, helping us to pin down the actual NP pattern hinted at by the b → s + − anomalies currently observed and to build accurate phenomenological models to be confirmed through other experimental probes such as direct production experiments.
Note added After the completion of this work, several global analyses have been performed to assess NP scenarios affecting b → s + − processes [14,28,36,37]. They agree well with our findings, with small differences stemming mainly from slightly different theoretical approaches as well as theoretical and experimental inputs. The improve-ment brought by RHC has been observed in Refs. [14,36], whereas the interest of LFU NP contributions is also identified in Refs. [14,28,38]. Most of the analyses observe that the slight deviation from B(B s → μ + μ − ) plays no specific role in the global fit [36,37], apart from Ref. [28]. In the latter analysis, the significance of a scenario with only C NP 10μ is much more important than in our case, and the hierarchies between the significances of 2D scenarios is different. After discussion with the authors of Ref. [28], this difference comes from their inclusion of B s -B s mixing and the assumption that F = 2 observables are purely governed by the SM, which helps them sharpening the prediction for B(B s → μ + μ − ) and increase the weight of this observable in the fit. Our present analysis does not rely on this strong hypothesis, which should be contrasted with the fact that most models invoked to explain b → s + − anomalies typically affect also F = 2 observables.
1 Correlation matrices of fits to LFUV NP First, we present the correlations between fit parameters of the NP scenarios defined in Tables 2 and 3. These are all NP solutions whose parameters assess LFUV NP.
By order of appearance in Table 2 The last two matrices correspond to Hyp. 1 and Hyp. 5 as defined in Table 2. Despite the high Pull SM of the 2D scenario {C NP 9μ , C 7 } (5.4σ ), its correlation matrix is not collected here due to the value of C 7 being negligible, with tiny errors.
Regarding the 6D fit of Table 3, where the columns are ordered as {C NP 7 , C NP 9μ , C NP 10μ , C 7 , C 9 μ , C 10 μ }. Interesting information can be extracted from Corr 6D . Most of the coefficients do not show particularly strong correlations with the others except for the pairs {C NP 10μ , C 9 μ }, {C NP 10μ , C 10 μ } and {C 9 μ , C 10 μ }, being the latter the highest in correlation. While C NP 9μ and C 9 μ show a non-negligible correlation in the fit to these coefficients only, in the 6D fit the aforementioned parameters are uncorrelated to a large extent. On the contrary, the correlation between C NP 9μ and C NP 10μ is very similar for both the global 6D and the 2D fit to these parameters alone.
2 Correlation matrices of fits to LFUV-LFU NP Second, the correlations between fit parameters of scenarios with both LFUV and LFU NP have also been considered.
Below one can find the correlation matrices of scenarios 5-11, in that order. No significant changes can be observed when comparing with the results in App. 2 of Ref. [2]. As expected, C V 9μ and C U 9 are highly anti-correlated, with its nominal value somewhat smaller than in [2]. Fit estimates of the parameters in scenario {C V 9μ = −C V 10μ , C U 9 = C U 10 } are now slightly correlated, while before their correlation was negligible. Interestingly, however, we find the parameters of the new scenario {C V 9μ , C U 10 } statistically independent up to a large extent.