On the braking index evolution of PSR B0540-69: wind braking of pulsars

Resumen Técnico: Sobre la evolución del índice de frenado de PSR B0540-69: frenado por viento de púlsares

Planteamiento del problema
El púlsar PSR B0540-69 (referido como PSR B0540-69 en el texto, aunque el resumen contiene un error tipográfico como PSR B0450-69) presenta un caso de estudio único en la física de púlsares porque posee tanto un índice de frenado medido como un cambio documentado en su estado de frenado (spin-down). Tras un cambio de estado donde su tasa de frenado aumentó un 36%, el púlsar exhibió una evolución dependiente del tiempo de su índice de frenado ($n$). Un trabajo teórico previo (Paper I, Kou et al. 2016) predijo que, si este cambio de estado fuera impulsado por un viento de partículas incrementado, el índice de frenado en el nuevo estado de alto frenado debería estabilizarse en un valor inferior a su valor de estado persistente ($n=2.13$), específicamente alrededor de $n \approx 1.79$ para un modelo de brecha de vacío. Sin embargo, los datos observacionales muestran que el índice de frenado evoluciona desde valores muy bajos ($n \approx 0.03$) hasta $n \approx 1.6$, aproximándose aparentemente a la predicción teórica. El problema central abordado es explicar la evolución transitoria del índice de frenado durante esta transición y determinar el valor del estado estacionario final para restringir el mecanismo de frenado subyacente.

Metodología
El artículo emplea el modelo de frenado por viento, un marco fenomenológico donde la pérdida de energía rotacional de la estrella de neutrones se atribuye tanto a la radiación de dipolo magnético como al flujo de viento de partículas. El modelo define un parámetro adimensional $\eta$ que combina el ángulo de inclinación ($\alpha$), la densidad numérica de partículas relativa a la densidad de Goldreich-Julian ($\kappa$) y el potencial de aceleración de partículas ($\Delta\phi$).

Los autores proponen que la evolución observada del índice de frenado no se debe a un cambio en la tasa media de frenado (que permanece en el estado alto), sino más bien a una pequeña variación transitoria en la densidad de partículas. Específicamente, el artículo asume que la densidad de partículas $\kappa$ en el nuevo estado contiene una perturbación decreciente:
$$\kappa(t) = \kappa_1(1 + \delta e^{-t/\tau})$$
donde $\delta$ es una amplitud pequeña y $\tau$ es una escala de tiempo de decaimiento. Al tomar la derivada temporal de la ecuación de frenado, los autores derivan cómo esta variación transitoria de la densidad de partículas contribuye a la segunda derivada de la velocidad angular ($\ddot{\Omega}$), lo que afecta directamente al índice de frenado medido. El índice de frenado total se modela como la suma del valor de estado estacionario ($n_{new}$) y un término exponencial transitorio:
$$n(t) = n_{new} + \Delta n e^{-t/\tau}$$

Contribuciones clave y resultados

1. Modelado de la evolución transitoria: El artículo deriva que una pequeña variación exponencialmente decreciente en la densidad de partículas ($\delta \ll 1$) puede producir una variación transitoria significativa en el índice de frenado si la relación entre la edad característica del púlsar y la escala de tiempo de decaimiento ($\tau_c/\tau$) es grande. Esto explica por qué la tasa de frenado permanece estable mientras que el índice de frenado evoluciona.
2. Ajuste cuantitativo: Ajustando este modelo a los datos observacionales (Ho et al. 2026), los autores encuentran que la evolución del índice de frenado sigue una forma de decaimiento exponencial. Los parámetros de mejor ajuste son una amplitud transitoria inicial de $\Delta n \approx -2.1$ y una escala de tiempo de decaimiento de $\tau \approx 1708$ días. Esto corresponde a una amplitud de variación de la densidad de partículas de $\delta \sim 10^{-3}$.
3. Convergencia a la predicción teórica: El modelo predice que el índice de frenado se estabilizará en un valor de estado estacionario constante en el futuro cercano. Este valor final se predice como $n_{new} \approx 1.79$ (basado en el modelo de brecha de vacío), el cual es inferior al valor de estado persistente ($n=2.13$) y es consistente con la hipótesis de frenado por viento.
4. Diferenciación de modelos alternativos: El artículo argumenta que las explicaciones alternativas para la evolución del índice de frenado —como la evolución de los campos magnéticos ($\dot{B}$), el momento de inercia ($\dot{I}$) o los ángulos de inclinación ($\dot{\alpha}$)— carecen de poder predictivo respecto al valor del estado estacionario final. En contraste, el modelo de frenado por viento proporciona una predicción específica y contrastable para el índice de frenado final basada en el potencial de aceleración de partículas.

Significancia
El artículo sostiene que la evolución del índice de frenado de PSR B0540-69 sirve como una herramienta de diagnóstico crítica para los mecanismos de frenado de púlsares. Al demostrar que la evolución es consistente con una variación transitoria en la densidad de partículas dentro del marco de frenado por viento, el estudio respalda la hipótesis de que el cambio de estado de frenado es causado por un viento de partículas incrementado, similar al mecanismo observado en púlsares intermitentes.

La significancia principal radica en la contrastabilidad futura del modelo. El artículo postula que el monitoreo continuo de PSR B0540-69 eventualmente revelará el índice de frenado de estado estacionario. Si el valor medido converge al valor predicho de $n \approx 1.79$ (o otro valor determinado por el potencial de aceleración específico), confirmará inequívocamente el mecanismo de frenado por viento y restringirá el potencial de aceleración de partículas en la magnetosfera del púlsar. Esto distinguiría al modelo de frenado por viento de las teorías competidoras que involucran la evolución del campo magnético o geométrica, las cuales no ofrecen tales restricciones específicas sobre el índice de frenado final.