Stellar-reddening-based Extinction Maps for Cosmological Applications

The cosmic infrared background (CIB) is a powerful probe of large-scale structure across a very large redshift range, and consists of unresolved redshifted infrared emission from dusty galaxies. It can be used to study the astrophysics of galaxies, the star formation history of the universe, and the connection between dark and luminous matter. It can furthermore be used as a tracer of the large-scale structure and thus assist in de-lensing of the cosmic microwave background. The major difficulty in its use lies in obtaining accurate and unbiased large-scale CIB images that are cleaned of the contamination by Galactic dust. We used data on neutral atomic hydrogen from the recently released HI4PI Survey to create template maps of Galactic dust, allowing us to remove this component from the Planck intensity maps from 353 to 857 GHz for approximately 25% of the sky. This allows us to constrain the CIB power spectrum down to ℓ ≳ 70. We present these CIB maps and the various processing and validation steps that we have performed to ensure their quality, as well as a comparison with previous studies. All our data products are made publicly available,4 thereby enabling the community to investigate a wide range of questions related to the universe’s large-scale structure.

Using the Planck 2015 data release (PR2) temperature maps, we separate Galactic thermal dust emission from cosmic infrared background (CIB) anisotropies. For this purpose, we implement a specifically tailored component-separation method, the so-called generalized needlet internal linear combination (GNILC) method, which uses spatial information (the angular power spectra) to disentangle the Galactic dust emission and CIB anisotropies. We produce significantly improved all-sky maps of Planck thermal dust emission, with reduced CIB contamination, at 353, 545, and 857 GHz. By reducing the CIB contamination of the thermal dust maps, we provide more accurate estimates of the local dust temperature and dust spectral index over the sky with reduced dispersion, especially at high Galactic latitudes above $b = \pm 20{\deg}$. We find that the dust temperature is $T = (19.4 \pm 1.3)$ K and the dust spectral index is $\beta = 1.6 \pm 0.1$ averaged over the whole sky, while $T = (19.4 \pm 1.5)$ K and $\beta = 1.6 \pm 0.2$ on 21 % of the sky at high latitudes. Moreover, subtracting the new CIB-removed thermal dust maps from the CMB-removed Planck maps gives access to the CIB anisotropies over 60 % of the sky at Galactic latitudes $|b| > 20{\deg}$. Because they are a significant improvement over previous Planck products, the GNILC maps are recommended for thermal dust science. The new CIB maps can be regarded as indirect tracers of the dark matter and they are recommended for exploring cross-correlations with lensing and large-scale structure optical surveys. The reconstructed GNILC thermal dust and CIB maps are delivered as Planck products.

The widely used Milky Way dust-reddening map, the Schlegel–Finkbeiner–Davis (SFD) map, was found to contain extragalactic large-scale structure (LSS) imprints. Such contamination is inherent in maps based on infrared emission, which pick up not only Galactic dust but also the cosmic infrared background (CIB). When SFD is used for extinction correction, overcorrection occurs in a spatially correlated and redshift-dependent manner, which could impact precision cosmology using galaxy clustering, lensing, and Type Ia supernova distances. Similarly, LSS imprints in other Galactic templates can affect intensity mapping and cosmic microwave background experiments. This paper presents a generic way to remove LSS traces in Galactic maps and applies it to SFD. First, we measure descriptive summary statistics of the CIB in SFD by cross-correlating the map with spectroscopic galaxies and quasars in the Sloan Digital Sky Survey tomographically as functions of redshift and angular scale. To reconstruct the LSS on the map level, however, additional information on the phases is needed. We build a large set of 180 overcomplete, full-sky basis template maps from the density fields of over 600 million galaxies in the Wide-field Infrared Survey Explorer and find a linear combination that reproduces all of the high-dimensional tomographic two-point statistics of the CIB in SFD. After subtracting this reconstructed LSS/CIB field, the end product is a full-sky Galactic dust-reddening map that supersedes SFD, carrying all Galactic features therein, with maximally suppressed CIB. We release this new dust map dubbed CSFD—the corrected SFD—at https://idv.sinica.edu.tw/ykchiang/CSFD.html and NASA’s LAMBDA archive.

