Are NICER and GW170817 constraints suggesting a compactified scenario for Neutron stars?

Sintesi Tecnica: I vincoli di NICER e GW170817 suggeriscono uno scenario di compattificazione per le stelle di neutroni?

Problema
L'equazione di stato (EoS) della materia a densità superiori a quella nucleare di saturazione rimane un problema fondamentale non ancora risolto della fisica. Mentre la teoria dell'efficacia chirale ($\chi$EFT) modella in modo affidabile le densità sub-nucleari e la QCD perturbativa descrive i regimi asintotici ad alta densità, la regione intermedia di densità all'interno dei nuclei delle stelle di neutroni (NS) è scarsamente compresa. Recenti osservazioni astrofisiche provenienti dal Neutron Star Interior Composition Explorer (NICER) e dall'evento di onde gravitazionali GW170817 hanno fornito vincoli precisi su masse, raggi e deformabilità tidale delle NS. Tuttavia, rimane poco chiaro se questi vincoli favoriscano una classe specifica di EoS (monotona, non monotona o discontinua) o se rivelino un modello strutturale universale della materia densa. Nello specifico, il documento indaga se i dati attuali suggeriscano uno scenario di "compattificazione" in cui la gravità favorisce naturalmente configurazioni che massimizzano la compattezza stellare piuttosto che quelle che massimizzano semplicemente la massa.

Metodologia
Gli autori utilizzano un quadro completo di inferenza bayesiana per analizzare tre distinti ensemble di EoS, indipendenti dalla classe:

1. Monotona: Caratterizzata da una velocità del suono ($c_s^2$) che aumenta monotonicamente con la densità (tipica delle EoS adroniche).
2. Non monotona: Caratterizzata da un picco in $c_s^2$ seguito da una diminuzione (associata a transizioni di fase fluide).
3. Discontinua: Caratterizzata da un salto di densità ($\Delta n$) in una transizione di fase (associata a transizioni di fase del primo ordine forti).

Costruzione dell'EoS:
L'EoS è costruita utilizzando un approccio ibrido. La materia a bassa densità ($< 0.5 n_0$) è definita dall'EoS tabulata BPS. La regione da $0.5 n_0$ a $1.1 n_0$ è modellata tramite una politrope con un indice adiabatico $\Gamma$ vincolato dalle bande di incertezza della $\chi$EFT. Per densità superiori a $1.1 n_0$, il quadrato della velocità del suono $c_s^2$ è parametrizzato come una funzione lineare a tratti del potenziale chimico $\mu$ attraverso quattro segmenti randomizzati, estendendosi fino al regime pQCD ($\mu \approx 2.6$ GeV).

Framework di Inferenza:
A differenza dei precedenti metodi di campionamento per rifiuto che generano milioni di campioni di EoS e li filtrano, questo studio utilizza l'algoritmo di campionamento annidato UltraNest per esplorare direttamente la distribuzione a posteriori. Il modello contiene 10 parametri liberi (potenziali chimici e valori del quadrato della velocità del suono per i segmenti, oltre a $\Gamma_{CET}$ e un parametro di salto di densità $\delta n$ per la classe discontinua).
L'analisi incorpora:

• Dati NICER: Distribuzioni a posteriori congiunte per massa e raggio di PSR J0030+0451, PSR J0437−4715, PSR J0740+6620 e del recentemente riportato PSR J0614−3329.
• Dati GW170817: Campioni a posteriori di masse componenti e deformabilità tidale dalla fusione di stelle di neutroni binarie.
• Vincoli $\chi$EFT: Verosimiglianze gaussiane per la pressione a basse densità.

Lo studio confronta le distribuzioni a posteriori dei parametri dell'EoS e delle proprietà stellari derivate (sequenze Massa-Raggio) sotto diverse combinazioni di vincoli osservativi.

Risultati Chiave

1. Potere di Vincolo: L'analisi bayesiana conferma che i dati osservativi vincolano efficacemente la regione a bassa densità dell'EoS (specificamente $\Gamma_{CET}$ e il primo segmento della velocità del suono $c_{s,1}^2$), che rimane robusto in tutte le classi di EoS. Al contrario, i parametri ad alta densità ($c_{s,2}^2, c_{s,3}^2$) esibiscono distribuzioni a posteriori più ampie, indicando una significativa incertezza rimanente.
2. Impatto dei Nuovi Dati: L'inclusione di PSR J0614−3329 e GW170817 sposta sistematicamente i parametri del potenziale chimico verso valori più elevati e riduce i limiti della massa massima. Tuttavia, l'evidenza bayesiana ($\ln Z$) non favorisce fortemente una classe di EoS rispetto a un'altra; le differenze nell'evidenza rientrano nel range "inconcludente o debole".
3. Massa Massima vs. Massima Compattezza: Viene tracciata una distinzione critica tra l'EoS che produce la massa possibile massima ($M_{TOV}$) e l'EoS che produce la massima compattezza ($C = GM/Rc^2$).
   • L'EoS a Massa Massima è rigida da densità molto basse, fornendo un supporto di pressione uniforme per sostenere le masse più grandi. Queste configurazioni occupano il bordo a grande raggio della distribuzione a posteriori M-R.
   • L'EoS a Massima Compattezza presenta un'EoS a bassa-media densità più morbida (minimizzando il raggio) seguita da un marcato irrigidimento ad alte densità (supportando la massa).
4. Preferenza Osservativa: L'EoS più favorita dai dati combinati NICER e GW170817 (massimo $\ln Z$) rispecchia strutturalmente l'EoS a Massima Compattezza, non l'EoS a Massa Massima. L'EoS osservazionalmente favorita presenta una relativa morbidezza a densità intermedie, transitando verso la rigidità ad alte densità.
5. Transizioni di Fase: Per la classe discontinua, i dati preferiscono piccole discontinuità di densità ($\delta n \approx 0.09-0.10$ fm$^{-3}$) e transizioni di fase tardive (che avvengono a densità più elevate). Sebbene una transizione di fase a bassa densità possa teoricamente produrre una stella altamente compatta con un grande nucleo di quark, le previsioni bayesiane favoriscono configurazioni in cui la transizione di fase avviene più tardi, confinando la materia di quark ai nuclei di stelle sufficientemente massicce.

Significato e Rivendicazioni
Il documento afferma che gli attuali vincoli astrofisici stanno selezionando configurazioni di EoS che producono stelle compatte piuttosto che massime in massa. La convergenza dell'EoS favorita dalle osservazioni con la sequenza di massima compattezza suggerisce una proprietà fondamentale della materia densa: l'EoS probabilmente subisce un significativo ammorbidimento a densità intermedie (potenzialmente segnalando l'inizio di nuovi gradi di libertà) prima di irrigidirsi nuovamente per sostenere le osservate stelle di neutroni da $\sim 2 M_\odot$.

Gli autori concludono che questa struttura "morbida-a-rigida" è ciò che la gravità favorisce naturalmente per raggiungere la compattificazione. Sebbene i dati non escludano definitivamente classi specifiche di EoS, essi puntano verso un modello universale dove l'equazione di stato transita dalla morbidezza alla rigidità, spingendo le stelle di neutroni verso i loro limiti di compattezza. Il documento nota con modestia che futuri vincoli sulle pulsar a bassa massa (per definire meglio l'EoS a bassa densità) e sulle pulsar ad alta massa (per vincolare l'EoS ad alta densità) saranno necessari per raffinare ulteriormente queste conclusioni e potenzialmente escludere scenari specifici, come le transizioni di fase precoci.