{"id":"b180b4f9-16b5-4f16-9651-60797e283c6f","arxiv_id":"2607.21161","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"In heavy-ion collisions the axis minimizing transverse spherocity aligns with the elliptic-flow plane, so S0 approximately equals (1 - 2/3 v2)^2 and low-S0 events are high-v2 events rather than jet-rich events.","lead":"Transverse spherocity is not an independent measure of jet-like topology in heavy-ion collisions: the direction it picks automatically lines up with the collective elliptic-flow plane, so low-spherocity events are just high-flow events. Because many analyses use spherocity to select jet-like events, this reframing affects how those published results should be read.","discovery_kind":"first_principles","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Equal-pT derivation of Eq. (21) is not shown to carry over to the pT-weighted spherocity actually measured.","rationale":"The reader's weakest assumption identifies the same gap that I find most load-bearing: the analytic derivation of Eq. (21) explicitly starts from pT_i = 1 and never justifies extension to the pT-weighted observable defined in Eq. (1). The central claim as stated includes the quantitative formula S0 ≈ (1 − 2/3 v2)^2 and the inference that low-S0 events are flow-selected. That inference depends on the direction alignment (ψmin = Ψ2) and on the monotonic anti-correlation; both are likely robust to pT weighting, as the AMPT result with fixed reaction plane supports alignment. However, the exact numerical relation may change because v2(pT) weighting modifies the effective harmonic content. This does not overturn the qualitative conclusion, but it makes the paper's headline quantitative prediction improperly qualified. The reader's CONDITIONAL verdict already captures this; my stress test does not move the verdict, so I recommend UNCHANGED. I also note the paper has independent support for the alignment mechanism via the AMPT fixed-reaction-plane test and the η-correlation study, which strengthens the qualitative claim even as the quantitative formula remains unvalidated.","tokens_in":8657,"tokens_out":11600,"duration_ms":117680,"concrete_test":"Run a toy Monte Carlo with a known unweighted v2 and a realistic pT(φ) correlation, e.g., draw pT_i from a distribution whose mean is a+b cos2(φ_i−Ψ2) with b>0. Compute S0 using the exact pT-weighted definition (Eq. 1) and record the standard v2 from the same events. Compare the event-averaged S0 versus v2 curve to Eq. (21) using both the input v2 and the pT-weighted v2_eff. If Eq. (21) only matches when v2 is replaced by v2_eff, the paper's stated relation to the usual v2 is unsupported; if the coefficient 2/3 changes by more than 20%, the quantitative claim requires revision.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's quantitative centerpiece, Eq. (21), is derived under the explicit assumption of equal particle weights ('we first consider equal particle weights i.e., pT,i = 1', Sec. III). The actual spherocity (Eq. 1) is pT-weighted: the ratio is Σ_i pT,i |sin(φ_i−ψ)| / Σ_i pT,i, so pT does not cancel. Because v2 depends on pT, the effective distribution probed by the minimization is P_eff(φ) ∝ pT(φ) P(φ), which has a second-harmonic coefficient v2_eff = ⟨pT cos2(φ−Ψ2)⟩/⟨pT⟩, generally different from the standard v2 (a simple model pT(φ)=a+b cos2(φ−Ψ2) gives v2_eff = (b+2a v2)/(2(a+b v2))). Hence Eq. (21) is not the relation for the measured observable; it could misstate the anti-correlation strength by an O(1) factor. The paper neither proves that ψmin = Ψ2 and Eq. (21) remain valid under pT weighting nor compares Eq. (21) against pT-weighted simulations. Since the abstract and conclusions explicitly invoke S0 ≈ (1 − 2/3 v2)^2 as the analytical evidence for the flow-selection interpretation, this is a load-bearing gap.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","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.","tokens_in":8997,"tokens_out":5189,"duration_ms":56082,"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":[{"comment":"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.","section":"§III, Eq. (5)–(21)"},{"comment":"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.","section":"§III, Eq. (12)–(21)"},{"comment":"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.","section":"§II.B and §III.A"}],"minor_comments":[{"comment":"The denominator is typeset as `Σ_i ⃗|pT,i|`; this should be `Σ_i |p_T,i|`.","section":"Eq. (1)"},{"comment":"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.","section":"§III"},{"comment":"The AMPT figure would benefit from error bars or a statement of statistical uncertainty, and the number of events used should be specified.","section":"Fig. 3"},{"comment":"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.","section":"Ref. [26]"}],"recommendation":"major_revision","confidential_remarks":"The equal-pT gap is fixable with additional derivations or numerical tests, and the core physical argument is plausible and timely. I would not reject the paper, but the current quantitative centerpiece overstates certainty. If the authors supply a p_T-weighted derivation or a direct p_T-weighted AMPT/toy validation, the paper could become acceptable for publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear colleague,\n\nHere's the quick take: the paper makes a point that needed making — in AA collisions, transverse spherocity minimization aligns with the event plane, so low S0 is mostly a flow selection effect, not a jet topology tag. The Fourier derivation (Eqs. 11–22) is clean and the toy/AMPT checks support the axis alignment. I buy the qualitative story.\n\nWhat's actually new is the analytic connection S0 ≈ (1 − 2/3 v2)^2 in the continuum equal-weight limit. The prior literature reported the anti-correlation empirically; nobody, as far as I know, derived it from the definition. The derivation is simple but correct, and it explains why the minimizing axis picks out Ψ2.\n\nThe soft spot is exactly what the stress test flags. The derivation starts with \"we first consider equal particle weights, pT,i = 1\" and never returns to the pT-weighted case. Eq. 1, the actual spherocity, is weighted by pT, so pT does not cancel. With v2(pT) nontrivial, the effective distribution probed by the minimization has a modified second harmonic; Eq. 21 could be off by an O(1) factor. The paper doesn't prove that ψmin = Ψ2 and Eq. 21 survive pT weighting, and doesn't compare Eq. 21 to pT-weighted simulations. Since the abstract and conclusions lean on Eq. 21 as the quantitative evidence, that's a load-bearing gap — but it's fixable, not fatal. A derivation with pT weight, even in a simple factorized form, or a plot of Eq. 21 against measured spherocity from AMPT would address it.\n\nMinor points: no code/data for the figures, which matters for the AMPT/PYTHIA panels but isn't disqualifying. The \"global nature\" section with η correlations is suggestive and consistent with the flow-selection picture. The sphericity appendix is fine but under the factorization approximation, which they state.\n\nWho should read this: anyone doing spherocity-based event classification in AA, and anyone citing Refs. [21–24,27] as evidence for jet-like topology in low-S0 events. The qualitative conclusion — that published spherocity-selected AA results need to be re-interpreted with v2 in mind — is likely robust even if the exact Eq. 21 changes.\n\nWorth refereeing? Yes. Send it to review; it's a serious, honest paper with a correct core and a well-defined technical gap. I'd ask for a pT-weighted derivation or explicit numerical check before acceptance. I'd cite it for the qualitative point.\n\nBest.","headline":"Qualitatively right and useful, but the quantitative centerpiece is only proven in an equal-pT limit that isn't the measured observable.","tokens_in":9442,"tokens_out":1671,"would_cite":true,"duration_ms":16765,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["25.75.-q","25.75.Ld"],"model":"deepseek-v4-flash","headline":"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.","keywords":["transverse spherocity","elliptic flow","event shape","anisotropic flow","event plane","heavy-ion collisions","event shape engineering","quark-gluon plasma"],"falsifier":"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.","tokens_in":8560,"feed_emoji":"🎯","tokens_out":7135,"duration_ms":64161,"temperature":0.7,"pith_summary":"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.","feed_headline":"Spherocity selects elliptic flow, not jets, in heavy-ion collisions","feed_subtitle":"Low-spherocity events are high-flow events, so previous jet-based interpretations of heavy-ion data need revisiting.","key_machinery":"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","core_discovery":"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","pith_inferences":["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."],"forward_implications":["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."],"fun_headline_variants":["Spherocity isn't jetty: it's flow in heavy ions","Low spherocity events are high-flow, not jet-rich","Spherocity and elliptic flow are intrinsically linked","Heavy-ion spherocity reads flow, not jets","Spherocity's jet tag is actually flow-driven"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"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.","fun_headline_variants_meta":{"raw":{"variants":["Spherocity isn't jetty: it's flow in heavy ions","Low spherocity events are high-flow, not jet-rich","Spherocity and elliptic flow are intrinsically linked","Heavy-ion spherocity reads flow, not jets","Spherocity's jet tag is actually flow-driven"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000255,"raw_usage":{"total_tokens":1463,"prompt_tokens":853,"completion_tokens":610,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":597,"completion_tokens_details":{"reasoning_tokens":537}},"tokens_in":597,"tokens_out":610,"duration_ms":5640,"temperature":1.0,"reasoning_tokens":537,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-01T08:16:21.547666+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"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.","supporting_citations":[],"review_version":1}