Quarkonia Spin Alignment

• Quarkonia spin alignment is defined by anisotropic spin distributions measured via dilepton decay angular patterns, directly probing heavy quark production mechanisms.
• Experimental analysis employs reference frames like the helicity and Collins–Soper to extract polarization parameters that distinguish between competing QCD models.
• Vorticity, magnetic fields, and QGP effects induce spin-dependent dissociation and regeneration, offering insights into non-equilibrium spin transport in extreme conditions.

Quarkonia spin alignment refers to the quantum statistical anisotropy in the spin orientation of heavy quarkonium states (bound states of heavy quark–antiquark pairs, such as $J/\psi$, $\Upsilon$) produced in high-energy collisions. The phenomenon manifests as deviations in the angular distribution of decay products, revealing underlying production mechanisms, in-medium effects, and QCD dynamics in both proton–proton and heavy-ion environments.

1. Definition and Experimental Observables

Quarkonium spin alignment is probed through the angular distribution of the dilepton pairs ($\mu^+\mu^-$ or $e^+e^-$) from vector quarkonium decays. In the rest frame of the parent state, the general two-dimensional angular distribution is given by

$$\frac{d^2N}{d\cos\theta\,d\phi} = N_0 \left[1 + \lambda_\theta\cos^2\theta + \lambda_\phi\sin^2\theta\cos2\phi + \lambda_{\theta\phi}\sin2\theta\cos\phi\right]$$

where $\theta$ and $\phi$ are the polar and azimuthal angles of the positive lepton with respect to a defined quantization axis, and $\lambda_\theta,\lambda_\phi,\lambda_{\theta\phi}$ are the spin alignment (polarization) parameters (Etzion, 2010). In practice, the one-dimensional projection

$$\frac{dN}{d\cos\theta^*} \propto 1 + \alpha \cos^2\theta^*$$

is frequently used, with $$\alpha \equiv (\sigma_T-2\sigma_L)/(\sigma_T+2\sigma_L)$$, relating transverse (T) and longitudinal (L) polarizations.

