REVIEW 3 major objections 4 minor 32 references
Relaxation time for quark spin and thermal vorticity alignment in heavy-ion collisions
T0 review · 3 major / 4 minor · reviewed 2026-08-14 · deepseek-v4-flash
Pith's one-line read The relaxation time for quark spin to align with thermal vorticity is about 10 fm or more for small plasma angular velocities and below 3 fm for large ones, with antiquarks aligning more slowly at finite chemical potential.
desk verdict Transparent but model-dominated estimates; the credible parts are the explicit HTL machinery and honest caveats, not the tau numbers. 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 object that carries the argument is the effective vertex $\lambda^\mu_a = g\,(\sigma^{\alpha\beta}/2)\,\omega_{\alpha\beta}\,\gamma^\mu t_a$ (Eq. (4)), which inserts the quark spin operator $\sigma^{\alpha\beta}/2$ into the elementary quark–gluon vertex and uses the thermal vorticity $\omega_{\alpha\beta}$ as the field strength. Inserting this vertex into the one-loop quark self-energy with hard-thermal-loop (HTL) gluon propagators, the interaction rate of Eq. (13) follows with an overall factor $(\omega/T)^2$, and the relaxation time is its inverse after phase-space integration (Eqs. (15)–(17)). The relevant scatterings are off space-like thermal gluons, where the HTL spectral densities have Landau-damping support, and the $\omega^2$ prefactor is what makes the alignment time so sensitive to the assumed angular velocity.
What would settle it
Measure the global polarization of $\Lambda$ and $\bar{\Lambda}$ in peripheral Au+Au collisions near $\sqrt{s_{NN}}\approx 10$ GeV, where the vorticity is expected to be small, and infer the effective spin–vorticity relaxation time from the observed polarization magnitude; if the inferred time is much shorter than the approximately 10 fm quoted here for small $\omega$, the assumed vertex underestimates the coupling. A simpler internal check is to recompute $\tau$ from a first-principles or lattice-derived spin–vorticity coupling and see whether $\tau$ at $\omega\approx 0.12~\mathrm{fm}^{-1}$ stays above 10 fm.
Extended reading notes
Core claim
The paper's central claim is that spin–vorticity equilibration in the quark–gluon plasma is a race between the rate $\Gamma$ and the fireball lifetime, and the rate is controlled by the ratio $\omega/T$. Concretely, the computed relaxation time $\tau \equiv 1/\Gamma$, obtained from Eqs. (13)–(17) with the effective vertex of Eq. (4) and hard-thermal-loop gluon propagators, is $\sim 10$ fm or larger for the small vorticity values inferred from event-generator calculations, while it drops below $\sim 3$ fm for the larger vorticity values of a scenario in which the plasma retains the initial angular momentum. Because the interaction rate is proportional to the quark/antiquark occupation number, replacing $\mu$ by $-\mu$ makes the antiquark relaxation time larger, so at finite chemical potential quark spins align faster than antiquark spins. These are model outputs of the assumed coupling, not measured constraints.
Load-bearing premise
The calculation assumes that the force aligning quark spin with the plasma's rotation is exactly the one written in Eq. (4); if the true force is weaker or has a different shape, all the relaxation times change.
Editorial extensions
If this is right
- If the small-vorticity scenario is realized in peripheral collisions, spin–vorticity alignment cannot be assumed to be complete, and observed hyperon polarization must be explained by other mechanisms or by a different vorticity distribution.
- If the large-vorticity scenario is realized, the relaxation time is short enough that thermal-vorticity-based models of global polarization are on solid ground.
- Because $\tau$ decreases as temperature and chemical potential increase, alignment is fastest in the hottest, densest part of the fireball and slowest at the edges.
- The slower antiquark relaxation at finite $\mu$ implies a measurable $\Lambda/\bar{\Lambda}$ polarization asymmetry if hadronization preserves the constituent quark polarization.
- The $(T/\omega)^2$ scaling means that extracting vorticity from polarization data requires folding in the relaxation efficiency, not assuming equilibrium alignment.
Reading between the lines
- An immediate extension the paper does not pursue is to feed this $\tau$ into spin-transport equations: if alignment is incomplete at small $\omega$, spin hydrodynamics needs a finite relaxation term rather than an equilibrium spin-vorticity relation.
- Because Eq. (4) is essentially the only free choice, the qualitative dichotomy (slow at low $\omega$, fast at high $\omega$) is likely robust to other Lorentz-invariant couplings sharing the $\omega^2$ factor, though the numerical boundary would shift.
- The massless-quark approximation likely changes the strange-quark relaxation time; including the strange mass could either shorten or lengthen the estimate, and since $\Lambda$ polarization is the main observable, this is the most relevant refinement.
- The model could be tested by comparing the predicted centrality dependence of the relaxation time (via the impact parameter dependence of $\omega$ and volume) against measured global polarization as a function of centrality.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript addresses whether quark spin can align with thermal vorticity within the lifetime of the quark-gluon plasma created in heavy-ion collisions. The authors model the spin-vorticity interaction by introducing a phenomenological modification of the quark-gluon vertex, Eq. (4), proportional to the spin operator contracted with the thermal vorticity. Using the imaginary-time formalism with a hard-thermal-loop gluon propagator and a bare quark propagator, they compute a one-loop self-energy and from it an interaction rate, Eqs. (13)-(15). The relaxation time is defined as the inverse of the phase-space-integrated total rate, Eq. (17). For a 'small' angular velocity scenario (ω ~ 0.10-0.12 fm^-1) they find τ ~ 3-10 fm or larger, implying that spin-vorticity alignment is not fully established, while for a 'large' angular velocity scenario (ω ~ 0.23-0.46 fm^-1, obtained from an explicitly artificial rigid-rotation assumption) they find τ ≲ 3 fm, implying efficient alignment. Antiquarks are found to relax more slowly at finite chemical potential, which is suggested as a possible source of hadron-antihadron polarization differences.
Significance. If the model and the identification of τ as a relaxation time were validated, this would provide a useful first estimate of the timescale for spin-vorticity equilibration and could inform interpretations of Λ and anti-Λ polarization measurements. The paper makes a clear, transparent computation using standard HTL thermal field theory techniques, and it openly acknowledges the phenomenological nature of its central input vertex. However, the quantitative results are entirely controlled by an uncalibrated effective vertex and by an interpretation of τ that is not derived from a kinetic or master-equation analysis. The explicit caveat that 'other modelings are possible' means that the quoted fm values are model outputs rather than firm QCD predictions. The paper's main qualitative message—that relaxation is slower at smaller ω and for antiquarks at finite μ—survives as a plausible tendency, but the specific thresholds (10 fm vs 3 fm) should be presented with appropriate uncertainty.
major comments (3)
- [§2, Eq. (15) and Eq. (17)] The definition τ ≡ 1/Γ, with Γ given by the phase-space-integrated total rate in Eq. (15), is not justified as the spin relaxation time. Equation (13) defines a momentum-dependent rate Γ(p0) for a quark of energy p0. Integrating this over all momenta and multiplying by the volume V yields a total number of interactions per unit time in the entire system, not the relaxation rate for an individual quark's spin or a thermally averaged single-particle rate. The time evolution of spin-vorticity alignment should be derived from a Boltzmann or master equation, producing a relaxation time as an inverse of an appropriately averaged transition rate. As it stands, the numerical values of τ in Figs. 2-6 are not demonstrably related to the physical equilibration time of spin alignment. The authors need to provide a kinetic-theory justification or at least compare their definition with the standard thermal-width approach.
- [§2, Eq. (4)] The effective vertex λ_a^μ = g (σ^{αβ}/2) ω_{αβ} γ^μ t_a is introduced purely phenomenologically, and the paper itself concedes that 'other modelings are possible.' Since Γ scales as the square of the vertex (explicitly as α_s (ω/T)^2 in Eq. (13)), any change in the overall normalization or Lorentz structure of this vertex changes every relaxation time in Figs. 2-6, not only by an overall factor but also by altering the kernel functions C_L and C_T in Eq. (14). There is no lattice calculation, no matching to a known effective theory, and no data constraint fixing the coupling strength g or the tensor structure. The authors should at least perform a sensitivity study, e.g., varying the overall coupling normalization and comparing different Lorentz structures (e.g., axial-vector versus vector), or provide a physical derivation from a more fundamental effective action. Without this, the quoted 10-fm versus 3-fm distinction is a property of the assumed vertex, not a prediction of QCD.
- [§3, Fig. 3 and accompanying text] The 'large angular velocity' scenario leading to the conclusion of efficient alignment is explicitly described in the text as 'an artificial way to describe the collision' because it assumes that all initial angular momentum is converted into rigid rotation, whereas Ref. [25] finds that most angular momentum manifests as local fluid shear. The abstract's statement that 'when the angular velocity created in the reaction is large, the alignment is efficient and well within the lifetime of the system' is therefore based on a limiting, not realistic, case. The conclusions should be reframed to distinguish a bounded, data-informed scenario (ω ~ 0.1 fm^-1) from an extreme upper-limit scenario, and the abstract should not present the efficient-alignment case as the typical outcome without this qualification.
minor comments (4)
- [§2, after Eq. (13)] The text states that 'for consistency of the approximation where we have considered massless quarks, we have also dropped terms proportional to the quark four-momentum components.' This step is not shown; please clarify which terms are dropped and whether that approximation is uniformly justified in the integrand of Eq. (13).
- [§3, Eq. (18)] The definition of ω in Eq. (18) uses the velocity along the beam axis, but the thermal vorticity in Eq. (1) is a four-dimensional tensor; the relation between this non-relativistic estimate and ω_{μν} used in Eq. (13) is not spelled out and should be stated explicitly.
- [Figures 2-5, captions] The captions of Figs. 2 and 4 say 'Notice that τ is of order ≲ 3 fm only for the largest T and μ considered' and similar, but the figures show a range of values and this phrasing is awkward; please clarify whether the conclusion applies to the plotted temperature interval or to the endpoints.
- [References] The comparison with Ref. [23] (Kapusta, Rrapaj, Rudaz) is only mentioned in the introduction; the present results should be quantitatively compared with that earlier estimate to help the reader judge the difference arising from the different interaction mechanism.
Circularity Check
No circularity: the spin-vorticity vertex is an explicit model input, and the relaxation times follow from an independent phase-space/HTL computation.
full rationale
The paper's central output, tau = 1/Gamma, is not derived from the quantity it is claimed to predict, nor is any fitted parameter renamed as a prediction. The interaction is introduced explicitly as a phenomenological ansatz in Eq. (4), lambda^mu_a = g (sigma^{alpha beta}/2) omega_{alpha beta} gamma^mu t_a, and the paper states plainly that this is a modeling choice and that other modelings are possible. The relaxation rate Gamma is then computed from the standard one-loop quark self-energy with HTL gluon propagators and spectral densities, Eqs. (3), (5)-(15), using only standard inputs: alpha_s = 0.3, T, mu, and omega values taken from external HIJING/AMPT and UrQMD simulations. No parameter is fitted to the resulting relaxation times, and the numerical threshold between the small-omega and large-omega scenarios is a genuine output of the phase-space integrals, not a restatement of the input vertex. The paper's own caveat that the description requires modeling is a limitation on external validity, not evidence of circular argumentation: a model-dependent prediction is still a prediction. The two self-citations, Refs. [30,31], appear only in the closing discussion about hadronization memory and are not load-bearing for the relaxation-time calculation. There is no self-defined quantity, no fitted-input-as-prediction step, no uniqueness claim imported from the authors' prior work, and no known result being renamed. The calculation is self-contained under its stated assumptions, so the circularity score is 0.
Assumptions & free parameters
free parameters (4)
- strong coupling alpha_s =
0.3
- angular velocity omega (small-vorticity scenario) =
0.12 fm^-1 at 10 GeV, 0.10 fm^-1 at 200 GeV
- angular velocity omega (large-vorticity scenario) =
0.23 fm^-1 at 10 GeV, 0.46 fm^-1 at 200 GeV
- spin-vorticity coupling strength in the effective vertex =
Implicitly 1 in Eq. (4)
assumptions (5)
- ad hoc to paper The effective vertex Eq. (4) correctly captures spin-vorticity alignment in the QGP.
- domain assumption The HTL gluon propagator (Eqs. 5-8) is adequate for the rate calculation.
- domain assumption Quarks can be treated as massless.
- domain assumption The vorticity values from Hijing/AMPT and UrQMD at Delta t = 0.4 fm represent the relevant QGP angular velocity.
- ad hoc to paper tau = 1/Gamma with Gamma defined by Eq. (15) is the spin relaxation time.
invented entities (1)
-
Phenomenological spin-vorticity quark-gluon vertex Eq. (4)
Cite this review
Pith. "Pith review of Relaxation time for quark spin and thermal vorticity alignment in heavy-ion collisions." pith.science (2026). https://pith.science/paper/7LTG7HPS
@misc{pith2026190900274,
author = {Pith},
title = {Pith review of: Relaxation time for quark spin and thermal vorticity alignment in heavy-ion collisions},
year = {2026},
howpublished = {\url{https://pith.science/paper/7LTG7HPS}},
note = {Machine review of arXiv:1909.00274}
}
read the original abstract
We compute the relaxation time for quark/antiquark spin and thermal vorticity alignment in a quark-gluon plasma at finite temperature and quark chemical potential. We model the interaction of quark/antiquark spin with thermal vorticity as driven by a phenomenological modification of the elementary quark interaction with gluons. We find that in a scenario where the angular velocity of the quark-gluon plasma produced in a peripheral heavy-ion collision is small, quarks/antiquarks take a long time to align their spin with the vorticity. However, when the angular velocity created in the reaction is large, the alignment is efficient and well within the lifetime of the system created in the reaction. The relaxation time is larger for antiquarks which points out to a difference for the polarization of hadrons and antihadrons when this alignment is preserved during hadronization.
Figures
Figures from the paper (1 more)
Reference graph
Works this paper leans on
-
[25]
Le Bellac, Thermal Field Theory , (Cambridge Uni- versity Press, Cambridge, 1996)
M. Le Bellac, Thermal Field Theory , (Cambridge Uni- versity Press, Cambridge, 1996)
work page 1996
-
[1]
The kinematical restrictions that Eq
= 2πϵ(p′ 0)δ(p′ 0 2−E2 p), (11) where Ep =|⃗ p−⃗k|. The kinematical restrictions that Eq. (11) imposes on Eq. (10), limit the integration over gluon energies to the space-like region, namely, |x| < 1, therefore, the part of the gluon spectral densities that contribute to the interaction rate are given by 3 ρL(k0,k ) = x 1−x2 2πm2θ(k2−k2 0) [k2 + 2m2 (1− (...
-
[2]
F. Becattini, G. Inghirami, V. Rolando, A. Beraudo, L. Del Zanna, A. De Pace, M. Nardi, G. Pagliara and V. Chandra, Eur. Phys. J. C 75, 406 (2015)
work page 2015
- [3]
-
[4]
D. E. Kharzeev and D. T. Son, Phys. Rev. Lett. 106, 062301 (2011)
2011
-
[5]
F. Becattini, I. Karpenko, M. A. Lisa, I. Upsal and S. A. Voloshin, Phys. Rev. C 95, 054902 (2017)
work page 2017
-
[6]
Y. Xie, R. C. Glastad and L. P. Csernai, Phys. Rev. C 92, 064901 (2015)
work page 2015
-
[7]
A. Sorin and O. Teryaev, Nucl. Part. Phys. Proc. 273- 275, 2587-2589 (2016)
work page 2016
Show all 32 references
-
[8]
Y. L. Xie, M. Bleicher, H. St¨ ocker, D. J. Wang, and L. P. Csernai, Phys. Rev. C 94, 054907 (2016)
2016
-
[9]
Sorin and O
A. Sorin and O. Teryaev, Phys. Rev. C95, 011902 (2017)
2017
-
[10]
B. I. Abelev, et al. (STAR Collaboration), Phys. Rev. C 76, 024915 (2007); Erratum, Phys. Rev. C 95, 039906 (2017)
2007
-
[11]
Adamczyk et al
L. Adamczyk et al. (STAR Collaboration), Nature 548, 62-65 (2017)
2017
-
[12]
Li, L.-G
H. Li, L.-G. Pang, Q. Wang and X.-L. Xia, Phys. Rev. C 96, 054908 (2017)
2017
-
[13]
Sun and C
Y. Sun and C. M. Ko, Phys. Rev. C 96, 024906 (2017)
2017
-
[14]
Han and J
Z.-Z. Han and J. Xu, Phys. Lett. B 786, 255-259 (2018)
2018
-
[15]
X.-L. Xia, H. Li, Z.-B. Tang and Q. Wang, Phys. Rev. C 98, 024905 (2018)
2018
-
[16]
Baznat, K
M. Baznat, K. Gudima, A. Sorin and O. Teryaev, Phys. Rev. C 97, 041902 (2018)
2018
-
[17]
Karpenko and F
I. Karpenko and F. Becattini, Nucl. Phys. A 982, 519- 522 (2019). 6
2019
-
[18]
Suvarieva, K
D. Suvarieva, K. Gudima and A. Zinchenko, Phys. Part. Nucl. Lett. 15, 182-188 (2018)
2018
-
[19]
E. E. Kolomeitsev, V. D. Toneev and V. Voronyuk, Phys. Rev. C 97, 064902 (2018)
2018
-
[20]
Y. Xie, D. Wang and L. P. Csernai, arXiv:1907.00773 [hep-ph]
1907 arXiv
-
[21]
Y. Guo, S. Shi, S. Feng and J. Liao, arXiv:1905.12613 [nucl-th]
1905 arXiv
- [22]
-
[23]
Adam et al
J. Adam et al. (STAR Collaboration), Phys. Rev. C 98, 014910 (2018)
2018
-
[24]
J. I. Kapusta, E. Rrapaj, S. Rudaz, arXiv:1907.10750 [nucl-th]
1907 arXiv
-
[26]
Deng and X.-G
W.-T. Deng and X.-G. Huang, Phys. Rev. C 93, 064907 (2016)
2016
-
[27]
Jiang, Z.-W
Y. Jiang, Z.-W. Lin and J. Liao, Phys. Rev. C 94, 044910 (2016); Erratum: Phys. Rev. C 95, 049904 (2017)
2016
-
[28]
S. A. Bass, M. Belkacem, M. Bleicher, M. Brandstet- ter, L. Bravina, C. Ernst, L. Gerland, M. Hofmann, S. Hofmann, J. Konopka, G. Mao, L. Neise, S. Soff, C. Spieles, H. Weber, L. A. Winckelmann, H. Stocker, W. Greiner, C. Hartnack, J. Aichelin and N. Amelin Prog. Part. Nucl. Ph...
1998
-
[29]
Kekelidze, A
V. Kekelidze, A. Kovalenko, R. Lednicky, V. Matveev, I. Meshkov, A. Sorin and G. Trubnikov, Nucl. Phys. A 967, 884-887 (2017)
2017
- [30]
-
[31]
Ayala, E
A. Ayala, E. Cuautle , G. Herrera, L. M. Montano, Phys. Rev. C 65, 024902 (2002)
2002
-
[32]
Ayala, E
A. Ayala, E. Cuautle, G. Herrera Corral, J. Magnin and L. M. Monta˜ no, Phys. Lett. B682, 408-412 (2010)
2010
Reviewed August 14, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.