Agile Observations of the 'Soft' Gamma-Ray Pulsar PSR B1509-58

We present the results of new Agile observations of PSR B1509-58 performed over a period of 2.5 years following the detection obtained with a subset of the present data. The modulation significance of the lightcurve above 30 MeV is at a 5$\sigma$ confidence level and the lightcurve is similar to those found earlier by Comptel up to 30 MeV: a broad asymmetric first peak reaching its maximum 0.39 +/- 0.02 cycles after the radio peak plus a second peak at 0.94 +/- 0.03. The gamma-ray spectral energy distribution of the pulsed flux detected by Comptel and Agile is well described by a power-law (photon index alpha=1.87+/-0.09) with a remarkable cutoff at E_c=81 +/- 20 MeV, representing the softest spectrum observed among gamma-ray pulsars so far. The pulsar luminosity at E > 1 MeV is $L_{\gamma}=4.2^{+0.5}_{-0.2} \times10^{35}$ erg/s, assuming a distance of 5.2 kpc, which implies a spin-down conversion efficiency to gamma-rays of $\sim 0.03$. The unusual soft break in the spectrum of PSR B1509-58 has been interpreted in the framework of polar cap models as a signature of the exotic photon splitting process in the strong magnetic field of this pulsar. In this interpretation our spectrum constrains the magnetic altitude of the emission point(s) at 3 km above the neutron star surface, implying that the attenuation may not be as strong as formerly suggested because pair production can substitute photon splitting in regions of the magnetosphere where the magnetic field becomes too low to sustain photon splitting. In the case of an outer-gap scenario, or the two pole caustic model, better constraints on the geometry of the emission would be needed from the radio band in order to establish whether the conditions required by the models to reproduce Agile lightcurves and spectra match the polarization measurements.

💡 Research Summary

The paper reports on a comprehensive 2.5‑year monitoring campaign of the young, high‑magnetic‑field pulsar PSR B1509‑58 with the AGILE satellite, extending the earlier detection obtained from a subset of the data. By accumulating more than one million seconds of exposure, the authors achieve a 5σ detection of pulsed emission above 30 MeV. The gamma‑ray light curve shows two distinct peaks: the first, broad and asymmetric, lags the radio main pulse by 0.39 ± 0.02 rotational cycles, while the second appears near phase 0.94 ± 0.03. This morphology reproduces the pattern previously seen with COMPTEL up to 30 MeV, but AGILE’s superior timing resolution refines the peak positions and widths.

Spectrally, the authors combine the COMPTEL (1–30 MeV) and AGILE (30–100 MeV) datasets to construct a broadband spectral energy distribution (SED). The pulsed flux follows a simple power law with photon index α = 1.87 ± 0.09, but exhibits a remarkably low exponential cutoff at E_c = 81 ± 20 MeV. This cutoff is the softest ever measured for a gamma‑ray pulsar, placing PSR B1509‑58 in a unique class of “soft” gamma‑ray emitters. Integrating the SED above 1 MeV yields a gamma‑ray luminosity L_γ = 4.2 × 10³⁵ erg s⁻¹ (assuming a distance of 5.2 kpc), which corresponds to a conversion efficiency of roughly 3 % of the spin‑down power (Ė ≈ 1.3 × 10³⁸ erg s⁻¹). Although this efficiency is modest compared with many GeV pulsars, it is significant given the unusually soft spectrum.

To interpret the low‑energy break, the authors invoke photon splitting, a third‑order quantum electrodynamics (QED) process that can dominate in magnetic fields exceeding the critical value B_cr ≈ 4.4 × 10¹³ G. In the polar‑cap scenario, photon splitting attenuates high‑energy photons before they can produce electron‑positron pairs, leading to a steep spectral turnover. By fitting the observed cutoff, the authors infer that the emission region lies roughly 3 km above the neutron‑star surface, where the magnetic field is still strong enough (∼10¹³ G) to sustain efficient splitting. They also note that, as the field weakens with altitude, pair production may take over, providing an additional attenuation channel that can soften the spectrum without requiring photon splitting to be the sole mechanism.

The paper also discusses alternative outer‑gap and two‑pole caustic models. In these frameworks, the observed light‑curve morphology could be reproduced if the geometry (magnetic inclination α and observer viewing angle ζ) is appropriately chosen. However, current radio polarization data for PSR B1509‑58 are insufficient to tightly constrain these angles. The authors therefore call for high‑precision radio polarimetry and X‑ray timing to determine the geometry and to test whether the outer‑gap or caustic models can simultaneously account for the light‑curve shape and the soft spectral cutoff.

In summary, the AGILE observations provide the first high‑significance detection of PSR B1509‑58 above 30 MeV, confirm the double‑peaked light curve, and reveal an unprecedentedly low spectral cutoff. The results support the presence of photon splitting in the magnetosphere of a high‑field pulsar, constrain the emission altitude to a few kilometres above the surface, and highlight the need for multi‑wavelength polarization measurements to discriminate between polar‑cap and outer‑gap emission scenarios. This work thus advances our understanding of how extreme magnetic fields shape high‑energy radiation processes in young pulsars.