Pith. sign in

REVIEW 3 major objections 5 minor 2 cited by

Systematic effects in the lattice calculation of inclusive semileptonic decays

T0 review · 3 major / 5 minor · reviewed 2026-08-12 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read By treating the D_s ground state exactly and applying the Chebyshev kernel approximation only to excited states, the dominant truncation error in lattice inclusive semileptonic decays is largely removed.

desk verdict The Chebyshev truncation error estimate is sound and useful; the finite-volume conclusion is a model-based plausibility argument, not yet a controlled systematic. read the letter →

arxiv 2411.18058 v2 pith:I3OSB5JC submitted 2024-11-27 hep-lat

Ryan Kellermann , Alessandro Barone , Ahmed Elgaziari , Shoji Hashimoto , Zhi Hu , Andreas Jüttner , Takashi Kaneko This is my paper · ORCID
classification hep-lat PACS 12.38.Gc13.20.Fc
keywords inclusivesemileptonicdecayslatticeQCDChebyshevpolynomialapproximationfinite-volumeeffectsD_smesonspectraldensitysystematicuncertainties
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

This paper reports progress on computing the inclusive semileptonic decay $D_s \to X_s\ell\nu_\ell$ from lattice QCD, where the full decay rate is obtained without enumerating exclusive final states. The authors focus on two systematic errors that stand between the lattice correlator and the physical rate: the error from approximating the integration kernel by a finite Chebyshev polynomial, and the error from the finite simulation volume. Their central result is that subtracting the ground-state contribution—treating it exactly in the energy integral—before applying the Chebyshev approximation to the remaining correlator keeps the truncation error small and stable even at large recoil momenta, where the full-data analysis suffers. For finite-volume effects, they build an explicit model of the dominant two-body $K\bar K$ final state and find the corrections small in the channel studied. If these strategies hold in other channels, they remove the two largest obstacles to a controlled first-principles inclusive rate.

What carries the argument

The central objects are the shifted Chebyshev polynomials $\tilde T_j(x)$ with $x=e^{-\omega}$, whose expansion coefficients $\tilde c_{\mu\nu,k}^{(l)}$ are known analytically and whose matrix elements are fit from the correlator via $\bar C(t)=\sum_j \tilde a_j^{(t)}\langle \tilde T_j\rangle$, together with the decomposition of the spectral density into an exact ground-state delta function plus an excited-state continuum, $\rho(\omega)=\rho_0\delta(\omega-m_X)+\rho_{Ex}(\omega)$. The Chebyshev order $N$ and the sigmoid smearing width are linked by $\sigma=1/N$, so the $\sigma\to 0$ and $N\to\infty$ limits are taken together, and the truncation error is bounded by drawing the uncomputed matrix elements uniformly from $[-1,1]$ and taking the standard deviation. The finite-volume model replaces the continuum two-body phase space by the finite-volume sum $\rho_V(\omega)=(\pi/V)\sum_{\mathbf q} q^2/(4(q^2+m_K^2))\,\delta(\omega-2\sqrt{q^2+m_K^2})$ and fits the correlator to $C(t)=A_0 e^{-E_0 t}+s(L)\sum_i A_i e^{-E_i t}F(E_i)$.

What would settle it

Run the same ground-state-subtracted analysis on a significantly larger volume (for example $L\simeq 4$--$5$ fm rather than $2.6$ fm) with the same action and currents, and compare the inclusive rate to the model extrapolation used here; if the result shifts by more than the model-predicted volume dependence, the finite-volume model is the wrong description.

Watch

Extended reading notes

Core claim

The paper establishes that the inclusive $D_s$ semileptonic decay rate can be computed on the lattice with controlled systematics using two techniques. First, the Chebyshev-polynomial reconstruction of the kernel $K_\sigma(\omega)$ is applied not to the full correlator but to the correlator with the ground state removed; the ground-state amplitude and energy are extracted from a single-exponential fit at large times and handled exactly in the energy integral. This reduces the truncation error of the polynomial approximation, which otherwise grows rapidly with recoil momentum as the kinematical phase space narrows. Second, finite-volume corrections are estimated by a model in which the excited-state spectral density is a non-interacting $K\bar K$ pair weighted by a vector-dominance form factor $F(E)=1/(E^2-m_\phi^2)$, with a volume-dependent normalization $s(L)$ fit to the lattice correlator; the model reproduces the lattice data and predicts negligible volume dependence in the spatial-current channel examined. The paper concludes that the truncation error is the more important systematic and that the ground-state-subtraction strategy largely removes it.

