Pith. sign in

REVIEW 3 major objections 5 minor 17 references

Hydrodynamic gradients alone cannot explain the measured Lambda spin polarization in p+Pb collisions.

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

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

Hydrodynamic spin contributions (thermal vorticity negative, thermal shear positive) reproduce Au+Au polarization but cannot match CMS p+Pb data, so new polarization mechanisms may be needed.

T0 review reviewed 2026-08-05 challenge →

load-bearing objection The p+Pb comparison is the only new piece, but the paper's conclusion ignores a hydrodynamic mechanism it cites itself, so the strong 'new mechanisms needed' claim doesn't follow. the 3 major comments →

arxiv 2509.00380 v1 pith:NTZWED53 submitted 2025-08-30 nucl-th

Hydrodynamic effects on spin polarization along the beam direction in Au+Au and p+Pb collisions

Cong Yi , Xiang-Yu Wu , Jie Zhu , Shi Pu , Guang-You Qin This is my paper
classification nucl-th
keywords spin polarizationLambda hyperonbeam-direction polarizationrelativistic hydrodynamicsthermal vorticitythermal-shear tensorAu+Au collisionsp+Pb collisions
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

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

The reading

This paper asks whether standard relativistic hydrodynamics can explain the measured spin alignment of Lambda hyperons along the beam direction, encoded in the Fourier coefficient . Using a (3+1)-dimensional hydrodynamic simulation with several equilibrium prescriptions, the authors find a consistent pattern: the thermal-shear tensor pushes this observable positive, while thermal vorticity pushes it negative, in both Au+Au and p+Pb collisions. In Au+Au at 200 GeV the combined effect can reproduce the measured centrality dependence for two of the three equilibrium scenarios, but in p+Pb at 8.16 TeV the net signal is negative for all three scenarios while measured data are positive. The conclusion is that current hydrodynamic mechanisms alone cannot explain spin polarization in small collision systems, so additional physics is needed.

Core claim

The paper's central claim is a sign competition and a failure. Decomposing the local-equilibrium spin polarization into a thermal-vorticity contribution (from fluid rotation) and a thermal-shear contribution (from fluid velocity and temperature gradients), the authors compute the second Fourier sine coefficient of longitudinal spin polarization, <Pz sin 2(phi_p - Psi_2)>, for Lambda hyperons using the CLVisc hydrodynamic model with AMPT initial conditions for Au+Au and Trento-3D initial conditions for p+Pb, followed by a hadronic afterburner. In both collision systems the shear term contributes positively and the vorticity term contributes negatively to this coefficient. In Au+Au at sqrt(s_N

What carries the argument

The central machinery is the decomposition of the spin-polarization vector into a thermal-vorticity piece (the antisymmetric gradient of velocity over temperature, i.e. local fluid rotation) and a thermal-shear piece (the symmetric gradient of velocity over temperature, i.e. local fluid deformation). Their signed contributions to <Pz sin 2(phi_p - Psi_2)> are computed from freeze-out hypersurface gradients through a modified Cooper-Frye freeze-out integral. Three equilibrium scenarios—full Lambda equilibrium, s-quark equilibrium, and isothermal equilibrium without temperature gradients—change which mass enters the formula and which gradients act, but they do not change the sign competition.

Load-bearing premise

That hadronic rescattering after freeze-out has little effect on the Lambda spin direction; if the afterburner substantially rotates spins, the predicted sign and size of the observable could change.

What would settle it

Recompute <Pz sin 2(phi_p - Psi_2)> with the hadronic afterburner effects fully included or removed; if the p+Pb total changes sign and matches the positive measured values, the claim that hydrodynamic gradients alone cannot explain the data would collapse. Alternatively, measure the same observable in p+Pb at higher multiplicities, where the hydrodynamic description is reliable, and check whether the sign stays positive.

Watch this falsifier. Get emailed when new claim-graph text bears on it.

