Pith. sign in

REVIEW 3 major objections 8 minor 1 cited by

The FAST Discovery of a Millisecond Pulsar M15O (PSR J2129+1210O) Hidden in the Harmonics of M15A (PSR J2129+1210A)

T0 review · 3 major / 8 minor · reviewed 2026-08-16 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read M15O is a real millisecond pulsar hidden in the 10th harmonic of M15A.

desk verdict New MSP in M15 is likely real, but the timing solution leans on a huge F2 and the paper doesn't demonstrate the alias search. read the letter →

arxiv 2504.16872 v2 pith:ZTZMIORN submitted 2025-04-23 astro-ph.HE

Yinfeng Dai , Zhichen Pan , Lei Qian , Liyun Zhang , Dejiang Yin , Baoda Li , Yaowei Li , Yuxiao Wu
show 1 more author
Yujie Lian
This is my paper · ORCID
classification astro-ph.HE
keywords millisecondpulsarglobularclusterM15M15OharmonicinterferencestacksearchtimingspinfrequencyderivativeFAST
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

The paper reports the discovery of M15O, a previously missed millisecond pulsar in the globular cluster M15. Its spin period, about 11.06686 ms, is so close to the tenth harmonic of the bright pulsar M15A that standard candidate screening mistook the two signals for one. The authors used a stack of 19 archival observations from the FAST telescope plus a phase-connected timing solution to show that M15O is a distinct, isolated pulsar with its own position, spin-down, and dispersion measure. If correct, the result means real pulsars can hide inside the harmonic structure of bright cluster pulsars, and that adding power spectra to candidate diagnostic plots would catch such cases.

What carries the argument

The load-bearing coincidence is the near-degeneracy between M15O's spin period ($\approx 11.06686$ ms) and the tenth harmonic of M15A ($\approx 11.06647$ ms), a difference the paper reports as roughly 0.004 ms. The mechanism that resolves the degeneracy is threefold: an incoherent stack of power spectra over many epochs raises M15O's faint peak above noise; a frequency-domain chop of M15A's harmonic removes the contaminating signal; and a phase-connected timing fit that includes $\ddot f$ fixes the integer cycle count, giving M15O a position distinct from M15A.

What would settle it

Fold a future FAST observation with the published ephemeris: if the predicted pulse phase drifts by much more than the $108\ \mu$s residual RMS, or if a re-fit of the existing TOAs finds another integer-cycle alias with comparable quality, then the claimed spin-down, position, and separation from M15A collapse.

Watch

Extended reading notes

Core claim

M15O (PSR J2129+1210O) is a genuine isolated millisecond pulsar in M15 with spin frequency $f = 90.3597874705(1)$ Hz, first derivative $\dot f = 1.79191(5)\times 10^{-14}$ Hz s$^{-1}$, second derivative $\ddot f = 3.3133(6)\times 10^{-23}$ Hz s$^{-2}$, and dispersion measure DM $= 67.44(1)$ pc cm$^{-3}$. Its projected position is $0.37''$ from the cluster center and $0.81''$ from M15A, making it the closest known pulsar to the M15 core. It was missed earlier because its period is nearly identical to the tenth harmonic of M15A, a difference the paper reports as roughly 0.004 ms. Timing of 65 arrival times with a phase-connected solution, including the second frequency derivative, is what finally separates it from that harmonic.

Load-bearing premise

The phase-connected timing solution is unique: the 65 arrival times have one correct set of integer pulse counts, and the fitted second frequency derivative is real rather than a stand-in for unmodeled orbital motion or phase wraps.

Editorial extensions

If this is right

  • M15O must be added to the census of M15 pulsars, and its position $0.37''$ from the cluster center means it is the closest known pulsar to the core.
  • Earlier M15 data can be re-examined: the paper notes M15O could already have been detected by January 2021 if the spectra had been inspected.
  • Stacking power spectra before folding is an effective way to uncover pulsars hidden under bright harmonics, so the same approach on other globular clusters may yield more pulsars.
  • Because M15O is the only M15 pulsar with a measured $\ddot f$, and its $\dot f$ is second only to M15D, the timing solution points to strong acceleration near the core; continued monitoring may reveal whether it is gravitational acceleration or orbital motion.
  • Adding the power spectrum to pulsar candidate diagnostic plots is a concrete change to search practice that can prevent similar misses.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • If harmonic proximity is not rare, some single-epoch 'harmonic' detections in other dense clusters may actually be unresolved pulsars; re-stacking archival spectra at matched dispersion measures could test this without new observations.
  • The large $\dot f$ and measured $\ddot f$ are also consistent with an unseen companion in a wide orbit; if future timing shows sinusoidal residuals on a decade timescale, M15O would become a probe of a massive central object in M15.
  • A practical extension would be to quantify the minimum period separation at which two same-DM pulsars can be disentangled by timing; that threshold would tell surveys how close to known harmonics they must search.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 8 minor

