Signs for new physics in the recent LHCb data?

We comment on some tensions with the Standard Model predictions in the recent LHCb data.


Introduction
With the first measurement of new angular observables in the exclusive decay B → K * µ + µ − based on the 1 fb −1 dataset, LHCb has found a 4.0σ local discrepancy in one of the q 2 bins for one of the angular observables [1], namely in the bin q 2 ∈ [4. 3,8.63] GeV 2 of the observable P 5 . The latter belongs to the set of optimised observables in which form factor dependence cancels out to first order. LHCb results are compared here with the theoretical predictions in Ref. [2]. Intriguingly, other smaller but consistent deviations are also present in other observables [1].
In the low-q 2 region, the up-to-date description of exclusive heavy-to-light B → K * µ + µ − decays is the method of QCD-improved Factorisation (QCDF) and its field-theoretical formulation of Soft-Collinear Effective Theory (SCET). In the combined limit of a heavy b-quark and of an energetic K * meson, the decay amplitude factorises to leading order in Λ/m b and to all orders in α s into process-independent non-perturbative quantities like B → K * form factors and light-cone distribution amplitudes (LCDAs) of the heavy (light) mesons and perturbatively calculable quantities, which are known to O(α 1 s ) [3,4]. Further, the seven a priori independent B → K * QCD form factors reduce to two universal soft form factors ξ ⊥, [5]. The factorisation formula applies well in the dilepton mass range 1 GeV 2 < q 2 < 6 GeV 2 .
Taking into account all these simplifications the various K * spin amplitudes at leading order in Λ QCD /m b and α S turn out to be linear in the soft form factors ξ ⊥, and also in the short-distance Wilson coefficients. As was explicitly shown in Refs. [6,7], these simplifications allow to design a set of optimised observables, in which any soft form factor dependence (and its corresponding uncertainty) cancels out for all low dilepton mass squared q 2 at leading order in α S and Λ QCD /m b . An optimised set of independent observables was constructed in Refs. [2,8], in which almost all observables are free from hadronic uncertainties which are related to the form factors.
However, the soft form factors are not the only source of hadronic uncertainties in these angular observables. It is well-known that within the QCDF/SCET approach, a general, quantitative method to estimate the important Λ QCD /m b corrections to the heavy quark limit is missing. It is clear that the interpretation of the LHC measurement strongly depends on the treatment of this problem as we discuss in the next section.
There is another issue: The validity of the theory predictions based on QCD factorisation approach within the region q 2 ∈ [4.3, 8.63] GeV 2 is highly questionable. The validity is commonly assumed up to 6 GeV 2 for two reasons. The perturbative description of the charm loops is valid in this region and also the kinematical assumptions about the large energy of the K * within the SCET/QCD factorisation approach are still reasonable. Thus, using the theory predictions up to 8.63 GeV 2 could induce larger hadronic corrections.
Leaving these two issues aside, it has been shown that the deviation in the observable P 5 and the small deviations in other observables in the low-q 2 area, can be consistently described by a smaller C 9 Wilson coefficient, together with a less significant contribution of a non-zero C 9 (see for example Ref. [9]). This is a challenge for the model-building as we will discuss below.
Thus, it is not clear if the anomaly is a sign for new physics beyond the SM, or a consequence of underestimated hadronic power corrections or non-perturbative charm effects or just a statistical fluctuation. The LHCb analysis based on the 3 fb −1 dataset is eagerly awaited to clarify the situation.
More recently, another small discrepancy occurred. The ratio R K = BR(B + → K + µ + µ − )/BR(B + → K + e + e − ) in the low-q 2 region has been measured by LHCb showing a 2.6σ deviation from the SM prediction [10]. In contrast to the anomaly in the rare decay B → K * µ + µ − which is affected by unknown power corrections, the ratio R K is theoretically rather clean. This might be a sign for lepton non-universality.

