Bayesian Model Comparison and Analysis of the Galactic Disk Population of Gamma-Ray Millisecond Pulsars

Abstract

Pulsed emission from almost one hundred millisecond pulsars (MSPs) has been detected in $\gamma$-rays by the Fermi Large-Area Telescope. The global properties of this population remain relatively unconstrained despite many attempts to model their spatial and luminosity distributions. We perform here a self-consistent Bayesian analysis of both the spatial distribution and luminosity function simultaneously. Distance uncertainties, arising from errors in the parallax measurement or Galactic electron-density model, are marginalized over. We provide a public Python package for calculating distance uncertainties to pulsars derived using the dispersion measure by accounting for the uncertainties in Galactic electron-density model YMW16. Finally, we use multiple parameterizations for the MSP population and perform Bayesian model comparison, finding that a broken power law luminosity function with Lorimer spatial profile are preferred over multiple other parameterizations used in the past. The best-fit spatial distribution and number of $\gamma$-ray MSPs is consistent with results for the radio population of MSPs.

arXiv comment13 pages, 8 figures, 3 tables + Appendix. Public code and source list available from http://github.com/tedwards2412/MSPDist

References in corpus (55)

• 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
• X-ray spectral modelling of the AGN obscuring region in the CDFS: Bayesian model selection and catalogue
• The ATNF Pulsar Catalogue
• Fermi Large Area Telescope Second Source Catalog
• Bayes in the sky: Bayesian inference and model selection in cosmology
• A New Electron Density Model for Estimation of Pulsar and FRB Distances
• The Second Fermi Large Area Telescope Catalog of Gamma-ray Pulsars
• The Parkes multibeam pulsar survey: VI. Discovery and timing of 142 pulsars and a Galactic population analysis
• The Characterization of the Gamma-Ray Signal from the Central Milky Way: A Compelling Case for Annihilating Dark Matter
• Improved Estimates of the Milky Way's Stellar Mass and Star Formation Rate from Hierarchical Bayesian Meta-Analysis
• NE2001.I. A New Model for the Galactic Distribution of Free Electrons and its Fluctuations
• Background model systematics for the Fermi GeV excess
• Dark Matter and Pulsar Model Constraints from Galactic Center Fermi-LAT Gamma Ray Observations
• Strong support for the millisecond pulsar origin of the Galactic center GeV excess
• Evidence for Unresolved Gamma-Ray Point Sources in the Inner Galaxy
• Possible Evidence For Dark Matter Annihilation In The Inner Milky Way From The Fermi Gamma Ray Space Telescope
• On Pulsar Distance Measurements and their Uncertainties
• The Consistency of Fermi-LAT Observations of the Galactic Center with a Millisecond Pulsar Population in the Central Stellar Cluster
• Fermi Large Area Telescope Third Source Catalog
• Neutron Star Population Dynamics.I: Millisecond Pulsars
• Galactic Bulge Preferred Over Dark Matter for the Galactic Center Gamma-Ray Excess
• Millisecond pulsars and the Galactic Center gamma-ray excess: the importance of luminosity function and secondary emission
• Combining cosmological datasets: hyperparameters and Bayesian evidence
• The Fermi-LAT GeV Excess Traces Stellar Mass in the Galactic Bulge
• Classification and Ranking of Fermi LAT Gamma-ray Sources from the 3FGL Catalog using Machine Learning Techniques
• Binary Millisecond Pulsar Discovery via Gamma-Ray Pulsations
• Population synthesis of radio and gamma-ray millisecond pulsars from the Galactic disk
• A 1.4-GHz Arecibo Survey for Pulsars in Globular Clusters
• Challenges in Explaining the Galactic Center Gamma-Ray Excess with Millisecond Pulsars
• The Pulsar Contribution to the Gamma-Ray Background
• Gamma-ray Pulsar Revolution
• Discovery of two millisecond pulsars in Fermi sources with the Nancay Radio Telescope
• The High Time Resolution Universe Pulsar Survey VIII: The Galactic millisecond pulsar population
• Diffuse gamma-ray emission from galactic pulsars
• Radio detection prospects for a bulge population of millisecond pulsars as suggested by Fermi LAT observations of the inner Galaxy
• Binary and Millisecond Pulsars at the New Millennium
• A Spitzer Five-Band Analysis of the Jupiter-Sized Planet TrES-1
• Angular Momentum Exchange by Gravitational Torques and Infall in the Circumbinary Disk of the Protostellar System L1551 NE
• Active-Region Tilt Angles: Magnetic Versus White-Light Determinations of Joy's Law
• The Gamma-Ray Luminosity Function of Millisecond Pulsars and Implications for the GeV Excess
• The Gamma-Ray Pulsar Population of Globular Clusters: Implications for the GeV Excess
• Einstein@Home discovers a radio-quiet gamma-ray millisecond pulsar
• Pulsed Gamma-rays from the millisecond pulsar J0030+0451 with the Fermi Large Area Telescope
• Consistency Between the Luminosity Function of Resolved Millisecond Pulsars and the Galactic Center Excess
• Anisotropies in the gamma-ray sky from millisecond pulsars
• The Parkes multibeam pulsar survey: VII. Timing of four millisecond pulsars and the underlying spin period distribution of the Galactic millisecond pulsar population
• New Black Widows and Redbacks in the Galactic Field
• Millisecond pulsar origin of the Galactic center excess and extended gamma-ray emission from Andromeda - a closer look
• Estimating the GeV Emission of Millisecond Pulsars in Dwarf Spheroidal Galaxies
• Constraining the Galactic millisecond pulsar population using Fermi Large Area Telescope
• Fermi LAT pulsed detection of PSR J0737-3039A in the double pulsar system
• A stacked analysis of 115 pulsars observed by the Fermi LAT
• A New Determination of the Spectra and Luminosity Function of Gamma-Ray Millisecond Pulsars

