Strongly interacting matter with criticality induced by modified excluded volume in core-collapse supernova simulations

Sintesi Tecnica: Materia Fortemente Interagente con Criticità Indotta da Volume Escluso Modificato nelle Simulazioni di Supernovae a Collasso Nucleare

Problema
Il meccanismo di esplosione delle supernovae a collasso nucleare (CCSN) rimane una delle questioni aperte più significative in astrofisica, fortemente dipendente dall'equazione di stato (EOS) ad alta densità. Uno scenario specifico ipotizza che una transizione di fase del primo ordine sufficientemente forte dalla materia adronica alla materia di quark deconfinata possa innescare un'esplosione riuscita. Tuttavia, studi precedenti su questo meccanismo si sono basati sull'approccio "a due fasi", dove modelli distinti di materia adronica e di quark vengono uniti tramite condizioni di Gibbs. Questo approccio introduce incertezze riguardo al salto di densità, alla densità di onset e alla pendenza della pressione tra le fasi, spesso mancando di un punto critico a densità e temperatura finite. Inoltre, i recenti vincoli dalle osservazioni di stelle di neutroni (masse $\approx 2 M_\odot$, raggi $\approx 9\text{--}13$ km) e dagli eventi di onde gravitazionali (GW170817) richiedono modelli di EOS che siano coerenti sia con la fisica nucleare che con i dati astrofisici.

Metodologia
Gli autori sviluppano una nuova classe di EOS multiuso, denominata DD2-MEV, basata sul framework del Campo Medio Relativistico (RMF) con accoppiamenti mesone-nucleone dipendenti dalla densità (specificamente il modello DD2). Invece di una tradizionale costruzione a due fasi, gli autori utilizzano un approccio a Volume Escluso Modificato (MEV).

• Meccanismo MEV: Il modello sostituisce l'interpretazione geometrica del volume escluso con un cambiamento nel numero effettivo di gradi di libertà. Un funzionale di volume escluso dipendente dal mezzo, $\Phi_N(T, n_B)$, modula i fattori di degenerazione dei nucleoni. Questo funzionale è progettato per mimare il comportamento di Van der Waals, presentando regimi sovra-critici e sotto-critici e un endpoint critico ad alte temperature.
• Consistenza Termodinamica: L'implementazione include termini di riarrangiamento derivanti dalla dipendenza in densità degli accoppiamenti e dal funzionale MEV per garantire la consistenza termodinamica.
• Modelli Comparativi: Lo studio confronta due versioni dell'EOS DD2-MEV:
  1. DD2-MEV (Gibbs): Una versione in cui la regione di instabilità di Van der Waals viene rimossa imponendo le condizioni di equilibrio di fase di Gibbs (pressione e potenziali chimici costanti).
  2. DD2-MEV (Van der Waals): Il modello continuo che mantiene la pendenza negativa della pressione e le regioni di instabilità caratteristiche di un'equazione di Van der Waals.
• Simulazioni: Gli autori eseguono simulazioni di idrodinamica delle radiazioni di neutrini in relatività generale in simmetria sferica utilizzando il codice AGILE-BOLTZTRAN. Utilizzano uno schema di trasporto di Boltzmann a sei specie e un progenitore da $40 M_\odot$ (s40a28).
• Analisi: Lo studio include un'analisi dettagliata del segnale dei neutrini e un'analisi del modo delle onde gravitazionali (GW) tramite il codice GREAT (astrosismologia in relatività generale) per tracciare i modi fondamentali ($f$) e di gravità ($g$).

Contributi Chiave e Risultati

1. Fallimento della Costruzione di Gibbs:
   Le simulazioni che utilizzano l'EOS DD2-MEV (Gibbs), che impone una standard costruzione di transizione di fase, non riescono a produrre un'esplosione di supernova. Il salto di densità tra le fasi adronica e di quark in questo modello è significativamente minore rispetto ai precedenti modelli ibridi RDF. Di conseguenza, la protoneutron star (PNS) subisce solo un lieve riassetto idrodinamico piuttosto che un collasso supersonico. La simulazione procede verso la formazione di un buco nero senza riattivazione dello shock.

