The Mass Function of Neutron Stars from Core-Collapse Supernova Simulations

Nel cuore profondo delle stelle massicce, si consuma un dramma finale e violento. Quando una stella molto più pesante del nostro Sole esaurisce il combustibile, il suo nucleo collassa sotto il proprio peso, innescando un'esplosione di supernova che può superare in luminosità un'intera galassia. In seguito, il nucleo spesso sopravvive come una stella di neutroni, un oggetto così denso che un singolo cucchiaino del suo materiale peserebbe quanto una montagna. Per decenni, gli astronomi hanno cercato di comprendere la nascita di questi residui cosmici. Hanno misurato le masse delle stelle di neutroni trovate in sistemi binari, dove due stelle orbitano l'una attorno all'altra, ma le ragioni alla base della specifica distribuzione di questi pesi sono rimaste elusive. Esiste un peso standard per una stella di neutroni neonata, o esse nascono con una grande varietà di masse? E come fa il caotico processo della morte di una stella a determinare se diventerà una stella di neutroni o se collasserà completamente in un buco nero?

Un nuovo studio condotto da un team di ricercatori mette insieme le più avanzate simulazioni al computer di queste esplosioni stellari con le più recenti osservazioni delle stelle di neutroni per rispondere a queste domande. Mappando la struttura interna della stella morente rispetto al peso del residuo che lascia dietro di sé, il team ha costruito un quadro teorico di come nascono le stelle di neutroni. Il loro lavoro suggerisce che la varietà di masse che vediamo oggi non è solo il risultato di stelle che acquistano peso più avanti nella loro vita, ma è scritta nel momento stesso della loro creazione. Le scoperte offrono una spiegazione coerente sul perché la maggior parte delle stelle di neutroni si raggruppi attorno a un peso specifico, mentre un gruppo più piccolo si estende verso masse molto più elevate, e forniscono una nuova stima di quanto spesso le stelle massicce non riescano affatto a diventare stelle di neutroni, collassando invece in buchi neri.

I ricercatori hanno iniziato esaminando l'intricata relazione tra la dimensione iniziale di una stella massiccia e la struttura del suo nucleo appena prima che esploda. Hanno utilizzato sofisticati modelli computerizzati tridimensionali che simulano la complessa fisica di una supernova, tracciando il comportamento del nucleo della stella negli ultimi momenti della sua vita. Queste simulazioni hanno rivelato un modello chiaro: la "compattezza" interna del nucleo della stella — una misura di quanto la materia sia densamente impacchettata in un punto specifico all'interno — determina direttamente la massa della stella di neutroni che si forma. Un nucleo meno compatto tende a lasciare dietro di sé una stella di neutroni più leggera, mentre un nucleo più compatto produce una stella di neutroni più pesante. Questa connessione ha permesso al team di prendere la distribuzione nota delle stelle massicce nell'universo e prevedere la distribuzione delle masse delle stelle di neutroni al momento della loro nascita.

Quando hanno confrontato questa distribuzione di nascita teorica con le masse reali delle stelle di neutroni osservate nel cielo, è emersa una discrepanza. Le simulazioni prevedevano che le stelle di neutroni neonate più comuni dovrebbero essere leggermente più leggere della massa media osservata, con un secondo gruppo di stelle più pesanti. Tuttavia, i dati osservati mostravano un picco a un peso leggermente superiore. I ricercatori si sono resi conto che questa differenza poteva essere spiegata da due fattori principali: le stelle di neutroni che acquistano massa dopo la loro nascita e il modo in cui la loro nascita influisce sulla permanenza in un sistema binario.

Le stelle di neutroni nei sistemi binari spesso attirano materiale dalle loro stelle compagne, un processo chiamato accrezione, che aggiunge peso alla stella di neutroni nel tempo. Il team ha applicato correzioni ai dati osservati per stimare quali sarebbero state le masse alla nascita, eliminando il peso guadagnato dalle loro compagne. Questo aggiustamento ha spostato il picco osservato verso il basso, allineandolo molto più strettamente alla previsione teorica. Inoltre, lo studio ha evidenziato il ruolo dei "calci" (kicks) impartiti alle stelle di neutroni alla nascita. Proprio come un razzo è spinto in avanti dal suo scarico, una stella di neutroni riceve una spinta improvvisa nella direzione opposta all'esplosione. Le simulazioni suggeriscono che le stelle di neutroni più pesanti ricevano in media calci più forti. Questi potenti calci possono spezzare il legame gravitazionale che tiene unito il sistema binario, espellendo le stelle di neutroni più pesanti nel vuoto, dove sono più difficili da trovare. Questa rimozione selettiva significa che i sistemi binari che osserviamo sono distorti verso le stelle di neutroni più leggere che hanno ricevuto calci più gentili e sono rimaste in orbita.

Combinando gli effetti della formazione di buchi neri, dell'accrezione di massa e di questi calci di nascita, il team ha creato un modello raffinato che si adatta sorprendentemente bene ai dati osservati. Il loro modello suggerisce che l'universo produca uno spettro ampio di masse di stelle di neutroni alla nascita, che varia da circa 1,2 a 2,0 volte la massa del nostro Sole. La massa di nascita più comune è di circa 1,35 volte la massa del Sole, il che spiega il prominente picco visto nelle osservazioni. La coda più pesante della distribuzione, che si estende fino a 2,0 masse solari, è un risultato naturale delle diverse strutture delle stelle progenitrici, piuttosto che solo del risultato dell'accrezione successiva.

Lo studio ha anche fornito una nuova stima della frazione di stelle massicce che terminano le loro vite come buchi neri invece che come stelle di neutroni. Analizzando quali tipi di stelle nelle loro simulazioni non sono riuscite a produrre una stella di neutroni, i ricercatori hanno calcolato che circa il 21 percento delle stelle massicce nell'intervallo di massa pertinente collassa direttamente in buchi neri. Questo numero non è una regola semplice basata sulla dimensione iniziale della stella; piuttosto, dipende dalla struttura complessa e non lineare del nucleo della stella, il che significa che stelle di dimensioni simili possono avere destini molto diversi.

Questo lavoro rappresenta un passo significativo nel connettere la teoria di come muoiano le stelle con la realtà di ciò che osserviamo. Supera la vecchia idea secondo cui tutte le stelle di neutroni nascono con quasi la stessa massa, sostituendola con un quadro dinamico in cui la struttura iniziale della stella detta una vasta gamma di risultati. L'accordo tra le complesse simulazioni e i dati osservati suggerisce che la nostra comprensione delle esplosioni di supernova stia maturando. Sebbene vi siano ancora incertezze, come l'influenza precisa della composizione chimica della stella o i dettagli delle interazioni binarie, lo studio dimostra che la moderna teoria delle supernovae può ora spiegare la popolazione di oggetti compatti che vediamo nel cielo. Apre una nuova era in cui la nascita di questi oggetti estremi può essere ricondotta alle specifiche condizioni delle stelle che li hanno creati, trasformando una collezione di misurazioni in una storia coerente di morte e rinascita stellare.