pith. machine review for the scientific record. sign in

arxiv: 2605.14348 · v1 · submitted 2026-05-14 · 🌌 astro-ph.IM · astro-ph.HE

Recognition: no theorem link

A Z₁² framework for rotational-parameter estimation and uncertainty quantification in high-energy pulsars

Akshat Singhal , Rohit Nair , Devendra Sahu , Gayathri Raman , Suman Bala
Authors on Pith no claims yet

Pith reviewed 2026-05-15 02:13 UTC · model grok-4.3

classification 🌌 astro-ph.IM astro-ph.HE
keywords Z_1^2 statisticpulsar spin estimationfrequency derivativeuncertainty quantificationhigh-energy pulsarsPoisson photon dataAstroSat observations
0
0 comments X

The pith

Referencing pulsar frequency to the observation midpoint removes its leading covariance with the frequency derivative in Z_1^2 analyses.

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

The paper develops a Z_1^2-based framework to recover spin frequency and frequency derivative from photon event lists in high-energy pulsars. Its central step is showing that writing the trial frequency at the exact midpoint of the observation interval cancels the dominant cross term between the two parameters in the local response surface. The resulting diagonal quadratic form then supplies closed-form uncertainty estimates that depend only on the fitted peak height and the measured widths along each axis. These expressions are tested across Monte Carlo ensembles that vary span, signal strength, and good-time-interval structure, and they match the observed run-to-run scatter without requiring a fresh simulation for each dataset. The same estimators are applied to AstroSat/LAXPC data on the Crab pulsar and two accreting sources.

Core claim

For sinusoidal signals the Z_1^2 statistic near its maximum is locally quadratic in frequency and frequency derivative. Shifting the frequency reference to the precise midpoint of the observation window eliminates the leading-order cross term, leaving a diagonal quadratic whose curvatures yield simple analytic variances for each parameter in terms of peak amplitude and local widths.

What carries the argument

The local quadratic form of the Z_1^2 response after frequency is expressed at the observation midpoint, which diagonalizes the covariance between frequency and frequency derivative.

If this is right

  • Uncertainty estimates follow directly from fitted peak amplitude and measured widths without exhaustive Monte Carlo runs for each observation.
  • The three tested estimators (segmented regression, coherent derivative scan, and localized two-dimensional fit) all yield consistent results under the midpoint re-referencing.
  • The predicted uncertainties reproduce the observed scatter across varied observing spans, signal strengths, and good-time-interval patterns.
  • Application to real AstroSat data returns stable rotational parameters for the Crab pulsar, Swift J0243.6+6124, and SAX J1808.4-3658.

Where Pith is reading between the lines

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

  • The midpoint re-referencing device may simplify covariance handling in other coherent search statistics that share the same two-parameter space.
  • For sources whose pulse profiles contain strong higher harmonics the local quadratic form would require a fresh derivation from the full multi-harmonic response.
  • The same analytic error propagation could be adapted to continuous-wave gravitational searches where frequency and spin-down parameters appear together.

Load-bearing premise

The incoming signal is purely sinusoidal, so that the Z_1^2 surface remains accurately quadratic near its peak.

What would settle it

Monte Carlo trials with non-sinusoidal pulse profiles (for example, those dominated by the second harmonic) would produce parameter scatter that deviates systematically from the analytic uncertainties predicted by peak amplitude and widths.

Figures

Figures reproduced from arXiv: 2605.14348 by Akshat Singhal, Devendra Sahu, Gayathri Raman, Rohit Nair, Suman Bala.

