Stochastic and secular anomalies in pulsar braking indices

Andrés F. Vargas,Andrew Melatos
2024-10-07
Abstract:Stochastic and secular variations in the spin frequency $\nu$ of a rotation-powered pulsar complicate the interpretation of the measured braking index, $n$, in terms of a power-law spin-down torque $\propto \nu^{n_{\rm pl}}$. Both categories of variation can lead to anomalous braking indices, with $\vert n \vert = \vert \nu \ddot{\nu} / \dot{\nu}^2 \vert \gg 1$, where the overdot symbolizes a derivative with respect to time. Here we quantify the combined effect of stochastic and secular deviations from pure power-law spin down on measurements of $n$. Through analytic calculations, Monte Carlo simulations involving synthetic data, and modern Bayesian timing techniques, it is shown that the variance of $n$ satisfies the predictive, falsifiable formula $\langle n^{2} \rangle = (n_{\rm pl}+\dot{K}_{\rm dim})^{2}+\sigma_{\rm dim}^{2}$, where $\dot{K}_{\rm dim}$ is inversely proportional to the time-scale $\tau_K$ over which the proportionality constant of the power-law spin-down torque varies, $\sigma_{\rm dim}$ is proportional to the timing noise amplitude and inversely proportional to the square root of the total observing time, and the average is over an ensemble of random realizations of the timing noise process. The anomalous regime $\langle n^2 \rangle \gg 1$ occurs for $\dot{K}_{\rm dim} \gg 1$, $\sigma_{\rm dim} \gg 1$, or both. The sign of $n$ depends in part on the sign of $\dot{K}_{\rm dim}$, so it is possible to measure unequal numbers of positive and negative $n$ values in a large sample of pulsars. The distinguishable impact of stochastic and secular anomalies on phase residuals is quantified to prepare for extending the analysis of synthetic data to real pulsars.
High Energy Astrophysical Phenomena
What problem does this paper attempt to address?
The problem that this paper attempts to solve is the random and long - term anomalies in the measurement of the braking index ($n$) of rotation - powered pulsars. Specifically, the paper focuses on how these anomalies affect the interpretation of the braking index $n$ measured based on the power - law spin - down torque ($\dot{\nu} \propto \nu^{n_{\text{pl}}}$). These problems include: 1. **Random Fluctuations**: The spin frequency $\nu$ of pulsars has random fluctuations, which may cause the measured braking index $n$ to have abnormal values, i.e., $|n| \gg 1$. 2. **Long - Term Variations**: The proportionality factor $K(t)$ changes slowly over time, resulting in $\dot{\nu} = K(t)\nu^{n_{\text{pl}}}$, and this change can also lead to an abnormal braking index $n$. 3. **Combined Effect**: How the combined action of random fluctuations and long - term variations affects the measurement results of the braking index $n$. Through analytical calculations, Monte Carlo simulations, and modern Bayesian timing techniques, the paper quantifies the effects of these random and long - term anomalies on the measurement of the braking index $n$, and proposes a predictive formula to describe the variance $\langle n^2 \rangle$ of $n$. The specific formula is: \[ \langle n^2 \rangle=(n_{\text{pl}}+\Delta K_{\text{dim}})^2+\sigma^2_{\text{dim}} \] where: - $\Delta K_{\text{dim}}$ is inversely proportional to the change time scale $\tau_K$ of the proportionality factor $K(t)$. - $\sigma_{\text{dim}}$ is proportional to the timing noise amplitude and inversely proportional to the square root of the total observation time. The paper also explores the effects of these anomalies on the phase residuals in order to prepare for extending the synthetic data analysis to real pulsars.