Search

The thermal Sunyaev-Zel'dovich effect (tSZ) is a sensitive probe of cosmology, as it traces the abundance of galaxy clusters and groups in the late-time Universe. Upcoming cosmic microwave background experiments such as the Simons Observatory (SO) and CMB-S4 will provide low-noise and high-resolution component-separated tSZ maps covering a large sky fraction. The tSZ signal is highly non-Gaussian; therefore, higher-order statistics are needed to optimally extract information from these maps. In this work, we study the cosmological constraining power of several tSZ statistics -- Minkowski functionals (MFs), peaks, minima, and moments -- that have yielded promising results in capturing non-Gaussian information from other cosmological data. Using a large suite of halo-model-based tSZ simulations with varying $\Omega_{c}$ and $\sigma_{8}$ (154 cosmologies and over $800, 000$ maps, each $10.5\times10.5$ deg$^{2}$), we show that by combining these observables, we can achieve $\approx 29\times$ tighter constraints compared to using the tSZ power spectrum alone in an idealized noiseless case, with the MFs dominating the constraints. We show that much of the MF constraining power arises from halos below the detection threshold of cluster surveys, suggesting promising synergies with cluster-count analyses. Finally, we demonstrate that these statistics have the potential to deliver tight constraints even in the presence of noise. For example, using post-component-separation tSZ noise expected for SO, we obtain $\approx1.6\times$ and $\approx1.8\times$ tighter constraints than the power spectrum with MFs and all statistics combined, respectively. We show that the constraints from MFs approach the noiseless case for white-noise levels $\lesssim 1 \,\, \mu$K-arcmin., Comment: 30 pages, 23 figures

The high electron temperature in galaxy clusters ($>1\,$keV or $>10^7\,$K) leads to corrections at the level of a few percent in their thermal Sunyaev-Zeldovich effect signatures. Both the size and frequency dependence of these corrections, which are known as relativistic temperature corrections, depend upon the temperature of the objects. In this work we exploit this effect to measure the average temperature of a stack of Compton-$y$ selected clusters. Specifically, we apply the "spectroscopic method" and search for the temperature that best fits the clusters' signal measured at frequencies from 30 to 545 GHz by the Atacama Cosmology Telescope and Planck satellite. We measure the average temperature of clusters detected in the ACT maps to be $8.5\pm 2.4\,$keV, with an additional systematic error of comparable amplitude dominated by passband uncertainty. Upcoming surveys, such as the Simons Observatory and CMB-S4, have the potential to dramatically improve upon these measurements and thereby enable precision studies of cluster temperatures with millimeter observations. The key challenge for future observations will be mitigating instrumental systematic effects, which already limit this analysis., Comment: 21 pages with 17 figures

The kinematic Sunyaev--Zel'dovich (kSZ) effect induces a non-zero density-density-temperature bispectrum, which we can use to reconstruct the large-scale velocity field from a combination of cosmic microwave background (CMB) and galaxy density measurements, in a procedure known as ``kSZ velocity reconstruction''. This method has been forecast to constrain large-scale modes with future galaxy and CMB surveys, improving their measurement beyond what is possible with the galaxy surveys alone. Such measurements will enable tighter constraints on large-scale signals such as primordial non-Gaussianity, deviations from homogeneity, and modified gravity. In this work, we demonstrate a statistically significant measurement of kSZ velocity reconstruction for the first time, by applying quadratic estimators to the combination of the ACT DR6 CMB+kSZ map and the DESI LRG galaxies (with photometric redshifts) in order to reconstruct the velocity field. We do so using a formalism appropriate for the 2-dimensional projected galaxy fields that we use, which naturally incorporates the curved-sky effects important on the largest scales. We find evidence for the signal by cross-correlating with an external estimate of the velocity field from the spectroscopic BOSS survey and rejecting the null (no-kSZ) hypothesis at $3.8\sigma$. Our work presents a first step towards the use of this observable for cosmological analyses., Comment: 16 pages (main)+5 pages (Appendix); 13 figures (main) + 8 figures (appendix)

Deviations of the cosmic microwave background (CMB) energy spectrum from a perfect blackbody uniquely probe a wide range of physics, ranging from fundamental physics in the primordial Universe ($\mu$-distortion) to late-time baryonic feedback processes (y-distortion). While the y-distortion can be detected with a moderate increase in sensitivity over that of COBE/FIRAS, the $\Lambda$CDM-predicted $\mu$-distortion is roughly two orders of magnitude smaller and requires substantial improvements, with foregrounds presenting a serious obstacle. Within the standard model, the dominant contribution to $\mu$ arises from energy injected via Silk damping, yielding sensitivity to the primordial power spectrum at wavenumbers $k \approx 1-10^{4}$ Mpc$^{-1}$. Here, we present a new instrument concept, SPECTER, with the goal of robustly detecting $\mu$. The instrument technology is similar to that of LiteBIRD, but with an absolute temperature calibration system. Using a Fisher approach, we optimize the instrument's configuration to target $\mu$ while robustly marginalizing over foreground contaminants. Unlike Fourier-transform-spectrometer-based designs, the specific bands and their individual sensitivities can be independently set in this instrument, allowing significant flexibility. We forecast SPECTER to observe the $\Lambda$CDM-predicted $\mu$-distortion at $\approx 5\sigma$ (10$\sigma$) assuming an observation time of 1 (4) year(s) (corresponding to mission duration of 2 (8) years), after foreground marginalization. Our optimized configuration includes 16 bands spanning 1-2000 GHz with degree-scale angular resolution at 150 GHz and 1046 total detectors. SPECTER will additionally measure the y-distortion at sub-percent precision and its relativistic correction at percent-level precision, yielding tight constraints on the total thermal energy and mean temperature of ionized gas., Comment: 26 pages, 12 figures

The resonant conversion of cosmic microwave background (CMB) photons into axions within large-scale structure induces an anisotropic spectral distortion in CMB temperature maps. Applying state-of-the-art foreground cleaning techniques to $\textit{Planck}$ CMB observations, we construct maps of axion-induced "patchy screening" of the CMB. We cross-correlate these maps with data from the $\textit{unWISE}$ galaxy survey and find no evidence of axions. We constrain the axion-photon coupling, $g_{a\gamma\gamma} \lesssim 2 \times 10^{-12}~{\rm GeV}^{-1}$, at the 95% confidence level for axion masses in the range $10^{-13}~{\rm eV} \lesssim m_a \lesssim 10^{-12}~{\rm eV}$. These constraints are competitive with the tightest astrophysical axion limits in this mass range and are inferred from robust population-level statistics, which makes them complementary to existing searches that rely on modeling of individual systems., Comment: 5+15 pages; 3+15 figures