Load-bearing premise

The estimate that finite-volume corrections are negligible rests on a model where the multi-hadron final state is a non-interacting $K\bar K$ pair with a vector-dominance form factor and a volume-dependent normalization fitted to the correlator; if interactions or other channels make the model miss the true volume dependence, the finite-volume error would be larger than estimated.

Editorial extensions

If this is right

  • The truncation error of the Chebyshev approximation can be drastically reduced by exact ground-state treatment, making the inclusive decay rate accessible for large recoil momenta where the full-data approach drifts.
  • Finite-volume corrections in the spatial-current channel are small and can be extrapolated with a simple two-body model, so large-volume simulations may not be needed for this channel.
  • The relation $\sigma=1/N$ ties the two limits (kernel smearing and polynomial truncation) into one, giving a single parameter to control the systematic error.
  • The combination of exact ground state plus Chebyshev-excited states gives results consistent with a ground-state-only estimate, indicating that excited-state contamination is small in this channel.
  • The finite-volume and truncation-error strategies must be repeated for other channels, as the small corrections seen here are not guaranteed to persist.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • The ground-state-subtraction strategy should apply equally to other kernel-reconstruction approaches, such as the Hansen-Lupo-Tantalo method, where the same ground-state-exact decomposition could reduce the smearing dependence.
  • The uniform-distribution bound on the uncomputed Chebyshev coefficients is likely conservative; a data-driven estimate of the omitted coefficients from the fit covariance could shrink the quoted truncation error without changing the central value.
  • The finite-volume model could be upgraded by including the $D_s$ dependence explicitly and by adding $K\bar K$ interactions through the L\"uscher formalism, moving beyond the non-interacting approximation and the single form-factor choice.
  • If the two-step procedure (ground-state subtraction plus finite-volume modeling) works in other channels, it could be applied to inclusive decays of $B$ or $B_s$ mesons, a central target for first-principles determinations of CKM matrix elements.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

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

Referee Report

3 major / 5 minor

Summary. This proceedings paper reports an update on the lattice calculation of the inclusive semileptonic decay D_s -> X_s l nu, focusing on two systematic effects. The authors model finite-volume effects by fitting the excited-state part of a four-point correlator with a non-interacting K Kbar spectrum, a vector-dominance form factor, and a per-volume normalization s(L) (Eq. (13)), and they use this model to compare V=48^3 with a 256^3 proxy for infinite volume. They also study the Chebyshev approximation error by relating polynomial order to smearing width (Eq. (10)) and by estimating the contribution of truncated Chebyshev matrix elements from a uniform distribution in [-1,1] (Eq. (11)). Their numerical results for the considered channel show that the model describes the correlator well, that the volume and smearing dependence is mild, and that treating the ground state exactly before applying the Chebyshev reconstruction reduces the truncation error and stabilizes the central value.

Significance. If the claims hold, the main value is methodological: the paper gives a concrete, falsifiable model for finite-volume corrections and a simple procedure (exact ground-state treatment) that visibly improves the convergence of the Chebyshev reconstruction. The authors are explicit about several limitations, including the neglect of the initial D_s in the form factor and the need to repeat the analysis for other channels. The paper is appropriate as a proceedings progress report, but the strength of the finite-volume conclusion currently exceeds what the model can support, and the Chebyshev truncation error estimate contains a distributional assumption that should be flagged as such.

