Pith. sign in

REVIEW 3 major objections 4 minor 27 references

Intrinsic coupling between transverse spherocity and elliptic flow in heavy-ion collisions

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

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

Pith's one-line read In heavy-ion collisions, transverse spherocity is set by elliptic flow: the minimizing axis aligns with the event plane, giving S0 ≈ (1 − 2/3 v2)^2, so low-S0 events are high-flow events, not jet topologies.

desk verdict Qualitatively right and useful, but the quantitative centerpiece is only proven in an equal-pT limit that isn't the measured observable. read the letter →

arxiv 2607.21161 v1 pith:4HMRYHFV submitted 2026-07-23 nucl-th hep-ph

Subikash Choudhury , Ekata Nandy This is my paper · ORCID
classification nucl-thhep-ph PACS 25.75.-q25.75.Ld
keywords transversespherocityellipticfloweventshapeanisotropicplaneheavy-ioncollisionsengineeringquark-gluonplasma
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

The paper sets out to overturn the standard reading of transverse spherocity in heavy-ion collisions. In proton collisions a low spherocity value marks a jet-like, collimated event, and the same classifier has been carried over to heavy-ion collisions to separate jetty from isotropic events. Using simple toy simulations, a transport-model calculation, and an analytic continuum derivation, the authors argue that in high-multiplicity heavy-ion events the axis that minimizes spherocity locks onto the second-order event plane, making the observable a function of elliptic flow: S0 ≈ (1 − 2/3 v2)^2. If this is right, low-spherocity events are simply events with large collective anisotropy, and several published observations about such events — enhanced v2, stronger radial flow, modified constituent-quark-number scaling — should be attributed to flow rather than to jet enrichment.

What carries the argument

The central object is transverse spherocity, an event-shape variable S0 = (π^2/4)[min_n Σ_i |p_{T,i} × n| / Σ_i |p_{T,i}|]^2 that measures how collimated an event's transverse momentum is. The analytical machinery is the exact Fourier expansion of |sin x|, which contains only even cosine harmonics; when the particle azimuthal distribution is written in the standard flow-harmonic series, orthogonality kills all odd harmonics, and the minimization functional F(ψ) depends only on even-order flow coefficients. Retaining the dominant v2 term, the minimum of F(ψ) occurs at ψ = Ψ2, yielding S0 ≈ (1 − 2/3 v2)^2. This identity is what converts spherocity from a jet-topology classifier into a flow obs

What would settle it

Measure, event-by-event, the angle between the axis that minimizes transverse spherocity and an independently reconstructed second-order event plane, in events where v2 is known to exceed several percent; if the distribution of that angle is not peaked at zero, the central claim fails. A complementary test is to compute S0 with pT weights turned off and on in the same simulated events: the formula (1 − 2/3 v2)^2 must hold in both cases for the paper's derivation to be complete.

Watch

Extended reading notes

Core claim

The central claim, on the paper's own terms, is that transverse spherocity is intrinsically coupled to elliptic flow in heavy-ion collisions: in the continuum limit the unit vector n that minimizes the spherocity functional aligns exactly with the second-order symmetry plane, ψ_min = Ψ2, so that F_min = (2/π)(1 − 2/3 v2) and S0 ≈ (1 − 2/3 v2)^2. The same mechanism creates an inherent anti-correlation between S0 and v2 — events with larger elliptic anisotropy automatically have smaller spherocity — independent of any jet activity. The authors conclude that, in heavy-ion collisions, spherocity should be interpreted as a probe of collective momentum-space anisotropy, and that earlier interpreta

Load-bearing premise

The derivation assigns equal weight to every particle and ignores the pT-dependence of flow; if the actual pT-weighted spherocity responds differently to high-momentum particles, the clean S0 ≈ (1 − 2/3 v2)^2 relation may shift.

Editorial extensions

