Lattice study of correlators for quarkonium decay

Antonio Vairo; Balbeer Singh; Debasish Banerjee; Julian Mayer-Steudte; Marc Janer; Nora Brambilla; Peter Petreczky; Saumen Datta; Viljami Leino

REVIEW 3 major objections 4 minor 14 references

Lattice study of correlators for quarkonium decay

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

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

Pith's one-line read The nonperturbative color-electric correlator driving quarkonium decay is asymmetric in Euclidean time and matches NLO at 10^4 Tc.

desk verdict A preliminary lattice report whose central asymmetry claim depends on an undefined renormalization convention; the NLO match is a good check, but the scheme ambiguity undercuts the headline result. read the letter →

arxiv 2506.22594 v1 pith:SH6CIW2L submitted 2025-06-27 hep-lat hep-phnucl-th

Saumen Datta , Debasish Banerjee , Nora Brambilla , Marc Janer , Viljami Leino , Julian Mayer-Steudte , Peter Petreczky , Balbeer Singh

show 1 more author

Antonio Vairo

This is my paper

classification hep-lathep-phnucl-th PACS 12.38.Gc12.38.Mh

keywords latticeQCDquarkoniumcolorelectriccorrelatoradjointWilsonlinequark-gluonplasmapotentialnonrelativisticheavyquarkdiffusionthermalfieldtheory

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 establishes, from lattice data in a gluonic plasma, that the color-electric correlator $G_E(\tau)$ governing singlet-to-octet decay of bottomonia is not symmetric about $\tau T=1/2$. That asymmetry separates it from the heavy-quark-diffusion correlator $G_F^E(\tau)$, whose fundamental Wilson-line orientation makes it symmetric. The antisymmetric part is an $O(g^4)$ effect in perturbation theory, so it should shrink as the temperature rises; the lattice results show this trend, and at $T=10^4 T_c$ the continuum-extrapolated correlator agrees with the NLO calculation. The result matters because $G_E(\tau)$ encodes the medium effect on quarkonium decay, so its true shape is necessary input for describing quarkonium suppression in heavy-ion collisions.

What carries the argument

The object carrying the argument is the adjoint color-electric correlator $G_E(\tau)$, together with its renormalized version $G_r^E(\tau)=(L_a^r(T)/L_a(a;T))^{\tau T} G_E(\tau;a)$, where $L_a$ is the adjoint Polyakov loop. This renormalization removes the exponential mass divergence of the adjoint Wilson line and leaves only mild cutoff effects between $\tau T=0.3$ and $0.8$. The comparison object is the heavy-quark-diffusion correlator $G_F^E(\tau)$ in Eq. (3), whose fundamental Wilson line and Polyakov-loop division make it symmetric; the difference in Wilson-line orientation is what produces the asymmetry. The leading-order lattice normalization $f(\tau;a)$ in Eq. (5) is used to expose the shape of the correlator and the size of cutoff effects.

What would settle it

Recompute $G_E(\tau)$ at $1.5 T_c$ using a renormalization scheme that does not rely on the adjoint Polyakov-loop ratio, for example a nonperturbative subtraction of the Wilson-line self-energy, and check whether the continuum extrapolation remains asymmetric about $\tau T=1/2$; if the asymmetry disappears under a different scheme, the claim of an intrinsic nonperturbative asymmetry fails.

Watch

Extended reading notes

Core claim

The central discovery is that the nonperturbative $G_E(\tau)$ in a gluonic plasma is asymmetric about $\tau T=1/2$. The correlator is defined as $G_E(\tau)=-\frac{1}{3}\sum_i \langle E_i^a(\tau) W^{ab}(\tau,0) E_i^b(0)\rangle_T$, where $W^{ab}$ is an adjoint Wilson line; the asymmetry originates from the orientation of this Wilson line, in contrast to the fundamental Wilson line used in $G_F^E(\tau)$. The paper shows that the asymmetry is mild at $1.5 T_c$, decreases with increasing temperature, and at $T=10^4 T_c$ the continuum limit of $G_E(\tau)$ agrees well with the NLO perturbative result. It also reports that within errors the octet-octet correlator $G_{\rm oct}(\tau)$ obeys the NLO scaling relative to $G_F^E(\tau)$.

Load-bearing premise

