Bayesian Inference of Neutron Star Properties in $f(Q)$ Gravity Using NICER Observations

Sintesi Tecnica: Inferenza Bayesiana delle Proprietà delle Stelle di Neutroni in Gravità $f(Q)$ utilizzando le Osservazioni NICER

Problema
Sebbene la Relatività Generale (GR) descriva con successo i fenomeni gravitazionali attraverso varie scale, essa affronta sfide concettuali riguardanti le singolarità, la materia oscura e l'energia oscura, motivando l'esplorazione delle teorie della Gravità Modificata (MG). Tra queste, la gravità teleparallea simmetrica (STEGR) e la sua estensione, la gravità $f(Q)$, offrono una distinta formulazione geometrica in cui la gravità è codificata nella non-metricità dello spaziotempo piuttosto che nella curvatura o nella torsione. Sebbene i modelli $f(Q)$ siano stati studiati in cosmologia, la loro applicazione nel regime di campo forte delle stelle di neutroni (NS) rimane ampiamente inesplorata. Nello specifico, vi è una mancanza di studi sistematici di inferenza bayesiana per vincolare i parametri dei modelli $f(Q)$ utilizzando dati osservativi, poiché le analisi precedenti si sono spesso basate su scelte fenomenologiche di parametri senza un rigoroso confronto statistico con i dati.

Metodologia
Gli autori investigano la struttura interna delle stelle di neutroni all'interno del quadro della gravità $f(Q)$ considerando tre forme funzionali rappresentative: modelli lineari, logaritmici ed esponenziali. Per isolare gli effetti della gravità modificata dalle incertezze della fisica della materia densa, lo studio adotta l'equazione di stato (EoS) realistica e dipendente dalla densità DDME2 del campo medio relativistico.

La metodologia procede come segue:

1. Quadro Teorico: Gli autori derivano le equazioni di Tolman–Oppenheimer–Volkoff (TOV) modificate per spazi-tempi statici e sfericamente simmetrici nella gravità $f(Q)$. Queste equazioni governano l'equilibrio stellare e vengono risolte numericamente utilizzando uno schema Runge-Kutta del quarto ordine (RK4).
2. Inferenza Bayesiana: Viene impiegato un quadro bayesiano per vincolare i parametri liberi ($\alpha, \beta, Q_0$) dei tre modelli $f(Q)$. L'analisi utilizza le distribuzioni a posteriori massa-raggio derivate dalle osservazioni NICER di quattro pulsar: PSR J0030+0451, PSR J0740+6620, PSR J0437+4715 e PSR J0614+3329.
3. Analisi Statistica: Lo studio calcola le distribuzioni a posteriori per i parametri del modello e per le proprietà derivate delle NS (massa, raggio, deformabilità tidale e compattezza). Il confronto tra i modelli è effettuato utilizzando il fattore di Bayes (BF) per determinare quale forma funzionale sia statisticamente preferita dai dati.

Contributi Chiave

• Prime Costrizioni Bayesiane in $f(Q)$: Questo lavoro rappresenta la prima analisi sistematica di inferenza bayesiana dei parametri dei modelli di gravità $f(Q)$ utilizzando osservabili di stelle di neutroni, andando oltre i test cosmologici verso il regime di campo forte.
• Analisi Comparativa dei Modelli: Il documento fornisce un confronto statistico rigoroso tra i modelli $f(Q)$ lineari, logaritmici ed esponenziali, valutandone la fattibilità rispetto ai dati ad alta precisione di NICER.
• Degenerazione e Correlazione dei Parametri: Lo studio mappa le correlazioni tra i parametri del modello e le proprietà delle NS, rivelando come diverse forme funzionali influenzino l'accoppiamento tra le modifiche gravitazionali e la struttura stellare.

Risultati

• Preferenza del Modello: L'analisi bayesiana, supportata dai confronti del fattore di Bayes, favorisce statisticamente il modello $f(Q)$ esponenziale rispetto ai casi lineare e logaritmico. Il modello esponenziale esibisce parametri ben vincolati con distribuzioni a posteriori quasi gaussiane e unimodali. Al contrario, il modello lineare mostra una pronunciata bimodalità e distribuzioni di parametri ampie, indicando forti degenerazioni, mentre il modello logaritmico presenta incertezze moderate.
• Proprietà delle Stelle di Neutroni: Per il modello esponenziale preferito, lo studio deriva un raggio e una deformabilità tidale a $1.4 M_\odot$ rispettivamente di $R_{1.4} = 11.27^{+0.53}_{-0.36}$ km e $\Lambda_{1.4} = 156.95^{+84.02}_{-41.73}$. Questi valori sono coerenti con i vincoli osservativi attuali.
• Massa Massima e il Mass Gap: Un risultato significativo è che tutti e tre i modelli vincolati prevedono masse massime delle stelle di neutroni che raggiungono $M_{\text{max}} \simeq 2.98 M_\odot$. Le regioni di confidenza al 95% per la relazione massa-raggio si estendono nel "lower mass gap" ($\sim 2.5\text{--}5 M_\odot$). Questa previsione emerge naturalmente dallo spazio dei parametri vincolato da Bayes senza essere imposta per costruzione.
• Correlazioni: Il modello esponenziale dimostra una forte interdipendenza strutturale tra la massa massima e altre proprietà delle NS, mentre il modello lineare mostra un'indipendenza statistica tra la massa massima e altre osservabili.

Significatività
L'articolo afferma che questi risultati evidenziano il potenziale delle stelle di neutroni come potenti sonde della gravità teleparallea simmetrica nel regime di campo forte. L'emergere di configurazioni stabili nel "lower mass gap" è presentato come una chiara previsione testabile del quadro $f(Q)$, distinguendolo dagli scenari standard della GR che tipicamente producono masse massime inferiori per la stessa EoS. Lo studio conclude che il modello $f(Q)$ esponenziale è la descrizione statisticamente più supportata della struttura delle stelle di neutroni tra i tre testati, suggerendo che gli effetti della gravità modificata possano significativamente potenziare il supporto gravitazionale ad alta densità. Gli autori osservano che future misurazioni precise della massa nella regione del "mass-gap" potrebbero servire come potente discriminatore tra la GR e gli scenari di gravità estesa. Suggeriscono inoltre che il lavoro futuro potrebbe coinvolgere l'inferenza bayesiana congiunta con parametrizzazioni EoS flessibili per affrontare ulteriormente le degenerazioni tra la gravità a campo forte e la microfisica.