Pith. sign in

REVIEW 3 major objections 3 minor 1 cited by

Static-light meson spectroscopy with optimal distillation profiles

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

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

Pith's one-line read A profile-weighted variant of lattice distillation suppresses excited-state contamination in static-light and static-charm meson correlation functions, producing cleaner effective-mass plateaus and S, P1/2, and P3/2 mass splittings at two…

desk verdict Useful application of profile-weighted distillation to static-light mesons, but the abstract's claim that the 'optimal profiles' are what improves overlap is not actually isolated from the GEVP benefit, and one splitting in Table 5 looks anomalous. read the letter →

arxiv 2501.12863 v1 pith:XYJOYAS6 submitted 2025-01-22 hep-lat

Laura Struckmeier , Roman Höllwieser , Francesco Knechtli , Tomasz Korzec , Michael Peardon , Juan Andrés Urrea-Niño This is my paper · ORCID
classification hep-lat PACS 12.38.Gc
keywords static-lightmesonslatticeQCDoptimaldistillationprofilesheavyquarkeffectivetheoryexcited-statecontaminationmasssplittingsgeneralizedeigenvalueproblemB*pistates
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 aims to establish that "optimal distillation profiles"—a profile function that weights the eigenmodes of the 3D gauge-covariant Laplacian in lattice-QCD smearing—improve the overlap with energy eigenstates for static-light and static-charm mesons compared with standard distillation. On the $m_\pi \approx 420\,\mathrm{MeV}$ ensemble, the improved method gives an effective ground-state mass of $a m = 0.30594(81)$ with a plateau that begins near $t/a = 4$, while standard distillation at the same eigenvector count still shows excited-state contamination ($a m = 0.30773(87)$). The same technique is used to extract S, P1/2, and P3/2 spectra and mass splittings at pion masses of about 420 and 800 MeV. If correct, the method provides a cleaner route to heavy-light spectroscopy and tests of heavy-quark effective theory, with a static-light spectrum that responds more strongly to the light-quark mass than the static-charm spectrum.

What carries the argument

The central object is the optimal distillation profile $\rho_i(t) = \rho(\lambda_i(t))$, a function of the eigenvalues $\lambda_i$ of the 3D gauge-covariant Laplacian that modulates each eigenvector's weight in the smeared quark field; standard distillation is the special case of a step-function profile that keeps only the $N_v$ lowest modes. Using $N = 7$ Gaussian profiles builds an $N \times N$ correlation matrix whose energy levels are extracted with the generalized eigenvalue problem (GEVP), the machinery of [13,14]. The profile functions are what suppress excited-state contamination, and the local plus derivative operators, projected onto fermionic irreducible representations of the doubled cubic group, provide access to the S, $P_{1/2}$, and $P_{3/2}$ channels.

What would settle it

Extend the GEVP basis on the A1 ensemble with $B^*\pi$ interpolating operators (pion momenta zero and one unit) and recompute the static-light spectrum. If the level quoted as $1P_{1/2}$ at $277.9(6.9)$ MeV above the ground state, or the $1P_{3/2}$ level at $408(10)$ MeV, shifts by more than the quoted error, or a third level appears below the quoted $2S$ state, then the single-meson labels in Tables 4 and 5 are falsified.

Watch

Extended reading notes

Core claim

The central claim, stated in the abstract, is that the use of optimal profiles improves the overlap with the energy states compared to standard distillation. Concretely, Figure 1 compares the static-light ground-state effective mass on the A1 ensemble: with $N_v = 100$ Laplacian eigenvectors, standard distillation yields $a m = 0.30773(87)$ with visible excited-state contamination, while improved distillation with seven Gaussian profiles yields $a m = 0.30594(81)$ and a plateau from about $t/a = 4$. The paper then applies the improved technique to measure the static-light and static-charm spectra in the $G_1^+$, $G_1^-$, and $H^-$ irreps, identified with continuum $S$, $P_{1/2}$, and $P_{3/2}$ states, and reports mass splittings (Tables 4 and 5) at two pion masses. Non-interacting $B^*\pi$ energies are included for comparison, and the authors state that a precise investigation of $B^*\pi$ excited-state contamination would require adding $B^*\pi$ operators to the basis. It also concludes that static-light meson splittings depend more strongly on the pion mass than static-charm splittings, which largely cancel the heavy-quark mass dependence when the ground-state mass is subtracted.

Load-bearing premise

