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

Resumen Técnico: Restricción de Inhomogeneidades y Asimetrías en Transitorios de Alta Energía a partir de Observaciones de Radio No Resueltas

Planteamiento del Problema
Los transitorios de alta energía, que incluyen supernovas (SNe), transitorios ópticos azules rápidos (FBOTs), eventos de disrupción de marea y brotes de rayos gamma, producen emisión de sincrotrón no térmico que rastrea las interacciones de choque con el medio circunestelar (CSM). Si bien la interferometría de muy larga base (VLBI) puede sondear la estructura de la fuente, estos objetos suelen estar no resueltos debido a sus distancias. En consecuencia, determinar la geometría y la estructura interna de los eyectos y del CSM depende de la interpretación de los espectros de radio.

Los modelos estándar de autoabsorción de sincrotrón (SSA) predicen una pendiente espectral de baja frecuencia de $F_\nu \propto \nu^{5/2}$ y una pendiente de alta frecuencia específica determinada por el índice de la ley de potencia de los electrones $p$. Sin embargo, las observaciones revelan frecuentemente desviaciones: pendientes de baja frecuencia más tenues que $5/2$, pendientes de alta frecuencia más planas de lo esperado, o transiciones más anchas cerca de la frecuencia de autoabsorción ($\nu_{sa}$). Las interpretaciones previas de estas desviaciones han sido a menudo fenomenológicas, careciendo de un marco riguroso para inferir la distribución espacial específica de las inhomogeneidades o asimetrías dentro de la región emisora a partir de datos no resueltos.

Metodología
Los autores proponen un marco general e independiente del modelo para inferir la estructura de la región emisora tratando el espectro observado como una superposición de emisiones de múltiples "parches" (subregiones) con propiedades físicas variables.

1. Modelado hacia adelante de parches: El flujo total observado $F_{\nu, \text{obs}}$ se modela como la integral de las contribuciones de parches caracterizados por su frecuencia de autoabsorción $\nu_{sa}$ y su flujo pico $F_{\nu_{sa}}$. La distribución de estos parches se describe mediante una función $n(\nu_{sa}, F_{\nu_{sa}})$.
2. Formulación del problema inverso: Los autores definen una función de distribución reducida $F(\nu_{sa})$, que representa el flujo pico total contribuido por todos los parches que comparten un $\nu_{sa}$ específico. El espectro observado es la convolución de $F(\nu_{sa})$ con la forma espectral estándar de SSA.
   • El problema se formula como un problema inverso (Ecuación 10): $F_{\nu, \text{obs}} = \int F(\nu_{sa}) f_\nu(\nu_{sa}) d\nu_{sa}$.
   • Debido a la naturaleza mal planteada de esta inversión (no unicidad e inestabilidad ante el ruido), los autores emplean la regularización de Tikhonov para resolver $F(\nu_{sa})$.
3. Análisis de degeneración: El artículo analiza rigurosamente las degeneraciones inherentes al sistema. Demuestra que, sin restricciones adicionales (por ejemplo, imágenes de alta resolución), no es posible recuperar de forma única la distribución conjunta del campo magnético ($B$) y el radio del choque ($R$). Sin embargo, la distribución $F(\nu_{sa})$ puede recuperarse si el espectro está densamente muestreado por debajo de $\nu_{sa}$.
4. Dependencia de la pendiente espectral: Los autores derivan relaciones analíticas que muestran cómo la correlación entre $F(\nu_{sa})$ y $\nu_{sa}$ (modelada como $F(\nu_{sa}) \propto \nu_{sa}^\alpha$) dicta la pendiente espectral observada. Específicamente:
   • Si $\alpha > 3/2$, la pendiente de baja frecuencia permanece en $\nu^{5/2}$, pero la región de transición se ensancha.
   • Si $-(p+1)/2 < \alpha < 3/2$, la pendiente de baja frecuencia cambia a $\nu^{\alpha+1}$.
   • Si $\alpha < -(p+1)/2$, la pendiente de alta frecuencia se ve modificada.

Contribuciones Clave y Resultados
El artículo aplica este marco a dos fuentes astrofísicas específicas para demostrar su eficacia:

• SN 2016coi (Supernova de envolvente desprendida): Utilizando datos de radio multiepoque (8.54 a 104.25 días), los autores invirtieron el espectro para recuperar $F(\nu_{sa})$.

  • Resultado: La inversión reveló que la región emisora no es una única cáscara homogénea. En su lugar, requiere múltiples parches con distintos valores de $\nu_{sa}$.
  • Interpretación Física: Los datos sugieren un medio de densidad irregular (clumpy). Por ejemplo, en épocas tempranas, la emisión está dominada por parches con campos magnéticos que difieren por factores de $\sim 1.7$ a $>3$, ocupando áreas superficiales que difieren por factores de $\sim 2.4$ a $7$. Los resultados apoyan la presencia de inhomogeneidades en el viento estelar a un radio de $\sim 1.1 \times 10^{16}$ cm.

• AT2018cow (FBOT): Debido al muestreo disperso, una solución de inversión completa se consideró degenerada. En su lugar, los autores utilizaron el modelado hacia adelante para restringir la evolución de la región emisora.

  • Resultado: La evolución espectral (el aplanamiento de la pendiente de baja frecuencia y el aumento de la pendiente de alta frecuencia) indica un cambio dramático en la estructura de la región emisora.
  • Interpretación Física: La región emisora transiciona de un estado homogéneo de alto campo magnético en tiempos tempranos ($t \sim 22$ días) a un estado altamente inhomogéneo y fragmentado (patchy) dominado por bajos campos magnéticos en tiempos tardíos ($t \sim 83-132$ días). El factor de cobertura de las regiones de alto campo cae por órdenes de magnitud ($\sim 100$ a $1000$ veces). Esta evolución se interpreta como evidencia sólida de una asimetría a gran escala, consistente con un flujo impulsado por un chorro (jet) interactuando con un CSM asimétrico y denso.

Significancia y Reivindicaciones
El artículo afirma proporcionar un método robusto e independiente del modelo para restringir el grado de inhomogeneidad y asimetría en transitorios de alta energía no resueltos.

• Aplicabilidad General: El método es aplicable a SNe, FBOTs, eventos de disrupción de marea y GRBs, siempre que se disponga de espectros bien muestreados.
• Resolución de la Ambigüedad: Los autores argumentan que las desviaciones de la pendiente estándar $\nu^{5/2}$ no son meramente anomalías fenomenológicas, sino que proporcionan evidencia directa de la distribución espacial de los campos magnéticos y las áreas emisoras.
• Corrección de Parámetros Físicos: El artículo demuestra que asumir homogeneidad conduce a una subestimación del radio real del choque (y la velocidad) por factores de hasta $\sim 2$ en fuentes altamente inhomogéneas. El formalismo derivado permite la corrección de estos parámetros.
• Limitaciones: Los autores señalan modestamente que, si bien el método recupera la distribución de los parámetros físicos, no puede distinguir entre asimetrías angulares (por ejemplo, un jet) e inhomogeneidades radiales (por ejemplo, turbulencia) sin restricciones adicionales (como VLBI o polarización). Además, la inversión está limitada por el número de puntos de datos observacionales.

El trabajo concluye que los futuros radiotelescopios (por ejemplo, SKA, ngVLA) con mayor sensibilidad y resolución de frecuencia más fina serán cruciales para obtener los espectros densamente muestreados necesarios para utilizar plenamente esta técnica de inversión y mapear las estructuras complejas de las regiones emisoras de los transitorios.