pith. sign in

arxiv: 2606.10049 · v1 · pith:VPTPRZHRnew · submitted 2026-06-08 · ✦ hep-lat · hep-ph· nucl-ex· nucl-th

Momentum Dependence of Heavy Quark Diffusion in a Thermal Gluonic Plasma on the Lattice

Pith reviewed 2026-06-27 13:55 UTC · model grok-4.3

classification ✦ hep-lat hep-phnucl-exnucl-th
keywords heavy quarksdrag coefficientdiffusion coefficientlattice QCDthermal plasmagluonsmomentum dependence
0
0 comments X

The pith

A lattice method extracts momentum dependence of heavy quark drag and diffusion in a gluonic plasma

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

This paper proposes a numerical strategy to simulate the real-time dynamics of a heavy quark with different initial momenta in a three-dimensional lattice discretization of an effective theory for high-temperature QCD. The plasma consists of non-perturbatively interacting soft gluons. The approach makes it possible to determine how the drag and diffusion coefficients vary with the heavy quark's momentum. Previously such quantities were only available at zero momentum or in perturbative approximations. Knowing the momentum dependence matters for understanding how heavy quarks lose energy while traversing the quark-gluon plasma created in heavy-ion collisions.

Core claim

The authors demonstrate a numerical strategy that simulates the real-time dynamics of a heavy quark with varying initial momenta in an effective three-dimensional lattice theory of thermal gluons, thereby determining the momentum dependence of the associated drag and diffusion coefficients for the first time in a non-perturbative setting.

What carries the argument

The numerical strategy for evolving the heavy quark in the lattice-discretized effective theory of thermal gluons while measuring its momentum loss.

If this is right

  • Drag and diffusion coefficients are obtained as functions of momentum.
  • The extraction is performed in a fully non-perturbative interacting plasma.
  • The method works in a non-Abelian gauge theory without Abelian approximations.
  • It provides input for transport models of heavy quarks in hot QCD matter.

Where Pith is reading between the lines

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

  • If the coefficients decrease at higher momenta, it would imply reduced energy loss for fast heavy quarks.
  • The technique might be adapted to study light quark or gluon transport on similar lattices.
  • Finite volume effects could be quantified by changing the lattice size in future runs.

Load-bearing premise

The lattice simulation accurately captures the heavy quark's momentum evolution without introducing large uncontrolled errors from discretization or the effective theory truncation.

What would settle it

If the computed drag coefficient at large momentum deviates significantly from the known perturbative result without matching the expected approach to that limit, the numerical extraction would be called into question.

Figures

Figures reproduced from arXiv: 2606.10049 by Harshit Pandey, Sayantan Sharma.

Figure 1
Figure 1. Figure 1: Variance of the momentum distribution (⟨p 2 x ⟩ − ⟨px⟩ 2 )/T 2 for a heavy quark with a zero and finite initial momenta |p|/T = 0.54 in a typical thermal configuration consisting of non-perturbatively interacting soft SU(3) gluons at T ∼ 480 MeV. We calculate the variance of the momentum distribution along the x-direction for a heavy quark with a zero initial mo￾mentum by solving the Dirac equation, but by… view at source ↗
Figure 2
Figure 2. Figure 2: The κ/ϵ3/4 for zero initial momentum of a heavy quark, obtained in an effective theory of soft (green points) SU(2) and (orange points) SU(3) gluons, shown as a function of temperature. For a comparison, we show the values of the same ratio in QCD without dynamical quarks as a function of temperature, where the parametric dependence of κ was taken from Ref. [18] and ϵ is simply the Stefan-Boltzmann value o… view at source ↗
Figure 3
Figure 3. Figure 3: The ratio of drag and diffusion coefficients ηD/T κL/T 3 × 2  M T  as a function of initial fermion momenta |p|/T, for a high temperature effective theory of soft (green) SU(2) and (orange) SU(3) gluons at T = 1.6 Tc. At T = 1.6 Tc, we calculate the MMOD for heavy quarks that have an initial momentum, both along and transverse to their direction. This allows us to extract the drag and diffusion coefficie… view at source ↗
Figure 5
Figure 5. Figure 5: The κT (p)/ϵ3/4 of a heavy quark as a function of its initial momenta |p|/T for different temperatures in a (top panel) SU(2) and (bottom panel) SU(3) plasma consisting of soft gluons. 5. Implications of our findings & Outlook In this work, we have developed a framework based on first￾principles lattice gauge theory to calculate the momentum de￾pendence of the motion of a heavy quark inside a non-Abelian p… view at source ↗
read the original abstract

