Accurate reconstruction of peculiar velocities from galaxy positions is important for probing the motion and evolution of large scale structure. They are sensitive to the cosmological effects of gravity and dark sector matter, and complement other velocity inference methods such as the kinematic Sunyaev-Zel'dovich (kSZ) imprinted in the CMB. We show that machine learning methods improve velocity reconstruction by capturing nonlinear contributions. Specifically, we train both a gradient boosting decision tree (GBDT) and a Transformer using multi-scale features to predict the residual between the actual velocity and the estimate from linear theory for both the line-of-sight and transverse components. We evaluate our approach in both periodic box and, more realistic, lightcone settings using mock galaxy catalogs from the AbacusSummit simulations tailored to DESI spectroscopic surveys of luminous red galaxies (LRGs) and emission line galaxies (ELGs). We also assess the impact of redshift uncertainties such as those in Rubin LSST photometry. Both models significantly outperform the linear theory with the Transformer achieving the best performance. They more accurately recover the velocity power spectrum and maintain a higher cross-correlation with the true velocities across a wider range of spatial scales. Finally, we demonstrate two applications relevant to kSZ analyses: estimating cluster pairwise velocity correlations and stacked cluster density profiles. This machine learning framework for nonlinear velocity reconstruction opens up powerful new applications of survey data from DESI, Rubin LSST, Euclid and the Roman Space Telescope.

We present a new method for measuring the EG statistic that combines two CMB secondaries-the kinematic Sunyaev-Zel'dovich (kSZ) effect and CMB lensing-for the first time to probe gravity on linear scales. The EG statistic is a discriminating tool for modified gravity theories, which leave imprints in lensing observables and peculiar velocities. Existing EG measurements rely on redshift space distortions (RSD) to infer the velocity field. Here, we employ kSZ velocity reconstruction instead of RSD, a complementary technique that constrains the largest-scale modes better than the galaxy survey it uses. We construct a novel VG estimator that involves a ratio between cross-correlations of a galaxy sample with a CMB convergence map and that with a 3D kSZ-reconstructed velocity field. We forecast for current and upcoming CMB maps from the Atacama Cosmology Telescope and the Simons Observatory (SO), respectively, in combination with three spectroscopic galaxy samples from the Dark Energy Spectroscopic Instrument (DESI). We find cumulative detection significances in the range S/N similar to 20-55, which can robustly test the scale-independent EG prediction under general relativity (GR) at different effective redshifts of the galaxy samples (z approximate to 0.73, 1.33, 1.84). In particular, the SO x DESI luminous red galaxies measurement would be able to distinguish between GR and certain modified gravity models, including Hu-Sawicki f(R) and Chameleon theories, with high confidence. The proposed VG estimator opens up a new avenue for stress-testing gravity and the ACDM & thorn; GR model at the largest observable scales.

The mean pairwise velocity of massive halos reflects the gravitational force law on cosmic scales. We combine cosmic microwave background intensity maps from the Atacama Cosmology Telescope and a galaxy catalog from the Sloan Digital Sky Survey to estimate the mean pairwise velocity using the kinematic Sunyaev-Zeldovich (kSZ) effect. On scales from 30 to 230 megaparsecs, we constrain the gravitational acceleration between pairs of halos at separation r to be g∝1/r^{n} with n=2.1±0.3, which is consistent with Newtonian gravity in an expanding spacetime (i.e., the standard ΛCDM model). This constraint shows agreement with an inverse quadratic radial dependence over the large distances separating galaxy halos, as expected in standard cosmology. Upcoming surveys have the potential to rule out n=1 at 10σ significance. Our results establish the kSZ effect as a powerful tool for testing gravity on cosmological scales.