2. Esplosione tramite Comportamento Continuo di Van der Waals:
   Al contrario, le simulazioni che utilizzano l'EOS continua DD2-MEV (Van der Waals) portano a un'esplosione riuscita. La transizione di fase induce un rapido aumento della densità centrale e della temperatura, portando a un collasso supersonico della PNS. Tuttavia, diversamente dai modelli precedenti dove un forte secondo shock si forma immediatamente al confine di fase, questo modello presenta una formazione dello shock ritardata. Un'onda di perturbazione sub-sonica si propaga verso l'esterno, raggiungendo lo shock di bounce stazionario dopo diversi millisecondi prima di trasformarsi in un'onda d'urto forte.

   • Ejecta: L'esplosione espelle una massa inizialmente maggiore ($\sim 0.2 M_\odot$) rispetto ai report precedenti, con un'energia diagnostica superiore a 10 Bethe ($10^{52}$ erg), sebbene un significativo fallback riduca la massa finale degli ejecta a $\sim 1 M_\odot$.

3. Segnature dei Neutrini Distinte:

   • Caso Gibbs: Non viene osservato alcun burst di millisecondi. Il segnale dei neutrini mostra solo fluttuazioni a breve termine corrispondenti al riassetto strutturale della PNS, seguite dal canonico comportamento di accrezione che porta alla formazione di un buco nero.
   • Caso Van der Waals: Viene rilasciata una segnatura di neutrini di tipo burst, ma è sostanzialmente più ampia (dura 5–10 ms) rispetto ai burst di 1–2 ms riportati in studi precedenti utilizzando modelli bag o EOS RDF. Questa allargamento è attribuito alla formazione e propagazione ritardata del secondo shock. Il burst è caratterizzato da una crescita in due stadi della luminosità: una salita iniziale lenta dovuta alla contrazione della PNS, seguita da un picco acuto mentre lo shock attraversa le neutrinosfere.

4. Analisi dei Modi delle Onde Gravitazionali:
   Lo studio esegue un'analisi dei modi delle oscillazioni della PNS ($f$ e $2g_1$).

   • Collasso delle Relazioni Standard: Gli autori riscontrano che le relazioni di fitting standard che collegano le frequenze dei modi alla densità media ($\bar{\rho}$) e alla compattezza ($\kappa$) della PNS crollano dopo l'apparizione dei gradi di libertà della materia di quark.
   • Comportamenti Distinti: Nel caso fallito (Gibbs), la pendenza del modo $f$ cambia bruscamente a causa dei cambiamenti strutturali, mentre il modo $2g_1$ diminuisce rapidamente. Nel caso riuscito (Van der Waals), entrambi i modi mostrano grandi fluttuazioni durante la transizione di fase e la formazione dello shock, indicando che la loro eccitazione è governata da fenomeni fisici distinti dalla semplice scalatura di densità/compattezza.

Significatività e Rivendicazioni
L'articolo sostiene di superare le criticità del tradizionale approccio a due fasi utilizzando un'EOS continua che imita una transizione di fase del primo ordine attraverso un meccanismo di volume escluso modificato. La principale significatività risiede nel dimostrare che i dettagli della costruzione della transizione di fase sono critici per il meccanismo di esplosione CCSN.

• Lo studio suggerisce che il "salto di densità" è un parametro cruciale; un piccolo salto (come nella costruzione di Gibbs di questo specifico modello MEV) impedisce l'esplosione, mentre il comportamento continuo di Van der Waals permette l'esplosione.
• Evidenzia come la durata del burst di neutrini e l'evoluzione del modo GW siano sonde sensibili capaci di distinguere tra diversi scenari di transizione di fase (fallito vs riuscito, costruzione continua vs Gibbs).
• Gli autori concludono con modestia che, sebbene questi risultati evidenzino grandi differenze rispetto ai precedenti modelli di esplosione guidati dalla QCD, essi si basano su una specifica parametrizzazione. Sottolineano la necessità di simulazioni multidimensionali per confermare questi risultati e per calcolare segnali GW diretti, notando che l'analisi corrente evita soluzioni spurie annullando la velocità del suono nelle regioni instabili.