Search

Because of their limited angular resolution, far-infrared telescopes are usually affected by confusion phenomenon. Since several galaxies can be located in the same instrumental beam, only the brightest objects emerge from the fluctuations caused by fainter sources. The probe far-infrared mission for astrophysics imager (PRIMAger) will observe the mid- and far-infrared (25-235 $\mu$m) sky both in intensity and polarization. We aim to provide predictions of the confusion level and its consequences for future surveys. We produced simulated PRIMAger maps affected only by the confusion noise using the simulated infrared extragalactic sky (SIDES) semi-empirical simulation. We then estimated the confusion limit in these maps and extracted the sources using a basic blind extractor. By comparing the input galaxy catalog and the extracted source catalog, we derived various performance metrics as completeness, purity, and the accuracy of various measurements. In intensity, we predict that the confusion limit increases rapidly with increasing wavelength. The confusion limit in polarization is more than 100x lower. The measured flux density is dominated by the brightest galaxy in the beam, but other objects also contribute at longer wavelength (~30% at 235 $\mu$m). We also show that galaxy clustering has a mild impact on confusion in intensity (up to 25%), while it is negligible in polarization. In intensity, a basic blind extraction will be sufficient to detect galaxies at the knee of the luminosity function up to z~3 and 10$^{11}$ M$_\odot$ main-sequence galaxies up to z~5. In polarization for a conservative sensitivity, we expect ~8 000 detections up to z=2.5 opening a totally new window on the high-z dust polarization. Finally, we show that intensity surveys at short wavelength and polarization surveys at long wavelength tend to reach confusion at similar depth. There is thus a strong synergy., Comment: 24 pages, 18 figures, 6 tables, Accepted by A&A