Power corrections
In spite of the fact that the power corrections cannot be calculated, the corresponding uncertainties should be made manifest within the theory predictions. Therefore, in Refs. [6,7] the effects of the Λ QCD /m b corrections has been parametrised for each of the K * 0 spin-amplitudes with some unknown linear correction. In case of CPconserving observables this just means where C i is the relative amplitude. In the case of CPviolating observables, a strong phase has to be included (see Ref. [7] for details). It is further assumed that these amplitudes (C i ) are not functions of q 2 , although in practice they may actually be, and any unknown correlations are also ignored.
Based on a simple dimensional estimate, one has chosen |C i | < 10%. There are also soft arguments for this choice: Under the assumption that the main part of the Λ QCD /m b corrections is included in the full form factors, the difference of the theoretical results using the full QCD form factors on one hand and the soft form factors on the other hand confirms this simple dimensional estimate. In fact, the comparison of the approaches leads to a 7% shift of the central value at the level of observables. Secondly, one can state that the chiral enhancement of Λ QCD /m b corrections in the case of hadronic B decays does not occur in the case of the semileptonic decay mode with a vector final state. Thus, it is not expected that they are as large as 20 − 30% as in the B → ππ decay.
This procedure was introduced to make the unknown Λ QCD /m b corrections manifest; this ansatz, put in by hand, was never meant as a real quantitative estimate. In fact assuming 10% corrections on the amplitude level one often finds power corrections of less than 5% on the observable level as one can read off from Table 5 and 6 of Ref. [2], in which the hand-made error due to Λ QCD /m b corrections is nicely separated from the error due to input parameters and scale dependence. So the errors due to the power corrections given in Ref. [2] are most probably an underestimation of the hadronic uncertainty if taken as a real quantitative estimate. However, if we assume 10% error due to the unknown power corrections -which corresponds to a naive dimension estimate of Λ/m b on the level of observables and is also backed up by some soft arguments (see above) -we find the pull in case of the third bin of the observable P 5 reduced from 4.0σ to 3.6σ what still represents a significant deviation. And even if one assumes 30% error then the pull in this case is still 2.2σ within the modelindependent analysis presented in Ref. [11].
In Ref. [12] a general parametrisation for the power corrections to the form factor terms (the factorisable piece in QCD factorisation) is given, There are two free parameters in the ansatz for each QCD form factor which have to be determined. In a more recent analysis [13], this procedure gets further developed. The authors analyse the factorisable power corrections by looking at the difference of QCD form factors and soft form factors, an argument that was already used to get an order of magnitude estimate of power corrections. Unfortunately, the errors of the light-cone sum rule calculations of QCD form factors are too large and the correlations between the various sum rule calculations are not commonly known in order to get reasonable results. Thus, in Ref. [13] only the central values get fixed by this procedure and a 10% error is attributed to the power corrections by hand again. Another argument is put forward that one can choose the renormalisation scheme of the soft form factors such that the power corrections to certain observables are reduced. In addition one has to estimate the non-factorisable power corrections. Here an observation is helpful [14], that non-factorisable corrections are not induced by the leading electromagnetic and semileptonic operators. However, also these contributions have to be estimated by a hand-made ansatz; this is done in Ref. [13] by a q 2 dependent correction of order 10% on the amplitude level. Clearly, an improved, but honest light-cone sum rule calculation of the various QCD form factors including their correlations is really needed to make progress on this issue. But the non-factorisable power contributions will then be still a source of uncertainty limiting the new physics sensitivity of these exclusive decay modes.
In principle, the light-cone sum rule (LCSR) approach allows to estimate these non-factorisable contributions (with the well-known uncertainties of QCD sum rules) as was demonstrated in the case of the decay B → K + − [15]. For the B → K * + − case, only softgluon contributions of the charm loop effects has been considered yet [16].

