Hadronic $J/ψ$ Regeneration in Pb+Pb Collisions

This study demonstrates that hadronic regeneration of $J/\psi$ mesons from $D$-meson collisions during the hadronic break-up phase in Pb+Pb collisions significantly contributes to the final observed yields, implying that regeneration must be accounted for in models and making it difficult to distinguish between regeneration occurring at hadronization versus during final-state interactions.

The Gist

Imagine a massive, high-speed collision between two lead nuclei (like smashing two giant, dense clouds of particles together). This happens at the Large Hadron Collider (LHC), creating a tiny, super-hot "soup" called the Quark-Gluon Plasma (QGP).

In this soup, particles called charm quarks are created in huge numbers. As the soup cools down, it "freezes" into ordinary matter (hadrons). One of the most interesting things that can form from these charm quarks is a particle called $J/\psi$ (pronounced "J-psi").

For a long time, scientists thought that almost all the $J/\psi$ particles they saw in the final results were formed immediately when the soup froze. However, the numbers they measured were surprisingly high. This led to a big question: Did the $J/\psi$ particles form right at the moment of freezing, or did they form later, from collisions between other particles as the system cooled down?

This paper investigates that second possibility: Regeneration.

The Analogy: The "Rebuilding" Workshop

Think of the collision aftermath like a chaotic construction site after a storm.

1. The Initial Freeze (Hadronization): When the hot soup cools, it instantly builds a bunch of "bricks" (D-mesons) and a few "houses" ($J/\psi$).
2. The Expansion: The construction site expands and cools down. The bricks are flying around, bumping into each other.
3. The Regeneration Process: Even after the initial "houses" are built, the flying bricks (D-mesons) keep crashing into each other. Sometimes, when two bricks hit just right, they smash together and rebuild a new house ($J/\psi$).

The authors of this paper wanted to know: How many of the final houses were built during the initial freeze, and how many were rebuilt later by the flying bricks?

How They Did the Math

The scientists used a "recipe" (a mathematical model) to calculate how often these bricks collide and rebuild houses.

• The Ingredients: They used real data from the LHC about how many "bricks" (D-mesons) were created.
• The Rules: They used physics rules (cross-sections) that tell us how likely it is for two bricks to smash together and make a house.
• The Simulation: They imagined the construction site expanding over time, calculating how many new houses get built as the bricks keep flying around and colliding.

The Big Discovery

The results were surprising and very important:

1. It's a Mix: The final number of $J/\psi$ particles we see is a mix of the ones built at the start and the ones rebuilt later.
2. The "Rebuild" is Huge: They found that the "rebuilding" process (regeneration) is so powerful that it could account for everything we see.
   • In their calculations, the final yield of $J/\psi$ could be explained if 0% to 100% of them were built at the very start.
   • Conversely, it could also mean that 28% to 113% of the final yield was built at the start, with the rest added later.

The "Invisible Hand" Problem

Here is the tricky part, explained with a metaphor:

Imagine you walk into a room and see 100 chairs. You want to know how many were there when you entered, and how many were assembled by people working in the room while you were watching.

The authors found that the "workers" (the D-meson collisions) are so efficient at building chairs that you cannot tell from the final count of 100 chairs how many were there at the start.

• If you started with 10 chairs, the workers could build 90 more.
• If you started with 90 chairs, the workers could build 10 more.
• The result looks exactly the same either way.

Why This Matters

This paper tells us that we cannot simply look at the final number of $J/\psi$ particles and say, "Aha! This proves the Quark-Gluon Plasma was saturated with charm quarks at the moment of freezing."

Why? Because the "rebuilding" phase (hadronic regeneration) is so effective that it masks the initial conditions. It's like trying to guess how much rain fell during a storm by looking at a puddle, but realizing that a giant hose was also spraying water into the puddle the whole time.

The Bottom Line

• Regeneration is real: The collisions between D-mesons do create a significant number of $J/\psi$ particles after the initial freeze.
• We can't be sure: Because of this regeneration, we cannot definitively say how many $J/\psi$ particles existed the instant the plasma turned into matter.
• Future models must change: Any computer simulation of these collisions must include this "rebuilding" phase, or the results will be wrong.

In short: The universe is a busy workshop where things are constantly being built and rebuilt. To understand the final product, we have to account for the work happening after the initial blueprint is drawn.

Technical Summary

1. Problem Statement

The paper addresses a critical ambiguity in interpreting the surprisingly high yield of $J/\psi$ mesons measured in Lead-Lead (Pb+Pb) collisions at the Large Hadron Collider (LHC).

• Context: At LHC energies ($\sqrt{s_{NN}} = 2.76$ and $5.02$ TeV), the initial collision produces an abundance of charm-anticharm ($c\bar{c}$) pairs, leading to an overpopulation of charm quarks in the Quark-Gluon Plasma (QGP) characterized by a fugacity factor $\lambda_c \approx 30$.
• The Issue: The final observed $J/\psi$ yield is the result of four potential mechanisms:
  1. Primordial formation (gluon-gluon fusion).
  2. Recombination within the QGP.
  3. Recombination at the QGP-hadron phase boundary (statistical hadronization).
  4. Regeneration via hadronic interactions in the late-stage hadron gas (specifically $D + \bar{D} \to J/\psi + h$).
• The Gap: Previous studies on hadronic regeneration predated the availability of precise LHC data for $D$-meson yields. It was unclear whether the measured $J/\psi$ yield could be entirely explained by regeneration from $D$-mesons in the hadronic phase, thereby "cloaking" the contributions from earlier QGP mechanisms. If regeneration alone can establish chemical equilibrium, the initial conditions at hadronization become impossible to deduce from final-state data alone.

