Pith. sign in

REVIEW 2 major objections 6 minor 1 cited by

NRQCD Bottomonium at non-zero temperature using time-derivative moments

T0 review · 2 major / 6 minor · reviewed 2026-08-09 · deepseek-v4-flash

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

Pith's one-line read Time-derivative moments of the lattice correlator, taken under a Gaussian spectral-function assumption, give the Upsilon and chi_b1 masses and thermal widths without reconstructing the full spectral function.

desk verdict A promising moment-based method for NRQCD quarkonia, but the fit forms drop terms that are not negligible at high T, so the reported mass and width shifts are not yet reliable. read the letter →

arxiv 2502.03951 v2 pith:6TXYVHDN submitted 2025-02-06 hep-lat

Rachel Horohan D'arcy , Gert Aarts , Chris Allton , Ryan Bignell , Timothy Burns , Benjamin Jaeger , Seyong Kim , Maria Paola Lombardo
show 3 more authors
Sinéad Ryan Antonio Smecca Jon-Ivar Skullerud
This is my paper · ORCID
classification hep-lat
keywords bottomoniumNRQCDlatticeQCDtime-derivativemomentsthermalwidthspectralfunctionquark-gluonplasmaheavyquarkonia
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 claims a new method for extracting the mass and thermal width of heavy quarkonium states from lattice QCD correlators without solving the ill-posed inverse problem of spectral reconstruction. The method assumes the spectral function is a sum of Gaussians, then reads the first and second time-derivative moments of the NRQCD correlator, which give the mass and width directly. On the Upsilon and the chi_b1 it reproduces the experimental zero-temperature masses, shows the thermal width rising with temperature, and — after a zero-temperature temporal-truncation control removes a kinematic effect — finds the mass decreasing with temperature. The result matters because quarkonium melting is a key diagnostic of the quark-gluon plasma produced in heavy-ion collisions.

What carries the argument

The central objects are the time-derivative moments of the Euclidean correlator: the effective mass $M_L(\tau)=\partial\log G(\tau)/\partial\tau$ and the squared thermal width $\Gamma^2_L(\tau)=\partial^2\log G(\tau)/\partial\tau^2$, together with their discrete-derivative partners $M_D$ and $\Gamma^2_D$ used to estimate systematic errors. They do the work because a Gaussian spectral function centred at $m_i$ with width $\Gamma_i$ maps the correlator to $\exp(-m_i \tau+\Gamma_i^2\tau^2/2)$; the first log-derivative then reads off the mass and the curvature of $\log G$ reads off the width, sidestepping the inverse problem entirely. The subsequent fit keeps one excited state, assumes the ground state is well separated, and uses the narrow-width condition $\Gamma^2\tau\ll m_0$ to extract $m_0$ and $\Gamma_0$ at each temperature.

What would settle it

Take the same lattice ensembles and run the moment method on synthetic correlators generated from a known non-Gaussian spectral function, for instance a skewed or two-peak shape, and from a $\chi_{b1}$-like large width; if the recovered mass and width deviate from the inputs beyond the quoted errors, the Gaussian Ansatz is biased. A second check is to compare the $\chi_{b1}$ width extracted here with the maximum-likelihood spectral widths obtained from the same correlators in refs. [8] and [16]; agreement would support the Ansatz, while disagreement would settle against it.

Watch

Extended reading notes

Core claim

On FASTSUM's Generation 2L anisotropic lattice ensembles, the authors replace the ill-posed reconstruction of the thermal spectral function $\rho(\omega;T)$ by assuming it is a finite sum of Gaussians centred at each state's mass with width set by its thermal broadening. Under that Ansatz the NRQCD correlator becomes a sum of exponentials $\exp(-m_i \tau + \Gamma_i^2 \tau^2/2)$, so the log-derivatives of the correlator give the ground-state mass $M_L(\tau)=\partial\log G/\partial\tau$ and squared width $\Gamma^2_L(\tau)=\partial^2\log G/\partial\tau^2$ without any full inversion. Fitting these moments with one excited state and the narrow-width approximation yields $M(\Upsilon)=9455(10)$ MeV and $M(\chi_{b1})=9965(79)$ MeV at zero temperature, consistent with the Particle Data Group values, and a thermal width that grows with temperature for both states. An apparent rise in the $\Upsilon$ mass with temperature is shown by the truncated zero-temperature ($T_0$) control to be a short-temporal-range kinematic effect; with that control the mass decreases with temperature.

Load-bearing premise