We present a joint analysis of the CMB lensing power spectra measured from the Data Release 6 of the Atacama Cosmology Telescope and Planck PR4, cross-correlations between the ACT and Planck lensing reconstruction and galaxy clustering from unWISE, and the unWISE clustering auto-spectrum. We obtain 1.5% constraints on the matter density fluctuations at late times parametrised by the best constrained parameter combination $S_8^{\rm 3x2pt}\equiv\sigma_8 (\Omega_m/0.3)^{0.4}=0.815\pm0.012$. The commonly used $S_8\equiv\sigma_8 (\Omega_m/0.3)^{0.5}$ parameter is constrained to $S_8=0.816\pm0.015$. In combination with baryon acoustic oscillation (BAO) measurements we find $\sigma_8=0.815\pm 0.012$. We also present sound-horizon-independent estimates of the present day Hubble rate of $H_0=66.4^{+3.2}_{-3.7} \,\mathrm{km}\,\mathrm{s}^{-1}\mathrm{Mpc}^{-1}$ from our large scale structure data alone and $H_0=64.3^{+2.1}_{-2.4}\,\mathrm{km}\,\mathrm{s}^{-1}\mathrm{Mpc}^{-1}$ in combination with uncalibrated supernovae from Pantheon+. Using parametric estimates of the evolution of matter density fluctuations, we place constraints on cosmic structure in a range of high redshifts typically inaccessible with cross-correlation analyses. Combining lensing cross- and auto-correlations, we derive a 3.3% constraint on the integrated matter density fluctuations above $z=2.4$, one of the tightest constraints in this redshift range and fully consistent with a $\Lambda$CDM model fit to the primary CMB from Planck. Combining with primary CMB observations and using the extended low redshift coverage of these combined data sets we derive constraints on a variety of extensions to the $\Lambda$CDM model including massive neutrinos, spatial curvature, and dark energy. We find in flat $\Lambda$CDM $\sum m_\nu<0.12$ eV at 95% confidence using the LSS data, BAO measurements from SDSS and primary CMB observations., Comment: 30 pages, 15 figures, to be submitted to PRD, comments welcome

Massive particles produced during inflation impact soft limits of primordial correlators. Searching for these signatures presents an exciting opportunity to uncover the particle spectrum in the inflationary epoch. We present non-perturbative methods to constrain intermediate-mass scalars ($0\leq m/H<3/2$, where $H$ is the inflationary Hubble scale) produced during inflation, which give rise to a power-law scaling in the squeezed primordial bispectrum. Exploiting the large-scale structure consistency relations and the separate universe approach, we derive models for the late-time squeezed matter bispectrum and collapsed matter trispectrum sourced by these fields. To validate our models, we run $N$-body simulations with the "Cosmological Collider" squeezed bispectrum for two different particle masses. Our models yield unbiased constraints on the amplitude of non-Gaussianity, $f_{\rm NL}^{\Delta}$, from the squeezed bispectrum and collapsed trispectrum deep into the non-linear regime ($k_{\rm max}\approx 2~h/{\rm Mpc}$ at $z=0$). We assess the information content of these summary statistics, emphasizing the importance of sample variance cancellation in the matter sector. We also study the scale-dependent halo bias in our simulations. For mass-selected halos, the non-Gaussian bias estimated from our simulations agrees with predictions based on (i) separate universe simulations and (ii) universal mass functions. With further work, these results can be used to search for inflationary massive particle production with upcoming galaxy surveys., Comment: 25 pages, 11 figures; comments welcome

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., Comment: 22+13 pages, 6+9 figures; companion paper to arXiv:2403.02261; code available at https://github.com/kmsurrao/NILC-Inference-Pipeline

Dark photons that kinetically mix with the Standard Model photon give rise to new spectral anisotropies (patchy dark screening) in the cosmic microwave background (CMB) due to conversion of photons to dark photons within large-scale structure. We utilize predictions for this patchy dark screening signal to provide the tightest constraints to date on the dark photon kinetic mixing parameter ($\varepsilon \lesssim 4\times 10^{-8}$ (95% confidence level)) over the mass range $10^{-13} \,\, {\rm eV} \lesssim m_{{A^\prime}} \lesssim 10^{-11}$ eV, almost an order of magnitude stronger than previous limits, by applying state-of-the-art component separation techniques to the cross-correlation of $\textit{Planck}$ CMB and $\textit{unWISE}$ galaxy survey data., Comment: 7+12 pages, 3+13 figures. Data products available at https://users.flatironinstitute.org/~fmccarthy/dark_photon_screening_maps/ V2 only has minor changes to these comments

Non-adiabatic production of massive particles is a generic feature of many inflationary mechanisms. If sufficiently massive, these particles can leave features in the cosmic microwave background (CMB) that are not well-captured by traditional correlation function analyses. We consider a scenario in which particle production occurs only in a narrow time-interval during inflation, eventually leading to CMB hot- or coldspots with characteristic shapes and sizes. Searching for such features in CMB data is analogous to searching for late-Universe hot- or coldspots, such as those due to the thermal Sunyaev-Zel'dovich (tSZ) effect. Exploiting this data-analysis parallel, we perform a search for particle-production hotspots in the Planck PR4 temperature dataset, which we implement via a matched-filter analysis. Our pipeline is validated on synthetic observations and found to yield unbiased constraints on sufficiently large hotspots across $\approx 60\%$ of the sky. After removing point sources and tSZ clusters, we find no evidence for new physics and place novel bounds on the coupling between the inflaton and massive particles. These bounds are strongest for larger hotspots, produced early in inflation, whilst sensitivity to smaller hotspots is limited by noise and beam effects. Through such methods we can constrain particles with masses $\mathcal{O}(100)$ times larger than the inflationary Hubble scale, which represents possibly the highest energies ever directly probed with observational data., Comment: 12 pages, 7 figures, submitted to Phys. Rev. D

Patchy reionization generates kinematic Sunyaev-Zeldovich (kSZ) anisotropies in the cosmic microwave background (CMB). Large-scale velocity perturbations along the line of sight modulate the small-scale kSZ power spectrum, leading to a trispectrum (or four-point function) in the CMB that depends on the physics of reionization. We investigate the challenges in detecting this trispectrum and use tools developed for CMB lensing, such as realization-dependent bias subtraction and cross-correlation based estimators, to counter uncertainties in the instrumental noise and assumed CMB power spectrum. We also find that both lensing and extragalactic foregrounds can impart larger trispectrum contributions than the reionization kSZ signal. We present a range of mitigation methods for both of these sources of contamination, validated on microwave-sky simulations. We use ACT DR6 and Planck data to calculate an upper limit on the reionization kSZ trispectrum from a measurement dominated by foregrounds. The upper limit is about 50 times the signal predicted from recent simulations., Comment: Measurements and covariances will be made public upon publication

