Next Article in Journal
Small-Signal Modeling and Frequency Support Capacity Analysis of Power Load Considering Voltage Variation Effect
Next Article in Special Issue
Short-Distance Physics with Rare Kaon Decays
Previous Article in Journal
RETRACTED: Sun et al. Application of Symmetry Law in Numerical Modeling of Hydraulic Fracturing by Finite Element Method. Symmetry 2020, 12, 1122
Previous Article in Special Issue
Lepton Flavour Universality in Rare B Decays
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

Current Status of the Standard Model Prediction for the Bsμ+μ Branching Ratio

Institute of Theoretical Physics, Faculty of Physics, University of Warsaw, PL-02-093 Warsaw, Poland
*
Author to whom correspondence should be addressed.
Symmetry 2024, 16(7), 917; https://doi.org/10.3390/sym16070917
Submission received: 29 June 2024 / Revised: 10 July 2024 / Accepted: 12 July 2024 / Published: 18 July 2024
(This article belongs to the Special Issue Symmetries and Anomalies in Flavour Physics)

Abstract

:
The rare decay B s μ + μ provides an important constraint on possible deviations from the Standard Model in b-s-- interactions. The present weighted average of its branching ratio measurements amounts to ( 3.34 ± 0.27 ) × 10 9 , which remains in good agreement with the theoretical prediction of ( 3.64 ± 0.12 ) × 10 9 within the Standard Model. In the present paper, we review calculations that have contributed to this prediction and discuss the associated uncertainties.

1. Introduction

The Higgs boson discovery [1,2] through direct production at the LHC completed the experimental search for the Standard Model (SM) particle content. Since then, no clear signal for Beyond-Standard Model (BSM) particle production has been seen at the high-energy frontier of experimental particle physics. Consequently, more and more focus is being shifted towards precise studies of rare processes that are sensitive to corrections from BSMs.
In particular, decays of the B meson mediated through Flavor-Changing Neutral Currents (FCNCs) have been a very active area of research. Theoretical studies of the B mesons are greatly aided by the framework of Heavy Quark Expansion which, in many cases, allows us to parameterize the effects of their hadronic structure through a series of non-perturbative matrix elements suppressed by powers of Λ Q C D / M B . Moreover, while the FCNC-mediated processes are loop-suppressed in the SM, they can receive sizable tree-level contributions from BSMs, which underlines them as primary candidates for observables where indirect signals from new particles may be detected [3]. The rich phenomenology of rare B meson decays is being actively studied by LHC experiments and Belle II that follow and extend past investigations at CLEO, Belle, and BABAR. To fully take advantage of the current and future experimental data, improvements in the precision of the SM predictions are often necessary.
One of the most interesting rare decays of the B meson and the focus of this review is the B s μ + μ channel that was first observed over a decade ago [4]. Since then, the experimental precision for its average time-integrated branching ratio B ¯ s μ has reached O ( 10 % ) [5,6,7,8]. The current world average reads [9]
B ¯ s μ exp = ( 3.34 ± 0.27 ) × 10 9 .
The similar B d μ + μ channel is suppressed by a factor | V t d / V t s | 2 0.04 , which has a negative effect on experimental precision due to significantly reduced statistics, as illustrated in Figure 1.
The SM description of B s μ + μ is greatly simplified by a factorization of long- and short-distance Quantum Chromodynamics (QCD) effects. In contrast to many other B decays, the SM decay amplitude of B s μ + μ depends, to a very good approximation, only on a single non-perturbative hadronic quantity, namely the B s -meson decay constant f B s . Its determinations from lattice simulations (see below) are mature and precise. Moreover, the hard QCD and electroweak (EW) corrections are, to a very good approximation, contained within a single Wilson coefficient C A , which can be calculated in a standard, perturbative matching procedure as a series in the strong and electromagnetic couplings α s and α e . Thanks to these properties,  SM calculations have reached the accuracy of a few percent. The current SM prediction for the branching ratio reads
B ¯ s μ S M = ( 3.64 ± 0.12 ) × 10 9 .
We describe its evaluation in the next sections. The numerical value in Equation (2) has been obtained by updating the input parameters in the semi-numerical expressions of Ref. [10] and including the power-enhanced QED correction from Refs. [11,12] that amounts to around 0.5 % . A difference with respect to ( 3.66 ± 0.14 ) × 10 9 in Ref. [12] is due to the parameter update. The above SM prediction is well in agreement with each individual measurement [5,6,7,8], as well as with their average in Equation (1).
In the SM, on top of the aforementioned FCNC loop suppression, B ¯ s μ receives a helicity suppression by the mass-squared ratio m μ 2 / m B s 2 . One or both of these suppressions may be lifted in models with additional Higgs doublets or with a  Z . In effect, one finds restrictions on allowed parameter spaces in, among others, Two Higgs Doublet Models (2HDMs) [13] and the Minimal Supersymmetric Standard Model [14]. In addition, several time-dependent observables in B s μ + μ can be used to study potential BSM CP -violation mechanisms [15].
The present review is organized as follows. In Section 2, we present the effective Lagrangian and the branching ratio formula for B s μ + μ in the SM. Section 3 and Section 4 are devoted to describing, respectively, the perturbative QCD and EW corrections to the Wilson coefficient C A . At the end of Section 4, the power-enhanced QED correction to B ¯ s μ is discussed. In Section 5, we summarize the current parameter update and evaluate the SM prediction for B ¯ s μ  (2) together with the corresponding uncertainty. We conclude in Section 6. In Appendix A, we recall the derivation of the branching ratio formula. In Appendix B, we present the current SM prediction for B ¯ d μ , i.e., the average time-integrated branching ratio of B d μ + μ .

2. The Effective Lagrangian and the Branching Ratio Formula

