Recognition: unknown
Quantum spin dynamics of heavy quarks and polarization observables in relativistic heavy-ion collisions
Pith reviewed 2026-05-08 08:24 UTC · model grok-4.3
The pith
A density-matrix framework evolves heavy-quark polarization from initial magnetic alignment to final hadron observables.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The central claim is that heavy-quark spin polarization obeys a density-matrix evolution equation whose analytic solution, when combined with a fragmentation-based hadronization prescription, reproduces the observed D* spin alignment and thereby determines an effective depolarization strength that sets the spin-relaxation timescale and yields benchmark predictions for Lambda_c polarization.
What carries the argument
The quantum spin-density-matrix evolution equation for the heavy-quark polarization vector, solved analytically and mapped to observables by coupling the heavy-quark spin to a light antiquark during fragmentation.
Load-bearing premise
The fragmentation process couples the heavy quark spin to the light antiquark spin in a way that directly determines the observed vector-meson rho_00 and baryon polarization.
What would settle it
A precision measurement of Lambda_c polarization that lies outside the range allowed by the depolarization strength fitted to ALICE D* data would falsify the spin-transfer assumption in the hadronization prescription.
Figures
read the original abstract
We develop a quantum spin-density-matrix framework for heavy-quark spin dynamics in relativistic heavy-ion collisions. Starting from an initial polarization induced along the magnetic-field direction, we derive the evolution equation for spin polarization within this framework and obtain analytic solutions. The evolved polarization is connected to open heavy-flavor observables via a fragmentation-based hadronization prescription. For vector mesons, the spin-alignment parameter $\rho_{00}$ is constructed by coupling the heavy-quark spin to that of the light antiquark produced during fragmentation. We confront our results with recent ALICE measurements of prompt $D^{*+}$ spin alignment in Pb--Pb collisions at $\sqrt{s_{\rm NN}}=5.02~{\rm TeV}$ and extract an effective depolarization strength that determines the spin-relaxation time scale. Using this fitted parameter, we provide benchmark estimates for $\Lambda_c^+$ and $\bar{\Lambda}_c^-$ polarization, up to an overall spin-transfer normalization. We further estimate the recently proposed elliptic polarization harmonic arising from path-length-dependent depolarization in an anisotropic fireball.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript develops a quantum spin-density-matrix framework for heavy-quark spin dynamics in relativistic heavy-ion collisions. Starting from initial polarization along the magnetic-field direction, it derives the evolution equation and obtains analytic solutions for the time-dependent spin polarization. These are connected to observables via a fragmentation-based hadronization prescription in which the heavy-quark spin density matrix is coupled to an assumed light-antiquark spin state to construct the vector-meson alignment parameter rho_00 and analogous expressions for Lambda_c baryon polarization. The model is confronted with ALICE prompt D*+ rho_00 data in Pb-Pb collisions at 5.02 TeV to extract an effective depolarization strength; this single fitted parameter is then used to provide benchmark estimates for Lambda_c^+ and anti-Lambda_c^- polarization (up to an overall normalization) and for the elliptic polarization harmonic arising from path-length-dependent depolarization in an anisotropic fireball.
Significance. If the fragmentation prescription holds, the work supplies a transparent analytic framework for extracting a spin-relaxation time scale from data and generating predictions for other heavy-flavor polarization observables. The derivation of closed-form analytic solutions for the density-matrix evolution is a genuine strength, as it makes the dependence on the depolarization parameter explicit and falsifiable. Timely comparison with recent ALICE measurements adds phenomenological relevance. However, the quantitative reach is constrained by the single-parameter fit and the untested hadronization assumption, limiting the result to a useful parametrization rather than a fully independent prediction.
major comments (2)
- [Hadronization and observables section] Hadronization and observables section: The fragmentation-based prescription that couples the heavy-quark spin density matrix to the light antiquark spin state to obtain rho_00 (and the analogous construction for Lambda_c polarization) is the load-bearing step for all quantitative results. No sensitivity tests, alternative coupling rules, or comparison to lattice or other hadronization models are provided to justify this specific spin-coupling assumption. If light-quark spin correlations or additional depolarization during fragmentation are present, the extracted depolarization strength ceases to represent purely medium-induced relaxation, and the subsequent Lambda_c estimates become uncontrolled extrapolations.
- [Results and data comparison section] Results and data comparison section: The central quantitative claims rest on fitting one effective depolarization parameter to the ALICE D*+ rho_00 data and directly inserting the same value into the expressions for Lambda_c polarization and the elliptic harmonic. No independent cross-check against other data sets, no error propagation from the fit, and no alternative fits are shown; by construction the Lambda_c and elliptic-harmonic results are therefore scaled versions of the fitted quantity rather than independent predictions.
minor comments (2)
- [Abstract] The abstract states that Lambda_c estimates are given 'up to an overall spin-transfer normalization'; this important caveat should be quantified with an estimated uncertainty range and discussed more explicitly in the main text when presenting the benchmark values.
- [Framework section] Notation for the depolarization strength and the spin-relaxation time scale should be introduced with a clear equation reference when first defined, to avoid ambiguity when the same parameter is reused across different observables.
Simulated Author's Rebuttal
We thank the referee for the careful reading of our manuscript and the constructive comments. We appreciate the recognition of the analytic solutions and the timely comparison with ALICE data. We address the two major comments point by point below, indicating the revisions we will incorporate.
read point-by-point responses
-
Referee: [Hadronization and observables section] The fragmentation-based prescription that couples the heavy-quark spin density matrix to the light antiquark spin state to obtain rho_00 (and the analogous construction for Lambda_c polarization) is the load-bearing step for all quantitative results. No sensitivity tests, alternative coupling rules, or comparison to lattice or other hadronization models are provided to justify this specific spin-coupling assumption. If light-quark spin correlations or additional depolarization during fragmentation are present, the extracted depolarization strength ceases to represent purely medium-induced relaxation, and the subsequent Lambda_c estimates become uncontrolled extrapolations.
Authors: We agree that the fragmentation prescription is a central assumption underlying the connection to observables and that the current manuscript provides no sensitivity tests or comparisons to alternative models. The prescription follows the standard fragmentation picture in which the heavy quark combines with a light antiquark whose spin state is taken to be uncorrelated on average. We acknowledge that light-quark spin correlations or additional depolarization during hadronization would render the extracted parameter effective rather than purely medium-induced. In the revised version we will expand the relevant section to state the assumption explicitly, discuss its uncertainties, and clarify that the depolarization strength should be interpreted as effective. We will also add a forward-looking remark on the desirability of future comparisons with more advanced hadronization models. revision: partial
-
Referee: [Results and data comparison section] The central quantitative claims rest on fitting one effective depolarization parameter to the ALICE D*+ rho_00 data and directly inserting the same value into the expressions for Lambda_c polarization and the elliptic harmonic. No independent cross-check against other data sets, no error propagation from the fit, and no alternative fits are shown; by construction the Lambda_c and elliptic-harmonic results are therefore scaled versions of the fitted quantity rather than independent predictions.
Authors: The referee is correct that the Lambda_c and elliptic-harmonic estimates are obtained by direct substitution of the single parameter fitted to the D*+ data. This construction is intentional: the manuscript presents a unified analytic framework in which the parameter extracted from the only currently available vector-meson polarization data is used to generate benchmark estimates for other observables. We do not present these estimates as independent predictions. In the revision we will add explicit error propagation from the fit, include a dedicated paragraph on the limitations arising from the single-parameter fit and the absence of independent cross-checks, and emphasize that the results remain conditional on the hadronization prescription. revision: partial
Circularity Check
Depolarization strength fitted to D*+ rho_00 data then inserted directly into Lambda_c polarization estimates
specific steps
-
fitted input called prediction
[Abstract and results section (confrontation with ALICE data)]
"We confront our results with recent ALICE measurements of prompt D*+ spin alignment in Pb-Pb collisions at sqrt(s_NN)=5.02 TeV and extract an effective depolarization strength that determines the spin-relaxation time scale. Using this fitted parameter, we provide benchmark estimates for Lambda_c+ and anti-Lambda_c- polarization, up to an overall spin-transfer normalization."
The depolarization strength is obtained by fitting the model's rho_00 prediction to the measured D*+ alignment; the same numerical value is then inserted into the Lambda_c polarization formula. The Lambda_c 'estimates' are therefore proportional to the fitted parameter by construction and do not constitute an independent test of the spin-evolution dynamics.
full rationale
The paper derives an evolution equation for the heavy-quark spin density matrix, obtains analytic solutions, and connects them to observables via a fragmentation prescription that couples heavy-quark spin to light antiquark spin. The central load-bearing step is the extraction of a single effective depolarization parameter from ALICE D*+ data; this fitted value is then substituted into the expressions for Lambda_c polarization (and the elliptic harmonic). By the paper's own construction the resulting 'benchmark estimates' are scaled versions of the fitted quantity rather than independent predictions. The fragmentation coupling rule itself is an ansatz, but the circularity flag is triggered specifically by the fit-then-predict structure, not by the ansatz alone. No self-citation chain or self-definitional loop is required for the reduction.
Axiom & Free-Parameter Ledger
free parameters (1)
- effective depolarization strength
axioms (1)
- domain assumption Fragmentation-based hadronization prescription that couples heavy-quark spin to light-antiquark spin
Forward citations
Cited by 1 Pith paper
-
Spin dynamics and polarization in relativistic systems: recent developments
The review summarizes developments in spin hydrodynamics, polarization from spin-vorticity coupling, pseudo-gauge freedom, and heavy-flavor spin dynamics in relativistic systems.
Reference graph
Works this paper leans on
-
[1]
E.V.Shuryak,TheoryofHadronicPlasma,Sov.Phys.JETP47(1978) 212–219
1978
-
[2]
Collective flow and viscosity in relativistic heavy-ion collisions.Ann
U. Heinz, R. Snellings, Collective flow and viscosity in relativistic heavy-ion collisions, Ann. Rev. Nucl. Part. Sci. 63 (2013) 123–151. arXiv:1301.2826,doi:10.1146/annurev-nucl-102212-170540
-
[3]
D. Everett, et al., Multisystem Bayesian constraints on the transport coefficients of QCD matter, Phys. Rev. C 103 (5) (2021) 054904. arXiv:2011.01430,doi:10.1103/PhysRevC.103.054904
-
[4]
B.Svetitsky,Diffusionofcharmedquarksinthequark-gluonplasma, Phys. Rev. D 37 (1988) 2484–2491.doi:10.1103/PhysRevD.37.2484
-
[5]
G. D. Moore, D. Teaney, How much do heavy quarks thermalize in a heavy ion collision?, Phys. Rev. C 71 (2005) 064904.arXiv: hep-ph/0412346,doi:10.1103/PhysRevC.71.064904
-
[6]
H. van Hees, V. Greco, R. Rapp, Heavy-quark probes of the quark- gluon plasma and interpretation of recent data taken at the BNL Relativistic Heavy Ion Collider, Phys. Rev. C 73 (2006) 034913. arXiv:nucl-th/0508055,doi:10.1103/PhysRevC.73.034913
-
[7]
arXiv:hep-th/0605182,doi:10.1103/PhysRevD.74.126005
S.S.Gubser,DragforceinAdS/CFT,Phys.Rev.D74(2006)126005. arXiv:hep-th/0605182,doi:10.1103/PhysRevD.74.126005
-
[8]
H. van Hees, M. Mannarelli, V. Greco, R. Rapp, Nonperturbative heavy-quark diffusion in the quark-gluon plasma, Phys. Rev. Lett. 100 (2008) 192301.arXiv:0709.2884,doi:10.1103/PhysRevLett.100. 192301
-
[9]
Y. Akamatsu, T. Hatsuda, T. Hirano, Heavy Quark Diffusion with RelativisticLangevinDynamicsintheQuark-GluonFluid,Phys.Rev. C 79 (2009) 054907.arXiv:0809.1499,doi:10.1103/PhysRevC.79. 054907
-
[10]
S. K. Das, J.-e. Alam, P. Mohanty, Dragging Heavy Quarks in Quark Gluon Plasma at the Large Hadron Collider, Phys. Rev. C 82 (2010) 014908.arXiv:1003.5508,doi:10.1103/PhysRevC.82.014908
-
[11]
D.Banerjee,S.Datta,R.Gavai,P.Majumdar,HeavyQuarkMomen- tum Diffusion Coefficient from Lattice QCD, Phys. Rev. D 85 (2012) 014510.arXiv:1109.5738,doi:10.1103/PhysRevD.85.014510
-
[12]
70 (2014) 00061.arXiv:1210.0292, doi:10.1051/epjconf/20147000061
H.T.Ding,A.Francis,O.Kaczmarek,F.Karsch,H.Satz,W.Söldner, Charmonium dissociation and heavy quark transport in hot quenched lattice QCD, EPJ Web Conf. 70 (2014) 00061.arXiv:1210.0292, doi:10.1051/epjconf/20147000061
-
[13]
S. K. Das, F. Scardina, S. Plumari, V. Greco, Toward a solution to the raa and v2 puzzle for heavy quarks, Physics Letters B 747 (2015) 260–264
2015
-
[14]
L. Adamczyk, et al., Measurement of𝐷 0 Azimuthal Anisotropy at Midrapidity in Au+Au Collisions at √𝑠𝑁𝑁=200 GeV, Phys. Rev. Lett. 118 (21) (2017) 212301.arXiv:1701.06060,doi:10.1103/ PhysRevLett.118.212301
-
[15]
X.Dong,V.Greco,Heavyquarkproductionandpropertiesofquark– gluon plasma, Progress in Particle and Nuclear Physics 104 (2019) 97–141
2019
-
[16]
S. Acharya, et al., Prompt D0, D+, and D∗+ production in Pb–Pb collisions at √𝑠NN = 5.02 TeV, JHEP 01 (2022) 174.arXiv:2110. 09420,doi:10.1007/JHEP01(2022)174
-
[17]
Braaten, M
E. Braaten, M. H. Thoma, Energy loss of a heavy quark in the quark - gluon plasma, Phys. Rev. D 44 (9) (1991) R2625.doi:10.1103/ PhysRevD.44.R2625
1991
-
[18]
M. Golam Mustafa, D. Pal, D. Kumar Srivastava, Propagation of charm quarks in equilibrating quark - gluon plasma, Phys. Rev. C 57 (1998) 889–898, [Erratum: Phys.Rev.C 57, 3499–3499 (1998)]. arXiv:nucl-th/9706001,doi:10.1103/PhysRevC.57.3499
-
[19]
P. Romatschke, M. Strickland, Collisional energy loss of a heavy quark in an anisotropic quark-gluon plasma, Phys. Rev. D 71 (2005) 125008.arXiv:hep-ph/0408275,doi:10.1103/PhysRevD.71.125008
-
[20]
H. van Hees, R. Rapp, Thermalization of heavy quarks in the quark- gluonplasma,Phys.Rev.C71(2005)034907.arXiv:nucl-th/0412015, doi:10.1103/PhysRevC.71.034907
-
[21]
P. P. Bhaduri, N. Borghini, A. Jaiswal, M. Strickland, Anisotropic escape mechanism and elliptic flow of bottomonia, Phys. Rev. C 100(5)(2019)051901.arXiv:1809.06235,doi:10.1103/PhysRevC.100. 051901
-
[22]
P. P. Bhaduri, M. Alqahtani, N. Borghini, A. Jaiswal, M. Strickland, Fireball tomography from bottomonia elliptic flow in relativistic heavy-ion collisions, Eur. Phys. J. C 81 (7) (2021) 585.arXiv: 2007.03939,doi:10.1140/epjc/s10052-021-09383-3
-
[23]
A. Andronic, P. Braun-Munzinger, M. K. Köhler, A. Mazeliauskas, K. Redlich, J. Stachel, V. Vislavicius, The multiple-charm hierarchy in the statistical hadronization model, JHEP 07 (2021) 035.arXiv: 2104.12754,doi:10.1007/JHEP07(2021)035
-
[24]
D. Kumar, N. Sarkar, P. P. Bhaduri, A. Jaiswal, Examination of thermalization of quarkonia at energies available at the CERN Large Hadron Collider, Phys. Rev. C 107 (6) (2023) 064906.arXiv:2303. 02900,doi:10.1103/PhysRevC.107.064906
-
[25]
A. Andronic, et al., Comparative study of quarkonium transport in hot QCD matter, Eur. Phys. J. A 60 (4) (2024) 88.arXiv:2402.04366, doi:10.1140/epja/s10050-024-01306-6
-
[26]
D. E. Kharzeev, L. D. McLerran, H. J. Warringa, The Effects of topological charge change in heavy ion collisions: ’Event by event P and CP violation’, Nucl. Phys. A803 (2008) 227–253.arXiv: 0711.0950,doi:10.1016/j.nuclphysa.2008.02.298
-
[27]
V.Skokov,A.Y.Illarionov,V.Toneev,Estimateofthemagneticfield strengthinheavy-ioncollisions,Int.J.Mod.Phys.A24(2009)5925– 5932.arXiv:0907.1396,doi:10.1142/S0217751X09047570
-
[28]
V. Voronyuk, V. D. Toneev, W. Cassing, E. L. Bratkovskaya, V. P. Konchakovski, S. A. Voloshin, (Electro-)Magnetic field evolution in relativistic heavy-ion collisions, Phys. Rev. C 83 (2011) 054911. arXiv:1103.4239,doi:10.1103/PhysRevC.83.054911
-
[29]
J. Bloczynski, X.-G. Huang, X. Zhang, J. Liao, Azimuthally fluc- tuating magnetic field and its impacts on observables in heavy-ion collisions, Phys. Lett. B 718 (2013) 1529–1535.arXiv:1209.6594, doi:10.1016/j.physletb.2012.12.030
-
[30]
K. Tuchin, Particle production in strong electromagnetic fields in relativisticheavy-ioncollisions,Adv.HighEnergyPhys.2013(2013) 490495.arXiv:1301.0099,doi:10.1155/2013/490495
-
[31]
L.McLerran,V.Skokov,CommentsAbouttheElectromagneticField in Heavy-Ion Collisions, Nucl. Phys. A 929 (2014) 184–190.arXiv: 1305.0774,doi:10.1016/j.nuclphysa.2014.05.008
-
[32]
A. Huang, Y. Jiang, S. Shi, J. Liao, P. Zhuang, Out-of-equilibrium chiral magnetic effect from chiral kinetic theory, Phys. Lett. B 777 (2018) 177–183.arXiv:1703.08856,doi:10.1016/j.physletb.2017.12. 025
-
[33]
Dynamical magnetic fields in heavy-ion collisions,
A. Huang, D. She, S. Shi, M. Huang, J. Liao, Dynamical magnetic fields in heavy-ion collisions, Phys. Rev. C 107 (3) (2023) 034901. arXiv:2212.08579,doi:10.1103/PhysRevC.107.034901
-
[34]
Z.-F. Jiang, S. Cao, W.-J. Xing, X.-Y. Wu, C. B. Yang, B.-W. Zhang, Probing the initial longitudinal density profile and electromagnetic field inultrarelativistic heavy-ioncollisions with heavyquarks, Phys. Rev. C 105 (5) (2022) 054907.arXiv:2202.13555,doi:10.1103/ PhysRevC.105.054907. :Preprint submitted to Elsevier Page 8 of 9 Quantum spin dynamics of h...
-
[35]
S. Dey, A. Jaiswal, Rotational Brownian motion and heavy quark polarization in QCD medium, Phys. Lett. B 873 (2026) 140202. arXiv:2502.20352,doi:10.1016/j.physletb.2026.140202
-
[36]
Heavy quark polarization anisotropy as a novel probe of fireball geometry,
A. Jaiswal, Heavy quark polarization anisotropy as a novel probe of fireball geometry (1 2026).arXiv:2601.22882
-
[37]
S. Acharya, et al., First measurement of D∗+ vector meson spin alignmentinPb–Pbcollisionsat √𝑠NN = 5.02TeV,JHEP10(2025) 094.arXiv:2504.00714,doi:10.1007/JHEP10(2025)094
-
[38]
J. J. Sakurai, J. Napolitano, Modern quantum mechanics., Addison- Wesley, 1994
1994
-
[39]
E.Leader,SpininParticlePhysics,CambridgeMonographsonParti- cle Physics, Nuclear Physics and Cosmology, Cambridge University Press, 2001
2001
-
[40]
Global Λ hyperon polarization in nuclear col- lisions: evidence for the most vortical fluid
L. Adamczyk, et al., GlobalΛhyperon polarization in nuclear colli- sions: evidence for the most vortical fluid, Nature 548 (2017) 62–65. arXiv:1701.06657,doi:10.1038/nature23004
-
[41]
Z.-T. Liang, X.-N. Wang, Spin alignment of vector mesons in non- central A+A collisions, Phys. Lett. B 629 (2005) 20–26.arXiv: nucl-th/0411101,doi:10.1016/j.physletb.2005.09.060
-
[42]
Y.-G. Yang, R.-H. Fang, Q. Wang, X.-N. Wang, Quark coalescence model for polarized vector mesons and baryons, Phys. Rev. C 97 (3) (2018) 034917.arXiv:1711.06008,doi:10.1103/PhysRevC.97.034917
-
[43]
A.V.Manohar,M.B.Wise,Heavyquarkphysics,Vol.10,Cambridge University Press, 2000.doi:10.1017/9781009402125
-
[44]
R. S. Bhalerao, Collectivity in large and small systems formed in ultrarelativistic collisions, Eur. Phys. J. ST 230 (3) (2021) 635–654. arXiv:2009.09586,doi:10.1140/epjs/s11734-021-00019-x
-
[45]
C.Loizides,J.Kamin,D.d’Enterria,ImprovedMonteCarloGlauber predictionsatpresentandfuturenuclearcolliders,Phys.Rev.C97(5) (2018) 054910, [Erratum: Phys.Rev.C 99, 019901 (2019)].arXiv: 1710.07098,doi:10.1103/PhysRevC.97.054910
-
[46]
:Preprint submitted to Elsevier Page 9 of 9
BayesianAnalysisofNuclearDynamics(BAND)Frameworkproject, 2020,https://bandframework.github.io/. :Preprint submitted to Elsevier Page 9 of 9
2020
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.