The early dark energy (EDE) extension to $\Lambda$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 $\Lambda$CDM beyond the EDE parameters by treating the total neutrino mass as a free parameter, finding $\sum m_\nu < 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., Comment: 8+7 pages, 3+6 figures

Analysis of microwave sky signals, such as the cosmic microwave background, often requires component separation with multi-frequency methods, where different signals are isolated by their frequency behaviors. Many so-called "blind" methods, such as the internal linear combination (ILC), make minimal assumptions about the spatial distribution of the signal or contaminants, and only assume knowledge of the frequency dependence of the signal. The ILC is a minimum-variance linear combination of the measured frequency maps. In the case of Gaussian, statistically isotropic fields, this is the optimal linear combination, as the variance is the only statistic of interest. However, in many cases the signal we wish to isolate, or the foregrounds we wish to remove, are non-Gaussian and/or statistically anisotropic (in particular for Galactic foregrounds). In such cases, it is possible that machine learning (ML) techniques can be used to exploit the non-Gaussian features of the foregrounds and thereby improve component separation. However, many ML techniques require the use of complex, difficult-to-interpret operations on the data. We propose a hybrid method whereby we train an ML model using only combinations of the data that $\textit{do not contain the signal}$, and combine the resulting ML-predicted foreground estimate with the ILC solution to reduce the error from the ILC. We demonstrate our methods on simulations of extragalactic temperature and Galactic polarization foregrounds, and show that our ML model can exploit non-Gaussian features, such as point sources and spatially-varying spectral indices, to produce lower-variance maps than ILC - eg, reducing the variance of the B-mode residual by factors of up to 5 - while preserving the signal of interest in an unbiased manner. Moreover, we often find improved performance when applying our model to foreground models on which it was not trained., Comment: 22 pages, 13 figures. v2 has some updated references

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., Comment: 19 pages, 3 figures; code available at https://github.com/kmsurrao/NILC-PS-Model

The predicted present-day amplitude of matter fluctuations based on cosmic microwave background (CMB) anisotropy data has sometimes been found discrepant with more direct measurements of late-time structure. This has motivated many extensions to the standard cosmological model, including kinetic interactions between dark matter and dark energy that introduce a drag force slowing the growth of structure at late times. Exploring this scenario, we develop a model for quasi-linear scales in the matter power spectrum by calculating the critical overdensity in the presence of this interaction and a varying dark energy equation of state. We explicitly avoid modeling or interpretation of data on non-linear scales in this model (such as use of $\Lambda$CDM-calibrated priors), which would require numerical simulations. We find that the presence of the drag force hinders halo formation, thus increasing the deviation from $\Lambda$CDM in the mildly non-linear regime. We use CMB lensing observations from the sixth data release of the Atacama Cosmology Telescope up to $L=1250$ (in combination with Planck, Sloan Digital Sky Survey, and 6dFGS data) to derive the strongest constraints to date on the amplitude of the drag term, finding the dimensionless interaction strength $\Gamma_\mathrm{DMDE}/(H_0\rho_\mathrm{c})<0.831\; (2.81)$ at the 68\% (95\%) confidence level. The inclusion of non-linear corrections improves our constraints by about 35\% compared to linear theory. Our results do not exclude the best-fit values of $\Gamma_\mathrm{DMDE}$ found in previous studies using information from galaxy weak lensing, though we find no statistical preference for the dark matter-dark energy kinetic interactions over $\Lambda$CDM. We implement our model in a publicly available fork of the Boltzmann code CLASS at https://github.com/fmccarthy/Class_DMDE., Comment: Comments welcome

We introduce GODMAX (Gas thermODynamics and Matter distribution using jAX), a novel code designed to calculate correlations between the cosmological matter distribution and various gas thermodynamic quantities. Utilizing the extensive ANTILLES suite of 200 hydrodynamical simulations with a diverse range of baryonic feedback strengths, we jointly fit the 3D profiles of total matter distribution, electron density, and electron pressure across various halo masses and redshifts. By accommodating significant variations in gas profiles expected due to baryonic feedback, solving exact hydrostatic equilibrium equation and offering flexible modeling of non-thermal pressure support, GODMAX has the capability to jointly fit all these profiles within the measurement uncertainties. This advancement enables, for the first time, robust joint analyses of multiple cosmic probes, including the kinetic and thermal Sunyaev-Zel'dovich effect, weak lensing, and X-ray observations. Furthermore, the model accurately captures correlations between the total matter power suppression due to baryonic feedback and local average thermodynamic quantities, such as the baryon fraction and integrated tSZ effect, in high-mass halos, aligning with observations from hydrodynamical simulations. Looking ahead, we forecast the expected constraints on cosmological and baryonic parameters from upcoming weak lensing catalogs from the LSST and tSZ maps from the Simons Observatory. This analysis underscores the importance of cross-correlations between weak lensing and tSZ in enhancing parameter constraints by resolving major systematic uncertainties due to baryonic physics. The GODMAX code leverages the JAX library, resulting in a fully differentiable halo model with native GPU compilation support., Comment: 17 pages, 9 figures

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 figures

We review and update constraints on the Early Dark Energy (EDE) model from cosmological data sets, in particular Planck PR3 and PR4 cosmic microwave background (CMB) data and large-scale structure (LSS) data sets including galaxy clustering and weak lensing data from the Dark Energy Survey, Subaru Hyper Suprime-Cam, and KiDS+VIKING-450, as well as BOSS/eBOSS galaxy clustering and Lyman-$\alpha$ forest data. We detail the fit to CMB data, and perform the first analyses of EDE using the CAMSPEC and Hillipop likelihoods for Planck CMB data, rather than Plik, both of which yield a tighter upper bound on the allowed EDE fraction than that found with Plik. We then supplement CMB data with large-scale structure data in a series of new analyses. All these analyses are concordant in their Bayesian preference for $\Lambda$CDM over EDE, as indicated by marginalized posterior distributions. We perform a series of tests of the impact of priors in these results, and compare with frequentist analyses based on the profile likelihood, finding qualitative agreement with the Bayesian results. All these tests suggest prior volume effects are not a determining factor in analyses of EDE. This work provides both a review of existing constraints and several new analyses., Comment: 59 pages, 23 figures, 11 tables, Invited review for International Journal of Modern Physics D

