Pith. sign in

REVIEW 4 minor 2 cited by

Gamma-ray pulsations are detected from the extreme-spin-down millisecond pulsar PSR J0435+3233, establishing it as a gamma-ray MSP with a very low gamma-ray efficiency.

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

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

T0 review · deepseek-v4-flash

2026-08-01 21:15 UTC pith:KVTG6XSJ

load-bearing objection A clean, externally-anchored Fermi-LAT detection of gamma-ray pulsations from an extreme MSP; the 6.8σ is probably real but the paper should quantify the two-interval trials correction.

arxiv 2607.16119 v1 pith:KVTG6XSJ submitted 2026-07-17 astro-ph.HE

Discovery of γ-Ray Pulsations from the Extreme-Spin-Down Millisecond Pulsar PSR J0435+3233

Mengqing Zhang , Shengbin Pei , Pengfei Zhang This is my paper
classification astro-ph.HE
keywords gamma-ray pulsarsmillisecond pulsarsspin-down luminositygamma-ray efficiencypulsar timingPSR J0435+32334FGL J0435.5+3232gamma-ray counterpart
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

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

The paper reports the detection of pulsed gamma-ray emission from PSR J0435+3233, a millisecond pulsar with an unusually large spin-down rate and a spin-down luminosity comparable to young pulsars. Using about 17.7 years of gamma-ray data but only the interval where the radio timing solution is valid, the authors find a signal at roughly 6.8 sigma significance. They identify the cataloged source 4FGL J0435.5+3232 as the gamma-ray counterpart based on positional coincidence and phase-dependent emission. The inferred gamma-ray luminosity is about 6.26e32 erg/s, which is roughly five orders of magnitude below the spin-down luminosity, corresponding to an apparent gamma-ray efficiency of about 1.1e-5. If correct, this makes the pulsar a striking outlier that tests how magnetospheric gamma-ray emission depends on spin-down power, magnetic field, and viewing geometry.

Core claim

The central claim is that PSR J0435+3233 is a gamma-ray millisecond pulsar. The evidence is the detection of gamma-ray pulsations at the radio spin period with a weighted H-test statistic of 63.9 (p = 7.6e-12, about 6.8 sigma) using events inside the ephemeris-valid interval, the positional coincidence with 4FGL J0435.5+3232, and the pulsed emission being concentrated in the rotational-phase interval 0.44-0.69. The paper further derives a gamma-ray luminosity and an apparent efficiency that are orders of magnitude lower than those of typical MSPs, a combination not seen before.

What carries the argument

The analysis hinges on the radio timing ephemeris of PSR J0435+3233, which provides the rotational phase for each gamma-ray photon. Photon weights from a likelihood source model are used in a weighted H-test to assess periodicity. The key diagnostics are the pulse profile, the phase-time diagram, and on-pulse/off-pulse TS maps, which together establish the association and the pulse's location in phase. The machinery is the phase-coherent folding of weighted photon arrival times, with the validity interval of the ephemeris defining the data window.

Load-bearing premise

The reported significance depends on the radio timing ephemeris being phase-coherent and accurate over MJD 59101-60898, and on the decision to analyze only this interval; if the ephemeris has unmodeled phase noise or the interval choice is counted as a search, the effective significance would be lower.

What would settle it

Reanalyze the same gamma-ray data with an updated or independently derived radio timing ephemeris covering the same dates; if the H-test value does not stay near 63.9 or the pulse phase shifts, the detection is not robust. Alternatively, recompute the p-value after including a trials factor for testing both the full 17.7-year interval and the ephemeris-valid interval.

Watch this falsifier — get emailed when new claim-graph text bears on it.

If this is right

  • PSR J0435+3233 joins the population of gamma-ray millisecond pulsars, adding an extreme point in period-derivative and spin-down power.
  • Its apparent gamma-ray efficiency of about 1e-5 is among the lowest seen in MSPs, implying that high spin-down power alone does not guarantee high GeV output.
  • The confined pulse profile and absence of off-pulse emission indicate that the pulsed magnetospheric component dominates the detected GeV emission.
  • Absolute phase alignment of radio and gamma-ray profiles in the future will constrain the emission geometry and may explain the low efficiency as a beaming effect.
  • The lack of significant orbital modulation leaves room for a possible intrabinary shock component, but the current data cannot discriminate it from magnetospheric emission.

