Lindblad-driven quarkonium production in heavy-ion collisions

Riepilogo Tecnico: Produzione di quarkonio guidata da Lindblad nelle collisioni di ioni pesanti

Enunciato del Problema
Gli stati di quarkonio pesante fungono da sonde di precisione del Plasma di Quark e Gluoni (QGP) formato nelle collisioni di ioni pesanti ultrarelativistiche. La loro produzione è governata dall'interplay tra lo schermo di Debye del potenziale dei quark pesanti e la decoerenza in mezzo, portando a una soppressione gerarchica basata sull'energia di legame. Sebbene il potenziale in mezzo sia noto per essere a valori complessi (con una parte reale che descrive lo schermo e una parte immaginaria che codifica lo smorzamento di Landau), è mancato un quadro teorico unificato che descriva simultaneamente, a partire dai primi principi, sia la soppressione degli stati preformati sia la ricombinazione (coalescenza) dei quark pesanti termalizzati. I precedenti approcci hanno spesso trattato soppressione e rigenerazione come meccanismi fenomenologici indipendenti.

Metodologia
Gli autori impiegano un quadro di sistema quantistico aperto basato sull'equazione di Lindblad per descrivere l'evoluzione della matrice di densità quark-antiquark pesante ($Q\bar{Q}$) accoppiata al bagno del QGP. La metodologia procede attraverso i seguenti passaggi:

1. Potenziale in Mezzo: Lo studio utilizza il modello della legge di Gauss (Lafferty e Rothkopf) per parametrizzare il potenziale complesso in mezzo $V(r, T) = \text{Re } V + i \text{Im } V$. Questo modello combina il potenziale Cornell nel vuoto con la teoria delle perturbazioni Hard-Thermal-Loop (HTL), governato da un singolo parametro dipendente dalla temperatura, la massa di Debye $m_D(T)$. La parte immaginaria è rigorosamente regolarizzata per rimuovere divergenze non fisiche nel termine della stringa.
2. Funzioni Spettrali e Proprietà di Dissociazione: Invece di risolvere l'equazione di Schrödinger dipendente dal tempo con un potenziale immaginario dipendente dalla posizione (il che è numericamente impegnativo), gli autori risolvono l'equazione di Schrödinger nello spazio delle frequenze. Estraggono le funzioni spettrali in mezzo calcolando la parte immaginaria della funzione di Green ritardata tracciata.
   • Energie di Legame ($E_n$): Determinate dalle posizioni dei picchi delle funzioni spettrali dopo la sottrazione delle masse dei quark costituenti.
   • Larghezze di Decadimento Termico ($\Gamma_n$): Estratte dalla larghezza a metà altezza (FWHM) dei picchi, utilizzando un profilo Breit-Wigner asimmetrico per tenere conto dell'asimmetria vicino alla soglia del continuo.
   • Temperature di Dissociazione ($T_d$): Identificate come la temperatura alla quale il picco spettrale si fonde nel continuo (l'energia di legame svanisce).
3. Probabilità di Sopravvivenza: Per un sistema che subisce un'espansione di Bjorken, la probabilità di sopravvivenza di uno stato preformato è calcolata integrando la larghezza di decadimento dipendente dalla temperatura $\Gamma(T(t))$ nel tempo.
4. Modello di Ricombinazione: Il quadro è esteso per includere la ricombinazione proiettando gli operatori di salto stocastici dell'equazione di Lindblad sullo spazio degli stati legati. Sotto un'approssimazione adiabatica (dove l'Hamiltoniana efficace evolve lentamente rispetto alle scale temporali dei salti quantistici), viene derivato un tasso di coalescenza. Questo tasso è proporzionale al quadrato del numero di quark pesanti liberi ($N^2$) e al volume del mezzo ($V$), modulato da un fattore di rilassamento per tenere conto del tempo finito richiesto affinché i quark pesanti si termalizzino.

Contributi Chiave

• Quadro Unificato: Il documento deriva sia i meccanismi di soppressione sia quelli di ricombinazione da un'unica equazione del moto (l'equazione di Lindblad) per la matrice di densità $Q\bar{Q}$, garantendo coerenza interna tra la parte immaginaria del potenziale (che guida il decadimento) e gli operatori di salto (che guidano la ricombinazione).
• Derivazione dai Primi Principi: L'approccio evita combinazioni ad hoc di modelli indipendenti. Il potenziale complesso fissa le larghezze di decadimento termico, che determinano direttamente sia la probabilità di sopravvivenza sia la normalizzazione del tasso di ricombinazione tramite il teorema di fluttuazione-dissipazione.
• Coalescenza Adiabatica: Gli autori derivano un modello specifico di coalescenza dall'equazione di Lindblad, incorporando un fattore di rilassamento per correggere la termalizzazione non istantanea dei quark pesanti nel mezzo.

Risultati
Lo studio applica questo quadro agli stati di charmonio ($J/\psi, \psi(2S)$) e bottomonio ($\Upsilon(nS)$) nelle collisioni Pb–Pb a $\sqrt{s_{NN}} = 5.02$ TeV.

• Temperature di Dissociazione:
  • Charmonio: $T_{J/\psi} \approx 372$ MeV; $T_{\psi(2S)} \lesssim 155$ MeV (si dissocia vicino o sotto la temperatura critica).
  • Bottomonio: $T_{\Upsilon(1S)} \approx 564$ MeV; $T_{\Upsilon(2S)} \approx 223$ MeV; $T_{\Upsilon(3S)} \approx 164$ MeV; $T_{\Upsilon(4S)} < 155$ MeV.
• Fattore di Modifica Nucleare ($R_{AA}$):
  • $J/\psi$: Il modello prevede un contributo significativo di ricombinazione nelle collisioni centrali ($N_{part} \gtrsim 200$), guidato dall'alta molteplicità di quark charm alle energie LHC. Questa rigenerazione compensa parzialmente la soppressione, risultando in un $R_{AA}$ totale vicino all'unità, coerente con i dati ALICE e CMS.
  • $\Upsilon(1S)$: Poiché la molteplicità di quark bottom è circa due ordini di grandezza inferiore a quella del charm, il termine di ricombinazione è trascurabile. L'$R_{AA}$ per $\Upsilon(1S)$ è dominato quasi interamente dalla soppressione, con i risultati del modello che concordano ragionevolmente bene con i dati CMS senza parametri aggiuntivi di adattamento.

Significato e Affermazioni
Gli autori affermano che questo lavoro fornisce una "descrizione unificata, ispirata ai primi principi" della produzione di quarkonio. Derivando sia la soppressione sia la ricombinazione dalle stesse dinamiche microscopiche (il potenziale complesso e gli operatori di salto di Lindblad), il quadro elimina la necessità di un aggiustamento fenomenologico indipendente di questi effetti concorrenti. Il documento afferma che una volta fissato il potenziale complesso (tramite il modello della legge di Gauss), tutti gli effetti del mezzo sulla resa del quarkonio seguono senza parametri liberi aggiuntivi. Questo rappresenta un passo significativo verso una descrizione veramente basata sui primi principi del quarkonio nelle collisioni di ioni pesanti, offrendo un'alternativa più fondata teoricamente ai modelli di trasporto tradizionali in cui soppressione e rigenerazione sono spesso combinate manualmente.