Search

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)

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

We present a high-significance cross-correlation of CMB lensing maps from the Atacama Cosmology Telescope (ACT) Data Release 6 (DR6) with spectroscopically calibrated luminous red galaxies (LRGs) from the Dark Energy Spectroscopic Instrument (DESI). We detect this cross-correlation at a significance of 38$\sigma$; combining our measurement with the Planck Public Release 4 (PR4) lensing map, we detect the cross-correlation at 50$\sigma$. Fitting this jointly with the galaxy auto-correlation power spectrum to break the galaxy bias degeneracy with $\sigma_8$, we perform a tomographic analysis in four LRG redshift bins spanning $0.4 \le z \le 1.0$ to constrain the amplitude of matter density fluctuations through the parameter combination $S_8^\times = \sigma_8 \left(\Omega_m / 0.3\right)^{0.4}$. Prior to unblinding, we confirm with extragalactic simulations that foreground biases are negligible and carry out a comprehensive suite of null and consistency tests. Using a hybrid effective field theory (HEFT) model that allows scales as small as $k_{\rm max}=0.6$ $h/{\rm Mpc}$, we obtain a 3.3% constraint on $S_8^\times = \sigma_8 \left(\Omega_m / 0.3\right)^{0.4} = 0.792^{+0.024}_{-0.028}$ from ACT data, as well as constraints on $S_8^\times(z)$ that probe structure formation over cosmic time. Our result is consistent with the early-universe extrapolation from primary CMB anisotropies measured by Planck PR4 within 1.2$\sigma$. Jointly fitting ACT and Planck lensing cross-correlations we obtain a 2.7% constraint of $S_8^\times = 0.776^{+0.019}_{-0.021}$, which is consistent with the Planck early-universe extrapolation within 2.1$\sigma$, with the lowest redshift bin showing the largest difference in mean. The latter may motivate further CMB lensing tomography analyses at $z<0.6$ to assess the impact of potential systematics or the consistency of the $\Lambda$CDM model over cosmic time., Comment: Prepared for submission to JCAP (47 pages, 13 figures)

We derive new constraints on the $E_G$ statistic as a test of gravity, combining the CMB lensing map estimated from Data Release 6 (DR6) of the Atacama Cosmology Telescope with SDSS BOSS CMASS and LOWZ galaxy data. We develop an analysis pipeline to measure the cross-correlation between CMB lensing maps and galaxy data, following a blinding policy and testing the approach through null and consistency checks. By testing the equivalence of the spatial and temporal gravitational potentials, the $E_G$ statistic can distinguish $\Lambda$CDM from alternative models of gravity. We find $E_G= 0.31^{+0.06}_{-0.05}$ for ACT and CMASS data at 68.28\% confidence level, and $E_G = 0.49^{+0.14}_{-0.11}$ for ACT and LOWZ. Systematic errors are estimated to be 3\% and 4\% respectively. Including CMB lensing information from Planck PR4 results in $E_G = 0.34^{+0.05}_{-0.05}$ with CMASS and $E_G= 0.43^{+0.11}_{-0.09}$ with LOWZ. These are consistent with predictions for the $\Lambda$CDM model that best fits the Planck CMB anisotropy and SDSS BOSS BAO, where $E_G^{\rm GR} (z_{\rm eff} = 0.555) = 0.401\pm 0.005$ for CMB lensing combined with CMASS and $E_G^{\rm GR} (z_{\rm eff} = 0.316) = 0.452\pm0.005$ combined with LOWZ. We also find $E_G$ to be scale independent, with PTE $>5\%$, as predicted by general relativity. The methods developed in this work are also applicable to improved future analyses with upcoming spectroscopic galaxy samples and CMB lensing measurements., Comment: 33 pages, 21 figures, prepared for submission to PRD

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

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

The $E_G$ statistic is a discriminating probe of gravity developed to test the prediction of general relativity (GR) for the relation between gravitational potential and clustering on the largest scales in the observable universe. We present a novel high-precision estimator for the $E_G$ statistic using CMB lensing and galaxy clustering correlations that carefully matches the effective redshifts across the different measurement components to minimize corrections. A suite of detailed tests is performed to characterize the estimator's accuracy, its sensitivity to assumptions and analysis choices and the non-Gaussianity of the estimator's uncertainty is characterized. After finalization of the estimator, it is applied to $\textit{Planck}$ CMB lensing and SDSS CMASS and LOWZ galaxy data. We report the first harmonic space measurement of $E_G$ using the LOWZ sample and CMB lensing and also updated constraints using the final CMASS sample and the latest $\textit{Planck}$ CMB lensing map. We find $E_G^{Planck+CMASS} = 0.36^{+0.06}_{-0.05}$ (68.27%) and $E_G^{\rm \textit{Planck}+LOWZ} = 0.40^{+0.11}_{-0.09} $ (68.27%), with additional subdominant systematic error budget estimates of 2% and 3% respectively. Using $\Omega_{\rm m,0}$ constraints from $\textit{Planck}$ and SDSS BAO observations, $\Lambda$CDM-GR predicts $E_G^{\rm GR} (z = 0.555) = 0.401 \pm 0.005$ and $E_G^{\rm GR} (z = 0.316) = 0.452 \pm 0.005$ at the effective redshifts of the CMASS and LOWZ based measurements. We report the measurement to be in good statistical agreement with the $\Lambda$CDM-GR prediction, and report that the measurement is also consistent with the more general GR prediction of scale-independence for $E_G$. This work provides a carefully constructed and calibrated statistic with which $E_G$ measurements can be confidently and accurately obtained with upcoming survey data., Comment: 36 pages, 20 figures, Accepted for publication in PRD

We use the integrated Sachs-Wolfe (ISW) effect, by now detectable at $\sim 5\sigma$ within the context of $\Lambda{}$CDM cosmologies, to place strong constraints on dynamical dark energy theories. Working within an effective field theory framework for dark energy we find that including ISW constraints from galaxy-CMB cross-correlations significantly strengthens existing large-scale structure constraints, yielding bounds consistent with $\Lambda{}$CDM and approximately reducing the viable parameter space by $\sim 70\%$. This is a direct consequence of ${\cal O}(1)$ changes induced in these cross-correlations by otherwise viable dark energy models, which we discuss in detail. We compute constraints by adapting the $\Lambda{}$CDM ISW likelihood from [1] for dynamical dark energy models using galaxy data from 2MASS, WISE $\times$ SuperCOSMOS, SDSS-DR12, QSOs and NVSS, CMB data from Planck 18, and BAO and RSD large scale structure measurements from BOSS and 6dF. We show constraints both in terms of EFT-inspired $\alpha_i$ and phenomenological $\mu/\Sigma$ parametrisations. Furthermore we discuss the approximations involved and related aspects of bias modelling in detail and highlight what these constraints imply for the underlying dark energy theories., Comment: 11 pages + appendices and references, 11 figures