We introduce a non-perturbative method to constrain the amplitude of local-type primordial non-Gaussianity ($f_{\rm NL}$) using squeezed configurations of the CMB lensing convergence and cosmic shear bispectra. First, we use cosmological consistency relations to derive a model for the squeezed limit of angular auto- and cross-bispectra of lensing convergence fields in the presence of $f_{\rm NL}$. Using this model, we perform a Fisher forecast with specifications expected for upcoming CMB lensing measurements from the Simons Observatory and CMB-S4, as well as cosmic shear measurements from a Rubin LSST/Euclid-like experiment. Assuming a minimum multipole $\ell_{\rm min}=10$ and maximum multipole $\ell_{\rm max}=1400$, we forecast $\sigma_{f_{\rm NL}}=175$ ($95$) for Simons Observatory (CMB-S4). Our forecasts improve considerably for an LSST/Euclid-like cosmic shear experiment with three tomographic bins and $\ell_{\rm min}=10$ and $\ell_{\rm max}=1400$ ($5000$) with $\sigma_{f_{\rm NL}}=31$ ($16$). A joint analysis of CMB-S4 lensing and LSST/Euclid-like shear yields little gain over the shear-only forecasts; however, we show that a joint analysis could be useful if the CMB lensing convergence can be reliably reconstructed at larger angular scales than the shear field. The method presented in this work is a novel and robust technique to constrain local primordial non-Gaussianity from upcoming large-scale structure surveys that is completely independent of the galaxy field (and therefore any nuisance parameters such as $b_\phi$), thus complementing existing techniques to constrain $f_{\rm NL}$ using the scale-dependent halo bias., Comment: 14 pages, 5 figures, comments welcome!

We present tomographic measurements of structure growth using cross-correlations of Atacama Cosmology Telescope (ACT) DR6 and Planck CMB lensing maps with the unWISE Blue and Green galaxy samples, which span the redshift ranges $0.2 \lesssim z \lesssim 1.1$ and $0.3 \lesssim z \lesssim 1.8$, respectively. We improve on prior unWISE cross-correlations not just by making use of the new, high-precision ACT DR6 lensing maps, but also by including additional spectroscopic data for redshift calibration and by analysing our measurements with a more flexible theoretical model. An extensive suite of systematic and null tests within a blind analysis framework ensures that our results are robust. We determine the amplitude of matter fluctuations at low redshifts ($z\simeq 0.2-1.6$), finding $S_8 \equiv \sigma_8 (\Omega_m / 0.3)^{0.5} = 0.813 \pm 0.021$ using the ACT cross-correlation alone and $S_8 = 0.810 \pm 0.015$ with a combination of Planck and ACT cross-correlations; these measurements are fully consistent with the predictions from primary CMB measurements assuming standard structure growth. The addition of Baryon Acoustic Oscillation data breaks the degeneracy between $\sigma_8$ and $\Omega_m$, allowing us to measure $\sigma_8 = 0.813 \pm 0.020$ from the cross-correlation of unWISE with ACT and $\sigma_8 = 0.813\pm 0.015$ from the combination of cross-correlations with ACT and Planck. These results also agree with the expectations from primary CMB extrapolations in $\Lambda$CDM cosmology; the consistency of $\sigma_8$ derived from our two redshift samples at $z \sim 0.6$ and $1.1$ provides a further check of our cosmological model. Our results suggest that structure formation on linear scales is well described by $\Lambda$CDM even down to low redshifts $z\lesssim 1$., Comment: 73 pages (incl. 31 pages of appendices), 52 figures, 16 tables, published in ApJ. Watch G. S. Farren and A. Krolewski discuss the analysis and results under https://cosmologytalks.com/2023/09/11/act-unwise

The recent detection of a stochastic gravitational wave background (SGWB) at nanohertz frequencies by pulsar timing arrays (PTAs) has sparked a flurry of interest. Beyond the standard interpretation that the progenitor is a network of supermassive black hole binaries, many exotic models have also been proposed, some of which can potentially offer a better fit to the data. We explore how the various connections between gravitational waves and CMB spectral distortions can be leveraged to help determine whether a SGWB was generated primordially or astrophysically. To this end, we present updated $k$-space window functions which can be used for distortion parameter estimation on enhancements to the primordial scalar power spectrum. These same enhancements can also source gravitational waves (GWs) directly at second order in perturbation theory, so-called scalar-induced GWs (SIGWs), and indirectly through the formation of primordial black holes (PBHs). We perform a mapping of scalar power spectrum constraints into limits on the GW parameter space of SIGWs for $\delta$-function features. We highlight that broader features in the scalar spectrum can explain the PTA results while simultaneously producing a spectral distortion (SD) within reach of future experiments. We additionally update PBH constraints from $\mu$- and $y$-type spectral distortions. Refined treatments of the distortion window functions widen existing SD constraints, and we find that a future CMB spectrometer could play a pivotal role in unraveling the origin of GWs imprinted at or below CMB anisotropy scales., Comment: 15 pages, 11 figures, comments welcome

We use the full-mission Planck PR4 data to measure the CMB lensing convergence ($\kappa$)--thermal Sunyaev-Zel'dovich (tSZ, $y$) cross-correlation, $C_\ell^{y\kappa}$. This is only the second measurement to date of this signal, following Hill \& Spergel (2014). We perform the measurement using foreground-cleaned tSZ maps built from the PR4 frequency maps via a tailored needlet internal linear combination (NILC) code in our companion paper, in combination with the Planck PR4 $\kappa$ maps and various systematic-mitigated PR3 $\kappa$ maps. A serious systematic is the residual cosmic infrared background (CIB) in the tSZ map, as the high CIB--$\kappa$ correlation can significantly bias the inferred tSZ--$\kappa$ cross-correlation. We mitigate this by deprojecting the CIB in our NILC, using a moment-deprojection approach to avoid leakage due to incorrect modelling of the CIB frequency dependence. We validate our method on mm-sky simulations. We fit a theoretical halo model to our measurement, finding a best-fit amplitude of $A=0.82\pm0.21$ (for the highest signal-to-noise PR4 $\kappa$ map) or $A=0.56\pm0.24$ (for a PR3 $\kappa$ map built from a tSZ-deprojected CMB map), indicating that the data are consistent with our model within $\sim 1$-$2\sigma$. Although our error bars are similar to those of the 2014 measurement, our method is significantly more robust to CIB contamination. Our moment-deprojection approach lays the foundation for future measurements of this signal with higher signal-to-noise maps from ground-based telescopes, which will precisely probe the astrophysics of the intracluster medium of galaxy groups and clusters in the intermediate-mass ($M\sim 10^{13} -10^{14} h^{-1} M_\odot$), high-$z$ ($z<\sim1.5$, c.f. $z<\sim0.8$ for the tSZ auto-power signal) regime, as well as CIB-decontaminated measurements of tSZ cross-correlations with other large-scale structure probes., Comment: 32 pages, 17 Figures. Companion paper to arXiv:2307.01043 . Public code is available to reproduce theory calculations at https://github.com/fmccarthy/tSZ_kappa_halomodel_classsz

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 public

