Search

Clusters of galaxies trace the most non-linear peaks in the cosmic density field. The weak gravitational lensing of background galaxies by clusters can allow us to infer their masses. However, galaxies associated with the local environment of the cluster can also be intrinsically aligned due to the local tidal gradient, contaminating any cosmology derived from the lensing signal. We measure this intrinsic alignment in Dark Energy Survey (DES) Year 1 REDMAPPER clusters. We find evidence of a non-zero mean radial alignment of galaxies within clusters between redshifts 0.1–0.7. We find a significant systematic in the measured ellipticities of cluster satellite galaxies that we attribute to the central galaxy flux and other intracluster light. We attempt to correct this signal, and fit a simple model for intrinsic alignment amplitude (AIA) to the measurement, finding AIA = 0.15 ± 0.04, when excluding data near the edge of the cluster. We find a significantly stronger alignment of the central galaxy with the cluster dark matter halo at low redshift and with higher richness and central galaxy absolute magnitude (proxies for cluster mass). This is an important demonstration of the ability of large photometric data sets like DES to provide direct constraints on the intrinsic alignment of galaxies within clusters. These measurements can inform improvements to small-scale modelling and simulation of the intrinsic alignment of galaxies to help improve the separation of the intrinsic alignment signal in weak lensing studies.

The fiducial cosmological analyses of imaging surveys like DES typically probe the Universe at redshifts z < 1. We present the selection and characterization of high-redshift galaxy samples using DES Year 3 data, and the analysis of their galaxy clustering measurements. In particular, we use galaxies that are fainter than those used in the previous DES Year 3 analyses and a Bayesian redshift scheme to define three tomographic bins with mean redshifts around z ∼0.9, 1.2, and 1.5, which extend the redshift coverage of the fiducial DES Year 3 analysis. These samples contain a total of about 9 million galaxies, and their galaxy density is more than 2 times higher than those in the DES Year 3 fiducial case. We characterize the redshift uncertainties of the samples, including the usage of various spectroscopic and high-quality redshift samples, and we develop a machine-learning method to correct for correlations between galaxy density and survey observing conditions. The analysis of galaxy clustering measurements, with a total signal to noise S/N ∼70 after scale cuts, yields robust cosmological constraints on a combination of the fraction of matter in the Universe

ABSTRACT: We present the first systematic follow-up of Planck Sunyaev–Zeldovich effect (SZE) selected candidates down to signal-to-noise (S/N) of 3 over the 5000 deg2 covered by the Dark Energy Survey. Using the MCMF cluster confirmation algorithm, we identify optical counterparts, determine photometric redshifts, and richnesses and assign a parameter, fcont, that reflects the probability that each SZE-optical pairing represents a random superposition of physically unassociated systems rather than a real cluster. The new PSZ-MCMF cluster catalogue consists of 853 MCMF confirmed clusters and has a purity of 90 per cent. We present the properties of subsamples of the PSZ-MCMF catalogue that have purities ranging from 90 per cent to 97.5 per cent, depending on the adopted fcont threshold. Halo mass estimates M500, redshifts, richnesses, and optical centres are presented for all PSZ-MCMF clusters. The PSZ-MCMF catalogue adds 589 previously unknown Planck identified clusters over the DES footprint and provides redshifts for an additional 50 previously published Planck-selected clusters with S/N>4.5. Using the subsample with spectroscopic redshifts, we demonstrate excellent cluster photo-z performance with an RMS scatter in Δz/(1 + z) of 0.47 per cent. Our MCMF based analysis allows us to infer the contamination fraction of the initial S/N>3 Planck-selected candidate list, which is ∼50 per cent. We present a method of estimating the completeness of the PSZ-MCMF cluster sample. In comparison to the previously published Planck cluster catalogues, this new S/N>3 MCMF confirmed cluster catalogue populates the lower mass regime at all redshifts and includes clusters up to z∼1.3.

We cross-match and compare characteristics of galaxy clusters identified in observations from two sky surveys using two completely different techniques. One sample is optically selected from the analysis of 3 years of Dark Energy Survey observations using the redMaPPer cluster detection algorithm. The second is X-ray selected from XMM observations analysed by the XMM Cluster Survey. The samples comprise a total area of 57.4 deg2, bounded by the area of four contiguous XMM survey regions that overlap the DES footprint. We find that the X-ray-selected sample is fully matched with entries in the redMaPPer catalogue, above λ > 20 and within 0.1

