Generative Models of Multi-channel Data from a Single Example -- Application to Dust Emission

Abstract

The quest for primordial $B$-modes in the cosmic microwave background has emphasized the need for refined models of the Galactic dust foreground. Here, we aim at building a realistic statistical model of the multi-frequency dust emission from a single example. We introduce a generic methodology relying on microcanonical gradient descent models conditioned by an extended family of wavelet phase harmonic (WPH) statistics. To tackle the multi-channel aspect of the data, we define cross-WPH statistics, quantifying non-Gaussian correlations between maps. Our data-driven methodology could apply to various contexts, and we have updated the software PyWPH, on which this work relies, accordingly. Applying this to dust emission maps built from a magnetohydrodynamics simulation, we construct and assess two generative models of: 1) a $(I,E,B)$ multi-observable input, 2) a $\{I_{\nu }{\}}_{\nu }$ multi-frequency input. The samples exhibit consistent features compared to the original maps. A statistical analysis of 1) shows that the power spectra, distributions of pixels, and Minkowski functionals are captured to a good extent. We analyze 2) by fitting the spectral energy distribution (SED) of both the synthetic and original maps with a modified blackbody (MBB) law. The maps are equally well fitted, and a comparison of the MBB parameters shows that our model succeeds in capturing the spatial variations of the SED from the data. Besides the perspectives of this work for dust emission modeling, the introduction of cross-WPH statistics opens a new avenue to characterize non-Gaussian interactions across different maps, which we believe will be fruitful for astrophysics.

18 pages, 7 figures, submitted to ApJ, code: https://github.com/bregaldo/dust_genmodels

References in corpus (16)

• Array Programming with NumPy
• A Joint Analysis of BICEP2/Keck Array and Planck Data
• A High Order Godunov Scheme with Constrained Transport and Adaptive Mesh Refinement for Astrophysical MHD
• Planck intermediate results. XX. Comparison of polarized thermal emission from Galactic dust with simulations of MHD turbulence
• 3D chemical structure of diffuse turbulent ISM. I. Statistics of the HI-to-H${}_2$ transition
• Going deep with Minkowski functionals of convergence maps
• Statistical description of dust polarized emission from the diffuse interstellar medium -- A RWST approach
• ForSE: a GAN based algorithm for extending CMB foreground models to sub-degree angular scales
• High precision modeling of polarized signals: Moment expansion method generalized to spin-2 fields
• A new approach for the statistical denoising of Planck interstellar dust polarization data
• Non-Gaussian modelling and statistical denoising of Planck dust polarization full-sky maps using scattering transforms
• Single frequency CMB B-mode inference with realistic foregrounds from a single training image
• Full-sky, arcminute-scale, 3D models of Galactic microwave foreground dust emission based on filaments
• A Generative Model of Galactic Dust Emission Using Variational Inference
• Characterizing line-of-sight variability of polarized dust emission with future CMB experiments
• Diagnosing Turbulence in the Neutral and Molecular Interstellar Medium of Galaxies

Cited by in corpus (6)

• Scattering Spectra Models for Physics
• Generative models of astrophysical fields with scattering transforms on the sphere
• Diffusion-HMC: Parameter Inference with Diffusion-model-driven Hamiltonian Monte Carlo
• PCAT-DE: Reconstructing point-like and diffuse signals in astronomical images using spatial and spectral information
• From few to many maps: A fast map-level emulator for extreme augmentation of CMB systematics datasets
• Mallat Scattering Transformation based surrogate for MagnetoHydroDynamics