If this is right

  • In Au+Au, the combined vorticity-plus-shear prediction reproduces the measured centrality dependence when the spin chemical potential is tied to s quarks or when temperature gradients are dropped; the Lambda-equilibrium choice does not work.
  • In p+Pb, the two hydrodynamic contributions have opposite signs, and the cancellation is strong enough that the total remains negative at all multiplicities and in all three equilibrium scenarios.
  • Because the sign mismatch is scenario-independent, choosing a different equilibrium prescription is not enough to cure the p+Pb discrepancy.
  • The result directly motivates adding non-hydrodynamic polarization mechanisms, such as hadronic final-state interactions or initial-state effects, to the standard framework.

Where Pith is reading between the lines

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

  • If the no-rescattering assumption is relaxed, hadronic final-state scattering could rotate Lambda spins enough to change the sign; a direct switch-on/switch-off test of the afterburner would settle this quickly.
  • The paper sets aside electromagnetic fields and baryon chemical potential; in a small, fast-expanding system those gradients might contribute positively to <Pz sin 2(phi_p - Psi_2)> and are a natural next test.
  • The sign competition itself could be used as a diagnostic: since shear and vorticity enter with opposite signs, observables that isolate one contribution, e.g. by particle species or by rapidity, would help distinguish the proposed new mechanisms.
  • A straightforward extension is to run the same calculation for other small systems, such as O+O collisions, to see whether the negative-total failure is generic to small systems or specific to p+Pb.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit.

Referee Report

3 major / 5 minor

Summary. The paper uses the (3+1)-D CLVisc hydrodynamic model with AMPT initial conditions for Au+Au at 200 GeV and Trento-3D for p+Pb at 8.16 TeV, followed by SMASH hadronic transport, to compute the second Fourier sine coefficient of the longitudinal spin polarization of Λ hyperons, <Pz sin(2(phi_p - Psi_2))>, under three local-equilibrium scenarios (Λ equilibrium, s-quark equilibrium, and iso-thermal equilibrium). The authors separate the contributions of thermal vorticity and the thermal-shear tensor. They report that in Au+Au the s-quark and iso-thermal scenarios reproduce the STAR centrality data, while the Λ-equilibrium scenario does not. In p+Pb, the thermal-shear contribution is positive and the vorticity contribution negative, but the total is negative across all scenarios, apparently contradicting the positive CMS data. The paper concludes that current hydrodynamic mechanisms cannot explain the p+Pb spin polarization and that new mechanisms are needed.

Significance. If the result were established, it would be a striking and important finding: it would indicate that the standard hydrodynamic spin-polarization framework, which successfully describes Au+Au, fails in small systems and that new mechanisms are required. The paper's strengths are its systematic comparison of three equilibrium scenarios, the explicit separation of vorticity and shear contributions, and the use of a consistent hydrodynamic + transport framework with initial conditions tuned to bulk observables. The authors are also transparent in noting that their low-multiplicity results are unreliable. However, because the conclusions are sign-level and the calculation omits at least one hydrodynamic mechanism identified in the paper's own introduction, the significance as it stands is limited.

major comments (3)
  1. [§2 (first paragraph) and §4 Summary] The decomposition S^mu = S^mu_thermal + S^mu_th-shear is introduced 'in absence of electromagnetic field and baryon chemical potential'. Section 1 explicitly lists 'gradient of chemical potential' as a hydrodynamic contribution [5-7]. The p+Pb measurement is at |eta| < 2.4, where net-baryon gradients are not guaranteed negligible, and the observable <Pz sin(2phi_p - 2Psi_2)> probes azimuthal modulations to which a mu_B/T gradient can contribute. The calculation therefore omits one of the hydrodynamic mechanisms cited in the paper's own introduction. The Summary's conclusion that 'current hydrodynamic effects alone cannot explain the spin polarization in p+Pb collisions' overreaches: at most the paper rules out the thermal-vorticity + thermal-shear combination. Please either include the mu_B-gradient contribution (or an estimate of its size) or restrict the conclusion accordingly.
  2. [§2 (numerical setup)] The calculation assumes 'hadronic scattering has little effect on spin polarization' and uses the modified Cooper-Frye formula at freeze-out. Since SMASH is used for the bulk, the same framework could in principle test whether rescattering rotates Lambda spins before decay. This assumption is load-bearing because the claimed sign mismatch in Fig. 2 is a sign-level statement; even a moderate change from rescattering could move the predicted total from negative to positive. No test of this assumption is presented. Please either quantify/justify the assumption with a comparison or treat it explicitly as a caveat limiting the central claim.
  3. [§3 (final paragraph) and Fig. 2] The authors state that 'our polarization results at lowest multiplicities are unreliable due to the limitations of hydrodynamic applicability,' yet Fig. 2 compares the full multiplicity range and the summary claims a generic contradiction. No numerical or chi-square comparison is shown for the high-multiplicity regime where the model is claimed to be reliable, and the detailed high-multiplicity discussion is deferred to Ref. [14]. Since the sign of the discrepancy is the central claim, the reader cannot assess whether it survives in the reliable region. Please show the comparison for the reliable high-multiplicity range and quantify the disagreement.
