skip to main content


Title: Cleaning our own dust: simulating and separating galactic dust foregrounds with neural networks
ABSTRACT Separating galactic foreground emission from maps of the cosmic microwave background (CMB) and quantifying the uncertainty in the CMB maps due to errors in foreground separation are important for avoiding biases in scientific conclusions. Our ability to quantify such uncertainty is limited by our lack of a model for the statistical distribution of the foreground emission. Here, we use a deep convolutional generative adversarial network (DCGAN) to create an effective non-Gaussian statistical model for intensity of emission by interstellar dust. For training data we use a set of dust maps inferred from observations by the Planck satellite. A DCGAN is uniquely suited for such unsupervised learning tasks as it can learn to model a complex non-Gaussian distribution directly from examples. We then use these simulations to train a second neural network to estimate the underlying CMB signal from dust-contaminated maps. We discuss other potential uses for the trained DCGAN, and the generalization to polarized emission from both dust and synchrotron.  more » « less
Award ID(s):
1812199 1560482 1852581
NSF-PAR ID:
10253855
Author(s) / Creator(s):
; ; ; ;
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Volume:
500
Issue:
3
ISSN:
0035-8711
Page Range / eLocation ID:
3889 to 3897
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    ABSTRACT Emission from the interstellar medium can be a significant contaminant of measurements of the intensity and polarization of the cosmic microwave background (CMB). For planning CMB observations, and for optimizing foreground-cleaning algorithms, a description of the statistical properties of such emission can be helpful. Here, we examine a machine learning approach to inferring the statistical properties of dust from observational data. In particular, we apply a type of neural network called a variational autoencoder (VAE) to maps of the intensity of emission from interstellar dust as inferred from Planck sky maps and demonstrate its ability to (i) simulate new samples with similar summary statistics as the training set, (ii) provide fits to emission maps withheld from the training set, and (iii) produce constrained realizations. We find VAEs are easier to train than another popular architecture: that of generative adversarial networks, and are better suited for use in Bayesian inference. 
    more » « less
  2. Abstract

    We characterize Galactic dust filaments by correlating BICEP/Keck and Planck data with polarization templates based on neutral hydrogen (Hi) observations. Dust polarization is important for both our understanding of astrophysical processes in the interstellar medium (ISM) and the search for primordial gravitational waves in the cosmic microwave background (CMB). In the diffuse ISM, Hiis strongly correlated with the dust and partly organized into filaments that are aligned with the local magnetic field. We analyze the deep BICEP/Keck data at 95, 150, and 220 GHz, over the low-column-density region of sky where BICEP/Keck has set the best limits on primordial gravitational waves. We separate the Hiemission into distinct velocity components and detect dust polarization correlated with the local Galactic Hibut not with the Hiassociated with Magellanic Streami. We present a robust, multifrequency detection of polarized dust emission correlated with the filamentary Himorphology template down to 95 GHz. For assessing its utility for foreground cleaning, we report that the Himorphology template correlates inBmodes at a ∼10%–65% level over the multipole range 20 << 200 with the BICEP/Keck maps, which contain contributions from dust, CMB, and noise components. We measure the spectral index of the filamentary dust component spectral energy distribution to beβ= 1.54 ± 0.13. We find no evidence for decorrelation in this region between the filaments and the rest of the dust field or from the inclusion of dust associated with the intermediate velocity Hi. Finally, we explore the morphological parameter space in the Hi-based filamentary model.

     
    more » « less
  3. Abstract

    We present the first linear polarization measurements from the 2015 long-duration balloon flight ofSpider, which is an experiment that is designed to map the polarization of the cosmic microwave background (CMB) on degree angular scales. The results from these measurements include maps and angular power spectra from observations of 4.8% of the sky at 95 and 150 GHz, along with the results of internal consistency tests on these data. While the polarized CMB anisotropy from primordial density perturbations is the dominant signal in this region of sky, Galactic dust emission is also detected with high significance. Galactic synchrotron emission is found to be negligible in theSpiderbands. We employ two independent foreground-removal techniques to explore the sensitivity of the cosmological result to the assumptions made by each. The primary method uses a dust template derived fromPlanckdata to subtract the Galactic dust signal. A second approach, which constitutes a joint analysis ofSpiderandPlanckdata in the harmonic domain, assumes a modified-blackbody model for the spectral energy distribution of the dust with no constraint on its spatial morphology. Using a likelihood that jointly samples the template amplitude andrparameter space, we derive 95% upper limits on the primordial tensor-to-scalar ratio from Feldman–Cousins and Bayesian constructions, findingr< 0.11 andr< 0.19, respectively. Roughly half the uncertainty inrderives from noise associated with the template subtraction. New data at 280 GHz fromSpider’s second flight will complement thePlanckpolarization maps, providing powerful measurements of the polarized Galactic dust emission.

     
    more » « less
  4. ABSTRACT

    In order to better analyse the polarization of the cosmic microwave background (CMB), which is dominated by emission from our Galaxy, we need tools that can detect residual foregrounds in cleaned CMB maps. Galactic foregrounds introduce statistical anisotropy and directionality to the polarization pseudo-vectors of the CMB, which can be investigated by using the $\mathcal {D}$ statistic of Bunn and Scott. This statistic is rapidly computable and capable of investigating a broad range of data products for directionality. We demonstrate the application of this statistic to detecting foregrounds in polarization maps by analysing the uncleaned Planck 2018 frequency maps. For the Planck 2018 CMB maps, we find no evidence for residual foreground contamination. In order to examine the sensitivity of the $\mathcal {D}$ statistic, we add a varying fraction of the polarized thermal dust and synchrotron foreground maps to the CMB maps and show the per cent-level foreground contamination that would be detected with 95 per cent confidence. We also demonstrate application of the $\mathcal {D}$ statistic to another data product by analysing the gradient of the minimum-variance CMB lensing potential map (i.e. the deflection angle) for directionality. We find no excess directionality in the lensing potential map when compared to the simulations provided by the Planck Collaboration.

     
    more » « less
  5. For the past decade, the BICEP/Keck collaboration has been operating a series of telescopes at the Amundsen-Scott South Pole Station measuring degree-scale B-mode polarization imprinted in the Cosmic Microwave Background (CMB) by primordial gravitational waves (PGWs). These telescopes are compact refracting polarimeters mapping about 2% of the sky, observing at a broad range of frequencies to account for the polarized foreground from Galactic synchrotron and thermal dust emission. Our latest publication "BK18" utilizes the data collected up to the 2018 observing season, in conjunction with the publicly available WMAP and Planck data, to constrain the tensor-to-scalar ratio r. It particularly includes (1) the 3-year BICEP3 data which is the current deepest CMB polarization map at the foreground-minimum 95 GHz; and (2) the Keck 220 GHz map with a higher signal-to-noise ratio on the dust foreground than the Planck 353 GHz map. We fit the auto- and cross-spectra of these maps to a multicomponent likelihood model (ΛCDM+dust+synchrotron+noise+r) and find it to be an adequate description of the data at the current noise level. The likelihood analysis yields σ(r)=0.009. The inference of r from our baseline model is tightened to r0.05=0.014+0.010−0.011 and r0.05<0.036 at 95% confidence, meaning that the BICEP/Keck B-mode data is the most powerful existing dataset for the constraint of PGWs. The up-coming BICEP Array telescope is projected to reach σ(r)≲0.003 using data up to 2027. 
    more » « less