Forming Double Neutron Stars using Detailed Binary Evolution Models with POSYDON: Comparison to the Galactic Systems

Sintesi Tecnica: Formazione di Doppie Stelle di Neutroni utilizzando Modelli di Evoluzione Binaria Dettagliati con POSYDON

Definizione del Problema
Le doppie stelle di neutroni (DNS) fungono da laboratori critici per la comprensione dell'evoluzione delle stelle massicce, dei meccanismi di supernova e delle sorgenti di onde gravitazionali. Sebbene siano state identificate oltre due dozzine di DNS nella Via Lattea, i loro canali di formazione rimangono complessi e oggetto di dibattito. Gli studi precedenti sulla popolazione sintetica si sono ampiamente basati su codici di sintesi della popolazione binaria (BPS) rapidi, che utilizzano formule di fitting o librerie di tracce di singole stelle. Tali approcci spesso mancano di un trattamento coerente delle interazioni binarie, in particolare durante la fase di involucro comune (CE) e le supernovae a collasso del nucleo (SN). Inoltre, la popolazione osservata di DNS galattiche esibisce una distinta biforcazione nei periodi orbitali e nelle eccentricità che i modelli standard faticano a riprodurre senza un raffinamento dei parametri (fine-tuning). Questo articolo affronta la necessità di vincolare le fasi evolutive dettagliate che portano alla formazione di DNS utilizzando un quadro più rigoroso e auto-coerente.

Metodologia
Gli autori utilizzano POSYDON (POpulation SYnthesis with Detailed binary-evolution simulatiONs), un codice di sintesi della popolazione binaria all'avanguardia che integra modelli dettagliati di evoluzione stellare e binaria basati sul codice 1D MESA. A differenza dei codici BPS rapidi, POSYDON impiega griglie di evoluzione stellare pre-calcolate e calcoli on-the-fly per gestire la fisica complessa, inclusi il frenamento magnetico, la radiazione di onde gravitazionali, l'evoluzione del CE e i calci di supernova (supernova kicks).

Lo studio simula $2 \times 10^6$ popolazioni binarie attraverso 25 distinte variazioni di modello (riassunte nella Tabella 2). Le variabili chiave manipolate includono:

