Search

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

At frequencies below 1\,Hz, fluctuations in atmospheric emission in the Chajnantor region in northern Chile are the primary source of interference for bolometric millimeter-wave observations. This paper focuses on the statistics of these fluctuations using measurements from the Atacama Cosmology Telescope (ACT) and the Atacama Pathfinder Experiment (APEX) water vapor radiometer. After introducing a method for separating atmospheric effects from other systematic effects, we present a direct measurement of the temporal outer scale of turbulence of $\tau_0\approx50$s corresponding to a spatial scale of $L_0\approx500$m. At smaller scales, the fluctuations are well described by the Kolmogorov 2/3 power law until, at yet smaller scales, the effects of beam smearing become important. As a part of this study, we present measurements of the atmosphere by the APEX radiometer over 20 years, focused on fluctuations in precipitable water vapor (PWV). We find that the 30-minute mean of the total PWV is not in general a robust estimator of the level of fluctuations. We show that the microwave frequency spectrum of these fluctuations is in good agreement with predictions by the \texttt{am} code for bands above 90~GHz. We then show that the variance of fluctuations in ACT's mm-wave bands correlates with the variance of fluctuations in PWV measured by APEX, even though the observatories are 6\,km apart and observe different lines of sight. We find that ACT's atmosphere-determined optical efficiencies are consistent with previous planet-based results., Comment: 13 pages plus appendix, 17 figures. Submitted to PRD

We present $\texttt{cunusht}$, a general-purpose Python package that wraps a highly efficient CUDA implementation of the nonuniform spin-$0$ spherical harmonic transform. The method is applicable to arbitrary pixelization schemes, including schemes constructed from equally-spaced iso-latitude rings as well as completely nonuniform ones. The algorithm has an asymptotic scaling of $\mathrm{O}{(\ell_{\rm max}^3)}$ for maximum multipole $\ell_{\rm max}$ and achieves machine precision accuracy. While $\texttt{cunusht}$ is developed for applications in cosmology in mind, it is applicable to various other interpolation problems on the sphere. We outperform the fastest available CPU algorithm by a factor of up to 5 for problems with a nonuniform pixelization and $\ell_{\rm max}>4\cdot10^3$ when comparing a single modern GPU to a modern 32-core CPU. This performance is achieved by utilizing the double Fourier sphere method in combination with the nonuniform fast Fourier transform and by avoiding transfers between the host and device. For scenarios without GPU availability, $\texttt{cunusht}$ wraps existing CPU libraries. $\texttt{cunusht}$ is publicly available and includes tests, documentation, and demonstrations.

Optimal extraction of cosmological information from observations of the Cosmic Microwave Background critically relies on our ability to accurately undo the distortions caused by weak gravitational lensing. In this work, we demonstrate the use of denoising diffusion models in performing Bayesian lensing reconstruction. We show that score-based generative models can produce accurate, uncorrelated samples from the CMB lensing convergence map posterior, given noisy CMB observations. To validate our approach, we compare the samples of our model to those obtained using established Hamiltonian Monte Carlo methods, which assume a Gaussian lensing potential. We then go beyond this assumption of Gaussianity, and train and validate our model on non-Gaussian lensing data, obtained by ray-tracing N-body simulations. We demonstrate that in this case, samples from our model have accurate non-Gaussian statistics beyond the power spectrum. The method provides an avenue towards more efficient and accurate lensing reconstruction, that does not rely on an approximate analytic description of the posterior probability. The reconstructed lensing maps can be used as an unbiased tracer of the matter distribution, and to improve delensing of the CMB, resulting in more precise cosmological parameter inference., Comment: 12 pages, 10 figures. v2: typo in one of the equations fixed, references added

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

We 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

