Point-like Off-pulse GeV Emission from the Millisecond Pulsar PSR J0437-4715

Ziwei Ou

arxiv: 2604.26442 · v1 · submitted 2026-04-29 · 🌌 astro-ph.HE

Point-like Off-pulse GeV Emission from the Millisecond Pulsar PSR J0437-4715

Ziwei Ou This is my paper

Pith reviewed 2026-05-07 12:37 UTC · model grok-4.3

classification 🌌 astro-ph.HE

keywords PSR J0437-4715millisecond pulsaroff-pulse emissionGeV gamma-raysFermi-LATtermination shockbow-shock PWNouter gap model

0 comments

The pith

Off-pulse GeV emission from PSR J0437-4715 appears point-like with an upper limit radius of 0.12 degrees.

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

The paper uses Fermi-LAT observations and the pulsar's timing ephemeris to separate on-pulse and off-pulse phases in PSR J0437-4715. A binned likelihood analysis finds no detectable spatial extension in the off-pulse GeV signal. The conversion efficiency from spin-down power to gamma-ray luminosity matches the outer gap model when the termination shock sits close to the pulsar. A reader would care because this constrains whether the steady off-pulse emission comes from the pulsar itself or from a larger nebula, and it links the signal to bow-shock interactions seen in X-rays.

Core claim

The off-pulse GeV emission is consistent with a point source, with no spatial extension detected above an upper limit of 0.12 degrees. Spectral analysis across phases shows that the gamma-ray luminosity and conversion efficiency from spin-down power align with an outer gap model in which the termination shock lies close to the pulsar. Comparison with a sample of bow-shock PWNe further suggests the off-pulse signal may connect to the pulsar's wind nebula rather than an unrelated background.

What carries the argument

Phase-resolved binned likelihood analysis of Fermi-LAT data, using pulsar timing ephemeris to isolate off-pulse intervals and test for spatial extension against outer gap model predictions.

If this is right

The off-pulse GeV emission originates close to the pulsar, consistent with a termination shock in the outer gap scenario.
Gamma-ray luminosity scales with X-ray luminosity in the same way seen for other bow-shock PWNe.
The efficiency of spin-down power conversion to GeV photons matches expectations only when the shock is nearby rather than distant.
Similar off-pulse components may appear in other millisecond pulsars that drive bow shocks.

Where Pith is reading between the lines

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

Higher-resolution instruments could test whether the 0.12-degree limit tightens or breaks, directly probing the shock location.
If the association holds, the total gamma-ray output of millisecond pulsars includes a steady off-pulse fraction that may contribute to the diffuse background.
Cross-checking the same phase separation method on other MSPs with known bow shocks would show how common this point-like off-pulse signature is.

Load-bearing premise

The timing ephemeris cleanly divides on-pulse and off-pulse phases, and the detected off-pulse GeV signal belongs to the pulsar instead of unrelated background.

What would settle it

Deeper Fermi-LAT data or higher-resolution observations that either detect spatial extension beyond 0.12 degrees or show the off-pulse flux disappearing when phase intervals are shifted by even a few percent.

Figures

Figures reproduced from arXiv: 2604.26442 by Ziwei Ou.

**Figure 1.** Figure 1: PSR J0437−4715 timing results from Tempo2 with the Fermi plug-in. Top-left panel: phase histogram of the analyzed Fermi-LAT data. Two full rotational phase are shown here. Bottom-left panel: H-test significance (TS) as a function of time. Right panel: pulse phase for each 𝛾-ray event vs. time. Phase Range Flux Γ 𝑑 TS log(L) (ph s−1 cm−2 ) whole (1.04 ± 0.02) × 10−08 1.82 ± 0.05 0.82 ± 0.05 7029 -141844.251… view at source ↗

**Figure 2.** Figure 2: Weighted pulse profile of PSR J0437−4715 at different energies. ing 𝐿𝛾 ∼ 𝑅 𝛼 bs with 𝛼 > 2 would point to an enhanced acceleration efficiency or reduced cooling losses in larger bow shocks. However, 𝛾-ray emission dominated by particles accelerated from the ISM can not be ruled out. The 𝛾-ray luminosity is also governed by other factors probably, such as the magnetic field strength, the injection spectral … view at source ↗

**Figure 3.** Figure 3: 𝛾-ray TS map (300 MeV to 300 GeV) of PSR J0437−4715 from the whole data (left) and the off-pulse data (right). PSR 𝑑 log𝐸¤ log 𝜏 𝐵11 𝑣⊥ 𝑟bs log 𝐿X log 𝐿𝛾 kpc erg/s yr 1011 G km/s 106 cm erg/s erg/s J1952+3251 3 36.57 5.03 4.86 460 < 12 33.02 ± 0.11 35.20 ± 0.08 J1826−1256 3.9 36.56 4.16 37 ... ... 33.38 ± 0.06 < 37.25 J1709−4429 2.6 36.53 4.24 31.2 ≲ 100 ∼ 70 32.60 ± 0.10 36.05 ± 0.03 J1801−2451 3.8 36.41 … view at source ↗