The correlation between the broad line region radius and continuum luminosity (R-L relation) of active galactic nuclei (AGNs) is critical for single-epoch mass estimates of supermassive black holes (SMBHs). At z ∼ 1-2, where AGN activity peaks, the R-L relation is constrained by the reverberation mapping (RM) lags of the Mg II line. We present 25 Mg II lags from the Australian Dark Energy Survey RM project based on 6 yr of monitoring. We define quantitative criteria to select good lag measurements and verify their reliability with simulations based on both the damped random walk stochastic model and the rescaled, resampled versions of the observed light curves of local, well-measured AGN. Our sample significantly increases the number of Mg II lags and extends the R-L relation to higher redshifts and luminosities. The relative iron line strength RFe has little impact on the R-L relation. The best-fitting Mg II R-L relation has a slope α = 0.39 ± 0.08 with an intrinsic scatter σrl = 0.15+−000203. The slope is consistent with previous measurements and shallower than the H β R-L relation. The intrinsic scatter of the new R-L relation is substantially smaller than previous studies and comparable to the intrinsic scatter of the H β R-L relation. Our new R-L relation will enable more precise single-epoch mass estimates and SMBH demographic studies at cosmic noon.

We cross-correlate positions of galaxies measured in data from the first three years of the Dark Energy Survey with Compton-y maps generated using data from the South Pole Telescope (SPT) and the Planck mission. We model this cross-correlation measurement together with the galaxy autocorrelation to constrain the distribution of gas in the Universe. We measure the hydrostatic mass bias or, equivalently, the mean halo bias-weighted electron pressure (bh Pe), using large-scale information. We find (bh Pe) to be [0.16+−000403, 0.28+−000504, 0.45+−001006, 0.54+−000708, 0.61+−000608, 0.63+−000807] meV cm−3 at redshifts z ∼ [0.30, 0.46, 0.62, 0.77, 0.89, 0.97]. These values are consistent with previous work where measurements exist in the redshift range. We also constrain the mean gas profile using small-scale information, enabled by the high-resolution of the SPT data. We compare our measurements to different parametrized profiles based on the cosmo-OWLS hydrodynamical simulations. We find that our data are consistent with the simulation that assumes an AGN heating temperature of 108.5 K but are incompatible with the model that assumes an AGN heating temperature of 108.0 K. These comparisons indicate that the data prefer a higher value of electron pressure than the simulations within r500c of the galaxies’ haloes.

Recent cosmological analyses with large-scale structure and weak lensing measurements, usually referred to as 3 × 2pt, had to discard a lot of signal to noise from small scales due to our inability to accurately model non-linearities and baryonic effects. Galaxy–galaxy lensing, or the position–shear correlation between lens and source galaxies, is one of the three two-point correlation functions that are included in such analyses, usually estimated with the mean tangential shear. However, tangential shear measurements at a given angular scale θ or physical scale R carry information from all scales below that, forcing the scale cuts applied in real data to be significantly larger than the scale at which theoretical uncertainties become problematic. Recently, there have been a few independent efforts that aim to mitigate the non-locality of the galaxy–galaxy lensing signal. Here, we perform a comparison of the different methods, including the Y-transformation, the point-mass marginalization methodology, and the annular differential surface density statistic. We do the comparison at the cosmological constraints level in a combined galaxy clustering and galaxy–galaxy lensing analysis. We find that all the estimators yield equivalent cosmological results assuming a simulated Rubin Observatory Legacy Survey of Space and Time (LSST) Year 1 like set-up and also when applied to DES Y3 data. With the LSST Y1 set-up, we find that the mitigation schemes yield ∼1.3 times more constraining S8 results than applying larger scale cuts without using any mitigation scheme.

We compare the two largest galaxy morphology catalogues, which separate early- and late-type galaxies at intermediate redshift. The two catalogues were built by applying supervised deep learning (convolutional neural networks, CNNs) to the Dark Energy Survey data down to a magnitude limit of ∼21 mag. The methodologies used for the construction of the catalogues include differences such as the cutout sizes, the labels used for training, and the input to the CNN - monochromatic images versus gri-band normalized images. In addition, one catalogue is trained using bright galaxies observed with DES (i < 18), while the other is trained with bright galaxies (r < 17.5) and 'emulated' galaxies up to r-band magnitude 22.5. Despite the different approaches, the agreement between the two catalogues is excellent up to i < 19, demonstrating that CNN predictions are reliable for samples at least one magnitude fainter than the training sample limit. It also shows that morphological classifications based on monochromatic images are comparable to those based on gri-band images, at least in the bright regime. At fainter magnitudes, i > 19, the overall agreement is good (∼95 per cent), but is mostly driven by the large spiral fraction in the two catalogues. In contrast, the agreement within the elliptical population is not as good, especially at faint magnitudes. By studying the mismatched cases, we are able to identify lenticular galaxies (at least up to i < 19), which are difficult to distinguish using standard classification approaches. The synergy of both catalogues provides an unique opportunity to select a population of unusual galaxies.

