Multi-Dimensional MHD simulations of young Core-Collapse Supernova Remnants

Resumen Técnico: Simulaciones MHD Multidimensionales de Restos de Supernovas de Colapso del Núcleo Jóvenes

Problema y Motivación
Los Restos de Supernovas (RS) de colapso del núcleo jóvenes son candidatos principales para la aceleración de rayos cósmicos (RC) galácticos, particularmente en lo que respecta a la característica "rodilla" en el espectro de RC a energías de peta-electrónvoltios (PeV). Una incertidumbre crítica en el modelado de estos sistemas radica en el medio interestelar circumestelar (CSM) moldeado por la pérdida de masa de la estrella progenitora antes de la explosión. Estudios anteriores a menudo se basaron en modelos unidimensionales (1D) o simplificados multidimensionales que asumían vientos estacionarios y estáticos, así como tasas promedio de pérdida de masa ($\dot{M}$) y velocidades terminales ($v_\infty$). Estas simplificaciones no logran capturar estructuras complejas del CSM que surgen de la evolución estelar dependiente del tiempo, como capas densas, burbujas rarefactadas y múltiples choques formados por interacciones entre diferentes fases evolutivas (por ejemplo, fases de Gigante Roja [RSG] frente a fases de Wolf-Rayet [WR]). Además, la configuración del campo magnético, crucial para la aceleración de partículas, depende de la rotación estelar y la dinámica del viento, que a menudo están desacopladas de la evolución hidrodinámica en los modelos existentes.

Metodología
Los autores emplean el código MHD pion para realizar simulaciones cilíndricas 2D y cartesianas 3D de RS jóvenes expandiéndose en CSMs generados de manera autoconsistente a partir de trayectorias detalladas de evolución estelar.

• Evolución Estelar: El CSM se deriva de una trayectoria de evolución estelar 1D calculada con MESA para una estrella primaria de $31.6 M_\odot$ en un sistema binario amplio (evolucionando como una estrella única). La trayectoria cubre las fases de Secuencia Principal (MS), RSG y WR. Las prescripciones de pérdida de masa varían según la fase (Vink et al. para MS; Nieuwenhuijzen & de Jager para RSG; Nugis & Lamers y Yoon para WR).
• Configuración MHD: Las simulaciones incluyen MHD ideal, refinamiento de malla estático y calentamiento/enfriamiento radiativo. Una suposición clave es la fotoionización completa del viento y el medio interestelar (ISM) (calentamiento a $\sim 8300$ K), lo que reduce los números de Mach del choque y altera la formación de capas. El campo magnético estelar se inyecta como un monopolo dividido, arrollado en una espiral por la rotación. Los autores asumen un campo magnético superficial modesto ($B_\star = 10$ G) y una tasa de rotación inicial baja ($20\%$ de la crítica), consistente con las observaciones de estrellas O galácticas.
• Escenarios: Se modelan dos escenarios canónicos de progenitores:
  1. RSG: Un viento lento y denso expandiéndose en el ISM.
  2. WR: Un viento rápido con alta pérdida de masa que barre el material RSG precedente, creando un CSM complejo y estructurado.
• Implementación de la Supernova: Una explosión esférica ($5 M_\odot$, $10^{51}$ erg) se inserta en el CSM pre-SN con perturbaciones gaussianas ($\pm 40\%$) para romper simetrías artificiales. El campo magnético del eyecta se establece para igualar el $\beta$ del plasma del CSM circundante.

Resultados Clave

1. Estructura del CSM y Fotoionización: La suposición de fotoionización altera significativamente la estructura del CSM en comparación con modelos anteriores no ionizados.
   • En el caso RSG, la burbuja de viento no forma una capa delgada y densa en el choque de terminación del viento (WTS). En su lugar, el viento fotoionizado y el ISM resultan en una región extendida de densidad creciente a lo largo de $\sim 3$ pc, con expansión subsónica hacia el ISM.
   • En el caso WR, el viento rápido barre el material RSG hacia una capa densa, radiativa e inestable por Rayleigh-Taylor. El viento WR crea un conjunto coherente de choques reflejados rápidos cuando el RS interactúa con esta capa.
