We perform a joint analysis of the counts and weak lensing signal of redMaPPer clusters selected from the Dark Energy Survey (DES) Year 1 dataset. Our analysis uses the same shear and source photometric redshifts estimates as were used in the DES combined probes analysis. Our analysis results in surprisingly low values for $S_8 =\sigma_8(\Omega_{\rm m}/0.3)^{0.5}= 0.65\pm 0.04$, driven by a low matter density parameter, $\Omega_{\rm m}=0.179^{+0.031}_{-0.038}$, with $\sigma_8-\Omega_{\rm m}$ posteriors in $2.4\sigma$ tension with the DES Y1 3x2pt results, and in $5.6\sigma$ with the Planck CMB analysis. These results include the impact of post-unblinding changes to the analysis, which did not improve the level of consistency with other data sets compared to the results obtained at the unblinding. The fact that multiple cosmological probes (supernovae, baryon acoustic oscillations, cosmic shear, galaxy clustering and CMB anisotropies), and other galaxy cluster analyses all favor significantly higher matter densities suggests the presence of systematic errors in the data or an incomplete modeling of the relevant physics. Cross checks with X-ray and microwave data, as well as independent constraints on the observable--mass relation from SZ selected clusters, suggest that the discrepancy resides in our modeling of the weak lensing signal rather than the cluster abundance. Repeating our analysis using a higher richness threshold ($\lambda \ge 30$) significantly reduces the tension with other probes, and points to one or more richness-dependent effects not captured by our model.

We characterize the luminosity distribution, halo mass dependence, and redshift evolution of red galaxies in galaxy clusters using the SDSS Data Release 8 RedMaPPer cluster sample. We propose a simple prescription for the relationship between the luminosity of both central and satellite galaxies and the mass of their host halos, and show that this model is well-fit by the data. Using a larger galaxy cluster sample than previously employed in the literature, we find that the luminosities of central galaxies scale as $\langle \log L \rangle \propto A_L \log (M_{200b})$, with $A_L=0.39\pm0.04$, and that the scatter of the central--galaxy luminosity at fixed $M_{200b}$ ( $\sigma_{\log L|M}$) is $0.23 ^{+0.05}_{-0.04}$ dex, with the error bar including systematics due to miscentering of the cluster finder, photometry, and photometric redshift estimation. Our data prefers a positive correlation between the luminosity of central galaxies and the observed richness of clusters at a fixed halo mass, with an effective correlation coefficient $d_{\rm{eff}}=0.36^{+0.17}_{-0.16}$. The characteristic luminosity of satellites becomes dimmer from $z=0.3$ to $z=0.1$ by $\sim 20\%$ after accounting for passive evolution. We estimate the fraction of galaxy clusters where the brightest galaxy is not the central to be $P_{\rm{BNC}} \sim 20\%$. We discuss implications of these findings in the context of galaxy evolution and the galaxy--halo connection.