We use the full-mission $\textit{Planck}$ PR4 data to construct maps of the thermal Sunyaev$--$Zel'dovich effect (Compton-$y$ parameter) in our Universe. To do so, we implement a custom needlet internal linear combination (NILC) pipeline in a Python package, $\texttt{pyilc}$, which we make publicly available. We publicly release our Compton-$y$ maps, which we construct using various constrained ILC ("deprojection") options in order to minimize contamination from the cosmic infrared background (CIB) in the reconstructed signal. In particular, we use a moment-based deprojection which minimizes sensitivity to the assumed frequency dependence of the CIB. Our code $\texttt{pyilc}$ performs needlet or harmonic ILC on mm-wave sky maps in a flexible manner, with options to deproject various components on all or some scales. We validate our maps and compare them to the official $\textit{Planck}$ 2015 $y$-map, finding that we obtain consistent results on large scales and 10-20$\%$ lower noise on small scales. We expect that these maps will be useful for many auto- and cross-correlation analyses; in a companion paper, we use them to measure the tSZ -- CMB lensing cross-correlation. We anticipate that $\texttt{pyilc}$ will be useful both for data analysis and for pipeline validation on simulations to understand the propagation of foreground components through a full NILC pipeline., Comment: 28 pages, 15 figures. Public maps and other data products are available at https://users.flatironinstitute.org/~fmccarthy/ymaps_PR4_McCH23/ ; public code is available at https://github.com/jcolinhill/pyilc . V2: some additional notes about the effective CIB SED parameters and the halo model, also removed discussion about apodization of our maps and Planck maps. V3: updated references

We present a measurement of the cross-correlation between the MagLim galaxies selected from the Dark Energy Survey (DES) first three years of observations (Y3) and cosmic microwave background (CMB) lensing from the Atacama Cosmology Telescope (ACT) Data Release 4 (DR4), reconstructed over $\sim 436$ sq.deg. of the sky. Our galaxy sample, which covers $\sim 4143$ sq.deg., is divided into six redshift bins spanning the redshift range of $0.20

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.13$ 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: 32 pages, 17 figures, replaced with version accepted in ApJ (Feb 2024). Cosmological likelihood data and mass maps are public 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

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 redshifts, Comment: 45+22 pages, 50 figures. v2 matches with published version in ApJ. Cosmological likelihood data and lensing maps are 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 Madhavacheril et al and MacCrann et al

We investigate the impact and mitigation of extragalactic foregrounds for the CMB lensing power spectrum analysis of Atacama Cosmology Telescope (ACT) data release 6 (DR6) data. Two independent microwave sky simulations are used to test a range of mitigation strategies. We demonstrate that finding and then subtracting point sources, finding and then subtracting models of clusters, and using a profile bias-hardened lensing estimator, together reduce the fractional biases to well below statistical uncertainties, with the inferred lensing amplitude, $A_{\mathrm{lens}}$, biased by less than $0.2\sigma$. We also show that another method where a model for the cosmic infrared background (CIB) contribution is deprojected and high frequency data from Planck is included has similar performance. Other frequency-cleaned options do not perform as well, incurring either a large noise cost, or resulting in biased recovery of the lensing spectrum. In addition to these simulation-based tests, we also present null tests performed on the ACT DR6 data which test for sensitivity of our lensing spectrum estimation to differences in foreground levels between the two ACT frequencies used, while nulling the CMB lensing signal. These tests pass whether the nulling is performed at the map or bandpower level. The CIB-deprojected measurement performed on the DR6 data is consistent with our baseline measurement, implying contamination from the CIB is unlikely to significantly bias the DR6 lensing spectrum. This collection of tests gives confidence that the ACT DR6 lensing measurements and cosmological constraints presented in companion papers to this work are robust to extragalactic foregrounds., Comment: Companion paper to Qu et al and Madhavacheril et al

A key assumption of the standard cosmological model is that the temperature of the cosmic microwave background (CMB) radiation scales with cosmological redshift $z$ as $T_{\rm CMB}(z) \propto (1+z)$ at all times after recombination at $z_\star \simeq 1090$. However, this assumption has only been precisely tested at $z \lesssim 3$. Here, we consider cosmological models with post-recombination reheating (PRR), in which the CMB monopole temperature abruptly increases due to energy injection after last scattering. Such a scenario can potentially resolve tensions between inferences of the current cosmic expansion rate (the Hubble constant, $H_0$). We consider an explicit model in which a metastable sub-component of dark matter (DM) decays to Standard Model photons, whose spectral energy distribution is assumed to be close to that of the CMB blackbody. A fit to Planck CMB anisotropy, COBE/FIRAS CMB monopole, and SH0ES distance-ladder measurements yields $H_0 = 71.2 \pm 1.1$ km/s/Mpc, matter fluctuation amplitude $S_8 = 0.774 \pm 0.018$, and CMB temperature increase $\delta T_{\rm CMB} = 0.109^{+0.033}_{-0.044}$ K, which is sourced by DM decay at $z \gtrsim 10$. However, matter density constraints from baryon acoustic oscillation and supernovae data highly constrain this scenario, with a joint fit to all datasets yielding $H_0 = 68.69 \pm 0.35$ km/s/Mpc, $S_8 = 0.8035 \pm 0.0081$, and $\delta T_{\rm CMB} < 0.0342$ K (95% CL upper limit). These bounds can be weakened if additional dark relativistic species are present in the early universe, yielding higher $H_0$. We conclude that current data disfavor models with significant PRR solely through its impact on background and linear-theory observables, completely independent of CMB spectral distortion constraints. However, a small amount of such energy injection could play a role in restoring cosmological concordance., Comment: 6+8 pages, 1+1 figures