If this is right

  • Earlier claims that low-spherocity heavy-ion events are jet-enriched — including enhanced elliptic flow, stronger radial expansion, and modified constituent-quark-number scaling — should be re-expressed as manifestations of large event-by-event v2.
  • S0-based event selection and Event Shape Engineering based on the reduced flow vector q2 should yield similar physics trends, because both are tied to the same symmetry plane.
  • The spherocity-minimizing axis itself can serve as an event-plane estimator in high-multiplicity events.
  • In small systems such as pp collisions, the same observable remains a topological classifier, so its physical meaning is genuinely system-dependent.

Reading between the lines

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

  • Inference: Because only even-order harmonics enter the Fourier expansion of |sin x|, odd harmonics such as v3 should leave S0 essentially untouched, suggesting S0 could select on v2 while filtering out triangular-flow fluctuations.
  • Inference: The derivation's equal-weight assumption could be tested directly by looking at the pT dependence: if S0 is computed with high-pT particles only, the relation may deviate; a data-driven parametrization of S0(v2, pT) would quantify how much jet-like hard particles can bias the axis.
  • Inference: If the relation holds, a natural follow-up is to check whether spherocity-selected events exhibit the same softening of the particle spectrum as q2-selected events; the paper implies they should.
  • Inference: One might exploit S0's global long-range correlation in heavy-ion collisions as a cheap way to correlate two η-separated detectors without a full event-plane reconstruction.
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.

Referee Report

3 major / 4 minor

Summary. The paper argues that transverse spherocity S0 in heavy-ion collisions is not primarily a jet-topology classifier but is intrinsically tied to elliptic flow v2. The authors support this with toy Monte Carlo events, AMPT calculations, and a continuum analytic derivation. In the equal-weight continuum limit, the minimizing axis of the spherocity functional aligns with the second-order symmetry plane Ψ2, and the minimal value gives S0 ≈ (1 - 2/3 v2)^2, so low-S0 events are naturally events with larger v2. The paper further shows a strong correlation between S0 values measured in two separated η intervals for AMPT Pb-Pb events, while PYTHIA pp events show no such correlation, supporting the global, flow-driven nature of S0 in heavy-ion collisions. The authors conclude that recent observations based on S0-selected events should be interpreted via collective anisotropic flow rather than jet enrichment.

Significance. If the central result holds, it provides a simple analytical explanation for the observed anti-correlation between S0 and v2 and reframes the interpretation of a widely used event-shape observable. The analytical derivation is transparent, the toy model is instructive, and the AMPT comparison is a useful independent check. The paper also makes a falsifiable prediction: for flow-dominated events, the minimizing spherocity axis should coincide with the event plane, which is consistent with the presented simulations. The main significance is conceptual: it warns the heavy-ion community against treating S0 as a clean jet/isotropy classifier in AA collisions. The quantitative relation, however, rests on an unvalidated equal-pT approximation and needs strengthening before the conclusions can be accepted at face value.

major comments (3)
  1. [§III, Eq. (5)–(21)] The derivation of the central result is performed under the explicit assumption p_T,i=1, while the measured observable (Eq. 1) is p_T-weighted. In the continuum limit the relevant distribution is W(φ) ∝ ∫ p_T dN/(dp_T dφ) dp_T, whose second harmonic v2^eff = ⟨p_T cos2(φ−Ψ2)⟩/⟨p_T⟩ is not equal to the usual v2 when v2(p_T) varies with p_T. Thus Eq. (21) is a relation between S0 and a p_T-weighted harmonic, not necessarily the inclusive v2 quoted in the abstract and conclusions. The paper does not prove that ψmin = Ψ2 remains valid for W(φ), nor does it compare Eq. (21) with p_T-weighted toy/AMPT calculations. This gap is load-bearing because Eq. (21) is the quantitative basis for the flow-selection interpretation. I recommend deriving the p_T-weighted version or benchmarking Eq. (21) directly against p_T-weighted simulations.
  2. [§III, Eq. (12)–(21)] The minimization in Eq. (12) retains only the m=1 term, and Eq. (21) additionally neglects v4, v6, ... in Eq. (18). Higher even harmonics contribute to F(ψ) with coefficients v_{2m}/(4m^2−1); when v4 is not negligible (e.g., in ultra-central collisions or low-multiplicity selections), the correction can shift the minimum away from Ψ2, especially if Ψ4 differs from Ψ2. A quantitative estimate of this bias, or a demonstration from the AMPT/toy models that these corrections are negligible for the kinematic range used, is needed before Eq. (21) is used to reinterpret published S0-based results.
  3. [§II.B and §III.A] The toy model and the AMPT test fix the reaction-plane angle at ΨRP = π/4 and show that the distribution of the minimizing direction peaks there. This establishes alignment on average, but not that ψmin tracks the per-event Ψ2 in realistic collisions with event-by-event fluctuating initial geometry. Section III.A uses the existence of a common per-event symmetry plane to explain the η-correlation; this would be strengthened by a direct AMPT check correlating ψmin with the per-event Ψ2 (or q2 vector) event-by-event, rather than only with a globally fixed input angle.
