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

Deze studie levert nauwkeurige atomaire gegevens voor telluur-ionen, berekend met de R-matrix- en MCDHF-methoden, om de spectra van kilonova's zoals AT2017gfo beter te modelleren en de mogelijke bijdrage van Te IV aan de emissie bij 1,08 μm te onderzoeken.

De kern

Titel: Het Verborgen Verhaal van Telluur in de Kosmische Explosies

Stel je voor dat het heelal een gigantische, donkere bioscoopzaal is. Soms, heel zelden, gebeurt er iets spectaculair: twee zware sterren botsen tegen elkaar. Dit is een kilonova. Het is als een kosmische vuurwerkshow die zo helder is dat hij tijdelijk de hele melkweg verlicht.

In 2017 zagen astronomen zo'n explosie, genaamd AT2017gfo. Maar toen ze naar de 'film' keken (het licht dat de explosie uitzond), zagen ze iets vreemds. Er was een specifiek lichtflitsje op een bepaalde golflengte (een soort kleur) dat ze niet helemaal konden verklaren. Het leek op een spook in de film.

De wetenschappers in dit artikel, Leo en zijn team, hebben zich afgevraagd: "Wie is die spook?"

1. De Puzzelstukjes (Atomaire Gegevens)

Om te begrijpen wat er in zo'n explosie gebeurt, moeten we weten hoe atomen zich gedragen. Atomen zijn als kleine LEGO-blokjes. Als ze botsen of licht opvangen, kunnen ze van kleur veranderen (ze stralen licht uit).

Voor de meeste atomen hebben we al een goede handleiding. Maar voor zware, zeldzame elementen zoals Telluur (een metaal dat we op aarde gebruiken in elektronica), hadden we geen goede handleiding. De oude handleidingen waren als het schatten van een afstand met je hand: "ongeveer zo ver". Dat is niet nauwkeurig genoeg voor de complexe natuurkunde van een sterrenexplosie.

Deze wetenschappers hebben dus een nieuwe, super-nauwkeurige handleiding geschreven voor Telluur. Ze hebben berekend hoe Telluur-atomen botsen met elektronen en hoe ze licht opnemen. Ze hebben dit gedaan voor verschillende versies van Telluur (waarvan sommige elektronen kwijt zijn geraakt, net als een auto die zijn deuren kwijtraakt).

2. De Detectivewerk (De 1,08 micron mysterie)

