A Three-Dimensional Exploration of Magnetic Fields, Rotation, and Shock Revival in a $39 M_\odot$ Core-Collapse Supernova Progenitor

Resumen Técnico: Una Exploración Tridimensional de Campos Magnéticos, Rotación y Reviviscencia de la Onda de Choque en una Progenitora de Supernova de Colapso del Núcleo de 39 M⊙

Problema y Motivación
La reviviscencia de la onda de choque estancada en las supernovas de colapso del núcleo (CCSNe) sigue siendo un problema complejo que implica la interacción entre el calentamiento por neutrinos, la rotación, los campos magnéticos y la hidrodinámica. Si bien el mecanismo de neutrinos es robusto en muchas progenitoras, particularmente aquellas de alta compacidad, el mecanismo magnetorrotacional (MR) ofrece una vía alternativa para la reviviscencia de la onda de choque y la generación de transitorios energéticos (por ejemplo, hipernovas, supernovas Tipo Ic-BL). Este mecanismo requiere rotación rápida y campos magnéticos intensos para convertir la energía libre rotacional en energía magnética mediante el enrollamiento y las inestabilidades, impulsando flujos de salida bipolares. Sin embargo, estudios anteriores han dependido a menudo de la axisimetría (2D), lo que puede estabilizar artificialmente los chorros y suprimir inestabilidades no axisimétricas como la inestabilidad de tipo kink. Además, el impacto específico de los campos magnéticos sobre las escalas de tiempo de reviviscencia de la onda de choque, la morfología del flujo de salida y la estabilidad de la estrella de neutrones proto (PNS) frente a la formación de agujeros negros en entornos 3D de alta compacidad y rápida rotación requiere mayor investigación.

Metodología
Los autores realizaron simulaciones hidrodinámicas (HD) y magnetohidrodinámicas (MHD) tridimensionales de una progenitora de secuencia principal de edad cero (ZAMS) de 39 M⊙ con alta compacidad ($\xi_{2.5} = 0.75$) y rápida rotación ($\Omega_c \approx 0.58$ rad s$^{-1}$). La progenitora fue evolucionada utilizando MESA con evolución cuasi-químicamente homogénea. Se simularon tres modelos distintos utilizando el código FLASH con el solucionador Spark, empleando refinamiento de malla adaptativo (AMR) y un esquema de transporte de neutrinos bidimensional, dependiente de la energía y de dos momentos con cierre M1:

1. B1R1: Un modelo rotante y magnetizado. Se inicializó un campo semilla poloidal ($B_0 \approx 5 \times 10^{11}$ G), alineado con el eje de rotación.
2. B0R1: Un modelo de control rotante y no magnetizado (solo HD).
3. B0R0: Un modelo base no rotante y no magnetizado (solo HD).

Las simulaciones evolucionaron desde el pre-colapso hasta el rebote del núcleo y hacia la fase temprana post-rebote (hasta $\sim$1 segundo), utilizando una cuadrícula cartesiana con una resolución más fina de $\sim$0.61 km cerca del núcleo. Los efectos de la relatividad general se aproximaron mediante un potencial efectivo (Caso A de Marek et al. 2006).

Resultados Clave
Los tres modelos lograron una expansión descontrolada de la onda de choque dentro del intervalo simulado, pero exhibieron morfologías, escalas de tiempo y destinos de remanente marcadamente diferentes:

• Reviviscencia de la Onda de Choque y Morfología:

  • B1R1 (Magnetizado/Rotante): Exhibió la reviviscencia de la onda de choque más temprana ($t_{rev} \approx 150$ ms post-rebote) y desarrolló el flujo de salida bipolar más claro. El flujo de salida fue inicialmente similar a un chorro, pero se ensanchó y distorsionó debido a inestabilidades no axisimétricas.
  • B0R0 (No rotante): Revivió más tarde ($t_{rev} \approx 250$ ms) mediante turbulencia impulsada por neutrinos. La onda de choque se expandió en grandes penachos desordenados sin una orientación polar preferida. Crucialmente, este modelo colapsó a un agujero negro aproximadamente 1 segundo después del rebote debido a la falta de soporte rotacional.
  • B0R1 (Rotante/No magnetizado): Experimentó la reviviscencia de la onda de choque más tardía ($t_{rev} \approx 370$ ms). La rotación suprimió la convección y redujo las luminosidades de neutrinos, conduciendo a una evolución más compacta en comparación con B0R0. No desarrolló un flujo de salida bipolar claro.

