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

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

Problem Statement
The pulsar PSR B0540-69 (referred to as PSR B0540-69 in the text, though the abstract contains a typo as PSR B0450-69) presents a unique case study in pulsar physics because it possesses both a measured braking index and a documented change in its spin-down state. Following a state change where its spin-down rate increased by 36%, the pulsar exhibited a time-dependent evolution of its braking index ($n$). Previous theoretical work (Paper I, Kou et al. 2016) predicted that if this state change were driven by an enhanced particle wind, the braking index in the new high spin-down state should settle at a value lower than its persistent state value ($n=2.13$), specifically around $n \approx 1.79$ for a vacuum gap model. However, observational data shows the braking index evolving from very low values ($n \approx 0.03$) up to $n \approx 1.6$, seemingly approaching the theoretical prediction. The central problem addressed is explaining the transient evolution of the braking index during this transition and determining the final steady-state value to constrain the underlying braking mechanism.

Methodology
The paper employs the wind braking model, a phenomenological framework where the neutron star's rotational energy loss is attributed to both magnetic dipole radiation and particle wind outflow. The model defines a dimensionless parameter $\eta$ that combines the inclination angle ($\alpha$), particle number density relative to the Goldreich-Julian density ($\kappa$), and the particle acceleration potential ($\Delta\phi$).

The authors propose that the observed evolution of the braking index is not due to a change in the mean spin-down rate (which remains in the high state) but rather a small, transient variation in the particle number density. Specifically, the paper assumes the particle density $\kappa$ in the new state contains a decaying perturbation:
$$\kappa(t) = \kappa_1(1 + \delta e^{-t/\tau})$$
where $\delta$ is a small amplitude and $\tau$ is a decay timescale. By taking the time derivative of the spin-down equation, the authors derive how this transient density variation contributes to the second derivative of the angular velocity ($\ddot{\Omega}$), which directly affects the measured braking index. The total braking index is modeled as the sum of the steady-state value ($n_{new}$) and a transient exponential term:
$$n(t) = n_{new} + \Delta n e^{-t/\tau}$$

Key Contributions and Results

1. Modeling the Transient Evolution: The paper derives that a small, exponentially decaying variation in particle density ($\delta \ll 1$) can produce a significant transient variation in the braking index if the ratio of the pulsar's characteristic age to the decay timescale ($\tau_c/\tau$) is large. This explains why the spin-down rate remains stable while the braking index evolves.
2. Quantitative Fit: Fitting this model to observational data (Ho et al. 2026), the authors find that the braking index evolution follows an exponential decay form. The best-fit parameters are an initial transient amplitude $\Delta n \approx -2.1$ and a decay timescale $\tau \approx 1708$ days. This corresponds to a particle density variation amplitude of $\delta \sim 10^{-3}$.
3. Convergence to Theoretical Prediction: The model predicts that the braking index will settle to a constant steady-state value in the near future. This final value is predicted to be $n_{new} \approx 1.79$ (based on the vacuum gap model), which is lower than the persistent state value ($n=2.13$) and consistent with the wind braking hypothesis.
4. Differentiation from Alternative Models: The paper argues that alternative explanations for the braking index evolution—such as evolving magnetic fields ($\dot{B}$), moment of inertia ($\dot{I}$), or inclination angles ($\dot{\alpha}$)—lack predictive power regarding the final steady-state value. In contrast, the wind braking model provides a specific, testable prediction for the final braking index based on the particle acceleration potential.

Significance
The paper claims that the braking index evolution of PSR B0540-69 serves as a critical diagnostic tool for pulsar braking mechanisms. By demonstrating that the evolution is consistent with a transient variation in particle density within the wind braking framework, the study supports the hypothesis that the spin-down state change is caused by an enhanced particle wind, similar to the mechanism observed in intermittent pulsars.

The primary significance lies in the future testability of the model. The paper posits that continued monitoring of PSR B0540-69 will eventually reveal the steady-state braking index. If the measured value converges to the predicted $n \approx 1.79$ (or another value determined by the specific acceleration potential), it will unambiguously confirm the wind braking mechanism and constrain the particle acceleration potential in the pulsar magnetosphere. This would distinguish the wind braking model from competing theories involving magnetic field or geometric evolution, which do not offer such specific constraints on the final braking index.