Summary. This manuscript reports the discovery of a new millisecond pulsar, M15O (PSR J2129+1210O), in the globular cluster M15, based on FAST L-band observations from 2019 to 2024. The pulsar has a spin period of about 11.06686 ms and a dispersion measure of 67.44(1) pc cm^-3, with a claimed phase-connected timing solution over 65 TOAs giving f = 90.3597874705(1) Hz, fdot = 1.79191(5) x 10^-14 Hz/s, fddot = 3.3133(6) x 10^-23 Hz/s^2, and a timing position separated by 0.81 arcsec from M15A and 0.37 arcsec from the cluster center. The authors argue that M15O was previously missed because its period is extremely close to the 10th harmonic of the bright pulsar M15A, and that spectral stacking plus careful visual inspection of power spectra can reveal such hidden signals. They also report flux-density estimates and discuss the large spin-frequency derivatives of M15O in the context of the dense cluster core.

Significance. If the timing solution is correct, M15O is a genuine addition to the M15 pulsar census and a useful demonstration that bright-pulsar harmonics can hide physically distinct MSP signals in globular cluster searches. The paper's strengths are that the candidate was confirmed by a phase-connected timing solution, that the comparison M15A parameters come from an independent published ephemeris, and that the frequency-domain subtraction of M15A provides a direct visual confirmation of the new signal. The measured second frequency derivative and the very small projected distance from the cluster center make the source potentially interesting for dynamical studies of the M15 core. However, the significance of the positional and dynamical claims rests entirely on the uniqueness of the phase connection, which the manuscript does not yet document.

major comments (3)
  1. [Section 4.2 and Table 2] The uniqueness of the phase-connected timing solution is not established. The paper states that no solution was found with F0 and F1 alone and that adding F2 yielded a single solution, but it does not report the Dracula search space, the number of integer-cycle aliases tested, the reduced chi-squared of the accepted solution, or the fit quality of rejected aliases. Because the quoted F2 = 3.3133(6) x 10^-23 Hz/s^2 produces a cubic phase of order 14 turns over the ~4.3-year span quoted in the text (about 27 turns if the Table 2 start MJD 58404.483 is used), an alternative cycle count with slightly different F1, F2, position, and DM could plausibly fit the 65 TOAs. Please provide the search details and a comparison with rejected alias solutions; this is load-bearing for the claims that M15O is distinct from M15A and located 0.37 arcsec from the cluster center.
  2. [Table 2 and Table 1] Table 2 lists the start of timing data as MJD 58404.483, which corresponds to 2018 October 13, but Table 1 and Section 2 state that the data are 19 FAST observations from 2019 November (MJD 58796) to 2024 February. This one-year discrepancy changes the length of the phase-connected timing span and appears to be related to the statement in Section 4.2 that the observations span six years. Please correct the MJD or explain which additional data were used in the timing solution.
  3. [Section 3 and Table 2] The timing fit is documented only by the residual RMS of 108 microseconds; the manuscript does not give the TOA uncertainties, the number of fitted parameters, the reduced chi-squared, or the per-epoch residual rms. These quantities are needed to judge whether the quoted F2 and timing position are statistically acceptable and to compare against alternative alias solutions. Please add them to the paper.