We present a first measurement of the galaxy-galaxy-CMB lensing bispectrum. The signal is detected at $26\sigma$ and $22\sigma$ significance using two samples from the unWISE galaxy catalog at mean redshifts $\bar{z}=0.6$ and $1.1$ and lensing reconstructions from Planck PR4. We employ a compressed bispectrum estimator based on the cross-correlation between the square of the galaxy overdensity field and CMB lensing reconstructions. We present a series of consistency tests to ensure the cosmological origin of our signal and rule out potential foreground contamination. We compare our results to model predictions from a halo model previously fit to only two-point spectra, finding reasonable agreement when restricting our analysis to large scales. Such measurements of the CMB lensing galaxy bispectrum will have several important cosmological applications, including constraining the uncertain higher-order bias parameters that currently limit lensing cross-correlation analyses., Comment: 10 pages, 4 figures, 1 table; v2: replaced HOD model from arXiv:2203.12583 with updated/improved model from arXiv:2303.08121, conclusions unchanged; to be submitted to PRL

A key challenge in the search for primordial B-modes is the presence of polarized Galactic foregrounds, especially thermal dust emission. Power-spectrum-based analysis methods generally assume the foregrounds to be Gaussian random fields when constructing a likelihood and computing the covariance matrix. In this paper, we investigate how non-Gaussianity in the dust field instead affects CMB and foreground parameter inference in the context of inflationary B-mode searches, capturing this effect via modifications to the dust power-spectrum covariance matrix. For upcoming experiments such as the Simons Observatory, we find no dependence of the tensor-to-scalar ratio uncertainty $\sigma(r)$ on the degree of dust non-Gaussianity or the nature of the dust covariance matrix. We provide an explanation of this result, noting that when frequency decorrelation is negligible, dust in mid-frequency channels is cleaned using high-frequency data in a way that is independent of the spatial statistics of dust. We show that our results hold also for non-zero levels of frequency decorrelation that are compatible with existing data. We find, however, that neglecting the impact of dust non-Gaussianity in the covariance matrix can lead to inaccuracies in goodness-of-fit metrics. Care must thus be taken when using such metrics to test B-mode spectra and models, although we show that any such problems can be mitigated by using only cleaned spectrum combinations when computing goodness-of-fit statistics., Comment: 16 pages, 8 figures, accepted version

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

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 propose a new probe of inflationary gravitational waves (IGWs): the cross-correlation of the lensing of inflationary $B$-mode polarization with a large-scale structure (LSS) tracer, which can also be a cosmic microwave background (CMB) lensing map. This is equivalent to measuring a three-point function of two CMB $B$-modes and an LSS tracer. We forecast expected $1\,\sigma$ constraints on the tensor-to-scalar ratio $r$, albeit with a simplistic foreground treatment, and find constraints of $\sigma_r \simeq 7 \times 10^{-3}$ from the correlation of CMB-S4-Deep $B$-mode lensing and LSST galaxies, $\sigma_r \simeq 5 \times 10^{-3}$ from the correlation of CMB-S4-Deep $B$-mode lensing and CMB-S4-Deep CMB lensing, and $\sigma_r \simeq 10^{-2}$ from the correlation of LiteBIRD $B$-mode lensing and CMB-S4-Wide lensing. Because this probe is inherently non-Gaussian, simple Gaussian foregrounds will not produce any biases to the measurement of $r$. While a detailed investigation of non-Gaussian foreground contamination for different cross-correlations will be essential, this observable has the potential to be a useful probe of IGWs, which, due to different sensitivity to many potential sources of systematic errors, can be complementary to standard methods for constraining $r$., Comment: 6 pages, 1 figure, accepted for publication in Physical Review Letters

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

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

Reconstruction of gravitational lensing effects in the CMB from current and upcoming surveys is still dominated by temperature anisotropies. Extragalactic foregrounds in temperature maps can induce significant biases in the lensing power spectrum obtained with the standard quadratic estimators. Techniques such as masking cannot remove these foregrounds fully, and the residuals can still lead to large biases if unaccounted for. In this paper, we study the "shear-only" estimator, an example of a class of geometric methods that suppress extragalactic foreground contamination while making only minimal assumptions about foreground properties. The shear-only estimator has only been formulated in the flat-sky limit and so is not easily applied to wide surveys. Here, we derive the full-sky version of the shear-only estimator and its generalisation to an $m=2$ multipole estimator that has improved performance for lensing reconstruction on smaller scales. The multipole estimator is generally not separable, and so is expensive to compute. We explore separable approximations based on a singular-value decomposition, which allow efficient evaluation of the estimator with real-space methods. Finally, we apply these estimators to simulations that include extragalactic foregrounds and verify their efficacy in suppressing foreground biases., Comment: v2 matches with published version

CMB lensing maps probe the mass distribution in projection out to high redshifts, but significant sensitivity to low-redshift structure remains. In this paper we discuss a method to remove the low-redshift contributions from CMB lensing mass maps by subtracting suitably scaled galaxy density maps, nulling the low redshift structure with a model-insensitive procedure that is similar to delensing. This results in a high-$z$-only mass map that can provide a probe of structure growth at uniquely high redshifts: if systematics can be controlled, we forecast that CMB-S4 lensing combined with a Rubin-LSST-like galaxy survey can probe the amplitude of structure at redshifts $z>3.75$ ($z>5$) to within $2.3\%$ ($3.3\%$). We then discuss other example applications of such high-$z$ CMB lensing maps. In standard analyses of CMB lensing, assuming the wrong dark energy model (or wrong model parametrization) can lead to biases in neutrino mass constraints. In contrast, we show with forecasts that a high-$z$ mass map constructed from CMB-S4 lensing and LSST galaxies can provide a nearly model-independent neutrino mass constraint, with only negligible sensitivity to the presence of non-standard dark energy models, irrespective of their parametrization., Comment: 19 pages, 11 figures

