Rotational evolution of the Vela pulsar during the 2016 glitch
Pith reviewed 2026-05-25 11:22 UTC · model grok-4.3
The pith
The 2016 Vela pulsar glitch rose in less than 12.6 seconds, exhibited a frequency overshoot with fast relaxation, and was preceded by a slowdown.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Using pulse-to-pulse timing data, the authors constrain the glitch rise time to less than 12.6 seconds with 90 percent . They report definitive evidence for a rotational-frequency overshoot and fast relaxation after the glitch. They also find evidence for a slowdown of the star's rotation immediately before the glitch. The overshoot matches predictions from some theoretical models of superfluid dynamics, while the preceding slowdown is interpreted as a possible trigger that builds a critical lag between the crust and superfluid.
What carries the argument
Parametric models of the star's rotation frequency evolution, including rise time, overshoot amplitude, and multiple relaxation timescales, fitted to individual pulse arrival times.
If this is right
- The short rise time restricts the allowed timescales for superfluid vortex motion or other internal transfer processes.
- The overshoot and fast relaxation support theoretical pictures in which angular momentum is exchanged in distinct stages.
- The pre-glitch slowdown indicates that a critical lag may build up and initiate the glitch.
- The measured relaxation timescales supply new constraints on how the neutron star returns to its steady spin-down rate.
Where Pith is reading between the lines
- Similar high-resolution monitoring of other glitching pulsars could test whether overshoots and pre-glitch slowdowns are common features.
- Continuous timing campaigns might eventually use the slowdown signature to anticipate when a glitch will occur.
- The approach could be extended to ask whether the same multi-timescale behavior appears in glitches of different sizes or in different pulsars.
Load-bearing premise
The observed timing residuals come solely from changes in the pulsar's rotation frequency and the chosen models capture all relevant pulse-to-pulse behavior without unmodeled noise or propagation effects.
What would settle it
A reanalysis or new observation that measures a glitch rise time longer than 12.6 seconds, or that finds no frequency overshoot when alternative noise models are used, would falsify the central claims.
Figures
read the original abstract
The 2016 Vela glitch observed by the Mt Pleasant radio telescope provides the first opportunity to study pulse-to-pulse dynamics of a pulsar glitch, opening up new possibilities to study the neutron star's interior. We fit models of the star's rotation frequency to the pulsar data, and present three new results. First, we constrain the glitch rise time to less than 12.6s with 90% confidence, almost three times shorter than the previous best constraint. Second, we find definitive evidence for a rotational-frequency overshoot and fast relaxation following the glitch. Third, we find evidence for a slow-down of the star's rotation immediately prior to the glitch. The overshoot is predicted theoretically by some models; we discuss implications of the glitch rise and overshoot decay times on internal neutron-star physics. The slow down preceding the glitch is unexpected; we propose the slow-down may trigger the glitch by causing a critical lag between crustal superfluid and the crust.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript analyzes high-cadence pulse arrival times from the 2016 Vela glitch observed at Mt Pleasant, fitting parametric models of rotational frequency evolution (instantaneous jump, exponential overshoot, power-law relaxation, and linear pre-glitch ramp) directly to the timing residuals. It reports three main results: an upper limit on glitch rise time of <12.6 s at 90% confidence, definitive evidence for a post-glitch frequency overshoot with fast relaxation, and evidence for a pre-glitch rotational slowdown, with discussion of implications for neutron-star superfluid dynamics.
Significance. If the modeling assumptions hold, the work supplies the tightest observational bound on glitch rise time to date and introduces new empirical features (overshoot decay and pre-glitch slowdown) that can be compared against theoretical predictions for vortex dynamics and crustal coupling. The pulse-to-pulse resolution of this event is a clear observational advance.
major comments (3)
- [Timing model and residual analysis] The three central claims rest on the assumption that all timing residuals arise solely from rotational frequency changes captured by the chosen parametric forms. The manuscript must demonstrate that post-fit residuals show no significant structure (e.g., via autocorrelation function or periodogram in the timing-analysis section) and that alternative functional families or unmodeled effects (DM variations, propagation) do not produce comparable improvements in fit quality.
- [Results on post-glitch evolution] The statement of 'definitive evidence' for the overshoot and fast relaxation requires quantitative model-selection metrics (Bayes factor or Δχ² with degrees-of-freedom correction) between the baseline glitch model and the model that includes the overshoot term; qualitative description alone is insufficient to support the quoted .
- [Pre-glitch analysis] The pre-glitch slowdown is modeled as a linear ramp; its statistical significance must be evaluated against a null hypothesis with no pre-glitch term, and the manuscript should test whether the feature persists under different data-selection windows or when timing noise is explicitly parameterized.
minor comments (2)
- [Abstract] The abstract uses 'definitive evidence'; this language should be moderated to reflect the model-dependent nature of the result.
- [Figures] Figure captions should explicitly state the number of free parameters in each plotted model and the data cadence used for the fits.
Simulated Author's Rebuttal
We thank the referee for the constructive report. The comments highlight important aspects of model validation and statistical rigor that we will address in revision. We respond to each major comment below.
read point-by-point responses
-
Referee: The three central claims rest on the assumption that all timing residuals arise solely from rotational frequency changes captured by the chosen parametric forms. The manuscript must demonstrate that post-fit residuals show no significant structure (e.g., via autocorrelation function or periodogram in the timing-analysis section) and that alternative functional families or unmodeled effects (DM variations, propagation) do not produce comparable improvements in fit quality.
Authors: We agree that explicit validation of the residuals is required. In the revised manuscript we will add an autocorrelation function and periodogram of the post-fit timing residuals in the timing-analysis section to demonstrate the absence of significant unmodeled structure. For unmodeled effects such as DM variations, the Mt Pleasant observations were made at 635 MHz and 990 MHz where dispersive delays are modest; we will nevertheless fit an additional DM parameter to the data and compare the resulting χ² improvement against the rotational models using the same information criteria. We expect this test to confirm that the parametric frequency models remain preferred. revision: yes
-
Referee: The statement of 'definitive evidence' for the overshoot and fast relaxation requires quantitative model-selection metrics (Bayes factor or Δχ² with degrees-of-freedom correction) between the baseline glitch model and the model that includes the overshoot term; qualitative description alone is insufficient to support the quoted confidence.
Authors: We accept that a quantitative model-selection statistic is necessary to support the strength of the claim. We will compute and report the Bayes factor (via nested sampling) and the Δχ² per degree of freedom between the baseline glitch model and the model that includes the overshoot-plus-fast-relaxation term. These metrics will be presented in the results section; if they exceed conventional thresholds for strong evidence we will retain the phrasing, otherwise we will adjust the language to reflect the precise statistical support. revision: yes
-
Referee: The pre-glitch slowdown is modeled as a linear ramp; its statistical significance must be evaluated against a null hypothesis with no pre-glitch term, and the manuscript should test whether the feature persists under different data-selection windows or when timing noise is explicitly parameterized.
Authors: We will add a formal model-comparison test (likelihood-ratio statistic or BIC difference) of the linear pre-glitch ramp against the null model that omits it. We will also repeat the fit using two alternative pre-glitch data windows and will introduce a simple power-law red-noise component to the timing model to check whether the ramp remains statistically significant. The outcomes of these robustness checks will be reported in the revised text. revision: yes
Circularity Check
No significant circularity; results from direct fits to new timing data
full rationale
The paper's central results (glitch rise-time upper limit, overshoot evidence, pre-glitch slowdown) are obtained by fitting parametric frequency-evolution models directly to the 2016 Mt Pleasant pulse arrival times. No quoted equations or claims reduce any reported quantity to a definition or fit performed in the authors' prior work; the derivation chain terminates at the external observational dataset rather than at any self-referential input. Self-citations, if present, are not load-bearing for the quoted constraints or evidence statements.
Axiom & Free-Parameter Ledger
Lean theorems connected to this paper
-
IndisputableMonolith/Foundation/RealityFromDistinction.leanreality_from_one_distinction unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
We fit models of the star's rotation frequency to the pulsar data... constrain the glitch rise time to less than 12.6 s... frequency overshoot and fast relaxation... slow-down... preceding the glitch
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
body-averaged model with two uncoupled spinning components... mutual-friction coefficient B
What do these tags mean?
- matches
- The paper's claim is directly supported by a theorem in the formal canon.
- supports
- The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
- extends
- The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
- uses
- The paper appears to rely on the theorem as machinery.
- contradicts
- The paper's claim conflicts with a theorem or certificate in the canon.
- unclear
- Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.
Reference graph
Works this paper leans on
-
[1]
Anderson, P . W. & Itoh, N. Pulsar glitches and restlessness as a hard superfluidity phenomenon. Nature 256, 25–27 (1975)
work page 1975
-
[2]
Palfreyman, J., Dickey, J. M., Hotan, A., Ellingsen, S. & van Straten, W. Alteration of the magnetosphere of the Vela pulsar during a glitch.Nature 556, 219–222 (2018)
work page 2018
-
[3]
Dodson, R. G., McCulloch, P . M. & Lewis, D. R. High Time Resolution Observations of the January 2000 Glitch in the Vela Pulsar. Astrophys. J. Lett. 564, L85–L88 (2002)
work page 2000
-
[4]
Dodson, R., Lewis, D. & McCulloch, P . Two decades of pulsar timing of Vela. Astrophys. Space Sci. 308, 585–589 (2007)
work page 2007
-
[5]
van Eysden, C. A. & Melatos, A. Pulsar glitch recovery and the superflu- idity coefficients of bulk nuclear matter. Mon. Notices Royal Astron. Soc. 409, 1253–1268 (2010)
work page 2010
-
[6]
Haskell, B., Pizzochero, P . M. & Sidery, T. Modelling pulsar glitches with realistic pinning forces: a hydrodynamical approach. Mon. Notices Royal Astron. Soc. 420, 658–671 (2012)
work page 2012
-
[7]
Antonelli, M. & Pizzochero, P . M. Axially symmetric equations for differ- ential pulsar rotation with superfluid entrainment. Mon. Notices Royal Astron. Soc. 464, 721–733 (2017)
work page 2017
-
[8]
Graber, V., Cumming, A. & Andersson, N. Glitch Rises as a Test for Rapid Superfluid Coupling in Neutron Stars. Astrophys. J. 865, 23 (2018)
work page 2018
-
[9]
Crust-breaking by neutron superfluids and the Vela pulsar glitches
Ruderman, M. Crust-breaking by neutron superfluids and the Vela pulsar glitches. Astrophys. J. 203, 213–222 (1976)
work page 1976
-
[10]
Cheng, K. S., Alpar, M. A., Pines, D. & Shaham, J. Spontaneous su- perfluid unpinning and the inhomogeneous distribution of vortex lines in neutron stars. Astrophys. J. 330, 835–846 (1988)
work page 1988
-
[11]
Alpar, M. A., Chau, H. F ., Cheng, K. S. & Pines, D. Postglitch Relaxation of the Crab Pulsar: Evidence for Crust Cracking. Astrophys. J. 427, L29 (1994)
work page 1994
-
[12]
Andersson, N., Comer, G. L. & Prix, R. Are Pulsar Glitches Triggered by a Superfluid Two-Stream Instability? Phys. Rev. Lett. 90, 091101 (2003)
work page 2003
- [14]
-
[15]
Warszawski, L. & Melatos, A. Knock-on processes in superfluid vortex avalanches and pulsar glitch statistics. Mon. Notices Royal Astron. Soc. 428, 1911–1926 (2013)
work page 1911
-
[16]
Andersson, N., Glampedakis, K. & Hogg, M. Superfluid instability of r- modes in “differentially rotating” neutron stars. Phys. Rev. D 87, 063007 (2013)
work page 2013
-
[17]
Hobbs, G. B., Edwards, R. T. & Manchester, R. N. TEMPO2, a new pulsar-timing package - I. An overview. Mon. Notices Royal Astron. Soc. 369, 655–672 (2006)
work page 2006
-
[18]
Edwards, R. T., Hobbs, G. B. & Manchester, R. N. TEMPO2, a new pulsar timing package - II. The timing model and precision estimates. Mon. Notices Royal Astron. Soc. 372, 1549–1574 (2006)
work page 2006
-
[19]
Ashton, G. et al. BILBY: A User-friendly Bayesian Inference Library for Gravitational-wave Astronomy. Astrophys. J. Supp. Ser. 241, 27 (2019)
work page 2019
-
[20]
Buchner, J. et al. X-ray spectral modelling of the AGN obscuring region in the CDFS: Bayesian model selection and catalogue. Astron. Astrophys. 564, A125 (2014)
work page 2014
-
[21]
Feroz, F . & Hobson, M. P . Multimodal nested sampling: an efficient and robust alternative to Markov Chain Monte Carlo methods for astronomical data analyses. Mon. Notices Royal Astron. Soc. 384, 449–463 (2008)
work page 2008
-
[22]
Feroz, F ., Hobson, M. P . & Bridges, M. MULTINEST: an efficient and robust Bayesian inference tool for cosmology and particle physics. Mon. Notices Royal Astron. Soc. 398, 1601–1614 (2009)
work page 2009
-
[23]
Sidery, T., Passamonti, A. & Andersson, N. The dynamics of pulsar glitches: contrasting phenomenology with numerical evolutions. Mon. Notices Royal Astron. Soc. 405, 1061–1074 (2010)
work page 2010
-
[24]
Alpar, M. A., Langer, S. A. & Sauls, J. A. Rapid postglitch spin-up of the superfluid core in pulsars. Astrophys. J. 282, 533–541 (1984)
work page 1984
-
[25]
Superfluid hydrodynamics in rotating neutron stars
Mendell, G. Superfluid hydrodynamics in rotating neutron stars. I - Nondissipative equations. II - Dissipative effects.Astrophys. J. 380, 515– 540 (1991)
work page 1991
-
[26]
Andersson, N., Sidery, T. & Comer, G. L. Mutual friction in superfluid neutron stars. Mon. Notices Royal Astron. Soc. 368, 162–170 (2006)
work page 2006
-
[27]
Kass, R. & Raftery, A. Bayes factors. Journal of the American Statistical Association 90, 773–795 (2015)
work page 2015
-
[28]
Cordes, J. M. & Shannon, R. M. A Measurement Model for Precision Pulsar Timing. arXiv e-prints (2010). 1010.3785
work page internal anchor Pith review Pith/arXiv arXiv 2010
-
[29]
Dib, R., Kaspi, V. M. & Gavriil, F . P . Glitches in Anomalous X-Ray Pulsars. Astrophys. J. 673, 1044–1061 (2008)
work page 2008
-
[30]
Weltevrede, P ., Johnston, S. & Espinoza, C. M. The glitch-induced iden- tity changes of PSR J1119-6127. Mon. Notices Royal Astron. Soc. 411, 1917–1934 (2011)
work page 1917
-
[31]
Archibald, R. F ., Kaspi, V. M., Tendulkar, S. P . & Scholz, P . A Magnetar- like Outburst from a High-B Radio Pulsar. Astrophys. J. 829, L21 (2016)
work page 2016
-
[32]
Alpar, M. A., Anderson, P . W., Pines, D. & Shaham, J. Giant glitches and pinned vorticity in the VELA and other pulsars. Astrophys. J. 249, L29–L33 (1981)
work page 1981
-
[33]
Link, B. K. & Epstein, R. I. Mechanics and Energetics of Vortex Unpinning in Neutron Stars. Astrophys. J. 373, 592 (1991)
work page 1991
-
[34]
Glampedakis, K. & Andersson, N. Hydrodynamical Trigger Mechanism for Pulsar Glitches. Phys. Rev. Lett. 102, 141101 (2009)
work page 2009
-
[35]
Pizzochero, P . M. Angular Momentum Transfer in Vela-like Pulsar Glitches. Astrophys. J. 743, L20 (2011)
work page 2011
-
[36]
Melatos, A., Peralta, C. & Wyithe, J. S. B. Avalanche Dynamics of Radio Pulsar Glitches. Astrophys. J. 672, 1103–1118 (2008)
work page 2008
-
[37]
Fuentes, J. R. et al. The glitch activity of neutron stars.Astron. Astrophys. 608, A131 (2017)
work page 2017
-
[38]
Howitt, G., Melatos, A. & Delaigle, A. Nonparametric Estimation of the Size and Waiting Time Distributions of Pulsar Glitches.Astrophys. J. 867, 60 (2018). Correspondence and requests for materialsshould we addressed to G.A, gregory.ashton@ligo.org. Acknowledgements We are grateful to Andrew Melatos, Ian Jones, and the anonymous reviewers for valuable com...
work page 2018
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.