pith. sign in

arxiv: 2605.14652 · v1 · pith:ARXEK4V3new · submitted 2026-05-14 · ✦ hep-lat

Extraction of spectral densities from lattice correlators: decoupling signal from noise

Alessandro Lupo , Nazario Tantalo This is my paper

Pith reviewed 2026-06-30 19:51 UTC · model grok-4.3

classification ✦ hep-lat
keywords spectral densitieslattice correlatorsinverse problemtruncationstatistical noiseBackus-GilbertEuclidean correlatorssmeared spectral density
0
0 comments X

The pith

The inverse problem solution for smeared spectral densities decomposes into terms where the largest noise contributors add the least to the central value, shifting control of systematics to optimal truncation.

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

The paper develops an alternative method to extract smeared spectral densities from Euclidean lattice correlators that avoids Backus-Gilbert regularization. It observes that the solution decomposes into a sum of terms ordered such that those dominating statistical noise contribute minimally to the central value. Systematics analysis therefore reduces to locating the truncation point at which the signal saturates before noise growth dominates. The procedure is tested both standalone and as a complement to stability checks in the regulated Backus-Gilbert approach. A reader would care because it offers a direct route to stable extractions from noisy correlators without introducing regularization bias.

Core claim

The solution can be decomposed into a sum of terms, in the spirit of the singular value decomposition, where those with the largest contribution to the statistical noise happen to contribute the least to the central value of the smeared spectral density. The analysis of the systematics of the inverse problem is then shifted to finding the optimal truncation of such summation, so that the signal is saturated before the noise explodes. The performance and systematics of this approach are scrutinized either as a standalone procedure or to complement the stability analysis required to extrapolate the unbiased result in the Backus-Gilbert regulated version of the solution.

What carries the argument

SVD-like decomposition of the inverse-problem solution into terms ordered by relative noise versus central-value contribution.

If this is right

  • The method supplies a regularization-free route to smeared spectral densities.
  • Systematics control reduces to locating a single truncation threshold.
  • The same decomposition supplies an independent stability diagnostic when used alongside Backus-Gilbert extrapolation.
  • Performance can be checked either in isolation or in hybrid mode on the same correlator ensemble.

Where Pith is reading between the lines

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

  • The observed ordering of noise versus signal contributions may appear in other linear inverse problems that arise from Euclidean-to-Minkowski transforms.
  • An automated criterion for choosing the truncation, such as monitoring the variance of successive partial sums, could be added without altering the underlying decomposition.
  • If the separation holds across different lattice actions or volumes, the truncation point might become a universal diagnostic rather than a fit parameter.

Load-bearing premise

An optimal truncation point exists and can be identified such that the signal saturates before the noise explodes.

What would settle it

A lattice correlator data set in which the central-value contribution from high-noise terms exceeds the low-noise terms before any clear saturation plateau appears in the truncated sum.

Figures

Figures reproduced from arXiv: 2605.14652 by Alessandro Lupo, Nazario Tantalo.

