Tomography of the Galactic free electron density with the Square Kilometer Array

Abstract

We present a new algorithm to reconstruct the Galactic free electron density from pulsar dispersion measures. The algorithm performs a nonparametric tomography for a density field with an arbitrary amount of degrees of freedom. It is based on approximating the Galactic free electron density as the product of a profile function with a statistically isotropic and homogeneous log-normal field. Under this approximation the algorithm generates a map of the free electron density as well as an uncertainty estimate without the need of information about the power spectrum. The uncertainties of the pulsar distances are treated consistently by an iterative procedure. We test the algorithm using the NE2001 model with modified fluctuations as a Galaxy model, pulsar populations generated from the Lorimer population model, and mock observations emulating the upcoming Square Kilometer Array. We show the quality of the reconstruction for mock data sets containing between 1000 and 10000 pulsars with distance uncertainties up to 25%. Our results show, that with the SKA nonparametric tomography of the Galactic free electron density becomes feasible, but the quality of the reconstruction is very sensitive to the distance uncertainties.

18 pages, 16 figures

References in corpus (8)

• The Parkes multibeam pulsar survey: VI. Discovery and timing of 142 pulsars and a Galactic population analysis
• Radio observational constraints on Galactic 3D-emission models
• The Vertical Structure of Warm Ionised Gas in the Milky Way
• PsrPopPy: An open-source package for pulsar population simulations
• Reconstruction of Gaussian and log-normal fields with spectral smoothness
• Reconstruction of signals with unknown spectra in information field theory with parameter uncertainty
• RESOLVE: A new algorithm for aperture synthesis imaging of extended emission in radio astronomy
• Three-dimensional extinction mapping using Gaussian random fields

Cited by in corpus (7)

• IMAGINE: A comprehensive view of the interstellar medium, Galactic magnetic fields and cosmic rays
• The Galactic Faraday depth sky revisited
• Bayesian weak lensing tomography: Reconstructing the 3D large-scale distribution of matter with a lognormal prior
• Disentangling the Faraday rotation sky
• The Galaxy in circular polarization: all-sky radio prediction, detection strategy, and the charge of the leptonic cosmic rays
• A new method of reconstructing Galactic three-dimensional structures using ultralong-wavelength radio observations
• D2O - a distributed data object for parallel high-performance computing in Python