We use imaging from the first three years of the Dark Energy Survey to characterize the dynamical state of 288 galaxy clusters at $0.1 \lesssim z \lesssim 0.9$ detected in the South Pole Telescope (SPT) Sunyaev-Zeldovich (SZ) effect survey (SPT-SZ). We examine spatial offsets between the position of the brightest cluster galaxy (BCG) and the center of the gas distribution as traced by the SPT-SZ centroid and by the X-ray centroid/peak position from Chandra and XMM data. We show that the radial distribution of offsets provides no evidence that SPT SZ-selected cluster samples include a higher fraction of mergers than X-ray-selected cluster samples. We use the offsets to classify the dynamical state of the clusters, selecting the 43 most disturbed clusters, with half of those at $z \gtrsim 0.5$, a region seldom explored previously. We find that Schechter function fits to the galaxy population in disturbed clusters and relaxed clusters differ at $z>0.55$ but not at lower redshifts. Disturbed clusters at $z>0.55$ have steeper faint-end slopes and brighter characteristic magnitudes. Within the same redshift range, we find that the BCGs in relaxed clusters tend to be brighter than the BCGs in disturbed samples, while in agreement in the lower redshift bin. Possible explanations includes a higher merger rate, and a more efficient dynamical friction at high redshift. The red-sequence population is less affected by the cluster dynamical state than the general galaxy population.

We introduce a galaxy cluster mass observable, $\mu_\star$, based on the stellar masses of cluster members, and we present results for the Dark Energy Survey (DES) Year 1 observations. Stellar masses are computed using a Bayesian Model Averaging method, and are validated for DES data using simulations and COSMOS data. We show that $\mu_\star$ works as a promising mass proxy by comparing our predictions to X-ray measurements. We measure the X-ray temperature-$\mu_\star$ relation for a total of 150 clusters matched between the wide-field DES Year 1 redMaPPer catalogue, and Chandra and XMM archival observations, spanning the redshift range $0.1<z<0.7$. For a scaling relation which is linear in logarithmic space, we find a slope of $\alpha = 0.488\pm0.043$ and a scatter in the X-ray temperature at fixed $\mu_\star$ of $\sigma_{{\rm ln} T_X|\mu_\star}=0.266^{+0.019}_{-0.020}$ for the joint sample. By using the halo mass scaling relations of the X-ray temperature from the Weighing the Giants program, we further derive the $\mu_\star$-conditioned scatter in mass, finding $\sigma_{{\rm ln} M|\mu_\star}=0.26^{+ 0.15}_{- 0.10}$. These results are competitive with well-established cluster mass proxies used for cosmological analyses, showing that $\mu_\star$ can be used as a reliable and physically motivated mass proxy to derive cosmological constraints.

Using $\sim 100$ X-ray selected clusters in the Dark Energy Survey Science Verification data, we constrain the luminosity function (LF) of cluster red sequence galaxies as a function of redshift. This is the first homogeneous optical/X-ray sample large enough to constrain the evolution of the luminosity function simultaneously in redshift ($0.1<z<1.05$) and cluster mass ($13.5 \le \rm{log_{10}}(M_{200crit}) \sim< 15.0$). We pay particular attention to completeness issues and the detection limit of the galaxy sample. We then apply a hierarchical Bayesian model to fit the cluster galaxy LFs via a Schecter function, including its characteristic break ($m^*$) to a faint end power-law slope ($\alpha$). Our method enables us to avoid known issues in similar analyses based on stacking or binning the clusters. We find weak and statistically insignificant ($\sim 1.9 \sigma$) evolution in the faint end slope $\alpha$ versus redshift. We also find no dependence in $\alpha$ or $m^*$ with the X-ray inferred cluster masses. However, the amplitude of the LF as a function of cluster mass is constrained to $\sim 20\%$ precision. As a by-product of our algorithm, we utilize the correlation between the LF and cluster mass to provide an improved estimate of the individual cluster masses as well as the scatter in true mass given the X-ray inferred masses. This technique can be applied to a larger sample of X-ray or optically selected clusters from the Dark Energy Survey, significantly improving the sensitivity of the analysis.

The center determination of a galaxy cluster from an optical cluster finding algorithm can be offset from theoretical prescriptions or $N$-body definitions of its host halo center. These offsets impact the recovered cluster statistics, affecting both richness measurements and the weak lensing shear profile around the clusters. This paper models the centering performance of the \RM~cluster finding algorithm using archival X-ray observations of \RM-selected clusters. Assuming the X-ray emission peaks as the fiducial halo centers, and through analyzing their offsets to the \RM~centers, we find that $\sim 75\pm 8 \%$ of the \RM~clusters are well centered and the mis-centered offset follows a Gamma distribution in normalized, projected distance. These mis-centering offsets cause a systematic underestimation of cluster richness relative to the well-centered clusters, for which we propose a descriptive model. Our results enable the DES Y1 cluster cosmology analysis by characterizing the necessary corrections to both the weak lensing and richness abundance functions of the DES Y1 redMaPPer cluster catalog.