We characterize the properties and evolution of bright central galaxies (BCGs) and the surrounding intracluster light (ICL) in galaxy clusters identified in the Dark Energy Survey and Atacama Cosmology Telescope Survey (DES-ACT) overlapping regions, covering the redshift range 0.20 < z < 0.80. Over this redshift range, we measure no change in the ICL’s stellar content (between 50 and 300 kpc) in clusters with log10(M200m,SZ/M☉) >14.4. We also measure the stellar mass–halo mass (SMHM) relation for the BCG+ICL system and find that the slope, β, which characterizes the dependence of M200m,SZ on the BCG+ICL stellar mass, increases with radius. The outskirts are more strongly correlated with the halo than the core, which supports that the BCG+ICL system follows a two-phase growth, where recent growth (z < 2) occurs beyond the BCG’s core. Additionally, we compare our observed SMHM relation results to the IllustrisTNG300-1 cosmological hydrodynamic simulations and find moderate qualitative agreement in the amount of diffuse light. However, the SMHM relation’s slope is steeper in TNG300-1 and the intrinsic scatter is lower, likely from the absence of projection effects in TNG300-1. Additionally, we find that the ICL exhibits a colour gradient such that the outskirts are bluer than the core. Moreover, for the lower halo mass clusters (log10(M200m,SZ/M☉) < 14.59), we detect a modest change in the colour gradient’s slope with lookback time, which combined with the absence of stellar mass growth may suggest that lower mass clusters have been involved in growth via tidal stripping more recently than their higher mass counterparts.

Recent cosmological analyses rely on the ability to accurately sample from high-dimensional posterior distributions. A variety of algorithms have been applied in the field, but justification of the particular sampler choice and settings is often lacking. Here, we investigate three such samplers to motivate and validate the algorithm and settings used for the Dark Energy Survey (DES) analyses of the first 3 yr (Y3) of data from combined measurements of weak lensing and galaxy clustering. We employ the full DES Year 1 likelihood alongside a much faster approximate likelihood, which enables us to assess the outcomes from each sampler choice and demonstrate the robustness of our full results. We find that the ellipsoidal nested sampling algorithm MULTINEST reports inconsistent estimates of the Bayesian evidence and somewhat narrower parameter credible intervals than the sliced nested sampling implemented in POLYCHORD. We compare the findings from MULTINEST and POLYCHORD with parameter inference from the Metropolis–Hastings algorithm, finding good agreement. We determine that POLYCHORD provides a good balance of speed and robustness for posterior and evidence estimation, and recommend different settings for testing purposes and final chains for analyses with DES Y3 data. Our methodology can readily be reproduced to obtain suitable sampler settings for future surveys.

Reverberation mapping measurements have been used to constrain the relationship between the size of the broad-line region and luminosity of active galactic nuclei (AGN). This R-L relation is used to estimate single-epoch virial black hole masses, and has been proposed to use to standardize AGN to determine cosmological distances. We present reverberation measurements made with Hβ from the 6-yr Australian Dark Energy Survey (OzDES) Reverberation Mapping Program. We successfully recover reverberation lags for eight AGN at 0.12 < z < 0.71, probing higher redshifts than the bulk of Hβ measurements made to date. Our fit to the R-L relation has a slope of α = 0.41 ± 0.03 and an intrinsic scatter of σ = 0.23 ± 0.02 dex. The results from our multi-object spectroscopic survey are consistent with previous measurements made by dedicated source-by-source campaigns, and with the observed dependence on accretion rate. Future surveys, including LSST, TiDES, and SDSS-V, which will be revisiting some of our observed fields, will be able to build on the results of our first-generation multi-object reverberation mapping survey.

