Cosmic Infrared Background Research Articles

Constraining the physics of star formation from CIB-cosmic shear cross-correlations

ABSTRACT Understanding the physics of star formation is one of the key problems facing modern astrophysics. The cosmic infrared background (CIB), sourced by the emission from all dusty star-forming galaxies since the epoch of reionization, is a complementary probe to study the star formation history, as well as an important extragalactic foreground for studies of the cosmic microwave background. In this paper, we make high signal-to-noise measurements of the cross-correlation between maps of the CIB from the Planck experiment, and cosmic shear measurements from the Dark Energy Survey and Kilo-Degree Survey. Cosmic shear is a direct tracer of the matter distribution and thus we can use its cross-correlation with the CIB to directly test our understanding of the link between the star formation rate (SFR) density and the matter density. We use our measurements to place constraints on a halo-based model of the SFR that parametrizes the efficiency with which gas is transformed into stars as a function of halo mass and redshift. These constraints are enhanced by using model-independent measurements of the bias-weighted SFR density extracted from the tomographic cross-correlation of galaxies and the CIB. We are able to place constraints on the peak efficiency at low redshifts, $\eta =0.445^{+0.055}_{-0.11}$, and on the halo mass at which this peak efficiency is achieved today log10(M1/M⊙) = 12.17 ± 0.25. Our constraints are in excellent agreement with direct measurements of the SFR density, as well as other CIB-based studies.

Importance of intracluster scattering and relativistic corrections from tSZ effect with cosmic infrared background

ABSTRACT The Sunyaev–Zeldovich effect towards clusters of galaxies has become a standard probe of cosmology. It is caused by the scattering of photons from the cosmic microwave background (CMB) by the hot cluster electron gas. In a similar manner, other photon backgrounds can be scattered when passing through the cluster medium. This problem has been recently considered for the radio and the cosmic infrared background. Here, we revisit the discussion of the cosmic infrared background (CIB) including several additional effects that were omitted before. We discuss the intracluster scattering of the CIB and the role of relativistic temperature corrections to the individual cluster and all-sky averaged signals. We show that the all-sky CIB distortion introduced by the scattering of the photon field was underestimated by a factor of ≃1.5 due to neglecting the intracluster scattering contribution. The CIB photons can scatter with the thermal electrons of both the parent halo or another halo, meaning that there are two ways to gain energy. Therefore, energy is essentially transferred twice from the thermal electrons to the CIB. We carefully clarify the origin of various effects in the calculation of the average CIB and also scattered signals. The single-cluster CIB scattering signal also exhibits a clear redshift dependence, which can be used in cosmological analyses, as we describe both analytically and numerically. This may open a new way for cosmological studies with future CMB experiments.

Herschel Optimized Tau and Temperature (HOTT) Maps: Uncertainty Analysis and Robust Parameter Extraction

We introduce the HOTT dust optical depth and temperature maps parameterizing thermal dust emission. Such maps have revolutionized studies of the distribution of matter in molecular clouds and processes relevant to star formation, including virial stability. HOTT maps for a suite of fields, including the Herschel Gould Belt Survey, are available online. The standardization of our robust pipeline for modified blackbody fitting of the spectral energy distribution (SED) of high-quality archival submillimeter data from the Herschel Space Observatory is based on a thorough analysis and quantification of the uncertainties of the data. This enables proper weighting in the SED fits. The uncertainties assessed fall into four main categories: instrument noise; the cosmic infrared background anisotropy, a contaminating sky signal; gradient-related noise arising because of dust signal morphology; and calibration uncertainty, scaling with the signal strength. Zero-level adjustments are important too. An analysis of residuals from the SED fits across many fields supports the overall appropriateness of the assumed modified blackbody model and points to where it breaks down. Finding χ 2 distributions close to the theoretical expectation boosts confidence in the pipeline and the optimized quality of the parameter maps and their estimated uncertainties. We compared our HOTT parameter maps to those from earlier studies to understand and quantify the potential for systematic differences.