We explore the use of the cosmic infrared background (CIB) as a tracer of the large scale structures for cross-correlating with the cosmic microwave background (CMB) and exploit the integrated Sachs–Wolfe (ISW) effect. We used an improved linear CIB model from our previous work and derived the theoretical CIB×ISW cross-correlation for different Planck HFI frequencies (217, 353, 545 and 857 GHz) and IRAS (3000 GHz). As expected, we predict a positive cross-correlation between the CIB and the CMB whose amplitude decreases rapidly at small scales. We perform a signal-to-noise ratio (S/N) analysis of the predicted cross-correlation. In the ideal case when the cross-correlation is obtained over 70% (40%) of the sky without residual contaminants (e.g. galactic dust) in maps, the S/N ranges from 4.2 to 5.6 (3.2 to 4.3); the highest S/N comes from 857 GHz. A Fisher matrix analysis shows that an ISW signal detected with a S/N this high on the 40% sky can considerably improve the constraints on the cosmological parameters; constraints on the equation of state of the dark energy especially are improved by 80%. We then performed a more realistic analysis considering the effect of residual galactic dust contamination in CIB maps. We calculated the dust power spectra for different frequencies and sky fractions that dominate the CIB power spectra at the lower multipoles we are interested in. Considering a conservative 10% residual level of galactic dust in the CIB power spectra, we observe that the S/N drops drastically, which makes it very challenging to detect the ISW. To determine the capability of current maps to detect the ISW effect through this method, we measured the cross-correlation of the CIB and the CMB Planck maps on the so-called GASS field, which covers an area of ∼11% in the southern hemisphere. We find that with such a small sky fraction and the dust residuals in the CIB maps, we do not detect any ISW signal, and the measured cross-correlation is consistent with zero. To avoid degrading the S/N for the ISW measurement by more than 10% on the 40% sky, we find that the dust needs to be cleaned up to the 0.01% level on the power spectrum.

The cosmic infrared background (CIB) sourced by infrared emission from dusty star-forming galaxies is a valuable source of information on the star formation history of the Universe. In measurements of the millimeter sky at frequencies higher than $\ensuremath{\sim}300\text{ }\text{ }\mathrm{GHz}$, the CIB and thermal emission from Galactic dust dominate. Insufficient understanding of the CIB contribution at lower frequencies can hinder efforts to measure the kinetic Sunyaev-Zeldovich spectrum on small scales as well as new physics that affects the damping tail of the cosmic microwave background (CMB). The Planck satellite has measured with high fidelity the CIB at 217, 353, 545 and 857 GHz. On very large scales, this measurement is limited by our ability to separate the CIB from Galactic dust, but on intermediate scales, the measurements are limited by sample variance in the underlying matter field traced by the CIB. We show how significant improvements (20--100%) can be obtained on parameters of star formation models by cross-correlating the CIB (as measured from existing Planck maps or upcoming CCAT-prime maps) with upcoming mass maps inferred from gravitational lensing of the CMB. This improvement comes from improved knowledge of the redshift distribution of star-forming galaxies as well as through the use of the unbiased matter density inferred from CMB lensing mass maps to cancel the sample variance in the CIB field. We also find that further improvements can be obtained on CIB model parameters if the cross-correlation of the CIB with CMB lensing is measured over a wider area while restricting the more challenging CIB autospectrum measurement to the cleanest 5% of the sky.

The measurement of the anisotropies in the cosmic infrared background (CIB) is a powerful mean of studying the evolution of galaxies and large-scale structures. These anisotropies have been measured by a number of experiments, from the far-infrared to the millimeter. One of the main impediments to an accurate measurement on large scales is the contamination of the foreground signal by Galactic dust emission. Our goal is to show that we can remove the Galactic cirrus contamination using HI data, and thus accurately measure the clustering of starburst galaxies in the CIB. We use observations of the ELAIS N1 field at far-infrared (100 and 160{\mu}m) and radio (21 cm) wavelengths. We compute the correlation between dust emission, traced by far-infrared observations, and HI gas traced by 21cm observations, and derive dust emissivities that enable us to subtract the cirrus emission from the far-infrared maps. We then derive the power spectrum of the CIB anisotropies, as well as its mean level at 100{\mu}m and 160{\mu}m. We also combine the HI data and Spitzer total power mode absolute measurements to determine the CIB mean level at 160{\mu}m. We find B160=0.77\pm0.04\pm0.12MJy/sr,where the first error is statistical and the second systematic. Combining this measurement with the B100/B160 color of the correlated anisotropies, we also derive the CIB mean at 100 {\mu}m, B100=0.24\pm0.08\pm0.04 MJy/sr. This measurement is in line with values obtained with recent models of infrared galaxy evolution and Herschel/PACS data, but is much smaller than the previous DIRBE measurements. The use of high-angular resolution Hi data is mandatory to accurately differentiate the cirrus from the CIB emission. The 100 {\mu}m IRAS map (and thus the map developed by Schlegel and collaborators) in such extragalactic fields is highly contaminated by the CIB anisotropies and hence cannot be used as a Galactic cirrus tracer.

