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

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.

18 pages, 9 figures

References in corpus (39)

• The ATNF Pulsar Catalogue
• A statistical study of 233 pulsar proper motions
• A New Electron Density Model for Estimation of Pulsar and FRB Distances
• Time-dependent force-free pulsar magnetospheres: axisymmetric and oblique rotators
• Birth and Evolution of Isolated Radio Pulsars
• The Parkes Multibeam Pulsar Survey: I. Observing and Data Analysis Systems, Discovery and Timing of 100 Pulsars
• The Parkes multibeam pulsar survey: VI. Discovery and timing of 142 pulsars and a Galactic population analysis
• Unifying the observational diversity of isolated neutron stars via magneto-thermal evolution models
• The High Time Resolution Universe Pulsar Survey I: System configuration and initial discoveries
• Real-time gravitational-wave science with neural posterior estimation
• Spectral properties of 441 radio pulsars
• Radio Pulsar Death Line Revisited: Is PSR J2144-3933 Anomalous?
• Time evolution of pulsar obliquity angle from 3D simulations of magnetospheres
• On the rate of core collapse supernovae in the Milky Way
• Population synthesis studies of isolated neutron stars with magnetic field decay
• The Swinburne Intermediate Latitude Pulsar Survey
• Nuisance hardened data compression for fast likelihood-free inference
• PsrPopPy: An open-source package for pulsar population simulations
• Population synthesis of isolated Neutron Stars with magneto--rotational evolution
• Supernova Fallback as Origin of Neutron Star Spins and Spin-kick Alignment
• The observed velocity distribution of young pulsars II: analysis of complete PSR$\pi$
• The Thousand-Pulsar-Array programme on MeerKAT I: Science objectives and first results
• Population Synthesis of Isolated Neutron Stars with magneto-rotational evolution II: from radio-pulsars to magnetars
• Kicks and Induced Spins of Neutron Stars at Birth
• Neural posterior estimation for exoplanetary atmospheric retrieval
• The Thousand-Pulsar-Array program on MeerKAT -- IX. The time-averaged properties of the observed pulsar population
• Magneto-thermal evolution of neutron stars with coupled Ohmic, Hall and ambipolar effects via accurate finite-volume simulations
• A neural simulation-based inference approach for characterizing the Galactic Center $\gamma$-ray excess
• The Thousand-Pulsar-Array programme on MeerKAT -- VIII. The subpulse modulation of 1198 pulsars
• Markov chain Monte Carlo population synthesis of single radio pulsars in the Galaxy
• Discovery of 37 new pulsars through GPU-accelerated reprocessing of archival data of the Parkes Multibeam Pulsar Survey
• Population statistics study of radio and gamma-ray pulsars in the Galactic plane
• Bayesian model comparison for simulation-based inference
• Prospects for Discovering Pulsars in Future Continuum Surveys Using Variance Imaging
• Isolated pulsar population synthesis with simulation-based inference
• Simulation-based inference of black hole ringdowns in the time domain
• Simulation-Based Inference of the sky-averaged 21-cm signal from CD-EoR with REACH
• Back to the Starting Point: on the Simulation of Initial Magnetic Fields and Spin Periods of Non-accretion Pulsars
• Pulsar Population Synthesis with Magnetorotational Evolution: Constraining the Decay of Magnetic field