Probing the Diversity of Type Ia Supernova Remnants in 3-D Hydrodynamic Simulations with X-ray Spectral Synthesis

Resumen Técnico: Sondeando la Diversidad de los Restos de Supernovas Tipo Ia en Simulaciones Hidrodinámicas 3-D con Síntesis Espectral de Rayos X

Planteamiento del Problema
Las supernovas Tipo Ia (SNe) son candelas estándar críticas para la cosmología, sin embargo, la diversidad de sus sistemas progenitores (degenerado simple vs. degenerado doble) y mecanismos de explosión (masa cercana a la de Chandrasekhar vs. masa sub-Chandrasekhar, deflagración vs. detonación) sigue siendo incierta. Si bien las observaciones de rayos X de los restos de supernovas (SNR) proporcionan diagnósticos para estos mecanismos —particularmente a través de las energías centroides y las luminosidades de las líneas Fe-K$\alpha$—, los estudios teóricos anteriores han sido limitados. La mayoría se basó en simulaciones hidrodinámicas 1-D que ignoran las asimetrías espaciales y de velocidad inherentes a las explosiones multidimensionales, o en simulaciones 3-D que se centraron en la morfología sin calcular espectros de rayos X observables. En consecuencia, ha faltado un marco autoconsistente que vincule los modelos de explosión 3-D con su emisión de rayos X evolucionada en el SNR, lo que dificulta la capacidad de interpretar la diversidad observada de los SNR Tipo Ia.

Metodología
Los autores realizaron un estudio exhaustivo que combina simulaciones hidrodinámicas 3-D con síntesis espectral de rayos X de alta resolución.

• Hidrodinámica: Se evolucionaron seis modelos de explosión Tipo Ia distintos desde la fase de supernova hasta 1000 años en un medio interestelar (ISM) uniforme con $n_{\text{ISM}} = 0.3 \text{ cm}^{-3}$. Los modelos incluyeron dos escenarios de deflagración pura, dos de detonación retardada (masa cercana a la de Chandrasekhar) y dos de doble detonación (masa sub-Chandrasekhar).
• Esquema Numérico: Las simulaciones utilizaron el código VH-1 con una malla cartesiana en movimiento para manejar el cambio dramático en las escalas físicas. La ionización fuera de equilibrio (NEI) se trató de manera autoconsistente utilizando partículas trazadoras lagrangianas (100,000 en total) para rastrear los estados locales del plasma ($n_e, T_e, T_{\text{ion}}$, abundancias) en cada paso de tiempo, evitando aproximaciones de posprocesamiento.
• Síntesis Espectral: Utilizando el software Simulated Observations of X-ray Sources (SOXS), los autores sintetizaron espectros de rayos X con una resolución de $\sim 1$ eV (superando las capacidades de XRISM/Resolve). Esto incluyó el ensanchamiento térmico, los desplazamientos Doppler debidos al movimiento global y la absorción interestelar, permitiendo tanto el análisis espectral integrado en volumen como el resuelto espacialmente.

Contribuciones Clave
Este trabajo presenta el primer cálculo de espectros de rayos X para SNR Tipo Ia derivados de simulaciones hidrodinámicas 3-D que siguen la evolución de manera autoconsistente desde la fase de SN hasta la fase de SNR. Las contribuciones técnicas clave incluyen:

1. Sondeo Sistemático de Parámetros: El esquema numérico eficiente permitió un sondeo completo 3-D de seis modelos de explosión diversos, aislando los efectos de las propiedades de la explosión en la evolución del SNR.
2. Diagnósticos Espectrales 3-D: El estudio avanza más allá de la morfología global para analizar cómo las distribuciones de eyección 3-D y las asimetrías se imprimen en los perfiles de líneas de rayos X, incluidos los efectos de la velocidad a lo largo de la línea de visión.
3. Integración de NEI: Al resolver la evolución de NEI de manera autoconsistente con la hidrodinámica mediante partículas trazadoras, el estudio proporciona una representación más precisa de los estados de ionización y las medidas de emisión que los estudios anteriores 1-D o 3-D posprocesados.