The effective Lagrangian used to describe B s l + l ,   l { e , μ , τ } in the SM is obtained through simultaneous integrating out of all fields heavier than the b quark at the scale μ 0 = O ( m t ) . It has the form
L = L QCD × QED ( fields lighter than W ) + N n C n Q n + h . c . ,
where
N V t b V t s G F 2 M W 2 π 2 ,
V i j are the Cabibbo–Kobayashi–Maskawa (CKM) matrix elements, G F is the Fermi constant, while M W is the W-boson mass defined in the on-shell renormalization scheme. The local operators Q n are polynomials in the light fields and their derivatives. The Wilson coefficients C n can be treated as real-valued (up to negligible corrections) once the global normalization factor N is set as in Equation (4). The operators Q n are of mass-dimension 5 or higher and have to be suppressed by powers of 1 / M W necessary to make the overall mass-dimension of L equal to 4. At the leading order in 1 / M W and α e , it is sufficient to consider the operators Q n , where a Δ B = Δ S = 1 flavor-changing quark current multiplies a lepton current. Moreover, the quark current must violate parity to annihilate the pseudoscalar B s meson. Once the Lorentz invariance is imposed, one is left with the following four operators only (no tensor current  b ¯ σ α β γ 5 s   needs to be considered because no antisymmetric tensor could be formed from p α alone in the corresponding matrix element, analogous to the one in Equation (7)):
Q A [ l ¯ γ α γ 5 l ] [ b ¯ γ α γ 5 s ] [ l ¯ γ α γ 5 l ] j A α , Q V [ l ¯ γ α l ] j A α , Q P [ l ¯ γ 5 l ] [ b ¯ γ 5 s ] [ l ¯ γ 5 l ] j P , Q S [ l ¯ l ] j P .
The Lagrangian (3) can be used to derive the formula for B ¯ s μ . A sketch of the derivation is given in Appendix A. While evaluating the contribution of Q A there, it becomes clear that Q V does not contribute at the leading order in α e . As far as Q P and Q S are concerned, their Wilson coefficients are computed in a matching to full SM amplitudes with the quark and lepton currents exchanging Higgs bosons. Such contributions are suppressed by m b 2 / M W 2 and can be neglected as being of the same order as dimension-8 operator effects. Hence, neglecting these tiny effects, B ¯ s μ in the SM depends on C A only. The explicit expression reads
B ¯ s μ S M = | N | 2 M B s 3 f B s 2 8 π Γ H s β r 2 | C A | 2 + O α e , m b 2 M W 2 ,
where r 2 m μ / M B s , β 1 r 2 , while Γ H s is the heavier mass eigenstate width in the B s B ¯ s system. Finally, the  B s decay constant f B s is defined through the relation
0 | j A α ( x ) | B s ( p ) i p α f B s e i p x .
It is calculated using lattice QCD methods with errors at a sub-percent level. The current world average based on 2 + 1 + 1 simulations [16,17,18,19] alone amounts to [20]
f B s = ( 230.3 ± 1.3 ) MeV .
In several popular BSMs, the Wilson coefficients C S and C P can become comparable in size to r C A . For example, in the 2HDM-II with large tan β H v 2 / v 1 (the ratio of the vacuum expectation values of the two Higgs doublets), one finds [21],
C S C P ln ρ ρ 1 m μ m b 4 M W 2 tan 2 β H ,
where ρ M H + 2 / m t 2 . The SM suppression factor of m b 2 / M W 2 can be compensated by values of tan β H = O ( 50 ) or larger. The branching ratio formula in the 2HDM becomes (see Appendix A)
B ¯ s μ 2 H D M = | N | 2 M B s 3 f B s 2 8 π Γ H s β | r C A u C P | 2 + Γ H s Γ L s | u β C S | 2 + O ( α e ) ,
where u M B s / ( m b + m s ) , and Γ L s is the lighter mass eigenstate width in the B s B ¯ s system.
Coming back to the SM expression for B ¯ s μ in Equation (6), several comments concerning the O ( α e ) terms are necessary. First, all the corrections of this order to C A are already known and included in the numerical result in Equation (2). They will be discussed in Section 4.1. Some of them are enhanced by 1 / sin 2 θ W , m t 2 / M W 2 or ln 2 M W 2 / m b 2 , where θ W is the weak mixing angle. Second, some of the remaining electromagnetic corrections in the O ( α e ) term in Equation (6) receive a power enhancement by M B s / Λ QCD [11]. They come from virtual photons emitted by the s quark and absorbed by the leptons. We shall comment on them in Section 4.2. Their evaluation requires extending the operator basis to include the effects of other dimension-6 operators, not present in Equation (5), such as
Q 2 [ b ¯ γ α P L c ] [ c ¯ γ α P L s ] or Q 7 e m b 16 π 2 [ b ¯ σ α β P L s ] F α β ,
where P L = ( 1 γ 5 ) / 2 . Third, the dependence of C A on the renormalization scale μ b (that arises due to QED effects only) induces around O ( 0.3 % ) uncertainty in B ¯ s μ S M [10]. Such a dependence must be compensated by the yet unknown O ( α e ) corrections that receive no extra enhancement factors.

3. The QCD Corrections to C A

As described in the previous section, at the leading order in α e and m b 2 / M W 2 , the only significant short-distance parameter entering the SM branching ratio formula is the C A ( μ b ) Wilson coefficient, where we specifically indicate its dependence on the low-energy renormalization scale μ b = O ( m b ) . It is evaluated by demanding the equality of corresponding Green’s functions in the SM and the effective theory (3) at the renormalization scale μ 0 = O ( m t ) at which the electroweak degrees of freedom are decoupled. Computationally, the simplest choice of a Green’s function for the extraction of C A is G [ b ¯ , s , l ¯ , l ] with vanishing external momenta:
G S M [ b ¯ , s , l ¯ , l ] ( μ 0 ) | p i = 0 = G eff [ b ¯ , s , l ¯ , l ] ( μ 0 ) | p i = 0 .
At the matching scale μ 0 , QCD on both sides is treated perturbatively. This equation is then solved for C A ( μ 0 ) order-by-order in α s and α e , resulting in a double series
C A ( μ 0 ) = m , n = 0 α ˜ s m α ˜ e n C A ( m , n ) ( μ 0 ) ,
where α i ( μ 0 ) 4 π α ˜ i are the running couplings renormalized in the MS ¯ scheme at the scale μ 0 . In this section, we focus on the leading terms in α e , namely C A ( m , 0 ) .
As we work at the leading order in 1 / M W , all light masses in the matching Equation (12) can be set to 0. In dimensional regularization, all scaleless loop integrals vanish. In consequence, on the effective theory side, one is left only with tree diagrams, including the UV-counterterm ones. On the SM side, partially massive tadpoles have to be calculated. Removing light masses from propagators in the SM Green’s function leads to spurious infrared divergences in loop integrals evaluated in d = 4 2 ϵ dimensions. The resulting additional ϵ -poles are not removed by the renormalization constants of the SM. Instead, they cancel in the matching Equation (12) against the tree-level UV counterterms in the effective theory.
Once such a procedure is applied in our case, one has to supplement the effective Lagrangian with additional operators that vanish when ϵ 0 due to the Dirac algebra identities. Such operators are called evanescent. The UV-counterterms with these operators cancel against some of the spurious infrared divergence effects on the SM side. Moreover, their Wilson coefficients evaluated at lower orders affect the physical operator Wilson coefficients at higher orders. For the C A ( m , 0 ) terms, it is sufficient to introduce only one evanescent operator [22]:
Q A E = [ b ¯ γ α γ β γ δ s ] [ l ¯ γ δ γ β γ α l ] 4 Q A .
The bare fields, couplings, and Wilson coefficients on the effective side are replaced by the renormalized ones, with mixing occurring between the physical and evanescent operators:
C A ( b ) Q A ( b ) + C A E ( b ) Q A E ( b ) = Z q Z l ( C A Z N N Q A + C A Z N E Q A E + C A E Z E N Q A + C A E Z E E Q A E ) ,
where Z q and Z l are the quark- and lepton-field MS ¯ renormalization constants, respectively, with  Z l = 1 + O ( α e ) . To fix the renormalization constants Z i j , one demands that in the renormalized Green’s functions, the terms proportional to C A are finite when ϵ 0 , while those proportional to C A E vanish in this limit [23]. Examples of diagrams contributing to Z i j at α s and α s 2 are shown in Figure 2. Results for the relevant Z i j up to O ( α e 0 α s 2 ) are given in Equation (15) of Ref. [24].
It is important to emphasize that Z N N = 1 to all orders in QCD, once higher-dimensional operators and O ( α e ) effects are neglected. Therefore, the Renormalization Group Equation (RGE) for C A is trivial at this level:
μ d d μ C A = O ( α e ) .
Consequently, there is no RG evolution of C A ( μ ) at the leading order in α e :
C A ( μ b ) = C A ( μ 0 ) + O ( α e ) .
The SM Green’s function in Equation (12) receives QCD corrections from two classes of diagrams: W-boxes and Z-penguins, shown in Figure 3. Such bare diagrams are renormalized using lower-loop SM diagrams with counterterms. The QCD coupling constant renormalization Z g has to be modified by an appropriate threshold correction to match the coupling constant of the effective theory with only five active flavors (see Section 3 in Ref. [25]). Similarly, the light-field wave-function renormalization constants have to be shifted to account for the decoupling threshold (see Section 4 in Ref. [26]). Heavy fields and masses are renormalized in the MS ¯ scheme.
Once both sides of the matching Equation (12) have been properly renormalized, the value of C A ( μ 0 ) can be extracted. All the Dirac structures appearing on the SM side are mapped onto either Q A or Q A E and their coefficients compared to the ones on the effective side.
The procedure described in this section was completed in Ref. [27] at the leading order, followed by Refs. [22,28,29] and [24] for the O ( α s 1 ) and O ( α s 2 ) corrections, respectively. In the latter case, high-order expansions in y M W / m t and w 1 M W 2 / m t 2 were computed, and their combination was used to obtain a numerical result at the physical value of M W / m t .

