Stellar Mass Fractions in z\sim 1 \madcows Clusters

The Massive and Distant Clusters of WISE Survey VI: Stellar Mass Fractions of a Sample of High-Redshift Infrared-selected Clusters

Bandon Decker11affiliation: Department of Physics and Astronomy, University of Missouri, 5110 Rockhill Road, Kansas City, MO 64110, USA , Mark Brodwin11affiliation: Department of Physics and Astronomy, University of Missouri, 5110 Rockhill Road, Kansas City, MO 64110, USA , Zubair Abdulla22affiliation: Kavli Institute for Cosmological Physics, University of Chicago, 5640 South Ellis Avenue, Chicago, IL 60637, USA 33affiliation: Department of Astronomy and Astrophysics, University of Chicago, 5640 South Ellis Avenue, Chicago, IL 60637, USA , Anthony H. Gonzalez44affiliation: Department of Astronomy, University of Florida, 211 Bryant Space Center, Gainesville, FL 32611, USA , Daniel P. Marrone55affiliation: Steward Observatory, University of Arizona, 933 North Cherry Avenue, Tucson, AZ 85721, USA , Christine O’Donnell55affiliation: Steward Observatory, University of Arizona, 933 North Cherry Avenue, Tucson, AZ 85721, USA , S. A. Stanford66affiliation: Department of Physics, University of California, One Shields Avenue, Davis, CA 95616, USA , Dominika Wylezalek77affiliation: European Southern Observatory, Karl Schwarzschild Straße 2, D-85748, Garching bei München, Germany , John E. Carlstrom22affiliation: Kavli Institute for Cosmological Physics, University of Chicago, 5640 South Ellis Avenue, Chicago, IL 60637, USA 33affiliation: Department of Astronomy and Astrophysics, University of Chicago, 5640 South Ellis Avenue, Chicago, IL 60637, USA , Peter R. M. Eisenhardt88affiliation: Jet Propulsion Laboratory, California Institute of Technology, Pasadena, CA 91109, USA , Adam Mantz99affiliation: Kavli Institute for Particle Astrophysics and Cosmology, Stanford University, 382 Via Pueblo Mall, Stanford, CA 94305-4060, USA 1010affiliation: Department of Physics, Stanford University, 382 Via Pueblo Mall, Stanford, CA 94305-4060, USA , Wenli Mo44affiliation: Department of Astronomy, University of Florida, 211 Bryant Space Center, Gainesville, FL 32611, USA , Emily Moravec44affiliation: Department of Astronomy, University of Florida, 211 Bryant Space Center, Gainesville, FL 32611, USA , Daniel Stern88affiliation: Jet Propulsion Laboratory, California Institute of Technology, Pasadena, CA 91109, USA , Greg Aldering1111affiliation: Lawrence Berkeley National Laboratory, 1 Cyclotron Road, MS 50B-4206, Berkeley, CA 94720, USA , Matthew L. N. Ashby1212affiliation: Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, MA 02138, USA , Kyle Boone1111affiliation: Lawrence Berkeley National Laboratory, 1 Cyclotron Road, MS 50B-4206, Berkeley, CA 94720, USA 1313affiliation: Department of Physics, University of California, 366 LeConte Hall MC 7300, Berkeley, CA 94720, USA , Brian Hayden1414affiliation: Space Telescope Science Institute, 3700 San Martin Drive, Baltimore, MD 21218, USA 1111affiliation: Lawrence Berkeley National Laboratory, 1 Cyclotron Road, MS 50B-4206, Berkeley, CA 94720, USA , Nikhel Gupta1515affiliation: School of Physics, University of Melbourne, Grattan Street, Parkville, VIC 3010, Australia , Michael A. McDonald1616affiliation: Kavli Institute for Astrophysics and Space Research, Massachusetts Institute of Technology, 77 Massachusetts Avenue, Cambridge, MA 02139, USA

We present measurements of the stellar mass fractions () for a sample of high-redshift () infrared-selected galaxy clusters from the Massive and Distant Clusters of WISE Survey (\madcows) and compare them to the stellar mass fractions of Sunyaev-Zel’dovich (SZ) effect-selected clusters in a similar mass and redshift range from the South Pole Telescope (SPT)-SZ Survey. We do not find a significant difference in mean  between the two selection methods, though we do find an unexpectedly large range in  for the SZ-selected clusters. In addition, we measure the luminosity function of the \madcows clusters and find , similar to other studies of clusters at or near our redshift range. Finally, we present SZ detections and masses for seven \madcows clusters and new spectroscopic redshifts for five \madcows clusters. One of these new clusters, MOO J1521+0452 at , is the most distant \madcows cluster confirmed to date.

Subject headings:
cosmology: observations — galaxies: clusters: general — galaxies: clusters: intracluster medium

1. Introduction

Galaxy clusters are the largest gravitationally-bound objects in the universe and a thorough knowledge of their composition, history and evolution is important for both cosmological abundance analyses and galaxy formation/evolution studies in the richest environments (e.g., Allen et al., 2011; Kravtsov & Borgani, 2012). It has been found in simulations (e.g., Ettori et al., 2006; Conroy et al., 2007) and suggested observationally (e.g., Lin et al., 2003) that the fraction of a cluster’s total mass that is in stars, , is related to the star formation history of that cluster. It follows that measuring  and , the fraction of mass in the intracluster medium (ICM), in clusters covering a range of masses and ages can constrain the growth and evolutionary history of clusters and the galaxies therein. A proper account of the total stellar mass of a cluster is also a necessary component of calculating the total baryon fraction in a cluster. The cluster baryon fraction is expected to be close to the total baryon fraction of the universe (White et al., 1993), but previous studies (e.g., Gonzalez et al., 2007, 2013; Lin et al., 2012) have found somewhat lower fractions. The size of this discrepancy and its relation to the total mass of the cluster is important cosmologically and can also provide clues to the baryon physics in clusters (He et al., 2005). Coupled with studies that show a cessation of star formation in the cores of large clusters since high redshifts (e.g., Brodwin et al., 2013), such measurements can shed light on the feedback processes involved in the partition of baryons into stars and gas in clusters.

Several studies have previously looked at the stellar mass fraction of clusters and generally find a trend of decreasing  with increasing halo mass. However, with the exception of van der Burg et al. (2014), all these studies were at (Lin et al., 2003; Gonzalez et al., 2007; Andreon, 2010; Zhang et al., 2011; Lin et al., 2012; Gonzalez et al., 2013) and/or used samples that selected clusters entirely on the strength of the signal from the ICM, either from X-ray observations or from the Sunyaev-Zel’dovich (SZ, Sunyaev & Zeldovich, 1970, 1972) decrement (Giodini et al., 2009; Hilton et al., 2013; Chiu et al., 2016, 2018). It is possible, for both SZ- and X-ray-selected samples, that selecting on an observable related to the ICM pressure or X-ray luminosity (approximately ICM density squared) could produce a sample with a bias toward a higher gas mass fraction, presumably at the expense of  (assuming a constant baryon fraction at fixed mass). Such a bias may also prevent the scatter in  from being fairly measured, though the measured scatter in  should be less affected, as the cluster selection does not have any intrinsic bias toward or against stellar mass.

To explore these issues, we use high-redshift infrared-selected clusters from the Massive and Distant Clusters of WISE Survey (\madcows, Gettings et al., 2012; Stanford et al., 2014; Brodwin et al., 2015; Gonzalez et al., 2015; Mo et al., 2018; Gonzalez et al., 2018). \madcows uses the Wide-field Infrared Survey Explorer (WISE, Wright et al., 2010) AllWISE data release (Cutri, 2013) and PanSTARRS (Chambers et al., 2016) optical data to identify overdensities of galaxies at across nearly the entire extragalactic sky. It therefore can provide a greater mass range at high-redshift than SZ surveys because it simultaneously has the area to find the rarest, most massive objects at high redshifts—such as MOO J1142+1527 ( = , ) reported in Gonzalez et al. (2015) and MOO J1521+0452 ( = , ) described herein—and the sensitivity to detect clusters to the same or lower mass limit of current SZ surveys.