ALMACAL IX: Multiband ALMA survey for dusty star-forming galaxies and the resolved fractions of the cosmic infrared background

ABSTRACT Wide, deep, blind continuum surveys at submillimetre/millimetre (submm/mm) wavelengths are required to provide a full inventory of the dusty, distant Universe. However, conducting such surveys to the necessary depth, with sub-arcsec angular resolution, is prohibitively time-consuming, even for the most advanced submm/mm telescopes. Here, we report the most recent results from the ALMACAL project, which exploits the ‘free’ calibration data from the Atacama Large Millimetre/submillimetre Array (ALMA) to map the lines of sight towards and beyond the ALMA calibrators. ALMACAL has now covered 1001 calibrators, with a total sky coverage around 0.3 deg2, distributed across the sky accessible from the Atacama desert, and has accumulated more than 1000 h of integration. The depth reached by combining multiple visits to each field makes ALMACAL capable of searching for faint, dusty, star-forming galaxies (DSFGs), with detections at multiple frequencies to constrain the emission mechanism. Based on the most up-to-date ALMACAL data base, we report the detection of 186 DSFGs with flux densities down to S870 µm ∼ 0.2 mJy, comparable with existing ALMA large surveys but less susceptible to cosmic variance. We report the number counts at five wavelengths between 870 μm and 3 mm, in ALMA bands 3, 4, 5, 6, and 7, providing a benchmark for models of galaxy formation and evolution. By integrating the observed number counts and the best-fitting functions, we also present the resolved fraction of the cosmic infrared background (CIB) and the CIB spectral shape. Combining existing surveys, ALMA has currently resolved about half of the CIB in the submm/mm regime.

CENN: A fully convolutional neural network for CMB recovery in realistic microwave sky simulations

Context. Component separation is the process with which emission sources in astrophysical maps are generally extracted by taking multi-frequency information into account. It is crucial to develop more reliable methods for component separation for future cosmic microwave background (CMB) experiments such as the Simons Observatory, the CMB-S4, or the LiteBIRD satellite. Aims. We aim to develop a machine learning method based on fully convolutional neural networks called the CMB extraction neural network (CENN) in order to extract the CMB signal in total intensity by training the network with realistic simulations. The frequencies we used are the Planck channels 143, 217, and 353 GHz, and we validated the neural network throughout the sky and at three latitude intervals: 0° < |b| < 5°, 5° < |b| < 30°, and 30° < |b| < 90°, Moreover, we used neither Galactic nor point-source (PS) masks. Methods. To train the neural network, we produced multi-frequency realistic simulations in the form of patches of 256 × 256 pixels that contained the CMB signal, the Galactic thermal dust, cosmic infrared background, and PS emissions, the thermal Sunyaev–Zel’dovich effect from galaxy clusters, and instrumental noise. After validating the network, we compared the power spectra from input and output maps. We analysed the power spectrum from the residuals at each latitude interval and throughout the sky, and we studied how our model handled high contamination at small scales. Results. We obtained a CMB power spectrum with a mean difference between input and output of 13 ± 113 µK2 for multipoles up to above 4000. We computed the residuals, obtaining 700 ± 60 µK2 for 0° < |b| < 5°, 80 ± 30 µK2 for 5° < |b| < 30°, and 30 ± 20 µK2 for 30° < |b| < 90° for multipoles up to above 4000. For the entire sky, we obtained 30 ± 10 µK2 for l ≤ 1000 and 20 ± 10 µK2 for l ≤ 4000. We validated the neural network in a single patch with strong contamination at small scales, obtaining a difference between input and output of 50 ± 120 µK2 and residuals of 40 ± 10 µK2 up to l ~ 2500. In all cases, the uncertainty of each measure was taken as the standard deviation. Conclusions. The results show that fully convolutional neural networks are promising methods for performing component separation in future CMB experiments. Moreover, we show that CENN is reliable against different levels of contamination from Galactic and PS foregrounds at both large and small scales.

