Lindblad-driven quarkonium production in heavy-ion collisions

This paper presents a unified open quantum system framework based on the Lindblad equation to describe both the suppression and thermal recombination of charmonium and bottomonium states in ultrarelativistic heavy-ion collisions undergoing Bjorken expansion.

The Gist

Imagine a heavy-ion collision (like smashing two lead atoms together at nearly the speed of light) as a massive, chaotic party where a new, super-hot substance called the "Quark-Gluon Plasma" (QGP) is created. This plasma is like a thick, boiling soup of tiny particles.

Inside this soup, scientists are looking for specific "guests" called quarkonium. These are pairs of heavy particles (like a charm and an anti-charm, or a bottom and an anti-bottom) that usually stick together tightly, like a couple holding hands.

The paper by Armesto and colleagues tries to answer two main questions about these couples in the hot soup:

1. How many get separated (dissociated) by the heat?
2. How many new couples form (recombine) from the loose particles floating around?

Here is the breakdown of their work using simple analogies:

1. The "Complex" Relationship (The Potential)

In a normal vacuum (empty space), these heavy particles have a clear, simple rule for how they stick together. But inside the hot plasma, the rules change. The authors use a mathematical tool called the Gauss law model to describe this new environment.

Think of the force holding the particles together as a rubber band.

• The Real Part: In the hot soup, the rubber band gets weaker and shorter (this is called "Debye screening"). It's like the heat is stretching the band until it's hard to hold on.
• The Imaginary Part: This is the tricky new part. The authors say the rubber band isn't just weak; it's also "jittery." The hot soup constantly bumps into the couple, shaking them apart. This shaking is called "Landau damping."

The paper calculates exactly how much the rubber band weakens and how much it shakes at different temperatures. This allows them to predict the dissociation temperature: the point where the heat is so intense that the couple can no longer hold hands and falls apart.

2. The "Survival" Game (Suppression)

Once the couple is formed, they have to survive the journey through the cooling plasma. The authors use a framework called the Lindblad equation.

Imagine the couple is walking through a crowded, noisy room (the plasma).

• The "Decay Width": This is like a measure of how likely the couple is to get bumped apart by the crowd. The hotter the room, the more likely they are to split up.
• The Survival Probability: The authors calculate the odds that a couple formed at the start of the party will still be holding hands when they reach the exit. If the room gets too hot (above the dissociation temperature), the couple never forms in the first place.

3. The "Reunion" Party (Recombination)

This is the most creative part of the paper. Usually, scientists only counted how many couples survived. But the authors realized that in a crowded room, loose particles might accidentally bump into each other and decide to hold hands, forming a new couple.

• The Analogy: Imagine the party is so crowded that even if the original couples get separated, there are so many loose "hands" (quarks) floating around that they randomly grab each other and form new pairs.
• The Math: They derived a formula to count these "newly formed" couples. However, this only happens if there are enough loose particles.
  • Charm Quarks (J/ψ): There are lots of these in the soup. It's like a room full of people; even if couples break up, new ones form constantly. The authors found that for these, the "reunion" effect is huge and actually helps the total number of couples stay high.
  • Bottom Quarks (Υ): There are very few of these in the soup (about 100 times fewer than charm). It's like a room with only a few people. Even if they break up, the chance of two random people finding each other to hold hands is tiny. So, for these, the "reunion" effect is almost zero.

4. The Results: What They Found

The authors combined the "Survival" (losing couples) and "Reunion" (gaining couples) math to predict what experiments should see.

• For the J/ψ (Charm couple): In the center of the collision (the hottest, densest part), the loss of couples is balanced out by the formation of new ones. The result is that the number of J/ψ particles stays surprisingly high, matching what experiments (ALICE and CMS) have actually measured.
• For the Υ (Bottom couple): Because there are so few bottom quarks, the "reunion" effect is negligible. The number of Υ particles is determined almost entirely by how many survived the heat. Since they are very tightly bound, the strongest ones (Υ(1S)) survive better than the weaker ones, but the overall count is much lower than the J/ψ.

The Big Picture

The paper's main achievement is that it uses one single set of rules (the Lindblad equation) to explain both the breaking up of couples and the forming of new ones.

Previously, scientists might have used one math model to explain why couples break up and a completely different, unrelated model to explain why new ones form. This paper shows that both processes are just two sides of the same coin: the interaction between the heavy particles and the hot plasma. It's a more unified, "first-principles" way of looking at the physics, requiring fewer guesswork parameters to get the right answer.

Technical Summary

Technical Summary: Lindblad-driven quarkonium production in heavy-ion collisions

Problem Statement
Heavy quarkonium states serve as precision probes of the Quark-Gluon Plasma (QGP) formed in ultrarelativistic heavy-ion collisions. Their production is governed by the interplay between Debye screening of the heavy-quark potential and in-medium decoherence, leading to a hierarchical suppression based on binding energy. While the in-medium potential is known to be complex-valued (with a real part describing screening and an imaginary part encoding Landau damping), a unified theoretical framework that simultaneously describes both the suppression of pre-formed states and the recombination (coalescence) of thermalized heavy quarks from first principles has been lacking. Previous approaches often treated suppression and regeneration as independent phenomenological mechanisms.

