Inclusive bar B_smapsto X_(bar sc) ell bar ν decays from lattice QCD: computational strategy and a first physical result
Pith reviewed 2026-07-02 01:25 UTC · model grok-4.3
The pith
The inclusive decay rate for anti-B_s to anti-charm lepton antineutrino is obtained by interpolating lattice QCD results at masses up to 4.3 GeV with operator product expansion predictions at infinite quark mass.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The physical decay rate is obtained from the interpolation of non-perturbative lattice data, obtained at lighter than physical heavy meson masses (M_Bs^max=4.3 GeV), with the Operator Product Expansion predictions, which become exact in the limit of infinitely heavy quarks. A new method for the computation of the required lattice four-point correlators is presented, and the strategy is shown to work on a subset of ETMC gauge ensembles, yielding a 7 percent result.
What carries the argument
Interpolation of lattice four-point correlators computed at finite heavy-meson masses up to 4.3 GeV onto operator product expansion results that hold exactly at infinite quark mass.
If this is right
- The inclusive rate can be computed from first principles once more ensembles and heavier masses become available.
- The new four-point correlator method lowers the cost of similar inclusive decay calculations.
- The 7 percent error is expected to drop substantially with increased statistics and fuller ensemble coverage.
- The same interpolation framework applies to other inclusive heavy-meson semileptonic processes.
Where Pith is reading between the lines
- If the interpolation holds, the result supplies a benchmark for testing the convergence of the heavy-quark expansion at physical masses.
- Direct lattice access to the physical B_s point would require overcoming the severe signal-to-noise degradation at heavier masses.
- Comparison of the interpolated rate with experimental data would test whether higher-order OPE terms remain small in the accessible mass window.
Load-bearing premise
Interpolating the lattice results between 4.3 GeV and the infinite-mass limit reproduces the physical decay rate without large uncontrolled higher-order corrections or bias from the limited mass range and ensemble size.
What would settle it
A direct lattice calculation performed at or very near the physical B_s mass that differs from the interpolated value by more than the quoted 7 percent uncertainty.
Figures
read the original abstract
We present a strategy to compute the inclusive decay rate for the process $\bar B_s \mapsto X_{\bar sc} \ell \bar{\nu}$ from first principles in lattice QCD. The physical decay rate is obtained from the interpolation of non-perturbative lattice data, obtained at lighter than physical heavy meson masses ($M_{\bar B_s}^\mathrm{max}=4.3$GeV), with the Operator Product Expansion predictions, which become exact in the limit of infinitely heavy quarks. We also present a new method for the computation of the required lattice four-point correlators, which represents a considerable improvement over the state-of-the-art on the subject. We show the effectiveness of the strategy by performing the calculation on a subset of the available $n_f=2+1+1$ physical-point Extended Twisted Mass Collaboration (ETMC) gauge ensembles. Our current determination of the inclusive decay rate has a 7% total error, that is dominated by uncertainties due to the relatively limited configuration ensembles considered herein, and can be significantly reduced in the near future.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript presents a strategy to compute the inclusive decay rate for ar B_s o X_{ar s c} ar u from lattice QCD. Non-perturbative lattice data at heavy-meson masses up to M_{ar B_s}^max = 4.3 GeV are combined via interpolation with Operator Product Expansion predictions that become exact at infinite mass; a new method for four-point correlators is introduced and demonstrated on a subset of ETMC n_f=2+1+1 ensembles, yielding a first physical-point result with a quoted 7% total uncertainty dominated by statistics.
Significance. If the interpolation procedure is shown to be robust, the work supplies a first-principles determination of an inclusive semileptonic rate that can be compared with exclusive channels and used in |V_cb| analyses. The new correlator technique is a concrete technical advance that could reduce computational cost for similar observables. The current result remains preliminary because of the limited ensemble set and the mass range covered.
major comments (2)
- [Abstract] Abstract: the physical result is obtained by interpolating lattice data at M ≤ 4.3 GeV with the exact OPE value at M o∞ and evaluating at the physical mass ≈5.37 GeV, yet no functional form, fit quality, or stability tests against alternative ansätze are reported. Because the highest simulated mass lies ~20 % below the physical point, any mismatch between the true 1/M expansion and the chosen interpolant directly affects the central value at a level comparable to the quoted 7 % error.
- [Abstract] Abstract: the total uncertainty is stated to be 7 % and statistics-dominated, but no quantitative breakdown or validation is given for the systematic component arising from the interpolation itself, nor for discretization or finite-volume effects on the four-point functions. This omission leaves the reliability of the central claim only partially supported.
minor comments (1)
- The abstract refers to “a subset of the available n_f=2+1+1 physical-point ETMC gauge ensembles” without listing the specific ensembles, their volumes, or lattice spacings; adding these parameters would improve reproducibility.
Simulated Author's Rebuttal
We thank the referee for the careful reading of our manuscript and for the constructive comments. We appreciate the positive assessment of the significance of the computational strategy and the new four-point correlator technique. We address the major comments point by point below.
read point-by-point responses
-
Referee: [Abstract] Abstract: the physical result is obtained by interpolating lattice data at M ≤ 4.3 GeV with the exact OPE value at M→∞ and evaluating at the physical mass ≈5.37 GeV, yet no functional form, fit quality, or stability tests against alternative ansätze are reported. Because the highest simulated mass lies ~20 % below the physical point, any mismatch between the true 1/M expansion and the chosen interpolant directly affects the central value at a level comparable to the quoted 7 % error.
Authors: We agree that the abstract would be improved by including a concise description of the interpolation procedure. The full manuscript employs a polynomial ansatz in inverse heavy-meson mass, constrained by the exact OPE result at infinite mass and motivated by the structure of the 1/M expansion. Fit quality is monitored through χ²/dof, and stability under changes in polynomial degree and included mass points is verified. To address the referee’s concern directly, we will revise the abstract to state the functional form used and note the outcome of the stability tests. The ~20 % mass gap is bridged by anchoring to the exact OPE limit; the associated systematic uncertainty from ansatz choice is estimated to be sub-dominant to statistics in our current error budget, but we will add an explicit quantification of this component in the revised version. revision: yes
-
Referee: [Abstract] Abstract: the total uncertainty is stated to be 7 % and statistics-dominated, but no quantitative breakdown or validation is given for the systematic component arising from the interpolation itself, nor for discretization or finite-volume effects on the four-point functions. This omission leaves the reliability of the central claim only partially supported.
Authors: We concur that a more detailed error breakdown would strengthen the presentation. The quoted 7 % total uncertainty is indeed statistics-dominated because of the limited number of gauge configurations employed in this first demonstration. The interpolation systematic is assessed by varying the polynomial degree and the mass range included in the fit; discretization effects are suppressed by the use of maximally twisted mass fermions, and finite-volume effects are expected to be small on the large ETMC volumes, but neither has been quantified numerically in the present work. We will add a dedicated error-budget table in the revised manuscript that separates the interpolation systematic from the statistical error and provides preliminary estimates for discretization and finite-volume contributions. A fuller quantification of all systematics will become possible once the full set of ensembles is analyzed. revision: yes
Circularity Check
No significant circularity: interpolation uses external OPE anchor independent of lattice inputs
full rationale
The derivation obtains the physical inclusive rate by interpolating non-perturbative lattice four-point correlators computed on ETMC ensembles at M_Bs ≤ 4.3 GeV to the physical point, using the standard OPE which is exact only at infinite heavy-quark mass. The OPE limit is an independent perturbative expansion whose validity does not depend on the present lattice data or any fitted parameter from this work. No equation in the abstract or described strategy reduces the target rate to a fit of itself, nor does any load-bearing step rely on a self-citation whose content is unverified or defined circularly. The 7 % error is dominated by ensemble statistics, not by a tautological construction. The method is therefore self-contained against external benchmarks.
Axiom & Free-Parameter Ledger
axioms (1)
- domain assumption Operator Product Expansion predictions become exact in the limit of infinitely heavy quarks
Reference graph
Works this paper leans on
-
[1]
The left, center and right panels correspond to the analysis orders CSI, SCI and SIC, respectively. The green band cor- responds to our final determination, ˆΓ(1) = 0.187(10)(8)[13], whose central value is obtained from the average of the six points and the error is the combination in quadrature of the statistical and systematic uncertainties. The systema...
-
[2]
We acknowledge EuroHPC Joint Undertaking for awarding the project ID EHPC-EXT-2023E02-052 access to MareNostrum5 hosted by the Barcelona Supercomput- ing Center, Spain
We acknowledge the Swiss National Supercomput- ing Centre (CSCS) and the EuroHPC Joint Undertaking for awarding this project access to the LUMI supercom- puter, owned by the EuroHPC Joint Undertaking, hosted by CSC (Finland) and the LUMI consortium through the Chronos programme under project IDs CH17-CSCS- CYP. We acknowledge EuroHPC Joint Undertaking for...
2022
-
[3]
A. J. Buras, Hunting New Animalcula with Flavour Changing Processes, Acta Phys. Polon. B57, 6 (2026), arXiv:2605.27525 [hep-ph]
work page internal anchor Pith review Pith/arXiv arXiv 2026
-
[4]
Gambinoet al., Challenges in semileptonicBdecays, Eur
P. Gambinoet al., Challenges in semileptonicBdecays, Eur. Phys. J. C80, 966 (2020), arXiv:2006.07287 [hep- ph]
-
[5]
Averages of $b$-hadron, $c$-hadron, and $\tau$-lepton properties as of 2023
S. Banerjeeet al.(Heavy Flavor Averaging Group (HFLAV)), Averages ofb-hadron,c-hadron, andτ-lepton properties as of 2023, Phys. Rev. D113, 012008 (2026), arXiv:2411.18639 [hep-ex]
work page internal anchor Pith review Pith/arXiv arXiv 2023
-
[6]
J. C. A. Barata and K. Fredenhagen, Particle scattering in Euclidean lattice field theories, Commun. Math. Phys. 138, 507 (1991)
1991
-
[7]
A. Patella and N. Tantalo, Scattering amplitudes from Euclidean correlators: Haag-Ruelle theory and approxi- mation formulae, JHEP01, 091, arXiv:2407.02069 [hep- lat]
-
[8]
M. T. Hansen, H. B. Meyer, and D. Robaina, From deep inelastic scattering to heavy-flavor semileptonic decays: Total rates into multihadron final states from lattice QCD, Phys. Rev. D96, 094513 (2017), arXiv:1704.08993 [hep-lat]
work page internal anchor Pith review Pith/arXiv arXiv 2017
-
[9]
S. Hashimoto, Inclusive semi-leptonicBmeson decay structure functions from lattice QCD, PTEP2017, 053B03 (2017), arXiv:1703.01881 [hep-lat]
-
[10]
On the extraction of spectral densities from lattice correlators
M. Hansen, A. Lupo, and N. Tantalo, Extraction of spec- tral densities from lattice correlators, Phys. Rev. D99, 094508 (2019), arXiv:1903.06476 [hep-lat]
work page internal anchor Pith review Pith/arXiv arXiv 2019
-
[11]
P. Gambino and S. Hashimoto, Inclusive Semileptonic Decays from Lattice QCD, Phys. Rev. Lett.125, 032001 (2020), arXiv:2005.13730 [hep-lat]
-
[12]
P. Gambino, S. Hashimoto, S. M¨ achler, M. Panero, F. Sanfilippo, S. Simula, A. Smecca, and N. Tantalo, Lat- tice QCD study of inclusive semileptonic decays of heavy mesons, JHEP07, 083, arXiv:2203.11762 [hep-lat]
-
[13]
A. De Santiset al., Inclusive Semileptonic Decays of the Ds Meson: Lattice QCD Confronts Experiments, Phys. Rev. Lett.135, 121901 (2025), arXiv:2504.06064 [hep- lat]
-
[14]
A. De Santiset al., Inclusive semileptonic decays of the Ds meson: A first-principles lattice QCD calculation, Phys. Rev. D112, 054503 (2025), arXiv:2504.06063 [hep- lat]
-
[15]
R. Kellermann, A. Barone, A. Elgaziari, S. Hashimoto, Z. Hu, A. J¨ uttner, and T. Kaneko, Inclusive semileptonic Ds →X sℓ¯νdecays from lattice QCD: continuum and chiral extrapolation, (2026), arXiv:2604.22201 [hep-lat]
work page internal anchor Pith review Pith/arXiv arXiv 2026
-
[16]
Y. Aokiet al.(Flavour Lattice Averaging Group (FLAG)), FLAG review 2024, Phys. Rev. D113, 014508 (2026), arXiv:2411.04268 [hep-lat]
work page internal anchor Pith review Pith/arXiv arXiv 2024
- [17]
-
[18]
A. X. El-Khadra, A. S. Kronfeld, and P. B. Mackenzie, Massive fermions in lattice gauge theory, Phys. Rev. D 55, 3933 (1997), arXiv:hep-lat/9604004
work page internal anchor Pith review Pith/arXiv arXiv 1997
-
[19]
N. H. Christ, M. Li, and H.-W. Lin, Relativistic Heavy Quark Effective Action, Phys. Rev. D76, 074505 (2007), arXiv:hep-lat/0608006
work page internal anchor Pith review Pith/arXiv arXiv 2007
-
[20]
Non-perturbatively Determined Relativistic Heavy Quark Action
H.-W. Lin and N. Christ, Non-perturbatively Deter- mined Relativistic Heavy Quark Action, Phys. Rev. D 76, 074506 (2007), arXiv:hep-lat/0608005
work page internal anchor Pith review Pith/arXiv arXiv 2007
-
[21]
A. V. Manohar and M. B. Wise, Inclusive semileptonic Band polarized Λ b decays from QCD, Phys. Rev. D49, 1310 (1994), arXiv:hep-ph/9308246. 29
work page internal anchor Pith review Pith/arXiv arXiv 1994
- [22]
-
[23]
I. I. Y. Bigi, N. G. Uraltsev, and A. I. Vainshtein, Non- perturbative corrections to inclusive beauty and charm decays: QCD versus phenomenological models, Phys. Lett. B293, 430 (1992), [Erratum: Phys.Lett.B 297, 477–477 (1992)], arXiv:hep-ph/9207214
work page internal anchor Pith review Pith/arXiv arXiv 1992
-
[24]
I. I. Y. Bigi, M. A. Shifman, N. G. Uraltsev, and A. I. Vainshtein, QCD predictions for lepton spectra in inclu- sive heavy flavor decays, Phys. Rev. Lett.71, 496 (1993), arXiv:hep-ph/9304225
work page internal anchor Pith review Pith/arXiv arXiv 1993
-
[25]
J. Chay, H. Georgi, and B. Grinstein, Lepton energy dis- tributions in heavy meson decays from QCD, Phys. Lett. B247, 399 (1990)
1990
- [26]
-
[27]
C. Alexandrouet al.(Extended Twisted Mass Collabo- ration (ETMC)), Probing the Energy-SmearedRRatio Using Lattice QCD, Phys. Rev. Lett.130, 241901 (2023), arXiv:2212.08467 [hep-lat]
-
[28]
A. Evangelista, R. Frezzotti, N. Tantalo, G. Gagliardi, F. Sanfilippo, S. Simula, and V. Lubicz (Extended Twisted Mass), Inclusive hadronic decay rate of theτlep- ton from lattice QCD, Phys. Rev. D108, 074513 (2023), arXiv:2308.03125 [hep-lat]
-
[29]
C. Bonanno, F. D’Angelo, M. D’Elia, L. Maio, and M. Naviglio, Sphaleron rate from a modified Backus- Gilbert inversion method, Phys. Rev. D108, 074515 (2023), arXiv:2305.17120 [hep-lat]
-
[30]
C. Bonanno, F. D’Angelo, M. D’Elia, L. Maio, and M. Naviglio, Sphaleron Rate ofN f = 2 + 1 QCD, Phys. Rev. Lett.132, 051903 (2024), arXiv:2308.01287 [hep- lat]
-
[31]
C. Alexandrouet al.(Extended Twisted Mass), Inclusive Hadronic Decay Rate of theτLepton from Lattice QCD: The ¯usFlavor Channel and the Cabibbo Angle, Phys. Rev. Lett.132, 261901 (2024), arXiv:2403.05404 [hep- lat]
- [32]
-
[33]
Sirlin, Largem W ,m Z Behavior of theO(α) Correc- tions to Semileptonic Processes Mediated byW, Nucl
A. Sirlin, Largem W ,m Z Behavior of theO(α) Correc- tions to Semileptonic Processes Mediated byW, Nucl. Phys. B196, 83 (1982)
1982
-
[34]
M. Neubert, Heavy quark symmetry, Phys. Rept.245, 259 (1994), arXiv:hep-ph/9306320
work page internal anchor Pith review Pith/arXiv arXiv 1994
- [35]
-
[36]
Nir, The Mass Ratiom c/mb in SemileptonicBDe- cays, Phys
Y. Nir, The Mass Ratiom c/mb in SemileptonicBDe- cays, Phys. Lett. B221, 184 (1989)
1989
-
[37]
Heavy-to-heavy quark decays at NNLO
A. Pak and A. Czarnecki, Heavy-to-heavy quark de- cays at NNLO, Phys. Rev. D78, 114015 (2008), arXiv:0808.3509 [hep-ph]
work page internal anchor Pith review Pith/arXiv arXiv 2008
-
[38]
Semileptonic decays in the limit of a heavy daughter quark
M. Dowling, J. H. Piclum, and A. Czarnecki, Semilep- tonic decays in the limit of a heavy daughter quark, Phys. Rev. D78, 074024 (2008), arXiv:0810.0543 [hep-ph]
work page internal anchor Pith review Pith/arXiv arXiv 2008
-
[39]
M. Beneke, Renormalons, Phys. Rept.317, 1 (1999), arXiv:hep-ph/9807443
work page internal anchor Pith review Pith/arXiv arXiv 1999
-
[40]
I. I. Y. Bigi, M. A. Shifman, N. Uraltsev, and A. I. Vainshtein, High powernofm b inb-flavored widths andn= 5→ ∞limit, Phys. Rev. D56, 4017 (1997), arXiv:hep-ph/9704245
work page internal anchor Pith review Pith/arXiv arXiv 1997
- [41]
-
[43]
Simulating twisted mass fermions at physical light, strange and charm quark masses
C. Alexandrouet al., Simulating twisted mass fermions at physical light, strange and charm quark masses, Phys. Rev. D98, 054518 (2018), arXiv:1807.00495 [hep-lat]
work page internal anchor Pith review Pith/arXiv arXiv 2018
-
[44]
G. Bergner, P. Dimopoulos, J. Finkenrath, E. Fiorenza, R. Frezzotti, M. Garofalo, B. Kostrzewa, F. Sanfilippo, S. Simula, and U. Wenger (Extended Twisted Mass), Quark masses and decay constants inN f = 2 + 1 + 1 isoQCD with Wilson clover twisted mass fermions, PoS LATTICE2019, 181 (2020), arXiv:2001.09116 [hep-lat]
-
[45]
C. Alexandrouet al.(Extended Twisted Mass), Ratio of kaon and pion leptonic decay constants withN f = 2+1+1 Wilson-clover twisted-mass fermions, Phys. Rev. D104, 074520 (2021), arXiv:2104.06747 [hep-lat]
-
[46]
J. Finkenrathet al., Twisted mass gauge ensem- bles at physical values of the light, strange and charm quark masses, PoSLATTICE2021, 284 (2022), arXiv:2201.02551 [hep-lat]
-
[47]
Lattice QCD with a chirally twisted mass term
R. Frezzotti, P. A. Grassi, S. Sint, and P. Weisz (Alpha), Lattice QCD with a chirally twisted mass term, JHEP 08, 058, arXiv:hep-lat/0101001
work page internal anchor Pith review Pith/arXiv arXiv
-
[48]
Twisted-mass lattice QCD with mass non-degenerate quarks
R. Frezzotti and G. C. Rossi, Twisted mass lattice QCD with mass nondegenerate quarks, Nucl. Phys. B Proc. Suppl.128, 193 (2004), arXiv:hep-lat/0311008
work page internal anchor Pith review Pith/arXiv arXiv 2004
-
[49]
Chirally improving Wilson fermions II. Four-quark operators
R. Frezzotti and G. C. Rossi, Chirally improving Wil- son fermions. II. Four-quark operators, JHEP10, 070, arXiv:hep-lat/0407002
work page internal anchor Pith review Pith/arXiv arXiv
-
[50]
C. Alexandrouet al.(Extended Twisted Mass), Strange and charm quark contributions to the muon anoma- lous magnetic moment in lattice QCD with twisted- mass fermions, Phys. Rev. D111, 054502 (2025), arXiv:2411.08852 [hep-lat]
-
[51]
Albaneseet al.(APE), Glueball Masses and String Tension in Lattice QCD, Phys
M. Albaneseet al.(APE), Glueball Masses and String Tension in Lattice QCD, Phys. Lett. B192, 163 (1987)
1987
- [52]
-
[53]
G. M. de Divitiis, R. Petronzio, and N. Tantalo, On the discretization of physical momenta in lattice QCD, Phys. Lett. B595, 408 (2004), arXiv:hep-lat/0405002
work page internal anchor Pith review Pith/arXiv arXiv 2004
-
[54]
Extraction of spectral densities from lattice correlators: decoupling signal from noise
A. Lupo and N. Tantalo, Extraction of spectral densities from lattice correlators: decoupling signal from noise, (2026), arXiv:2605.14652 [hep-lat]
work page internal anchor Pith review Pith/arXiv arXiv 2026
-
[55]
A. P. Valentine and M. Sambridge, Gaussian pro- cess models–I. A framework for probabilistic continu- ous inverse theory, Geophysical Journal International 220, 1632 (2020), https://academic.oup.com/gji/article- pdf/220/3/1632/31578341/ggz520.pdf
2020
-
[56]
C. K. Williams and C. E. Rasmussen,Gaussian processes 30 for machine learning, Vol. 2 (MIT press Cambridge, MA, 2006)
2006
-
[57]
L. Del Debbio, A. Lupo, M. Panero, and N. Tantalo, Bayesian solution to the inverse problem and its rela- tion to Backus–Gilbert methods, Eur. Phys. J. C85, 185 (2025), arXiv:2409.04413 [hep-lat]
-
[58]
G. Martinelli, M. Naviglio, S. Simula, and L. Vitto- rio,|V cb|, lepton flavor universality andSU(3) F symme- try breaking inB s →D (∗) s ℓνℓ decays through unitar- ity and lattice QCD, Phys. Rev. D106, 093002 (2022), arXiv:2204.05925 [hep-ph]
-
[59]
C. Oswaldet al.(Belle), Measurement of the inclu- sive semileptonic branching fractionB(B 0 s →X −ℓ+νℓ) at Belle, Phys. Rev. D87, 072008 (2013), [Erratum: Phys.Rev.D 90, 119901 (2014)], arXiv:1212.6400 [hep- ex]
work page internal anchor Pith review Pith/arXiv arXiv 2013
- [60]
-
[61]
G. Martinelli, S. Simula, and L. Vittorio, Updates on the determination of|V cb|, R(D∗) and|V ub|/|Vcb|, Eur. Phys. J. C84, 400 (2024), arXiv:2310.03680 [hep-ph]
-
[62]
G. Martinelli, S. Simula, and L. Vittorio, Semileptonic B→D ∗ decays from light toτleptons: the extraction of the form factorF 2 from data, Eur. Phys. J. C85, 242 (2025), arXiv:2410.17974 [hep-ph]
-
[63]
M. Bordone and A. Juttner, New strategies for probing B→D ∗ℓ¯νℓ lattice and experimental data, Eur. Phys. J. C85, 129 (2025), arXiv:2406.10074 [hep-ph]
-
[64]
A. Carvunis, G. Finauri, P. Gambino, M. Jung, and S. M¨ achler, New Physics in inclusive semileptonicBde- cays, JHEP01, 037, arXiv:2507.22123 [hep-ph]
-
[65]
J¨ ulich Supercomputing Centre, JUWELS Cluster and Booster: Exascale Pathfinder with Modular Supercom- puting Architecture at Juelich Supercomputing Centre, Journal of large-scale research facilities7, 10.17815/jlsrf- 7-183 (2021)
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.