Lindblad-driven quarkonium production in heavy-ion collisions

Technisches Fazit: Lindblad-getriebene Quarkonium-Produktion in Schwerionenkollisionen

Problemstellung
Schwere Quarkonium-Zustände dienen als Präzisionssonden für das in ultrarelativistischen Schwerionenkollisionen gebildete Quark-Gluon-Plasma (QGP). Ihre Produktion wird durch das Zusammenspiel zwischen der Debye-Abschirmung des Schwerquark-Potenzials und der In-Medium-Dekohärenz bestimmt, was zu einer hierarchischen Unterdrückung basierend auf der Bindungsenergie führt. Obwohl bekannt ist, dass das In-Medium-Potenzial komplexwertig ist (mit einem reellen Teil, der die Abschirmung beschreibt, und einem imaginären Teil, der die Landau-Dämpfung kodiert), fehlte bisher ein vereinheitlichtes theoretisches Rahmenwerk, das gleichzeitig die Unterdrückung vorformierter Zustände und die Rekombination (Koaleszenz) thermisierter Schwerquarks aus ersten Prinzipien beschreibt. Frühere Ansätze behandelten Unterdrückung und Regeneration oft als unabhängige phänomenologische Mechanismen.

Methodik
Die Autoren verwenden ein Rahmenwerk offener Quantensysteme, das auf der Lindblad-Gleichung basiert, um die Entwicklung der Dichtematrix des Schwerquark-Antiquark-Paares ($Q\bar{Q}$) gekoppelt an das QGP-Bad zu beschreiben. Die Methodik verläuft in folgenden Schritten:

1. In-Medium-Potenzial: Die Studie nutzt das Gauss-Gesetz-Modell (Lafferty und Rothkopf), um das komplexe In-Medium-Potenzial $V(r, T) = \text{Re } V + i \text{Im } V$ zu parametrisieren. Dieses Modell kombiniert das Vakuum-Cornell-Potenzial mit der Störungstheorie der Hart-Thermal-Loop (HTL) und wird durch einen einzigen temperaturabhängigen Parameter, die Debye-Masse $m_D(T)$, gesteuert. Der imaginäre Teil wird rigoros regularisiert, um nichtphysikalische Divergenzen im String-Term zu entfernen.
2. Spektralfunktionen und Dissoziations-Eigenschaften: Anstatt die zeitabhängige Schrödinger-Gleichung mit einem ortsabhängigen imaginären Potenzial zu lösen (was numerisch herausfordernd ist), lösen die Autoren die Schrödinger-Gleichung im Frequenzraum. Sie extrahieren In-Medium-Spektralfunktionen durch Berechnung des imaginären Teils der getrackten retardierten Green-Funktion.
   • Bindungsenergien ($E_n$): Bestimmt aus den Peak-Positionen der Spektralfunktionen nach Abzug der Konstituentenquark-Massen.
   • Thermische Zerfallsbreiten ($\Gamma_n$): Extrahiert aus der Halbwertsbreite (FWHM) der Peaks unter Verwendung eines schiefen Breit-Wigner-Profils, um die Asymmetrie nahe der Kontinuumsschwelle zu berücksichtigen.
   • Dissoziations-Temperaturen ($T_d$): Identifiziert als die Temperatur, bei der der Spektralpeak in das Kontinuum übergeht (die Bindungsenergie verschwindet).
3. Überlebenswahrscheinlichkeit: Für ein System, das sich einer Bjorken-Expansion unterzieht, wird die Überlebenswahrscheinlichkeit eines vorformierten Zustands durch Integration der temperaturabhängigen Zerfallsbreite $\Gamma(T(t))$ über die Zeit berechnet.
4. Rekombinations-Modell: Das Rahmenwerk wird erweitert, um Rekombination einzubeziehen, indem die stochastischen Sprungoperatoren der Lindblad-Gleichung auf den gebundenen Zustandsraum projiziert werden. Unter einer adiabatischen Näherung (bei der sich der effektive Hamilton-Operator langsam im Vergleich zu den Zeitskalen quantenmechanischer Sprünge entwickelt) wird eine Koaleszenzrate abgeleitet. Diese Rate ist proportional zum Quadrat der Anzahl freier Schwerquarks ($N^2$) und zum Mediumvolumen ($V$), moduliert durch einen Relaxationsfaktor, der die endliche Zeit berücksichtigt, die Schwerquarks zur Thermalisierung benötigen.

