Chandra X-Ray Imaging and Spatially Resolved Spectroscopy of SN 1987A: Energy-Dependent Morphology of the Equatorial Ring

Resumen Técnico: Imágenes de Rayos X de Chandra y Espectroscopía con Resolución Espacial de SN 1987A

Planteamiento del Problema
La SN 1987A, ubicada en la Gran Nube de Magallanes, sirve como un punto de referencia crítico para estudiar la interacción entre los eyectos de la supernova y el material circunestelar (CSM). Si bien la evolución del remanente ha sido monitoreada extensamente, comprender sistemáticamente sus cambios morfológicos y espectrales a largo plazo a través de diferentes bandas de energía sigue siendo un desafío. Específicamente, estudios previos han observado que la morfología de rayos X varía según la energía, donde la emisión de banda suave se asocia con protuberancias del anillo ecuatorial (ER) denso y la emisión de banda dura traza choques más rápidos. Sin embargo, observaciones recientes sugieren una transición donde la emisión de banda suave se extiende hacia el interior mientras que la de banda dura se extiende hacia el exterior. Además, la línea de emisión de Fe K (~6.7 keV), una sonda del plasma más caliente, ha sido detectada en datos de XMM-Newton desde 2010, pero su evolución espacial y su aparición temprana han sido difíciles de constreñir debido a la limitada estadística de fotones en los datos de banda dura de Chandra. Se requiere un análisis unificado de la estructura radial del remanente y su evolución espectral durante todo el periodo de Chandra (1999–2025) para caracterizar esta transición de una fase dominada por el ER a una que involucra una emisión interior significativa.

Metodología
Los autores presentan un estudio sistemático utilizando observaciones de archivo de Chandra desde 1999 hasta 2025, que comprenden 21 observaciones ACIS-S "NONE" y 79 observaciones de la Red de Transmisión de Alta Energía (HETG). Para abordar los desafíos de fusionar datos multi-época con estructuras de sub-segundo de arco y campos de visión limitados (que imposibilitan la astrometría estándar de fuente de referencia), los autores desarrollaron una estrategia de alineación interna especializada:

1. Reducción de Datos y Agrupamiento: Los datos fueron reducidos utilizando CIAO v4.17 y agrupados en intervalos de cuatro años para equilibrar la resolución temporal con la estadística de fotones.
2. Determinación del Centro y Alineación: Para cada observación, el centro del remanente se determinó aplicando el reposicionamiento de eventos subpíxel (EDSER), simulando una función de dispersión de punto (PSF) mediante MARX, aplicando la deconvolución de Richardson-Lucy (RL) y ajustando un modelo de toro elíptico. Este centro sirvió como referencia para alinear todas las épocas.
3. Análisis de Imágenes: Se co-añadieron imágenes reconstruidas por subpíxel y mapas de exposición. Las PSF fusionadas se generaron mediante el promedio ponderado por exposición. Se aplicó la deconvolución RL (30 iteraciones) a las imágenes fusionadas en tres bandas: total (0.5–7.0 keV), suave (0.5–1.5 keV) y dura (1.5–7.0 keV).
4. Modelado Morfológico: La estructura radial se cuantificó ajustando un modelo de toro proyectado en dos dimensiones a las imágenes deconvolucionadas, derivando el radio del eje mayor ($r_{ell}$) y el ancho ($\sigma_r$) en el plano de la imagen.
5. Análisis Espectral: Se extrajeron espectros con resolución espacial de la región completa y de cuatro sectores (Este, Oeste, Norte, Sur). Los ajustes se realizaron en la banda de 5.0–8.0 keV utilizando un modelo de ley de potencia absorbida más Gaussiana ($tbabs*(powerlaw+gauss)$) para caracterizar la línea de Fe K, con el ancho de la Gaussiana fijado según los resultados de XRISM.
6. Contexto Multi-longitud de Onda: Se analizaron las observaciones de 2022 de ATCA (radio), JWST (infrarrojo cercano) y HST (óptico) y se compararon con los datos contemporáneos de Chandra para contextualizar la morfología de rayos X.

Resultados Clave

• Morfología Dependiente de la Energía: El estudio confirma que, si bien las bandas suave y dura exhibieron estructuras similares en épocas tempranas, surgió una divergencia después de principios de la década de 2010. La emisión de banda suave (0.5–1.5 keV) se volvió sistemáticamente más ancha que la emisión de banda dura (1.5–7.0 keV).
• Extensión hacia el Interior: El análisis del radio y el ancho combinados revela que la emisión de banda suave ha experimentado una extensión hacia el interior. Para la década de 2020, el límite interno de la distribución de la banda suave alcanza un radio proyectado de $\sim0.5''$, comparable al radio observado alrededor del año 2000. En contraste, el límite externo permanece similar entre las bandas, lo que indica que la dependencia de la energía surge principalmente de la región interior.
• Evolución de la Emisión de Fe K: La evolución del flujo de la línea de Fe K es consistente con los resultados previos de XMM-Newton. Sin embargo, la combinación de espectros de Chandra permite la detección de la línea de Fe K en épocas anteriores (específicamente 2007–2009, y potencialmente 2004–2007 en el sector este), lo cual era difícil de constreñir anteriormente.
• Asimetría Espacial: El análisis con resolución espacial indica que la emisión de Fe K está aumentada en la región este. Aunque esta asimetría se notó claramente después de ~2018, los autores encuentran evidencia que sugiere que esta tendencia podría haber estado presente en épocas anteriores, aunque con menor significancia estadística.
• Comparación Multi-longitud de Onda: En 2022, la emisión de radio parece más simétrica azimutalmente en comparación con épocas anteriores, mientras que las imágenes ópticas y de rayos X muestran una emisión más brillante en el lado oeste. El ordenamiento radial sugiere que el componente de rayos X suaves se origina en regiones interiores al anillo óptico.

Significancia y Reivindicaciones
El artículo afirma proporcionar una visión unificada de la evolución morfológica y espectral a largo plazo de la SN 1987A utilizando un conjunto de datos de Chandra de varias décadas reducido de manera uniforme. La principal significancia radica en la identificación de la estructura radial dependiente de la energía como una característica clave de la evolución tardía del remanente.

Los autores interpretan la extensión hacia el interior de la emisión de rayos X suaves y el ensanchamiento de la banda suave como evidencia de una contribución creciente de eyectos calentados por el choque inverso y/o componentes de emisión de alta latitud. Esta observación se alinea con simulaciones hidrodinámicas y magnetohidrodinámicas tridimensionales recientes que predicen que, aproximadamente 30–35 años después de la explosión, la emisión de rayos X debería reflejar cada vez más los eyectos con choque en lugar de solo la interacción con el anillo ecuatorial. El estudio destaca así la transición de la SN 1987A de una fase dominada por el anillo ecuatorial a una donde los componentes interiores juegan un papel progresivamente importante.