We describe the Outer Rim cosmological simulation, one of the largest high-resolution N-body simulations performed to date, aimed at promoting science to be carried out with large-scale structure surveys. The simulation covers a volume of (4.225Gpc)^3 and evolves more than one trillion particles. It was executed on Mira, a BlueGene/Q system at the Argonne Leadership Computing Facility. We discuss some of the computational challenges posed by a system like Mira, a many-core supercomputer, and how the simulation code, HACC, has been designed to overcome these challenges. We have carried out a large range of analyses on the simulation data and we report on the results as well as the data products that have been generated. The full data set generated by the simulation totals more than 5PB of data, making data curation and data handling a large challenge in of itself. The simulation results have been used to generate synthetic catalogs for large-scale structure surveys, including DESI and eBOSS, as well as CMB experiments. A detailed catalog for the LSST DESC data challenges has been created as well. We publicly release some of the Outer Rim halo catalogs, downsampled particle information, and lightcone data.

Recent detections of the cross-correlation of the thermal Sunyaev-Zel'dovich (tSZ) effect and weak gravitational lensing (WL) enable unique studies of cluster astrophysics and cosmology. In this work, we present constraints on the amplitude of the non-thermal pressure fraction in galaxy clusters, $\alpha_0$, and the amplitude of the matter power spectrum, $\sigma_8$, using measurements of the tSZ power spectrum from Planck, and the tSZ-WL cross-correlation from Planck and the Red Cluster Sequence Lensing Survey. We fit the data to a semi-analytic model with the covariance matrix using $N$-body simulations. We find that the tSZ power spectrum alone prefers $\sigma_8 \sim 0.85$ and a large fraction of non-thermal pressure ($\alpha_0 \sim 0.2$-$0.3$). The tSZ-WL cross-correlation on the other hand prefers a significantly lower $\sigma_8 \sim 0.6$, and low $\alpha_0 \sim 0.05$. We show that this tension can be mitigated by allowing for a steep slope in the stellar-mass-halo-mass relation, which would cause a reduction of the gas in low-mass halos. In such a model, the combined data prefer $\sigma_8 \sim 0.7$ and $\alpha_0 \sim 0.2$, consistent with predictions from hydrodynamical simulations.

Future data from galaxy redshift surveys, combined with high-resolutions maps of the cosmic microwave background, will enable measurements of the pairwise kinematic Sunyaev-Zel'dovich (kSZ) signal with unprecedented statistical significance. This signal probes the matter-velocity correlation function, scaled by the average optical depth ($\tau$) of the galaxy groups and clusters in the sample, and is thus of fundamental importance for cosmology. However, in order to translate pairwise kSZ measurements into cosmological constraints, external constraints on $\tau$ are necessary. In this work, we present a new model for the intra-cluster medium, which takes into account star-formation, feedback, non-thermal pressure, and gas cooling. Our semi-analytic model is computationally efficient and can reproduce results of recent hydrodynamical simulations of galaxy cluster formation. By calibrating the model using recent X-ray measurements of gas density profiles of clusters and $M_{\mathrm{gas}}-M$ relations of groups and clusters, we show that the integrated $\tau$ is constrained to better than 6% (with a confidence level of 95%), which sets the tightest constraints on the external prior on the optical depth of groups and clusters reported in the literature to date. Our strong prior on $\tau$ should significantly improve cosmological constraints derived from future pairwise kSZ measurements and shed new insights into the nature of dark energy, modified gravity, and neutrino mass.