4. The Electroweak Corrections

4.1. The C A ( 0 , 1 ) Correction

The next-to-leading order EW correction C A ( 0 , 1 ) ( μ 0 ) was first obtained in Ref. [30] (Let us note that it was called c 10 ( 2 , 2 ) in that paper due to different notational conventions). Its calculation follows a similar pattern to the one described in the previous section. The main difference comes from additional Dirac structures appearing in the O ( α e ) corrections on the SM side of the matching equation, which necessitate including additional operators in the effective Lagrangian. For the complete matching at the scale μ 0 , at the leading order in 1 / M W and including O ( α e ) terms, one has to retain all the operators from Equation (5), and supplement them with Q 2 defined in Equation (11). The cancellation of spurious infrared divergence effects requires also the inclusion of additional evanescent operators (see Appendix A of Ref. [30]). The renormalization constants of all resulting Wilson coefficients are then calculated in a way analogous to the C A ( m , 0 ) calculation.
Examples of diagrams contributing to C A ( 0 , 1 ) ( μ 0 ) on the SM side are shown in Figure 4. Before this correction can be extracted, the UV divergences on the SM side have to be renormalized. For a discussion of different renormalization schemes for the electroweak boson and top-quark masses, we refer to the original article [30]. Here, we will continue the analysis in the OS-2 scheme defined therein, as it was used in the subsequent phenomenological analysis [10]. In this scheme, all the QCD corrections to mass renormalization constants are defined in the MS ¯ , but the O ( α e ) corrections to Z m t , Z M W , and Z M Z are defined on-shell. In practice, the calculation is first performed fully in the MS ¯ scheme, and the finite terms in these three renormalization constants are subsequently added to the renormalized results.
The value of C A extracted at the scale μ 0 must then be RG-evolved down to the μ b scale of the B s μ + μ decay. The renormalized Wilson coefficients of the extended set of effective operators are related to the bare ones through a relation similar to Equation (15):
C j ( b ) Q j ( b ) = Z q Z l k C k Z k j Q j C j ( b ) = k C k Z k j .
The one-loop RGE for Wilson coefficients C has a general form
μ d d μ C = α ˜ s γ s ( 0 ) + α ˜ e γ e ( 0 ) T C ,
where the Anomalous Dimension Matrices (ADMs) γ s ( 0 ) and γ e ( 0 ) are obtained from the renormalization constants of Wilson coefficients. Once we restrict to non-evanescent operators only, the necessary MS-scheme relation reads
Z k j = δ k j + 1 2 ϵ α ˜ s γ s ( 0 ) + α ˜ e γ e ( 0 ) k j + O α ˜ s 2 , α ˜ e 2 , α ˜ s α ˜ e .
For the RGE (19) to close, one has to extend the set of operators in the effective Lagrangian further (see the discussion under Equation (15) in Ref. [30]). Analytical expressions for the evolution operator associated with Equation (19) can be found, e.g., in Ref. [31]. The details and results of the numerical solution can be found in Appendix B of Ref. [30].

4.2. Power-Enhanced QED Corrections

Within the effective theory framework, the S-matrix element corresponding to the B s μ + μ decay receives O ( α e ) contributions from diagrams with tree-level Wilson coefficients and a virtual photon exchanged between the fermions. In particular, there is a class of diagrams with a photon emitted from the spectator s quark and absorbed by one of the outgoing muons (see Figure 5). The leading terms of the power series in { m s , m μ , Λ QCD } / m b of this correction were computed in Ref. [11]. It was observed that its helicity suppression by the factor of r 2 in the tree-level branching ratio (6) is partially lifted due to a relative enhancement by M B s / Λ QCD .
Effects of hard virtual photons with energies and virtualities above the μ b scale are contained in the Wilson coefficients C i . The remaining virtual photons need to be taken into account in the physical matrix elements that are evaluated at the scale μ b . In particular, virtual photons exchanged in the diagrams in Figure 5 probe the inner structure of the B s meson, smearing the annihilation point of the valence quarks over a distance O ( ( M B s Λ QCD ) 1 2 ) , which corresponds to inverse virtuality of the off-shell s quark. Such long-distance QCD effects cannot be parameterized solely by f B s , as at the leading order. Instead, one has to estimate the effects of matrix elements like
0 | d d x T [ [ μ ¯ γ α μ ] ( x ) Q A ( 0 ) ] | B s ,
where T is the time-ordering operator. They involve the B meson light-cone distribution amplitude [32] and its first logarithmic moments.
Virtual photon exchanges leading to power-enhanced QED corrections were thoroughly studied in the formalism of the Soft-Collinear Effective Theory by Beneke, Bobeth, and Szafron in Ref. [12]. Their effect can be included in the SM prediction through a replacement
B ¯ s μ η B B S B ¯ s μ ,
with
η B B S = 0 . 995 0.005 + 0.003 .
The main uncertainty in η B B S comes from poorly known values of hadronic parameters [33]. We have extracted the numerical value of η B B S as well as its uncertainty in Equation (23) from Equations (8.8) and (8.10) of Ref. [12]. One can observe there that the relatively small central value of the correction ( 0.5 % ) arises as an effect of partial cancellation between potentially unrelated contributions from the Q V and Q 7 operators. Thanks to this cancellation, the overall effect remains below the ± 1.5 % non-parametric uncertainty estimated in Ref. [10]. That uncertainty was primarily due to unknown O ( α e ) effects, although the possibility of their power enhancement remained unknown at the time.
The contribution of Q 7 to the considered corrections was reanalyzed in Ref. [34]. Despite several differences in the analytical expressions with respect to the earlier analysis of Ref. [12], the numerical results of the two papers remain in qualitative agreement, and no modification of the factor η B B S in Equation (23) is necessary.