1. Trattamenti dell'Involucro Comune (CE): Variazioni dell'efficienza di espulsione del CE ($\alpha_{CE}$), la definizione del confine nucleo-involucro (frazioni di massa di idrogeno dell'1%, 10% e 30%) e il trattamento della perdita di massa post-CE (trasferimento di massa stabile vs perdita di vento vs fase singola).
2. Modelli di Supernova (SN): Diverse prescrizioni per gli esiti del collasso del nucleo (Sukhbold et al. 2016; Patton & Sukhbold 2020; Fryer et al. 2012), distinguendo tra collasso del nucleo di ferro (CCSN) e supernovae a cattura elettronica (ECSN).
3. Prescrizioni del Calcio Natale (Natal Kick): Confronto tra la distribuzione standard di Hobbs-Maxwellian ($\sigma \approx 265$ km/s) contro modelli basati sulla fisica:
   • AsymEj: Calci derivati da meccanismi di espulsione asimmetrica (Janka 2017), che scalano con la massa dell'espulso.
   • Linear: Calci che scalano linearmente con la massa dell'espulso (Richards et al. 2023).
   • Low Maxwellian: Un modello a dispersione ridotta ($\sigma \approx 20$ km/s).

I modelli assumono una metallicità solare e una storia di formazione stellare costante. Le popolazioni sintetiche risultanti sono confrontate qualitativamente con il campione osservato di 25 DNS galattiche (Tabella 1), concentrandosi sulle distribuzioni del periodo orbitale ($P_{orb}$) ed eccentricità ($e$), nonché sui tassi di fusione (merger rates).

Contributi Chiave e Risultati

• Biforcazione dei Canali di Formazione: Lo studio identifica una fondamentale biforcazione nel canale di formazione CE, determinata dallo stato evolutivo della stella donatrice all'inizio del trasferimento di massa di Roche (RLO) con la prima stella di neutroni nata:

  • Canale del Nucleo di Elio (RLO Caso B): La donatrice avvia il trasferimento di massa mentre possiede un nucleo di elio. Questi sistemi hanno tipicamente periodi orbitali più brevi e energie di legame dell'involucro più elevate.
  • Canale del Nucleo C/O (RLO Caso C): La donatrice evolve ulteriormente sviluppando un nucleo di carbonio-ossigeno prima del RLO. Questi sistemi possiedono involucri debolmente legati e risultano in periodi orbitali più ampi.
  • Significato: Questa biforcazione spiega naturalmente la divisione nella popolazione di DNS galattiche tra sistemi a periodo breve che fondono (Sotto-popolazione i) e sistemi ampi che non fondono (S sotto-popolazione ii).

• Vincoli sulla Fisica dell'Involucro Comune (CE): La formazione del canale del nucleo di elio (che produce le DNS in fusione) è altamente sensibile all'efficienza del CE. I modelli rivelano che questo canale sopravvive solo se:

  • Viene utilizzata una definizione del nucleo generosa (frazione di H $\simeq 30\%$), o
  • Viene adottata un'alta efficienza di espulsione del CE ($\alpha_{CE} \gtrsim 1.2$).
    Senza queste condizioni, i sistemi con nucleo di elio si fondono all'interno del CE. Il canale del nucleo C/O è più robusto e sopravvive anche con valori di $\alpha_{CE}$ inferiori.

• Vincoli sul Calcio Natale (Natal Kick): La distribuzione standard di Hobbs-Maxwellian ($\sigma \approx 265$ km/s) è esclusa poiché interrompe la maggior parte delle binarie durante la prima supernova, risultando in tassi di fusione di ordini di grandezza inferiori a quelli osservati. La popolazione di DNS galattiche osservata è meglio riprodotta da prescrizioni di calcio a bassa velocità ($\lesssim 50$ km/s). I modelli "basati sulla fisica" AsymEj e Linear, che producono naturalmente calci più bassi per i progenitori strippati, riproducono con successo sia la distribuzione orbitale che la densità del tasso di fusione locale ($O(10) - O(100)$ Gpc$^{-3}$ yr$^{-1}$).

• Riciclo dei Pulsar e Accrezione di Massa: Lo studio trova che il riciclo delle pulsar (l'accelerazione della rotazione della prima stella di neutroni nata) avviene prevalentemente durante la fase di trasferimento di massa stabile di tipo BB con la compagna di stella di elio, piuttosto che durante la fase CE. I sistemi formati tramite il canale del nucleo di elio subiscono un riciclo più efficiente, coerentemente con i periodi di rotazione più brevi osservati nella Sotto-popolazione (i).

• Tipi di Supernova e Progenitori:

  • La prima supernova è prevalentemente un collasso del nucleo di ferro (CCSN).
  • La seconda supernova coinvolge spesso un progenitore ultra-strippato. Circa il 60% dei progenitori di DNS in fusione sta superando il limite del lobo di Roche al momento della seconda SN.
  • Le DNS in fusione nascono da progenitori con masse medie dell'involucro pre-SN di $\sim 0.2 M_\odot$. Ciò suggerisce che una frazione significativa di DNS in fusione potrebbe essere osservata come supernovae di Tipo Ib o Tipo Ic.

• Tassi di Fusione e Sorgenti di GW: I modelli prevedono che solo le DNS formate attraverso il canale del nucleo di elio si fondano entro un tempo di Hubble. Di conseguenza, si inferisce che eventi di onde gravitazionali come GW170817 e GW190425 originino da questo canale. I tassi di fusione locali previsti per i modelli con calci basati sulla fisica si allineano con i vincoli di LIGO-Virgo-Kagra e del campione di DNS galattiche.

Significato e Rivendicazioni
L'articolo sostiene che i modelli evolutivi dettagliati sono essenziali per risolvere la storia della formazione delle DNS. Andando oltre le formule di fitting verso griglie auto-coerenti, gli autori dimostrano che la diversità osservata delle DNS galattiche è una conseguenza naturale dello stato evolutivo della stella donatrice all'inizio della fase CE.

Gli autori notano con modestia che le loro conclusioni si basano su confronti qualitativi a causa dell'attuale assenza di un modello di evoluzione delle pulsar all'interno di POSYDON (che terrebbe conto degli effetti di selezione e del rallentamento della rotazione). Essi riconoscono inoltre che i loro risultati sono specifici per la metallicità solare. Tuttavia, lo studio stabilisce che:

1. Una biforcazione nei canali di formazione (nucleo di He vs nucleo di C/O) è una caratteristica robusta dell'evoluzione binaria che spiega la divisione del periodo orbitale osservata.
2. Un'efficiente espulsione del CE (tramite definizioni del nucleo generose o un alto $\alpha_{CE}$) è richiesta per formare la popolazione in fusione.
3. Calci natali bassi sono una condizione necessaria per la sopravvivenza delle binarie DNS, supportando i recenti spostamenti teorici lontano dalle distribuzioni di calcio ad alta velocità.
4. Il trasferimento di massa instabile di tipo BB non produce DNS in fusione in questi modelli dettagliati, sfidando le ipotesi precedenti riguardanti la formazione di sistemi a fusione ultra-rapida.

Il lavoro fornisce una base per future analisi bayesiane quantitative una volta che i modelli di evoluzione delle pulsar saranno integrati nel framework di POSYDON.