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

Mid-infrared emission features probe the properties of ionized gas, and hot or warm molecular gas. The Orion Bar is a frequently studied photodissociation region (PDR) containing large amounts of gas under these conditions, and was observed with the MIRI IFU aboard JWST as part of the "PDRs4All" program. The resulting IR spectroscopic images of high angular resolution (0.2") reveal a rich observational inventory of mid-IR emission lines, and spatially resolve the substructure of the PDR, with a mosaic cutting perpendicularly across the ionization front and three dissociation fronts. We extracted five spectra that represent the ionized, atomic, and molecular gas layers, and measured the most prominent gas emission lines. An initial analysis summarizes the physical conditions of the gas and the potential of these data. We identified around 100 lines, report an additional 18 lines that remain unidentified, and measured the line intensities and central wavelengths. The H I recombination lines originating from the ionized gas layer bordering the PDR, have intensity ratios that are well matched by emissivity coefficients from H recombination theory, but deviate up to 10% due contamination by He I lines. We report the observed emission lines of various ionization stages of Ne, P, S, Cl, Ar, Fe, and Ni, and show how certain line ratios vary between the five regions. We observe the pure-rotational H$_2$ lines in the vibrational ground state from 0-0 S(1) to 0-0 S(8), and in the first vibrationally excited state from 1-1 S(5) to 1-1 S(9). We derive H$_2$ excitation diagrams, and approximate the excitation with one thermal (~700 K) component representative of an average gas temperature, and one non-thermal component (~2700 K) probing the effect of UV pumping. We compare these results to an existing model for the Orion Bar PDR and highlight the differences with the observations., Comment: 26 pages, 12 figures, 3 tables. Submitted to A&A, under review (1st revision)

Most low-mass stars form in stellar clusters that also contain massive stars, which are sources of far-ultraviolet (FUV) radiation. Theoretical models predict that this FUV radiation produces photo-dissociation regions (PDRs) on the surfaces of protoplanetary disks around low-mass stars, impacting planet formation within the disks. We report JWST and Atacama Large Millimetere Array observations of a FUV-irradiated protoplanetary disk in the Orion Nebula. Emission lines are detected from the PDR; modelling their kinematics and excitation allows us to constrain the physical conditions within the gas. We quantify the mass-loss rate induced by the FUV irradiation, finding it is sufficient to remove gas from the disk in less than a million years. This is rapid enough to affect giant planet formation in the disk.

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

(Abridged) We investigate the impact of radiative feedback from massive stars on their natal cloud and focus on the transition from the HII region to the atomic PDR (crossing the ionisation front (IF)), and the subsequent transition to the molecular PDR (crossing the dissociation front (DF)). We use high-resolution near-IR integral field spectroscopic data from NIRSpec on JWST to observe the Orion Bar PDR as part of the PDRs4All JWST Early Release Science Program. The NIRSpec data reveal a forest of lines including, but not limited to, HeI, HI, and CI recombination lines, ionic lines, OI and NI fluorescence lines, Aromatic Infrared Bands (AIBs including aromatic CH, aliphatic CH, and their CD counterparts), CO2 ice, pure rotational and ro-vibrational lines from H2, and ro-vibrational lines HD, CO, and CH+, most of them detected for the first time towards a PDR. Their spatial distribution resolves the H and He ionisation structure in the Huygens region, gives insight into the geometry of the Bar, and confirms the large-scale stratification of PDRs. We observe numerous smaller scale structures whose typical size decreases with distance from Ori C and IR lines from CI, if solely arising from radiative recombination and cascade, reveal very high gas temperatures consistent with the hot irradiated surface of small-scale dense clumps deep inside the PDR. The H2 lines reveal multiple, prominent filaments which exhibit different characteristics. This leaves the impression of a "terraced" transition from the predominantly atomic surface region to the CO-rich molecular zone deeper in. This study showcases the discovery space created by JWST to further our understanding of the impact radiation from young stars has on their natal molecular cloud and proto-planetary disk, which touches on star- and planet formation as well as galaxy evolution., Comment: 52 pages, 30 figures, submitted to A&A

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 JWST has captured the most detailed and sharpest infrared images ever taken of the inner region of the Orion Nebula, the nearest massive star formation region, and a prototypical highly irradiated dense photo-dissociation region (PDR). We investigate the fundamental interaction of far-ultraviolet photons with molecular clouds. The transitions across the ionization front (IF), dissociation front (DF), and the molecular cloud are studied at high-angular resolution. These transitions are relevant to understanding the effects of radiative feedback from massive stars and the dominant physical and chemical processes that lead to the IR emission that JWST will detect in many Galactic and extragalactic environments. Due to the proximity of the Orion Nebula and the unprecedented angular resolution of JWST, these data reveal that the molecular cloud borders are hyper structured at small angular scales of 0.1-1" (0.0002-0.002 pc or 40-400 au at 414 pc). A diverse set of features are observed such as ridges, waves, globules and photoevaporated protoplanetary disks. At the PDR atomic to molecular transition, several bright features are detected that are associated with the highly irradiated surroundings of the dense molecular condensations and embedded young star. Toward the Orion Bar PDR, a highly sculpted interface is detected with sharp edges and density increases near the IF and DF. This was predicted by previous modeling studies, but the fronts were unresolved in most tracers. A complex, structured, and folded DF surface was traced by the H2 lines. This dataset was used to revisit the commonly adopted 2D PDR structure of the Orion Bar. JWST provides us with a complete view of the PDR, all the way from the PDR edge to the substructured dense region, and this allowed us to determine, in detail, where the emission of the atomic and molecular lines, aromatic bands, and dust originate.

(Abridged) Mid-infrared observations of photodissociation regions (PDRs) are dominated by strong emission features called aromatic infrared bands (AIBs). The most prominent AIBs are found at 3.3, 6.2, 7.7, 8.6, and 11.2 $\mu$m. The most sensitive, highest-resolution infrared spectral imaging data ever taken of the prototypical PDR, the Orion Bar, have been captured by JWST. We provide an inventory of the AIBs found in the Orion Bar, along with mid-IR template spectra from five distinct regions in the Bar: the molecular PDR, the atomic PDR, and the HII region. We use JWST NIRSpec IFU and MIRI MRS observations of the Orion Bar from the JWST Early Release Science Program, PDRs4All (ID: 1288). We extract five template spectra to represent the morphology and environment of the Orion Bar PDR. The superb sensitivity and the spectral and spatial resolution of these JWST observations reveal many details of the AIB emission and enable an improved characterization of their detailed profile shapes and sub-components. While the spectra are dominated by the well-known AIBs at 3.3, 6.2, 7.7, 8.6, 11.2, and 12.7 $\mu$m, a wealth of weaker features and sub-components are present. We report trends in the widths and relative strengths of AIBs across the five template spectra. These trends yield valuable insight into the photochemical evolution of PAHs, such as the evolution responsible for the shift of 11.2 $\mu$m AIB emission from class B$_{11.2}$ in the molecular PDR to class A$_{11.2}$ in the PDR surface layers. This photochemical evolution is driven by the increased importance of FUV processing in the PDR surface layers, resulting in a "weeding out" of the weakest links of the PAH family in these layers. For now, these JWST observations are consistent with a model in which the underlying PAH family is composed of a few species: the so-called 'grandPAHs'., Comment: 25 pages, 10 figures, to appear in A&A

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

Massive stars disrupt their natal molecular cloud material through radiative and mechanical feedback processes. These processes have profound effects on the evolution of interstellar matter in our Galaxy and throughout the Universe, from the era of vigorous star formation at redshifts of 1-3 to the present day. The dominant feedback processes can be probed by observations of the Photo-Dissociation Regions (PDRs) where the far-ultraviolet photons of massive stars create warm regions of gas and dust in the neutral atomic and molecular gas. PDR emission provides a unique tool to study in detail the physical and chemical processes that are relevant for most of the mass in inter- and circumstellar media including diffuse clouds, proto-planetary disks and molecular cloud surfaces, globules, planetary nebulae, and star-forming regions. PDR emission dominates the infrared (IR) spectra of star-forming galaxies. Most of the Galactic and extragalactic observations obtained with the James Webb Space Telescope (JWST) will therefore arise in PDR emission. In this paper we present an Early Release Science program using the MIRI, NIRSpec, and NIRCam instruments dedicated to the observations of an emblematic and nearby PDR: the Orion Bar. These early JWST observations will provide template datasets designed to identify key PDR characteristics in JWST observations. These data will serve to benchmark PDR models and extend them into the JWST era. We also present the Science-Enabling products that we will provide to the community. These template datasets and Science-Enabling products will guide the preparation of future proposals on star-forming regions in our Galaxy and beyond and will facilitate data analysis and interpretation of forthcoming JWST observations., Comment: Submitted to PASP

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

Context. Planck observations demonstrated that the grain alignment efficiency is almost constant in the diffuse ISM. Aims. We test if the Radiative Torque (RAT) theory is compatible with observational constraints on grain alignment. Methods. We combine a numerical simulation with the radiative transfer code POLARIS that incorporates a physical dust model and the detailed grain alignment physics of RATs. A dust model is designed to reproduce the spectral dependence of extinction of the ISM. From a RAMSES simulation of interstellar turbulence, we extract a cube representative of the diffuse ISM. We post-process the cube with POLARIS to get the grain temperature and RATs to simulate synthetic dust polarization maps. Results. In our simulation the grain alignment efficiency is correlated with gas pressure, but not with the RAT intensity. Because of the low dust extinction, the magnitude of RATs varies little, decreasing only for high column densities $N_H$. Comparing our maps with a uniform alignment efficiency, we find no systematic difference. The dependence of polarization fraction $p$ with $N_H$ or polarization dispersion $S$ is similar. The drop of RATs in dense regions barely affects the polarization pattern, the signal being dominated by the LOS and magnetic field geometry. If a star is inserted, the polarization increases, with no specific pattern around the star. The angle-dependence of RATs is not observed in the maps, and is weak using a uniform magnetic field. Conclusions. RATs are compatible with Planck data for the diffuse ISM such that both uniform alignment and RAT alignment lead to similar observations. To further test the predictions of RATs where an important drop of grain alignment is expected, polarization observations of dense regions must be confronted to numerical simulations sampling high column densities through dense clouds, with enough statistics., Comment: 22 pages, 22 figures

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 primary source of systematic uncertainty in the quest for the B-mode polarization of the Cosmic Microwave Background (CMB) introduced by primordial gravitational waves is polarized thermal emission from Galactic dust. Therefore, accurate characterization and separation of the polarized thermal dust emission is an essential step in distinguishing such a faint CMB B-mode signal. We provide a modelling framework to simulate polarized thermal dust emission based on the model described in Ghosh et al. (2017), making use of both the Planck dust and Effelsberg-Bonn HI surveys over the northern Galactic cap. Our seven-parameter dust model, incorporating both HI gas in three different column density templates as a proxy for spatially variable dust intensity and a phenomenological model of Galactic magnetic field, is able to reproduce both 1- and 2-point statistics of the observed dust polarization maps seen by Planck at 353 GHz over a selected low-column density region in the northern Galactic cap. This work has important applications in assessing the accuracy of component separation methods and in quantifying the confidence level of separating polarized Galactic emission and the CMB B-mode signal, as is needed for ongoing and future CMB missions., Comment: 12 pages, 12 figures, Accepted in A&A

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

Magnetic fields have important or dominant effects in many areas of astrophysics, but have been very difficult to quantify. Spectropolarimetry from Ground State Alignment (GSA) has been suggested as a direct tracer of magnetic field in interstellar diffuse medium. Owing to the long life of the atoms on ground states the Larmor precession in an external magnetic field imprints the direction of the field onto the polarization of absorbing species. This provides a unique tool for studies of sub-gauss magnetic fields using polarimetry of UV, optical and radio lines. Many spectral lines with strong signals from GSA are in the UV band. By discerning magnetic fields in gas with different dynamical properties, high spectral resolution measurement of spectral polarization will allow the study of 3D magnetic field distribution and interstellar turbulence. GSA provides also a unique chance to map 3D direction of magnetic field on small scales, e.g., disks, where grain alignment is unreliable. The range of objects suitable for studies is extremely wide and includes magnetic fields in the interplanetary medium, in the interstellar medium, and in circumstellar regions as well as diffuse media in extragalactic objects., Comment: 5 pages, 1 figure; science white paper submitted to the Astro2020 decadal survey. arXiv admin note: substantial text overlap with arXiv:1203.5571, arXiv:1302.3264

Understanding the physics of how stars form is a highly-prioritized goal of modern Astrophysics, in part because star formation is linked to both galactic dynamics on large scales and to the formation of planets on small scales. It is well-known that stars form from the gravitational collapse of molecular clouds, which are in turn formed out of the turbulent interstellar medium. Star formation is highly inefficient, with one of the likely culprits being the regulation against gravitational collapse provided by magnetic fields. Measurement of the polarized emission from interstellar dust grains, which are partially aligned with the magnetic field, provides a key tool for understanding the role these fields play in the star formation process. Over the past decade, much progress has been made by the most recent generation of polarimeters operating over a range of wavelengths (from the far-infrared through the millimeter part of the spectrum) and over a range of angular resolutions (from less than an arcsecond through fractions of a degree). Future developments in instrument sensitivity for ground-based, airborne, and space-borne polarimeters operating over range of spatial scales are critical for enabling revolutionary steps forward in our understanding of the magnetized turbulence from which stars are formed., Comment: Submitted to the Astro2020 Decadal Survey (with a few updated citations)

Galaxy formation depends on a complex interplay between gravitational collapse, gas accretion, merging, and feedback processes. Yet, after many decades of investigation, these concepts are poorly understood. This paper presents the argument that warm H$_2$ can be used as a tool to unlock some of these mysteries. Turbulence, shocks and outflows, driven by star formation, AGN activity or inflows, may prevent the rapid buildup of star formation in galaxies. Central to our understanding of how gas is converted into stars is the process by which gas can dissipate its mechanical energy through turbulence and shocks in order to cool. H$_2$ lines provide direct quantitative measurements of kinetic energy dissipation in molecular gas in galaxies throughout the Universe. Based on the detection of very powerful H$_2$ lines from z = 2 galaxies and proto-clusters at the detection limits of {\it Spitzer}, we are confident that future far-IR and UV H$_2$ observations will provide a wealth of new information and insight into galaxy evolution to high-z. Finally, at the very earliest epoch of star and galaxy formation, warm H$_2$ may also provide a unique glimpse of molecular gas collapse at 7 $<$ z $<$ 12 in massive dark matter (DM) halos on their way to forming the very first galaxies. Such measurements are beyond the reach of existing and planned observatories., Comment: Submitted as a science White Paper to the Astronomy and Astrophysics Astro 2020 Decadal Survey call issued by the National Academies of Sciences, Engineering and Medicine (March 11 2019)

A comprehensive study of interstellar medium phases and the nature of their boundaries or connections requires comparing abundances and velocity profiles from tracers of the different phases. In the UV, studies of a wealth of absorption features appearing in the spectra of hot stars yield fundamental insights into the composition and physical characteristics of all phases of the ISM along with the processes that influence them. They also inform us on the nature of boundaries between them. However no single instrument has as yet given access to species in all ISM phases at the same high spectral resolution: from the molecular bands of CO and H2 in the far-UV, to the cold neutral medium tracers C I and S I and the warm medium tracers like C II, NI, OI, Mg II, Fe II, Si II etc, and to the high ions of the hot ionized medium C IV, Si IV in the UV, as well as O VI in the far UV. We have yet to design the spectrometer that will enable observing the full UV domain at resolving power R>200 000 and signal-to-noise ratio SNR>500. The line FWHM being governed by turbulence, temperature, and species mass, such a resolution is necessary to resolve lines from both the cold molecular hydrogen and the warm metal ions with a turbulence of ~1 km/s, and to differentiate distinct velocity components, typically separated by less than 2 km/s., Comment: Science white paper submitted to the Astro2020 Decadal Survey

The IMAGINE Consortium aims to bring modeling of the magnetic field of the Milky Way to a next level, by using Bayesian inference. IMAGINE includes an open-source modular software pipeline that optimizes parameters in a user-defined Galactic magnetic field model against various selected observational datasets. Bayesian priors can be added as external probabilistic constraints of the model parameters. These conference proceedings describe the science goals of the IMAGINE Consortium, the software pipeline and its inputs, viz observational data sets, Galactic magnetic field models, and Bayesian priors., Comment: Accepted to Galaxies Special Issue "The Power of Faraday Tomography", proceedings of workshop of the same name in Miyazaki, Japan, May 2018

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

In this white paper we introduce the IMAGINE Consortium and its scientific background, goals and structure. Our purpose is to coordinate and facilitate the efforts of a diverse group of researchers in the broad areas of the interstellar medium, Galactic magnetic fields and cosmic rays, and our goal is to develop more comprehensive insights into the structures and roles of interstellar magnetic fields and their interactions with cosmic rays. To achieve a higher level of self-consistency, depth and rigour can only be achieved by the coordinated efforts of experts in diverse areas of astrophysics involved in observational, theoretical and numerical work. We present our view of the present status of this topic, identify its key unsolved problems and suggest a strategy that will underpin our work. The backbone of the consortium is the Interstellar MAGnetic field INference Engine, a publicly available Bayesian platform that employs robust statistical methods to explore the multi-dimensional likelihood space using any number of modular inputs. It provides an interpretation and modelling framework that has the power and flexibility to include a variety of observational, theoretical and numerical lines of evidence into a self-consistent and comprehensive picture of the thermal and non-thermal interstellar media. An important innovation is that a consistent understanding of the phenomena that are directly or indirectly influenced by the Galactic magnetic field, such as the deflection of ultra-high energy cosmic rays or extragalactic backgrounds, is made an integral part of the modelling. The IMAGINE Consortium, which is informal by nature and open to new participants, hereby presents a methodological framework for the modelling and understanding of Galactic magnetic fields that is available to all communities whose research relies on a state-of-the-art solution to this problem. (Abridged.), Comment: 56 pages, 15 figures, minor updates to match published JCAP version

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

Polarized extinction and emission from dust in the interstellar medium (ISM) are hard to interpret, as they have a complex dependence on dust optical properties, grain alignment and magnetic field orientation. This is particularly true in molecular clouds. The data available today are not yet used to their full potential. The combination of emission and extinction, in particular, provides information not available from either of them alone. We combine data from the scientific literature on polarized dust extinction with Planck data on polarized emission and we use them to constrain the possible variations in dust and environmental conditions inside molecular clouds, and especially translucent lines of sight, taking into account magnetic field orientation. We focus on the dependence between \lambda_max -- the wavelength of maximum polarization in extinction -- and other observables such as the extinction polarization, the emission polarization and the ratio of the two. We set out to reproduce these correlations using Monte-Carlo simulations where the relevant quantities in a dust model -- grain alignment, size distribution and magnetic field orientation -- vary to mimic the diverse conditions expected inside molecular clouds. None of the quantities chosen can explain the observational data on its own: the best results are obtained when all quantities vary significantly across and within clouds. However, some of the data -- most notably the stars with low emission-to-extinction polarization ratio -- are not reproduced by our simulation. Our results suggest not only that dust evolution is necessary to explain polarization in molecular clouds, but that a simple change in size distribution is not sufficient to explain the data, and point the way for future and more sophisticated models.

Super star clusters (SSCs), likely the progenitors of globular clusters, are one of the most extreme forms of star formation. Understanding how SSCs form is an observational challenge. Theoretical studies establish that, to form such clusters, the dynamical timescale of their parent clouds has to be shorter than the timescale of the disruption of their parent clouds by stellar feedback. However, due to insufficient observational support, it is still unclear how feedback from SSCs acts on their surrounding matter. Studying feedback in SSCs is essential to understand how such clusters form. Based on ALMA and VLT observations, we study this process in a SSC in the overlap region of the Antennae galaxies. We analyze a unique massive (~10^7 Msun) and young (1-3.5 Myr) SSC, still associated with compact molecular and ionized gas emission. The cluster has two CO velocity components, a low velocity one spatially associated with the cluster and a high velocity one distributed in a bubble-like shape around the cluster. Our results on the low velocity component suggest that this gas did not participate in the formation of the SSC. We propose that most of the parent cloud has already been blown away, accelerated at the early stages of the SSC evolution by radiation pressure, in a timescale ~1 Myr. The high velocity component may trace outflowing molecular gas from the parent cloud. Supporting evidence is found in shock heated H2 gas. The low velocity component may be gas that was near the SSC when it formed but not part of its parent cloud or clumps that migrated from the SGMC environment. This gas would be dispersed by stellar winds and supernova explosions. Within ~100 pc from the cluster, we estimate SFE>17%, smaller than the theoretical limit of 30% needed to form a bound cluster. Further higher spatial resolution observations are needed to test and quantify our proposed scenario., Comment: 15 pages, accepted to A&A

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

Simulations of cosmological filamentary accretion reveal flows ("streams") of warm gas, ~$10^4$ K, which are efficient in bringing gas into galaxies. We present a phenomenological scenario where gas in such flows -- if it is shocked as it enters the halo as we assume -- become biphasic and, as a result, turbulent. We consider a collimated stream of warm gas that flows into a halo from an over dense filament of the cosmic web. The post-shock streaming gas expands because it has a higher pressure than the ambient halo gas, and fragments as it cools. The fragmented stream forms a two phase medium: a warm cloudy phase embedded in hot post-shock gas. We argue that the hot phase sustains the accretion shock. A fraction of the initial kinetic energy of the infalling gas is converted into turbulence among and within the warm clouds. The thermodynamic evolution of the post-shock gas is largely determined by the relative timescales of several processes -- the cooling, the expansion of the post-shock gas, the amount of turbulence in the clouds, and the halo dynamics. We expect the gas to become multiphase when the cooling and dynamical times are of the same order-of-magnitude. In this framework, we show that this occurs in the important mass range of ~$10^{11}$ to $10^{13}$ M$_\odot$ , where the bulk of stars have formed in galaxies. Gas accreting along cosmic web filaments may eventually lose coherence and mix with the ambient halo gas. Through both the phase separation and "disruption" of the stream, the accretion efficiency onto a galaxy in a halo dynamical time is lowered. De-collimating flows make the direct interaction between galaxy feedback and accretion streams more likely, thereby further reducing the overall accretion efficiency. Moderating the gas accretion efficiency through these mechanisms may help to alleviate a number of significant challenges in theoretical galaxy formation., Comment: 13 pages, 8 figures, Submitted to A&A, New version includes new figures

This study aims to elucidate a possible link between chemical properties of ices in star-forming regions and environmental characteristics of the host galaxy. We performed 3--4 micron spectroscopic observations toward nine embedded high-mass YSOs in the Large Magellanic Cloud (LMC) with the ISAAC at the VLT. Additionally, we analyzed archival ISAAC data of two LMC YSOs. As a result, we detected absorption bands due to solid H2O and CH3OH as well as the 3.47 micron absorption band. The 3.53 micron CH3OH ice absorption band for the LMC YSOs is found to be absent or very weak compared to that seen toward Galactic sources. The result suggests the low abundance of CH3OH ice in the LMC. The 3.47 micron absorption band is detected toward six out of eleven LMC YSOs. We found that the 3.47 micron band and the H2O ice band correlate similarly between the LMC and Galactic samples, but the LMC sources seem to require a slightly higher H2O ice threshold for the appearance of the 3.47 micron band. For the LMC sources with relatively large H2O ice optical depths, we found that the strength ratio of the 3.47 micron band relative to the water ice band is only marginally lower than those of the Galactic sources. We propose that grain surface reactions at a relatively high dust temperature (warm ice chemistry) are responsible for the observed characteristics of ice chemical compositions in the LMC. We suggest that this warm ice chemistry is one of the important characteristics of interstellar and circumstellar chemistry in low metallicity environments. The low abundance of CH3OH in the solid phase implies that formation of complex organic molecules from methanol-derived species is less efficient in the LMC. For the 3.47 micron band, the observed difference in the water ice threshold may suggest that a more shielded environment is necessary for the formation of the 3.47 micron band carrier in the LMC., Comment: 21 pages, 13 figures, Accepted for publication in A&A

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

We present the first X-ray observations of the Taffy galaxies (UGC 12914/5) with the Chandra observatory, and detect soft X-ray emission in the region of the gas-rich, radio-continuum-emitting Taffy bridge. The results are compared to Herschel observations of dust and diffuse [CII] line-emitting gas. The diffuse component of the Taffy bridge has an X-ray luminosity of L(X) (0.5-8keV) =5.4 x 10^39 erg s^-1, which accounts for 19% of the luminosity of the sum for the two galaxies. The total mass in hot gas is (0.8--1.3) x 10^8 M_sun, which is approximately 1% of the total (HI~+~H2) gas mass in the bridge, and ~11% of the mass of warm molecular hydrogen discovered by Spitzer. The soft X-ray and dense CO-emitting gas in the bridge have offset distributions, with the X-rays peaking along the north-western side of the bridge in the region where stronger far-IR dust and diffuse [CII] gas is observed by Herschel. We detect nine Ultra Luminous X-ray sources (ULXs) in the system, the brightest of which is found in the bridge, associated with an extragalactic HII region. We suggest that the X-ray--emitting gas has been shocked--heated to high temperatures and "splashed" into the bridge by the collision. The large amount of gas stripped from the galaxies into the bridge and its very long gas depletion timescale (>10 Gyr) may explain why this system, unlike most major mergers, is not a powerful IR emitter., Comment: Accepted for Publication in ApJ, Sept 15 2015

We compare the performance of several dust models in reproducing the dust spectral energy distribution (SED) per unit extinction in the diffuse interstellar medium (ISM). We use our results to constrain the variability of the optical properties of big grains in the diffuse ISM, as published by the Planck collaboration. We use two different techniques to compare the predictions of dust models to data from the Planck HFI, IRAS and SDSS surveys. First, we fit the far-infrared emission spectrum to recover the dust extinction and the intensity of the interstellar radiation field (ISRF). Second, we infer the ISRF intensity from the total power emitted by dust per unit extinction, and then predict the emission spectrum. In both cases, we test the ability of the models to reproduce dust emission and extinction at the same time. We identify two issues. Not all models can reproduce the average dust emission per unit extinction: there are differences of up to a factor $\sim2$ between models, and the best accord between model and observation is obtained with the more emissive grains derived from recent laboratory data on silicates and amorphous carbons. All models fail to reproduce the variations in the emission per unit extinction if the only variable parameter is the ISRF intensity: this confirms that the optical properties of dust are indeed variable in the diffuse ISM. Diffuse ISM observations are consistent with a scenario where both ISRF intensity and dust optical properties vary. The ratio of the far-infrared opacity to the $V$ band extinction cross-section presents variations of the order of $\sim20\%$ ($40-50\%$ in extreme cases), while ISRF intensity varies by $\sim30\%$ ($\sim60\%$ in extreme cases). This must be accounted for in future modelling., Comment: A&A, in press

Given the current and past star-formation in the Milky Way in combination with the limited gas supply, the re-fuelling of the reservoir of cool gas is an important aspect of Galactic astrophysics. The infall of \ion{H}{i} halo clouds can, among other mechanisms, contribute to solving this problem. We study the intermediate-velocity cloud IVC135+54 and its spatially associated high-velocity counterpart to look for signs of a past or ongoing interaction. Using the Effelsberg-Bonn \ion{H}{i} Survey data, we investigated the interplay of gas at different velocities. In combination with far-infrared Planck and IRIS data, we extended this study to interstellar dust and used the correlation of the data sets to infer information on the dark gas. The velocity structure indicates a strong compression and deceleration of the infalling high-velocity cloud (HVC), associated with far-infrared excess emission in the intermediate-velocity cloud. This excess emission traces molecular hydrogen, confirming that IVC135+54 is one of the very few molecular halo clouds. The high dust emissivity of IVC135+54 with respect to the local gas implies that it consists of disk material and does not, unlike the HVC, have an extragalactic origin. Based on the velocity structure of the HVC and the dust content of the IVC, a physical connection between them appears to be the logical conclusion. Since this is not compatible with the distance difference between the two objects, we conclude that this particular HVC might be much closer to us than complex C. Alternatively, the indicators for an interaction are misleading and have another origin., Comment: 11 pages, 10 figures, accepted for publication in A&A

Towards the high galactic latitude sky, the far-infrared (FIR) intensity is tightly correlated to the total hydrogen column density which is made up of atomic (HI) and molecular hydrogen (H$_{2})$. Above a certain column density threshold, atomic hydrogen turns molecular. We analyse gas and dust properties of intermediate-velocity clouds (IVCs) in the lower galactic halo to explore their transition from the atomic to the molecular phase. Driven by observations, we investigate the physical processes that transform a purely atomic IVC into a molecular one. Data from the Effelsberg-Bonn HI-Survey (EBHIS) are correlated to FIR wavebands of the Planck satellite and IRIS. Modified black-body emission spectra are fitted to deduce dust optical depths and grain temperatures. We remove the contribution of atomic hydrogen to the FIR intensity to estimate molecular hydrogen column densities. Two IVCs show different FIR properties, despite their similarity in HI, such as narrow spectral lines and large column densities. One FIR bright IVC is associated with H$_{2}$, confirmed by $^{12}$CO $(1\rightarrow0)$ emission; the other IVC is FIR dim and shows no FIR excess, which indicates the absence of molecular hydrogen. We propose that the FIR dim and bright IVCs probe the transition between the atomic and molecular gas phase. Triggered by dynamical processes, this transition happens during the descent of IVCs onto the galactic disk. The most natural driver is ram pressure exerted onto the cloud by the increasing halo density. Because of the enhanced pressure, the formation timescale of H$_{2}$ is reduced, allowing the formation of large amounts of H$_{2}$ within a few Myr., Comment: 13 pages, 14 figures, accepted for publication by A&A

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 examine the chemical and emission properties of mildly irradiated (G0=1) magnetised shocks in diffuse media (nH=10^2 to 10^4 /cm3) at low to moderate velocities (from 3 to 40 km/s). Results: The formation of some molecules relies on endoergic reactions. In J-shocks, their abundances are enhanced by several orders of magnitude for shock velocities as low as 7 km/s. Otherwise most chemical properties of J-type shocks vary over less than an order of magnitude between velocities from about 7 to about 30 km/s, where H2 dissociation sets in. C-type shocks display a more gradual molecular enhancement as the shock velocity increases. We quantify the energy flux budget (fluxes of kinetic, radiated and magnetic energies) with emphasis on the main cooling lines of the cold interstellar medium. Their sensitivity to shock velocity is such that it allows observations to constrain statistical distributions of shock velocities. We fit various probability distribution functions (PDFs) of shock velocities to spectroscopic observations of the galaxy-wide shock in Stephan's Quintet (SQ) and of a Galactic line of sight sampling diffuse molecular gas in Chamaeleon. In both cases, low velocities bear the greatest statistical weight and the PDF is consistent with a bimodal distribution. In the very low velocity shocks (below 5 km/s), dissipation is due to ion-neutral friction which powers H2 low energy transitions and atomic lines. In moderate velocity shocks (20 km/s and above), the dissipation is due to viscous heating and accounts for most of the molecular emission. In our interpretation a significant fraction of the gas on the line of sight is shocked (from 4% to 66%). For example, C+ emission may trace shocks in UV irradiated gas where C+ is the dominant carbon species., Comment: In press, version prior to language edition, 15 pages, 19 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