5. Numerical Analysis

In this section, we update the SM prediction for B ¯ s μ based on Equation (6), including the complete O ( α s 2 ) and O ( α e ) corrections to C A , as well as the QED correction factor in Equation (23). In practice, it is sufficient to use the semi-numerical expressions from Equation (6) of Ref. [10], which gives
B ¯ s μ S M × 10 9 = ( 3.65 ± 0.06 ) M t [ GeV ] 173.1 3.06 α s ( M Z ) 0.1184 0.18 R s η B B S ,
where
R s = f B s [ MeV ] 227.7 2 | V c b | 0.0424 2 | V t b V t s / V c b | 0.980 2 τ H s [ ps ] 1.615 .
In the above expressions, all the explicitly displayed input parameters are normalized to their 2013 central values. In Table 1, we update these central values together with the corresponding uncertainties. All the remaining parameters are retained unaltered with respect to Table I of Ref. [10], as their update would not affect B ¯ s μ in a noticeable manner—they are either very precisely measured or have little effect on B ¯ s μ .
As already mentioned in Equation (2), we find B ¯ s μ S M = ( 3.64 ± 0.12 ) × 10 9 . The overall uncertainty is now almost a factor of two smaller than that found in Ref. [10], while the central value remains almost unchanged. The latter fact can be attributed to an approximate cancellation of shifts stemming from the parameter updates and η B B S , as in the following sum:
+ 2.3 % ( f B s ) 1.6 % ( C K M ) + 0.4 % ( τ H s ) 0.9 % ( M t ) + 0.1 % ( α s ) 0.5 % ( η B B S ) 0.3 % .
As far as the uncertainty breakdown is concerned, its current version is compared to the 2013 one [10] in Table 2. In its last column, the uncertainties are combined in quadrature.
One can observe a significant improvement in the first four columns where the dominant parametric uncertainties originate from, in particular in the case of f B s , which is determined on the lattice. As far as the top-quark mass M t is concerned, let us recall that the PDG [9] value is treated as the on-shell mass, which is not strictly correct. However, the overall non-parametric uncertainty of ± 1.5 % is understood to contain a contribution from such an approximation. As is evident from Table 1 and Table 2, a  300 MeV shift in M t implies a  0.5 % shift in B ¯ s μ S M .
At present, the most important uncertainty originates from | V c b | , in which case we use the inclusive determination only [35]. A combination with exclusive determinations would not lead to an improvement, given the persistent discrepancy between the inclusive and exclusive results [9]. Our preference is the same as in Ref. [10], i.e., we treat the inclusive determination as theoretically cleaner and more reliable.
As already mentioned in Ref. [10], one can get rid of | V c b | in the ratio of B ¯ s μ to the measured B s ( H ) - B s ( L ) mass difference, at the cost of introducing an extra uncertainty from lattice determinations of the “bag parameter” B B s . Such an approach is likely to become relevant once the experimental accuracy in B ¯ s μ (currently 8.1 % in Equation (1)) becomes closer to the theoretical one.

6. Conclusions

In this review, we analyzed the current SM prediction for the branching ratio of the rare B s μ + μ decay. This channel continues to be among the most promising candidates for detecting BSM physics without the direct production of new particles, due to its SM suppression and possible BSM enhancements.
The SM analysis is, to a very good approximation, contained in the perturbative calculation of the Wilson coefficient C A and the lattice calculation of the long-distance QCD parameter f B s . The perturbative part was already complete up to and including next-to-next-to-leading QCD and next-to-leading EW effects.
The currently dominant theoretical uncertainty originates from the CKM-matrix element | V c b | . The next-to-dominant uncertainty is already non-parametric, stemming mainly from the unknown higher-order electromagnetic corrections at the scale μ b . They depend on non-perturbative effects that are not contained in f B s . The current result for B ¯ s μ S M changes by around 0.3 % when μ b is varied between m b / 2 and 2 m b . However, since it provides a lower bound only on the possible size of unknown electromagnetic effects, the actual uncertainty estimate should be somewhat more conservative. Here, we have retained the overall non-parametric uncertainty at the same level as in Ref. [10], namely ± 1.5 % .
The experimental error is currently much larger, around 8%, which sets a limit on the power of B s μ + μ as a means for testing various BSM theories. The situation is expected to improve with time, when higher statistics are collected at the LHC and future experiments.
A final message that we would like to share with the reader is as follows. Our numerical result in Equation (2) will become outdated as soon as any of the input parameters are determined in a new analysis. Performing another update of B ¯ s μ S M does not require being an expert. It is sufficient to substitute new inputs into the simple formula (24).

Author Contributions

Both authors have equally contributed. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Science Center, Poland, under the research project 2020/37/B/ST2/02746.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. The Branching Ratio Formula