In de explosie AT2017gfo was er een mysterieus lichtflitsje bij een golflengte van 1,08 micron (een kleur die we niet kunnen zien met het blote oog, maar wel met speciale camera's).

• De oude theorie: Iedereen dacht dat dit licht kwam van Strontium (een ander element), net als een bekende acteur die een rol speelt.
• De nieuwe theorie: De onderzoekers zeggen: "Wacht even. Als we kijken naar onze nieuwe, super-nauwkeurige handleiding voor Telluur, past Telluur (in een specifieke, zware versie genaamd Te IV) perfect bij dit lichtflitsje!"

Het is alsof je een spook in de bioscoop zag en dacht dat het een bekende acteur was, maar toen je de camera vergrootte, zag je dat het eigenlijk een heel goed vermomde andere acteur was.

3. De Uitdaging: Is het mogelijk?

Er is een probleem. Telluur in deze vorm (Te IV) is als een zeer energieke versie van het element. Normaal gesproken heeft het heel veel hitte nodig om zo energiek te worden. De explosie was op dat moment misschien niet heet genoeg.

Maar, de onderzoekers zeggen: "Misschien is er een onzichtbare kracht die het element toch zo heet maakt." Denk aan radioactief verval (het 'verval' van atomen in de explosie) als een onzichtbare verwarming die de atomen extra energie geeft. Als dat zo is, dan is het heel goed mogelijk dat Telluur die mysterieuze lichtflits veroorzaakt.

4. Waarom is dit belangrijk?

Stel je voor dat je een recept probeert te volgen om een taart te bakken, maar je hebt geen exacte maten voor de suiker. Je taart wordt misschien wel goed, maar je weet niet precies waarom.

• Vroeger: Wetenschappers gebruikten 'geschatte' maten (ongeveer) om de spectra van deze explosies te modelleren.
• Nu: Deze paper geeft de exacte maten voor Telluur.

Dit betekent dat we in de toekomst veel beter kunnen begrijpen:

1. Wat er precies gebeurt in deze kosmische explosies.
2. Hoe zware elementen (zoals goud en platina, die ontstaan bij deze explosies) in het heelal worden verspreid.
3. Of we de 'spookachtige' lichtflitsjes in de toekomst beter kunnen voorspellen.

Samenvatting in één zin

Deze wetenschappers hebben een nieuwe, super-nauwkeurige 'gebruiksaanwijzing' gemaakt voor het element Telluur, en daarmee bewezen dat Telluur waarschijnlijk de schuldige is van een mysterieus lichtflitsje in een sterrenexplosie, wat ons helpt om de geboorte van zware elementen in het heelal beter te begrijpen.

Het is alsof ze de sleutel hebben gevonden om een deur te openen die we al jaren probeerden te forceren, en nu zien we eindelijk wat er achter zit.

Technische samenvatting

Probleemstelling

De spectroscopische modellering van kilonovae (zoals AT2017gfo en AT2023vfi) vereist uitgebreide atomaire data voor zware elementen, met name voor laag geïoniseerde soorten (neutraal, enkel en dubbel geïoniseerd) die ontstaan tijdens de $r$-proces nucleosynthese. Helaas is de beschikbare data voor telluur (Te, $Z=52$), een element dat zich bij de tweede piek van de $r$-proces overvloed bevindt, vaak beperkt of gebaseerd op benaderingen.

• Huidige beperkingen: Veel bestaande modellen gebruiken semi-empirische formules (zoals die van van Regemorter en Axelrod) voor excitatieprocessen of waterstofachtige benaderingen voor foto-ionisatie.
• Specifiek probleem: Er is een gebrek aan nauwkeurige $R$-matrix data voor hogere ionisatietoestanden van telluur (Te IV en V), die mogelijk een rol spelen in de late fasen van kilonovae. Bovendien is er onzekerheid over de oorsprong van de emissiepiek bij 1,08 $\mu$m in AT2017gfo, die eerder toegeschreven werd aan Sr II, maar mogelijk een bijdrage bevat van Te IV.

Methodologie

De auteurs hebben een robuuste berekeningsframework toegepast om atomaire data te genereren voor Te I tot en met Te V:

1. Atomaire Structuur (MCDHF):

   • Gebruik van de Multi-Configuration-Dirac-Hartree-Fock (MCDHF) methode via de grasp0 code.
   • Er werden geoptimaliseerde single-electron orbitalen gegenereerd op basis van configuratie-state-functies (CSFs).
   • Voor Te IV en V werden specifieke configuraties geselecteerd om een representatief model te creëren dat klein genoeg is voor grote $R$-matrix berekeningen, maar groot genoeg om de belangrijkste energieniveaus en stralingsovergangen nauwkeurig weer te geven.
   • Einstein A-coëfficiënten werden berekend en waar mogelijk gecorrigeerd met spectroscopische golflengten.

2. Elektron-impact Excitatie ($R$-matrix):

   • Berekening van botsingssterktes ($\Omega$) en effectieve botsingssterktes ($\Upsilon$) voor Te IV en V.
   • Gebruik van de $R$-matrix methode met een close-coupling expansie (100 doelniveaus).
   • Inclusie van resonantiestructuur door fijne energie-meshes en een "top-up" regel voor hoge partiële golven.
   • Thermische middeling volgens een Maxwelliaanse verdeling voor verschillende elektrontemperaturen.

3. Foto-ionisatie:

   • Berekening van niveau-opgeloste foto-ionisatie doorsneden voor Te I, II, III en IV.
   • Gebruik van de darc code (een implementatie van de $R$-matrix methode voor foto-ionisatie).
   • De berekeningen omvatten resonanties die in benaderende methoden vaak worden gemist.

4. Stralingsmodelling (Collisional-Radiative):

   • De berekende excitatie- en stralingsdata werden ingevoerd in het ColRadPy pakket.
   • Synthetische spectra werden gegenereerd en vergeleken met waarnemingen van AT2017gfo, specifiek gericht op de emissie bij 1,08 $\mu$m.

Belangrijkste Bijdragen

• Nieuwe Atomaire Data: Voor het eerst zijn nauwkeurige $R$-matrix data geproduceerd voor Te IV en V (excitatie) en voor Te I-IV (foto-ionisatie).
• Validatie: De berekende energieniveaus en stralingsovergangen (A-waarden) tonen goede overeenkomst met experimentele waarden en eerdere theoretische berekeningen (o.a. met de fac code en grasp2k).
• Toepassing op AT2017gfo: Een gedetailleerd onderzoek naar de mogelijke bijdrage van Te IV aan de 1,08 $\mu$m emissiepiek, naast de reeds bekende Sr II P-Cygni kenmerken.
• Open Data: De effectieve botsingssterktes en foto-ionisatie doorsneden zijn beschikbaar gesteld via OPEN-ADAS en supplementair materiaal, wat de standaard in de gemeenschap verhoogt.

Resultaten

1. Atomaire Structuur:

   • Voor Te IV werden 21 configuraties gebruikt; voor Te V werden 27 configuraties gebruikt.
   • De berekende energieniveaus voor Te IV vertonen een afwijking van maximaal 0,09 Ry ten opzichte van experimentele waarden.
   • De verboden overgang $5s^2 5p \ ^2P^o_{3/2} \rightarrow \ ^2P^o_{1/2}$ in Te IV (golflengte 1,08 $\mu$m) heeft een berekende A-waarde van 7,03 s$^{-1}$, wat uitstekend overeenkomt met de literatuurwaarde van 7,09 s$^{-1}$.

2. Excitatie en Foto-ionisatie:

   • Effectieve botsingssterktes zijn berekend voor een breed temperatuurbereik.
   • Foto-ionisatie doorsneden tonen duidelijke resonantiestructuren, met name een brede resonantie rond 3,0-3,05 Ry die correleert met foto-excitatie van de 5s innerschil.

3. Modellering van AT2017gfo:

   • 1,08 $\mu$m Emissie: De synthetische spectra suggereren dat Te IV een significante bijdrage kan leveren aan de emissie-dominante fase van de 1,08 $\mu$m piek rond dag 7 na de explosie.
   • Massa en Temperatuur: Om de waargenomen lichtkracht ($\sim 2,0 \times 10^{39}$ erg/s) te reproduceren, is een Te IV-massa van ongeveer $2,5 \times 10^{-3} M_{\odot}$ nodig bij een temperatuur van 3000 K en een dichtheid van $2,0 \times 10^7$ cm$^{-3}$.
   • Ionisatieprobleem: In lokale thermodynamische evenwicht (LTE) is de aanwezigheid van Te IV bij deze temperaturen onwaarschijnlijk (vereist ~8500 K). Echter, onder Niet-LTE (NLTE) omstandigheden, waarbij niet-thermische ionisatie door radioactief verval een rol speelt, is de aanwezigheid van Te IV haalbaar.
   • Vergelijking met Sr II: Terwijl Sr II waarschijnlijk de P-Cygni structuur in de vroege fasen veroorzaakt, kan Te IV de zuivere emissie in de latere fasen verklaren zonder de "vervuiling" van andere emissielijnen die bij Sr II modellen vaak optreedt.

Significantie

Deze studie is van cruciaal belang voor de astrofysica van zware elementen en kilonovae:

• Verbeterde Modellering: Het vervangt onnauwkeurige benaderingen door gepreciseerde $R$-matrix data, wat essentieel is voor het interpreteren van waarnemingen van zware elementen in neutronenster-samensmeltingen.
• Nieuwe Interpretatie: Het biedt een plausibel alternatief (of aanvulling) voor de interpretatie van de 1,08 $\mu$m emissie in AT2017gfo, wat de chemische samenstelling van de ejecta beter kan verklaren.
• Toekomstgericht: De data vormt de basis voor toekomstige radiative transfer simulaties die rekening houden met niet-thermische ionisatieprocessen, wat essentieel is voor het begrijpen van de evolutie van kilonovae.
• Uniek: Dit is de eerste $R$-matrix foto-ionisatie berekening specifiek gericht op Te-ionen en hun toepassing op kilonovae-modellering.