Born to be recycled: a comprehensive population synthesis of the Galactic millisecond pulsars

1. Planteamiento del Problema

Los púlsares de milisegundos (MSPs) son los objetos más antiguos pero también los más rápidos conocidos, con periodos de rotación menores a 30 ms. La hipótesis predominante sugiere que se forman mediante un mecanismo de "reciclaje": un púlsar antiguo y lento en un sistema binario acreta materia de su compañera, lo que aumenta su momento angular y lo acelera hasta frecuencias de rotación extremas.

Sin embargo, existen incertidumbres significativas sobre los detalles de este proceso evolutivo, incluyendo:

• La tasa de nacimiento real de los MSPs.
• La evolución del campo magnético durante la acreción y el tiempo libre.
• La distribución de masas finales de los púlsares reciclados.
• La contribución de los MSPs al "exceso de GeV" observado en el centro galáctico (¿son de origen en cúmulos globulares migrados o nacidos en el bulbo galáctico?).
• La existencia de canales de formación alternativos, como el colapso inducido por acreción (AIC) de enanas blancas de oxígeno-neón (ONeWD).

El objetivo del trabajo es desarrollar una síntesis de población completa para verificar si el mecanismo de reciclaje por acreción puede explicar la mayoría de los púlsares observados y derivar estadísticas sobre sus propiedades que son difíciles de obtener solo mediante observaciones.

2. Metodología

Los autores desarrollaron un código de síntesis de población que combina la evolución estelar binaria con la evolución física del púlsar y modelos de detección observacional.

• Evolución Binaria (SEVN): Utilizaron el código Stellar EVolution for N-body (SEVN) para simular la evolución de sistemas binarios desde estrellas de la secuencia principal hasta la formación de objetos compactos. El código maneja procesos clave como el desbordamiento del lóbulo de Roche (RLO), transferencia de masa (estable e inestable) y la evolución orbital.
• Evolución del Púlsar (Autoconsistente): A diferencia de estudios anteriores que asumen campos magnéticos constantes o modelos simplificados, este trabajo evoluciona cada púlsar de manera autoconsistente considerando:
  • Frenado y aceleración: Torques de acreción, frenado magnético (interacción disco-campo), pérdida de radiación y frenado gravitacional (ondas gravitatorias) si la elipticidad es significativa.
  • Decaimiento del campo magnético: Se adopta una ley de potencia para el decaimiento del campo magnético ($B(t)$) con tres escalas de tiempo características ($\tau_d$) seleccionadas aleatoriamente (0.05, 1.4 y 2.3 Myr) para reproducir la distribución observada. También se incluye el decaimiento inducido por acreción (enterramiento del campo).
  • Geometría de emisión: Se modelan las emisiones de radio (polo cap) y rayos gamma (viento rayado o striped wind) para determinar la visibilidad.
• Condiciones Iniciales: Las estrellas de neutrones nacen en binarias con compañeras de secuencia principal. Se utilizan distribuciones específicas para la masa inicial, la velocidad de kick (ajustada a 70 km/s para mantener la binaria), la excentricidad y la orientación de los ejes.
• Modelo de Detección: Se aplican criterios de detección realistas basados en dos encuestas de radio (FAST GPPS y PMPS) y el catálogo de rayos gamma del telescopio Fermi-LAT (3PC y 4FGL). Se calcula el flujo de radio y gamma considerando la distancia, la geometría del haz, el fondo de cielo y la sensibilidad instrumental.

3. Contribuciones Clave

• Simulación de Evolución Completa: Es el primer estudio que sigue la evolución completa de la binaria antes y durante la fase de reciclaje, permitiendo la formación de púlsares "ligeramente reciclados" (30 ms < P < 1000 ms) y totalmente reciclados.
• Modelado Autoconsistente: Integra la evolución del campo magnético, la acreción y la geometría de emisión en un solo marco, superando las limitaciones de estudios previos que trataban estas fases por separado.
• Análisis de Canales Alternativos: Incluye una simulación específica para evaluar la contribución del colapso inducido por acreción (AIC) de enanas blancas ONe a la población de púlsares reciclados.
• Predicciones para Futuros Encuestas: Proporciona estimaciones cuantitativas sobre el número de púlsares detectables con el Square Kilometre Array (SKA) y telescopios de rayos gamma de próxima generación.

4. Resultados Principales

• Reproducción de la Población Observada: Las simulaciones reproducen bien la distribución en el diagrama $P-\dot{P}$ de los púlsares de radio y rayos gamma observados, validado mediante pruebas de Kolmogorov-Smirnov (KS).
• Púlsares No Identificados: Se estima que debe haber menos de 220 púlsares no identificados en el catálogo 4FGL de Fermi-LAT.
• Duración de la Acreción: La mayoría de los púlsares acreta materia durante duraciones típicas de $10^8 - 10^9$ años.
• Alineación de Ejes: El tiempo de alineación entre el eje de rotación y el normal al plano orbital es mucho más corto que la duración de la acreción. Como resultado, aproximadamente el 80% de la población de púlsares binarios tiene un ángulo de alineación ($\alpha$) cercano a 0°.
• Distribución de Masas: La mayoría de los púlsares reciclados tienen masas alrededor de 1.8 $M_\odot$, con una cola hacia masas más altas de hasta 2.7 $M_\odot$. Esto concuerda con las observaciones recientes de púlsares "araña" (spiders) masivos.
• Origen del Exceso de GeV: Se estima que solo ~190 MSPs nacidos en los brazos espirales galácticos contribuyen a aproximadamente el 5% del exceso de GeV en el centro galáctico. Esto refuerza la hipótesis de que si el exceso es totalmente de origen en púlsares, la gran mayoría debe provenir de cúmulos globulares que migraron hacia el centro o de púlsares nacidos en el bulbo galáctico, los cuales son actualmente indetectables.
• Geometría de Detección: Es difícil observar púlsares reciclados con un eje de rotación y magnético alineados ($\chi \le 10^\circ$). Se requieren ángulos de visión ($\zeta$) altos para detectar la mayoría de los púlsares reciclados.
• Colapso Inducido por Acreción (AIC): De las enanas blancas ONe en sistemas binarios, solo una fracción muy pequeña (~0.0076%) parece contribuir a la población de púlsares ligeramente reciclados a través del AIC. Estos púlsares tienden a tener periodos más largos y no alcanzan la categoría de MSPs completos.
• Población de "Arañas" (Spiders): El modelo reproduce la existencia de sistemas tipo "black widow" y "redback", pero subestima el número total observado, sugiriendo que faltan canales de formación adicionales o ajustes en la eficiencia de acreción/evaporación.

5. Significado e Implicaciones

Este trabajo valida robustamente el escenario de reciclaje por acreción como el mecanismo principal para la formación de la mayoría de los púlsares de milisegundos observados. Al integrar la evolución binaria completa con la física del púlsar, el estudio resuelve discrepancias anteriores sobre las tasas de nacimiento y las distribuciones de masa.

La conclusión de que los púlsares del campo galáctico contribuyen mínimamente al exceso de GeV del centro galáctico tiene implicaciones profundas para la astrofísica de altas energías, apoyando fuertemente la idea de que el exceso es causado por una población antigua y densa de púlsares en el bulbo o en cúmulos globulares migrados, en lugar de ser un signo de materia oscura.

Finalmente, las predicciones de que el SKA aumentará las detecciones de radio en un factor de ~3 y que instrumentos de rayos gamma más sensibles podrían aumentar las detecciones en un factor de ~5, proporcionan una hoja de ruta clara para las futuras campañas observacionales y la búsqueda de ondas gravitacionales continuas provenientes de esta población.