In this appendix, we sketch the derivation of the branching ratio formula (10) that holds in models with SM-like CP -violation, including the SM and the 2HDM. We work at the leading order in QED throughout, i.e., the final-state muons are understood to be emitted directly from the operator vertex ( Q A , Q P , or Q S ).
Let | B s and | B ¯ s denote the meson flavor eigenstates with valence quarks b ¯ s and b s ¯ , respectively. We fix conventions in their overall phases by demanding that CP   | B s = | B ¯ s , and  CPT | B s = | B ¯ s . Once this is done, the heavier (H) and lighter (L) mass eigenstates (see Section 13 of Ref. [9]) can, respectively, be written as
| B s ( H ) = 1 2 | N | N | B s N | B ¯ s , | B s ( L ) = 1 2 | N | N | B s + N | B ¯ s ,
where N has been defined in Equation (4). In the limit of no CP -violation ( N = | N | ), B s ( H ) and B s ( L ) are CP -odd and CP -even, respectively.
From the form of the lepton currents in Equation (5), we observe that the Q A and Q P interactions can lead to production of CP -odd lepton pairs only (in the CM frame), while Q S can lead to production of CP -even pairs only. In the Q A case, it follows from the fact that the timelike component of the lepton current is CP -odd (see, e.g., the table below Equation (3.150) of Ref. [38]), while the spacelike components play no role, as they get contracted with the vanishing spacelike components of the meson momentum (see Equation (7)).
Let us now show that even in the presence of SM-like CP -violation, ( N | N | but real Wilson coefficients), the operators Q A and Q P have no effect on dimuonic decays of B s ( L ) , while Q S has no effect on such decays of B s ( H ) . When the “ h . c . ” terms in the Lagrangian (3) are taken into account, the leading Q A contribution to the B s ( L ) decay amplitude is proportional to
μ + μ | N Q A + N Q A | B s ( L ) = | N | 2 μ + μ | Q A | B s + μ + μ | Q A | B ¯ s ,
where the matrix elements on the r.h.s. are the only two that do not vanish due to flavor conservation in QCD. Next, we observe that
μ + μ | Q A | B ¯ s = μ + μ | ( CP ) Q A CP | B ¯ s = μ + μ | Q A | B s ,
where we have taken advantage of the fact that the dimuon state is CP -odd, as argued above. Consequently, the sum of the two matrix elements on the r.h.s. of Equation (A2) vanishes. An identical reasoning holds for Q P . As far as Q S and B s ( H ) are concerned, we proceed by analogy:
μ + μ | N Q S + N Q S | B s ( H ) = | N | 2 μ + μ | Q S | B s μ + μ | Q S | B ¯ s ,
μ + μ | Q S | B ¯ s = μ + μ | ( CP ) Q S CP | B ¯ s = μ + μ | Q S | B s ,
where, this time, the dimuon state is CP -even. Consequently, the difference in the two matrix elements on the r.h.s. of Equation (A4) vanishes. As a by-product of the above reasoning, we can simplify the non-vanishing matrix elements as follows:
μ + μ | N Q A , P + N Q A , P | B s ( H ) = | N | 2 μ + μ | Q A , P | B s , μ + μ | N Q S + N Q S | B s ( L ) = | N | 2 μ + μ | Q S | B s .
Both at the LHC and at e + e machines, the production rates of B s and B ¯ s are practically equal. Thus, to a very good approximation, the heavy and light mass eigenstates are produced in the same quantities, as can be seen by inverting the relation (A1). Since the decay products of B s ( H ) ( CP -odd dimuons) and B s ( L ) ( CP -even dimuons) do not interfere, the average time-integrated branching ratio is simply given by
B ¯ s μ = 1 2 Γ [ B s ( H ) μ + μ ] Γ H s + Γ [ B s ( L ) μ + μ ] Γ L s .
In each of the two cases, the decay rate is given by the well-known formula
Γ [ B s ( H , L ) μ + μ ] = 1 2 M B s d P S 2 M ¯ ( H , L ) 2 = β 16 π M B s M ¯ ( H , L ) 2 ,
with
d P S 2 = 1 4 π 2 d 4 k + d 4 k δ ( k + 2 m μ 2 ) θ ( k + 0 ) δ ( k 2 m μ 2 ) θ ( k 0 ) δ ( 4 ) ( p k + k ) .
Here, M ( H , L ) are the corresponding invariant matrix elements, β has been defined below Equation (6), and we have neglected the tiny mass splitting between B s ( H ) and B s ( L ) . Summing over spins of the final-state muons is understood in Equation (A8). The two-body phase-space integral is trivial (see Section 3.2 of Ref. [39]), as M ¯ ( H , L ) 2 is constant in the integration domain due to rotational symmetry in the decaying scalar rest frame.
The Q A (and Q A ) contribution to M ( H ) reads
M A ( H ) = i C A | N | 2 e i p x 0 | j A α ( x ) | B s ( p ) u ¯ ( k ) γ α γ 5 v ( k + ) = f B s C A | N | 2 u ¯ ( k ) p γ 5 v ( k + ) ,
where the identity (A6) has already been taken into account. Moreover, the  | B s state normalization has been adjusted to the one that is conventionally used in the decay-constant definition (7), namely
B s ( q ) | B s ( p ) = 2 p 0 ( 2 π ) 3 δ ( 3 ) p q .
To verify that the global normalization in Equation (A10) is correct, one can begin with the relevant S matrix element
S   out μ + μ | B s ( H , L ) in = l + l | T [ exp ( i d 4 x L int ) ] | B s ( H , L ) .
In its evaluation, the  B s meson hadronic structure cannot be treated in a perturbative manner. Therefore, in the definition of the in/out states in the above equation, the interaction part of the Lagrangian L int that is switched off at timelike infinities consists only of the weak part N n C n Q n + h . c . . It is assumed that non-perturbative QCD effects have been solved beforehand and are contained in | B s ( H , L ) in that is treated as an asymptotic state from the weak interaction perspective.
In practice, the right normalization in Equation (A10) can be determined with possibly least effort by considering an analogy with a certain purely perturbative theory. Suppose the muons interact with a massive real pseudoscalar ϕ via a dimension-five operator L ˜ int = λ M ( α ϕ ) μ ¯ γ α γ 5 μ . Once the tree-level invariant matrix element M for the decay ϕ μ + μ is quickly determined, it should be expressed in terms of the free-field matrix element  e i p x 0 | λ M α ϕ ( x ) | ϕ ( p ) , with the free-particle state | ϕ ( p ) normalized as in Equation (A11). Finally, replacing
e i p x 0 | λ M α ϕ ( x ) | ϕ ( p ) by C A | N | 2 e i p x 0 | j A α ( x ) | B s ( p ) ,
(see Equation (A6)), one obtains Equation (A10) with the proper overall normalization.
The r.h.s. of Equation (A10) can be further simplified by noticing that
u ¯ p γ 5 v = u ¯ ( k + k + ) γ 5 v = u ¯ ( k γ 5 γ 5 k + ) v = 2 m μ u ¯ γ 5 v ,
where the identities u ¯ ( k ) k = m μ u ¯ ( k )   and  k + v ( k + ) = m μ v ( k + )   have been used in the last step. Let us note that our expression would vanish in the absence of γ 5 , which is related to the vanishing divergence of the vectorlike current, the same one that enters the QED interactions of muons and photons. This is precisely the reason why we were allowed to skip Q V in our initial operator list in Equation (5).
Let us now turn to the operators Q P and Q S . The pseudoscalar current j P = b ¯ γ 5 s that is present in both of them can be eliminated in favor of the axial current j A μ using Equations of Motion (EoMs) for the s- and b ¯ -quark operator fields:
[ i g A a T a e A m s ] s = EoM O ( M W 2 ) , b ¯ [ i + g A a T a + e A + m b ] = EoM O ( M W 2 ) ,
where g and e are the QCD and QED couplings, respectively, A μ a is the gluon field, A μ is the photon field, and  T a are the SU(3) generators. The  O ( M W 2 ) effects on the r.h.s. of the above equations stand for higher-order weak-interaction effects that stem from the operators Q n in Equation (). Such effects will be neglected below, as we are going to use the EoM to transform the weak operators Q P and Q S , while we work at the leading order in weak interactions. Operators that vanish by the EoM will commonly be denoted by E . They can be skipped in evaluation of observables at the leading order in weak interactions, as their physical matrix elements vanish.
To proceed, we multiply the first equation in (A14) by ( b ¯ γ 5 ) from the left, the second one by ( γ 5 s ) from the right, take their difference, rearrange, and obtain
α j A α = i ( m b + m s ) j P + E .
Using the above identity, one can express the following total derivative as
α [ ( l ¯ l ) j A α ] = j A α α [ l ¯ l ] + i ( m b + m s ) Q S + E .
Consequently,
Q S = ( i α [ l ¯ l ] ) j A α m b + m s + E + T ,
where T commonly denotes total derivatives of operators that are invariant under QCD and QED gauge transformations. Similarly,
Q P = ( i α [ l ¯ γ 5 l ] ) j A α m b + m s + E + T .
Since T in the Lagrangian has no effect on physical observables, we are allowed to replace Q S and Q P by the first (explicit) terms in the above equations. After such replacements, the same quark current shows up in Q A , Q S , and Q P , which means that only a single non-perturbative quantity, namely f B s , is sufficient to describe their physical matrix elements at the leading order in QED.
Using Equations (A17) and (A18), and performing calculations along the same way as in the case of Q A , one obtains
M P ( H ) = M B s 2 f B s m b + m s C P | N | 2 u ¯ ( k ) γ 5 v ( k + ) ,
and
M S ( L ) = M B s 2 f B s m b + m s C S | N | 2 u ¯ ( k ) v ( k + ) .
The total invariant matrix element for B s ( H ) reads
M ( H ) = M P ( H ) + M A ( H ) = M B s C P m b + m s 2 m μ C A M B s M B s f B s | N | 2 u ¯ ( k ) γ 5 v ( k + ) .
It remains to take the moduli squared and perform the sums over spins:
spins | u ¯ γ 5 v | 2 = Tr ( k + m μ ) γ 5 ( k + m μ ) γ 5 = 2 M B s 2 , spins | u ¯ v | 2 = Tr ( k + m μ ) ( k + m μ ) = 2 M B s 2 β 2 .
Finally, after substitution to Equation (A8) and then (A7), we end up with the branching ratio formula quoted in Equation (10).
In BSM models with generic (not SM-like) CP -violation, the Wilson coefficients C A , C P , and C S are not necessarily real. Moreover, the mass eigenstates are not necessarily given by Equation (A1), as the phase factors may differ from N / | N | and N / | N | . The latter effect can be described by introducing an extra complex phase ϕ s B S M ϕ s c c ¯ s arg [ ( V t s V t b ) 2 ] , see Section 2.2 of Ref. [15]. In such models, the branching ratio formula generalizes to
B ¯ s μ = | N | 2 M B s 3 f B s 2 8 π Γ H s β | r C A u C P | 2 F P + | u β C S | 2 F S + O ( α e ) ,
with r, β , and u that have been defined below Equations (6) and (10), as well as
F P 1 Γ L s Γ H s Γ L s sin 2 1 2 ϕ s BSM + arg ( r C A u C P ) , F S 1 Γ L s Γ H s Γ L s cos 2 1 2 ϕ s BSM + arg ( r C S ) .
The above expressions for F P and F S have been derived from the results of Refs. [40,41] where an interesting discussion concerning time-dependent observables can be found.

