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 →
Hunting for new physics in B meson b to u transitions: A fresh look with Belle II data
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
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.
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
- 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.
Referee Report
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)
- 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.
- 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)
- 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.
- 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.
- 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.
- 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
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
free parameters (3)
- Re(ε_μ^{1}), Im(ε_μ^{1}), Re(ε_μ^{2}), Im(ε_μ^{2}), Re(ε_μ^{s}), Im(ε_μ^{s}), Re(ε_μ^{p}), Im(ε_μ^{p}), Re(ε_μᵀ), Im(ε_μ
- |V_ub|
- BCL and LCSR form-factor coefficients for B→π and B→ρ
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.
- domain assumption Isospin symmetry relates B̄⁰→ρ^{+} form factors to the B⁻→ρ⁰ amplitudes used in the fit.
- 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.
- ad hoc to paper Flat priors on each Wilson coefficient inside (−1,1) and the Metropolis-Hastings proposal widths used for MCMC sampling.
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
Reference graph
Works this paper leans on
-
[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]
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]
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]
= 0 is disfavoured at 2σ level for all values of|V ub|, and that Re(ϵ µ
-
[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]
for all benchmark values of|V ub|, with Re(ϵµ
-
[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]
= 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]
Cabibbo, Phys
N. Cabibbo, Phys. Rev. Lett.10, 531 (1963)
1963
-
[10]
Kobayashi and T
M. Kobayashi and T. Maskawa, Prog. Theor. Phys. 49, 652 (1973)
1973
-
[11]
P. Gambinoet al., Eur. Phys. J. C80, 966 (2020), arXiv:2006.07287 [hep-ph]
Pith/arXiv arXiv 2020
-
[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]
Pith/arXiv arXiv 2026
-
[13]
P. Colangelo and F. De Fazio, Phys. Rev. D95, 011701 (2017), arXiv:1611.07387 [hep-ph]
Pith/arXiv arXiv 2017
-
[14]
I. Adachiet al.(Belle-II), Phys. Rev. D111, 112009 (2025), arXiv:2407.17403 [hep-ex]
Pith/arXiv arXiv 2025
-
[15]
D. Beˇ cirevi´ c, M. Martines, S. Rosauro-Alcaraz, and O. Sumensari, Phys. Lett. B880, 140713 (2026), arXiv:2603.25837 [hep-ph]
Pith/arXiv arXiv 2026
-
[16]
B. Grzadkowski, M. Iskrzynski, M. Misiak, and J. Rosiek, JHEP10, 085 (2010), arXiv:1008.4884 [hep-ph]
Pith/arXiv arXiv 2010
-
[17]
C. Cornella, M. Ferr´ e, M. K¨ onig, and M. Neubert, (2026), arXiv:2601.14361 [hep-ph]
Pith/arXiv arXiv 2026
-
[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]
Pith/arXiv arXiv 2026
-
[19]
P. Colangelo, F. De Fazio, and F. Loparco, Phys. Rev. D100, 075037 (2019), arXiv:1906.07068 [hep- ph]
Pith/arXiv arXiv 2019
-
[20]
D. Leljak, B. Meli´ c, and D. van Dyk, JHEP07, 036 (2021), arXiv:2102.07233 [hep-ph]
Pith/arXiv arXiv 2021
-
[21]
C. Bourrely, I. Caprini, and L. Lellouch, Phys. Rev. D79, 013008 (2009), [Erratum: Phys.Rev.D 82, 099902 (2010)], arXiv:0807.2722 [hep-ph]
Pith/arXiv arXiv 2009
-
[22]
P. Colangelo, F. De Fazio, F. Loparco, and N. Losacco, Phys. Rev. D109, 075047 (2024), arXiv:2401.12304 [hep-ph]
Pith/arXiv arXiv 2024
-
[23]
A. Bharucha, D. M. Straub, and R. Zwicky, JHEP 08, 098 (2016), arXiv:1503.05534 [hep-ph]
Pith/arXiv arXiv 2016
-
[24]
H. Haet al.(Belle), Phys. Rev. D83, 071101 (2011), arXiv:1012.0090 [hep-ex]
Pith/arXiv arXiv 2011
-
[25]
del Amo Sanchezet al.(BaBar), Phys
P. del Amo Sanchezet al.(BaBar), Phys. Rev. D83, 032007 (2011), arXiv:1005.3288 [hep-ex]
Pith/arXiv arXiv 2011
-
[26]
J. P. Leeset al.(BaBar), Phys. Rev. D86, 092004 (2012), arXiv:1208.1253 [hep-ex]
Pith/arXiv arXiv 2012
-
[27]
A. Sibidanovet al.(Belle), Phys. Rev. D88, 032005 (2013), arXiv:1306.2781 [hep-ex]
Pith/arXiv arXiv 2013
-
[28]
B. H. Behrenset al.(CLEO), Phys. Rev. D61, 052001 (2000), arXiv:hep-ex/9905056
Pith/arXiv arXiv 2000
-
[29]
N. E. Adamet al.(CLEO), Phys. Rev. Lett.99, 041802 (2007), arXiv:hep-ex/0703041
Pith/arXiv arXiv 2007
-
[30]
T. Hokuueet al.(Belle), Phys. Lett. B648, 139 (2007), arXiv:hep-ex/0604024
Pith/arXiv arXiv 2007
-
[31]
Trotta, Contemporary Physics49, 71–104 (2008)
R. Trotta, Contemporary Physics49, 71–104 (2008)
2008
-
[32]
Navaset al.(Particle Data Group), Phys
S. Navaset al.(Particle Data Group), Phys. Rev. D 110, 030001 (2024)
2024
-
[33]
G. Roberts and J. Rosenthal, Statistical Science16 (2001), 10.1214/ss/1015346320
-
[34]
Wolfram Language Documentation,
SmoothKernelDistribution, “Wolfram Language Documentation,” [reference.wolfram.com]
-
[35]
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...
Pith/arXiv arXiv 2018
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.