Constraining inhomogeneities and asymmetries in SNe, FBOTs, and other high-energy transients from unresolved radio observations

Technische Samenvatting: Beperken van Inhomogeniteiten en Asymmetrieën in Hoogenergetische Transiënten vanuit Onopgeloste Radio-observaties

Probleemstelling
Hoogenergetische transiënten, waaronder supernova's (SNe), fast blue optical transients (FBOTs), tidal disruption events en gamma-ray bursts, produceren niet-thermische synchrotronstraling die de interactie van schokken met het circumstellaire medium (CSM) traceert. Hoewel Very Long Baseline Interferometry (VLBI) de bronstructuur kan onderzoeken, zijn deze bronnen doorgaat onopgelost vanwege hun afstand. Bijgevolg berust het bepalen van de geometrie en de interne structuur van het ejecta en het CSM op het interpreteren van radiospectra.

Standaard synchrotron zelfabsorptie (SSA) modellen voorspellen een laagfrequente spectrale helling van $F_\nu \propto \nu^{5/2}$ en een specifieke hoogfrequente helling bepaald door de elektronische machtswetindex $p$. Echter, observaties onthullen frequent afwijkingen: laagfrequente hellingen die minder steil zijn dan $5/2$, hoogfrequente hellingen die platter zijn dan verwacht, of bredere overgangen nabij de zelfabsorptiefrequentie ($\nu_{sa}$). Eerdere interpretaties van deze afwijkingen waren vaak fenomeenologisch en misten een rigoureus kader om de specifieke ruimtelijke verdeling van inhomogeniteiten of asymmetrieën binnen het emitterende gebied af te leiden uit onopgeloste data.