The load-bearing assumption is that each level extracted from the GEVP basis is an S, P1/2, or P3/2 single-meson state, even though the basis contains no $B^*\pi$ two-particle operators; if one of the extracted levels is actually a $B^*\pi$ scattering state, the corresponding splitting in Tables 4 and 5 mislabels a meson state.

Editorial extensions

If this is right

  • Using the same number of Laplacian eigenvectors, optimal profiles give a longer, cleaner plateau than standard distillation, so static-light and static-charm ground-state masses can be extracted with smaller systematic error.
  • The improved method makes higher radial and orbital excitations (2S, P1/2, P3/2) accessible enough that splittings such as $1P_{1/2}-1S$, $1P_{3/2}-1S$, and $2S-1S$ can be quoted at two pion masses.
  • The static-light splittings quoted on the $m_\pi \approx 420$ MeV ensemble, for example $1P_{3/2}-1S = 408(10)$ MeV, are the values the authors put forward for comparison with experimental B-meson splittings after extrapolation to the physical pion mass.
  • The heavier-pion A1h result shows $B^*\pi$ energies closer to the measured levels than the A1 result, so the paper's proposed next step of adding $B^*\pi$ operators is needed before any level can be labeled unambiguously as a meson excitation rather than a scattering state.

Reading between the lines

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

  • Beyond the paper, if the profile improvement persists at the physical pion mass and with more eigenmodes, optimal-profile distillation could reduce the need to construct explicit multi-particle interpolators for the lower heavy-light spectrum, since the profiles automatically emphasize the relevant wavefunction components.
  • Because the A1 and A1h ensembles differ in lattice spacing as well as pion mass, the attributed pion-mass dependence of the static-light splittings, such as $1P_{1/2}-1S$ changing from $277.9(6.9)$ to $369.9(5.6)$ MeV, is not isolated from discretization effects; a comparison at fixed lattice spacing would separate the two.
  • The planned $H$/$G_2$ comparison for the putative $5/2$ state doubles as a lattice-artifact test: if the two irreps that should form the continuum $5/2$ level do not align, the $P_{3/2}$ assignments carry residual symmetry-breaking contamination.
  • If the $B^*\pi$ check leaves the A1 excited levels unchanged, the first excited states in the $G_1^-$ and $H^-$ channels are predominantly single-meson radial excitations, which would make them usable inputs for heavy-meson chiral perturbation theory.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

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

Referee Report

3 major / 3 minor

Summary. The paper applies improved (profile-weighted) distillation, using Laplacian eigenmodes modulated by Gaussian profiles, to static-light and static-charm mesons on two N_f = 3+1 ensembles with pion masses of about 420 and 800 MeV. Operators are constructed in the relevant fermionic irreps of the doubled cubic group, and a 7-profile GEVP is used to extract energy levels. The central claim is that optimal profiles improve overlap with energy states compared to standard distillation, supported by effective-mass comparisons in Figure 1. The paper also reports S, P1/2 and P3/2 spectra and mass splittings in Tables 4 and 5, compares them with non-interacting B*pi thresholds, and quotes PDG values.

Significance. If the central claim holds, profile-weighted distillation is a practical improvement for heavy-light spectroscopy, and the two-ensemble dataset provides useful static-light and static-charm splittings. The paper deserves credit for the direct same-data comparison in Figure 1, for testing several N_v values for standard distillation, and for openly acknowledging that B*pi operators are missing from the basis. However, the improvement claim is not yet isolated from the GEVP, and the spectroscopic identifications are provisional because scattering-state operators are absent.

major comments (3)
  1. [Section 3, Figure 1] The improved-distillation effective mass is extracted from a 7x7 GEVP, while the standard-distillation curves shown in Figure 1 are not reported as GEVP results. The earlier plateau and reduced excited-state contamination of the improved curve could therefore be due to the variational benefit of the GEVP rather than to the profile weighting. To support the abstract and title claim, add a control in which standard hard-cutoff distillation operators are used in a GEVP with comparable basis size (for example N_v = 10, 30, 60, 100), or rephrase the claim as a combined effect of profile weighting plus GEVP.
  2. [Section 4, Tables 4 and 5] The levels labeled S, P1/2 and P3/2 are assumed to be single-meson states, but the operator basis contains no B*pi scattering operators. The manuscript itself states at the end of Section 4 that a precise investigation of this contamination requires including B*pi operators. Since the non-interacting B*pi energies lie close to some measured splittings, especially on A1h, the quoted splittings may misidentify scattering states as radial or orbital excitations. The statement that the A1 results "appear to be more likely radial excitations" is not a quantitative criterion; the tables should be qualified accordingly, or B*pi operators should be added before assigning these quantum numbers.
  3. [Section 4, Figures 3 and 4, and Section 5] The claimed dependence of the spectrum on the pion mass is not isolated, because the A1 and A1h ensembles differ in lattice spacing (0.05359 fm versus 0.0690 fm) as well as in pion mass. The observed differences in splittings therefore include discretization effects, and without a third ensemble or a continuum extrapolation the differences cannot be uniquely attributed to the light-quark mass. This limitation should be stated explicitly in the discussion of Tables 4 and 5.
