Pith. sign in

REVIEW 2 major objections 4 minor 35 references

Combined Belle II data on three B o u decays disfavors a vanishing right-minus-left vector Wilson coefficient at 2σ for every tested |V_ub|.

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

T0 review · grok-4.5

2026-07-13 04:36 UTC pith:V4IWYENY

load-bearing objection Solid dual-method SMEFT fit to the new Belle II binned π and ρ spectra that mildly prefers Re(ε₂^μ) ≠ 0 at 2σ; form-factor systematics remain the main limit. the 2 major comments →

arxiv 2607.09216 v1 pith:V4IWYENY submitted 2026-07-10 hep-ph hep-exhep-lat

Hunting for new physics in B meson b to u transitions: A fresh look with Belle II data

classification hep-ph hep-exhep-lat
keywords b to u transitionsBelle IIWilson coefficientsSMEFT|V_ub|semileptonic B decaysnew physicsform factors
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

This paper asks whether the newest Belle II q^{2}-binned branching fractions for B̄^{0} oπ^{+}μν̄ and B^{-} oρ^{0}μν̄, together with the B^{-} oμν̄ upper limit, already prefer non-standard couplings in the most general dimension-six Hamiltonian with left-handed neutrinos. The three modes are sensitive to different linear combinations of the five complex Wilson coefficients, so their joint analysis can isolate effects that a single channel would miss. Working with three representative values of |V_ub| and then floating |V_ub| itself, the authors map the allowed regions and find that a pure Standard-Model vector coupling is disfavored at the two-sigma level for every choice of |V_ub|. The global fit still returns a |V_ub| close to the exclusive average, showing that new-physics parameters and the CKM element can be extracted together without tension.

Core claim

When the full set of complex Wilson coefficients is floated against Belle II binned spectra plus the leptonic bound, Re(ε_μ^{2})=0 lies outside the 2σ highest-posterior-density region for every benchmark |V_ub|, while a simultaneous fit yields |V_ub|=(3.79±0.17)×10^{-3}.

What carries the argument

The generalized low-energy Hamiltonian (Eq. 1) that encodes five complex, lepton-flavor-dependent Wilson coefficients multiplying the full set of dimension-six operators with left-handed neutrinos; its different linear combinations are isolated by the three complementary decay modes.

Load-bearing premise

The external B oπ and B oρ form-factor parametrizations and their full covariance matrices are taken as free of residual systematic bias large enough to move the preferred Wilson-coefficient regions.

What would settle it

A new lattice or LCSR determination of the B oπ and B oρ form factors whose central values shift by more than the present uncertainties would move the 2σ contours so that Re(ε_μ^{2})=0 re-enters the allowed region for all three |V_ub| benchmarks.

Watch this falsifier — get emailed when new claim-graph text bears on it.

If this is right

  • Future fully differential B oρ(→ππ)μν̄ angular analyses will further tighten the tensor and pseudoscalar bounds already set by the q^{2} spectra.
  • The same SMEFT coefficient set can be re-fitted once Belle II luminosity doubles, testing whether the 2σ preference for non-zero Re(ε_μ^{2}) strengthens or disappears.
  • A simultaneous extraction of |V_ub| and all five Wilson coefficients is already statistically viable and will become the default procedure for exclusive |V_ub| averages.
  • Any ultraviolet model that generates a right-handed vector current for b o uμν must confront the present 2σ exclusion of Re(ε_μ^{2})=0.

Where Pith is reading between the lines

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

  • The same complementarity of leptonic, pion and rho modes can be repeated for the electron and tau channels once comparable binned data appear, testing lepton-flavor universality of the preferred non-zero Re(ε^{2}).
  • If the 2σ preference persists with higher statistics, it would favor ultraviolet models that generate a right-handed charged current at the weak scale rather than purely scalar or tensor new physics.
  • The mild downward shift of |V_ub| when NP is floated suggests that part of the exclusive–inclusive tension may be re-absorbed into Wilson coefficients rather than into form-factor systematics alone.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit.

Referee Report

2 major / 4 minor

