Constraints on tensor and scalar couplings from [Formula: see text] and [Formula: see text].

The angular distribution of [Formula: see text] ([Formula: see text]) depends on two parameters, the lepton forward-backward asymmetry, [Formula: see text], and the flat term, [Formula: see text]. Both are strongly suppressed in the standard model and constitute sensitive probes of tensor and scalar contributions. We use the latest experimental results for [Formula: see text] in combination with the branching ratio of [Formula: see text] to derive the strongest model-independent bounds on tensor and scalar effective couplings to date. The measurement of [Formula: see text] provides a complementary constraint to that of the branching ratio of [Formula: see text] and allows us - for the first time - to constrain all complex-valued (pseudo-)scalar couplings and their chirality-flipped counterparts in one fit. Based on Bayesian fits of various scenarios, we find that our bounds even become tighter when vector couplings are allowed to deviate from the standard model and that specific combinations of angular observables in [Formula: see text] are still allowed to be up to two orders of magnitude larger than in the standard model, which would place them in the region of LHCb's sensitivity.

Introduction

With the analysis of the data collected by the LHCb Collaboration during run I at the Large Hadron Collider (LHC), we now have access to rather large samples of rare B-meson decays with branching ratios below . As a consequence, angular analyses of three- and four-body final states can be used to measure a larger number of observables than previously possible at the B factories BaBar and Belle. In this work we focus on rare B decays driven at the parton level by the flavor-changing neutral-current (FCNC) transition that constitutes a valuable probe of the standard model (SM) and provides constraints on its extensions.

The angular distribution of – normalized to the width – in the angle between B and as measured in the dilepton rest frame isLHCb analyzed their full run 1 data set of 3 fb and measured the angular distribution of the mode , i.e.  [1], with unprecedented precision. They provide the lepton forward–backward asymmetry and the flat term in CP-averaged form and integrated over several bins in the dilepton invariant mass . Similarly, the CP-averaged branching ratios, , [2] and the rate CP asymmetry [3] are also available from 3 fb.

Both angular observables, and , exhibit strong suppression factors for vector and dipole couplings present in the SM, thereby enhancing their sensitivity to tensor and scalar couplings [4, 5]. A similar enhancement of scalar couplings compared to helicity-suppressed vector couplings of the SM is well known from . Unfortunately the limited data set of from LHCb [6] did not yet allow one to perform a full angular analysis without the assumption of vanishing scalar and tensor couplings in this decay mode. In the future with more data or special-purpose analysis techniques like the method of moments [7], certain angular observables in will provide additional constraints on such couplings, as for example [8] and the linear combinations and [5, 9] as well as the experimental test of the relations and  [5] at low hadronic recoil.

Here we exploit current data from and to derive stronger constraints than before on tensor and scalar couplings in various model-independent scenarios and study their impact on the not-yet-measured sensitive observables in . In Sect. 2, we specify the effective theory of decays on which our model-independent fits are based. Within this theory, we discuss the dependence of observables in and on the tensor and scalar couplings in Sect. 3 and specify also the experimental input used in the fits. The constraints on tensor and scalar couplings from the data are presented for several model-independent scenarios in Sect. 4. Technical details of the angular observables in , the branching fraction of , the treatment of theory uncertainties, and the Monte Carlo methods used are relegated to appendices.

Effective theory

In the framework of the effective theorythe most general dimension-six flavor-changing operators mediating and are classified according to their chiral structure. There are dipole () and vector () operators,further scalar () operators,and tensor () operators,where the notation is also used frequently in the literature. The respective short-distance couplings, the Wilson coefficients , are evaluated at a scale of the order of the b-quark mass and can be modified from SM predictions in the presence of new physics.

The SM values are obtained at next-to-next-to leading order (NNLO) [10, 11] and depend on the fundamental parameters of the top-quark and W-boson masses, as well as on the sine of the weak mixing angle. Moreover, they are universal for the three lepton flavors . All other Wilson coefficients are numerically suppressed or zero: , , and . The Wilson coefficients of the four-quark current–current and QCD-penguin operators as well as of the chromomagnetic dipole operators are set to their NNLO SM values at [10, 11].

For the rest of this article, we will suppress the lepton-flavor index on the Wilson coefficients and operators . In Sect. 4 we exploit data with only, hence all derived constraints apply in principle only to the muonic case but can be carried over to the other lepton flavors for NP models that do not violate lepton flavor. In general, the Wilson coefficients are decomposed into SM and NP contributions but often we will use for Wilson coefficients with zero (or suppressed) SM contributions synonymously with .