minor comments (8)
  1. [Abstract and Section 3] The period difference between M15O and the 10th harmonic of M15A is stated as about 0.004 ms, but the quoted frequencies imply a difference of about 0.0004 ms; please correct the numerical value.
  2. [Table 1] The flux densities in Table 1 are quoted to three significant figures with no uncertainties, and they are derived from the radiometer equation with an assumed pulse width W_obs = 0.15 P_spin; please state the uncertainty budget or clearly label these values as estimates.
  3. [Section 4.1] The statement that 'at least 50% of experienced persons' did not see the signal in the diagnostic plot is informal; please replace it with a quantitative description of the test or remove the anecdotal percentage.
  4. [Section 4.3] The comparison with Galactic-field pulsars uses a 5-arcmin diameter and 5 cm^-3 DM difference criterion, but M15 pulsars are concentrated within a single FAST beam and have a DM spread of only a few tenths of pc cm^-3; the logical connection between the field-pulsar statistics and the cluster case should be clarified.
  5. [Section 2.2] The stack search method is essential to the discovery narrative, but its details are deferred to Dai et al. (2025, in preparation); please include a concise description of the stacking procedure and the candidate threshold in the present paper.
  6. [Abstract and Keywords] The pulsar name appears as J2129+1210O in the title and text but as J2129-1210O in the keywords; please use a consistent sign convention.
  7. [Figure 2 caption] Please clarify what exactly is meant by 'chopping the M15A signal' in the frequency domain and how the removed component was identified and validated.
  8. [Section 4.2] The sentence suggesting an orbital period on the order of decades based on F1 and F2 needs a shown calculation or should be explicitly labeled as speculative; if a wide orbit is possible, a direct binary or acceleration-model fit would be the appropriate test.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the pulsar discovery and timing solution are empirical outputs of independent search and fitting procedures.

full rationale

The paper reports an empirical discovery rather than a derivation. The central claim that M15O is an isolated millisecond pulsar rests on (1) a candidate detection in the power spectrum with SNR ~11, (2) confirmation by phase-coherent timing of 65 TOAs using TEMPO and the Dracula algorithm, and (3) a phase-connected timing solution whose parameters (spin frequency, F1, F2, DM, position) are fitted outputs from those TOAs, not pre-supplied inputs. The close proximity of M15O's period to the 10th harmonic of M15A is a comparison to an independently published timing solution for M15A (Wu et al. 2024), and none of the M15O parameters are defined in terms of that comparison. The need to include F2 to obtain a phase-connected solution is a modeling choice and a potential validity concern about aliases or alternative fits, but that is a correctness issue, not circularity: no equation in the paper reduces the claimed result to an input or to a fitted parameter renamed as a prediction. Self-citations (Pan et al. 2021, Wu et al. 2024, Yin et al. 2024) provide contextual or comparative material and are not load-bearing for the existence, position, or spin parameters of M15O. The paper is self-contained against external benchmarks in the sense that the timing fit uses standard software and procedures, and the discovery claim does not borrow its conclusion from the cited works.

Assumptions & free parameters 0 free parameters · 3 assumptions · 0 invented entities

The paper is an observational discovery rather than a parametric derivation, so there are no free constants tuned to make the central claim work. It relies on standard pulsar-timing assumptions: DE421 ephemeris, a cold-plasma dispersion relation, and a phase model with F0, F1, and F2. The assumption that M15O is isolated is not independently verified and could bias F1 and F2 if the object is actually in a long-period binary.

assumptions (3)
  • standard math DE421 solar system ephemeris and barycentric time are correct.
    Table 2 states DE421 and Barycentric Dynamic Time units; all TOAs are barycentered using this ephemeris, so any error in it would propagate into all spin and position parameters.
  • domain assumption Pulse phase over 2019 to 2024 is fully described by a constant spin frequency plus first and second derivatives, without glitches, timing noise, or unmodeled orbital motion.
    Section 4.2 says a fit with F0 and F1 alone could not connect all epochs and adding F2 yielded a single solution. If unmodeled processes contribute, F1, F2, and the derived position would be biased.
  • domain assumption The new signal follows the standard cold-plasma dispersion delay with one dispersion measure, so dedispersing at DM 67.44 separates it from M15A.
    The search dedispersed over DM 63 to 70 pc/cm^3 and the timing solution fits DM 67.44(1) pc/cm^3. This assumption underlies the interpretation that the candidate is an independent pulsar rather than a spurious harmonic artifact.

how reviews work

0 comments
Cite this review

Pith. "Pith review of The FAST Discovery of a Millisecond Pulsar M15O (PSR J2129+1210O) Hidden in the Harmonics of M15A (PSR J2129+1210A)." pith.science (2026). https://pith.science/paper/ZTZMIORN