minor comments (4)
  1. [Eq. (1)] The denominator is typeset as `Σ_i ⃗|pT,i|`; this should be `Σ_i |p_T,i|`.
  2. [§III] The phrase 'we first consider equal particle weights' suggests the equal-pT assumption will be lifted later, but the paper never returns to the p_T-weighted case. This should be stated explicitly at the end of Sec. III, and the limitation should be acknowledged.
  3. [Fig. 3] The AMPT figure would benefit from error bars or a statement of statistical uncertainty, and the number of events used should be specified.
  4. [Ref. [26]] The reference to the PYTHIA8 online manual is not ideal; the standard PYTHIA8 paper (Sjöstrand et al., Comput. Phys. Commun. 191, 159 (2015)) should be cited instead.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: Eq. (21) is a derived kinematic identity, not an input masquerading as a prediction.

full rationale

The paper's central relation, S0 ≈ (1 − (2/3)v2)^2, is obtained by substituting two standard Fourier expansions — the azimuthal flow distribution (Eq. 7) and |sin x| (Eq. 8) — into the continuum equal-weight spherocity functional (Eq. 6). The elliptic flow coefficient v2 is not fitted to reproduce S0; rather, both S0 and v2 are moments of the same underlying distribution, and their relation is a derived mathematical consequence of the definitions. The toy Monte Carlo and AMPT comparisons serve as external illustrations and independent model checks, not as inputs that force Eq. (21). No load-bearing self-citations are present: references [21-24,27] are other groups' measurements and are reinterpreted rather than invoked as proof, and refs [25,26] are standard public codes. The only notable weakness is that the analytic derivation explicitly assumes equal particle weights (pT,i = 1), whereas the experimentally measured spherocity in Eq. (1) is pT-weighted; this is an approximation/scope limitation and a possible quantitative gap, but it is not circular, because the paper does not define S0 to equal the v2 formula. The derivation is self-contained and does not reduce to its inputs by construction.

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

The derivation rests on standard Fourier/flow formalism plus two unverified approximations: the equal-weight continuum limit and v2 dominance. No new entities are postulated and no values are fitted to data; the toy-model numbers are illustrative.

assumptions (6)
  • domain assumption The azimuthal distribution can be treated as a continuous normalized probability density (high-multiplicity continuum limit).
    Used in Section III Eq. 4 to convert the discrete spherocity sum to an integral; requires large multiplicity and neglects finite-N fluctuations.
  • ad hoc to paper All particles carry equal transverse momentum (pT,i = 1) in the analytical derivation.
    Introduced in Section III 'we first consider equal particle weights i.e, pT,i=1' to isolate geometry; never generalized to the pT-weighted spherocity of Eq. 1.
  • domain assumption Elliptic flow dominates higher harmonics (v2 >> v4 > v6 ...).
    Used after Eq. 12 to drop m >= 2 terms; reasonable in non-central AA collisions but not quantified for central or strongly fluctuating events.
  • standard math The flow Fourier expansion P(phi) = (1/2π)[1 + 2 Σ vn cos n(phi - Ψn)] is the correct description of the final-state distribution.
    Definition of harmonic flow used in Eq. 7; standard in the field.
  • domain assumption The AMPT transport model adequately represents the relevant dynamics of Pb-Pb collisions.
    Used in Section II to validate the toy-model conclusion (Fig. 3).
  • domain assumption pT and azimuthal anisotropy factorize (⟨pT cos 2(φ - Ψ2)⟩ ≈ ⟨pT⟩⟨cos 2(φ - Ψ2)⟩).
    Invoked in Appendix A Eq. A5 for the analogous sphericity relation; not needed for the main spherocity derivation.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Intrinsic coupling between transverse spherocity and elliptic flow in heavy-ion collisions." pith.science (2026). https://pith.science/paper/4HMRYHFV