Observables and experimental input

The full dependence of and on tensor and scalar couplings has been presented in [4, 5], adopting the effective theory (2.1), i.e. neglecting higher-dimensional operators with . These results imply that for SM values of the effective couplingsHence, for both observables are quasi-null tests. The flat term is strongly suppressed by small lepton masses for the considered kinematic region GeV [4, 12, 13]. Nonzero values of can be induced by higher-order QED corrections, which will modify the simple dependence of the angular distribution (1.1); however, currently there is no solid estimate available for this source of SM background.1 This picture does not change in the presence of new-physics contributions to vector and dipole operators . On the other hand, nonvanishing tensor or scalar contributions are enhanced unless the dynamics of NP implies similar suppression factors, i.e., lepton–Yukawa couplings for or in the case of . In particular, is very sensitive to tensor couplings (see (3.3) below) and is sensitive to the interference of tensor and scalar couplings (see (3.6) below).

There are some angular observables in with the same properties; i.e. tensor and scalar contributions are kinematically enhanced by a factor over vector ones present in the SM or their respective interference terms. These are and the two linear combinations and with explicit formulas given in Appendix A. In our fits and predictions we include all kinematically suppressed terms. But for the purpose of illustration, we now consider the analytical dependence for vanishing lepton mass. In this limit,is sensitive to the interference of tensor and scalar operators, complementary to in i.e., with an interchange of tensor contributions . We note also that contributes to the lepton forward–backward asymmetry of being . Since it has to compete with in this observable, a separate measurement of and is necessary.

Only tensor contributions enterwhere the dots indicate different kinematic and form-factor dependencies. But tensor and scalar contributions enterwhich is similar to the dependence of in Concerning , the involved kinematic factors – see [4, 5] – are such that tensor and scalar couplings contribute only constructively/cumulatively, apart from cancellations among and . Interference terms in the numerator of of the form and are suppressed by . They become numerically relevant in case or where the smallness of is of the same level as the suppression factor accompanying the large vectorial SM Wilson coefficients . This implies, however, no large enhancement of over the SM prediction.

On the one hand, the observables (3.6) and (3.3) are measured in the angular distribution (1.1) of normalized to the decay width such that uncertainties due to form factors can cancel in part [4, 5]. On the other hand, , , and appear in the unnormalized angular distribution of . “Optimized” versions , , and for the low- region for which form factors cancel in the limit of have been identified in [9]. For the high- region, potential normalizations are discussed in Appendix A, which could serve to form optimized observables for special scenarios of either vanishing chirality-flipped vector or tensor or scalar couplings. In the most general case, however, there are no optimized observables at high . Although form factors do not cancel in this case, it might still be preferable to use normalizations, for example when the overall normalization of form factors constitutes a major theoretical uncertainty.

To illustrate the sensitivity of to tensor couplings, we compare it in Fig. 1 to the branching ratios of for , integrated over one low- and one high- bin. The details of the numerical input and the uncertainty propagation can be found in Appendices B and C. In the light of the hint of new physics in from recent global analyses of data [26-30], we show predictions for in addition to .

Fig. 1

The sensitivity to the tensor coupling of and in as well as , , and in . Angular observables are rescaled by the lifetime of the B meson, . The bands represent the theory uncertainties at 68 and 95 % probability of the prior predictive. Two sets of bands are shown for (blue) and (red). If available, the gray band indicates the latest 68 % confidence interval reported by LHCb. All observables are integrated over bins denoted as to match LHCb

From Fig. 1, the highest sensitivity to tensor couplings of any observable is attained by at high due to a partial cancellation of form factors [5]. If the experimental uncertainty could be reduced further, would give a very strong constraint on a simultaneous negative shift in and . The prediction of is essentially insensitive to but sensitive to . A stronger impact on global fits, however, would require a reduced theory uncertainty.

The observable shows moderate dependence on at least at low and has some impact on the constraints on tensor couplings as will be discussed in Sect. 4. At the moment, theory and experimental uncertainty are of similar size.

is sensitive to in both regimes. At low , it is mildly affected by , whereas at high it is unaffected. Regarding , the situation is reversed: here the strong dependence on appears at low . Overall, , at high , and are sensitive to and theoretically very clean around .

