Sintesi Tecnica: Studio della materia QCD nelle collisioni tra ioni pesanti utilizzando il motore di calcolo MUSES

Definizione del Problema

L'equazione di stato (EoS) della materia QCD calda e densa è fondamentale per descrivere le collisioni tra ioni pesanti a tutte le energie di collisione. Tuttavia, determinare l'EoS presenta due sfide primarie:

1. Limitazioni Teoriche: I calcoli della Lattice QCD (LQCD) basati sui primi principi sono limitati a potenziali chimici barionici reali ($\mu_B$) nulli o piccoli a causa del problema del segno dei fermioni. Sebbene la LQCD fornisca dati affidabili ad alte temperature e bassi $\mu_B$, sono necessari modelli fenomenologici per densità più elevate. Garantire che questi modelli corrispondano ai vincoli della LQCD in modo termodinamicamente coerente rimane un ostacolo significativo.
2. Complessità Dinamica: Le collisioni tra ioni pesanti sono sistemi intrinsecamente dinamici e fuori dall'equilibrio. Il confronto tra i dati sperimentali e la teoria richiede simulazioni idrodinamiche che evolvano molteplici cariche conservate (numero barionico, strangeness e carica elettrica). Queste simulazioni richiedono input di EoS che siano non solo termodinamicamente coerenti, ma anche invertibili (conversione da $(T, \mu_X)$ alle variabili naturali dell'idrodinamica come densità di entropia $s$ e densità di carica $n_X$) e capaci di gestire transizioni di fase, inclusi potenziali punti critici (CP) e transizioni del primo ordine.

Inoltre, connettere i vincoli dalle collisioni tra ioni pesanti alle stelle di neutroni richiede un quadro unificato capace di gestire il diagramma di fase quadridimensionale $(T, \mu_B, \mu_S, \mu_Q)$, poiché i sistemi esplorano regioni diverse di questo spazio pur avendo potenzialmente densità barioniche simili.

Metodologia

Gli autori utilizzano il motore di calcolo MUSES (Modular Unified Solver of the Equation of State), specificamente la nuova release Calliope. Questo framework fornisce un sistema di gestione del workflow composibile che orchestra una collezione di moduli software per calcolare EoS e osservabili.

Moduli Chiave e Workflow

Il documento si concentra su moduli progettati per la fenomenologia degli ioni pesanti, espandendo l'insieme precedente di 10 moduli a 17:

• Fase Adronica: Il modulo HRG (Hadron Resonance Gas), basato su Thermal-FIST, calcola la termodinamica nell'intero spazio 4D $(T, \mu_B, \mu_S, \mu_Q)$. Supporta interazioni ideali, a volume escluso e interazioni Quantum van der Waals (QvdW).
• Espansioni Lattice QCD:
  • 4D-Taylor (BQS): Espande la pressione in potenze di $\hat{\mu}_X = \mu_X/T$ utilizzando suscettibilità stimate nel continuo.
  • 4D-TExS: Generalizza lo schema di espansione $T'$ a 4D, mappando i potenziali chimici cartesiani in coordinate sferiche per estendere l'intervallo di validità oltre le standard espansioni di Taylor.
• Modelli Fenomenologici:
  • Ising 2D-TExS: Combina l'espansione $T'$ del reticolo 2D con lo scaling critico di Ising 3D per modellare un punto critico e una linea del primo ordine nel piano $(T, \mu_B)$.
  • Olografia (NumRelHolo): Utilizza un modello di corrispondenza gauge/gravità Einstein-Maxwell-Dilaton (EMD) per descrivere la materia fortemente accoppiata, calibrato sulla LQCD a $\mu_B=0$ ed esteso a densità finite. Include coefficienti di trasporto.
• Modulo di Sintesi: Una nuova procedura basata sulla miscelazione statistica e sui criteri di stabilità permette di fondere due diverse EoS (ad esempio, Olografia e HRG) in una singola EoS termodinamicamente coerente che copre un intervallo più ampio di $(T, \mu_X)$. Questa procedura può generare transizioni di crossover, transizioni del primo ordine o punti critici.
• Invertitore di EoS: Un modulo che inverte numericamente le tabelle EoS dalle coordinate $(T, \mu_X)$ alle variabili naturali dell'idrodinamica $(s, n_X)$ o $(\epsilon, n_X)$, essenziale per i solver idrodinamici lagrangiani come CCAKE.

Framework di Simulazione

Gli autori eseguono simulazioni idrodinamiche viscose relativistiche utilizzando il framework NuclearConfectionary (specificamente AMPT per le condizioni iniziali, CCAKE 2.0 per l'idrodinamica e SMASH per l'afterburner adronico).

• Collisioni: Collisioni centrali Au+Au a $\sqrt{s_{NN}} = 7.7, 19.6,$ e $39$ GeV.
• Variabili: Le simulazioni tracciano le singole particelle SPH (Smoothed Particle Hydrodynamics) attraverso il diagramma di fase, permettendo l'analisi di traiettorie dinamiche, mappe di calore e distribuzioni di freeze-out.
• Coefficienti di Trasporto: La viscosità di taglio è impostata come costante $\eta T/(\epsilon+P) = 0.08$. La viscosità di bulk include una componente di scaling critico mimata tramite la suscettibilità barionica $\chi_B^2$.

Contributi Chiave