Impact of extragalactic foregrounds on internal delensing of the CMB B -mode polarization

The search for primordial $B$-mode polarization of the CMB is limited by the sample variance of $B$-modes produced at later times by gravitational lensing. Constraints can be improved by `delensing': using some proxy of the matter distribution to partially remove the lensing-induced $B$-modes. Current and soon-upcoming experiments will infer a matter map -- at least in part -- from the temperature anisotropies of the CMB. These reconstructions are contaminated by extragalactic foregrounds: radio-emitting galaxies, the cosmic infrared background, or the Sunyaev--Zel'dovich effects. Using the Websky simulations, we show that the foregrounds add spurious power to the angular auto-spectrum of delensed $B$-modes via non-Gaussian higher-point functions, biasing constraints on the tensor-to-scalar ratio, $r$. We consider an idealized experiment similar to the Simons Observatory, with no Galactic or atmospheric foregrounds. After removing point sources detectable at 143 GHz and reconstructing lensing from CMB temperature modes $l<3500$ using a Hu-Okamoto quadratic estimator (QE), we infer a value of $r$ that is $1.5\,\sigma$ higher than the true $r=0$. Reconstructing instead from a minimum-variance ILC map only exacerbates the problem, bringing the bias above $3\,\sigma$. When the $TT$ estimator is co-added with other QEs or with external matter tracers, new couplings ensue which partially cancel the diluted bias from $TT$. We provide a simple and effective prescription to model these effects. In addition, we demonstrate that the point-source-hardened or shear-only QEs can not only mitigate the biases to acceptable levels, but also lead to lower power than the Hu-Okamoto QE after delensing. Thus, temperature-based reconstructions remain powerful tools in the quest to measure $r$.

Forecasting ground-based sensitivity to the Rayleigh scattering of the CMB in the presence of astrophysical foregrounds

The Rayleigh scattering of cosmic microwave background (CMB) photons off the neutral hydrogen produced during recombination effectively creates an additional scattering surface after recombination that encodes new cosmological information, including the expansion and ionization history of the universe. A first detection of Rayleigh scattering is a tantalizing target for next-generation CMB experiments. We have developed a Rayleigh scattering forecasting pipeline that includes instrumental effects, atmospheric noise, and astrophysical foregrounds (e.g., Galactic dust, cosmic infrared background, or CIB, and the thermal Sunyaev-Zel'dovich effect). We forecast the Rayleigh scattering detection significance for several upcoming ground-based experiments, including SPT-3G+, Simons Observatory, CCAT-prime, and CMB-S4, and examine the limitations from atmospheric and astrophysical foregrounds as well as potential mitigation strategies. When combined with Planck data, we estimate that the ground-based experiments will detect Rayleigh scattering with a significance between 1.6 and 3.7, primarily limited by atmospheric noise and the CIB.

Cosmic star formation history with tomographic cosmic infrared background-galaxy cross-correlation

In this work we present a new method for probing the star formation history of the Universe, namely tomographic cross-correlation between the cosmic infrared background (CIB) and galaxy samples. The galaxy samples are from the Kilo-Degree Survey (KiDS), while the CIB maps are made from Planck sky maps at 353, 545, and 857 GHz. We measure the cross-correlation in harmonic space within 100 &lt; ℓ &lt; 2000 with a significance of 43σ. We model the cross-correlation with a halo model, which links CIB anisotropies to star formation rates (SFRs) and galaxy abundance. We assume that the SFR has a lognormal dependence on halo mass and that the galaxy abundance follows the halo occupation distribution (HOD) model. The cross-correlations give a best-fit maximum star formation efficiency of ηmax = 0.41−0.14+0.09 at a halo mass log10(Mpeak/M⊙) = 12.14 ± 0.36. The derived star formation rate density (SFRD) is well constrained up to z ∼ 1.5. The constraining power at high redshift is mainly limited by the KiDS survey depth. We also show that the constraint is robust to uncertainties in the estimated redshift distributions of the galaxy sample. A combination with external SFRD measurements from previous studies gives log10(Mpeak/M⊙) = 12.42−0.19+0.35. This tightens the SFRD constraint up to z = 4, yielding a peak SFRD of 0.09−0.004+0.003 M⊙ yr−1 Mpc−3 at z = 1.74−0.02+0.06, corresponding to a lookback time of 10.05−0.03+0.12 Gyr. Both constraints are consistent, and the derived SFRD agrees with previous studies and simulations. This validates the use of CIB tomography as an independent probe of the star formation history of the Universe. Additionally, we estimate the galaxy bias, b, of KiDS galaxies from the constrained HOD parameters and obtain an increasing bias from b = 1.1−0.31+0.17 at z = 0 to b = 1.96−0.64+0.18 at z = 1.5, which highlights the potential of this method as a probe of galaxy abundance. Finally, we provide a forecast for future galaxy surveys and conclude that, due to their considerable depth, future surveys will yield a much tighter constraint on the evolution of the SFRD.

First Identification of a CMB Lensing Signal Produced by 1.5Million Galaxies at z∼4: Constraints on Matter Density Fluctuations at High Redshift.

We report the first detection of the dark matter distribution around Lyman break galaxies (LBGs) at high redshift through the cosmic microwave background (CMB) lensing measurements with the public Planck PR3 κ map. The LBG sample consists of 1 473 106 objects with the median redshift of z∼4 that are identified in a total area of 305 deg^{2} observed by the Hyper Suprime-Cam Strategic Survey Program survey. After careful investigations of systematic uncertainties, such as contamination from foreground galaxies and cosmic infrared background, we obtain the significant detection of the CMB lensing signal at 5.1σ that is dominated by 2-halo term signals of the LBGs. Fitting a simple model consisting of the Navarro-Frenk-White profile and the linear-bias model, we obtain the typical halo mass of M_{h}=2.9_{-2.5}^{+9.5}×10^{11} h^{-1} M_{⊙}. Combining the CMB lensing and galaxy-galaxy clustering signals on the large scales, we demonstrate the first cosmological analysis at z∼4 that constrains (Ω_{m0},σ_{8}). We find that our constraint on σ_{8} is roughly consistent with the Planck cosmology, while this σ_{8} constraint is lower than the Planck cosmology over the 1σ level. This study opens up a new window for constraining cosmological parameters at high redshift by the combination of CMB and high-z galaxies, as well as studying the interplay between galaxy evolution and large-scale structure at such high redshift, by upcoming CMB and optical and near-infrared imaging surveys.

Simulated catalogs and maps of radio galaxies at millimeter wavelengths in Websky

We present simulated millimeter-wavelength maps and catalogs of radio galaxies across the full sky that trace the nonlinear clustering and evolution of dark matter halos from the Websky simulation at z < 4.6 and M halo > 1012 m ⊙/h, and the accompanying framework for generating a new sample of radio galaxies from any halo catalog of positions, redshifts, and masses. Object fluxes are generated using a hybrid approach that combines (1) existing astrophysical halo models of radio galaxies from the literature to determine the positions and rank-ordering of the observed fluxes with (2) empirical models from the literature based on fits to the observed distribution of flux densities and (3) spectral indices drawn from an empirically-calibrated frequency-dependent distribution. The resulting population of radio galaxies is in excellent agreement with the number counts, polarization fractions, and distribution of spectral slopes from the data from observations at millimeter wavelengths from 20-200 GHz, including Planck, ALMA, SPT, and ACT. Since the radio galaxies are correlated with the existing cosmic infrared background (CIB), Compton-y (tSZ), and CMB lensing maps from Websky, our model makes new predictions for the cross-correlation power spectra and stacked profiles of radio galaxies and these other components. These simulations will be important for unbiased analysis of a wide variety of observables that are correlated with large-scale structure, such as gravitational lensing and SZ clusters.

Inverse-Compton scattering of the cosmic infrared background

The thermal Sunyaev-Zel'dovich (tSZ) effect is the distortion generated in the cosmic microwave background (CMB) spectrum by the inverse-Compton scattering of CMB photons off free, energetic electrons, primarily located in the intracluster medium. Cosmic infrared background (CIB) photons from thermal dust emission in star-forming galaxies are expected to undergo the same process. In this work, we perform the first calculation of the resulting tSZ-like distortion in the CIB. Focusing on the CIB monopole, we use a halo model approach to calculate both the CIB signal and the Compton-$y$ field that generates the distortion. We self-consistently account for the redshift coevolution of the CIB and Compton-$y$ fields: they are (partially) sourced by the same dark matter halos, which introduces new aspects to the calculation as compared to the CMB case. We find that the inverse-Compton distortion to the CIB monopole spectrum has a positive (negative) peak amplitude of $\ensuremath{\approx}4\text{ }\text{ }\mathrm{Jy}/\mathrm{sr}$ ($\ensuremath{\approx}\ensuremath{-}5\text{ }\text{ }\mathrm{Jy}/\mathrm{sr}$) at 2260 GHz (940 GHz). In contrast to the usual tSZ effect, the distortion to the CIB spectrum has two null frequencies, at approximately 196 and 1490 GHz. We perform a Fisher matrix calculation to forecast the detectability of this new distortion signal by future experiments. PIXIE would have sufficient instrumental sensitivity to detect the signal at $4\ensuremath{\sigma}$, but foreground contamination reduces the projected signal-to-noise by a factor of $\ensuremath{\approx}70$. A future ESA Voyage 2050 spectrometer could detect the CIB distortion at $\ensuremath{\approx}5\ensuremath{\sigma}$ significance, even after marginalizing over foregrounds. A measurement of this signal would provide new information on the star formation history of the Universe, and the distortion anisotropies may be accessible by near-future ground-based experiments.

CMB lensing reconstruction biases from masking extragalactic sources

Observed cosmic microwave background (CMB) temperature and polarization maps can be powerful cosmological probes and used for CMB lensing reconstruction. However, the CMB maps are inevitably contaminated by foregrounds, some of which are usually masked to perform the analysis. If this mask is correlated to the lensing signal, measurements over the unmasked sky may give biased estimates and hence biased cosmological inferences. For example, masking extragalactic astrophysical emissions associated with objects located in dark matter halos will systematically remove parts of the sky that have a mass density higher than average. This can lead to a modified lensed CMB power spectrum over the unmasked area and biased measurements of lensing reconstruction auto- and cross-correlation power spectra with external matter tracers (both from the direct impact on the lensing power and via modifications to the lensing-dependent reconstruction power spectrum corrections, ${N}_{L}^{(0)}$, ${N}_{L}^{(1)}$, and ${N}_{L}^{(3/2)}$). For direct masking of the CMB lensing field, we give an approximate halo-model prediction of the size of the effect and derive simple analytic models for point sources and threshold masks constructed on a correlated Gaussian foreground field. We show that biases are significantly reduced by optimal filtering of the CMB maps in the lensing reconstruction, which effectively fills back some of the information in small mask holes. We test the resulting lensing power spectrum biases on numerical simulations, masking radio sources, and peaks of thermal Sunyaev-Zeldovich (tSZ) and cosmic infrared background (CIB) emission. For radio point sources, the remaining bias is negligible, but a temperature lensing reconstruction power spectrum bias remains at the 0.5%--5% level if clusters with a mass larger than $\ensuremath{\sim}{10}^{14}\text{ }\text{ }{M}_{\ensuremath{\bigodot}}/h$ are masked or if peaks in the CIB and tSZ maps are removed with ${f}_{\mathrm{sky}}^{\mathrm{mask}}\ensuremath{\simeq}0.5--7%$. In any case, these biases can only be measured with a statistical significance $\ensuremath{\lesssim}2\ensuremath{\sigma}$ for future datasets. Moreover, we quantified the impact of the mask biases in the cross-correlation power spectrum between CMB lensing and tSZ and CIB and found them to be larger (up to $\ensuremath{\sim}30%$). We found that masking tSZ-selected galaxy clusters leads to the largest mask biases, potentially detectable with high significance, and should therefore be avoided as much as possible. For the most realistic masks we considered, masking biases can only be measured with marginal significance. We found that the calibration of cluster masses using CMB lensing, in particular for objects at $z\ensuremath{\lesssim}0.6$, might be significantly affected by mask biases for near-future observations if the lensing signal recovered inside the mask holes is used without further corrections. Conversely, mass calibration of high redshift objects will still deliver unbiased results.

Delensing the CMB with the cosmic infrared background: the impact of foregrounds

ABSTRACT 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.

Cross-correlation of Dark Energy Survey Year 3 lensing data with ACT and Planck thermal Sunyaev-Zel’dovich effect observations. I. Measurements, systematics tests, and feedback model constraints

We present a tomographic measurement of the cross-correlation between thermal Sunyaev-Zeldovich (tSZ) maps from ${\it Planck}$ and the Atacama Cosmology Telescope (ACT) and weak galaxy lensing shears measured during the first three years of observations of the Dark Energy Survey (DES Y3). This correlation is sensitive to the thermal energy in baryons over a wide redshift range, and is therefore a powerful probe of astrophysical feedback. We detect the correlation at a statistical significance of $21\sigma$, the highest significance to date. We examine the tSZ maps for potential contaminants, including cosmic infrared background (CIB) and radio sources, finding that CIB has a substantial impact on our measurements and must be taken into account in our analysis. We use the cross-correlation measurements to test different feedback models. In particular, we model the tSZ using several different pressure profile models calibrated against hydrodynamical simulations. Our analysis marginalises over redshift uncertainties, shear calibration biases, and intrinsic alignment effects. We also marginalise over $\Omega_{\rm m}$ and $\sigma_8$ using ${\it Planck}$ or DES priors. We find that the data prefers the model with a low amplitude of the pressure profile at small scales, compatible with a scenario with strong AGN feedback and ejection of gas from the inner part of the halos. When using a more flexible model for the shear profile, constraints are weaker, and the data cannot discriminate between different baryonic prescriptions.

Probing the rest-frame of the Universe with the near-IR cosmic infrared background

ABSTRACT While the cosmic microwave background (CMB) dipole is largely assumed to be entirely kinematic, there is evidence that part of it is primordial. Such a possibility arises in models implying a tilt, interpreted as a dark flow, across the observable Universe. The kinematic nature of the entire CMB dipole can be probed using the dipole of cosmic backgrounds from galaxies after the last scattering. The near-infrared (near-IR) cosmic infrared background (CIB) spectral energy distribution leads to an amplified dipole compared with the CMB. The CIB dipole is affected by galaxy clustering, decreasing with fainter, more distant galaxies, and by Solar System emissions and Galactic dust, which dominate the net CIB cosmological dipole in the optical/near-IR. We propose a technique that enables an accurate measurement of the kinematic near-IR CIB dipole. The CIB, effectively the integrated galaxy light (IGL), would be reconstructed from resolved galaxies in forthcoming space-borne wide surveys covering four bands, 0.9–2.5 μm. The galaxies will be subselected from the identified magnitude range where the dipole component from galaxy clustering is below the expected kinematic dipole. Using this technique, the dipole can be measured in each of the bands at the statistical signal-to-noise ratio S/N ≳50–100 with the forthcoming Euclid and Roman surveys, isolating the CMB dipole’s kinematic nature.

Candidate high-redshift protoclusters and lensed galaxies in the Planck list of high-z sources overlapping with Herschel-SPIRE imaging

ABSTRACT The Planck list of high-redshift source candidates (the PHz catalogue) contains 2151 peaks in the cosmic infrared background, unresolved by Planck’s 5 arcmin beam. Follow-up spectroscopic observations have revealed that some of these objects are $z\, {\approx }\, 2$ protoclusters and strong gravitational lenses but an unbiased survey has not yet been carried out. To this end, we have used archival Herschel-SPIRE observations to study a uniformly selected sample of 187 PHz sources. In contrast with follow-up studies that were biased towards bright, compact sources, we find that only one of our PHz sources is a bright gravitationally lensed galaxy (peak flux ${\gtrsim }\, 300$ mJy), indicating that such objects are rarer in the PHz catalogue than previously believed (&amp;lt;1 per cent). The majority of our PHz sources consist of many red, star-forming galaxies, demonstrating that typical PHz sources are candidate protoclusters. However, our new PHz sources are significantly less bright than found in previous studies and differ in colour, suggesting possible differences in redshift and star formation rate. None the less, 40 of our PHz sources contain ${\gt }\, 3\, \sigma$ galaxy overdensities, comparable to the fraction of ${\gt }\, 3\, \sigma$ overdensities found in earlier biased studies. We additionally use a machine-learning approach to identify less extreme (peak flux ${\sim }\, 100$ mJy) gravitationally lensed galaxies among Herschel-SPIRE observations of PHz sources, finding a total of seven candidates in our unbiased sample, and 13 amongst previous biased samples. Our new uniformly selected catalogues of ${\gt }\, 3\, \sigma$ candidate protoclusters and strong gravitational lenses provide interesting targets for follow up with higher resolution facilities, such as ALMA and JWST.

Measurement of the Relativistic Sunyaev–Zeldovich Correction in RX J1347.5-1145

We present a measurement of the relativistic corrections to the thermal Sunyaev–Zel’dovich (SZ) effect spectrum, the rSZ effect, toward the massive galaxy cluster RX J1347.5-1145 by combining submillimeter images from Herschel-SPIRE with millimeter wavelength Bolocam maps. Our analysis simultaneously models the SZ effect signal, the population of cosmic infrared background galaxies, and the galactic cirrus dust emission in a manner that fully accounts for their spatial and frequency-dependent correlations. Gravitational lensing of background galaxies by RX J1347.5-1145 is included in our methodology based on a mass model derived from the Hubble Space Telescope observations. Utilizing a set of realistic mock observations, we employ a forward modeling approach that accounts for the non-Gaussian covariances between the observed astrophysical components to determine the posterior distribution of SZ effect brightness values consistent with the observed data. We determine a maximum a posteriori (MAP) value of the average Comptonization parameter of the intracluster medium (ICM) within R 2500 to be 〈y〉2500 = 1.56 × 10−4, with corresponding 68% credible interval [1.42, 1.63] × 10−4, and a MAP ICM electron temperature of 〈T sz〉2500 = 22.4 keV with 68% credible interval spanning [10.4, 33.0] keV. This is in good agreement with the pressure-weighted temperature obtained from Chandra X-ray observations, 〈T x,pw〉2500 = 17.4 ± 2.3 keV. We aim to apply this methodology to comparable existing data for a sample of 39 galaxy clusters, with an estimated uncertainty on the ensemble mean 〈T sz〉2500 at the ≃ 1 keV level, sufficiently precise to probe ICM physics and to inform X-ray temperature calibration.

Joint constraints on cosmology and the impact of baryon feedback: Combining KiDS-1000 lensing with the thermal Sunyaev–Zeldovich effect from Planck and ACT

We conduct a pseudo-Cℓ analysis of the tomographic cross-correlation between 1000 deg2 of weak-lensing data from the Kilo-Degree Survey (KiDS-1000) and the thermal Sunyaev–Zeldovich (tSZ) effect measured by Planck and the Atacama Cosmology Telescope (ACT). Using HMX, a halo-model-based approach that consistently models the gas, star, and dark matter components, we are able to derive constraints on both cosmology and baryon feedback for the first time from these data, marginalising over redshift uncertainties, intrinsic alignment of galaxies, and contamination by the cosmic infrared background (CIB). We find our results to be insensitive to the CIB, while intrinsic alignment provides a small but significant contribution to the lensing–tSZ cross-correlation. The cosmological constraints are consistent with those of other low-redshift probes and prefer strong baryon feedback. The inferred amplitude of the lensing–tSZ cross-correlation signal, which scales as σ8(Ωm/0.3)0.2, is low by ∼2 σ compared to the primary cosmic microwave background constraints by Planck. The lensing–tSZ measurements are then combined with pseudo-Cℓ measurements of KiDS-1000 cosmic shear into a novel joint analysis, accounting for the full cross-covariance between the probes, providing tight cosmological constraints by breaking parameter degeneracies inherent to both probes. The joint analysis gives an improvement of 40% on the constraint of S8 = σ8Ωm/0.3 over cosmic shear alone, while providing constraints on baryon feedback consistent with hydrodynamical simulations, demonstrating the potential of such joint analyses with baryonic tracers such as the tSZ effect. We discuss remaining modelling challenges that need to be addressed if these baryonic probes are to be included in future precision-cosmology analyses.

Retrieving cosmological information from small-scale CMB foregrounds

We propose a new analysis of small-scale cosmic microwave background (CMB) data by introducing the cosmological dependency of the foreground signals, focussing first on the thermal Sunyaev-Zel’dovich (tSZ) power spectrum, derived from the halo model. We analyse the latest observations by the South Pole Telescope (SPT) of the high-ℓ power (cross) spectra at 95, 150, and 220 GHz, as the sum of CMB and tSZ signals, both depending on cosmological parameters and remaining contaminants. In order to perform faster analyses, we propose a new tSZ modelling based on machine learning algorithms (namely Random Forest). We show that the additional information contained in the tSZ power spectrum tightens constraints on cosmological and tSZ scaling relation parameters. We combined for the first time the Planck tSZ data with SPT high-ℓ to derive new constraints. Finally, we show how the amplitude of the remaining kinetic SZ power spectrum varies depending on the assumptions made on both tSZ and cosmological parameters. These results show the importance of a thorough modelling of foregrounds in the cosmological analysis of small-scale CMB data. Reliable constraints on cosmological parameters can only be achieved once other significant foregrounds, such as the kinetic SZ and the cosmic infrared background (CIB), are also properly accounted for.

The interstellar dust emission spectrum

Context. Most of the modelling of the interstellar dust infrared emission spectrum is done by assuming some variations around a single-temperature grey-body approximation. For example, the foreground modelling of Planck mission maps involves a single dust temperature, T, along a given line-of-sight with a single emissivity index, β. The two parameters are then fitted and therefore variable from one line-of-sight to the other. Aims. Our aim is to go beyond that modelling in an economical way. Methods. We model the dust spectrum with a temperature distribution around the mean value and show that only the second temperature moment matters. We advocate the use of the temperature logarithm as the proper variable. Results. If the interstellar medium is not too heterogeneous, there is a universal analytical spectrum, which is derived here, that goes beyond the grey-body assumption. We show how the cosmic microwave background radiatively interacts with the dust spectrum (a non-negligible corrective term at millimetre wavelengths). Finally, we construct a universal ladder of discrete temperatures, which gives a minimal and fast description of dust emission spectra as measured by photometric mapping instruments that lends itself to an almost linear fitting. This data modelling can include contributions from the cosmic infrared background fluctuations.