The load-bearing premise is that the true spectral function is a finite sum of well-separated Gaussians and that the ground-state width obeys the narrow-width condition $\Gamma^2\tau\ll m_0$; if the real spectral shape is non-Gaussian or the $\chi_{b1}$ width is too large for that limit, the extracted mass and width are biased, and the authors single out the $\chi_{b1}$ case as an outstanding question.

Editorial extensions

If this is right

  • The method extracts ground-state masses and thermal widths for $\Upsilon$ and $\chi_{b1}$ using only derivatives of the correlator, avoiding the ill-posed inverse problem that dominates spectral-function studies.
  • Zero-temperature masses agree with experiment, so the Gaussian-moment approach reproduces established bottomonium spectroscopy from the same ensembles.
  • The thermal width $\Gamma$ increases with temperature for both states, consistent with quarkonium suppression in the quark-gluon plasma.
  • The apparent upward shift of the $\Upsilon$ mass with temperature is a kinematic truncation effect, not a thermal one; after the $T_0$ control the mass falls with temperature.
  • Relaxing the single-excited-state and narrow-width assumptions, by fitting Eqs. (15) and (16) with $N>1$, is a direct route to reduce systematic uncertainties in future applications.

Reading between the lines

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

  • The same derivative-moment machinery should transfer to charmonium or other heavy-quark bound states, where the Gaussian Ansatz's bias could be quantified by comparing with full spectral reconstructions on identical ensembles.
  • A synthetic-data test, generating correlators from known non-Gaussian spectral shapes and running the method on them, would map exactly where the Gaussian-plus-narrow-width assumption distorts the extracted mass and width.
  • The $T_0$ truncation control is a transferable subtraction scheme: any lattice spectral study whose high-temperature correlators suffer from short temporal extent could apply the same zero-temperature truncation to isolate genuine thermal effects.
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

2 major / 6 minor

Summary. The manuscript proposes a time-derivative-moment method for NRQCD bottomonium correlators: assuming the spectral function is a sum of Gaussians, the authors define logarithmic first and second derivatives of the correlator and fit them to extract the ground-state mass and thermal width of the Upsilon and chi_b1 for temperatures from 47 to 380 MeV. They report zero-temperature masses consistent with experiment after adding the NRQCD energy shift, a thermal broadening of both states, and a mass decrease once a zero-temperature control ('T0 analysis') is applied to separate short-time kinematic effects from genuine thermal effects. The central technical claim is that Eqs. (17)-(18) follow from the Gaussian ansatz under an N=1 narrow-width simplification.

Significance. The method is attractive because it avoids full spectral reconstruction and is computationally cheap, and the paper includes useful control checks: bootstrap errors, comparison of two derivative discretizations (M_L/M_D and Gamma_L/Gamma_D), and a T0 analysis that separates kinematic from thermal effects. If the derivation of the fit forms were consistent, the approach would be a valuable complementary tool to MEM and Bayesian reconstructions for quarkonium at nonzero temperature. However, the manuscript contains a load-bearing gap: the fit forms used for the reported masses and widths do not follow from the stated Gaussian ansatz. The concern is confirmed by comparing Eqs. (15)-(18). The paper does not provide machine-checked code or a full model-systematic error budget for the dropped terms, so the numerical results should be regarded as preliminary until the fit forms are justified or replaced.

major comments (2)
  1. [Sec. 3, Eqs. (15) and (17)] The fit form (17) does not follow from the Gaussian ansatz. For N=1, Eq. (15) gives M_L(tau) = -m0 + Gamma0^2 tau + (A1/A0)(-Delta m + DeltaGamma^2 tau) exp(-Delta m tau + DeltaGamma^2 tau^2/2). Equation (17) instead fits M(tau) = m0 + A exp(-Delta m tau), dropping the Gamma0^2 tau term and the tau-dependent prefactor. The stated narrow-width condition Gamma^2 tau << m0 only makes Gamma0^2 tau small relative to m0 ~ 9.5 GeV; it does not make it small relative to the quoted mass uncertainties (about 10 MeV) or the thermal mass shifts reported in Fig. 2 (tens of MeV). For the chi_b1, with Gamma of order 150-400 MeV as shown in Fig. 4 (right), Gamma^2 tau is of order 100 MeV over the fitted window, which is comparable to the magnitude of the reported effects. The linear large-tau trend in the chi_b1 panel of Fig. 1 is the signature of the neglected term. The authors should fit the full expression (15) or quantify the bias from omitting it; otherwise the extracted masses are biased by an uncontrolled approximation.
  2. [Sec. 3, Eqs. (16) and (18)] The width fit is internally inconsistent. With N=1, Eq. (16) gives Gamma^2_L(tau) = Gamma0^2 + (A1/A0)[(-Delta m + DeltaGamma^2 tau)^2 + DeltaGamma^2] exp(-Delta m tau + DeltaGamma^2 tau^2/2). Equation (18) retains DeltaGamma^2 in the exponent but replaces the entire prefactor by a constant A. This is justified only if DeltaGamma^2 = 0, in which case the exponent should be -Delta m tau without the DeltaGamma^2 tau^2/2 term. If DeltaGamma^2 is retained because it is needed to represent the data, the tau-dependence of the prefactor is generically as important as the exponential factor over the fitted range. The Gamma(T) values in Fig. 4 are therefore obtained from a functional form that is not a consequence of the model stated in Eq. (2), and the reported thermal widths inherit an unquantified systematic error.