We present a cross-correlation analysis between $1'$ resolution total intensity and polarization observations from the Atacama Cosmology Telescope (ACT) at 150 and 220 GHz and 15$''$ mid-infrared photometry from the Wide-field Infrared Survey Explorer (WISE) over 107 12.5$^\circ\times$12.5$^\circ$ patches of sky. We detect a spatially isotropic signal in the WISE$\times$ACT $TT$ cross power spectrum at 30$\sigma$ significance that we interpret as the correlation between the cosmic infrared background at ACT frequencies and polycyclic aromatic hydrocarbon (PAH) emission from galaxies in WISE, i.e., the cosmic PAH background. Within the Milky Way, the Galactic dust $TT$ spectra are generally well-described by power laws in $\ell$ over the range 10$^3 < \ell < $10$^4$, but there is evidence both for variability in the power law index and for non-power law behavior in some regions. We measure a positive correlation between WISE total intensity and ACT $E$-mode polarization at 1000$ < \ell \lesssim $6000 at $>$3$\sigma$ in each of 35 distinct $\sim$100 deg$^2$ regions of the sky, suggesting alignment between Galactic density structures and the local magnetic field persists to sub-parsec physical scales in these regions. The distribution of $TE$ amplitudes in this $\ell$ range across all 107 regions is biased to positive values, while there is no evidence for such a bias in the $TB$ spectra. This work constitutes the highest-$\ell$ measurements of the Galactic dust $TE$ spectrum to date and indicates that cross-correlation with high-resolution mid-infrared measurements of dust emission is a promising tool for constraining the spatial statistics of dust emission at millimeter wavelengths., Comment: 20 pages, 14 figures, submitted to ApJ

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

We use time-domain simulations of Jupiter observations to test and develop a beam reconstruction pipeline for the Simons Observatory Small Aperture Telescopes. The method relies on a map maker that estimates and subtracts correlated atmospheric noise and a beam fitting code designed to compensate for the bias caused by the map maker. We test our reconstruction performance for four different frequency bands against various algorithmic parameters, atmospheric conditions and input beams. We additionally show the reconstruction quality as function of the number of available observations and investigate how different calibration strategies affect the beam uncertainty. For all of the cases considered, we find good agreement between the fitted results and the input beam model within a ~1.5% error for a multipole range l = 30 - 700 and an ~0.5% error for a multipole range l = 50 - 200. We conclude by using a harmonic-domain component separation algorithm to verify that the beam reconstruction errors and biases observed in our analysis do not significantly bias the Simons Observatory r-measurement., Comment: 23 pages, 21 figures, published in ApJ

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

The increasing statistical power of cosmic microwave background (CMB) datasets requires a commensurate effort in understanding their noise properties. The noise in maps from ground-based instruments is dominated by large-scale correlations, which poses a modeling challenge. This paper develops novel models of the complex noise covariance structure in the Atacama Cosmology Telescope Data Release 6 (ACT DR6) maps. We first enumerate the noise properties that arise from the combination of the atmosphere and the ACT scan strategy. We then prescribe a class of Gaussian, map-based noise models, including a new wavelet-based approach that uses directional wavelet kernels for modeling correlated instrumental noise. The models are empirical, whose only inputs are a small number of independent realizations of the same region of sky. We evaluate the performance of these models against the ACT DR6 data by drawing ensembles of noise realizations. Applying these simulations to the ACT DR6 power spectrum pipeline reveals a $\sim 20\%$ excess in the covariance matrix diagonal when compared to an analytic expression that assumes noise properties are uniquely described by their power spectrum. Along with our public code, $\mathtt{mnms}$, this work establishes a necessary element in the science pipelines of both ACT DR6 and future ground-based CMB experiments such as the Simons Observatory (SO)., Comment: 41 pages (+10 appendix), 22 figures (+5 appendix), submitted to JCAP

Polarization of the cosmic microwave background (CMB) can probe new parity-violating physics such as cosmic birefringence (CB), which requires exquisite control over instrumental systematics. The non-idealities of the half-wave plate (HWP) represent a source of systematics when used as a polarization modulator. We study their impact on the CMB angular power spectra, which is partially degenerate with CB and miscalibration of the polarization angle. We use full-sky beam convolution simulations including HWP to generate mock noiseless time-ordered data, process them through a bin averaging map-maker, and calculate the power spectra including $TB$ and $EB$ correlations. We also derive analytical formulae which accurately model the observed spectra. For our choice of HWP parameters, the HWP-induced angle amounts to a few degrees, which could be misinterpreted as CB. Accurate knowledge of the HWP is required to mitigate this. Our simulation and analytical formulae will be useful for deriving requirements for the accuracy of HWP calibration., Comment: 12 pages + appendices and bibliography, 5 figures, 1 table; accepted for publication in JCAP