Using archival X-ray observations and a log-normal population model, we estimate constraints on the intrinsic scatter in halo mass at fixed optical richness for a galaxy cluster sample identified in Dark Energy Survey Year-One (DES-Y1) data with the redMaPPer algorithm. We examine the scaling behavior of X-ray temperatures, $T_X$, with optical richness, $\lambda_{RM}$, for clusters in the redshift range $0.2<z<0.7$. X-ray temperatures are obtained from Chandra and XMM observations for 58 and 110 redMaPPer systems, respectively. Despite non-uniform sky coverage, the $T_X$ measurements are $> 50\%$ complete for clusters with $\lambda_{RM} > 130$. Regression analysis on the two samples produces consistent posterior scaling parameters, from which we derive a combined constraint on the residual scatter, $\sigma_{\ln Tx | \lambda} = 0.275 \pm 0.019$. Joined with constraints for $T_X$ scaling with halo mass from the Weighing the Giants program and richness--temperature covariance estimates from the LoCuSS sample, we derive the richness-conditioned scatter in mass, $\sigma_{\ln M | \lambda} = 0.30 \pm 0.04\, _{({\rm stat})} \pm 0.09\, _{({\rm sys})}$, at an optical richness of approximately 70. Uncertainties in external parameters, particularly the slope and variance of the $T_X$--mass relation and the covariance of $T_X$ and $\lambda_{RM}$ at fixed mass, dominate the systematic error. The $95\%$ confidence region from joint sample analysis is relatively broad, $\sigma_{\ln M | \lambda} \in [0.14, \, 0.55]$, or a factor ten in variance.

We constrain the mass--richness scaling relation of redMaPPer galaxy clusters identified in the Dark Energy Survey Year 1 data using weak gravitational lensing. We split clusters into $4\times3$ bins of richness $\lambda$ and redshift $z$ for $\lambda\geq20$ and $0.2 \leq z \leq 0.65$ and measure the mean masses of these bins using their stacked weak lensing signal. By modeling the scaling relation as $\langle M_{\rm 200m}|\lambda,z\rangle = M_0 (\lambda/40)^F ((1+z)/1.35)^G$, we constrain the normalization of the scaling relation at the 5.0 per cent level as $M_0 = [3.081 \pm 0.075 ({\rm stat}) \pm 0.133 ({\rm sys})] \cdot 10^{14}\ {\rm M}_\odot$ at $\lambda=40$ and $z=0.35$. The richness scaling index is constrained to be $F=1.356 \pm 0.051\ ({\rm stat})\pm 0.008\ ({\rm sys})$ and the redshift scaling index $G=-0.30\pm 0.30\ ({\rm stat})\pm 0.06\ ({\rm sys})$. These are the tightest measurements of the normalization and richness scaling index made to date. We use a semi-analytic covariance matrix to characterize the statistical errors in the recovered weak lensing profiles. Our analysis accounts for the following sources of systematic error: shear and photometric redshift errors, cluster miscentering, cluster member dilution of the source sample, systematic uncertainties in the modeling of the halo--mass correlation function, halo triaxiality, and projection effects. We discuss prospects for reducing this systematic error budget, which dominates the uncertainty on $M_0$. Our result is in excellent agreement with, but has significantly smaller uncertainties than, previous measurements in the literature, and augurs well for the power of the DES cluster survey as a tool for precision cosmology and upcoming galaxy surveys such as LSST, Euclid and WFIRST.

We describe updates to the \redmapper{} algorithm, a photometric red-sequence cluster finder specifically designed for large photometric surveys. The updated algorithm is applied to $150\,\mathrm{deg}^2$ of Science Verification (SV) data from the Dark Energy Survey (DES), and to the Sloan Digital Sky Survey (SDSS) DR8 photometric data set. The DES SV catalog is locally volume limited, and contains 786 clusters with richness $\lambda>20$ (roughly equivalent to $M_{\rm{500c}}\gtrsim10^{14}\,h_{70}^{-1}\,M_{\odot}$) and $0.2<z<0.9$. The DR8 catalog consists of 26311 clusters with $0.08<z<0.6$, with a sharply increasing richness threshold as a function of redshift for $z\gtrsim 0.35$. The photometric redshift performance of both catalogs is shown to be excellent, with photometric redshift uncertainties controlled at the $\sigma_z/(1+z)\sim 0.01$ level for $z\lesssim0.7$, rising to $\sim0.02$ at $z\sim0.9$ in DES SV. We make use of \emph{Chandra} and \emph{XMM} X-ray and South Pole Telescope Sunyaev-Zeldovich data to show that the centering performance and mass--richness scatter are consistent with expectations based on prior runs of \redmapper{} on SDSS data. We also show how the \redmapper{} \photoz{} and richness estimates are relatively insensitive to imperfect star/galaxy separation and small-scale star masks.