Evolution of proto-neutron stars to pulsars, magnetars and central compact objects

Resumen Técnico: Evolución de las Proto-estrellas de Neutrones a Púlsares, Magnetares y Objetos Compactos Centrales

Planteamiento del Problema
El origen de los campos magnéticos ultrafuertes observados en los magnetares ($B \sim 10^{14}–10^{15}$ G) sigue siendo objeto de debate. Mientras que la hipótesis del "campo fósil" (conservación del flujo de las estrellas progenitoras) puede explicar los campos de las estrellas de neutrones convencionales y de los Objetos Compactos Centrales (CCOs), tiene dificultades para explicar los campos de los magnetares sin invocar condiciones progenitoras extremas y raras. Por el contrario, la hipótesis del dínamo, que requiere una rotación rápida y una convección vigorosa durante la fase de la proto-estrella de neutrones (PNS), es ampliamente aceptada para los magnetares, pero enfrenta un desafío: la tasa de nacimiento inferida de los magnetares es comparable a la tasa de supernovas de colapso de núcleo, lo que sugiere que la formación de magnetares no es tan rara como implicarían los requisitos de parámetros extremos. Además, los cálculos estándar de conservación de flujo sugieren que incluso con un núcleo de hierro colapsando, los campos resultantes ($B \sim 10^{10}$ G) son insuficientes para explicar los campos típicos de los púlsares ($B \sim 10^{12}$ G), lo que implica la necesidad de un mecanismo de dínamo incluso para los púlsares ordinarios. El problema central es reconciliar las altas tasas de nacimiento de diversos tipos de estrellas de neutrones con sus diversos grados de intensidad de campo magnético utilizando un marco evolutivo unificado.

Metodología
Los autores emplean un modelo de dínamo $\alpha-\Omega$ reducido, unidimensional y dirigido por cizalladura, adaptado de Wickramasinghe et al. (2014), quienes originalmente lo aplicaron a fusiones de enanas blancas. El modelo rastrea la evolución de los campos magnéticos poloidal ($B_p$) y toroidal ($B_\phi$) durante la fase de la PNS. Las características clave de la metodología incluyen:

• Ecuaciones Gobernantes: El modelo resuelve ecuaciones diferenciales acopladas para $B_p$ y $B_\phi$. El campo toroidal es generado mediante el efecto $\Omega$ (rotación diferencial actuando sobre el campo poloidal), mientras que el campo poloidal es regenerado mediante el efecto $\alpha$ (convección turbulenta actuando sobre el campo toroidal).
• Evolución de la PNS: A diferencia de los modelos estáticos, este marco incorpora la contracción dinámica de la PNS desde un radio inicial $R_0 \sim 40$ km hasta un radio final $R_\infty \approx 12$ km. La evolución del radio está parametrizada según Camelio et al. (2017) para diversas masas y ecuaciones de estado (EoS).
• Masa y Momento Angular: El modelo tiene en cuenta la disminución de la masa gravitacional debido al enfriamiento por neutrinos y la conservación del momento angular, lo que impulsa el aumento de la velocidad de rotación (spin-up) de la PNS.
• Viscosidad y Amortiguación de Cizalladura: El amortiguamiento de la rotación diferencial (tasa de cizalladura $q$) se modela mediante una combinación de viscosidad impulsada por la inestabilidad magnetorrotacional (MRI) y viscosidad convectiva. La tasa de cizalladura decae con el tiempo, terminando eventualmente el dínamo.
• Pérdida de Flujo y Estabilidad: El modelo incluye un término de reconexión para contabilizar la pérdida de flujo magnético (establecido inicialmente en 20%) e incorpora criterios de estabilidad (Braithwaite 2009) donde las inestabilidades de campo determinan las escalas de tiempo de decaimiento de $B_p$ y $B_\phi$.
• Escaneo del Espacio de Parámetros: Los autores integran numéricamente las ecuaciones utilizando el método Runge-Kutta, variando sistemáticamente los parámetros iniciales, incluyendo la tasa de cizalladura inicial ($q_0$), la masa inicial ($M_0$), el radio inicial ($R_0$), los campos semilla iniciales y el periodo de rotación final ($P_\infty$).

