Beyond the Simple Power Law: A Bayesian Analysis of 897 Pulsar Spectra

Abstract

We present a comprehensive re-evaluation of pulsar radio spectra using the largest curated dataset of calibrated flux densities to date, comprising 897 pulsars, and employing a robust Bayesian framework for model comparison alongside frequentist methods. Contrary to the established consensus that pulsar spectra are predominantly simple power laws, our analysis reveals that complex spectral shapes with curvature or breaks are in fact the norm. The broken power law emerges as the most common spectral shape, accounting for 60\% of pulsars, while the simple power law describes only 13.5\%, with 68.8\% of pulsars decisively favoring curved or broken models. We further identify 74 confident gigahertz-peaked spectrum pulsars, and demonstrate that millisecond pulsars frequently exhibit spectral curvature. A key finding is that the previously reported dominance of the simple power law was largely a statistical artifact of the frequentist method used in earlier work. These findings substantially revise the prevailing view of pulsar spectra and establish a critical, model-classified foundation for future theoretical work.

Accepted for publication in Physical Review D; 13 pages, 8 figures, 7 tables

References in corpus (13)

• A Quick Look at the $3$GHz Radio Sky I. Source Statistics from the Very Large Array Sky Survey
• Polarimetry of 600 pulsars from observations at 1.4 GHz with the Parkes radio telescope
• The impact of glitches on young pulsar rotational evolution
• The MeerTime Pulsar Timing Array -- A Census of Emission Properties and Timing Potential
• The Thousand-Pulsar-Array program on MeerKAT -- IX. The time-averaged properties of the observed pulsar population
• Turn-over in pulsar spectra around 1 GHz
• The Galactic population and properties of young, highly-energetic pulsars
• Flux-density spectral analysis for several pulsars and two newly-identified gigahertz-peaked spectra
• pulsar_spectra: A pulsar flux density catalogue and spectrum fitting repository
• Discovery of a steep-spectrum low-luminosity pulsar with the Murchison Widefield Array
• Flux density monitoring of 89 millisecond pulsars with MeerKAT
• The birth mass function of neutron stars
• Radio spectra of pulsars fitted with the spectral distribution function of the emission from their current sheet