We present a 9.3 sigma detection of the pairwise kinematic Sunyaev-Zel'dovich (kSZ) effect by combining a sample of 913,286 Luminous Red Galaxies (LRGs) from the Dark Energy Spectroscopic Instrument Data Release 1 (DESI DR1) catalog and coadded Atacama Cosmology Telescope (ACT DR6) and Planck cosmic microwave background (CMB) temperature maps. This represents the highest-significance pairwise kSZ measurement to date. The analysis uses three ACT CMB temperature maps: coadded 150 GHz, total frequency maps, and a component separated Internal Linear Combination (ILC) map, all of which cover 19,000 square degrees of the sky from Advanced ACTPol observations conducted between 2017 and 2022. Comparison of the results of these three maps serves as a consistency check for potential foreground contamination that may depend on the observation frequency. An estimate of the best-fit mass-averaged optical depth is obtained by comparing the pairwise kSZ curve with the linear theory prediction of the pairwise velocity under the best-fit Planck cosmology and is compared with predictions from simulations. This estimate serves as a reference point for future comparisons with thermal SZ-derived optical depth measurements for the same DESI cluster samples, which will be presented in a companion paper. Finally, we employ a machine learning approach, trained on simulations to estimate the optical depth for 456,803 DESI LRG-identified clusters within the simulated mass range (greater than or similar to 1013M circle dot). These are combined with the measured kSZ signal to infer the individual cluster peculiar velocities, providing the opportunity to constrain the behavior of gravity and the dark sector over a range of cosmic scales and epochs.

We present power spectra of the cosmic microwave background (CMB) anisotropy in temperature and polarization, measured from the Data Release 6 maps made from Atacama Cosmology Telescope (ACT) data. These cover 19,000 deg(2) of sky in bands centered at 98, 150 and 220 GHz, with white noise levels three times lower than Planck in polarization. We find that the ACT angular power spectra estimated over 10,000 deg(2), and measured to arcminute scales in TT, TE and EE, are well fit by the sum of CMB and foregrounds, where the CMB spectra are described by the Lambda CDM model. Combining ACT with larger-scale Planck data, the joint P-ACT dataset provides tight limits on the ingredients, expansion rate, and initial conditions of the universe. We find similar constraining power, and consistent results, from either the Planck power spectra or from ACT combined with WMAP data, as well as from either temperature or polarization in the joint P-ACT dataset. When combined with CMB lensing from ACT and Planck, and baryon acoustic oscillation data from the Dark Energy Spectroscopic Instrument (DESI DR1), we measure a baryon density of Omega(b)h(2) = 0.0226 +/- 0.0001, a cold dark matter density of Omega(c)h(2) = 0.118 +/- 0.001, a Hubble constant of H-0 = 68.22 +/- 0.36 km/s/Mpc, a spectral index of n(s) = 0.974 +/- 0.003, and an amplitude of density fluctuations of sigma(8) = 0.813 +/- 0.005. Including the DESI DR2 data tightens the Hubble constant to H-0 = 68.43 +/- 0.27 km/s/Mpc; Lambda CDM parameters agree between the P-ACT and DESI DR2 data at the 1.6 sigma level. We find no evidence for excess lensing in the power spectrum, and no departure from spatial flatness. The contribution from Sunyaev-Zel'dovich (SZ) anisotropy is detected at high significance; we find evidence for a tilt with suppressed small-scale power compared to our baseline SZ template spectrum, consistent with hydrodynamical simulations with feedback.

