REVIEW 3 major objections 4 minor 8 references
Investigating quarkonium collectivity in heavy-ion collisions
T0 review · 3 major / 4 minor · reviewed 2026-08-15 · deepseek-v4-flash
Pith's one-line read This paper reports the first quarkonium elliptic-flow measurements from 5.36 TeV Pb–Pb collisions, including a 2.7σ deficit of the four-particle cumulant that may signal event-by-event J/ψ flow fluctuations.
desk verdict The new Run 3 quarkonium v2 data are genuinely useful, but the paper's only new physics claim — a 2.7σ sensitivity to J/ψ flow fluctuations — is load-bearing and not yet supported by the analysis details. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The central objects are Fourier flow harmonics $v_2\{EP\}$, $v_2\{SP\}$, $v_2\{2\}$, and $v_2\{4\}$, extracted from azimuthal correlations between dimuons at forward rapidity and charged particles or the event plane. The multi-particle cumulant technique, applied here for the first time to $J/\psi$ flow, provides the new $v_2\{2\}$ and $v_2\{4\}$ observables. The interpretative machinery is the Gaussian-fluctuation assumption that relates $v_2\{SP\}^2$ and $v_2\{4\}^2$ to the mean flow $\langle v_2\rangle$ and its event-by-event variance $\sigma_v^2$; the observed $v_2\{4\}<v_2\{SP\}$ deficit is what carries the claim of sensitivity to $J/\psi$ flow fluctuations.
What would settle it
A measurement with larger statistics that does not reproduce the $v_2\{4\}$ deficit in a given centrality or $p_T$ bin, or a demonstration that residual non-flow survives the $|\Delta\eta|>1.7$ gap, would falsify the claim of strong sensitivity to $J/\psi$ flow fluctuations.
Extended reading notes
Core claim
The central claim is that quarkonia participate in the collective flow of the quark-gluon plasma to a degree that can now be quantified with Run 3 statistics. The $J/\psi$ elliptic flow measured by the event-plane, scalar-product, and two-particle cumulant methods agree with one another and with Run 2 results, establishing a positive anisotropy from $p_T=0$ up to about 10 GeV/$c$ in several centrality classes. The new four-particle cumulant result lies below $v_2\{SP\}$ with $2.7\sigma$ significance, which the paper takes, under the assumption of Gaussian flow fluctuations, as evidence for event-by-event fluctuations in $J/\psi$ flow. The $\Upsilon(1S)$ measurement is compatible with zero and is compared with two theoretical models; the data agree qualitatively but cannot yet distinguish the models.
Load-bearing premise
The interpretation of the observed $v_2\{4\}<v_2\{SP\}$ deficit as flow fluctuations depends on the assumption that event-by-event elliptic flow is Gaussian-distributed and that the $|\Delta\eta|>1.7$ gap removes all non-flow contributions; if either fails, the central claim collapses.
Editorial extensions
If this is right
- If the $2.7\sigma$ deficit is real, higher-statistics Run 3 data will map $J/\psi$ flow fluctuations as a function of centrality and $p_T$.
- The cross-method agreement among $v_2\{EP\}$, $v_2\{SP\}$, and $v_2\{2\}$ indicates that the observed $J/\psi$ anisotropy is not an artifact of a single extraction method.
- The positive low-$p_T$ $J/\psi$ $v_2$ strengthens the regeneration scenario for charmonium in the quark-gluon plasma and will tighten model constraints on charm-quark diffusion and recombination.
- The $\Upsilon(1S)$ $v_2$ being consistent with zero constrains beauty-quark thermalization and bottomonium regeneration, even if current uncertainties do not allow firm conclusions.
Reading between the lines
- The same cumulant comparison could be made for charged hadrons in the same events; the ratio of the inferred flow variance for $J/\psi$ to that for light hadrons would test whether heavy quarks inherit the full medium fluctuation spectrum or are diluted by their hard-production vertex.
- If the $v_2\{4\}$ deficit grows with centrality in the same way as initial-state eccentricity fluctuations, it would support a hydrodynamic origin; if it is centrality-independent, it would favour initial-state correlation mechanisms of the kind the paper cites.
- A future $\Upsilon$ $v_2\{4\}$ measurement with the same technique could sharpen the contrast between the near-zero $\Upsilon$ $v_2$ and the positive $J/\psi$ signal, testing whether bottom quarks are merely less thermalized or also less fluctuating.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. This ALICE preliminary proceedings paper reports new Run 3 Pb-Pb measurements at sqrt(s_NN)=5.36 TeV of J/psi elliptic flow at forward rapidity (2.5<y<4) using event-plane, scalar-product, and, for the first time, multiparticle cumulant methods, together with a first Upsilon(1S) v2 measurement. The J/psi v2{SP} centrality and pT dependence is found consistent with Run 2. The new quantitative claim is a 2.7 sigma deficit of v2{4} relative to v2{SP} in 10-30% centrality, which the paper interprets as the first strong sensitivity to J/psi flow fluctuations under a Gaussian fluctuation assumption. The Upsilon(1S) v2 is consistent with zero and is compared qualitatively with the TAMU and BBJS models.
Significance. If the central fluctuation claim survives scrutiny, the paper would provide the first experimental indication of event-by-event fluctuations in quarkonium elliptic flow, a genuinely new constraint on quarkonium production and thermalization in the quark-gluon plasma. The comparison of J/psi, Upsilon(1S), and open heavy-flavour v2 values across rapidity also adds useful model discrimination. The paper is strengthened by the use of multiple independent flow methods, explicit statistical and systematic uncertainty bars on the figures, and direct comparison to published theoretical models without introducing free parameters or fitting model parameters.
major comments (3)
- [Results, Fig. 3] The paper's central new claim, the first sensitivity to J/psi flow fluctuations, rests entirely on the 2.7 sigma deficit of v2{4} relative to v2{SP} in the 10-30% centrality class. This interpretation requires two unverified assumptions: (i) that event-by-event elliptic flow fluctuations are Gaussian so that Eqs. (1)-(2) apply, and (ii) that the |Delta eta| > 1.7 gap removes all relevant non-flow correlations. The paper provides no test of Gaussianity and no quantification of residual non-flow. The reference particles for v2{SP} and v2{EP} (FT0 forward detectors) differ from the midrapidity charged particles used for v2{2} and v2{4}, so the two sets can carry different non-flow and resolution corrections. Since the quoted 2.7 sigma is statistical, the observed deficit could be a method difference rather than a flow-fluctuation signal. I ask the authors to quantify or otherwise bound these effects before making the fluctuation claim.
- [Results, Fig. 3] The systematic uncertainties shown in the figures are never defined. The paper does not describe how the systematic uncertainties were estimated, how the event-mixing normalization was performed, or how the event-plane resolution correction was applied to the v2{EP} and v2{SP} results. More importantly, no systematic uncertainty is quoted for the 2.7 sigma significance of v2{4} < v2{SP}, and no correlation between the v2{4} and v2{SP} systematic uncertainties is given. Without this information the reader cannot evaluate whether the central deficit is robust. A short description of the systematic procedure and a table of results with combined uncertainties would resolve this.
- [Results, Eq. (2)] Equation (2) as written is dimensionally inconsistent: it states v2{4} = <v2>^2 - sigma_v^2, whereas the left-hand side is a flow coefficient and the right-hand side has dimensions of squared flow. The standard relation is v2{4}^2 = <v2>^2 - sigma_v^2. If the authors intend v2{4} to denote the squared quantity, this should be stated explicitly; otherwise the formula is a technical error that obscures the fluctuation interpretation.
minor comments (4)
- [Results, Fig. 2] The left panel legend lists centrality classes 10-30% and 30-50%, while the caption says the panel shows 0-10%, 10-30%, and 30-50%; please make the legend and caption consistent.
- [Results, Fig. 4] The right panel caption contains an garbled pT range ('c < 15 GeV/Tp') and the left panel caption also has formatting issues; please correct these to the intended 'pT < 15 GeV/c' and '2.5 < y < 4'.
- [Results, p. 3] The phrase 'a strong sensitivity to J/psi flow fluctuations' is stronger than a 2.7 sigma statistical deficit warrants; the authors may wish to soften this to 'first evidence' or 'first indication' unless a full systematic treatment changes the significance.
- [References] Reference 2 combines a 2008 JINST detector paper with a 2023 arXiv identifier that appears unrelated; please verify the citation and provide the correct arXiv number for the ALICE detector description.
Circularity Check
No circularity: the paper reports new measurements using standard flow methods and cites external relations; no fitted input is renamed as a prediction.
full rationale
The paper is an experimental conference report presenting Run 3 quarkonium v2 measurements. No parameter is fitted to the data and then presented as a prediction; the flow coefficients are extracted with standard, separately established methods (event plane, scalar product, multi-particle cumulants). The Gaussian-fluctuation relation used to interpret the difference v2{4} < v2{SP} is quoted from CMS (Eq. 1 and Eq. 2) and is not derived from or fitted to the present data; the 2.7σ deficit is a measured difference, not a constructed output. The agreement of v2{EP}, v2{SP}, and v2{2} is an empirical consistency check, not an input to the cumulant extraction. The auxiliary assumptions (Gaussian fluctuations, non-flow suppression by |Δη| > 1.7) are explicitly stated conditions for the physics interpretation; if they fail, the conclusion would weaken, but the derivation is not circular because the observable values are measured independently of those assumptions. The ALICE self-references are to the detector description and to preliminary plots of the present data, not to a uniqueness theorem, a fitted parameter, or a prior result that is load-bearing for the new claim. Theoretical model comparisons (TAMU, BBJS) are external and are not used to tune the measured values. The paper also explicitly acknowledges limitations: no definitive conclusions beyond 10 GeV/c due to statistical uncertainties, and no firm conclusions for Υ(1S) v2. No step in the derivation chain reduces by construction to its own input, and no self-citation is invoked to forbid alternatives. Overall circularity score is 0.
Assumptions & free parameters
assumptions (4)
- domain assumption Event mixing reproduces the combinatorial dimuon background in Pb-Pb collisions.
- domain assumption Non-flow correlations are suppressed by the pseudorapidity gap |Delta eta| > 1.7.
- domain assumption Event-by-event elliptic flow fluctuations are Gaussian, so v2{2} and v2{4} are related through the variance sigma_v^2.
- domain assumption Charged particles at midrapidity and dimuons at forward rapidity share the same underlying flow plane.
Cite this review
Pith. "Pith review of Investigating quarkonium collectivity in heavy-ion collisions." pith.science (2026). https://pith.science/paper/VGLSINXN
@misc{pith2026250513892,
author = {Pith},
title = {Pith review of: Investigating quarkonium collectivity in heavy-ion collisions},
year = {2026},
howpublished = {\url{https://pith.science/paper/VGLSINXN}},
note = {Machine review of arXiv:2505.13892}
}
abstract
The quark-gluon plasma (QGP) produced in ultrarelativistic heavy-ion collisions has exhibited properties of a mostly perfect fluid. These properties can be observed through the hydrodynamic expansion of the QGP. Experimentally, this was established by measuring azimuthal anisotropies in the final state, known as elliptic flow ($v_2$) or higher order harmonics such as triangular flow ($v_3$). These Fourier harmonic coefficients have been extensively measured in past experiments using inclusive charged particles or identified particles in the soft sector. Interestingly, measuring such coefficients using hard probes, such as quarkonia, brings additional information about heavy-quarks production and thermalization in the QGP. In this study, we investigate quarkonia collectivity using Run 3 data collected in 2023, presenting new flow measurements. We employ different experimental methods to extract flow coefficients, including the scalar product, event plane, and cumulant methods. These results will impose new constraints on theoretical models, enhancing our understanding of quarkonia behavior in heavy-ion collisions.
Reference graph
Works this paper leans on
-
[1]
Jean-Yves Ollitrault, Phys. Rev. D , 1992, 46, 229–245. https://doi.org/10.1103/ PhysRevD.46.229
work page 1992
-
[2]
ALICE Collaboration, JINST 3,S08002, (2008) https://doi.org/10.48550/arXiv. 2307.14084
-
[3]
Xiaojian Du, Ralf Rapp https://doi.org/10.48550/arXiv.1504.00670
-
[4]
Ante Bilandzic et al., https://doi.org/10.1103/PhysRevC.83.044913
-
[5]
CMS collaboration, JHEP05(2023)007 https://doi.org/10.48550/arXiv.2205.00080
-
[6]
Xiaojian Du, et al. Phys. Rev. C. 96, 054901 (2017), https://doi.org/10.1103/ PhysRevC.96.054901
work page 2017
-
[7]
Partha Pratim Bhaduri, Nicolas Borghini, Amaresh Jaiswal, and Michael Strickland, Phys. Rev. C. 100, 051901(R) (2019), https://doi.org/10.1103/PhysRevC.100.051901
- [8]
Reviewed August 15, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.