The cosmic microwave background (CMB) energy spectrum is a near-perfect blackbody. The standard model of cosmology predicts small spectral distortions to this form, but no such distortion of the sky-averaged CMB spectrum has yet been measured. We calculate the largest expected distortion, which arises from the inverse Compton scattering of CMB photons off hot, free electrons, known as the thermal Sunyaev-Zel'dovich (TSZ) effect. We show that the predicted signal is roughly one order of magnitude below the current bound from the COBE-FIRAS experiment, but it can be detected at enormous significance (≳1000σ) by the proposed Primordial Inflation Explorer (PIXIE). Although cosmic variance reduces the effective signal-to-noise ratio to 230σ, this measurement will still yield a subpercent constraint on the total thermal energy of electrons in the observable Universe. Furthermore, we show that PIXIE can detect subtle relativistic effects in the sky-averaged TSZ signal at 30σ, which directly probe moments of the optical depth-weighted intracluster medium electron temperature distribution. These effects break the degeneracy between the electron density and the temperature in the mean TSZ signal, allowing a direct inference of the mean baryon density at low redshift. Future spectral distortion probes will thus determine the global thermodynamic properties of ionized gas in the Universe with unprecedented precision. These measurements will impose a fundamental "integral constraint" on models of galaxy formation and the injection of feedback energy over cosmic time.

The kinematic Sunyaev-Zel'dovich (kSZ) signal is a powerful probe of the cosmic baryon distribution. The kSZ signal is proportional to the integrated free electron momentum rather than the electron pressure (which sources the thermal SZ signal). Since velocities should be unbiased on large scales, the kSZ signal is an unbiased tracer of the large-scale electron distribution, and thus can be used to detect the "missing baryons" that evade most observational techniques. While most current methods for kSZ extraction rely on the availability of very accurate redshifts, we revisit a method that allows measurements even in the absence of redshift information for individual objects. It involves cross-correlating the square of an appropriately filtered cosmic microwave background (CMB) temperature map with a projected density map constructed from a sample of large-scale structure tracers. We show that this method will achieve high signal-to-noise when applied to the next generation of high-resolution CMB experiments, provided that component separation is sufficiently effective at removing foreground contamination. Considering statistical errors only, we forecast that this estimator can yield S/N≈3, 120 and over 150 for Planck, Advanced ACTPol, and a hypothetical Stage IV CMB experiment, respectively, in combination with a galaxy catalog from WISE, and about 20% larger S/N for a galaxy catalog from the proposed SPHEREx experiment. We show that the basic estimator receives a contribution due to leakage from CMB lensing, but that this term can be effectively removed by either direct measurement or marginalization, with little effect on the kSZ significance. We discuss possible sources of systematic contamination and propose mitigation strategies for future surveys. We compare the theoretical predictions to numerical simulations and validate the approximations in our analytic approach. This work serves as a companion paper to the first kSZ measurement with this method, where we used CMB temperature maps constructed from Planck and WMAP data, together with galaxies from the WISE survey, to obtain a 3.8-4.5σ detection of the kSZ2 amplitude.