In this work we use SZ observations and follow-up Spitzer Space Telescope data on twelve \madcows clusters to calculate  for this high-redshift, infrared-selected sample. We also analyze a comparable sample of SZ-selected clusters from the South Pole Telescope (SPT)-SZ survey (Bleem et al., 2015) using the same methodology and compare these to the same quantities measured for our infrared-selected \madcows clusters. Because the SPT sample is SZ-selected, it fairly measures the average value and scatter in .

The cluster samples and data we use are described in §2 and the analysis thereof is described in §3. We present the results of our  measurements in §4 and discuss them in §5. Our conclusions are in §6. Throughout we use AB magnitudes and a concordance CDM cosmology of , and . We define  as the radius inside which a cluster has an average density 500 times the critical density of the universe and  as the mass inside that radius.

2. Cluster Sample and Data

For our infrared-selected sample, we use twelve \madcows clusters with halo masses calculated from SZ detections from the Combined Array for Research in Millimeter-wave Astronomy (CARMA). SZ observations of four of these (MOO J0319-0025, MOO J1014+0038, MOO J1155+3901 and MOO J1514+1346) are described in Brodwin et al. (2015). A fifth, MOO J1142+1527, the most massive cluster yet found by any method at , is reported in Gonzalez et al. (2015). Here we report new SZ detections for the other seven clusters, along with total masses and radii determined from those data as well as new masses and radii of the previously-reported clusters derived from an updated reduction of the CARMA data, described in §3.1. All twelve clusters have imaging with the Infrared Array Camera (IRAC, Fazio et al., 2004) on Spitzer, which enables us to determine the stellar mass of the clusters as described in §4.1.

The SZ-selected clusters we use for comparison are drawn from the SPT-SZ survey described in Bleem et al. (2015). To ensure we are making a fair comparison between the infrared- and SZ-selected samples, we only use the 33 SPT clusters that lie in a similar range of mass and redshift as the \madcows clusters, specifically and  , and for which comparable IRAC data exist. We do not impose a lower limit on the mass for the SPT subsample as the SPT-SZ catalog has a higher mass threshold than \madcows at these redshifts. A plot of mass versus redshift for both samples is shown in Figure 1.

Figure 1.— Plot showing the distribution of the \madcows clusters (red diamonds) and the comparison SPT clusters (blue circles) in mass and redshift. The open diamonds denote the \madcows clusters in this analysis that currently lack spectroscopic redshifts.

2.1. CARMA Data

Before its closure in early 2015, CARMA was a heterogenous 23-element interferometer with six 10.4 m antennae, nine 6.1 m antennae and eight 3.5 m antennae. All of the antennae were equipped with 30 and 90 GHz receivers and the 10.4 and 6.1 m antennae had additional 230 GHz receivers. CARMA had two correlators, a wide-band (WB) and spectral-line (SL) correlator, and the 3.5 m antennae could operate as an independent array (CARMA-8 mode) or alongside the other 15 antennae in CARMA-23 mode. In its most compact ‘E’ configuration, the shortest CARMA baselines provided an appropriate beam size for SZ observations while the longer baselines enabled point source identification and subtraction.

The CARMA data for the seven new clusters were taken in the summer and autumn of 2014 and the observation dates of all twelve of our \madcows clusters, as well as the on-source observation times excluding observations of the gain and flux calibrators, are given in Table 1. Point source-subtracted SZ maps of the seven clusters newly reported here are shown in Figure 2. The maps are in units of signal-to-noise with negative signal to denote the SZ effect being a decrement at 30 GHz. A 4k taper was applied to the uv data to produce an illustrative beam size and the maps were CLEANed (Högbom, 1974) in a box on a side and centered on the SZ centroid.

Cluster ID RA Dec. UT Dates Exp. TimeaaOn-source, unflagged. S/N
(J2000) (J2000) (hr) () () (Mpc) ()
MOO J00373306 00:37:45.8 33:06:51 2014 Sep 12,27-28 6.0 3.3 1.133
MOO J01051323eeIdentified as a merging cluster from follow-up Chandra imaging (see Gonzalez et al., 2018). 01:05:31.5 13:23:55 2014 Jul 6; Oct 11 7.3 8.1 1.143
MOO J01232545 01:23:50.3 25:45:31 2014 Sep 27 1.9 4.4 1.224
MOO J03190025bbBrodwin et al. (2015), with a mass and radius re-calculated from an improved CARMA reduction pipeline. 03:19:24.4 00:25:21 2013 Sep 30 1.0 5.7 1.194
MOO J10140038bbBrodwin et al. (2015), with a mass and radius re-calculated from an improved CARMA reduction pipeline. 10:14:08.4 00:38:26 2013 Oct 6-7 2.2 8.0 1.231
MOO J11111503 11:11:42.6 15:03:44 2014 Jul 23,25 4.4 5.0 1.32ddPhotometric redshift estimated from IRAC  and  images, with error .
MOO J11421527ccGonzalez et al. (2015), with a mass and radius re-calculated from an improved CARMA reduction pipeline and using a slightly different cosmology. 11:42:45.1 15:27:05 2014 Jul 3 5.0 10.4 1.189
MOO J11553901bbBrodwin et al. (2015), with a mass and radius re-calculated from an improved CARMA reduction pipeline. 11:55:45.6 39:01:15 2012 May 11-12 7.2 2.9 1.009
MOO J12316533 12:31:14.8 65:33:29 2014 Sep 7-8 1.5 4.3 0.99ddPhotometric redshift estimated from IRAC  and  images, with error .
MOO J15141346b,eb,efootnotemark: 15:14:42.7 13:46:31 2013 Jun 1,3,5-7,9,11 8.4 2.8 1.059
MOO J15210452 15:21:04.6 04:52:08 2014 Sep 23 2.5 2.7 1.312
MOO J22060906eeIdentified as a merging cluster from follow-up Chandra imaging (see Gonzalez et al., 2018). 22:06:28.6 09:06:32 2014 Jul 5,8 5.7 3.1 0.926
Table 1Summary of CARMA observations and properties of \madcows clusters
Figure 2.— CARMA maps of the seven new \madcows clusters presented here. Each map is , centered on the centroid of the SZ decrement and in units of signal-to-noise. Emissive point sources have been subtracted out of all the maps and they have all been CLEANed around the decrement. A representative beam pattern is shown in the lower left-hand corner of the map of MOO J1111+1503.

2.2. Spitzer Data

Eight of the \madcows clusters were observed in Spitzer Cycle 9 (Program ID 90177; PI Gonzalez) and have depth in the IRAC  and  channels, while the remaining four were observed to the same depth as part of a Cycle 11-12 snapshot program (PID 11080; PI Gonzalez). This depth allows us to detect objects down to one magnitude fainter than the characteristic magnitude () on all of our clusters with high (%) completeness. The SPT clusters were observed with Spitzer over four Cycles (PID 60099, 70053, 80012, 10101; PI Brodwin) to a depth of in  and in .

2.3. Optical Data

Five of the \madcows clusters have follow-up - and -band imaging with the Gemini Multi-Object Spectrograph (GMOS, Hook et al., 2004) on Gemini-North with five 180 s exposures in the -band and twelve 80 s exposures in the -band from programs GN-2013A-Q-44 and GN-2013B-Q-8 (PI Brodwin). The data were taken between 2013 February and 2015 July.

2.4. New Spectroscopic Redshifts