We study the dynamics of a heavy quark in a thermal plasma consisting of non-perturbatively interacting soft momentum gluons at high temperatures, described in terms of an effective theory of QCD discretized on a three-dimensional lattice. We propose a numerical strategy that allows us to simulate the dynamics of a heavy quark for different values of initial momenta in this thermalized plasma. This allows, for the first time, to extract the momentum dependence of the heavy quark drag and diffusion coefficients in a non-perturbatively interacting thermal, non-Abelian plasma.

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

1 major / 0 minor

Summary. The manuscript proposes a numerical strategy to simulate the dynamics of a heavy quark with varying initial momenta in a thermal plasma of non-perturbatively interacting soft gluons, described via an effective theory of QCD discretized on a three-dimensional lattice. This is claimed to enable, for the first time, extraction of the momentum dependence of the heavy quark drag and diffusion coefficients in a non-perturbative, non-Abelian thermal plasma.

Significance. If the proposed strategy is validated without uncontrolled systematics, the work would provide the first non-perturbative lattice results on momentum-dependent heavy-quark transport coefficients in a thermal gluonic plasma. This is relevant for heavy-flavor phenomenology in the quark-gluon plasma, where momentum dependence enters energy-loss and flow calculations. The lattice discretization of the effective theory is a methodological strength for accessing non-perturbative effects beyond weak-coupling expansions.

major comments (1)
  1. [Numerical strategy description] The manuscript presents the numerical strategy but does not report explicit validation against known limits (e.g., perturbative drag coefficient or free-theory diffusion) or quantitative error estimation for the extracted coefficients. This is load-bearing for the central claim that the approach enables reliable extraction of momentum dependence.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for their positive evaluation of the potential impact of our work and for the constructive comment. We respond to the major comment below.

read point-by-point responses
  1. Referee: [Numerical strategy description] The manuscript presents the numerical strategy but does not report explicit validation against known limits (e.g., perturbative drag coefficient or free-theory diffusion) or quantitative error estimation for the extracted coefficients. This is load-bearing for the central claim that the approach enables reliable extraction of momentum dependence.

    Authors: We agree that explicit validation against known limits and quantitative error estimation are essential to substantiate the reliability of the extracted coefficients. The current manuscript is a proposal of the numerical strategy and focuses on its description, without performing the full set of simulations or reporting numerical extractions. In a revised version we will add a dedicated section outlining validation procedures (including comparisons to the perturbative drag coefficient at high momentum and to free-theory diffusion) together with a discussion of how systematic uncertainties can be quantified. This directly addresses the load-bearing aspect of the central claim. revision: yes

Circularity Check

0 steps flagged

No significant circularity detected

full rationale

The paper proposes a numerical strategy for simulating heavy quark dynamics on a 3D lattice effective theory to extract momentum-dependent drag and diffusion coefficients for the first time. No load-bearing derivation steps, self-definitions, fitted inputs renamed as predictions, or self-citation chains are identifiable in the abstract or described approach. The central claim is a methodological proposal that remains independent of its own outputs and does not reduce to any input by construction.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

Review based on abstract; full details on parameters and assumptions not accessible.

axioms (1)
  • domain assumption The effective three-dimensional theory captures the relevant dynamics of heavy quarks in high-temperature QCD.
    Extracted from abstract description of the setup.

pith-pipeline@v0.9.1-grok · 5616 in / 1017 out tokens · 21578 ms · 2026-06-27T13:55:20.904612+00:00 · methodology

discussion (0)

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

Reference graph

Works this paper leans on