The kinematic Sunyaev-Zel'dovich (KSZ) effect-the Doppler boosting of cosmic microwave background (CMB) photons due to Compton scattering off free electrons with nonzero bulk velocity-probes the abundance and the distribution of baryons in the Universe. All KSZ measurements to date have explicitly required spectroscopic redshifts. Here, we implement a novel estimator for the KSZ-large-scale structure cross-correlation based on projected fields: it does not require redshift estimates for individual objects, allowing KSZ measurements from large-scale imaging surveys. We apply this estimator to cleaned CMB temperature maps constructed from Planck and WMAP data and a galaxy sample from the Wide-field Infrared Survey Explorer (WISE). We measure the KSZ effect at 3.8σ-4.5σ significance, depending on the use of additional WISE galaxy bias constraints. We verify that our measurements are robust to possible dust emission from the WISE galaxies. Assuming the standard Λ cold dark matter cosmology, we directly constrain (f_{b}/0.158)(f_{free}/1.0)=1.48±0.19 (statistical error only) at redshift z≈0.4, where f_{b} is the fraction of matter in baryonic form and f_{free} is the free electron fraction. This is the tightest KSZ-derived constraint reported to date on these parameters. Astronomers have long known that baryons do not trace dark matter on ∼ kiloparsec scales and there has been strong evidence that galaxies are baryon poor. The consistency between the f_{b} value found here and the values inferred from analyses of the primordial CMB and big bang nucleosynthesis verifies that baryons approximately trace the dark matter distribution down to ∼ megaparsec scales. While our projected-field estimator is already competitive with other KSZ approaches when applied to current data sets (because we are able to use the full-sky WISE photometric survey), it will yield enormous signal-to-noise ratios when applied to upcoming high-resolution, multifrequency CMB surveys.

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 per year, depending on the distance (r). 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 predicted aphelion and most likely location of the planet corresponds to the shallower end of these ranges. The search covers the whole 18,000 square degrees of the ACT survey. No significant detections are found, which is used to place limits on the millimeter-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. These bounds approach those of a recent INPOP19a ephemeris-based analysis, but do not exceed it. 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 au < r < 600 au and heliocentric angular velocity , corresponding to low-to-moderate eccentricities. These limits worsen gradually beyond 600 au, reaching 5–15 mJy by 1500 au.

SPHEREx is a proposed SMEX mission selected for Phase A. SPHEREx will carry out the first all-sky spectral survey and provide for every 6.2" pixel a spectra between 0.75 and 4.18 $\mu$m [with R$\sim$41.4] and 4.18 and 5.00 $\mu$m [with R$\sim$135]. The SPHEREx team has proposed three specific science investigations to be carried out with this unique data set: cosmic inflation, interstellar and circumstellar ices, and the extra-galactic background light. It is readily apparent, however, that many other questions in astrophysics and planetary sciences could be addressed with the SPHEREx data. The SPHEREx team convened a community workshop in February 2016, with the intent of enlisting the aid of a larger group of scientists in defining these questions. This paper summarizes the rich and varied menu of investigations that was laid out. It includes studies of the composition of main belt and Trojan/Greek asteroids; mapping the zodiacal light with unprecedented spatial and spectral resolution; identifying and studying very low-metallicity stars; improving stellar parameters in order to better characterize transiting exoplanets; studying aliphatic and aromatic carbon-bearing molecules in the interstellar medium; mapping star formation rates in nearby galaxies; determining the redshift of clusters of galaxies; identifying high redshift quasars over the full sky; and providing a NIR spectrum for most eROSITA X-ray sources. All of these investigations, and others not listed here, can be carried out with the nominal all-sky spectra to be produced by SPHEREx. In addition, the workshop defined enhanced data products and user tools which would facilitate some of these scientific studies. Finally, the workshop noted the high degrees of synergy between SPHEREx and a number of other current or forthcoming programs, including JWST, WFIRST, Euclid, GAIA, K2/Kepler, TESS, eROSITA and LSST.

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σ, 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 will be made available on LAMBDA no later than three months after the journal publication of this article, along with an interactive sky atlas.