The early dark energy (EDE) solution to the Hubble tension comes at the cost of an increased clustering amplitude that has been argued to worsen the fit to galaxy clustering data. We explore whether freeing the total neutrino mass $M_{\nu}$, which can suppress small-scale structure growth, improves EDE's fit to galaxy clustering. Using Planck Cosmic Microwave Background and BOSS galaxy clustering data, a Bayesian analysis shows that freeing $M_{\nu}$ does not appreciably increase the inferred EDE fraction $f_{\rm EDE}$: we find the 95% C.L. upper limits $f_{\rm EDE}<0.092$ and $M_{\nu}<0.15\,{\rm eV}$. Similarly, in a frequentist profile likelihood setting (where our results support previous findings that prior volume effects are important), we find that the baseline EDE model (with $M_{\nu}=0.06\,{\rm eV}$) provides the overall best fit. For instance, compared to baseline EDE, a model with $M_\nu=0.24\,{\rm eV}$ maintains the same $H_0$(km/s/Mpc)=(70.08, 70.11, respectively) whilst decreasing $S_8$=(0.837, 0.826) to the $\Lambda$CDM level, but worsening the fit significantly by $\Delta \chi^2=7.5$. For the datasets used, these results are driven not by the clustering amplitude, but by background modifications to the late-time expansion rate due to massive neutrinos, which worsen the fit to measurements of the BAO scale., Comment: 9 pages, 6 figures. v3: included additional discussions and 3 figures which further strengthen our results and their interpretation, added further references, slight format change. Accepted for publication in MNRAS. and associated Jupyter notebook available at https://github.com/alexreevesy/ede_mnu

Many theoretical resolutions to the so-called "Hubble tension" rely on modifying the sound horizon at recombination, $r_s$, and thus the acoustic scale used as a standard ruler in the cosmic microwave background (CMB) and large scale structure (LSS) datasets. As shown in a number of recent works, these observables can also be used to compute $r_s$-independent constraints on $H_0$ by making use of the horizon scale at matter-radiation equality, $k_{\rm eq}$, which has different sensitivity to high redshift physics than $r_s$. As such, $r_s$- and $k_{\rm eq}$-based measurements of $H_0$ (within a $\Lambda$CDM framework) may differ if there is new physics present pre-recombination. In this work, we present the tightest constraints on the latter from current data, finding $H_0=64.8^{+2.2}_{-2.5}$ at 68% CL (in $\mathrm{km}\,\mathrm{s}^{-1}\mathrm{Mpc}^{-1}$ units) from a combination of BOSS galaxy power spectra, Planck CMB lensing, and the newly released Pantheon+ supernova constraints, as well as physical priors on the baryon density, neutrino mass, and spectral index. The BOSS and Planck measurements have different degeneracy directions, leading to the improved combined constraints, with a bound of $H_0 = 67.1^{+2.5}_{-2.9}$ ($63.6^{+2.9}_{-3.6}$) from BOSS (Planck) alone. The results show some dependence on the neutrino mass bounds, with the constraint broadening to $H_0 = 68.0^{+2.9}_{-3.2}$ if we instead impose a weak prior on $\sum m_\nu$ from terrestrial experiments, or shifting to $H_0 = 64.6\pm2.4$ if the neutrino mass is fixed to its minimal value. Even without dependence on the sound horizon, our results are in $\approx 3\sigma$ tension with those obtained from the Cepheid-calibrated distance ladder, which begins to cause problems for new physics models that vary $H_0$ by changing acoustic physics or the expansion history immediately prior to recombination., Comment: 11 pages, 3 figures, accepted by Phys. Rev. D

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

$H_0$ constraints from galaxy surveys are sourced by the geometric properties of two standardisable rulers: the sound horizon scale, $r_s$, and the matter-radiation equality scale, $k_{\rm eq}$. While most analyses over the last decade have focused on the first scale, recent work has emphasised that the second can provide an independent source of information about the expansion rate of the universe. Recent approaches to obtain a sound-horizon-independent measurement of $H_0$ from the equality scale have avoided $r_s$-based information by removing the sound-horizon-calibrating prior on the baryon density. We present a new method to marginalise over $r_s$; this allows baryon information to be retained enabling tighter parameter constraints. For a Euclid-like spectroscopic survey, we forecast sound-horizon-independent $H_0$ constraints of $\sigma_{H_0} = 0.7\rm{\ km\ s^{-1}\ Mpc^{-1}}$ for our method using the equality scale, compared with $\sigma_{H_0} = 0.5\rm{\ km\ s^{-1}\ Mpc^{-1}}$ from the sound horizon. Upcoming equality scale $H_0$ measurements thus can be highly competitive, although we caution that the impact of observational systematics on such measurements still needs to be investigated in detail. Applying our new approach to the BOSS power spectrum gives $H_0 = 69.5^{+3.0}_{-3.5}\rm{\ km\ s^{-1}\ Mpc^{-1}}$ from equality alone, somewhat tighter than previous constraints. Consistency of $r_s$- and $k_{\rm eq}$-based $H_0$ measurements can provide a valuable internal consistency test of the cosmological model; as an example, we consider the change in $H_0$ created by early dark energy. Assuming the \textit{Planck}+SH0ES best-fit EDE model we find a $2.6\sigma$ shift ($\Delta H_0 = 2.6\rm{\ km\ s^{-1}\ Mpc^{-1}}$) between the two measurements for Euclid; if we instead assume the ACT best-fit model, this increases to $9.0\sigma$ ($\Delta H_0 = 7.8\rm{\ km\ s^{-1}\ Mpc^{-1}}$)., Comment: 19 Pages, 9 Figures, 4 Tables, comments welcome, published in PRD

A key challenge for current and upcoming CMB lensing measurements is their sensitivity to biases from extragalactic foregrounds, such as Sunyaev-Zeldovich (SZ) signals or cosmic infrared background emission. Several methods have been developed to mitigate these lensing foreground biases, dividing broadly into multi-frequency cleaning approaches and modifications to the estimator geometry, but how to optimally combine these methods has not yet been explored in detail. In this paper, we examine which combination of lensing foreground mitigation strategies is best able to reduce the impact of foreground contamination for a Simons-Observatory-like experiment while preserving maximal signal-to-noise. Although the optimal combination obtained depends on whether bias- or variance-reduction are prioritized and on whether polarization data is used, generally, we find that combinations involving both geometric (profile hardening, source hardening or shear) and multifrequency (symmetric cleaning) methods perform best. For lensing power spectrum measurements from temperature (polarization and temperature), our combined estimator methods are able to reduce the bias to the lensing amplitude below $\sigma/4$ or 0.3% (0.1%), a factor of 16 (30) lower than the standard QE bias, at a modest signal-to-noise cost of only 18% (12%). In contrast, single-method foreground-mitigation approaches struggle to reduce the bias to a negligible level below $\sigma/2$ without incurring a large noise penalty. For upcoming and current experiments, our combined methods therefore represent a promising approach for making lensing measurements with negligible foreground bias.