Appendix B. Numerical Update for B ¯ d μ SM

Here, we present the current SM prediction for B ¯ d μ , i.e., the average time-integrated branching ratio of B d μ + μ . The necessary formula is analogous to Equation (24), obtained by combining the semi-numerical expression from Equation (7) of Ref. [10] with η B B S (23) derived from Ref. [12]. It reads
B ¯ d μ S M × 10 10 = ( 1.06 ± 0.02 ) M t [ GeV ] 173.1 3.06 α s ( M Z ) 0.1184 0.18 R d η B B S ,
where
R d = f B d [ MeV ] 190.5 2 | V t b V t d | 0.0088 2 τ av d [ ps ] 1.519 .
As input, we need three more parameters in addition to those already listed in Table 1, namely the decay constant f B d of the B d meson, the average lifetime of this meson τ av d , and the relevant CKM factor | V t b V t d | . Their current values are listed in Table A1. Our use of the well-measured τ av d instead of τ H d is a good approximation thanks to the very small width difference predicted in the SM for the B d - B ¯ d system: Δ Γ d / ( 2 Γ av d ) = 0.00172 ( 46 ) [42].
Table A1. Numerical values of the extra input parameters that matter for B ¯ d μ .
Table A1. Numerical values of the extra input parameters that matter for B ¯ d μ .
ParameterValueUnitRefs.
f B d 190.0 (1.3)MeV[16,17,18,19,20]
| V t b V t d | 0.00851 (10)-derived from Ref. [36]
τ av d 1.517 (4)ps[37]
After substituting all the numerical inputs to Equation (A25), we find
B ¯ d μ S M = ( 9.71 ± 0.33 ) × 10 11 .
The uncertainty breakdown is presented in Table A2, in the same manner as it was done for B ¯ s μ S M in Table 2.
Table A2. The current uncertainty breakdown in B ¯ d μ S M , as compared to the 2013 one.
Table A2. The current uncertainty breakdown in B ¯ d μ S M , as compared to the 2013 one.
f B d CKM τ av d M t α s η BBS OtherNon-Parametric
2024 [this paper] 1.4 % 2.4 % 0.3 % 0.5 % 0.1 % 0.5 %<0.1% 1.5 % 3.4 %
2013 [10] 4.5 % 6.9 % 0.5 % 1.6 % 0.1 % 0.0 %<0.1% 1.5 % 8.5 %