minor comments (6)
  1. [Sec. 3, Eq. (5)] The printed formula contains a malformed integral and a stray minus sign in the numerator; please correct the typesetting so that the definition of M_D is unambiguous.
  2. [Sec. 4] The sentence describing the bootstrap and systematic errors is duplicated with slightly different wording; please consolidate it and specify which systematic estimate is shown in each figure.
  3. [Fig. 3 (right)] The legend omits the T=47 MeV entry that appears for the Upsilon; please state whether chi_b1 data at this temperature are unavailable and why.
  4. [Fig. 4 (right)] The text says systematic errors from Gamma_D are estimated, but only Gamma_L and Gamma_L(T0) are plotted for chi_b1; please clarify whether Gamma_D was omitted because it is too noisy, and show it or state the reason.
  5. [Sec. 4] The zero-temperature Upsilon agreement with experiment is not an independent validation because the additive shift E0 is fixed by the Upsilon(1S) mass; please state this explicitly and present the chi_b1 mass as the genuine prediction.
  6. [Sec. 4] For reproducibility, please report the fit ranges, the number of points, and the fit quality (chi^2 or comparable) for the fits to Eqs. (17) and (18) at each temperature.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the extraction fits explicit Gaussian-model parameters to lattice correlator data; the one calibrating input (NRQCD additive shift set by the Upsilon mass) is disclosed and does not drive the thermal or chi_b1 results.

full rationale

The central claim is a parameter extraction rather than a derivation from an external principle: the paper explicitly assumes a sum-of-Gaussians spectral function (Eq. (2), 'We approximate the spectral function... by a finite sum of Gaussian functions'), computes the resulting correlator and its log-derivatives (Eqs. (12)-(16)), and then fits lattice data with the forms (17)-(18). The fit parameters m0 and Gamma0 are independent outputs determined by the correlator data; they are not defined to equal any input quantity, and the Gaussian ansatz is disclosed as an approximation rather than imported as a theorem. The citations to [8,9] (one of which is a prior FASTSUM paper) support the choice of ansatz but do not supply a forced or uniqueness result; the in-preparation review [7] is not load-bearing. The paper's own Section 5 flags the limit Gamma^2 -> 0 for chi_b1 as an outstanding question and defers fits to the full Eqs. (15)-(16), which is a validity caveat, not a hidden circular step. The only point where an input controls a reported number is the additive NRQCD shift E0 = 7463 MeV, stated to be 'set by the Upsilon(1S) mass'; consequently the quoted zero-temperature Upsilon agreement is a calibration consistency check rather than an independent prediction. That calibration does not determine the chi_b1 mass, the thermal widths, or the temperature dependence, and the T0 (truncated zero-temperature correlator) control provides an independent cross-check of kinematic vs thermal effects. Although the approximations leading to (17)-(18) may bias the fit if the narrow-width assumption fails, that is a model-misspecification issue, not a circular reduction of the result to its inputs. I therefore find no load-bearing circularity.

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

The central extraction rests on a Gaussian model for the spectral function and a set of truncation approximations. The only calibrated input is the NRQCD shift E0; all other numbers come from fits to the correlator. No new physical entities are introduced.

free parameters (2)
  • Additive NRQCD energy shift E0 = 7463 MeV
    Set to reproduce the Upsilon(1S) mass at zero temperature, so the Upsilon mass comparison is not independent.
  • Fit parameters m0, Gamma0, A, Delta m, DeltaGamma^2 in Eqs (17)-(18) = not tabulated in paper
    For each state and temperature, the mass and width are extracted from fits with these parameters. Values appear only in figures; no tables or fit ranges are given.
