Bayesian inference methodology to characterize the dust emissivity at far-infrared and submillimeter frequencies

Abstract

We present a Bayesian inference method to characterise the dust emission properties using the well-known dust-HI correlation in the diffuse interstellar medium at Planck frequencies $\nu \ge 217$ GHz. We use the Galactic HI map from the Galactic All-Sky Survey (GASS) as a template to trace the Galactic dust emission. We jointly infer the pixel-dependent dust emissivity and the zero level present in the Planck intensity maps. We use the Hamiltonian Monte Carlo technique to sample the high dimensional parameter space ($D\sim 10^3$). We demonstrate that the methodology leads to unbiased recovery of dust emissivity per pixel and the zero level when applied to realistic Planck sky simulations over a 6300 deg${}^2$ area around the Southern Galactic pole. As an application on data, we analyse the Planck intensity map at 353 GHz to jointly infer the pixel-dependent dust emissivity at Nside=32 resolution (1.8° pixel size) and the global offset. We find that the spatially varying dust emissivity has a mean of 0.031 MJysr${}^{-1}(10^{20}\mathrm{c}{\mathrm{m}}^{-2}{)}^{-1}$ and $1\sigma$ standard deviation of 0.007 MJysr${}^{-1}(10^{20}\mathrm{c}{\mathrm{m}}^{-2}{)}^{-1}$. The mean dust emissivity increases monotonically with increasing mean HI column density. We find that the inferred global offset is consistent with the expected level of Cosmic Infrared Background (CIB) monopole added to the Planck data at 353 GHz. This method is useful in studying the line-of-sight variations of dust spectral energy distribution in the multi-phase interstellar medium.

19 pages, 14 figures, Accepted for publication in MNRAS

References in corpus (48)

• HEALPix -- a Framework for High Resolution Discretization, and Fast Analysis of Data Distributed on the Sphere
• The Leiden/Argentine/Bonn (LAB) Survey of Galactic HI: Final data release of the combined LDS and IAR surveys with improved stray-radiation corrections
• MCMC using Hamiltonian dynamics
• Rank-normalization, folding, and localization: An improved $\hat{R}$ for assessing convergence of MCMC
• Planck 2018 results. I. Overview and the cosmological legacy of Planck
• HI4PI: A full-sky HI survey based on EBHIS and GASS
• Planck 2013 results. XI. All-sky model of thermal dust emission
• Planck 2018 results. IV. Diffuse component separation
• Planck 2018 results. III. High Frequency Instrument data processing and frequency maps
• GASS: The Parkes Galactic All-Sky Survey. I. Survey Description, Goals, and Initial Data Release
• Planck intermediate results. LVII. Joint Planck LFI and HFI data processing
• Planck 2013 results. XXX. Cosmic infrared background measurements and implications for star formation
• Dusty Infrared Galaxies: Sources of the Cosmic Infrared Background
• The Origin of the Magellanic Stream and Its Leading Arm
• Planck 2015 results. VIII. High Frequency Instrument data processing: Calibration and maps
• Planck 2013 results. XXVII. Doppler boosting of the CMB: Eppur si muove
• A unified empirical model for infrared galaxy counts based on observed physical evolution of distant galaxies
• Global properties of the HI distribution in the outer Milky Way
• Planck Early Results XVIII: The power spectrum of cosmic infrared background anisotropies
• Planck Early Results XXIV: Dust in the diffuse interstellar medium and the Galactic halo
• Planck 2018 results. XI. Polarized dust foregrounds
• GASS: The Parkes Galactic All-Sky Survey. II. Stray-Radiation Correction and Second Data Release
• Planck intermediate results. XXII. Frequency dependence of thermal emission from Galactic dust in intensity and polarization
• Planck 2013 results. IX. HFI spectral response
• The Magellanic Stream: Circumnavigating the Galaxy
• The 200 Degree-Long Magellanic Stream System
• Planck intermediate results. XVII. Emission of dust in the diffuse interstellar medium from the far-infrared to microwave frequencies
• Planck 2013 results. XVIII. Gravitational lensing-infrared background correlation
• Planck 2013 results. VIII. HFI photometric calibration and mapmaking
• Planck 2015 results. XXIII. The thermal Sunyaev-Zeldovich effect--cosmic infrared background correlation
• Bayesian power-spectrum inference for Large Scale Structure data
• Randomized Hamiltonian Monte Carlo
• The Atacama Cosmology Telescope: Lensing of CMB Temperature and Polarization Derived from Cosmic Infrared Background Cross-Correlation
• Fast optimal CMB power spectrum estimation with Hamiltonian sampling
• Demonstration of cosmic microwave background delensing using the cosmic infrared background
• Efficient Cosmological Parameter Estimation with Hamiltonian Monte Carlo
• Methods for Bayesian power spectrum inference with galaxy surveys
• Large-scale Maps of the Cosmic Infrared Background from Planck
• Modelling and simulation of large-scale polarized dust emission over the southern Galactic cap using the GASS HI data
• Impact of polarised extragalactic sources on the measurement of CMB B-mode anisotropies
• Bayesian inference of CMB gravitational lensing
• The four leading arms of the Magellanic Cloud system
• Determining thermal dust emission from Planck HFI data using a sparse, parametric technique
• Constraining cosmology with the cosmic microwave and infrared backgrounds correlation
• Dust polarization modeling at large-scale over the Northern Galactic cap using EBHIS and Planck data
• Hierarchical Bayesian CMB Component Separation with the No-U-Turn Sampler
• Constraining the physics of star formation from CIB-cosmic shear cross-correlations
• Polarization of the Cosmic Infrared Background Fluctuations