Search for the lepton-flavor-violating τ⁻ rightarrow e^(mp) ell^(pm) ell^(mp) decays at Belle II
Pith reviewed 2026-05-19 02:52 UTC · model grok-4.3
The pith
Belle II sets upper limits between 1.3 and 2.5 × 10^{-8} on branching fractions for lepton-flavor-violating tau decays at 90% confidence level.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
We present the result of a search for the charged-lepton-flavor violating decays τ⁻ → e∓ ℓ± ℓ∓, where ℓ is a muon or an electron, using a data sample with an integrated luminosity of 428 fb^{-1} recorded by the Belle II experiment at the SuperKEKB e⁺e⁻ collider. The selection of e⁺e⁻ → τ⁺τ⁻ events containing a signal candidate is based on an inclusive-tagging reconstruction and on a boosted decision tree to suppress background. Upper limits on the branching fractions between 1.3 and 2.5 × 10^{-8} are set at the 90% confidence level. These results are the most stringent bounds to date for four of the modes.
What carries the argument
Inclusive-tagging reconstruction of tau-pair events combined with a boosted decision tree that suppresses background while preserving signal efficiency.
If this is right
- The derived upper limits constitute the most stringent constraints available for four of the six decay modes examined.
- The results improve upon all previously published bounds for those channels.
- The analysis establishes a new reference sensitivity for future searches of the same decays with larger Belle II data sets.
Where Pith is reading between the lines
- These experimental bounds can be translated into limits on the strength of hypothetical new interactions that mediate the decays.
- The same reconstruction and classification strategy can be adapted to additional rare tau decay modes not included in the present search.
- Accumulation of more integrated luminosity at SuperKEKB will directly tighten the same limits without requiring changes to the core method.
Load-bearing premise
The boosted decision tree and inclusive-tagging selection suppress background without significant bias or loss of signal efficiency, relying on accurate simulation of detector response and background processes.
What would settle it
Observation of one or more events inside the signal region after all selection cuts that cannot be accounted for by the estimated background would indicate a branching fraction above the reported upper limit.
read the original abstract
We present the result of a search for the charged-lepton-flavor violating decays $\tau^- \rightarrow e^\mp \ell^\pm \ell^-$, where $\ell$ is a muon or an electron, using a data sample with an integrated luminosity of 428 fb$^{-1}$ recorded by the Belle II experiment at the SuperKEKB $e^+e^-$ collider. The selection of $e^+e^- \to\tau^+\tau^-$ events containing a signal candidate is based on an inclusive-tagging reconstruction and on a boosted decision tree to suppress background. Upper limits on the branching fractions between 1.3 and 2.5 $\times 10^{-8}$ are set at the 90% confidence level. These results are the most stringent bounds to date for four of the modes.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript reports a search for the six lepton-flavor-violating tau decays of the form τ⁻ → e∓ ℓ± ℓ∓ (ℓ = e or μ) in 428 fb⁻¹ of Belle II data. Events are selected via inclusive tagging of the opposite-sign tau and a boosted decision tree trained to suppress standard-model backgrounds; 90% CL upper limits on the branching fractions are extracted in the range 1.3–2.5 × 10^{-8}, stated to be the most stringent limits for four of the modes.
Significance. If the central results hold, the work improves the existing experimental bounds on charged-lepton-flavor violation in the tau sector by roughly a factor of two for several channels, using a large data set and standard multivariate techniques. The analysis follows established high-energy-physics practices for limit setting and background modeling; the explicit statement that four modes now have the world’s tightest constraints is a clear, falsifiable claim.
major comments (1)
- [Event selection and BDT description] The quoted upper limits depend on the product of signal efficiency and background normalization after the BDT cut and inclusive tagging. Both quantities are taken from Monte Carlo; the manuscript does not describe an independent data-driven validation of the BDT response or PID performance in a signal-like kinematic region. Any systematic mismatch in low-momentum lepton modeling would rescale the expected background or efficiency and directly affect the reported limits at the 10–20 % level typical for these analyses.
minor comments (2)
- List the six explicit final states (e.g., τ⁻ → e⁻e⁺e⁻, τ⁻ → e⁻μ⁺μ⁻, etc.) in the abstract so the reader immediately sees which four modes receive the new world-best limits.
- [Systematics section] Add a short paragraph or table summarizing the main sources of systematic uncertainty on efficiency and background yield, with their magnitudes.
Simulated Author's Rebuttal
We thank the referee for the careful reading of our manuscript and the constructive comment. We address the point below and have updated the manuscript to improve the description of our validation procedures.
read point-by-point responses
-
Referee: [Event selection and BDT description] The quoted upper limits depend on the product of signal efficiency and background normalization after the BDT cut and inclusive tagging. Both quantities are taken from Monte Carlo; the manuscript does not describe an independent data-driven validation of the BDT response or PID performance in a signal-like kinematic region. Any systematic mismatch in low-momentum lepton modeling would rescale the expected background or efficiency and directly affect the reported limits at the 10–20 % level typical for these analyses.
Authors: We thank the referee for highlighting this important aspect of the analysis. The efficiencies and background yields are indeed obtained from Monte Carlo simulation after the inclusive tagging and BDT selection. While the manuscript already contains a brief description of the BDT training and the PID criteria, we agree that an explicit discussion of data-driven cross-checks would strengthen the presentation. We have therefore added a new paragraph in the event-selection section that details the validation of PID efficiencies using control samples in data (e.g., radiative Bhabha and two-photon events for electrons, and J/ψ → μμ decays for muons) in the low-momentum regime relevant to the signal. We also compare the BDT output distribution in background-enriched sidebands between data and simulation, finding agreement within the assigned systematic uncertainties. These additions make the reliance on Monte Carlo modeling more transparent and quantify the associated systematic effects on the final limits. revision: yes
Circularity Check
No significant circularity in upper-limit derivation from data
full rationale
The paper reports a search for rare tau decays using 428 fb^{-1} of Belle II data. Event selection employs inclusive tagging plus a boosted decision tree, after which observed yields, background estimates, and signal efficiencies are combined via standard statistical methods to extract 90% CL upper limits on branching fractions. None of the load-bearing steps reduce by construction to the final limits themselves: efficiencies and backgrounds are obtained from independent Monte Carlo simulation (with data-driven validation where described), and the limits are not used to define or fit any input parameter. No self-citation chain, ansatz smuggling, or renaming of known results is invoked to justify the central result. The derivation is therefore self-contained against external benchmarks and receives the default non-circularity finding.
Axiom & Free-Parameter Ledger
axioms (1)
- domain assumption Background processes and detector response are accurately modeled in Monte Carlo simulation for the purpose of training the boosted decision tree and estimating efficiencies.
Reference graph
Works this paper leans on
- [1]
-
[2]
G. Hern´ andez-Tom´ e, G. L´ opez Castro and P. Roig,Flavor violating leptonic decays of τ and µ leptons in the standard model with massive neutrinos , Eur. Phys. J. C 79 (2019) 84 [1807.06050]
-
[3]
P. Blackstone, M. Fael and E. Passemar, τ → µµµ at a rate of one out of 1014 tau decays?, Eur. Phys. J. C 80 (2020) 506 [ 1912.09862]
- [4]
- [5]
-
[6]
Flavour physics of leptons and dipole moments
M. Raidal et al., Flavour physics of leptons and dipole moments , Eur. Phys. J. C 57 (2008) 13 [0801.1826]
work page internal anchor Pith review Pith/arXiv arXiv 2008
-
[7]
Theoretical aspects of charged Lepton Flavour Violation
A.M. Teixeira, Theoretical aspects of charged lepton flavour violation , J. Phys. Conf. Ser. 888 (2017) 012029 [ 1612.05561]
work page internal anchor Pith review Pith/arXiv arXiv 2017
-
[8]
On the Importance of Electroweak Corrections for B Anomalies
F. Feruglio, P. Paradisi and A. Pattori, On the importance of electroweak corrections for B anomalies, JHEP 09 (2017) 061 [ 1705.00929]
work page internal anchor Pith review Pith/arXiv arXiv 2017
-
[9]
Model-independent Analysis of Lepton Flavour Violating Tau Decays
B.M. Dassinger, T. Feldmann, T. Mannel and S. Turczyk, Model-independent analysis of lepton flavour violating tau decays , JHEP 10 (2007) 039 [ 0707.0988]
work page internal anchor Pith review Pith/arXiv arXiv 2007
-
[10]
Higgs-Mediated tau --> mu and tau --> e transitions in II Higgs doublet Model and Supersymmetry
P. Paradisi, Higgs-mediated τ → µ and τ → e transitions in II Higgs doublet model and supersymmetry, JHEP 02 (2006) 050 [ hep-ph/0508054]
work page internal anchor Pith review Pith/arXiv arXiv 2006
- [11]
- [12]
- [13]
-
[14]
Belle collaboration, Search for Lepton Flavor Violating Tau Decays into Three Leptons with 719 Million Produced τ +τ − Pairs, Phys. Lett. B 687 (2010) 139 [ 1001.3221]. – 15 –
work page internal anchor Pith review Pith/arXiv arXiv 2010
-
[15]
BaBar collaboration, Limits on tau lepton-flavor violating decays in three charged leptons , Phys. Rev. D 81 (2010) 111101 [ 1002.4550]
work page internal anchor Pith review Pith/arXiv arXiv 2010
-
[16]
Belle II collaboration, Belle II technical design report , 1011.0352
work page internal anchor Pith review Pith/arXiv arXiv
- [17]
-
[18]
S. Banerjee, B. Pietrzyk, J.M. Roney and Z. W¸ as, Tau and muon pair production cross-sections in electron-positron annihilations at √s = 10.58 GeV, Phys. Rev. D 77 (2008) 054012 [0706.3235]
work page internal anchor Pith review Pith/arXiv arXiv 2008
-
[19]
Particle Data Group collaboration, Review of Particle Physics , PTEP 2020 (2020) 083C01
work page 2020
-
[20]
The Precision Monte Carlo Event Generator KK For Two-Fermion Final States In e+e- Collisions
S. Jadach, B.F.L. Ward and Z. W¸ as,The precision Monte Carlo event generator KK for two-fermion final states in e+e− collisions, Comput. Phys. Commun. 130 (2000) 260 [hep-ph/9912214]
work page internal anchor Pith review Pith/arXiv arXiv 2000
-
[21]
S. Jadach, J.H. Kuhn and Z. W¸ as, TAUOLA: A library of Monte Carlo programs to simulate decays of polarized tau leptons , Comput. Phys. Commun. 64 (1990) 275
work page 1990
-
[22]
E. Barberio, B. van Eijk and Z. W¸ as, PHOTOS: A universal Monte Carlo for QED radiative corrections in decays, Comput. Phys. Commun. 66 (1991) 115
work page 1991
-
[23]
T. Sj¨ ostrand, S. Ask, J.R. Christiansen, R. Corke, N. Desai, P. Ilten et al., An Introduction to PYTHIA 8.2 , Comput. Phys. Commun. 191 (2015) 159 [ 1410.3012]
work page internal anchor Pith review Pith/arXiv arXiv 2015
-
[24]
Lange, The EvtGen particle decay simulation package , Nucl
D.J. Lange, The EvtGen particle decay simulation package , Nucl. Instrum. Meth.A 462 (2001) 152
work page 2001
-
[25]
G. Balossini, C.M. Carloni Calame, G. Montagna, O. Nicrosini and F. Piccinini, Matching perturbative and parton shower corrections to Bhabha process at flavour factories , Nucl. Phys.B 758 (2006) 227
work page 2006
-
[26]
G. Balossini, C. Bignamini, C.M.C. Calame, G. Montagna, O. Nicrosini and F. Piccinini, Photon pair production at flavour factories with per mille accuracy , Phys. Lett.B 663 (2008) 209
work page 2008
-
[27]
C.M. Carloni Calame, G. Montagna, O. Nicrosini and F. Piccinini, The BABAYAGA event generator, Nucl. Phys. B Proc. Suppl. 131 (2004) 48
work page 2004
-
[28]
Carloni Calame, An improved parton shower algorithm in QED , Phys
C.M. Carloni Calame, An improved parton shower algorithm in QED , Phys. Lett.B 520 (2001) 16
work page 2001
-
[29]
C.M. Carloni Calame, C. Lunardini, G. Montagna, O. Nicrosini and F. Piccinini, Large angle Bhabha scattering and luminosity at flavor factories , Nucl. Phys.B 584 (2000) 459
work page 2000
-
[30]
F. Berends, P. Daverveldt and R. Kleiss, Radiative corrections to the process e+e− → e+e−µ+µ−, Nucl. Phys.B 253 (1985) 421
work page 1985
-
[31]
F. Berends, P. Daverveldt and R. Kleiss, Complete lowest-order calculations for four-lepton final states in electron-positron collisions , Nucl. Phys.B 253 (1985) 441
work page 1985
-
[32]
F. Berends, P. Daverveldt and R. Kleiss, Monte Carlo simulation of two-photon processes: II: Complete lowest order calculations for four-lepton production processes in electron-positron collisions, Comp. Phys. Commun. 40 (1986) 285
work page 1986
-
[33]
S. Uehara, TREPS: A Monte-Carlo event generator for two-photon processes at e+e− colliders using an equivalent photon approximation , 1310.0157. – 16 –
work page internal anchor Pith review Pith/arXiv arXiv
-
[34]
Belle II Framework Software Group, The Belle II Core Software , Comput. Softw. Big Sci. 3 (2019) 1 [ 1809.04299]
work page internal anchor Pith review Pith/arXiv arXiv 2019
-
[35]
Belle II Analysis Software Framework (basf2)
Belle II collaboration, “Belle II Analysis Software Framework (basf2).” https://doi.org/10.5281/zenodo.5574115
-
[36]
GEANT4 collaboration, GEANT4: A simulation toolkit , Nucl.Instrum.Meth.A 506 (2003) 250
work page 2003
- [37]
-
[38]
Farhi, A QCD Test for Jets , Phys
E. Farhi, A QCD Test for Jets , Phys. Rev. Lett. 39 (1977) 1587
work page 1977
- [39]
- [40]
-
[41]
Particle Data Group collaboration, Review of particle physics , Phys. Rev. D 110 (2024) 030001
work page 2024
- [42]
-
[43]
CLEO collaboration, Search for exclusive charmless hadronic B decays , Phys. Rev. D 53 (1996) 1039 [ hep-ex/9508004]
work page internal anchor Pith review Pith/arXiv arXiv 1996
-
[44]
T. Chen and C. Guestrin, XGBoost: A scalable tree boosting system , in Proceedings of the 22nd ACM SIGKDD International Conference on Knowledge Discovery and Data Mining , KDD ’16, ACM, Aug., 2016, DOI
work page 2016
-
[45]
Optuna: A Next-generation Hyperparameter Optimization Framework
T. Akiba, S. Sano, T. Yanase, T. Ohta and M. Koyama, Optuna: A Next-generation Hyperparameter Optimization Framework, 1907.10902
work page internal anchor Pith review Pith/arXiv arXiv 1907
-
[46]
Skwarnicki, A study of the radiative CASCADE transitions between the Υ′ and Υ resonances, Ph.D
T. Skwarnicki, A study of the radiative CASCADE transitions between the Υ′ and Υ resonances, Ph.D. thesis, Cracow, INP, 1986
work page 1986
-
[47]
Belle II collaboration, Test of light-lepton universality in τ decays with the Belle II experiment, JHEP 08 (2024) 205
work page 2024
- [48]
-
[49]
Confidence Level Computation for Combining Searches with Small Statistics
T. Junk, Confidence level computation for combining searches with small statistics , Nucl. Instrum. Meth. A 434 (1999) 435 [ hep-ex/9902006]
work page internal anchor Pith review Pith/arXiv arXiv 1999
-
[50]
Read, Presentation of search results: The CL s technique, J
A.L. Read, Presentation of search results: The CL s technique, J. Phys. G 28 (2002) 2693. – 17 – Additional Material The ( Meℓℓ,∆Eeℓℓ) plane for the signal, the lepton pair invariant mass and visible energy of the event, and the Meℓℓ and ∆ Eeℓℓ distributions, for the modes omitted throughout the paper, are shown in Fig. 6, Fig. 7, and Fig. 8, respectively...
work page 2002
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.