minor comments (5)
  1. [Figure 1] The label 'CLVsic' should be 'CLVisc'.
  2. [Figure 2 caption/panels] The caption order of the three scenarios appears inconsistent with the panel labels: panel (b) is labeled 'iso_th' but the caption calls it s-quark equilibrium, and panel (c) is labeled 's_th' but the caption calls it iso-thermal. Please unify the caption and panel labels.
  3. [§2 (definitions)] The thermal-shear tensor is defined without an explicit factor 1/2 and without trace subtractions; clarify the convention so that the sign and normalization of the shear contribution are unambiguous.
  4. [General] The text alternates 'iso-thermal' and 'isothermal'; use one spelling consistently.
  5. [§2] All formulas in Section 2 are inline; adding numbered equations would make the results easier to check and cite.

Circularity Check

0 steps flagged

No significant circularity: the spin-polarization signal is a genuine prediction from a bulk-tuned hydrodynamic background.

full rationale

The derivation chain is self-contained in the relevant sense. Hydrodynamic parameters are tuned to reproduce bulk multiplicity and v2, while the beam-direction spin polarization <Pz sin 2(phi_p - Psi_2)> is then computed from temperature/flow gradients via the modified Cooper-Frye formula. No polarization datum is used in tuning, so the p+Pb sign disagreement is a genuine prediction rather than a fitted input renamed as a prediction. The spin-polarization decomposition S^mu = S^mu_thermal + S^mu_th-shear is imported from Refs. [8,9,10,14], including the authors' own previous work [10,14]; this is a provenance self-citation for a standard Kubo/kinetic-theory formula, but it is parameter-free and does not assume the Au+Au or p+Pb sign of the observable, so it is not load-bearing circularity. The paper explicitly self-limits its result ('our polarization results at lowest multiplicities are unreliable', 'Assuming that the hadronic scattering has little effect on spin polarization') and separately notes the chemical-potential gradient contribution in the Introduction, which it omits; those are completeness/assumption caveats about the physical input set, not circular reductions of the output to the input. The central claim is therefore an independent numerical result, and no step in the chain is equivalent by construction to its inputs.

Axiom & Free-Parameter Ledger

2 free parameters · 5 axioms · 0 invented entities

The model imports spin polarization formulas from previous literature, assumes local equilibrium and negligible afterburner effects, and uses hydrodynamic parameters tuned to bulk observables. No new particles, fields, or conserved quantities are introduced.

free parameters (2)
  • hydrodynamic model parameters (shear viscosity, initial condition normalization, freeze-out details) = not stated in the paper
    Section 2 says the simulation parameters can successfully reproduce multiplicity and v2; these are calibrated to bulk observables and inherited from Refs [10,14,15,16].
  • strange quark mass m_s = 0.3 GeV
    Used in the Cooper-Frye formula for the s-quark equilibrium scenario; a phenomenological constituent mass, not derived in the paper or fitted to spin data.
