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

Sintesi Tecnica: Vincolare le Inomogeneità e le Asimmetrie nei Transienti ad Alta Energia da Osservazioni Radio Non Risolte

Problematica
I transienti ad alta energia, tra cui supernovae (SNe), fast blue optical transients (FBOTs), eventi di distruzione mareale e gamma-ray bursts, producono emissione sincrotrone non termica che traccia le interazioni d'urto con il mezzo circumstellare (CSM). Sebbene l'interferometria a lunghissima base (VLBI) possa sondare la struttura della sorgente, questi oggetti sono tipicamente non risolti a causa della loro distanza. Di conseguenza, determinare la geometria e la struttura interna dell'ejecta e del CSM dipende dall'interpretazione degli spettri radio.

I modelli standard di auto-assorbimento sincrotrone (SSA) prevedono una pendenza spettrale a bassa frequenza di $F_\nu \propto \nu^{5/2}$ e una specifica pendenza ad alta frequenza determinata dall'indice della legge di potenza degli elettroni $p$. Tuttavia, le osservazioni rivelano frequentemente deviazioni: pendenze a bassa frequenza più piatte del valore $5/2$, pendenze ad alta frequenza più piatte del previsto, o transizioni più ampie vicino alla frequenza di auto-assorbimento ($\nu_{sa}$). Le precedenti interpretazioni di queste deviazioni sono state spesso fenomenologiche, prive di un quadro rigoroso per inferire la specifica distribuzione spaziale delle inomogeneità o delle asimmetrie all'interno della regione emissiva partendo da dati non risolti.

Metodologia
Gli autori propongono un quadro generale, indipendente dal modello, per inferire la struttura della regione emissiva trattando lo spettro osservato come una sovrapposizione di emissioni provenienti da molteplici "patch" (sotto-regioni) con proprietà fisiche variabili.

1. Modellazione Forward delle Patch: L'emissione totale osservata $F_{\nu, \text{obs}}$ è modellata come l'integrale dei contributi di patch caratterizzate dalla loro frequenza di auto-assorbimento $\nu_{sa}$ e dal flusso di picco $F_{\nu_{sa}}$. La distribuzione di queste patch è descritta da una funzione $n(\nu_{sa}, F_{\nu_{sa}})$.
2. Formulazione del Problema Inverso: Gli autori definiscono una funzione di distribuzione ridotta $F(\nu_{sa})$, che rappresenta il flusso di picco totale contribuito da tutte le patch che condividono lo stesso $\nu_{sa}$. Lo spettro osservato è la convoluzione di $F(\nu_{sa})$ con la forma spettrale SSA standard.
   • Il problema è formulato come un problema inverso (Equazione 10): $F_{\nu, \text{obs}} = \int F(\nu_{sa}) f_\nu(\nu_{sa}) d\nu_{sa}$.
   • A causa della natura mal posta dell'inversione (non unicità e instabilità con il rumore), gli autori impiegano la regolarizzazione di Tikhonov per risolvere $F(\nu_{sa})$.
3. Analisi della Degenerazione: Il documento analizza rigorosamente le degenerazioni inerenti al sistema. Dimostra che, senza vincoli aggiuntivi (ad esempio, imaging ad alta risoluzione), non è possibile recuperare univocamente la distribuzione congiunta del campo magnetico ($B$) e del raggio d'urto ($R$). Tuttavia, la distribuzione $F(\nu_{sa})$ può essere recuperata se lo spettro è densamente campionato al di sotto di $\nu_{sa}$.
4. Dipendenza dalla Pendenza Spettrale: Gli autori derivano relazioni analitiche mostrando come la correlazione tra $F(\nu_{sa})$ e $\nu_{sa}$ (modellata come $F(\nu_{sa}) \propto \nu_{sa}^\alpha$) determini la pendenza spettrale osservata. Nello specifico:
   • Se $\alpha > 3/2$, la pendenza a bassa frequenza rimane $\nu^{5/2}$, ma la regione di transizione si allarga.
   • Se $-(p+1)/2 < \alpha < 3/2$, la pendenza a bassa frequenza cambia in $\nu^{\alpha+1}$.
   • Se $\alpha < -(p+1)/2$, la pendenza ad alta frequenza viene modificata.