Hauptbeiträge

• Vereinheitlichtes Rahmenwerk: Die Arbeit leitet sowohl Unterdrückungs- als auch Rekombinationsmechanismen aus einer einzigen Bewegungsgleichung (der Lindblad-Gleichung) für die $Q\bar{Q}$-Dichtematrix ab und gewährleistet so eine innere Konsistenz zwischen dem imaginären Teil des Potenzials (der den Zerfall antreibt) und den Sprungoperatoren (die die Rekombination antreiben).
• Ableitung aus ersten Prinzipien: Der Ansatz vermeidet ad-hoc-Kombinationen unabhängiger Modelle. Das komplexe Potenzial fixiert die thermischen Zerfallsbreiten, die direkt sowohl die Überlebenswahrscheinlichkeit als auch die Normierung der Rekombinationsrate über das Fluktuations-Dissipations-Theorem bestimmen.
• Adiabatische Koaleszenz: Die Autoren leiten ein spezifisches Koaleszenzmodell aus der Lindblad-Gleichung ab, das einen Relaxationsfaktor integriert, um die nicht-instantane Thermalisierung von Schwerquarks im Medium zu korrigieren.

Ergebnisse
Die Studie wendet dieses Rahmenwerk auf Charmonium ($J/\psi, \psi(2S)$) und Bottomonium ($\Upsilon(nS)$) in Pb–Pb-Kollisionen bei $\sqrt{s_{NN}} = 5.02$ TeV an.

• Dissoziations-Temperaturen:
  • Charmonium: $T_{J/\psi} \approx 372$ MeV; $T_{\psi(2S)} \lesssim 155$ MeV (dissoziiert nahe oder unterhalb der kritischen Temperatur).
  • Bottomonium: $T_{\Upsilon(1S)} \approx 564$ MeV; $T_{\Upsilon(2S)} \approx 223$ MeV; $T_{\Upsilon(3S)} \approx 164$ MeV; $T_{\Upsilon(4S)} < 155$ MeV.
• Nuklearer Modifikationsfaktor ($R_{AA}$):
  • $J/\psi$: Das Modell sagt einen signifikanten Rekombinationsbeitrag in zentralen Kollisionen ($N_{part} \gtrsim 200$) voraus, getrieben durch die hohe Charm-Quark-Multiplizität bei LHC-Energien. Diese Regeneration kompensiert die Unterdrückung teilweise, was zu einem Gesamt-$R_{AA}$ nahe Eins führt, konsistent mit ALICE- und CMS-Daten.
  • $\Upsilon(1S)$: Da die Bottom-Quark-Multiplizität etwa zwei Größenordnungen kleiner ist als die der Charm-Quarks, ist der Rekombinationsterm vernachlässigbar. Der $R_{AA}$ für $\Upsilon(1S)$ wird fast ausschließlich durch Unterdrückung dominiert; die Modellergebnisse stimmen ohne zusätzliche Anpassungsparameter vernünftig gut mit den CMS-Daten überein.

Bedeutung und Behauptungen
Die Autoren behaupten, dass diese Arbeit eine „vereinheitlichte, aus ersten Prinzipien inspirierte Beschreibung" der Quarkonium-Produktion liefert. Durch die Ableitung sowohl von Unterdrückung als auch von Rekombination aus derselben mikroskopischen Dynamik (dem komplexen Potenzial und den Lindblad-Sprungoperatoren) eliminiert das Rahmenwerk die Notwendigkeit einer unabhängigen phänomenologischen Anpassung dieser konkurrierenden Effekte. Die Arbeit stellt fest, dass einmal das komplexe Potenzial fixiert ist (via des Gauss-Gesetz-Modells), alle Medium-Effekte auf die Quarkonium-Ausbeute ohne zusätzliche freie Parameter folgen. Dies stellt einen bedeutenden Schritt hin zu einer wirklich aus ersten Prinzipien abgeleiteten Beschreibung von Quarkonium in Schwerionenkollisionen dar und bietet eine theoretisch fundiertere Alternative zu traditionellen Transportmodellen, bei denen Unterdrückung und Regeneration oft manuell kombiniert werden.