The Primordial Inflation Explorer (PIXIE) is an Explorer-class mission concept to measure the spectrum and polarization of the cosmic microwave background. Cosmological signals are small compared to the instantaneous instrument noise, requiring strict control of instrumental signals. The instrument design provides multiple levels of null operation, signal modulation, and signal differences, with only few-percent systematic error suppression required at each level. Jackknife tests based on discrete instrument symmetries provide an independent means to identify, model, and remove remaining instrumental signals. We use detailed time-ordered simulations, including realistic performance and tolerance parameters, to evaluate the instrument response to broad classes of systematic errors for both spectral distortions and polarization. The largest systematic errors contribute additional white noise at the few-percent level compared to the dominant photon noise. Coherent instrumental effects which do not integrate down are smaller still, and remain several orders of magnitude below the targeted cosmological signals., Comment: 37 pages including 23 figures

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: 23+29 pages, 6+18 figures; matches the version accepted for publication in PRD; code is available at https://github.com/olakusiak/deCIBing

We use the emulation framework CosmoPower to construct and publicly release neural network emulators of cosmological observables, including the Cosmic Microwave Background (CMB) temperature and polarization power spectra, matter power spectrum, distance-redshift relation, baryon acoustic oscillation (BAO) and redshift-space distortion (RSD) observables, and derived parameters. We train our emulators on Einstein-Boltzmann calculations obtained with high-precision numerical convergence settings, for a wide range of cosmological models including $\Lambda$CDM, $w$CDM, $\Lambda$CDM+$N_\mathrm{eff}$, and $\Lambda$CDM+$\Sigma m_\nu$. Our CMB emulators are accurate to better than 0.5% out to $\ell=10^4$ which is sufficient for Stage-IV data analysis, and our $P(k)$ emulators reach the same accuracy level out to $k=50 \,\, \mathrm{Mpc}^{-1}$, which is sufficient for Stage-III data analysis. We release the emulators via an online repository CosmoPower Organisation, which will be continually updated with additional extended cosmological models. Our emulators accelerate cosmological data analysis by orders of magnitude, enabling cosmological parameter extraction analyses, using current survey data, to be performed on a laptop. We validate our emulators by comparing them to CLASS and CAMB and by reproducing cosmological parameter constraints derived from Planck TT, TE, EE, and CMB lensing data, as well as from the Atacama Cosmology Telescope Data Release 4 CMB data, Dark Energy Survey Year-1 galaxy lensing and clustering data, and Baryon Oscillation Spectroscopic Survey Data Release 12 BAO and RSD data.

Diffusion damping of the cosmic microwave background (CMB) power spectrum results from imperfect photon-baryon coupling in the pre-recombination plasma. At redshift $5 \times 10^4 < z < 2 \times 10^6$, the plasma acquires an effective chemical potential, and energy injections from acoustic damping in this era create $\mu$-type spectral distortions of the CMB. These $\mu$ distortions trace the underlying photon density fluctuations, probing the primordial power spectrum in short-wavelength modes $k_\mathrm{S}$ over the range $50 \ \mathrm{Mpc}^{-1} \lesssim k \lesssim 10^4 \ \mathrm{Mpc}^{-1}$. Small-scale power modulated by long-wavelength modes $k_\mathrm{L}$ from squeezed-limit non-Gaussianities introduces cross-correlations between CMB temperature anisotropies and $\mu$ distortions. Under single-field inflation models, $\mu \times T$ correlations measured from an observer in an inertial frame should vanish up to a factor of $(k_\mathrm{L}/k_\mathrm{S})^2 \ll 1$. Thus, any measurable correlation rules out single-field inflation models. We forecast how well the next-generation ground-based CMB experiment CMB-S4 will be able to constrain primordial squeezed-limit non-Gaussianity, parameterized by $f_\mathrm{NL}$, using measurements of $C_{\ell}^{\mu T}$ as well as $C_{\ell}^{\mu E}$ from CMB $E$ modes. Using current experimental specifications and foreground modeling, we expect $\sigma(f_\mathrm{NL}) \lesssim 1000$. This is roughly four times better than the current limit on $f_\mathrm{NL}$ using $\mu \times T$ and $\mu \times E$ correlations from Planck and is comparable to what is achievable with LiteBIRD, demonstrating the power of the CMB-S4 experiment. This measurement is at an effective scale of $k \simeq 740 \ \text{Mpc}^{-1}$ and is thus highly complementary to measurements at larger scales from primary CMB and large-scale structure., Comment: 19 pages, 6 figures