The load-bearing premise is that the entire mass divergence of the adjoint Wilson line correlator is exponential and is exactly removed by the ratio of renormalized to bare adjoint Polyakov loops; if any cutoff-dependent piece survives that subtraction, the continuum extrapolation and the reported asymmetry could be artifacts.

Editorial extensions

If this is right

A non-symmetric $G_E(\tau)$ means the real-time correlator $G_{EE}^>(\omega)$ has an antisymmetric component, so quarkonium decay widths and the transport coefficient $\kappa$ cannot be extracted from fits that assume $G_E(\tau)$ is symmetric.
Agreement with NLO at $T=10^4 T_c$ validates the dipole effective-field-theory description of quarkonium-medium interaction at asymptotically high temperatures.
Because the adjoint and fundamental Wilson-line correlators have different symmetry properties, heavy-quark-diffusion results cannot be directly reused for quarkonium decay; each process needs its own nonperturbative correlator.
The small cutoff effects after the Eq. (4) renormalization at $1.5 T_c$ make continuum extrapolation feasible at phenomenologically relevant temperatures, allowing direct input for Lindblad and Boltzmann descriptions of $\Upsilon$ suppression.

Reading between the lines

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

If the antisymmetric part grows as the coupling increases, the $\omega\to 0$ limit of $G_{EE}^>$ may differ substantially from the NLO estimate near $T_c$, which would sharpen the predicted temperature dependence of quarkonium dissociation; this is a testable extension of the present data.
The same adjoint-Polyakov-loop renormalization applied to $G_{\rm oct}(\tau)$, which the paper finds symmetric, provides a control: a scheme-induced asymmetry would likely affect both correlators, so comparing the two could separate physical from cutoff effects.
For spectral-function reconstruction, the observed asymmetry implies standard methods that assume $G(-i\omega)$ symmetry would need to be generalized to include an odd part, or they would mis-model the decay width.

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.

Desk Editor's Note

read the letter

Quick read: this is a conference-style report from a strong collaboration, presenting lattice results for the color-electric correlator relevant to quarkonium decay. The genuinely new pieces relative to the companion paper [11] are the temperature scan (1.4 to 10^4 Tc) and the NLO comparison at 10^4 Tc. The data handling looks careful: multiple lattice spacings at each temperature, a continuum extrapolation at 1.5 Tc, and a clean NLO match at 10^4 Tc. The observation that G_E is not symmetric about τT=1/2, unlike the heavy-quark diffusion correlator, is interesting and visually striking.

The soft spot is the renormalization. Equation (4) defines G_E^r via the ratio (L_r_a/L_a)^{τT}, but L_r_a(T) is never defined. The stress-test note is right: any finite renormalization choice, say a temperature-independent Z_L multiplying the loop, changes G_E^r by e^{ln Z_L · τT}. That is not an overall constant; it is an exponential tilt that manufactures an antisymmetric component. Unless the paper fixes L_r_a(T) by an explicit condition—for example, matching NLO at 10^4 Tc and then using the same convention at all T—the plotted asymmetry is not a scheme-independent observable. The paper simply asserts that "the major part of the cutoff effect" is removed, which is not enough.

Also, the continuum extrapolated curve in Fig. 2 has no error band, a minor issue for a proceedings but worth noting. The overlap with [11] is substantial; the new content is the temperature dependence and the NLO comparison rather than the correlator itself.

Is this worth engaging? Yes, as a proceedings contribution it is a useful status report, and the NLO match at high T is a good sanity check. But if the asymmetry claim is meant to carry phenomenological weight, the renormalization scheme has to be nailed down. I would send it to peer review if the venue wants to enforce reproducibility, but the referee should ask for a definition of L_r_a(T) and an error band on the continuum result. If this is only a proceedings note, the pressure is lower, but the claim as stated is premature.

Bottom line: the paper is honest and technically careful in its own terms, but the central observation is not yet fully defined. A serious referee should engage with it, mainly to get the renormalization scheme on the table.

Referee Report

3 major / 4 minor

Summary. This paper reports a lattice study of the Euclidean color-electric correlator GE(τ) with an adjoint Wilson line, which enters the pNRQCD description of quarkonium decay in a quark-gluon plasma. The authors propose a renormalization of the bare lattice correlator via the adjoint Polyakov loop in Eq. (4), present results for a gluonic plasma at several temperatures and lattice spacings, and argue that, unlike the heavy-quark-diffusion correlator GEF(τ), the nonperturbative GE(τ) is not symmetric about τT=1/2. They also report agreement with the NLO perturbative result at T=10^4 Tc and discuss the octet-octet correlator GEoct(τ). The manuscript is explicitly preliminary, but the central asymmetry claim is stated without qualification.