48 extracted references · 23 linked inside Pith

  1. [1]

    Andronicet al., Eur

    A. Andronicet al., Eur. Phys. J. C76, 107 (2016), arXiv:1506.03981 [nucl-ex]

  2. [2]

    Acharyaet al.(ALICE), Phys

    S. Acharyaet al.(ALICE), Phys. Lett. B813, 136054 (2021), arXiv:2005.11131 [nucl-ex]

  3. [3]

    N. J. Abdulameeret al.(PHENIX), Phys. Rev. C112, 014904 (2025), arXiv:2409.12756 [nucl-ex]

  4. [4]

    Acharyaet al.(ALICE), Eur

    S. Acharyaet al.(ALICE), Eur. Phys. J. C83, 1123 (2023), arXiv:2307.14084 [nucl-ex]

  5. [5]

    Svetitsky, Phys

    B. Svetitsky, Phys. Rev. D37, 2484 (1988)

  6. [6]

    Casalderrey-Solana and D

    J. Casalderrey-Solana and D. Teaney, Phys. Rev. D74, 085012 (2006), arXiv:hep-ph/0605199

  7. [7]

    Beraudoet al., Nucl

    A. Beraudoet al., Nucl. Phys. A979, 21 (2018), arXiv:1803.03824 [nucl-th]

  8. [8]

    M. He, H. van Hees, and R. Rapp, Prog. Part. Nucl. Phys. 130, 104020 (2023), arXiv:2204.09299 [hep-ph]

  9. [9]

    Caron-Huot and G

    S. Caron-Huot and G. D. Moore, JHEP02, 081, arXiv:0801.2173 [hep-ph]

  10. [10]

    Braaten and M

    E. Braaten and M. H. Thoma, Phys. Rev. D44, 1298 (1991)

  11. [11]

    Braaten and M

    E. Braaten and M. H. Thoma, Phys. Rev. D44, R2625 (1991). 6

  12. [12]

    Altenkort, D

    L. Altenkort, D. de la Cruz, O. Kaczmarek, R. Larsen, G. D. Moore, S. Mukherjee, P. Petreczky, H.-T. Shu, and S. Stendebach (HotQCD), Phys. Rev. Lett.132, 051902 (2024), arXiv:2311.01525 [hep-lat]

  13. [13]

    Laine, G

    M. Laine, G. D. Moore, O. Philipsen, and M. Tassler, JHEP05, 014, arXiv:0902.2856 [hep-ph]

  14. [14]

    Banerjee, S

    D. Banerjee, S. Datta, R. Gavai, and P. Majumdar, Phys. Rev. D85, 014510 (2012), arXiv:1109.5738 [hep-lat]

  15. [15]

    Francis, O

    A. Francis, O. Kaczmarek, M. Laine, and J. Langelage, PoSLATTICE2011, 202 (2011), arXiv:1109.3941 [hep- lat]

  16. [16]

    Francis, O

    A. Francis, O. Kaczmarek, M. Laine, T. Neuhaus, and H. Ohno, Phys. Rev. D92, 116003 (2015), arXiv:1508.04543 [hep-lat]

  17. [17]

    Altenkort, A

    L. Altenkort, A. M. Eller, O. Kaczmarek, L. Mazur, G. D. Moore, and H.-T. Shu, Phys. Rev. D103, 014511 (2021), arXiv:2009.13553 [hep-lat]

  18. [18]

    Brambilla, V

    N. Brambilla, V . Leino, P. Petreczky, and A. Vairo, Phys. Rev. D102, 074503 (2020), arXiv:2007.10078 [hep-lat]

  19. [19]

    Banerjee, S

    D. Banerjee, S. Datta, R. V . Gavai, and P. Majumdar, Nucl. Phys. A1038, 122721 (2023), arXiv:2206.15471 [hep-ph]

  20. [20]

    Brambilla, V

    N. Brambilla, V . Leino, J. Mayer-Steudte, and P. Pe- treczky (TUMQCD), Phys. Rev. D107, 054508 (2023), arXiv:2206.02861 [hep-lat]

  21. [21]

    Boguslavski, A

    K. Boguslavski, A. Kurkela, T. Lappi, and J. Peuron, JHEP09, 077, arXiv:2005.02418 [hep-ph]

  22. [22]

    Pandey, S

    H. Pandey, S. Schlichting, and S. Sharma, Phys. Rev. Lett. 132, 222301 (2024), arXiv:2312.12280 [hep-lat]

  23. [23]

    Altenkort, O

    L. Altenkort, O. Kaczmarek, R. Larsen, S. Mukherjee, P. Petreczky, H.-T. Shu, and S. Stendebach (HotQCD), Phys. Rev. Lett.130, 231902 (2023), arXiv:2302.08501 [hep-lat]

  24. [24]

    Bouttefeux and M

    A. Bouttefeux and M. Laine, JHEP12, 150, arXiv:2010.07316 [hep-ph]

  25. [25]

    Altenkort, D

    L. Altenkort, D. de la Cruz, O. Kaczmarek, G. D. Moore, and H.-T. Shu, Phys. Rev. D109, 114505 (2024), arXiv:2402.09337 [hep-lat]

  26. [26]

    de la Cruz, A

    D. de la Cruz, A. M. Eller, and G. D. Moore, Phys. Rev. D110, 094057 (2024), arXiv:2410.01578 [hep-ph]

  27. [27]

    Banerjee, S

    D. Banerjee, S. Datta, and M. Laine, JHEP08, 128, arXiv:2204.14075 [hep-lat]

  28. [28]

    Altenkort, A

    L. Altenkort, A. M. Eller, O. Kaczmarek, L. Mazur, G. D. Moore, and H.-T. Shu, PoSLATTICE2021, 367 (2022), arXiv:2111.12462 [hep-lat]

  29. [29]

    Bollweg, J

    D. Bollweg, J. L. Dasilva Golán, O. Kaczmarek, R. N. Larsen, G. D. Moore, S. Mukherjee, P. Petreczky, H.-T. Shu, S. Stendebach, and J. H. Weber (HotQCD), JHEP 09, 180, arXiv:2506.11958 [hep-lat]

  30. [30]

    G. D. Moore and D. Teaney, Phys. Rev. C71, 064904 (2005), arXiv:hep-ph/0412346

  31. [31]

    Bodeker, Phys

    D. Bodeker, Phys. Lett. B426, 351 (1998), arXiv:hep- ph/9801430

  32. [32]

    van Hees and R

    H. van Hees and R. Rapp, Phys. Rev. C71, 034907 (2005), arXiv:nucl-th/0412015

  33. [33]

    M. G. Mustafa, Phys. Rev. C72, 014905 (2005), arXiv:hep-ph/0412402

  34. [34]

    Golam Mustafa, D

    M. Golam Mustafa, D. Pal, and D. Kumar Srivastava, Phys. Rev. C57, 889 (1998), [Erratum: Phys.Rev.C 57, 3499–3499 (1998)], arXiv:nucl-th/9706001

  35. [35]

    Engels, F

    J. Engels, F. Karsch, and K. Redlich, Nucl. Phys. B435, 295 (1995), arXiv:hep-lat/9408009

  36. [36]

    Francis, O

    A. Francis, O. Kaczmarek, M. Laine, T. Neuhaus, and H. Ohno, Phys. Rev. D91, 096002 (2015), arXiv:1503.05652 [hep-lat]

  37. [37]

    P. B. Arnold and L. G. Yaffe, Phys. Rev. D62, 125014 (2000), arXiv:hep-ph/9912306

  38. [38]

    S. Guin, H. Pandey, and S. Sharma, Phys. Lett. B866, 139490 (2025), arXiv:2501.04397 [hep-lat]

  39. [39]

    Caron-Huot and G

    S. Caron-Huot and G. D. Moore, Phys. Rev. Lett.100, 052301 (2008), arXiv:0708.4232 [hep-ph]

  40. [40]

    S. S. Gubser, Nucl. Phys. B790, 175 (2008), arXiv:hep- th/0612143

  41. [41]

    Rajagopal, B

    K. Rajagopal, B. Scheihing-Hitschfeld, and U. A. Wiedemann, Phys. Rev. Lett.135, 242301 (2025), arXiv:2504.21139 [hep-ph]

  42. [42]

    Rajagopal, B

    K. Rajagopal, B. Scheihing-Hitschfeld, and U. A. Wiede- mann, (2026), arXiv:2606.02693 [hep-ph]

  43. [43]

    S. K. Das, F. Scardina, S. Plumari, and V . Greco, Phys. Rev. C90, 044901 (2014), arXiv:1312.6857 [nucl-th]

  44. [44]

    Caoet al., Phys

    S. Caoet al., Phys. Rev. C99, 054907 (2019), arXiv:1809.07894 [nucl-th]

  45. [45]

    Y . Xu, J. E. Bernhard, S. A. Bass, M. Nahrgang, and S. Cao, Phys. Rev. C97, 014907 (2018), arXiv:1710.00807 [nucl-th]

  46. [46]

    W. Ke, Y . Xu, and S. A. Bass, Phys. Rev. C98, 064901 (2018), arXiv:1806.08848 [nucl-th]

  47. [47]

    M. He, R. J. Fries, and R. Rapp, Phys. Rev. C86, 014903 (2012), arXiv:1106.6006 [nucl-th]

  48. [48]

    S. K. Das, F. Scardina, S. Plumari, and V . Greco, Phys. Lett. B747, 260 (2015), arXiv:1502.03757 [nucl-th] . 7