ABSTRACT: We present a method for mapping variations between probability distribution functions and apply this method within the context of measuring galaxy redshift distributions from imaging survey data. This method, which we name PITPZ for the probability integral transformations it relies on, uses a difference in curves between distribution functions in an ensemble as a transformation to apply to another distribution function, thus transferring the variation in the ensemble to the latter distribution function. This procedure is broadly applicable to the problem of uncertainty propagation. In the context of redshift distributions, for example, the uncertainty contribution due to certain effects can be studied effectively only in simulations, thus necessitating a transfer of variation measured in simulations to the redshift distributions measured from data. We illustrate the use of PITPZ by using the method to propagate photometric calibration uncertainty to redshift distributions of the Dark Energy Survey Year 3 weak lensing source galaxies. For this test case, we find that PITPZ yields a lensing amplitude uncertainty estimate due to photometric calibration error within 1 per cent of the truth, compared to as much as a 30 per cent underestimate when using traditional methods.

ABSTRACT: Recent analyses have found intriguing correlations between the colour (c) of type Ia supernovae (SNe Ia) and the size of their ‘mass-step’, the relationship between SN Ia host galaxy stellar mass (Mstellar) and SN Ia Hubble residual, and suggest that the cause of this relationship is dust. Using 675 photometrically classified SNe Ia from the Dark Energy Survey 5-yr sample, we study the differences in Hubble residual for a variety of global host galaxy and local environmental properties for SN Ia subsamples split by their colour. We find a 3σ difference in the mass-step when comparing blue (c < 0) and red (c > 0) SNe. We observe the lowest r.m.s. scatter (∼0.14 mag) in the Hubble residual for blue SNe in low mass/blue environments, suggesting that this is the most homogeneous sample for cosmological analyses. By fitting for c-dependent relationships between Hubble residuals and Mstellar, approximating existing dust models, we remove the mass-step from the data and find tentative ∼2σ residual steps in rest-frame galaxy U − R colour. This indicates that dust modelling based on Mstellar may not fully explain the remaining dispersion in SN Ia luminosity. Instead, accounting for a c-dependent relationship between Hubble residuals and global U − R, results in ≤1σ residual steps in Mstellar and local U − R, suggesting that U − R provides different information about the environment of SNe Ia compared to Mstellar, and motivating the inclusion of galaxy U − R colour in SN Ia distance bias correction.

Cosmological analyses with type Ia supernovae (SNe Ia) often assume a single empirical relation between colour and luminosity (β) and do not account for varying host-galaxy dust properties. However, from studies of dust in large samples of galaxies, it is known that dust attenuation can vary significantly. Here, we take advantage of state-of-the-art modelling of galaxy properties to characterize dust parameters (dust attenuation AV, and a parameter describing the dust law slope RV) for 1100 Dark Energy Survey (DES) SN host galaxies. Utilizing optical and infrared data of the hosts alone, we find three key aspects of host dust that impact SN cosmology: (1) there exists a large range (∼1-6) of host RV; (2) high-stellar mass hosts have RV on average ∼0.7 lower than that of low-mass hosts; (3) for a subsample of 81 spectroscopically classified SNe there is a significant (>3σ) correlation between the Hubble diagram residuals of red SNe Ia and the host RV that when corrected for reduces scatter by ∼ 13 per cent and the significance of the 'mass step' to ∼1σ. These represent independent confirmations of recent predictions based on dust that attempted to explain the puzzling 'mass step' and intrinsic scatter (σint) in SN Ia analyses.

We evaluate the consistency between lensing and clustering based on measurements from Baryon Oscillation Spectroscopic Survey combined with galaxy-galaxy lensing from Dark Energy Survey (DES) Year 3, Hyper Suprime-Cam Subaru Strategic Program (HSC) Year 1, and Kilo-Degree Survey (KiDS)-1000. We find good agreement between these lensing data sets. We model the observations using the Dark Emulator and fit the data at two fixed cosmologies: Planck (S8 = 0.83), and a Lensing cosmology (S8 = 0.76). For a joint analysis limited to large scales, we find that both cosmologies provide an acceptable fit to the data. Full utilization of the higher signal-to-noise small-scale measurements is hindered by uncertainty in the impact of baryon feedback and assembly bias, which we account for with a reasoned theoretical error budget. We incorporate a systematic inconsistency parameter for each redshift bin, A, that decouples the lensing and clustering. With a wide range of scales, we find different results for the consistency between the two cosmologies. Limiting the analysis to the bins for which the impact of the lens sample selection is expected to be minimal, for the Lensing cosmology, the measurements are consistent with A = 1; A = 0.91 ± 0.04 (A = 0.97 ± 0.06) using DES+KiDS (HSC). For the Planck case, we find a discrepancy: A = 0.79 ± 0.03 (A = 0.84 ± 0.05) using DES+KiDS (HSC). We demonstrate that a kinematic Sunyaev-Zeldovich-based estimate for baryonic effects alleviates some of the discrepancy in the Planck cosmology. This analysis demonstrates the statistical power of small-scale measurements; however, caution is still warranted given modelling uncertainties and foreground sample selection effects.

