Electron impact excitation of Te IV and V and Level Resolved R-matrix Photoionization of Te I - IV with application to modelling of AT2017gfo

Questo studio presenta dati atomici calcolati con il metodo R-matrix per ioni di tellurio (Te I-IV) e risultati specifici per Te IV e V, utilizzandoli per modellare spettri di kilonove e indagare il contributo di Te IV alla caratteristica emissione a 1,08 µm osservata in AT2017gfo.

L'essenziale

Immagina di essere un detective che sta cercando di risolvere il mistero di un'esplosione cosmica gigantesca, avvenuta miliardi di anni fa. Questo "crimine" è stato un kilonova, un evento creato quando due stelle di neutroni (le spoglie ultra-dense di stelle morte) si sono scontrate.

Quando queste stelle si scontrano, espellono una nuvola di detriti incredibilmente calda e luminosa. Per capire cosa c'è dentro questa nuvola, gli astronomi guardano la luce che ne esce, come se fosse un codice a barre. Ogni elemento chimico lascia un'impronta digitale unica nella luce.

Il problema? Per leggere questo codice a barre, abbiamo bisogno di una "chiave" molto precisa: i dati su come questi atomi si comportano quando vengono colpiti da elettroni o da luce. Fino a poco tempo fa, per elementi pesanti e rari come il Tellurio (Te), usavamo delle "chiavi" fatte in casa, approssimative, come se provassimo a chiudere una serratura complessa con un pezzo di legno.

Ecco cosa hanno fatto gli autori di questo studio, in termini semplici:

1. La Fabbrica di Chiavi Precise (I Dati Atomici)

Gli scienziati hanno costruito delle "chiavi" matematiche perfette per il Tellurio. Hanno usato un metodo sofisticato chiamato metodo R-matrix.

• L'analogia: Immagina di voler prevedere come si comporta una palla da biliardo quando colpisce un'altra palla. Se usi una formula semplice, potresti sbagliare di poco. Se usi un supercomputer che simula ogni singola forza tra le molecole, ottieni una previsione perfetta.
• Hanno calcolato esattamente come gli ioni di Tellurio (atomi che hanno perso alcuni elettroni, come Te IV e Te V) si eccitano quando vengono colpiti da elettroni e come si ionizzano quando vengono colpiti dalla luce.

2. Il Mistero della Luce Blu (La riga a 1,08 micron)

C'è un particolare mistero nella luce di una famosa kilonova chiamata AT2017gfo. Circa una settimana dopo l'esplosione, gli astronomi hanno visto un picco di luce molto luminoso a una lunghezza d'onda specifica (1,08 micron, che è nel vicino infrarosso).

• Il sospettato n. 1: Per molto tempo si è pensato che fosse lo Stronzio (Sr), un elemento più leggero. È come se avessimo accusato un ladro comune per un crimine molto specifico.
• Il nuovo sospettato: Gli autori di questo studio dicono: "Aspetta, guarda qui!". Hanno scoperto che il Tellurio IV (una versione del Tellurio con tre elettroni in meno) ha una "impronta digitale" che corrisponde perfettamente a quel picco di luce.
• L'analogia: È come se avessimo trovato un'impronta di scarpa che corrisponde esattamente a quella lasciata sulla scena del crimine, mentre quella dello Stronzio era solo simile ma non perfetta.

3. È possibile che il Tellurio sia lì? (Il problema della temperatura)

C'è un ostacolo. Per avere il Tellurio IV, l'ambiente deve essere molto caldo e "ionizzato" (gli atomi devono perdere elettroni). Ma la kilonova, in quella fase, sembra essere piuttosto fredda.