Cited by in corpus (30)

• Dark Matter Strikes Back at the Galactic Center
• Strong Evidence that the Galactic Bulge is Shining in Gamma Rays
• Foreground Mismodeling and the Point Source Explanation of the Fermi Galactic Center Excess
• Pseudo-Goldstone dark matter confronts cosmic ray and collider anomalies
• The Enigmatic Galactic Center Excess: Spurious Point Sources and Signal Mismodeling
• Transient Radio Signatures from Neutron Star Encounters with QCD Axion Miniclusters
• Characterizing the Nature of the Unresolved Point Sources in the Galactic Center
• The GCE in a New Light: Disentangling the $\gamma$-ray Sky with Bayesian Graph Convolutional Neural Networks
• The discovery and scientific potential of fast radio bursts
• A neural simulation-based inference approach for characterizing the Galactic Center $\gamma$-ray excess
• Comparing the Galactic Bulge and Galactic Disk Millisecond Pulsars
• Maximum Entropy Estimation of the Galactic Bulge Morphology via the VVV Red Clump
• Luminosity functions consistent with a pulsar-dominated Galactic Center Excess
• Pevatron at the Galactic Center: Multi-Wavelength Signatures from Millisecond Pulsars
• Anti-Deuterons and Anti-Helium Nuclei from Annihilating Dark Matter
• The Galactic bulge millisecond pulsars shining in X rays: A gamma-ray perspective
• Millisecond Pulsars Modify the Radio-SFR Correlation in Quiescent Galaxies
• Inverse Compton emission from millisecond pulsars in the Galactic bulge
• Reconciling hints on axion-like-particles from high-energy gamma rays with stellar bounds
• Probing the Fermi-LAT GeV excess with gravitational waves
• Mind the gap: The discrepancy between simulation and reality drives interpretations of the Galactic Center Excess
• Comment on "Understanding the $\gamma$-ray emission from the globular cluster 47 Tuc: evidence for dark matter?"
• Harnessing the Population Statistics of Subhalos to Search for Annihilating Dark Matter
• A New Determination of the Millisecond Pulsar Gamma-Ray Luminosity Function and Implications for the Galactic Center Gamma-Ray Excess
• Towards resolving the Galactic center GeV excess with millisecond-pulsar-like sources using machine learning
• Role of Fan Beam Model in Population Synthesis of Isolated Radio Pulsars
• Applying Simulation-Based Inference to Spectral and Spatial Information from the Galactic Center Gamma-Ray Excess
• Characterizing the Expected Behavior of Non-Poissonian Template Fitting
• Constraining Gamma-ray Lines from Dark Matter Annihilation using Fermi-LAT and H.E.S.S. data
• Millisecond Pulsars in Globular Clusters and Implications for the Galactic Center Gamma-Ray Excess