minor comments (3)
  1. [Section 3] The seven Gaussian profiles used for the GEVP are not specified: neither the functional form nor the shape parameters are given in the text. Please provide the explicit definition or point to the relevant equations in refs. [10, 11].
  2. [Table 4] The PDG entry for the B_s(5840) splitting is garbled; it should read m_{B_s2^*(5840)^0} - m_{B_s^0}. Also, the comparison of the static-light 1P_{3/2} - 1S splitting with both B and B_s PDG values should be commented on, since these states have different light-quark content.
  3. [Figure 1] Panel (a) does not identify the standard-distillation curve in the legend; please add the symbol and N_v value used. The caption should also define the shaded bands and state the fit range used for the plateaus.

Circularity Check

0 steps flagged · score 2.0 of 10

No significant circularity: the improved-overlap claim is supported by a direct same-data comparison, and the self-citations supply the method rather than the empirical result.

full rationale

The central claim is that optimal distillation profiles improve overlap with the energy states compared with standard distillation. The paper supports this with a direct numerical comparison on the same A1 ensemble and the same number of Laplacian eigenvectors: improved distillation with seven Gaussian profiles plus a GEVP is compared against standard distillation at N_v = 100, 60, 30 and 10 (Fig. 1). This is an independent measurement, not a re-statement of an input. The improved-distillation technique is taken from the authors' prior works [10,11], so self-citations are present, but the load-bearing empirical comparison is performed in this paper and does not reduce to the citations. The spectra and mass splittings in Tables 4 and 5 are new lattice measurements; they are not derived from the PDG values with which they are compared, and the paper explicitly reports disagreement with those values. Two limitations are real but are not circularity: the comparison in Fig. 1 does not include a standard-distillation-plus-GEVP control, so the earlier plateau could be attributed partly to the variational GEVP rather than uniquely to the profile weighting; and the operator basis lacks B*pi operators, so some extracted levels could be misidentified scattering states. Both are experimental-design and interpretation concerns, not reductions of the result to its inputs. No fitted parameter is renamed as a prediction, and no equation defines the claimed improvement into existence. Accordingly, the circularity score is low, reflecting only the minor self-citation for method provenance.

Assumptions & free parameters 3 free parameters · 5 assumptions · 0 invented entities

The splittings in Tables 4 and 5 are direct GEVP measurements with no fitted physics parameters; the B*pi thresholds are computed from measured inputs through the lattice dispersion relation, not fitted. The free choices are analysis-level: the 7-profile Gaussian basis, the undocumented plateau fit windows, and the eigenvector counts N_v. The method itself is prior work by the same collaboration (refs [10] and [11] share authors with this paper), so the operator construction is self-referential, though the improvement claim is tested here against standard distillation on the same data. No new particles, forces, or dimensions are introduced; the 'optimal profiles' are weight functions on existing Laplacian eigenmodes, not new physical entities.

free parameters (3)
  • Profile basis: 7 Gaussian profiles and their shape parameters = not quoted
    The variational profiles rho(lambda) in Eq. (1) are built from 7 Gaussian basis functions with parameters chosen following refs [10, 11]. This is an operator-construction choice that affects GEVP quality, not a physics constant.
  • Effective-mass plateau fit ranges = not specified
    The quoted masses come from a 'correlated linear fit to the logarithm of the effective mass' (Figure 1 caption, Section 3); the fit windows are not documented, and different windows would shift the central values and errors.
  • Number of Laplacian eigenvectors N_v = 100 (light), 200 (charm), per Table 2
    Convergence parameter for the distillation subspace, chosen separately per ensemble and quark mass. Figure 1 shows how the standard method depends on this choice.
