Lindblad-driven quarkonium production in heavy-ion collisions

Résumé technique : Production de quarkonium pilotée par Lindblad dans les collisions d'ions lourds

Énoncé du problème
Les états de quarkonium lourd servent de sondes de précision du Plasma de Quarks et de Gluons (QGP) formé lors de collisions d'ions lourds ultrarelativistes. Leur production est régie par l'interplay entre l'écranage de Debye du potentiel des quarks lourds et la décohérence en milieu, conduisant à une suppression hiérarchique basée sur l'énergie de liaison. Bien que le potentiel en milieu soit connu pour être à valeurs complexes (avec une partie réelle décrivant l'écranage et une partie imaginaire codant l'amortissement de Landau), un cadre théorique unifié décrivant simultanément, à partir des premiers principes, à la fois la suppression des états préformés et la recombinaison (coalescence) des quarks lourds thermalisés a fait défaut. Les approches précédentes traitaient souvent la suppression et la régénération comme des mécanismes phénoménologiques indépendants.

Méthodologie
Les auteurs emploient un cadre de système quantique ouvert basé sur l'équation de Lindblad pour décrire l'évolution de la matrice de densité quark-antiquark lourd ($Q\bar{Q}$) couplée au bain QGP. La méthodologie procède par les étapes suivantes :

1. Potentiel en milieu : L'étude utilise le modèle de la loi de Gauss (Lafferty et Rothkopf) pour paramétrer le potentiel complexe en milieu $V(r, T) = \text{Re } V + i \text{Im } V$. Ce modèle combine le potentiel Cornell du vide avec la théorie des perturbations Hard-Thermal-Loop (HTL), régie par un seul paramètre dépendant de la température, la masse de Debye $m_D(T)$. La partie imaginaire est rigoureusement régularisée pour éliminer les divergences non physiques dans le terme de la corde.
2. Fonctions spectrales et propriétés de dissociation : Au lieu de résoudre l'équation de Schrödinger dépendante du temps avec un potentiel imaginaire dépendant de la position (ce qui est numériquement difficile), les auteurs résolvent l'équation de Schrödinger dans l'espace des fréquences. Ils extraient les fonctions spectrales en milieu en calculant la partie imaginaire de la fonction de Green retardée tracée.
   • Énergies de liaison ($E_n$) : Déterminées à partir des positions des pics des fonctions spectrales après soustraction des masses des quarks constituants.
   • Largeurs de désintégration thermique ($\Gamma_n$) : Extraites de la largeur à mi-hauteur (FWHM) des pics, utilisant un profil de Breit-Wigner asymétrique pour tenir compte de l'asymétrie près du seuil du continuum.
   • Températures de dissociation ($T_d$) : Identifiées comme la température à laquelle le pic spectral se fond dans le continuum (l'énergie de liaison s'annule).
3. Probabilité de survie : Pour un système subissant une expansion de Bjorken, la probabilité de survie d'un état préformé est calculée en intégrant la largeur de désintégration dépendante de la température $\Gamma(T(t))$ sur le temps.
4. Modèle de recombinaison : Le cadre est étendu pour inclure la recombinaison en projetant les opérateurs de saut stochastiques de l'équation de Lindblad sur le sous-espace des états liés. Sous une approximation adiabatique (où l'hamiltonien effectif évolue lentement par rapport aux échelles de temps des sauts quantiques), un taux de coalescence est dérivé. Ce taux est proportionnel au carré du nombre de quarks lourds libres ($N^2$) et au volume du milieu ($V$), modulé par un facteur de relaxation pour tenir compte du temps fini requis pour que les quarks lourds se thermalisent.

Contributions clés

• Cadre unifié : L'article dérive à la fois les mécanismes de suppression et de recombinaison à partir d'une seule équation du mouvement (l'équation de Lindblad) pour la matrice de densité $Q\bar{Q}$, assurant une cohérence interne entre la partie imaginaire du potentiel (pilotant la désintégration) et les opérateurs de saut (pilotant la recombinaison).
• Dérivation à partir des premiers principes : L'approche évite les combinaisons ad hoc de modèles indépendants. Le potentiel complexe fixe les largeurs de désintégration thermique, qui déterminent directement à la fois la probabilité de survie et la normalisation du taux de recombinaison via le théorème fluctuation-dissipation.
• Coalescence adiabatique : Les auteurs dérivent un modèle de coalescence spécifique à partir de l'équation de Lindblad, incorporant un facteur de relaxation pour corriger la thermalisation non instantanée des quarks lourds dans le milieu.

Résultats
L'étude applique ce cadre aux états de charmonium ($J/\psi, \psi(2S)$) et de bottomonium ($\Upsilon(nS)$) dans les collisions Pb–Pb à $\sqrt{s_{NN}} = 5.02$ TeV.

• Températures de dissociation :
  • Charmonium : $T_{J/\psi} \approx 372$ MeV ; $T_{\psi(2S)} \lesssim 155$ MeV (dissocie près ou en dessous de la température critique).
  • Bottomonium : $T_{\Upsilon(1S)} \approx 564$ MeV ; $T_{\Upsilon(2S)} \approx 223$ MeV ; $T_{\Upsilon(3S)} \approx 164$ MeV ; $T_{\Upsilon(4S)} < 155$ MeV.
• Facteur de modification nucléaire ($R_{AA}$) :
  • $J/\psi$ : Le modèle prédit une contribution de recombinaison significative dans les collisions centrales ($N_{part} \gtrsim 200$), pilotée par la multiplicité élevée de quarks charmés aux énergies du LHC. Cette régénération compense partiellement la suppression, résultant en un $R_{AA}$ total proche de l'unité, cohérent avec les données ALICE et CMS.
  • $\Upsilon(1S)$ : En raison de la multiplicité de quarks bottom étant environ deux ordres de grandeur inférieure à celle des quarks charmés, le terme de recombinaison est négligeable. Le $R_{AA}$ pour $\Upsilon(1S)$ est dominé presque entièrement par la suppression, les résultats du modèle s'accordant raisonnablement bien avec les données CMS sans paramètres d'ajustement supplémentaires.

Signification et affirmations
Les auteurs affirment que ce travail fournit une « description unifiée, inspirée des premiers principes » de la production de quarkonium. En dérivant à la fois la suppression et la recombinaison à partir des mêmes dynamiques microscopiques (le potentiel complexe et les opérateurs de saut de Lindblad), le cadre élimine le besoin d'un ajustement phénoménologique indépendant de ces effets concurrents. L'article affirme que, une fois le potentiel complexe fixé (via le modèle de la loi de Gauss), tous les effets du milieu sur le rendement du quarkonium suivent sans paramètres libres supplémentaires. Cela représente une étape significative vers une description véritablement à partir des premiers principes du quarkonium dans les collisions d'ions lourds, offrant une alternative plus théoriquement fondée aux modèles de transport traditionnels où la suppression et la régénération sont souvent combinées manuellement.