We search for potential galactic and extragalactic dust contamination in thermal Sunyaev–Zeldovich maps derived from the Planck data. To test for contamination, we apply a variety of galactic dust and cosmic infrared background (CIB) models to the data as part of the y map reconstruction process. We evaluate the level of contamination by cross-correlating these y maps with mass tracers based on weak lensing data. The lensing data we use are the convergence map, κ, from the Red Sequence Cluster Lensing survey, and the cosmic microwave background (CMB) lensing potential map, ϕ, from the Planck Collaboration. We make a CIB-subtracted y map and measure the cross-correlation between it and the lensing data. By comparing it with CIB-contaminated cross-correlation, we find that the cross-correlation between κ and y is only slightly contaminated by CIB signal, at the level of 6.8 ± 3.5%, which implies that previous detections of κ × y are robust to CIB contamination. However, we find that ϕ × y is more significantly contaminated, by 16.7 ± 3.5%, because the CMB lensing potential probes higher redshift sources that overlap more with the CIB sources. We find that Galactic dust does not significantly contaminate either cross-correlation signal.

Foreground emission and scattered light from interplanetary dust (IPD) particles and emission from Galactic stellar sources are the greatest obstacles to determining the cosmic infrared background (CIB) from diffuse sky measurements in the ~1-5 μm range. We use ground-based observational limits on the K-band intensity of the CIB in conjunction with sky maps obtained by the Diffuse Infrared Background Experiment (DIRBE) on the Cosmic Background Explorer satellite to reexamine the limits on the CIB at 1.25, 3.5, and 4.9 μm. Adopting a CIB intensity of 7.4 nW m−2 sr−1 at 2.2 μm, and using the 2.2 μm DIRBE sky map from which the emission from the IPD cloud has been subtracted, we create a spatial template of the Galactic stellar contribution to the diffuse infrared sky. This template is then used to subtract the contribution of the diffuse Galactic stellar emission from the IPD emission-subtracted DIRBE sky maps at 1.25, 3.5, and 4.9 μm. The DIRBE 100 μm data are used to estimate the small contribution of emission from interstellar dust at 3.5 and 4.9 μm. Our method significantly reduces the errors associated with the subtraction of Galactic starlight, leaving only the IPD emission component as the primary obstacle to the detection of the CIB at these wavelengths. The analysis leads to a tentative detection of the CIB at 3.5 μm with an intensity of νIν={9.9+0.312[ν0ICIB(ν0)-7.4]} ± 2.9 nW m−2 sr−1, where ν0ICIB(ν0) is the CIB intensity at 2.2 μm in units of nW m−2 sr−1. The analysis also yields new upper limits (95% confidence limit) on the CIB at 1.25 and 4.9 μm of 68 and 36 nW m−2 sr−1, respectively. The cosmological implications of these results are discussed in this Letter.

We present a full-sky 100 μm map that is a reprocessed composite of the COBE/DIRBE and IRAS/ISSA maps, with the zodiacal foreground and confirmed point sources removed. Before using the ISSA maps, we remove the remaining artifacts from the IRAS scan pattern. Using the DIRBE 100 and 240 μm data, we have constructed a map of the dust temperature so that the 100 μm map may be converted to a map proportional to dust column density. The dust temperature varies from 17 to 21 K, which is modest but does modify the estimate of the dust column by a factor of 5. The result of these manipulations is a map with DIRBE quality calibration and IRAS resolution. A wealth of filamentary detail is apparent on many different scales at all Galactic latitudes. In high-latitude regions, the dust map correlates well with maps of H I emission, but deviations are coherent in the sky and are especially conspicuous in regions of saturation of H I emission toward denser clouds and of formation of H2 in molecular clouds. In contrast, high-velocity H I clouds are deficient in dust emission, as expected. To generate the full-sky dust maps, we must first remove zodiacal light contamination, as well as a possible cosmic infrared background (CIB). This is done via a regression analysis of the 100 μm DIRBE map against the Leiden-Dwingeloo map of H I emission, with corrections for the zodiacal light via a suitable expansion of the DIRBE 25 μm flux. This procedure removes virtually all traces of the zodiacal foreground. For the 100 μm map no significant CIB is detected. At longer wavelengths, where the zodiacal contamination is weaker, we detect the CIB at surprisingly high flux levels of 32 ± 13 nW m-2 sr-1 at 140 μm and of 17 ± 4 nW m-2 sr-1 at 240 μm (95% confidence). This integrated flux ~2 times that extrapolated from optical galaxies in the Hubble Deep Field. The primary use of these maps is likely to be as a new estimator of Galactic extinction. To calibrate our maps, we assume a standard reddening law and use the colors of elliptical galaxies to measure the reddening per unit flux density of 100 μm emission. We find consistent calibration using the B-R color distribution of a sample of the 106 brightest cluster ellipticals, as well as a sample of 384 ellipticals with B-V and Mg line strength measurements. For the latter sample, we use the correlation of intrinsic B-V versus Mg2 index to tighten the power of the test greatly. We demonstrate that the new maps are twice as accurate as the older Burstein-Heiles reddening estimates in regions of low and moderate reddening. The maps are expected to be significantly more accurate in regions of high reddening. These dust maps will also be useful for estimating millimeter emission that contaminates cosmic microwave background radiation experiments and for estimating soft X-ray absorption. We describe how to access our maps readily for general use.

