pith. sign in

arxiv: 2605.21750 · v1 · pith:ARD4P3EYnew · submitted 2026-05-20 · 🌌 astro-ph.HE

A Generalized Template Matching Algorithm for Correcting Jitter Noise in Pulsar Timing

Ross J. Jennings , James M. Cordes , Shami Chatterjee , Maura A. McLaughlin This is my paper

Pith reviewed 2026-05-22 08:36 UTC · model grok-4.3

classification 🌌 astro-ph.HE
keywords pulsar timingjitter noisetemplate matchingprincipal component analysisgravitational wavestime of arrivalpulse shape variation
0
0 comments X

The pith

Generalized template matching with principal components corrects jitter noise in pulsar timing

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

Pulsar timing arrays rely on precise pulse arrival times to detect gravitational waves and other high-energy phenomena. Standard template matching assumes a fixed average pulse shape, yet intrinsic jitter produces short-term changes in amplitude and shape that add noise to the measured times of arrival. The paper generalizes the template matching approach by folding in principal component analysis so that the dominant modes of shape variation can be included in the fit and removed. The authors compare the new procedure against other jitter-mitigation ideas, with explicit checks that the corrections do not absorb the astrophysical signals being sought, and they verify performance on simulated data.

Core claim

By incorporating principal component analysis into template matching, variations in pulse shape can be modeled and subtracted, leading to cleaner time-of-arrival estimates that are less contaminated by jitter noise.

What carries the argument

Principal component analysis applied to pulse profiles to derive basis functions that represent the main modes of shape variation, which are then included in the template matching fit.

If this is right

  • Time-of-arrival measurements become more precise for pulsars that exhibit jitter.
  • Timing residuals exhibit lower noise floors, improving the sensitivity of pulsar timing arrays to gravitational waves.
  • The correction can be adjusted so that it does not remove long-term astrophysical effects.
  • Simulations show advantages over some existing jitter-mitigation proposals.

Where Pith is reading between the lines

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

  • Application to real pulsar-timing-array data sets could increase the significance of any existing gravitational-wave background.
  • Analogous dimensionality-reduction steps might reduce shape-related noise in other periodic signals observed in radio or X-ray astronomy.
  • The same principal-component basis could be updated continuously as new observations arrive, allowing the method to track slow changes in average pulse shape.

Load-bearing premise

The principal components capture the relevant jitter variations without absorbing or removing other astrophysical signals of interest such as gravitational wave signatures.

What would settle it

A simulation in which a known gravitational-wave signal is injected into data containing realistic jitter; after correction the wave signal should be recovered at higher significance than in the uncorrected data.

Figures

Figures reproduced from arXiv: 2605.21750 by James M. Cordes, Maura A. McLaughlin, Ross J. Jennings, Shami Chatterjee.

