{"id":"ce53e091-42af-47f6-9d23-7bb713cf1aba","arxiv_id":"2412.16470","paper_version":1,"verdict":"ACCEPT","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"low","formal_verification":"none","parameter_count":7,"one_line_summary":"First search for B0 to K_S^0 tau± l∓ shows no signal and sets 90% CL upper limits on branching fractions of (0.8 to 3.6) times 10^-5.","lead":"Physicists searched for two never-before-seen decays in which a B meson turns into a kaon, a tau lepton, and a muon or electron, a process the Standard Model forbids. No such decays were found, and the new limits are the best yet for the electron channels and competitive for the muon channels.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The headline limits inherit the uniform-phase-space signal model; non-uniform BSM kinematics could shift the Table I ULs by up to ~2-4x, but the disclosed efficiency maps keep the ACCEPT verdict intact.","rationale":"I read the paper in good faith. It is a carefully executed first search: the hadronic B-tagging chain is calibrated with B0→D−π+ (RFEI), the signal PDF and BDT selections are calibrated with B0→Ds+D−, and the measured B(B0→Ds+D−) = (10.1±1.2)x10^-3 agrees with the world average within 2σ, providing a closure test. The internal arithmetic checks out: the quoted central values follow from the Nsig and efficiencies via Eq. 2 (e.g., OSμ: -1.8/(1.7x10^-4 x 2 x 1.159x10^9 x 0.487) ≈ -0.94x10^-5), and the ULs are consistent with the pseudo-experiment procedure including negative fluctuations. I found no internal inconsistency or statistical-treatment error. The genuinely load-bearing assumption is the uniform-phase-space signal model, exactly as the reader identified; the efficiency maps in the supplemental material are a credible mitigation but are in reconstructed variables and do not cover the Mτ signal-shape dependence. Since the assumption is stated explicitly and the field-standard mitigation (maps) is provided, this does not change the reader's ACCEPT verdict.","tokens_in":14395,"tokens_out":23704,"duration_ms":215068,"concrete_test":"Recompute the 90% CL UL for one channel, e.g., OSμ, under a non-uniform signal Dalitz distribution: reweight the phase-space signal MC to a benchmark BSM operator from Refs. [3-5] (or, as a bounding exercise, weight events to the lowest-efficiency 25% of the supplemental (M²τℓ, M²KS0ℓ) map), regenerate the Mτ signal PDF, and rerun the pseudo-experiment limit procedure described in the text. If the resulting UL exceeds the Table I value by more than ~40%, the headline limit range is materially model-dependent and the phase-space caveat should be raised to the abstract; otherwise the concern does not land.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim - 90% CL upper limits of [0.8, 3.6]x10^-5 on B(B0→KS0 τ±ℓ∓) - rests on the signal model stated in the simulation section: 'The B0→KS0 τ±ℓ∓ signal channels are modeled using an uniform three-body phase space model.' The supplemental efficiency maps show the per-event acceptance varies from near zero to about 8x10^-4 across the (M²τℓ, M²KS0ℓ) plane, whereas the phase-space-averaged efficiencies in Table I are only 1.7-2.1x10^-4. A plausible BSM operator (e.g., the b→sτℓ models in Refs. [3-5]) that populates a different Dalitz region would therefore change the signal efficiency, and hence the UL, by up to roughly a factor of two to four, larger than the 22-24% total systematic uncertainty quoted. In the non-conservative direction (signal concentrated in low-efficiency regions), the true 90% CL limit could be several times weaker than the headline numbers, so the claim that true branching fractions lie below [0.8, 3.6]x10^-5 is strictly valid only for the phase-space model. The paper mitigates this by providing efficiency maps, which is a real mitigation: the limit for any model can in principle be obtained by reweighting. However, the maps are binned in reconstructed/inferred variables (the τ momentum is inferred via the Btag), so a reinterpretation requires detector-level simulation, and the Mτ signal-PDF shape used in the profile fit stays fixed to the phase-space sample, an effect the maps do not encode. This is a disclosed, standard model-dependence rather than an internal inconsistency; it is the correct target for a stress test but does not invalidate the analysis as presented.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper reports the first search for the lepton flavor-violating decays B0→KS0τ±ℓ∓ (ℓ=μ,e), using 711 fb−1 of Belle and 365 fb−1 of Belle II data. Events are selected with a hadronic FEI B-tag; the signal is identified via the recoil mass Mτ reconstructed from the Btag, the KS0, and a primary lepton, with τ candidates reconstructed in the e, μ, π, and ρ decay modes. An unbinned maximum-likelihood fit to Mτ is performed on the combined data sample; no significant signal is found, and 90% CL upper limits on the branching fractions in [0.8, 3.6]×10−5 are set. The analysis is calibrated on control channels: B0→D−π+ for the tag-efficiency scale factors R_FEI, B0→Ds−D+ for the signal PDF and BDT efficiency ratios R_BDT, with a closure test of B(B0→Ds−D+) consistent with the world average. Limit-setting uses pseudo-experiments with systematic smearing.","tokens_in":14807,"tokens_out":22076,"duration_ms":193283,"significance":"If it holds, this result provides the first constraints on B0→KS0τ±ℓ∓; the electron modes are the most stringent limits on b→sτe transitions, and the muon modes are competitive with the best existing b→sτμ bounds. The limits approach the O(10−6) enhancements suggested by the BSM models cited in the paper. The strengths of the analysis are its data-driven calibration chain with an independent closure test, pseudo-experiment validation of the fitting procedure, and the Supplemental efficiency maps as functions of (M2τℓ, M2KS0ℓ), which provide a genuine route toward model reinterpretation. The principal caveat is the uniform-phase-space signal model underlying the headline limits, assessed in the major comments.","major_comments":[{"comment":"The quoted 90% CL upper limits in Table I and in the abstract and conclusion are derived from a signal model with a uniform three-body phase-space distribution, as stated in the Simulation section ('The B0→KS0τ±ℓ∓ signal channels are modeled using an uniform three-body phase space model'). This assumption is load-bearing: the Supplemental efficiency maps show that the per-event acceptance varies from near zero to about 8×10−4 across the (M2τℓ, M2KS0ℓ) plane, whereas the phase-space-averaged efficiencies in Table I are only 1.7–2.1×10−4. A BSM signal with a different Dalitz distribution (such as the models of Refs. [3–5]) could therefore change the signal efficiency, and hence the limits, by roughly a factor of two to four, which is larger than the quoted 22–24% total systematic uncertainty. The efficiency maps are a useful mitigation, but they are binned in reconstructed variables (the τ momentum is inferred from the Btag), so a reinterpretation requires detector-level simulation, and the Mτ signal-PDF shape used in the profile fit is fixed to the phase-space sample, an effect the maps do not encode. I request: (i) an explicit qualification in the abstract and conclusion that the limits assume a uniform three-body phase-space model, and (ii) preferably, a quantitative robustness study in which the signal MC is reweighted to a few representative non-uniform Dalitz models and the efficiencies (or, if feasible, the full fit and ULs) are recomputed, so that the reader can see the range of possible limit shifts.","section":"Simulation section; Table I; Supplemental Fig. 1"}],"minor_comments":[{"comment":"The phrase 'an uniform three-body phase space model' should read 'a uniform three-body phase-space model'; please also use 'phase-space' consistently throughout.","section":"Simulation"},{"comment":"The text quotes the R_FEI uncertainty as '4%', but the quoted values 0.74±0.04 and 0.81±0.04 correspond to relative uncertainties of 5.4% and 4.9%; the sentence should be corrected.","section":"Systematic uncertainties"},{"comment":"The horizontal-axis label of the lower-right panel in Figure 2 appears garbled ('M2 −0' ...); please verify the typesetting.","section":"Fig. 2"},{"comment":"Since the best-fit branching-fraction central values are negative in three of the four channels, the column label 'central value' is confusing; 'best-fit value' would be more accurate.","section":"Table I"},{"comment":"Reference [10] combines two distinct publications (Kurokawa and Kikutani; Abe et al.); they should be cited separately.","section":"Introduction, Ref. [10]"},{"comment":"The caption of Supplemental Fig. 1 should state explicitly that the τ four-momentum used for M2τℓ is the value inferred from the Btag reconstruction, as this information currently appears only in the body text.","section":"Supplemental material"},{"comment":"The definitions of 'same-sign' (SSℓ) and 'opposite-sign' (OSℓ) should be clarified in one sentence, since τ±ℓ∓ pairs have opposite charges in both cases.","section":"Signal categories"},{"comment":"The expected ULs are quoted only as the range [2.1, 2.2]×10−5; listing the expected limit for each channel in Table I would help the reader assess the observed SSe limit (0.8×10−5), which is a downward fluctuation relative to expectation.","section":"Results"}],"recommendation":"major_revision","confidential_remarks":"The analysis is careful and the internal consistency checks are convincing; my main substantive request is the qualification or robustness study of the phase-space signal model, which affects the literal reading of the headline limits. I also note that the figures are marked 'preliminary' and the collaboration's approval status should be confirmed before publication. The Supplemental efficiency maps are a strength and should be retained."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The short version: this is a genuine first search, carefully done, and the headline limits stand as stated under the phase-space signal model they use. The phase-space assumption is the only real soft spot, and it is disclosed and mitigated.\n\nWhat's new: no one had looked for B0 -> KS0 tau± l∓ before, and this is also the first direct LFV search in B decays using Belle II data. Adding the tau -> rho nu channel is a nice technical step that previous B -> K tau l analyses hadn't used.\n\nWhat it does well: the calibration chain is solid. They use B0 -> D- pi+ for the tagging efficiency, B0 -> D_s^- D+ to fix the signal PDF and BDT selection, and get a closure test consistent with the world average. The fit validation with pseudo-experiments is standard and fine. The systematic budget is honest, with the dominant terms (BDT and signal PDF) estimated from control samples. The supplemental efficiency maps in (M^2_tau_l, M^2_KS_l) are exactly the right thing to provide.\n\nThe soft spot: the signal is generated with a uniform three-body phase space. The maps show per-event acceptance varies by an order of magnitude over the Dalitz plane, so a non-uniform BSM signal could shift the 90% CL ULs by maybe a factor of 2-4 relative to Table I. That is larger than the quoted ~22-24% systematic uncertainties. The maps let you reweight, but the tau momentum is inferred via the Btag, so reinterpretation requires detector-level simulation, and the M_tau signal shape in the fit stays tied to the phase-space sample. This is a disclosed and fairly standard model-dependence, not an internal inconsistency. The paper should make the caveat more prominent in the abstract or conclusion.\n\nOverall: the central argument holds up. I'd take the limits at face value as phase-space-model limits, which is how they are framed.","headline":"First search for B0 -> KS0 tau l is cleanly executed and the null-result limits are credible as stated, with a disclosed phase-space model dependence that the efficiency maps help mitigate.","tokens_in":17776,"tokens_out":2319,"would_cite":true,"duration_ms":21773,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"The first search for $B^0\\to K_S^0\\tau^\\pm\\ell^\\mp$ decays finds no signal and sets 90% confidence-level upper limits on the branching fractions in the range $[0.8,3.6]\\times10^{-5}$.","keywords":["lepton flavor violation","B meson decays","tau lepton","hadronic B tagging","recoil mass","upper limits","Belle II","b to s tau ell transitions"],"falsifier":"A future dataset with substantially more integrated luminosity would settle the claim: if the reconstructed $\\tau$ mass distribution shows a peak at the $\\tau$ mass in any of the four channels with a yield incompatible with background, the no-signal claim is wrong; alternatively, recomputing the limits with a specific non-uniform signal model using the published efficiency maps would show how much the limits shift.","tokens_in":14215,"feed_emoji":"🔍","tokens_out":8932,"duration_ms":69254,"temperature":0.7,"pith_summary":"The paper reports the first search for the lepton-flavor-violating decays $B^0\\to K_S^0\\tau^\\pm\\ell^\\mp$ with $\\ell=\\mu,e$, processes that the Standard Model forbids but that several beyond-Standard-Model scenarios, motivated by anomalies in semileptonic $B$ decays, predict could be enhanced to observable rates. The analysis combines 711 fb$^{-1}$ of data from Belle with 365 fb$^{-1}$ from Belle II, fully reconstructing one $B$ meson and searching for the signal in the recoiling system through the reconstructed $\\tau$ mass. No evidence for any of the four charge-flavor channels is found, and 90% confidence-level upper limits on the branching fractions are set in the range $[0.8,3.6]\\times10^{-5}$. If the true branching fractions lie below these limits, models of lepton flavor violation in $b\\to s\\tau\\ell$ transitions become more constrained, and the limits are within an order of magnitude of the $10^{-6}$ level some models predict.","feed_headline":"No signal in first search for B0 to KS tau-lepton decays","feed_subtitle":"Limits of (0.8–3.6)×10^-5 on b→sτℓ branching fractions narrow the room for beyond-Standard-Model lepton flavor violation.","key_machinery":"The central object is the reconstructed recoil mass $M_\\tau$, obtained from the beam energy and the momentum of the fully reconstructed other $B$ meson: since the signal $B$ is not fully visible (the $\\tau$ decays with neutrinos), its momentum is taken as $\\vec p_\\tau = -\\vec p_{B_{\\rm tag}} - \\vec p_{K_S^0} - \\vec p_\\ell$ and its energy as $E_\\tau = E_{\\rm beam} - E_{K_S^0} - E_\\ell$. Signal events accumulate at the known $\\tau$ mass while backgrounds are smooth. The tagging itself uses the full-event-interpretation algorithm, a machine-learning $B$-tagging tool with an average efficiency of 0.59% and 44% purity, calibrated on control samples ($B^0\\to D^-\\pi^+$ for the tag efficiency and $B^0\\to D_s^+D^-$ for the signal PDF and the boosted-decision-tree selection). Signal decays are modeled with a uniform three-body phase-space distribution.","core_discovery":"On its own terms, the paper establishes that $B^0\\to K_S^0\\tau^+\\mu^-$, $B^0\\to K_S^0\\tau^-\\mu^+$, $B^0\\to K_S^0\\tau^+e^-$, and $B^0\\to K_S^0\\tau^-e^+$ do not appear in the combined Belle and Belle II data at rates above a few times $10^{-5}$. The signal is searched for in the mass recoiling against a hadronically tagged $B$ meson; the reconstructed quantity $M_\\tau$ peaks at the known $\\tau$ mass for signal and is smooth for background. With no significant signal in any channel, the authors set observed upper limits at 90% CL of $1.1\\times10^{-5}$, $3.6\\times10^{-5}$, $1.5\\times10^{-5}$, and $0.8\\times10^{-5}$, respectively, and report expected limits of $[2.1,2.2]\\times10^{-5}$. They also supply efficiency maps in ($M_{\\tau\\ell}^2$, $M_{K_S^0\\ell}^2$) so the limits can be reinterpreted for signal models that do not follow a uniform three-body phase space.","pith_inferences":["If the Belle II $B^+\\to K^+\\nu\\bar\\nu$ excess is due to new physics with off-diagonal lepton-flavor couplings, predictions in that framework place $B(B\\to K\\tau^\\pm\\mu^\\mp)$ near $[2,3]\\times10^{-6}$; the limits here are roughly a factor of 5-10 above that, so with several times more Belle II data the same technique would either observe the decay or exclude the simplest version of that explanation","The uniform phase-space assumption is likely conservative for signals that populate the high-efficiency region of phase space, but could understate limits for signals concentrated elsewhere; the published efficiency maps make this testable, and a recast for vector, scalar, or tensor operators would quantify the shift.","The inclusion of the $\\tau\\to\\rho\\nu$ decay mode, which carries more than 20% of the $\\tau$ branching fraction and had not been used in earlier $B\\to K\\tau\\ell$ searches, suggests that similar searches for related modes such as $B^0\\to K^{*0}\\tau^\\pm\\ell^\\mp$ or $B_s$ decays could gain sensitivity by the same route."],"forward_implications":["These are the first upper limits on $B^0\\to K_S^0\\tau^\\pm\\ell^\\mp$; no other experiment has searched for these modes.","The limits on $B^0\\to K_S^0\\tau^\\pm e^\\mp$ are the most stringent on $b\\to s\\tau e$ transitions to date.","The $B^0\\to K_S^0\\tau^\\pm\\mu^\\mp$ limits are among the best on $b\\to s\\tau\\mu$ transitions achieved so far.","The results approach the $\\mathcal{O}(10^{-6})$ branching-fraction level predicted by some models motivated by the $B^+\\to K^+\\nu\\bar\\nu$ excess, so the search is beginning to test those scenarios.","The efficiency maps provided in the supplemental material let the limits be recast for arbitrary signal kinematics, not just the uniform phase-space model."],"supporting_citations":[{"why":"Motivates the search by reporting anomalies in semileptonic $B$ decays such as $b\\to c\\tau\\nu$ that hint at lepton flavor universality violation.","marker":"[1]"},{"why":"Reports the Belle II $B^+\\to K^+\\nu\\bar\\nu$ excess that, if confirmed, could reflect off-diagonal lepton-flavor couplings.","marker":"[2]"},{"why":"Predicts enhanced $b\\to s\\tau\\ell$ branching fractions in BSM models with couplings to third-generation quarks and leptons.","marker":"[3]"},{"why":"Predicts $B(B\\to K\\tau^\\pm\\mu^\\mp)$ at $[2,3]\\times10^{-6}$ from the $B^+\\to K^+\\nu\\bar\\nu$ excess, setting the target sensitivity.","marker":"[4]"},{"why":"Reports the first LFV searches in $B^+\\to K^+\\tau^\\pm\\ell^\\mp$, providing the prior limits this analysis extends.","marker":"[6]"},{"why":"Gives the most stringent previous limit on $B^+\\to K^+\\tau^+\\mu^-$, the benchmark for the muon channels.","marker":"[7]"},{"why":"Sets limits on the related mode $B^0\\to K^{*0}\\tau^\\pm\\mu^\\mp$, a comparison point for $b\\to s\\tau\\mu$ constraints.","marker":"[9]"},{"why":"Supplies the full-event-interpretation algorithm used for hadronic $B$-tagging, the backbone of the recoil-mass method.","marker":"[22]"},{"why":"Provides the world-average $\\tau$ mass, $B^0$ mass, and branching fractions used to define the signal PDF and efficiency.","marker":"[23]"},{"why":"Provides the $f_{+-}/f_{00}$ ratio of $\\Upsilon(4S)$ decays used to convert yields to branching fractions.","marker":"[34]"}],"fun_headline_variants":["No lepton flavor violation in B0→KS τ decays, say Belle & Belle II","First search for B0→KS τℓ decays finds nothing, sets limits","B0→KS τℓ flavor-violating decays: null, limits 0.8–3.6×10^-5","Belle and Belle II see no B0→KS τℓ, restrict new physics"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The quoted upper limits assume that the signal decays follow a uniform three-body phase-space distribution; if a real new-physics signal were distributed differently in phase space, the efficiency and therefore the limits would change.","fun_headline_variants_meta":{"raw":{"variants":["No lepton flavor violation in B0→KS τ decays, say Belle & Belle II","First search for B0→KS τℓ decays finds nothing, sets limits","B0→KS τℓ flavor-violating decays: null, limits 0.8–3.6×10^-5","Belle and Belle II see no B0→KS τℓ, restrict new physics"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000344,"raw_usage":{"total_tokens":1913,"prompt_tokens":989,"completion_tokens":924,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":605,"completion_tokens_details":{"reasoning_tokens":825}},"tokens_in":605,"tokens_out":924,"duration_ms":8528,"temperature":1.0,"reasoning_tokens":825,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T10:32:55.400298+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A future dataset with substantially more integrated luminosity would settle the claim: if the reconstructed $\\tau$ mass distribution shows a peak at the $\\tau$ mass in any of the four channels with a yield incompatible with background, the no-signal claim is wrong; alternatively, recomputing the limits with a specific non-uniform signal model using the published efficiency maps would show how much the limits shift.","supporting_citations":[{"cited_title":"Florian et al ., Rev","cited_arxiv_id":null,"evidence_quote":"Motivates the search by reporting anomalies in semileptonic $B$ decays such as $b\\to c\\tau\\nu$ that hint at lepton flavor universality violation."},{"cited_title":"Adachi et al","cited_arxiv_id":null,"evidence_quote":"Reports the Belle II $B^+\\to K^+\\nu\\bar\\nu$ excess that, if confirmed, could reflect off-diagonal lepton-flavor couplings."},{"cited_title":"Becirevi´ c, O","cited_arxiv_id":null,"evidence_quote":"Predicts enhanced $b\\to s\\tau\\ell$ branching fractions in BSM models with couplings to third-generation quarks and leptons."},{"cited_title":"Allwicher et al ., Phys","cited_arxiv_id":null,"evidence_quote":"Predicts $B(B\\to K\\tau^\\pm\\mu^\\mp)$ at $[2,3]\\times10^{-6}$ from the $B^+\\to K^+\\nu\\bar\\nu$ excess, setting the target sensitivity."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Reports the first LFV searches in $B^+\\to K^+\\tau^\\pm\\ell^\\mp$, providing the prior limits this analysis extends."},{"cited_title":"Watanuki et al","cited_arxiv_id":null,"evidence_quote":"Gives the most stringent previous limit on $B^+\\to K^+\\tau^+\\mu^-$, the benchmark for the muon channels."},{"cited_title":"Aaij et al","cited_arxiv_id":null,"evidence_quote":"Sets limits on the related mode $B^0\\to K^{*0}\\tau^\\pm\\mu^\\mp$, a comparison point for $b\\to s\\tau\\mu$ constraints."},{"cited_title":"Keck et al ., Comput","cited_arxiv_id":null,"evidence_quote":"Supplies the full-event-interpretation algorithm used for hadronic $B$-tagging, the backbone of the recoil-mass method."}],"review_version":1}