Resultos Clave
Las simulaciones revelan que la configuración final del campo magnético es altamente sensible a la tasa de rotación inicial y a la cizalladura, lo que permite que un solo mecanismo de dínamo produzca todo el espectro de tipos de estrellas de neutrones observados:

1. Umbral de Tasa de Cizalladura: El valor de saturación del campo toroidal ($B_{\phi, \infty}$) depende fuertemente de la tasa de cizalladura inicial $q_0$. Existe un umbral crítico en $q_0 \simeq 0.23$.

   • Por debajo del Umbral ($q_0 < 0.23$): El efecto $\alpha$ no se activa suficientemente. El campo poloidal permanece cerca del valor dictado por la conservación del flujo ($B_p \approx 2.5 \times 10^{10}$ G), y el campo toroidal no crece significativamente. Este régimen corresponde a los Objetos Compactos Centrales (CCOs).
   • Por encima del Umbral: El dínamo se activa. Tasas de cizalladura más altas conducen a campos toroidales más fuertes.

2. Formación de Diferentes Poblaciones:

   • Magnetares: Para los periodos iniciales más cortos (rotación rápida, por ejemplo, $P_\infty \sim 2.5$ ms), el modelo produce intensidades tipo magnetar con $B_p \simeq 10^{15}$ G y $B_\phi \simeq 10^{16}$ G.
   • Púlsares Ordinarios: Para rotadores iniciales más lentos, el modelo arroja campos típicos de púlsares ($B_p \simeq 10^{12}$ G, $B_\phi \simeq 10^{14}$ G).
   • Magnetares de Bajo-B / Púlsares de Alto-B: Las tasas de rotación intermedias producen campos en el rango de los magnetares de bajo-B o púlsares de alto campo.

3. Papel de la Conservación del Flujo: La intensidad del campo poloidal inicial está determinada en gran medida por la conservación del flujo durante el colapso. Sin embargo, el proceso de dínamo (específicamente el efecto $\alpha$) es necesario para amplificar estos campos a los niveles observados en púlsares y magnetares. El campo toroidal está genéricamente dominado por el efecto $\Omega$.

4. Objetos Compactos Centrales (CCOs): El modelo sugiere que los CCOs no experimentan un proceso de dínamo vigoroso debido a su rotación inicial lenta ($P_0 \sim 100$ ms). Sus campos poloidales están determinados por la conservación del flujo, pero aún pueden adquirir campos toroidales amplificados únicamente por el efecto $\Omega$, lo que explica sus pequeños campos dipolares pero potencialmente significativos componentes toroidales internos.

Significación y Reivindicaciones
El artículo afirma reconciliar la aparente discrepancia entre las altas tasas de nacimiento de los magnetares y las condiciones extremas que a menudo se requieren en los escenarios de formación. Al demostrar que un proceso de dínamo que opera durante la fase de la PNS puede producir naturalmente un continuo de intensidades de campo magnético dependiendo de la tasa de rotación inicial, los autores proponen que:

• Tanto los púlsares ordinarios como los magnetares probablemente pasan por una fase de dínamo; la diferencia radica en la eficiencia de este proceso impulsado por la rotación inicial.
• La hipótesis del "campo fósil" es insuficiente para explicar los campos típicos de los púlsares, lo que requiere un dínamo incluso para las estrellas de neutrones que no son magnetares.
• Los CCOs representan el extremo de baja rotación de este espectro donde el dínamo se inhibe, dejando campos dominados por la conservación del flujo.

Los autores enfatizan que su modelo unidimensional, aunque simplificado, captura las dependencias esenciales de la rotación, la cizalladura y la flotabilidad, y se alinea cualitativamente con simulaciones tridimensionales más complejas (por ejemplo, Raynaud et al. 2020). Concluyen que el dínamo $\alpha-\Omega$ proporciona un marco unificado para comprender la diversidad de los campos magnéticos de las estrellas de neutrones jóvenes sin requerir condiciones progenitoras extremas y raras para los magnetares.