Isolated pulsar population synthesis with simulation-based inference

Abstract

We combine pulsar population synthesis with simulation-based inference (SBI) to constrain the magnetorotational properties of isolated Galactic radio pulsars. We first develop a framework to model neutron star birth properties and their dynamical and magnetorotational evolution. We specifically sample initial magnetic field strengths, $B$, and spin periods, $P$, from lognormal distributions and capture the late-time magnetic field decay with a power law. Each lognormal is described by a mean, ${\mu }_{\log B},{\mu }_{\log P}$, and standard deviation, ${\sigma }_{\log B},{\sigma }_{\log P}$, while the power law is characterized by the index, $a_{\mathrm{late}}$. We subsequently model the stars' radio emission and observational biases to mimic detections with three radio surveys, and we produce a large database of synthetic $P$--$\dot{P}$ diagrams by varying our five magnetorotational input parameters. We then follow an SBI approach that focuses on neural posterior estimation and train deep neural networks to infer the parameters' posterior distributions. After successfully validating these individual neural density estimators on simulated data, we use an ensemble of networks to infer the posterior distributions for the observed pulsar population. We obtain ${\mu }_{\log B}=13.10_{-0.10}^{+0.08}$, ${\sigma }_{\log B}=0.45_{-0.05}^{+0.05}$ and ${\mu }_{\log P}=-1.00_{-0.21}^{+0.26}$, ${\sigma }_{\log P}=0.38_{-0.18}^{+0.33}$ for the lognormal distributions and $a_{\mathrm{late}}=-1.80_{-0.61}^{+0.65}$ for the power law at the $95\%$ credible interval. We contrast our results with previous studies and highlight uncertainties of the inferred $a_{\mathrm{late}}$ value. Our approach represents a crucial step toward robust statistical inference for complex population synthesis frameworks and forms the basis for future multiwavelength analyses of Galactic pulsars.

arXiv comment31 pages, 16 figures, 5 tables, 2 appendices; published version

References in corpus (80)

• SciPy 1.0--Fundamental Algorithms for Scientific Computing in Python
• Array Programming with NumPy
• Astropy: A Community Python Package for Astronomy
• emcee: The MCMC Hammer
• The NumPy array: a structure for efficient numerical computation
• GW170817: Observation of Gravitational Waves from a Binary Neutron Star Inspiral
• The Structure of Cold Dark Matter Halos
• The Astropy Project: Building an inclusive, open-science project and status of the v2.0 core package
• HEALPix -- a Framework for High Resolution Discretization, and Fast Analysis of Data Distributed on the Sphere
• The ATNF Pulsar Catalogue
• MultiNest: an efficient and robust Bayesian inference tool for cosmology and particle physics
• dynesty: A Dynamic Nested Sampling Package for Estimating Bayesian Posteriors and Evidences
• galpy: A Python Library for Galactic Dynamics
• A statistical study of 233 pulsar proper motions
• A unified equation of state of dense matter and neutron star structure
• Bilby: A user-friendly Bayesian inference library for gravitational-wave astronomy
• Short-Duration Gamma-Ray Bursts
• A New Electron Density Model for Estimation of Pulsar and FRB Distances
• Time-dependent force-free pulsar magnetospheres: axisymmetric and oblique rotators
• An Ultraluminous X-ray Source Powered by An Accreting Neutron Star
• Birth and Evolution of Isolated Radio Pulsars
• Physics of Neutron Star Crusts
• The Parkes Multibeam Pulsar Survey: I. Observing and Data Analysis Systems, Discovery and Timing of 100 Pulsars
• Binary and Millisecond 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
• Searches for Fast Radio Transients
• Neutron stars - cooling and transport
• Markov Chain Monte Carlo Methods for Bayesian Data Analysis in Astronomy
• Real-time gravitational-wave science with neural posterior estimation
• A time-resolved picture of our Milky Way's early formation history
• Spectral properties of 441 radio pulsars
• Fast likelihood-free cosmology with neural density estimators and active learning
• The observed velocity distribution of young pulsars
• On the birthrates of Galactic neutron stars
• Radio Pulsar Death Line Revisited: Is PSR J2144-3933 Anomalous?
• A three-dimensional map of the Milky Way using classical Cepheid variable stars
• 2D Cooling of Magnetized Neutron Stars
• A highly resistive layer within the crust of X-ray pulsars limits their spin periods
• Revisiting the radial distribution of pulsars in the Galaxy
• Understanding the radio emission geometry of PSR B0329+54
• 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
• Ionization Break-Out from Millisecond Pulsar Wind Nebulae: an X-ray Probe of the Origin of Superluminous Supernovae
• The Swinburne Intermediate Latitude Pulsar Survey
• Tests of the nuclear equation of state and superfluid and superconducting gaps using the Cassiopeia A neutron star
• PsrPopPy: An open-source package for pulsar population simulations
• Magnetic, thermal and rotational evolution of isolated neutron stars
• Population synthesis of isolated Neutron Stars with magneto--rotational evolution
• Gaia DR2 in 6D: Searching for the fastest stars in the Galaxy
• Unveiling the Engines of Fast Radio Bursts, Super-Luminous Supernovae, and Gamma-Ray Bursts
• Magnetic, thermal and rotational evolution of isolated neutron stars
• Scatter broadening measurements of 124 pulsars at 327 MHz
• 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$
• A guided map to the Spiral arms in the galactic disk of the Milky Way
• Population Synthesis of Isolated Neutron Stars with magneto-rotational evolution II: from radio-pulsars to magnetars
• The Sun was not born in M 67
• Kicks and Induced Spins of Neutron Stars at Birth
• On the pulse--width statistics in radio pulsars. I. Importance of the interpulse emission
• 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
• On the pulse-width statistics in radio pulsars. II. Importance of the core profile components
• The Galactic population and properties of young, highly-energetic pulsars
• Robust Simulation-Based Inference in Cosmology with Bayesian Neural Networks
• 3D code for MAgneto-Thermal evolution in Isolated Neutron Stars, MATINS: The Magnetic Field Formalism
• 3D Magneto-thermal Simulations of Tangled Crustal Magnetic Field in Central Compact Objects
• X-ray Emission from Isolated Neutron Stars revisited: 3D magnetothermal simulations
• Meterwavelength Single-pulse Polarimetric Emission Survey IV: The Period Dependence of Component Widths of Pulsars
• Markov chain Monte Carlo population synthesis of single radio pulsars in the Galaxy
• Population statistics study of radio and gamma-ray pulsars in the Galactic plane
• Accurate X-ray Timing in the Presence of Systematic Biases With Simulation-Based Inference
• A simulation-based inference pipeline for cosmic shear with the Kilo-Degree Survey
• Analyzing the Galactic pulsar distribution with machine learning
• Testing the robustness of simulation-based gravitational-wave population inference
• Back to the Starting Point: on the Simulation of Initial Magnetic Fields and Spin Periods of Non-accretion Pulsars

Cited by in corpus (4)

• The Galactic population of magnetars : a simulation-based inference study
• Radio pulsar population synthesis with consistent flux measurements using simulation-based inference
• Pulsar Population Synthesis with Magnetorotational Evolution: Constraining the Decay of Magnetic field
• Pulsars identified in the LOFAR Two-metre Sky Survey at 144 MHz