Methodology
The authors employ an open quantum system framework based on the Lindblad equation to describe the evolution of the heavy quark-antiquark ($Q\bar{Q}$) density matrix coupled to the QGP bath. The methodology proceeds through the following steps:

1. In-Medium Potential: The study utilizes the Gauss law model (Lafferty and Rothkopf) to parametrize the complex in-medium potential $V(r, T) = \text{Re } V + i \text{Im } V$. This model combines the vacuum Cornell potential with Hard-Thermal-Loop (HTL) perturbation theory, governed by a single temperature-dependent parameter, the Debye mass $m_D(T)$. The imaginary part is rigorously regularized to remove nonphysical divergences in the string term.
2. Spectral Functions and Dissociation Properties: Instead of solving the time-dependent Schrödinger equation with a position-dependent imaginary potential (which is numerically challenging), the authors solve the frequency-space Schrödinger equation. They extract in-medium spectral functions by computing the imaginary part of the traced retarded Green's function.
   • Binding Energies ($E_n$): Determined from the peak positions of the spectral functions after subtracting constituent quark masses.
   • Thermal Decay Widths ($\Gamma_n$): Extracted from the full width at half maximum (FWHM) of the peaks, utilizing a skewed Breit-Wigner profile to account for asymmetry near the continuum threshold.
   • Dissociation Temperatures ($T_d$): Identified as the temperature at which the spectral peak merges into the continuum (binding energy vanishes).
3. Survival Probability: For a system undergoing Bjorken expansion, the survival probability of a pre-formed state is calculated by integrating the temperature-dependent decay width $\Gamma(T(t))$ over time.
4. Recombination Model: The framework is extended to include recombination by projecting the stochastic jump operators of the Lindblad equation onto the bound-state subspace. Under an adiabatic approximation (where the effective Hamiltonian evolves slowly compared to quantum jump timescales), a coalescence rate is derived. This rate is proportional to the square of the number of free heavy quarks ($N^2$) and the medium volume ($V$), modulated by a relaxation factor to account for the finite time required for heavy quarks to thermalize.

Key Contributions

• Unified Framework: The paper derives both suppression and recombination mechanisms from a single equation of motion (the Lindblad equation) for the $Q\bar{Q}$ density matrix, ensuring internal consistency between the imaginary part of the potential (driving decay) and the jump operators (driving recombination).
• First-Principles Derivation: The approach avoids ad-hoc combinations of independent models. The complex potential fixes the thermal decay widths, which directly determine both the survival probability and the normalization of the recombination rate via the fluctuation-dissipation theorem.
• Adiabatic Coalescence: The authors derive a specific coalescence model from the Lindblad equation, incorporating a relaxation factor to correct for the non-instantaneous thermalization of heavy quarks in the medium.

Results
The study applies this framework to charmonium ($J/\psi, \psi(2S)$) and bottomonium ($\Upsilon(nS)$) states in Pb–Pb collisions at $\sqrt{s_{NN}} = 5.02$ TeV.

• Dissociation Temperatures:
  • Charmonium: $T_{J/\psi} \approx 372$ MeV; $T_{\psi(2S)} \lesssim 155$ MeV (dissociates near or below the critical temperature).
  • Bottomonium: $T_{\Upsilon(1S)} \approx 564$ MeV; $T_{\Upsilon(2S)} \approx 223$ MeV; $T_{\Upsilon(3S)} \approx 164$ MeV; $T_{\Upsilon(4S)} < 155$ MeV.
• Nuclear Modification Factor ($R_{AA}$):
  • $J/\psi$: The model predicts a significant recombination contribution in central collisions ($N_{part} \gtrsim 200$), driven by the high charm-quark multiplicity at LHC energies. This regeneration partially compensates for suppression, resulting in a total $R_{AA}$ close to unity, consistent with ALICE and CMS data.
  • $\Upsilon(1S)$: Due to the bottom-quark multiplicity being roughly two orders of magnitude smaller than that of charm, the recombination term is negligible. The $R_{AA}$ for $\Upsilon(1S)$ is dominated almost entirely by suppression, with the model results agreeing reasonably well with CMS data without additional fitting parameters.

Significance and Claims
The authors claim that this work provides a "unified, first-principles-inspired description" of quarkonium production. By deriving both suppression and recombination from the same microscopic dynamics (the complex potential and Lindblad jump operators), the framework eliminates the need for independent phenomenological tuning of these competing effects. The paper asserts that once the complex potential is fixed (via the Gauss law model), all medium effects on the quarkonium yield follow without additional free parameters. This represents a significant step toward a truly first-principles description of quarkonium in heavy-ion collisions, offering a more theoretically grounded alternative to traditional transport models where suppression and regeneration are often combined manually.