We present a new method for measuring the $E_G$ statistic that combines two CMB secondaries -- the kinematic Sunyaev-Zeldovich (kSZ) effect and CMB lensing -- for the first time to probe gravity on linear scales. The $E_G$ statistic is a discriminating tool for modified gravity theories, which leave imprints in lensing observables and peculiar velocities. Existing $E_G$ measurements rely on redshift space distortions (RSD) to infer the velocity field. Here, we employ kSZ velocity-reconstruction instead of RSD, a complementary technique that constrains the largest-scale modes better than the galaxy survey it uses. We construct a novel $\widehat{V}_G$ estimator that involves a ratio between cross-correlations of a galaxy sample with a CMB convergence map and that with a 3D kSZ-reconstructed velocity field. We forecast for current and upcoming CMB maps from the Atacama Cosmology Telescope (ACT) and the Simons Observatory (SO), respectively, in combination with three spectroscopic galaxy samples from the Dark Energy Spectroscopic Instrument (DESI). We find cumulative detection significances in the range $S/N \sim 20-55$, which can robustly test the scale-independent $E_G$ prediction under general relativity (GR) at different effective redshifts of the galaxy samples ($z\approx 0.73, 1.33, 1.84$). In particular, the SO$\times$DESI LRG measurement would be able to distinguish between GR and certain modified gravity models, including Hu-Sawicki $f(R)$ and Chameleon theories, with high confidence. The proposed $\widehat{V}_G$ estimator opens up a new avenue for stress-testing gravity and the $Λ$CDM+GR model at the largest observable scales.

We systematically analyze the implications of statistical noise within numerical derivatives on simulation-based Fisher forecasts for large scale structure surveys. Noisy numerical derivatives resulting from a finite number of simulations, Nsims, act to bias the associated Fisher forecast such that the resulting marginalized constraints can be significantly tighter than the noise-free limit. We show the source of this effect can be traced to the influence of the noise on the marginalization process. Parameters such as the neutrino mass, M nu, for which higher-order forward differentiation schemes are commonly used, are more prone to noise; the predicted constraints can be akin to those purely from a random instance of statistical noise even using (1 Gpc/h)3 simulations with Nsims = 500 realizations. We demonstrate how derivative noise can artificially reduce parameter degeneracies and seemingly null the effects of adding nuisance parameters to the forecast, such as halo occupation distribution fitting parameters. We mathematically characterize these effects through a full statistical analysis, and demonstrate how confidence intervals for the true noise-free, Nsims -> infinity, Fisher constraints can be recovered even when noise comprises a consequential component of the measured signal. The findings and approaches developed here are important for ensuring simulation-based analyses can be accurately used to assess upcoming survey capabilities.

We describe updated scientific goals for the wide-field, millimeter-wave survey that will be produced by the Simons Observatory (SO). Significant upgrades to the 6-meter SO Large Aperture Telescope (LAT) are expected to be complete by 2028, and will include a doubled mapping speed with 30,000 new detectors and an automated data reduction pipeline. In addition, a new photovoltaic array will supply most of the observatory's power. The LAT survey will cover about 60% of the sky at a regular observing cadence, with five times the angular resolution and ten times the map depth of the Planck satellite. The science goals are to: (1) determine the physical conditions in the early universe and constrain the existence of new light particles; (2) measure the integrated distribution of mass, electron pressure, and electron momentum in the late-time universe, and, in combination with optical surveys, determine the neutrino mass and the effects of dark energy via tomographic measurements of the growth of structure at redshifts z ≲ 3; (3) measure the distribution of electron density and pressure around galaxy groups and clusters, and calibrate the effects of energy input from galaxy formation on the surrounding environment; (4) produce a sample of more than 30,000 galaxy clusters, and more than 100,000 extragalactic millimeter sources, including regularly sampled AGN light-curves, to study these sources and their emission physics; (5) measure the polarized emission from magnetically aligned dust grains in our Galaxy, to study the properties of dust and the role of magnetic fields in star formation; (6) constrain asteroid regoliths, search for Trans-Neptunian Objects, and either detect or eliminate large portions of the phase space in the search for Planet 9; and (7) provide a powerful new window into the transient universe on time scales of minutes to years, concurrent with observations from the Vera C. Rubin Observatory of overlapping sky.

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.