From the available measurements, at high currently provides the most stringent constraints on the size of tensor couplings. Moreover, the dependence on vector couplings is such that leads to stronger constraints on  than .

Important additional constraints on scalar couplings come from the branching ratio of as given in (A.8). It provides the most stringent constraints on the moduli and and further depends only on . Thus it is complementary to in ; see Eq. (3.6).

Eventually we also explore the effect of interference with NP contributions in the vector couplings on the bounds on tensor and scalar couplings. For this purpose we include also the branching ratio, the lepton forward–backward asymmetry, and the rate CP asymmetry of as they provide additional constraints on the real and imaginary parts of . The experimental input of all observables entering our fits is listed in Table 1 together with input for the form factors. More details on the latter can be found in Appendix B.

Table 1

List of all observables of the various decays entering the fits with the respective kinematics and experiments that provide the measurements. LCSR and lattice results of form factors are used to constrain a -dependent form-factor parametrization. For more details see Sect. 3 and Appendix B

Channel	Constraints	Kinematics	Source
B_s→μ¯μ	B¯≡∫dτB(τ)	–	[17–19]
B^+→K^+μ¯μ	B_μ	q^2∈[1,6],[14.18,16],[>16] GeV^2	[20]
		q^2∈[1.1,6.0],[15.0,22.0] GeV^2	[2]
	AFBμ	q^2∈[1.1,6.0],[15.0,22.0] GeV^2	[1, 20]
	FHμ	q^2∈[1.1,6.0],[15.0,22.0] GeV^2	[1]
	ACPμ	q^2∈[1.1,6.0],[15.0,22.0] GeV^2	[3]
B^0→K^∗0μ¯μ	B_μ	q^2∈[1,6],[14.18,16],[>16] GeV^2	[20–22]
	AFBμ	q^2∈[1,6],[14.18,16],[>16] GeV^2	[20–22]
	ACPμ	q^2∈[1.1,6.0],[15.0,22.0] GeV^2	[3]
B→K form factors	f_0,+,T	q^2=17,20,23 GeV^2	[13]
B→K^∗ form factors	V,A_0,1,2,T_1,2,3	q^2=0.1,4.1,8.1,12.1 GeV^2	[23]
	V,A_0,1,2,T_1,2,3	q^2∈[11.9,17.8] GeV^2	[24, 25]

Fits and constraints

There are no discrepancies between the latest measurements for (throughout this section) of and in and their tiny SM predictions; cf. Fig. 1. Thus our main objective is to derive constraints on tensor and scalar couplings through the enhanced sensitivity of both observables to these couplings compared to vector couplings. For this purpose, we will consider several model-independent scenarios, progressing from rather restricted to more general ones in order to asses the effect of cancellations due to interference of various contributions.

For each coupling that we vary in a fit, we remain as general as possible, treat it as a complex number and use the Cartesian parametrization assuming uniform priors for ease of comparison with previous studies. Specifically, we setand the same for the imaginary parts. The priors of the nuisance parameters are given in Appendix B.

We start with the scenario of only tensor couplings and see that they are well constrained by alone. In a second scenario we consider only scalar couplings in order to investigate the complementarity of and . Here we find that – for the first time – all complex-valued scalar couplings can be bounded simultaneously by the combination of both measurements. Finally we consider as a special scenario the SM augmented by dimension-six operators as an effective theory of new physics below some high scale assumed much larger than the typical scale of electroweak symmetry breaking. In addition, the model contains one scalar doublet under as in the SM. For each scenario, we also investigate interference effects with new physics in vector couplings . Finally, we conclude this section with posterior predictions – conditional on all experimental constraints – of the probable ranges of the not-yet-measured angular observables , and ) in and in .

Tensor couplings

In a scenario with only complex-valued tensor couplings , the experimental measurement of constrains the combination , up to some small interference of with vector couplings ; cf. (3.6). The according 68 % (95 %) 1D-marginalized probability intervals are listed in the second column of Table 2. From the third column, it is seen that the constraints become tighter when utilizing all observables in Table 1, mainly due to the sensitivity of the branching ratio of to tensor couplings (see also Fig. 1). The latter stronger bounds are driven by the new lattice results of form factors that predict values above the measured ones [31]. Since tensor couplings contribute constructively to , large values are better constrained. In this scenario with vanishing scalar couplings, current measurements of barely provide any constraint; cf. (3.3).

Table 2