Context. Gaia data and stellar surveys open the way to the construction of detailed 3D maps of the Galactic interstellar (IS) dust based on the synthesis of star distances and extinctions. Dust maps are tools of broad use, also for Gaia-related Milky Way studies. Aims. Reliable extinction measurements require very accurate photometric calibrations. We show the first step of an iterative process linking 3D dust maps and photometric calibrations, and improving them simultaneously. Methods. Our previous 3D map of nearby IS dust was used to select low-reddening SDSS/APOGEE-DR14 red giants, and this database served for an empirical effective temperature- and metallicity-dependent photometric calibration in the Gaia G and 2MASS Ks bands. This calibration has been combined with Gaia G-band empirical extinction coefficients recently published, G, J, and Ks photometry and APOGEE atmospheric parameters to derive the extinction of a large fraction of the survey targets. Distances were estimated independently using isochrones and the magnitude-independent extinction KJ−Ks. This new dataset has been merged with the one used for the earlier version of dust map. A new Bayesian inversion of distance-extinction pairs has been performed to produce an updated 3D map. Results. We present several properties of the new map. A comparison with 2D dust emission reveals that all large dust shells seen in emission at middle and high latitudes are closer than 300 pc. The updated distribution constrains the well-debated, X-ray bright North Polar Spur to originate beyond 800 pc. We use the Orion region to illustrate additional details and distant clouds. On the large scale the map reveals a complex structure of the Local Arm. Chains of clouds of 2–3 kpc in length appear in planes tilted by ≃15° with respect to the Galactic plane. A series of cavities oriented along a l ≃ 60–240° axis crosses the Arm. Conclusions. The results illustrate the ongoing synergy between 3D mapping of IS dust and stellar calibrations in the context of Gaia. Dust maps provide prior foregrounds for future calibrations appropriate to different target characteristics or ranges of extinction, allowing us in turn to increase extinction data and produce more detailed and extended maps.

Aims. Three-dimensional (3D) maps of Galactic interstellar dust are a tool for a wide range of uses. We aim to construct 3D maps of dust extinction in the Local Arm and surrounding regions. Methods. To do this, Gaia EDR3 photometric data were combined with 2MASS measurements to derive extinction towards stars with accurate photometry and relative uncertainties on EDR3 parallaxes of less than 20%. We applied our hierarchical inversion algorithm adapted to inhomogeneous spatial distributions of target stars to this catalogue of individual extinctions. Results. We present the updated 3D dust extinction distribution and provide an estimate of the error on integrated extinctions from the Sun to each area in the 3D map. The full computational area is similar to the one of the previous DR2 map, that is to say with a 6 × 6 × 0.8 kpc3 volume around the Sun. Due to the addition of fainter target stars, the volume in which the clouds can be reconstructed has increased. Due to the improved accuracy of the parallaxes and photometric data in EDR3, extinctions among neighbouring targets are more consistent, allowing one to reach an increased contrast in the dense areas, while cavity contours are more regular. We show several comparisons with recent results on dust and star distributions. The wavy pattern around the Plane of the dust concentrations is better seen and exists over large regions. Its mean vertical peak-to-peak amplitude is of the order of 300 pc; interestingly, it is similar to the vertical period of the spectacular snail-shaped stellar kinematical pattern discovered in Gaia data. Conclusions. The Gaia EDR3 catalogue allows for a significant improvement of the extinction maps to be made, both in extent and quality. The hierarchical technique confirms its efficiency in the inversion of massive datasets. Future comparisons between 3D maps of interstellar matter and stellar distributions may help to understand which mergers or internal perturbations have shaped the Galaxy within the first 3 kpc.