Methodologie
De auteurs stellen een algemeen, model-onafhankelijk kader voor om de structuur van het emitterende gebied te achterhalen door het geobserveerde spectrum te behandelen als een superpositie van emissies van meerdere "patches" (subregio's) met variërende fysische eigenschappen.

1. Forward Modeling van Patches: De totale geobserveerde flux $F_{\nu, \text{obs}}$ wordt gemodelleerd als de integraal van bijdragen van patches die worden gekenmerkt door hun zelfabsorptiefrequentie $\nu_{sa}$ en piekflux $F_{\nu_{sa}}$. De distributie van deze patches wordt beschreven door een functie $n(\nu_{sa}, F_{\nu_{sa}})$.
2. Formulering van het Inverse Probleem: De auteurs definiëren een gereduceerde distributiefunctie $F(\nu_{sa})$, die de totale piekflux vertegenwoordigt die door alle patches met een specifieke $\nu_{sa}$ wordt bijgedragen. Het geobserveerde spectrum is de convolutie van $F(\nu_{sa})$ met de standaard SSA-spectrale vorm.
   • Het probleem wordt geformuleerd als een invers probleem (Vergelijking 10): $F_{\nu, \text{obs}} = \int F(\nu_{sa}) f_\nu(\nu_{sa}) d\nu_{sa}$.
   • Vanwege de slecht gestelde aard van deze inversie (niet-uniciteit en instabiliteit met ruis), gebruiken de auteurs Tikhonov-regularisatie om $F(\nu_{sa})$ op te lossen.
3. Degeneratieanalyse: Het artikel analyseert rigoureus de degeneraties die inherent zijn aan het systeem. Het demonstreert dat men zonder aanvullende beperkingen (bijv. hoge-resolutie beeldvorming) niet uniek de gezamenlijke distributie van het magnetisch veld ($B$) en de schokradius ($R$) kan herstellen. Echter, de distributie $F(\nu_{sa})$ kan worden hersteld als het spectrum dicht bemonsterd is onder $\nu_{sa}$.
4. Afhankelijkheid van de Spectrale Helling: De auteurs leiden analytische relaties af die laten zien hoe de correlatie tussen $F(\nu_{sa})$ en $\nu_{sa}$ (gemodelleerd als $F(\nu_{sa}) \propto \nu_{sa}^\alpha$) de geobserveerde spectrale helling dicteert. Specifiek:
   • Als $\alpha > 3/2$, blijft de laagfrequente helling $\nu^{5/2}$, maar verbreedt de transitiezone zich.
   • Als $-(p+1)/2 < \alpha < 3/2$, verandert de laagfrequente helling naar $\nu^{\alpha+1}$.
   • Als $\alpha < -(p+1)/2$, wordt de hoogfrequente helling gewijzigd.

Belangrijkste Bijdragen en Resultaten
Het artikel past dit kader toe op twee specifieke astrofysische bronnen om de effectiviteit ervan te demonstreren:

• SN 2016coi (Stripped-Envelope SN): Gebruikmakend van multi-epoch radio-data (8,54 tot 104,25 dagen), hebben de auteurs het spectrum geïnverteerd om $F(\nu_{sa})$ te herstellen.

  • Resultaat: De inversie onthulde dat het emitterende gebied geen enkele homogene schil is. In plaats daarvan zijn meerdere patches vereist met verschillende $\nu_{sa}$-waarden.
  • Fysische Interpretatie: De data suggereren een klonterig (clumpy) densiteitsmedium. Zo wordt de emissie bij vroege tijdstippen gedomineerd door patches met magnetische velden die verschillen door factoren van $\sim 1,7$ tot $>3$, die oppervlaktegebieden beslaan die verschillen door factoren van $\sim 2,4$ tot $7$. De resultaten ondersteunen de aanwezigheid van inhomogeniteiten in de stellaire wind op een radius van $\sim 1,1 \times 10^{16}$ cm.

• AT2018cow (FBOT): Vanwege de schaarse bemonstering werd een volledige inverse oplossing als degeneratief beschouwd. In plaats daarvan gebruikten de auteurs forward modeling om de evolutie van het emitterende gebied te beperken.

  • Resultaat: De spectrale evolutie (het afvlakken van de laagfrequente helling en het steiler worden van de hoogfrequente helling) wijst op een dramatische verandering in de structuur van het emitterende gebied.
  • Fysische Interpretatie: Het emitterende gebied gaat over van een homogene, hoogmagnetische toestand op vroege tijdstippen ($t \sim 22$ dagen) naar een zeer inhomogene, patchy toestand die gedomineerd wordt door lage magnetische velden op latere tijdstippen ($t \sim 83-132$ dagen). De dekkingsfactor (covering factor) van hoogveld-regio's daalt met ordes van grootte ($\sim 100$ tot $1000$ keer). Deze evolutie wordt geïnterpreteerd als sterk bewijs voor grootschalige asymmetrie, consistent met een jet-gedreven outflow die interacteert met een dicht, asymmetrisch CSM.

Betekenis en Claims
Het artikel beweert een robuuste, model-onafhankelijke methode te bieden om de mate van inhomogeniteit en asymmetrie in onopgeloste hoogenergetische transiënten te beperken.

• Algemene Toepasbaarheid: De methode is toepasbaar op SNe, FBOTs, tidal disruption events en GRBs, mits er dicht bemonsterde spectra beschikbaar zijn.
• Oplossing van Ambiguïteit: De auteurs stellen dat afwijkingen van de standaard $\nu^{5/2}$ helling niet louter fenomeenologische anomalieën zijn, maar direct bewijs leveren voor de ruimtelijke verdeling van magnetische velden en emitterende oppervlakken.
• Correctie van Fysische Parameters: Het artikel demonstreert dat het aannemen van homogeniteit leidt tot een onderschatting van de ware schokradius (en snelheid) met factoren tot wel $\sim 2$ in zeer inhomogene bronnen. De afgeleide formalisme maakt het mogelijk om deze parameters te corrigeren.
• Beperkingen: De auteurs merken bescheiden op dat hoewel de methode de distributie van fysische parameters herstelt, zij niet tussen angulaire asymmetrieën (bijv. een jet) en radiale inhomogeniteiten (bijv. turbulentie) kan onderscheiden zonder aanvullende beperkingen (zoals VLBI of polarisatie). Bovendien wordt de inversie beperkt door het aantal observationele datapunten.

Het werk concludeert dat toekomstige radiotelescopen (bijv. SKA, ngVLA) met hogere gevoeligheid en fijnere frequentieresolutie cruciaal zullen zijn voor het verkrijgen van de dicht bemonsterde spectra die nodig zijn om deze inversietechniek volledig te benutten en de complexe structuren van transiënte emitterende regio's in kaart te brengen.