Relativistic jets from millisecond proto-magnetars

Resumen Técnico: Jets Relativistas de Proto-magnetares de Milisegundos

Planteamiento del Problema
Los proto-magnetares de milisegundos —estrellas de neutrones que rotan rápidamente ($P \sim 1$ ms) y están fuertemente magnetizadas ($B \gtrsim 10^{15}$ G), formadas en supernovas de colapso de núcleo, fusiones de estrellas de neutrones o colapso inducido por acreción de enanas blancas (AIC)— han sido propuestos durante mucho tiempo como motores centrales de los brotes de rayos gamma (GRB). Un desafío teórico crítico persiste: durante los primeros segundos tras su nacimiento, estos objetos emiten flujos intensos de neutrinos que calientan la superficie e impulsan vientos ricos en bariones. Este "viento de neutrinos" carga el flujo con materia, lo que potencialmente suprime la magnetización ($\sigma$) y limita el factor de Lorentz de volumen ($\Gamma_\infty$) alcanzable. Los modelos semianalíticos previos, calibrados para simulaciones unidimensionales y asumiendo geometrías de monopolo dividido promediadas angularmente, predijeron $\sigma \lesssim 1$ en los tiempos tempranos, lo que sugiere que los proto-magnetares no podrían lanzar los jets ultra-relativistas ($\Gamma \gtrsim 10-100$) requeridos para alimentar los GRB de corta duración. No se había abordado plenamente la cuestión de si los efectos multidimensionales podrían superar esta carga de bariones para producir jets relativistas dentro de los segundos de su formación.

Metodología
Los autores presentan las primeras simulaciones de magnetohidrodinámica de relatividad general (GRMHD) en tres dimensiones de vientos de proto-magnetares calentados por neutrinos. El estudio utiliza una versión modificada del código de código abierto GRHydro (parte de la Einstein Toolkit) empleando un esquema de volumen finito con reconstrucción WENO-Z y un resolvedor de Riemann HLLE aproximado. Características numéricas clave incluyen:

• Transporte de Neutrinos: Se utiliza una aproximación de un momento ($M_0$) de la ecuación de Boltzmann de relatividad general para modelar de forma autoconsistente el calentamiento y enfriamiento por neutrinos, incluyendo procesos $\beta$ de corriente cargada, aniquilación de pares y dispersión.
• Evolución del Campo Magnético: El campo magnético se evoluciona mediante el potencial electromagnético de 4 vectores en el gauge de Lorenz generalizado utilizando un esquema de transporte de constricción upwind para mantener la restricción solenoidal con precisión de máquina.
• Condiciones Iniciales: Los modelos asumen una estrella de neutrones de $1.4 M_\odot$ con una temperatura central de $20$ MeV. Las simulaciones abarcan un rango de periodos de rotación iniciales ($P \approx 1 - 18$ ms) y fuerzas de campo magnético superficial polar ($B_P \approx 10^{13} - 3 \times 10^{15}$ G).
• Configuración: Las simulaciones emplean simetría de reflexión a través del plano ecuatorial y utilizan refinamiento de malla adaptativa (Berger-Oliger) con la rejilla más interna resolviendo la superficie del proto-magnetar. Los modelos se evolucionan primero hacia un estado estacionario hidrodinámico impulsado por neutrinos antes de inicializar un campo dipolar poloidal.

Contribuciones Clave y Resultados
Las simulaciones revelan que los efectos multidimensionales alteran fundamentalmente la estructura del viento en comparación con los modelos 1D o promediados angularmente, particularmente para los modelos de rotación rápida ($P \approx 1$ ms) y alta magnetización ($B_P \approx 3 \times 10^{15}$ G):

1. Estratificación Angular del Flujo: El flujo no es uniforme, sino que se divide en tres regiones distintas:

   • Cinturón Ecuatorial ($\theta \gtrsim 30^\circ$): Las fuerzas centrífugas aumentan significativamente la pérdida de masa cerca del ecuador rotacional. Esto crea un viento sub-relativista denso ($v \sim 0.1c$) con una alta carga de bariones.
   • Jet Polar ($\theta \lesssim 10^\circ$): Las líneas de campo polar experimentan un aumento mínimo de la carga de masa por fuerza centrífuga. En consecuencia, esta región permanece "limpia" (baja carga de bariones) y altamente magnetizada.
   • Región de Transición: Latitudes intermedias donde el viento se acelera más eficientemente.

2. Colimación Intrínseca del Jet: El viento ecuatorial denso y sub-relativista actúa como un medio de confinamiento, colimando naturalmente el flujo polar con menor carga de bariones. Esta auto-colimación ocurre a radios pequeños ($\lesssim 100$ km) sin requerir el confinamiento por parte de eyecciones de supernova o de fusión externas.

3. Logro de Condiciones Ultra-Relativistas: En los modelos de rotación y magnetización más rápidas (por ejemplo, el modelo B3P1), el jet polar alcanza valores de magnetización pico de $\sigma \sim 30 - 100$. Esto es suficiente para alcanzar factores de Lorentz de volumen $\Gamma_\infty \sim 100$ en escalas mayores. El ángulo de apertura medio del jet es de aproximadamente $8^\circ$, consistente con las estimaciones observacionales para GRB cortos y largos.

4. Partición de Energía: La energía total del viento se particiona entre un jet polar ultra-relativista y un viento ecuatorial de ángulo ancho y sub-relativista. Este último transporta la mayor parte de la energía cinética y la masa, potencialmente alimentando la eyección de la supernova o kilonova, mientras que el primero alimenta el GRB.

Significancia y Reivindicaciones
El artículo afirma demostrar que los proto-magnetares de milisegundos pueden lanzar jets relativistas dentro de los segundos de su formación, resolviendo la tensión de larga data entre la carga de bariones impulsada por neutrinos y los requisitos de los motores de GRB. Los autores argumentan que:

• La Multidimensionalidad es Esencial: La carga de bariones no impide la formación de jets relativistas cuando se tiene en cuenta la dependencia latitudinal de las fuerzas centrífugas y el calentamiento por neutrinos.
• Escenario de Motor Unificado: Un solo motor central puede alimentar simultáneamente un GRB de corta duración (a través del jet polar) y la supernova/kilonova asociada (a través del viento ecuatorial), explicando la partición observada entre las energías de GRB corregidas por el haz y las energías cinéticas de la eyección.
• Consistencia Temporal: El mecanismo opera dentro de los primeros segundos de la vida del proto-magnetar, consistente con la fase de emisión prompt de los GRB cortos.

Los autores señalan que, aunque sus modelos asumen campos dipolares ordenados y luminosidades de neutrinos fijas, los resultados sugieren que la emergencia de jets relativistas es una consecuencia genérica de la rotación y magnetización extremas en las etapas iniciales de la evolución de la proto-estrella de neutrones. Reconocen que el trabajo futuro debe abordar los campos magnéticos generados por dinamos de forma autoconsistente, el enfriamiento/desaceleración dependiente del tiempo y la interacción con discos de acreción masivos que rodean los remanentes de fusión.