• Momento Angular y Transporte de Energía:

  • En B1R1, las tensiones de Maxwell redistribuyeron eficientemente el momento angular hacia afuera desde la PNS interna, reduciendo el giro del núcleo interno y extrayendo energía libre rotacional. Este proceso aprovechó el reservorio de rotación diferencial para impulsar el flujo de salida polar.
  • El presupuesto de momento angular reveló que, si bien la advección suministró momento angular al núcleo interno, los torques magnéticos proporcionaron el sumidero dominante, compensando el aumento de giro. La emisión de neutrinos actuó como un sumidero secundario y menor.
  • La energía diagnóstica ($E_{diag}$) en B1R1 creció más temprano y más rápidamente, consistente con la extracción de energía rotacional. El presupuesto de energía magnética indicó que una fracción sustancial de la energía extraída pasó a través de un término de disipación efectivo ($\dot{Q}_{diss}$), probablemente asociado con reconexión no resuelta y transferencia subcuadrícula, en lugar de retenerse como energía magnética resuelta.

• Inestabilidades y Dinámica 3D:

  • El flujo de salida magnetizado en B1R1 estuvo sujeto a una inestabilidad de tipo kink no axisimétrica (modo $m=1$). Los diagnósticos mostraron un crecimiento exponencial temprano del desplazamiento del baricentro de la energía magnética (tiempo de e-folding $\sim$1.5 ms), consistente con las escalas de tiempo teóricas del kink.
  • A diferencia de las simulaciones axisimétricas que podrían predecir un chorro estable y estrecho, la simulación 3D mostró que el embudo del flujo de salida se inclinó (varios grados) y se ensanchó con el tiempo. El flujo de salida sobrevivió y se propagó a grandes radios, pero perdió su colimación estrecha, lo que sugiere que los efectos 3D son críticos para modelar con precisión la estabilidad de los chorros.

• Nucleosíntesis:

  • El post-procesamiento de los eyectas mostró diferencias composicionales distintas. B0R0 produjo eyectas ricos en protones ($Y_e > 0.5$) típicos de explosiones impulsadas por neutrinos, sintetizando hasta elementos del grupo del hierro ($A \sim 56$).
  • B1R1 produjo un componente significativo rico en neutrones ($Y_e$ hasta $\sim 0.33$) debido al desarrollo rápido del flujo de salida, permitiendo un proceso r débil hasta $A \sim 90$. Sin embargo, el campo semilla inicial y el ensanchamiento inducido por el kink impidieron la producción de eyectas extremadamente ricos en neutrones ($Y_e < 0.32$) observados en modelos con campos magnéticos más extremos.

Significado y Afirmaciones
El artículo afirma que en una progenitora extrema de alta compacidad, la rotación rápida combinada con campos magnéticos intensos puede abrir una vía temprana, asistida magnéticamente, hacia la reviviscencia de la onda de choque que es intrínsecamente tridimensional. Las conclusiones clave incluyen:

1. Campos Magnéticos como Impulsores: Las tensiones magnéticas no solo dan forma al flujo de salida; extraen activamente energía libre rotacional, frenan el giro del núcleo interno y facilitan la reviviscencia de la onda de choque antes que en los casos rotantes no magnetizados.
2. Necesidad de 3D: La evolución del flujo de salida impulsado magnéticamente es inherentemente no axisimétrica. La inestabilidad de tipo kink ensancha e inclina el flujo de salida, evitando la formación de un chorro estrecho persistente. Las simulaciones axisimétricas pueden sobreestimar tanto la estabilidad como la colimación de dichos flujos de salida.
3. Estabilidad del Remanente: La rotación proporciona un soporte significativo contra la formación de agujeros negros inmediatos en progenitoras de alta compacidad, aunque la estabilidad a largo plazo del remanente permanece incierta más allá del intervalo simulado.
4. Vía de Energía: La transferencia de energía desde la rotación diferencial hacia el flujo de salida implica una disipación no resuelta significativa, destacando la importancia de la física subcuadrícula en las simulaciones MHD de CCSN.

Los autores señalan limitaciones, incluido el uso de un campo semilla fuerte para compensar la amplificación no resuelta del dínamo a pequeña escala, la sensibilidad a la resolución en los diagnósticos de disipación y la necesidad de simulaciones más largas para determinar el destino final de los modelos rotantes. Enfatizan que, si bien las simulaciones apoyan un escenario de reviviscencia de la onda de choque magnetorrotacional, los resultados específicos dependen en gran medida de la dinámica 3D completa del flujo de salida.