DESI (Dark Energy Spectropic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. The DESI instrument is a robotically-actuated, fiber-fed spectrograph capable of taking up to 5,000 simultaneous spectra over a wavelength range from 360 nm to 980 nm. The fibers feed ten three-arm spectrographs with resolution $R= \lambda/\Delta\lambda$ between 2000 and 5500, depending on wavelength. The DESI instrument will be used to conduct a five-year survey designed to cover 14,000 deg$^2$. This powerful instrument will be installed at prime focus on the 4-m Mayall telescope in Kitt Peak, Arizona, along with a new optical corrector, which will provide a three-degree diameter field of view. The DESI collaboration will also deliver a spectroscopic pipeline and data management system to reduce and archive all data for eventual public use.

DESI (Dark Energy Spectroscopic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations (BAO) and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. To trace the underlying dark matter distribution, spectroscopic targets will be selected in four classes from imaging data. We will measure luminous red galaxies up to $z=1.0$. To probe the Universe out to even higher redshift, DESI will target bright [O II] emission line galaxies up to $z=1.7$. Quasars will be targeted both as direct tracers of the underlying dark matter distribution and, at higher redshifts ($ 2.1 < z < 3.5$), for the Ly-$\alpha$ forest absorption features in their spectra, which will be used to trace the distribution of neutral hydrogen. When moonlight prevents efficient observations of the faint targets of the baseline survey, DESI will conduct a magnitude-limited Bright Galaxy Survey comprising approximately 10 million galaxies with a median $z\approx 0.2$. In total, more than 30 million galaxy and quasar redshifts will be obtained to measure the BAO feature and determine the matter power spectrum, including redshift space distortions.

We detect the kinematic Sunyaev-Zel'dovich (kSZ) effect with a statistical significance of $4.2 \sigma$ by combining a cluster catalogue derived from the first year data of the Dark Energy Survey (DES) with CMB temperature maps from the South Pole Telescope Sunyaev-Zel'dovich (SPT-SZ) Survey. This measurement is performed with a differential statistic that isolates the pairwise kSZ signal, providing the first detection of the large-scale, pairwise motion of clusters using redshifts derived from photometric data. By fitting the pairwise kSZ signal to a theoretical template we measure the average central optical depth of the cluster sample, $\bar{\tau}_e = (3.75 \pm 0.89)\cdot 10^{-3}$. We compare the extracted signal to realistic simulations and find good agreement with respect to the signal-to-noise, the constraint on $\bar{\tau}_e$, and the corresponding gas fraction. High-precision measurements of the pairwise kSZ signal with future data will be able to place constraints on the baryonic physics of galaxy clusters, and could be used to probe gravity on scales $ \gtrsim 100$ Mpc.

The pairwise kinematic Sunyaev-Zel'dovich (kSZ) signal from galaxy clusters is a probe of their line-of-sight momenta, and thus a potentially valuable source of cosmological information. In addition to the momenta, the amplitude of the measured signal depends on the properties of the intra-cluster gas and observational limitations such as errors in determining cluster centers and redshifts. In this work we simulate the pairwise kSZ signal of clusters at z<1, using the output from a cosmological N-body simulation and including the properties of the intra-cluster gas via a model that can be varied in post-processing. We find that modifications to the gas profile due to star formation and feedback reduce the pairwise kSZ amplitude of clusters by ~50%, relative to the naive 'gas traces mass' assumption. We demonstrate that mis-centering can reduce the overall amplitude of the pairwise kSZ signal by up to 10%, while redshift errors can lead to an almost complete suppression of the signal at small separations. We confirm that a high-significance detection is expected from the combination of data from current-generation, high-resolution CMB experiments, such as the South Pole Telescope, and cluster samples from optical photometric surveys, such as the Dark Energy Survey. Furthermore, we forecast that future experiments such as Advanced ACTPol in conjunction with data from the Dark Energy Spectroscopic Instrument will yield detection significances of at least 20{\sigma}, and up to 57{\sigma} in an optimistic scenario. Our simulated maps are publicly available at: http://www.hep.anl.gov/cosmology/ksz.html

We investigate the hemispherical power asymmetry in the cosmic microwave background on small angular scales. We find an anomalously high asymmetry in the multipole range l=601-2048, with a naive statistical significance of 6.5 sigma. However, we show that this extreme anomaly is simply a coincidence of three other effects, relativistic power modulation, edge effects from the mask applied, and inter-scale correlations. After correcting for all of these effects, the significance level drops to ~1 sigma, i.e., there is no anomalous intrinsic asymmetry in the small angular scales. Using this null result, we derive a constraint on a potential dipolar modulation amplitude, A(k)<0.0045 on the ~10 Mpc-scale, at 95% C.L. This new constraint must be satisfied by any theoretical model attempting to explain the hemispherical asymmetry at large angular scales.

A detection of the stacked integrated Sachs-Wolfe (ISW) signal in the CMB of rare superstructures identified in the SDSS Luminous Red Galaxy catalogue has been reported at very high statistical significance. The magnitude of the observed signal has previously been argued to be more than 3 standard deviations larger than the theoretical \Lambda CDM expectation. However, this calculation was made in the linear approximation, and relied on assumptions that may potentially have caused the \Lambda CDM expectation to be underestimated. Here we update the theoretical model calculation and compare it with an analysis of ISW maps obtained from N-body simulations of a \Lambda CDM universe. The differences between model predictions and the map analyses are found to be small and cannot explain the discrepancy with observation, which remains at >3 s.d. significance. We discuss the cosmological significance of this anomaly and speculate on the potential of alternative models to explain it.

We show how the relativistic matter and velocity power spectra behave in different gauges. We construct a new gauge where both spectra coincide with Newtonian theory on all scales. However, in this gauge there are geometric quantities present which do not exist in Newtonian theory, for example the local variation of the Hubble parameter. Comparing this quantity to second order Newtonian quantities, we find that Newtonian theory is inaccurate on scales larger than 10 Mpc. This stresses the importance of relativistic corrections to Newtonian cosmological N-body simulations on these scales.