REVIEW 4 major objections 6 minor 27 references
Investigation into Pulsar Glitch Parameters
T0 review · 4 major / 6 minor · reviewed 2026-08-12 · deepseek-v4-flash
Pith's one-line read Glitch rate correlates weakly with pulsar spin and falls with age.
desk verdict A modest but honest statistical update of the pulsar glitch-rate correlations; worth engaging with once the exposure proxy and missing error estimates are fixed. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The carrying tool is Pearson's correlation coefficient ($r$, a standard measure of linear association from -1 to 1), applied to the quantity glitch rate $N/T_{\rm obs}$, where $N$ is the number of glitches recorded for a pulsar and $T_{\rm obs}$ is the time from the pulsar's reported discovery year to November 8, 2024. Each relation is computed twice, once on raw values and once after logarithmic transformation, and the log-log versions spread out the small-value tail that dominates the sample. The two-scale comparison is what lets the paper claim that the relationships are weak on a linear scale while still showing a clearer inverse age trend on a log-log scale.
What would settle it
If the same correlations are recomputed using each pulsar's actual monitored span (first to last observation) instead of discovery-to-2024, and the six Pearson coefficients move far from 0.03, 0.13, 0.12, 0.40, -0.13, and -0.42—for instance into the strong range above 0.7—then the paper's weak-correlation conclusion would be disproved.
Extended reading notes
Core claim
The central claim is that, on a sample of 215 pulsars with 677 glitch events, glitch rate is only weakly linearly related to the pulsar's rotational frequency and absolute spin-down rate, and weakly inversely related to characteristic age. In the paper's own numbers: against rotational frequency, Pearson's $r$ is 0.03 for glitch rate versus log frequency and 0.13 for log glitch rate versus log frequency; against absolute spin-down rate, $r$ is 0.12 and 0.40; against characteristic age, $r$ is -0.13 and -0.42. The paper concludes that changes in rotational frequency or spin-down rate do not strongly predict changes in glitch rate, and that glitch rate tends to decrease as characteristic age increases. This is presented as an updated confirmation of the earlier result that middle-aged pulsars glitch most frequently.
Load-bearing premise
The analysis assumes that every pulsar in the sample was monitored continuously and completely from its discovery to November 8, 2024, so every glitch (a sudden spin-up) that occurred was recorded; if coverage was patchy or sensitivity varied, the glitch rates and all six correlation coefficients are biased.
Editorial extensions
If this is right
- If the weak correlations hold, glitch rate is not predictable from rotational frequency alone, so single-parameter spin-based glitch models are insufficient.
- The modest log-log spin-down relation ($r=0.40$) allows a mild tendency for faster-slowing pulsars to glitch more often, but with large scatter.
- The inverse age relation supports the earlier picture that middle-aged pulsars glitch most and older pulsars glitch less.
- Extending the analysis to the full 2024 glitch catalog preserves the earlier age trend, so the middle-aged peak is not a small-sample artifact.
- Near-zero raw correlations imply future predictive work should search for thresholds, subpopulations, or nonlinear laws rather than simple linear scaling.
Reading between the lines
- A coverage-weighted reanalysis that replaces discovery-to-2024 time with each pulsar's actual monitored span could separate a true ageing trend from observing-history effects, because early monitoring of older pulsars was sparser.
- The weak spin-down correlations might sharpen inside subpopulations such as young pulsars with high spin-down, so splitting the sample by characteristic age or glitch size is a test the paper does not run.
- If the inverse age relation is physical, glitch rate could become a rough age indicator for isolated neutron stars, and the paper's suggestion of higher-order polynomial fits is a natural first step toward that clock.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper compiles an updated sample of 215 glitching pulsars and 677 glitches from the ATNF pulsar catalog (retrieved 8 November 2024), computes glitch rates as R = N/T_obs with T_obs = 2024.8552 minus the discovery year, and plots these rates against rotational frequency, absolute spin-down rate, and characteristic age. Pearson correlation coefficients are reported on both linear and log-transformed scales: r = 0.03/0.13 for frequency, 0.12/0.40 for absolute spin-down rate, and -0.13/-0.42 for characteristic age. The paper concludes that glitch rate is only weakly correlated with rotational frequency and spin-down rate and that it decreases weakly with characteristic age, consistent with Espinoza et al. (2011) and Millhouse et al. (2022).
Significance. If the correlations are robust, the paper provides a useful descriptive update on a larger sample than previous glitch-rate studies, and it explicitly places its results in the context of Espinoza et al. (2011) and Millhouse et al. (2022). The main strength is the larger, more recent catalog, which allows the earlier qualitative trends to be re-examined. However, the absence of uncertainties, the exposure definition, and the selection of only glitching pulsars currently limit the strength of the conclusions. With proper sensitivity analysis and error estimation, the paper could be a compact research note; in its present form it does not establish the quantitative claims beyond the raw descriptive correlations.
major comments (4)
- [Section 2 (Methodology)] The glitch rate is defined as R = N/T_obs with T_obs = 2024.8552 - discovery_year. This equates observational exposure with calendar time since discovery, and assumes that every pulsar was monitored continuously and completely since its discovery. In practice monitoring cadence, sensitivity, and glitch-detection methods vary across telescopes and epochs, and early data are often sparser. For a pulsar discovered in 1980 but regularly timed only from 2000, this definition overstates the exposure while N undercounts missed glitches, biasing R low. Because older pulsars tend to have longer and sparser monitoring histories, this bias can produce a spurious negative correlation with characteristic age, which is exactly the paper's headline result (r = -0.42 for log R vs log tau_c). The authors should either use actual monitoring intervals or demonstrate robustness by repeating the analysis on a uniformly observed subsample; at a minimum the assumption should be stated as a limitation and tested.
- [Section 4 (Discussions)] No uncertainties, p-values, or confidence intervals are reported for any Pearson r value. With n = 215, the difference between r = 0.03 and r = 0.13 is not qualitatively meaningful without error bars, and the rate variable is a count divided by a noisy exposure, so the error distribution is heteroscedastic and non-normal. The paper should provide bootstrap confidence intervals or a method that accounts for count noise (e.g., Poisson regression or a nonparametric rank correlation). Without this, the central distinction between 'weak' and 'no' correlation is not quantitatively established.
- [Section 2 and Section 4] The sample is restricted to pulsars with at least one recorded glitch, so the computed rate is conditional on having glitched. This selection can distort correlations, especially for pulsars with short exposure times, where a single glitch yields a very high inferred rate. The paper claims consistency with Millhouse et al. (2022), who included non-glitching pulsars, but the different sample selection means the comparison is not direct. The authors should quantify this selection effect or include non-glitching pulsars with appropriate upper limits on their rates.
- [Data Availability and Figures] The manuscript references 'Table 2' and states that a spreadsheet 'can be accessed here', but the version I inspected contains no table, URL, or access link, and Figures 1-3 appear only as captions. The exact ATNF catalog query parameters and version are not specified, and no code is provided. This makes the analysis unreproducible. Please include the full dataset (or a stable DOI) and the exact steps used for querying and processing the data.
minor comments (6)
- [Abstract] The phrase 'each glitch at least once' is ungrammatical; it should read 'each pulsar having glitched at least once', and 'inversely weak' would be clearer as 'weakly inversely correlated'.
- [Section 1 (Introduction)] The sentence '102 pulsars that glitched at least, including radio, X-ray' is missing the word 'once' and should be completed.
- [Section 2 (Methodology)] The spin-down rate is denoted F1 with units s^-2 in the query list, but in the results and discussion it is called 'absolute spin-down rate' without consistently showing that the absolute value is used; this should be stated explicitly.
- [Figures] The figure captions do not state that error bars are absent, and the axis labels should include units for both variables; additionally, the caption style switches between 'Figure 1' and 'Figures 1'.
- [Section 4 (Discussions)] The interpretation of r = 0.40 as 'weak' is subjective; the paper should either reference an effect-size convention or provide context from previous studies comparing similar correlation magnitudes.
- [Conclusion and Acknowledgments] The recommendation to 'explore polynomials of higher ranks' is vague and is not tied to any diagnostic evidence in the paper; it should be removed or made specific. The acknowledgment thanking 'the anonymous referee' is inappropriate before the review process and should be deleted.
Circularity Check
No circularity: glitch rates and correlation coefficients are computed directly from catalog data, with no fitted parameter or self-cited premise standing in for a result.
full rationale
The paper's load-bearing claim is that Pearson correlations between glitch rate and rotational frequency, spin-down rate, and characteristic age are weak. These correlations are computed by a direct formula: glitch rate equals N divided by Tobs, with N taken from the ATNF glitch counts and Tobs defined as the interval from the discovery year to 2024.8552. No parameter in this calculation is fitted from the target correlations, and no result is assumed in the inputs. The only external methodological choice, the Tobs definition, is attributed to Millhouse et al. (2022), which is not the present author's work and is not used to assert a conclusion. The conclusion that spin-down rate and frequency do not strongly predict glitch rate is an arithmetic summary of the computed r values, not a quantity that was inserted into the data before the analysis. The weak negative age correlation is also read directly from the plotted log-transformed data. There is no self-definitional step, no fitted input relabeled as prediction, no load-bearing self-citation, no imported uniqueness theorem, no ansatz smuggled in via citation, and no renaming of a known result. The possible observational bias in using discovery-to-now time as an exposure proxy is a data-quality and selection-effect concern, not circular reasoning, because the paper does not assume its conclusion when defining Tobs. The sample is restricted to pulsars that glitched at least once, which may bias rate estimates, but that is a sampling limitation rather than a circular derivation. Therefore no circular step is identified and the score is 0.
Assumptions & free parameters
free parameters (1)
- Observation end date =
2024.8552 (MJD 57789.5)
assumptions (4)
- domain assumption ATNF catalog glitch counts and discovery years are accurate and complete for the selected pulsars.
- domain assumption Time since discovery is a valid monitoring time for detecting glitches.
- standard math Pearson correlation is an appropriate measure of relationship strength in this context.
- domain assumption Only pulsars that have glitched at least once form the sample.
Cite this review
Pith. "Pith review of Investigation into Pulsar Glitch Parameters." pith.science (2026). https://pith.science/paper/3KHJZZV4
@misc{pith2026241112139,
author = {Pith},
title = {Pith review of: Investigation into Pulsar Glitch Parameters},
year = {2026},
howpublished = {\url{https://pith.science/paper/3KHJZZV4}},
note = {Machine review of arXiv:2411.12139}
}
read the original abstract
An updated analysis of pulsar glitch parameters was conducted using a sample of 215 pulsars, encompassing 677 recorded glitch events, each glitch at least once between 1968 and 2024. The glitch rates were estimated and plotted against various pulsar parameters, including rotational frequency, spin-down rate, and characteristic age. These relationships were analyzed using Pearson's correlation coefficient. The results indicate linear and weak relationships between glitch rate and the pulsar parameters. However, the relationship between glitch rate and characteristic age is inversely weak.
Reference graph
Works this paper leans on
-
[1]
Andersson, N., Glampedakis, K., Ho, W. C., & Espinoza, C. M. (2012). Pulsar glitches: The crust is not enough. Physical review letters, 109(24), 241103, https://journals.aps.org/prl/abstract/10.1103/PhysRevLett.109.241103
-
[2]
Baade, W., & Zwicky, F. (1934). On super-novae. Proceedings of the National Academy of Sciences, 20(5), 254-259, https://www.pnas.org/doi/abs/10.1073/pnas.20.5.254
-
[3]
Basu, A., Shaw, B., Antonopoulou, D., Keith, M. J., Lyne, A. G., Mickaliger, M. B., ... & Jordan, C. A. (2022). The Jodrell bank glitch catalogue: 106 new rotational glitches in 70 pulsars. Monthly Notices of the Royal Astronomical Society, 510(3), 4049-4062, https://academic.oup.com/mnras/article/510/3/4049/6440177
work page 2022
-
[4]
Bhattacharya, D., & van den Heuvel, E. P. J. (1991). Formation and evolution of binary and millisecond radio pulsars. Physics Reports, 203(1-2), 1-124, https://www.sciencedirect.com/science/article/abs/pii/037015739190064S
arXiv 1991
-
[5]
Burnell, J. B., & Bell, J. (1977). Little green men, white dwarfs or pulsars. Annals of the New York Academy of Science, 302, 685-689, http://www.bigear.org/CSMO/PDF/CS01/cs01p16.pdf
work page 1977
-
[6]
Burnell, S. J. (2004). An introduction to the sun and stars. Cambridge University Press, https://books.google.com.ng/books?hl=en&lr=&id=lb5owLGIQGsC&oi=fnd&pg=PA1&dq=Gre en,+Simon+F.%3B+Jones,+Mark+H.%3B+Burnell,+S.+Jocelyn+(2004).+An+Introduction+to+t he+Sun+and+Stars+(illustrated+ed.).+Cambridge+University+Press.+p.+322.+ISBN+978-0-521- 54622-5.+Archive...
work page 2004
-
[7]
Espinoza, C. M., Lyne, A. G., Stappers, B. W., & Kramer, M. (2011). A study of 315 glitches in the rotation of 102 pulsars. Monthly Notices of the Royal Astronomical Society, 414, pp. 1679-1704, https://academic.oup.com/mnras/article/414/2/1679/978848
work page 2011
-
[8]
Eya, I. O., Alhassan, J. A., Iyida, E. U., Chukwude, A. E., & Urama, J. O. (2022). On the pulsar spin frequency derivatives and the glitch activity. Astrophysics and Space Science, 367(3), 28, https://link.springer.com/article/10.1007/s10509-022-04059-0
Show all 27 references
-
[9]
O., Urama, J
Eya, I. O., Urama, J. O., & Chukwude, A. E. (2019). Entrainment factor of individual glitch fractional moment of inertia. Revista mexicana de astronomía y astrofísica, 55(1), 3-9. https://www.scielo.org.mx/scielo.php?pid=S0185- 11012019000100003&script=sci_arttext&tlng=en
2019
-
[10]
O., & Urama, J
Eya, I. O., & Urama, J. O. (2014). Statistical study of neutron star glitches, https://www.researchgate.net/profile/Innocent- Eya/publication/262637319_Statistical_Study_of_Neutron_Star_Glitches/links/02e7e5384a624b1 49f000000/Statistical-Study-of-Neutron-Star-Glitches.pdf
2014
-
[11]
Kutschera, M. (1998). Neutron stars: Formation and structure. arXiv preprint astro-ph/9801235
1998 arXiv
-
[12]
Lattimer, J. M. (2015, February). Introduction to neutron stars. In AIP Conference Proceedings (Vol. 1645, No. 1, pp. 61-78). American Institute of Physics, https://pubs.aip.org/aip/acp/article/1645/1/61/852737/Introduction-to-neutron-stars
2015
-
[13]
I., & Lattimer, J
Link, B., Epstein, R. I., & Lattimer, J. M. (1999). Pulsar constraints on neutron star structure and equation of state. Physical Review Letters, 83(17), 3362, https://journals.aps.org/prl/abstract/10.1103/PhysRevLett.83.3362
1999 doi
-
[14]
G., Shemar, S
Lyne, A. G., Shemar, S. L., & Graham Smith, F. (2000). Statistical studies of pulsar glitches. Monthly Notices of the Royal Astronomical Society, 315(3), 534-542, https://academic.oup.com/mnras/article/315/3/534/972014
2000
-
[15]
N., Hobbs, G
Manchester, R. N., Hobbs, G. B., Teoh, A., & Hobbs, M. (2005). The Australia telescope national facility pulsar catalogue. The Astronomical Journal, 129(4), 1993, https://iopscience.iop.org/article/10.1086/428488/meta
2005 doi
-
[16]
Negreiros, R., Bernal, C., Dexheimer, V., & Troconis, O. (2018). Many Aspects of Magnetic Fields in Neutron Stars. Universe, 4, pp. 43, https://www.mdpi.com/2218-1997/4/3/43
2018
-
[17]
Potekhin, A. Y. (2010). The physics of neutron stars. Physics-Uspekhi, 53(12), 1235, https://iopscience.iop.org/article/10.3367/UFNe.0180.201012c.1279/meta
2010
-
[18]
Reisenegger, A. (2003). Origin and evolution of neutron star magnetic fields. arXiv preprint astro-ph/0307133
2003 arXiv
-
[19]
& Cohen, I
Cohen, I., Huang, Y., Chen, J., Benesty, J., Benesty, J., Chen, J., ... & Cohen, I. (2009). Pearson correlation coefficient. Noise reduction in speech processing, 1-4, https://link.springer.com/chapter/10.1007/978-3-642-00296-0_5
2009 doi
-
[20]
L., Teukolsky, S
Shapiro, S. L., Teukolsky, S. A., & Lightman, A. P. (1983). Black Holes, White Dwarfs, and Neutron Stars: The Physics of Compact Objects. Physics Today, 36, pp. 89-90, https://books.google.com/books?hl=en&lr=&id=d1CRQIcP1zoC&oi=fnd&pg=PP2&dq=Shapiro, +S.+L.,+Teukolsky,+S.+A.,+...
1983
-
[21]
O., & Okeke, P
Urama, J. O., & Okeke, P. N. (1999). Vela-size glitch rates in youthful pulsars. Monthly Notices of the Royal Astronomical Society, 310(2), 313-316, https://academic.oup.com/mnras/article/310/2/313/1046930?login=false Van Paradijs, J., Taam, R. E., & Van Den Heuvel, E. P. J. (...
1999
-
[22]
N., Zhang, J., Wu, X
Wang, N., Manchester, R. N., Zhang, J., Wu, X. J., Yusup, A., Lyne, A. G., & Chen, M. Z. (2001). Pulsar Timing at Urumqi Astronomical Observatory: Observing System and Results. Monthly Notices of the Royal Astronomical Society, 328, pp. 855-866, https://academic.oup.com/mnras/...
2001
-
[23]
N., Pace, R
Wang, N., Manchester, R. N., Pace, R. T., Bailes, M., Kaspi, V. M., Stappers, B. W., & Lyne, A. G. (2000). Glitches in southern pulsars. Monthly Notices of the Royal Astronomical Society, 317(4), 843-860. https://academic.oup.com/mnras/article/317/4/843/1037128
2000
-
[24]
Demorest, P., Freire, P., Hessels, J., Margueron, J., Oertel, M., Patruno, A., Possenti, A., Ransom, S., Stairs, I., & Stappers, B. (2014). Probing the Neutron Star Interior and the Equation of State of Cold Dense Matter with the SKA, http://arxiv.org/abs/1501.00042
2014 arXiv
-
[25]
M., Shannon, R., & Xu, R
Ravi, V., Sarkissian, J. M., Shannon, R., & Xu, R. X. (2013). Detection of 107 glitches in 36 southern pulsars. Monthly Notices of the Royal Astronomical Society, 429, pp. 688-724, https://academic.oup.com/mnras/article/429/1/688/1025310
2013
-
[26]
Zwicky, F. (1958). Nuclear Goblins and Flare Stars. Publications of the Astronomical Society of the Pacific, 70, pp. 506, https://www.jstor.org/stable/40676817
1958
-
[27]
Zwicky, F. (1974). Nuclear goblins and Cosmic Gamma Ray Bursts. Astrophysics and Space Science, 28, pp. 111-114, https://link.springer.com/article/10.1007/BF00642241
1974 doi
Reviewed August 12, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.