The spin-density matrix $\rho_{mm'}$ (with $m=\pm1,0$ in the spin-1 basis) parameterizes the ensemble. The longitudinal alignment observable is $\rho_{00}$ (population in $m=0$ state), directly related to $\lambda_\theta$ via

$$\rho_{00} = \frac{1 - \lambda_\theta}{3 + \lambda_\theta}$$

By default, $\rho_{00}=1/3$ signals no alignment. Deviations quantify the net alignment along the quantization axis.

2. Production Mechanisms and Reference Frames

Quarkonium production models directly impact expected spin alignment. The helicity and Collins–Soper frames are standard reference axes: in the helicity frame, the $z$-axis aligns with the quarkonium momentum in the laboratory, while the Collins–Soper frame uses the bisector of the colliding beam directions (Etzion, 2010).

• Color Singlet Model (CSM): Predicts weak or slightly longitudinal polarization, i.e., $\alpha\lesssim 0$ at all $p_T$.
• NRQCD Color Octet Model (COM): Anticipates strong transverse polarization at high $p_T$, with $\alpha \to +1$. Experimental data, however, deviates from COM expectations, often exhibiting near-zero or negative $\alpha$ even up to substantial $p_T$ (Etzion, 2010).

Frame choice is critical for unambiguous interpretation. Nontrivial differences in $\lambda_\theta,\lambda_\phi$ between frames can elucidate the QCD production dynamics.

3. Medium-Induced Spin Alignment: Vorticity and Magnetic Fields

Heavy-ion collisions create extreme conditions where rotation (fluid vorticity) and intense electromagnetic fields imprint further spin dynamics.

a. Vortical Effects

Quarkonium in a rotating QGP experiences a Hamiltonian: $H = \frac{p^2}{2\mu} + V(r) - \omega (L_z + S_z)$ where $\omega$ is the vorticity along the reaction plane normal, and $S_z$ the spin projection (Sahoo et al., 11 Jun 2025, Moura et al., 2023, Sahoo et al., 21 Dec 2025).

• The $-\omega S_z$ term splits the triplet into $m_j=\pm1,0$ substates, favoring alignment along the vorticity axis.
• The equilibrium spin density matrix is then

$$\rho_{mm} = \frac{e^{\beta m \omega}}{\sum_{m'=-1}^1 e^{\beta m' \omega}}$$

leading to $\rho_{00}\ne1/3$.

• Calculations show $\rho_{00}>1/3$ for $J/\psi$ and $<1/3$ for excited states like $\psi(2S)$ under strong vorticity (Sahoo et al., 11 Jun 2025, Sahoo et al., 21 Dec 2025).

b. Magnetic Field Effects

Strong magnetic fields induce both orbital deformation (diamagnetic/Stark) and, dominantly, spin-state mixing (Zeeman). The Zeeman effect mixes spin-triplet and singlet states, producing: $$\lambda_\theta^{\rm spin} \simeq \frac{1}{2}\left[\frac{Q B}{m_Q \Delta E}\right]^2$$ with $Q$ the heavy-quark charge and $\Delta E$ the triplet-singlet splitting. For representative $eB\sim 0.2$ GeV$^2$, $\lambda_\theta^{\rm spin}\sim 10^{-3}$, dominating over the orbital contributions ($\lesssim10^{-5}$) (Yan et al., 30 Jul 2025).

4. Dissociation and Spin-Dependent Dynamics in QGP

The QGP environment modifies quarkonium via dissociation, with recent kinetic and quantum open-systems approaches elucidating the spin dependence.

a. Polarized Dissociation

• Spin-dependent dissociation rates arise from chromomagnetic couplings in pNRQCD, especially the magnetic-dipole (M1) transitions.
• Both leading-order gluo-dissociation and inelastic Coulomb scattering yield suppressed dissociation for $m=0$ compared to $m=\pm1$ states in moving quarkonia, leading to $\rho_{00} > 1/3$ for $J/\psi$ at moderate $p_T$. Regeneration effects from recombination counteract, possibly inverting the sign at low $p_T$ (Chen et al., 28 Jan 2025).
• In a vortical QGP, light quark/gluon polarization induces further spin splitting in dissociation rates, generally enhancing the survival probability of $m=0$ and resulting in $\rho_{00}<1/3$, consistent with ALICE and CMS observations (Liang et al., 9 Feb 2025, Sahoo et al., 21 Dec 2025).

b. Global Spin Correlations and Spin-Transport Kinetics

• The open quantum systems + pNRQCD framework yields Boltzmann and Lindblad equations with explicit chromomagnetic correlators governing polarization dynamics (Yang et al., 26 Feb 2025, Yang et al., 2024).
• The spin-1 density matrix elements $f_{mm'}$ evolve according to spin-dependent dissociation and recombination kernels. Chromoelectric field correlators control energy-momentum relaxation; chromomagnetic correlators drive polarization evolution.
• In the high-temperature limit, Lindblad evolution generates exponential relaxation toward $\rho_{00}=1/3$, but anisotropic correlators or initial conditions can sustain $\rho_{00}\neq1/3$ at freeze-out.

5. Experimental Extraction and Impact on QCD Diagnostics

Spin alignment is extracted via fits to the corrected angular decay distributions, applying kinematic acceptance and efficiency corrections derived from data and MC. In ATLAS, di-muon triggers, vertexing, and prompt selection suppress backgrounds, allowing fits for $\alpha$ in $p_T$ bins; present uncertainties on $\alpha$ are typically of order $0.03$–$0.17$ per bin for $J/\psi$ in early LHC datasets (Etzion, 2010).

Spin alignment parameters are sensitive to production model discrimination: e.g., no significant rise of $\alpha$ with $p_T$ up to 20 GeV disfavors NRQCD color-octet dominance (Etzion, 2010). In heavy-ion collisions, quarkonium spin alignment becomes a probe of QGP vorticity, magnetic field, and anisotropic structure. The sign and magnitude of $\rho_{00}-1/3$—positive for velocity/dissociation-driven, negative for vorticity-driven—distinguish mechanisms (Chen et al., 28 Jan 2025, Liang et al., 9 Feb 2025, Wei et al., 2023).

Combined measurements of $\rho_{00}$ for ground and excited states, as functions of $p_T$, centrality, and event-plane angle, map the emergent QGP spin structure. Off-diagonal elements of the spin-density matrix, accessible via full angular analyses, decisively test local equilibrium between vorticity and spin, and can signal hydrodynamic spin–orbit non-equilibrium (Moura et al., 2023, Gonçalves et al., 2022).

6. Theoretical Extensions and Current Limitations

Realistic modeling requires integrating the full kinetic evolution, including both dissociation and regeneration, and accounting for the spatial and temporal profile of vorticity and electromagnetic fields in QGP. Statistical freeze-out models predict "maximally impure" matrices with $\rho_{00}=1/3$ and vanishing $\lambda_\phi$, $\lambda_{\theta\phi}$; deviations probe dynamical QCD spin transport (Gonçalves et al., 2022).

Innovative formulations now incorporate both local and long-range quark spin correlations, with effective quark spin correlators mapping directly onto observable hadronic spin alignment, enabling the discrimination of dynamical and induced spin correlations in heavy-ion experiments (Lv et al., 2024). Open-quantum-systems approaches provide a unified, gauge-invariant framework crucial for strong and weak coupling studies (Yang et al., 26 Feb 2025, Yang et al., 2024), yet quantitative progress depends on determining nonperturbative chromomagnetic correlators, e.g., via lattice QCD.

7. Outlook and Significance

Quarkonia spin alignment has become a precision tool for dissecting both the initial production mechanisms and emergent collective dynamics of quark-gluon plasma in high-energy collisions. Measurements of $\rho_{00}(p_T)$ for various states, correlation of spin-alignment observables with kinematic and geometric variables, and mapping of off-diagonal spin-matrix elements stand as central priorities for ongoing and future heavy-ion physics programs.

Spin alignment thus underpins the study of QCD under extreme conditions, enabling direct access to the color-magnetic structure, vorticity, and possible spin-hydrodynamic regime of deconfined matter (Sahoo et al., 21 Dec 2025, Sahoo et al., 11 Jun 2025, Yang et al., 26 Feb 2025). The theoretical advancements, anchored by kinetic equations, open quantum systems, and non-equilibrium hydrodynamics, continue to sharpen predictive control and experimental interpretability in this domain.