Five of the \madcows clusters presented here have previously unreported spectroscopic redshifts. We obtained spectroscopic observations of these five clusters from 2015 through 2017, primarily with the Low Resolution Imaging Spectrometer (LRIS, Oke et al., 1995) at the W. M. Keck Observatory, the details of which are given in Table 2. The mask used for each cluster was designed from the Spitzer imaging and focused on the IRAC sequence members identified in a  color-magnitude diagram.

One of the clusters with new spectroscopic redshifts reported here, MOO J1521+0452, is the highest-redshift \madcows cluster with spectroscopy, and with , it is the third-most massive cluster to be found at by any method. The spectroscopy confirmed six cluster members and established as the cluster redshift. Representative spectra of two of the confirmed members are shown in Figure 3.

Figure 3.— Spectra of two of the six confirmed members of MOO J1521+0452 at . The vertical lines show , left-to-right the locations of [O II] , Ca II K & H lines and D4000.

Four cluster members were confirmed for MOO J00373306, establishing the cluster redshift of . MOO J01051323 and MOO J01232545 each had five identified members, placing their redshifts at and , respectively.

In addition to the newly reported clusters above, we also present a new spectroscopic redshift for MOO J1014+0038, previously reported at a photometric redshift of (Brodwin et al., 2015). In addition to LRIS spectroscopy, we also observed this cluster with the Multi-Object Spectrometer For Infra-Red Exploration (MOSFIRE, McLean et al., 2010, 2012) at Keck on 2016 February 01. These new spectra identified seven members and established the redshift for MOO J1014+0038 as . Spectra for two of these members are shown in Figure 4.

Figure 4.— Spectra of confirmed MOO J1014+0038 cluster members from LRIS (top) and MOSFIRE (bottom), establishing a cluster redshift of . Left-to-right, the vertical lines of the top spectrum show the Ca II K&H lines and the 4000 Å break and the vertical lines of the lower plot show the H and [O III] emission features.
RA Dec. Instrument UT Date QualityaaQualities ‘A’ and ‘B’ denote redshifts of high and reasonable certainty, respectively (Stanford et al., 2014). Features
MOO J00373306
00:37:45.77 33:07:50.9 LRIS 2016 August 05 1.131 A D4000
00:37:46.18 33:07:28.2 LRIS 2016 August 05 1.123 A Ca HK
00:37:48.82 33:07:08.4 LRIS 2016 August 05 1.15 B D4000
00:37:47.03 33:06:45.7 LRIS 2016 August 05 1.13 B D4000
MOO J01051323
01:05:26.64 13:23:36.9 LRIS 2015 December 04 1.13 B D4000
01:05:26.20 13:23:53.7 LRIS 2015 December 04 1.14 B D4000
01:05:29.95 13:23:54.6 LRIS 2015 December 04 1.15 A Ca HK,D4000
01:05:35.27 13:23:10.4 LRIS 2015 December 04 1.144 B [O II]3727,D4000
MOO J01232545
01:23:50.95 25:45:47.19 LRIS 2017 July 20 1.20 B D4000
01:23:57.16 25:44:16.67 LRIS 2017 July 20 1.22 B D4000
01:23:47.37 25:46:50.65 LRIS 2017 July 20 1.2214 A [O II]3727
01:23:41.53 25:47:32.78 LRIS 2017 July 20 1.2196 B [O II]3727
MOO J10140038
10:14:07.31 00:38:27.1 LRIS 2015 February 21 1.231 B Ca HK
10:14:10.51 00:37:56.2 LRIS 2015 February 21 1.23 B D4000
10:14:08.11 00:37:36.6 LRIS 2015 February 21 1.239 A Ca HK
10:14:00.32 00:36:43.7 LRIS 2015 February 21 1.22 B [O II]3727
10:14:08.13 00:38:21.3 LRIS 2015 December 06 1.23 B Ca HK,D4000
10:14:12.80 00:38:12.2 MOSFIRE 2016 February 01 1.2318 A H,[O III]4959,5007
10:14:09.71 00:41:11.1 LRIS 2016 March 06 1.226 B [O II]3727
MOO J15210452
15:21:13.66 04:53:28.0 LRIS 2016 July 05 1.308 B Ca HK
15:21:12.10 04:51:16.9 LRIS 2016 July 05 1.317 B Ca HK
15:21:06.79 04:52:09.1 LRIS 2016 July 05 1.312 B Ca HK,D4000
15:21:04.90 04:51:59.8 LRIS 2016 July 05 1.302 B Ca HK,D4000
15:21:04.15 04:52:12.4 LRIS 2016 July 05 1.32 B Ca HK,D4000
15:20:59.35 04:51:40.7 LRIS 2016 July 05 1.314 A Ca HK,D4000
Foreground/Background Objects
00:37:51.56 33:10:07.0 LRIS 2016 August 05 1.453 A [O II]3727,D4000
01:05:22.72 13:23:55.2 LRIS 2015 December 04 0.229 A [O II]3727
01:05:35.14 13:23:36.9 LRIS 2015 December 04 0.248 A [O II]3727,H,H
01:23:48.16 25:46:01.2 LRIS 2017 July 20 0.2120 A H,H
01:23:42.28 25:46:31.4 LRIS 2017 July 20 0.4659 A H,[O III]4959,5007
01:23:42.32 25:47:17.5 LRIS 2017 July 20 0.4364 A H
01:23:56.71 25:46:31.7 LRIS 2017 July 20 1.4781 A [O II]3727,Ca HK
10:14:11.57 00:38:39.3 LRIS 2015 February 21 1.158 A Ca HK,D4000
10:14:02.48 00:34:53.0 LRIS 2015 February 21 0.326 A H
10:14:13.36 00:39:57.8 LRIS 2015 December 06 0.966 A Ca HK,D4000,G
10:14:04.15 00:41:03.5 LRIS 2015 December 06 0.981 A [O II]3727, Ca HK,D4000
10:14:00.76 00:40:23.2 LRIS 2015 December 06 0.3283 A H,[N II],Na D
15:21:08.78 04:52:59.5 LRIS 2016 July 05 0.514 A H,[N II]
15:20:52.34 04:51:32.0 LRIS 2016 July 05 0.489 A H,[N II]
Table 2Spectroscopic Cluster Members

3. Analysis

3.1. Total Cluster Mass

Details of the CARMA observations are given in Table 1. The data, including those for clusters previously reported in Brodwin et al. (2015) and Gonzalez et al. (2015), were re-reduced using a new MATLAB pipeline designed specifically for CARMA data. Mars was used as the flux calibrator for each cluster with the Rudy et al. (1987) flux model and observations of a bright monochromatic quasar were interleaved with the cluster observations for gain calibration. The cluster Comptonization () was calculated by using a Monte Carlo Markov Chain to fit an Arnaud et al. (2010) pressure profile and point source models (where indicated by the long baseline data) to the CARMA data in uv space. The significance of the detection was calculated by comparing for the fit to the Arnaud model and point source(s) to for a fit to just the point source(s) with no cluster model.  and  were calculated from by forcing consistency with the scaling relation from Andersson et al. (2011). The resulting masses, radii and values are shown in Table 1. Updated masses and radii, based on the new pipeline, are reported for the clusters reported in Brodwin et al. (2015) and Gonzalez et al. (2015). These are all consistent within one sigma with the originally reported quantities. The total masses for the SPT-SZ sample are from the Bleem et al. (2015) catalog.

3.2. Catalogs

For each cluster we ran SExtractor (Bertin & Arnouts, 1996) in dual-image mode on the  and  images, selecting on the  image. We used the IRAC coverage maps as weights and SExtractor parameters similar to those in Lacy et al. (2005). These parameters are optimized for IRAC, but we changed DEBLEND_NTHRESH to 64 and DEBLEND_MINCONT to 0.00005 to better deblend sources in the cluster cores. Magnitudes were measured in diameter apertures and corrected to diameter apertures using the corrections from Ashby et al. (2009). Catalogs for the optical images were produced with the same SExtractor parameters, but with MAG_AUTO magnitudes instead of corrected aperture magnitudes. The optical and infrared catalogs were then matched using the IRAC astrometry to produce combined catalogs for each cluster. All of the catalogs have IRAC  and  fluxes that are complete down to magnitudes of 21.0 and 22.5, respectively. The clusters with optical data have additional - and -band data similarly complete to depths of 25.5 and 24.5 magnitudes.

