Non-LTE Ionization Modeling for Helium and Strontium in Neutron Star Merger Ejecta

Titolo: Modellazione dell'ionizzazione non-LTE per Elio e Stronzio nei ejecta di fusioni di stelle di neutroni

Autore: Koya Chiba et al.
Data: 8 aprile 2026 (Draft)

---

1. Il Problema Scientifico

Le fusioni di stelle di neutroni binarie (BNS) producono "kilonove", emissioni elettromagnetiche alimentate dal decadimento radioattivo di nuclei sintetizzati tramite il processo rapido di cattura neutronica (r-processo). Lo spettro della kilonova AT2017gfo, associata all'evento gravitazionale GW170817, mostra una forte caratteristica di assorbimento intorno a 1 µm nella fase iniziale (primi giorni).

Fino ad ora, l'origine di questa caratteristica è stata attribuita principalmente allo Stronzio (Sr II) o, alternativamente, all'Elio (He I). Tuttavia, esistono due limitazioni fondamentali negli studi precedenti:

1. Assunzione di LTE: La maggior parte dei modelli per la fase iniziale ha assunto l'equilibrio termodinamico locale (LTE). Tuttavia, negli ejecta delle kilonove, la densità è bassa e l'ionizzazione da parte di elettroni ad alta energia (provenienti dal decadimento radioattivo) è significativa, rendendo l'assunzione LTE inadeguata.
2. Mancanza di modelli congiunti: Non esistevano modelli coerenti che considerassero simultaneamente Elio e Stronzio. Studi precedenti hanno trattato gli elementi separatamente o hanno fissato le abbondanze senza esplorare l'intero spazio dei parametri, ignorando l'interazione tra i due elementi e gli effetti non-LTE.

2. Metodologia

Gli autori hanno sviluppato un modello di ionizzazione non-LTE (Non-Local Thermodynamic Equilibrium) auto-consistente per Elio e Stronzio, tenendo conto dell'ionizzazione indotta da elettroni ad alta energia.

• Ionizzazione Non-Termica: Il modello calcola il tasso di ionizzazione ($\Gamma_i$) dovuto agli elettroni non termici generati dal decadimento radioattivo ($\dot{q}$), utilizzando un potenziale di ionizzazione efficace ($w_i$). Questo è cruciale perché gli elettroni ad alta energia (keV-MeV) ionizzano gli atomi molto più efficientemente della radiazione termica (0.1-1 eV).
• Modellazione dell'Elio (He): È stato utilizzato un modello atomico dettagliato con 21 stati (19 stati legati di He I fino a $n=4$, più stati fondamentali di He II e He III). Sono stati inclusi processi radiativi e collisionali, nonché l'ionizzazione non termica. La transizione responsabile della caratteristica a 1 µm (1.08 µm) avviene tra stati altamente eccitati (2 $^3$S – 2 $^3$P) che non sono popolati efficacemente in condizioni LTE.
• Modellazione dello Stronzio (Sr): A causa della complessità degli elementi pesanti, è stato adottato un modello semplificato che considera gli stati di ionizzazione fino a Sr XX (con focus su Sr VII). Il modello risolve l'equazione di bilancio per gli stati di ionizzazione, includendo l'ionizzazione termica (fotoionizzazione) e non termica. Si assume una distribuzione di Boltzmann per i livelli eccitati di Sr II.
• Confronto con i Dati Osservativi: I modelli sono stati applicati agli spettri osservati di AT2017gfo a quattro epoche diverse (1.43, 2.42, 3.41, 4.40 giorni). È stato utilizzato un profilo di fitting spettrale con un profilo P-Cygni e un corpo nero relativistico per stimare le condizioni fisiche (temperatura, velocità, profondità ottica $\tau$) appena sopra la fotosfera.

3. Contributi Chiave

• Sviluppo di un modello non-LTE auto-consistente: Per la prima volta, un modello calcola simultaneamente le popolazioni dei livelli di Elio e gli stati di ionizzazione dello Stronzio, tenendo conto dell'ionizzazione da elettroni non termici in un ambiente di kilonova.
• Quantificazione degli effetti non-LTE: Dimostrazione che l'assunzione LTE porta a sottostimare drasticamente le abbondanze richieste per spiegare le caratteristiche spettrali, a causa della sovr-ionizzazione degli elementi.
• Vincoli sulle abbondanze: Determinazione delle frazioni di massa necessarie di He e Sr per riprodurre la caratteristica a 1 µm, esplorando l'intero spazio dei parametri possibili.

