The B^+ to K^+ ν bar ν decay as a search for the QCD axion
Pith reviewed 2026-05-18 04:36 UTC · model grok-4.3
The pith
A model-independent method reinterprets Belle II public data on B+ to K+ nu nubar to set the strongest limit on the decay to a QCD axion.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
By deriving an analytic reconstruction of the mapping between true and reconstructed kinematic variables solely from public data in the Inclusive Tagging Analysis, the B+ to K+ nu nubar decay can be reinterpreted as a search for the QCD axion, yielding the strongest bound on the branching fraction for B+ to K+ a and constraining the axion's flavour-changing coupling to b and s quarks.
What carries the argument
The analytic reconstruction of the mapping between true and reconstructed kinematic variables within the Inclusive Tagging Analysis, which allows model-independent limits on invisible particles without internal simulations.
If this is right
- The B+ to K+ nu nubar decay acts as a dual probe for short-distance new physics and light invisible states, with the two tests operating independently.
- The axion's fundamental flavour-changing coupling between b and s quarks receives a tighter constraint than before.
- The method supplies a general strategy for model-independent reinterpretation of other rare-decay measurements at colliders.
- Existing experimental limits on the branching fraction for B+ to K+ a are improved by roughly a factor of nine.
Where Pith is reading between the lines
- The same analytic reconstruction technique could be adapted to reinterpret other public datasets from Belle II or LHCb for axion-like particle searches in additional channels.
- Releasing detailed kinematic information in public data releases enables independent groups to perform new physics searches beyond the original analysis goals.
- Higher-statistics future measurements of B+ to K+ nu nubar would allow even stronger limits on the axion coupling through the same method.
Load-bearing premise
The analytic reconstruction of the kinematic mapping from public data accurately captures the dominant effects in the Inclusive Tagging Analysis without significant unaccounted detector or background biases.
What would settle it
A direct comparison of the extracted limit on B+ to K+ a using the analytic mapping against the limit from the experiment's internal full simulation or dedicated detector study would confirm or refute the accuracy of the reconstruction.
Figures
read the original abstract
We introduce a model-independent framework to reinterpret Belle II results using only public data, analytically reconstructing the mapping between true and reconstructed kinematic variables within the statistically dominant Inclusive Tagging Analysis. This enables rare-decay measurements to probe light invisible particles -- such as the QCD axion or axion-like particles, collectively denoted $a$ -- without relying on internal simulations. Applying the method to $B^+ \! \to \! K^+ \nu \bar\nu$ yields the strongest bound on the branching fraction for $B^+ \! \to \! K^+ a$, improving existing limits by about a factor of nine and constraining the axion's fundamental flavour-changing coupling to $b$ and $s$ quarks. The approach establishes $B^+ \! \to \! K^+ \nu \bar\nu$ as a dual probe -- simultaneously testing short-distance new physics and light invisible states, the two probes working independently to an excellent approximation -- and provides a general strategy for model-independent reinterpretation of collider data.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript introduces a model-independent framework to reinterpret Belle II results using only public data, analytically reconstructing the mapping between true and reconstructed kinematic variables within the statistically dominant Inclusive Tagging Analysis. This is used to probe light invisible particles such as the QCD axion in the decay B⁺ → K⁺ a. The paper claims that applying this method to B⁺ → K⁺ νν̄ yields the strongest bound on the branching fraction for B⁺ → K⁺ a, improving existing limits by about a factor of nine, and constrains the axion's fundamental flavour-changing coupling to b and s quarks. It also positions B⁺ → K⁺ νν̄ as a dual probe for short-distance new physics and light invisible states.
Significance. If the analytic mapping is shown to be robust, this result would offer a significant advancement by enabling searches for axions and ALPs using public collider data without internal simulations. The approach's strength is its model-independence and reproducibility from public information, providing a general strategy for reinterpretation of rare decay measurements. This could impact axion phenomenology and encourage similar analyses in other channels.
major comments (2)
- [§3] The analytic reconstruction of the mapping between true and reconstructed kinematic variables is the load-bearing element for the claimed improvement. The paper should include a quantitative validation, for example by comparing the reconstructed distributions or efficiency for the fixed-mass axion signal against expectations, to address potential unaccounted detector or background biases that could affect the extracted limit on BR(B⁺ → K⁺ a).
- [Results] The factor-of-nine improvement in the bound should be accompanied by a detailed breakdown of how the analytic mapping contributes to the uncertainty, including any assumptions about the Inclusive Tagging Analysis.
minor comments (2)
- Clarify the definition of the fundamental flavour-changing coupling in the context of the axion model.
- [Abstract] The abstract could specify the mass range for the axion considered in the bound.
Simulated Author's Rebuttal
We thank the referee for their careful reading of the manuscript and for the constructive comments, which have helped us improve the presentation and robustness of our results. We address each major comment below and have revised the manuscript accordingly.
read point-by-point responses
-
Referee: [§3] The analytic reconstruction of the mapping between true and reconstructed kinematic variables is the load-bearing element for the claimed improvement. The paper should include a quantitative validation, for example by comparing the reconstructed distributions or efficiency for the fixed-mass axion signal against expectations, to address potential unaccounted detector or background biases that could affect the extracted limit on BR(B⁺ → K⁺ a).
Authors: We agree that explicit quantitative validation strengthens the central claim. In the revised manuscript we have added a new subsection to §3 that presents a direct comparison of the analytically reconstructed kinematic distributions and signal efficiencies for a fixed-mass axion against the publicly available Belle II kinematic maps and efficiency parametrizations. The comparison shows agreement at the level of a few percent across the relevant mass range, with a brief discussion of how residual detector effects are absorbed into the conservative systematic uncertainty assigned to the mapping. This addition directly addresses the concern about unaccounted biases while remaining fully within the model-independent, public-data framework. revision: yes
-
Referee: [Results] The factor-of-nine improvement in the bound should be accompanied by a detailed breakdown of how the analytic mapping contributes to the uncertainty, including any assumptions about the Inclusive Tagging Analysis.
Authors: We appreciate the request for greater transparency on the uncertainty budget. The revised Results section now contains an expanded paragraph and an accompanying table that decomposes the total uncertainty on the extracted BR(B⁺ → K⁺ a) limit. The table isolates the contribution arising from the analytic mapping (propagated through the kinematic reweighting), the statistical uncertainty of the original Belle II measurement, and the dominant assumptions of the Inclusive Tagging Analysis (primarily the tagging efficiency and the modeling of the dominant background components). We explicitly state that the mapping itself is derived solely from public information and that its uncertainty is evaluated by varying the input kinematic parametrizations within their reported uncertainties. This breakdown makes the origin of the factor-of-nine improvement fully traceable. revision: yes
Circularity Check
No circularity: analytic reinterpretation from public data is self-contained
full rationale
The paper introduces an analytic reconstruction of the true-to-reconstructed kinematic mapping derived solely from public Belle II data for the Inclusive Tagging Analysis. This mapping is then applied to reinterpret the B+ → K+ νν̄ measurement as a limit on B+ → K+ a without fitting any parameters to the axion signal itself. The derivation chain therefore remains independent of the final bound and the claimed factor-of-nine improvement; it does not reduce by construction to a fitted input, self-definition, or load-bearing self-citation. The central result is obtained by applying the externally derived mapping to existing data, satisfying the requirement for a self-contained, non-circular derivation.
Axiom & Free-Parameter Ledger
axioms (1)
- domain assumption The statistically dominant Inclusive Tagging Analysis kinematics can be mapped analytically from public data without material bias from unpublicized detector effects.
Lean theorems connected to this paper
-
IndisputableMonolith/Foundation/AbsoluteFloorClosure.leanabsolute_floor_iff_bare_distinguishability unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
the mapping between q2 and q2_rec can be accurately derived from kinematic arguments alone... q2_rec ≡ s/4 + mK^2 - sqrt(s) E*_K = q2 + ... -2 |p*_B| gamma** (beta** E**_K(q2) + |p**_K(q2)| cos theta**)
What do these tags mean?
- matches
- The paper's claim is directly supported by a theorem in the formal canon.
- supports
- The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
- extends
- The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
- uses
- The paper appears to rely on the theorem as machinery.
- contradicts
- The paper's claim conflicts with a theorem or certificate in the canon.
- unclear
- Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.
Forward citations
Cited by 2 Pith papers
-
Braking protons at the EIC: from invisible meson decay to new physics searches
The EIC can probe invisible pseudoscalar meson decays down to branching ratios of 10^{-8} and invisibly decaying ALPs with couplings up to 10^5 GeV for masses 0.1-2 GeV.
-
Dark Matter emission at Belle II and NA62 in Minimal Flavor Violation framework
A single nearly degenerate dark matter multiplet in the MFV framework can accommodate either the K+ to pi+ nu nubar or B+ to K+ nu nubar excess but not both simultaneously.
Reference graph
Works this paper leans on
-
[1]
I. G. Irastorza and J. Redondo, New experimental ap- proaches in the search for axion-like particles, Prog. Part. Nucl. Phys. 102, 89 (2018), arXiv:1801.08127 [hep-ph]
work page internal anchor Pith review Pith/arXiv arXiv 2018
-
[2]
G. Lanfranchi, M. Pospelov, and P. Schuster, The Search for Feebly Interacting Particles, Ann. Rev. Nucl. Part. Sci. 71, 279 (2021), arXiv:2011.02157 [hep-ph]
-
[3]
Invisible Axion Search Methods,
P. Sikivie, Invisible Axion Search Methods, Rev. Mod. Phys. 93, 015004 (2021), arXiv:2003.02206 [hep-ph]
-
[4]
G. G. Raffelt, Astrophysical methods to constrain axions and other novel particle phenomena, Phys. Rept.198, 1 (1990)
work page 1990
-
[5]
G. G. Raffelt, Astrophysical axion bounds, Lect. Notes Phys. 741, 51 (2008), arXiv:hep-ph/0611350
work page internal anchor Pith review Pith/arXiv arXiv 2008
-
[6]
A. Caputo and G. Raffelt, Astrophysical Axion Bounds: The 2024 Edition, in 1st Training School of the COST Action COSMIC WISPers (CA21106) (2024) arXiv:2401.13728 [hep-ph]
-
[7]
D. J. E. Marsh, Axion Cosmology, Phys. Rept.643, 1 (2016), arXiv:1510.07633 [astro-ph.CO]
work page internal anchor Pith review Pith/arXiv arXiv 2016
- [8]
-
[9]
The landscape of QCD axion models
L. Di Luzio, M. Giannotti, E. Nardi, and L. Visinelli, The landscape of QCD axion models, Phys. Rept.870, 1 (2020), arXiv:2003.01100 [hep-ph]
work page internal anchor Pith review Pith/arXiv arXiv 2020
- [10]
-
[11]
J. E. Kim, Weak Interaction Singlet and Strong CP In- variance, Phys. Rev. Lett.43, 103 (1979)
work page 1979
-
[12]
M. A. Shifman, A. I. Vainshtein, and V. I. Zakharov, Can Confinement Ensure Natural CP Invariance of Strong In- teractions?, Nucl. Phys. B166, 493 (1980)
work page 1980
-
[13]
M. Dine, W. Fischler, and M. Srednicki, A Simple Solu- tion to the Strong CP Problem with a Harmless Axion, Phys. Lett. B104, 199 (1981)
work page 1981
-
[14]
A. R. Zhitnitsky, On Possible Suppression of the Axion Hadron Interactions. (In Russian), Sov. J. Nucl. Phys. 31, 260 (1980)
work page 1980
-
[15]
R. D. Peccei and H. R. Quinn, CP Conservation in the Presence of Instantons, Phys. Rev. Lett.38, 1440 (1977)
work page 1977
-
[16]
R. D. Peccei and H. R. Quinn, Constraints Imposed by CP Conservation in the Presence of Instantons, Phys. Rev. D 16, 1791 (1977)
work page 1977
-
[17]
Weinberg, A New Light Boson?, Phys
S. Weinberg, A New Light Boson?, Phys. Rev. Lett.40, 223 (1978)
work page 1978
-
[18]
Wilczek, Problem of StrongP andT Invariance in the Presence of Instantons, Phys
F. Wilczek, Problem of StrongP andT Invariance in the Presence of Instantons, Phys. Rev. Lett.40, 279 (1978). 6
work page 1978
-
[19]
L.Ubaldi,Effectsofthetaonthedeuteronbindingenergy and the triple-alpha process, Phys. Rev. D81, 025011 (2010), arXiv:0811.1599 [hep-ph]
work page internal anchor Pith review Pith/arXiv arXiv 2010
-
[20]
M. Dine, L. Stephenson Haskins, L. Ubaldi, and D. Xu, Some Remarks on Anthropic Approaches to the Strong CP Problem, JHEP05, 171, arXiv:1801.03466 [hep-th]
work page internal anchor Pith review Pith/arXiv arXiv
-
[21]
J. Preskill, M. B. Wise, and F. Wilczek, Cosmology of the Invisible Axion, Phys. Lett. B120, 127 (1983)
work page 1983
-
[22]
L. F. Abbott and P. Sikivie, A Cosmological Bound on the Invisible Axion, Phys. Lett. B120, 133 (1983)
work page 1983
-
[23]
M. Dine and W. Fischler, The Not So Harmless Axion, Phys. Lett. B120, 137 (1983)
work page 1983
-
[24]
Belle II Technical Design Report
T. Abeet al. (Belle-II), Belle II Technical Design Report, (2010), arXiv:1011.0352 [physics.ins-det]
work page internal anchor Pith review Pith/arXiv arXiv 2010
-
[25]
The spread, or bite, around these momenta is tiny and completely negligible for the purposes of this paper
-
[26]
The conjugated mode with negatively-charged mesons is assumed to be included throughout our study
- [27]
-
[28]
The sign difference in thefa-definingaG ˜G term in Eq.1 with respect to Ref.[9] arises from our use of opposite ϵ-tensor conventions
-
[29]
N. Gubernari, M. Reboud, D. van Dyk, and J. Virto, Dispersive analysis ofB→K(∗) and Bs→ϕform fac- tors, JHEP 12, 153, [Erratum: JHEP 01, 125 (2025)], arXiv:2305.06301 [hep-ph]
-
[30]
S. Navaset al. (Particle Data Group), Review of particle physics, Phys. Rev. D110, 030001 (2024)
work page 2024
- [31]
-
[32]
Adachiet al., Evidence for B+ → K+ν¯ν decays, Phys
I. Adachiet al. (Belle-II), Evidence for B+→K+νν¯ de- cays, Phys. Rev. D109, 112006 (2024), arXiv:2311.14647 [hep-ex]
- [33]
-
[34]
Abudinén et al., Search for B+ → K+ν¯ν decays using an inclusive tagging method at Belle II , Phys
F. Abudinén et al. (Belle-II), Search for B+→K+νν¯ Decays Using an Inclusive Tagging Method at Belle II, Phys. Rev. Lett. 127, 181802 (2021), arXiv:2104.12624 [hep-ex]
- [35]
-
[36]
D. Bečirević, G. Piazza, and O. Sumensari, Revisiting B→K(∗)ν¯νdecays in the Standard Model and beyond, Eur. Phys. J. C83, 252 (2023), arXiv:2301.06990 [hep- ph]
-
[37]
A. J. Buras, J. Girrbach-Noe, C. Niehoff, and D. M. Straub, B→K(∗)ννdecays in the Standard Model and beyond, JHEP 02, 184, arXiv:1409.4557 [hep-ph]
work page internal anchor Pith review Pith/arXiv arXiv
-
[38]
J. Brod, M. Gorbahn, and E. Stamou, Two-Loop Elec- troweak Corrections for theK→πν¯νDecays, Phys. Rev. D 83, 034030 (2011), arXiv:1009.0947 [hep-ph]
work page internal anchor Pith review Pith/arXiv arXiv 2011
-
[39]
V. Bertacchiet al. (Belle II Tracking Group), Track find- ing at Belle II, Comput. Phys. Commun. 259, 107610 (2021), arXiv:2003.12466 [physics.ins-det]
-
[40]
The below-2%NB error has a negligible impact on our analysis, becauseni ν¯ν,bkg have an error approaching 10%
-
[41]
K. Fridell, M. Ghosh, T. Okui, and K. Tobioka, De- coding the B→Kννexcess at Belle II: Kinematics, op- erators, and masses, Phys. Rev. D109, 115006 (2024), arXiv:2312.12507 [hep-ph]
- [42]
-
[43]
Praz, Search for B → Kvv decays with a machine learning method at the Belle II experiment , Ph.D
C. Praz, Search for B → Kvv decays with a machine learning method at the Belle II experiment , Ph.D. thesis, U. Hamburg (main), Hamburg U., Hamburg U., Dept. Math., Hamburg (2022)
work page 2022
-
[44]
This tail originates from radiative decays with unde- tected photons. Strictly speaking, B →Kν¯ν(+nγ) is a distinct process fromB →Kν¯νand should be sub- tracted if sizeable, as inBs→µ+µ−[50]. Here, however, the effect is negligible (well below 1% in Ref. [43]), since the lightest emitter is theK+, which makes the exponent b in B(B→Kν¯ν(+nγ)) = ( 2∆E mB )...
-
[45]
The Full Event Interpretation -- An exclusive tagging algorithm for the Belle II experiment
T. Keck et al. , The Full Event Interpretation: An Ex- clusive Tagging Algorithm for the Belle II Experiment, Comput. Softw. Big Sci. 3, 6 (2019), arXiv:1807.08680 [hep-ex]
work page internal anchor Pith review Pith/arXiv arXiv 2019
-
[46]
In principle, the same smearing should also be applied to the di-neutrino channel, Eq. (12). In the axion case this treatment is essential, as it broadens the Dirac-delta shape of theq2 distribution. For the di-neutrino mode, however, the impact is much milder because itsq2 spec- trum is already smooth. Such an effect would become relevant only once a mea...
-
[47]
W. Altmannshofer, A. Crivellin, H. Haigh, G. In- guglia, and J. Martin Camalich, Light new physics in B→K(*)νν¯?, Phys. Rev. D 109, 075008 (2024), arXiv:2311.14629 [hep-ph]
-
[48]
M. Cavan-Piton, D. Guadagnoli, A. Iohner, D. Mar- tinez Santos, and L. Vittorio, Probing QCD Axions or Axion-like Particles in three-body K Decays, (2024), arXiv:2411.04170 [hep-ph]
-
[49]
J. Martin Camalich, M. Pospelov, P. N. H. Vuong, R. Ziegler, and J. Zupan, Quark Flavor Phenomenology of the QCD Axion, Phys. Rev. D 102, 015023 (2020), arXiv:2002.04623 [hep-ph]
-
[50]
A. J. Buras, J. Girrbach, D. Guadagnoli, and G. Isidori, On the Standard Model prediction for BR(Bs,d to mu+ mu-), Eur. Phys. J. C72, 2172 (2012), arXiv:1208.0934 [hep-ph]. 7 ITA HTA Combined 0 2 4 6 80 2.5 5 7.5 10 12.5 15 μνν _ - 2 log (ℒ/ℒ max) FIG. 3: Twice the negative profile log-likelihood ratio as a function ofµν¯νfor the ITA (red), HTA (blue) or ...
work page internal anchor Pith review Pith/arXiv arXiv 2012
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.