Current cosmological data exhibit discordance between indirect and some direct inferences of the present-day expansion rate, $H_0$. Early dark energy (EDE), which briefly increases the cosmic expansion rate prior to recombination, is a leading scenario for resolving this "Hubble tension" while preserving a good fit to cosmic microwave background (CMB) data. However, this comes at the cost of changes in parameters that affect structure formation in the late-time universe, including the spectral index of scalar perturbations, $n_s$. Here, we present the first constraints on axion-like EDE using data from the Lyman-$\alpha$ forest, i.e., absorption lines imprinted in background quasar spectra by neutral hydrogen gas along the line of sight. We consider two independent measurements of the one-dimensional Ly$\alpha$ forest flux power spectrum, from the Sloan Digital Sky Survey (SDSS eBOSS) and from the MIKE/HIRES and X-Shooter spectrographs. We combine these with a baseline dataset comprised of Planck CMB data and baryon acoustic oscillation (BAO) measurements. Combining the eBOSS Ly$\alpha$ data with the CMB and BAO dataset reduces the 95% confidence level (CL) upper bound on the maximum fractional contribution of EDE to the cosmic energy budget, $f_{\rm EDE}$, from 0.07 to 0.03 and constrains $H_0=67.9_{-0.4}^{+0.4}$ km/s/Mpc (68% CL), with maximum a posteriori value $H_0=67.9$ km/s/Mpc. Similar results are obtained for the MIKE/HIRES and X-Shooter Ly$\alpha$ data. Our Ly$\alpha$-based EDE constraints yield $H_0$ values that are in $>4\sigma$ tension with the SH0ES distance-ladder measurement and are driven by the preference of the Ly$\alpha$ forest data for $n_s$ values lower than those required by EDE cosmologies that fit Planck CMB data. Taken at face value, the Ly$\alpha$ forest severely constrains canonical EDE models that could resolve the Hubble tension., Comment: 7+9 pages, 2+9 figures, accepted by Phys. Rev. Lett

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., Comment: 14+9 pages, 6+6 figures; matches the version accepted for publication in PRD (https://journals.aps.org/prd/abstract/10.1103/PhysRevD.107.083521); code available at at https://github.com/kmsurrao/reMASTERed

Early dark energy (EDE), whose cosmological role is localized in time around the epoch of matter-radiation equality in order to resolve the Hubble tension, introduces a new coincidence problem: why should the EDE dynamics occur near equality if EDE is decoupled from both matter and radiation? The resolution of this problem may lie in an {\it early dark sector} (EDS), wherein the dark matter mass is dependent on the EDE scalar field. Concretely, we consider a Planck-suppressed coupling of EDE to dark matter, as would naturally arise from breaking of the global $U(1)$ shift symmetry of the former by quantum gravity effects. With a sufficiently flat potential, the rise to dominance of dark matter at matter-radiation equality itself triggers the rolling and subsequent decay of the EDE. We show that this {\it trigger} EDS (tEDS) model can naturally resolve the EDE coincidence problem at the background level without any fine tuning of the coupling to dark matter or of the initial conditions. When fitting to current cosmological data, including that from the local distance ladder and the low-redshift amplitude of fluctuations, the tEDS maximum-likelihood model performs comparably to EDE for resolving the Hubble tension, achieving $H_0 =71.2$ km/s/Mpc. However, fitting the \emph{Planck} cosmic microwave background data requires a specific range of initial field positions to balance the scalar field fluctuations that drive acoustic oscillations, providing testable differences with other EDE models., Comment: 12 pages, 6 figures

PICO is a concept for a NASA probe-scale mission aiming to detect or constrain the tensor to scalar ratio $r$, a parameter that quantifies the amplitude of inflationary gravity waves. We carry out map-based component separation on simulations with five foreground models and input $r$ values $r_{in}=0$ and $r_{in} = 0.003$. We forecast $r$ determinations using a Gaussian likelihood assuming either no delensing or a residual lensing factor $A_{\rm lens}$ = 27%. By implementing the first full-sky, post component-separation, map-domain delensing, we show that PICO should be able to achieve $A_{\rm lens}$ = 22% - 24%. For four of the five foreground models we find that PICO would be able to set the constraints $r < 1.3 \times 10^{-4} \,\, \mbox{to} \,\, r <2.7 \times 10^{-4}\, (95\%)$ if $r_{in}=0$, the strongest constraints of any foreseeable instrument. For these models, $r=0.003$ is recovered with confidence levels between $18\sigma$ and $27\sigma$. We find weaker, and in some cases significantly biased, upper limits when removing few low or high frequency bands. The fifth model gives a $3\sigma$ detection when $r_{in}=0$ and a $3\sigma$ bias with $r_{in} = 0.003$. However, by correlating $r$ determinations from many small 2.5% sky areas with the mission's 555 GHz data we identify and mitigate the bias. This analysis underscores the importance of large sky coverage. We show that when only low multipoles $\ell \leq 12$ are used, the non-Gaussian shape of the true likelihood gives uncertainties that are on average 30% larger than a Gaussian approximation., Comment: 34 pages, 13 figures, published in JCAP

Tensions between cosmological parameters (in particular the local expansion rate $H_0$ and the amplitude of matter clustering $S_8$) inferred from low-redshift data and data from the cosmic microwave background (CMB) and large-scale structure (LSS) experiments have inspired many extensions to the standard cosmological model, $\Lambda$CDM. Models which simultaneously lessen both tensions are of particular interest. We consider one scenario with the potential for such a resolution, in which some fraction of the dark matter has converted into dark radiation since the release of the CMB. Such a scenario encompasses and generalizes the more standard "decaying dark matter" model, allowing additional flexibility in the rate and time at which the dark matter converts into dark radiation. In this paper, we constrain this scenario with a focus on exploring whether it can solve (or reduce) these tensions. We find that such a model is effectively ruled out by low-$\ell$ CMB data, in particular by the reduced peak-smearing due to CMB lensing and the excess Integrated Sachs--Wolfe (ISW) signal caused by the additional dark energy density required to preserve flatness after dark matter conversion into dark radiation. Thus, such a model does not have the power to reduce these tensions without further modifications. This conclusion extends and generalizes related conclusions derived for the standard decaying dark matter model., Comment: 19 pages, 8 figures. Our modified Boltzmann code is available at https://github.com/fmccarthy/class_DMDR Re-uploaded with some more details on the failure of the DMDR model (v2); Reuploaded with the version published by PRD (v3)

We show how consistency relations can be used to robustly extract the amplitude of local primordial non-Gaussianity ($f_{\rm NL}$) from the squeezed limit of the matter bispectrum, well into the non-linear regime. First, we derive a non-perturbative relation between primordial non-Gaussianity and the leading term in the squeezed bispectrum, revising some results present in the literature. This relation is then used to successfully measure $f_{\rm NL}$ from $N$-body simulations. We discuss the dependence of our results on different scale cuts and redshifts. Specifically, the analysis is strongly dependent on the choice of the smallest soft momentum, $q_{\rm min}$, which is the most sensitive to primordial bispectrum contributions, but is largely independent of the choice of the largest hard momentum, $k_{\rm max}$, due to the non-Gaussian nature of the covariance. We also show how the constraints on $f_{\rm NL}$ improve at higher redshift, due to a reduced off-diagonal covariance. In particular, for a simulation with $f_{\rm NL} = 100$ and a volume of $(2.4 \text{ Gpc}/h)^3$, we measure $f_{\rm NL} = 98 \pm 12$ at redshift $z=0$ and $f_{\rm NL} = 97 \pm 8$ at $z=0.97$. Finally, we compare our results with a Fisher forecast, showing that the current version of the analysis is satisfactorily close to the Fisher error. We regard this as a first step towards the realistic application of consistency relations to constrain primordial non-Gaussianity using observations., Comment: 17 pages, 8 figures. Minor changes. Matches version published in PRD

Feedback from active galactic nuclei (AGN) and supernovae can affect measurements of integrated SZ flux of halos ($Y_\mathrm{SZ}$) from CMB surveys, and cause its relation with the halo mass ($Y_\mathrm{SZ}-M$) to deviate from the self-similar power-law prediction of the virial theorem. We perform a comprehensive study of such deviations using CAMELS, a suite of hydrodynamic simulations with extensive variations in feedback prescriptions. We use a combination of two machine learning tools (random forest and symbolic regression) to search for analogues of the $Y-M$ relation which are more robust to feedback processes for low masses ($M\lesssim 10^{14}\, h^{-1} \, M_\odot$); we find that simply replacing $Y\rightarrow Y(1+M_*/M_\mathrm{gas})$ in the relation makes it remarkably self-similar. This could serve as a robust multiwavelength mass proxy for low-mass clusters and galaxy groups. Our methodology can also be generally useful to improve the domain of validity of other astrophysical scaling relations. We also forecast that measurements of the $Y-M$ relation could provide percent-level constraints on certain combinations of feedback parameters and/or rule out a major part of the parameter space of supernova and AGN feedback models used in current state-of-the-art hydrodynamic simulations. Our results can be useful for using upcoming SZ surveys (e.g., SO, CMB-S4) and galaxy surveys (e.g., DESI and Rubin) to constrain the nature of baryonic feedback. Finally, we find that the an alternative relation, $Y-M_*$, provides complementary information on feedback than $Y-M$, Comment: Version appearing in MNRAS. Minor change to Fig.6 and added Fig. A5 compared to the previous version. 7+5 pages. The code and data associated with this paper are available at https://github.com/JayWadekar/ScalingRelations_ML

The kinetic Sunyaev-Zel'dovich (kSZ) effect, i.e., the Doppler boost of cosmic microwave background (CMB) photons caused by their scattering off free electrons in galaxy clusters and groups with non-zero bulk velocity, is a powerful window on baryons in the universe. We present the first halo-model computation of the cross-power spectrum of the ``projected-field'' kSZ signal with large-scale structure (LSS) tracers. We compare and validate our calculations against previous studies, which relied on $N$-body-calibrated effective formulas rather than the halo model. We forecast results for CMB maps from the Atacama Cosmology Telescope (AdvACT), Simons Observatory (SO), and CMB-S4, and LSS survey data from the Dark Energy Survey, the Vera C.~Rubin Observatory (VRO), and \textit{Euclid}. In cross-correlation with galaxy number density, for AdvACT $\times$ \textit{unWISE} we forecast an 18$\sigma$ projected-field kSZ detection using data already in hand. Combining SO CMB maps and \textit{unWISE} galaxy catalogs, we expect a $62\sigma$ detection, yielding precise measurements of the gas density profile radial slopes. Additionally, we forecast first detections of the kSZ -- galaxy weak lensing cross-correlation with AdvACT $\times$ VRO/\textit{Euclid} (at 6$\sigma$) and of the kSZ -- CMB weak lensing cross-correlation with SO (at 16$\sigma$). Finally, $\approx 10-20$\% precision measurements of the shape of the gas density profile should be possible with CMB-S4 kSZ -- CMB lensing cross-correlation without using any external datasets.

We use data from the Atacama Cosmology Telescope (ACT) DR4 to search for the presence of neutrino self-interaction in the cosmic microwave background. Consistent with prior works, the posterior distributions we find are bimodal, with one mode consistent with $\Lambda$CDM and one where neutrinos strongly self-interact. By combining ACT data with large-scale information from WMAP, we find that a delayed onset of neutrino free streaming caused by significantly strong neutrino self-interaction is compatible with these data at the $2-3\sigma$ level. As seen in the past, the preference shifts to $\Lambda$CDM with the inclusion of Planck data. We determine that the preference for strong neutrino self-interaction is largely driven by angular scales corresponding to $700 \lesssim \ell \lesssim 1000$ in the ACT E-mode polarization data. This region is expected to be key to discriminate between neutrino self-interacting modes and will soon be probed with more sensitive data., Comment: added additional references, 9+10 pages, 4+8 figures

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.

This is a solicited whitepaper for the Snowmass 2021 community planning exercise. The paper focuses on measurements and science with the Cosmic Microwave Background (CMB). The CMB is foundational to our understanding of modern physics and continues to be a powerful tool driving our understanding of cosmology and particle physics. In this paper, we outline the broad and unique impact of CMB science for the High Energy Cosmic Frontier in the upcoming decade. We also describe the progression of ground-based CMB experiments, which shows that the community is prepared to develop the key capabilities and facilities needed to achieve these transformative CMB measurements., Comment: contribution to Snowmass 2021

This Snowmass 2021 White Paper describes the Cosmic Microwave Background Stage 4 project CMB-S4, which is designed to cross critical thresholds in our understanding of the origin and evolution of the Universe, from the highest energies at the dawn of time through the growth of structure to the present day. We provide an overview of the science case, the technical design, and project plan., Comment: Contribution to Snowmass 2021. arXiv admin note: substantial text overlap with arXiv:1908.01062, arXiv:1907.04473

Cosmological data in the next decade will be characterized by high-precision, multi-wavelength measurements of thousands of square degrees of the same patches of sky. By performing multi-survey analyses that harness the correlated nature of these datasets, we will gain access to new science, and increase the precision and robustness of science being pursued by each individual survey. However, effective application of such analyses requires a qualitatively new level of investment in cross-survey infrastructure, including simulations, associated modeling, coordination of data sharing, and survey strategy. The scientific gains from this new level of investment are multiplicative, as the benefits can be reaped by even present-day instruments, and can be applied to new instruments as they come online., Comment: Contribution to Snowmass 2021

In this paper we will list a few important goals that need to be addressed in the next decade, also taking into account the current discordances between the different cosmological probes, such as the disagreement in the value of the Hubble constant $H_0$, the $\sigma_8$--$S_8$ tension, and other less statistically significant anomalies. While these discordances can still be in part the result of systematic errors, their persistence after several years of accurate analysis strongly hints at cracks in the standard cosmological scenario and the necessity for new physics or generalisations beyond the standard model. In this paper, we focus on the $5.0\,\sigma$ tension between the {\it Planck} CMB estimate of the Hubble constant $H_0$ and the SH0ES collaboration measurements. After showing the $H_0$ evaluations made from different teams using different methods and geometric calibrations, we list a few interesting new physics models that could alleviate this tension and discuss how the next decade's experiments will be crucial. Moreover, we focus on the tension of the {\it Planck} CMB data with weak lensing measurements and redshift surveys, about the value of the matter energy density $\Omega_m$, and the amplitude or rate of the growth of structure ($\sigma_8,f\sigma_8$). We list a few interesting models proposed for alleviating this tension, and we discuss the importance of trying to fit a full array of data with a single model and not just one parameter at a time. Additionally, we present a wide range of other less discussed anomalies at a statistical significance level lower than the $H_0$--$S_8$ tensions which may also constitute hints towards new physics, and we discuss possible generic theoretical approaches that can collectively explain the non-standard nature of these signals.[Abridged], Comment: Contribution to Snowmass 2021. 224 pages, 27 figures. Accepted for publication in JHEAp