References

  1. Aad, G.; Abajyan, T.; Abbott, B.; Abdallah, J.; Khalek, S.A.; Abdelalim, A.; Abdinov, O.; Aben, R.; Abi, B.; Abolins, M.; et al. Observation of a new particle in the search for the Standard Model Higgs boson with the ATLAS detector at the LHC. Phys. Lett. B 2012, 716, 16565–16590. [Google Scholar] [CrossRef]
  2. Chatrchyan, S.; Khachatryan, V.; Sirunyan, A.; Tumasyan, A.; Adam, W.; Aguilo, E.; Bergauer, T.; Dragicevic, M.; Erö, J.; Fabjan, C.; et al. Observation of a New Boson at a Mass of 125 GeV with the CMS Experiment at the LHC. Phys. Lett. B 2012, 716, 30–61. [Google Scholar] [CrossRef]
  3. Altmannshofer, W.; Stangl, P. New physics in rare B decays after Moriond 2021. Eur. Phys. J. C 2021, 81, 952. [Google Scholar] [CrossRef] [PubMed]
  4. Aaij, R.; Abellan Beteta, C.; Adametz, A.; Adeva, B.; Adinolfi, M.; Adrover, C.; Affolder, A.; Ajaltouni, Z.; Albrecht, J.; Alessio, F.; et al. First Evidence for the Decay B s 0 μ+μ. Phys. Rev. Lett. 2013, 110, 021801. [Google Scholar] [CrossRef] [PubMed]
  5. Sirunyan, A.M.; Tumasyan, A.; Adam, W.; Ambrogi, F.; Bergauer, T.; Brandstetter, J.; Dragicevic, M.; Erö, J.; Escalante Del Valle, A.; Flechl, M.; et al. Measurement of properties of B s 0 μ+μ decays and search for B0μ+μ with the CMS experiment. J. High Energy Phys. 2020, 2020, 188. [Google Scholar] [CrossRef]
  6. Tumasyan, A.; Adam, W.; Andrejkovic, J.W.; Bergauer, T.; Chatterjee, S.; Damanakis, K.; Dragicevic, M.; Del Valle, A.E.; Hussain, P.S.; Jeitler, M.; et al. Measurement of the B S 0 μ+μ decay properties and search for the B0μ+μ decay in proton-proton collisions at s = 13 TeV. Phys. Lett. B 2023, 842, 137955. [Google Scholar] [CrossRef]
  7. Aaij, R.; Beteta, C.A.; Ackernley, T.; Adeva, B.; Adinolfi, M.; Afsharnia, H.; Aidala, C.A.; Aiola, S.; Ajaltouni, Z.; Akar, S.; et al. Analysis of Neutral B-Meson Decays into Two Muons. Phys. Rev. Lett. 2022, 128, 041801. [Google Scholar] [CrossRef] [PubMed]
  8. Aaboud, M.; Aad, G.; Abbott, B.; Abbott, D.C.; Abdinov, O.; Abeloos, B.; Abhayasinghe, D.K.; Abidi, S.H.; AbouZeid, O.S.; Abraham, N.L.; et al. Study of the rare decays of B s 0 and B0 mesons into muon pairs using data collected during 2015 and 2016 with the ATLAS detector. J. High Energy Phys. 2019, 2019, 098. [Google Scholar] [CrossRef]
  9. Navas, S.; Amsler, C.; Gutsche, T.; Hanhart, C.; Hernández-Rey, J.J.; Lourenço, C.; Masoni, A.; Mikhasenko, M.; Mitchell, R.E.; Patrignani, C.; et al. Review of Particle Physics. Phys. Rev. D 2024, 110, 030001. Available online: https://pdg.lbl.gov (accessed on 1 June 2024).
  10. Bobeth, C.; Gorbahn, M.; Hermann, T.; Misiak, M.; Stamou, E.; Steinhauser, M. Bs,dl+l in the Standard Model with Reduced Theoretical Uncertainty. Phys. Rev. Lett. 2014, 112, 101801. [Google Scholar] [CrossRef]
  11. Beneke, M.; Bobeth, C.; Szafron, R. Enhanced electromagnetic correction to the rare B-meson decay Bs,dμ+μ. Phys. Rev. Lett. 2018, 120, 011801. [Google Scholar] [CrossRef]
  12. Beneke, M.; Bobeth, C.; Szafron, R. Power-enhanced leading-logarithmic QED corrections to Bqμ+μ. J. High Energy Phys. 2019, 10, 232, Erratum in J. High Energy Phys. 2022, 11, 099. [Google Scholar] [CrossRef]
  13. Arbey, A.; Mahmoudi, F.; Stal, O.; Stefaniak, T. Status of the Charged Higgs Boson in Two Higgs Doublet Models. Eur. Phys. J. C 2018, 78, 182. [Google Scholar] [CrossRef]
  14. Arbey, A.; Battaglia, M.; Mahmoudi, F. Constraints on the MSSM from the Higgs Sector: A pMSSM Study of Higgs Searches, B s 0 μ+μ and Dark Matter Direct Detection. Eur. Phys. J. C 2012, 72, 1906. [Google Scholar] [CrossRef]
  15. Buras, A.J.; Fleischer, R.; Girrbach, J.; Knegjens, R. Probing New Physics with the Bsμ+μ Time-Dependent Rate. J. High Energy Phys. 2013, 2013, 077. [Google Scholar] [CrossRef]
  16. Bazavov, A.; Bernard, C.; Brown, N.; DeTar, C.; El-Khadra, X.; Gámiz, E.; Gottlieb, S.; Heller, M.; Komijani, J.; Kronfeld, S.; et al. B- and D-meson leptonic decay constants from four-flavor lattice QCD. Phys. Rev. D 2018, 98, 074512. [Google Scholar] [CrossRef]
  17. Bussone, A.; Carrasco, N.; Dimopoulos, P.; Frezzotti, R.; Lami, P.; Lubicz, V.; Picca, E.; Riggio, L.; Rossi, G.C.; Simula, S.; et al. Mass of the b quark and B -meson decay constants from Nf=2+1+1 twisted-mass lattice QCD. Phys. Rev. D 2016, 93, 114505. [Google Scholar] [CrossRef]
  18. Dowdall, R.J.; Davies, C.T.H.; Horgan, R.R.; Monahan, C.J.; Shigemitsu, J. B-Meson Decay Constants from Improved Lattice Nonrelativistic QCD with Physical u, d, s, and c Quarks. Phys. Rev. Lett. 2013, 110, 222003. [Google Scholar] [CrossRef]
  19. Hughes, C.; Davies, C.T.H.; Monahan, C.J. New methods for B meson decay constants and form factors from lattice NRQCD. Phys. Rev. D 2018, 97, 054509. [Google Scholar] [CrossRef]
  20. Aoki, Y.; Blum, T.; Colangelo, G.; Collins, S.; Della Morte, M.; Dimopoulos, P.; Dürr, S.; Feng, X.; Fukaya, H.; Golterman, M.; et al. FLAG Review 2021. Eur. Phys. J. C 2022, 82, 869. [Google Scholar] [CrossRef]
  21. Logan, H.E.; Nierste, U. Bs,d+ in a two Higgs doublet model. Nucl. Phys. B 2000, 586, 39–55. [Google Scholar] [CrossRef]
  22. Misiak, M.; Urban, J. QCD corrections to FCNC decays mediated by Z-penguins and W-boxes. Phys. Lett. B 1999, 451, 161–169. [Google Scholar] [CrossRef]
  23. Dugan, M.J.; Grinstein, B. On the vanishing of evanescent operators. Phys. Lett. B 1991, 256, 239–244. [Google Scholar] [CrossRef]
  24. Hermann, T.; Misiak, M.; Steinhauser, M. Three-loop QCD corrections to Bsμ+μ. J. High Energy Phys. 2013, 12, 097. [Google Scholar] [CrossRef]
  25. Steinhauser, M. Results and techniques of multiloop calculations. Phys. Rept. 2002, 364, 247–357. [Google Scholar] [CrossRef]
  26. Misiak, M.; Steinhauser, M. Three loop matching of the dipole operators for b and bsg. Nucl. Phys. B 2004, 683, 277–305. [Google Scholar] [CrossRef]
  27. Inami, T.; Lim, C.S. Effects of Superheavy Quarks and Leptons in Low-Energy Weak Processes KLμ μ ¯ , K+π+ν ν ¯ and K0 K ¯ 0. Prog. Theor. Phys. 1981, 65, 297–314. [Google Scholar] [CrossRef]
  28. Buchalla, G.; Buras, A.J. QCD corrections to rare K- and B-decays for arbitrary top quark mass. Nucl. Phys. B 1993, 400, 225–239. [Google Scholar] [CrossRef]
  29. Buchalla, G.; Buras, A.J. QCD corrections to the sdZ vertex for arbitrary top quark mass. Nucl. Phys. B 1993, 398, 285–300. [Google Scholar] [CrossRef]
  30. Bobeth, C.; Gorbahn, M.; Stamou, E. Electroweak Corrections to Bs,d+. Phys. Rev. D 2014, 89, 034023. [Google Scholar] [CrossRef]
  31. Huber, T.; Lunghi, E.; Misiak, M.; Wyler, D. Electromagnetic logarithms in B ¯ Xsl+l. Nucl. Phys. B 2006, 740, 105–137. [Google Scholar] [CrossRef]
  32. Beneke, M.; Buchalla, G.; Neubert, M.; Sachrajda, C.T. QCD factorization for Bππ decays: Strong phases and CP violation in the heavy quark limit. Phys. Rev. Lett. 1999, 83, 1914–1917. [Google Scholar] [CrossRef]
  33. Beneke, M.; Rohrwild, J. B meson distribution amplitude from Bγℓ ν ¯ . Eur. Phys. J. C 2011, 71, 1818. [Google Scholar] [CrossRef]
  34. Feldmann, T.; Gubernari, N.; Huber, T.; Seitz, N. Contribution of the electromagnetic dipole operator O7 to the B ¯ sμ+μ decay amplitude. Phys. Rev. D 2023, 107, 013007. [Google Scholar] [CrossRef]
  35. Finauri, G.; Gambino, P. The q2 moments in inclusive semileptonic B decays. J. High Energy Phys. 2024, 2024, 206. [Google Scholar] [CrossRef]
  36. Charles, J.; Hocker, A.; Lacker, H.; Laplace, S.; Le Diberder, F.R.; Malcles, J.; Ocariz, J.; Pivk, M.; Roos, L. CP violation and the CKM matrix: Assessing the impact of the asymmetric B factories. Eur. Phys. J. C 2005, 41, 1–131. [Google Scholar] [CrossRef]
  37. Amhis, Y.S.; Banerjee, S.; Ben-Haim, E.; Bertholet, E.; Bernlochner, F.U.; Bona, M.; Bozek, A.; Bozzi, C.; Brodzicka, J.; Chobanova, V.; et al. Averages of b-hadron, c-hadron, and τ-lepton properties as of 2021. Phys. Rev. D 2023, 107, 052008. [Google Scholar] [CrossRef]
  38. Peskin, M.E.; Schroeder, D.V. An Introduction to Quantum Field Theory; Addison-Wesley: Reading, PA, USA, 1995. [Google Scholar] [CrossRef]
  39. Lang, M. Leptonic Decays of Neutral B Mesons in the Three-Spurion Two-Higgs-Doublet Model. Ph.D. Thesis, KIT, Karlsruhe, Germany, 2023. [Google Scholar] [CrossRef]
  40. De Bruyn, K.; Fleischer, R.; Knegjens, R.; Koppenburg, P.; Merk, M.; Pellegrino, A.; Tuning, N. Probing New Physics via the B s 0 μ+μ Effective Lifetime. Phys. Rev. Lett. 2012, 109, 041801. [Google Scholar] [CrossRef]
  41. De Bruyn, K.; Fleischer, R.; Knegjens, R.; Koppenburg, P.; Merk, M.; Tuning, N. Branching Ratio Measurements of Bs Decays. Phys. Rev. D 2012, 86, 014027. [Google Scholar] [CrossRef]
  42. Asatrian, H.M.; Asatryan, H.H.; Hovhannisyan, A.; Nierste, U.; Tumasyan, S.; Yeghiazaryan, A. Penguin contribution to the width difference and CP asymmetry in Bq- B ¯ q mixing at order α s 2 Nf. Phys. Rev. D 2020, 102, 033007. [Google Scholar] [CrossRef]