Microwave telescopes require an ever-increasing control of experimental systematics in their quest to measure the Cosmic Microwave Background (CMB) to exquisite levels of precision. One important systematic for ground and balloon-borne experiments is ground pickup, where beam sidelobes detect the thermal emission of the much warmer ground while the main beam is scanning the sky. This generates scan-synchronous noise in experiment timestreams, which is difficult to filter out without also deleting some of the signal from the sky. Therefore, efficient modelling of pickup can help guide the design of experiments and of analysis pipelines. In this work, we present an extension to the beamconv algorithm that enables us to generate time-ordered data (TOD) from beam-convolved sky and ground maps simultaneously. We simulate ground pickup for both a ground-based experiment and a telescope attached to a stratospheric balloon. Ground templates for the balloon experiment are obtained by re-projecting satellite maps of the Earth's microwave emission., Comment: 10 pages, 6 figures, paper 12190-165 for the "SPIE AT+I 2022 mm, sub-mmm, and Far-Infrared Detectors and Instrumentation for Astronomy XI" conference

We describe the measurement and treatment of the telescope beams for the Atacama Cosmology Telescope's fourth data release, DR4. Observations of Uranus are used to measure the central portion (<12') of the beams to roughly -40 dB of the peak. Such planet maps in intensity are used to construct azimuthally averaged beam profiles, which are fit with a physically motivated model before being transformed into Fourier space. We investigate and quantify a number of percent-level corrections to the beams, all of which are important for precision cosmology. Uranus maps in polarization are used to measure the temperature-to-polarization leakage in the main part of the beams, which is <1% (2.5%) at 150 GHz (98 GHz). The beams also have polarized sidelobes, which are measured with observations of Saturn and deprojected from the ACT time-ordered data. Notable changes relative to past ACT beam analyses include an improved subtraction of the atmospheric effects from Uranus calibration maps, incorporation of a scattering term in the beam profile model, and refinements to the beam model uncertainties and the main temperature-to-polarization leakage terms in the ACT power spectrum analysis., Comment: 32 pages, 21 figures

Fluctuating atmospheric emission is a dominant source of noise for ground-based millimeter-wave observations of the CMB temperature anisotropy at angular scales $\gtrsim 0.5^{\circ}$. We present a model of the atmosphere as a discrete set of emissive turbulent layers that move with respect to the observer with a horizontal wind velocity. After introducing a statistic derived from the time-lag dependent correlation function for detector pairs in an array, referred to as the pair-lag, we use this model to estimate the aggregate angular motion of the atmosphere derived from time-ordered data from the Atacama Cosmology Telescope (ACT). We find that estimates derived from ACT's CMB observations alone agree with those derived from satellite weather data that additionally include a height-dependent horizontal wind velocity and water vapor density. We also explore the dependence of the measured atmospheric noise spectrum on the relative angle between the wind velocity and the telescope scan direction. In particular, we find that varying the scan velocity changes the noise spectrum in a predictable way. Computing the pair-lag statistic opens up new avenues for understanding how atmospheric fluctuations impact measurements of the CMB anisotropy.

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

