pith. sign in

arxiv: 2510.18953 · v2 · submitted 2025-10-21 · ✦ hep-ph · hep-ex

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

classification ✦ hep-ph hep-ex
keywords Belle IIB meson decaysQCD axionrare decaysmodel-independent analysisinvisible particlesflavour-changing couplingsaxion-like particles
0
0 comments X

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.

The paper develops a framework that reinterprets existing measurements of rare B meson decays using only public data to search for light invisible particles such as the QCD axion. It does so by analytically reconstructing the mapping between the true kinematic variables and the reconstructed ones used in the statistically dominant Inclusive Tagging Analysis. Applying this to the B+ to K+ nu nubar decay produces an improved upper bound on the branching fraction for B+ to K+ a. The approach shows that this single measurement can test both short-distance new physics and light invisible states to a good approximation.

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

These are editorial extensions of the paper, not claims the author makes directly.

  • 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

Figures reproduced from arXiv: 2510.18953 by Axel Iohner, Claudio Toni, Diego Guadagnoli, Giulio Dujany, Merna Abumusabh.

Figure 1
Figure 1. Figure 1: FIG. 1: Correlation between the true ( [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 5
Figure 5. Figure 5: figure 5.1 of Ref. [ [PITH_FULL_IMAGE:figures/full_fig_p003_5.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2: Upper limit on [PITH_FULL_IMAGE:figures/full_fig_p005_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3: Twice the negative profile log-likelihood ratio as a function of [PITH_FULL_IMAGE:figures/full_fig_p007_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4: Regions [PITH_FULL_IMAGE:figures/full_fig_p008_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5: Upper limit at 90% CL on [PITH_FULL_IMAGE:figures/full_fig_p008_5.png] view at source ↗
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.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

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)
  1. [§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).
  2. [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)
  1. Clarify the definition of the fundamental flavour-changing coupling in the context of the axion model.
  2. [Abstract] The abstract could specify the mass range for the axion considered in the bound.

Simulated Author's Rebuttal

2 responses · 0 unresolved

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
  1. 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

  2. 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

0 steps flagged

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

0 free parameters · 1 axioms · 0 invented entities

The central claim rests on the accuracy of an analytic kinematic mapping derived from public information and on standard assumptions about the dominance of the Inclusive Tagging Analysis; no new particles or forces are postulated.

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.
    Invoked when the paper states that the reconstruction enables rare-decay measurements to probe light invisible particles using only public data.

pith-pipeline@v0.9.0 · 5729 in / 1394 out tokens · 36478 ms · 2026-05-18T04:36:55.477792+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

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

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

  1. Braking protons at the EIC: from invisible meson decay to new physics searches

    hep-ph 2025-12 unverdicted novelty 6.0

    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.

  2. Dark Matter emission at Belle II and NA62 in Minimal Flavor Violation framework

    hep-ph 2026-01 unverdicted novelty 5.0

    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

50 extracted references · 50 canonical work pages · cited by 2 Pith papers · 11 internal anchors

  1. [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]

  2. [2]

    Lanfranchi, M

    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. [3]

    Invisible Axion Search Methods,

    P. Sikivie, Invisible Axion Search Methods, Rev. Mod. Phys. 93, 015004 (2021), arXiv:2003.02206 [hep-ph]

  4. [4]

    G. G. Raffelt, Astrophysical methods to constrain axions and other novel particle phenomena, Phys. Rept.198, 1 (1990)

  5. [5]

    G. G. Raffelt, Astrophysical axion bounds, Lect. Notes Phys. 741, 51 (2008), arXiv:hep-ph/0611350

  6. [6]

    Caputo and G

    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. [7]

    D. J. E. Marsh, Axion Cosmology, Phys. Rept.643, 1 (2016), arXiv:1510.07633 [astro-ph.CO]

  8. [8]

    J. E. Kim and G. Carosi, Axions and the Strong CP Problem, Rev. Mod. Phys. 82, 557 (2010), [Erratum: Rev.Mod.Phys.91, 049902(2019)],arXiv:0807.3125[hep- ph]

  9. [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]

  10. [10]

    K. Choi, S. H. Im, and C. Sub Shin, Recent Progress in the Physics of Axions and Axion-Like Particles, Ann. Rev. Nucl. Part. Sci. 71, 225 (2021), arXiv:2012.05029 [hep-ph]

  11. [11]

    J. E. Kim, Weak Interaction Singlet and Strong CP In- variance, Phys. Rev. Lett.43, 103 (1979)

  12. [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)

  13. [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)

  14. [14]

    A. R. Zhitnitsky, On Possible Suppression of the Axion Hadron Interactions. (In Russian), Sov. J. Nucl. Phys. 31, 260 (1980)

  15. [15]

    R. D. Peccei and H. R. Quinn, CP Conservation in the Presence of Instantons, Phys. Rev. Lett.38, 1440 (1977)

  16. [16]

    R. D. Peccei and H. R. Quinn, Constraints Imposed by CP Conservation in the Presence of Instantons, Phys. Rev. D 16, 1791 (1977)

  17. [17]

    Weinberg, A New Light Boson?, Phys

    S. Weinberg, A New Light Boson?, Phys. Rev. Lett.40, 223 (1978)

  18. [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

  19. [19]

    L.Ubaldi,Effectsofthetaonthedeuteronbindingenergy and the triple-alpha process, Phys. Rev. D81, 025011 (2010), arXiv:0811.1599 [hep-ph]

  20. [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]

  21. [21]

    Preskill, M

    J. Preskill, M. B. Wise, and F. Wilczek, Cosmology of the Invisible Axion, Phys. Lett. B120, 127 (1983)

  22. [22]

    L. F. Abbott and P. Sikivie, A Cosmological Bound on the Invisible Axion, Phys. Lett. B120, 133 (1983)

  23. [23]

    Dine and W

    M. Dine and W. Fischler, The Not So Harmless Axion, Phys. Lett. B120, 137 (1983)

  24. [24]

    Belle II Technical Design Report

    T. Abeet al. (Belle-II), Belle II Technical Design Report, (2010), arXiv:1011.0352 [physics.ins-det]

  25. [25]

    The spread, or bite, around these momenta is tiny and completely negligible for the purposes of this paper

  26. [26]

    The conjugated mode with negatively-charged mesons is assumed to be included throughout our study

  27. [27]

    Georgi, D

    H. Georgi, D. B. Kaplan, and L. Randall, Manifesting the Invisible Axion at Low-energies, Phys. Lett. B169, 73 (1986)

  28. [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. [29]

    Gubernari, M

    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. [30]

    Navaset al

    S. Navaset al. (Particle Data Group), Review of particle physics, Phys. Rev. D110, 030001 (2024)

  31. [31]

    P. D. Bolton, S. Fajfer, J. F. Kamenik, and M. Novoa-Brunet, Signatures of light new particles in B→K(*)Emiss, Phys. Rev. D110, 055001 (2024), [Erra- tum: Phys.Rev.D 111, 039903 (2025)], arXiv:2403.13887 [hep-ph]

  32. [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. [33]

    M. T. Prim et al. (Belle), Search forB+→µ+νµand B+→µ+N with inclusive tagging, Phys. Rev. D101, 032007 (2020), arXiv:1911.03186 [hep-ex]

  34. [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. [35]

    W. G. Parrott, C. Bouchard, and C. T. H. Davies (HPQCD), Standard Model predictions for B→Kℓ+ℓ- , B→Kℓ1-ℓ2+ and B→Kνν¯ using form factors from Nf=2+1+1 lattice QCD, Phys. Rev. D 107, 014511 (2023), [Erratum: Phys.Rev.D 107, 119903 (2023)], arXiv:2207.13371 [hep-ph]

  36. [36]

    Beˇ cirevi´ c, G

    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. [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]

  38. [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]

  39. [39]

    Bertacchiet al

    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. [40]

    The below-2%NB error has a negligible impact on our analysis, becauseni ν¯ν,bkg have an error approaching 10%

  41. [41]

    Fridell, M

    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. [42]

    P. D. Bolton, S. Fajfer, J. F. Kamenik, and M. Novoa- Brunet, Impact of new invisible particles on B → K(∗)Emiss observables, (2025), arXiv:2503.19025 [hep- ph]

  43. [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)

  44. [44]

    theoretical

    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. [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]

  46. [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. [47]

    Altmannshofer, A

    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. [48]

    Cavan-Piton, D

    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. [49]

    Martin Camalich, M

    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. [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 ...