REVIEW 5 minor 9 cited by
Test of lepton flavor universality and search for lepton flavor violation in $B \to K \ell\ell$ decays
T0 review · 0 major / 5 minor · reviewed 2026-08-14 · deepseek-v4-flash
Pith's one-line read This paper finds all measured $R_K$ values consistent with the Standard Model, with $R_K = 1.03^{+0.28}_{-0.24} \pm 0.01$ in the $1.0 < q^2 < 6.0$ GeV$^2$/c$^4$ bin, and bounds lepton-flavor-violating $B \to K\mu^{\pm}e^{\mp}$ branching…
desk verdict Belle's full-data R_K measurement is a careful, null-result paper that anchors the LFU picture; the analysis is solid and the few soft spots are not load-bearing. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The machinery is a three-dimensional unbinned extended maximum-likelihood fit, a simultaneous fit that extracts signal and background yields from three discriminating variables: the beam-energy-constrained mass $M_{bc}$, the energy difference $\Delta E$, and a transformed topology/output variable $O'$. Signal probability density functions come from Monte Carlo simulation and are calibrated on high-statistics $B \to J/\psi K$ control samples, which also validate the method by recovering $R_K(J/\psi)$ consistent with unity. $R_K$ is formed from efficiency-corrected signal yields, so most systematic uncertainties cancel in the double ratio; $A_I$ uses the same yields together with the $B^+$/$B^0$ lifetime ratio and production fractions, and the lepton-flavor-violating limits come from a frequentist pseudoexperiment procedure on the fitted yields.
What would settle it
Compare the $q^2 \in (1.0, 6.0)$ GeV$^2$/c$^4$ $R_K$ value and the 2.6$\sigma$ negative $A_I$ obtained here with the same measurement on an independent, larger data set; if both return to unity and zero, the deviations were statistical. More directly, measure the $K_S^0$ reconstruction efficiency in data versus simulation as a function of $K_S^0$ momentum from the $B \to J/\psi K_S^0$ control sample: a momentum-dependent bias of about 5-10% at low momentum would be enough to reproduce the negative muon $A_I$, and would invalidate the neutral-mode branching fractions.
Extended reading notes
Core claim
The measurement establishes that, at the sensitivity of a 711 fb$^{-1}$ sample, the ratio of $B \to K\mu^{+}\mu^{-}$ to $B \to Ke^{+}e^{-}$ branching fractions is consistent with the Standard Model prediction of lepton flavor universality in every $q^2$ bin. The value most comparable to earlier hints is $R_K = 1.03^{+0.28}_{-0.24} \pm 0.01$ for $1.0 < q^2 < 6.0$ GeV$^2$/c$^4$, and the whole-range value is $1.10^{+0.16}_{-0.15} \pm 0.02$. The same fits yield CP-averaged isospin asymmetries that are null within uncertainties, the single largest deviation being a 2.6$\sigma$ negative asymmetry in the muon mode in the (1.0, 6.0) bin. Differential branching fractions for charged $B$ decays agree with theoretical predictions, while neutral $B$ modes sit below them. In the lepton-flavor-violating search, no signal is established; 90% confidence-level upper limits are set in the $10^{-8}$ range, improving the neutral-mode limit by a factor of 7.1.
Load-bearing premise
The measurement assumes that the Monte Carlo signal shapes, calibrated with $B \to J/\psi K$ control samples, and the efficiency corrections describe the real $B \to K\ell\ell$ events equally well in every $q^2$ bin, especially in the low-yield neutral-$B$ modes where a small shape or efficiency bias would shift $R_K$ and the isospin asymmetry directly.
Editorial extensions
If this is right
- The $R_K$ value in the (1.0, 6.0) GeV$^2$/c$^4$ bin is above the previously reported low value by 1.6$\sigma$; if it holds, new-physics models predicting a large muon deficit in this bin are disfavored.
- The 2.6$\sigma$ negative isospin asymmetry for muons in the same bin, together with neutral-$B$ differential branching fractions below theory, points to a possible isospin-dependent effect that needs confirmation with more data.
- The 90% confidence-level upper limits of $8.5 \times 10^{-8}$, $3.0 \times 10^{-8}$ and $3.8 \times 10^{-8}$ for $B^+ \to K^+\mu^{+}e^{-}$, $B^+ \to K^+\mu^{-}e^{+}$ and $B^0 \to K^0\mu^{\pm}e^{\mp}$, respectively, improve the neutral-mode constraint by a factor of 7.1 and tighten the room for lepton-flavor-violating new physics.
- The control-sample result $R_K(J/\psi) = 0.994 \pm 0.011 \pm 0.010$ shows the analysis procedure recovers lepton flavor universality in a channel where the Standard Model is certain, supporting the reliability of the $R_K$ result.
Reading between the lines
- An implicit consequence is that next-generation B-factory data, with several ab$^{-1}$, will decide whether the (1.0, 6.0) $R_K$ central value drifts toward or away from unity, since the uncertainty here is still statistics-dominated.
- A testable explanation of the 2.6$\sigma$ negative $A_I$ is that it traces to a momentum-dependent $K_S^0$ efficiency rather than new physics; comparing the neutral-to-charged $B \to J/\psi K$ yield ratio in fine momentum bins would expose such a bias.
- The lepton-flavor-violation limits reached here imply that other $b \to s\ell\ell'$ modes, such as $B \to K^*\mu e$, should be within reach of next-generation B-factory data if any lepton-flavor violation accompanies lepton-flavor-universality violation.
- If the neutral-$B$ differential branching fractions remain below theory at higher statistics, the more likely explanation would be an isospin-breaking hadronic effect or a $K_S^0$ reconstruction normalization issue rather than new physics, given that the charged modes agree with predictions.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The Belle Collaboration presents measurements of the branching fractions of $B \to K \mu^+\mu^-$ and $B \to K e^+e^-$, the ratio $R_K$, the $CP$-averaged isospin asymmetry $A_I$, and a search for lepton-flavor-violating $B \to K \mu^\pm e^\mp$ decays, using 711 fb$^{-1}$ of $\Upsilon(4S)$ data containing $772 \times 10^6$ $B\bar{B}$ events. The analysis employs a three-dimensional unbinned extended maximum-likelihood fit to $M_{\rm bc}$, $\Delta E$, and $O'$ in five $q^2$ bins and the full region. The key results are $R_K = 1.03^{+0.28}_{-0.24} \pm 0.01$ for $q^2 \in (1.0, 6.0)$ GeV$^2/c^4$, with all $R_K$ values consistent with the Standard Model; $A_I$ values consistent with zero, with the largest deviation being $2.6\sigma$ in the muon mode in the $(1.0, 6.0)$ bin; and 90% CL upper limits on the LFV branching fractions of $8.5 \times 10^{-8}$, $3.0 \times 10^{-8}$, and $3.8 \times 10^{-8}$ for $B^+ \to K^+ \mu^+ e^-$, $B^+ \to K^+ \mu^- e^+$, and $B^0 \to K^0 \mu^\pm e^\mp$, respectively. The measurement is blinded, signal shapes are calibrated with $B \to J/\psi K$ control samples, and systematic uncertainties are evaluated in detail.
Significance. The measurement provides the most precise Belle determination of $R_K$ and the first Belle result in the $(1.0, 6.0)$ GeV$^2/c^4$ bin, where LHCb reported a $2.5\sigma$ deviation from the Standard Model. A result consistent with the SM is important for the global picture of $b \to s \ell\ell$ anomalies. The internal consistency checks, including $R_K(J/\psi) = 0.994 \pm 0.011 \pm 0.010$ and $A_I(B \to J/\psi K) = -0.002 \pm 0.006 \pm 0.014$, support the reliability of the shape and efficiency transfer. The LFV limits improve the neutral-mode limit by a factor of 7.1 over the previous BaBar result. The paper is carefully worded, with statistical and systematic uncertainties appropriately separated, and the $2.6\sigma$ isospin asymmetry is correctly framed as a consistency check rather than as evidence for new physics. If the central claims hold, the results constitute a meaningful constraint on lepton flavor universality and on new-physics scenarios that violate lepton flavor.
minor comments (5)
- [Section 4] The sentence "The results are listed in Table 1 andRK andAI are also shown in Figs. 3 and 4, respectively" is missing spaces after "and" in the preprint text; please correct the typesetting in the final version.
- [References] References [13] and [41] are identical (both cite R. Aaij et al., JHEP 06, 133 (2014)); they should be consolidated or given distinct contextual identifiers to avoid duplication.
- [Table 1] Table 1 is dense and the formatting of the $A_I$ and $R_K$ columns makes it difficult to distinguish per-mode values from combined values; consider reformatting so that each row is unambiguous.
- [Section 5] The text states which systematic uncertainties cancel in $R_K$ and $A_I$ but does not explicitly say that the lepton-identification uncertainties for the signal $R_K$ measurement are included in the quoted total systematics; please state this explicitly, since otherwise a reader may infer that the $\pm 0.97\%$ lepton-ID uncertainty in Table 4 is absent from the $R_K$ results.
- [Section 4] The fit validation for the low-yield neutral modes, such as $B^0 \to K_S^0 \mu^+\mu^-$ with approximately 3.9 signal events in the $(1.0, 6.0)$ bin, is not described; a brief report of pseudoexperiment checks on fit bias and pull distributions would increase confidence in the quoted statistical uncertainties and in the $2.6\sigma$ isospin asymmetry.
Circularity Check
No significant circularity: the central R_K, A_I, and LFV results are measured yields divided by efficiencies and compared with external SM predictions, with no fitted input renamed as a prediction.
full rationale
The paper's derivation chain is self-contained and data-driven. R_K is defined by Eq. (1.1) as a ratio of measured partial widths, and the quoted values in Table 1 are obtained from three-dimensional unbinned maximum-likelihood fits to Mbc, Delta E, and O' with signal yields floated in each q^2 bin. The Standard Model comparison is external: differential branching fractions are compared with theoretical predictions from light-cone sum rule and lattice QCD calculations in Refs. [38,39], and no SM value is used to define or constrain the fit. The J/psi K control samples calibrate signal PDF shapes and the O>O_min efficiency and provide cross-checks such as R_K(J/psi)=0.994 +/- 0.011 +/- 0.010 and A_I(B->J/psi K)=-0.002 +/- 0.006 +/- 0.014, but these calibrations are not the target observables and do not force R_K or A_I toward any particular value. The LFV upper limits are obtained from fitted yields converted with efficiencies and the number of B-Bbar pairs, with no assumed signal branching fraction. The only self-citations (Refs. [5,7]) are prior Belle measurements that this analysis supersedes; they are not used as load-bearing justification. I find no step in which a prediction is equivalent to an input by construction, no fitted parameter renamed as a prediction, and no self-citation chain that imports the conclusion.
Assumptions & free parameters
assumptions (4)
- domain assumption Monte Carlo simulation of signal and background shapes and efficiencies is accurate after control-sample calibration.
- domain assumption External Standard Model predictions for differential branching fractions and R_K in references 38 and 39 are correct within their quoted uncertainties.
- domain assumption Production fractions f_plus-minus/f_00 and the lifetime ratio tau_B+/tau_B0 from PDG and HFLAV are accurate.
- domain assumption The ARGUS, Gaussian, and Crystal Ball parameterizations adequately describe the fit variables M_bc, Delta E, and O'.
Cite this review
Pith. "Pith review of Test of lepton flavor universality and search for lepton flavor violation in $B \to K \ell\ell$ decays." pith.science (2026). https://pith.science/paper/7UTHT23N
@misc{pith2026190801848,
author = {Pith},
title = {Pith review of: Test of lepton flavor universality and search for lepton flavor violation in $B \to K \ell\ell$ decays},
year = {2026},
howpublished = {\url{https://pith.science/paper/7UTHT23N}},
note = {Machine review of arXiv:1908.01848}
}
abstract
We present measurements of the branching fractions for the decays $B\to K \mu^{+}\mu^{-}$ and $B\to K e^{+}e^{-}$, and their ratio ($R_{K}$), using a data sample of 711 $fb^{-1}$ that contains $772 \times 10^{6}$ $B\bar{B}$ events. The data were collected at the $\Upsilon(4S)$ resonance with the Belle detector at the KEKB asymmetric-energy $e^{+}e^{-}$ collider. The ratio $R_{K}$ is measured in five bins of dilepton invariant-mass-squared ($q^{2}$): $q^{2} \in (0.1, 4.0), (4.0, 8.12), (1.0, 6.0)$, $(10.2, 12.8)$ and ($>14.18) GeV^{2}/c^{4}$, along with the whole $q^2$ region. The $R_{K}$ value for $q^{2} \in (1.0, 6.0) GeV^{2}/c^{4}$ is $1.03^{+0.28}_{-0.24} \pm 0.01$. The first and second uncertainties listed are statistical and systematic, respectively. All results for $R_{K}$ are consistent with Standard Model predictions. We also measure $C\!P$-averaged isospin asymmetries in the same $q^{2}$ bins. The results are consistent with a null asymmetry, with the largest difference of 2.6 standard deviations occurring for the $q^{2}\in(1.0,6.0) GeV^{2}/c^{4}$ bin in the mode with muon final states. The measured differential branching fractions, ${d\cal B}/{dq^{2}}$, are consistent with theoretical predictions for charged $B$ decays, while the corresponding values are below the expectations for neutral $B$ decays. We have also searched for lepton-flavor-violating $B \rightarrow K\mu^{\pm}e^{\mp}$ decays and set $90\%$ confidence-level upper limits on the branching fraction in the range of $10^{-8}$ for $B^{+} \rightarrow K^{+}\mu^{\pm}e^{\mp}$, and $B^{0} \rightarrow K^{0}\mu^{\pm}e^{\mp}$ modes.
Forward citations
Cited by 9 Pith papers
-
Gauging the accidental symmetries of the Standard Model, and implications for the flavour anomalies
Gauging an almost arbitrary combination of baryon number, lepton numbers and hypercharge, with three new singlet fermions, yields a Z' framework that can explain the B decay flavour anomalies with a universal axial co...
-
Long distance effects in inclusive rare $B$ decays and phenomenology of $\bar{B}\to X_d \ell^+\ell^-$
A new Standard Model prediction for the inclusive \bar B -> X_d \ell^+\ell^- decay that incorporates charmonium resonances, resolved power corrections, QED logarithms, and five-body contributions.
-
Exhaustive Model Selection in $b \to s \ell \ell$ Decays: Pitting Cross-Validation against AIC$_c$
An exhaustive model selection over 511 Wilson-coefficient combinations for b to s l l decays finds that all surviving scenarios contain O9 (Delta C9 near -1.1 to -1.4), the only surviving one-operator scenario.
-
The effects of a scalar singlet Leptoquark at the $Z$ factory
A scalar singlet leptoquark that explains B-meson anomalies produces a ~0.7% decrease in Z→τ+τ−, which future Z-factory measurements could detect, while Z→μ+μ− is essentially unchanged.
-
Resolving the $(g-2)_{\mu}$ and $B$ anomalies with leptoquarks and a dark Higgs boson
A model with a light dark Higgs and TeV-scale leptoquarks can simultaneously explain the muon g-2 and B decay anomalies, and predicts rare B, K, and Higgs decays within reach of current experiments.
-
Features of Charged Lepton Flavor Violation in an $A_4$ Symmetric Neutrino Mass Model
A reanalysis of an A4 neutrino mass model concludes meson CLFV decays can be near current bounds while radiative CLFV vanishes, but the radiative amplitude is computed incorrectly.
-
Constraints on Dark Photon and Dark $Z$ Model Parameters in the $B$ and $K$ Meson Decays
Only a fine-tuned dark Z model with cancelled electron couplings survives the combined constraints, but its muon g-2 contribution is orders of magnitude too large.
-
The rare decay $B^+ \to K^+\ell^+\ell^-(\nu\bar{\nu})$ under the QCD sum rules approach
A QCD light-cone sum-rule calculation with the authors' kaon distribution amplitudes predicts B(B+ to K+ nu nubar) = 4.14 x 10^-6 and B(B+ to K+ l+l-) around 6.6 x 10^-7, consistent with other SM estimates.
-
Fitting in or odd one out? Pulls vs residual responses in $b\to s \ell^+\ell^-$
The updated LHCb R_K measurement shifts the best-fit new physics in b→sμμ transitions toward smaller C10mu and reduces its uncertainty.
Reference graph
Works this paper leans on
- [1]
- [2]
- [3]
- [4]
- [5]
- [6]
-
[7]
J. T. Wei et al. (Belle Collaboration), Phys. Rev. Lett. 103, 171801 (2009)
work page 2009
- [8]
Show all 41 references
-
[9]
Zyla et al
P. Zyla et al. (Particle Data Group), Prog. Theor. Exp. Phys. 2020, 083C01 (2020)
2020
-
[10]
Amhis, et al., Online update at http://www.slac.sta nford.edu/xorg/hfag, arXiv:1909.12524
Y. Amhis, et al., Online update at http://www.slac.sta nford.edu/xorg/hfag, arXiv:1909.12524
1909 arXiv
-
[11]
Lyon and R
J. Lyon and R. Zwicky, Phys. Rev. D 88, 094004 (2013)
2013
-
[12]
J. P. Lees et al. (BaBar Collaboration), Phys. Rev. D 86, 032012 (2012)
2012
-
[13]
Aaij et al
R. Aaij et al. (LHCb Collaboration), J. High Energy Phys. 06, 133 (2014)
2014
-
[14]
S. L. Glashow, D. Guadagnoli, and K. Lane, Phys. Rev. Let t. 114, 091801 (2015)
2015
-
[15]
J. C. Helo, S. Kovalenko, and I. Schmidt, Nucl. Phys. B853, 80 (2011)
2011
-
[16]
Aaij et al
R. Aaij et al. (LHCb Collaboration), Phys. Rev. Lett 123, 241802 (2019)
2019
-
[17]
Aubert et al
B. Aubert et al. (BaBar Collaboration), Phys. Rev. D 73, 092001 (2006)
2006
-
[18]
Kurokawa and E
S. Kurokawa and E. Kikutani, Nucl. Instrum. Methods Phy s. Res., Sec. A 499, 1 (2003), and other papers included in this volume; T. Abe et al. , Prog. Theor. Exp. Phys. 2013, 03A001 (2013) and following articles up to 03A011
2003
-
[19]
Abashian et al
A. Abashian et al. (Belle Collaboration), Nucl. Instrum. Methods Phys. Res., Sec. A 479, 117 (2002); also, see the detector section in J. Brodzicka et al. , Prog. Theor. Exp. Phys. 2012, 04D001 (2012)
2002
-
[20]
Natkaniec et al
Z. Natkaniec et al. (Belle SVD2 Group), Nucl. Instrum. Methods Phys. Res., Sec. A 560, 1 (2006)
2006
-
[21]
D. J. Lange, Nucl. Instrum. Methods Phys. Res., Sec. A 462, 152 (2001)
2001
-
[22]
A. Ali, P. Ball, L. T. Handoko, and G. Hiller, Phys. Rev. D 61, 074024 (2000)
2000
-
[23]
Barberio and Z
E. Barberio and Z. W¸ as, Comput. Phys. Commun. 79, 291 (1994); P. Golonka and Z. W¸ as, Eur. Phys. J. C 45, 97 (2006); P. Golonka and Z. W¸ as, Eur. Phys. J. C 50, 53 (2007)
1994
-
[24]
Brun et al
R. Brun et al. , CERN Report No. DD/EE/84-1 (1984)
1984
-
[25]
The inclusion of the charge-conjugate decay mode is imp lied unless otherwise stated
-
[26]
Abashian et al
A. Abashian et al. , Nucl. Instrum. Methods Phys. Res., Sec. A 491, 69 (2002). – 17 –
2002
-
[27]
Nakano, Nucl
E. Nakano, Nucl. Instrum. Methods Phys. Res., Sec. A 494, 402 (2002)
2002
-
[28]
Hanagaki, H
K. Hanagaki, H. Kakuno, H. Ikeda, T. Iijima, and T. Tsuka moto, Nucl. Instrum. Methods Phys. Res., Sec. A 485, 490 (2002)
2002
-
[29]
S. H. Lee et al. (Belle Collaboration), Phys. Rev. Lett. 91, 261801 (2003)
2003
-
[30]
G. C. Fox and S. Wolfram, Phys. Rev. Lett. 41, 1581 (1978)
1978
-
[31]
Kakuno et al
H. Kakuno et al. (Belle Collaboration), Nucl. Instrum. Methods Phys. Res., Sec. A 533, 516 (2004)
2004
-
[32]
D. M. Asner et al. (CLEO Collaboration), Phys. Rev. D 53, 1039 (1996)
1996
-
[33]
Skwarnicki, Ph.D
T. Skwarnicki, Ph.D. thesis, Institute for Nuclear Phy sics, Krakow; DESY Internal Report No. DESY F31-86-02, 1986
1986
-
[34]
Albrecht et al
H. Albrecht et al. (ARGUS Collaboration), Phys. Lett. B 241, 278 (1990)
1990
-
[35]
Dash et al
N. Dash et al. (Belle Collaboration), Phys. Rev. Lett. 119, 171801 (2017)
2017
-
[36]
Melikhov et al
D. Melikhov et al. , Phys. Lett. B 410, 290 (1997)
1997
-
[37]
Colangelo et al
P. Colangelo et al. , Phys. Rev. D 53, 3672 (1996)
1996
-
[38]
Bobeth, G
C. Bobeth, G. Hiller, and D. van Dyk, J. High Energy Phys. 07, 067 (2011)
2011
-
[39]
Bobeth, G
C. Bobeth, G. Hiller, D. van Dyk, and C. Wacker, J. High En ergy Phys. 01, 107 (2012)
2012
-
[40]
Aaij et al
R. Aaij et al. (LHCb Collaboration), Phys. Rev. Lett. 113, 151601 (2014)
2014
-
[41]
Aaij et al
R. Aaij et al. (LHCb Collaboration), J. High Energy Phys. 06, 133 (2014). – 18 –
2014
Reviewed August 14, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.