Where Pith is reading between the lines

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

  • No trials factor is applied for the two-interval comparison (full data versus ephemeris-valid interval); treating the interval choice as a search would reduce the formal significance, though it likely remains above 5 sigma.
  • If the low efficiency is intrinsic rather than geometric, existing scaling relations between gamma-ray luminosity and spin-down power for MSPs may need a secondary dependence on magnetic field or evolutionary history.
  • The extreme spin-down luminosity combined with a millisecond period and possible non-standard formation suggests PSR J0435+3233 could be a rare transition object; gamma-ray pulsation detection supports magnetospheric emission operating in such systems.
  • A longer timing campaign that provides a continuous phase-coherent ephemeris over the whole data set would test whether the pulsation exists beyond the current ephemeris-valid interval, clarifying whether the low signal outside that interval is due to ephemeris decoherence or a genuine cessation of pulsed emission.

Editorial analysis

A structured set of objections, weighed in public.

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

Referee Report

0 major / 4 minor

Summary. The paper analyzes ~17.7 yr of Fermi-LAT data toward the extremely fast-spinning-down millisecond pulsar PSR J0435+3233. Using the radio timing ephemeris of Wu et al. (2026) and restricting to its validity interval (MJD 59101–60898), the authors detect gamma-ray pulsations at a weighted H-test statistic H=63.9, corresponding to p=7.6e-12 (~6.8σ). They identify the cataloged source 4FGL J0435.5+3232 as the gamma-ray counterpart based on positional coincidence, pulsations at the spin period, and phase-dependent emission. The pulse profile shows a prominent peak at phases 0.44–0.69; phase-resolved analysis finds the emission concentrated there. The phase-averaged EVI flux yields L_gamma = 6.26e32 erg/s (d=1.2 kpc, isotropic) and an apparent efficiency of 1.1e-5, exceptionally low for a pulsar with a young-pulsar-like spin-down power. The paper discusses physical implications, including beaming and intrinsic emission models, and explicitly notes limitations in phase alignment and orbital variability.

Significance. If the detection is correct, it establishes PSR J0435+3233 as a gamma-ray MSP with extreme spin-down properties, adding a unique data point to the P–Pdot plane and to studies of magnetospheric efficiency in recycled pulsars. The very low apparent efficiency (~1e-5) compared to typical MSPs (~0.1–1%) is a striking result that, even with substantial distance or beaming uncertainties, robustly indicates suppressed gamma-ray output. The analysis uses standard, well-established methodology: weighted H-test, gtsrcprob weights, and externally motivated ephemeris validity interval. The paper also provides a phase–time diagram, cumulative H-test curve, and on/off-pulse TS maps that support the detection and the phase-resolved conclusions.

minor comments (4)
  1. [§2.2] The quoted p-value (p=7.6e-12) conditions on the EVI being the only interval analyzed, but the text reports an earlier exploratory whole-data search (H~9.4). Though the EVI is externally fixed by the radio ephemeris, which mitigates the look-elsewhere effect, the paper should explicitly state that the EVI selection is a priori with respect to the LAT data, or provide a simple trial-corrected p-value. As written, a reader could perceive a post-hoc interval selection. This is a presentational but important statistical clarification.
  2. [§2.1, Eq. (1)] The PLEC4 model equation appears to contain a typographical error: the power-law index is written as Γ0 + d/b, but standard PLEC4 (as implemented in Fermitools) uses Γ0 alone in the power-law factor, with the exponential cutoff term carrying the d/b^2 dependence. Please correct the equation or explicitly define the parameterization if the exponent is intentional.
  3. [§2.2 and §2.3] The on-pulse interval (ϕ=0.4375–0.6875) is defined from the same pulse profile used for the phase-resolved analysis. This is a common and accepted practice, but the statement that 'the pulsed emission is concentrated predominantly' in this interval is partly a restatement of the definition. The claim is strengthened by the off-pulse TS=0 result, which is independent. Please add a sentence acknowledging this descriptive circularity.
  4. [Abstract and §1] The abstract says '~17 yr' of observations while the text says 17.7 yr; please harmonize. Also, the abstract and intro mention the efficiency as an 'exceptionally low' value; it would be useful to state explicitly that the quoted L_gamma and efficiency assume f_Omega=1, a caveat already present in §3.

Circularity Check

0 steps flagged

No significant circularity: detection rests on an external radio ephemeris and the 4FGL catalog, with significance from a full-profile weighted H-test.

full rationale