1. Rilascio di MUSES Calliope: Il documento introduce un nuovo set di moduli e workflow specificamente adattati per le collisioni tra ioni pesanti, completando l'esistente set di moduli per le stelle di neutroni.
2. Copertura EoS Estesa: Gli autori dimostrano la costruzione di EoS che coprono lo spazio di fase rilevante per il RHIC Beam Energy Scan (BES-II), incluse le regioni con alto $\mu_B$ e basso $T$ dove le espansioni LQCD tipicamente falliscono.
3. Fusione Termodinamica: L'implementazione del modulo Synthesis consente di fondere senza soluzione di continuità modelli EoS disparati (ad esempio, Olografia e HRG) con punti critici regolabili, garantendo la coerenza termodinamica attraverso i confini di fase.
4. Inversione e Stabilità: Il modulo EoS Inverter è testato per la sua capacità di gestire la mappatura non lineare tra $(T, \mu_X)$ e $(s, n_X)$, incluse le strategie per gestire transizioni di fase del primo ordine e regimi metastabili/instabili.
5. Analisi delle Traiettorie Dinamiche: Utilizzando un solver idrodinamico lagrangiano, gli autori mappano le specifiche traiettorie delle celle di fluido attraverso il diagramma di fase QCD, identificando quali regioni vengono sondate a diverse energie di fascio e rapidità.

Risultati

• Validità e Copertura EoS:
  • L'EoS 4D-TExS fornisce una copertura significativamente migliore nel piano $(s, n_B)$ rispetto alla standard espansione 4D-Taylor, particolarmente a temperature più basse e densità barioniche più elevate.
  • La EoS fusa Holography+QvdW-HRG offre la copertura più ampia, coprendo efficacemente sia le fasi adroniche che quelle deconfinate.
  • Le simulazioni a $\sqrt{s_{NN}} = 19.6$ GeV usando l'EoS 4D-Taylor mostrano che solo il ~50-70% delle celle di fluido a rapidità centrale rimane entro la tabella EoS valida, rendendo necessario l'uso di EoS conformi di backup. Al contrario, i modelli 4D-TExS e Holography+HRG mantengono una copertura >90% a rapidità centrale per $\sqrt{s_{NN}} \ge 19.6$ GeV.
• Sensibilità al Punto Critico:
  • La posizione del punto critico influenza significativamente quali energie di fascio sono sensibili alla sua presenza.
  • Per un punto critico a $\mu_B^c = 598$ MeV (coerente con recenti previsioni teoriche), solo l'energia più bassa ($\sqrt{s_{NN}} = 7.7$ GeV) mostra una frazione significativa di celle di fluido ($\sim$40%) che passa vicino alla regione critica. Le energie superiori sondano la regione critica principalmente a grandi rapidità pro o retro.
  • Non è stata osservata alcuna significativa "lente critica" (raggruppamento di traiettorie isentropiche) nelle simulazioni, suggerendo che la regione critica possa essere orientata orizzontalmente lungo l'asse $\mu_B$ piuttosto che l'asse $T$, o che gli effetti fuori dall'equilibrio ne smorzino il segnale.
• Transizioni del Primo Ordine:
  • Le simulazioni che includono regimi metastabili e instabili attraverso una transizione del primo ordine non hanno causato crash numerici nel solver CCAKE, poiché la regione di velocità del suono al quadrato negativa ($c_s^2 < 0$) era estremamente piccola e di breve durata ($\langle \Delta \tau \rangle \approx 1.4 \times 10^{-5}$ fm/c).
  • Tuttavia, gli autori notano che l'attuale trattamento della linea del primo ordine è un'esplorazione preliminare e non una teoria dinamica completa della conversione di fase.
• Osservabili:
  • La scelta dell'EoS (e della posizione del punto critico) ha avuto un effetto minimo sulla molteplicità totale dei pioni ($dN/dy$), suggerendo che i parametri di normalizzazione dominano questi osservabili.
  • Si è scoperto che la larghezza della distribuzione di rapidità ($dN/d\eta$) dipende fortemente dalla lunghezza di smearing dello stato iniziale in AMPT piuttosto che dai dettagli dell'EoS.

Significato e Prospettive

Il documento stabilisce l'infrastruttura necessaria per studi sistematici della materia QCD densa all'interno di un quadro unificato e open-source. Fornendo EoS modulari e termodinamicamente coerenti che possono essere fuse e invertite, il motore MUSES Calliope consente confronti più robusti tra i dati delle collisioni tra ioni pesanti e i modelli teorici.

Gli autori sottolineano che il loro lavoro è un "primo passo" e un' "esplorazione preliminare", in particolare riguardo al trattamento delle transizioni di fase del primo ordine e all'inclusione dello scaling critico nei coefficienti di trasporto. Dichiarano esplicitamente che un trattamento sistematico delle transizioni del primo ordine nelle simulazioni dinamiche richiede ulteriori lavori teorici e numerici. Il significato primario risiede nella capacità di tracciare quali regioni del diagramma di fase QCD vengono campionate durante l'evoluzione del fireball e di quantificare come la validità limitata delle EoS basate su espansioni influenzi le predizioni idrodinamiche. Questo framework è posizionato per facilitare future analisi congiunte che colleghino i vincoli delle collisioni tra ioni pesanti e delle stelle di neutroni.