Future Space-based Gamma-ray Pulsar Timing Arrays

Radio pulsar timing array (PTA) experiments using millisecond pulsars (MSPs) are beginning to detect nHz gravitational waves (GWs). MSPs are bright GeV gamma-ray emitters, and all-sky monitoring of about 100 MSPs with the Fermi Large Area Telescope (LAT) has enabled a gamma-ray Pulsar Timing Array. The GPTA provides a complementary view of nHz GWs because its MSP sample is different, and because the gamma-ray data are immune to plasma propagation effects, have minimal data gaps, and rely on homogeneous instrumentation. To assess GPTA performance for future gamma-ray observatories, we simulated the population of Galactic MSPs and developed a high-fidelity method to predict their gamma-ray spectra. This combination reproduces the properties of the LAT MSP sample, validating it for future population studies. We determined the expected signal from the simulated gamma-ray MSPs for instrument concepts with a wide range of capabilities. We found that the optimal GPTA energy range runs about 0.1 to 5 GeV, but we also examined Compton/MeV instruments. With the caveat that the MSP spectra models are extrapolated beyond observational constraints, we found low signal-to-background ratios, yielding few MSP detections. GeV-band concepts would detect 10$^3$ to 10$^4$ MSPs and achieve GW sensitivity on par with and surpassing the current generation of radio PTAs, reaching the GW self-noise regime. When considering two possible scenarios for the formation of MSPs in the Galactic bulge, the collective signal from which is a potential source of an excess GeV signal observed towards the Galactic center, we find that most of the concepts can both detect this bulge population and distinguish the production channel. In summary, the high discovery potential, strong GW performance, and tremendous synergy with radio PTAs all argue for the pursuit of next-generation gamma-ray pulsar timing.

💡 Research Summary

This paper evaluates the potential of future space‑based gamma‑ray observatories to build a Gamma‑ray Pulsar Timing Array (GPTA) that can rival or surpass current radio pulsar timing arrays (PTAs) in detecting nanohertz gravitational waves (GWs). The authors begin by synthesizing realistic Galactic millisecond pulsar (MSP) populations using the PSRPOPPy framework. Three scenarios are considered: (S1) a pure disk MSP population, and (S2, S3) disk plus a bulge component. The bulge component is motivated by the Galactic Center Excess (GCE) and is modeled under two formation channels – accretion‑induced collapse (AIC) of white dwarfs (S2) and a population with the same spin‑down properties as the disk (S3). Spatially, the bulge follows a r⁻² law consistent with the GCE morphology.

For each simulated MSP, the authors attach a gamma‑ray spectral model that scales with spin‑down power (Ė) and magnetic field, using a power‑law with an exponential cutoff (Γ≈1.5, E_cut≈3 GeV). This prescription reproduces the observed Fermi‑LAT MSP sample, validating the approach for extrapolation to fainter, as‑yet‑undetected pulsars.

The next step is to assess a suite of hypothetical gamma‑ray instruments. Key parameters include effective area, energy range, angular resolution, and background rejection. The authors find that the optimal energy band for timing is roughly 0.1–5 GeV; the lower bound is set by the need for good point‑source localization, while the upper bound captures the bulk of MSP emission before the exponential cutoff. Instruments operating in the MeV (Compton) regime suffer from very low signal‑to‑background ratios because MSP spectra fall off sharply and the diffuse Galactic background dominates.

In the optimal GeV band, a next‑generation instrument with an effective area of a few square meters and sub‑0.1° angular resolution would detect between 1,000 and 10,000 MSPs, a factor of ten to a hundred more than the current LAT sample. Many of these new detections would be radio‑quiet or too faint for existing radio surveys, thereby expanding the PTA sky coverage. The timing precision achievable from gamma‑ray photon arrival times is estimated at ~0.5 µs for the brightest sources, comparable to the best radio MSPs.

Using the Hellings‑Downs correlation framework, the authors translate the enlarged MSP set into a GW sensitivity forecast. They show that such a GPTA could detect a stochastic GW background with characteristic strain amplitude A_gwb ≈ 2 × 10⁻¹⁵ (the level currently hinted at by radio PTAs) within roughly five years of observation, reaching the “self‑noise” regime where the array’s own timing noise limits sensitivity. Inclusion of a bulge MSP population (S2 or S3) further enhances sensitivity to GW signals originating near the Galactic Center, and the GPTA could distinguish between the two bulge formation scenarios based on the spatial distribution of the detected pulsars.

The paper also highlights the complementary nature of gamma‑ray and radio PTAs. Gamma‑ray timing is immune to interstellar dispersion measure variations, suffers no data gaps, and uses a homogeneous instrument, eliminating many systematic uncertainties that plague radio PTAs (e.g., ionospheric effects, RFI, heterogeneous telescopes). Consequently, a combined analysis could improve overall GW detection significance by ~20 % and provide cross‑validation of individual pulsar timing solutions.

In summary, the study demonstrates that a future space‑based gamma‑ray mission optimized for the 0.1–5 GeV band could discover thousands of new MSPs, achieve timing precision on par with the best radio pulsars, and deliver GW sensitivity that matches or exceeds the current generation of radio PTAs. This would not only accelerate the detection of the nanohertz GW background but also shed light on the origin of the Galactic Center Excess and deepen our understanding of MSP formation and evolution.