3.3. Cluster Membership

Because our cluster masses are measured at an overdensity of , we only consider galaxies projected within  (as determined from the SZ data) from the centroid of the SZ decrement in our measurement of stellar masses and fractions (e.g., Figure 5). To ensure our choice of center does not significantly impact our results, we also ran our analysis using the centroid of the galaxy distribution and using the BCG as the center. We find no appreciable differences in our results. Within , we also reject objects that likely lie in the foreground by not including any source with an apparent magnitude brighter than at the redshift of our cluster. The effects of this choice of cutoff are discussed in §5.4. The characteristic magnitude was calculated using the same model as was used for our K-corrections (described in §4.1). To limit the effect of incompleteness at the faint-end, we reject objects more than one magnitude fainter than .

Figure 5.— A Spitzer/IRAC  image of MOO J1521+0452 showing an radius circle corresponding to the cluster  of 0.67 Mpc at . Only galaxies inside the red circle were included in the analysis. Cluster members with spectroscopic redshifts are marked with cyan squares.

We used the available optical data for five of the \madcows clusters to identify stars in color-color space. Following Eisenhardt et al. (2004), we plot versus  colors for each of our possible cluster members. To the limit where our optical data are complete for all clusters, we characterize as stars objects falling above the line shown in Figure 6 that separates objects with the colors of stars from objects that are likely galaxies. Only objects bright enough to be clearly detected in even the shallowest of our optical images are so characterized to ensure a consistent cut across all clusters. We also match our catalogs to objects in the Gaia DR2 catalog (Gaia Collaboration et al., 2016, 2018) with greater than parallax, to confirm that objects known to be stars are the objects being rejected by this approach. We cannot do this for the SPT clusters due to a lack of comparable optical data.

Figure 6.— Combined color-color plot of all the \madcows clusters for which there are GMOS data, showing color plotted against . Objects above the blue dashed line have colors consistent with being stars. Objects with parallax in the Gaia DR2 catalog are plotted as red stars.

Although the bulk of galaxies within  are cluster members, there is still a line-of-sight interloper contribution that must be subtracted. To account for this, we determine the expected contribution to the total flux density from field galaxies within the projected  area and subtract it off the flux density calculated from our cluster. We use the Spitzer Deep Wide-Field Survey (SDWFS, Ashby et al., 2009) to do this, applying the same brightness cuts to reject non-cluster members as we apply to our cluster catalogs. For the clusters with optical data allowing the rejection of stars, we use optical photometry from the NOAO Deep Wide-Field Survey (NDWFS, Jannuzi & Dey, 1999) to make the analogous stellar rejection in our background. For each cluster, we treat all remaining objects in the SDWFS catalog as though they were at the redshift of that cluster and calculate how much spurious luminosity they would add. We use the SDWFS field to determine our background because the IRAC imaging is deeper than that of our clusters and because SDWFS is large enough to smooth out small-scale variations in the background level. This background selection methodology does produce an appreciable systematic uncertainty in our results, as discussed in §5.4.

3.4. Completeness

To correct for incompleteness in our IRAC catalogs, we ran completeness simulations over the range of magnitudes at which we were looking using IRAF’s mkobjects task in the noao/artdata package. For each cluster we added ten random point sources in each half magnitude bin to the IRAC  image, and ran SExtractor to see how many were recovered. We repeated this process 1,000 times in each magnitude bin. This was done for both the \madcows and SPT clusters and we performed a similar analysis on the SDWFS  image and on the optical images of the clusters. The average completeness curve for \madcows and SPT are shown in Figure 7. At , the faint-end limit of our analysis, the catalogs of both surveys are approximately 70% complete, depending slightly on cluster redshift. Because our clusters have slightly different (depending on redshift), our faint-end cutoff varies slightly, as shown in the figure.

Figure 7.— Average completeness curves for the \madcows (red) and SPT (blue) IRAC  images in AB magnitudes. The shaded region represents the range of maximum depths to which our analysis extends.

4. Results

4.1. Stellar Mass

We calculate the stellar mass of the galaxies selected as possible cluster members using their rest-frame -band luminosity. The rest frame -band is centered at the peak of the emission from the old, red stars that dominate the stellar mass of the galaxy. It is therefore a relatively low-scatter proxy for total stellar mass (e.g., Hainline et al., 2011) with a relatively small dependence on the overall SED. At this is easily probed by the IRAC  band. To determine the K-correction from observed IRAC  to rest-frame -band, we use EZGal (Mancone & Gonzalez, 2012). We construct a synthetic galaxy SED with a Bruzual & Charlot (2003) tau model, formation redshift , solar metallicity and a Chabrier (2003) IMF. From this SED we derive a K-correction to the absolute magnitude in the -band, from which we calculate . We statistically correct our luminosities for incompleteness using the simulations described above. We use the same EZGal model to determine the stellar mass-to-light ratio in the -band at the cluster redshift. This ratio is different for each cluster, depending on the redshift, but is close to 0.34 on average. We apply the stellar ratio to the sum of the luminosities of all the objects along the line of sight minus the background contribution estimated from SDWFS to get our final cluster stellar mass. Both the \madcows and SPT clusters were analyzed in the same way and to the same depth to allow for direct comparison of the two samples.

4.2. Estimating Stellar Corrections with Luminosity Functions

Before calculating total stellar mass fractions, we need to account for foreground stars along the line-of-sight to our clusters. We do this by combining the optical stellar identification discussed above with cluster luminosity functions to estimate and correct the total impact from stars on our clusters that lack optical data.

The mean IRAC  luminosity function (LF) for the five \madcows clusters with optical data for stellar rejection is shown in Figure 8. To make this LF we applied the membership cuts from §3.3, including stellar rejection from the optical data, to each cluster and evolution-corrected the members to . The galaxies from all the clusters were then binned in quarter-magnitude wide bins and the appropriate completeness and statistical background corrections were applied. The uncertainties on the values are Poisson errors.

We fit to the data a parameterized Schechter function of the form

(Schechter, 1976) and we fix as our data are not deep enough to constrain the faint-end slope. This choice is consistent with Mancone et al. (2010) and is a reasonable value for our data. The best-fit value is and we take this LF as representative of clusters independent of selection. The error on the fit is calculated from the range of , and is the same as the error calculated from bootstrap resampling. This value for is slightly lower than, but close to that of Muzzin et al. (2008) who found (in AB magnitudes) for IRAC  at and Mancone et al. (2010) who found at . It is also consistent with the value of found for infrared-selected clusters in a higher redhift bin () by Wylezalek et al. (2014).

Figure 8.— The IRAC  luminosity function for the five \madcows clusters with optical data for stellar rejection. The solid circles are background subtracted number per magnitude in each bin and the error bars are from Poisson noise. The dashed line is a best-fit Schechter function with a fixed and the best-fit value of is shown.

We used a similar approach to make luminosity functions for the full sample of twelve \madcows clusters and for the SPT clusters, shown in Figures 9 and 10, respectively. The stellar contamination is more extensive for the SPT clusters because the sample extends to a lower galactic latitude, where there is more line-of-sight contamination, than does the \madcows sample. We do not have adequate optical data for all of these clusters and thus do not attempt stellar corrections on a per-cluster basis. Rather, we construct a statistical stellar correction as follows. We fit the Schechter function determined above, allowing only to vary (i.e., with fixed and as for the clusters without stellar contamination), to the points at the faint end of \madcows and SPT LFs that show no evidence of stellar contamination (as determined by the SPT LF). These are the points plotted with filled circles in Figures 9 and 10; the unfit portion of the LFs, where there appears to be significant stellar contamination in the SPT LF, is plotted with black crosses. The ratios between the areas under these ‘no-stars’ fits for each sample, over the full magnitude range in this work, to the area under their respective observed LFs is the statistical stellar correction factor for that sample. We multiply the measured luminosity of each cluster by the correction factor of the sample to get the true luminosity for that cluster absent stellar contamination.