We present the sensitivity calculation results for observing the Cosmic Infrared Background (CIRB) by the Multi-purpose IR Imaging System (MIRIS), which will be launched in 2010 as a main payload of the Science and Technology Satellite 3 (STSAT-3). MIRIS will observe in I () and H () band with a degree field of view to obtain the large scale structure ( degree) of the CIRB. With the given specifications of the MIRIS, our sensitivity calculation results show that the MIRIS has a detection limit of (I band) and (H band), which is appropriate to observe the large scale structure of CIRB.

Context. The Planck data releases have provided the community with submillimetre and full-sky radio observations at unprecedented resolutions. We make use of the Planck 353, 545, and 857 GHz maps alongside the IRAS 3000 GHz map. These maps contain information on the cosmic microwave background (CMB), cosmic infrared background (CIB), extragalactic point sources, and diffuse thermal dust emission. Aims. We aim to determine the modified black-body (MBB) model parameters of thermal dust emission in total intensity and produce all-sky maps of pure thermal dust, having separated this Galactic component from the CMB and CIB. Methods. This separation is completed using a new, sparsity-based, parametric method, Parameter Recovery Exploiting Model Informed Sparse Estimates (premise). The method is comprised of three main stages: 1) filtering the raw data to reduce the effect of the CIB on the MBB fit; 2) fitting an MBB model to the filtered data across super-pixels of various sizes determined by the algorithm itself; and 3) refining these super-pixel estimates into full-resolution maps of the MBB parameters. Results. We present our maps of MBB temperature, spectral index, and optical depth at 5 arcmin resolution and compare our estimates to those of GNILC and to the two-step MBB fit presented by the Planck Collaboration in 2013. Conclusions. By exploiting sparsity we avoid the need for smoothing, enabling us to produce the first full-resolution MBB parameter maps from intensity measurements of thermal dust emission. We consider the premise parameter estimates to be competitive with the existing state-of-the-art solutions, outperforming these methods within low signal-to-noise regions as we account for the CIB without removing thermal dust emission through oversmoothing.

We use analytic computations to predict the power spectrum as well as the bispectrum of Cosmic Infrared Background (CIB) anisotropies. Our approach is based on the halo model and takes into account the mean luminosity-mass relation. The model is used to forecast the possibility to simultaneously constrain cosmological, CIB and halo occupation distribution (HOD) parameters in the presence of foregrounds. For the analysis we use wavelengths in eight frequency channels between 200 and 900 GHz with survey specifications given by Planck and LiteBird. We explore the sensitivity to the model parameters up to multipoles of ℓ = 1000 using auto- and cross-correlations between the different frequency bands. With this setting, cosmological, HOD and CIB parameters can be constrained to a few percent. Galactic dust is modeled by a power law and the shot noise contribution as a frequency dependent amplitude which are marginalized over. We find that dust residuals in the CIB maps only marginally influence constraints on standard cosmological parameters. Furthermore, the bispectrum yields tighter constraints (by a factor four in 1σ errors) on almost all model parameters while the degeneracy directions are very similar to the ones of the power spectrum. The increase in sensitivity is most pronounced for the sum of the neutrino masses. Due to the similarity of degeneracies a combination of both analysis is not needed for most parameters. This, however, might be due to the simplified bias description generally adopted in such halo model approaches.

The most promising avenue for detecting primordial gravitational waves from cosmic inflation is through measurements of degree-scale cosmic microwave background (CMB) B-mode polarization. This approach must face the challenge posed by gravitational lensing of the CMB, which obscures the signal of interest. Fortunately, the lensing effects can be partially removed by combining high-resolution E-mode measurements with an estimate of the projected matter distribution. For near-future experiments, the best estimate of the latter will arise from co-adding internal reconstructions (derived from the CMB itself) with external tracers such as the cosmic infrared background (CIB). In this work, we characterize how foregrounds impact the delensing procedure when CIB intensity, I, is used as the matter tracer. We find that higher point functions of the CIB and Galactic dust such as 〈BEI〉c and 〈EIEI〉c can, in principle, bias the power spectrum of delensed B-modes. To quantify these, we first estimate the dust residuals in currently available CIB maps and upcoming, foreground-cleaned Simons Observatory CMB data. Then, using non-Gaussian simulations of Galactic dust – extrapolated to the relevant frequencies, assuming the spectral index of polarized dust emission to be fixed at the value determined by Planck – we show that the bias to any primordial signal is small compared to statistical errors for ground-based experiments, but might be significant for space-based experiments probing very large angular scales. However, mitigation techniques based on multifrequency cleaning appear to be very effective. We also show, by means of an analytical model, that the bias arising from the higher point functions of the CIB itself ought to be negligible.