Figure 1
Figure 1. Figure 1: FIG. 1. Eigenvalues [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Values of the coefficients used to approximate a Gaus [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. An instance of synthetic correlator used in this work [PITH_FULL_IMAGE:figures/full_fig_p005_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. Partial sums leading to the smeared spectral density [PITH_FULL_IMAGE:figures/full_fig_p006_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6. Top: example of stability analysis to remove the [PITH_FULL_IMAGE:figures/full_fig_p007_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7. Smeared spectral density from the sum in Eq. ( [PITH_FULL_IMAGE:figures/full_fig_p008_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: FIG. 8. Top: single example of mock vector-vector correlator [PITH_FULL_IMAGE:figures/full_fig_p009_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: FIG. 9. Histogram for the pull defined in Eq. ( [PITH_FULL_IMAGE:figures/full_fig_p009_9.png] view at source ↗
Figure 11
Figure 11. Figure 11: FIG. 11. Closure tests on the same datasets, and with the [PITH_FULL_IMAGE:figures/full_fig_p010_11.png] view at source ↗
Figure 10
Figure 10. Figure 10: FIG. 10. Companion plot of Fig [PITH_FULL_IMAGE:figures/full_fig_p010_10.png] view at source ↗
Figure 13
Figure 13. Figure 13: FIG. 13. The figure updates Figs [PITH_FULL_IMAGE:figures/full_fig_p011_13.png] view at source ↗
read the original abstract

We expand the treatment of the problem of the extraction of smeared spectral densities from Euclidean correlators introduced in [Phys. Rev. D 99, 094508], providing an alternative which does not rely on the Backus-Gilbert regularization. This is possible due to the observation that the solution can be decomposed into a sum of terms, in the spirit of the singular value decomposition, where those with the largest contribution to the statistical noise happen to contribute the least to the central value of the smeared spectral density. The analysis of the systematics of the inverse problem is then shifted to finding the optimal truncation of such summation, so that the signal is saturated before the noise explodes. We scrutinise the performance and systematics of this approach either as a standalone procedure, or to complement the stability analysis required to extrapolate the unbiased result in the Backus-Gilbert regulated version of the solution.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

0 major / 2 minor

Summary. The manuscript expands the treatment of extracting smeared spectral densities from Euclidean lattice correlators, offering an alternative to Backus-Gilbert regularization. It observes that the solution decomposes into SVD-like terms where those dominating statistical noise contribute least to the central value of the smeared spectral density. Systematics analysis is shifted to identifying an optimal truncation that saturates the signal before noise dominates. Performance is scrutinized both as a standalone procedure and as a complement to Backus-Gilbert for stability and extrapolation.

Significance. If the decomposition and truncation strategy holds under numerical validation, the work supplies a practical alternative for the inverse problem in lattice spectral density extraction, reducing dependence on regularization while shifting focus to truncation systematics. Explicit testing both standalone and in tandem with Backus-Gilbert is a clear strength, as is the emphasis on separating noise-dominant from signal-dominant contributions. This could improve error control and reproducibility in non-perturbative QCD calculations of smeared spectral functions.

minor comments (2)
  1. [Abstract] Abstract: the description of the decomposition would benefit from a brief reference to the defining equation or matrix whose SVD is performed, to make the central observation immediately traceable.
  2. The manuscript should include a dedicated section or subsection detailing the lattice ensembles, correlator sources, and smearing kernels used in the numerical tests, together with quantitative measures (e.g., bias vs. variance trade-off curves) that demonstrate the existence and identifiability of the optimal truncation point.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for their positive summary and recommendation of minor revision. No specific major comments were listed in the report, so we have no individual points requiring detailed response or manuscript changes at this stage. We appreciate the recognition that the SVD-style decomposition and truncation approach offers a practical alternative for the inverse problem.

Circularity Check

0 steps flagged

No significant circularity

full rationale

The paper's central construction rests on the stated observation that the inverse-problem solution admits an SVD-like decomposition in which noise-dominant terms contribute least to the smeared spectral density; truncation is then introduced as the control parameter. This decomposition is presented as an independent mathematical property of the correlator-to-density map rather than a quantity fitted from the target result or defined by the truncation itself. The cited prior work (Phys. Rev. D 99, 094508) supplies only the problem statement and the BG baseline; the new truncation procedure is scrutinized both standalone and as a complement, without any load-bearing step reducing to a self-citation chain or to a parameter that is statistically forced by construction. No equation or claim in the abstract or described method equates the output to its inputs by definition.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

Abstract provides no explicit free parameters, axioms, or invented entities; all details are deferred to the full text which was unavailable.

pith-pipeline@v0.9.1-grok · 5675 in / 934 out tokens · 39963 ms · 2026-06-30T19:51:35.897852+00:00 · methodology

discussion (0)

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

Forward citations

Cited by 3 Pith papers

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

  1. Inclusive $\bar B_s\mapsto X_{\bar sc} \ell \bar \nu$ decays from lattice QCD: computational strategy and a first physical result

    hep-lat 2026-07 unverdicted novelty 7.0

    Lattice QCD strategy interpolates sub-physical mass data with OPE to obtain the inclusive anti-B_s to X_sc l nu decay rate, yielding 7% precision on limited ETMC ensembles via a new four-point correlator method.

  2. The Causal Bootstrap: Bounding Smeared Spectral Functions from Non-Perturbative Euclidean Data

    hep-lat 2026-05 unverdicted novelty 7.0

    The causal bootstrap computes rigorous bounds on smeared spectral functions from non-perturbative Euclidean data by optimizing over the convex set of compatible positive spectral densities and reducing dual problems t...

  3. Kernel transformations and bounds for smeared spectral functions

    hep-lat 2026-06 unverdicted novelty 5.0

    Framework for exact and approximate kernel transformations between smeared spectral functions, including systematic error bounds computable from input data.

Reference graph

Works this paper leans on

26 extracted references · 21 canonical work pages · cited by 3 Pith papers · 4 internal anchors

  1. [1]

    On the extraction of spectral densities from lattice correlators

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

    work page internal anchor Pith review Pith/arXiv arXiv 2019

  2. [2]

    Barata and K

    J.C.A. Barata and K. Fredenhagen,Particle scattering in Euclidean lattice field theories,Commun. Math. Phys.138(1991) 507

    1991

  3. [3]

    Bertero, C.D

    M. Bertero, C.D. Mol and E.R. Pike,Linear inverse problems with discrete data: Ii. stability and regularisation,Inverse Problems4(1988) 573

    1988

  4. [4]

    Furmanski and R

    W. Furmanski and R. Petronzio,A Method of Analyzing the Scaling Violation of Inclusive Spectra in Hard Processes,Nucl. Phys. B195(1982) 237

    1982

  5. [5]

    Backus and F

    G. Backus and F. Gilbert,The Resolving Power of Gross Earth Data,Geophysical Journal International 16(1968) 169

    1968

  6. [6]

    Pijpers and M.J

    F.P. Pijpers and M.J. Thompson,The sola method for helioseismic inversion,Astronomy and Astrophysics (ISSN 0004-6361), vol. 281, no. 1, p. 231-240281 (1994) 231

    1994

  7. [7]

    Bulava, M

    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]

    work page arXiv 2022

  8. [8]

    Del Debbio, A

    L. Del Debbio, A. Lupo, M. Panero and N. Tantalo, Bayesian solution to the inverse problem and its relation to Backus–Gilbert methods,Eur. Phys. J. C85 (2025) 185 [2409.04413]

    work page arXiv 2025

  9. [9]

    Buzzicotti, A

    M. Buzzicotti, A. De Santis and N. Tantalo,Teaching to extract spectral densities from lattice correlators to a broad audience of learning-machines,Eur. Phys. J. C 84(2024) 32 [2307.00808]

    work page arXiv 2024

  10. [10]

    Del Debbio, A

    L. Del Debbio, A. Lupo, M. Panero and N. Tantalo, Multi-representation dynamics of SU(4) composite Higgs models: chiral limit and spectral reconstructions, Eur. Phys. J. C83(2023) 220 [2211.09581]. 5 The approach of Ref. [27] considers a restricted class of smearing kernels that can be represented asS(E) = R ∞ 0 dt g(t)e−tE, i.e. as integrals that corresp...

    work page arXiv 2023

  11. [11]

    Bennett et al.,Meson spectroscopy from spectral densities in lattice gauge theories,Phys

    E. Bennett et al.,Meson spectroscopy from spectral densities in lattice gauge theories,Phys. Rev. D110 (2024) 074509 [2405.01388]. [12]Extended Twisted Masscollaboration,Probing the Energy-Smeared R Ratio Using Lattice QCD,Phys. Rev. Lett.130(2023) 241901 [2212.08467]. [13]Extended Twisted Masscollaboration,Inclusive Hadronic Decay Rate of theτLepton from...

    work page arXiv 2024

  12. [12]

    De Santis et al.,Inclusive semileptonic decays of the Ds meson: A first-principles lattice QCD calculation,Phys

    A. De Santis et al.,Inclusive semileptonic decays of the Ds meson: A first-principles lattice QCD calculation, Phys. Rev. D112(2025) 054503 [2504.06063]

    work page arXiv 2025

  13. [13]

    De Santis et al.,Inclusive Semileptonic Decays of theD s Meson: Lattice QCD Confronts Experiments,Phys

    A. De Santis et al.,Inclusive Semileptonic Decays of the Ds Meson: Lattice QCD Confronts Experiments,Phys. Rev. Lett.135(2025) 121901 [2504.06064]

    work page arXiv 2025

  14. [14]

    Bonanno, F

    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 (2023) 074515 [2305.17120]

    work page arXiv 2023

  15. [15]

    Bonanno, F

    C. Bonanno, F. D’Angelo, M. D’Elia, L. Maio and M. Naviglio,Sphaleron Rate of Nf=2+1 QCD,Phys. Rev. Lett.132(2024) 051903 [2308.01287]

    work page arXiv 2024

  16. [16]

    A novel Bayesian approach to spectral function reconstruction

    Y. Burnier and A. Rothkopf,Bayesian Approach to Spectral Function Reconstruction for Euclidean Quantum Field Theories,Phys. Rev. Lett.111(2013) 182003 [1307.6106]

    work page internal anchor Pith review Pith/arXiv arXiv 2013

  17. [17]

    From deep inelastic scattering to heavy-flavor semi-leptonic decays: Total rates into multi-hadron final states from lattice QCD

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

    work page internal anchor Pith review Pith/arXiv arXiv 2017

  18. [18]

    Del Debbio, T

    L. Del Debbio, T. Giani, J. Karpie, K. Orginos, A. Radyushkin and S. Zafeiropoulos,Neural-network analysis of Parton Distribution Functions from Ioffe-time pseudodistributions,JHEP02(2021) 138 [2010.03996]

    work page arXiv 2021

  19. [19]

    Del Debbio, T

    L. Del Debbio, T. Giani and M. Wilson,Bayesian approach to inverse problems: an application to NNPDF closure testing,Eur. Phys. J. C82(2022) 330 [2111.05787]. 14

    work page arXiv 2022

  20. [20]

    Medrano, H

    Y.C. Medrano, H. Dutrieux, J. Karpie, K. Orginos and S. Zafeiropoulos,Gaussian Processes for Inferring Parton Distributions,2510.21041

    work page arXiv

  21. [21]

    Barone, S

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

    work page arXiv 2023

  22. [22]

    Kellermann, Z

    R. Kellermann, Z. Hu, A. Barone, A. Elgaziari, S. Hashimoto, T. Kaneko et al.,Inclusive semileptonic decays from lattice QCD: Analysis of systematic effects, Phys. Rev. D112(2025) 014501 [2504.03358]

    work page arXiv 2025

  23. [23]

    Bergamaschi, W

    T. Bergamaschi, W.I. Jay and P.R. Oare,Hadronic structure, conformal maps, and analytic continuation, Phys. Rev. D108(2023) 074516 [2305.16190]

    work page arXiv 2023

  24. [24]

    Bruno, L

    M. Bruno, L. Giusti and M. Saccardi,Spectral densities from Euclidean lattice correlators via the Mellin transform,2407.04141

    work page arXiv

  25. [25]

    Patella and N

    A. Patella and N. Tantalo,Scattering Amplitudes from Euclidean Correlators: Haag-Ruelle theory and approximation formulae,2407.02069

    work page arXiv

  26. [26]

    The anomalous magnetic moment of the muon in the Standard Model: an update

    R. Aliberti et al.,The anomalous magnetic moment of the muon in the Standard Model: an update,Phys. Rept.1143(2025) 1 [2505.21476]

    work page internal anchor Pith review Pith/arXiv arXiv 2025