Radio pulsar population synthesis with consistent flux measurements using simulation-based inference

TL;DR

This paper introduces a simulation-based inference method using TSNPE to efficiently infer neutron star properties from pulsar population data, incorporating new flux measurements from the TPA program to improve luminosity estimates.

Contribution

It applies TSNPE to pulsar population synthesis, demonstrating improved efficiency and incorporating new flux data to better constrain neutron star luminosity.

Findings

TSNPE outperforms standard NPE in efficiency.

Inclusion of flux data enhances luminosity parameter constraints.

Results are consistent with previous studies.

Abstract

The properties of the entire neutron star population can be inferred by modeling their evolution, from birth to the present, through pulsar population synthesis. This involves simulating a mock population, applying observational filters, and comparing the resulting sources to the limited subset of detected pulsars. We specifically focus on the magneto-rotational properties of Galactic isolated neutron stars and provide new insights into the intrinsic radio luminosity law by combining pulsar population synthesis with a simulation-based inference (SBI) technique called truncated sequential neural posterior estimation (TSNPE). We employ TSNPE to train a neural density estimator on simulated pulsar populations to approximate the posterior distribution of the underlying parameters. This technique efficiently explores the parameter space by concentrating on regions that are most likely to match the observed data thus allowing a significant reduction in training dataset size. We demonstrate the efficiency of TSNPE over standard neural posterior estimation (NPE), achieving robust inferences of magneto-rotational parameters consistent with previous studies using only around 4% of the simulations required by NPE approaches. Moreover, for the first time, we incorporate data from the Thousand Pulsar Array (TPA) program on MeerKAT, the largest unified sample of neutron stars with consistent fluxes measurement to date, to help constrain the stars' intrinsic radio luminosity. We find that adding flux information as an input to the neural network largely improves the constraints on the pulsars' radio luminosity, as well as improving the estimates on other input parameters.