New physics interpretations of the anomaly
Consistent SM and new physics interpretations of the measured deviation in the B → K * + − mode have been discussed in a large number of references [11,[17][18][19][20][21][22][23][24][25][26][27]. In a model-independent analysis, the anomaly can be consistently described by smaller C 9 and C 9 Wilson coefficients. The usual suspects like the MSSM, warped extra dimension scenarios, or models with partial compositeness, cannot accommodate the deviation at the 1σ level, but Z models may do this [18].
In the MSSM, we have no means of generating any sizeable contribution to the coefficient C 9 , but also any significant contribution to C 9 is correlated to contributions to other Wilson coefficients affecting the other observables. Nevertheless, combining all the observables in a fit one can check the global agreement of the model with the available data [26]. This is shown in Fig. 1 for the relevant Wilson coefficients in the constrained MSSM (CMSSM) and in the more general setup of the phenomenological MSSM (pMSSM), where all the generated points are shown (grey), indicating the values for the Wilson coefficients reachable within the MSSM, as well as the points satisfying the global 1,2,3 σ constraints (red, yellow, green). The Higgs mass constraint has also been imposed. As can be seen, the overall agreement is fairly good, with regions in SUSY parameter space where the absolute χ 2 is sufficiently small and an agreement at the 1σ level is reached.

Signs for lepton non-universality
Besides this known anomaly in the angular analysis of B → K * µ + µ − decay, another small discrepancy recently occurred. The ratio R K = BR(B + → K + µ + µ − )/BR(B + → K + e + e − ) in the low-q 2 region has been measured by LHCb showing a 2.6σ deviation from the SM prediction [10]. This discrepancy has been addressed in a few recent studies [28][29][30][31][32][33][34]. A global fit to all observables considering separately new physics contributions to the electron and muon semileptonic Wilson coefficients C e 9,10 and C µ 9,10 (and the corresponding chirality flipped coefficients) is shown in Figs. 2 and 3 (see [33] for more detail). Fig. 2 shows the fit results for C µ 9 , C e 9 , C µ 10 , C e 10 , while Fig. 3 presents the results for C µ 9 , C µ 9 , C e 9 , C e 9 . We see that the SM is disfavoured at the 2σ level. Yet there is tension in the muon sector for C 9 . In order to compare the results with the scan for two operators O µ 9 , O e 9 only, the contours corresponding to the 1 and 2σ fit result for C µ 9 , C e 9 are overlaid in the figures. This shows that considering arbitrarily only two operators can be too restrictive and even misleading since a large area of new physics parameter space might be unjustifiably overlooked. For example, in the two-operator fit lepton-universality, δC µ 9 = δC e 9 , is disfavoured by 2σ, while within the four-operator fit the agreement is improved.

Cross-checks with the inclusive mode
The inclusive mode B → X s + − can only be measured at e + e − machines and is theoretically cleaner than the exclusive modes [35,36]. The theoretical accuracy in the low-q 2 region is of the order of 10% [37]. But the branching fraction has been measured by Belle and BaBar using the sum-of-exclusive technique only.
In Refs. [25,33] we checked the compatibility of the present datasets of the exclusive and inclusive modes. It is remarkable result that the sets of exclusion plots are nicely compatible with each other. This is a nontrivial consistency check. At the moment, the measurements of the B → K * + − observables are the most powerful ones. However, the latest published measurement of Belle [38] is based on a sample of 152 × 10 6 BB events only, which corresponds to less than 30% of the dataset available at the end of the Belle experiment while BaBar has just recently presented an analysis based on the whole dataset using a sample of 471 × 10 6 BB events [39] overwriting the previous measurement from 2004 based on 89 × 10 6 BB events [40]. Moreover, there will be a Super-B factory Belle-II with a final integrated luminosity of 50 ab −1 [41]. There is a recent analysis [42] of the expected total uncertainty on the partial decay width and the forward-backward asymmetry in several bins of dilepton mass-squared for the fully inclusive B → X s + − decays assuming a 50 ab −1 total integrated luminosity (for details see Ref. [25]). One finds a relative fractional uncertainty of 2.9% (4.1%) for the branching fraction in the low-(high-)q 2 region and a total absolute uncertainty of 0.050 in the low-q 2 bin 1 (1 < q 2 < 3.5 GeV 2 ), 0.054 in the low-q 2 bin 2 (3.5 < q 2 < 6 GeV 2 ) and 0.058 in the high-q 2 interval (q 2 > 14.4 GeV 2 ) for the normalised A FB . So the inclusive mode will lead to very strong constraints on the Wilson coefficients We illustrate the usefulness of these future measurements of the inclusive mode at Belle-II in the following way [33]. We make a model independent fit for the coefficients C 7 , C 8 , C 9 , C 10 and C l (for notation see Ref. [25]). In addition to all the b → s + − observables, we consider the inclusive branching ratio of B → X s γ as well as the isospin asymmetry in B → K * γ decay which are relevant to constrain C 7 and C 8 . Based on our model-independent analysis we predict the branching ratio at low-and high-q 2 . In Fig. 4, we show the 1, 2, and 3σ ranges for these observables. In addition, we add the future measurements at Belle-II assuming the best fit solution of our model-independent analysis as central value. These measurements are indicated by the black error bars. They should be compared with the theoretical SM predictions given by the red error bars. Fig. 4 indicates that the future measurement of the inclusive branching ratios separates nicely from the SM prediction as the model-independent fit. And also the future measurement of the forward-backward asymmetry at Belle-II will allow us to separate the potential new physics measurement from the SM prediction in a significant way as shown in Fig. 5.

Global fit in MFV
Assuming the tensions seen in the recent LHCb results are hints for new physics, it is important to con- Figure 5: 1, 2 and 3σ ranges for the unnormalised forward-backward asymmetry in bin 1 (1 < q 2 < 3.5 GeV 2 ) and in bin 2 (3.5 < q 2 < 6 GeV 2 ) within the model-independent analysis. Future measurement at the high-luminosity Belle-II Super-B-Factory assuming the best-fit point of the model-independent analysis as central value (black) and the SM predictions (red/grey). sider whether new flavour structures are needed to explain the data. The hypothesis of MFV [43][44][45][46][47] implies that flavour and CP symmetries are broken as in the SM. Thus, it requires that all flavour-and CP-violating interactions be linked to the known structure of Yukawa couplings. The MFV hypothesis represents an important benchmark in the sense that any measurement which is inconsistent with the general constraints and relations induced by the MFV hypothesis unambiguously indicates the existence of new flavour structures.
We study the results of the global fit for the new physics contributions to the Wilson coefficients C 7 , C 8 , C 9 , C 10 and C l , and update the results of Refs. [25,47] based on the latest experimental results. As can be seen in Fig. 6, while the 1 and 2σ allowed regions are squeezed compared to those of [25,47] which shows the impact of the new measurements, the overall agreement of the MFV solutions with the data is still very good, and no new flavour structure is needed to explain the experimental results. Figure 6: Global fit to the NP contributions δC i in the MFV effective theory, at 1 (red) and 2σ (green).