Kick bimodality of neutron stars and mode dependence of their parameters

Planteamiento del Problema
Se observa que las estrellas de neutrones (EN) poseen velocidades espaciales que exceden significativamente las de sus estrellas progenitoras, un fenómeno atribuido a los patadas natales (kicks) recibidas durante las explosiones de supernovas. Si bien se ha respaldado una distribución bimodal de estas velocidades de patada (que comprende modos de baja y alta velocidad) mediante el modelado teórico y la síntesis de poblaciones (por ejemplo, Igoshev 2020), la confirmación observacional sigue siendo un desafío. Derivar la distribución de la patada natal a partir de las velocidades transversales 2D observadas no es trivial debido al desconocimiento de los componentes de velocidad radial, las incertidumbres en los lugares de nacimiento y la evolución cinemática dentro del potencial Galáctico. Además, estudios existentes no han investigado exhaustivamente si parámetros físicos específicos de los púlsares de radio se correlacionan con estos distintos modos de velocidad de patada, particularmente distinguiendo entre efectos causados por el sesgo de selección y diferencias físicas intrínsecas.

Metodología
Los autores analizaron una muestra de 202 púlsares de radio normales e aislados seleccionados del Catálogo de Púlsares ATNF (v2.7.0). Los criterios de selección requirieron mediciones robustas de paralaje o distancias de medida de dispersión < 10 kpc, movimientos propios disponibles, parámetros de rotación conocidos ($P, \dot{P}$), y la exclusión de púlsares reciclados, binarios y fuentes exóticas.

Para clasificar estos púlsares en modos de baja o alta velocidad, los autores emplearon un modelo de integración de trayectoria:

1. Integración de Trayectoria: Para cada púlsar, se muestrearon 100 valores de distancia (basados en el paralaje o la medida de dispersión). Utilizando coordenadas galácticas y movimientos propios, se generaron 100 posiciones actuales.
2. Reconstrucción de Velocidad: Dado que las velocidades radiales ($v_r$) no son observadas, los autores asumieron una orientación de velocidad isotrópica respecto al Estándar Local de Reposo (LSR). Se muestreó el ángulo $\theta$ entre la línea de visión y el vector de velocidad 3D a partir de una distribución uniforme truncada (Uniform(0.05$\pi$, 0.95$\pi$)) para evitar valores atípicos extremos, generando $10^3$ vectores de velocidad por distancia. Esto resultó en $10^5$ combinaciones de posición-velocidad por púlsar.
3. Integración hacia Atrás: Utilizando el potencial Galáctico (McMillan 2017), las trayectorias se integraron hacia atrás en el tiempo hasta $1.5 \times \tau_c$ (edad característica). Los lugares de nacimiento se identificaron ya sea como intersecciones con el plano Galáctico o en $t = -\tau_c$ si no ocurría tal intersección.
4. Cálculo de la Velocidad de Patada: En cada posible lugar de nacimiento, se calculó la velocidad de la patada ($v_{kick} = |\vec{v} - \vec{v}_{LSR}|$).
5. Categorización: Utilizando los parámetros de la distribución bimodal de patada de Igoshev (2020) (una suma ponderada de dos Maxwellianas: $\sigma_1 = 56$ km s$^{-1}$, $\sigma_2 = 336$ km s$^{-1}$, con pesos $w_1=0.2, w_2=0.8$), los autores calcularon la probabilidad de que cada púlsar perteneciera al modo de baja o alta velocidad. Los púlsares se categorizaron como pertenecientes "inequívocamente" a un modo si la probabilidad superaba 0.5 en todos los posibles lugares de nacimiento; de lo contrario, se categorizaron según la mayoría de los lugares de nacimiento.

Resultados Clave

• Distribución de Modos: Aproximadamente el 30% (60 de 202) de la muestra pertenece al modo de baja velocidad, lo cual es consistente con estimaciones previas para púlsares jóvenes. El ~70% restante pertenece al modo de alta velocidad.
• Efectos de Sesgo de Selección: Como era de esperar, los púlsares de alta velocidad fueron encontrados, en promedio, más jóvenes (edades de frenado más cortas) y más distantes que los de baja velocidad. Estas diferencias fueron confirmadas mediante pruebas de Kolmogorov-Smirnov (K-S) ($p < 0.01$) y se atribuyen a efectos de selección observacional (dificultad para medir velocidades pequeñas a grandes distancias y objetos de alta velocidad saliendo del volumen de sondeo).
• Ancho de Pulso: No se encontró una diferencia significativa en el ancho de pulso ($W_{10}$ y $W_{50}$) entre los dos modos ($p \approx 0.5$), lo que sugiere que este parámetro no está sujeto a un fuerte sesgo de selección o dependencia intrínseca del modo.
• Correlación del Campo Magnético: Se detectó una diferencia estadísticamente significativa ($p \approx 0.02$) en las distribuciones del campo magnético superficial. Los púlsares en el modo de baja velocidad presentan una sobreabundancia de campos magnéticos más bajos ($B \lesssim 10^{12}$ G). Notablemente, no se identificaron púlsares con $B \lesssim 10^{11}$ G en el modo de alta velocidad. Esta discrepancia persiste incluso al contabilizar la distancia y la edad, aunque los autores señalan que puede estar parcialmente influenciada por el decaimiento del campo magnético a lo largo del tiempo.

Significancia y Reivindicaciones
El artículo afirma proporcionar una clasificación detallada de los púlsares de radio normales en grupos de velocidad de patada bimodal y analizar las distribuciones de parámetros resultantes. La principal significancia reside en la identificación de una correlación entre los modos de velocidad de patada y la intensidad del campo magnético. Los autores afirman que, si bien el sesgo de selección y el decaimiento del campo magnético pueden contribuir a la diferencia observada, el origen de la discrepancia no está completamente claro. Sugieren que, si la estructura de la patada bimodal surge de diferencias en la física de la explosión de la supernova, la diferencia del campo magnético podría compartir un origen similar, vinculando potencialmente los mecanismos de patada de las EN con el campo magnético dipolar externo. Sin embargo, los autores se mantienen modestos, reconociendo que se requiere un modelado de síntesis de poblaciones detallado para confirmar estos orígenes y que no se pueden descartar por completo los efectos de selección. El estudio no propone nuevas aplicaciones experimentales, sino que ofrece un análisis estadístico de datos existentes para refinar la comprensión de la cinemática y las propiedades magnéticas de las EN.