Charmonium properties at high temperatures from lattice QCD
Pith reviewed 2026-05-20 03:19 UTC · model grok-4.3
The pith
Lattice QCD finds charmonium states survive below open charm threshold up to 305 MeV with growing thermal widths.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Our lattice QCD results are consistent with the existence of all charmonium states below the open charm threshold in this temperature region. However, charmonium states acquire sizable thermal width, which increases with increasing temperature. The size of the thermal width follows the hierarchy of charmonium sizes, i.e. the smaller ground state charmonium has a smaller thermal width than the larger excited charmonia.
What carries the argument
Correlation functions of extended meson operators on lattices with HISQ dynamical quarks and Wilson clover valence charm quarks, from which thermal widths and state existence are extracted.
If this is right
- All charmonium states below the open charm threshold remain present between 153 MeV and 305 MeV.
- Thermal widths become sizable and continue to grow as temperature rises.
- The thermal width scales with state size, smaller for the ground state than for excited states.
- This pattern holds for states that stay below the open charm threshold.
Where Pith is reading between the lines
- The size-dependent widths suggest screening in the medium affects larger states first.
- These results could guide hydrodynamic models of heavy-ion collisions by providing temperature-dependent widths.
- Similar calculations for bottomonium would test whether the hierarchy is universal across heavy quark flavors.
Load-bearing premise
The correlation functions from the chosen extended operators and lattice actions can be interpreted as directly indicating the continued existence of distinct states and the magnitude of their thermal widths without dominant contamination from lattice artifacts or excited-state mixing.
What would settle it
Higher-statistics lattice runs or calculations with alternate operators at the same temperatures showing zero thermal width or sudden disappearance of states below the open charm threshold.
read the original abstract
We study charmonium properties at non-zero temperature in the temperature range 153 MeV $<T<$ 305 MeV using lattice QCD. We use HISQ action for dynamical quarks and Wilson clover action for valence charm quarks and calculate the correlation function of extended meson operators. Our lattice QCD results are consistent with the existence of all charmonium states below the open charm threshold in this temperature region. However, charmonium states acquire sizable thermal width, which increases with increasing temperature. The size of the thermal width follows the hierarchy of charmonium sizes, i.e. the smaller ground state charmonium has a smaller thermal width than the larger excited charmonia.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper reports lattice QCD calculations of charmonium correlation functions in the temperature range 153 MeV < T < 305 MeV, employing the HISQ action for dynamical quarks and the Wilson clover action for valence charm quarks together with extended meson operators. The central claims are that the results remain consistent with the existence of all charmonium states below the open-charm threshold while these states develop sizable thermal widths that grow with temperature and obey the hierarchy of the states' spatial sizes.
Significance. If the thermal-width extractions are shown to be robust against systematics, the results would supply useful non-perturbative information on the in-medium behavior of heavy quarkonia, with potential relevance to quarkonium suppression phenomenology in heavy-ion collisions. The reported size-dependent width hierarchy would be a distinctive observation worth further theoretical comparison.
major comments (2)
- Abstract: the claim that the lattice results are 'consistent with the existence of all charmonium states below the open charm threshold' is not accompanied by any description of the fitting procedure, number of exponentials, fit-range stability, or spectral reconstruction method used to extract the thermal widths, leaving the central conclusion without visible support.
- Methods (correlation-function analysis): at T = 305 MeV the temporal extent Nt is necessarily small; the manuscript gives no indication that excited-state contamination or operator-dependent artifacts were controlled through smearing-parameter variation, multi-state fits, or cross-checks with point sources, which directly affects the reliability of the reported widths and size hierarchy.
minor comments (2)
- Specify the lattice spacings, volumes, and number of configurations employed so that discretization and finite-volume effects can be assessed.
- Clarify how the 'extended meson operators' were constructed and whether their smearing parameters were varied as a systematic check.
Simulated Author's Rebuttal
We thank the referee for the careful reading of our manuscript and the constructive comments. We address each major point below and have made revisions to improve clarity and robustness of the presentation.
read point-by-point responses
-
Referee: Abstract: the claim that the lattice results are 'consistent with the existence of all charmonium states below the open charm threshold' is not accompanied by any description of the fitting procedure, number of exponentials, fit-range stability, or spectral reconstruction method used to extract the thermal widths, leaving the central conclusion without visible support.
Authors: We agree that the abstract would be strengthened by a concise reference to the analysis methods supporting the central claim. In the revised version we have added a brief clause noting that the results are obtained from multi-exponential fits to correlation functions constructed with extended operators, with stability verified across fit ranges. The full technical details remain in the Methods section. revision: yes
-
Referee: Methods (correlation-function analysis): at T = 305 MeV the temporal extent Nt is necessarily small; the manuscript gives no indication that excited-state contamination or operator-dependent artifacts were controlled through smearing-parameter variation, multi-state fits, or cross-checks with point sources, which directly affects the reliability of the reported widths and size hierarchy.
Authors: We thank the referee for emphasizing this important control. The manuscript already describes the use of extended meson operators, which are chosen precisely to improve ground-state overlap and reduce excited-state contamination. Multi-state fits were performed and fit-range stability was checked; these procedures are documented in the correlation-function analysis subsection. To make the controls more explicit we have added a short paragraph summarizing the smearing-parameter variations and consistency checks with the extracted widths and size hierarchy. revision: yes
Circularity Check
Lattice QCD simulation results are direct computations with no circular derivation chain
full rationale
The paper reports outcomes from explicit lattice QCD simulations: HISQ action for sea quarks, Wilson clover for valence charm, and correlation functions computed with extended meson operators. Claims of consistency with charmonium states below threshold and temperature-dependent thermal widths are extracted from these correlators via standard analysis. No derivation reduces to self-definition, fitted inputs renamed as predictions, or load-bearing self-citations; the work is self-contained against external lattice benchmarks and does not invoke uniqueness theorems or ansatze from prior author work to force results.
Axiom & Free-Parameter Ledger
axioms (1)
- domain assumption Lattice QCD with HISQ dynamical quarks and Wilson-clover valence charm quarks yields correlation functions that can be interpreted in terms of physical charmonium states and their thermal widths.
Lean theorems connected to this paper
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
We assume that in-medium S-wave charmonia masses and widths are described by a Lorentzian form... σ_med_α(ω,T)=A_α(T)/π Γ_α(ω,T)/((ω−M_α(T))²+Γ²_α(ω,T))+z_α ω δ(ω)
-
IndisputableMonolith/Foundation/RealityFromDistinction.leanreality_from_one_distinction unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
Our lattice QCD results are consistent with the existence of all charmonium states below the open charm threshold... thermal width follows the hierarchy of charmonium sizes
What do these tags mean?
- matches
- The paper's claim is directly supported by a theorem in the formal canon.
- supports
- The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
- extends
- The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
- uses
- The paper appears to rely on the theorem as machinery.
- contradicts
- The paper's claim conflicts with a theorem or certificate in the canon.
- unclear
- Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.
Reference graph
Works this paper leans on
- [1]
-
[2]
Heavy-flavor production and medium properties in high-energy nuclear collisions - What next?
G. Aartset al., Eur. Phys. J. A53, 93 (2017), arXiv:1612.08032 [nucl-th]
work page internal anchor Pith review Pith/arXiv arXiv 2017
- [3]
-
[4]
A. Andronicet al., Eur. Phys. J. A60, 88 (2024), arXiv:2402.04366 [nucl-th]
-
[5]
Color screening in (2+1)-flavor QCD
A. Bazavov, N. Brambilla, P. Petreczky, A. Vairo, and J. H. Weber (TUMQCD), Phys. Rev. D98, 054511 (2018), arXiv:1804.10600 [hep-lat]
work page internal anchor Pith review Pith/arXiv arXiv 2018
-
[6]
Potential NRQCD: an effective theory for heavy quarkonium
N. Brambilla, A. Pineda, J. Soto, and A. Vairo, Nucl. Phys. B566, 275 (2000), arXiv:hep-ph/9907240
work page internal anchor Pith review Pith/arXiv arXiv 2000
-
[7]
Effective field theories for heavy quarkonium
N. Brambilla, A. Pineda, J. Soto, and A. Vairo, Rev. Mod. Phys.77, 1423 (2005), arXiv:hep-ph/0410047
work page internal anchor Pith review Pith/arXiv arXiv 2005
-
[8]
Real-time static potential in hot QCD
M. Laine, O. Philipsen, P. Romatschke, and M. Tassler, JHEP03, 054 (2007), arXiv:hep-ph/0611300
work page internal anchor Pith review Pith/arXiv arXiv 2007
-
[9]
Static quark-antiquark pairs at finite temperature
N. Brambilla, J. Ghiglieri, A. Vairo, and P. Petreczky, Phys. Rev. D78, 014017 (2008), arXiv:0804.0993 [hep-ph]
work page internal anchor Pith review Pith/arXiv arXiv 2008
-
[10]
M. A. Escobedo and J. Soto, Phys. Rev. A78, 032520 (2008), arXiv:0804.0691 [hep-ph]
work page internal anchor Pith review Pith/arXiv arXiv 2008
-
[11]
How to compute the thermal quarkonium spectral function from first principles?
M. Laine, Nucl. Phys. A820, 25C (2009), arXiv:0810.1112 [hep-ph]
work page internal anchor Pith review Pith/arXiv arXiv 2009
-
[12]
Charmonium at finite temperature in quenched lattice QCD
T. Umeda, K. Nomura, and H. Matsufuru, Eur. Phys. J. C39S1, 9 (2005), arXiv:hep-lat/0211003
work page internal anchor Pith review Pith/arXiv arXiv 2005
-
[13]
A Study of Charmonium Systems across the Deconfinement Transition
S. Datta, F. Karsch, P. Petreczky, and I. Wetzorke, Nucl. Phys. B Proc. Suppl.119, 487 (2003), arXiv:hep-lat/0208012
work page internal anchor Pith review Pith/arXiv arXiv 2003
-
[14]
Hadron correlators, spectral functions and thermal dilepton rates from lattice QCD
F. Karsch, S. Datta, E. Laermann, P. Petreczky, S. Stickan, and I. Wetzorke, Nucl. Phys. A 715, 701 (2003), arXiv:hep-ph/0209028 . – 35 –
work page internal anchor Pith review Pith/arXiv arXiv 2003
-
[15]
Behavior of Charmonium Systems after Deconfinement
S. Datta, F. Karsch, P. Petreczky, and I. Wetzorke, Phys. Rev. D69, 094507 (2004), arXiv:hep-lat/0312037
work page internal anchor Pith review Pith/arXiv arXiv 2004
-
[16]
$J/\psi$ and $\eta_c$ in the Deconfined Plasma from Lattice QCD
M. Asakawa and T. Hatsuda, Phys. Rev. Lett.92, 012001 (2004), arXiv:hep-lat/0308034
work page internal anchor Pith review Pith/arXiv arXiv 2004
-
[17]
Quarkonium correlators and spectral functions at zero and finite temperature
A. Jakovac, P. Petreczky, K. Petrov, and A. Velytsky, Phys. Rev. D75, 014506 (2007), arXiv:hep-lat/0611017
work page internal anchor Pith review Pith/arXiv arXiv 2007
-
[18]
H. Ohno, S. Aoki, S. Ejiri, K. Kanaya, Y. Maezawa, H. Saito, and T. Umeda (WHOT-QCD), Phys. Rev. D84, 094504 (2011), arXiv:1104.3384 [hep-lat]
work page internal anchor Pith review Pith/arXiv arXiv 2011
-
[19]
H. T. Ding, A. Francis, O. Kaczmarek, F. Karsch, H. Satz, and W. Soeldner, Phys. Rev. D 86, 014509 (2012), arXiv:1204.4945 [hep-lat]
work page internal anchor Pith review Pith/arXiv arXiv 2012
-
[20]
H.-T. Ding, O. Kaczmarek, S. Mukherjee, H. Ohno, and H. T. Shu, Phys. Rev. D97, 094503 (2018), arXiv:1712.03341 [hep-lat]
work page internal anchor Pith review Pith/arXiv arXiv 2018
- [21]
-
[22]
Bottomonium above deconfinement in lattice nonrelativistic QCD
G. Aarts, S. Kim, M. P. Lombardo, M. B. Oktay, S. M. Ryan, D. K. Sinclair, and J. I. Skullerud, Phys. Rev. Lett.106, 061602 (2011), arXiv:1010.3725 [hep-lat]
work page internal anchor Pith review Pith/arXiv arXiv 2011
-
[23]
G. Aarts, C. Allton, S. Kim, M. P. Lombardo, M. B. Oktay, S. M. Ryan, D. K. Sinclair, and J. I. Skullerud, JHEP11, 103 (2011), arXiv:1109.4496 [hep-lat]
work page internal anchor Pith review Pith/arXiv arXiv 2011
-
[24]
S wave bottomonium states moving in a quark-gluon plasma from lattice NRQCD
G. Aarts, C. Allton, S. Kim, M. P. Lombardo, M. B. Oktay, S. M. Ryan, D. K. Sinclair, and J.-I. Skullerud, JHEP03, 084 (2013), arXiv:1210.2903 [hep-lat]
work page internal anchor Pith review Pith/arXiv arXiv 2013
-
[25]
Melting of P wave bottomonium states in the quark-gluon plasma from lattice NRQCD
G. Aarts, C. Allton, S. Kim, M. P. Lombardo, S. M. Ryan, and J. I. Skullerud, JHEP12, 064 (2013), arXiv:1310.5467 [hep-lat]
work page internal anchor Pith review Pith/arXiv arXiv 2013
-
[26]
The bottomonium spectrum at finite temperature from $N_f=2+1$ lattice QCD
G. Aarts, C. Allton, T. Harris, S. Kim, M. P. Lombardo, S. M. Ryan, and J.-I. Skullerud, JHEP07, 097 (2014), arXiv:1402.6210 [hep-lat]
work page internal anchor Pith review Pith/arXiv arXiv 2014
-
[27]
S. Kim, P. Petreczky, and A. Rothkopf, Phys. Rev. D91, 054511 (2015), arXiv:1409.3630 [hep-lat]
work page internal anchor Pith review Pith/arXiv arXiv 2015
-
[28]
S. Kim, P. Petreczky, and A. Rothkopf, JHEP11, 088 (2018), arXiv:1808.08781 [hep-lat]
work page internal anchor Pith review Pith/arXiv arXiv 2018
-
[29]
Can quarkonia survive deconfinement ?
A. Mocsy and P. Petreczky, Phys. Rev. D77, 014501 (2008), arXiv:0705.2559 [hep-ph]
work page internal anchor Pith review Pith/arXiv arXiv 2008
-
[30]
Quarkonium spectral functions with complex potential
P. Petreczky, C. Miao, and A. Mocsy, Nucl. Phys. A855, 125 (2011), arXiv:1012.4433 [hep-ph]
work page internal anchor Pith review Pith/arXiv arXiv 2011
-
[31]
Quarkonium at finite temperature: Towards realistic phenomenology from first principles
Y. Burnier, O. Kaczmarek, and A. Rothkopf, JHEP12, 101 (2015), arXiv:1509.07366 [hep-ph]
work page internal anchor Pith review Pith/arXiv arXiv 2015
-
[32]
H. Ohno, T. Umeda, and K. Kanaya (WHOT-QCD), PoSLATTICE2008, 203 (2008), arXiv:0810.3066 [hep-lat]
work page internal anchor Pith review Pith/arXiv arXiv 2008
- [33]
- [34]
- [35]
-
[36]
H.-T. Ding, W.-P. Huang, R. Larsen, S. Meinel, S. Mukherjee, P. Petreczky, and Z. Tang, JHEP05, 149 (2025), arXiv:2501.11257 [hep-lat]
-
[37]
C. T. H. Davies, K. Hornbostel, A. Langnau, G. P. Lepage, A. Lidsey, J. Shigemitsu, and J. H. Sloan, Phys. Rev. D50, 6963 (1994), arXiv:hep-lat/9406017
work page internal anchor Pith review Pith/arXiv arXiv 1994
-
[38]
The bottomonium spectrum from lattice QCD with 2+1 flavors of domain wall fermions
S. Meinel, Phys. Rev. D79, 094501 (2009), arXiv:0903.3224 [hep-lat]
work page internal anchor Pith review Pith/arXiv arXiv 2009
-
[39]
S. Meinel, Phys. Rev. D82, 114502 (2010), arXiv:1007.3966 [hep-lat]
work page internal anchor Pith review Pith/arXiv arXiv 2010
-
[40]
T. C. Hammant, A. G. Hart, G. M. von Hippel, R. R. Horgan, and C. J. Monahan, Phys. Rev. Lett.107, 112002 (2011), [Erratum: Phys.Rev.Lett. 115, 039901 (2015)], arXiv:1105.5309 [hep-lat]
work page internal anchor Pith review Pith/arXiv arXiv 2011
-
[41]
R. J. Dowdallet al.(HPQCD), Phys. Rev. D85, 054509 (2012), arXiv:1110.6887 [hep-lat]
work page internal anchor Pith review Pith/arXiv arXiv 2012
-
[42]
J. O. Daldrop, C. T. H. Davies, and R. J. Dowdall (HPQCD), Phys. Rev. Lett.108, 102003 (2012), arXiv:1112.2590 [hep-lat]
work page internal anchor Pith review Pith/arXiv arXiv 2012
-
[43]
A. Bazavov, D. Hoying, R. N. Larsen, S. Mukherjee, P. Petreczky, A. Rothkopf, and J. H. Weber (HotQCD), Phys. Rev. D109, 074504 (2024), arXiv:2308.16587 [hep-lat]
-
[44]
D. Bollweg, J. L. Dasilva Gol´ an, O. Kaczmarek, R. N. Larsen, G. D. Moore, S. Mukherjee, P. Petreczky, H.-T. Shu, S. Stendebach, and J. H. Weber (HotQCD), JHEP09, 180 (2025), arXiv:2506.11958 [hep-lat]
-
[45]
Results for light pseudoscalar mesons
A. Bazavovet al.(MILC), PoSLATTICE2010, 074 (2010), arXiv:1012.0868 [hep-lat]
work page internal anchor Pith review Pith/arXiv arXiv 2010
- [46]
-
[47]
Flavor Symmetry and the Static Potential with Hypercubic Blocking
A. Hasenfratz and F. Knechtli, Phys. Rev. D64, 034504 (2001), arXiv:hep-lat/0103029
work page internal anchor Pith review Pith/arXiv arXiv 2001
-
[48]
L. Mazuret al.(HotQCD), Comput. Phys. Commun.300, 109164 (2024), arXiv:2306.01098 [hep-lat]
-
[49]
T. Izubuchi, L. Jin, C. Kallidonis, N. Karthik, S. Mukherjee, P. Petreczky, C. Shugert, and S. Syritsyn, Phys. Rev. D100, 034516 (2019), arXiv:1905.06349 [hep-lat]
- [50]
-
[51]
On the generalized eigenvalue method for energies and matrix elements in lattice field theory
B. Blossier, M. Della Morte, G. von Hippel, T. Mendes, and R. Sommer, JHEP04, 094 (2009), arXiv:0902.1265 [hep-lat]
work page internal anchor Pith review Pith/arXiv arXiv 2009
-
[52]
E. Eichten, K. Gottfried, T. Kinoshita, K. D. Lane, and T.-M. Yan, Phys. Rev. D21, 203 (1980)
work page 1980
-
[53]
Navaset al.(Particle Data Group), Phys
S. Navaset al.(Particle Data Group), Phys. Rev. D110, 030001 (2024)
work page 2024
-
[54]
Constant contribution in meson correlators at finite temperature
T. Umeda, Phys. Rev. D75, 094502 (2007), arXiv:hep-lat/0701005
work page internal anchor Pith review Pith/arXiv arXiv 2007
-
[55]
On temperature dependence of quarkonium correlators
P. Petreczky, Eur. Phys. J. C62, 85 (2009), arXiv:0810.0258 [hep-lat]
work page internal anchor Pith review Pith/arXiv arXiv 2009
-
[56]
Continuum and lattice meson spectral functions at nonzero momentum and high temperature
G. Aarts and J. M. Martinez Resco, Nucl. Phys. B726, 93 (2005), arXiv:hep-lat/0507004
work page internal anchor Pith review Pith/arXiv arXiv 2005
-
[57]
Heavy Quark Diffusion from the Lattice
P. Petreczky and D. Teaney, Phys. Rev. D73, 014508 (2006), arXiv:hep-ph/0507318
work page internal anchor Pith review Pith/arXiv arXiv 2006
- [58]
- [59]
-
[60]
Signatures of charmonium modification in spatial correlation functions
F. Karsch, E. Laermann, S. Mukherjee, and P. Petreczky, Phys. Rev. D85, 114501 (2012), arXiv:1203.3770 [hep-lat]
work page internal anchor Pith review Pith/arXiv arXiv 2012
-
[61]
In-medium modifications of open and hidden strange-charm mesons from spatial correlation functions
A. Bazavov, F. Karsch, Y. Maezawa, S. Mukherjee, and P. Petreczky, Phys. Rev. D91, 054503 (2015), arXiv:1411.3018 [hep-lat]
work page internal anchor Pith review Pith/arXiv arXiv 2015
-
[62]
P. Petreczky, S. Sharma, and J. H. Weber, Phys. Rev. D104, 054511 (2021), arXiv:2107.11368 [hep-lat]
- [63]
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.