REVIEW
Heavy-quark thermalization in the quark-gluon plasma via a fractional Langevin equation
Nuclear Theory · Submitted 19 Oct 2026 · 295 references
Pith's verdict: The fractional-FDT idea is genuinely new and the framework is salvageable, but the printed central formula has a sign error that makes D negative for every alpha in (1,2) and fails the stated alpha -> 1 limit.
T0 · AI review
Why Pith says this
Highlight any text to comment on it.
Referee report
The paper
The paper claims
This paper derives a fractional fluctuation–dissipation theorem for heavy quarks in the quark-gluon plasma: charm quarks thermalize to the same momentum for every fractional order α, while energy loss, suppression, and flow all grow with α.
The reading
What carries the argument
The load-bearing object is the Caputo fractional derivative of order $\alpha$, ${}^C D^\alpha_{0+} u(t) = \frac{1}{\Gamma(n-\alpha)} \int_0^t \frac{u^{(n)}(s)}{(t-s)^{1+\alpha-n}} \, ds$ for $1 < \alpha \le 2$, which installs a memory kernel into the equation of motion; paired with the characteristic time $\tau_p \approx 1/\gamma$, it replaces the first-order time derivative of the standard Langevin equation. The identity carrying the argument is the fractional fluctuation-dissipation theorem, obtained by taking the Fourier transform of the fractional Langevin equation, computing the equilibrium momentum variance from the transfer function $H(\omega) = [\gamma + \tau_p^{\alpha-1}(i\omega)^\alpha]^{-1}$ through the Wiener-Khintchine relation $S_p(\omega) = |H(\omega)|^2 S_\xi(\omega)$ with white-noise spectrum $S_\xi(\omega) = 2D$, and imposing $\langle p^2 \rangle_{\rm eq} = MT$. Numerically, the Caputo operator is discretized with a three-level L2 scheme whose history coefficients $b_j$ retain the entire preceding trajectory, because fractional evolution is nonlocal in time.
Extended reading notes
Core claim
Working from a fractional Langevin equation in which the Caputo derivative of order $\alpha$ ($1 < \alpha \le 2$) replaces the first-order time derivative, with a characteristic time $\tau_p$ restoring dimensional consistency, the paper imposes stationarity and the equipartition condition $\langle p^2 \rangle_{\rm eq} = MT$ and derives the fractional fluctuation-dissipation theorem $D(\alpha) = K_\alpha M T \gamma^{2-1/\alpha} \tau_p^{(\alpha-1)/\alpha}$, with $K_\alpha = \frac{\alpha}{2} \frac{\sin(\pi\alpha/2)\,\sin(\pi\alpha)}{\sin(\pi(\alpha-1)/2)}$, which reduces to the standard relation $D = \gamma M T$ as $\alpha \to 1$. Under this relation, $\langle p^2(t) \rangle$ reaches the same $\alpha$-independent plateau at late times, while $\langle x^2(t) \rangle$ grows as $t^{2\beta-1}$ and is therefore superdiffusive. The momentum distribution, mean energy loss, longitudinal momentum loss, transverse momentum broadening, nuclear suppression factor $R_{AA}$, and a path-length estimate of $v_2$ all acquire systematic $\alpha$-dependence: energy and longitudinal-momentum loss grow rapidly with $\alpha$ and their separation widens toward high $p_T$, high-$p_T$ suppression deepens, and for $\alpha \gtrsim 1.4$ the relaxation becomes nonmonotonic and oscillatory, with partial energy recovery after rapid degradation. The near-Markovian cases $\alpha = 1.05$ and $1.001$ reproduce the standard Langevin results for the same input parameters, which the authors use as the internal consistency check of the framework.
Load-bearing premise
The derivation assumes the thermal kicks on the heavy quark are delta-correlated white noise, $\langle \xi(t) \xi(t') \rangle = 2D \, \delta(t-t')$, with all memory carried by the fractional derivative kernel; if the plasma's noise is itself correlated in time, the derived drag-diffusion relation and every $\alpha$-dependent prediction built on it would change.
Editorial extensions
If this is right
- The standard drag-diffusion relation used across heavy-quark phenomenology is recovered as the $\alpha \to 1$ limit, so the conventional Langevin framework becomes a special case of a one-parameter family rather than a competing model.
- Thermal equilibrium cannot by itself reveal fractional dynamics: $\langle p^2(t) \rangle$ reaches the same plateau for every $\alpha$, so the discriminating signals are dynamical, namely the superdiffusive $\langle x^2(t) \rangle$, the early-time thermalization rate, and the shape of the momentum distribution.
- Energetic charm quarks are the most sensitive probe, since the separation between different $\alpha$ in energy loss and longitudinal momentum loss widens with initial $p_T$; high-$p_T$ suppression measurements therefore carry the strongest constraint on $\alpha$.
- For $\alpha \gtrsim 1.4$ the model predicts nonmonotonic, oscillatory relaxation in which energy can partially recover after rapid degradation, a qualitative signature that no Markovian description produces.
- After matching the overall suppression, the path-length estimate of the elliptic flow at intermediate $p_T$ rises by roughly a factor of two between the near-Markovian case and $\alpha = 1.8$, so flow and suppression data together constrain the fractional order more tightly than either alone.
Reading between the lines
- The pairing of a Caputo derivative with delta-correlated white noise is not the standard fractional-Langevin construction: canonical treatments the paper cites pair the fractional kernel with time-correlated fractional-Gaussian noise precisely so that the fluctuation-dissipation relation is consistent. The paper lists time-correlated noise as future work but does not assess whether the derived $D(
- The F-FDT also depends on $\tau_p$ whenever $\alpha \neq 1$, and the paper fixes $\tau_p$ from lattice and phenomenological inputs; extracting $\alpha$ from data will therefore be entangled with the value and temperature dependence of $\tau_p$, so a joint fit over both parameters is the natural measurement strategy.
- A testable corollary of the invariance result is that single-time equilibrium measurements cannot identify fractional dynamics, so the cleanest observables are transient ones: the early-time slope of $\langle p^2(t) \rangle$, the mean-squared-displacement exponent $2\beta-1$, and the nonmonotonic energy-loss pattern predicted for large $\alpha$.
- Because the near-Markovian limit is shown to reproduce existing standard-Langevin results, the practical question is not whether the fractional framework fits current $R_{AA}$ data but whether any observable can distinguish it from a standard model with a rescaled drag coefficient; the paper's own calibration step, which renormalizes the transport strength at each $\alpha$ before comparing $v_2$,
Original abstract
The dynamics of heavy quarks in the quark-gluon plasma is commonly described within a standard Langevin equation framework. In this work, we extend the conventional framework by using a fractional Langevin equation based on the Caputo derivative of order $\alpha$ with $1<\alpha\leq 2$. We solve the fractional Langevin equation and derive the corresponding {\textit {fractional fluctuation-dissipation theorem}} for the heavy quark thermalization in the quark-gluon plasma medium. We calculate the time evolution of the mean-squared momentum, mean-squared displacement, momentum distribution, energy loss, longitudinal momentum loss, transverse momentum broadening, and the nuclear suppression factor of the heavy quarks. The mean-squared momentum thermalizes to the same value for all $\alpha$ at late times, in agreement with the fractional fluctuation-dissipation theorem, and the nuclear suppression factor shows a systematic $\alpha$-dependence. The energy loss of an energetic heavy quark grows strongly with $\alpha$. The time evolution of the heavy quark momentum distribution provides additional sensitivity to the fractional dynamics of the heavy quarks in the QGP medium.
The compiled record
The paper's own sentences, arranged, kept beside the review and separate from it. Each label is drawn from a closed vocabulary; the number beside it is the probability the classifier gave that label. Numbers in the tables are the authors', checked by code. No sentence in this record was written by a model.
The abstract, sentence by sentence
The dynamics of heavy quarks in the quark-gluon plasma is commonly described within a standard Langevin equation framework. In this work, we extend the conventional framework by using a fractional Langevin equation based on the Caputo derivative of order $\alpha$ with $1<\alpha\leq 2$. We solve the fractional Langevin equation and derive the corresponding {\textit {fractional fluctuation-dissipation theorem}} for the heavy quark thermalization in the quark-gluon plasma medium. We calculate the time evolution of the mean-squared momentum, mean-squared displacement, momentum distribution, energy loss, longitudinal momentum loss, transverse momentum broadening, and the nuclear suppression factor of the heavy quarks. The mean-squared momentum thermalizes to the same value for all $\alpha$ at late times, in agreement with the fractional fluctuation-dissipation theorem, and the nuclear suppression factor shows a systematic $\alpha$-dependence. The energy loss of an energetic heavy quark grows strongly with $\alpha$. The time evolution of the heavy quark momentum distribution provides additional sensitivity to the fractional dynamics of the heavy quarks in the QGP medium.
background contribution method result interpretation main claim
| # | sentence | role | the authors' own work | states the main conclusion | strength |
|---|---|---|---|---|---|
| s0 | The dynamics of heavy quarks in the quark-gluon plasma is commonly described within a standard Langevin equation framework. | background 1.00 | 0.04 | 0.01 | asserted 0.99 |
| s1 | In this work, we extend the conventional framework by using a fractional Langevin equation based on the Caputo derivative of order $\alpha$ with $1<\alpha\leq 2$. | method 0.51 | 0.87 | 0.05 | asserted 0.98 |
| s2 | We solve the fractional Langevin equation and derive the corresponding {\textit {fractional fluctuation-dissipation theorem}} for the heavy quark thermalization in the quark-gluon plasma medium. | contribution 0.50 | 0.97 | 0.12 | asserted 1.00 |
| s3 | We calculate the time evolution of the mean-squared momentum, mean-squared displacement, momentum distribution, energy loss, longitudinal momentum loss, transverse momentum broadening, and the nuclear suppression factor of the heavy quarks. | contribution 0.66 | 0.96 | 0.06 | asserted 0.99 |
| s4 | The mean-squared momentum thermalizes to the same value for all $\alpha$ at late times, in agreement with the fractional fluctuation-dissipation theorem, and the nuclear suppression factor shows a systematic $\alpha$-dependence. | result 0.99 | 0.94 | 0.41 | asserted 1.00 |
| s5 | The energy loss of an energetic heavy quark grows strongly with $\alpha$. | result 1.00 | 0.94 | 0.21 | asserted 1.00 |
| s6 | The time evolution of the heavy quark momentum distribution provides additional sensitivity to the fractional dynamics of the heavy quarks in the QGP medium. | interpretation 0.85 | 0.84 | 0.13 | asserted 0.95 |
What was examined, and what was not
- arXiv listing abstract v1
- full text
- figures and tables
- references and citing sentences
- relations to other papers
- P abstract
- S sections
- T body sentences
- C citations
- R cross-paper relations
Figures
Showing 7 of the paper's 8 figures. See all figures in the PDF ↗ Figures reproduced from arXiv:2610.26086 by the authors.
References
Discussion
Disagree with the review, or know something it missed? Leave a comment, highlight any passage above to comment on it, or challenge this review formally.