Is the coexistence of strange quark stars and hadronic stars favored by astrophysical data? A Bayesian analysis

Sintesi Tecnica: Coesistenza di Stelle di Quark Strane e Stelle Adroniche

Definizione del Problema
La composizione microscopica delle stelle compatte e l'equazione di stato (EOS) della materia fortemente interagente rimangono problemi centrali nella fisica nucleare e nell'astrofisica. Esiste una specifica tensione tra i vincoli osservativi che richiedono un'EOS "morbida" a densità intermedie (per spiegare i piccoli raggi degli oggetti di massa ordinaria e la deformabilità tidale da GW170817) e la necessità di un'EOS "rigida" ad alte densità (per sostenere pulsar con masse superiori a $2 M_\odot$). Inoltre, l'inclusione di iperoni e risonanze $\Delta$ nella materia adronica tipicamente ammorbidisce l'EOS, portando al "problema degli iperoni" dove le masse massime scendono al di sotto dei valori osservati.

Lo scenario a "due famiglie" (2F) propone una risoluzione: la coesistenza di due rami disconnessi di stelle compatte. In questo quadro, le stelle adroniche (HS) sono metastabili e popolano un ramo con raggi minori, mentre la materia di quark strana (SQM) assolutamente stabile forma un ramo auto-legato di stelle di quark (QS) capace di sostenere masse più elevate. Sebbene questo scenario sia stato esplorato teoricamente, non era stato precedentemente sottoposto a un'analisi bayesiana rigorosa che confrontasse direttamente l'ipotesi standard a "una sola famiglia" (1F) (dove tutte le stelle sono adroniche) utilizzando un insieme unificato di dati astrofisici e di laboratorio.

Metodologia
Gli autori eseguono la prima analisi bayesiana dettagliata per vincolare i parametri dell'EOS per entrambi gli scenari 1F e 2F simultaneamente.

• Modelli:
  • Settore Adronico: Un modello di campo medio relativistico non lineare (RMF) che include nucleoni, iperoni ($\Lambda, \Sigma, \Xi$) e risonanze $\Delta$. Il modello utilizza parametri empirici nucleari alla saturazione e accoppiamenti dipendenti dalla densità per il mesone $\rho$ per controllare la pendenza dell'energia di simmetria ($L_{sym}$).
  • Settore di Quark (solo 2F): Un'EOS del modello a sacchetto (bag-model) con color-flavor-locking (CFL) che descrive la SQM assolutamente stabile. Questo settore include i parametri per la costante del sacchetto ($B$), il gap di accoppiamento CFL ($\Delta_{CFL}$) e le correzioni perturbative della QCD ($a_4$).
• Dati e Verosimiglianze: L'analisi incorpora:
  • Vincoli Massa-Raggio (M-R): Posteri joint da NICER (PSR J0030+0451, PSR J0740+6620, PSR J0614–3329) e osservazioni a raggi X di HESS J1731–347.
  • Vincoli basati solo sulla Massa: PSR J0952–0607 e la componente secondaria di GW190814 (trattata come un potenziale oggetto compatto).
  • Onde Gravitazionali: Il posteriore joint completo per GW170817 (masse e deformabilità tidali), trattato con marginalizzazione sulle etichette di famiglia (HS-HS, HS-QS, QS-HS, QS-QS) nello scenario 2F.
  • Vincoli di Laboratorio/Teorici: Bande della teoria dell'campo efficace chirale ($\chi$EFT) per la materia neutronica pura e dati di flusso di collisioni di ioni pesanti (HIC) per la materia nucleare simmetrica.
• Framework Statistico:
  • Lo studio utilizza il campionamento annidato dinamico (tramite il pacchetto dynesty) per calcolare le distribuzioni posteriori e le evidenze bayesiane ($Z$).
  • Il fattore di Bayes ($B_{2F,1F} = Z_{2F}/Z_{1F}$) viene utilizzato per confrontare i modelli, tenendo conto delle diverse dimensionalità dello spazio dei parametri (11 parametri per 1F contro 14 per 2F).
  • Nello scenario 2F, l'assegnazione della famiglia di ciascuna sorgia osservata è trattata come una variabile latente e marginalizzata, con prior definiti in base alla stabilità strutturale e ai canali di formazione astrofisica.