We present clustering redshift measurements for Dark Energy Survey (DES) lens sample galaxies used in weak gravitational lensing and galaxy clustering studies. To perform these measurements, we cross-correlate with spectroscopic galaxies from the Baryon Acoustic Oscillation Survey (BOSS) and its extension, eBOSS. We validate our methodology in simulations, including a new technique to calibrate systematic errors that result from the galaxy clustering bias, and we find that our method is generally unbiased in calibrating the mean redshift. We apply our method to the data, and estimate the redshift distribution for 11 different photometrically selected bins. We find general agreement between clustering redshift and photometric redshift estimates, with differences on the inferred mean redshift found to be below |Δz| = 0.01 in most of the bins. We also test a method to calibrate a width parameter for redshift distributions, which we found necessary to use for some of our samples. Our typical uncertainties on the mean redshift ranged from 0.003 to 0.008, while our uncertainties on the width ranged from 4 to 9 per cent. We discuss how these results calibrate the photometric redshift distributions used in companion papers for DES Year 3 results.

In this work, we present the galaxy clustering measurements of the two DES lens galaxy samples: a magnitude-limited sample optimized for the measurement of cosmological parameters, maglim, and a sample of luminous red galaxies selected with the redmagic algorithm. maglim/redmagic sample contains over 10 million/2.5 million galaxies and is divided into six/five photometric redshift bins spanning the range z [0.20, 1.05]/z [0.15, 0.90]. Both samples cover 4143 °2 over which we perform our analysis blind, measuring the angular correlation function with an S/N ∼63 for both samples. In a companion paper, these measurements of galaxy clustering are combined with the correlation functions of cosmic shear and galaxy-galaxy lensing of each sample to place cosmological constraints with a 3 × 2pt analysis. We conduct a thorough study of the mitigation of systematic effects caused by the spatially varying survey properties and we correct the measurements to remove artificial clustering signals. We employ several decontamination methods with different configurations to ensure the robustness of our corrections and to determine the systematic uncertainty that needs to be considered for the final cosmology analyses. We validate our fiducial methodology using lognormal mocks, showing that our decontamination procedure induces biases no greater than 0.5σ in the (ωm, b) plane, where b is the galaxy bias.

Lensing without borders is a cross-survey collaboration created to assess the consistency of galaxy-galaxy lensing signals (Δς) across different data sets and to carry out end-to-end tests of systematic errors. We perform a blind comparison of the amplitude of Δς using lens samples from BOSS and six independent lensing surveys. We find good agreement between empirically estimated and reported systematic errors which agree to better than 2.3σ in four lens bins and three radial ranges. For lenses with zL > 0.43 and considering statistical errors, we detect a 3-4σ correlation between lensing amplitude and survey depth. This correlation could arise from the increasing impact at higher redshift of unrecognized galaxy blends on shear calibration and imperfections in photometric redshift calibration. At zL > 0.54, amplitudes may additionally correlate with foreground stellar density. The amplitude of these trends is within survey-defined systematic error budgets that are designed to include known shear and redshift calibration uncertainty. Using a fully empirical and conservative method, we do not find evidence for large unknown systematics. Systematic errors greater than 15 per cent (25 per cent) ruled out in three lens bins at 68 per cent (95 per cent) confidence at z < 0.54. Differences with respect to predictions based on clustering are observed to be at the 20-30 per cent level. Our results therefore suggest that lensing systematics alone are unlikely to fully explain the 'lensing is low' effect at z < 0.54. This analysis demonstrates the power of cross-survey comparisons and provides a promising path for identifying and reducing systematics in future lensing analyses.