Figure 9.— The mean IRAC  luminosity function for the full sample of twelve \madcows clusters with no optical rejection of stars. All of the points are background-subtracted number per magnitude in each bin and the error bars are from Poisson noise. The black crosses on the bright end are points with potential stellar contamination that we did not include when fitting the Schechter function, which is represented by the dashed red line. For the Schechter function, we fixed and to match the luminosity function derived using optical stellar rejection.
Figure 10.— The mean IRAC  luminosity function for the 33 comparison SPT clusters in this work. All of the points are background-subtracted number per magnitude in each bin and the error bars are from Poisson noise. The black crosses on the bright end are points with likely stellar contamination that we did not include when fitting the Schechter function, which is represented by the dashed blue line. For the Schechter function, we fixed and to match the \madcows LF. Note the stellar contamination at the bright end.

4.3. Stellar Mass Fraction

Figure 11.— Upper: Stellar mass fraction versus total mass for the \madcows (red) and SPT (blue) clusters. The size of the systematic error in , which varies with , is represented by the black error bars on either end of the figure. The green dashed line is the low-redshift relation from G13. Lower: The stellar mass fractions of each cluster normalized by the G13 relation versus total mass. The error about the mean normalized  for both samples is calculated from bootstrap resampling and for \madcows (SPT) is plotted in pink (cyan) across the figure. The scatter in normalized  is shown by the thick, vertical red (blue) error bars.

To calculate , we divide the stellar mass of the cluster by the total mass calculated from the SZ decrement described above. The stellar mass that we use is calculated by summing the completeness- and K-corrected -band luminosity of every object projected within  of the cluster SZ centroid and subtracting the average background calculated from SDWFS. We then multiply this luminosity by the M/L ratio from our EZGal model for the cluster redshift and the average stellar correction for either the \madcows or SPT subsample calculated above. The systematic uncertainties inherent in this method are discussed in §5.4.

A plot of  versus  is shown in the upper panel of Figure 11, in which the red diamonds represent the infrared-selected \madcows clusters and the blue circles represent the SZ-selected SPT clusters. The dashed green line is the low-redshift relation found by Gonzalez et al. (2013, hereafter G13) and the black error bars on either side of the plot indicate the systematic error introduced by the background subtraction. For each cluster in both samples the stellar mass fraction was calculated without any stellar rejection and then the average stellar correction for the appropriate sample, as described in §4.2, was applied in order to achieve a consistent stellar correction for all the clusters in each sample.

On average, the \madcows clusters do not have significantly higher stellar mass fractions than the SPT clusters. There is a sizable systematic error, largely from the background subtraction, which is both larger than the statistical error and mass dependent, but it should affect both samples to the same degree and thus does not affect the direct comparison. This is dicussed further in §5.4. To ensure that this comparison of is unrelated to the trend of  with mass seen at low redshift, we also divide out the G13 trend line, as shown in the lower half of Figure 11. The errors on the resulting G13-normalized means for each sample are calculated from bootstrap resampling and shown as horizontal pink and cyan bars across the data. This normalization still does not show a significant difference between the mean of the twelve \madcows clusters and the 33 SPT clusters, though there is still a relatively large error on the individual  errors for both sets of clusters. Stellar masses and stellar mass fractions for the \madcows clusters are given in Table 3.

MOO J00373306 1.139
MOO J01051323 1.143
MOO J01232545 1.224
MOO J03190025 1.194
MOO J10140038 1.229
MOO J11111503 1.32
MOO J11421527 1.189
MOO J11553901 1.009
MOO J12316533 0.99
MOO J15141346 1.059
MOO J15210452 1.312
MOO J22060906 0.951
Table 3\madcows Stellar Mass Fractions

As the vertical red and blue error bars in the lower panel of Figure 11 show, the scatter in the SPT stellar mass fractions is larger than that of the \madcows clusters. There is also a much larger range in the SPT stellar mass fractions, with an order of magnitude separating the highest  clusters from the lowest. The scatter in  seen in the \madcows clusters is lower, but may not be representative of the general cluster population because of two selection biases. First, \madcows is a stellar mass-selected cluster sample. As such, it may be biased toward systems with higher-than-average  values. Second, this particular subset of \madcows clusters consists of the most significant detections from the first stage of the study, so may not be representative of the sample or of clusters as a whole. We do not expect the different redshift distributions to introduce a bias, however, as we find no evidence that  evolves with redshift. The SPT clusters, however, should provide a fair sample of the mean value and scatter of the stellar mass fraction at the redshift of those SZ-selected clusters because they are selected independently of those components. We compared the stellar mass fractions of the \madcows and SPT samples using a Kolmogorov-Smirnov test and found they were consistent with being drawn from the same underlying distribution.

The \madcows sample contains three clusters known to be merging from high-resolution Chandra X-ray Observatory follow-up observations (Gonzalez et al., 2018). Previous studies of the effect merging has on the inferred mass of a cluster have produced mixed conclusions, with some (e.g., Poole et al., 2007; Krause et al., 2012) finding that major mergers bias the inferred mass of a system low for most of the observed timescale and others, (e.g., Marrone et al., 2012) finding the mass of merging clusters was overestimated. We do not expect merging to affect the observed richness of a cluster in the same way as the mass, however, so any effect on the inferred mass will bias our measurement of . We do not have X-ray data for the full \madcows sample or the comparison SPT sample, so we cannot fully remove mergers from our current analysis. However the effect of excluding these clusters, for which we know our  measurement is likely to be wrong, is shown in Figure 12. The clusters are plotted in the same manner as the lower part of Figure 11, however the three clusters known to be mergers are now plotted as open red diamonds and the mean is recalculated to exclude them. Although they are not large outliers, the three merging systems do have the highest normalized stellar mass fractions of the \madcows sample. When they are excluded, the mean-normalized  for \madcows decreases to , still higher than that of the SPT clusters, but now consistent within 1 . We also removed two clusters from the SPT sample identified as mergers in Nurgaliev et al. (2017, shown as open circles) which did not affect the mean of the SPT clusters.

Figure 12.— The same as the lower panel of Figure 11, with the merging \madcows clusters (now shown as open diamonds) removed from the calculation of the mean normalized . The effect of removing these clusters for which the total mass is known to be underestimated relative to the stellar mass is to drop the G13-normalized mean to , higher than the unchanged SPT mean.

5. Discussion

5.1. Comparison of Stellar Mass Fractions

As discussed above, Figures 11 and 12 show that the average stellar mass fraction in the \madcows sample is not significantly higher than that of the SPT sample, though there is considerable scatter. To confirm that this is not an artifact of the trend of  with mass we normalized all the  measurements relative to the G13 relation and measured the normalized mean  for both samples, shown in the lower panel of Figure 11. While the mean normalized  for \madcows, , is higher than the corresponding mean for the SPT sample, , these are consistent within .

5.2. Scatter in the Stellar Mass Fraction

The SZ-selected SPT clusters are best-suited to measure the scatter in  at high-redshift as they are selected independently of stellar content and thus should represent an unbiased sampling of the stellar mass fraction in the full cluster population. The large range in  seen in this sample, approximately an order of magnitude (see Figure 11), is perhaps surprising. As Figure 13 shows, however, this variation is clearly apparent in a visual inspection of the richnesses of two clusters with the same halo mass. Although both clusters in this figure have an SZ mass of   (Bleem et al., 2015), SPT-CL J0154-4824 (left) has a stellar mass fraction of   whereas SPT-CL J2148-4843 (right) has a stellar mass fraction of  , an order of magnitude higher.