assumptions (5)
  • domain assumption NRQCD expansion is valid for bottomonium
    The calculation uses the NRQCD kernel K=e^{-omega tau} throughout (Section 1).
  • ad hoc to paper Spectral function is a finite sum of Gaussians
    Eq (2) is the central modeling assumption; it is not derived from QCD, and the authors flag its limitations for chi_b1 in Section 5.
  • ad hoc to paper Only one excited state contributes (N=1)
    Stated in Section 3: 'we will be assuming only a single excited state contributes, i.e., N=1'; this relies on ground-state dominance at large tau.
  • ad hoc to paper Narrow-width approximation Gamma^2 tau << m0
    Used to simplify the mass and width fit functions; the authors question it for chi_b1 in Section 5.
  • domain assumption FASTSUM Gen2L ensembles represent the QCD thermal medium
    The results depend on these specific ensembles with m_pi=240 MeV and T_c=167 MeV (Section 2).

how reviews work

0 comments
Cite this review

Pith. "Pith review of NRQCD Bottomonium at non-zero temperature using time-derivative moments." pith.science (2026). https://pith.science/paper/6TXYVHDN

@misc{pith2026250203951,
  author       = {Pith},
  title        = {Pith review of: NRQCD Bottomonium at non-zero temperature using time-derivative moments},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/6TXYVHDN}},
  note         = {Machine review of arXiv:2502.03951}
}
abstract

A well-known challenge for the lattice community is calculating the spectral function from the Euclidean correlator. We have approximated the spectral function and derived the mass and thermal width of particles through the time derivatives of the lattice correlator moments. We have focused on extracting the properties of bottomonium states, specifically $\Upsilon$ and $\chi_{b1}$. We will give an overview of the time-derivative moments approach and present results for the temperature dependence of the mass and width of both bottomonium states. The zero temperature results are consistent with experimental values, while results at higher temperatures are similar to those obtained using other methods.

Figures

Figures reproduced from arXiv: 2502.03951 by the authors.

Figure 1
Figure 1. Left: 𝑀𝐿 (𝜏) of the Υ at each temperature. Right: 𝑀𝐿 (𝜏) of the 𝜒𝑏1 at each temperature. 100 150 200 250 300 350 T (MeV) 9490 9500 9510 9520 9530 9540 M(T) (M e V) MD ML ML(T0) 50 100 150 200 250 300 350 T (MeV) 9850 9900 9950 10000 10050 10100 M(T) (M e V) MD ML ML(T0) [PITH_FULL_IMAGE:figures/full_fig_p006_1.png] view at source ↗
Figure 2
Figure 2. Left: 𝑀𝐿 (blue) from the fit for the Υ at each temperature and the corresponding 𝑇0 (red) result. Right: 𝑀𝐿 (blue) from the fit for the 𝜒𝑏1 at each temperature and the corresponding 𝑇0 (red) result. although there is a lot of noise in the data at large 𝜏. We fit Γ 2 (𝜏) for each temperature to Equation (18), shown in [PITH_FULL_IMAGE:figures/full_fig_p006_2.png] view at source ↗
Figure 3
Figure 3. Left: Γ 2 𝐿 (𝜏) of the Υ at each temperature. Right: Γ 2 𝐿 (𝜏) of the 𝜒𝑏1 at each temperature. 50 100 150 200 250 300 350 T (MeV) 50 100 150 200 250 300 (M e V) L L(T0) D 100 150 200 250 300 T (MeV) 100 150 200 250 300 350 400 450 (M e V) L L(T0) D [PITH_FULL_IMAGE:figures/full_fig_p007_3.png] view at source ↗
Figures from the paper (1 more)
Figure 4
Figure 4. Figure 4: Left: Γ𝐿 (blue) and Γ𝐷 (green) from the fit for the Υ at each temperature and the corresponding 𝑇0 (red) result. Right: Γ𝐿 (blue) from the fit for the 𝜒𝑏1 at each temperature and the corresponding 𝑇0 (red) result. Points offset for clarity due to the limited number of …

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. Finite temperature hadronic spectral properties

    hep-lat 2025-05 conditional novelty 5.0 of 10

    Doubly charmed baryons stay stable and the bottomonium ground state loses 30 to 40 MeV as quark-gluon plasma temperature rises, according to new lattice QCD ensembles.

Reference graph

