Quarkonium $p_{\rm T}$ spectra in heavy--ion collisions at LHC energies within a hydrodynamic core--corona framework

This paper presents a unified analytical hydrodynamic core-corona framework that successfully describes the transverse momentum spectra and yield ratios of charmonium and bottomonium states in Pb-Pb collisions at LHC energies, demonstrating its effectiveness in modeling quarkonium production across a broad kinematic range.

The Gist

Imagine smashing two giant, heavy balls (lead nuclei) together at nearly the speed of light. When they collide, they create a tiny, super-hot, super-dense soup of particles called a "quark-gluon plasma" (QGP). This soup is so hot that the usual rules of physics change; the particles that normally stick together to form atoms (like protons and neutrons) melt into a free-flowing fluid.

The scientists in this paper are trying to understand how specific "heavy" particles, called quarkonia, behave inside this soup. Think of quarkonia as heavy-duty couples: a heavy quark and its anti-quark partner holding hands. In normal conditions, they stay together. But in this hot soup, the heat tries to pull them apart.

Here is a simple breakdown of what the researchers did and found:

1. The Two-Part Model: The "Core" and the "Corona"

To explain how these heavy couples survive the crash, the authors used a clever two-part recipe, like a core-and-crust model for a pizza or a core-and-corona for a star.

• The Core (The Hot Soup): This is the middle of the collision where the density is highest. Here, the soup is so thick and hot that it acts like a fluid. The researchers used a mathematical "hydrodynamic" framework (think of it like a weather model for fluids) to describe how this soup expands and cools down. They assumed the soup expands like a balloon being blown up, but in a specific, symmetrical way.
• The Corona (The Outer Edge): Not every part of the collision is a perfect fluid. On the very edges, the density is lower, like the thin outer crust of a pizza. Here, the particles don't melt into a soup; they just bounce off each other like billiard balls. The researchers modeled this part using data from simpler collisions (proton-on-proton) to represent these "hard" interactions.

By combining the fluid-like Core and the billiard-ball-like Corona, they created a complete picture of what happens to the heavy particles.

2. The Experiment: Catching the Particles

The team looked at data from the Large Hadron Collider (LHC), specifically from collisions of lead nuclei. They focused on two types of heavy couples:

• Charmonium (J/ψ and ψ(2S)): Made of "charm" quarks. These are like lighter couples in the heavy-quark world.
• Bottomonium (ϒ(1S), ϒ(2S), ϒ(3S)): Made of "bottom" quarks. These are much heavier and tighter couples.

They measured how much "sideways" energy (transverse momentum, or $p_T$) these particles had when they finally escaped the collision.

3. The Results: Different Couples, Different Stories

The paper found that these two types of couples tell different stories about the soup:

• The Bottomonium Story (The Early Bird):
  The heavy bottom couples are so tightly bound that they can survive the hottest, earliest moments of the collision. The model showed they "freeze out" (stop interacting with the soup) at a very high temperature (around 224 MeV) and don't get pushed around by the fluid flow as much.

  • The Analogy: Imagine a heavy anchor dropped in a river. It sinks quickly and stays put, feeling the current only for a short time. The bottom couples are like that anchor; they probe the very hottest, earliest stage of the soup.
  • The Pattern: The model successfully predicted that the looser bottom couples (like ϒ(2S) and ϒ(3S)) get melted more easily than the tightest one (ϒ(1S)). This is called "sequential suppression," and the model got it right.

• The Charmonium Story (The Latecomer):
  The charm couples are lighter and looser. They seem to survive longer and get swept up by the expanding fluid flow more than the bottom couples. They "freeze out" at a lower temperature (around 160 MeV) and have more sideways push.

  • The Analogy: Imagine a leaf floating on that same river. It gets carried along by the current for a long time, feeling the flow of the water. The charm couples are like that leaf; they interact with the soup for a longer time and are more influenced by its movement.
  • The Twist: The model worked great for low and medium speeds, but at very high speeds, it slightly underestimated the number of particles. This suggests there are some other "hard" mechanisms (like high-energy collisions) happening that the fluid model doesn't fully capture yet.

4. The Big Picture

The main takeaway is that this Core-Corona approach, combined with a fluid dynamics model, works very well to explain the data.

• It successfully describes how the heavy particles move and how many of them survive.
• It confirms that bottomonium acts as a thermometer for the very early, hottest moments of the collision.
• It confirms that charmonium is more influenced by the later stages of the collision, where the fluid flow is stronger.

In short, the paper shows that by treating the collision as a mix of a hot, expanding fluid (the core) and some leftover hard collisions (the corona), scientists can get a clear, unified view of how heavy particles behave in the extreme conditions created at the LHC.

Technical Summary

1. Problem Statement

The production and suppression of quarkonium states (charmonium $c\bar{c}$ and bottomonium $b\bar{b}$) in relativistic heavy-ion collisions are primary probes for studying the Quark-Gluon Plasma (QGP). While the nuclear modification factor ($R_{AA}$) is commonly used to quantify suppression, it does not fully capture the thermal properties or collective behavior of the medium.

• The Challenge: Existing models often struggle to provide a unified description of transverse momentum ($p_T$) spectra across the entire kinematic range. Pure hydrodynamic models describe the low-$p_T$ thermalized region well but fail to reproduce the high-$p_T$ tail dominated by hard scattering. Conversely, pure perturbative QCD approaches miss the collective flow effects observed at low $p_T$.
• Specific Questions:
  • How do charmonium and bottomonium states decouple from the medium at different stages?
  • Can an analytical hydrodynamic framework, combined with a non-thermal component, simultaneously describe the spectra and yield ratios of $J/\psi$, $\psi(2S)$, and $\Upsilon(nS)$ states at LHC energies ($\sqrt{s_{NN}} = 5.02$ TeV)?
  • What are the distinct freeze-out parameters (temperature and flow) for charm vs. bottom quarks?