The constraints on complex-valued when using measurements of only , and other data in Table 1 (except ), and finally allowing for complex-valued new-physics contributions to

Data set	Only FHμ	FHμ + other	FHμ + other
Set of couplings	CT,T5	CT,T5	CT,T5,9,10
Credibility level	68 %, 95 %	68 %, 95 %	68 %, 95 %
ReCT	[-0.32,0.16], [-0.52,0.35]	[-0.29,0.09], [-0.43,0.27]	[-0.23,0.14], [-0.39,0.30]
ImCT	[-0.25,0.24], [-0.44,0.44]	[-0.19,0.21], [-0.37,0.35]	[-0.17,0.22], [-0.33,0.37]
ReCT5	[-0.25,0.24], [-0.44,0.44]	[-0.19,0.17], [-0.33,0.33]	[-0.18,0.16], [-0.32,0.32]
ImCT5	[-0.24,0.25], [-0.44,0.45]	[-0.21,0.19], [-0.37,0.35]	[-0.20,0.18], [-0.36,0.35]

We also perform a fit with nonzero in order to assess the robustness of the bounds with respect to interference. Note that appears in linear combinations with such that its interference with tensor and scalar couplings is captured implicitly by allowing for new physics in . Thus we fix to the SM value without loss of generality. In this case, by itself still provides bounds on that are weakened by a factor of 2, since does not pose constraints on (in the chosen prior range). Once additional experimental measurements of Table 1 are taken into account, the potential destructive effects of new physics in become reduced and almost the same constraints on are recovered, as shown in the last column in Table 2. If in addition we allow for (not shown in Table 2), the credible regions further shrink by about 10%, which we attribute to the cumulative effect of in ; cf. (3.6). In summary, the measurement [1] of LHCb with 3 fb shrinks the previous bounds [5] on by roughly 50 %.

Scalar couplings

Scalar couplings enter without kinematic suppression – see (3.6) – as the sum whereas in the time-integrated branching ratio they appear as the difference . Since the existing measurement of constrains the sum, the combination of and allows us – for the first time – to bound the real and imaginary parts of all four couplings. The corresponding 2D-marginalized regions in the () planes are shown in Fig. 2. The corresponding plots for are very similar to those shown and thus omitted. These bounds do not change when including all other data in Table 1, since requires interference of scalar with tensor couplings and the other observables are not very sensitive to scalar couplings. Quantitatively, the constraint from on is about a factor 4 to 5 stronger than the one of on .

Fig. 2

The constraints on complex-valued from only  (gray dotted), only (blue dashed), and the combination with all other data in Table 1 as well as nonzero (red solid) at 68 % (darker) and 95% (lighter) probability. The constraints on Re are identical to Re, apart from a small translation of the contours by . The SM prediction is indicated by the black diamond

Interference terms of with vector couplings might weaken these bounds. For , the relevant term is (see (A.7)) and for it is [4]; both are suppressed by the factors and , respectively. Nevertheless, these terms become important for small due to the large SM value of . We compile bounds on complex-valued scalar couplings in Table 3 for only , only , and their combination with all the other observables in Table 1. Neither nor alone can bound all four complex-valued scalar couplings, however, their combination is capable to do so and moreover, the bounds are stable against destructive interference with vector couplings.

Table 3

The 1D-marginalized constraints on complex-valued at 68% (95%) probability from measurements of only , only , and all the data in Table 1 and additional new-physics contributions to

Data set	Only FHμ	Only B¯(B_s→μ¯μ)	All
Set of couplings	CS,S^′,P,P^′	CS,S^′,P,P^′	CS,S^′,P,P^′,9,9^′,10,10^′
Credibility level	68 %, 95 %	68 %, 95 %	68 %, 95 %
Re(CS-CS^′)	–	[-0.10,0.08], [-0.14,0.13]	[-0.08,0.07], [-0.13,0.13]
Im(CS-CS^′)	–	[-0.07,0.07], [-0.11,0.12]	[-0.07,0.07], [-0.12,0.11]
Re(CS+CS^′)	[-0.36,0.39], [-0.69,0.68]	-	[-0.32,0.32], [-0.59,0.62]
Im(CS+CS^′)	[-0.37,0.35], [-0.68,0.66]	-	[-0.25,0.41], [-0.57,0.64]
Re(CP-CP^′)	-	[0.05,0.20], [0.01,0.26]	[0.00,0.16], [-0.07,0.22]
Im(CP-CP^′)	-	[-0.07,0.08], [-0.12,0.12]	[-0.07,0.09], [-0.14,0.16]
Re(CP+CP^′)	[-0.24,0.51], [-0.51,0.82]	-	[-0.12,0.52], [-0.42,0.78]
Im(CP+CP^′)	[-0.36,0.37], [-0.67,0.67]	-	[-0.37,0.29], [-0.68,0.57]