major comments (3)
  1. [§4.1, Eq. (13) and Fig. 2] The conclusion in Sec. 5 that finite-volume corrections are insignificant in this setup is not yet a controlled estimate. The model in Eq. (13) contains a free per-volume normalization s(L) fitted to the V=48^3 lattice correlator, and the V=256^3 curve is generated from the same model rather than from lattice data at a second physical volume. The agreement shown in Fig. 2 therefore validates the model's ability to describe the correlator, but it does not independently constrain the infinite-volume extrapolation. Together with the acknowledged neglect of the initial D_s in F(E) and the restriction to one channel, this means the finite-volume statement should be presented as a model-dependent estimate, or supplemented by a second-volume or multi-channel cross-check.
  2. [§3.2, Eq. (11) and Fig. 3] The truncation error estimate for the Chebyshev approximation is based on an additional assumption that is not stated as such: the unreconstructed matrix elements <T_j> are drawn from a uniform distribution in [-1,+1]. Boundedness of the Chebyshev polynomials only gives |<T_j>| <= 1; the standard deviation of the uniform distribution is a modeling choice, and the relation sigma = 1/(alpha N) with alpha=1 is also an assumption. The numerical size of the error bars in Fig. 3, and hence the quantitative strength of the improvement in Sec. 4.2, depends on these choices. Please state explicitly that this is a prior/model assumption and test the sensitivity to the assumed distribution (e.g., uniform versus a Gaussian or a worst-case bound).
  3. [§4.1, Eq. (9)] The finite-volume model uses the non-interacting K Kbar spectrum, but the channel considered contains the phi resonance. Because the discretized levels in an interacting two-body system are shifted by the scattering phase shift (cf. Ref. [9]), the free-spectrum approximation may not be reliable for the volume dependence of smeared spectral densities. The conclusion that volume effects are insignificant should either be restricted to the non-interacting model or accompanied by an estimate of the sensitivity to the interaction, since the current text cites Luscher's work as motivation for the importance of scattering states.
minor comments (5)
  1. [Eq. (11)] In the second sum of Eq. (11), the last term uses T_j(e^{-omega}) while the summation index is k; it should be T_k(e^{-omega}).
  2. [Fig. 2 (right panel)] The horizontal axis label 'th' should be written as omega_th, and the legend entries such as 'K^(2), V = 48^3' are ambiguous; please define the symbol K^(2) in the caption.
  3. [Fig. 3] The caption repeatedly says 'Gound state contribution'; this should be 'Ground state contribution'.
  4. [Secs. 4.1 and 4.2] The text in Sec. 4.1 refers to X^parallel_AA for spatial axial-current insertions, while Sec. 4.2 and Fig. 3 use X^parallel_VV for what appears to be the same or closely related quantity; please harmonize the notation.
  5. [Eq. (3)] The sigmoid function and the smearing parameter sigma are used in Eq. (3) but not defined in the text; a brief definition would help the reader.

Circularity Check

0 steps flagged · score 2.0 of 10

No significant circularity: the finite-volume conclusion is model-based but not a circular reduction, and the Chebyshev error estimate is self-contained.

full rationale

The paper does not exhibit a circular reduction. The Chebyshev truncation-error estimate (Sec. 3.2) is constructed from the explicit bounded-polynomial property |T~_j|≤1 and a random uniform prior on unresolved coefficients; it does not use the target X values as input. The finite-volume analysis (Sec. 4.1) fits Eq. (13) to the V=48^3 correlator with a free normalization s(L), then uses the same model to generate the V=256^3 curve. This is a model-based extrapolation, not a parameter-free prediction, but the ω_th-shape and the volume dependence are not determined by the single fitted normalization s(L); they come from the stated non-interacting KK spectrum and vector-dominance form factor. The paper itself scopes the result: 'This form factor does not take into account the initial D_s' and 'This study needs to be repeated for other channels.' These are acknowledged limitations rather than circular reasoning. The 'nice agreement' between lattice data and model after fitting s(L) is partly a fit-quality statement, but the volume-extrapolation conclusion is not statistically forced by that parameter. Self-citations [1,2,8] are used for methodological details, but the relevant equations are reproduced in this paper, so the conclusions do not rest on an unverified self-citation. Score 2 reflects minor caveats only.

Assumptions & free parameters 4 free parameters · 7 assumptions · 0 invented entities

The central systematic conclusions rest on the fitted or chosen parameters listed above, plus the modeling assumptions in the axioms. The Chebyshev truncation error estimate is cheap and conservative. The finite-volume significance claim is the most assumption-heavy part, since its template normalization is fitted to the data it is used to interpret.