We propose a new approach to improve the precision of astrophysical parameter constraints for the 21cm signal from the epoch of reionisation (EoR). Our method introduces new sets of summary statistics, hereafter `evolution compressed' statistics, which quantify the spectral evolution of the 2D spatial statistics computed at fixed redshift. We defined such compressed statistics for power spectrum (PS), wavelet scattering transforms (WST), and wavelet moments (WM), which also characterise non-Gaussian features. To compare these different statistics, along with the 3D power spectrum, we estimated their Fisher information on three cosmological parameters from an ensemble of simulations of 21cm EoR data, both in noiseless and noisy scenarios using Square Kilometre Array (SKA) noise levels equivalent to 100 and 1000 hours of observations. We also compare wavelet statistics, in particular WST, built from standard directional Morlet wavelets, as well as from a set of isotropic wavelets derived from the binning window function of the 2D power spectrum. For the noiseless case, the compressed wavelet statistics give constraints that are up to five times more precise than those obtained from the 3D isotropic power spectrum. At the same time, for 100h SKA noise, from which it is difficult to extract non-Gaussian features, compressed wavelet statistics still give over 30\% tighter constraints. We find that the wavelet statistics with wavelets derived from the power-spectrum binning window function provide the tightest constraints of all the statistics, with the WSTs seemingly performing better than the WMs, in particular when working with noisy data. The findings of this study demonstrate that evolution-compressed statistics extract more information than usual 3D isotropic power-spectra approaches and that our wavelet-based statistics can consistently outmatch power-spectrum-based statistics., Comment: 17 pages, 8 figures. Accepted to A&A. Difference to version 1 includes the addition of the same wavelet used for the Wavelet Moments application within the Wavelet Scattering Transform formalism, as well as better normalising of the scattering coefficients

We present a Bayesian inference method to characterise the dust emission properties using the well-known dust-HI correlation in the diffuse interstellar medium at Planck frequencies $\nu \ge 217$ GHz. We use the Galactic HI map from the Galactic All-Sky Survey (GASS) as a template to trace the Galactic dust emission. We jointly infer the pixel-dependent dust emissivity and the zero level present in the Planck intensity maps. We use the Hamiltonian Monte Carlo technique to sample the high dimensional parameter space ($D \sim 10^3$). We demonstrate that the methodology leads to unbiased recovery of dust emissivity per pixel and the zero level when applied to realistic Planck sky simulations over a 6300 deg$^2$ area around the Southern Galactic pole. As an application on data, we analyse the Planck intensity map at 353 GHz to jointly infer the pixel-dependent dust emissivity at Nside=32 resolution (1.8\deg\ pixel size) and the global offset. We find that the spatially varying dust emissivity has a mean of 0.031 MJysr$^{-1} (10^{20} \mathrm{cm^{-2}})^{-1}$ and $1\sigma$ standard deviation of 0.007 MJysr$^{-1} (10^{20} \mathrm{cm^{-2}})^{-1}$. The mean dust emissivity increases monotonically with increasing mean HI column density. We find that the inferred global offset is consistent with the expected level of Cosmic Infrared Background (CIB) monopole added to the Planck data at 353 GHz. This method is useful in studying the line-of-sight variations of dust spectral energy distribution in the multi-phase interstellar medium., Comment: 19 pages, 14 figures, Accepted for publication in MNRAS

The low brightness dust emission at high Galactic latitude is of interest to study the interplay between physical processes in shaping the structure of the interstellar medium (ISM), as well as to statistically characterize dust emission as a foreground to the Cosmic Microwave Background (CMB). Progress in this avenue of research have been hampered by the difficulty of separating the dust emission from the Cosmic Infrared Background (CIB). We demonstrate that dust and CIB may be effectively separated based on their different structure on the sky and use the separation to characterize the structure of diffuse dust emission on angular scales where CIB is a significant component in terms of power. We use scattering transform statistics, the Wavelet Phase Harmonics (WPH), to perform a statistical component separation using Herschel SPIRE observations. This component separation is done only from observational data using non-Gaussian properties as a lever arm, and is done at a single 250 microns frequency. This method, that we validate on mock data, gives us access to non-Gaussian statistics of the interstellar dust and an output dust map essentially free from CIB contamination. Our statistical modelling characterizes the non-Gaussian structure of the diffuse ISM down to the smallest scales observed by Herschel. We recover the power-law shape of the dust power spectrum up to a wavenumber of 2 arcmin$^{-1}$ where the dust signal represents 2 percent of the total power. The output dust map reveals coherent structures at the smallest scales which were hidden by the CIB anisotropies. It opens new observational perspectives on the formation of structure in the diffuse ISM which we discuss with reference to past work. We have succeeded to perform a statistical separation from observational data only at a single frequency by using non-Gaussian statistics., Comment: Accepted in A&A on October 23, 2023

The change of physical conditions across the turbulent and magnetized interstellar medium (ISM) induces a 3D spatial variation of the properties of Galactic polarized emission. The observed signal results from the averaging of different spectral energy distributions (SED) and polarization angles, along and between lines of sight. As a consequence, the total Stokes parameters $Q$ and $U$ will have different distorted SEDs, so that the polarization angle becomes frequency dependent. In the present work, we show how this phenomenon similarly induces a different distorted SED for the three polarized angular power spectra $EE$, $BB$ and $EB$, implying a variation of the $EE/BB$ ratio with frequency. We demonstrate how the previously introduced spin-moment formalism provides a natural framework to grasp these effects, allowing us to derive analytical predictions for the spectral behaviors of the polarized spectra, focusing here on the example of thermal dust polarized emission. After a quantitative discussion based on a model combining emission from a filament with its background, we further reveal that the spectral complexity implemented in the dust models commonly used by the cosmic microwave background (CMB) community produce such effects. This new understanding is crucial for CMB component separation, in which an extreme accuracy is required in the modeling of the dust signal to allow for the search of the primordial imprints of inflation or cosmic birefringence. For the latter, as long as the dust $EB$ signal is not measured accurately, great caution is required about the assumptions made to model its spectral behavior, as it may not simply follow from the other dust angular power spectra., Comment: Published in A&A

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., Comment: 18 pages, 7 figures, submitted to ApJ, code: https://github.com/bregaldo/dust_genmodels

The search for the primordial B-modes of the cosmic microwave background (CMB) relies on the separation from the brighter foreground dust signal. In this context, the characterisation of the spectral energy distribution (SED) of thermal dust in polarization has become a critical subject of study. We present a power-spectra analysis of Planck data, which improves on previous studies by using the newly released SRoll2 maps that correct residual data systematics, and by extending the analysis to regions near the Galactic plane. Our analysis focuses on the lowest multipoles between l=4 and 32, and three sky areas with sky fractions of fsky = 80%, 90%, and 97%. The mean dust SED for polarization and the 353 GHz Q and U maps are used to compute residual maps at 100, 143 and 217 GHz, highlighting spatial variations of the dust polarization SED. Residuals are detected at the three frequencies for the three sky areas. We show that models based on total intensity data are underestimating by a significant factor the complexity of dust polarized CMB foreground. Our analysis emphasizes the need to include variations of polarization angles of the dust polarized CMB foreground. The frequency dependence of the EE and BB power spectra of the residual maps yields further insight. We find that the moments expansion to the first order of the modified black-body (MBB) spectrum provides a good fit to the EE power-spectra. This result suggests that the residuals could follow mainly from variations of dust MBB spectral parameters. However, this conclusion is challenged by cross-spectra showing that the residuals maps at the three frequencies are not fully correlated, and the fact that the BB power-spectra do not match the first order moment expansion of a MBB SED. This work sets new requirements for simulations of the dust polarized foreground and component separation methods (abridged), Comment: 13 pages, 13 figures, Accepted for publication in the Cosmology section of A&A journal

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., Comment: Accepted by Monthly Notices of the Royal Astronomical Society Letters. 5 pages with 3 figures (plus 1 page of Supporting Materials with 2 figures)

Dust emission is the main foreground for cosmic microwave background (CMB) polarization. Its statistical characterization must be derived from the analysis of observational data because the precision required for a reliable component separation is far greater than what is currently achievable with physical models of the turbulent magnetized interstellar medium. This letter takes a significant step toward this goal by proposing a method that retrieves non-Gaussian statistical characteristics of dust emission from noisy Planck polarization observations at 353 GHz. We devised a statistical denoising method based on wavelet phase harmonics (WPH) statistics, which characterize the coherent structures in non-Gaussian random fields and define a generative model of the data. The method was validated on mock data combining a dust map from a magnetohydrodynamic simulation and Planck noise maps. The denoised map reproduces the true power spectrum down to scales where the noise power is an order of magnitude larger than that of the signal. It remains highly correlated to the true emission and retrieves some of its non-Gaussian properties. Applied to Planck data, the method provides a new approach to building a generative model of dust polarization that will characterize the full complexity of the dust emission. We also release PyWPH, a public Python package, to perform GPU-accelerated WPH analyses on images., Comment: 10 pages, 9 figures, accepted by A&A

The statistical characterization of the diffuse magnetized ISM and Galactic foregrounds to the CMB poses a major challenge. To account for their non-Gaussian statistics, we need a data analysis approach capable of efficiently quantifying statistical couplings across scales. This information is encoded in the data, but most of it is lost when using conventional tools, such as one-point statistics and power spectra. The wavelet scattering transform (WST), a low-variance statistical descriptor of non-Gaussian processes introduced in data science, opens a path towards this goal. We applied the WST to noise-free maps of dust polarized thermal emission computed from a numerical simulation of MHD turbulence. We analyzed normalized complex Stokes maps and maps of the polarization fraction and polarization angle. The WST yields a few thousand coefficients; some of them measure the amplitude of the signal at a given scale, and the others characterize the couplings between scales and orientations. The dependence on orientation can be fitted with the reduced WST (RWST), an angular model introduced in previous works. The RWST provides a statistical description of the polarization maps, quantifying their multiscale properties in terms of isotropic and anisotropic contributions. It allowed us to exhibit the dependence of the map structure on the orientation of the mean magnetic field and to quantify the non-Gaussianity of the data. We also used RWST coefficients, complemented by additional constraints, to generate random synthetic maps with similar statistics. Their agreement with the original maps demonstrates the comprehensiveness of the statistical description provided by the RWST. This work is a step forward in the analysis of observational data and the modeling of CMB foregrounds. We also release PyWST, a Python package to perform WST/RWST analyses at: https://github.com/bregaldo/pywst., Comment: 20 pages, 15 figures, accepted by Astronomy & Astrophysics

The characterization of the spectral energy distribution (SED) of dust emission has become a critical issue in the quest for primordial B-modes. The dust SED is often approximated by a modified black body (MBB) emission law but the extent to which this is accurate is unclear. This paper addresses this question, expanding the dust SED at the power spectrum level. The expansion is performed by means of moments around the MBB law, related to derivatives with respect to the dust spectral index. We present the mathematical formalism and apply it to simulations and Planck total intensity data, from 143 to 857 GHz, because no polarized data are yet available that provide the required sensitivity to perform this analysis. With simulations, we demonstrate the ability of high-order moments to account for spatial variations in MBB parameters. Neglecting these moments leads to poor fits and a bias in the recovered dust spectral index. We identify the main moments that are required to fit the Planck data. The comparison with simulations helps us to disentangle the respective contributions from dust and the cosmic infrared background to the high-order moments, but the simulations give an insufficient description of the actual Planck data. Extending our model to cosmic microwave background B-mode analyses within a simplified framework, we find that ignoring the dust SED distortions, or trying to model them with a single decorrelation parameter, could lead to biases that are larger than the targeted sensitivity for the next generation of CMB B-mode experiments., Comment: 18 pages, 19 figures, accepted in A&A: A&A 647, A52 (2021)

The analysis of the Planck polarization E and B mode power spectra of interstellar dust emission at 353 GHz recently raised new questions on the impact of Galactic foregrounds to the detection of the polarization of the Cosmic Microwave Background (CMB) and on the physical properties of the interstellar medium (ISM). In the diffuse ISM a clear E-B asymmetry is observed, with twice as much power in E modes than in B modes; as well as a positive correlation between the total power, T, and both E and B modes, presently interpreted in terms of the link between the structure of interstellar matter and that of the Galactic magnetic field. In this paper we aim at extending the Planck analysis of the high-latitude sky to low Galactic latitude, investigating the correlation between the TEB power spectra with the gas column density from the diffuse ISM to molecular clouds. We divide the sky between Galactic latitude |b|>5 deg and |b|<60 deg in 552 circular patches and we study the cross-correlations between the TEB power spectra and the column density of each patch using the latest release of the Planck polarization data. We find that the B-to-E power ratio (BB/EE) and the TE correlation ratio (rTE) depend on column density. While the former increases going from the diffuse ISM to molecular clouds in the Gould Belt, the latter decreases. This systematic variation must be related to actual changes in ISM properties. The data show significant scatter about this mean trend. The variations of BB/EE and rTE are observed to be anti-correlated for all column densities. In the diffuse ISM, the variance of these two ratios is consistent with a stochastic non-Gaussian model in which the values of BB/EE and rTE are fixed. We finally discuss the dependencies of TB and EB with column density, which are however hampered by instrumental noise., Comment: Submitted to A&A on May 24, 2019. Accepted on October 14, 2019

The Probe of Inflation and Cosmic Origins (PICO) is an imaging polarimeter that will scan the sky for 5 years in 21 frequency bands spread between 21 and 799 GHz. It will produce full-sky surveys of intensity and polarization with a final combined-map noise level of 0.87 $\mu$K arcmin for the required specifications, equivalent to 3300 Planck missions, and with our current best-estimate would have a noise level of 0.61 $\mu$K arcmin (6400 Planck missions). PICO will either determine the energy scale of inflation by detecting the tensor to scalar ratio at a level $r=5\times 10^{-4}~(5\sigma)$, or will rule out with more than $5\sigma$ all inflation models for which the characteristic scale in the potential is the Planck scale. With LSST's data it could rule out all models of slow-roll inflation. PICO will detect the sum of neutrino masses at $>4\sigma$, constrain the effective number of light particle species with $\Delta N_{\rm eff}<0.06~(2\sigma)$, and elucidate processes affecting the evolution of cosmic structures by measuring the optical depth to reionization with errors limited by cosmic variance and by constraining the evolution of the amplitude of linear fluctuations $\sigma_{8}(z)$ with sub-percent accuracy. Cross-correlating PICO's map of the thermal Sunyaev-Zeldovich effect with LSST's gold sample of galaxies will precisely trace the evolution of thermal pressure with $z$. PICO's maps of the Milky Way will be used to determine the make up of galactic dust and the role of magnetic fields in star formation efficiency. With 21 full sky legacy maps in intensity and polarization, which cannot be obtained in any other way, the mission will enrich many areas of astrophysics. PICO is the only single-platform instrument with the combination of sensitivity, angular resolution, frequency bands, and control of systematic effects that can deliver this compelling, timely, and broad science., Comment: Probe class mission study submitted to NASA and 2020 Decadal Panel. Executive summary: 2.5 pages; Science: 28 pages; Total: 50 pages, 36 figures

BFORE is a high-altitude ultra-long-duration balloon mission to map the cosmic microwave background (CMB). During a 28-day mid-latitude flight launched from Wanaka, New Zealand, the instrument will map half the sky to improve measurements of the optical depth to reionization tau. This will break parameter degeneracies needed to detect neutrino mass. BFORE will also hunt for the gravitational wave B-mode signal, and map Galactic dust foregrounds. The mission will be the first near-space use of TES/mSQUID multichroic detectors (150/217 GHz and 280/353 GHz bands) with low-power readout electronics., Comment: 11 pages, 2 figures, Submitted to SPIE Astronomical Telescopes and Instrumentation

Lensing of the CMB is now a well-developed probe of large-scale clustering over a broad range of redshifts. By exploiting the non-Gaussian imprints of lensing in the polarization of the CMB, the CORE mission can produce a clean map of the lensing deflections over nearly the full-sky. The number of high-S/N modes in this map will exceed current CMB lensing maps by a factor of 40, and the measurement will be sample-variance limited on all scales where linear theory is valid. Here, we summarise this mission product and discuss the science that it will enable. For example, the summed mass of neutrinos will be determined to an accuracy of 17 meV combining CORE lensing and CMB two-point information with contemporaneous BAO measurements, three times smaller than the minimum total mass allowed by neutrino oscillations. In the search for B-mode polarization from primordial gravitational waves with CORE, lens-induced B-modes will dominate over instrument noise, limiting constraints on the gravitational wave power spectrum amplitude. With lensing reconstructed by CORE, one can "delens" the observed polarization internally, reducing the lensing B-mode power by 60%. This improves to 70% by combining lensing and CIB measurements from CORE, reducing the error on the gravitational wave amplitude by 2.5 compared to no delensing (in the null hypothesis). Lensing measurements from CORE will allow calibration of the halo masses of the 40000 galaxy clusters that it will find, with constraints dominated by the clean polarization-based estimators. CORE can accurately remove Galactic emission from CMB maps with its 19 frequency channels. We present initial findings that show that residual Galactic foreground contamination will not be a significant source of bias for lensing power spectrum measurements with CORE. [abridged], Comment: 44 pages, 12 figures

BFORE is a NASA high-altitude ultra-long-duration balloon mission proposed to measure the cosmic microwave background (CMB) across half the sky during a 28-day mid-latitude flight launched from Wanaka, New Zealand. With the unique access to large angular scales and high frequencies provided by the balloon platform, BFORE will significantly improve measurements of the optical depth to reionization tau, breaking parameter degeneracies needed for a measurement of neutrino mass with the CMB. The large angular scale data will enable BFORE to hunt for the large-scale gravitational wave B-mode signal, as well as the degree-scale signal, each at the r~0.01 level. The balloon platform allows BFORE to map Galactic dust foregrounds at frequencies where they dominate, in order to robustly separate them from CMB signals measured by BFORE, in addition to complementing data from ground-based telescopes. The combination of frequencies will also lead to velocity measurements for thousands of galaxy clusters, as well as probing how star-forming galaxies populate dark matter halos. The mission will be the first near-space use of TES multichroic detectors (150/217 GHz and 280/353 GHz bands) using highly-multiplexed mSQUID microwave readout, raising the technical readiness level of both technologies., Comment: 8 pages, 4 figures, to appear in Proceedings of Low Temperature Detectors 17

We forecast the scientific capabilities to improve our understanding of cosmic inflation of CORE, a proposed CMB space satellite submitted in response to the ESA fifth call for a medium-size mission opportunity. The CORE satellite will map the CMB anisotropies in temperature and polarization in 19 frequency channels spanning the range 60-600 GHz. CORE will have an aggregate noise sensitivity of $1.7 \mu$K$\cdot \,$arcmin and an angular resolution of 5' at 200 GHz. We explore the impact of telescope size and noise sensitivity on the inflation science return by making forecasts for several instrumental configurations. This study assumes that the lower and higher frequency channels suffice to remove foreground contaminations and complements other related studies of component separation and systematic effects, which will be reported in other papers of the series "Exploring Cosmic Origins with CORE." We forecast the capability to determine key inflationary parameters, to lower the detection limit for the tensor-to-scalar ratio down to the $10^{-3}$ level, to chart the landscape of single field slow-roll inflationary models, to constrain the epoch of reheating, thus connecting inflation to the standard radiation-matter dominated Big Bang era, to reconstruct the primordial power spectrum, to constrain the contribution from isocurvature perturbations to the $10^{-3}$ level, to improve constraints on the cosmic string tension to a level below the presumptive GUT scale, and to improve the current measurements of primordial non-Gaussianities down to the $f_{NL}^{\rm local} < 1$ level. For all the models explored, CORE alone will improve significantly on the present constraints on the physics of inflation. Its capabilities will be further enhanced by combining with complementary future cosmological observations., Comment: Latex 107 pages, revised with updated author list and minor modifications

We forecast the main cosmological parameter constraints achievable with the CORE space mission which is dedicated to mapping the polarisation of the Cosmic Microwave Background (CMB). CORE was recently submitted in response to ESA's fifth call for medium-sized mission proposals (M5). Here we report the results from our pre-submission study of the impact of various instrumental options, in particular the telescope size and sensitivity level, and review the great, transformative potential of the mission as proposed. Specifically, we assess the impact on a broad range of fundamental parameters of our Universe as a function of the expected CMB characteristics, with other papers in the series focusing on controlling astrophysical and instrumental residual systematics. In this paper, we assume that only a few central CORE frequency channels are usable for our purpose, all others being devoted to the cleaning of astrophysical contaminants. On the theoretical side, we assume LCDM as our general framework and quantify the improvement provided by CORE over the current constraints from the Planck 2015 release. We also study the joint sensitivity of CORE and of future Baryon Acoustic Oscillation and Large Scale Structure experiments like DESI and Euclid. Specific constraints on the physics of inflation are presented in another paper of the series. In addition to the six parameters of the base LCDM, which describe the matter content of a spatially flat universe with adiabatic and scalar primordial fluctuations from inflation, we derive the precision achievable on parameters like those describing curvature, neutrino physics, extra light relics, primordial helium abundance, dark matter annihilation, recombination physics, variation of fundamental constants, dark energy, modified gravity, reionization and cosmic birefringence. (ABRIDGED), Comment: 90 pages, 25 Figures. Revised version with new authors list and references

The characterization of the dust polarization foreground to the Cosmic Microwave Background (CMB) is a necessary step towards the detection of the B-mode signal associated with primordial gravitational waves. We present a method to simulate maps of polarized dust emission on the sphere, similarly to what is done for the CMB anisotropies. This method builds on the understanding of Galactic polarization stemming from the analysis of Planck data. It relates the dust polarization sky to the structure of the Galactic magnetic field and its coupling with interstellar matter and turbulence. The Galactic magnetic field is modelled as a superposition of a mean uniform field and a random component with a power-law power spectrum of exponent $\alpha_{\rm M}$. The model parameters are constrained to fit the power spectra of dust polarization EE, BB and TE measured using Planck data. We find that the slopes of the E and B power spectra of dust polarization are matched for $\alpha_{\rm M} = -2.5$. The model allows us to compute multiple realizations of the Stokes Q and U maps for different realizations of the random component of the magnetic field, and to quantify the variance of dust polarization spectra for any given sky area outside of the Galactic plane. The simulations reproduce the scaling relation between the dust polarization power and the mean total dust intensity including the observed dispersion around the mean relation. We also propose a method to carry out multi-frequency simulations including the decorrelation measured recently by Planck, using a given covariance matrix of the polarization maps. These simulations are well suited to optimize component separation methods and to quantify the confidence with which the dust and CMB B-modes can be separated in present and future experiments. We also provide an astrophysical perspective on our modeling of the dust polarization spectra., Comment: 14 pages, 8 figures

The Planck survey has quantified polarized Galactic foregrounds and established that they are a main limiting factor in the quest for the cosmic microwave background (CMB) B-mode signal induced by primordial gravitational waves during cosmic inflation. Accurate separation of the Galactic foregrounds therefore binds this quest to our understanding of the magnetized interstellar medium (ISM). The two most relevant empirical results from analysis of Planck data are line of sight depolarization arising from fluctuations of the Galactic magnetic field orientation and alignment of filamentary dust structures with the magnetic field at high Galactic latitude. Furthermore, Planck and HI emission data in combination indicate that most of the filamentary dust structures are in the cold neutral medium. The goal of this paper is to test whether these salient observational results, taken together, can account fully for the statistical properties of the dust polarization over a selected low column density region comprising 34 % of the southern Galactic cap ($b \le -30\deg$). To do this, we construct a dust model that incorporates HI column density maps as tracers of the dust intensity structures and a phenomenological description of the Galactic magnetic field. By adjusting the parameters of the dust model, we were able to reproduce the Planck dust observations at 353 GHz in the selected region. Realistic simulations of the polarized dust emission enabled by such a dust model are useful for testing the accuracy of component separation methods, studying non-Gaussianity, and constraining the amount of decorrelation with frequency., Comment: 19 pages, 11 figures. A&A accepted

The B-mode Foreground Experiment (BFORE) is a proposed NASA balloon project designed to make optimal use of the sub-orbital platform by concentrating on three dust foreground bands (270, 350, and 600 GHz) that complement ground-based cosmic microwave background (CMB) programs. BFORE will survey ~1/4 of the sky with 1.7 - 3.7 arcminute resolution, enabling precise characterization of the Galactic dust that now limits constraints on inflation from CMB B-mode polarization measurements. In addition, BFORE's combination of frequency coverage, large survey area, and angular resolution enables science far beyond the critical goal of measuring foregrounds. BFORE will constrain the velocities of thousands of galaxy clusters, provide a new window on the cosmic infrared background, and probe magnetic fields in the interstellar medium. We review the BFORE science case, timeline, and instrument design, which is based on a compact off-axis telescope coupled to >10,000 superconducting detectors., Comment: 7 pages, 4 figures, conference proceedings published in Journal of Low Temperature Physics

PRISM (Polarized Radiation Imaging and Spectroscopy Mission) was proposed to ESA in May 2013 as a large-class mission for investigating within the framework of the ESA Cosmic Vision program a set of important scientific questions that require high resolution, high sensitivity, full-sky observations of the sky emission at wavelengths ranging from millimeter-wave to the far-infrared. PRISM's main objective is to explore the distant universe, probing cosmic history from very early times until now as well as the structures, distribution of matter, and velocity flows throughout our Hubble volume. PRISM will survey the full sky in a large number of frequency bands in both intensity and polarization and will measure the absolute spectrum of sky emission more than three orders of magnitude better than COBE FIRAS. The aim of this Extended White Paper is to provide a more detailed overview of the highlights of the new science that will be made possible by PRISM, Comment: Figure 8 replaced by the correct version, affiliations of authors added, and minor changes to match the published JCAP version

PRISM (Polarized Radiation Imaging and Spectroscopy Mission) was proposed to ESA in response to the Call for White Papers for the definition of the L2 and L3 Missions in the ESA Science Programme. PRISM would have two instruments: (1) an imager with a 3.5m mirror (cooled to 4K for high performance in the far-infrared---that is, in the Wien part of the CMB blackbody spectrum), and (2) an Fourier Transform Spectrometer (FTS) somewhat like the COBE FIRAS instrument but over three orders of magnitude more sensitive. Highlights of the new science (beyond the obvious target of B-modes from gravity waves generated during inflation) made possible by these two instruments working in tandem include: (1) the ultimate galaxy cluster survey gathering 10e6 clusters extending to large redshift and measuring their peculiar velocities and temperatures (through the kSZ effect and relativistic corrections to the classic y-distortion spectrum, respectively) (2) a detailed investigation into the nature of the cosmic infrared background (CIB) consisting of at present unresolved dusty high-z galaxies, where most of the star formation in the universe took place, (3) searching for distortions from the perfect CMB blackbody spectrum, which will probe a large number of otherwise inaccessible effects (e.g., energy release through decaying dark matter, the primordial power spectrum on very small scales where measurements today are impossible due to erasure from Silk damping and contamination from non-linear cascading of power from larger length scales). These are but a few of the highlights of the new science that will be made possible with PRISM., Comment: 20 pages Latex

CO and dust emission at millimeter wavelengths are independent tracers of cold interstellar matter, which have seldom been compared on the scale of GMCs in other galaxies. In this study, and for the first time in the Large Magellanic Cloud, we compute the molecular cloud masses from the mm emission of the dust and compare them with the masses derived from their CO luminosity and virial theorem. We present CO (J=1-0,2-1) and 1.2 mm continuum observations of the N11 star forming region in the LMC obtained with the SEST telescope and the SIMBA bolometer, respectively. We use the CO data to identify individual molecular clouds and measure their physical properties. The correlations between the properties of the N11 clouds are in agreement with those found in earlier studies in the LMC that sample a larger set of clouds and a larger range of cloud masses. For the N11 molecular clouds, we compare the masses estimated from the CO luminosity (Xco\Lco), the virial theorem (Mvir) and the millimeter dust luminosity (L_d). The measured ratios Lco/Mvir and L_d/Mvir constrain the Xco and K_d (dust emissivity at 1.2 mm per unit gas mass) parameters as a function of the virial parameter a_vir. The comparison between the different mass estimates yields a Xco-factor of 8.8x10^20 cm-2/(K km s-1) x a_vir and a K_d parameter of 1.5x10^-3 cm2/g x a_vir. We compare our N11 results with a similar analysis for molecular clouds in the Gould's Belt. We do not find in N11 a large discrepancy between the dust mm and virial masses as reported in earlier studies of molecular clouds in the SMC. The ratio between L_d and Mvir in N11 is half of that measured for Gould's Belt clouds, which can be accounted for by a factor of two lower gas-to-dust mass ratio, as the difference in gas metallicities. If the two samples have similar a_vir values, this result implies that their dust far-IR properties are also similar., Comment: A&A accepted, 13 pages, 10 figures

We report on single-dish radio CO observations towards the inter-galactic medium (IGM) of the Stephan's Quintet (SQ) group of galaxies. Extremely bright mid-IR H2 rotational line emission from warm molecular gas has been detected by Spitzer in the kpc-scale shock created by a galaxy collision. We detect in the IGM CO(1-0), (2-1) and (3-2) line emission with complex profiles, spanning a velocity range of 1000 km/s. The spectra exhibit the pre-shock recession velocities of the two colliding gas systems (5700 and 6700 km/s), but also intermediate velocities. This shows that much of the molecular gas has formed out of diffuse gas accelerated by the galaxy-tidal arm collision. A total H2 mass of 5x10^9 Msun is detected in the shock. The molecular gas carries a large fraction of the gas kinetic energy involved in the collision, meaning that this energy has not been thermalized yet. The turbulent kinetic energy of the H2 gas is at least a factor of 5 greater than the thermal energy of the hot plasma heated by the collision. The ratio between the warm H2 mass derived from Spitzer IRS spectroscopy and the H2 mass derived from CO fluxes is ~0.3 in the IGM of SQ, which is 10-100 times higher than in star-forming galaxies. In the shocked region, the ratio of the PAH-to-CO surface luminosities, commonly used to measure the star formation efficiency of the H2 gas, is lower (up to a factor 75) than the observed values in star-forming galaxies. We suggest that turbulence fed by the galaxy-tidal arm collision maintains a high heating rate within the H2 gas. This interpretation implies that the velocity dispersion on the scale of giant molecular clouds in SQ is one order of magnitude larger than the Galactic value. The high amplitude of turbulence may explain why this gas is not forming stars efficiently. [abridged version], Comment: Revised abstract and small editing to match published version. 15 pages, 5 figures. Accepted for publication in ApJ

We present an analysis of super-giant molecular complexes (SGMCs) in the overlap region of the Antennae galaxy merger, based on ALMA CO(3-2) interferometry and VLT/SINFONI imaging spectroscopy of H2 1-0 S(1) at angular resolutions of 0.9" and 0.7", respectively. All but one SGMC have multiple velocity components offset from each other by up to 150 km/s. H2 line emission is found in all SGMCs and the kinematics of H2 and CO are well matched. H2/CO line ratios vary by up to a factor of 10 among SGMCs and different velocity components of the same SGMCs. We also identify the CO counterpart of a bright, compact source of near-IR H2 line emission, which shows no Brgamma, and was first identified with SINFONI. This source has the highest H2/CO line ratio, and coincides with the steepest CO velocity gradient of the entire overlap region. With a size of 50 pc and a virial mass of a few 10^7 Msun it is perhaps a pre-cluster cloud that has not yet formed significant numbers of massive stars. We present observational evidence that the H2 emission is powered by shocks, and demonstrate how the H2 1-0 S(1) and the CO(3-2) lines can be used as tracers of energy dissipation and gas mass, respectively. The variations in the H2/CO line ratio may indicate that the SGMCs are dissipating their turbulent kinetic energy at different rates. The compact source could represent a short (~ 1 Myr) evolutionary stage in the early formation of super-star clusters., Comment: 4 pages, 3 figures, A&A accepted: december 28th 2011

We report the discovery of a 3 kpc disk of few 10^9 Ms of dense, warm H_2 in the nearby radio galaxy 3C326 N, which shows no signs of on-going or recent star formation and falls a factor 60 below the Schmidt-Kennicutt law. VLT/SINFONI imaging spectroscopy shows broad (FWHM \sim 500 km/s) ro-vibrational H_2 lines across all of the disk, with irregular profiles and line ratios consistent with shocks. The ratio of turbulent and gravitational energy suggests that the gas is highly turbulent and not gravitationally bound. In absence of the driving by the jet, short turbulent dissipation times suggest the gas should collapse rapidly and form stars, at odds with the recent star-formation history. Motivated by hydrodynamic models of rapid H_2 formation boosted by turbulent compression, we propose that the molecules formed from diffuse atomic gas in the turbulent jet cocoon. Since the gas is not self-gravitating, it cannot form molecular clouds or stars while the jet is active, and is likely to disperse and become atomic again after the nuclear activity ceases. We speculate that very low star-formation rates are to be expected under such conditions, provided that the large-scale turbulence sets the gas dynamics in molecular clouds. Our results illustrate that jets may create large molecular reservoirs as expected in 'positive feedback' scenarios of AGN-triggered star formation, but that this alone is not sufficient to trigger star formation., Comment: A&A accepted

H2 pure-rotational emission lines are detected from warm (100-1500 K) molecular gas in 17/55 (31% of) radio galaxies at redshift z<0.22 observed with the Spitzer IR Spectrograph. The summed H2 0-0 S(0)-S(3) line luminosities are L(H2)=7E38-2E42 erg/s, yielding warm H2 masses up to 2E10 Msun. These radio galaxies, of both FR radio morphological types, help to firmly establish the new class of radio-selected molecular hydrogen emission galaxies (radio MOHEGs). MOHEGs have extremely large H2 to 7.7 micron PAH emission ratios: L(H2)/L(PAH7.7) = 0.04-4, up to a factor 300 greater than the median value for normal star-forming galaxies. In spite of large H2 masses, MOHEGs appear to be inefficient at forming stars, perhaps because the molecular gas is kinematically unsettled and turbulent. Low-luminosity mid-IR continuum emission together with low-ionization emission line spectra indicate low-luminosity AGNs in all but 3 radio MOHEGs. The AGN X-ray emission measured with Chandra is not luminous enough to power the H2 emission from MOHEGs. Nearly all radio MOHEGs belong to clusters or close pairs, including 4 cool core clusters (Perseus, Hydra, A 2052, and A 2199). We suggest that the H2 in radio MOHEGs is delivered in galaxy collisions or cooling flows, then heated by radio jet feedback in the form of kinetic energy dissipation by shocks or cosmic rays., Comment: ApJ in press, 40 pages, 18 figures, 14 tables

We want to study the mid-infrared properties and the starburst and AGN contributions, of 24um sources at z~2, through analysis of mid-infrared spectra combined with millimeter, radio, and infrared photometry. Mid-infrared spectroscopy allows us to recover accurate redshifts. A complete sample of 16 Spitzer-selected sources (ULIRGs) believed to be starbursts at z~2 ("5.8um-peakers") was selected in the (0.5 sq.deg.) J1064+56 SWIRE Lockman Hole field. These sources have S(24um)>0.5mJy, a stellar emission peak redshifted to 5.8um, and r'(Vega)>23. The entire sample was observed with the low resolution units of the Spitzer/IRS infrared spectrograph. These sources have 1.2mm observations with IRAM 30m/MAMBO and very deep 20cm observations from the VLA. Nine of our sources also benefit from 350um observation and detection from CSO/SHARC-II. The entire sample shows good quality IRS spectra dominated by strong PAH features. The main PAH features at 6.2, 7.7, 8.6, and 11.3um have high S/N average luminosities of 2.90, 10.38, 3.62, and 2.29x10^{10}Lsun, respectively. We derived accurate redshifts spanning from 1.75 to 2.28. The average of these redshifts is 2.017. This result confirms that the selection criteria of "5.8um-peakers" associated with a strong detection at 24um are reliable to select sources at z~2. We have analyzed the different correlations between PAH emission and infrared, millimeter, and radio emission. Practically all our sources are strongly dominated by starburst emission. We have also defined two subsamples based on the equivalent width at 7.7um to investigate AGN contributions. Our sample contains strong starbursts and represents a particularly 24um-bright class of SMGs. The very good correlation between PAH and far-IR luminosities is now confirmed in high-z starburst ULIRGs. These sources show a small AGN contribution to the mid-IR, around ~20% in most cases., Comment: 19 pages, 14 figures + 2 figure in appendix Accepted in A&A on 26/08/201

The amount of molecular gas is a key for understanding the future star formation in a galaxy. Because H2 is difficult to observe directly in dense and cold clouds, tracers like CO are used. However, at low metallicities especially, CO only traces the shielded interiors of the clouds. mm dust emission can be used as a tracer to unveil the total dense gas masses. The comparison of masses deduced from the continuum SIMBA 1.2 mm emission and virial masses in a sample of giant molecular clouds (GMCs), in the SW region of the Small Magellanic Cloud (SMC), showed a discrepancy that is in need of an explanation. This study aims at better assessing possible uncertainties on the dust emission observed in the sample of GMCs from the SMC and focuses on the densest parts of the GMCs where CO is detected. New observations were obtained with the LABOCA camera on the APEX telescope. All GMCs previously observed in CO are detected and their emission at 870microns is compared to ancillary data. The different contributions to the sub-mmm emission are estimated, as well as dust properties, in order to deduce molecular cloud masses precisely. The (sub-)mm emission observed in the GMCs in the SW region of the SMC is dominated by dust emission and masses are deduced for the part of each cloud where CO is detected and compared to the virial masses. The mass discrepancy between both methods is confirmed at 870microns with the LABOCA observations: the virial masses are on average 4 times smaller than the masses of dense gas deduced from dust emission, contrary to what is observed for equivalent clouds in our Galaxy. The origin of this mass discrepancy in the SMC remains unkown. The direct interpretation of this effect is that the CO linewidth used to compute virial masses do not measure the full velocity distribution of the gas. Geometrical effects and uncertainties on the dust properties are also discussed., Comment: 12 pages, 4 figures, Accepted in A&A