In the special case of real-valued couplings, would lead to rings [32] instead of circles in Fig. 2. Our results improve and extend previous bounds in the literature to the most general case of complex-valued couplings. For example they are a factor 2 to 5 more stringent than [33] and comparable to [34] once restricting to the simpler scenarios considered there.

SM-EFT-constrained scalar couplings

In the following we consider a scenario in which it is assumed that there is a sizable hierarchy between the electroweak scale and the new-physics scale, , and that the SM gauge symmetries are only broken at the electroweak scale. This results in the augmentation of the SM by dimension-six operators that respect the SM gauge group and are composed of SM fields only. Such a scenario becomes more and more viable for two reasons. The first is the discovery of a scalar resonance at the LHC in agreement with all requirements of the Higgs particle in the SM. The second is the steadily rising lower bound on the mass of new particles reported by ATLAS and CMS in various more or less specific models.

A nonredundant set of dimension-six operators of this effective theory (SM-EFT) that requires a linear realization of the electroweak symmetry was given in [35]. The matching of the SM-EFT to the effective theory of decays (2.1) at the scale of the order of the W-boson mass was performed for vector couplings in [36]. The matching of tensor (2.4) and scalar (2.3) operators [32] shows that SM gauge groups in conjunction with the linear representation impose the relationson scalar couplings and require tensor couplings to be suppressed to the level of dimension-eight operators. In consequence only two scalar couplings arise that scale as where denotes the scale of electroweak symmetry breaking.

It must be noted that the relations (4.2) are a consequence of embedding the Higgs in a weak doublet along with the Goldstone bosons. For example, choosing a nonlinear representation of the scalar sector allows for additional dimension-six operators in the according effective theory, such that the couplings are all independent and tensor operators have nonvanishing couplings already at dimension six [37].

Omitting for the sake of simplicity terms of order and , the couplings can be bound from [32]In a similar spirit, dropping terms of order and givesIn the SM-EFT no relations between and arise, so they are in general additional independent parameters. Here we find that destructive interference with contributions involving does not significantly alter the bounds on . The results of two fits are shown in Fig. 3. In the first fit, we set and include all constraints on and . In the second fit, we allow for and further include all constraints from Table 1. For both fits, all six 2D marginals of real and imaginary parts of vs. have nearly circular contours of equal size that contain the SM point at the 68 % level except for vs. where it is within the 95 % credible region. The regions hardly vary between the two fits.

Fig. 3

68 and 95 % contours of the 2D-marginalized distributions of scalar couplings in the scenario SM-EFT with all constraints in Table 1 (red) when marginalizing over nonzero . For comparison, we superimpose the corresponding contours from using only and all constraints (dashed blue contours) with fixed . The SM is indicated by the black diamond

Since we consider here complex-valued couplings the allowed regions are circles rather than rings as for the case of real-valued couplings [32]. Compared to those rings, the circles are smaller because the probability moves from the ring toward the center of the circle.

Tensor, scalar, and vector couplings

The most general fit of complex-valued tensor and scalar couplings in combination with vector couplings – the combination of Sect. 4.1 and Sect. 4.2 – yields bounds very similar to those in Tables 2 and 3. The changes are only small and in fact the bounds tend to be even more stringent by (10–20) % because tensor and scalar couplings can contribute to only constructively – see (3.6). As before, branching-ratio measurements of help to improve the constraints on . This demonstrates that even in the case of complex-valued couplings there is enough information in the data to bound all 20 real and imaginary parts.

Angular observables in

Now we discuss what the fits tell us about likely values of observables that have not been measured yet but have sensitivity to tensor and scalar couplings. In the decay, we again consider and as in Sect. 3 and additionally . We compute the posterior predictive distribution (see Appendix C) for each observable integrated over the low- bin  GeV and high- bin  GeV matching LHCb’s range. The distributions resemble Gaussians, thus we summarize them by their modes and smallest 68 % intervals in Table 4 comparing the SM (prior predictive, ) to three NP scenarios. In each, we allow for interference with the vector couplings and additionally vary only (Sect. 4.1), only (Sect. 4.2), and finally both tensor and scalar couplings (Sect. 4.4).