We present the statistical methods that have been developed to analyse the OzDES reverberation mapping sample. To perform this statistical analysis we have created a suite of customizable simulations that mimic the characteristics of each source in the OzDES sample. These characteristics include: the variability in the photometric and spectroscopic light curves, the measurement uncertainties, and the observational cadence. By simulating the sources in the OzDES sample that contain the C iv emission line, we developed a set of criteria that rank the reliability of a recovered time-lag depending on the agreement between different recovery methods, the magnitude of the uncertainties, and the rate at which false positives were found in the simulations. These criteria were applied to simulated light curves and these results used to estimate the quality of the resulting Radius-Luminosity relation. We grade the results using three quality levels (gold, silver, and bronze). The input slope of the R-L relation was recovered within 1σ for each of the three quality samples, with the gold standard having the lowest dispersion with a recovered a R-L relation slope of 0.454 ± 0.016 with an input slope of 0.47. Future work will apply these methods to the entire OzDES sample of 771 AGN.

Galaxy-galaxy lensing is a powerful probe of the connection between galaxies and their host dark matter haloes, which is important both for galaxy evolution and cosmology. We extend the measurement and modelling of the galaxy-galaxy lensing signal in the recent Dark Energy Survey Year 3 cosmology analysis to the highly non-linear scales (100 kpc). This extension enables us to study the galaxy-halo connection via a Halo Occupation Distribution (HOD) framework for the two lens samples used in the cosmology analysis: a luminous red galaxy sample (redmagic) and a magnitude-limited galaxy sample (maglim). We find that redmagic (maglim) galaxies typically live in dark matter haloes of mass log10(Mh/M) ≈ 13.7 which is roughly constant over redshift (13.3-13.5 depending on redshift). We constrain these masses to 15 per cent, approximately 1.5 times improvement over the previous work. We also constrain the linear galaxy bias more than five times better than what is inferred by the cosmological scales only. We find the satellite fraction for redmagic (maglim) to be 0.1-0.2 (0.1-0.3) with no clear trend in redshift. Our constraints on these halo properties are broadly consistent with other available estimates from previous work, large-scale constraints, and simulations. The framework built in this paper will be used for future HOD studies with other galaxy samples and extensions for cosmological analyses.

In this paper, we present and validate the galaxy sample used for the analysis of the baryon acoustic oscillation (BAO) signal in the Dark Energy Survey (DES) Y3 data. The definition is based on a colour and redshift-dependent magnitude cut optimized to select galaxies at redshifts higher than 0.5, while ensuring a high-quality determination. The sample covers ~4100 deg2 to a depth of i = 22.3 (AB) at 10s. It contains 7031 993 galaxies in the redshift range from z = 0.6 to 1.1, with a mean effective redshift of 0.835. Redshifts are estimated with the machine learning algorithm DNF, and are validated using the VIPERS PDR2 sample. We find a mean redshift bias of zbias~0.01 and a mean uncertainty, in units of 1 + z, of σ68~0.03. We evaluate the galaxy population of the sample, showing it is mostly built upon Elliptical to Sbc types. Furthermore, we find a low level of stellar contamination of ≤ 4 per cent. We present the method used to mitigate the effect of spurious clustering coming from observing conditions and other large-scale systematics.We apply it to the BAO sample and calculate weights that are used to get a robust estimate of the galaxy clustering signal. This paper is one of a series dedicated to the analysis of the BAO signal in DES Y3. In the companion papers, we present the galaxy mock catalogues used to calibrate the analysis and the angular diameter distance constraints obtained through the fitting to the BAO scale.

We present measurements of the radial profiles of the mass and galaxy number density around Sunyaev–Zel’dovich (SZ)-selected clusters using both weak lensing and galaxy counts. The clusters are selected from the Atacama Cosmology Telescope Data Release 5 and the galaxies from the Dark Energy Survey Year 3 data set. With signal-to-noise ratio of 62 (45) for galaxy (weak lensing) profiles over scales of about 0.2–20 h-1 Mpc, these are the highest precision measurements for SZ-selected clusters to date. Because SZ selection closely approximates mass selection, these measurements enable several tests of theoretical models of the mass and light distribution around clusters. Our main findings are: (1) The splashback feature is detected at a consistent location in both the mass and galaxy profiles and its location is consistent with predictions of cold dark matter N-body simulations. (2) The full mass profile is also consistent with the simulations. (3) The shapes of the galaxy and lensing profiles are remarkably similar for our sample over the entire range of scales, from well inside the cluster halo to the quasilinear regime. We measure the dependence of the profile shapes on the galaxy sample, redshift, and cluster mass. We extend the Diemer & Kravtsov model for the cluster profiles to the linear regime using perturbation theory and show that it provides a good match to the measured profiles. We also compare the measured profiles to predictions of the standard halo model and simulations that include hydrodynamics. Applications of these results to cluster mass estimation, cosmology, and astrophysics are discussed.