Resultados

• Evolución del Plasma: Las escalas de tiempo de ionización y las temperaturas de los electrones varían significativamente entre los modelos. Los modelos de deflagración, con una expansión más lenta, retienen densidades más altas en el interior de los ejecta, lo que conduce a una ionización más eficiente y estados de ionización más altos en comparación con los modelos de detonación retardada. Con el tiempo, la emisión pasa de estar dominada por los ejecta chocados al ISM choquado; esta transición ocurre antes en los modelos de expansión rápida (detonación retardada) que en los más lentos (deflagración).
• Distribución Elemental: Las fracciones de masa no chocadas de los elementos evolucionan de manera diferente según sus distribuciones radiales 3-D iniciales. Los elementos pesados (Fe, Ni) concentrados en el núcleo permanecen sin chocar más tiempo que los IME exteriores. La relación de masas chocadas para elementos específicos (por ejemplo, alto Ti/bajo Ni en modelos de doble detonación) sirve como una firma distintiva del mecanismo de explosión.
• Diagnósticos Fe-K$\alpha$: Las energías centroides y las luminosidades de la línea Fe-K$\alpha$ sintetizadas reproducen gran parte de la diversidad observada de los SNR Tipo Ia. Los modelos se agrupan en tres regiones basadas en el mecanismo:
  • Detonación Retardada (DDT): Menor ionización, mayor luminosidad.
  • Deflagración (DEF): Mayor ionización, mayor luminosidad.
  • Doble Detonación (DD): Menor ionización, menor luminosidad.
    Sin embargo, el estudio señala que los restos observados más altamente ionizados y luminosos (por ejemplo, 0519–69.0, N103B, 3C 397) no pueden reproducirse bajo la suposición de un ISM uniforme, lo que sugiere que estos objetos interactúan con medios circumestelares (CSM) densos e inhomogéneos.
• Perfiles de Línea y Asimetría: Los efectos 3-D introducen asimetrías significativas en los perfiles de línea. Si bien el ensanchamiento térmico domina en épocas tempranas ($\sim 400$ años), los desplazamientos Doppler se vuelven cada vez más visibles, particularmente para Fe-K$\alpha$ debido a su mayor energía. Los modelos asimétricos exhiben perfiles de línea sesgados o desplazados dependiendo del ángulo de visión. En los modelos de doble detonación, la presencia de una enana blanca compañera crea una región de "sombra" (un cono de baja densidad), lo que lleva a firmas espectrales distintas en regiones específicas fuera del centro que difieren entre los escenarios de explosión secundaria simple y doble.

Significado y Afirmaciones
Los autores afirman que su marco de modelado 3-D reproduce con éxito la diversidad observada de los SNR Tipo Ia y proporciona restricciones cualitativas sobre los sistemas progenitores y los mecanismos de explosión.

• Poder Diagnóstico: El estudio demuestra que las propiedades de la línea Fe-K$\alpha$ (centroide y luminosidad) y los perfiles de línea resueltos espacialmente pueden distinguir entre mecanismos de explosión (por ejemplo, cerca de MCh vs. sub-MCh, deflagración vs. detonación).
• Limitaciones del 1-D: La comparación con modelos promediados 1-D destaca que ignorar las distribuciones espaciales y de velocidad 3-D conduce a errores sistemáticos en las temperaturas de electrones inferidas y en las formas de línea, subestimando la contribución del material choquado en ejecta asimétricos.
• Perspectiva Futura: Los autores posicionan este trabajo como un paso necesario hacia la interpretación de datos de alta resolución de XRISM. Declaran explícitamente que, si bien sus modelos de ISM uniforme explican una amplia gama de observaciones, la diversidad completa (específicamente los restos de alta luminosidad/alta ionización) requiere extensiones futuras para incluir entornos 3-D de CSM. El artículo concluye que confrontar estos espectros sintéticos con observaciones de restos específicos como Tycho (SN 1572) ofrece una ruta prometedora para discriminar entre modelos de progenitores.