We apply two machine learning methods, a CNN deep-learning model and a gradient-boosting decision tree, to estimate individual cluster optical depths from observed properties derived from multiple complementary datasets. The models are trained and tested with simulated N-body derived halo catalogs and synthetic full-sky cosmic microwave background (CMB) maps designed to mirror data from the DESI and Simons Observatory experiments. Specifically, the thermal Sunyaev-Zel'dovich (tSZ) and CMB lensing convergence along with cluster virial mass estimates are used as features to train the machine learning models. The predicted optical depths are combined with kinematic Sunyaev-Zel'dovich (kSZ) measurements to estimate individual cluster radial peculiar velocities. The method is shown to recover an unbiased estimate of the pairwise velocity statistics of the simulated cluster sample. The models are demonstrated to be effective for halos with mass range 10(13)M(circle dot) < M-200 < 10(15)M(circle dot) over a redshift range 0 < z < 1, and validated in the presence of primary CMB, instrument noise, lensing convergence noise, and potential uncertainties in halo virial mass estimates. We apply the method to Atacama Cosmology Telescope (ACT) CMB data, using ACT DR4 component-separated maps for tSZ and CMB lensing and ACT DR5 maps for kSZ, in conjunction with galaxy clusters observed in the SDSS DR15 spectroscopic survey. We demonstrate that the machine learning approach is an effective one to analyze data from current and upcoming CMB experiments such as Simons Observatory and Cerro Chajnantor Atacama Telescope, and galaxy surveys, such as DESI and Roman, for which the pairwise velocity statistics can provide valuable insights into the properties of neutrinos and gravity on cosmic scales.

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σ. Our work presents a first step towards the use of this observable for cosmological analyses.

We present Atacama Cosmology Telescope (ACT) Data Release 6 (DR6) maps of the Cosmic Microwave Background temperature and polarization anisotropy at arcminute resolution over three frequency bands centered on 98, 150 and 220 GHz. The maps are based on data collected with the AdvancedACT camera over the period 2017-2022 and cover 19,000 square degrees with a median combined depth of 10 mu K arcmin. We describe the instrument, mapmaking and map properties and illustrate them with a number of figures and tables. The ACT DR6 maps and derived products are available on LAMBDA at https://lambda.gsfc.nasa.gov/product/act/actadv_prod_table.html. We also provide an interactive web atlas at https://phy-act1.princeton.edu/public/snaess/actpol/dr6/atlas and HiPS data sets in Aladin (e.g. https://alasky.cds.unistra.fr/ACT/DR4DR6/color_CMB).

We use new cosmic microwave background (CMB) primary temperature and polarization anisotropy measurements from the Atacama Cosmology Telescope (ACT) Data Release 6 (DR6) to test foundational assumptions of the standard cosmological model, Lambda CDM, and set constraints on extensions to it. We derive constraints from the ACT DR6 power spectra alone, as well as in combination with legacy data from the Planck mission. To break geometric degeneracies, we include ACT and Planck CMB lensing data and baryon acoustic oscillation data from DESI Year-1. To test the dependence of our results on non-ACT data, we also explore combinations replacing Planck with WMAP and DESI with BOSS, and further add supernovae measurements from Pantheon+ for models that affect the late-time expansion history. We verify the near-scale-invariance (running of the spectral index dn(s)/d ln k = 0.0062 +/- 0.0052) and adiabaticity of the primordial perturbations. Neutrino properties are consistent with Standard Model predictions: we find no evidence for new light, relativistic species that are free-streaming (N-eff = 2.86 +/- 0.13, which combined with astrophysical measurements of primordial helium and deuterium abundances becomes N-eff = 2.89 +/- 0.11), for non-zero neutrino masses (Sigma m(nu) < 0.089 eV at 95% CL), or for neutrino self-interactions. We also find no evidence for self-interacting dark radiation (N-idr < 0.134), or for early-universe variation of fundamental constants, including the fine-structure constant (alpha(EM)/alpha(EM,0) = 1.0043 +/- 0.0017) and the electron mass (m(e)/m(e,0) = 1.0063 +/- 0.0056). Our data are consistent with standard big bang nucleosynthesis (we find Y-p = 0.2312 +/- 0.0092), the COBE/FIRAS-inferred CMB temperature (we find T-CMB = 2.698 +/- 0.016 K), a dark matter component that is collisionless and with only a small fraction allowed as axion-like particles, a cosmological constant (w = -0.986 +/- 0.025), and the late-time growth rate predicted by general relativity (gamma = 0.663 +/- 0.052). We find no statistically significant preference for a departure from the baseline Lambda CDM model. In fits to models invoking early dark energy, primordial magnetic fields, or an arbitrary modified recombination history, we find H-0 = 69.9(-1.5)(+0.8), 69.1 +/- 0.5, or 69.6 +/- 1.0 km/s/Mpc, respectively; using BOSS instead of DESI BAO data reduces the central values of these constraints by 1-1.5 km/s/Mpc while only slightly increasing the error bars. In general, models introduced to increase the Hubble constant or to decrease the amplitude of density fluctuations inferred from the primary CMB are not favored over Lambda CDM by our data.

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.