Summary. The paper analyzes the decays B⁻ → ℓ⁻ ν̄_ℓ, B̄⁰ → π⁺ ℓ⁻ ν̄_ℓ and B⁻ → ρ⁰ ℓ⁻ ν̄_ℓ in a model-independent SMEFT-consistent low-energy Hamiltonian containing all dimension-six operators with left-handed neutrinos and complex, lepton-flavor-dependent Wilson coefficients. Using Belle II q²-binned differential branching fractions for the semileptonic modes (together with the B⁻ → μ⁻ ν̄_μ upper limit and the B⁻ → τ⁻ ν̄_τ branching fraction), the authors first extract |V_ub| in the SM, then constrain the Wilson-coefficient space for three fixed benchmark values of |V_ub|, and finally perform a global fit that floats |V_ub| simultaneously with all new-physics couplings. Complementary sensitivity of the three channels is exploited, and two independent statistical pipelines (MCMC highest-posterior-density regions and profile-likelihood Δχ² accelerated by a Random-Forest emulator) are shown to yield consistent best-fit points and overlapping confidence regions. The combined analysis mildly disfavors Re(ε_μ^{2}) = 0 at the 2σ level for all benchmarks, while the global fit returns |V_ub| = (3.79 ± 0.17) × 10^{-3}.

Significance. If the mild preference for non-vanishing Re(ε_μ^{2}) survives improved form-factor inputs and larger Belle II samples, the work supplies a concrete, falsifiable target for SMEFT analyses of b → u transitions and demonstrates that exclusive |V_ub| extractions can be performed simultaneously with new-physics Wilson coefficients. The dual statistical construction, the explicit inclusion of theoretical form-factor covariances into the total covariance matrix, and the transparent use of three |V_ub| benchmarks constitute a solid methodological template that can be reused as more precise data become available. The result is therefore of immediate interest to both the flavor-physics and EFT communities, even though the present 2σ tension is not claimed as evidence for new physics.

major comments (2)
  1. Section III C–E and Eqs. (5), (7): the 2σ preference for Re(ε_μ^{2}) ≠ 0 is driven primarily by the B → ρ channel and therefore inherits any residual systematic bias present in the Bharucha et al. (2016) LCSR+lattice form-factor set (and its covariance). While C_th is correctly added, the paper does not quantify how the preferred regions shift under alternative modern form-factor determinations (e.g., more recent lattice-only or combined LCSR+lattice updates). A short robustness check with at least one independent set would strengthen the central claim that the preference is not an artifact of the chosen hadronic inputs.
  2. Section III E and Table I: the posterior for ε_μ^{1} and ε_μ^{2} exhibits a continuous ridge of essentially degenerate maxima (semicircle centered at (−1,0)). The tabulated “best-fit” points therefore change sign with |V_ub| merely because of Monte-Carlo sampling fluctuations. The text should state more explicitly that only the ridge itself (not any single point on it) is robust, and that the quoted coordinates are representative rather than unique global maxima.
minor comments (4)
  1. Page 3 and footnote 3: the decision to average electron and muon samples and treat them as a single muon-like data set is stated but not quantified. A one-sentence estimate of the size of possible lepton-flavor-dependent NP effects that are thereby averaged away would be useful.
  2. Figures 2–4: the 1σ/2σ/3σ shading is clear, yet the location of the SM point (0,0) is not marked on every panel; adding a small cross or open circle would aid the reader.
  3. Appendix B, Table II: the interference terms I_INT,PT vanish identically; a brief remark why this occurs (parity or helicity selection) would improve readability.
  4. References: the arXiv numbers of the Belle II measurement [6] and of the form-factor papers [12,15] should be checked for final journal citations once available.

Circularity Check

0 steps flagged

No significant circularity: free Wilson coefficients and |V_ub| are fitted to external Belle II/PDG data using independent form-factor inputs; self-citations are only methodological.

full rationale

The paper's central results (allowed regions for complex Wilson coefficients ε_μ^{1,2,S,P,T} and the global |V_ub| extraction) are obtained by maximizing a likelihood constructed from Belle II binned branching fractions (external data) plus the PDG B→τν rate, with theoretical predictions built from the standard SMEFT Hamiltonian (Eq. 1) and external form-factor sets (Leljak et al. 2021 for B→π, Bharucha et al. 2016 for B→ρ) whose covariances enter C_th. No equation equates a claimed prediction to a quantity defined by the same fit; the mild 2σ preference for Re(ε_μ^{2})≠0 is simply the outcome of that fit. Self-citations ([11],[14],[27]) supply only the kinematic expressions for the differential widths and helicity amplitudes already derived in prior work by the same group; they do not supply uniqueness theorems, force an ansatz, or close a logical loop that would make the fit results tautological. The dual MCMC/profile-likelihood procedure and the three |V_ub| benchmarks are independent statistical constructions applied to the same external inputs. Hence the derivation chain is self-contained against external benchmarks and exhibits only the ordinary, non-load-bearing self-citation common in sequential phenomenological papers.