Table 4

The posterior predictive 68 % probability intervals of not-yet-measured angular observables for several new-physics scenarios given all the considered experimental constraints. The corresponding values for the SM (prior predictive) are given, too, where “” indicates zero in the considered approximation, see text for details

Observable	q^2-bin [GeV^2]	SM	T(5),9,10	S^(^′),P^(^′),9^(^′),10^(^′)	S^(^′),P^(^′),T(5),9^(^′),10^(^′)
τ_B^0×J_6c	[1.1,6]	≃0	(0.6-1.9+1.8)×10^-9	(-0.1-1.9+2.3)×10^-10	(0.2-1.1+1.1)×10^-9
	[15,19]	≃0	(2.1-6.2+5.1)×10^-10	(0.7-6.6+5.4)×10^-11	(-0.2-2.9+3.5)×10^-11
τ_B^0×(J_1c+J_2c)	[1.1,6]	(3.30-0.56+0.65)×10^-9	(4.8-1.9+1.7)×10^-9	(2.6-0.4+0.5)×10^-9	(2.8-1.1+1.5)×10^-9
	[15,19]	(1.72-0.16+0.16)×10^-10	(4.7-2.6+3.3)×10^-9	(1.9-0.7+0.3)×10^-10	(3.1-2.3+2.5)×10^-9
τ_B^0×(J_1s-3J_2s)	[1.1,6]	(2.44-0.48+0.46)×10^-10	(1.9-1.1+1.4)×10^-8	(4.0-1.2+1.6)×10^-10	(1.5-1.1+1.0)×10^-8
	[15,19]	(1.12-0.10+0.10)×10^-10	(7.1-3.8+5.6)×10^-9	(1.4-0.4+0.3)×10^-10	(5.3-3.8+3.8)×10^-9
A_ΔΓ(B_s→μ¯μ)	–	1	1	[-1,1]	[-1,1]

We rescale and combinations by the -meson life time  ps [38] to judge the experimental sensitivity in the near future by comparing to current measurements of the branching ratio,In the SM, the typical magnitude of the branching ratio of is approximately equal to for both the and the GeV bins. For comparison, the predicted ranges for and in the SM are suppressed by orders of magnitude down to ; cf. Table 4.

The angular observable is strictly zero in the absence of tensor and scalar couplings. Nonzero contributions can be generated in the SM by QED corrections or potentially from higher-dimensional () operators, leading to parametric suppression by or . These factors should be compared to the potential suppression present for tensor and scalar contributions in particular NP models in order to gauge their relevance. Our model-independent fits are still in a regime where such considerations are insignificant since current experimental measurements, in combination with theory uncertainties, do not yet impose sufficiently stringent constraints on tensor and scalar couplings.

Beyond the SM, can become of order in scenarios involving tensor couplings only and about in the presence of scalar couplings only. Both effects are due to interference with vector couplings. Schematically, is a function of , , and . The largest interval is obtained for the scenario without scalar couplings because then the uncertainty on the (tensor) couplings is largest. But even then, it seems that the experimental sensitivity will not be high enough to have an impact in global fits.

Concerning and , substantial deviations from the SM prediction are again only possible in the presence of nonzero tensor couplings. In this case, an enhancement by two orders of magnitude is possible up to at high and also at low in the case of . We want to stress again that we make these statements conditional on all included experimental constraints, the scenario, and our prior. In view of the current experimental precision of on the branching ratio at LHCb [21] with only 1 fb, corresponding to the , one can indeed hope for some sensitivity to such large effects in and for the not-yet-published 3 fb data set. At least, we can hope for some measurement if the method of moments [7] is applied.

For the mass-eigenstate rate asymmetry induced by the nonvanishing width of the meson (cf. Appendix A), we find a rather uniform distribution in scenarios with nonzero scalar couplings. So any value in the range is plausible whereas the SM and the scenario with only tensor couplings predict a value of precisely one [39]; cf. the last row in Table 4. Hence any deviation from one would unambiguously hint at the presence of scalar operators.

Conclusions

We have derived the most stringent constraints to date on tensor and scalar couplings that mediate transitions. They are based on the latest measurements of angular observables and the lepton forward–backward asymmetry in from LHCb [1], supplemented by measurements of the branching ratios of and .