We introduce and validate a delensing framework for the Simons Observatory (SO), which will be used to improve constraints on inflationary gravitational waves (IGWs) by reducing the lensing noise in measurements of the $B$-modes in CMB polarization. SO will initially observe CMB by using three small aperture telescopes and one large-aperture telescope. While polarization maps from small-aperture telescopes will be used to constrain IGWs, the internal CMB lensing maps used to delens will be reconstructed from data from the large-aperture telescope. Since lensing maps obtained from the SO data will be noise-dominated on sub-degree scales, the SO lensing framework constructs a template for lensing-induced $B$-modes by combining internal CMB lensing maps with maps of the cosmic infrared background from Planck as well as galaxy density maps from the LSST survey. We construct a likelihood for constraining the tensor-to-scalar ratio $r$ that contains auto- and cross-spectra between observed $B$-modes and the lensing $B$-mode template. We test our delensing analysis pipeline on map-based simulations containing survey non-idealities, but that, for this initial exploration, does not include contamination from Galactic and extragalactic foregrounds. We find that the SO survey masking and inhomogeneous and atmospheric noise have very little impact on the delensing performance, and the $r$ constraint becomes $\sigma(r)\approx 0.0015$ which is close to that obtained from the idealized forecasts in the absence of the Galactic foreground and is nearly a factor of two tighter than without delensing. We also find that uncertainties in the external large-scale structure tracers used in our multi-tracer delensing pipeline lead to bias much smaller than the $1\,\sigma$ statistical uncertainties., Comment: 22 pages, 14 figures

The early dark energy (EDE) scenario aims to increase the value of the Hubble constant ($H_0$) inferred from cosmic microwave background (CMB) data over that found in $\Lambda$CDM, via the introduction of a new form of energy density in the early universe. The EDE component briefly accelerates cosmic expansion just prior to recombination, which reduces the physical size of the sound horizon imprinted in the CMB. Previous work has found that non-zero EDE is not preferred by Planck CMB power spectrum data alone, which yield a 95% confidence level (CL) upper limit $f_{\rm EDE} < 0.087$ on the maximal fractional contribution of the EDE field to the cosmic energy budget. In this paper, we fit the EDE model to CMB data from the Atacama Cosmology Telescope (ACT) Data Release 4. We find that a combination of ACT, large-scale Planck TT (similar to WMAP), Planck CMB lensing, and BAO data prefers the existence of EDE at $>99.7$% CL: $f_{\rm EDE} = 0.091^{+0.020}_{-0.036}$, with $H_0 = 70.9^{+1.0}_{-2.0}$ km/s/Mpc (both 68% CL). From a model-selection standpoint, we find that EDE is favored over $\Lambda$CDM by these data at roughly $3\sigma$ significance. In contrast, a joint analysis of the full Planck and ACT data yields no evidence for EDE, as previously found for Planck alone. We show that the preference for EDE in ACT alone is driven by its TE and EE power spectrum data. The tight constraint on EDE from Planck alone is driven by its high-$\ell$ TT power spectrum data. Understanding whether these differing constraints are physical in nature, due to systematics, or simply a rare statistical fluctuation is of high priority. The best-fit EDE models to ACT and Planck exhibit coherent differences across a wide range of multipoles in TE and EE, indicating that a powerful test of this scenario is anticipated with near-future data from ACT and other ground-based experiments., Comment: v1: 25+20 pages, 6+18 figures, submitted to PRD; v2: 26+20 pages, 7+18 figures, results unchanged, matches PRD accepted version. MCMC chains available at https://lambda.gsfc.nasa.gov/product/act/actpol_mcmc_chains_get.html

The search for cosmic polarization rotation or birefringence in the CMB is well-motivated because it can provide powerful constraints on parity-violating new physics, such as axion-like particles. In this paper we point out that since the CMB polarization is produced at two very different redshifts - it is generated at both reionization and recombination - new parity-violating physics can generically rotate the polarization signals from these different sources by different amounts. We explore two implications of this. First, measurements of CMB birefringence are challenging because the effect is degenerate with a miscalibration of CMB polarization angles; however, by taking the difference of the reionization and recombination birefringence angles (measured from different CMB angular scales), we can obtain a cosmological signal that is immune to instrumental angle miscalibration. Second, we note that the combination with other methods for probing birefringence can give tomographic information, constraining the redshift origin of any physics producing birefringence. We forecast that the difference of the reionization and recombination birefringence angles can be competitively determined to within ~0.05 degrees for future CMB satellites such as LiteBIRD. Although much further work is needed, we argue that foreground mitigation for this measurement should be less challenging than for inflationary B-mode searches on similar scales due to larger signals and lower foregrounds., Comment: 7 pages, 3 figures; comments welcome

We seek to clarify the origin of constraints on the dark energy equation of state parameter from CMB lensing tomography, that is the combination of galaxy clustering and the cross-correlation of galaxies with CMB lensing in a number of redshift bins. In particular, we consider the two-point correlation functions which can be formed with a catalog of galaxy locations and photometric redshifts from the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST) and CMB lensing maps from the CMB-S4 experiment. We focus on the analytic understanding of the origin of the constraints. Dark energy information in these data arises from the influence of three primary relationships: distance as a function of redshift (geometry), the amplitude of the power spectrum as a function of redshift (growth), and the power spectrum as a function of wavenumber (shape). We find that the effects from geometry and growth play a significant role and partially cancel each other out, while the shape effect is unimportant. We also show that Dark Energy Task Force (DETF) Figure of Merit (FoM) forecasts from the combination of LSST galaxies and CMB-S4 lensing are comparable to the forecasts from cosmic shear in the absence of the CMB lensing map, thus providing an important independent check. Compared to the forecasts with the LSST galaxies alone, combining CMB lensing and LSST clustering information (together with the primary CMB spectra) increases the FoM by roughly a factor of 3-4 in the optimistic scenario where systematics are fully under control. We caution that achieving these forecasts will likely require a full analysis of higher-order biasing, photometric redshift uncertainties, and stringent control of other systematic limitations, which are outside the scope of this work, whose primary purpose is to elucidate the physical origin of the constraints., Comment: 9 pages, 6 figures