**Figure 4.** Figure 4: 𝛾-ray spectra (300 MeV to 300 GeV) of PSR J0437−4715 in the on-pulse data and off-pulse data. The on-pulse data was modeled by a PLSC model while the off-pulse data was modeled by a PL model. Statistical relations from a bow-shock PWN sample imply 𝐿𝛾 ∼ 𝐿 0.61±0.21 𝑋 and 𝐿𝛾 ∼ 𝑅 1.70±0.73 bs . The former suggests that IC scattering of external photon fields likely dominates the 𝛾-ray emission, although SSC … view at source ↗

**Figure 5.** Figure 5: Left panel: 𝛾-ray luminosity 𝐿𝛾 vs. X-ray luminosity 𝐿𝑋. Right panel: 𝛾-ray luminosity 𝐿𝛾 vs. bow-shock radius 𝑅bs. The dashed lines and the gray shadows provide the best-fit results of each panel view at source ↗

**Figure 6.** Figure 6: Spectral cutoff energy 𝐸cut values (red points) as a function of pulse phase. The blue line show the weighted counts of PSR J0437−4715. Off-pulse regions are shown as gray shadow view at source ↗

**Figure 7.** Figure 7: Linear correlation between 𝛾-ray luminosity log10 𝐸cut and cutoff energy log10 𝐿 in different phase bin. The line and shadow represent the best-fit results with error. The spectral cutoff energy values are provided in view at source ↗

read the original abstract

PSR J0437-4715 is a gamma-ray millisecond pulsar, which has been detected by Fermi-LAT. For understanding the nature, we analyze the GeV gamma-ray data obtained with Fermi-LAT around the pulsar region. Based on the pulsar timing ephemeris, we derived the gamma-ray pulse profile and defined on-pulse and off-pulse phase intervals. A binned likelihood analysis was performed to investigate the spectral properties of the pulsar across different phase ranges. No spatial extension was detected for off-pulse, with an upper limit radius of 0.12 degree. We further investigate the relationship between gamma-ray luminosity, X-ray luminosity, and bow-shock radius for a sample of pulsars with detected bow-shock PWN. The relationship between gamma-ray luminosity and X-ray luminosity is explored. The conversion efficiency from spin-down power to GeV emission of outer gap model is consistent with a termination shock located close to the pulsar. We discuss the potential nature of off-pulse GeV emission and the connection to bow shock pulsar wind nebulae.

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.

Desk Editor's Note private letter to a colleague

The paper reports a 0.12-degree upper limit on off-pulse GeV extension for PSR J0437-4715 with an outer-gap consistency check, but the off-pulse excess significance is not quantified.

read the letter

This paper gives a new 0.12-degree upper limit on the spatial size of off-pulse GeV emission from PSR J0437-4715 and checks that the outer-gap efficiency is consistent with a termination shock close to the pulsar. The analysis is a standard phase-resolved binned likelihood fit using Fermi-LAT data and the pulsar's timing ephemeris. It does a decent job of applying familiar techniques to this well-studied millisecond pulsar and adding a bow-shock PWN sample comparison for context on luminosities. The specific limit is tighter than some earlier work and could help constrain emission models for this object. The main weakness is the lack of reported details on the off-pulse signal strength. There is no test statistic, no residual map, and no clear statement on how the Galactic diffuse model and other sources were handled in the off-pulse phase. Without those, it's hard to know if the point-like excess is convincingly detected or could be background. The stress-test note flags this correctly, and it affects how much weight to give the model consistency claim. The phase intervals seem reasonably defined, but the spatial and spectral separation step is where the evidence is thinnest. This work is for pulsar emission modelers who track individual sources and bow-shock connections. It is incremental rather than transformative. I would recommend sending it to peer review because the full methods might resolve the gaps and the limit itself is a concrete addition worth checking.

Referee Report

2 major / 1 minor

Summary. The manuscript analyzes Fermi-LAT GeV gamma-ray data for the millisecond pulsar PSR J0437-4715. Using the pulsar timing ephemeris to define on-pulse and off-pulse phase intervals, a binned likelihood analysis is performed to examine spectral properties across phase ranges. The central results are that no spatial extension is detected for the off-pulse emission (upper limit radius 0.12 degrees) and that the spin-down to GeV conversion efficiency in the outer gap model is consistent with a termination shock located close to the pulsar. The paper additionally explores gamma-ray versus X-ray luminosity relations for a sample of bow-shock PWN pulsars and discusses possible connections between the off-pulse emission and bow-shock nebulae.

Significance. If the off-pulse GeV excess is statistically robust and spatially coincident with the pulsar, the result would add a useful data point to studies of millisecond pulsar emission mechanisms, particularly by constraining the termination shock location via the reported efficiency and extension limit. The sample comparison with other bow-shock PWNe could help contextualize the findings. However, the current lack of quantitative detection metrics limits the immediate impact.

major comments (2)