Works this paper leans on

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

  1. [1]

    Brambilla et al.,Heavy Quarkonium: Progress, Puzzles, and Opportunities, Eur

    N. Brambilla et al.,Heavy Quarkonium: Progress, Puzzles, and Opportunities, Eur. Phys. J. C 71(2011) 1534 [1010.5827]

    arXiv 2011

  2. [2]

    G.P.Lepage,L.Magnea,C.Nakhleh,U.MagneaandK.Hornbostel, Improvednonrelativistic QCD for heavy quark physics, Phys. Rev. D46(1992) 4052 [hep-lat/9205007]

    arXiv 1992

  3. [3]

    Brambilla, A

    N. Brambilla, A. Pineda, J. Soto and A. Vairo,Effective Field Theories for Heavy Quarkonium,Rev. Mod. Phys.77(2005) 1423 [hep-ph/0410047]

    arXiv 2005

  4. [4]

    Aarts et al.,Properties of the QCD thermal transition with Nf=2+1 flavors of Wilson quark, Phys

    G. Aarts et al.,Properties of the QCD thermal transition with Nf=2+1 flavors of Wilson quark, Phys. Rev. D105 (2022) 034504 [2007.04188]

    arXiv 2022

  5. [5]

    Aarts, C

    G. Aarts, C. Allton, R. Bignell, T.J. Burns, S.C. García-Mascaraque, S. Hands et al.,Open charm mesons at nonzero temperature: results in the hadronic phase from lattice QCD, 2209.14681

    arXiv

  6. [6]

    A comparison of spectral reconstruction methods applied to non-zero temperature NRQCD meson correlation functions

    T. Spriggs et al.,A comparison of spectral reconstruction methods applied to non-zero temperature NRQCD meson correlation functions,EPJ Web Conf.258 (2022) 05011 [2112.04201]

    work page Pith review arXiv 2022

  7. [7]

    Aarts, C

    G. Aarts, C. Allton, R. Bignell, T.J. Burns, B. Jäger, R. Horohan Darcy et al.,Systematic study of the temperature dependence of the mass and width of quarkonia in lattice nonrelativistic QCD,(in preparation)

    work page

  8. [8]

    Spriggs et al.,Bottomonium spectral widths at nonzero temperature using maximum likelihood, PoS LATTICE2021(2022) 077 [2112.01599]

    T. Spriggs et al.,Bottomonium spectral widths at nonzero temperature using maximum likelihood, PoS LATTICE2021(2022) 077 [2112.01599]

    arXiv 2022

Show all 17 references
  1. [9]

    Larsen, S

    R. Larsen, S. Meinel, S. Mukherjee and P. Petreczky,Thermal broadening of bottomonia: Lattice nonrelativistic QCD with extended operators, Phys. Rev. D100 (2019) 074506 [1908.08437]

    2019 arXiv

  2. [10]

    Baer,findiff software package, 2018

    M. Baer,findiff software package, 2018

    2018

  3. [11]

    Aarts, C

    G. Aarts, C. Allton, T. Harris, S. Kim, M.P. Lombardo, S.M. Ryan et al.,The bottomonium spectrum at finite temperature from N𝑓 = 2 + 1 lattice QCD, JHEP 07(2014) 097 [1402.6210]

    2014 arXiv

  4. [12]

    Offler, G

    S. Offler, G. Aarts, C. Allton, J. Glesaaen, B. Jäger, S. Kim et al.,News from bottomonium spectral functions in thermal QCD, PoS LATTICE2019(2019) 076 [1912.12900]

    2019 arXiv

  5. [13]

    Particle Data Groupcollaboration, Review of particle physics,Phys. Rev. D110 (2024) 030001. 8 NRQCD Bottomonium at non-zero temperature using time-derivative momentsRachel Horohan D’arcy

    2024

  6. [14]

    Aarts, C

    G. Aarts, C. Allton, S. Kim, M.P. Lombardo, M.B. Oktay, S.M. Ryan et al.,What happens to theΥ and𝜂𝑏 in the quark-gluon plasma? Bottomonium spectral functions from lattice QCD, JHEP 11 (2011) 103 [1109.4496]

    2011 arXiv

  7. [15]

    Kelly, A

    A. Kelly, A. Rothkopf and J.-I. Skullerud,Bayesian study of relativistic open and hidden charm in anisotropic lattice QCD,Phys. Rev. D97 (2018) 114509 [1802.00667]

    2018 arXiv

  8. [16]

    S. Kim, P. Petreczky and A. Rothkopf,Quarkonium in-medium properties from realistic lattice NRQCD, JHEP 11(2018) 088 [1808.08781]

    2018 arXiv

  9. [17]

    Bignell, G

    R. Bignell, G. Aarts, C. Allton, M.N. Anwar, T.J. Burns, R. Horohan D’arcy et al., Anisotropic excited bottomonia from a basis of smeared operators, PoS LATTICE2024 (2025) 202. 9

    2025

Pith tools

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