@misc{pith2026250416872,
  author       = {Pith},
  title        = {Pith review of: The FAST Discovery of a Millisecond Pulsar M15O (PSR J2129+1210O) Hidden in the Harmonics of M15A (PSR J2129+1210A)},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/ZTZMIORN}},
  note         = {Machine review of arXiv:2504.16872}
}
abstract

We report the discovery of an isolated millisecond pulsar M15O (J2129+1210O) from the globular cluster M15 (NGC 7078) with a period of $\sim$11.06686 ms and a dispersion measure of $\sim$67.44 cm$^{-3}$ pc. Its spin period is so close to the 10th harmonic of the bright pulsar M15A ($\sim$11.06647 ms) that it was missed in previous pulsar search. We suggest adding the spectrum in the pulsar candidate diagnostic plot to identify new signals near the harmonics. M15O has the first spin frequency derivative and the second spin frequency derivative,being 1.79191(5) $\times$ $10^{-14}$ Hz $s^{-1}$ and 3.3133(6)$\times$ $10^{-23}$ Hz $s^{-2}$, respectively. Its projected distance from the optical center of M15 is the closest among all the pulsars in M15. The origin can be something from the center of the massive and core-collapsed globular cluster M15.

Discussion (0). Continue with ORCID to comment.

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score. Full citation record

  1. The FAST Globular Cluster Pulsar Survey (GC FANS)

    astro-ph.HE 2025-06 conditional novelty 6.0 of 10

    The FAST globular cluster pulsar survey reports 60 new pulsars and identifies M71D as a likely double neutron star system formed through standard massive binary evolution.

Reference graph

Works this paper leans on

19 extracted references · 3 canonical work pages · cited by 1 Pith paper

  1. [1]

    B., Gorham P

    Anderson S. B., Gorham P. W., Kulkarni S. R., Prince T. A., Wolszczan A., 1990, Natur, 346,

    work page 1990

  2. [29]

    V ., Ransom S

    doi:10.1088/0004-637X/795/1/29 2 Padmanabh P. V ., Ransom S. M., Freire P. C. C., Ridolfi A., Taylor J. D., Choza C., Clark C. J., et al., 2024, A&A, 686, A166. doi:10.1051/0004-6361/202449303 2 Pan Z., Hobbs G., Li D., Ridolfi A., Wang P., Freire P., 2016, MNRAS, 459, L26. doi:10.1093/mnrasl/slw037 2 10 Dai et al. Pan Z., Ma X.-Y ., Qian L., Wang L., Yan...

    doi:10.1088/0004-637x/795/1/29 2024

  3. [42]

    B., 1993, PhDT 2 Backer D

    doi:10.1038/346042a0 2 Anderson S. B., 1993, PhDT 2 Backer D. C., Kulkarni S. R., Heiles C., Davis M. M., Goss W. M., 1982, Natur, 300,

    doi:10.1038/346042a0 1993

  4. [43]

    R., Kramer M., 2012, hpa..book 7 Morello V ., Barr E

    doi:10.3847/1538-4357/ad5a82 2 Lorimer D. R., Kramer M., 2012, hpa..book 7 Morello V ., Barr E. D., Stappers B. W., Keane E. F., Lyne A. G., 2020, MNRAS, 497,

    doi:10.3847/1538-4357/ad5a82 2012

  5. [64]

    doi:10.1007/s11433-018-9376-1 3 Li B., Zhang L.-

    doi:10.1088/1674-4527/20/5/64 7 Jiang P., Yue Y ., Gan H., Yao R., Li H., Pan G., Sun J., et al., 2019, SCPMA, 62, 959502. doi:10.1007/s11433-018-9376-1 3 Li B., Zhang L.-. yun ., Yao J., Yin D., Eatough R. P., Li M., Li Y ., et al., 2024, ApJ, 972,

    doi:10.1088/1674-4527/20/5/64 2019

  6. [91]

    doi:10.1088/1674-4527/20/6/91 7

    work page doi:10.1088/1674-4527/20/6/91

  7. [125]

    B., Pennucci T., Ransom S

    doi:10.3847/1538-4357/aaac2a 2 Demorest P. B., Pennucci T., Ransom S. M., Roberts M. S. E., Hessels J. W. T., 2010, Natur, 467,

    doi:10.3847/1538-4357/aaac2a 2010

  8. [143]

    doi:10.3847/2041- 8213/ac0bbd 6 Wu Y ., Pan Z., Qian L., Ransom S

    doi:10.1088/1674-4527/21/6/143 3 Pan Z., Qian L., Ma X., Liu K., Wang L., Luo J., Yan Z., et al., 2021, ApJL, 915, L28. doi:10.3847/2041- 8213/ac0bbd 6 Wu Y ., Pan Z., Qian L., Ransom S. M., Eatough R. P., Wang B., Freire P. C. C., et al., 2024, ApJL, 974, L23. doi:10.3847/2041-8213/ad7b9e 4, 8 Ransom S. M., Eikenberry S. S., Middleditch J., 2002, AJ, 124,

    doi:10.1088/1674-4527/21/6/143 2021