axioms (5)
  • domain assumption The spin polarization vector can be decomposed into thermal vorticity and thermal-shear contributions in the absence of electromagnetic field and baryon chemical potential.
    Section 2 invokes this decomposition from Refs [8,9,10] without re-deriving it; the entire calculation depends on it.
  • domain assumption Local thermal equilibrium at freeze-out and the modified Cooper-Frye formula apply to spin polarization.
    Section 2 uses the modified Cooper-Frye formula to convert hydrodynamic gradients on the freeze-out hypersurface into Lambda spin polarization.
  • ad hoc to paper Hadronic scattering has little effect on spin polarization.
    Section 2 states this assumption explicitly; if afterburner rescattering changes Lambda spin, the predicted observable could shift.
  • domain assumption Hydrodynamics is reliable only at high multiplicity.
    Section 3 says polarization results at lowest multiplicities are unreliable due to limitations of hydrodynamic applicability, so the comparison to CMS data rests on the high-multiplicity regime.
  • domain assumption AMPT and TRENTo-3D initial conditions provide realistic bulk profiles for the two collision systems.
    Section 2 chooses these initial conditions for Au+Au and p+Pb; the spin prediction inherits any uncertainty in the initial state.

reviewed 2026-08-05 · how reviews work

0 comments
Cite this review

Pith. "Pith review of Hydrodynamic effects on spin polarization along the beam direction in Au+Au and p+Pb collisions." pith.science (2026). https://pith.science/paper/NTZWED53

@misc{pith2026250900380,
  author       = {Pith},
  title        = {Pith review of: Hydrodynamic effects on spin polarization along the beam direction in Au+Au and p+Pb collisions},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/NTZWED53}},
  note         = {Machine review of arXiv:2509.00380}
}
Share X Bluesky LinkedIn Reddit HN
read the original abstract

We investigate hydrodynamic effects on the spin polarization of $\Lambda$ hyperons in Au+Au collisions at $\sqrt{s_{NN}} = 200$ GeV and p+Pb collisions at $\sqrt{s_{NN}} = 8.16$ TeV using the CLVisc hydrodynamic framework. We present numerical results for the second Fourier sine coefficient of the longitudinal spin polarization, $\langle P_{z} \sin 2(\phi_{p} - \Psi_{2}) \rangle$, as a function of multiplicity (centrality) under three equilibrium scenarios: $\Lambda$ equilibrium, $s$-quark equilibrium, and isothermal equilibrium. We highlight the respective roles of thermal vorticity and the thermal-shear tensor in generating $\langle P_{z} \sin 2(\phi_{p} - \Psi_{2}) \rangle$ across collision systems and scenarios.

Figures

Figures reproduced from arXiv: 2509.00380 by Cong Yi, Guang-You Qin, Jie Zhu, Shi Pu, Xiang-Yu Wu.

Figure 1
Figure 1. Figure 1: The centrality dependence of ⟨Pz sin(2ϕp − 2Ψ2)⟩ for Λ hyperons in √ sNN = 200 GeV Au+Au collisions. Results are shown for three scenarios: Λ equilibrium (blue dashed line), s quark equilibrium (orange dash-dotted line), and iso-thermal equilibrium (green solid line). Red markers represent experimental data from Ref. [3]. polarization induced by the thermal vorticity ϖαβ = 1 2 h ∂α  uβ T  − ∂β  uα T i … view at source ↗
Figure 2
Figure 2. Figure 2: The multiplicity ⟨Nch⟩ dependence of ⟨Pz sin(2ϕp − 2Ψ2)⟩ for Λ hyperons at √ sNN = 8.16 TeV p+Pb collisions. Results are shown for three scenarios: Λ equilibrium (a) , s quark equilibrium (b), and iso-thermal equilibrium (c). The contributions from thermal vorticity (blue triangles), thermal￾shear tensor (green diamonds), and their combined effects (black circles) are presented. Red points correspond to CM… view at source ↗

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

