pith. sign in

arxiv: 2502.12333 · v5 · submitted 2025-02-17 · ✦ hep-ph

Impact of momentum-dependent drag coefficient on energy loss of charm and bottom quarks in QGP

Pith reviewed 2026-05-23 02:29 UTC · model grok-4.3

classification ✦ hep-ph
keywords charm quarksbottom quarksquark-gluon plasmadrag coefficientenergy lossnuclear modification factorFokker-Planckheavy-ion collisions
0
0 comments X

The pith

Expressing energy loss coefficients as polynomial expansions in momentum provides a flexible framework to test sensitivity to momentum dependence for charm and bottom quarks in the quark-gluon plasma.

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

This paper establishes that a phenomenological extension using polynomial expansions of the drag and energy loss coefficients in momentum offers a way to incorporate additional momentum dependence beyond simplified treatments. Numerical evolution of initial distribution functions via the Fokker-Planck equation then yields the nuclear modification factor R_AA for charm and bottom quarks in Pb-Pb collisions at 5.02 TeV, which is compared to recent ALICE and ATLAS data. A sympathetic reader would care because this tests how sensitive heavy-quark observables are to the momentum dependence in transport coefficients and enables more accurate separation of radiative and collisional energy loss.

Core claim

The central claim is that the phenomenological extension of the drag coefficient by polynomial expansions in momentum allows the effects of particle momentum on radiative and collisional energy loss to be determined more accurately, with the resulting R_AA values for charm and bottom quarks compared against experimental data.

What carries the argument

Polynomial expansions of the energy loss coefficients as functions of momentum, which provide a flexible phenomenological framework to test the sensitivity of heavy-quark observables to additional momentum dependence in transport coefficients.

If this is right

  • The nuclear modification factor R_AA of charm and bottom quarks can be computed with better accounting for momentum effects.
  • Radiative and collisional energy losses are determined more accurately.
  • Comparisons with ALICE and ATLAS data from Pb-Pb collisions at 5.02 TeV become feasible under this extended model.
  • Initial distributions are evolved using the Fokker-Planck equation under the momentum-dependent coefficients.

Where Pith is reading between the lines

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

  • This framework could be used to explore the minimal order of the polynomial needed to fit the data.
  • Application to other heavy flavors or energies might reveal patterns in the momentum dependence.
  • The approach may help bridge between phenomenological models and microscopic calculations of QGP transport.

Load-bearing premise

The polynomial expansion of the drag and energy-loss coefficients accurately captures the true momentum dependence of heavy-quark interactions in the QGP without introducing uncontrolled artifacts.

What would settle it

An observation that the R_AA predictions remain insensitive to the choice of polynomial order or functional form for the momentum dependence, or that they fail to improve agreement with data compared to constant-coefficient models, would falsify the utility of this extension.

Figures

Figures reproduced from arXiv: 2502.12333 by Fatemeh Taghavi-Shahri, Kurosh Javidan, Marjan Rahimi Nezhad.