We measure the projected number density profiles of galaxies and the splashback feature in clusters selected by the Sunyaev--Zeldovich (SZ) effect from the Advanced Atacama Cosmology Telescope (AdvACT) survey using galaxies observed by the Dark Energy Survey (DES). The splashback radius for the complete galaxy sample is consistent with theoretical measurements from CDM-only simulations, and is located at $2.4^{+0.3}_{-0.4}$ Mpc $h^{-1}$. We split the sample based on galaxy color and find significant differences in the profile shapes. Red galaxies and those in the green valley show a splashback-like minimum in their slope profile consistent with theoretical predictions, while the bluest galaxies show a weak feature that appears at a smaller radius. We develop a mapping of galaxies to subhalos in $N$-body simulations by splitting subhalos based on infall time onto the cluster halos. We find that the location of the steepest slope and differences in the shapes of the profiles can be mapped to differences in the average time of infall of galaxies of different colors. The minima of the slope in the galaxy profiles trace a discontinuity in the phase space of dark matter halos. By relating spatial profiles to infall time for galaxies of different colours, we can use splashback as a clock to understand galaxy quenching. We find that red galaxies have on average been in their clusters for over $3.2 ~\rm Gyrs$, green galaxies about $2.2 ~\rm Gyrs$, while blue galaxies have been accreted most recently and have not reached apocenter. Using the information from the complete radial profiles, we fit a simple quenching model and find that the onset of galaxy quenching in clusters occurs after a delay of about a gigayear, and that galaxies quench rapidly thereafter with an exponential timescale of $0.6$ Gyr.

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 deg2, the deepest 600 deg2 with noise levels below 10µK-arcmin. We use the power spectrum derived from almost 6,000 deg2 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, H0. By combining ACT data with large-scale information from WMAP we measure H0 = 67.6±1.1 km/s/Mpc, at 68% confidence, in excellent agreement with the independently-measured Planck satellite estimate (from ACT alone we find H0 = 67.9 ± 1.5 km/s/Mpc). The Λ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σ; the number of relativistic species, the primordial Helium fraction, and the running of the spectral index are consistent with ΛCDM predictions to within 1.5–2.2σ. 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.

We present the temperature and polarization angular power spectra of the CMB measured by the Atacama Cosmology Telescope (ACT) from 5400 deg2 of the 2013–2016 survey, which covers >15000 deg2 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 ` = 4000. At large angular scales, foreground emission at 150 GHz is ∼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 ΛCDM for the ACT data alone with a prior on the optical depth of τ = 0.065 ± 0.015. ΛCDM is a good fit. The best-fit model has a reduced χ2 of 1.07 (PTE = 0.07) with H0 = 67.9 ± 1.5 km/s/Mpc. We show that the lensing BB signal is consistent with ΛCDM and limit the celestial EB polarization angle to ψP = −0.07◦ ±0.09◦. 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.

The Simons Observatory (SO) is a new cosmic microwave background experiment being built on Cerro Toco in Chile, due to begin observations in the early 2020s. We describe the scientific goals of the experiment, motivate the design, and forecast its performance. SO will measure the temperature and polarization anisotropy of the cosmic microwave background in six frequency bands centered at: 27, 39, 93, 145, 225 and 280 GHz. The initial configuration of SO will have three small-aperture 0.5-m telescopes and one large-aperture 6-m telescope, with a total of 60,000 cryogenic bolometers. Our key science goals are to characterize the primordial perturbations, measure the number of relativistic species and the mass of neutrinos, test for deviations from a cosmological constant, improve our understanding of galaxy evolution, and constrain the duration of reionization. The small aperture telescopes will target the largest angular scales observable from Chile, mapping 10% of the sky to a white noise level of 2 μK-arcmin in combined 93 and 145 GHz bands, to measure the primordial tensor-to-scalar ratio, r, at a target level of σ(r)=0.003. The large aperture telescope will map 40% of the sky at arcminute angular resolution to an expected white noise level of 6 μK-arcmin in combined 93 and 145 GHz bands, overlapping with the majority of the Large Synoptic Survey Telescope sky region and partially with the Dark Energy Spectroscopic Instrument. With up to an order of magnitude lower polarization noise than maps from the Planck satellite, the high-resolution sky maps will constrain cosmological parameters derived from the damping tail, gravitational lensing of the microwave background, the primordial bispectrum, and the thermal and kinematic Sunyaev-Zel'dovich effects, and will aid in delensing the large-angle polarization signal to measure the tensor-to-scalar ratio. The survey will also provide a legacy catalog of 16,000 galaxy clusters and more than 20,000 extragalactic sources.