assumptions (5)
  • domain assumption The N_f = 3+1 ensembles of ref. [15] with the parameters of Table 2 correctly represent QCD at the quoted pion masses and lattice spacings, with the scale a set as in ref. [9].
    All physical units in Tables 4 and 5 come from the scale setting of the authors' earlier paper [9], which is not re-derived here.
  • domain assumption Heavy quark spin symmetry holds in the static limit, so the static quark is a color source without spin and the states are classified by the double cover of the cubic group OD_h.
    Used in Section 2 for the operator construction and in Section 4 for the B*pi quantum-number combinations.
  • domain assumption The improved-distillation profile construction of refs [10, 11] works as stated for static-light correlation functions, including the derivative-operator generalization in Section 3.
    The paper relies on the group's own method papers for the profile optimization and does not re-derive the algorithm in this proceedings text.
  • domain assumption The generalized eigenvalue problem built from 7 Gaussian profiles reliably extracts the lowest energy levels from the correlation matrix.
    Section 3: 'By choosing N different profiles, in this work N = 7 Gaussian profiles, a GEVP can be solved to extract the energy eigenstates.'
  • standard math The pion energies follow the relativistic lattice dispersion relation, and the subduction of pion irreps given in Table 3 is complete.
    Used in Section 4 to build the non-interacting B*pi threshold energies for the contamination comparison.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Static-light meson spectroscopy with optimal distillation profiles." pith.science (2026). https://pith.science/paper/XYJOYAS6

@misc{pith2026250112863,
  author       = {Pith},
  title        = {Pith review of: Static-light meson spectroscopy with optimal distillation profiles},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/XYJOYAS6}},
  note         = {Machine review of arXiv:2501.12863}
}
abstract

The spectrum of static-light and static-charm mesons is studied using optimized distillation in two different $N_{\rm{f}} = 3 + 1$ QCD ensembles with pion masses of $m_{\pi} \approx 800 \, \text{MeV}$ and $m_{\pi} \approx 420 \,\text{MeV}$ and a heavy (charm) quark. Local and derivative-based operators are used to access states of multiple quantum numbers. The use of optimal profiles is shown to improve the overlap with the energy states compared to standard distillation.

Figures

Figures reproduced from arXiv: 2501.12863 by the authors.

Figure 1
Figure 1. Comparison of the effective ground state mass between standard distillation (colored x’s) and improved distillation (black dots) obtained using 𝑁𝑣 = 100 eigenvectors for the ground state of the static￾light meson on the A1 ensemble. The mass plateaus and their errors (light shaded bands) are obtained by a correlated linear fit to the logarithm of the effective mass. using standard distillation with fewer eigenvector… view at source ↗
Figure 2
Figure 2. Optimal meson distillation profiles ˜𝑓 (𝜆) for the ground state (black) and first excited state (red) of the static-light meson. static-light spectrum, the pion has to have momentum in certain cases. The subduction of pion representations for the smallest lattice momenta | ®𝑝| is given in table 3 [5]. 𝑝® Irreducible content (0, 0, 0) 𝐴 − 1 (1, 0, 0) 𝐴 − 1 ⊕ 𝐸 − ⊕ 𝑇 + 1 (1, 1, 0) 𝐴 − 1 ⊕ 𝐸 − ⊕ 𝑇 + 1 ⊕ 𝑇 + 2 ⊕ 𝑇 − 2 … view at source ↗
Figure 3
Figure 3. Static-light meson spectrum with accessible radial and orbital excitations, obtained using improved distillation and non-interacting 𝐵 ∗𝜋-states (magenta). S P1/2 P3/2 0 200 400 600 800 1000 1200 1400 1600 E E ( 0 ) Q c [ M e V ] (a) A1 ensemble. S P1/2 P3/2 0 200 400 600 800 1000 1200 1400 1600 E E ( 0 ) Q c [ M e V ] (b) A1h ensemble [PITH_FULL_IMAGE:figures/full_fig_p006_3.png] view at source ↗
Figures from the paper (1 more)
Figure 4
Figure 4. Figure 4: Static-charm meson spectrum with accessible radial and orbital excitations, obtained using improved distillation. The magenta-colored bars in figure 3 correspond to these energies in all given channels. Although the pion mass is heavier in the A1h ensemble, the 𝐵 ∗𝜋 st…

Discussion (0). Continue with ORCID to comment.

Forward citations

Cited by 1 Pith paper

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

  1. Hybrid static potentials and gluelumps on $N_f=3+1$ ensembles

    hep-lat 2025-01 conditional novelty 6.0 of 10

    New lattice QCD measurements of hybrid static potentials, static-light thresholds, and gluelump masses on N_f=3+1 ensembles with pions near 420 MeV, using Laplace trial states.

Reference graph

Works this paper leans on