2. Methodology

The authors employ a hybrid theoretical framework combining analytical relativistic hydrodynamics with a core–corona approach.

A. Hydrodynamic Core (Thermal Component)

• Framework: The medium evolution is modeled using an exact analytical solution of relativistic hydrodynamics (based on Ref. [55]).
• Assumptions:
  • Cylindrical symmetry.
  • Boost-invariant longitudinal expansion.
  • Hubble-like transverse flow.
  • Ideal hydrodynamics with a constant speed of sound ($c_s$).
• Particle Emission: Quarkonium spectra are calculated on a constant-temperature freeze-out hypersurface using the Cooper–Frye prescription. This yields a closed-form analytical expression for the transverse momentum distribution (Eq. 1 in the paper), dependent on the freeze-out temperature ($T$) and the expansion parameter ($\nu$, related to average transverse velocity $\langle v_T \rangle$).

B. Core–Corona Approach (Non-Thermal Component)

• Rationale: To account for high-$p_T$ particles produced in initial hard scatterings that do not thermalize.
• Definition:
  • Core: The high-density central region where thermalization occurs and collective expansion is dominant.
  • Corona: The dilute peripheral region (where local density drops to $\approx 10\%$ of the central density) where interactions resemble proton-proton ($pp$) collisions.
• Implementation:
  • The corona contribution is modeled using a power-law functional form fitted to $pp$ collision data at $\sqrt{s} = 5.02$ TeV.
  • The relative fractions of core and corona contributions are determined using a Glauber Monte Carlo approach based on the nuclear charge density distribution.
  • The total spectrum is the sum of the hydrodynamic core and the scaled corona contribution.

C. Data Comparison

The model is tested against experimental data from:

• ALICE: $J/\psi$ and $\psi(2S)$ spectra at mid-rapidity ($|y| < 0.9$) and forward rapidity ($2.5 < y < 4$).
• CMS: $\Upsilon(1S)$, $\Upsilon(2S)$, and $\Upsilon(3S)$ spectra at mid-rapidity ($|y| < 2.4$).
• Kinematics: Pb–Pb collisions at $\sqrt{s_{NN}} = 5.02$ TeV.

3. Key Contributions

1. Unified Analytical Framework: The paper demonstrates that an analytical hydrodynamic solution, often considered a simplification, can effectively describe complex heavy-ion data when combined with a core–corona mechanism.
2. Distinct Freeze-out Scenarios: The study quantitatively establishes a clear distinction between the freeze-out conditions of charmonium and bottomonium:
   • Charmonium: Freezes out at a lower temperature ($T \approx 160$ MeV) with higher transverse flow ($\langle v_T \rangle \approx 0.60$), indicating sensitivity to later stages of the medium evolution and potential regeneration.
   • Bottomonium: Freezes out at a significantly higher temperature ($T \approx 224$ MeV) with lower transverse flow ($\langle v_T \rangle \approx 0.56$), confirming that bottomonium probes the earlier, hotter stages of the QGP where color screening dominates and regeneration is negligible.
3. Sequential Suppression Validation: The model successfully reproduces the sequential suppression pattern of $\Upsilon(nS)$ states, consistent with their binding energies.

4. Results

• Charmonium ($J/\psi, \psi(2S)$):
  • The model provides an excellent description of the $p_T$ spectra for $p_T < 8$ GeV.
  • The $\psi(2S)/J/\psi$ yield ratio is well reproduced across the measured $p_T$ range, validating the model's ability to capture the interplay between suppression and regeneration.
  • Deviation: At high $p_T$ ($> 8$ GeV), the model slightly underestimates the $J/\psi$ yield. This suggests the presence of additional hard production mechanisms (e.g., gluon fragmentation, viscous corrections) not included in the ideal hydrodynamic core.
• Bottomonium ($\Upsilon(nS)$):
  • The spectra for $\Upsilon(1S)$, $\Upsilon(2S)$, and $\Upsilon(3S)$ are well described up to $p_T \approx 15$ GeV.
  • The yield ratios $\Upsilon(2S)/\Upsilon(1S)$ and $\Upsilon(3S)/\Upsilon(1S)$ match the experimental data, correctly reflecting the expected sequential melting pattern (loosely bound states melt more easily).
  • The results are consistent with the assumption that bottomonium production is governed almost entirely by suppression, with negligible regeneration.

5. Significance

• Physical Insight: The work confirms that heavy quarkonia decouple from the medium at different temperatures, providing a "thermometer" for different stages of the QGP evolution. Bottomonium acts as a probe for the initial hot phase, while charmonium reflects the later, cooler, and more collective phase.
• Methodological Value: It validates the use of analytical hydrodynamic solutions as a transparent and efficient benchmark for more complex numerical simulations. The core–corona approach successfully bridges the gap between soft (thermal) and hard (perturbative) physics without requiring full transport model complexity.
• Future Directions: The deviations observed at high $p_T$ highlight the necessity of incorporating viscous corrections and microscopic parton fragmentation mechanisms in future refinements of the model to achieve a complete description of the full $p_T$ spectrum.

In conclusion, the paper presents a robust, unified description of quarkonium production in LHC heavy-ion collisions, demonstrating that the combination of analytical hydrodynamics and a core–corona framework effectively captures the essential thermal and non-thermal dynamics of the QGP.