pith. sign in

arxiv: 1907.01124 · v1 · pith:OS4YF6RBnew · submitted 2019-07-02 · 🌌 astro-ph.HE · gr-qc

Rotational evolution of the Vela pulsar during the 2016 glitch

Pith reviewed 2026-05-25 11:22 UTC · model grok-4.3

classification 🌌 astro-ph.HE gr-qc
keywords Vela pulsarpulsar glitchneutron starglitch rise timerotational overshootsuperfluidtiming residuals
0
0 comments X

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.

The paper uses high time-resolution radio observations of the Vela pulsar to fit models directly to the pulse arrival times during its 2016 glitch. It reports a firm upper bound on the time taken for the frequency jump and identifies an overshoot in rotation rate followed by quick recovery. The data also show a brief slowdown in rotation just before the main event. These measurements test how angular momentum is redistributed inside the neutron star on short timescales. A reader would care because the results directly limit the possible physical processes occurring in the star's interior during glitches.

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

These are editorial extensions of the paper, not claims the author makes directly.

  • 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

Figures reproduced from arXiv: 1907.01124 by Gregory Ashton, Jim Palfreyman, Paul D. Lasky, Vanessa Graber.

Figure 1
Figure 1. Figure 1: Furthermore, it is one of many simple phenomenological mod [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
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.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

3 major / 2 minor

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)
  1. [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.
  2. [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 .
  3. [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)
  1. [Abstract] The abstract uses 'definitive evidence'; this language should be moderated to reflect the model-dependent nature of the result.
  2. [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

3 responses · 0 unresolved

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
  1. 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

  2. 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

  3. 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

0 steps flagged

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

0 free parameters · 0 axioms · 0 invented entities

The central claims rest on the assumption that the observed timing variations are produced by rotational-frequency changes that can be captured by a small number of parametric glitch models; no free parameters, axioms, or invented entities are enumerated in the abstract.

pith-pipeline@v0.9.0 · 5705 in / 1182 out tokens · 21170 ms · 2026-05-25T11:22:00.072674+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

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

37 extracted references · 37 canonical work pages · 1 internal anchor

  1. [1]

    Anderson, P . W. & Itoh, N. Pulsar glitches and restlessness as a hard superfluidity phenomenon. Nature 256, 25–27 (1975)

  2. [2]

    M., Hotan, A., Ellingsen, S

    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)

  3. [3]

    G., McCulloch, P

    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)

  4. [4]

    & McCulloch, P

    Dodson, R., Lewis, D. & McCulloch, P . Two decades of pulsar timing of Vela. Astrophys. Space Sci. 308, 585–589 (2007)

  5. [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)

  6. [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)

  7. [7]

    & Pizzochero, P

    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)

  8. [8]

    & Andersson, N

    Graber, V., Cumming, A. & Andersson, N. Glitch Rises as a Test for Rapid Superfluid Coupling in Neutron Stars. Astrophys. J. 865, 23 (2018)

  9. [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)

  10. [10]

    S., Alpar, M

    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)

  11. [11]

    A., Chau, H

    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)

  12. [12]

    Andersson, N., Comer, G. L. & Prix, R. Are Pulsar Glitches Triggered by a Superfluid Two-Stream Instability? Phys. Rev. Lett. 90, 091101 (2003)

  13. [14]

    & Ooi, A

    Peralta, C., Melatos, A., Giacobello, M. & Ooi, A. Transitions between Turbulent and Laminar Superfluid Vorticity States in the Outer Core of a Neutron Star. Astrophys. J. 651, 1079–1091 (2006)

  14. [15]

    & Melatos, A

    Warszawski, L. & Melatos, A. Knock-on processes in superfluid vortex avalanches and pulsar glitch statistics. Mon. Notices Royal Astron. Soc. 428, 1911–1926 (2013)

  15. [16]

    differentially rotating

    Andersson, N., Glampedakis, K. & Hogg, M. Superfluid instability of r- modes in “differentially rotating” neutron stars. Phys. Rev. D 87, 063007 (2013)

  16. [17]

    B., Edwards, R

    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)

  17. [18]

    T., Hobbs, G

    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)

  18. [19]

    Ashton, G. et al. BILBY: A User-friendly Bayesian Inference Library for Gravitational-wave Astronomy. Astrophys. J. Supp. Ser. 241, 27 (2019)

  19. [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)

  20. [21]

    & Hobson, M

    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)

  21. [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)

  22. [23]

    & Andersson, N

    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)

  23. [24]

    A., Langer, S

    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)

  24. [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)

  25. [26]

    & Comer, G

    Andersson, N., Sidery, T. & Comer, G. L. Mutual friction in superfluid neutron stars. Mon. Notices Royal Astron. Soc. 368, 162–170 (2006)

  26. [27]

    & Raftery, A

    Kass, R. & Raftery, A. Bayes factors. Journal of the American Statistical Association 90, 773–795 (2015)

  27. [28]

    Cordes, J. M. & Shannon, R. M. A Measurement Model for Precision Pulsar Timing. arXiv e-prints (2010). 1010.3785

  28. [29]

    Dib, R., Kaspi, V. M. & Gavriil, F . P . Glitches in Anomalous X-Ray Pulsars. Astrophys. J. 673, 1044–1061 (2008)

  29. [30]

    & Espinoza, C

    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)

  30. [31]

    F ., Kaspi, V

    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)

  31. [32]

    A., Anderson, P

    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)

  32. [33]

    Link, B. K. & Epstein, R. I. Mechanics and Energetics of Vortex Unpinning in Neutron Stars. Astrophys. J. 373, 592 (1991)

  33. [34]

    & Andersson, N

    Glampedakis, K. & Andersson, N. Hydrodynamical Trigger Mechanism for Pulsar Glitches. Phys. Rev. Lett. 102, 141101 (2009)

  34. [35]

    Pizzochero, P . M. Angular Momentum Transfer in Vela-like Pulsar Glitches. Astrophys. J. 743, L20 (2011)

  35. [36]

    & Wyithe, J

    Melatos, A., Peralta, C. & Wyithe, J. S. B. Avalanche Dynamics of Radio Pulsar Glitches. Astrophys. J. 672, 1103–1118 (2008)

  36. [37]

    Fuentes, J. R. et al. The glitch activity of neutron stars.Astron. Astrophys. 608, A131 (2017)

  37. [38]

    & Delaigle, A

    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...