Figure 13.— IRAC  images of SPT-CL J0154-4824 (left) and SPT-CL J2148-4843 (right) showing the large difference in richness between clusters of the same halo mass. The projected  of each cluster is shown as a red circle. The difference in the angular size of the two circles is due to the redshift difference, which boosts the richness by 28% in the nearer cluster, but the comparison is relatively unaffected by the differential K-correction between the clusters as in Spitzer  is not significantly different between the two redshifts.

The \madcows clusters in this work do not exhibit the same wide peak-to-trough range of stellar mass fractions nor as large a scatter, measured by the standard deviation of , presumably because they represent the high-richness end of an infrared-selected sample rather than a fair cross-section of all clusters. We attempt to quantify the intrinsic scatter in  of both samples about their respective means, independent of our measurement errors, by assuming that the reduced chi-squared will be equal to unity when all the errors are included in the error budget. We therefore set the reduced chi-squared for each sample to unity and solve for the intrinsic scatter term. We find a significant intrinsic scatter, dex for the SPT and dex for \madcows. This discrepancy supports the idea described in §4.3, that the \madcows clusters may not provide a fair measurement of the scatter in  due to their selection, but the SPT clusters should. By the same token, the \madcows clusters should provide fair measurements of the scatter in  that the SZ-selected surveys may not; this is a topic for future analyses with \madcows. The SPT clusters show a larger intrinsic scatter in  than is predicted in simulations, such as those of Kravtsov et al. (2005), Ettori et al. (2006) and Planelles et al. (2013). Very recently, IllustrisTNG (Pillepich et al., 2017) directly measured the scatter in the stellar-total mass relationship in simulated clusters at and and found a very low scatter in the relationship, only 0.07 dex. Some of the low values and high scatter in the SPT  measurements may be due to the masses of low signal-to-noise clusters being overestimated. The clusters we use go to the low signal-to-noise limit of the SPT-SZ catalog and it is possible that some of these are lower mass clusters that scattered up above the cutoff. If we exclude these clusters, the intrinsic scatter of the SPT sample becomes consistent with that of the \madcows clusters. This effect notwithstanding, understanding the baryonic processes causing the remaining large intrinsic scatter in stellar mass fraction, for which the \madcows measurement of dex may be considered a lower limit, is a challenge for the next generation of cosmological simulations.

5.3. Comparison to Other Works

Given the systematic uncertainties described above, it is difficult to make direct comparisons to other works with different systematics. Nevertheless, other works with similar methodologies provide good external checks on our results, and in particular, allow us to test the effect of infrared- versus ICM-selection.

Chiu et al. (2018) also measured  for 84 clusters from the SPT-SZ survey, some of which overlap with our SPT comparison clusters. We do not expect to find the same  values for these clusters, as they use a slightly different cluster mass estimation (from de Haan et al., 2016) and an SED-fitting method to calculate stellar mass. Nevertheless, their average value for  is consistent with ours for the clusters in the same range of mass and redshift.

Hilton et al. (2013) reported stellar and total masses for a sample of 14 SZ-selected clusters from the Atacama Cosmology Telescope (ACT) in a redshift range of . They have a mean stellar mass fraction of  = , which is larger than what we find for our SZ-selected clusters. However, we use a Chabrier (2003) IMF to calculate stellar mass-to-light ratios which results in lower stellar masses than the Salpeter (1955) IMF Hilton et al. (2013) used. Accounting for the difference in stellar mass resulting from the choice of IMFs (0.24 dex), our results are consistent with theirs.

Similarly, van der Burg et al. (2014) reported stellar and halo masses for ten red sequence-selected clusters in a similar redshift range as ours. Using SED-fitting to determine the stellar mass of each galaxy, they find a mean stellar mass fraction for their IR-selected clusters of  . This is consistent with our \madcows mean of   , however their method of calculating stellar mass has different systematics to ours. Correcting for these, as described below, shifts their average stellar mass fraction higher than the \madcows value, but it remains consistent with the G13 trend due to their lower mass range. When we divide out the G13 line in the same manner as in Figure 11, we find they have an average normalized stellar mass fraction of /G13 , consistent with what we find for \madcows.

Figure 14 shows  versus  for our \madcows and SPT clusters plotted alongside the values found by the studies described above. To make a meaningful comparison, we corrected the Hilton et al. (2013) and van der Burg et al. (2014) results to a Chabrier IMF. We further corrected the latter for the offset between SED-fitted and -based stellar masses reported in that work. The infrared-selected \madcows and van der Burg et al. (2014) clusters are plotted as red and violet diamonds, and the SZ-selected SPT clusters in this work, the Chiu et al. (2018) SPT clusters and the Hilton et al. (2013) ACT clusters are plotted as blue, green and cyan circles, respectively. The SZ-selected studies again find broadly similar stellar mass fractions to the infrared-selected studies, consistent with what we find here. The G13 relation is plotted as a dashed line and for each sample error bars are plotted for three representative clusters.

Figure 14.— Comparison of the  measured in this work (red diamonds, blue circles) to  measured by Chiu et al. (2016, green circles), Hilton et al. (2013, cyan circles) and van der Burg et al. (2014, violet diamonds). All samples have been adjusted to be consistent with our methodology. Error bars are plotted for three representative clusters for each sample.

5.4. Systematics

There are three main sources of systematic error in our analysis. The largest is due to our background subtraction; this error is represented by the black error bars in Figure 11. We quantify the size of this uncertainty by measuring the background luminosity from the SDWFS field in radius cutouts across the field and measure the scatter in this background to estimate small-scale variation due to clustering. We add this scatter in quadrature with the field-to-field scatter derived by comparing SDWFS to similar measurements in the EGS (Davis et al., 2007) and COSMOS fields (Scoville et al., 2007). Since this is an error in the luminosity—and therefore the stellar mass—of each cluster, the size of the systematic error in  decreases with increasing . This systematic error is a uniform shift affecting both the \madcows and SPT clusters equally, so it does not affect our comparison of the infrared and ICM selection methods.

The second source of systematic uncertainty in the absolute value of  for our clusters is our choice of stellar mass-to-light ratio. There are two components to this systematic. The first is the choice of tau model described in §4.1, but this is a small effect. The 1.6 bump is largely insensitive to the star formation history of the galaxy, so varying tau does not have a large effect on the M/L ratio. The second component is the choice of IMF. We use a Chabrier (2003) IMF, but other choices, such as the Salpeter (1955) IMF, are also common. This has a large effect on our M/L ratio, almost doubling it for a 1 Gyr tau model. However, since this is easily corrected for and does not affect any comparisons we make, we do not include it in our systematic error bar in Figure 11.

A final possible source of systematic uncertainty stems from our rejection of cluster non-members using magnitude cuts. Our choice of as a brightness threshold strikes a balance between maximizing the bright members included and minimizing the inclusion of bright foreground interlopers. Although this choice is a somewhat arbitrary threshold, changing it has only a small effect on our values for  since we already statistically correct for non-member contamination, and one that is quite consistent from cluster-to-cluster. It does not make an appreciable difference to our analysis.

Our faint-end cutoff leads to a modest underestimate of the total stellar mass. Integrating a luminosity function with beyond suggests we could be missing of the stellar mass from fainter galaxies. If we correct our stellar masses for this, the result is a simple multiplicative increase of all our  values, but by an amount less than both the scatter and the existing systematic error. Since this offset affects all clusters equally, it does not affect the scatter in either sample, or our comparison between the \madcows and SPT stellar mass fractions. As a practical matter, the large uncertainties in and make it difficult to accurately quantify the size of this uncertainty, and thus we choose not to include it in our analysis.

6. Conclusions

