Lindblad-driven quarkonium production in heavy-ion collisions

Technische Samenvatting: Lindblad-gedreven quarkoniumproductie in zware-ionenbotsingen

Probleemstelling
Zware quarkoniumtoestanden fungeren als precisieprobes van het Quark-Gluon Plasma (QGP) dat wordt gevormd in ultrarelativistische zware-ionenbotsingen. Hun productie wordt bepaald door het samenspel tussen Debye-scherming van het zware-quarkpotentiaal en decoherentie in het medium, wat leidt tot een hiërarchische onderdrukking gebaseerd op bindingsenergie. Hoewel bekend is dat het potentiaal in het medium complexwaardig is (met een reëel deel dat scherming beschrijft en een imaginair deel dat Landau-demping codeert), ontbrak er tot nu toe een verenigd theoretisch kader dat zowel de onderdrukking van vooraf gevormde toestanden als de recombinatie (coalescentie) van gethermaliseerde zware quarks vanuit eerste principes tegelijkertijd beschrijft. Eerdere benaderingen behandelden onderdrukking en regeneratie vaak als onafhankelijke fenomenologische mechanismen.

Methodologie
De auteurs maken gebruik van een kader voor open kwantumsystemen gebaseerd op de Lindblad-vergelijking om de evolutie van de dichtheidsmatrix van het zware quark-antiquark-paar ($Q\bar{Q}$) gekoppeld aan het QGP-bad te beschrijven. De methodologie verloopt via de volgende stappen:

1. Potentiaal in het medium: De studie maakt gebruik van het Gauss-wetmodel (Lafferty en Rothkopf) om de complexe potentiaal in het medium $V(r, T) = \text{Re } V + i \text{Im } V$ te parametriseren. Dit model combineert het vacuüm-Cornellpotentiaal met Hard-Thermal-Loop (HTL) verstoringstheorie, die wordt bestuurd door één temperatuurafhankelijke parameter, de Debye-massa $m_D(T)$. Het imaginaire deel wordt strikt geregulariseerd om niet-fysische divergenties in het snaarterm te verwijderen.
2. Spectrale functies en dissociatie-eigenschappen: In plaats van de tijdsafhankelijke Schrödingervergelijking op te lossen met een positie-afhankelijk imaginair potentiaal (wat numeriek uitdagend is), lossen de auteurs de Schrödingervergelijking in frequentieruimte op. Ze extraheren spectrale functies in het medium door het imaginaire deel van de getraceerde vertraagde Green-functie te berekenen.
   • Bindingsenergieën ($E_n$): Bepaald uit de piekposities van de spectrale functies na aftrek van de massa's van de samenstellende quarks.
   • Thermische vervalbreedten ($\Gamma_n$): Geëxtraheerd uit de volle breedte op halve maximum (FWHM) van de pieken, waarbij een scheefgetrokken Breit-Wigner-profiel wordt gebruikt om rekening te houden met asymmetrie nabij de continuüm-drempel.
   • Dissociatietemperaturen ($T_d$): Geïdentificeerd als de temperatuur waarbij de spectrale piek samensmelt met het continuüm (bindingsenergie verdwijnt).
3. Overlevingskans: Voor een systeem dat een Bjorken-expansie ondergaat, wordt de overlevingskans van een vooraf gevormde toestand berekend door de temperatuurafhankelijke vervalbreedte $\Gamma(T(t))$ over de tijd te integreren.
4. Recombinatiemodel: Het kader wordt uitgebreid om recombinatie op te nemen door de stochastische jump-operatoren van de Lindblad-vergelijking te projecteren op de gebonden-toestandsruimte. Onder een adiabatische benadering (waarbij de effectieve Hamiltoniaan langzaam evolueert ten opzichte van kwantumjump-tijdschalen) wordt een coalescentiesnelheid afgeleid. Deze snelheid is evenredig met het kwadraat van het aantal vrije zware quarks ($N^2$) en het mediumvolume ($V$), gemoduleerd door een relaxatiefactor om rekening te houden met de eindige tijd die nodig is voor zware quarks om te thermaliseren.

Belangrijkste bijdragen

• Verenigd kader: Het artikel leidt zowel onderdrukkings- als recombinatiemechanismen af uit één bewegingsvergelijking (de Lindblad-vergelijking) voor de $Q\bar{Q}$-dichtheidsmatrix, waardoor interne consistentie wordt gewaarborgd tussen het imaginaire deel van het potentiaal (dat verval aandrijft) en de jump-operatoren (die recombinatie aandrijven).
• Afleiding vanuit eerste principes: De aanpak vermijdt ad-hoc combinaties van onafhankelijke modellen. Het complexe potentiaal bepaalt de thermische vervalbreedten, die direct zowel de overlevingskans als de normalisatie van de recombinatiesnelheid bepalen via de fluctuatie-dissipatietheorema.
• Adiabatische coalescentie: De auteurs leiden een specifiek coalescentiemodel af uit de Lindblad-vergelijking, waarbij een relaxatiefactor wordt opgenomen om de niet-instantane thermalisatie van zware quarks in het medium te corrigeren.

Resultaten
De studie past dit kader toe op charmonium ($J/\psi, \psi(2S)$) en bottomonium ($\Upsilon(nS)$) toestanden in Pb–Pb-botsingen bij $\sqrt{s_{NN}} = 5.02$ TeV.

• Dissociatietemperaturen:
  • Charmonium: $T_{J/\psi} \approx 372$ MeV; $T_{\psi(2S)} \lesssim 155$ MeV (dissocieert nabij of onder de kritische temperatuur).
  • Bottomonium: $T_{\Upsilon(1S)} \approx 564$ MeV; $T_{\Upsilon(2S)} \approx 223$ MeV; $T_{\Upsilon(3S)} \approx 164$ MeV; $T_{\Upsilon(4S)} < 155$ MeV.
• Nucleaire modificatiefactor ($R_{AA}$):
  • $J/\psi$: Het model voorspelt een aanzienlijke recombinatiebijdrage in centrale botsingen ($N_{part} \gtrsim 200$), gedreven door de hoge multipliciteit van charm-quarks bij LHC-energieën. Deze regeneratie compenseert gedeeltelijk de onderdrukking, wat resulteert in een totale $R_{AA}$ dicht bij één, consistent met ALICE- en CMS-data.
  • $\Upsilon(1S)$: Omdat de multipliciteit van bottom-quarks ongeveer twee ordes van grootte kleiner is dan die van charm, is de recombinatieterm verwaarloosbaar. De $R_{AA}$ voor $\Upsilon(1S)$ wordt bijna uitsluitend gedomineerd door onderdrukking, waarbij de modelresultaten redelijk goed overeenkomen met CMS-data zonder extra aanpassingsparameters.

Betekenis en claims
De auteurs claimen dat dit werk een "verenigde, vanuit eerste principes geïnspireerde beschrijving" van quarkoniumproductie biedt. Door zowel onderdrukking als recombinatie af te leiden uit dezelfde microscopische dynamica (het complexe potentiaal en Lindblad-jump-operatoren), elimineert het kader de noodzaak voor onafhankelijke fenomenologische afstemming van deze concurrerende effecten. Het artikel stelt dat zodra het complexe potentiaal is vastgesteld (via het Gauss-wetmodel), alle mediumeffecten op de quarkoniumopbrengst volgen zonder extra vrije parameters. Dit vormt een aanzienlijke stap naar een echt eerste-principes-beschrijving van quarkonium in zware-ionenbotsingen, en biedt een meer theoretisch onderbouwde alternatief voor traditionele transportmodellen waarbij onderdrukking en regeneratie vaak handmatig worden gecombineerd.