Show all 19 references
  1. [584]

    A., Taylor J

    doi:10.1088/0004-637X/692/1/584 6 Hulse R. A., Taylor J. H., 1975, ApJL, 195, L51. doi:10.1086/181708 1 Harris W. E., 1996, AJ, 112,

    1975 doi

  2. [615]

    M., Freire P

    doi:10.1038/300615a0 1 B´egin S., 2006, MsT, 5 1 Cadelano M., Ransom S. M., Freire P. C. C., Ferraro F. R., Hessels J. W. T., Lanzoni B., Pallanca C., et al., 2018, ApJ, 855,

    2006 doi

  3. [709]

    O., Cohn H

    doi:10.1038/217709a0 1 Heinke C. O., Cohn H. N., Lugger P. M., 2009, ApJ, 692,

    2009 doi

  4. [857]

    doi:10.1093/mnras/stx1533 2 Freire P. C. C., Ridolfi A., 2018, MNRAS, 476,

    2018 doi

  5. [892]

    yun ., Qian L., Eatough R

    doi:10.1126/science.1108632 2, 5 Yin D., Zhang L.-. yun ., Qian L., Eatough R. P., Li B., Lorimer D. R., Dai Y ., et al., 2024, ApJL, 969, L7. doi:10.3847/2041-8213/ad534e 2 Yu Q.-Y ., Pan Z.-C., Qian L., Wang S., Yue Y .-L., Huang M.-L., Hao Q.-L., et al., 2020, RAA, 20,

    2024 doi

  6. [989]

    M., Ferraro F

    doi:10.1142/S0218271811019335 2 Pallanca C., Ransom S. M., Ferraro F. R., Dalessandro E., Lanzoni B., Hessels J. W. T., Stairs I., et al., 2014, ApJ, 795,

    2014 doi

  7. [1081]

    doi:10.1038/nature09466 1 Freire P. C. C., Ridolfi A., Kramer M., Jordan C., Manchester R. N., Torne P., Sarkissian J., et al., 2017, MNRAS, 471,

    2017 doi

  8. [1487]

    doi:10.1086/118116 2 Jiang P., Tang N.-Y ., Hou L.-G., Liu M.-T., Krˇco M., Qian L., Sun J.-H., et al., 2020, RAA, 20,

    2020 doi

  9. [1788]

    M., Hessels J

    doi:10.1086/342285 2, 3 Ransom S. M., Hessels J. W. T., Stairs I. H., Freire P. C. C., Camilo F., Kaspi V . M., Kaplan D. L., 2005, Sci, 307,

    2005 doi

  10. [4654]

    ascl:1509.002 2 Nan R., Li D., Jin C., Wang Q., Zhu L., Zhu W., Zhang H., et al., 2011, IJMPD, 20,

    doi:10.1093/mnras/staa2291 2 Nice D., Demorest P., Stairs I., Manchester R., Taylor J., Peters W., Weisberg J., et al., 2015, ascl.soft. ascl:1509.002 2 Nan R., Li D., Jin C., Wang Q., Zhu L., Zhu W., Zhang H., et al., 2011, IJMPD, 20,

    2015 doi

  11. [4794]

    J., Pilkington J

    doi:10.1093/mnras/sty524 5 Hewish A., Bell S. J., Pilkington J. D. H., Scott P. F., Collins R. A., 1968, Natur, 217,

    1968 doi

Pith tools

Reviewed August 16, 2026 · model on record in the stance chip above.