We report a significant detection of the hot intergalactic medium in the filamentary bridge connecting the galaxy clusters Abell 399 and Abell 401. This result is enabled by a low-noise, high-resolution map of the thermal Sunyaev-Zeldovich signal from the Atacama Cosmology Telescope (ACT) and Planck satellite. The ACT data provide the $1.65'$ resolution that allows us to clearly separate the profiles of the clusters, whose centres are separated by $37'$, from the gas associated with the filament. A model that fits for only the two clusters is ruled out compared to one that includes a bridge component at $>5\sigma$. Using a gas temperature determined from Suzaku X-ray data, we infer a total mass of $(3.3\pm0.7)\times10^{14}\,\mathrm{M}_{\odot}$ associated with the filament, comprising about $8\%$ of the entire Abell 399-Abell 401 system. We fit two phenomenological models to the filamentary structure; the favoured model has a width transverse to the axis joining the clusters of ${\sim}1.9\,\mathrm{Mpc}$. When combined with the Suzaku data, we find a gas density of $(0.88\pm0.24)\times10^{-4}\,\mathrm{cm}^{-3}$, considerably lower than previously reported. We show that this can be fully explained by a geometry in which the axis joining Abell 399 and Abell 401 has a large component along the line of sight, such that the distance between the clusters is significantly greater than the $3.2\,\mathrm{Mpc}$ projected separation on the plane of the sky. Finally, we present initial results from higher resolution ($12.7"$ effective) imaging of the bridge with the MUSTANG-2 receiver on the Green Bank Telescope., Comment: 23 pages, 10 figures, 3 tables. This is a pre-copyedited, author-produced PDF of an article accepted for publication in the Monthly Notices of the Royal Astronomical Society following peer review. The version of record is available online at: https://academic.oup.com/mnras/advance-article/doi/10.1093/mnras/stab3391/6442294

We present arcminute-resolution intensity and polarization maps of the Galactic center made with the Atacama Cosmology Telescope (ACT). The maps cover a 32 deg$^2$ field at 98, 150, and 224 GHz with $\vert l\vert\le4^\circ$, $\vert b\vert\le2^\circ$. We combine these data with Planck observations at similar frequencies to create coadded maps with increased sensitivity at large angular scales. With the coadded maps, we are able to resolve many known features of the Central Molecular Zone (CMZ) in both total intensity and polarization. We map the orientation of the plane-of-sky component of the Galactic magnetic field inferred from the polarization angle in the CMZ, finding significant changes in morphology in the three frequency bands as the underlying dominant emission mechanism changes from synchrotron to dust emission. Selected Galactic center sources, including Sgr A*, the Brick molecular cloud (G0.253+0.016), the Mouse pulsar wind nebula (G359.23-0.82), and the Tornado supernova remnant candidate (G357.7-0.1), are examined in detail. These data illustrate the potential for leveraging ground-based Cosmic Microwave Background polarization experiments for Galactic science., Comment: 26 pages, 14 figures, accepted for publication in ApJ

We use Atacama Cosmology Telescope (ACT) observations at 98 GHz (2015--2019), 150 GHz (2013--2019) and 229 GHz (2017--2019) to perform a blind shift-and-stack search for Planet 9. The search explores distances from 300 AU to 2000 AU and velocities up to 6.3 arcmin per year, depending on the distance. For a 5 Earth-mass Planet 9 the detection limit varies from 325 AU to 625 AU, depending on the sky location. For a 10 Earth-mass planet the corresponding range is 425 AU to 775 AU. The search covers the whole 18,000 square degrees of the ACT survey, though a slightly deeper search is performed for the parts of the sky consistent with Planet 9's expected orbital inclination. No significant detections are found, which is used to place limits on the mm-wave flux density of Planet 9 over much of its orbit. Overall we eliminate roughly 17% and 9% of the parameter space for a 5 and 10 Earth-mass Planet 9 respectively. We also provide a list of the 10 strongest candidates from the search for possible follow-up. More generally, we exclude (at 95% confidence) the presence of an unknown Solar system object within our survey area brighter than 4--12 mJy (depending on position) at 150 GHz with current distance $300 \text{ AU} < r < 600 \text{ AU}$ and heliocentric angular velocity $1.5'/\text{yr} < v \cdot \frac{500 \text{ AU}}{r} < 2.3'\text{yr}$, corresponding to low-to-moderate eccentricities. These limits worsen gradually beyond 600 AU, reaching 5--15 mJy by 1500 AU., Comment: 23 pages, 10 figures, 5 tables, submitted to ApJ

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/

We present a 5.4$\sigma$ detection of the pairwise kinematic Sunyaev-Zel'dovich (kSZ) effect using Atacama Cosmology Telescope (ACT) and $\it{Planck}$ CMB observations in combination with Luminous Red Galaxy samples from the Sloan Digital Sky Survey (SDSS) DR15 catalog. Results are obtained using three ACT CMB maps: co-added 150 GHz and 98 GHz maps, combining observations from 2008-2018 (ACT DR5), which overlap with SDSS DR15 over 3,700 sq. deg., and a component-separated map using night-time only observations from 2014-2015 (ACT DR4), overlapping with SDSS DR15 over 2,089 sq. deg. Comparisons of the results from these three maps provide consistency checks in relation to potential frequency-dependent foreground contamination. A total of 343,647 galaxies are used as tracers to identify and locate galaxy groups and clusters from which the kSZ signal is extracted using aperture photometry. We consider the impact of various aperture photometry assumptions and covariance estimation methods on the signal extraction. Theoretical predictions of the pairwise velocities are used to obtain best-fit, mass-averaged, optical depth estimates for each of five luminosity-selected tracer samples. A comparison of the kSZ-derived optical depth measurements obtained here to those derived from the thermal SZ effect for the same sample is presented in a companion paper., Comment: 17 pages, 10 figures. Updated to match published version in PRD

We present high signal-to-noise measurements (up to 12$\sigma$) of the average thermal Sunyaev Zel'dovich (tSZ) effect from optically selected galaxy groups and clusters and estimate their baryon content within a 2.1$^\prime$ radius aperture. Sources from the Sloan Digital Sky Survey (SDSS) Baryon Oscillation Spectroscopic Survey (BOSS) DR15 catalog overlap with 3,700 sq. deg. of sky observed by the Atacama Cosmology Telescope (ACT) from 2008 to 2018 at 150 and 98 GHz (ACT DR5), and 2,089 sq. deg. of internal linear combination component-separated maps combining ACT and $\it{Planck}$ data (ACT DR4). The corresponding optical depths, $\bar{\tau}$, which depend on the baryon content of the halos, are estimated using results from cosmological hydrodynamic simulations assuming an AGN feedback radiative cooling model. We estimate the mean mass of the halos in multiple luminosity bins, and compare the tSZ-based $\bar{\tau}$ estimates to theoretical predictions of the baryon content for a Navarro-Frenk-White profile. We do the same for $\bar{\tau}$ estimates extracted from fits to pairwise baryon momentum measurements of the kinematic Sunyaev-Zel'dovich effect (kSZ) for the same data set obtained in a companion paper. We find that the $\bar{\tau}$ estimates from the tSZ measurements in this work and the kSZ measurements in the companion paper agree within $1\sigma$ for two out of the three disjoint luminosity bins studied, while they differ by 2-3$\sigma$ in the highest luminosity bin. The optical depth estimates account for one third to all of the theoretically predicted baryon content in the halos across luminosity bins. Potential systematic uncertainties are discussed. The tSZ and kSZ measurements provide a step towards empirical Compton-$\bar{y}$-$\bar{\tau}$ relationships to provide new tests of cluster formation and evolution models., Comment: 19 pages, 9 figures

We report on the serendipitous discovery of three transient mm-wave sources using data from the Atacama Cosmology Telescope. The first, detected at RA = 273.8138, dec = -49.4628 at ${\sim}50\sigma$ total, brightened from less than 5 mJy to at least 1100 mJy at 150 GHz with an unknown rise time shorter than thirteen days, during which the increase from 250 mJy to 1100 mJy took only 8 minutes. Maximum flux was observed on 2019-11-8. The source's spectral index in flux between 90 and 150 GHz was positive, $\alpha = 1.5\pm0.2$. The second, detected at RA = 105.1584, dec = -11.2434 at ${\sim}20\sigma$ total, brightened from less than 20 mJy to at least 300 mJy at 150 GHz with an unknown rise time shorter than eight days. Maximum flux was observed on 2019-12-15. Its spectral index was also positive, $\alpha = 1.8\pm0.2$. The third, detected at RA = 301.9952, dec = 16.1652 at ${\sim}40\sigma$ total, brightened from less than 8 mJy to at least 300 mJy at 150 GHz over a day or less but decayed over a few days. Maximum flux was observed on 2018-9-11. Its spectrum was approximately flat, with a spectral index of $\alpha = -0.2\pm0.1$. None of the sources were polarized to the limits of these measurements. The two rising-spectrum sources are coincident in position with M and K stars, while the third is coincident with a G star., Comment: 8 pages, 4 figures, 1 table

We study systematic effects from half-wave plates (HWPs) for cosmic microwave background (CMB) experiments using full-sky time-domain beam convolution simulations. Using an optical model for a fiducial spaceborne two-lens refractor telescope, we investigate how different HWP configurations optimized for dichroic detectors centred at 95 and 150 GHz impact the reconstruction of primordial B-mode polarization. We pay particular attention to possible biases arising from the interaction of frequency dependent HWP non-idealities with polarized Galactic dust emission and the interaction between the HWP and the instrumental beam. To produce these simulations, we have extended the capabilities of the publicly available beamconv code. To our knowledge, we produce the first time-domain simulations that include both HWP non-idealities and realistic full-sky beam convolution. Our analysis shows how certain achromatic HWP configurations produce significant systematic polarization angle offsets that vary for sky components with different frequency dependence. Our analysis also demonstrates that once we account for interactions with HWPs, realistic beam models with non-negligible cross-polarization and sidelobes will cause significant B-mode residuals that will have to be extensively modelled in some cases.

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/

The thermal and kinematic Sunyaev-Zel'dovich effects (tSZ, kSZ) probe the thermodynamic properties of the circumgalactic and intracluster medium (CGM and ICM) of galaxies, groups, and clusters, since they are proportional, respectively, to the integrated electron pressure and momentum along the line-of-sight. We present constraints on the gas thermodynamics of CMASS galaxies in the Baryon Oscillation Spectroscopic Survey (BOSS) using new measurements of the kSZ and tSZ signals obtained in a companion paper. Combining kSZ and tSZ measurements, we measure within our model the amplitude of energy injection $\epsilon M_\star c^2$, where $M_\star$ is the stellar mass, to be $\epsilon=(40\pm9)\times10^{-6}$, and the amplitude of the non-thermal pressure profile to be $\alpha_{\rm Nth}<0.2$ (2$\sigma$), indicating that less than 20% of the total pressure within the virial radius is due to a non-thermal component. We estimate the effects of including baryons in the modeling of weak-lensing galaxy cross-correlation measurements using the best-fit density profile from the kSZ measurement. Our estimate reduces the difference between the original theoretical model and the weak-lensing galaxy cross-correlation measurements in arXiv:1611.08606 by half but does not fully reconcile it. Comparing the tSZ measurements to cosmological simulations, we find that simulations underestimate the CGM pressure at large radii while they fare better in comparison with the kSZ measurements. This suggests that the energy injected via feedback models in the simulations that we compared against does not sufficiently heat the gas at these radii. We do not find significant disagreement at smaller radii. These measurements provide novel tests of current and future simulations. This work demonstrates the power of joint, high signal-to-noise kSZ and tSZ observations, upon which future cross-correlation studies will improve., Comment: Corrected error in the algorithm that calculates the kSZ temperature profile for a given GNFW density model. The value of $\log_{\rm10} \rho_0$ changed by 0.75$\sigma$ in Tab.II and Fig.2, and affected the results in the left panels of Fig. 6. Conclusions are unchanged. Erratum published at https://link.aps.org/doi/10.1103/PhysRevD.107.049903

The scattering of cosmic microwave background (CMB) photons off the free-electron gas in galaxies and clusters leaves detectable imprints on high resolution CMB maps: the thermal and kinematic Sunyaev-Zel'dovich effects (tSZ and kSZ respectively). We use combined microwave maps from the Atacama Cosmology Telescope (ACT) DR5 and Planck in combination with the CMASS and LOWZ galaxy catalogs from the Baryon Oscillation Spectroscopic Survey (BOSS DR10 and DR12), to study the gas associated with these galaxy groups. Using individual reconstructed velocities, we perform a stacking analysis and reject the no-kSZ hypothesis at 6.5$\sigma$, the highest significance to date. This directly translates into a measurement of the electron number density profile, and thus of the gas density profile. Despite the limited signal to noise, the measurement shows at high significance that the gas density profile is more extended than the dark matter density profile, for any reasonable baryon abundance (formally $>90\sigma$ for the cosmic baryon abundance). We simultaneously measure the tSZ signal, i.e. the electron thermal pressure profile of the same CMASS objects, and reject the no-tSZ hypothesis at 10$\sigma$. We combine tSZ and kSZ measurements to estimate the electron temperature to 20% precision in several aperture bins, and find it comparable to the virial temperature. In a companion paper, we analyze these measurements to constrain the gas thermodynamics and the properties of feedback inside galaxy groups. We present the corresponding LOWZ measurements in this paper, ruling out a null kSZ (tSZ) signal at 2.9 (13.9)$\sigma$, and leave their interpretation to future work. Our stacking software ThumbStack is publicly available at https://github.com/EmmanuelSchaan/ThumbStack and directly applicable to future Simons Observatory and CMB-S4 data., Comment: Accepted in Physical Review D, Editors' Suggestion

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

The next generation of ground-based CMB experiments aim to measure temperature and polarization fluctuations up to $\ell_{\rm max} \approx 5000$ over half of the sky. Combined with Planck data on large scales, this will provide improved constraints on primordial non-Gaussianity. However, the impressive resolution of these experiments will come at a price. Besides signal confusion from galactic foregrounds, extra-galactic foregrounds and late-time gravitational effects, gravitational lensing will introduce large non-Gaussianity that can become the leading contribution to the bispectrum covariance through the connected 4-point function. Here, we compute this effect analytically for the first time on the full sky for both temperature and polarization. We compare our analytical results with those obtained directly from map-based simulations of the CMB sky for several levels of instrumental noise. Of the standard shapes considered in the literature, the local shape is most affected, resulting in a 35\% increase of the estimator standard deviation for an experiment like the Simons Observatory (SO) and a 110\% increase for a cosmic-variance limited experiment, including both temperature and polarization modes up to $\ell_{\rm max} = 3800$. Because of the nature of the lensing 4-point function, the impact on other shapes is reduced while still non negligible for the orthogonal shape. Two possible avenues to reduce the non-Gaussian contribution to the covariance are proposed. First by marginalizing over lensing contributions, such as the ISW-lensing 3pt function in temperature, and second by delensing the CMB. We show the latter method can remove almost all extra covariance, reducing the effect to below $<$5\% for local bispectra. At the same time, delensing would remove signal biases from secondaries induced by lensing, such as ISW-lensing., Comment: 22 pages, 9 figures. Abstract abridged

Upcoming cosmic microwave background (CMB) data can be used to explore harmonic 3-point functions that involve the B-mode component of the CMB polarization signal. We focus on bispectra describing the non-Gaussian correlation of the B-mode field and the CMB temperature anisotropies (T) and/or E-mode polarization, i.e. , , and . Such bispectra probe violations of the tensor consistency relation: the model-independent behavior of cosmological correlation functions that involve a large-wavelength tensor mode (gravitational wave). An observed violation of the tensor consistency relation would exclude a large number of inflation models. We describe a generalization of the Komatsu-Spergel-Wandelt (KSW) bispectrum estimator that allows statistical inference on this type of primordial non-Gaussianity with data of the CMB temperature and polarization anisotropies. The generalized estimator shares its statistical properties with the existing KSW estimator and retains the favorable numerical scaling with angular resolution. In this paper we derive the estimator and present a set of Fisher forecasts. We show how the forecasts scale with various experimental parameters such as lower and upper angular band-limit, relevant for e.g. the upcoming ground-based Simons Observatory experiment and proposed LiteBIRD satellite experiment. We comment on possible contaminants due to secondary cosmological and astrophysical sources., Comment: 40 pages, 7 figures

We introduce a publicly available full-sky beam convolution code library intended to inform the design of future cosmic microwave background (CMB) instruments and help current experiments probe potential systematic effects. The code can be used to assess the impact of optical systematics on all stages of data reduction for a realistic experiment, including analyses beyond power spectrum estimation, by generating signal timelines that may serve as input to full analysis pipelines. The design and mathematical framework of the Python code is discussed along with a few simple benchmarking results. We present a simple two-lens refracting telescope design and use it together with the code to simulate a year-long dataset for 400 detectors scanning the sky on a satellite instrument. The simulation results identify a number of sub-leading optical non-idealities and demonstrate significant B-mode residuals caused by extended sidelobes that are sensitive to polarized radiation from the Galaxy. For the proposed design and satellite scanning strategy, we show that a full physical optics beam model generates B-mode systematics that differ significantly from the simpler elliptical Gaussian model. The code is available at https://github.com/adrijd/beamconv