17 extracted references · 6 canonical work pages · cited by 1 Pith paper

  1. [1]

    Introduction to Non-perturbative Heavy Quark Effective Theory

    R. Sommer,Introduction to Non-perturbative Heavy Quark Effective Theory, 1008.0710

    work page Pith review arXiv

  2. [2]

    Morte, A

    M.D. Morte, A. Shindler and R. Sommer,On lattice actions for static quarks,JHEP 2005 (2005) 51 [hep-lat/0506008]

    arXiv 2005

  3. [3]

    O. Bar, A. Broll and R. Sommer,𝐵𝜋 excited-state contamination in lattice calculations of B-meson correlation functions,Eur. Phys. J. C83(2023) 757 [2306.02703]

    work page Pith review arXiv 2023

  4. [4]

    Michael and J

    C. Michael and J. Peisa,Maximal variance reduction for stochastic propagators with applications to the static quark spectrum, Phys. Rev. D58 (1998) 034506 [hep-lat/9802015]

    arXiv 1998

  5. [5]

    Radial and orbital excitations of static-light mesons

    J. Foley, A. O’Cais, M. Peardon and S.M. Ryan,Radial and orbital excitations of static-light mesons,Phys. Rev. D75 (2007) 094503 [hep-lat/0702010]

    work page Pith review arXiv 2007

  6. [6]

    Altmann and P

    S. Altmann and P. Herzig,Point-group Theory Tables, Oxford science publications, Clarendon Press (1994)

    work page 1994

  7. [7]

    Peskin and D.V

    M.E. Peskin and D.V. Schroeder,An Introduction to quantum field theory, Addison-Wesley, Reading, USA (1995), 10.1201/9780429503559

    doi:10.1201/9780429503559 1995

  8. [8]

    The static-light meson spectrum from twisted mass lattice QCD

    K. Jansen, C. Michael, A. Shindler and M. Wagner,The Static-light meson spectrum from twisted mass lattice QCD, JHEP 12(2008) 058 [0810.1843]. 8 Static-light meson spectroscopy with optimal distillation profiles Laura Struckmeier

    work page Pith review arXiv 2008

Show all 17 references
  1. [9]

    Höllwieser, F

    R. Höllwieser, F. Knechtli and T. Korzec,Scale setting for𝑁𝑓 = 3+ 1 QCD,Eur. Phys. J. C 80 (2020) 349 [1907.04309v2]

    2020 arXiv

  2. [10]

    Knechtli, T

    F. Knechtli, T. Korzec, M. Peardon and J.A. Urrea-Niño,Optimizing creation operators for charmonium spectroscopy on the lattice, Phys. Rev. D106 (2022) 034501 [2205.11564]

    2022 arXiv

  3. [11]

    Höllwieser, F

    R. Höllwieser, F. Knechtli, T. Korzec, M. Peardon and J.A. Urrea-Niño,Constructing static quark-antiquark creation operators from Laplacian eigenmodes, Phys. Rev. D107 (2023) 034511 [2212.08485]

    2023 arXiv

  4. [12]

    Peardon, J

    M. Peardon, J. Bulava, J. Foley, C. Morningstar, J. Dudek, R.G. Edwards et al.,A Novel quark-field creation operator construction for hadronic physics in lattice QCD,Phys. Rev. D 80 (2009) 054506 [0905.2160]

    2009 arXiv

  5. [13]

    M.LuscherandU.Wolff, HowtoCalculatetheElasticScatteringMatrixinTwo-dimensional Quantum Field Theories by Numerical Simulation,Nucl. Phys. B339 (1990) 222

    1990

  6. [14]

    Blossier, M.D

    B. Blossier, M.D. Morte, G. von Hippel, T. Mendes and R. Sommer,On the generalized eigenvalue method for energies and matrix elements in lattice field theory, JHEP 2009 (2009) 94 [0902.1265]

    2009 arXiv

  7. [15]

    Fritzsch, R

    P. Fritzsch, R. Sommer, F. Stollenwerk and U. Wolff,Symanzik improvement with dynamical charm: a 3+1 scheme for Wilson quarks, JHEP 6 (2018) 25 [1805.01661]

    2018 arXiv

  8. [16]

    Joswig, S

    F. Joswig, S. Kuberski, J.T. Kuhlmann and J. Neuendorf,pyerrors: A python framework for error analysis of Monte Carlo data,Comput. Phys. Commun.288(2023) 108750 [2209.14371]

    2023 arXiv

  9. [17]

    Navas et al.,Review of particle physics, Phys

    S. Navas et al.,Review of particle physics, Phys. Rev. D110 (2024) 3 030001. 9

    2024

Pith tools

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