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_{\rm 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.08}_{-0.10}$, $\sigma_{\log B} = 0.45^{+0.05}_{-0.05}$ and $\mu_{\log P} = -1.00^{+0.26}_{-0.21}$, $\sigma_{\log P} = 0.38^{+0.33}_{-0.18}$ for the lognormal distributions and $a_{\rm late} = -1.80^{+0.65}_{-0.61}$ for the power law at the $95\%$ credible interval. We contrast our results with previous studies and highlight uncertainties of the inferred $a_{\rm 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.

What problem does this paper attempt to address?

The problem that this paper attempts to solve is to constrain the magneto - rotational characteristics of isolated Galactic radio pulsars by combining pulsar population synthesis with Simulation - Based Inference (SBI). Specifically, the authors developed a framework for modeling the initial properties of neutron stars, as well as their dynamical and magneto - rotational evolution, and used the SBI method to infer the posterior distributions of these properties. This helps to understand the formation and evolution processes of pulsars more accurately. The following are the main problems and solutions in the paper: ### 1. **Research Background** - **Problem**: How to better understand the magneto - rotational characteristics of isolated Galactic radio pulsars? - **Challenge**: Traditional Bayesian Markov Chain Monte Carlo (MCMC) or nested sampling methods become infeasible when dealing with complex pulsar population models because explicit likelihood functions cannot be written. - **Solution**: Introduce Simulation - Based Inference (SBI) and use deep neural networks to infer the posterior distributions of parameters. ### 2. **Research Objectives** - **Objective**: Through the SBI method, infer the parameters required to describe the initial magnetic field of pulsars (such as the mean and standard deviation of the log - normal distribution) and the power - law index of late - time magnetic decay. - **Specific Parameters**: - Log - normal distribution parameters of the initial magnetic field strength $ B $ and rotation period $ P $: - Mean: $ \mu_{\log B} = 13.10^{+0.08}_{-0.10} $, $ \mu_{\log P} = - 1.00^{+0.26}_{-0.21} $ - Standard deviation: $ \sigma_{\log B} = 0.45^{+0.05}_{-0.05} $, $ \sigma_{\log P} = 0.38^{+0.33}_{-0.18} $ - Power - law index of late - time magnetic decay: $ a_{\text{late}} = - 1.80^{+0.65}_{-0.61} $ ### 3. **Research Methods** - **Steps**: 1. **Modeling**: Establish a framework for modeling the initial properties of neutron stars (including the initial magnetic field strength $ B $ and rotation period $ P $) and their dynamical and magneto - rotational evolution. 2. **Simulation**: Generate a large number of synthetic $ P-\dot{P} $ diagrams, where $ P $ is the rotation period of the pulsar and $ \dot{P} $ is its time derivative. 3. **Training and Validation**: Use deep neural networks for training and verify the effectiveness of these networks through simulated data. 4. **Application**: Apply the optimized deep - learning framework to the actually observed pulsar population to infer the posterior distributions of parameters. ### 4. **Results and Significance** - **Results**: Through the SBI method, the parameters required for the initial magnetic field of pulsars and the power - law index of their late - time magnetic decay were successfully inferred. - **Significance**: This research provides an important step for robust statistical inference in complex pulsar population synthesis frameworks and lays the foundation for further multi - band analysis. ### Summary This paper successfully solves the problem of the magneto - rotational characteristics of isolated Galactic radio pulsars by combining pulsar population synthesis with Simulation - Based Inference, providing a new perspective and tool for understanding the formation and evolution of pulsars.

Isolated pulsar population synthesis with simulation-based inference

Isolated Pulsar Population Synthesis with Simulation-based Inference

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

Population Synthesis of Isolated Neutron Stars with magneto-rotational evolution II: from radio-pulsars to magnetars

Markov Chain Monte Carlo population synthesis of single radio pulsars in the Galaxy

Population synthesis of Galactic pulsars with machine learning

The Galactic population of canonical pulsars II

Population Synthesis of Radio Pulsars in the Fermi Era

Pulsar Population Synthesis with Magnetorotational Evolution: Constraining the Decay of Magnetic field

Pulsar Population Synthesis with Magnetorotational Evolution: Constraining the Decay of the Magnetic Field

40 Years of Pulsars: The Birth and Evolution of Isolated Radio Pulsars

Population Synthesis of Radio and Gamma-ray Pulsars in the Fermi Era

Gaussian mixture model and population synthesis of radio pulsars

Binary neutron star populations in the Milky Way

Binary population synthesis of the Galactic canonical pulsar population

Populating the Galaxy with pulsars I: stellar & binary evolution

Furnishing the Galaxy with Pulsars

A NICER View of the Nearest and Brightest Millisecond Pulsar: PSR J0437$\unicode{x2013}$4715

Multi-wavelength pulse profiles from the force-free neutron star magnetosphere

Progress on Radio Pulsar Survey and Estimate the Neighbor Pulsar Population

The radio pulsar population of the Small Magellanic Cloud