Two recent large data releases for the Atacama Cosmology Telescope (ACT), called DR4 and DR5, are available for public access. These data include temperature and polarization maps that cover nearly half the sky at arcminute resolution in three frequency bands; lensing maps and component-separated maps covering ~ 2,100 deg^2 of sky; derived power spectra and cosmological likelihoods; a catalog of over 4,000 galaxy clusters; and supporting ancillary products including beam functions and masks. The data and products are described in a suite of ACT papers; here we provide a summary. In order to facilitate ease of access to these data we present a set of Jupyter IPython notebooks developed to introduce users to DR4, DR5, and the tools needed to analyze these data. The data products (excluding simulations) and the set of notebooks are publicly available on the NASA Legacy Archive for Microwave Background Data Analysis (LAMBDA); simulation products are available on the National Energy Research Scientific Computing Center (NERSC)., Comment: Accepted to ApJS. 21 pages, 8 figures. Data and notebooks available on LAMBDA https://lambda.gsfc.nasa.gov/product/act/

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

We propose a new reionization probe that uses cosmic microwave background (CMB) observations; the cross-correlation between fluctuations in the CMB optical depth which probes the integrated electron density, $\delta\tau$, and the Compton $y$-map which probes the integrated electron pressure. This cross-correlation is much less contaminated than the $y$-map power spectrum by late-time cluster contributions. In addition, this cross-correlation can constrain the temperature of ionized bubbles while the optical-depth fluctuations and kinetic SZ effect can not. We measure this new observable using a Planck $y$-map as well as a map of optical-depth fluctuations that we reconstruct from Planck CMB temperature data. We use our measurements to derive a first CMB-only upper limit on the temperature inside ionized bubbles, $T_{\rm b}\lesssim 7.0\times10^5\,$K ($2\,\sigma$). We also present future forecasts, assuming a fiducial model with characteristic reionization bubble size $R_{\rm b}=5\,$Mpc and $T_{\rm b}=5\times10^4\,$K. The signal-to-noise ratio of the fiducial cross-correlation using a signal dominated PICO-like $y$-map becomes $\simeq7$ with CMB-S4 $\delta\tau$ and $\simeq13$ with CMB-HD $\delta\tau$. For the fiducial model, we predict that the CMB-HD $-$ PICO cross-correlation should achieve an accurate measurement of the reionization parameters; $T_{\rm b}\simeq 49800^{+4500}_{-5100}\,$K and $R_{\rm b}\simeq 5.09^{+0.66}_{-0.79}\,$Mpc. Since the power spectrum of the electron density fluctuations is constrained by the $\delta\tau$ auto spectrum, the temperature constraints should be only weakly model-dependent on the details of the electron distributions and should be statistically representative of the temperature in ionized bubbles during reionization. This cross-correlation could, therefore, become an important observable for future CMB experiments., Comment: 13 pages, 9 figures, Accepted for publication in PRD

Internal Linear Combination (ILC) methods are some of the most widely used multi-frequency cleaning techniques employed in CMB data analysis. These methods reduce foregrounds by minimizing the total variance in the coadded map (subject to a signal-preservation constraint), although often significant foreground residuals or biases remain. A modification to the ILC method is the constrained ILC (cILC), which explicitly nulls certain foreground components; however, this foreground nulling often comes at a high price for ground-based CMB datasets, with the map noise increasing significantly on small scales. In this paper we explore a new method, the partially constrained ILC (pcILC), which allows us to optimize the tradeoff between foreground bias and variance in ILC methods. In particular, this method allows us to minimize the variance subject to an inequality constraint requiring that the constrained foregrounds are reduced by at least a fixed factor, which can be chosen based on the foreground sensitivity of the intended application. We test our method on simulated sky maps for a Simons Observatory-like experiment; we find that for cleaning thermal Sunyaev-Zel'dovich (tSZ) contamination at $\ell \in [3000,4800]$, if a small tSZ residual of 20% of the standard ILC residual can be tolerated, the variance of the CMB temperature map is reduced by at least 50% over the cILC value. We also demonstrate an application of this method to reduce noise in CMB lensing reconstruction.