This paper introduces cosmoDC2, a large synthetic galaxy catalog designed to support precision dark energy science with the Large Synoptic Survey Telescope (LSST). CosmoDC2 is the starting point for the second data challenge (DC2) carried out by the LSST Dark Energy Science Collaboration (LSST DESC). The catalog is based on a trillion-particle, 4.225 Gpc^3 box cosmological N-body simulation, the `Outer Rim' run. It covers 440 deg^2 of sky area to a redshift of z=3 and is complete to a magnitude depth of 28 in the r-band. Each galaxy is characterized by a multitude of properties including stellar mass, morphology, spectral energy distributions, broadband filter magnitudes, host halo information and weak lensing shear. The size and complexity of cosmoDC2 requires an efficient catalog generation methodology; our approach is based on a new hybrid technique that combines data-driven empirical approaches with semi-analytic galaxy modeling. A wide range of observation-based validation tests has been implemented to ensure that cosmoDC2 enables the science goals of the planned LSST DESC DC2 analyses. This paper also represents the official release of the cosmoDC2 data set, including an efficient reader that facilitates interaction with the data.

We present the UNIT $N$-body cosmological simulations project, designed to provide precise predictions for nonlinear statistics of the galaxy distribution. We focus on characterizing statistics relevant to emission line and luminous red galaxies in the current and upcoming generation of galaxy surveys. We use a suite of precise particle mesh simulations (FastPM) as well as with full $N$-body calculations with a mass resolution of $\sim 1.2\times10^9\,h^{-1}$M$_{\odot}$ to investigate the recently suggested technique of Angulo & Pontzen 2016 to suppress the variance of cosmological simulations We study redshift space distortions, cosmic voids, higher order statistics from $z=2$ down to $z=0$. We find that both two- and three-point statistics are unbiased. Over the scales of interest for baryon acoustic oscillations and redshift-space distortions, we find that the variance is greatly reduced in the two-point statistics and in the cross correlation between halos and cosmic voids, but is not reduced significantly for the three-point statistics. We demonstrate that the accuracy of the two-point correlation function for a galaxy survey with effective volume of 20 ($h^{-1}$Gpc)$^3$ is improved by about a factor of 40, indicating that two pairs of simulations with a volume of 1 ($h^{-1}$Gpc)$^3$ lead to the equivalent variance of $\sim$150 such simulations. The $N$-body simulations presented here thus provide an effective survey volume of about seven times the effective survey volume of DESI or Euclid. The data from this project, including dark matter fields, halo catalogues, and their clustering statistics, are publicly available at http://www.unitsims.org.

Using data collected by the Dark Energy Survey (DES), we report the detection of intracluster light (ICL) with $\sim300$ galaxy clusters in the redshift range of 0.2-0.3. We design methods to mask detected galaxies and stars in the images and stack the cluster light profiles, while accounting for several systematic effects (sky subtraction, instrumental point-spread function, cluster selection effects and residual light in the ICL raw detection from background and cluster galaxies). The methods allow us to acquire high signal-to-noise measurements of the ICL and central galaxies (CGs), which we separate with radial cuts. The ICL appears as faint and diffuse light extending to at least 1 Mpc from the cluster center, reaching a surface brightness level of 30 mag arcsec$^{-2}$. The ICL and the cluster CG contribute to $44\%\pm17$\% of the total cluster stellar luminosity within 1 Mpc. The ICL color is overall consistent with that of the cluster red sequence galaxies, but displays the trend of becoming bluer with increasing radius. The ICL demonstrates an interesting self-similarity feature -- for clusters in different richness ranges, their ICL radial profiles are similar after scaling with cluster $R_\mathrm{200m}$, and the ICL brightness appears to be a good tracer of the cluster radial mass distribution. These analyses are based on the DES redMaPPer cluster sample identified in the first year of observations.

We present a statistical detection of 1.5 GHz radio continuum emission from a sample of faint z~4 Lyman-break galaxies (LBGs). LBGs are key tracers of the high-redshift star formation history and important sources of UV photons that ionized the intergalactic medium in the early universe. In order to better constrain the extinction and intrinsic star formation rate (SFR) of high-redshift LBGs, we combine the latest ultradeep Karl G. Jansky Very Large Array 1.5 GHz radio image and the Hubble Space Telescope Advance Camera for Surveys (ACS) optical images in the Great Observatories Origins Deep Survey-North. We select a large sample of 1771 z~4 LBGs from the ACS catalogue using $\bband$-dropout color criteria. Our LBG samples have $\iband$~25-28 (AB), ~0-3 magnitudes fainter than M*_UV at z~4. In our stacked radio images, we find the LBGs to be point-like under our 2" angular resolution. We measure their mean 1.5 GHz flux by stacking the measurements on the individual objects. We achieve a statistical detection of $S_{1.5GHz}$=0.210+-0.075 uJy at ~3 sigma, first time on such a faint LBG population at z~4. The measurement takes into account the effects of source size and blending of multiple objects. The detection is visually confirmed by stacking the radio images of the LBGs, and the uncertainty is quantified with Monte Carlo simulations on the radio image. The stacked radio flux corresponds to an intrinsic SFR of 16.0+-5.7 M/yr, which is 2.8X the SFR derived from the rest-frame UV continuum luminosity. This factor of 2.8 is in excellent agreement with the extinction correction derived from the observed UV continuum spectral slope, using the local calibration of meurer99. This result supports the use of the local calibration on high-redshift LBGs for deriving the extinction correction and SFR, and also disfavors a steep reddening curve such as that of the Small Magellanic Cloud.