skip to main content

Explore Research Products in the NSF-PAR

Title: A generative model of galactic dust emission using variational autoencoders
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
Award ID(s):
1852617
NSF-PAR ID:
10232682
Author(s) / Creator(s):
; ;
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Volume:
504
Issue:
2
ISSN:
0035-8711
Page Range / eLocation ID:
2603 to 2613
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ; ; ; ; ( , Monthly Notices of the Royal Astronomical Society)
    null (Ed.)
    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
  2. ; ; ; ; ; ; ; ; ; ; et al ( , The Astrophysical Journal)
    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. ; ; ; ; ; ; ( , Astronomy & Astrophysics)
    Context. The formation of molecular gas in interstellar clouds is a slow process, but can be enhanced by gas compression. Magneto-hydrodynamic (MHD) waves can create compressed quasi-periodic linear structures, referred to as striations. Striations are observed at the column densities at which the transition from atomic to molecular gas takes place. Aims. We explore the role of MHD waves in the CO chemistry in regions with striations within molecular clouds. Methods. We targeted a region with striations in the Polaris Flare cloud. We conducted a CO J = 2−1 survey in order to probe the molecular gas properties. We used archival starlight polarization data and dust emission maps in order to probe the magnetic field properties and compare against the CO morphological and kinematic properties. We assessed the interaction of compressible MHD wave modes with CO chemistry by comparing their characteristic timescales. Results. The estimated magnetic field is 38–76 µG. In the CO integrated intensity map, we observe a dominant quasiperiodic intensity structure that tends to be parallel to the magnetic field orientation and has a wavelength of approximately one parsec. The periodicity axis is ~17° off from the mean magnetic field orientation and is also observed in the dust intensity map. The contrast in the CO integrated intensity map is ~2.4 times higher than the contrast of the column density map, indicating that CO formation is enhanced locally. We suggest that a dominant slow magnetosonic mode with an estimated period of 2.1–3.4 Myr and a propagation speed of 0.30–0.45 km s −1 is likely to have enhanced the formation of CO, hence created the observed periodic pattern. We also suggest that within uncertainties, a fast magnetosonic mode with a period of 0.48 Myr and a velocity of 2.0 km s −1 could have played some role in increasing the CO abundance. Conclusions. Quasiperiodic CO structures observed in striation regions may be the imprint of MHD wave modes. The Alfvénic speed sets the dynamical timescales of the compressible MHD modes and determines which wave modes are involved in the CO chemistry. 
    more » « less
  4. ; ; ; ( , Monthly Notices of the Royal Astronomical Society)
    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. ; ; ; ; ; ; ; ( , Monthly Notices of the Royal Astronomical Society)
    ABSTRACT Line intensity mapping (LIM) is rapidly emerging as a powerful technique to study galaxy formation and cosmology in the high-redshift Universe. We present LIM estimates of select spectral lines originating from the interstellar medium (ISM) of galaxies and 21 cm emission from neutral hydrogen gas in the Universe using the large volume, high resolution thesan reionization simulations. A combination of subresolution photoionization modelling for H ii regions and Monte Carlo radiative transfer calculations is employed to estimate the dust-attenuated spectral energy distributions (SEDs) of high-redshift galaxies (z ≳ 5.5). We show that the derived photometric properties such as the ultraviolet (UV) luminosity function and the UV continuum slopes match observationally inferred values, demonstrating the accuracy of the SED modelling. We provide fits to the luminosity–star formation rate relation (L–SFR) for the brightest emission lines and find that important differences exist between the derived scaling relations and the widely used low-z ones because the ISM of reionization era galaxies is generally less metal enriched than in their low-redshift counterparts. We use these relations to construct line intensity maps of nebular emission lines and cross-correlate with the 21 cm emission. Interestingly, the wavenumber at which the correlation switches sign (ktransition) depends heavily on the reionization model and to a lesser extent on the targeted emission line, which is consistent with the picture that ktransition probes the typical sizes of ionized regions. The derived scaling relations and intensity maps represent a timely state-of-the-art framework for forecasting and interpreting results from current and upcoming LIM experiments. 
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email