Significance. If the asymmetry claim survives scrutiny, it is a genuinely useful nonperturbative input for quarkonium-in-QGP phenomenology: it distinguishes the decay correlator from the diffusion correlator and provides a target for perturbative calculations. The strength of the paper is that the lattice calculation is an independent measurement with no fitted parameters, uses multiple lattice spacings, includes a continuum extrapolation at 1.5 Tc, and benchmarks against an external NLO result at very high temperature. The main weakness is that the renormalization condition underlying Eq. (4) is not defined, so the central observable is not yet scheme-independent.

major comments (3)

[Eq. (4)] The renormalization prescription in Eq. (4) is not complete. The factor (Lr_a(T)/La(a;T))^{τT} removes the power divergence of the adjoint Wilson line, but the finite normalization of Lr_a(T) is never specified. Any multiplicative renormalization of Lr_a(T), even a temperature-independent constant Z, changes GEr(τ) by Z^{τT}, which is a τ-dependent reweighting rather than an overall constant. Such a factor itself generates an antisymmetric component around τT=1/2, so the plotted asymmetry is not yet a scheme-independent observable. The text must state the condition that fixes Lr_a(T) and should show that the asymmetry persists for a range of admissible conventions, or the central claim is a subtraction artifact.
[Fig. 2 and surrounding text] The continuum extrapolated curve in the left panel of Fig. 2 is shown without an error band, uncertainties, or the extrapolation ansatz. Since the asymmetry claim is based on this curve, the reader cannot assess whether the deviation from symmetry is statistically significant. Please provide the continuum values with errors, the fit form (e.g., linear in a^2), the fit range, and the chi^2/dof. In addition, quantify the asymmetry directly, for example by plotting or tabulating GEr(τ) - GEr(1/T-τ), rather than relying on visual inspection.
[Fig. 3 and high-temperature comparison] The agreement with NLO at T=10^4 Tc is presented as a validation, but the continuum result in Fig. 3 has no error band and no goodness-of-fit is reported, so the strength of the agreement is unclear. Moreover, because the NLO calculation in Ref. [10] shares authors with the present work, the text should state explicitly that no parameters from that calculation are adjusted in the comparison and that the same renormalization convention for Lr_a(T) is used at T=10^4 Tc and at lower temperatures. Otherwise the high-temperature match cannot serve to fix the scheme for the asymmetry claim.

minor comments (4)

[Eq. (5)] The function f(τ;a) is called both the normalization and the leading-order lattice correlator; please clarify that the plotted quantity is GEr(τ)/f(τ;a) and that f is evaluated with the same lattice discretization.
[Fig. 1 caption] The caption contains a typo ('pN RQCDdue') and does not clearly explain how the three diagrams correspond to Eqs. (2), (3), and (6); a few words connecting each diagram to its equation would help.
[Text near Eq. (4)] The exponential factor e^{δm(a)τ} is introduced but δm(a) is not defined, and the equality with the Polyakov-loop ratio is stated without derivation or reference; please define δm(a) and give the argument for the factorization.
[General presentation] The manuscript is explicitly preliminary; for a journal publication, the lattice setup should be documented in more detail, including the gauge action, the number of configurations, the physical scales, and the determination of the lattice spacing.

Circularity Check

0 steps flagged · score 1.0 of 10

No significant circularity: the lattice correlator is an independent measurement, and the NLO comparison is an external benchmark; the only caveats are routine self-citations and a scheme-dependence in Eq. (4) that is a correctness issue, not a circular derivation.

full rationale

The paper's central result is a direct lattice measurement of the color-electric correlator G_E(τ) after multiplicative renormalization. The renormalization in Eq. (4) uses the adjoint Polyakov loop ratio as a separately measured lattice quantity; it is not fitted to the target correlator, and the continuum extrapolation is shown for several lattice spacings. The observed lack of symmetry about τT=1/2 is a property of the renormalized lattice data, not an output of a formula that explicitly contains the asymmetry as an input. The high-temperature comparison with the NLO calculation of Ref. [10] is a benchmark check, not an input: the NLO result is an independent perturbative calculation and is not used to fix parameters in the lattice measurement. Self-citations to Refs. [10] and [11] are from the same collaboration and provide context, the renormalization convention, and earlier results, but they do not force the present conclusion; the central lattice calculation is self-contained. A legitimate caveat, noted in the text only implicitly, is that the finite normalization of the renormalized adjoint Polyakov loop L_r^a(T) is not specified, and a different finite renormalization would multiply G_E^r(τ) by a τ-dependent factor, altering the asymmetry. This is a scheme-dependence/correctness concern rather than a circularity: it does not make the lattice derivation equivalent to its inputs by construction.