The paper's derivation chain is anchored to external inputs: the radio timing solution of Wu et al. (2026) defines the phase model, the 4FGL-DR4 catalog provides the candidate γ-ray counterpart position, and the spin-down luminosity and distance come from the radio discovery and DM estimate. The pulsation significance is evaluated with the weighted H-test over the full pulse profile (H = 63.9, p = 7.6e-12), not against a region chosen after seeing the data. The on-pulse region is indeed defined from the observed profile ('The γ-ray profile exhibits a prominent peak within the phase interval ϕ = 0.4375–0.6875. We therefore defined this interval as the on-pulse region'), so describing the emission as concentrated there is to some extent descriptive tautology, but the pulsed detection itself does not depend on this definition. The PLEC4 spectral model is adopted from standard pulsar spectral fitting, not constructed to force pulsations. The restriction to the ephemeris-valid interval is externally motivated by the stated validity of the radio solution; any concern about a trials factor for first folding the whole dataset is a statistical robustness issue, not circularity. No load-bearing step reduces to the paper's own conclusions, and no self-citation chain is used to justify the central claim. The paper explicitly acknowledges its limitations (unabsolute phase alignment, distance/beaming dependence of efficiency, inconclusive orbital modulation), further supporting that the derivation is not circular.

Axiom & Free-Parameter Ledger

6 free parameters · 6 axioms · 0 invented entities

The paper introduces no new physical entities. It relies on externally measured inputs (radio ephemeris, distance, spin-down luminosity) and on standard Fermi-LAT analysis assumptions. The only data-defined choice is the on-pulse phase interval, which is descriptive rather than a fitted model parameter.

free parameters (6)
  • Distance d = 1.2 kpc
    DM-inferred distance taken from Wu et al. (2026); enters L_gamma and eta_gamma linearly. No uncertainty propagated.
  • Beaming factor f_Omega = 1
    Isotropic emission assumed; the paper acknowledges f_Omega could be much less than 1, which would lower the already low efficiency.
  • On-pulse phase boundaries = 0.4375-0.6875
    Defined from the pulse profile of the same dataset; used for phase-resolved spectral fits and TS maps.
  • PLEC4 spectral parameters Gamma, d = Gamma=1.93(20), d=0.17(9) (whole); Gamma=2.12(33), d=0.13(9) (EVI)
    Fitted to LAT data. These shape the photon weights used in the pulsation search, though small changes do not drive the detection.
  • PLEC4 sub-exponential index b = 2/3
    Fixed to a canonical pulsar spectral value, not fitted.
  • Spin-down luminosity Edot and Pdot = Edot=5.89e37 erg/s, Pdot=4.88e-17 s/s
    Taken from Wu et al. (2026) as inputs; used for the apparent efficiency. Not fitted in this paper.
axioms (6)
  • domain assumption Radio timing ephemeris of Wu et al. (2026) is phase-coherent only within MJD 59101-60898; extrapolation loses phase coherence.
    Used to restrict the pulsation search to the EVI. If false, the 6.8 sigma significance is not valid.
  • domain assumption The gamma-ray source 4FGL J0435.5+3232 is the counterpart of PSR J0435+3233.
    Supported by 0.01 degree positional coincidence and pulsations; used to assign photon weights.
  • domain assumption Distance to the pulsar is 1.2 kpc (DM-inferred).
    Used for L_gamma and eta_gamma; unquantified uncertainty.
  • domain assumption Emission is isotropic (f_Omega=1).
    Used for luminosity; the paper acknowledges it may be wrong.
  • standard math Weighted H-test null distribution and gtsrcprob weights are valid under the Fermi-LAT source model.
    Standard Fermi-LAT methodology; the H-test significance depends on the assumed null distribution.
  • domain assumption PLEC4 is the appropriate spectral form for pulsar magnetospheric emission.
    Used for the source model and photon weights; b is fixed to 2/3.

pith-pipeline@v1.3.0-alltime-deepseek · 7615 in / 12497 out tokens · 101590 ms · 2026-08-01T21:15:52.541571+00:00 · methodology

0 comments
read the original abstract