17 extracted references · 2 canonical work pages

  1. [1]

    Z. T. Liang and X. N. Wang, Phys. Rev. Lett. 94, 102301 (2005) [erratum: Phys. Rev. Lett. 96, 039901 (2006)] doi:10.1103 /PhysRevLett.94.102301 [arXiv:nucl-th /0410079 [nucl-th]]

    work page 2005

  2. [2]

    Adamczyk et al

    L. Adamczyk et al. [STAR], Nature 548, 62-65 (2017) doi:10.1038 /nature23004 [arXiv:1701.06657 [nucl-ex]]

    Pith/arXiv arXiv 2017

  3. [3]

    Adam et al

    J. Adam et al. [STAR], Phys. Rev. Lett. 123, no.13, 132301 (2019) doi:10.1103/PhysRevLett.123.132301 [arXiv:1905.11917 [nucl-ex]]

    arXiv 2019

  4. [4]

    Becattini and I

    F. Becattini and I. Karpenko, Phys. Rev. Lett. 120, no.1, 012302 (2018) doi:10.1103/PhysRevLett.120.012302 [arXiv:1707.07984 [nucl-th]]

    Pith/arXiv arXiv 2018

  5. [5]

    Hidaka, S

    Y . Hidaka, S. Pu and D. L. Yang, Phys. Rev. D 97, no.1, 016004 (2018) doi:10.1103/PhysRevD.97.016004 [arXiv:1710.00278 [hep-th]]

    Pith/arXiv arXiv 2018

  6. [6]

    Becattini, M

    F. Becattini, M. Buzzegoli and A. Palermo, Phys. Lett. B 820, 136519 (2021) doi:10.1016/j.physletb.2021.136519 [arXiv:2103.10917 [nucl-th]]

    work page arXiv 2021

  7. [7]

    S. Y . F. Liu and Y . Yin, JHEP 07, 188 (2021) doi:10.1007 /JHEP07(2021)188 [arXiv:2103.09200 [hep-ph]]

    Pith/arXiv arXiv 2021

  8. [8]

    B. Fu, S. Y . F. Liu, L. Pang, H. Song and Y . Yin, Phys. Rev. Lett. 127, no.14, 142301 (2021) doi:10.1103/PhysRevLett.127.142301 [arXiv:2103.10403 [hep-ph]]

    Pith/arXiv arXiv 2021

  9. [9]

    Becattini, M

    F. Becattini, M. Buzzegoli, G. Inghirami, I. Karpenko and A. Palermo, Phys. Rev. Lett. 127, no.27, 272302 (2021) doi:10.1103/PhysRevLett.127.272302 [arXiv:2103.14621 [nucl-th]]

    Pith/arXiv arXiv 2021

  10. [10]

    C. Yi, S. Pu and D. L. Yang, Phys. Rev. C 104, no.6, 064901 (2021) doi:10.1103/PhysRevC.104.064901 [arXiv:2106.00238 [hep-ph]]

    Pith/arXiv arXiv 2021

  11. [11]

    Alzhrani, S

    S. Alzhrani, S. Ryu and C. Shen, Phys. Rev. C 106, no.1, 014905 (2022) doi:10.1103/PhysRevC.106.014905 [arXiv:2203.15718 [nucl-th]]

    Pith/arXiv arXiv 2022

  12. [12]

    X. Y . Wu, C. Yi, G. Y . Qin and S. Pu, Phys. Rev. C 105, no.6, 064909 (2022) doi:10.1103/PhysRevC.105.064909 [arXiv:2204.02218 [hep-ph]]

    Pith/arXiv arXiv 2022

  13. [13]

    Hayrapetyan et al

    A. Hayrapetyan et al. [CMS], [arXiv:2502.07898 [nucl-ex]]

    arXiv

  14. [14]

    C. Yi, X. Y . Wu, J. Zhu, S. Pu and G. Y . Qin, Phys. Rev. C 111, no.4, 044901 (2025) doi:10.1103/PhysRevC.111.044901 [arXiv:2408.04296 [hep-ph]]

    Pith/arXiv arXiv 2025

  15. [15]

    L. Pang, Q. Wang and X. N. Wang, Phys. Rev. C 86, 024911 (2012) doi:10.1103/PhysRevC.86.024911 [arXiv:1205.5019 [nucl-th]]

    Pith/arXiv arXiv 2012

  16. [16]

    X. Y . Wu, G. Y . Qin, L. G. Pang and X. N. Wang, Phys. Rev. C105, no.3, 034909 (2022) doi:10.1103/PhysRevC.105.034909 [arXiv:2107.04949 [hep-ph]]

    Pith/arXiv arXiv 2022

  17. [17]

    Weil et al

    J. Weil et al. [SMASH], Phys. Rev. C 94, no.5, 054905 (2016) doi:10.1103/PhysRevC.94.054905 [arXiv:1606.06642 [nucl-th]]

    Pith/arXiv arXiv 2016

This paper was first reviewed by deepseek-v4-flash on August 5, 2026.