Assumptions & free parameters 0 free parameters · 4 assumptions · 0 invented entities

The central lattice result relies on the pNRQCD dipole approximation, the exponential mass-divergence renormalization, and the analytic continuation assumption; no free parameters are fitted. The only input quantities are the temperature in units of Tc and the lattice spacing, which are set by the simulation.

assumptions (4)

domain assumption The dipole approximation r << 1/T, where r is quarkonium size and T the temperature, so the medium interaction is a color dipole g r.E.
Invoked in the abstract and Section 1 to justify the pNRQCD treatment and the correlator definition; valid for ground state bottomonia at LHC energies but restricts the regime of applicability.
domain assumption The mass divergence of the adjoint Wilson line correlator is purely exponential and is removed by the adjoint Polyakov loop ratio in Eq. (4).
Eq. (4) defines GEr = (Lr_a/La)^{τT} GE; the paper states this removes the major cutoff effect but does not prove it removes all divergences.
domain assumption The Euclidean correlator GE(τ) can be analytically continued to real time and its spectral function yields the decay width via Eq. (1).
Standard analytic continuation assumption in thermal field theory, the basis for connecting the lattice observable to quarkonium decay.
domain assumption The NLO perturbative result of Ref. [10] is a reliable benchmark at T=10^4 Tc.
Used for the comparison in Fig. 3; assumes that at this high temperature perturbation theory is valid and the lattice discretization effects are under control.

how reviews work

0 comments

Cite this review

Pith. "Pith review of Lattice study of correlators for quarkonium decay." pith.science (2026). https://pith.science/paper/SH6CIW2L

@misc{pith2026250622594,
  author       = {Pith},
  title        = {Pith review of: Lattice study of correlators for quarkonium decay},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/SH6CIW2L}},
  note         = {Machine review of arXiv:2506.22594}
}

read the original abstract

While there has been a lot of progress in developing a formalism for the study of quarkonia in QGP, a nonperturbative study is still difficult. For bottomonia, where the system size is much less than the inverse temperature, the interaction of the system with the medium can be approximated by a dipole interaction with the color electric field. The decay of the quarkonia can be connected to a correlation function of the color electric field. We present preliminary results from a lattice study of the relevant color electric field correlator. The structure of the correlator, and its difference from the corresponding correlator studied for heavy quark diffusion, is discussed.

Figures

Figures reproduced from arXiv: 2506.22594 by the authors.

**Figure 2.** (left) GE r (τ ; a) at different lattice spacings and the continuum extrapolation, at 1.5 Tc. (Right) GE r (τ ; a)/f(τ ) at different temperatures. At each temperature results for two different lattice spacings are shown. been calculated nonperturbatively using lattice QCD [11]. In this report we present lattice results for GE (τ ) at various temperatures in a gluonic plasma. GE (τ ; a) has a mass divergence coming … view at source ↗