free parameters (4)
  • s(L) = not specified in paper
    Free normalization of the excited-state sum in the finite-volume model, Eq. (13), fitted per volume to the lattice correlator; it absorbs the overall strength of the K Kbar contribution and therefore shapes the finite-volume correction estimate.
  • alpha in sigma = 1/(alpha N) = 1
    Introduced in Eq. (10) to combine the smearing and truncation limits; set equal to one by hand with no derivation or sensitivity study. It determines the mapping between Chebyshev order and smearing width used for the error estimate.
  • N_Cut = not specified
    Highest Chebyshev order that can be reconstructed from the data; used in Eq. (11) to separate reliably fitted coefficients from the randomized higher-order tail. The choice affects the size of the truncation error estimate.
  • t0 = 1/2
    Chosen to maximize the available data for the Chebyshev fit (Sec. 2); this choice defines the kernel and the correlator normalization, but it is not tuned to make the result come out a particular way.
assumptions (7)
  • standard math The lattice four-point correlator C_mu nu(q,t) is the Laplace transform of the hadronic tensor W_mu nu(q,omega), Eq. (4)
    Spectral representation of the correlator, standard in lattice QCD; the basis of the Chebyshev reconstruction method.
  • domain assumption The physical inclusive rate is obtained by the ordered limit lim_{sigma->0} lim_{V->infinity} Xbar_sigma(q^2), Eq. (7)
    Assumed in the smeared-spectral-density framework; requires that smearing be removed after the infinite-volume limit.
  • domain assumption The dominant multi-hadron final state for the finite-volume study is a non-interacting K Kbar pair, Eq. (8)
    The model ignores interactions between K and Kbar as well as other multi-hadron channels; the finite-volume correction estimate depends on this dominance.
  • ad hoc to paper The vector-dominance form factor F(E)=1/(E^2-m_phi^2) models the transition weight and does not depend on the initial D_s
    Introduced in Sec. 3.1 to make the model more realistic; the paper explicitly notes it does not take into account the initial D_s.
  • ad hoc to paper Chebyshev order N and smearing width sigma are related by sigma = 1/(alpha N) with alpha=1
    Chosen in Sec. 3.2 to combine the N->infinity and sigma->0 limits; no derivation is given for this identification.
  • ad hoc to paper Higher-order Chebyshev matrix elements in the unreconstructed tail follow a uniform distribution in [-1,+1]
    Used in Sec. 3.2 to estimate the truncation error; this is a conservative prior exploiting boundedness, not a measured distribution.
  • domain assumption The ground state of the channel can be described by a single exponential and is separable from excited states
    Standard assumption in lattice spectroscopy when a plateau is present; used as the basis of the ground-state subtraction strategy.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Systematic effects in the lattice calculation of inclusive semileptonic decays." pith.science (2026). https://pith.science/paper/I3OSB5JC

@misc{pith2026241118058,
  author       = {Pith},
  title        = {Pith review of: Systematic effects in the lattice calculation of inclusive semileptonic decays},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/I3OSB5JC}},
  note         = {Machine review of arXiv:2411.18058}
}
abstract

We report on the calculation of the inclusive semileptonic decay of the $D_s$ meson on the lattice. We simulate the $D_s \rightarrow X_s\ell\nu_\ell$ process with M\"obius domain-wall charm and strange quarks, whose masses are approximately tuned to their physical values. Our simulations cover the whole kinematical region. The focus of this work is to present updates on our strategies towards estimating the systematic uncertainties in the determination of the inclusive decay rate. We specifically focus on the systematic errors due to the choice of our approximation strategy and finite-volume effects.

Figures

Figures reproduced from arXiv: 2411.18058 by the authors.

Figure 1
Figure 1. Comparison to highlight the differences in the kernel function depending on the choice of the Heaviside or sigmoid function. The Heaviside function in both plots is represented by the solid line, while the sigmoid function uses the dashed line. The smearing width of the sigmoid function is 𝜎 = 0.1. 𝑁 → ∞. In [PITH_FULL_IMAGE:figures/full_fig_p005_1.png] view at source ↗
Figure 2
Figure 2. Left: Four-point correlation function for the spatial components⟨𝐴 † 𝑖 (𝑡)𝐴𝑖(0)⟩ with zero momentum insertion. The dashed line represents the single-exponential fit to determine the ground state. The black dash-dotted line represents the fit results using the prescription (13). Right: A comparison of the model prediction and the lattice data for the finite-volume corrections. The model prediction is obtained for the… view at source ↗
Figure 3
Figure 3. Errors of 𝑋¯ ∥ 𝑉𝑉 (𝑞 2 ) as a function of 𝜎 = 1/𝑁. Estimates are plotted for all values of 𝑞 used in our simulations. The gray band represents the expected ground state contribution. The empty symbols are obtained from applying the Chebyshev approximation on the full data set, while the filled symbols are obtained treating the ground-state contribution exactly and using the Chebyshev approximation only on the remain… view at source ↗