We measure the cross-correlation between galaxy weak lensing data from the Kilo Degree Survey (KiDS-1000, DR4) and cosmic microwave background (CMB) lensing data from the Atacama Cosmology Telescope (ACT, DR4) and the Planck Legacy survey. We use two samples of source galaxies, selected with photometric redshifts, $(0.1

The power spectrum of cosmic microwave background (CMB) lensing will be measured to sub-percent precision with upcoming surveys, enabling tight constraints on the sum of neutrino masses and other cosmological parameters. Measuring the lensing power spectrum involves the estimation of the connected trispectrum of the four-point function of the CMB map, which requires the subtraction of a large Gaussian disconnected noise bias. This reconstruction noise bias receives contributions both from CMB and foreground fluctuations as well as instrument noise (both detector and atmospheric noise for ground-based surveys). The debiasing procedure therefore relies on the quality of simulations of the instrument noise which may be expensive or inaccurate. We propose a new estimator that makes use of at least four splits of the CMB maps with independent instrument noise. This estimator makes the CMB lensing power spectrum completely insensitive to any assumptions made in modeling or simulating the instrument noise. We show that this estimator, in many practical situations, leads to no substantial loss in signal-to-noise. We provide an efficient algorithm for its computation that scales with the number of splits $m$ as $\mathcal{O}(m^2)$ as opposed to a naive $\mathcal{O}(m^4)$ expectation., Comment: 18 pages, 4 figures, minor changes in v2 to match published version

We use gravitational lensing of the cosmic microwave background (CMB) to measure the mass of the most distant blindly-selected sample of galaxy clusters on which a lensing measurement has been performed to date. In CMB data from the the Atacama Cosmology Telescope (ACT) and the Planck satellite, we detect the stacked lensing effect from 677 near-infrared-selected galaxy clusters from the Massive and Distant Clusters of WISE Survey (MaDCoWS), which have a mean redshift of $ \langle z \rangle = 1.08$. There are no current optical weak lensing measurements of clusters that match the distance and average mass of this sample. We detect the lensing signal with a significance of $4.2 \sigma$. We model the signal with a halo model framework to find the mean mass of the population from which these clusters are drawn. Assuming that the clusters follow Navarro-Frenk-White density profiles, we infer a mean mass of $\langle M_{500c}\rangle = \left(1.7 \pm 0.4 \right)\times10^{14}\,\mathrm{M}_\odot$. We consider systematic uncertainties from cluster redshift errors, centering errors, and the shape of the NFW profile. These are all smaller than 30% of our reported uncertainty. This work highlights the potential of CMB lensing to enable cosmological constraints from the abundance of distant clusters populating ever larger volumes of the observable Universe, beyond the capabilities of optical weak lensing measurements., Comment: 14 pages, 3 figures, matches version accepted in ApJL, code available at https://github.com/ACTCollaboration/madcows_lensing/

CMB-S4---the next-generation ground-based cosmic microwave background (CMB) experiment---is set to significantly advance the sensitivity of CMB measurements and enhance 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. Among the science cases pursued with CMB-S4, the quest for detecting primordial gravitational waves is a central driver of the experimental design. This work details the development of a forecasting framework that includes a power-spectrum-based semi-analytic projection tool, targeted explicitly towards optimizing constraints on the tensor-to-scalar ratio, $r$, in the presence of Galactic foregrounds and gravitational lensing of the CMB. This framework is unique in its direct use of information from the achieved performance of current Stage 2--3 CMB experiments to robustly forecast the science reach of upcoming CMB-polarization endeavors. The methodology allows for rapid iteration over experimental configurations and offers a flexible way to optimize the design of future experiments given a desired scientific goal. To form a closed-loop process, we couple this semi-analytic tool with map-based validation studies, which allow for the injection of additional complexity and verification of our forecasts with several independent analysis methods. We document multiple rounds of forecasts for CMB-S4 using this process and the resulting establishment of the current reference design of the primordial gravitational-wave component of the Stage-4 experiment, optimized to achieve our science goals of detecting primordial gravitational waves for $r > 0.003$ at greater than $5\sigma$, or, in the absence of a detection, of reaching an upper limit of $r < 0.001$ at $95\%$ CL., Comment: 24 pages, 8 figures, 9 tables, submitted to ApJ. arXiv admin note: text overlap with arXiv:1907.04473

Two sources of geometric information are encoded in the galaxy power spectrum: the sound horizon at recombination and the horizon at matter-radiation equality. Analyzing the BOSS DR12 galaxy power spectra using perturbation theory with $\Omega_m$ priors from Pantheon supernovae but no priors on $\Omega_b$, we obtain constraints on $H_0$ from the second scale, finding $H_0 = 65.1^{+3.0}_{-5.4}\,\mathrm{km}\,\mathrm{s}^{-1}\mathrm{Mpc}^{-1}$; this differs from the best-fit of SH0ES at 95\% confidence. Similar results are obtained if $\Omega_m$ is constrained from uncalibrated BAO: $H_0 = 65.6^{+3.4}_{-5.5}\,\mathrm{km}\,\mathrm{s}^{-1}\mathrm{Mpc}^{-1}$. Adding the analogous lensing results from Baxter \& Sherwin 2020, the posterior shifts to $70.6^{+3.7}_{-5.0}\,\mathrm{km}\,\mathrm{s}^{-1}\mathrm{Mpc}^{-1}$. Using mock data, Fisher analyses, and scale-cuts, we demonstrate that our constraints do not receive significant information from the sound horizon scale. Since many models resolve the $H_0$ controversy by adding new physics to alter the sound horizon, our measurements are a consistency test for standard cosmology before recombination. A simple forecast indicates that such constraints could reach $\sigma_{H_0} \simeq 1.6\,\mathrm{km}\,\mathrm{s}^{-1}\mathrm{Mpc}^{-1}$ in the era of Euclid., Comment: 8 pages, 4 figures, accepted by Phys. Rev. D

We present LCDM cosmological parameter constraints obtained from delensed microwave background power spectra. Lensing maps from a subset of DR4 data from the Atacama Cosmology Telescope (ACT) are used to undo the lensing effect in ACT spectra observed at 150 and 98 GHz. At 150 GHz, we remove the lensing distortion with an effective efficiency of 30% (TT), 30% (EE), 26% (TE) and 20% (BB); this results in detections of the delensing effect at 8.7 sigma (TT), 5.1 sigma (EE), 2.6 sigma (TE), and 2.4 sigma (BB) significance. The combination of 150 and 98 GHz TT, EE, and TE delensed spectra is well fit by a standard LCDM model. We also measure the shift in best-fit parameters when fitting delensed versus lensed spectra; while this shift does not inform our ability to measure cosmological parameters, it does provide a three-way consistency check among the lensing inferred from the best-fit parameters, the lensing in the CMB power spectrum, and the reconstructed lensing map. This shift is predicted to be zero when fitting with the correct model since both lensed and delensed spectra originate from the same region of sky. Fitting with a LCDM model and marginalizing over foregrounds, we find that the shift in cosmological parameters is consistent with zero. Our results show that gravitational lensing of the microwave background is internally consistent within the framework of the standard cosmological model., Comment: 29 pages, 17 figures, version matches that accepted by JCAP

Large-scale Fourier modes of the cosmic density field are of great value for learning about cosmology because of their well-understood relationship to fluctuations in the early universe. However, cosmic variance generally limits the statistical precision that can be achieved when constraining model parameters using these modes as measured in galaxy surveys, and moreover, these modes are sometimes inaccessible due to observational systematics or foregrounds. For some applications, both limitations can be circumvented by reconstructing large-scale modes using the correlations they induce between smaller-scale modes of an observed tracer (such as galaxy positions). In this paper, we further develop a formalism for this reconstruction, using a quadratic estimator similar to the one used for lensing of the cosmic microwave background. We incorporate nonlinearities from gravity, nonlinear biasing, and local-type primordial non-Gaussianity, and verify that the estimator gives the expected results when applied to N-body simulations. We then carry out forecasts for several upcoming surveys, demonstrating that, when reconstructed modes are included alongside directly-observed tracer density modes, constraints on local primordial non-Gaussianity are generically tightened by tens of percents compared to standard single-tracer analyses. In certain cases, these improvements arise from cosmic variance cancellation, with reconstructed modes taking the place of modes of a separate tracer, thus enabling an effective "multitracer" approach with single-tracer observations., Comment: 30 pages plus 14 pages appendices, 19 figures

This paper presents a maximum-likelihood algorithm for combining sky maps with disparate sky coverage, angular resolution and spatially varying anisotropic noise into a single map of the sky. We use this to merge hundreds of individual maps covering the 2008-2018 ACT observing seasons, resulting in by far the deepest ACT maps released so far. We also combine the maps with the full Planck maps, resulting in maps that have the best features of both Planck and ACT: Planck's nearly white noise on intermediate and large angular scales and ACT's high-resolution and sensitivity on small angular scales. The maps cover over 18,000 square degrees, nearly half the full sky, at 100, 150 and 220 GHz. They reveal 4,000 optically-confirmed clusters through the Sunyaev Zel'dovich effect (SZ) and 18,500 point source candidates at $> 5\sigma$, the largest single collection of SZ clusters and millimeter wave sources to date. The multi-frequency maps provide millimeter images of nearby galaxies and individual Milky Way nebulae, and even clear detections of several nearby stars. Other anticipated uses of these maps include, for example, thermal SZ and kinematic SZ cluster stacking, CMB cluster lensing and galactic dust science. The method itself has negligible bias. However, due to the preliminary nature of some of the component data sets, we caution that these maps should not be used for precision cosmological analysis. The maps are part of ACT DR5, and are available on LAMBDA at https://lambda.gsfc.nasa.gov/product/act/actpol_prod_table.cfm. There is also a web atlas at https://phy-act1.princeton.edu/public/snaess/actpol/dr5/atlas., Comment: 38 pages, 29 figures, data release on lambda. Published in JCAP

We present new arcminute-resolution maps of the Cosmic Microwave Background temperature and polarization anisotropy from the Atacama Cosmology Telescope, using data taken from 2013-2016 at 98 and 150 GHz. The maps cover more than 17,000 deg$^2$, the deepest 600 deg$^2$ with noise levels below 10 $\mu$K-arcmin. We use the power spectrum derived from almost 6,000 deg$^2$ of these maps to constrain cosmology. The ACT data enable a measurement of the angular scale of features in both the divergence-like polarization and the temperature anisotropy, tracing both the velocity and density at last-scattering. From these one can derive the distance to the last-scattering surface and thus infer the local expansion rate, $H_0$. By combining ACT data with large-scale information from WMAP we measure $H_0 = 67.6 \pm 1.1$ km/s/Mpc, at 68% confidence, in excellent agreement with the independently-measured Planck satellite estimate (from ACT alone we find $H_0 = 67.9 \pm 1.5$ km/s/Mpc). The $\Lambda$CDM model provides a good fit to the ACT data, and we find no evidence for deviations: both the spatial curvature, and the departure from the standard lensing signal in the spectrum, are zero to within 1$\sigma$; the number of relativistic species, the primordial Helium fraction, and the running of the spectral index are consistent with $\Lambda$CDM predictions to within $1.5 - 2.2\sigma$. We compare ACT, WMAP, and Planck at the parameter level and find good consistency; we investigate how the constraints on the correlated spectral index and baryon density parameters readjust when adding CMB large-scale information that ACT does not measure. The DR4 products presented here will be publicly released on the NASA Legacy Archive for Microwave Background Data Analysis., Comment: 33 pages, 24 figures, products available on the NASA LAMBDA website, version accepted for publication in JCAP

We present the temperature and polarization angular power spectra of the CMB measured by the Atacama Cosmology Telescope (ACT) from 5400 deg$^2$ of the 2013-2016 survey, which covers $>$15000 deg$^2$ at 98 and 150 GHz. For this analysis we adopt a blinding strategy to help avoid confirmation bias and, related to this, show numerous checks for systematic error done before unblinding. Using the likelihood for the cosmological analysis we constrain secondary sources of anisotropy and foreground emission, and derive a "CMB-only" spectrum that extends to $\ell=4000$. At large angular scales, foreground emission at 150 GHz is $\sim$1% of TT and EE within our selected regions and consistent with that found by Planck. Using the same likelihood, we obtain the cosmological parameters for $\Lambda$CDM for the ACT data alone with a prior on the optical depth of $\tau=0.065\pm0.015$. $\Lambda$CDM is a good fit. The best-fit model has a reduced $\chi^2$ of 1.07 (PTE=0.07) with $H_0=67.9\pm1.5$ km/s/Mpc. We show that the lensing BB signal is consistent with $\Lambda$CDM and limit the celestial EB polarization angle to $\psi_P =-0.07^{\circ}\pm0.09^{\circ}$. We directly cross correlate ACT with Planck and observe generally good agreement but with some discrepancies in TE. All data on which this analysis is based will be publicly released., Comment: 44 pages, 27 figures, products available on the NASA LAMBDA website, version accepted for publication in JCAP

Measurements of the Hubble constant, $H_0$, from the cosmic distance ladder are currently in tension with the value inferred from Planck observations of the CMB and other high redshift datasets if a flat $\Lambda$CDM cosmological model is assumed. One of the few promising theoretical resolutions of this tension is to invoke new physics that changes the sound horizon scale in the early universe; this can bring CMB and BAO constraints on $H_0$ into better agreement with local measurements. In this paper, we discuss how a measurement of the Hubble constant can be made from the CMB without using information from the sound horizon scale, $r_s$. In particular, we show how measurements of the CMB lensing power spectrum can be used to place interesting constraints on $H_0$ when combined with measurements of either supernovae or galaxy weak lensing, which constrain the matter density parameter. The constraints arise from the sensitivity of the CMB lensing power spectrum to the horizon scale at matter-radiation equality (in projection); this scale could have a different dependence on new physics than the sound horizon. From an analysis of current CMB lensing data from Planck and Pantheon supernovae with conservative external priors, we derive an $r_s$-independent constraint of $H_0 = 73.5\pm 5.3$ km/s/Mpc. Forecasts for future CMB surveys indicate that improving constraints beyond an error of $\sigma(H_0) = 3$ km/s/Mpc will be difficult with CMB lensing, although applying similar methods to the galaxy power spectrum may allow for further improvements., Comment: 14 pages, 6 figures; minor revision to match version accepted by MNRAS

We present new constraints on anisotropic birefringence of the cosmic microwave background polarization using two seasons of data from the Atacama Cosmology Telescope covering $456$ square degrees of sky. The birefringence power spectrum, measured using a curved-sky quadratic estimator, is consistent with zero. Our results provide the tightest current constraint on birefringence over a range of angular scales between $5$ arcminutes and $9$ degrees. We improve previous upper limits on the amplitude of a scale-invariant birefringence power spectrum by a factor of between $2$ and $3$. Assuming a nearly-massless axion field during inflation, our result is equivalent to a $2\,\sigma$ upper limit on the Chern-Simons coupling constant between axions and photons of $g_{\alpha\gamma}<4.0\times 10^{-2}/H_I$ where $H_I$ is the inflationary Hubble scale., Comment: 18 pages, 3 figures, Accepted for publication in PRD

We present a measurement of the gravitational lensing deflection power spectrum reconstructed with two seasons cosmic microwave background polarization data from the POLARBEAR experiment. Observations were taken at 150 GHz from 2012 to 2014 which survey three patches of sky totaling 30 square degrees. We test the consistency of the lensing spectrum with a Cold Dark Matter (CDM) cosmology and reject the no-lensing hypothesis at a confidence of 10.9 sigma including statistical and systematic uncertainties. We observe a value of A_L = 1.33 +/- 0.32 (statistical) +/- 0.02 (systematic) +/- 0.07 (foreground) using all polarization lensing estimators, which corresponds to a 24% accurate measurement of the lensing amplitude. Compared to the analysis of the first year data, we have improved the breadth of both the suite of null tests and the error terms included in the estimation of systematic contamination.

We report a $5.3\sigma$ detection of the gravitational lensing effect of cosmic voids from the Baryon Oscillation Spectroscopic (BOSS) Data Release 12 seen in the $Planck$ 2018 cosmic microwave background (CMB) lensing convergence map. To make this detection, we introduce new optimal techniques for void stacking and filtering of the CMB maps, such as binning voids by a combination of their observed galaxy density and size to separate those with distinctive lensing signatures. We calibrate theoretical expectations for the void-lensing signal using mock catalogs generated in a suite of 108 full-sky lensing simulations from Takahashi et al. (2017). Relative to these templates, we measure the lensing amplitude parameter in the data to be $A_L=1.10 \pm 0.21$ using a matched-filter stacking technique, and confirm it using an alternative Wiener filtering method. We demonstrate that the result is robust against thermal Sunyaev-Zel'dovich contamination and other sources of systematics. We use the lensing measurements to test the relationship between the matter and galaxy distributions within voids, and show that the assumption of linear bias with a value consistent with galaxy clustering results is discrepant with observation at $\sim 3\sigma$; we explain why such a result is consistent with simulations and previous results, and is expected as a consequence of void selection effects. We forecast the potential for void-CMB lensing measurements in future data from the Advanced ACT, Simons Observatory and CMB-S4 experiments, showing that, for the same number of voids, the achievable precision improves by a factor of more than two compared to $Planck$., Comment: 20 pages, 8 figures, 3 tables; addressed minor comments from the reviewer; accepted for publication in ApJ

Optimal analyses of many signals in the cosmic microwave background (CMB) require map-level extraction of individual components in the microwave sky, rather than measurements at the power spectrum level alone. To date, nearly all map-level component separation in CMB analyses has been performed exclusively using satellite data. In this paper, we implement a component separation method based on the internal linear combination (ILC) approach which we have designed to optimally account for the anisotropic noise (in the 2D Fourier domain) often found in ground-based CMB experiments. Using this method, we combine multi-frequency data from the Planck satellite and the Atacama Cosmology Telescope Polarimeter (ACTPol) to construct the first wide-area, arcminute-resolution component-separated maps (covering approximately 2100 sq. deg.) of the CMB temperature anisotropy and the thermal Sunyaev-Zel'dovich (tSZ) effect sourced by the inverse-Compton scattering of CMB photons off hot, ionized gas. Our ILC pipeline allows for explicit deprojection of various contaminating signals, including a modified blackbody approximation of the cosmic infrared background (CIB) spectral energy distribution. The cleaned CMB maps will be a useful resource for CMB lensing reconstruction, kinematic SZ cross-correlations, and primordial non-Gaussianity studies. The tSZ maps will be used to study the pressure profiles of galaxies, groups, and clusters through cross-correlations with halo catalogs, with dust contamination controlled via CIB deprojection. The data products described in this paper are available on LAMBDA., Comment: 24 pages, 11 figures, matches version accepted by PRD. Data products are available on LAMBDA at https://lambda.gsfc.nasa.gov/product/act/act_dr4_derived_maps_get.cfm

Over the past decade, the gravitational lensing of the Cosmic Microwave Background (CMB) has become a powerful tool for probing the matter distribution in the Universe. The standard technique used to reconstruct the CMB lensing signal employs the quadratic estimator (QE) method, which has recently been shown to be suboptimal for lensing measurements on very small scales in temperature and polarization data. We implement a simple, more optimal method for the small-scale regime, which involves taking the direct inverse of the background gradient. We derive new techniques to make continuous maps of lensing using this "Gradient-Inversion" (GI) method and validate our method with simulated data, finding good agreement with predictions. For idealized simulations of lensing cross- and autospectra that neglect foregrounds, we demonstrate that our method performs significantly better than previous quadratic estimator methods in temperature; at $L=5000-9000$, it reduces errors on the lensing auto-power spectrum by a factor of $\sim 4$ for both idealized CMB-S4 and Simons Observatory-like experiments and by a factor of $\sim 2.6$ for cross-correlations of CMB-S4-like lensing reconstruction and the true lensing field. We caution that the level of the neglected small-scale foreground power, while low in polarization, is very high in temperature; though we briefly outline foreground mitigation methods, further work on this topic is required. Nevertheless, our results show the future potential for improved small-scale CMB lensing measurements, which could provide stronger constraints on cosmological parameters and astrophysics at high redshifts., Comment: 14 pages, 4 figures

We constrain the rotational kinematic Sunyaev-Zel'dovich (rkSZ) effect in Planck data using a sample of rotating galaxy clusters identified in the Sloan Digital Sky Survey (SDSS). We extract cluster-centered cutouts from Planck cosmic microwave background (CMB) maps that have been cleaned of thermal SZ signal. Using previous constraints on the cluster rotation vectors determined from the motions of galaxies, we fit for the amplitude of the rkSZ effect in the CMB cutouts, marginalizing over parameters describing the cluster electron distribution. We also employ an alternative, less model-dependent approach to measuring the rkSZ signal that involves measuring the dipole induced by the rkSZ in rotation velocity-oriented CMB stacks. In both cases, we find roughly $2\sigma$ evidence for a rkSZ signal consistent with the expected amplitude and morphology. We comment on future directions for measurements of the rkSZ signal., Comment: 10 pages, 4 figures; comments welcome!