Both and belong to a class of observables in which vector and dipole couplings – present in the standard model (SM) – are suppressed (mostly kinematically by ) with respect to tensor and scalar couplings. We provide predictions for the equivalent but not-yet-measured angular observables , , and in .

In a Bayesian analysis of the complex-valued couplings of the effective theory, we find that the measurement of , especially at high-,Our updated bounds on complex-valued tensor and scalar couplings are summarized in Tables 2 and 3, accounting also for interference effects with vector couplings. These bounds hold even in the most general scenario of complex-valued tensor, scalar, and vector couplings, showing that the data are good enough to bound the real and imaginary parts of all Wilson coefficients simultaneously.

imposes by itself constraints on real and imaginary parts of the tensor couplings such that the upper bound of the smallest 68 % (95 %) credibility interval is , superseding previous bounds. In combination with current data from and lattice predictions of form factors, the bounds are lowered to , even in the presence of nonstandard contributions in vector couplings.

for the first time allows one to simultaneously bound all four scalar couplings due to its complementarity to . Even when taking into account destructive interference with vector couplings, all credibility intervals contain the SM. The corresponding 1D intervals of real and imaginary parts have widths of 0.6 (1.2) for and 0.15  (0.3) for for . Currently, the bounds from are weaker than those from by about a factor of 4. Future measurements of at LHCb and Belle II will further tighten the bounds. Moreover, measurements of () will provide constraints on scalar couplings in the electron channel in the absence of a direct determination of the branching ratio of .

As a special case, we consider the scenario arising from the SM augmented by dimension-six operators generalizing existing studies to the case of complex-valued couplings. In this scenario, tensor couplings are absent and additional relations between scalar couplings are enforced by the linear realization of the electroweak symmetry group.

Our study of the yet unmeasured angular observables , , and in (see Table 4) shows that despite the current bounds on tensor couplings, enhancements of up to two orders of magnitude over the SM predictions are allowed for and , placing them in the reach of the LHCb analysis of the full run I data set. Our bounds on scalar couplings from and , however, are already quite restrictive, permitting only small deviations from SM predictions in and . Notably, the mass-eigenstate rate asymmetry given nonzero scalar couplings can take on any value in the range .

Table 5

Prior distributions of the nuisance parameters for hadronic quantities

Quantity	Prior	References
B_s decay constant
f_B_s	(227.7±4.5) MeV	[45]
B→K form factors
f_+(0)=f_0(0)	0.34±0.05	[51, 52]
f_T(0)	0.38±0.06	[51, 52]
b10	-4.3-0.9+0.8	[51]
b1+	-2.1-1.6+0.9	[51]
b1T	-2.2-2.0+1.0	[51]
B→K^∗ form factors
α0A_0	0.35-0.03+0.02	[23, 25]
α1A_0	-1.21-0.08+0.08	[23, 25]
α2A_0	0.77-0.48+0.45	[23, 25]
α0A_1	0.2625-0.015+0.015	[23, 25]
α1A_1	0.07-0.08+0.08	[23, 25]
α2A_1	0.045-0.15+0.12	[23, 25]
α1A_12	0.53-0.14+0.10	[23, 25]
α2A_12	0.32-0.39+0.33	[23, 25]
α0V	0.35-0.03+0.02	[23, 25]
α1V	-1.19-0.08+0.08	[23, 25]
α2V	1.55-0.33+0.33	[23, 25]
α0T_1	0.31-0.02+0.02	[23, 25]
α1T_1	-1.07-0.08+0.08	[23, 25]
α2T_1	1.40-0.30+0.30	[23, 25]
α1T_2	0.33-0.10+0.08	[23, 25]
α2T_2	0.26-0.24+0.21	[23, 25]
α0T_23	0.67-0.05+0.04	[23, 25]
α1T_23	0.93-0.22+0.28	[23, 25]
α2T_23	-0.15-0.71+0.86	[23, 25]

Table 6

Mean values, standard deviations (top) and correlation coefficients (bottom) of lattice points [13] for the form factors

q^2 (GeV ^2)	17	20	23
f_0(q^2)	0.62±0.03	0.72±0.03	0.87±0.04
f_+(q^2)	1.13±0.05	1.63±0.07	2.68±0.13
f_T(q^2)	1.02±0.06	1.47±0.08	2.42±0.18