Constraining pulsar gap models with the light-curve and flux properties of the gamma-ray pulsar population

We compare population synthesis results for inner and outer magnetosphere emission models with the various characteristics measured in the first LAT pulsar catalogue for both the radio-loud and radio-weak or radio-quiet gamma-ray pulsars. We show that all models fail to reproduce the observations: for each model there is a lack of luminous and energetic objects that suggest a non dipolar magnetic field structure or spin-down evolution. The large dispersion that we find in the simulated gamma-ray luminosity versus spin-down power relation does not allow to use the present trend seen in the Fermi data to distinguish among models. For each model and each Fermi detected pulsar, we have generated light curves as a function of obliquity and inclination angles. The theoretical curves were fitted to the observed one, using a maximum-likelihood approach, to derive the best-fit orientations and to compare how well each model can reproduce the data. Including the radio light-curve gives an additional key constraint to restrict the orientation space

💡 Research Summary

The paper presents a comprehensive assessment of pulsar gap models by confronting population‑synthesis predictions with the observational properties of the first Fermi‑LAT pulsar catalogue (1PC). Four widely used magnetospheric emission scenarios are examined: the Polar‑Cap (PC) model, the Slot‑Gap (SG) model, the Outer‑Gap (OG) model, and the One‑Pole Caustic (OPC) model, which is a variant of the OG geometry. For each model the authors generate synthetic pulsar populations by drawing initial spin periods, magnetic field strengths, magnetic‑axis inclination angles (α) and observer viewing angles (ζ) from prescribed distributions. The spin‑down power (Ė) of each synthetic pulsar is evolved according to the standard dipole braking law, and the corresponding γ‑ray luminosity (Lγ) is calculated using the model‑specific radiation efficiency prescriptions.

The simulated populations are then compared with the measured distributions of Ė, Lγ, and the ratio of radio‑loud to radio‑quiet γ‑ray pulsars in the 1PC sample. Three principal discrepancies emerge. First, none of the models can reproduce the high‑luminosity, high‑Ė tail of the observed distribution. The synthetic pulsars with Ė > 10^36 erg s⁻¹ systematically under‑predict Lγ, yielding a deficit of objects with Lγ ≈ 10^35 erg s⁻¹ that are abundant in the LAT data. The PC and OG models show the steepest decline in luminosity at high Ė, while the SG model fares slightly better in the intermediate Ė regime but still fails at the extreme end.

Second, the simulated Lγ–Ė relation exhibits a large intrinsic scatter, reflecting the strong dependence of the predicted γ‑ray flux on the geometry (α, ζ) and on model‑specific gap width parameters. Even after averaging over the full angular distribution and applying observational selection effects (e.g., LAT sensitivity, radio‑beam visibility), the synthetic Lγ–Ė trend remains much broader than the relatively tight correlation seen in the LAT data. This suggests that the current gap prescriptions do not capture the true scaling of γ‑ray efficiency with spin‑down power.

Third, the authors construct model light curves for each synthetic pulsar over a dense grid of (α, ζ) values and fit them to the observed γ‑ray pulse profiles using a maximum‑likelihood approach. When the radio pulse profile is added as an additional constraint, the allowed region in (α, ζ) space shrinks dramatically and differs markedly among the models. The OG and OPC geometries permit only a narrow set of orientations that produce both a radio and a γ‑ray beam intersecting the line of sight, leading to an under‑prediction of the observed radio‑loud fraction. The SG geometry allows a broader range of orientations, improving the match to the radio‑loud/quiet statistics, but it does not resolve the high‑Ė luminosity shortfall.

Collectively, these findings indicate that the four canonical gap models, each built on the assumption of a pure dipolar magnetic field and a simple spin‑down law, are insufficient to describe the full γ‑ray pulsar population. The inability to generate enough luminous, energetic pulsars points to missing physics, such as non‑dipolar (multipolar) magnetic field components that can enhance particle acceleration, or more complex spin‑down evolution (e.g., magnetic field decay, torque variations) that modifies the Ė distribution over time. Moreover, the large dispersion in the simulated Lγ–Ė relation undermines the use of the observed trend as a discriminant among models.

The paper concludes by outlining future directions: incorporating multipolar field configurations into gap electrodynamics, allowing the γ‑ray efficiency to evolve with Ė in a physically motivated way, and refining the treatment of observational biases (both γ‑ray and radio) in population synthesis. By addressing these aspects, subsequent studies may bridge the gap between theory and the rich dataset now provided by Fermi‑LAT, ultimately leading to a more accurate physical picture of particle acceleration and high‑energy emission in rotation‑powered pulsars.