Axiom & Free-Parameter Ledger

3 free parameters · 4 axioms · 0 invented entities

The central claims rest on the SMEFT dim-6 left-handed-neutrino Hamiltonian, external form-factor determinations, the Belle II covariance matrices, and flat priors on the Wilson coefficients. No new dynamical entities are postulated; the free parameters are the standard Wilson coefficients and |V_ub| itself.

free parameters (3)
  • Re(ε_μ^{1}), Im(ε_μ^{1}), Re(ε_μ^{2}), Im(ε_μ^{2}), Re(ε_μ^{s}), Im(ε_μ^{s}), Re(ε_μ^{p}), Im(ε_μ^{p}), Re(ε_μᵀ), Im(ε_μ
    Ten real parameters floated in the global MCMC/profile-likelihood fit; best-fit values listed in Table I for three |V_ub| benchmarks.
  • |V_ub|
    Treated both as three fixed benchmarks and as a free parameter in the final global fit, yielding (3.79±0.17)×10^{-3}.
  • BCL and LCSR form-factor coefficients for B→π and B→ρ
    Taken from external references [12,15] with their published covariance matrices; they enter every theoretical bin and dominate C_th.
axioms (4)
  • domain assumption The low-energy Hamiltonian contains only dimension-six operators with left-handed neutrinos (Eq. 1), consistent with SMEFT after integrating out heavy degrees of freedom.
    Stated in Sec. II; right-handed neutrinos and higher-dimensional operators are excluded by construction.
  • domain assumption Isospin symmetry relates B̄⁰→ρ^{+} form factors to the B⁻→ρ⁰ amplitudes used in the fit.
    Invoked in Sec. II C when the Bharucha et al. form factors are adopted.
  • ad hoc to paper Electron and muon samples may be averaged and treated as a single muon-like data set for the purpose of the present analysis.
    Explicitly stated in footnote 3; lepton-flavor universality is assumed at the level of the experimental bins.
  • ad hoc to paper Flat priors on each Wilson coefficient inside (−1,1) and the Metropolis-Hastings proposal widths used for MCMC sampling.
    Sec. III, Eqs. (21) and surrounding text; the prior volume directly shapes the reported HPD regions.

pith-pipeline@v1.1.0-grok45 · 18965 in / 2991 out tokens · 37797 ms · 2026-07-13T04:36:37.908871+00:00 · methodology

0 comments
read the original abstract

We present an analysis of the decays $B^- \to \ell^- {\bar \nu}_\ell$, $\bar B^0 \to \pi^+ \ell^- {\bar \nu}_\ell$, and $B^- \to \rho^0 \ell^- {\bar \nu}_\ell$ leveraging the precise measurements of binned differential branching fractions recently reported by the Belle II Collaboration for the semileptonic modes. We adopt a generalized low-energy Hamiltonian that incorporates all dimension-six operators with left-handed neutrinos -- consistent with the Standard Model effective field theory (SMEFT) -- featuring complex, lepton-flavor-dependent Wilson coefficients. For three representative values of $|V_{ub}|$, we perform a fit to the measured $q^2$-binned observables to constrain the coupling parameter space and determine the corresponding confidence regions. Finally, we perform a global fit for the three decay channels, determining $|V_{ub}|$ together with the new physics couplings.

Figures

Figures reproduced from arXiv: 2607.09216 by Davide Milillo, Fulvia De Fazio, Maria Letizia Di Cuia, Pietro Colangelo.

Figure 1
Figure 1. Figure 1: FIG. 1. Real and imaginary parts of the Wilson coefficient [PITH_FULL_IMAGE:figures/full_fig_p005_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Allowed regions for real and imaginary parts of the NP coefficients [PITH_FULL_IMAGE:figures/full_fig_p007_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. Allowed regions for the couplings [PITH_FULL_IMAGE:figures/full_fig_p008_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. Allowed regions for the couplings [PITH_FULL_IMAGE:figures/full_fig_p009_4.png] view at source ↗

discussion (0)

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

Reference graph

Works this paper leans on

35 extracted references · 21 linked inside Pith

  1. [1]

    In the profile-likelihood analysis we set ∆χ 2 1σ(6) = 7.04

    = 0: the 1σallowed region favors a nonzero value of Re(ϵ µ 1), with the Standard Model point sit- ting outside this region. In the profile-likelihood analysis we set ∆χ 2 1σ(6) = 7.04. The obtained result confirms the same preferred NP solution, while yielding moderately broader con- fidence regions. D.B − →ρ 0µ−¯νµ In this decay mode the involved NP para...

  2. [2]

    This parameter is compatible with zero only at the 2σlevel for|V ub|min

    is favoured at 1σC.L.. This parameter is compatible with zero only at the 2σlevel for|V ub|min. This result is not unexpected: actually, in the analysis performed by the Belle II Collaboration when this channel is separately con- sidered not allowing for NP contributions, the value |Vub|= (3.19±0.12±0.17±0.26)×10 −3 has been ob- tained (with the errors re...

  3. [3]

    4, in which the position of the maximum probability is indicated by dark dots

    This is illustrated in Fig. 4, in which the position of the maximum probability is indicated by dark dots. This also explains the change of sign of the coordinates of some best-fit points in Table I varying|V ub|. The remarkable result of the com- bined analysis is that Re(ϵµ

  4. [4]

    = 0 is disfavoured at 2σ level for all values of|V ub|, and that Re(ϵ µ

  5. [5]

    = 0 is disfavoured at 1σfor|V ub|max. F. Global fit and determination of|V ub| We finally perform a global fit in which also|V ub|is included in the set of parameters, together with the NP couplings. The result |Vub|= (3.79±0.17)×10 −3 ,(26) is slightly smaller than the outcome (24) obtained in the SM. The selected regions of the NP couplings are similar ...

  6. [6]

    for all benchmark values of|V ub|, with Re(ϵµ

  7. [7]

    More- over, for the largest benchmark value of|V ub|, the point Re(ϵ µ

    = 0 lying outside the 2σallowed region. More- over, for the largest benchmark value of|V ub|, the point Re(ϵ µ

  8. [8]

    = 0 is also excluded at the 1σlevel. While these results cannot be interpreted as evidence for physics beyond the Standard Model, they indi- cate regions of the parameter space that deserve fur- ther scrutiny as more precise measurements become available. Finally, we have determined|V ub|allowing simultaneously for generic new-physics contributions. The e...

  9. [9]

    Cabibbo, Phys

    N. Cabibbo, Phys. Rev. Lett.10, 531 (1963)

  10. [10]

    Kobayashi and T

    M. Kobayashi and T. Maskawa, Prog. Theor. Phys. 49, 652 (1973)

  11. [11]

    Gambinoet al., Eur

    P. Gambinoet al., Eur. Phys. J. C80, 966 (2020), arXiv:2006.07287 [hep-ph]

  12. [12]

    Banerjeeet al.(Heavy Flavor Averaging Group (HFLAV)), Phys

    S. Banerjeeet al.(Heavy Flavor Averaging Group (HFLAV)), Phys. Rev. D113, 012008 (2026), arXiv:2411.18639 [hep-ex]

  13. [13]

    Colangelo and F

    P. Colangelo and F. De Fazio, Phys. Rev. D95, 011701 (2017), arXiv:1611.07387 [hep-ph]

  14. [14]

    Adachiet al.(Belle-II), Phys

    I. Adachiet al.(Belle-II), Phys. Rev. D111, 112009 (2025), arXiv:2407.17403 [hep-ex]

  15. [15]

    Beˇ cirevi´ c, M

    D. Beˇ cirevi´ c, M. Martines, S. Rosauro-Alcaraz, and O. Sumensari, Phys. Lett. B880, 140713 (2026), arXiv:2603.25837 [hep-ph]

  16. [16]

    Grzadkowski, M

    B. Grzadkowski, M. Iskrzynski, M. Misiak, and J. Rosiek, JHEP10, 085 (2010), arXiv:1008.4884 [hep-ph]

  17. [17]

    Cornella, M

    C. Cornella, M. Ferr´ e, M. K¨ onig, and M. Neubert, (2026), arXiv:2601.14361 [hep-ph]

  18. [18]

    Aokiet al.(Flavour Lattice Averaging Group (FLAG)), Phys

    Y. Aokiet al.(Flavour Lattice Averaging Group (FLAG)), Phys. Rev. D113, 014508 (2026), arXiv:2411.04268 [hep-lat]

  19. [19]

    Colangelo, F

    P. Colangelo, F. De Fazio, and F. Loparco, Phys. Rev. D100, 075037 (2019), arXiv:1906.07068 [hep- ph]

  20. [20]

    Leljak, B

    D. Leljak, B. Meli´ c, and D. van Dyk, JHEP07, 036 (2021), arXiv:2102.07233 [hep-ph]

  21. [21]

    Bourrely, I

    C. Bourrely, I. Caprini, and L. Lellouch, Phys. Rev. D79, 013008 (2009), [Erratum: Phys.Rev.D 82, 099902 (2010)], arXiv:0807.2722 [hep-ph]

  22. [22]

    Colangelo, F

    P. Colangelo, F. De Fazio, F. Loparco, and N. Losacco, Phys. Rev. D109, 075047 (2024), arXiv:2401.12304 [hep-ph]

  23. [23]

    Bharucha, D

    A. Bharucha, D. M. Straub, and R. Zwicky, JHEP 08, 098 (2016), arXiv:1503.05534 [hep-ph]

  24. [24]

    Haet al.(Belle), Phys

    H. Haet al.(Belle), Phys. Rev. D83, 071101 (2011), arXiv:1012.0090 [hep-ex]

  25. [25]

    del Amo Sanchezet al.(BaBar), Phys

    P. del Amo Sanchezet al.(BaBar), Phys. Rev. D83, 032007 (2011), arXiv:1005.3288 [hep-ex]

  26. [26]

    J. P. Leeset al.(BaBar), Phys. Rev. D86, 092004 (2012), arXiv:1208.1253 [hep-ex]

  27. [27]

    Sibidanovet al.(Belle), Phys

    A. Sibidanovet al.(Belle), Phys. Rev. D88, 032005 (2013), arXiv:1306.2781 [hep-ex]

  28. [28]

    B. H. Behrenset al.(CLEO), Phys. Rev. D61, 052001 (2000), arXiv:hep-ex/9905056

  29. [29]

    N. E. Adamet al.(CLEO), Phys. Rev. Lett.99, 041802 (2007), arXiv:hep-ex/0703041

  30. [30]

    Hokuueet al.(Belle), Phys

    T. Hokuueet al.(Belle), Phys. Lett. B648, 139 (2007), arXiv:hep-ex/0604024

  31. [31]

    Trotta, Contemporary Physics49, 71–104 (2008)

    R. Trotta, Contemporary Physics49, 71–104 (2008)

  32. [32]

    Navaset al.(Particle Data Group), Phys

    S. Navaset al.(Particle Data Group), Phys. Rev. D 110, 030001 (2024)

  33. [33]

    Roberts and J

    G. Roberts and J. Rosenthal, Statistical Science16 (2001), 10.1214/ss/1015346320

  34. [34]

    Wolfram Language Documentation,

    SmoothKernelDistribution, “Wolfram Language Documentation,” [reference.wolfram.com]

  35. [35]

    Colangelo and F

    P. Colangelo and F. De Fazio, JHEP06, 082 (2018), arXiv:1801.10468 [hep-ph]. 12 |Vub| ×10 3 lnLmax Re(ϵµ 1 ) Im(ϵ µ 1 ) Re(ϵ µ 2 ) Im(ϵ µ 2 ) Re(ϵ µ S) Im(ϵ µ S) Re(ϵ µ P ) Im(ϵ µ P ) Re(ϵ µ T ) Im(ϵ µ T ) 3.51 288.790 -0.190 0.771 -0.856 -0.715 -0.057 0.020 -0.046 0.031 -0.090 -0.082 3.85 288.789 -0.582 -0.937 -0.884 0.705 0.017 0.072 -0.036 -0.026 -0.00...