We have measured the stellar mass fractions of twelve infrared-selected clusters from \madcows and 33 SZ-selected clusters from the SPT-SZ survey and found little difference in average  between the two selection methods. We measured  using IRAC  images of the clusters at as a proxy for stellar mass along with total masses derived from SZ measurements. We found that when accounting for mergers in the \madcows sample and normalizing over the trend of stellar mass fraction with total mass, the infrared-selected \madcows clusters have an average stellar mass fraction of , higher than the average stellar mass fraction of  for the SPT, but not significantly so.

We also compare our results to those of Hilton et al. (2013), van der Burg et al. (2014) and Chiu et al. (2016) who also looked at stellar mass fractions in cluster samples of comparable mass and redshift to ours. When we correct for the differences between our methodologies and those of the other studies, we find our results are consistent with all three and they support our result that infrared-selected clusters do not have an appreciably higher mean than SZ-selected clusters. We also compare the value we calculate for of the IRAC  luminosity function to that found by Muzzin et al. (2008), Mancone et al. (2010) and Wylezalek et al. (2014) and find similar results.

We found an unexpectedly large range in the stellar mass fractions of individual clusters in the SPT sample and a larger range and scatter in  than in our \madcows clusters. It is possible that the SZ-selected SPT clusters give a fairer sample of the full range of  than the infrared-selected \madcows clusters do. Future work with \madcows will compare  measurements in infrared- and SZ-selected cluster samples to look for a comparable selection effect in the latter.

Finally, we have presented SZ observations of seven new \madcows clusters and new spectroscopic redshifts for five clusters. Among the SZ observations of the seven new \madcows clusters is MOO J1521+0452, which at is one of the most massive clusters yet found at . Along with the previous discovery of a cluster of   at , reported in Gonzalez et al. (2015), this further demonstrates the ability of \madcows’ nearly all-sky infrared selection to find the most massive clusters at high redshifts.