Figure 1
Figure 1. Figure 1: Schematic simulated photon arrival series for a sinusoidal source with an intentionally large spin-down term. The underlying rate is 𝑦(𝑡) = 100+80 sin[2𝜋(20𝑡 −𝑡 2 ) ], corresponding to 𝑓0 = 20 Hz and ¤𝑓0 = −2 Hz s−1 at the reference epoch 𝑡 = 0. The event times are binned at 0.05 s over an observation time of 𝑇 = 9.5 s. These parameters are unphysical for pulsars and are used only to make the effect of the… view at source ↗
Figure 3
Figure 3. Figure 3: Variation of the 𝑍 2 1 profiles across independent Poisson realisa￾tions of the same faint simulated source. All realisations use identical injected parameters, 𝑎 = 100, 𝑏 = 8, 𝑇 = 500 s, 𝑓0 = 1 Hz, and ¤𝑓0 = −10−6 Hz s−1 . Poisson fluctuations change both the peak amplitude and the recovered peak location. The scatter in peak location illustrates why the macroscopic peak width alone is not a statistical u… view at source ↗
Figure 4
Figure 4. Figure 4: Illustration of segmented frequency regression for a simulated spin-down signal with data gaps. Left: binned event counts over the observation, with the three selected intervals highlighted. Middle: independent 𝑍 2 1 ( 𝑓 ) profiles computed in the selected intervals. Right: interval-wise recovered frequencies, 𝑓 (int) 𝑟 , plotted against interval midpoint time. The weighted linear fit gives the recovered f… view at source ↗
Figure 5
Figure 5. Figure 5: Effect of segment duration on the segmented frequency-regression method. Rows show three choices of 𝑡int: short segments (top), an intermediate choice (middle), and long segments (bottom). Left: segment selection, with the analysed time intervals shaded. For very short segments, only a subset of intervals is shaded for visual clarity. Middle: representative 𝑍 2 1 ( 𝑓 ) profiles from independent Poisson rea… view at source ↗
Figure 6
Figure 6. Figure 6: One-dimensional 𝑍 2 1 ( 𝑓 ) profiles obtained from the same simulated event list for three selected trial values of the frequency derivative, ¤𝑓trial. The simulated source parameters are 𝑎 = 100, 𝑏 = 80, 𝑓0 = 2 Hz, ¤𝑓0 = −3×10−5 Hz s−1 , and 𝑇 = 3000 s. The blue, red, and green curves correspond approximately to ¤𝑓trial = −3.025 × 10−5 , −3.002 × 10−5 , and −2.986 × 10−5 Hz s−1 , respectively. The profile … view at source ↗
Figure 7
Figure 7. Figure 7: Fitted peak properties obtained from the coherent derivative scan. Each point is obtained by fitting the local 𝑍 2 1 ( 𝑓 ) peak at a fixed trial value of ¤𝑓trial. The coloured markers correspond to the same three trial derivatives shown in [PITH_FULL_IMAGE:figures/full_fig_p008_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: Localized two-dimensional coherent fitting of a simulated 𝑍 2 1 ( 𝑓 , ¤𝑓 ) surface using the local model in Eq. (21). The simulation uses 𝑎 = 40, 𝑏 = 30, 𝑓0 = 5 Hz, ¤𝑓0 = −2 × 10−5 Hz s−1 , and 𝑇 = 1500 s. The first panel shows the raw 𝑍 2 1 ( 𝑓 , ¤𝑓 ) surface evaluated on a coarse grid around the injected parameters, with the shaded rectangle marking the region selected for the local fit. The second panel… view at source ↗
Figure 10
Figure 10. Figure 10: Monte Carlo validation of the segment-duration choice in seg￾mented frequency regression. The recovered spin-down scatter, 𝜎( ¤𝑓𝑟 ), is shown as a function of the dimensionless intra-segment phase-drift parame￾ter | ¤𝑓 |𝑡 2 int. Very short segments are photon-noise limited because each local 𝑍 2 1 ( 𝑓 ) peak is weak. Very long segments are affected by intra-segment phase evolution, which violates the loca… view at source ↗
Figure 9
Figure 9. Figure 9: Convergence of the Monte Carlo scatter estimates as a function of the number of independent Poisson realizations included in the ensemble. The example shown uses the coherent derivative scan for a fixed simulated source. Curves show the running standard deviation, the three-sigma-clipped standard deviation, and the median absolute deviation from the injected pa￾rameter value for representative recovered pa… view at source ↗
Figure 11
Figure 11. Figure 11: shows that the uncertainty decreases gradually as ad￾ditional segments are included. The dependence is not as steep as a simple independent-measurement average because the spin-down estimate also depends on the temporal placement of the segments and on the quality of each local frequency recovery. Nevertheless, the trend demonstrates that the method is not dominated by a single interval once several usabl… view at source ↗
Figure 12
Figure 12. Figure 12: Monte Carlo scatter of the recovered rotational parameters for the three estimators. Top: recovery in the start-epoch frequency 𝑓start and frequency derivative ¤𝑓 . Bottom: the same realisations expressed in terms of the midpoint frequency 𝑓mid and ¤𝑓 . The coherent estimators show substantially smaller scatter than the segmented regression baseline. The midpoint-frequency parametrization also reduces the… view at source ↗
Figure 13
Figure 13. Figure 13: Dependence of the coherent uncertainty calibration on the total observing span 𝑇. The panels compare the Monte Carlo scatter of 𝑓mid,𝑟 and ¤𝑓𝑟 with the uncertainty predicted from the fitted local peak amplitude and widths. For fixed source parameters and uninterrupted observations, the expected scalings are approximately 𝜎( 𝑓mid,𝑟 ) ∝ 𝑇 −3/2 and 𝜎( ¤𝑓𝑟 ) ∝ 𝑇 −5/2 . As discussed in Section 4.3, this coordi… view at source ↗
Figure 14
Figure 14. Figure 14: Uncertainty calibration of the two coherent estimators as a function of signal strength. The horizontal axis shows 𝑏 2 /𝑎, which sets the expected coherent peak power for fixed observing span. The vertical axis shows the ratio between the Monte Carlo scatter of the recovered parameter and the fitted uncertainty prediction. Values close to unity indicate that the fitted peak quantities provide a calibrated… view at source ↗
Figure 15
Figure 15. Figure 15: Representative raw 𝑍 2 1 ( 𝑓 , ¤𝑓 ) surface evaluated on an 11 × 11 grid. The white marker denotes the injected parameters and the black marker denotes the recovered parameters after fitting the local coherent response. The figure illustrates the role of the raw grid in Method 3: it is used to identify and sample the coherent ridge, while the final parameter estimate is obtained from the local two-dimensi… view at source ↗
Figure 16
Figure 16. Figure 16: Dependence of the Method 3 uncertainty calibration on the raw two-dimensional grid resolution. The horizontal axis gives the grid size 𝑛 for an 𝑛 × 𝑛 search in ( 𝑓mid, ¤𝑓 ). The vertical axis shows the ratio of the Monte Carlo scatter to the fitted uncertainty prediction. The recovery remains stable once the local coherent peak is adequately sampled. Nevertheless, the fitted local response absorbs much of… view at source ↗
Figure 18
Figure 18. Figure 18: Robustness of the coherent uncertainty estimates to observational gaps. The horizontal axis shows the fraction of removed observing time. For each gap fraction, several random gap patterns are tested. The vertical axis shows the ratio between the Monte Carlo scatter of the recovered parameter and the median fitted uncertainty prediction from the individual realisations. The dotted horizontal line marks pe… view at source ↗
Figure 19
Figure 19. Figure 19: Segmented frequency recovery for the three AstroSat event lists using Method I. Left: Crab pulsar, analysed using six independent intervals of duration 20000 s from a longer event list of duration ∼ 300000 s. Middle: SAX J1808.4−3658, analysed using intervals of duration 10000 s from an event list of duration ∼ 100000 s. Right: Swift J0243.6+6124, analysed using intervals of duration 1000 s from an event … view at source ↗
Figure 20
Figure 20. Figure 20: Coherent derivative scan for the Crab pulsar, Swift J0243.6+6124, and SAX J1808.4−3658 using Method II. For each trial frequency derivative ¤𝑓trial, a local one-dimensional 𝑍 2 1 ( 𝑓 ) peak is fitted and its recovered amplitude 𝐴𝑟 is recorded. The resulting amplitude envelope is then fitted with a local sinc2 response to estimate the preferred ¤𝑓 . The fitted peak location gives the recovered frequency de… view at source ↗
Figure 21
Figure 21. Figure 21: Two-dimensional coherent recovery for the Crab pulsar, Swift J0243.6+6124, and SAX J1808.4−3658 using Method III. Each row corresponds to one source. The left panels show the raw two-dimensional 𝑍 2 1 surface in the local ( 𝑓 , ¤𝑓 ) search region. The middle panels show the local region used for fitting the central coherent response. The right panels show the corresponding fitted two-dimensional response … view at source ↗
Figure 24
Figure 24. Figure 24: Convergence of the recovered-parameter scatter for the simulated Crab-like signal shown in [PITH_FULL_IMAGE:figures/full_fig_p016_24.png] view at source ↗
Figure 25
Figure 25. Figure 25: Illustration of the midpoint-frequency reparametrization for a single simulated pulsar event list. The events are generated using a sinusoidal count-rate model with mean rate 𝑎 = 80 counts s−1 , pulse amplitude 𝑏 = 40 counts s−1 , start-epoch frequency 𝑓start,0 = 350 Hz, frequency derivative ¤𝑓0 = −2 × 10−5 Hz s−1 , and exposure time 𝑇 = 120 s. Top: the 𝑍 2 1 search surface in the original ( 𝑓start, ¤𝑓 ) … view at source ↗
Figure 26
Figure 26. Figure 26: Source-matched simulated event lists used for additional validation of the timing-uncertainty estimates. The left column shows the input simulated light curves, with the corresponding binned simulated event data overlaid, for the Crab pulsar, SAX J1808.4−3658, and Swift J0243.6+6124 from top to bottom. The right column compares the binned AstroSat/LAXPC event lists with the corresponding simplified simula… view at source ↗
read the original abstract

We present a $Z_1^2$-based framework for estimating the spin frequency and frequency derivative of high-energy pulsars from Poisson-limited photon event lists. The key point is that the width of a coherent detection peak is not, by itself, the statistical uncertainty on the recovered rotational parameters. We develop and compare three computationally efficient estimators: segmented frequency regression, a coherent derivative scan, and a localized two-dimensional coherent fit. For sinusoidal signals, we derive the local form of the Z-squared response as a function of frequency and frequency derivative, and show that expressing the frequency at the midpoint of the observation removes the leading-order covariance between the two parameters. This gives simple uncertainty estimates in terms of the fitted peak amplitude and local widths, without requiring an exhaustive Monte Carlo simulation for each observation. We test these estimates with Monte Carlo simulations over a range of observing spans, signal strengths, grid resolutions, and good-time-interval structures, and show that the predicted uncertainties reproduce the run-to-run scatter of the recovered parameters in the tested regimes. We then apply the framework to AstroSat/LAXPC event lists for the Crab pulsar, Swift J0243.6+6124, and SAX J1808.4-3658. The results provide a practical and statistically motivated route to rotational-parameter estimation for targeted high-energy pulsar searches.

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 manuscript presents a Z_1^2-based framework for estimating spin frequency and frequency derivative from Poisson-limited photon event lists in high-energy pulsars. It compares three estimators (segmented frequency regression, coherent derivative scan, and localized 2D coherent fit), derives the local quadratic form of the Z_1^2 response for sinusoidal signals, demonstrates that evaluating frequency at the observation midpoint removes leading-order covariance between frequency and derivative, and obtains analytic uncertainty estimates from fitted peak amplitude and local widths. These estimates are validated against Monte Carlo simulations across observing spans, signal strengths, grid resolutions, and GTI structures, showing that predicted uncertainties reproduce run-to-run parameter scatter. The framework is applied to AstroSat/LAXPC data for the Crab pulsar, Swift J0243.6+6124, and SAX J1808.4-3658.

Significance. If the analytic uncertainties remain accurate, the work supplies a computationally efficient route to rotational-parameter estimation and uncertainty quantification for targeted high-energy pulsar searches, avoiding exhaustive per-observation Monte Carlo runs while remaining grounded in the shape of the Z_1^2 statistic itself. The explicit validation against independent simulations and the practical demonstration on real event lists strengthen its utility for missions producing Poisson-limited timing data.

major comments (2)
  1. [Derivation of local Z_1^2 response and Monte Carlo validation section] The derivation of the local quadratic Z_1^2 form and the midpoint covariance removal (described in the abstract and the section on the coherent 2D fit) is stated to hold for sinusoidal signals. The Monte Carlo validation reproduces scatter in the tested regimes, but the manuscript does not report whether the simulated signals included higher harmonics; real targets such as the Crab and SAX J1808.4-3658 exhibit non-sinusoidal profiles, so cross-harmonic terms could modify the curvature matrix and residual covariance, undermining direct application of the analytic widths.
  2. [Applications to real data] In the applications to Crab, Swift J0243.6+6124, and SAX J1808.4-3658, the reported uncertainties rely on the sinusoidal-derived formulas. No explicit test is shown confirming that the local quadratic approximation and midpoint decorrelation still reproduce the observed scatter when the Z_1^2 statistic sums multiple harmonics, which is load-bearing for the claim that the framework provides statistically motivated uncertainties for these sources.
minor comments (2)
  1. [Abstract and methods overview] The abstract lists three estimators but the main text should more clearly indicate which estimator supplies the primary uncertainty formulas used in the applications.
  2. [Uncertainty estimation subsection] Notation for the local widths and peak amplitude should be defined once with consistent symbols before the uncertainty expressions are introduced.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their thorough and constructive review. We appreciate the recognition of the framework's potential for efficient uncertainty quantification in high-energy pulsar timing. We address the major comments below and outline the revisions we will make to the manuscript.

read point-by-point responses

  1. Referee: [Derivation of local Z_1^2 response and Monte Carlo validation section] The derivation of the local quadratic Z_1^2 form and the midpoint covariance removal (described in the abstract and the section on the coherent 2D fit) is stated to hold for sinusoidal signals. The Monte Carlo validation reproduces scatter in the tested regimes, but the manuscript does not report whether the simulated signals included higher harmonics; real targets such as the Crab and SAX J1808.4-3658 exhibit non-sinusoidal profiles, so cross-harmonic terms could modify the curvature matrix and residual covariance, undermining direct application of the analytic widths.

    Authors: We thank the referee for highlighting this point. The Z_1^2 statistic is defined exclusively for the fundamental harmonic and does not incorporate or sum higher harmonics; therefore, cross-harmonic terms do not arise within the statistic itself. The local quadratic form of the Z_1^2 surface follows directly from the linear phase model applied to the fundamental Fourier components (sums of cos(2π f t_i) and sin(2π f t_i)), whose curvature in (f, ḟ) space is determined solely by the observation span and photon count, independent of pulse shape. Higher harmonics appear at integer multiples of the trial frequency and do not contribute to the Z_1^2 value or its local curvature at the fundamental. The fitted amplitude A is measured directly from the observed peak height, automatically reflecting the strength of the fundamental component present in any given profile. The Monte Carlo simulations used pure sinusoids to isolate this fundamental response, but the analytic expressions are expected to remain valid for non-sinusoidal signals. In the revised manuscript we will add Monte Carlo experiments using simulated non-sinusoidal profiles (with significant power in the second and third harmonics) to explicitly confirm that the predicted uncertainties continue to match the observed parameter scatter. revision: yes

  2. Referee: [Applications to real data] In the applications to Crab, Swift J0243.6+6124, and SAX J1808.4-3658, the reported uncertainties rely on the sinusoidal-derived formulas. No explicit test is shown confirming that the local quadratic approximation and midpoint decorrelation still reproduce the observed scatter when the Z_1^2 statistic sums multiple harmonics, which is load-bearing for the claim that the framework provides statistically motivated uncertainties for these sources.

    Authors: We agree that an explicit check on real data strengthens the claim. As noted above, the Z_1^2 statistic does not sum multiple harmonics, so the concern about cross-harmonic terms within the statistic does not apply. In the revised manuscript we will augment the applications section with a direct validation: we will subdivide each real event list (Crab, Swift J0243.6+6124, SAX J1808.4-3658) into independent segments of comparable length, recompute the rotational parameters on each segment, and compare the observed run-to-run scatter against the analytic uncertainties derived from the full-observation peak amplitude and widths. We will also report the results of a bootstrap resampling exercise on the full lists. These additions will demonstrate that the midpoint-decorrelated analytic formulas reproduce the empirical scatter for the actual (non-sinusoidal) profiles of these sources. revision: yes

Circularity Check

0 steps flagged

No circularity: analytic derivation of local Z_1^2 quadratic form and midpoint covariance removal stands on the statistic definition and is validated by independent Monte Carlo

full rationale

The paper derives the local quadratic form of the Z_1^2 response surface for sinusoidal signals directly from the definition of the statistic, then shows that shifting the frequency reference to the observation midpoint removes leading-order covariance between frequency and frequency derivative. The resulting analytic uncertainty expressions (in terms of peak amplitude and local widths) are presented as consequences of that derivation. These expressions are subsequently tested against fresh Monte Carlo realizations spanning signal strengths, spans, grids, and GTI structures; the simulations are not used to fit or tune the formulas. No load-bearing step reduces to a fitted input renamed as prediction, no self-citation chain is invoked for uniqueness, and the central claim does not equate to its inputs by construction. The framework is therefore self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The framework rests on the standard assumption that photon arrivals are Poisson and that the signal is sinusoidal; no new free parameters are introduced beyond the usual grid resolution and signal amplitude, and no new entities are postulated.

axioms (2)
  • domain assumption Photon arrival times follow a Poisson process with constant background rate.
    Invoked throughout the derivation of the Z_1^2 response and the Monte Carlo tests.
  • domain assumption The pulsed signal is purely sinusoidal at the fundamental frequency.
    Required for the local quadratic form of the Z-squared surface derived in the paper.

pith-pipeline@v0.9.0 · 5559 in / 1464 out tokens · 33007 ms · 2026-05-15T02:13:26.426017+00:00 · methodology

discussion (0)

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

Reference graph

Works this paper leans on

89 extracted references · 89 canonical work pages · 5 internal anchors

  1. [1]

    Nature , volume=

    A millisecond pulsar in an X-ray binary system , author=. Nature , volume=. 1998 , publisher=

    work page 1998

  2. [2]

    Nature , volume=

    The two-hour orbit of a binary millisecond X-ray pulsar , author=. Nature , volume=. 1998 , publisher=

    work page 1998

  3. [3]

    arXiv e-prints, art , author=

    Accretion powered X-ray millisecond pulsars. arXiv e-prints, art , author=. arXiv preprint arXiv:2010.09005 , year=

    work page arXiv 2010

  4. [4]

    The Astrophysical Journal Letters , volume=

    NICER Discovery of Millisecond X-ray Pulsations and an Ultracompact Orbit in IGR J17494-3030 , author=. The Astrophysical Journal Letters , volume=. 2021 , publisher=

    work page 2021

  5. [5]

    The Astrophysical Journal Letters , volume=

    The discovery of the 528.6 Hz accreting millisecond X-ray pulsar MAXI J1816--195 , author=. The Astrophysical Journal Letters , volume=. 2022 , publisher=

    work page 2022

  6. [6]

    Monthly Notices of the Royal Astronomical Society: Letters , volume=

    MAXI J1957+ 032: a new accreting millisecond X-ray pulsar in an ultra-compact binary , author=. Monthly Notices of the Royal Astronomical Society: Letters , volume=. 2022 , publisher=

    work page 2022

  7. [7]

    4--3658 during its 2019 outburst , author=

    AstroSat observation of the accreting millisecond X-ray pulsar SAX J1808. 4--3658 during its 2019 outburst , author=. Monthly Notices of the Royal Astronomical Society , volume=. 2023 , publisher=

    work page 2019

  8. [8]

    arXiv e‑prints , eprint =

    Hu, Chin‑Ping and Chou, Yi and Chung, Yi‑Ying , title =. arXiv e‑prints , eprint =. 2007 , archivePrefix=

    work page 2007

  9. [9]

    arXiv e‑prints , eprint =

    Chou, Yi , title =. arXiv e‑prints , eprint =. 2014 , archivePrefix=

    work page 2014

  10. [10]

    Annual Review of Astronomy and Astrophysics , volume=

    Ultraluminous X-ray sources , author=. Annual Review of Astronomy and Astrophysics , volume=. 2017 , publisher=

    work page 2017

  11. [11]

    6+ 6124 , author=

    AstroSat observations of the first Galactic ULX pulsar Swift J0243. 6+ 6124 , author=. Monthly Notices of the Royal Astronomical Society , volume=. 2021 , publisher=

    work page 2021

  12. [12]

    6+ 6124: Swift discovery of an accreting NS transient , author=

    Swift J0243. 6+ 6124: Swift discovery of an accreting NS transient , author=. The Astronomer's Telegram , volume=

    work page

  13. [13]

    , author=

    GRB 171003A: Swift detection of a burst or a Galactic Transient. , author=. GRB Coordinates Network , volume=

    work page

  14. [14]

    6+ 6124 , author=

    Fermi GBM detects pulsations from Swift J0243. 6+ 6124 , author=. The Astronomer's Telegram , volume=

    work page

  15. [15]

    , author=

    GRB 171004A: Swift-XRT refined Analysis. , author=. GRB Coordinates Network , volume=

    work page

  16. [16]

    6+ 6124 , author=

    Swift and NuSTAR observations of Swift J0243. 6+ 6124 , author=. The Astronomer's Telegram , volume=

    work page

  17. [17]

    6+ 6124 , author=

    Understanding the spectral and timing behaviour of a newly discovered transient X-ray pulsar Swift J0243. 6+ 6124 , author=. Monthly Notices of the Royal Astronomical Society , volume=. 2018 , publisher=

    work page 2018

  18. [18]

    and Harrison, F

    Bachetti, M. and Harrison, F. A. and Walton, D. J. and others , title =. Nature , year =

    work page

  19. [19]

    Lyne, A. G. and Pritchard, R. S. and Smith, F. G. , title =. Monthly Notices of the Royal Astronomical Society , volume =. 1993 , doi =

    work page 1993

  20. [20]

    Weisskopf, M. C. and O'Dell, S. L. and Paerels, F. and Elsner, R. F. and Becker, W. and Tennant, A. F. and Swartz, D. A. , title =. The Astrophysical Journal , volume =. 2004 , doi =

    work page 2004

  21. [21]

    , title =

    van Straten, W. , title =. arXiv e‑prints , eprint =

    work page

  22. [22]

    and Verbiest, J

    Liu, K. and Verbiest, J. P. W. and Kramer, M. and Stappers, B. W. and van Straten, W. and Cordes, J. M. , title =. arXiv e‑prints , eprint =

    work page

  23. [23]

    Hellings, R. W. and Downs, G. S. , title =. The Astrophysical Journal , volume =

    work page

  24. [24]

    and Lu, F.-J

    Ge, M.-Y. and Lu, F.-J. and Qu, J.-L. and Zheng, S.-J. and Chen, Y. and Han, D.-W. , title =. arXiv e‑prints , eprint =

    work page

  25. [25]

    Lutovinov, A. A. and Tsygankov, S. S. and Suleimanov, V. F. and Mushtukov, A. A. and Doroshenko, V. and Nagirner, D. I. and Poutanen, J. , title =. Monthly Notices of the Royal Astronomical Society , volume =

    work page

  26. [26]

    arXiv preprint arXiv:2209.07954 , year=

    Fourier Domain , author=. arXiv preprint arXiv:2209.07954 , year=

    work page arXiv

  27. [27]

    Monthly Notices of the Royal Astronomical Society , volume=

    An improved test for periodicity , author=. Monthly Notices of the Royal Astronomical Society , volume=

    work page

  28. [28]

    Astronomy and Astrophysics , volume=

    Search for pulsed gamma-ray emission from radio pulsars in the COS-B data , author=. Astronomy and Astrophysics , volume=

    work page

  29. [29]

    The Astronomical Journal , volume=

    Fourier techniques for very long astrophysical time-series analysis , author=. The Astronomical Journal , volume=. 2002 , publisher=

    work page 2002

  30. [30]

    Verbiest, J. P. W. and Bailes, M. and van Straten, W. and Hobbs, G. B. and others , title =. Monthly Notices of the Royal Astronomical Society , volume =. 2009 , doi =

    work page 2009

  31. [31]

    The Astrophysical Journal , volume=

    Precise timing and phase-resolved spectroscopy of the young pulsar J1617--5055 with NuSTAR , author=. The Astrophysical Journal , volume=. 2021 , publisher=

    work page 2021

  32. [32]

    The Astrophysical Journal , volume=

    Stingray: a modern Python library for spectral timing , author=. The Astrophysical Journal , volume=. 2019 , publisher=

    work page 2019

  33. [33]

    Space Telescopes and Instrumentation 2016: Ultraviolet to Gamma Ray , year = 2016, editor =

    Space Telescopes and Instrumentation 2016: Ultraviolet to Gamma Ray. Space Telescopes and Instrumentation 2016: Ultraviolet to Gamma Ray , year = 2016, editor =

    work page 2016

  34. [34]

    , year = 2017, month = jul, volume =

    Calibration of the Large Area X-Ray Proportional Counter (LAXPC) Instrument on board AstroSat. , year = 2017, month = jul, volume =. doi:10.3847/1538-4365/aa7a0e , archivePrefix =

    work page doi:10.3847/1538-4365/aa7a0e 2017

  35. [35]

    , year = 2019, month = feb, volume =

    Thermonuclear X-ray bursts in rapid succession in 4U 1636-536 with AstroSat-LAXPC. , year = 2019, month = feb, volume =. doi:10.1093/mnras/sty2975 , archivePrefix =

    work page doi:10.1093/mnras/sty2975 2019

  36. [36]

    , year = 2020, month = jul, volume =

    Timing the Pulsations of the Accreting Millisecond Pulsar SAX J1808.4-3658 during Its 2019 Outburst. , year = 2020, month = jul, volume =. doi:10.3847/1538-4357/ab9827 , archivePrefix =

    work page doi:10.3847/1538-4357/ab9827 2019

  37. [37]

    Space Telescopes and Instrumentation 2014: Ultraviolet to Gamma Ray , year = 2014, editor =

    ASTROSAT mission. Space Telescopes and Instrumentation 2014: Ultraviolet to Gamma Ray , year = 2014, editor =

    work page 2014

  38. [38]

    , year = 2010, month = jan, volume =

    Absolute Timing of the Crab Pulsar with the INTEGRAL/SPI Telescope. , year = 2010, month = jan, volume =. doi:10.1088/0004-637X/708/1/403 , archivePrefix =

    work page doi:10.1088/0004-637x/708/1/403 2010

  39. [39]

    Science , year = 2010, month = jun, volume =

    Switched Magnetospheric Regulation of Pulsar Spin-Down. Science , year = 2010, month = jun, volume =

    work page 2010

  40. [40]

    , year = 1983, month = sep, volume =

    On searches for periodic pulsed emission - The Rayleigh test compared to epoch folding. , year = 1983, month = sep, volume =

    work page 1983

  41. [41]

    , year = 1975, month = jun, volume =

    The period derivatives of pulsars. , year = 1975, month = jun, volume =

    work page 1975

  42. [42]

    , year = 2025, month = feb, volume =

    A Method for Estimating X-Ray Pulsar Period and Pulse Time Delay: Applying the Improved <inline-formula> </inline-formula>-test to Complex Profiles. , year = 2025, month = feb, volume =

    work page 2025

  43. [43]

    , year = 2016, month = may, volume =

    On More Sensitive Periodogram Statistics. , year = 2016, month = may, volume =. doi:10.3847/0004-637X/822/1/14 , archivePrefix =

    work page doi:10.3847/0004-637x/822/1/14 2016

  44. [44]

    Frontiers of Information Technology and Electronic Engineering , year = 2015, volume =

    Shen, Li-Rong and Li, Xiao-Ping and Sun, Hai-Feng and Fang, Hai-Yan and Xue, Meng-Fan , title = ". Frontiers of Information Technology and Electronic Engineering , year = 2015, volume =

    work page 2015

  45. [45]

    , year = 1976, month = feb, volume =

    Least-Squares Frequency Analysis of Unequally Spaced Data. , year = 1976, month = feb, volume =

    work page 1976

  46. [46]

    Studies in astronomical time series analysis. II. Statistical aspects of spectral analysis of unevenly spaced data. , year = 1982, month = dec, volume =

    work page 1982

  47. [47]

    , year = 1989, month = mar, volume =

    Fast Algorithm for Spectral Analysis of Unevenly Sampled Data. , year = 1989, month = mar, volume =

    work page 1989

  48. [48]

    Six signals in the radial velocity data for GJ876

    Improved signal detection algorithms for unevenly sampled data. Six signals in the radial velocity data for GJ876. , year = 2014, month = jul, volume =. doi:10.1093/mnras/stu683 , archivePrefix =

    work page doi:10.1093/mnras/stu683 2014

  49. [49]

    Espinoza and L

    D.Antonopoulou and C.M. Espinoza and L. kuiper and N. Andersson , title = ". , year = 2017, month = sep, doi =

    work page 2017

  50. [50]

    Advances in Space Research , year = 2006, month = jan, volume =

    A broad spectral band Indian Astronomy satellite Astrosat. Advances in Space Research , year = 2006, month = jan, volume =

    work page 2006

  51. [51]

    and Bell, S

    Hewish, A. and Bell, S. J. and Pilkington, J. D. H. and Scott, P. F. and Collins, R. A. , title =. Nature , year =

    work page

  52. [52]

    Lattimer, J. M. and Prakash, M. , title =. Science , year =

    work page

  53. [53]

    and Stairs, I

    Kramer, M. and Stairs, I. H. and Manchester, R. N. and others , title =. Science , year =

    work page

  54. [54]

    Manchester, R. N. and Taylor, J. H. , title =. 1977 , address =

    work page 1977

  55. [55]

    Lorimer, D. R. and Kramer, M. , title =. 2005 , address =

    work page 2005

  56. [56]

    and others , title =

    Agazie, G. and others , title =. The Astrophysical Journal Letters , year =

    work page

  57. [57]

    and others , title =

    Antoniadis, J. and others , title =. Astronomy & Astrophysics , year =

    work page

  58. [58]

    Taylor, J. H. , title =. Philosophical Transactions of the Royal Society A , year =

    work page

  59. [59]

    Hobbs, G. B. and Edwards, R. T. and Manchester, R. N. , title =. Monthly Notices of the Royal Astronomical Society , year =

    work page

  60. [60]

    and others , title =

    Luo, J. and others , title =. The Astrophysical Journal , year =

    work page

  61. [61]

    Precise Gamma-Ray Timing and Radio Observations of 17 Fermi Gamma-Ray Pulsars

    Precise -ray Timing and Radio Observations of 17 Fermi -ray Pulsars. , keywords =. doi:10.1088/0067-0049/194/2/17 , archivePrefix =. 1011.2468 , primaryClass =

    work page internal anchor Pith review Pith/arXiv arXiv doi:10.1088/0067-0049/194/2/17

  62. [62]

    , title =

    Kerr, M. , title =. The Astrophysical Journal , year =

    work page

  63. [63]

    Abdo, A. A. and others , title =. The Astrophysical Journal Supplement Series , year =

    work page

  64. [64]

    Smith, D. A. and others , title =. The Astrophysical Journal , year =

    work page

  65. [65]

    , year = 2009, volume = 697, pages =

    The Large Area Telescope on the Fermi Gamma-ray Space Telescope Mission. , year = 2009, volume = 697, pages =

    work page 2009

  66. [66]

    Space Telescopes and Instrumentation 2016: Ultraviolet to Gamma Ray , series =

    The Neutron star Interior Composition Explorer (NICER): design and development. Space Telescopes and Instrumentation 2016: Ultraviolet to Gamma Ray , series =

    work page 2016

  67. [67]

    , year = 1978, month = sep, volume =

    Period determination using phase dispersion minimization. , year = 1978, month = sep, volume =

    work page 1978

  68. [68]

    On the problem of random vibrations, and of random flights in one, two, or three dimensions , journal =

    XXXI. On the problem of random vibrations, and of random flights in one, two, or three dimensions , journal =. 1919 , volume =

    work page 1919

  69. [69]

    , year = 2018, month = may, volume =

    Orbit and intrinsic spin-up of the newly discovered transient X-ray pulsar Swift J0243.6+6124. , year = 2018, month = may, volume =. doi:10.1051/0004-6361/201732208 , archivePrefix =

    work page doi:10.1051/0004-6361/201732208 2018

  70. [70]

    , year = 2012, month = feb, volume =

    Accelerated Orbital Expansion and Secular Spin-down of the Accreting Millisecond Pulsar SAX J1808.4-3658. , year = 2012, month = feb, volume =. doi:10.1088/2041-8205/746/2/L27 , archivePrefix =

    work page doi:10.1088/2041-8205/746/2/l27 2012

  71. [71]

    AAS/High Energy Astrophysics Division \#10 , year = 2008, series =

    The Long-term Timing of Accreting Millisecond Pulsar SAX J1808.4-3658. AAS/High Energy Astrophysics Division \#10 , year = 2008, series =

    work page 2008

  72. [72]

    , year = 1987, month = jun, volume =

    Searches for pulsed emission with application to four globular cluster X-ray sources. , year = 1987, month = jun, volume =

    work page 1987

  73. [73]

    , year = 2005, month = jan, volume =

    X-ray timing of the accretion-powered millisecond pulsar XTE J1814-338. , year = 2005, month = jan, volume =

    work page 2005

  74. [74]

    , year = 2017, month = oct, volume =

    A timing and spectral study of the transient X-ray pulsar. , year = 2017, month = oct, volume =

    work page 2017

  75. [75]

    Fourier Analysis of Time Series: An Introduction

    work page

  76. [76]

    , year = 1981, month = aug, volume =

    Fast period-searching in unevenly sampled observational data. , year = 1981, month = aug, volume =

    work page 1981

  77. [77]

    , year = 1996, month = jun, volume =

    Parameter estimation in periodograms. , year = 1996, month = jun, volume =

    work page 1996

  78. [78]

    The timing model and precision estimates

    TEMPO2, a new pulsar-timing package - II. The timing model and precision estimates. , year = 2006, month = nov, volume =

    work page 2006

  79. [79]

    Bayesian Spectrum Analysis and Parameter Estimation

    work page

  80. [80]

    The Young Pulsar J1357-6429 and Its Pulsar Wind Nebula

    X-Ray Observations of the Young Pulsar J1357 6429 and Its Pulsar Wind Nebula. , keywords =. doi:10.1088/0004-637X/744/2/81 , archivePrefix =. 1107.1819 , primaryClass =

    work page internal anchor Pith review Pith/arXiv arXiv doi:10.1088/0004-637x/744/2/81

Showing first 80 references.