@misc{pith2026260721161,
  author       = {Pith},
  title        = {Pith review of: Intrinsic coupling between transverse spherocity and elliptic flow in heavy-ion collisions},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/4HMRYHFV}},
  note         = {Machine review of arXiv:2607.21161}
}
abstract

Transverse spherocity ($S_{0}$) is an event-shape observable widely used to classify collision events according to their topology, particularly to distinguish jet-like from isotropic events. Low-spherocity events are generally interpreted as being associated with enhanced jet activity. This event classification has also been applied to heavy-ion collisions to investigate the influence of event topology on several observables, including elliptic flow and constituent-quark-number scaling. In this work, we demonstrate that such an interpretation requires careful reconsideration. Using toy Monte Carlo simulations, A Multiphase Transport (AMPT) model calculations, and an analytical formulation of the spherocity observable, we show that transverse spherocity is intrinsically related to the elliptic flow coefficient, $v_{2}$. This connection arises because the axis that minimizes the spherocity aligns with the event symmetry plane, causing events with larger elliptic anisotropy to naturally exhibit smaller spherocity values even in the absence of genuine jet-like topologies. We further show that this intrinsic relation gives rise to an inherent anti-correlation between transverse spherocity and $v_{2}$, implying that several characteristics previously attributed to the enhanced jet-like nature of low-spherocity events in heavy-ion collisions can instead be understood as consequences of collective anisotropic flow. Our results indicate that, in heavy-ion collisions, transverse spherocity should be interpreted primarily as a probe of the collective momentum-space anisotropy rather than as a direct measure of jetty event topology. Consequently, physics conclusions drawn from spherocity-selected events should explicitly account for its intrinsic correlation with elliptic flow.

Figures

Figures reproduced from arXiv: 2607.21161 by the authors.