Figure 1
Figure 1. Figure 1: FIG. 1 [PITH_FULL_IMAGE:figures/full_fig_p007_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2 [PITH_FULL_IMAGE:figures/full_fig_p009_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3 [PITH_FULL_IMAGE:figures/full_fig_p010_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4 [PITH_FULL_IMAGE:figures/full_fig_p011_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5 [PITH_FULL_IMAGE:figures/full_fig_p012_5.png] view at source ↗
read the original abstract

This paper investigates the influence of heavy-quark momentum on their interaction rate and the resulting drag coefficient in a quark-gluon plasma. To go beyond simplified treatments, we introduce a phenomenological extension of the drag coefficient by expressing the energy loss coefficients as polynomial expansions of momentum, thereby providing a flexible framework to test the sensitivity of heavy-quark observables to additional momentum dependence in transport coefficients. Furthermore, the effects of particle momentum on radiative and collisional energy loss are determined more accurately. The study focuses on calculating the nuclear modification factor ($R_{AA}$) of charm and bottom quarks in Pb-Pb collisions at $\sqrt{S_{NN}} = 5.02 \: TeV$. The initial distribution functions are evolved numerically using the Fokker-Planck equation. The results are compared with the latest experimental data from ALICE and ATLAS, collected in 2021 and 2022.

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

2 major / 0 minor

Summary. The manuscript investigates the momentum dependence of the drag coefficient for charm and bottom quarks in the QGP by introducing a phenomenological polynomial expansion of the energy-loss coefficients. It evolves initial heavy-quark distributions via the Fokker-Planck equation to compute the nuclear modification factor R_AA in 5.02 TeV Pb-Pb collisions and compares the results to ALICE and ATLAS data. The central claim is that the momentum-dependent extension supplies a flexible framework for testing sensitivity to additional momentum dependence in transport coefficients and thereby yields a more accurate determination of radiative versus collisional energy loss.

Significance. If the polynomial form is shown to capture genuine momentum dependence without uncontrolled artifacts, the approach would supply a practical tool for exploring how momentum-dependent transport coefficients affect heavy-quark suppression observables. The numerical Fokker-Planck evolution and direct comparison to recent experimental data are standard and reproducible elements of the analysis.

major comments (2)
  1. [Abstract] Abstract (phenomenological extension paragraph): the claim that the polynomial extension yields 'more accurate' energy loss is not accompanied by any explicit functional form, polynomial degree, basis choice, or fitting procedure. Without these specifications it is impossible to determine whether reported improvements arise from physical momentum dependence or from additional free parameters.
  2. [Abstract] Abstract (phenomenological extension paragraph): the drag and energy-loss coefficients are defined via polynomials whose coefficients are free parameters; any improvement in R_AA agreement is therefore achieved by fitting those coefficients to the same data used for validation, rendering the reported 'prediction' and 'more accurate determination' circular by construction.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the detailed review and constructive feedback on our manuscript. We address the major comments point by point below, agreeing that the abstract requires clarification and revision to avoid ambiguity.

read point-by-point responses
  1. Referee: [Abstract] Abstract (phenomenological extension paragraph): the claim that the polynomial extension yields 'more accurate' energy loss is not accompanied by any explicit functional form, polynomial degree, basis choice, or fitting procedure. Without these specifications it is impossible to determine whether reported improvements arise from physical momentum dependence or from additional free parameters.

    Authors: We agree that the abstract as written does not provide the necessary technical details. In the revised version we will specify the polynomial form (e.g., a quadratic expansion in momentum p), the degree chosen, the basis functions employed, and the procedure used to fix the coefficients (comparison with measured R_AA supplemented by theoretical guidance on the relative strength of radiative and collisional contributions). This will make clear that the extension is a controlled phenomenological parameterization rather than an ad-hoc addition of free parameters. revision: yes

  2. Referee: [Abstract] Abstract (phenomenological extension paragraph): the drag and energy-loss coefficients are defined via polynomials whose coefficients are free parameters; any improvement in R_AA agreement is therefore achieved by fitting those coefficients to the same data used for validation, rendering the reported 'prediction' and 'more accurate determination' circular by construction.

    Authors: The referee correctly identifies a potential circularity if the abstract is read as claiming an independent prediction. The manuscript's intent is to introduce a flexible parameterization that allows systematic exploration of momentum dependence, not to predict R_AA from first principles. We will revise the abstract to remove the phrasing 'more accurate determination' and 'prediction', replacing it with language that emphasizes the framework's utility for testing sensitivity to momentum-dependent transport coefficients. The body of the paper already compares constant versus momentum-dependent cases to illustrate the effect; the revised abstract will align with that exploratory scope. revision: yes

Circularity Check

0 steps flagged

No significant circularity detected; derivation remains self-contained.

full rationale

The paper introduces a phenomenological polynomial extension of the drag and energy-loss coefficients as a flexible framework, evolves initial distributions numerically via the Fokker-Planck equation to obtain R_AA for charm and bottom quarks, and compares the results to ALICE/ATLAS data. No quoted step in the provided text shows the polynomial coefficients being fitted to the same R_AA data used for validation, nor any self-definitional reduction, self-citation load-bearing premise, or renaming of a known result. The numerical evolution step is independent of the specific momentum-dependent form chosen, so the comparison to external data constitutes a genuine test rather than a tautology by construction.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 0 invented entities

The model rests on a domain assumption that the Fokker-Planck equation suffices and on free parameters introduced by the polynomial expansion; no new entities are postulated.

free parameters (1)
  • coefficients of the polynomial expansion of drag/energy-loss terms
    Introduced explicitly to make the drag coefficient momentum-dependent; their values are not fixed by first principles and must be chosen or fitted.
axioms (1)
  • domain assumption Fokker-Planck equation accurately describes the evolution of heavy-quark momentum distributions in an expanding QGP
    Invoked when the initial distributions are evolved numerically to obtain final spectra and R_AA.

pith-pipeline@v0.9.0 · 5693 in / 1386 out tokens · 29188 ms · 2026-05-23T02:29:05.587013+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 · 48 canonical work pages · 11 internal anchors

  1. [1]

    02 and β = 0

    for the c quark are found to be α = 0 . 02 and β = 0 . 05. These values indicate that both collisional and radiative e nergy losses of particles increase effec- tively at higher momentum. Specifically, the impact of parti cle’s momentum on radiative energy loss becomes more significant and must be taken into account. We also obtain k2 ≈ 1. 4 k1 through 10 /s...

  2. [2]

    Adams [STAR Collaboration], et al., Nucl

    J. Adams [STAR Collaboration], et al., Nucl. Phys. A 757 (2005) 102-183

  3. [3]

    Adcox [PHENIX Collaboration] et al., Nucl

    K. Adcox [PHENIX Collaboration] et al., Nucl. Phys. A 757 (2005) 184-283

  4. [4]

    Superdense matter: Neutrons or asym ptotically free quarks?

    J.C. Collins, M.J. Perry, “Superdense matter: Neutrons or asym ptotically free quarks?”, Phys. Rev. Lett. 34.1353 (1975)

  5. [5]

    Relativistic hydrodynamics for heavy-ion co llisions

    Jean-Yves Ollitrault,“Relativistic hydrodynamics for heavy-ion co llisions”, Eur. J. Phys. 29 (2008) 275-302

  6. [6]

    The Quark Gluon Plasma and relativistic hea vy ion collisions in the LHC era

    Francesco Becattini, “The Quark Gluon Plasma and relativistic hea vy ion collisions in the LHC era”, J. Phys. Conf. Ser. 527 (2014) 012012

  7. [7]

    Heavy-quark probes of the q uark-gluon plasma at RHIC

    H. van Hees, V. Greco, R. Rapp, “Heavy-quark probes of the q uark-gluon plasma at RHIC”, Phys. Rev. C 73 (2006) 034913

  8. [8]

    Heavy Quarks in the Quark–Gluon Plasma

    R. Rapp, H. van Hees, “Heavy Quarks in the Quark–Gluon Plasma” , Quark–Gluon Plasma 4, 111–206 (2010)

  9. [9]

    Mazumder, and T

    S. Mazumder, and T. Bhattacharyya, and J. Alam, and S. K Das, Phys. Rev. C 84 (2011) 044901

  10. [10]

    Jai Prakash, Mohammad Yousuf Jamal, Eur. Phys. J. Plus 139 (2024) 9, 778

  11. [11]

    Juhee Hong, Phys. Rev. C 109 (2024) 2, 024913 , e-Print: 2308.14530 [hep-ph]

  12. [12]

    Shuang Li, Jinfeng Liao , Eur. Phys. J. C 80 (2020) 7, 671 , e-Print: 1912.08965 [hep-ph]

  13. [13]

    The QCD collisional energy loss revised

    Peshier, A. ,“The QCD collisional energy loss revised”, Phys. Rev. Lett. 97 (2006) 212301, e-Print: hep-ph/0605294 [h ep-ph]

  14. [14]

    Energy loss of a heavy fe rmion in a hot plasma

    Braaten, Eric and Thoma, Markus H. “Energy loss of a heavy fe rmion in a hot plasma”, Phys. Rev. D 44 (1991) 1298-1310

  15. [15]

    Collisional energy loss of a fast heavy quark in a quark-gluon plasma

    Peigne, Stephane and Peshier, Andre, “Collisional energy loss o f a fast heavy quark in a quark-gluon plasma”, Phys. Rev. D 77 (2008) 114017 , e-Print: 0802.4364 [hep-ph]

  16. [16]

    J. D. Bjorken, Fermilab preprint Pub-82/59-THY (1982)

  17. [17]

    Braaten and Thoma, Markus H., Phys. Rev. D 44 (1991) 9, R2625

  18. [18]

    Armesto, C

    N. Armesto, C. A. Salgado and U. A. Wiedemann, Phys. Rev. D 69 (2004) 114003 , e-Print: hep- ph/0312106 [hep-ph]

  19. [19]

    Armesto, C

    N. Armesto, C. A. Salgado and U. A. Wiedemann, Phys. Rev. C 72 (2005) 064910 , e-Print: hep- ph/0411341 [hep-ph]

  20. [20]

    Non-Abelian Energy Loss at Finite Opacity

    M. Gyulassy, P. Levai and I. Vitev, Phys. Rev. Lett. 85 (2000) 5535 [nucl-th/0005032]

  21. [21]

    Heavy Quark Radiative Energy Loss in QCD Matter

    M. Djordjevic and M. Gyulassy, Nucl. Phys. A733 (2004) 265 [nucl-th/0310076]

  22. [22]

    Nuclear suppression of heavy quark production at forward rapidities in relativistic heavy ion collisions

    U. Jamil and D. K. Srivastava, J. Phys. G37 (2010) 085106, [arXiv: 1005.1208 [nucl-th]]

  23. [23]

    G. D. Moore and D. Teaney, Phys. Rev. C 71, 064904 (2005)

  24. [24]

    J. Alam, S. Raha and B. Sinha, Phys. Rev. Lett.73, 1895 (1994)

  25. [25]

    Svetitsky, Phys

    B. Svetitsky, Phys. Rev. D 37 (1988) 2484-2491

  26. [26]

    Heavy Quarks in the Quark-Gluon Plasma

    Ralf Rapp, Hendrik van Hees,“Heavy Quarks in the Quark-Gluon P lasma”, e-Print: 0903.1096 [hep-ph]

  27. [27]

    Impact of EMC effect on D meson modification factor in 15 equilibrating QGP

    Sheibani, J. and Javidan, K. and Mirjalili, “Impact of EMC effect on D meson modification factor in 15 equilibrating QGP”, Eur. Phys. J. Plus 137 (2022) 807

  28. [28]

    and Taghavi, R

    Modarres, M. and Taghavi, R. and Nik, R. Aminzadeh and Valesha badi, R. Kord. , Phys. Rev. D 104 (2021) 5, 056005

  29. [29]

    Olanj, M

    N. Olanj, M. Modarres, Nucl. Phys. A 998 (2020) 121735

  30. [30]

    On higher order and anisotropic hydrodynamics for Bjorken and Gubser flows

    C. Chattopadhyay, U. Heinz, S. Pal, G. Vujanovic, “Higher ord er and anisotropic hydrodynamics for Bjorken and Gubser flows”, Phys. Rev. C 97 (2018) 6, 064909 , arXiv:1801.07755 [nucl-th]

  31. [31]

    Constructing higher-order hydrodynamics: The third order

    S. Grozdanov, N. Kaplis, “Constructing higher-order hydrod ynamics: The third order”, Phys. Rev. D 93 (2016) 6, 066012 . , arXiv:1507.02461 [hep-th]

  32. [32]

    Ruggieri, F

    M. Ruggieri, F. Scardina, S. Plumari and V. Greco, Phys. Rev. C 89 (2014) 5, 054914

  33. [33]

    Resummation and Gauge Invariance of the Gluon Damping Rate in Hot QCD

    E. Braaten, R.D. Pisarski, “Resummation and Gauge Invariance of the Gluon Damping Rate in Hot QCD”, Phys. Rev. Lett. 64 (1990) 1338

  34. [34]

    Dragging Heavy Quarks in Quark Gluon Plasma at the Large Hadron Collider

    S.K. Das, J. Alam, P. Mohanty, “Drag of heavy quarks in Quark g luon plasma at the large hadron collider”, Phys. Rev. C 82 (2010) 014908 . , arXiv:1003.5508

  35. [35]

    Drag and Diffusion of Heavy Qua rks in a hot and anisotropic QCD medium

    P. K. Srivastava, B. K. Patra, “Drag and Diffusion of Heavy Qua rks in a hot and anisotropic QCD medium”, Eur. Phys. J. A 53, no. 6, 116 (2017)

  36. [36]

    Heavy Quark Diffusion with Relativistic Langevin Dynamics in the Quark-Gluon Fluid

    Y. Akamatsu, T. Hatsuda and T. Hirano, “Heavy Quark Diffusion with Relativistic Langevin Dynamics in the Quark-Gluon Fluid”, Phys. Rev. C 79 (2009) 054907

  37. [37]

    Dragging heavy quarks in qu ark gluon plasma at the large hadron collider

    S. K. Das, J.-e. Alam, P. Mohanty, “Dragging heavy quarks in qu ark gluon plasma at the large hadron collider”, Phys. Rev. C 82 (2010) 014908

  38. [38]

    Thermalization of heavy quarks in the qu ark-gluon plasma

    H. van Hees, R. Rapp, “Thermalization of heavy quarks in the qu ark-gluon plasma”, Phys. Rev. C 71 (2005) 034907

  39. [39]

    Rahimi Nezhad, and F

    M. Rahimi Nezhad, and F. Taghavi Shahri, and K. Javidan, Eur. Phys. J. Plus 139 (2024) 1, 92 , e- Print: 2307.06168 [hep-ph]

  40. [40]

    Wicks, W

    S. Wicks, W. Horowitz, M. Djordjevic and M. Gyulassy, Nucl. Phys. A784 (2007) 426 [nucl- th/0512076]

  41. [41]

    Energy loss of $B$ and $D$ mesons in PbPb collisions at $\sqrt{s_{NN}}$ = 2.76 TeV

    K. Saraswat, P. Shukla, V. Singh, PoS ICPAQGP2015 (2015) 082 , e-Print: 1511.05349 [hep-ph]

  42. [42]

    Energy loss o f heavy quarks and B and D meson spectra in PbPb collisions at LHC energies

    K. Saraswat, P. Shukla, V. Kumar, and V. Singh, “Energy loss o f heavy quarks and B and D meson spectra in PbPb collisions at LHC energies”, Nucl. Phys. A 961, 169-182 (2017)

  43. [43]

    ALICE Collaboration, Phys. Lett. B 820 (2021) 136558 , Report number: CERN-EP-2020-186

  44. [44]

    ALICE Collaboration, JHEP 01 (2022) 174 , Report number: CERN-EP-2021-213

  45. [45]

    ATLAS Collaboration, Phys. Lett. B 829 (2022) 137077 , e-Print: 2109.00411 [nucl-ex]

  46. [46]

    James and M

    F. James and M. Roos, Comput. Phys. Commun. 10, 343-367 (1975)

  47. [47]

    Shu-Qing Li, Wen-Jing Xing, Shanshan Cao, Guang-You Qin, Chin. Phys. C 48 (2024) 8, 084106

  48. [48]

    Y. L. Dokshitzer and D. E. Kharzeev, Phys. Lett. B 519 (2001) 199-206 , e-Print: hep-ph/0106202 [hep- ph]