• La soluzione: Gli scienziati spiegano che, anche se la temperatura media è bassa, ci sono particelle radioattive che decadono e creano "elettroni non termici" (elettroni veloci e arrabbiati). Questi elettroni possono strappare via gli elettroni dal Tellurio anche se l'aria circostante è fresca.
• L'analogia: Immagina di essere in una stanza fredda (la kilonova), ma qualcuno sta lanciando proiettili caldi (elettroni radioattivi) contro una bottiglia d'acqua (l'atomo di Tellurio). Anche se la stanza è fredda, la bottiglia si rompe (si ionizza) a causa dei proiettili.

4. Perché è importante?

Prima di questo studio, se volevamo simulare queste esplosioni al computer, dovevamo usare dati "fatti a mano" per il Tellurio, che potevano essere sbagliati.

• Il risultato: Ora abbiamo dati precisi, calcolati con la massima cura possibile. Questo permette agli astronomi di dire con più sicurezza: "Sì, quel bagliore è Tellurio, non Stronzio".
• Il futuro: Questo aiuta a capire meglio come si formano gli elementi pesanti nell'universo (come l'oro e l'uranio) durante queste esplosioni. È come passare da una mappa disegnata a mano con errori a una mappa satellitare ad alta definizione.

In sintesi

Gli autori hanno creato una mappa di precisione per il comportamento del Tellurio nello spazio. Hanno scoperto che questo elemento potrebbe essere il "colpevole" principale di un particolare bagliore osservato dopo una collisione di stelle di neutroni, risolvendo un mistero che gli astronomi stavano cercando di decifrare da anni. Hanno dimostrato che, con le giuste condizioni (anche un po' di calore extra dai decadimenti radioattivi), il Tellurio può essere lì e brillare proprio come previsto.

È un po' come se avessimo finalmente trovato la lente giusta per guardare attraverso la nebbia cosmica e vedere chiaramente cosa c'è davvero dentro.

Sommario tecnico

Titolo: Eccitazione per impatto elettronico di Te IV e V e fotoionizzazione risolta per livelli di Te I-IV con applicazione alla modellazione di AT2017gfo

1. Il Problema

La modellazione spettroscopica delle kilonove (KNe), eventi risultanti dalla fusione di stelle di neutroni come AT2017gfo, richiede grandi quantità di dati atomici precisi per specie pesanti a basso grado di ionizzazione (neutri, singolarmente e doppiamente ionizzate). Attualmente, molti studi si basano su risultati idrogenoidi approssimati o formule semi-empiriche (es. van Regemorter, Axelrod) per i processi di eccitazione e fotoionizzazione. Questi approcci, sebbene utili per la scalabilità, introducono incertezze significative nella sintesi spettrale, specialmente quando si trattano elementi complessi come il Tellurio (Te, Z=52), che si trova al secondo picco del processo di cattura neutronica rapida (r-process) ed è previsto essere abbondante nel materiale espulso. Inoltre, c'è una carenza di dati calcolati rigorosamente per gli stati di ionizzazione più elevati (Te IV e V) e per la fotoionizzazione risolta per livelli, necessari per modellare correttamente le fasi iniziali (dominate dalla fotoionizzazione) e le fasi successive (dominate dalle collisioni) delle kilonove.

2. Metodologia

Gli autori hanno adottato un approccio basato sulla teoria R-matrix, considerata lo standard aureo per il calcolo di processi atomici in plasmi astrofisici, combinato con il metodo Multi-Configuration-Dirac-Hartree-Fock (MCDHF).

• Struttura Atomica: Sono stati calcolati i modelli strutturali per Te IV e Te V utilizzando il pacchetto grasp0. Sono state incluse configurazioni di valenza non relativistiche selezionate attentamente per bilanciare accuratezza e complessità computazionale (21 configurazioni per Te IV, 27 per Te V). I livelli energetici sono stati confrontati con dati sperimentali e calcoli teorici precedenti (es. grasp2k, fac).
• Eccitazione per Impatto Elettronico: Utilizzando i modelli strutturali, sono stati calcolati i parametri di collisione (forze di collisione $\Omega$) per Te IV e V. I calcoli sono stati eseguiti con il codice darc, mantenendo 100 livelli target e utilizzando un'espansione di onde parziali fino a $2J=62$ (o superiore con regole di "top-up"). Le forze di collisione sono state mediate termicamente secondo una distribuzione di Maxwell-Boltzmann per ottenere le velocità di eccitazione/de-eccitazione efficaci ($\Upsilon$).
• Fotoionizzazione: Sono stati calcolati i sezioni d'urto di fotoionizzazione risolti per livello per gli ioni Te I, II, III e IV. Questi calcoli utilizzano i modelli target di Te II, III, IV e V rispettivamente, trattando esplicitamente le risonanze che i metodi approssimati trascurano.
• Modellazione Collisionale-Radiativa: I nuovi dati sono stati integrati nel codice ColRadPy per generare spettri sintetici e confrontarli con le osservazioni di AT2017gfo, investigando in particolare l'origine della caratteristica emissiva a 1.08 $\mu$m.

3. Contributi Chiave

• Nuovi Dati Atomici: Fornitura di dati calcolati con R-matrix per la prima volta per gli ioni Te IV e V (eccitazione) e per la fotoionizzazione risolta per livelli di Te I-IV.
• Validazione dei Modelli: Confronto dettagliato dei livelli energetici e dei coefficienti di Einstein ($A$-values) con dati sperimentali e codici indipendenti (come fac), dimostrando un buon accordo, specialmente per le transizioni proibite e intercombinate.
• Applicazione a AT2017gfo: Investigazione specifica del contributo di Te IV alla caratteristica emissiva a 1.08 $\mu$m osservata nelle fasi centrali (circa 7 giorni dopo l'esplosione) di AT2017gfo, un'alternativa o complemento all'identificazione precedente basata su Sr II.
• Disponibilità dei Dati: I dati sono resi disponibili nel formato standard ADAS (OPEN-ADAS) e nel materiale supplementare, permettendo l'uso diretto nei codici di sintesi spettrale su larga scala.

4. Risultati Principali

• Struttura Atomica: I livelli calcolati per Te IV e V mostrano un accordo eccellente con i dati sperimentali per i livelli più bassi (differenze $\le 0.09$ Ry). Per Te IV, il coefficiente di transizione proibita per il decadimento $5s^2 5p \ ^2P^o_{3/2} \to \ ^2P^o_{1/2}$ (lunghezza d'onda 1.08 $\mu$m) è stato calcolato come $7.03 \ s^{-1}$, in ottimo accordo con il valore sperimentale di $7.09 \ s^{-1}$.
• Eccitazione: Sono stati forniti i parametri di collisione termicamente mediati per tutte le transizioni tra i livelli inclusi. I dati mostrano che le forze di collisione convergono bene grazie all'alto numero di onde parziali e alla corretta gestione del comportamento asintotico.
• Fotoionizzazione: Le sezioni d'urto calcolate mostrano risonanze ampie (es. intorno a 3.0-3.05 Ry) attribuite a transizioni di eccitazione del nucleo interno (shell 5s). I dati coprono un intervallo di energie da 0 a 6 Ry.
• Modellazione di AT2017gfo:
  • L'analisi suggerisce che Te IV potrebbe essere il principale responsabile della forte emissione a 1.08 $\mu$m osservata tra i 6 e gli 8 giorni dopo l'esplosione.
  • Per riprodurre la luminosità osservata ($\sim 2.0 \times 10^{39}$ erg/s) a una temperatura di 3000 K e densità di $2.0 \times 10^7$ cm$^{-3}$, è necessaria una massa di Te IV di circa $2.5 \times 10^{-3} M_{\odot}$.
  • A differenza di Sr II, che produce caratteristiche P-Cygni e linee spurie verso il blu, Te IV produce una linea di emissione "pulita" a 1.08 $\mu$m senza contaminanti significativi, anche a temperature più elevate (fino a 10.000 K).
  • Nota Critica: La presenza di Te IV (stato triplamente ionizzato) a temperature relativamente basse (3000 K) è improbabile in condizioni di equilibrio termodinamico locale (LTE). Tuttavia, la presenza di elettroni non termici derivanti dal decadimento radioattivo (come simulato da Pognan et al.) può mantenere popolati stati di ionizzazione più alti, rendendo plausibile l'identificazione.

5. Significato e Prospettive Future

Questo lavoro rappresenta un passo fondamentale verso la sostituzione di dati atomici approssimati con dati di alta precisione basati sulla teoria R-matrix per gli elementi pesanti nelle kilonove.

• Impatto Scientifico: Migliora l'affidabilità della sintesi spettrale, permettendo di vincolare meglio la composizione chimica, la massa espulsa e le condizioni fisiche (temperatura, densità) delle kilonove.
• Identificazione di Te IV: Offre una spiegazione plausibile per l'evoluzione temporale della caratteristica a 1.08 $\mu$m, suggerendo un ruolo di Te IV nelle fasi intermedie, mentre Sr II domina la fase iniziale (P-Cygni).
• Futuro: I dati forniti sono essenziali per modelli radiativi più sofisticati che includano bilanci di ionizzazione non termici. Il lavoro apre la strada a studi simili per altri elementi del processo r e sottolinea la necessità di includere effetti di ionizzazione non termica per spiegare la presenza di stati di ionizzazione elevati in ambienti freddi come le kilonove.