Interpolation techniques for reconstructing Galactic Faraday rotation

Abstract

The line-of-sight structure of the Galactic magnetic field (GMF) can be studied using Faraday rotation measure (RM) grids. We analyze how the choice of interpolation kernel can affect the accuracy and reliability of reconstructed RM maps. We test the following kernels: inverse distance weighting (IDW), natural neighbour interpolation (NNI), inverse multiquadric interpolation (IM), thin-plate spline interpolation (TPS), and a Bayesian rotation measure sky (BRMS); all techniques were tested on two simulated Galactic foreground RMs (one assuming the GMF has patchy structures and the other assuming it has filamentary structures) using magnetohydrodynamic simulations. Both foregrounds were sampled to form RM grids with densities of $\sim$40 sources deg${}^{-2}$ and area $\sim$144 deg${}^2$. The techniques were tested on data sets with different noise levels and Gaussian random extragalactic RM contributions. The data set that most closely emulates expected data from current surveys, such as the POlarization Sky Survey of the Universe's Magnetism (POSSUM), had extragalactic contributions and a noise standard deviation of $\sim 1.5$ rad m${}^{-2}$. For this data set, the accuracy of the techniques for the patchy structures from best to worst was: BRMS, NNI, TPS, IDW and IM; while in the filamentary simulate foreground it was: BRMS, NNI, TPS, and IDW. IDW is the most computationally expensive technique, while TPS and IM are the least expensive. BRMS and NNI have the same, intermediate computational cost. This analysis lays the groundwork for Galactic RM studies with large radio polarization sky surveys, such as POSSUM.

32 pages, 22 figures, 8 tables

References in corpus (33)

• SciPy 1.0--Fundamental Algorithms for Scientific Computing in Python
• Array Programming with NumPy
• The Interstellar Environment of our Galaxy
• Comparing the statistics of interstellar turbulence in simulations and observations: Solenoidal versus compressive turbulence forcing
• Faraday Rotation Measure Synthesis
• A New Model of the Galactic Magnetic Field
• An improved map of the Galactic Faraday sky
• Mach Number Dependence of Turbulent Magnetic Field Amplification: Solenoidal versus Compressive Flows
• A Bayesian view on Faraday rotation maps - Seeing the magnetic power spectra in galaxy clusters
• Estimating extragalactic Faraday rotation
• Complex Faraday depth structure of Active Galactic Nuclei as revealed by broadband radio polarimetry
• Integrated Polarization of Sources at lambda ~1m and New Rotation Measure Ambiguities
• Structures in the RM sky
• Variable structures in M87* from space, time and frequency resolved interferometry
• The Coherent Magnetic Field of the Milky Way
• The Galactic magnetic field as spectrograph for ultra-high energy cosmic rays
• Magnetic Fields and Ionized Gas in the inner Galaxy: An Outer Scale for Turbulence and the Possible Role of HII Regions
• Plasma turbulence in the interstellar medium
• Cosmic Rays in Intermittent Magnetic Fields
• The Galactic Faraday depth sky revisited
• The correct sense of Faraday rotation
• The Faraday Rotation Measure Grid of the LOFAR Two-metre Sky Survey: Data Release 2
• Information theory for fields
• Geometric variational inference
• Magnetic fields in the Milky Way from pulsar observations: effect of the correlation between thermal electrons and magnetic fields
• A physical approach to modelling large-scale galactic magnetic fields
• Rotation measure structure functions with higher-order stencils as a probe of small-scale magnetic fluctuations and its application to the Small and Large Magellanic Clouds
• RMTable2023 and PolSpectra2023: standards for reporting polarization and Faraday rotation measurements of radio sources
• A Faraday Rotation Template for the Galactic Sky
• Re-Envisioning Numerical Information Field Theory (NIFTy.re): A Library for Gaussian Processes and Variational Inference
• Disentangling the Faraday rotation sky
• Alignment of dense molecular core morphology and velocity gradients with ambient magnetic fields
• A method for reconstructing the Galactic magnetic field using dispersion of fast radio bursts and Faraday rotation of radio galaxies