Motivated by the recent discovery of PSR~J0435+3233, a millisecond pulsar with an exceptionally large period derivative of $\dot{P}=4.9\times10^{-17}\ {\rm s\,s^{-1}}$ and a high spin-down luminosity of $\dot{E}=5.89\times10^{37}\ {\rm erg\,s^{-1}}$, we analyze $\sim$17~yr of observations obtained with the \textit{Fermi} Large Area Telescope for the pulsar. We identify the cataloged source 4FGL~J0435.5+3232, located only $0.01^{\circ}$ from the radio timing position, as the $\gamma$-ray counterpart of PSR~J0435+3233. Using $\gamma$-ray events collected within the validity interval of the radio timing ephemeris, we detect its $\gamma$-ray pulsations at a $\sim6.8\sigma$ confidence level. In this interval, the phase-resolved analysis shows that the pulsed emission is concentrated predominantly within the rotational-phase of $\phi\sim0.44$--$0.69$. We also derive its $\gamma$-ray luminosity of $L_{\gamma}=6.26\times10^{32}\ {\rm erg\,s^{-1}}$, assuming a distance of 1.2~kpc and isotropic emission. This luminosity corresponds to an apparent $\gamma$-ray efficiency of $\eta_{\gamma}\sim1.1\times10^{-5}$, revealing an exceptionally low $\gamma$-ray output despite the pulsar's young-pulsar-like rotational-energy budget. Our detection establishes PSR~J0435+3233 as a $\gamma$-ray MSP, and the striking combination of its high spin-down power and low apparent $\gamma$-ray efficiency provides a new probe of particle acceleration, radiation beaming, and viewing geometry in the magnetospheres of millisecond pulsars with extreme rotational properties.

Figures

Figures reproduced from arXiv: 2607.16119 by Mengqing Zhang, Pengfei Zhang, Shengbin Pei.

Figure 1
Figure 1. Figure 1: Fermi-LAT timing analysis of target over the EVI. (A) Integrated weighted γ-ray pulse profile. The horizontal dashed line indicates the estimated background level. (B) Two-dimensional phase–time diagram of the weighted photons. The color scale represents the probability weight assigned to each photon, with a maximum weight of ∼0.85. (C) Cumulative weighted H-test statistic as a function of time. (D) Schema… view at source ↗
Figure 2
Figure 2. Figure 2: γ-ray SEDs of target obtained from the whole data, the EVI, and the on-pulse selection within the EVI, shown in green, blue, and black, respectively. The best￾fit PLEC4 and LP models are represented by the solid and dashed curves, respectively. Flux measurements are shown for energy bins with TS ⩾ 4, whereas 95% upper limits are given for bins with TS < 4. stantial fraction of the rotational energy could i… view at source ↗
Figure 3
Figure 3. Figure 3: TS maps of the 5 ◦ ×5 ◦ region centered on PSR J0435+3233 in the 0.1–500.0 GeV, constructed using events collected within the EVI. Left and right panels show the on-pulse and off-pulse TS maps, respectively. Both maps have a pixel size of 0.1 ◦ × 0.1 ◦ . In each panel, the green circle indicates the 1σ positional uncertainty region of 4FGL J0435.5+3232, and the white cross marks the radio timing position o… view at source ↗

discussion (0)

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

Forward citations

Cited by 2 Pith papers

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

  1. On the triple nature of the PSR J0435+3233 system

    astro-ph.HE 2026-07 conditional novelty 7.0

    PSR J0435+3233 is a hierarchical triple: its anomalous spin-down is an acceleration artifact from a ~1.2 M⊙ outer companion on a ~70-year eccentric orbit.

  2. Search for continuous gravitational waves from the pulsar J0435+3233

    gr-qc 2026-07 accept novelty 6.0

    A LIGO O4a search for continuous gravitational waves from millisecond pulsar J0435+3233 finds no signal, setting h0<5.8×10^-27 at 95% confidence and an ellipticity limit of 1.6×10^-8.

Reference graph

Works this paper leans on