Support for CARMA construction was derived from the Gordon and Betty Moore Foundation, the Kenneth T. and Eileen L. Norris Foundation, the James S. McDonnell Foundation, the Associates of the California Institute of Technology, the University of Chicago, the states of California, Illinois, and Maryland, and the National Science Foundation. CARMA development and operations are supported by the National Science Foundation under a cooperative agreement, and by the CARMA partner universities; the work at Chicago was supported by NSF grant AST- 1140019. Additional support was provided by PHY- 0114422. This work is based in part on observations made with the Spitzer Space Telescope, which is operated by the Jet Propulsion Laboratory, California Institute of Technology under a contract with NASA. Support for this work was provided by NASA through an award issued by JPL/Caltech. This publication makes use of data products from the Wide-field Infrared Survey Explorer, which is a joint project of the University of California, Los Angeles, and the Jet Propulsion Laboratory/California Institute of Technology, funded by the National Aeronautics and Space Administration. Some of the data presented herein were obtained at the W.M. Keck Observatory, which is operated as a scientific partnership among the California Institute of Technology, the University of California and the National Aeronautics and Space Administration. The Observatory was made possible by the generous financial support of the W.M. Keck Foundation. This work was based in part on observations obtained at the Gemini Observatory, which is operated by the Association of Universities for Research in Astronomy, Inc., under a cooperative agreement with the NSF on behalf of the Gemini partnership: the National Science Foundation (United States), the National Research Council (Canada), CONICYT (Chile), Ministerio de Ciencia, Tecnología e Innovación Productiva (Argentina), and Ministério da Ciência, Tecnologia e Inovação (Brazil). This work has made use of data from the European Space Agency (ESA) mission Gaia (, processed by the Gaia Data Processing and Analysis Consortium (DPAC, Funding for the DPAC has been provided by national institutions, in particular the institutions participating in the Gaia Multilateral Agreement. Funding for this program is provided by NASA through the NASA Astrophysical Data Analysis Program, award NNX12AE15GB. Support for this work was provided by the National Aeronautics and Space Administration through Chandra Award Number GO6-17130X issued by the Chandra X-ray Center, which is operated by the Smithsonian Astrophysical Observatory for and on behalf of the National Aeronautics Space Administration under contract NAS8-03060.


  • Allen et al. (2011) Allen, S. W., Evrard, A. E., & Mantz, A. B. 2011, ARA&A, 49, 409, doi: 10.1146/annurev-astro-081710-102514
  • Andersson et al. (2011) Andersson, K., Benson, B. A., Ade, P. A. R., et al. 2011, ApJ, 738, 48, doi: 10.1088/0004-637X/738/1/48
  • Andreon (2010) Andreon, S. 2010, MNRAS, 407, 263, doi: 10.1111/j.1365-2966.2010.16856.x
  • Arnaud et al. (2010) Arnaud, M., Pratt, G. W., Piffaretti, R., et al. 2010, A&A, 517, A92, doi: 10.1051/0004-6361/200913416
  • Ashby et al. (2009) Ashby, M. L. N., Stern, D., Brodwin, M., et al. 2009, ApJ, 701, 428, doi: 10.1088/0004-637X/701/1/428
  • Bertin & Arnouts (1996) Bertin, E., & Arnouts, S. 1996, A&AS, 117, 393, doi: 10.1051/aas:1996164
  • Bleem et al. (2015) Bleem, L. E., Stalder, B., de Haan, T., et al. 2015, ApJS, 216, 27, doi: 10.1088/0067-0049/216/2/27
  • Brodwin et al. (2013) Brodwin, M., Stanford, S. A., Gonzalez, A. H., et al. 2013, ApJ, 779, 138, doi: 10.1088/0004-637X/779/2/138
  • Brodwin et al. (2015) Brodwin, M., Greer, C. H., Leitch, E. M., et al. 2015, ApJ, 806, 26, doi: 10.1088/0004-637X/806/1/26
  • Bruzual & Charlot (2003) Bruzual, G., & Charlot, S. 2003, MNRAS, 344, 1000, doi: 10.1046/j.1365-8711.2003.06897.x
  • Chabrier (2003) Chabrier, G. 2003, PASP, 115, 763, doi: 10.1086/376392
  • Chambers et al. (2016) Chambers, K. C., Magnier, E. A., Metcalfe, N., et al. 2016, ArXiv e-prints.
  • Chiu et al. (2016) Chiu, I., Mohr, J., McDonald, M., et al. 2016, MNRAS, 455, 258, doi: 10.1093/mnras/stv2303
  • Chiu et al. (2018) Chiu, I., Mohr, J. J., McDonald, M., et al. 2018, MNRAS, 478, 3072, doi: 10.1093/mnras/sty1284
  • Conroy et al. (2007) Conroy, C., Wechsler, R. H., & Kravtsov, A. V. 2007, ApJ, 668, 826, doi: 10.1086/521425
  • Cutri (2013) Cutri, R. M. 2013, VizieR Online Data Catalog, 2328, 0
  • Davis et al. (2007) Davis, M., Guhathakurta, P., Konidaris, N. P., et al. 2007, ApJ, 660, L1, doi: 10.1086/517931
  • de Haan et al. (2016) de Haan, T., Benson, B. A., Bleem, L. E., et al. 2016, ApJ, 832, 95, doi: 10.3847/0004-637X/832/1/95
  • Eisenhardt et al. (2004) Eisenhardt, P. R., Stern, D., Brodwin, M., et al. 2004, ApJS, 154, 48, doi: 10.1086/423180
  • Ettori et al. (2006) Ettori, S., Dolag, K., Borgani, S., & Murante, G. 2006, MNRAS, 365, 1021, doi: 10.1111/j.1365-2966.2005.09788.x
  • Fazio et al. (2004) Fazio, G. G., Hora, J. L., Allen, L. E., et al. 2004, ApJS, 154, 10, doi: 10.1086/422843
  • Gaia Collaboration et al. (2016) Gaia Collaboration, Prusti, T., de Bruijne, J. H. J., et al. 2016, A&A, 595, A1, doi: 10.1051/0004-6361/201629272
  • Gaia Collaboration et al. (2018) Gaia Collaboration, Brown, A. G. A., Vallenari, A., et al. 2018, A&A, 616, A1, doi: 10.1051/0004-6361/201833051
  • Gettings et al. (2012) Gettings, D. P., Gonzalez, A. H., Stanford, S. A., et al. 2012, ApJ, 759, L23, doi: 10.1088/2041-8205/759/1/L23
  • Giodini et al. (2009) Giodini, S., Pierini, D., Finoguenov, A., et al. 2009, ApJ, 703, 982, doi: 10.1088/0004-637X/703/1/982
  • Gonzalez et al. (2013) Gonzalez, A. H., Sivanandam, S., Zabludoff, A. I., & Zaritsky, D. 2013, ApJ, 778, 14, doi: 10.1088/0004-637X/778/1/14
  • Gonzalez et al. (2007) Gonzalez, A. H., Zaritsky, D., & Zabludoff, A. I. 2007, ApJ, 666, 147, doi: 10.1086/519729
  • Gonzalez et al. (2015) Gonzalez, A. H., Decker, B., Brodwin, M., et al. 2015, ApJ, 812, L40, doi: 10.1088/2041-8205/812/2/L40
  • Gonzalez et al. (2018) Gonzalez, A. H., Gettings, D. P., Brodwin, M., et al. 2018, ArXiv e-prints.
  • Hainline et al. (2011) Hainline, L. J., Blain, A. W., Smail, I., et al. 2011, ApJ, 740, 96, doi: 10.1088/0004-637X/740/2/96
  • He et al. (2005) He, P., Feng, L.-L., & Fang, L.-Z. 2005, ApJ, 623, 601, doi: 10.1086/428708
  • Hilton et al. (2013) Hilton, M., Hasselfield, M., Sifón, C., et al. 2013, MNRAS, 435, 3469, doi: 10.1093/mnras/stt1535
  • Högbom (1974) Högbom, J. A. 1974, A&AS, 15, 417
  • Hook et al. (2004) Hook, I. M., Jørgensen, I., Allington-Smith, J. R., et al. 2004, PASP, 116, 425, doi: 10.1086/383624
  • Jannuzi & Dey (1999) Jannuzi, B. T., & Dey, A. 1999, in Astronomical Society of the Pacific Conference Series, Vol. 191, Photometric Redshifts and the Detection of High Redshift Galaxies, ed. R. Weymann, L. Storrie-Lombardi, M. Sawicki, & R. Brunner, 111
  • Krause et al. (2012) Krause, E., Pierpaoli, E., Dolag, K., & Borgani, S. 2012, MNRAS, 419, 1766, doi: 10.1111/j.1365-2966.2011.19844.x
  • Kravtsov & Borgani (2012) Kravtsov, A. V., & Borgani, S. 2012, ARA&A, 50, 353, doi: 10.1146/annurev-astro-081811-125502
  • Kravtsov et al. (2005) Kravtsov, A. V., Nagai, D., & Vikhlinin, A. A. 2005, ApJ, 625, 588, doi: 10.1086/429796
  • Lacy et al. (2005) Lacy, M., Wilson, G., Masci, F., et al. 2005, ApJS, 161, 41, doi: 10.1086/432894
  • Lin et al. (2003) Lin, Y.-T., Mohr, J. J., & Stanford, S. A. 2003, ApJ, 591, 749, doi: 10.1086/375513
  • Lin et al. (2012) Lin, Y.-T., Stanford, S. A., Eisenhardt, P. R. M., et al. 2012, ApJ, 745, L3, doi: 10.1088/2041-8205/745/1/L3
  • Mancone & Gonzalez (2012) Mancone, C. L., & Gonzalez, A. H. 2012, PASP, 124, 606, doi: 10.1086/666502
  • Mancone et al. (2010) Mancone, C. L., Gonzalez, A. H., Brodwin, M., et al. 2010, ApJ, 720, 284, doi: 10.1088/0004-637X/720/1/284
  • Marrone et al. (2012) Marrone, D. P., Smith, G. P., Okabe, N., et al. 2012, ApJ, 754, 119, doi: 10.1088/0004-637X/754/2/119
  • McLean et al. (2010) McLean, I. S., et al. 2010, in Proc. SPIE, Vol. 7735, Ground-based and Airborne Instrumentation for Astronomy III, 77351E–77351E–12
  • McLean et al. (2012) McLean, I. S., et al. 2012, in Proc. SPIE, Vol. 8446, Ground-based and Airborne Instrumentation for Astronomy IV, 84460J
  • Mo et al. (2018) Mo, W., Gonzalez, A., Stern, D., et al. 2018, ArXiv e-prints.
  • Muzzin et al. (2008) Muzzin, A., Wilson, G., Lacy, M., Yee, H. K. C., & Stanford, S. A. 2008, ApJ, 686, 966, doi: 10.1086/591542
  • Nurgaliev et al. (2017) Nurgaliev, D., McDonald, M., Benson, B. A., et al. 2017, ApJ, 841, 5, doi: 10.3847/1538-4357/aa6db4
  • Oke et al. (1995) Oke, J. B., Cohen, J. G., Carr, M., et al. 1995, PASP, 107, 375, doi: 10.1086/133562
  • Pillepich et al. (2017) Pillepich, A., Nelson, D., Hernquist, L., et al. 2017, ArXiv e-prints.
  • Planelles et al. (2013) Planelles, S., Borgani, S., Dolag, K., et al. 2013, MNRAS, 431, 1487, doi: 10.1093/mnras/stt265
  • Poole et al. (2007) Poole, G. B., Babul, A., McCarthy, I. G., et al. 2007, MNRAS, 380, 437, doi: 10.1111/j.1365-2966.2007.12107.x
  • Rudy et al. (1987) Rudy, D. J., Muhleman, D. O., Berge, G. L., Jakosky, B. M., & Christensen, P. R. 1987, Icarus, 71, 159, doi: 10.1016/0019-1035(87)90170-9
  • Salpeter (1955) Salpeter, E. E. 1955, ApJ, 121, 161, doi: 10.1086/145971
  • Schechter (1976) Schechter, P. 1976, ApJ, 203, 297, doi: 10.1086/154079
  • Scoville et al. (2007) Scoville, N., Aussel, H., Brusa, M., et al. 2007, ApJS, 172, 1, doi: 10.1086/516585
  • Stanford et al. (2014) Stanford, S. A., Gonzalez, A. H., Brodwin, M., et al. 2014, ApJS, 213, 25, doi: 10.1088/0067-0049/213/2/25
  • Sunyaev & Zeldovich (1970) Sunyaev, R. A., & Zeldovich, Y. B. 1970, Comments on Astrophysics and Space Physics, 2, 66
  • Sunyaev & Zeldovich (1972) —. 1972, Comments on Astrophysics and Space Physics, 4, 173
  • van der Burg et al. (2014) van der Burg, R. F. J., Muzzin, A., Hoekstra, H., et al. 2014, A&A, 561, A79, doi: 10.1051/0004-6361/201322771
  • White et al. (1993) White, S. D. M., Navarro, J. F., Evrard, A. E., & Frenk, C. S. 1993, Nature, 366, 429, doi: 10.1038/366429a0
  • Wright et al. (2010) Wright, E. L., Eisenhardt, P. R. M., Mainzer, A. K., et al. 2010, AJ, 140, 1868, doi: 10.1088/0004-6256/140/6/1868
  • Wylezalek et al. (2014) Wylezalek, D., Vernet, J., De Breuck, C., et al. 2014, ApJ, 786, 17, doi: 10.1088/0004-637X/786/1/17
  • Zhang et al. (2011) Zhang, Y.-Y., Laganá, T. F., Pierini, D., et al. 2011, A&A, 535, A78, doi: 10.1051/0004-6361/201116803
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description