Contributi Chiave

1. Primo Confronto Bayesiano: Questo lavoro fornisce il primo confronto bayesiano rigoroso tra gli scenari 1F e 2F utilizzando un dataset completo che include sia vincoli nucleari a bassa densità che osservazioni astrofisiche ad alta massa.
2. Framework Unificato: Esso vincola simultaneamente l'EOS adronica metastabile e l'EOS della SQM assolutamente stabile all'interno dell'ipotesi 2F, invece di assumere a priori la natura delle sorgenti osservate.
3. Risoluzione della Tensione: L'analisi dimostra come lo scenario 2F risolva naturalmente la tensione tra la morbidezza a densità intermedia e la rigidità ad alta densità separando questi requisiti in due distinti rami stellari.

Risultati

• Preferenza del Modello: Lo scenario 2F è favorito rispetto allo scenario 1F in tutte le combinazioni di dataset. Per il dataset completo (inclusi HIC e tutti i vincoli astrofisici), il log-fattore di Bayes è $\ln B_{2F,1F} = 11.9$, indicando un supporto "decisivo" per lo scenario a due famiglie.
• Caratteristiche dell'EOS:
  • Scenario 1F: Per soddisfare sia i vincoli sul piccolo raggio che il limite di $2 M_\odot$, l'EOS 1F richiede un compromesso finemente bilanciato. Ciò è ottenuto attraverso un esordio precoce delle risonanze $\Delta$ ($\sim 1.8 n_{sat}$) che ammorbidisce l'EOS a densità intermedie, seguito da un forte irrigidimento ad alte densità. La massa massima è vincolata a $M_{max} \approx 2.05 M_\odot$.
  • Scenario 2F: La tensione viene alleviata. Il ramo adronico rimane morbido (supportato da risonanze $\Delta$ e iperoni) con una massa massima di $M_{HS,max} \approx 1.54 M_\odot$ e un raggio di $R_{1.4} \approx 11.1$ km. Le pulsar massicce sono sostenute dal più rigido ramo QS, raggiungendo $M_{QS,max} \approx 2.33 M_\odot$.
• Vincoli dei Parametri: Il posteriore 2F preferisce una massa efficace di Dirac ($m^*/m$) maggiore e un tuning meno estremo degli accoppiamenti $\Delta$ rispetto al caso 1F. I parametri del settore di quark ($B, \Delta_{CFL}, a_4$) sono moderatamente vincolati, con una correlazione positiva tra la costante del sacchetto e il gap di accoppiamento.
• Interpretazione di GW170817: Nello scenario 2F, l'evento GW170817 è molto probabilmente una fusione HS-HS, ma i canali misti (HS-QS) contribuiscono significativamente (~35% della verosimiglianza), mentre il canale QS-QS è disfavorito per la migliore EOS.

Significatività e Rivendicazioni
L'articolo afferma che i dati astrofisici e di laboratorio attuali sono più naturalmente organizzati da due famiglie di stelle compatte disconnesse piuttosto che da una singola sequenza puramente adronica. Lo scenario 2F offre una spiegazione naturale per la coesistenza di oggetti molto compatti di massa ordinaria (come HESS J1731–347) e pulsar molto massicce, rimuovendo efficacementamente il problema degli iperoni.

Gli autori sottolineano che la conclusione è robusta rispetto alle variazioni dei prior (potenziali degli iperoni, intervalli dei parametri SQM) e alle assunzioni sulla massa delle sorgenti, purché venga utilizzato l'intero dataset. Tuttavia, notano che la preferenza statistica è meno definitiva se si considerano dataset ristretti (escludendo HESS J1731–347, PSR J0614–3329 e i dati HIC).

Lo studio esclude esplicitamente le stelle ibride (stelle con nuclei di quark) dal confronto attuale, notando che una completa selezione del modello dovrebbe includere le stelle ibride come una terza ipotesi in lavori futuri. Gli autori concludono che future rilevazioni di oggetti compatti più massicci, oggetti di massa ordinaria molto compatti o misurazioni precise di due masse distinte con lo stesso raggio fornirebbero indicazioni ancora più forti a favore dello scenario a due famiglie.