[Abstract] The description of the binned likelihood analysis (Abstract) does not report the test statistic (TS) value, detection significance, or any residual maps for the off-pulse phase interval. Without these, it is not possible to verify that the claimed point-like excess is pulsar-associated rather than residual diffuse or unmodeled background emission, which directly undermines the reliability of the 0.12-degree extension upper limit.
[Abstract] No information is given on the energy range, the Galactic diffuse emission model, the 4FGL source catalog version, or systematic uncertainties in the likelihood fits (Abstract). These details are required to assess whether the background is fully accounted for in the off-pulse phase range and to support the spectral and spatial conclusions.

minor comments (1)

[Abstract] The abstract would be clearer if it specified the exact phase intervals used for on-pulse and off-pulse gating and the energy range of the analysis.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful and constructive review of our manuscript. We have addressed the concerns about the abstract by incorporating the requested quantitative details and methodological information into a revised version of the abstract.

read point-by-point responses

Referee: [Abstract] The description of the binned likelihood analysis (Abstract) does not report the test statistic (TS) value, detection significance, or any residual maps for the off-pulse phase interval. Without these, it is not possible to verify that the claimed point-like excess is pulsar-associated rather than residual diffuse or unmodeled background emission, which directly undermines the reliability of the 0.12-degree extension upper limit.

Authors: We agree that the abstract should include these metrics to strengthen the presentation of the results. In the revised manuscript we have added the test statistic value and corresponding detection significance for the off-pulse emission. We have also clarified that residual maps were examined as part of the analysis and show no significant unmodeled emission, supporting the association with the pulsar. The 0.12-degree extension upper limit is obtained from a likelihood-ratio test comparing point-like and spatially extended source models; we have expanded the abstract text to make this procedure explicit. revision: yes
Referee: [Abstract] No information is given on the energy range, the Galactic diffuse emission model, the 4FGL source catalog version, or systematic uncertainties in the likelihood fits (Abstract). These details are required to assess whether the background is fully accounted for in the off-pulse phase range and to support the spectral and spatial conclusions.

Authors: We thank the referee for highlighting these omissions. The revised abstract now states the energy range used for the likelihood analysis, the Galactic diffuse emission model adopted, the version of the 4FGL catalog employed, and the approach taken to evaluate systematic uncertainties (by varying the diffuse model and source parameters). These additions provide the necessary context for readers to assess the background treatment and the robustness of the reported spectral and spatial results. revision: yes

Circularity Check

0 steps flagged

No circularity: standard observational analysis with external model comparison

full rationale

The paper conducts a Fermi-LAT binned likelihood analysis on phase-gated data using an external pulsar timing ephemeris, reports a non-detection of spatial extension (upper limit 0.12°), and compares derived luminosities to an independent sample of bow-shock PWNe. No equations, fitted parameters, or derivations are presented that reduce the claimed efficiency, location, or point-like nature to quantities defined by the same dataset. The outer-gap model consistency statement invokes an external model and sample statistics rather than any self-referential fit or ansatz. This is a self-contained observational report; the central claims rest on data reduction and external benchmarks, not internal redefinition.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

No free parameters, axioms, or invented entities are introduced; the work consists of observational data analysis and comparison to an existing outer-gap model.

pith-pipeline@v0.9.0 · 5480 in / 1108 out tokens · 37755 ms · 2026-05-07T12:37:24.427046+00:00 · methodology

discussion (0)

Reference graph

Works this paper leans on

2 extracted references · 1 canonical work pages

[1]

2019, Science with the Cherenkov Telescope Array, doi: 10.1142/10986

Abeysekara A. U., et al., 2017, Science, 358, 911 Ackermann M., et al., 2011, ApJ, 726, 35 Aldcroft T. L., Romani R. W., Cordes J. M., 1992, ApJ, 400, 638 Auchettl K., et al., 2015, ApJ, 802, 68 BalletJ.,BruelP.,BurnettT.H.,LottB.,TheFermi-LATcollaboration2023, arXiv e-prints, p. arXiv:2307.12546 Bednarek W., Sitarek J., 2013, MNRAS, 430, 2951 Brownsberge...

work page doi:10.1142/10986 2017
[2]

MNRAS000, 1–6 (2026)

2026

[1] [1]

2019, Science with the Cherenkov Telescope Array, doi: 10.1142/10986

Abeysekara A. U., et al., 2017, Science, 358, 911 Ackermann M., et al., 2011, ApJ, 726, 35 Aldcroft T. L., Romani R. W., Cordes J. M., 1992, ApJ, 400, 638 Auchettl K., et al., 2015, ApJ, 802, 68 BalletJ.,BruelP.,BurnettT.H.,LottB.,TheFermi-LATcollaboration2023, arXiv e-prints, p. arXiv:2307.12546 Bednarek W., Sitarek J., 2013, MNRAS, 430, 2951 Brownsberge...

work page doi:10.1142/10986 2017

[2] [2]

MNRAS000, 1–6 (2026)

2026