4. Risultati Principali

Effetti Non-LTE

• Elio: La popolazione degli stati eccitati di He I (necessari per l'assorbimento a 1.08 µm) è fortemente dipendente dalla temperatura e dall'ionizzazione non termica. In condizioni non-LTE, la popolazione di questi stati è significativamente diversa rispetto al caso LTE.
• Stronzio: L'ionizzazione non termica porta a una sovr-ionizzazione dello Sr II (trasformandolo in stati di ionizzazione più alti come Sr III, IV, ecc.). Di conseguenza, la frazione di Sr II disponibile per l'assorbimento è ridotta di diversi ordini di grandezza rispetto alle previsioni LTE. Per ottenere la stessa profondità della linea di assorbimento, è necessaria una massa di Stronzio molto maggiore nel caso non-LTE.

Vincoli sulle Abbondanze

Per riprodurre la caratteristica a 1 µm osservata negli ejecta che si muovono a $v \sim 0.15 c$:

• Caso Elio: È necessaria una frazione di massa di Elio di circa 1% ($X_{He} \sim 10^{-2}$).
• Caso Stronzio: È necessaria una frazione di massa di Stronzio compresa tra 1% e 10% ($X_{Sr} \sim 10^{-2} - 10^{-1}$).
• Confronto con l'Abbondanza Solare: La frazione di massa di Stronzio richiesta (1-10%) è in ottimo accordo con l'abbondanza di r-processo nel Sole, supportando l'ipotesi che le fusioni di stelle di neutroni siano la fonte principale degli elementi r-processo nell'universo.

Condizioni di Nucleosintesi

Confrontando i vincoli osservativi con calcoli di nucleosintesi parametrici:

• Le condizioni migliori per spiegare i dati sono bassa frazione di elettroni ($Y_e \lesssim 0.35$) e bassa entropia ($s \lesssim 30 k_B$/nucleone).
• Scenario $Y_e \lesssim 0.15$: Se la frazione di elettroni è molto bassa, la caratteristica a 1 µm è probabilmente dovuta all'Elio prodotto dal decadimento $\alpha$ di nuclei trans-Piombo. Questo fornirebbe una "firma indiretta" della produzione di elementi oltre il terzo picco del r-processo.
• Scenario $Y_e \sim 0.15 - 0.35$: La caratteristica è probabilmente dovuta allo Stronzio.

Implicazioni per il Meccanismo di Eiezione

Il confronto con simulazioni idrodinamiche viscose (modello DD2-135-135, che prevede un residuo di stella di neutroni a lunga vita) mostra che tali modelli tendono a sovraprodurre sia Elio che Stronzio (a causa di $Y_e$ ed entropia troppo alti negli ejecta post-fusione).

• Questo suggerisce che, per GW170817, l'eiezione di massa significativa deve avvenire in condizioni di bassa $Y_e$ e bassa entropia, o in un tempo scala più breve rispetto a quanto previsto dalle simulazioni viscose attuali, prima che la frazione di elettroni aumenti troppo a causa del raffreddamento e delle interazioni con i neutrini.

5. Significato e Conclusioni

Questo studio dimostra che la modellazione spettrale non-LTE è essenziale per interpretare correttamente le kilonove. Ignorare l'ionizzazione non termica porta a conclusioni errate sulle abbondanze chimiche e, di conseguenza, sulle condizioni fisiche della nucleosintesi e sui meccanismi di espulsione della materia.

I risultati supportano fortemente il ruolo delle fusioni di stelle di neutroni come siti principali del processo r, con abbondanze di Stronzio coerenti con quelle solari. Inoltre, offrono nuove prospettive sulla fisica del residuo della fusione (durata di vita della stella di neutroni) e sulla possibile produzione di elementi super-pesanti (oltre il terzo picco del r-processo) attraverso il canale dell'Elio.

Limitazioni e Lavori Futuri:
Il modello attuale utilizza alcune semplificazioni, come la trattazione approssimata della ricombinazione dielettrica per gli elementi pesanti e l'assunzione di uno stato stazionario. Futuri lavori dovranno includere modelli non-LTE più realistici che tengano conto degli effetti temporali e delle interazioni radiative più complesse per risolvere la degenerazione tra i contributi di Elio e Stronzio.