2. Methodology

The authors employ a kinetic rate equation approach combined with published experimental data and schematic expansion models.

• Reaction Model:

  • They utilize a broken $SU(4)$ flavor symmetry effective interaction model (based on Abreu et al.) to describe interactions between pseudoscalar ($P$) and vector ($V$) mesons.
  • Cross Sections: They calculate thermal, isospin-averaged cross-sections for $J/\psi$ production processes (e.g., $D\bar{D} \to J/\psi + \pi$, $D^*\bar{D} \to J/\psi + \rho$). The thermal averages $\langle \sigma v \rangle$ are found to be nearly temperature-independent in the relevant range ($90 \text{ MeV} \leq T \leq 155 \text{ MeV}$).
  • Key Processes: The study focuses on reactions involving open charm mesons ($D, D^*, D_s, D_s^*$) and light hadrons ($\pi, \rho, K, K^*$).

• Rate Equation:

  • The time evolution of the $J/\psi$ number ($N_{J/\psi}$) is governed by a balance between formation and absorption:
    $$\frac{dN_{J/\psi}}{d\tau} = C_{J/\psi} - A_{J/\psi}(\tau) \frac{N_{J/\psi}(\tau)}{V(\tau)}$$
  • $C_{J/\psi}$ represents the formation rate (driven by $D$-meson abundances).
  • $A_{J/\psi}$ represents the absorption rate (driven by light hadron densities).
  • The equation is solved for volume-integrated numbers per unit rapidity.

• Initial Conditions & Inputs:

  • $D$-meson Yields: They use measured yields for $D^0, D^+, D_s$ and their excited states from ALICE (0-10% central Pb+Pb collisions at 5.02 TeV). They account for feed-down from short-lived states.
  • Expansion Models: Instead of a full transport simulation, they test schematic expansion models for the hadronic fireball:
    1. Abreu et al. Model: A cylindrical fireball with a specific transverse radius evolution.
    2. Andronic et al. / Mazeliauskas et al. Models: Models based on freeze-out contours with transverse flow profiles ($\beta(r)$).
  • Time Evolution: The integration is performed from the hadronization time ($\tau_H \approx 10\text{--}16$ fm/c) to $\tau \to \infty$. The authors argue that stopping at the nominal kinetic freeze-out ($\sim 20$ fm/c) underestimates the yield because $D$-meson collisions remain effective due to the exothermic nature of the reactions.

3. Key Contributions

• Updated Calculation: This is the first study to incorporate the latest LHC $D$-meson yield data into the calculation of hadronic $J/\psi$ regeneration.
• Quantification of Regeneration: The study quantifies the maximum possible contribution of the hadronic phase to the final $J/\psi$ yield.
• Robustness Analysis: The authors demonstrate that their conclusions hold regardless of specific assumptions regarding the expansion profile (isentropic vs. non-isentropic) or exact cross-section magnitudes (within a factor of 2).

4. Results

• Regeneration Capacity: Starting with zero $J/\psi$ at the moment of hadronization ($N_{J/\psi}(\tau_H) = 0$), the hadronic regeneration mechanism alone produces a final yield of $0.047 < N_{J/\psi} < 0.073$.
• Comparison to Equilibrium: This regenerated yield constitutes roughly 40% to 60% of the measured chemical equilibrium yield ($N_{eq} \approx 0.12$).
• Bounding Initial Abundance: By varying the initial $J/\psi$ abundance at hadronization, the authors determine the range of initial conditions consistent with the final measured yield. They find that the fractional abundance of $J/\psi$ immediately after hadronization ($f_{init}$) relative to the equilibrium yield ($N_{eq}$) must satisfy:
  $$0.28 \leq \frac{dN_{J/\psi}^0/dy}{dN_{J/\psi}^{eq}/dy} \leq 1.13$$
• Implication of Bounds: This range implies that the initial $J/\psi$ yield could be as low as 28% of the equilibrium value (meaning ~72% was regenerated) or as high as 113% (meaning regeneration contributed little or the system was slightly over-saturated).
• Sensitivity: The results are robust against variations in the expansion model parameters ($n, \beta_H$) and the temperature dependence of cross-sections.

5. Significance and Conclusion

• Indistinguishability: The study concludes that it is principally impossible to distinguish the contribution of early-stage mechanisms (QGP recombination) from late-stage hadronic regeneration using only the final $J/\psi$ yield data. If regeneration is efficient enough to drive the system toward equilibrium, the "memory" of the initial production mechanism is lost.
• Modeling Requirement: Any microscopic transport model of heavy-ion collisions at the LHC must include hadronic regeneration via $D$-meson collisions to accurately describe $J/\psi$ yields. Ignoring this phase leads to a significant underestimation of the total yield.
• Underestimation Warning: The authors note their calculation likely underestimates the true regeneration yield because they neglected certain channels (e.g., $D_s + \bar{D}_s \to J/\psi + \phi$) and feed-down from excited charmonium states ($\chi_c$). If these are included, the regeneration contribution could be even larger, potentially allowing for 100% regeneration of the observed yield.
• Final Takeaway: The observed equilibrium-like yield of $J/\psi$ does not provide unambiguous evidence that all observed $J/\psi$ mesons were formed during the hadronization transition; a substantial portion (or potentially all) could be generated by final-state hadronic interactions.