Discussion (0). Continue with ORCID to comment.

Forward citations

Cited by 2 Pith papers

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

  1. The NRQCD $\Upsilon$ spectrum at non-zero temperature using Backus-Gilbert regularisations

    hep-lat 2025-02 conditional novelty 4.0 of 10

    On FASTSUM Generation 2L ensembles, Backus-Gilbert and Tikhonov reconstructions agree on the Upsilon ground-state mass, but the decay width is dominated by reconstruction smearing and high-temperature broadening is un...

  2. Approaching the Inverse Problem: Toward Lattice QCD Calculations of Inclusive Hadronic Quantities

    hep-lat 2025-01 conditional novelty 4.0 of 10

    For lattice QCD spectral reconstruction, the Wertevorrat from Nevanlinna-Pick interpolation bounds the analytic-continuation uncertainty, and in a model test this bound falls roughly exponentially with the number of d...

Reference graph

Works this paper leans on

10 extracted references · 3 canonical work pages · cited by 2 Pith papers

  1. [9]

    Luscher,Volume Dependence of the Energy Spectrum in Massive Quantum Field Theories

    M. Luscher,Volume Dependence of the Energy Spectrum in Massive Quantum Field Theories

    work page

  2. [1]

    Inclusive semi-leptonic decays of charmed mesons with M\"obius domain wall fermions

    R. Kellermann, A. Barone, S. Hashimoto, A. Jüttner and T. Kaneko,Inclusive semi-leptonic decays of charmed mesons with Möbius domain wall fermions, PoSLATTICE2022(2023) 414 [2211.16830]

    work page Pith review arXiv 2023

  3. [2]

    Kellermann, A

    R. Kellermann, A. Barone, S. Hashimoto, A. Jüttnerc and T. Kanekoa, Studies on finite-volume effects in the inclusive semileptonic decays of charmed mesons, PoS LATTICE2023(2024) 272 [2312.16442]

    arXiv 2024

  4. [3]

    Gambino and S

    P. Gambino and S. Hashimoto,Inclusive Semileptonic Decays from Lattice QCD, Phys. Rev. Lett.125 (2020) 032001 [2005.13730]. 8 Systematic effects in the lattice calculation of inclusive semileptonic decays R. Kellermann

    arXiv 2020

  5. [4]

    P.Gambino,S.Hashimoto,S.Mächler,M.Panero,F.Sanfilippo,S.Simulaetal., LatticeQCD study of inclusive semileptonic decays of heavy mesons,JHEP07 (2022) 083 [2203.11762]

    arXiv 2022

  6. [5]

    Hansen, H.B

    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. D 96(2017) 094513 [1704.08993]

    arXiv 2017

  7. [6]

    Hansen, A

    M. Hansen, A. Lupo and N. Tantalo,Extraction of spectral densities from lattice correlators, Phys. Rev. D99(2019) 094508 [1903.06476]

    arXiv 2019

  8. [7]

    Bulava, M.T

    J. Bulava, M.T. Hansen, M.W. Hansen, A. Patella and N. Tantalo,Inclusive rates from smeared spectral densities in the two-dimensional O(3) non-linear𝜎-model,JHEP07 (2022) 034 [2111.12774]

    arXiv 2022

Show all 10 references
  1. [8]

    Barone, S

    A. Barone, S. Hashimoto, A. Jüttner, T. Kaneko and R. Kellermann,Approaches to inclusive semileptonic B(𝑠)-meson decays from Lattice QCD, JHEP07(2023) 145 [2305.14092]

    2023 arXiv

  2. [10]

    Scattering States, Commun. Math. Phys.105 (1986) 153. 9

    1986

Pith tools

Reviewed August 12, 2026 · model on record in the stance chip above.