Figure 1
Figure 1. FIG. 1: [Color online] (Left) Azimuthal ( [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. FIG. 2: [Color online] (Left panel) Flow modulated [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗
Figure 3
Figure 3. FIG. 3: Phase angle distribution of unit vector obtained from events simulated with AMPT [25] model for Pb-Pb collisions in [PITH_FULL_IMAGE:figures/full_fig_p004_3.png] view at source ↗
Figures from the paper (1 more)
Figure 4
Figure 4. Figure 4: FIG. 4: Correlation between self normalized spherocity calculated in two disjoint [PITH_FULL_IMAGE:figures/full_fig_p007_4.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

27 extracted references · 2 linked inside Pith

  1. [1]

    Takahashi et al., Phys

    J. Takahashi et al., Phys. Rev. Lett. 103, 242301 (2009)

    2009

  2. [2]

    Alver and G

    B. Alver and G. Roland, Phys. Rev. C 81, 054905 (2010) [Erratum-ibid. C 82, 039903 (2010)]

    2010

  3. [3]

    Pang, G.-Y

    L.-G. Pang, G.-Y. Qin, V. Roy, X.-N. Wang, and G.-L. Ma, Phys. Rev. C 91, 044904 (2015)

    2015

  4. [4]

    Aamodt et al

    K. Aamodt et al. [ALICE Collaboration], Phys. Rev. Lett. 107, 032301 (2011)

    2011

  5. [5]

    Adare et al

    A. Adare et al. [PHENIX Collaboration], Phys. Rev. Lett. 107, 252301 (2011)

    2011

  6. [6]

    Aad et al

    G. Aad et al. [ATLAS Collaboration], Phys. Rev. C 86, 014907 (2012)

    2012

  7. [7]

    Aaboud et al

    M. Aaboud et al. [ATLAS Collaboration], Eur. Phys. J. C 78, 142 (2018)

    2018

  8. [8]

    Voloshin, Phys

    Jurgen Schukraft, Anthony Timmins, Sergei A. Voloshin, Phys. Lett. B 719, 394 (2013)

    2013

Show all 27 references
  1. [9]

    Aad et al

    G. Aad et al. [ATLAS Collaboration], Phys. Rev. C 92, 034903 (2015)

    2015

  2. [10]

    Adam et al

    J. Adam et al. [ALICE Collaboration], Phys. Rev. C 93, 034916 (2016)

    2016

  3. [11]

    Banfi, G

    A. Banfi, G. P. Salam, G. Zanderighi, J. High Energ. Phys. 2010, 38 (2010). 9

    2010

  4. [12]

    Movilla Fern´ andez, O

    P.A. Movilla Fern´ andez, O. Biebel, S. Bethke, S. Kluth, P. Pfeifenschneider and the JADE Collaboration, Eur. Phys. J. C 1, 461 (1998)

    1998

  5. [13]

    Dasgupta, G

    M. Dasgupta, G. P. Salam, J. Phys. G 30 R143 (2004)

    2004

  6. [14]

    The OPAL Collaboration, Eur. Phys. J. 40, 287 (2005)

    2005

  7. [15]

    Bethke, S

    S. Bethke, S. Kluth, C. Pahl, J. Schieck and the JADE Collaboration, Eur. Phys. J. C 64, 351 (2009)

    2009

  8. [16]

    Abbiendi et al

    G. Abbiendi et al. [The OPAL Collaboration], Eur. Phys. J. C 71, 1733 (2011)

    2011

  9. [17]

    Cuautle, R

    E. Cuautle, R. Jimenez, I. Maldonado, A. Ortiz, G. Paic and E. Perez, arXiv:1404.2372 [hep-ph]

    arXiv

  10. [18]

    Ortiz, G

    A. Ortiz, G. Paic and E. Cuautle, Nucl. Phys. A 941, 78 (2015)

    2015

  11. [19]

    Abelev et al

    B. Abelev et al. [ALICE Collaboration], Eur. Phys. J. C 72, 2124 (2012)

    2012

  12. [20]

    Acharya et al

    S. Acharya et al. [ALICE Collaboration], Eur. Phys. J. C 79, 857 (2019)

    2019

  13. [21]

    Mallick, R

    N. Mallick, R. Sahoo, S. Tripathy and A. Ortiz, J. Phys. G 48, 045104 (2021)

    2021

  14. [22]

    Prasad, N

    S. Prasad, N. Mallick, S. Tripathy and R. Sahoo, [arXiv:2207.12133 [hep-ph]]

    arXiv

  15. [23]

    Mallick, S

    N. Mallick, S. Tripathy and R. Sahoo, Eur. Phys. J. C 82, 524 (2022)

    2022

  16. [24]

    Sarkar, P

    S. Sarkar, P. Mali, A. Mukhopadhyay, Eur. Phys. J. A 58, 139 (2022)

    2022

  17. [25]

    Z. W. Lin, C. M. Ko, B. A. Li, B. Zhang and S. Pal, Phys. Rev. C 72, 064901 (2005)

    2005

  18. [26]

    Pythia8 Online Manual https://pythia.org/manuals/pythia8215/Welcome.html

  19. [27]

    Prasad, N

    S. Prasad, N. Mallick, D. Behera, R. Sahoo and S. Tripathy, Sci. Rep. 12, 3917 (2022). Appendix A: Relation between transverse sphericity and elliptic flow In this appendix, we derive the relation between transverse sphericity and elliptic flow under the factorization approxim...

    2022

Pith tools

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