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.

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

References in corpus (44)

• SciPy 1.0--Fundamental Algorithms for Scientific Computing in Python
• Array Programming with NumPy
• Cosmological Hydrodynamics with Adaptive Mesh Refinement: a new high resolution code called RAMSES
• Planck 2018 results. I. Overview and the cosmological legacy of Planck
• The Simons Observatory: Science goals and forecasts
• A Joint Analysis of BICEP2/Keck Array and Planck Data
• A Probe of Primordial Gravity Waves and Vorticity
• Planck 2013 results. XI. All-sky model of thermal dust emission
• Understanding Deep Convolutional Networks
• Planck 2018 results. VII. Isotropy and Statistics of the CMB
• A High Order Godunov Scheme with Constrained Transport and Adaptive Mesh Refinement for Astrophysical MHD
• Planck 2018 results. IV. Diffuse component separation
• The Quest for B Modes from Inflationary Gravitational Waves
• The Python Sky Model: software for simulating the Galactic microwave sky
• The pre-launch Planck Sky Model: a model of sky emission at submillimetre to centimetre wavelengths
• Planck 2018 results. XI. Polarized dust foregrounds
• The Atacama Cosmology Telescope: CMB Polarization at $200<\ell <9000$
• Planck intermediate results. XX. Comparison of polarized thermal emission from Galactic dust with simulations of MHD turbulence
• Rethinking CMB foregrounds: systematic extension of foreground parameterizations
• New Interpretable Statistics for Large Scale Structure Analysis and Generation
• SPIDER: Probing the Early Universe with a Suborbital Polarimeter
• Modeling the magnetized Local Bubble from dust data
• The nature of the E-B decomposition of CMB polarization
• Non-Gaussian Minkowski functionals & extrema counts in redshift space
• The RWST, a comprehensive statistical description of the non-Gaussian structures in the ISM
• Mitigating Complex Dust Foregrounds in Future CMB Polarization Experiments
• 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
• Local Bubble contribution to the 353-GHz dust polarized emission
• High precision modeling of polarized signals: Moment expansion method generalized to spin-2 fields
• ForSE: a GAN based algorithm for extending CMB foreground models to sub-degree angular scales
• Dust moments: towards a new modelling of the galactic dust emission for CMB B-modes analysis
• A new approach for the statistical denoising of Planck interstellar dust polarization data
• A systematic study of radiative torque grain alignment in the diffuse interstellar medium
• Maximum Entropy Models from Phase Harmonic Covariances
• Cleaning our own Dust: Simulating and Separating Galactic Dust Foregrounds with Neural Networks
• 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
• On the statistics of the polarized submillimetre emission maps from thermal dust in the turbulent, magnetized, diffuse ISM
• 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