We demonstrate that pairwise peculiar velocity correlations for galaxy clusters can be directly reconstructed from the kinematic Sunyaev-Zel'dovich (kSZ) signature imprinted in the cosmic microwave background (CMB) using a machine learning model with a gradient boosting algorithm trained on highfidelity kSZ simulations. The machine learning model is trained using six to seven cluster features that are directly related to observables from CMB and large-scale structure surveys. We validate the capabilities of the approach in light of the presence of primary CMB, detector noise, and potential uncertainties in the cluster mass estimate and cluster center location. The pairwise velocity statistics extracted using the techniques developed here have the potential to elicit valuable cosmological constraints on dark energy, modified gravity models, and massive neutrinos with kSZ measurements from upcoming CMB surveys, including the Simons Observatory, CMB-S4 and CCAT, and the DESI and SDSS galaxy surveys.

We systematically analyze the implications of statistical noise within numerical derivatives on simulation-based Fisher forecasts for large scale structure surveys. Noisy numerical derivatives resulting from a finite number of simulations, N_sims, act to bias the associated Fisher forecast such that the resulting marginalized constraints can be significantly tighter than the noise-free limit. We show the source of this effect can be traced to the influence of the noise on the marginalization process. Parameters such as the neutrino mass, , for which higher-order forward differentiation schemes are commonly used, are more prone to noise; the predicted constraints can be akin to those purely from a random instance of statistical noise even using (1Gpc/h)^3 simulations with N_sims=500 realizations. We demonstrate how derivative noise can artificially reduce parameter degeneracies and seemingly null the effects of adding nuisance parameters to the forecast, such as HOD fitting parameters. We mathematically characterize these effects through a full statistical analysis, and demonstrate how confidence intervals for the true noise-free, N_sims→∞, Fisher constraints can be recovered even when noise comprises a consequential component of the measured signal. The findings and approaches developed here are important for ensuring simulation-based analyses can be accurately used to assess upcoming survey capabilities.

We determine the efficacy of the kinematic Sunyaev-Zel'dovich (kSZ) signal extraction pipeline, using pairwise kSZ measurements, in recovering unbiased estimates of the signal and inference of the associated optical depth. We consider the impact of cluster coalignments along the line of sight, the modeling of baryonic clustering, and the presence of diffuse gas, as well as instrument beam convolution and noise. We demonstrate that two complementary approaches, aperture photometry, and a matched filter, can be used to recover an unbiased estimate of the cluster kSZ signal and the associated optical depth. Aperture photometry requires a correction factor accounting for the subtraction of signal in the annulus while the matched filter requires a tuning of the signal template profile. We show that both of these can be calibrated from simulated survey data. The optical depth estimates are also consistent with those inferred from stacked thermal SZ measurements. We apply the approaches to the publicly available Atacama Cosmology Telescope (ACT) data. The techniques developed here provide a promising method to leverage upcoming kSZ measurements, from ACT, Simons Observatory, CCAT, and CMB-S4 with spectroscopic galaxy surveys from DESI, Euclid, and Roman, to constrain cosmological properties of the dark energy, gravity, and neutrino masses.

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.

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

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