2. Dinámica de Choques:
   • Choques Directos: En ambos escenarios, las velocidades de los choques directos son más rápidas de lo predicho por teorías analíticas (por ejemplo, Truelove & McKee 1999) que asumen CSMs estáticos. Para RSGs, esto se debe a la aceleración adicional del viento por presión térmica (fotoionización). Para WRs, se debe a la velocidad de expansión del propio CSM.
   • Choques Reflejados: En la simulación WR, la interacción con la capa RSG barrida genera choques reflejados que se propagan hacia adentro, interactuando con el choque inverso y procesando energía significativa ($\sim 3 \times 10^{50}$ erg).
3. Campos Magnéticos:
   • Las intensidades del campo magnético alcanzan valores máximos de decenas de $\mu$G en regiones comprimidas y de enfriamiento (capas).
   • El efecto Axford-Cranfill (acumulación de campo azimutal $B_\phi \propto r$) no se reproduce significativamente; el campo permanece dominado radialmente debido a las bajas tasas de rotación estelar ($v_{rot}/v_{crit} \sim 1\%$).
   • La estructura del campo magnético es altamente asimétrica, particularmente en el caso RSG debido a efectos de movimiento propio.
4. Potencial de Aceleración de Partículas:
   • Choques Directos: Las energías máximas de partículas se estiman en $\sim 100$ TeV para RSGs y $\sim 30$ TeV para WRs en la primera década post-explosión. Estos valores son significativamente más bajos que las energías PeV requeridas para explicar la rodilla de los rayos cósmicos.
   • Choques Reflejados: Aunque el choque reflejado en el caso WR procesa energía sustancial, el campo magnético interno es débil ($B \ll \mu$G). En consecuencia, la energía máxima de electrones se limita a decenas de TeV. Sin embargo, los autores predicen emisión de rayos gamma detectable por efecto Compton inverso de estos electrones si están confinados dentro del RS.

Significancia y Conclusiones
El artículo demuestra que un tratamiento MHD detallado y multidimensional acoplado con evolución estelar autoconsistente es esencial para modelar con precisión la evolución temprana de los RS.

• Más allá de las Aproximaciones Analíticas: Los modelos analíticos que asumen CSMs estacionarios y estáticos subestiman las velocidades de los choques directos y no logran capturar las estructuras de densidad complejas (por ejemplo, gradientes de densidad extendidos frente a capas delgadas) que surgen de la pérdida de masa dependiente del tiempo y la fotoionización.
• Rotación Estelar y Magnetismo: La falta de rotación rápida en los modelos de progenitores conduce a campos magnéticos dominados radialmente que no se acumulan significativamente aguas abajo del WTS, limitando el potencial de aceleración de partículas a PeV en RS aislados.
• PeVatrones: Los resultados sugieren que los RS jóvenes aislados de progenitores estándar RSG o WR son poco propensos a acelerar protones a energías PeV. Lograr tales energías probablemente requiere condiciones extremas, como supernovas ocurriendo en cúmulos estelares densos, dentro de superburbujas, o interactuando con material denso expulsado inmediatamente antes de la explosión.
• Firmas Observacionales: El estudio destaca que los choques reflejados en progenitores WR, aunque no son aceleradores de PeV, pueden producir emisión de rayos gamma detectable mediante dispersión Compton inversa, ofreciendo una firma observacional potencial para estos canales evolutivos específicos.

Los autores concluyen que los modelos futuros deben tener en cuenta las inestabilidades dinámicas y la naturaleza dependiente del tiempo del CSM para predecir correctamente las propiedades de los choques y la emisión no térmica de los RS jóvenes.