21 extracted references · 1 canonical work pages · cited by 2 Pith papers

  1. [1]

    A., Ackermann , M., Ajello , M., et al

    Abdo , A. A., Ackermann , M., Ajello , M., et al. 2010, , 708, 1254, 10.1088/0004-637X/708/2/1254

    doi:10.1088/0004-637x/708/2/1254 2010

  2. [2]

    A., Ajello , M., Allafort , A., et al

    Abdo , A. A., Ajello , M., Allafort , A., et al. 2013, , 208, 17, 10.1088/0067-0049/208/2/17

    doi:10.1088/0067-0049/208/2/17 2013

  3. [3]

    M., Stairs , I

    Archibald , A. M., Stairs , I. H., Ransom , S. M., et al. 2009, Science, 324, 1411, 10.1126/science.1172740

    doi:10.1126/science.1172740 2009

  4. [4]

    C., Kulkarni , S

    Backer , D. C., Kulkarni , S. R., Heiles , C., Davis , M. M., & Goss , W. M. 1982, , 300, 615, 10.1038/300615a0

    doi:10.1038/300615a0 1982

  5. [5]

    H., Lott , B., & The Fermi-LAT collaboration

    Ballet , J., Bruel , P., Burnett , T. H., Lott , B., & The Fermi-LAT collaboration . 2023, arXiv e-prints, arXiv:2307.12546, 10.48550/arXiv.2307.12546

    Pith/arXiv arXiv doi:10.48550/arxiv.2307.12546 2023

  6. [6]

    Bhattacharya , D., & van den Heuvel , E. P. J. 1991, , 203, 1, 10.1016/0370-1573(91)90064-S

    doi:10.1016/0370-1573(91)90064-s 1991

  7. [7]

    C., & B \"u sching , I

    de Jager , O. C., & B \"u sching , I. 2010, , 517, L9, 10.1051/0004-6361/201014362

    doi:10.1051/0004-6361/201014362 2010

  8. [8]

    Z., Chang , J., Guo , J

    Fan , Y. Z., Chang , J., Guo , J. H., et al. 2022, Acta Astronomica Sinica, 63, 27

    2022

  9. [9]

    2006, Chinese Journal of Astronomy and Astrophysics Supplement, 6, 189

    Hobbs , G., Edwards , R., & Manchester , R. 2006, Chinese Journal of Astronomy and Astrophysics Supplement, 6, 189

    2006

  10. [10]

    M., Hessels , J

    Jaodand , A., Archibald , A. M., Hessels , J. W. T., et al. 2016, , 830, 122, 10.3847/0004-637X/830/2/122

    work page doi:10.3847/0004-637x/830/2/122 2016

  11. [11]

    2011, , 732, 38, 10.1088/0004-637X/732/1/38

    Kerr , M. 2011, , 732, 38, 10.1088/0004-637X/732/1/38

    doi:10.1088/0004-637x/732/1/38 2011

  12. [12]

    Lomb , N. R. 1976, , 39, 447, 10.1007/BF00648343

    doi:10.1007/bf00648343 1976

  13. [13]

    Lorimer , D. R. 2008, Living Reviews in Relativity, 11, 8, 10.12942/lrr-2008-8

    doi:10.12942/lrr-2008-8 2008

  14. [14]

    N., Hobbs , G

    Manchester , R. N., Hobbs , G. B., Teoh , A., & Hobbs , M. 2005, , 129, 1993, 10.1086/428488

    doi:10.1086/428488 2005

  15. [15]

    2013, , 501, 517, 10.1038/nature12470

    Papitto , A., Ferrigno , C., Bozzo , E., et al. 2013, , 501, 517, 10.1038/nature12470

    doi:10.1038/nature12470 2013

  16. [16]

    Patruno , A., & Watts , A. L. 2021, in Astrophysics and Space Science Library, Vol. 461, Timing Neutron Stars: Pulsations, Oscillations and Explosions, ed. T. M. Belloni , M. M \'e ndez , & C. Zhang , 143--208, 10.1007/978-3-662-62110-3_4

    doi:10.1007/978-3-662-62110-3_4 2021

  17. [17]

    S., Kerr , M., Parent , D., et al

    Ray , P. S., Kerr , M., Parent , D., et al. 2011, , 194, 17, 10.1088/0067-0049/194/2/17

    doi:10.1088/0067-0049/194/2/17 2011

  18. [18]

    Scargle , J. D. 1982, , 263, 835, 10.1086/160554

    doi:10.1086/160554 1982

  19. [19]

    A., Abdollahi , S., Ajello , M., et al

    Smith , D. A., Abdollahi , S., Ajello , M., et al. 2023, , 958, 191, 10.3847/1538-4357/acee67

    doi:10.3847/1538-4357/acee67 2023

  20. [20]

    1998, , 394, 344, 10.1038/28557

    Wijnands , R., & van der Klis , M. 1998, , 394, 344, 10.1038/28557

    doi:10.1038/28557 1998

  21. [21]

    2026, Nature Astronomy, 10.1038/s41550-026-02836-3

    Wu , Q., Wang , N., Yuan , J., et al. 2026, Nature Astronomy, 10.1038/s41550-026-02836-3

    doi:10.1038/s41550-026-02836-3 2026