Search

With a single training image and using wavelet phase harmonic augmentation, we present polarized Cosmic Microwave Background (CMB) foreground marginalization in a high-dimensional likelihood-free (Bayesian) framework. We demonstrate robust foreground removal using only a single frequency of simulated data for a BICEP-like sky patch. Using Moment Networks, we estimate the pixel-level posterior probability for the underlying { E, B } signal and validate the statistical model with a quantile-type test using the estimated marginal posterior moments. The Moment Networks use a hierarchy of U-Net convolutional neural networks. This work validates such an approach in the most difficult limiting case: pixel-level, noise-free, highly non-Gaussian dust foregrounds with a single training image at a single frequency. For a real CMB experiment, a small number of representative sky patches would provide the training data required for full cosmological inference. These results enable robust likelihood-free , simulation-based parameter, and model inference for primordial B-mode detection using observed CMB polarization data. [ABSTRACT FROM AUTHOR]

A foreword of the journal "Interstellar Dust from Astronomical Observations to Fundamental Studies," is presented.

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 multifrequency 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 multichannel 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: (1) a (I, E, B) multi-observable input, and (2) a { I ν } ν multifrequency input. The samples exhibit consistent features compared to the original maps. A statistical analysis of model 1 shows that the power spectra, distributions of pixels, and Minkowski functionals are captured to a good extent. We analyze model 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. [ABSTRACT FROM AUTHOR]

The zodiacal light is the dominant source of the mid-infrared sky brightness seen from Earth, and exozodiacal light is the dominant emission from planetary and debris systems around other stars. We observed the zodiacal light spectrum with the mid-infrared camera ISOCAM over the wavelength range 5–16 μm and a wide range of orientations relative to the Sun (solar elongations 68°–113°) and the ecliptic (plane to pole). The temperature in the ecliptic ranged from 269 K at solar elongation 68° to 244 K at 113°, and the polar temperature, characteristic of dust 1 AU from the Sun, is 274 K. The observed temperature is exactly as expected for large (>10 μm radius), low-albedo (<0.08), rapidly-rotating, gray particles 1 AU from the Sun. Smaller particles (<10 μm radius) radiate inefficiently in the infrared and are warmer than observed. We present theoretical models for a wide range of particle size distributions and compositions; it is evident that the zodiacal light is produced by particles in the 10–100 μm radius range. In addition to the continuum, we detect a weak excess in the 9–11 μm range, with an amplitude of 6% of the continuum. The shape of the feature can be matched by a mixture of silicates: amorphous forsterite/olivine provides most of the continuum and some of the 9–11 μm silicate feature, dirty crystalline olivine provides the red wing of the silicate feature (and a bump at 11.35 μm), and a hydrous silicate (montmorillonite) provides the blue wing of the silicate feature. The presence of hydrous silicate suggests the parent bodies of those particles were formed in the inner solar nebula. Large particles dominate the size distribution, but at least some small particles (radii ∼1 μm) are required to produce the silicate emission feature. The strength of the feature may vary spatially, with the strongest features being at the lowest solar elongations as well as at high ecliptic latitudes; if confirmed, this would imply that the dust properties change such that dust further from the Sun has a weaker silicate feature. To compare the properties of zodiacal dust to dust around other main sequence stars, we reanalyzed the exozodiacal light spectrum for β Pic to derive the shape of its silicate feature. The zodiacal and exozodiacal spectra are very different. The exozodiacal spectra are dominated by cold dust, with emission peaking in the far-infrared, while the zodiacal spectrum peaks around 20 μm. We removed the debris disk continuum from the spectra by fitting a blackbody with a different temperature for each aperture (ranging from 3.7″ to 27″); the resulting silicate spectra for β Pic are identical for all apertures, indicating that the silicate feature arises close to the star. The shape of the silicate feature from β Pic is nearly identical to that derived from the ISO spectrum of 51 Oph; both exozodiacal features are very different from that of the zodiacal light. The exozodiacal features are roughly triangular, peaking at 10.3 μm, while the zodiacal feature is more boxy, indicating a different mineralogy. [Copyright &y& Elsevier]