Enhancing Measurements of the CMB Blackbody Temperature Power Spectrum by Removing CIB and Thermal Sunyaev-Zel'dovich Contamination Using External Galaxy Catalogs

Extracting the CMB blackbody temperature power spectrum -- which is dominated by the primary CMB signal and the kinematic Sunyaev-Zel'dovich (kSZ) effect -- from mm-wave sky maps requires cleaning other sky components. In this work, we develop new methods to use large-scale structure (LSS) tracers to remove cosmic infrared background (CIB) and thermal Sunyaev-Zel'dovich (tSZ) contamination in such measurements. Our methods rely on the fact that LSS tracers are correlated with the CIB and tSZ signals, but their two-point correlations with the CMB and kSZ signals vanish on small scales, thus leaving the CMB blackbody power spectrum unbiased after cleaning. We develop methods analogous to delensing ($\textit{de-CIB}$ or $\textit{de-(CIB+tSZ)}$) to clean CIB and tSZ contaminants using these tracers. We compare these methods to internal linear combination (ILC) methods, including novel approaches that incorporate the tracer maps in the ILC procedure itself, without requiring exact assumptions about the CIB SED. As a concrete example, we use the $\textit{unWISE}$ galaxy samples as tracers. We provide calculations for a combined Simons Observatory and $\textit{Planck}$-like experiment, with our simulated sky model comprising eight frequencies from 93 to 353 GHz. Using $\textit{unWISE}$ tracers, improvements with our methods over current approaches are already non-negligible: we find improvements up to 20% in the kSZ power spectrum signal-to-noise ratio (SNR) when applying the de-CIB method to a tSZ-deprojected ILC map. These gains could be more significant when using additional LSS tracers from current surveys, and will become even larger with future LSS surveys, with improvements in the kSZ power spectrum SNR up to 50%. For the total CMB blackbody power spectrum, these improvements stand at 4% and 7%, respectively. Our code is publicly available at https://github.com/olakusiak/deCIBing.Comment: 35+21 pages, 20+11 figures; code is available at https://github.com/olakusiak/deCIBin

Constraining the galaxy-halo connection of infrared-selected unWISE galaxies with galaxy clustering and galaxy-CMB lensing power spectra

We present the first detailed analysis of the connection between galaxies and their dark matter halos for the unWISE galaxy catalog -- a full-sky, infrared-selected sample built from WISE data, containing over 500 million galaxies. Using unWISE galaxy-galaxy auto-correlation and Planck CMB lensing-galaxy cross-correlation measurements down to 10 arcmin angular scales, we constrain the halo occupation distribution (HOD), a model describing how central and satellite galaxies are distributed within dark matter halos, for three unWISE} galaxy samples at mean redshifts $\bar{z} \approx 0.6$, $1.1$, and $1.5$. We constrain the characteristic minimum halo mass to host a central galaxy, $M_\mathrm{min}^\mathrm{HOD} = 1.83^{+0.41}_{-1.63} \times 10^{12} M_\odot/h$, $5.22^{+0.34}_{-4.80} \times 10^{12} M_\odot/h$, $6.60 ^{+0.30}_{-1.11} \times 10^{13} M_\odot/h$ for the unWISE samples at $\bar{z}\approx 0.6$, $1.1$, and $1.5$, respectively. We find that all three samples are dominated by central galaxies, rather than satellites. Using our constrained HOD models, we infer the effective linear galaxy bias for each unWISE sample, and find that it does not evolve as steeply with redshift as found in previous perturbation-theory-based analyses of these galaxies. We discuss possible sources of systematic uncertainty in our results, the most significant of which is the uncertainty on the galaxy redshift distribution. Our HOD constraints provide a detailed, quantitative understanding of how the unWISE galaxies populate the underlying dark matter halo distribution. These constraints will have a direct impact on future studies employing the unWISE galaxies as a cosmological and astrophysical probe, including measurements of ionized gas thermodynamics and dark matter profiles via Sunyaev-Zel'dovich and lensing cross-correlations

Constraining the Baryon Abundance with the Kinematic Sunyaev-Zel'dovich Effect: Projected-Field Detection Using Planck, WMAP, and unWISE

The kinematic Sunyaev-Zel'dovich (kSZ) effect -- the Doppler boosting of cosmic microwave background (CMB) photons scattering off free electrons with non-zero line-of-sight velocity -- is an excellent probe of the distribution of baryons in the Universe. In this paper, we measure the kSZ effect due to ionized gas traced by infrared-selected galaxies from the \emph{unWISE} catalog. We employ the "projected-field" kSZ estimator, which does not require spectroscopic galaxy redshifts. To suppress non-kSZ foreground signals associated with the galaxies (e.g., dust emission and thermal SZ), this estimator requires cleaned CMB maps, which we obtain from \emph{Planck} and \emph{WMAP} data. Using a new "asymmetric" estimator that combines different foreground-cleaned CMB maps to maximize the signal-to-noise, we measure the kSZ$^2$-galaxy cross-power spectrum for three subsamples of the \emph{unWISE} galaxy catalog, which peak at mean redshifts $z \approx$ 0.6, 1.1, and 1.5, have average halo mass $\sim 1$-$5\times 10^{13}$ $h^{-1} M_{\odot}$, and in total contain over 500 million galaxies. After marginalizing over CMB lensing contributions, we measure the amplitude of the kSZ signal $A_{\rm kSZ^2} = 0.42 \pm 0.31(stat.) \pm 0.14(sys.)$, $5.02 \pm 1.01(stat.) \pm 0.47(sys.)$, and $8.23 \pm 3.23(stat.) \pm 1.60(sys.)$, for the three subsamples, where $A_{\rm kSZ^2} = 1$ corresponds to our fiducial model. The combined kSZ detection S/N $>$ 5. We discuss possible explanations for the excess kSZ signal associated with the $z \approx 1.1$ sample, and show that foreground contamination in the CMB maps is very unlikely to be the cause. Our measurements illustrate clearly that no baryons are missing on large scales at low redshifts.Comment: Matches the version accepted for publication in PRD (https://journals.aps.org/prd/accepted/d907bQ13Oae1fd21e4ff4f13177300acb1c23e88e

The Atacama Cosmology Telescope: Detection of Patchy Screening of the Cosmic Microwave Background

Spatial variations in the cosmic electron density after reionization generate cosmic microwave background anisotropies via Thomson scattering, a process known as the ``patchy screening" effect. In this paper, we propose a new estimator for the patchy screening effect that is designed to mitigate biases from the dominant foreground signals. We use it to measure the cross-correlation between \textit{unWISE} galaxies and patchy screening, the latter measured by the Atacama Cosmology Telescope and \textit{Planck} satellite. We report the first detection of the patchy screening effect, with the statistical significance of the cross-correlation exceeding $7\sigma$. This measurement directly probes the distribution of electrons around these galaxies and provides strong evidence that gas is more extended than the underlying dark matter. By comparing our measurements to electron profiles extracted from simulations, we demonstrate the power of these observations to constrain galaxy evolution models. Requiring only the 2D positions of objects and no individual redshifts or velocity estimates, this approach is complementary to existing gas probes, such as those based on the kinetic Sunyaev-Zeldovich effect.Comment: See Schutt et al for a detailed comparison of patchy screening estimators. 17 pages with 8 figure

The Atacama Cosmology Telescope: A Measurement of the DR6 CMB Lensing Power Spectrum and its Implications for Structure Growth

We present new measurements of cosmic microwave background (CMB) lensing over $9400$ sq. deg. of the sky. These lensing measurements are derived from the Atacama Cosmology Telescope (ACT) Data Release 6 (DR6) CMB dataset, which consists of five seasons of ACT CMB temperature and polarization observations. We determine the amplitude of the CMB lensing power spectrum at $2.3\%$ precision ($43\sigma$ significance) using a novel pipeline that minimizes sensitivity to foregrounds and to noise properties. To ensure our results are robust, we analyze an extensive set of null tests, consistency tests, and systematic error estimates and employ a blinded analysis framework. The baseline spectrum is well fit by a lensing amplitude of $A_{\mathrm{lens}}=1.013\pm0.023$ relative to the Planck 2018 CMB power spectra best-fit $\Lambda$CDM model and $A_{\mathrm{lens}}=1.005\pm0.023$ relative to the $\text{ACT DR4} + \text{WMAP}$ best-fit model. From our lensing power spectrum measurement, we derive constraints on the parameter combination $S^{\mathrm{CMBL}}_8 \equiv \sigma_8 \left({\Omega_m}/{0.3}\right)^{0.25}$ of $S^{\mathrm{CMBL}}_8= 0.818\pm0.022$ from ACT DR6 CMB lensing alone and $S^{\mathrm{CMBL}}_8= 0.813\pm0.018$ when combining ACT DR6 and Planck NPIPE CMB lensing power spectra. These results are in excellent agreement with $\Lambda$CDM model constraints from Planck or $\text{ACT DR4} + \text{WMAP}$ CMB power spectrum measurements. Our lensing measurements from redshifts $z\sim0.5$--$5$ are thus fully consistent with $\Lambda$CDM structure growth predictions based on CMB anisotropies probing primarily $z\sim1100$. We find no evidence for a suppression of the amplitude of cosmic structure at low redshiftsComment: 45+21 pages, 50 figures. Prepared for submission to ApJ. Also see companion papers Madhavacheril et al and MacCrann et a

The Atacama Cosmology Telescope: High-resolution component-separated maps across one-third of the sky

Observations of the millimeter sky contain valuable information on a number of signals, including the blackbody cosmic microwave background (CMB), Galactic emissions, and the Compton-$y$ distortion due to the thermal Sunyaev-Zel'dovich (tSZ) effect. Extracting new insight into cosmological and astrophysical questions often requires combining multi-wavelength observations to spectrally isolate one component. In this work, we present a new arcminute-resolution Compton-$y$ map, which traces out the line-of-sight-integrated electron pressure, as well as maps of the CMB in intensity and E-mode polarization, across a third of the sky (around 13,000 sq.~deg.). We produce these through a joint analysis of data from the Atacama Cosmology Telescope (ACT) Data Release 4 and 6 at frequencies of roughly 93, 148, and 225 GHz, together with data from the \textit{Planck} satellite at frequencies between 30 GHz and 545 GHz. We present detailed verification of an internal linear combination pipeline implemented in a needlet frame that allows us to efficiently suppress Galactic contamination and account for spatial variations in the ACT instrument noise. These maps provide a significant advance, in noise levels and resolution, over the existing \textit{Planck} component-separated maps and will enable a host of science goals including studies of cluster and galaxy astrophysics, inferences of the cosmic velocity field, primordial non-Gaussianity searches, and gravitational lensing reconstruction of the CMB.Comment: The Compton-y map and associated products will be made publicly available upon publication of the paper. The CMB T and E mode maps will be made available when the DR6 maps are made publi

The Atacama Cosmology Telescope: DR6 Gravitational Lensing Map and Cosmological Parameters

We present cosmological constraints from a gravitational lensing mass map covering 9400 sq. deg. reconstructed from CMB measurements made by the Atacama Cosmology Telescope (ACT) from 2017 to 2021. In combination with BAO measurements (from SDSS and 6dF), we obtain the amplitude of matter fluctuations $\sigma_8 = 0.819 \pm 0.015$ at 1.8% precision, $S_8\equiv\sigma_8({\Omega_{\rm m}}/0.3)^{0.5}=0.840\pm0.028$ and the Hubble constant $H_0= (68.3 \pm 1.1)\, \text{km}\,\text{s}^{-1}\,\text{Mpc}^{-1}$ at 1.6% precision. A joint constraint with CMB lensing measured by the Planck satellite yields even more precise values: $\sigma_8 = 0.812 \pm 0.013$, $S_8\equiv\sigma_8({\Omega_{\rm m}}/0.3)^{0.5}=0.831\pm0.023$ and $H_0= (68.1 \pm 1.0)\, \text{km}\,\text{s}^{-1}\,\text{Mpc}^{-1}$. These measurements agree well with $\Lambda$CDM-model extrapolations from the CMB anisotropies measured by Planck. To compare these constraints to those from the KiDS, DES, and HSC galaxy surveys, we revisit those data sets with a uniform set of assumptions, and find $S_8$ from all three surveys are lower than that from ACT+Planck lensing by varying levels ranging from 1.7-2.1$\sigma$. These results motivate further measurements and comparison, not just between the CMB anisotropies and galaxy lensing, but also between CMB lensing probing $z\sim 0.5-5$ on mostly-linear scales and galaxy lensing at $z\sim 0.5$ on smaller scales. We combine our CMB lensing measurements with CMB anisotropies to constrain extensions of $\Lambda$CDM, limiting the sum of the neutrino masses to $\sum m_{\nu} < 0.12$ eV (95% c.l.), for example. Our results provide independent confirmation that the universe is spatially flat, conforms with general relativity, and is described remarkably well by the $\Lambda$CDM model, while paving a promising path for neutrino physics with gravitational lensing from upcoming ground-based CMB surveys.Comment: 30 pages, 16 figures, prepared for submission to ApJ. Cosmological likelihood data is here: https://lambda.gsfc.nasa.gov/product/act/actadv_prod_table.html ; likelihood software is here: https://github.com/ACTCollaboration/act_dr6_lenslike . Also see companion papers Qu et al and MacCrann et al. Mass maps will be released when papers are publishe

The Atacama Cosmology Telescope: A measurement of the DR6 CMB lensing power spectrum and its implications for structure growth

We present new measurements of cosmic microwave background (CMB) lensing over 9400 deg2 of the sky. These lensing measurements are derived from the Atacama Cosmology Telescope (ACT) Data Release 6 (DR6) CMB data set, which consists of five seasons of ACT CMB temperature and polarization observations. We determine the amplitude of the CMB lensing power spectrum at 2.3% precision (43σ significance) using a novel pipeline that minimizes sensitivity to foregrounds and to noise properties. To ensure that our results are robust, we analyze an extensive set of null tests, consistency tests, and systematic error estimates and employ a blinded analysis framework. Our CMB lensing power spectrum measurement provides constraints on the amplitude of cosmic structure that do not depend on Planck or galaxy survey data, thus giving independent information about large-scale structure growth and potential tensions in structure measurements. The baseline spectrum is well fit by a lensing amplitude of A lens = 1.013 ± 0.023 relative to the Planck 2018 CMB power spectra best-fit ΛCDM model and A lens = 1.005 ± 0.023 relative to the ACT DR4 + WMAP best-fit model. From our lensing power spectrum measurement, we derive constraints on the parameter combination S8CMBL≡σ8Ωm/0.30.25 of S8CMBL=0.818±0.022 from ACT DR6 CMB lensing alone and S8CMBL=0.813±0.018 when combining ACT DR6 and Planck NPIPE CMB lensing power spectra. These results are in excellent agreement with ΛCDM model constraints from Planck or ACT DR4 + WMAP CMB power spectrum measurements. Our lensing measurements from redshifts z ∼ 0.5–5 are thus fully consistent with ΛCDM structure growth predictions based on CMB anisotropies probing primarily z ∼ 1100. We find no evidence for a suppression of the amplitude of cosmic structure at low redshifts

The Atacama Cosmology Telescope: DR6 gravitational lensing map and cosmological parameters

We present cosmological constraints from a gravitational lensing mass map covering 9400 deg2 reconstructed from measurements of the cosmic microwave background (CMB) made by the Atacama Cosmology Telescope (ACT) from 2017 to 2021. In combination with measurements of baryon acoustic oscillations and big bang nucleosynthesis, we obtain the clustering amplitude σ 8 = 0.819 ± 0.015 at 1.8% precision, S8≡σ8(Ωm/0.3)0.5=0.840±0.028 , and the Hubble constant H 0 = (68.3 ± 1.1) km s−1 Mpc−1 at 1.6% precision. A joint constraint with Planck CMB lensing yields σ 8 = 0.812 ± 0.013, S8≡σ8(Ωm/0.3)0.5=0.831±0.023 , and H 0 = (68.1 ± 1.0) km s−1 Mpc−1. These measurements agree with ΛCDM extrapolations from the CMB anisotropies measured by Planck. We revisit constraints from the KiDS, DES, and HSC galaxy surveys with a uniform set of assumptions and find that S 8 from all three are lower than that from ACT+Planck lensing by levels ranging from 1.7σ to 2.1σ. This motivates further measurements and comparison, not just between the CMB anisotropies and galaxy lensing but also between CMB lensing probing z ∼ 0.5–5 on mostly linear scales and galaxy lensing at z ∼ 0.5 on smaller scales. We combine with CMB anisotropies to constrain extensions of ΛCDM, limiting neutrino masses to ∑m ν < 0.13 eV (95% c.l.), for example. We describe the mass map and related data products that will enable a wide array of cross-correlation science. Our results provide independent confirmation that the universe is spatially flat, conforms with general relativity, and is described remarkably well by the ΛCDM model, while paving a promising path for neutrino physics with lensing from upcoming ground-based CMB surveys