Figure 1
Figure 1. Figure 1: Average profiles for three millisecond pulsars (PSRs J1909−3744, J2145−0750, and J1012+5307), demonstrating the variety of pulse shapes exhibited by MSPs. All three profiles were taken from the NANOGrav 12.5-year data release (Alam et al. 2021), and are derived from observations made using the 1–2 GHz receiver and GUPPI backend on the Green Bank Telescope between 2010 and 2017. More sophisticated template-… view at source ↗
Figure 2
Figure 2. Figure 2: A comparison of the results of different TOA correction methods on simulated pulse profiles. The results shown here are for simulated profiles with two identical Gaussian components, with the mean separation between the profile components shown on the horizontal axis. The standard deviation of the TOA estimates resulting from each method is shown on the vertical axis. The left panel shows the results for p… view at source ↗
Figure 3
Figure 3. Figure 3: A comparison of TOA correction methods on simulated data, as in [PITH_FULL_IMAGE:figures/full_fig_p011_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: The standard deviation of TOA estimates in simulated profiles with two identical components of width w = 0.0425 and varying spacing. Crosses indicate results of simulated cases, and solid lines indicate predictions based on equations (8) and (B36). Four cases are shown: one in which the pulses have amplitude and phase variations (m = 1.0, f = 0.3 for each component), one in which the pulses have only phase… view at source ↗
Figure 5
Figure 5. Figure 5: The standard deviation of TOA estimates in simulated profiles with several identical components of width w = 0.0425 and varying spacing. Results with amplitude variations only (m = 1.0) is shown in the upper left, results with phase variations only (fj = 0.3) in the upper right, and results with both amplitude and phase variations (m = 1.0, fj = 0.3) in the lower left. As in [PITH_FULL_IMAGE:figures/full_… view at source ↗
read the original abstract

Pulsar timing is a valuable source of high-precision astrophysical measurements which can be used to probe gravitational physics, including by detecting gravitational waves. An important factor limiting the precision of these measurements is pulse jitter; i.e., intrinsic, short-timescale variation in the amplitude and shape of pulses from a given pulsar. Because conventional pulse time-of-arrival (TOA) measurement relies on template matching, which assumes the average pulse shape is stable, such variation gives rise to jitter noise in TOA measurements. Here we introduce a generalization of the template matching technique, making use of principal component analysis, which can account for variations in pulse shape. We compare this technique to other proposals for mitigating jitter noise in pulsar timing, paying particular attention to the possibility of corrections absorbing other astrophysical signals of interest, and demonstrate its effectiveness using simulated data.

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

2 major / 2 minor

Summary. The paper introduces a generalization of template matching for pulsar timing that incorporates principal component analysis (PCA) to account for intrinsic pulse shape variations and thereby reduce jitter noise in time-of-arrival measurements. It compares the approach to existing jitter-mitigation methods while emphasizing the need to avoid absorbing other astrophysical signals such as gravitational waves, and demonstrates performance on simulated data in which jitter is injected independently of the timing model.

Significance. If the PCA-augmented template matching successfully isolates jitter-induced shape changes without projecting out deterministic signals in the timing model, the technique could improve the precision of pulsar timing arrays and enhance their sensitivity to nanohertz gravitational waves. The paper's attention to the absorption risk and its use of simulations are positive features, though translation to real data remains to be shown.

major comments (2)
  1. [Abstract / Method] Abstract and method description: the paper notes the importance of ensuring corrections do not absorb gravitational-wave or other timing-model signals, yet provides no analytic demonstration or explicit cross-validation that the PCA projection operator commutes with the least-squares timing solution (i.e., that retained principal components remain orthogonal to the design matrix containing spin-down, position, and GW basis functions). This is load-bearing for the central claim.
  2. [Simulations] Simulations section: while jitter is injected independently of the timing model, the tests do not appear to include cases in which jitter variations share frequency content with GW signals or in which the PCA basis is estimated from data already containing timing residuals; such checks would be required to confirm the method does not subtract part of the target signal.
minor comments (2)
  1. [Abstract] The abstract states that quantitative comparisons are performed but does not report specific metrics (e.g., reduction in TOA uncertainty or residual RMS); adding these numbers would improve clarity.
  2. [Method] Notation for the PCA components and the projection operator should be defined explicitly early in the text to aid readers unfamiliar with the generalization.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We are grateful to the referee for their constructive feedback on our manuscript. We have carefully considered each major comment and provide point-by-point responses below. Where appropriate, we have outlined revisions to address the concerns raised.

read point-by-point responses

  1. Referee: [Abstract / Method] Abstract and method description: the paper notes the importance of ensuring corrections do not absorb gravitational-wave or other timing-model signals, yet provides no analytic demonstration or explicit cross-validation that the PCA projection operator commutes with the least-squares timing solution (i.e., that retained principal components remain orthogonal to the design matrix containing spin-down, position, and GW basis functions). This is load-bearing for the central claim.

    Authors: We agree that an explicit demonstration of the orthogonality between the PCA components and the timing model design matrix is crucial for validating the method's safety with respect to gravitational wave signals. In the revised manuscript, we will add a new subsection in the Methods section providing an analytic argument based on the independence of pulse shape variations from the timing parameters. We will also include results from cross-validation simulations where no jitter is injected, showing that the timing solution remains unaffected by the PCA correction. This addresses the load-bearing aspect of the claim. revision: yes

  2. Referee: [Simulations] Simulations section: while jitter is injected independently of the timing model, the tests do not appear to include cases in which jitter variations share frequency content with GW signals or in which the PCA basis is estimated from data already containing timing residuals; such checks would be required to confirm the method does not subtract part of the target signal.

    Authors: The referee correctly identifies potential limitations in the current simulation setup. To strengthen the validation, we will expand the Simulations section to incorporate additional test cases. Specifically, we will simulate jitter variations that include frequency components overlapping with those of injected gravitational wave signals, and we will also derive the PCA basis from datasets that include realistic timing residuals. These additions will provide further evidence that the method does not inadvertently subtract astrophysical signals of interest. revision: yes

Circularity Check

0 steps flagged

No circularity detected in the derivation chain

full rationale

The paper presents a generalization of template matching that incorporates principal component analysis to model pulse shape variations for jitter noise mitigation in pulsar timing. It applies standard PCA to pulse profile data and evaluates performance on simulated datasets with independently injected jitter. The abstract and described approach contain no equations, self-definitions, or reductions where a claimed prediction or result is constructed directly from fitted parameters or prior self-citations in a load-bearing way. The method is framed as an extension of conventional techniques with attention to signal absorption risks, but without any quoted steps that equate outputs to inputs by construction. This constitutes a self-contained empirical and algorithmic proposal rather than a circular derivation.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

Based solely on the abstract; the central method rests on the assumption that PCA can model pulse variations without affecting other signals.

axioms (1)
  • domain assumption Pulse shape variations are adequately represented by a linear combination of principal components derived from the data.
    This underpins the generalization of template matching described in the abstract.

pith-pipeline@v0.9.0 · 5685 in / 1127 out tokens · 57746 ms · 2026-05-22T08:36:25.366134+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

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

  1. [1]

    2023, ApJL, 951, L11, doi: 10.3847/2041-8213/acdc91

    Afzal, A., Agazie, G., Anumarlapudi, A., et al. 2023, ApJL, 951, L11, doi: 10.3847/2041-8213/acdc91

    work page doi:10.3847/2041-8213/acdc91 2023

  2. [2]

    doi:10.3847/2041-8213/acdac6 , keywords =

    Agazie, G., Anumarlapudi, A., Archibald, A. M., et al. 2023a, ApJL, 951, L8, doi: 10.3847/2041-8213/acdac6 —. 2023b, ApJL, 952, L37, doi: 10.3847/2041-8213/ace18b

    work page doi:10.3847/2041-8213/acdac6 2041

  3. [3]

    F., Arzoumanian, Z., Baker, P

    Alam, M. F., Arzoumanian, Z., Baker, P. T., et al. 2021, ApJS, 252, 4, doi: 10.3847/1538-4365/abc6a0

    work page doi:10.3847/1538-4365/abc6a0 2021

  4. [4]

    M., & Downs, G

    Cordes, J. M., & Downs, G. S. 1985, ApJS, 59, 343, doi: 10.1086/191076

    work page doi:10.1086/191076 1985

  5. [5]

    NE2001.I. A New Model for the Galactic Distribution of Free Electrons and its Fluctuations

    Cordes, J. M., & Lazio, T. J. W. 2002, arXiv e-prints, astro. https://arxiv.org/abs/astro-ph/0207156

    work page internal anchor Pith review Pith/arXiv arXiv 2002

  6. [6]

    Craft, Jr., H. D. 1970, PhD thesis, Cornell University, New York

    work page 1970

  7. [7]

    T., Fonseca , E., Ransom , S

    Cromartie, H. T., Fonseca, E., Ransom, S. M., et al. 2020, Nature Astronomy, 4, 72, doi: 10.1038/s41550-019-0880-2

    work page doi:10.1038/s41550-019-0880-2 2020

  8. [8]

    Demorest, P. B. 2007, PhD thesis, University of California, Berkeley. https://www.cv.nrao.edu/∼pdemores/thesis.pdf

    work page 2007

  9. [9]

    B., Ferdman, R

    Demorest, P. B., Ferdman, R. D., Gonzalez, M. E., et al. 2013, ApJ, 762, 94, doi: 10.1088/0004-637X/762/2/94

    work page doi:10.1088/0004-637x/762/2/94 2013

  10. [10]

    Desvignes and R

    Desvignes, G., Caballero, R. N., Lentati, L., et al. 2016, MNRAS, 458, 3341, doi: 10.1093/mnras/stw483

    work page doi:10.1093/mnras/stw483 2016

  11. [11]

    S., & Krause-Polstorff, J

    Downs, G. S., & Krause-Polstorff, J. 1986, ApJS, 62, 81, doi: 10.1086/191134

    work page doi:10.1086/191134 1986

  12. [12]

    D., & Craft, Jr., H

    Drake, F. D., & Craft, Jr., H. D. 1968, Science, 160, 758, doi: 10.1126/science.160.3829.758

    work page doi:10.1126/science.160.3829.758 1968

  13. [13]

    T., & Stappers, B

    Edwards, R. T., & Stappers, B. W. 2003, A&A, 407, 273, doi: 10.1051/0004-6361:20030716

    work page doi:10.1051/0004-6361:20030716 2003

  14. [14]

    J., Manchester, R

    Helfand, D. J., Manchester, R. N., & Taylor, J. H. 1975, ApJ, 198, 661, doi: 10.1086/153644

    work page doi:10.1086/153644 1975

  15. [15]

    2013, Classical and Quantum Gravity, 30, 224007, doi: 10.1088/0264-9381/30/22/224007

    Hobbs, G. 2013, Classical and Quantum Gravity, 30, 224007, doi: 10.1088/0264-9381/30/22/224007

    work page doi:10.1088/0264-9381/30/22/224007 2013

  16. [16]

    A., Anderson, S

    Jenet, F. A., Anderson, S. B., Kaspi, V. M., Prince, T. A., & Unwin, S. C. 1998, ApJ, 498, 365, doi: 10.1086/305529

    work page doi:10.1086/305529 1998

  17. [17]

    L., McLaughlin, M

    Jones, M. L., McLaughlin, M. A., Roy, J., et al. 2021, ApJ, 915, 15, doi: 10.3847/1538-4357/abfdc5

    work page doi:10.3847/1538-4357/abfdc5 2021

  18. [18]

    2015, MNRAS, 452, 607, doi: 10.1093/mnras/stv1296

    Kerr, M. 2015, MNRAS, 452, 607, doi: 10.1093/mnras/stv1296

    work page doi:10.1093/mnras/stv1296 2015

  19. [19]

    Lentati, L., Alexander, P., & Hobson, M. P. 2015, MNRAS, 447, 2159, doi: 10.1093/mnras/stu2611

    work page doi:10.1093/mnras/stu2611 2015

  20. [20]

    A., et al., 2011, @doi [ ] 10.1111/j.1365-2966.2011.18706.x , http://adsabs.harvard.edu/abs/2011MNRAS.417.1621D 417, 1621

    Liu, K., Keane, E. F., Lee, K. J., et al. 2012, MNRAS, 420, 361, doi: 10.1111/j.1365-2966.2011.20041.x

    work page doi:10.1111/j.1365-2966.2011.20041.x 2012

  21. [21]

    S., Miles, M

    Nathan, R. S., Miles, M. T., Ashton, G., et al. 2023, MNRAS, 523, 4405, doi: 10.1093/mnras/stad1660 Os lowski, S., van Straten, W., Hobbs, G. B., Bailes, M., &

    work page doi:10.1093/mnras/stad1660 2023

  22. [22]

    A., et al., 2011, @doi [ ] 10.1111/j.1365-2966.2011.18706.x , http://adsabs.harvard.edu/abs/2011MNRAS.417.1621D 417, 1621

    Demorest, P. 2011, MNRAS, 418, 1258, doi: 10.1111/j.1365-2966.2011.19578.x

    work page doi:10.1111/j.1365-2966.2011.19578.x 2011

  23. [23]

    T., Demorest, P

    Pennucci, T. T., Demorest, P. B., & Ransom, S. M. 2014, ApJ, 790, 93, doi: 10.1088/0004-637X/790/2/93 19

    work page doi:10.1088/0004-637x/790/2/93 2014

  24. [24]

    Rankin, J. M. 1986, ApJ, 301, 901, doi: 10.1086/163955

    work page doi:10.1086/163955 1986

  25. [25]

    M., Os lowski, S., Dai, S., et al

    Shannon, R. M., Os lowski, S., Dai, S., et al. 2014, MNRAS, 443, 1463, doi: 10.1093/mnras/stu1213

    work page doi:10.1093/mnras/stu1213 2014

  26. [26]

    2023, MNRAS, 520, 4562, doi: 10.1093/mnras/stad135 Sosa Fiscella, S

    Song, X., Weltevrede, P., Szary, A., et al. 2023, MNRAS, 520, 4562, doi: 10.1093/mnras/stad135 Sosa Fiscella, S. V., Lam, M. T., & McLaughlin, M. A. 2025, ApJ, 984, 111, doi: 10.3847/1538-4357/adc1c2

    work page doi:10.1093/mnras/stad135 2023

  27. [27]

    Stairs, I. H. 2003, Living Reviews in Relativity, 6, 5, doi: 10.12942/lrr-2003-5

    work page doi:10.12942/lrr-2003-5 2003

  28. [28]

    Stinebring, D., & Cordes, J. M. 1981, ApJ, 249, 704, doi: 10.1086/159332

    work page doi:10.1086/159332 1981

  29. [29]

    Taylor, J. H. 1992, Philosophical Transactions of the Royal Society of London Series A, 341, 117, doi: 10.1098/rsta.1992.0088

    work page doi:10.1098/rsta.1992.0088 1992

  30. [30]

    H., & Weisberg, J

    Taylor, J. H., & Weisberg, J. M. 1982, ApJ, 253, 908, doi: 10.1086/159690 van Straten, W. 2006, ApJ, 642, 1004, doi: 10.1086/501001

    work page doi:10.1086/159690 1982

  31. [31]

    1994, Science, 264, 538, doi: 10.1126/science.264.5158.538

    Wolszczan, A. 1994, Science, 264, 538, doi: 10.1126/science.264.5158.538

    work page doi:10.1126/science.264.5158.538 1994