**Figure 3.** GE r (τ ) at 104Tc. (Left) Results at different lattice spacings, and the continuum limit. (Right) Comparison of the continuum GE r (τ ; T = 104Tc) with the NLO result of Ref. [10]. where the leading order calculation is done for the lattice discretized correlator. In the range τT ∈ [0.3, 0.8] we find mild cutoff effect in the renormalized correlator. The continuum extrapolated correlator is also shown in [PITH_FUL… view at source ↗

Discussion (0). Sign in to comment.

Reference graph

Works this paper leans on

14 extracted references · 12 canonical work pages

[11]

Lattice study of correlators of chromoelectric fields for heavy quarkonium dynamics in the quark-gluon plasma

N. Brambilla, S. Datta, M. Janer, V. Leino, J. Mayer-Steudte, P. Petreczky and A. Vairo, arXiv:2505.16603

work page Pith review arXiv
[10]

The chromoelectric adjoint correlators in Euclidean space at next-to-leading order

N. Brambilla, P. Panayiotou, S. S¨ appi and A. Vairo, arXiv:2505.16604

work page Pith review arXiv
[1]

Matsui and H

T. Matsui and H. Satz, Phys. Lett. B 178 (1986) 416

1986
[2]

Brambilla, J

N. Brambilla, J. Ghiglieri, P. Petreczky and A. Vairo, Phys. Rev. D78 (2008) 014017

work page 2008
[3]

Akamatsu, Phys

Y. Akamatsu, Phys. Rev. D 87 (2013) 045016; Phys. Rev. D 91 (2015) 056002

work page 2013
[4]

Pineda and J

A. Pineda and J. Soto, Nucl. Phys. B Proc. Suppl. 64 (1998) 428. N. Brambilla, A. Pineda, J. Soto and A. Vairo, Nucl. Phys. B 566 (2000) 275

work page 1998
[5]

Brambilla, M

N. Brambilla, M. Escobedo, J. Soto and A. Vairo, Phys. Rev. D 96 (2017) 034201; Phys. Rev. D97 (2018) 074009

work page 2017
[6]

Brambilla, T

N. Brambilla, T. Magorsch, M. Strickland, A. Vairo and P.V. Griend, Phys. Rev. D 109 (2024) 114016

work page 2024

Show all 14 references

[7]

and Rishi Sharma, arXiv:2504.19348

Vyshakh B.R. and Rishi Sharma, arXiv:2504.19348
[8]

Yao and T

X. Yao and T. Mehen, JHEP 02 (2021) 062

2021
[9]

Scheihing-Hitschfeld and X

B. Scheihing-Hitschfeld and X. Yao, Phys. Rev. D 108 (2023) 054024

2023
[12]

Caron-Huot, M

S. Caron-Huot, M. Laine and G. Moore, JHEP 04 (2009) 053

2009
[13]

Banerjee, S

D. Banerjee, S. Datta, R. Gavai and P. Majumdar, Phys. Rev. D 85 (2012) 014510; Nucl. Phys. A 1038 (2023) 122721. L. Altenkort, A.M. Eller, O. Kaczmarek, L. Mazur, G.D. Moore and H-T.Shu, Phys. Rev. D 103 (2021) 014511. N. Brambilla, V. Leino, J. Mayer-Steudte and P. Petreczky...

2012
[14]

Altenkort, O

L. Altenkort, O. Kaczmarek, R. Larsen, S. Mukherjee, P. Petreczky, H-T. Shu and S. Stendebach, Phys. Rev. Lett. 130 (2023) 231902. 5

2023

Pith tools

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

[1] [11]

Lattice study of correlators of chromoelectric fields for heavy quarkonium dynamics in the quark-gluon plasma

N. Brambilla, S. Datta, M. Janer, V. Leino, J. Mayer-Steudte, P. Petreczky and A. Vairo, arXiv:2505.16603

work page Pith review arXiv

[2] [10]

The chromoelectric adjoint correlators in Euclidean space at next-to-leading order

N. Brambilla, P. Panayiotou, S. S¨ appi and A. Vairo, arXiv:2505.16604

work page Pith review arXiv

[3] [1]

Matsui and H

T. Matsui and H. Satz, Phys. Lett. B 178 (1986) 416

1986

[4] [2]

Brambilla, J

N. Brambilla, J. Ghiglieri, P. Petreczky and A. Vairo, Phys. Rev. D78 (2008) 014017

work page 2008

[5] [3]

Akamatsu, Phys

Y. Akamatsu, Phys. Rev. D 87 (2013) 045016; Phys. Rev. D 91 (2015) 056002

work page 2013

[6] [4]

Pineda and J

A. Pineda and J. Soto, Nucl. Phys. B Proc. Suppl. 64 (1998) 428. N. Brambilla, A. Pineda, J. Soto and A. Vairo, Nucl. Phys. B 566 (2000) 275

work page 1998

[7] [5]

Brambilla, M

N. Brambilla, M. Escobedo, J. Soto and A. Vairo, Phys. Rev. D 96 (2017) 034201; Phys. Rev. D97 (2018) 074009

work page 2017

[8] [6]

Brambilla, T

N. Brambilla, T. Magorsch, M. Strickland, A. Vairo and P.V. Griend, Phys. Rev. D 109 (2024) 114016

work page 2024