Contributi Chiave e Risultati
Il documento applica questo quadro a due sorgenti astrofisiche specifiche per dimostrarne l'efficacia:

• SN 2016coi (Supernova a involucro rimosso): Utilizzando dati radio multi-epoca (da 8.54 a 104.25 giorni), gli autori hanno invertito lo spettro per recuperare $F(\nu_{sa})$.

  • Risultato: L'inversione ha rivelato che la regione emissiva non è un singolo guscio omogeneo. Al contrario, richiede molteplici patch con distinti valori di $\nu_{sa}$.
  • Interpretazione Fisica: I dati suggeriscono un mezzo a densità clumpy (a grumi). Ad esempio, nelle fasi iniziali, l'emissione è dominata da patch con campi magnetici differenti per fattori di $\sim 1.7$ a $>3$, che occupano aree superficiali differenti per fattori di $\sim 2.4$ a $7$. I risultati supportano la presenza di inomogeneità nel vento stellare a un raggio di $\sim 1.1 \times 10^{16}$ cm.

• AT2018cow (FBOT): A causa del campionamento rado, una soluzione inversa completa è stata ritenuta degenera. Invece, gli autori hanno utilizzato la modellazione forward per vincolare l'evoluzione della regione emissiva.

  • Risultato: L'evoluzione spettrale (appiattimento della pendenza a bassa frequenza e inasprimento della pendenza ad alta frequenza) indica un cambiamento drammatico nella struttura della regione emissiva.
  • Interpretazione Fisica: La regione emissiva transita da uno stato omogeneo ad alto campo magnetico nelle fasi precoci ($t \sim 22$ giorni) a uno stato altamente inomogeneo e a patch dominato da bassi campi magnetici nelle fasi tardive ($t \sim 83-132$ giorni). Il fattore di copertura delle regioni ad alto campo scende di ordini di grandezza ($\sim 100$ a $1000$ volte). Questa evoluzione è interpretata come prova significativa di una asimmetria su larga scala, coerente con un outflow guidato da un jet che interagisce con un CSM asimmetrico e denso.

Significatività e Rivendicazioni
Il documento sostiene di fornire un metodo robusto e indipendente dal modello per vincolare il grado di inomogeneità e asimmetria nei transienti ad alta energia non risolti.

• Generalità dell'Applicazione: Il metodo è applicabile a SNe, FBOTs, eventi di distruzione mareale e GRB, purché siano disponibili spettri densamente campionati.
• Risoluzione dell'Ambiguità: Gli autori sostengono che le deviazioni dalla pendenza standard $\nu^{5/2}$ non sono semplici anomalie fenomenologiche, ma forniscono prova diretta della distribuzione spaziale dei campi magnetici e delle aree emissive superficiali.
• Correzione dei Parametri Fisici: Il documento dimostra che assumere l'omogeneità porta a sottostimare il vero raggio d'urto (e la velocità) fino a fattori di $\sim 2$ in sorgenti altamente inomogenee. Il formalismo derivato permette di correggere tali parametri.
• Limitazioni: Gli autori notano con modestia che, sebbene il metodo recuperi la distribuzione dei parametri fisici, non può distinguere tra asimmetrie angolari (es. un jet) e inomogeneità radiali (es. turbolenza) senza vincoli aggiuntivi (come VLBI o polarizzazione). Inoltre, l'inversione è limitata dal numero di punti dati osservativi.

Il lavoro conclude che i futuri radiotelescopi (es. SKA, ngVLA) con maggiore sensibilità e una risoluzione in frequenza più fine saranno cruciali per ottenere gli spettri densamente campionati necessari per utilizzare appieno questa tecnica di inversione e mappare le strutture complesse delle regioni emissive dei transienti.