Figure 1. LHCb measurements of B ¯ s μ and B ¯ d μ at 4.4 and 9 fb 1 integrated luminosity, compared to the current SM predictions. Plot from Figure 2 of Ref. [7].
Figure 1. LHCb measurements of B ¯ s μ and B ¯ d μ at 4.4 and 9 fb 1 integrated luminosity, compared to the current SM predictions. Plot from Figure 2 of Ref. [7].
Symmetry 16 00917 g001
Figure 2. Examples of the Feynman diagrams appearing in the calculation of Z i j . The double square in the middle denotes an insertion of either Q A or Q A E . Diagrams from Figure 2 of Ref. [24].
Figure 2. Examples of the Feynman diagrams appearing in the calculation of Z i j . The double square in the middle denotes an insertion of either Q A or Q A E . Diagrams from Figure 2 of Ref. [24].
Symmetry 16 00917 g002
Figure 3. Examples of diagrams entering the SM side of the matching Equation (12) at various orders in QCD. The W-boxes and Z-penguins are shown in the top and bottom rows, respectively. Contributions of the orders α s 0 , α s 1 , and α s 2 are arranged from left to right. Diagrams from Figures 1 and 4 of Ref. [24].
Figure 3. Examples of diagrams entering the SM side of the matching Equation (12) at various orders in QCD. The W-boxes and Z-penguins are shown in the top and bottom rows, respectively. Contributions of the orders α s 0 , α s 1 , and α s 2 are arranged from left to right. Diagrams from Figures 1 and 4 of Ref. [24].
Symmetry 16 00917 g003
Figure 4. Examples of diagrams contributing to the O ( α e ) EW correction on the SM side of the matching Equation (12). Diagrams from Figure 1 of Ref. [30].
Figure 4. Examples of diagrams contributing to the O ( α e ) EW correction on the SM side of the matching Equation (12). Diagrams from Figure 1 of Ref. [30].
Symmetry 16 00917 g004
Figure 5. Examples of diagrams with virtual photon exchanges that are responsible for the power-enhanced QED correction in Equation (23). Diagrams from Figure 1 of Ref. [11].
Figure 5. Examples of diagrams with virtual photon exchanges that are responsible for the power-enhanced QED correction in Equation (23). Diagrams from Figure 1 of Ref. [11].
Symmetry 16 00917 g005
Table 1. Numerical values of the updated input parameters.
Table 1. Numerical values of the updated input parameters.
ParameterValueUnitRefs.
f B s 230.3 (1.3)MeV[16,17,18,19,20]
| V c b |   ×   10 3 41.97 (48)-[35]
| V t b V t s / V c b | 0.9820 (4)-derived from Ref. [36]
τ H s 1.622 (8)ps[37]
M t 172.57 (29)GeV[9]
α s ( M Z ) 0.1180 (9)-[9]
Table 2. The current uncertainty breakdown in B ¯ s μ S M , as compared to the 2013 one.
Table 2. The current uncertainty breakdown in B ¯ s μ S M , as compared to the 2013 one.
f B s CKM τ H s M t α s η BBS OtherNon-Parametric
2024 [this paper] 1.1 % 2.3 % 0.5 % 0.5 % 0.1 % 0.5 %<0.1% 1.5 % 3.2 %
2013 [10] 4.0 % 4.3 % 1.3 % 1.6 % 0.1 % 0.0 %<0.1% 1.5 % 6.4 %
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Czaja, M.; Misiak, M. Current Status of the Standard Model Prediction for the Bsμ+μ Branching Ratio. Symmetry 2024, 16, 917. https://doi.org/10.3390/sym16070917

AMA Style

Czaja M, Misiak M. Current Status of the Standard Model Prediction for the Bsμ+μ Branching Ratio. Symmetry. 2024; 16(7):917. https://doi.org/10.3390/sym16070917

Chicago/Turabian Style

Czaja, Mateusz, and Mikołaj Misiak. 2024. "Current Status of the Standard Model Prediction for the Bsμ+μ Branching Ratio" Symmetry 16, no. 7: 917. https://doi.org/10.3390/sym16070917

APA Style

Czaja, M., & Misiak, M. (2024). Current Status of the Standard Model Prediction for the Bsμ+μ Branching Ratio. Symmetry, 16(7), 917. https://doi.org/10.3390/sym16070917

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop