-
Constraining Cosmological Parameters with Needlet Internal Linear Combination Maps II: Likelihood-Free Inference on NILC Power Spectra
Authors:
Kristen M. Surrao,
J. Colin Hill
Abstract:
Standard cosmic microwave background (CMB) analyses constrain cosmological and astrophysical parameters by fitting parametric models to multifrequency power spectra (MFPS). However, such methods do not optimally weight maps in power spectrum (PS) measurements for non-Gaussian CMB foregrounds. We propose needlet internal linear combination (NILC), operating on wavelets with compact support in pixel…
▽ More
Standard cosmic microwave background (CMB) analyses constrain cosmological and astrophysical parameters by fitting parametric models to multifrequency power spectra (MFPS). However, such methods do not optimally weight maps in power spectrum (PS) measurements for non-Gaussian CMB foregrounds. We propose needlet internal linear combination (NILC), operating on wavelets with compact support in pixel and harmonic space, as a weighting scheme to yield more optimal parameter constraints. In a companion paper, we derived an analytic formula for NILC map PS, which is physically insightful but computationally difficult to use in parameter inference pipelines. In this work, we analytically show that fitting parametric templates to MFPS and harmonic ILC PS yields identical parameter constraints when the number of sky components equals or exceeds the number of frequency channels. We numerically show that, for Gaussian random fields, the same holds for NILC PS. This suggests that NILC can reduce parameter error bars in the presence of non-Gaussian fields since it uses non-Gaussian information. As Gaussian likelihoods may be inaccurate, we use likelihood-free inference (LFI) with neural posterior estimation. We show that performing inference with auto- and cross-PS of NILC component maps as summary statistics yields smaller parameter error bars than inference with MFPS. For a model with CMB, an amplified thermal Sunyaev--Zel'dovich (tSZ) signal, and noise, we find a 60% reduction in the area of the 2D 68% confidence region for component amplitude parameters inferred from NILC PS, as compared to inference from MFPS. Primordial $B$-mode searches are a promising application for our new method, as the amplitude of the non-Gaussian dust foreground is known to be larger than a potential signal. Our code is available in https://github.com/kmsurrao/NILC-Inference-Pipeline.
△ Less
Submitted 24 June, 2024;
originally announced June 2024.
-
Accelerated inference on accelerated cosmic expansion: New constraints on axion-like early dark energy with DESI BAO and ACT DR6 CMB lensing
Authors:
Frank J. Qu,
Kristen M. Surrao,
Boris Bolliet,
J. Colin Hill,
Blake D. Sherwin,
Hidde T. Jense
Abstract:
The early dark energy (EDE) extension to $Λ$CDM has been proposed as a candidate scenario to resolve the "Hubble tension". We present new constraints on the EDE model by incorporating new data from the Dark Energy Spectroscopic Instrument (DESI) Baryon Acoustic Oscillation (BAO) survey and CMB lensing measurements from the Atacama Cosmology Telescope (ACT) DR6 and \textit{Planck} NPIPE data. We do…
▽ More
The early dark energy (EDE) extension to $Λ$CDM has been proposed as a candidate scenario to resolve the "Hubble tension". We present new constraints on the EDE model by incorporating new data from the Dark Energy Spectroscopic Instrument (DESI) Baryon Acoustic Oscillation (BAO) survey and CMB lensing measurements from the Atacama Cosmology Telescope (ACT) DR6 and \textit{Planck} NPIPE data. We do not find evidence for EDE. The maximum fractional contribution of EDE to the total energy density is $f_\mathrm{EDE}< 0.091 \; (95\% \; \mathrm{CL} )$ from our baseline combination of \textit{Planck} CMB, CMB lensing, and DESI BAO. Our strongest constraints on EDE come from the combination of \textit{Planck} CMB and CMB lensing alone, yielding $f_\mathrm{EDE}< 0.070 \; (95\% \; \mathrm{CL} )$. We also explore extensions of $Λ$CDM beyond the EDE parameters by treating the total neutrino mass as a free parameter, finding $\sum m_ν< 0.096 \,\, {\rm eV} \; (95\% \; \mathrm{CL} )$ and $f_\mathrm{EDE}< 0.087 \; (95\% \; \mathrm{CL} )$. For the first time in EDE analyses, we perform Bayesian parameter estimation using neural network emulators of cosmological observables, which are on the order of a hundred times faster than full Boltzmann solutions.
△ Less
Submitted 25 April, 2024;
originally announced April 2024.
-
Constraining Cosmological Parameters with Needlet Internal Linear Combination Maps I: Analytic Power Spectrum Formalism
Authors:
Kristen M. Surrao,
J. Colin Hill
Abstract:
The internal linear combination (ILC) method is a popular approach for constructing component-separated maps in cosmic microwave background (CMB) analyses. It optimally combines observed maps at different frequencies to produce an unbiased minimum-variance map of a component. When performed in harmonic space, it is straightforward to analytically compute the contributions of individual sky compone…
▽ More
The internal linear combination (ILC) method is a popular approach for constructing component-separated maps in cosmic microwave background (CMB) analyses. It optimally combines observed maps at different frequencies to produce an unbiased minimum-variance map of a component. When performed in harmonic space, it is straightforward to analytically compute the contributions of individual sky components to the power spectrum of the resulting ILC map. ILC can also be performed on a basis of needlets, spherical wavelets that have compact support in both pixel and harmonic space, capturing both scale-dependent and spatially varying information. However, an analytic understanding of the power spectra of needlet ILC (NILC) component-separated maps, as needed to enable their use in cosmological parameter inference, has remained an outstanding problem. In this paper, we derive the first analytic expression for the power spectra of NILC maps, as well as an expression for the cross-spectrum of a NILC map with an arbitrary second map, in terms of contributions from individual sky components. We validate our result using simulations, finding that it is exact. These results contain useful insights: we explicitly see how NILC power spectra contain information from contaminant fields beyond the two-point level, and we obtain a formalism with which to parameterize NILC power spectra. However, because this parameter dependence is complicated by correlations and higher-point functions of the component maps and weight maps, we find that it is intractable to perform parameter inference using these analytic expressions. Instead, numerical techniques are needed to estimate parameters using NILC maps -- we explore the use of likelihood-free inference with neural posterior estimation in a companion paper. Our code to produce the results in this paper is available in https://github.com/kmsurrao/NILC-PS-Model.
△ Less
Submitted 4 March, 2024;
originally announced March 2024.
-
class_sz I: Overview
Authors:
B. Bolliet,
A. Kusiak,
F. McCarthy,
A. Sabyr,
K. Surrao,
J. C. Hill,
J. Chluba,
S. Ferraro,
B. Hadzhiyska,
D. Han,
J. F. Macías-Pérez,
M. Madhavacheril,
A. Maniyar,
Y. Mehta,
S. Pandey,
E. Schaan,
B. Sherwin,
A. Spurio Mancini,
Í. Zubeldia
Abstract:
class_sz is a versatile and robust code in C and Python that can compute theoretical predictions for a wide range of observables relevant to cross-survey science in the Stage IV era. The code is public at https://github.com/CLASS-SZ/class_sz along with a series of tutorial notebooks (https://github.com/CLASS-SZ/notebooks). It will be presented in full detail in paper II. Here we give a brief overv…
▽ More
class_sz is a versatile and robust code in C and Python that can compute theoretical predictions for a wide range of observables relevant to cross-survey science in the Stage IV era. The code is public at https://github.com/CLASS-SZ/class_sz along with a series of tutorial notebooks (https://github.com/CLASS-SZ/notebooks). It will be presented in full detail in paper II. Here we give a brief overview of key features and usage.
△ Less
Submitted 27 October, 2023;
originally announced October 2023.
-
Enhancing Measurements of the CMB Blackbody Temperature Power Spectrum by Removing CIB and Thermal Sunyaev-Zel'dovich Contamination Using External Galaxy Catalogs
Authors:
Aleksandra Kusiak,
Kristen M. Surrao,
J. Colin Hill
Abstract:
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…
▽ More
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.
△ Less
Submitted 1 November, 2023; v1 submitted 14 March, 2023;
originally announced March 2023.
-
Accurate estimation of angular power spectra for maps with correlated masks
Authors:
Kristen M. Surrao,
Oliver H. E. Philcox,
J. Colin Hill
Abstract:
The widely used MASTER approach for angular power spectrum estimation was developed as a fast $C_{\ell}$ estimator on limited regions of the sky. This method expresses the power spectrum of a masked map ("pseudo-$C_\ell$") in terms of the power spectrum of the unmasked map (the true $C_\ell$) and that of the mask or weight map. However, it is often the case that the map and mask are correlated in…
▽ More
The widely used MASTER approach for angular power spectrum estimation was developed as a fast $C_{\ell}$ estimator on limited regions of the sky. This method expresses the power spectrum of a masked map ("pseudo-$C_\ell$") in terms of the power spectrum of the unmasked map (the true $C_\ell$) and that of the mask or weight map. However, it is often the case that the map and mask are correlated in some way, such as point source masks used in cosmic microwave background (CMB) analyses, which have nonzero correlation with CMB secondary anisotropy fields and other mm-wave sky signals. In such situations, the MASTER approach gives biased results, as it assumes that the unmasked map and mask have zero correlation. While such effects have been discussed before with regard to specific physical models, here we derive a completely general formalism for any case where the map and mask are correlated. We show that our result ("reMASTERed") reconstructs ensemble-averaged pseudo-$C_\ell$ to effectively exact precision, with significant improvements over traditional estimators for cases where the map and mask are correlated. In particular, we obtain an improvement in the mean absolute percent error from 30% with the MASTER result to essentially no error with the reMASTERed result for an integrated Sachs-Wolfe (ISW) field map with a mask built from the thresholded ISW field, and 10% to effectively zero for a Compton-$y$ map combined with an infrared source mask (the latter being directly relevant to actual data analysis). An important consequence of our result is that for maps with correlated masks it is no longer possible to invert a simple equation to obtain the true $C_\ell$ from the pseudo-$C_\ell$. Instead, our result necessitates the use of forward modeling from theory space into the observable domain of the pseudo-$C_\ell$. Our code is publicly available at https://github.com/kmsurrao/reMASTERed.
△ Less
Submitted 25 April, 2023; v1 submitted 10 February, 2023;
originally announced February 2023.
