{"id":"361a8962-feab-4ed2-ad43-4e3353caaf51","arxiv_id":"2411.18058","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"A lattice QCD study shows that subtracting the ground state reduces the Chebyshev truncation error in inclusive D_s semileptonic decays, and that modeled finite-volume corrections are small in the tested channel.","lead":"This paper tests two sources of error in lattice calculations of a quark decay process (D_s meson decaying inclusively to strange hadrons plus leptons and neutrinos). It finds that modeling the two-particle final state gives small finite-volume corrections, and that subtracting the known ground state before approximating shrinks the truncation error of the Chebyshev polynomials.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The finite-volume conclusion rests on a model with a free per-volume normalization fitted to one physical volume, so the claim that corrections are insignificant is not yet a controlled estimate.","rationale":"The reader's weakest_assumption identified the same finite-volume model and its free per-volume normalization s(L), so I agree. The paper's stronger contribution, the Chebyshev truncation-error analysis with ground-state subtraction, is supported by the data shown: the filled symbols are stable and consistent with the ground-state band, and the error prescription does not move central values. The finite-volume part is the weaker link because it is a single-channel model with a fitted normalization and an acknowledged missing D_s dependence. This does not make the paper unacceptable for a proceedings contribution, but it does justify the reader's CONDITIONAL verdict. No additional independent concern about the Chebyshev error estimate is strong enough to change the verdict, although a sensitivity study of alpha in sigma = 1/(alpha N) would be welcome. The verdict should remain CONDITIONAL, with the condition that the finite-volume study be repeated in other channels and on a second physical volume.","tokens_in":7561,"tokens_out":4687,"duration_ms":47506,"concrete_test":"Run the same ground-state-subtracted Chebyshev analysis on a second physical volume (e.g., L ~ 3.6 fm) with the same action and quark masses, and compare X_VV(q^2) with the L ~ 2.6 fm result without invoking the Eq. (13) model. If the two volumes disagree by more than the assigned finite-volume systematic, the model-based 'insignificant' claim is not reliable. A cheaper analytical check: test whether s(L) fitted at V = 48^3 and V = 256^3 follows the 1/L^3 scaling of Eq. (9) without an extra per-volume factor; if s(L) must absorb a different volume dependence, then the free normalization is masking the true finite-volume effect.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim that finite-volume corrections are insignificant (Sec. 5) depends entirely on the model in Eq. (13). Its excited-state term uses a non-interacting K Kbar spectrum with a vector-dominance form factor F(E) = 1/(E^2 - m_phi^2) and an overall per-volume normalization s(L) fitted to the lattice correlator (Sec. 4.1). Because s(L) is free and fitted, the model is not an ab initio prediction of the volume dependence. The comparison between V = 48^3 and V = 256^3 in Fig. 2 is dominated by the model's assumed spectrum and normalization, not by lattice data at a second physical volume; the 256^3 result is generated from the same model as a proxy for infinite volume. The paper explicitly states that the form factor neglects the initial D_s and that the study must be repeated for other channels (Sec. 4.1). These limitations are acknowledged, but they mean the conclusion 'finite-volume corrections are insignificant in our setup' is not yet a controlled systematic estimate: it is a plausible model-based statement that has not been validated against a second lattice volume or against channels with different angular-momentum kernels. Since the paper's stated purpose is to estimate systematic uncertainties, this is the most load-bearing weak point.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","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.","tokens_in":7862,"tokens_out":7241,"duration_ms":68697,"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":[{"comment":"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.","section":"§4.1, Eq. (13) and Fig. 2"},{"comment":"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).","section":"§3.2, Eq. (11) and Fig. 3"},{"comment":"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.","section":"§4.1, Eq. (9)"}],"minor_comments":[{"comment":"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}).","section":"Eq. (11)"},{"comment":"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.","section":"Fig. 2 (right panel)"},{"comment":"The caption repeatedly says 'Gound state contribution'; this should be 'Ground state contribution'.","section":"Fig. 3"},{"comment":"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.","section":"Secs. 4.1 and 4.2"},{"comment":"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.","section":"Eq. (3)"}],"recommendation":"major_revision","confidential_remarks":"This is a proceedings contribution, and the authors are transparent about the limitations of their model. My recommendation is driven by the gap between the strength of the finite-volume conclusion and the model evidence behind it; this is fixable by rephrasing the claims and adding sensitivity checks, so I do not regard it as grounds for rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: the Chebyshev truncation analysis is the real content here, and it's in good shape; the finite-volume piece is a plausible model but not a controlled estimate. Read it for the ground-state subtraction, and treat the FV conclusion as provisional.\n\nWhat's new: the error estimate that uses the boundedness of Chebyshev polynomials and random sampling of higher-order coefficients, combined with treating the ground state exactly. The plots in Fig. 3 make the case: without subtraction, the truncation error explodes at larger recoil; with subtraction it stays flat and consistent with a ground-state-only estimate. That is a practical step forward for the inclusive program, and the estimate is conservative by construction—it only adds to the error, it does not move central values. Good, and the math behind it is solid.\n\nThe finite-volume part is weaker. The model in Eq. (13) uses a non-interacting K Kbar spectrum, a vector-dominance form factor that the authors admit neglects the initial D_s, and a per-volume normalization s(L) fitted to the correlator. The comparison between V = 48^3 and 256^3 is generated entirely within the model; the 256^3 curve is not a second lattice volume. So the statement that finite-volume corrections are insignificant is a model-based plausibility argument, not a controlled systematic estimate. The authors are upfront about this limitation—they say the study needs to be repeated for other channels—but the conclusion in Sec. 5 reads stronger than the evidence supports. Also, the sigma = 1/(alpha N) mapping with alpha = 1 is an ansatz, and there is no sensitivity study for alpha. These are real soft spots, though they are acknowledged rather than hidden.\n\nFor a proceedings, this is acceptable as a status report. No code or data is shipped, which is normal for this venue, but it means the lattice results cannot be independently checked from the manuscript alone. Still, the Chebyshev error estimate is a genuine methodological contribution that deserves referee attention.\n\nWho gets value: lattice QCD practitioners working on inclusive decays and spectral reconstruction. A serious referee should engage with it. My recommendation: engage, but require the authors to soften the finite-volume claim and add a brief sensitivity note on alpha.","headline":"The Chebyshev truncation error estimate is sound and useful; the finite-volume conclusion is a model-based plausibility argument, not yet a controlled systematic.","tokens_in":8361,"tokens_out":2070,"would_cite":false,"duration_ms":19200,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["12.38.Gc","13.20.Fc"],"model":"deepseek-v4-flash","headline":"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.","keywords":["inclusive semileptonic decays","lattice QCD","Chebyshev polynomial approximation","finite-volume effects","D_s meson","spectral density","systematic uncertainties"],"falsifier":"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.","tokens_in":7339,"feed_emoji":"⚛️","tokens_out":5497,"duration_ms":46373,"temperature":0.7,"pith_summary":"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.","feed_headline":"Ground-state subtraction stabilizes lattice decay-rate errors","feed_subtitle":"A lattice QCD path to inclusive D_s decays now controls its main systematic uncertainties.","key_machinery":"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)$.","core_discovery":"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.","pith_inferences":["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."],"forward_implications":["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."],"supporting_citations":[{"why":"Previous work by this group that first estimated the Chebyshev approximation systematic in inclusive semileptonic decays; the present paper extends that error estimate.","marker":"[1]"},{"why":"Prior study of finite-volume effects in inclusive semileptonic decays of charmed mesons, on which the current finite-volume model builds.","marker":"[2]"},{"why":"Introduced the Chebyshev-polynomial kernel approximation for inclusive semileptonic decays from lattice QCD; this is the method whose truncation error is analyzed here.","marker":"[3]"},{"why":"Provides the practical implementation details, including the definitions of the coefficients $\\tilde a_j^{(t)}$ and the Chebyshev matrix-element fitting procedure used in the analysis.","marker":"[8]"},{"why":"Foundational result on the volume dependence of the energy spectrum in quantum field theory, which underlies the finite-volume modeling of the present paper.","marker":"[9]"}],"fun_headline_variants":["Ground-state subtraction tames lattice D_s decay errors","Lattice D_s: truncation and volume systematics under leash","Inclusive D_s decays: subtraction shrinks systematic errors","Lattice QCD method pins down inclusive D_s rate systematics","Controlled systematics for semileptonic D_s on the lattice"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"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.","fun_headline_variants_meta":{"raw":{"variants":["Ground-state subtraction tames lattice D_s decay errors","Lattice D_s: truncation and volume systematics under leash","Inclusive D_s decays: subtraction shrinks systematic errors","Lattice QCD method pins down inclusive D_s rate systematics","Controlled systematics for semileptonic D_s on the lattice"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000186,"raw_usage":{"total_tokens":1281,"prompt_tokens":860,"completion_tokens":421,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":476,"completion_tokens_details":{"reasoning_tokens":334}},"tokens_in":476,"tokens_out":421,"duration_ms":4932,"temperature":1.0,"reasoning_tokens":334,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T11:33:54.032996+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"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.","supporting_citations":[{"cited_title":"Inclusive semi-leptonic decays of charmed mesons with M\\\"obius domain wall fermions","cited_arxiv_id":"2211.16830","evidence_quote":"Previous work by this group that first estimated the Chebyshev approximation systematic in inclusive semileptonic decays; the present paper extends that error estimate."},{"cited_title":"Luscher,Volume Dependence of the Energy Spectrum in Massive Quantum Field Theories","cited_arxiv_id":null,"evidence_quote":"Foundational result on the volume dependence of the energy spectrum in quantum field theory, which underlies the finite-volume modeling of the present paper."}],"review_version":1}