REVIEW 4 major objections 6 minor 49 references
The measured J/ψ spin alignment at forward rapidity is reproduced by a pT-dependent mixture of coalescence and primordial production, with low-pT charmonia inheriting vorticity-induced charm-quark polarization.
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 →
T0 review · deepseek-v4-flash
2026-08-02 07:04 UTC pith:RYJLQMQT
load-bearing objection The forward-rapidity curve is interpolated between two fitted endpoints; the mid-rapidity prediction is the real content and deserves a test. the 4 major comments →
Investigating J/psi spin alignment in heavy-ion collisions within a two-component transport model
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Using the relativistic spin Boltzmann equation for vector mesons, the authors derive a non-relativistic expression for ρ00 of heavy quarkonia and then construct the observed ρ00 as the pT-weighted sum ρ_obs^00(pT) = r_coal(pT) ρ_coal^00(pT) + r_init(pT) ρ_init^00. In the coalescence channel, the deviation of ρ00 from 1/3 is driven by the squared thermal-vorticity component ⟨ω_y^2⟩; in the primordial channel, ρ_init^00 is fixed to 0.36 by the highest-pT ALICE point. With the two fractions taken from a transport model, the forward-rapidity data are reproduced, and the calculation yields a distinctive mid-rapidity prediction: the Lorentz transformation between the J/ψ rest frame and the lab fra
What carries the argument
The central object is the two-component superposition formula (Eq. 32), ρ_obs^00(pT) = r_coal(pT) ρ_coal^00(pT) + r_init(pT) ρ_init^00. The coalescence piece ρ_coal^00 is computed from the non-relativistic limit of the vector-meson spin density matrix derived from the spin Boltzmann equation (Eq. 29), where the only production-relevant terms are the squared thermal-vorticity component ⟨ω_y^2⟩ and the vector-field fluctuation parameters f_T^2, f_z^2; the functions F and G encode the Lorentz boost that produces the rapidity dependence. The primordial baseline ρ_init^00=0.36 represents the effective spin state of high-pT charmonium and is not predicted from first principles.
Load-bearing premise
The result stands or falls on the two-component decomposition: the measured inclusive J/ψ sample is assumed to be dominated by exactly two prompt channels—coalescence and primordial production—with their pT-dependent fractions given by a transport model and the primordial baseline fixed empirically to 0.36 by the highest-pT data point.
What would settle it
Measure the inclusive J/ψ ρ00(pT) at mid-rapidity (|Y| < 0.9) in 30–50% Pb-Pb at 5.02 TeV. The model predicts that the vorticity contribution is kinematically suppressed there, so ρ00−1/3 should be near zero at low pT in the vorticity-only case; a clearly positive low-pT deviation larger than the vector-field sensitivity would falsify the kinematic suppression and, with it, the two-component explanation of the forward-rapidity data. A forward-rapidity high-pT point that does not approach ρ00 = 0.36 would likewise falsify the assumed primordial baseline.
If this is right
- If the two-component mechanism is correct, the non-monotonic forward-rapidity ρ00(pT) curve is understood as a transition from vorticity-polarized coalescence at low pT to primordial production at high pT, with no exotic spin physics needed.
- The same thermal vorticity produces a much weaker spin-alignment signal at mid-rapidity, so J/ψ spin alignment should be rapidity-dependent in a specific, kinematically predictable way.
- The high-pT plateau of ρ00 directly measures the primordial J/ψ spin state; its value can be compared across collision systems and energies as a probe of initial-production spin alignment.
- A positive low-pT deviation at mid-rapidity would signal an additional polarization source beyond thermal vorticity, such as the effective vector-field fluctuations introduced here.
Where Pith is reading between the lines
- The calibration scheme (⟨ω_y^2⟩ from the lowest-pT forward point, ρ_init=0.36 from the highest-pT point) means the forward-rapidity agreement is partly a reproduction; the genuinely testable prediction is the mid-rapidity shape and the sign of (ρ00−1/3) at low pT.
- If non-prompt J/ψ from B-hadron decays contribute non-negligibly in the measured inclusive sample, their different spin-alignment pattern would dilute the two-component interpretation; separate prompt and non-prompt measurements would settle this.
- The framework could be applied to other quarkonium states or to φ mesons: the rapidity-dependent suppression of vorticity-induced alignment is a generic kinematic effect, not specific to J/ψ, so the same signature should appear in other vector mesons.
- A stronger test would use a hydrodynamic model to compute ⟨ω_y^2(pT, Y)⟩ event-by-event instead of calibrating it to one data point; the resulting ρ00(pT) could then be compared with the data at both rapidities simultaneously.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper derives a non-relativistic approximation to the J/ψ spin density matrix element ρ00 from the relativistic spin Boltzmann equation, retaining thermal-vorticity and effective vector-field fluctuation contributions to the charm-quark polarization. The observed ρ00 is then written as a pT-weighted sum of a coalescence and an initial-production component (Eq. (32)), with the production fractions r_coal and r_init taken from a transport model. ⟨ω_y²⟩ is calibrated to the lowest-pT ALICE forward-rapidity point, and ρ_init^00 = 0.36 is fixed from the highest-pT point. With these inputs the model gives a non-monotonic forward-rapidity curve that follows the ALICE inclusive J/ψ data. The same framework is used to predict a strongly suppressed mid-rapidity effect and to illustrate the effect of vector-field fluctuations.
Significance. The theoretical construction is a useful extension of the MVSD/spin-Boltzmann formalism to heavy quarkonia, and the two-component production picture is physically motivated: coalescence-dominated low pT and primordial-dominated high pT naturally produce a non-monotonic ρ00 if the coalescence component inherits a spin-dependent polarization. The mid-rapidity prediction is a clear, falsifiable consequence. The paper is also transparent that ⟨ω_y²⟩ and ρ_init^00 are calibrated and that the vector-field parameters are illustrative. However, because both endpoints of the forward-rapidity curve are fixed by the same ALICE points that the paper claims to reproduce, and because the inclusive sample is modeled without non-prompt feed-down, the current comparison is significantly weaker than the abstract suggests. If these issues are addressed, the framework would be a valuable contribution.
major comments (4)
- [Sec. IV, Eq. (34) and the following paragraph] The forward-rapidity comparison is not a parameter-free reproduction. ⟨ω_y²⟩ is calibrated with the lowest-pT ALICE point, and ρ_init^00 is fixed to 0.36 from the highest-pT ALICE point. The observed curve is therefore an interpolation between two fitted endpoints; only the intermediate shape is a nontrivial consequence of Eq. (32) and r_coal/r_init. The abstract and Sec. I should be reworded from "well reproduced" to a statement that the model can accommodate the data, or the authors should provide a fit statistic and demonstrate that the shape is robust.
- [Abstract vs. Sec. IV and Fig. 2] There is a direct contradiction about the high-pT limit. The abstract states that primordial production causes ρ00 to approach 1/3, while Sec. IV and Fig. 2 show that the curve approaches the fitted input ρ_init^00 = 0.36. The value 0.36 is read off the data, not derived. The abstract must be corrected; otherwise the result misrepresents the paper's central claim.
- [Sec. IV, Eq. (32)] The ALICE sample is inclusive, but Eq. (32) contains only the two prompt components "assumed to provide the dominant prompt contribution." Non-prompt J/ψ from B-hadron decays are not included. At forward rapidity and pT ≳ 5 GeV, the non-prompt fraction is non-negligible. If non-prompt J/ψ have ρ00 ≈ 1/3, their admixture changes the high-pT inclusive ρ00 and the extracted ρ_init^00 (and indirectly ⟨ω_y²⟩). The authors should either quantify the feed-down fraction and extend Eq. (32), use a prompt-selected data set, or argue quantitatively that the effect is below the ~0.01 scale visible in Fig. 2.
- [Sec. IV, Fig. 1] The pT-dependent fractions r_coal and r_init are taken from the transport model of Ref. [48], and the intermediate non-monotonic shape is generated largely by these fractions. Only the nuclear-shadowing uncertainty is propagated; no uncertainty in the transport model itself, or in its applicability to ALICE 30–50% Pb-Pb at 5.02 TeV, is shown. This is load-bearing because the claimed agreement depends on these fractions. A sensitivity study or at least a discussion of the model's range of validity should be added.
minor comments (6)
- [Sec. IV] "data poin" is a typo for "data point."
- [Eq. (10)] The notation E_q^{p_T} is undefined; it should be E_q^p or explicitly defined as the on-shell energy.
- [Fig. 2, lower panel] The vertical axis label "(ρ00 − 1/3) × 10" should be clarified; the reader has to infer that the plotted quantity is dimensionless and amplified.
- [Eq. (29)] The large brackets in the f_T² and f_z² terms are hard to parse; please re-check the parentheses and define all symbols (F, G, Y) immediately before use.
- [Sec. III.B] The symbol g_V is reused for two different couplings (quark–meson and quark–effective-vector-field). Please use distinct notation to avoid confusion.
- [Sec. IV, Fig. 1] Ref. [48] is treated as a black box for r_coal/r_init. Provide enough detail or a short validation so that a reader can assess whether the fractions are appropriate for the ALICE centrality and rapidity window.
Circularity Check
Forward-rapidity 'reproduction' is anchored at both ends to the ALICE data it claims to explain: ⟨ω_y²⟩ is calibrated to the lowest-pT point and ρ_init^00 = 0.36 is fixed from the highest-pT point; the mid-rapidity curve is the only unanchored output.
specific steps
-
fitted input called prediction
[Sec. IV, Eq. (32)-(34) and paragraph after Eq. (34)]
"Since the averaged squared thermal vorticity is not determined independently in the present framework, we calibrate it using the lowest-pT forward-rapidity ALICE data poin."
In Eq. (32) the forward-rapidity observable is ρobs00(pT) = r_coal(pT) ρcoal00(pT) + r_init(pT) ρinit00 with r_coal ≈ 1 at low pT. The coalescence term is computed from Eq. (34) as ρcoal00 ≈ 1/3 − (1/9)⟨ω_y²⟩[F²+½G²+2FG v2(pT)]. Calibrating the single free parameter ⟨ω_y²⟩ to the lowest-pT ALICE point forces the model's low-pT value to equal that data point by construction; the claimed reproduction of the low-pT dip is therefore not an independent prediction.
-
fitted input called prediction
[Sec. IV, paragraph defining ρ_init^00 before Eq. (32) application; Abstract]
"we fix ρinit 00 = 0.36, motivated by the highest-pT ALICE data point, where the initial-production fraction is expected to dominate. This value should be regarded as an effective phenomenological input rather than a first-principles prediction of the primordial J/ψ spin alignment."
Since r_init(pT) → 1 at high pT, Eq. (32) makes ρobs00 → ρinit00 = 0.36, a value chosen from the highest-pT ALICE point; the high-pT plateau is an input, not a prediction. The Abstract's claim that ρ00 'approach[es] 1/3' at high pT is moreover not what the calculation does: Sec. IV states 'all curves approach the same limiting value ρinit00 = 0.36.' Only the intermediate non-monotonic shape and the mid-rapidity curve are unconstrained outputs.
full rationale
The spin-Boltzmann derivation (Secs. II-III) is an independent framework taken from external references [41, 29, 44], none of which are authored by the present authors, so no self-citation load-bearing circularity is present; the coalescence/initial fractions come from the external transport model [48]. The mid-rapidity result is a genuine prediction: the same calibrated ⟨ω_y²⟩ and external fractions are fed through the Lorentz-transformed Eq. (29), and the rapidity-dependent suppression is not anchored to mid-rapidity data. The vector-field parameter f is scanned for sensitivity, not fitted, so it does not add circularity. The circularity is confined to the forward-rapidity claim: Eq. (32) plus the two calibrations (⟨ω_y²⟩ tied to the lowest-pT ALICE point; ρ_init^00 = 0.36 tied to the highest-pT ALICE point) mean the endpoint behavior of the 'reproduced' curve is imposed by the same data set. The intermediate shape retains independent content from the transport-model fractions, and the paper is transparent that these are inputs, which keeps the score at partial (6) rather than full. The omission of non-prompt B-hadron feed-down in modeling the inclusive sample is a correctness risk, not a circularity, and is not scored here.
Axiom & Free-Parameter Ledger
free parameters (3)
- ⟨ω_y²⟩ (spacetime-averaged squared thermal vorticity along OAM) =
upper boundary 0.00466 used in Fig. 3; central calibrated value not quoted
- ρ_init^00 (primordial-production baseline) =
0.36
- f_T, f_z (vector-field fluctuation strengths) =
varied: f=0.01–0.05 (forward), f=0.03–0.30 (mid-rapidity)
axioms (8)
- domain assumption Vector-meson spin Boltzmann equation (Eq. 2) with coalescence/dissociation rates of Ref. [41] and formal solution f ≈ R_coal Δt (Eq. 8)
- domain assumption Quark polarization four-vector form, Eq. (10): P^μ ∝ (ω̃^μν ± …) p_ν [1−f]
- domain assumption Non-relativistic limit: m_q = m_qbar, m_V ≈ 2m_q, p^μ ≈ (m_V, 0) (Eq. 14)
- domain assumption Neglect of the electric part of thermal vorticity ε and of transverse vorticity components ω_x, ω_z (Sec. III, after Eq. 29)
- ad hoc to paper Two-component superposition, Eq. (32): ρ_obs^00 = r_coal ρ_coal^00 + r_init ρ_init^00, r_coal + r_init = 1
- ad hoc to paper Effective Abelian color-singlet vector field for gluon-field fluctuations (Sec. III.B, Sec. IV)
- domain assumption Only the y-component (global OAM direction) of quark polarization is retained (Sec. V)
- standard math Small-polarization expansion of Eq. (19) to Eq. (20): 1/(3+P·P̄) ≈ (1/3)(1 − P·P̄/3)
invented entities (1)
-
Effective Abelian-like vector field representing gluon-field fluctuations
no independent evidence
read the original abstract
We investigate the spin alignment of $J/\psi$ mesons in relativistic heavy-ion collisions within a two-component Boltzmann transport model. Starting from the relativistic spin Boltzmann equation, we derive the spin density matrix element $\rho_{00}$ under a non-relativistic approximation for heavy quarks. To interpret the recent ALICE measurements in Pb+Pb collisions, the observed $\rho_{00}$ is described as a $p_T$-dependent mixture of contributions from primordial production and the coalescence process. At forward rapidity, the $p_T$ dependence of charmonium $\rho_{00}$ is well reproduced by this two-component mechanism: at low $p_T$, charmonium production is dominated by the coalescence of partially polarized charm quarks induced by thermal vorticity; with increasing $p_T$, primordially produced charmonia become dominant, causing $\rho_{00}$ to approach $1/3$. To further test this spin alignment mechanism, we provide predictions for the $J/\psi$ $\rho_{00}$ in the mid-rapidity region, which exhibits a distinct $p_T$ trend due to the kinematic suppression of the thermal vorticity contribution. This study elucidates the underlying mechanism of $J/\psi$ spin alignment and advances our understanding of heavy quarkonium spin dynamics in strongly interacting matter.
Figures
Reference graph
Works this paper leans on
-
[1]
Deng and X.-G
W.-T. Deng and X.-G. Huang, Phys. Rev. C93, 064907 (2016)
2016
-
[2]
Wei, W.-T
D.-X. Wei, W.-T. Deng, and X.-G. Huang, Phys. Rev. C 99, 014905 (2019)
2019
-
[3]
Huang et al., in Strongly Interacting Matter under Rotation (Springer International Publishing, Cham, 2021) pp
X.-G. Huang et al., in Strongly Interacting Matter under Rotation (Springer International Publishing, Cham, 2021) pp. 281–308
2021
-
[4]
Becattini et al., Eur
F. Becattini et al., Eur. Phys. J. C75, 406 (2015)
2015
-
[5]
Deng and X.-G
W.-T. Deng and X.-G. Huang, Phys. Rev. C85, 044907 (2012)
2012
-
[6]
V. V. Skokov, A. Yu. Illarionov, and V. D. Toneev, Int. J. Modern Phys. A24, 5925 (2009)
2009
-
[7]
Tuchin, Adv
K. Tuchin, Adv. High Energy Phys.2013, 1 (2013)
2013
-
[8]
Voronyuk et al., Phys
V. Voronyuk et al., Phys. Rev. C83, 054911 (2011)
2011
-
[9]
Buzzegoli, Nucl
M. Buzzegoli, Nucl. Phys. A1036, 122674 (2023)
2023
-
[10]
Becattini, V
F. Becattini, V. Chandra, L. Del Zanna, et al., Ann. Physics338, 32 (2013)
2013
-
[11]
Becattini, I
F. Becattini, I. Karpenko, M. A. Lisa, et al., Phys. Rev. C95, 054902 (2017)
2017
-
[12]
Becattini and M
F. Becattini and M. A. Lisa, Annu. Rev. Nucl. Part. Sci. 70, 395 (2020)
2020
-
[13]
B. Betz, M. Gyulassy, and G. Torrieri, Phys. Rev. C: Nucl. Phys.76, 44901 (2007), arXiv:0708.0035 [nucl-th]
Pith/arXiv arXiv 2007
-
[14]
Fang, L.-g
R.-h. Fang, L.-g. Pang, Q. Wang, et al., Phys. Rev. C94, 024904 (2016)
2016
-
[15]
K. G. Wilson, Phys. Rev. D10, 2445 (1974)
1974
-
[16]
The STAR Collaboration, Nature548, 62 (2017)
2017
-
[17]
Adam et al., Phys
J. Adam et al., Phys. Rev. C98, 014910 (2018)
2018
-
[18]
Adam et al., Phys
J. Adam et al., Phys. Rev. Lett.126, 162301 (2021)
2021
-
[19]
Sheng, L
X.-L. Sheng, L. Oliva, and Q. Wang, Phys. Rev. D101, 096005 (2020)
2020
-
[20]
Liang and X.-N
Z.-T. Liang and X.-N. Wang, Phys. Rev. Lett.94, 102301 (2005)
2005
-
[21]
Liang and X.-N
Z.-T. Liang and X.-N. Wang, Phys. Lett. B629, 20 (2005)
2005
-
[22]
Lv, Z.-h
J.-p. Lv, Z.-h. Yu, Z.-t. Liang, et al., Phys. Rev. D109, 114003 (2024)
2024
-
[23]
Xu and M
K. Xu and M. Huang, Phys. Rev. D110, 094034 (2024)
2024
-
[24]
Sheng, Q
X.-L. Sheng, Q. Wang, and X.-N. Wang, Phys. Rev. D 102, 056013 (2020)
2020
-
[25]
Yang, R.-h
Y.-g. Yang, R.-h. Fang, Q. Wang, et al., Phys. Rev. C 97, 034917 (2018)
2018
-
[26]
Chen, W.-j
H.-L. Chen, W.-j. Fu, X.-G. Huang, et al., Phys. Rev. Lett.135, 032302 (2025)
2025
-
[27]
STAR Collaboration et al., Nature614, 244 (2023)
2023
-
[28]
B. I. Abelev et al., Phys. Rev. Lett.99, 112301 (2007)
2007
-
[29]
Sheng, L
X.-L. Sheng, L. Oliva, Z.-T. Liang, et al., Phys. Rev. Lett.131, 042304 (2023)
2023
-
[30]
P. H. De Moura, K. J. Gon¸ calves, and G. Torrieri, Phys. Rev. D108, 34032 (2023)
2023
-
[31]
J. Zhao, K. Zhou, S. Chen, et al., Prog. Part. Nucl. Phys. 114, 103801 (2020)
2020
-
[32]
H. Satz, J. Phys. G: Nucl. Part. Phys.32, R25 (2006)
2006
-
[33]
Rothkopf, Phys
A. Rothkopf, Phys. Rep.858, 1 (2020)
2020
-
[34]
ALICE Collaboration et al., Eur. Phys. J. C78, 562 (2018)
2018
-
[35]
Andronic et al., Eur
A. Andronic et al., Eur. Phys. J. C76, 107 (2016)
2016
-
[36]
Adam et al
J. Adam et al. (STAR Collaboration), Phys. Rev. D102, 092009 (2020)
2020
-
[37]
Acharya et al., Phys
S. Acharya et al., Phys. Rev. Lett.131, 042303 (2023)
2023
-
[38]
High Energy Phys
The ALICE collaboration et al., J. High Energy Phys. 2025(10), 94
2025
-
[39]
Chou, Z.-b
K.-c. Chou, Z.-b. Su, B.-l. Hao, et al., Phys. Rep.118, 1 (1985)
1985
-
[40]
Cassing, Eur
W. Cassing, Eur. Phys. J. Spec. Top.168, 3 (2009)
2009
-
[41]
Sheng, L
X.-L. Sheng, L. Oliva, Z.-T. Liang, et al., Phys. Rev. D 109, 036004 (2024)
2024
-
[42]
Sheng, N
X.-L. Sheng, N. Weickgenannt, E. Speranza, et al., Phys. Rev. D104, 016029 (2021)
2021
-
[43]
Xu et al., Eur
Y.-Z. Xu et al., Eur. Phys. J. C81, 895 (2021)
2021
-
[44]
Sheng, S
X.-L. Sheng, S. Pu, and Q. Wang, Phys. Rev. C108, 054902 (2023)
2023
-
[45]
High Energy Phys
The ALICE collaboration et al., J. High Energy Phys. 2020(10), 141
2020
-
[46]
L.-G. Pang, H. Petersen, and X.-N. Wang, Phys. Rev. C 97, 064918 (2018)
2018
-
[47]
L. Pang, Q. Wang, and X.-N. Wang, Phys. Rev. C86, 024911 (2012)
2012
-
[48]
K. Zhou, N. Xu, Z. Xu, et al., Phys. Rev. C89, 54911 (2014)
2014
-
[49]
Eskola, H
K.J. Eskola, H. Paukkunen, and C.A. Salgado, J. High Energy Phys.2009(4), 65
2009
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.