Advanced
MASSIVE STRUCTURES OF GALAXIES AT HIGH REDSHIFTS IN THE GREAT OBSERVATORIES ORIGINS DEEP SURVEY FIELDS
MASSIVE STRUCTURES OF GALAXIES AT HIGH REDSHIFTS IN THE GREAT OBSERVATORIES ORIGINS DEEP SURVEY FIELDS
Journal of The Korean Astronomical Society. 2015. Feb, 48(1): 21-55
Copyright © 2015, The Korean Astronomical Soiety
  • Received : August 25, 2014
  • Accepted : December 15, 2014
  • Published : February 28, 2015
Download
PDF
e-PUB
PubReader
PPT
Export by style
Share
Article
Author
Metrics
Cited by
TagCloud
About the Authors
EUGENE KANG
Korea Educational Broadcasting System, 2748, Nambusunhwanro, Gangnam-Gu, Seoul 135-854, Korea
MYUNGSHIN IM
CEOU, Astronomy Program, Department of Physics and Astronomy, Seoul National University, 1 Gwanak-ro, Gwanak-gu, Seoul 151-742, Korea;mim@astro.snu.ac.kr

Abstract
If the Universe is dominated by cold dark matter and dark energy as in the currently popular CDM cosmology, it is expected that large scale structures form gradually, with galaxy clusters of mass M 10 14 M appearing at around 6 Gyrs after the Big Bang ( z ~ 1). Here, we report the discovery of 59 massive structures of galaxies with masses greater than a few times 10 13 M at redshifts between z = 0.6 and 4.5 in the Great Observatories Origins Deep Survey fields. The massive structures are identified by running top-hat filters on the two dimensional spatial distribution of magnitude-limited samples of galaxies using a combination of spectroscopic and photometric redshifts. We analyze the Millennium simulation data in a similar way to the analysis of the observational data in order to test the CDM cosmology. We find that there are too many massive structures ( M > 7×10 13 M ) observed at z > 2 in comparison with the simulation predictions by a factor of a few, giving a probability of < 1/2500 of the observed data being consistent with the simulation. Our result suggests that massive structures have emerged early, but the reason for the discrepancy with the simulation is unclear. It could be due to the limitation of the simulation such as the lack of key, unrecognized ingredients (strong non-Gaussianity or other baryonic physics), or simply a difficulty in the halo mass estimation from observation, or a fundamental problem of the CDM cosmology. On the other hand, the over-abundance of massive structures at high redshifts does not favor heavy neutrino mass of ~ 0.3 eV or larger, as heavy neutrinos make the discrepancy between the observation and the simulation more pronounced by a factor of 3 or more.
Keywords
1. INTRODUCTION
Currently, the energy density of the Universe is believed to be dominated by dark energy (73%), and cold dark matter (23%), based on various pieces of observational evidence ( Im et al. 1997 ; Chiba & Yoshii 1997 ; Riess et al. 1998 ; Perlmutter et al. 1999 ; Chae et al. 2002 ; Eisenstein et al. 2005 ; Mantz et al. 2010 ; Komatsu et al. 2011 ; Cay´on et al. 2011 ). In such a Universe, called the ΛCDM cosmology, the formation of galaxies and large scale structures of galaxies proceeds from less massive objects to more massive objects through hierarchical mergers. Thus, based on the ΛCDM model, it is expected that massive galaxies and most groups and clusters of galaxies appear relatively late, at z < 1 (e.g., Allen et al. 2011 ; Kravtsov et al. 2012 ; Park et al. 2012 ; De Lucia & Blaizot 2007 ; Im et al. 2002 ). Analytical calculations of halo mass functions and numerical simulations of galaxy formation support this expectation. Therefore, the evolution of the number density of massive structures can serve as a useful test of the ΛCDM cosmology ( Bahcall et al. 1999 ; Haiman et al. 2001 ). Furthermore, the epoch of the emergence of massive structures depends on other physical conditions that are poorly constrained so far such as the initial density fluctuations ( Reed et al. 2007 ) and the neutrino mas ( Bond et al. 1980 ; reviewed in Lesgourgues & Pastor 2006 ), and hence its study can provide new insights on these physical conditions as well.
Current searches for high redshift clusters are unveiling massive clusters and proto-clusters at high redshifts, making it possible to investigate the evolution of massive structures of galaxies throughout the age of the Universe and to test the predictions of the ΛCDM models. Various cluster-search techniques exist, such as searches using extended X-ray sources, redsequence galaxies, the Sunyaev-Zeldovich (SZ) effect, and narrow-band imaging around high redshift radio galaxies, AGNs and Lyman break galaxies (LBGs). The cluster searches using X-ray emissions are based on the fact that a massive, gravitationally bound object can be identified as an extended X-ray source since deep gravitational potential wells of galaxy clusters provide the necessary conditions for the X-ray emitting hot gas to exist. Spectroscopic observation of galaxies associated with extended X-ray sources has revealed clusters of galaxies out to z ~ 2 (e.g., Stanford et al. 2006 ; Hilton et al. 2007 ; Fassbender et al. 2011 ). The redsequence technique, on the other hand, utilizes the fact that galaxies (mostly early-type) in a cluster form a tight sequence in the color-magnitude parameter space. If early-type galaxies in clusters formed the bulk of their stars at a similar epoch at high redshift ( z > 2), one would expect to find a red sequence in high redshift clusters too. Therefore, one can identify cluster/protoclusters at high redshift by searching for areas that are populated by red galaxies with a certain color that is expected for passively evolving galaxies at high redshift. Groups/clusters/proto-clusters of galaxies are identified in this way out to z ~ 2 ( Gladders & Yee 2005 ; McCarthy et al. 2007 ; Goto et al. 2008 ; Muzzin et al. 2009 ; Bielby et al. 2010 ; Papovich et al. 2010 ), although in most cases, the identified structures are at z < 1 ( Muzzin 2008 ; Gilbank et al. 2011 ; Li et al. 2012 ). Going out to higher redshifts, radio galaxies or quasars are being used as signposts for massive structures. Radio galaxies at low redshift are often massive early-type galaxies that are in cluster regions, and if a similar trend holds at high redshift, one can expect to see massive structures around radio galaxies. A similar argument applies to quasars too. Overdense areas are searched within a certain projected radius around high redshift radio galaxies or quasars by looking for excessive number of extremely red objects (EROs), BzK galaxies, LBGs, or emission line galaxies through optical/NIR broad band imaging or narrow-band imaging centered at the wavelength where emission lines (e.g., Ly α , H α ) show up at the radio galaxy’s redshift. Such studies reveal possible large scale structures as far as z ~ 6 ( Le Fevre et al. 1996 ; Kurk et al. 2000 ; Pentericci et al. 2000 ; Wold et al. 2003 ; Kajisawa et al. 2006 ; Kodama et al. 2007 ; Venemans et al. 2007 ; Doherty et al. 2010 ; Capak et al. 2011 ; Kuiper et al. 2011 ; Toshikawa et al. 2012 ). Recently, the Sunyaev-Zeldovich effect has been utilized to uncover clusters. These searches are now revealing clusters out to z ~ 1.9, but mostly at z < 1, with reliable estimates on cluster mass (Foley et al. 2011; Sehgal et al. 2011; Benson et al. 2013 ; Williamson et al. 2011 ; Menanteau et al. 2012 ; Stalder et al. 2013 ; Bayliss et al. 2014 ).
Using the above techniques, a handful of massive clusters of galaxies have been identified at z ~ 1.5, and as several previous works emphasize, the implication of the existence of massive clusters at high redshift is quite tantalizing. While the ΛCDM models can successfully explain the abundance of galaxy clusters at z < 1 (e.g., Benson et al. 2013 ; Williamson et al. 2011 ) the situation is more controversial at higher redshift. Several studies suggest that there are too many massive clusters at z ~ 1.5 compared to theoretical expectations based on the ΛCDM cosmology ( Jee et al. 2009 ; Papovich et al. 2010 ; Fassbender et al. 2011 ; Nastasi et al. 2011 ; Holz et al. 2012 ; Gonzalez et al. 2012 ; Menanteau et al. 2012 ). This does not necessarily mean the failure of the ΛCDM model, but could reflect either our poor understanding of the physics involved in the growth of massive structures and the galaxies within them, or limitations of current observational studies. For example, the studies disfavoring the ΛCDM cosmology using clusters mainly focus on a small number of the most massive structures. Thus, they could be subject to bias from small number statistics ( Biviano et al. 2006 ). Other explanations exist which try to match the over-abundance of massive clusters at high redshift using exotic forms of dark energy ( Baldi 2012 ; Carlesi et al. 2012 ). On the other hand, different studies suggest that there is no significant lack of high redshift massive clusters ( Capak et al. 2011 ; Coupon et al. 2012 ; Williamson et al. 2011 ; Stalder et al. 2013 ; Bayliss et al. 2014 ), and that the discovery of massive high redshift clusters so far is consistent with the ΛCDM cosmology ( Mortonson et al. 2011 ; Hotchkiss 2011 ; Waizmann et al. 2012 , 2013 ; Harrison et al. 2012 , 2013 ). Extending these kinds of studies to higher redshifts can place stronger constraints on the ΛCDM cosmology, but there is a general lack of unbiased searches for high redshift clusters/proto-clusters at z > 2. Searches based on radio galaxies or quasars could miss massive structures where radio galaxies or quasars are not found.
One promising way to uncover high redshift massive structures in an unbiased way is to utilize spectroscopic/photometric redshift information from multiwavelength data (e.g., Kang & Im 2009 ; Cucciati et al. 2014 ). Mixing spectroscopic redshift information with photometric redshifts has disadvantages, such as the dilution of signals from over-dense areas and the confusion of foreground and background structures into a single structure. However, such disadvantages can be overcome if one also applies a similar observational conditions to simulation data in order to interpret the observational results.
In order to place constraints on the formation of large scale structures at high redshift and see how well ΛCDM models can reproduce the observed trend, we embarked on a blind search of massive structures of galaxies (hereafter, MSGs) out to z ~ 4.5, extending the redshift range of MSG searches well beyond the red-shift limit of existing works ( z = 2). Here, we use the term “massive structures of galaxies” instead of terms such as clusters, massive halos, or proto-clusters that are commonly used in the literature, since it is dif- ficult to separate each halo or cluster when they are undergoing hierarchical build-up. The search is performed on the Great Observatories Origins Deep Survey (GOODS) fields, which include two fields, each covering 160 arcmin 2 . These fields offer a wealth of ultra-deep, multi-wavelength imaging and spectroscopic datasets that are necessary to identify and physically characterize distant objects. The area coverage is large enough to mitigate the cosmic variance problem, if MSGs are searched over a sufficiently large depths along the redshift direction (~ a few hundred Mpc). The GOODSSouth is the field where we previously identified an overdensity of galaxies at z ~ 3.7 ( Kang & Im 2009 ), by using photometric and spectroscopic redshifts derived from the multi-wavelength data ranging from U −band through 8 µ m band of the Spitzer. Our work basically builds on the previous work of Kang & Im (2009) , extending both the redshift range and the mass range of the search for MSGs. We will derive the number density of MSGs out to z ~ 4.5, and compare them to the Millennium Simulation ( Springel et al. 2005 ) where we impose the same selection criteria as those used in the analysis of the observed data. The use of the simulation data is essential, to minimize selection effects associated with the MSG identification process. In the analysis of the observed data, we will adopt the cosmological parameter values identical to those adopted by the Millennium Simulation, i.e., Ω m = Ω dm +Ω b = 0.25, Ω b = 0.045, h = 0.73, and Ω Λ = 0.75 ( Colless et al. 2001 ; Spergel et al. 2003 ; Seljak et al. 2005 ). Note that magnitudes are in AB units throughout the paper.
2. DATA
- 2.1. Imaging Data and Catalog
The GOODS fields cover two separate regions in the northern and the southern hemispheres centered on 12 h 36 m 55 s , +62°14 m 15 s (J2000) of the Hubble Deep Field-North (HDF-N) and 3 h 32 m 30 s , −27°48 m 20 s (J2000) of the Chandra Deep Field-South (CDF-S), respectively. Each field provides an area coverage of 160 arcmin 2 . A wealth of multi-wavelength imaging data is available in the GOODS fields, and we utilize these for our study.
For the optical imaging data, we take the version 1.1 of the HST Advanced Camera for Surveys (ACS) multi-band source catalog (F435W, F606W, F775W, F850LP) released by the GOODS team ( Giavalisco et al. 2004a ). The HST data are augmented by groundbased, deep, U -band data taken with VIMOS on the VLT with a depth of 28.0 mag at 5 σ within a 1′′ radius aperture ( Nonino et al. 2009 ), or with the MOSAIC prime focus camera of the Kitt Peak National Observatory (KPNO) 4m telescope to the depth of U = 27.1 mag ( Capak et al. 2004 ).
For the GOODS-South area, we use the publicly available deep J -, H - and Ks -band images ( Retzlaff et al. 2010 ). These images were taken by the ESOGOODS team using Infrared Spectrometer And Array Camera (ISAAC) on the Very Large Telescope (VLT). The depth of the VLT/ISAAC data reaches to J = 25.2 mag, H and Ks = 24.4 mags at 5 σ . We use the version 2.0 released images which are made of 24 fields in the J -, H -bands and 26 fields in Ks -band.
For the GOODS-North area, we use the NIR source catalog ( J , Ks ; e.g., Wang et al. 2010 ) that is based on images taken by the Wide Field Infrared Cam- era (WIRCam) on the Canada-France-Hawaii Tele- scope (CFHT). The depth of the CFHT/WIRCam data reaches to J = 24.6 mag and Ks = 24.0 mag at 5 σ , which are comparable to the depths of the GOODS-South NIR images.
The imaging data at 3.6, 4,5, 5.8, and 8.0 µ m are taken by the Spitzer Infrared Array Camera (IRAC), and we use the photometry in the catalog created by R. Chary (private communication) in both the northern and the southern fields. Table 1 summarizes the available multi-wavelength imaging data, along with their depths and angular resolutions.
Summary of the multiwavelength imaging data
PPT Slide
Lager Image
a: Filter central wavelength; b: Capak et al. (2004); c: Nonino et al. (2009); d: Giavalisco et al. (2004a) and references therein.; e: Retzlaff et al. (2010); f: Wang et al. (2010); g: M. Dickinson et al. in preparation, R. Chary (private communication)
To derive photometry of objects in U , J , H , and Ks , we ran SExtractor ( Bertin & Arnouts 1996 ). The source detection threshold was set to a minimum contiguous area of 16 pixels, with each pixel value above 1.0 σ . For deblending parameter, we use 0.0075 of minimum contrast, and the background mesh size is set to 128 pixels. We used Ks -band detections as a reference, and derived photometry in the other bands using SEx-tractor in dual-mode at the Ks -band detection position. Auto-magnitudes, generally representing total magnitudes of sources, are used for the photometry, since the multi-wavelength images have a wide range of spatial resolutions.
The cross-matching of the objects in the multiwavelength data is done by matching U , Ks and IRAC-band catalogs against the objects in the HST ACS catalog using
PPT Slide
Lager Image
and
PPT Slide
Lager Image
matching radii, respectively. The matching radius is defined as the radius where the ACS sources start to be multiply matched beyond ~ 50% of cases. We find that our adopted cross-match radius provides counterparts in 98% and 87% of the Ks -band limited objects against the ACS-band and IRAC-band sources, respectively. Multiple matches are examined visually on the image to determine which of the multiple matches should be chosen as a good match. The objects that have the smallest positional difference between the different wavelength images are generally chosen to be good matches. In some cases, multiple objects in an ACS image are blended together in Ks or IRAC images, making it impossible to construct a reliable spectral energy distribution (SED). Such cases are excluded from the analysis, which represent 4% and 5% of the Ks -band and IRAC-band detected objects, respectively, for a K -band limited sample at Ks < 24 mag.
- 2.2. Spectroscopic Redshifts
Galaxies in both the GOODS-South and the GOODS-North fields have been subject to extensive spectroscopic follow-up, and a fair number of spectroscopic redshifts are available. The observations were carried out on a heterogeneous conglomeration of samples, including magnitude-limited samples in specific bands, or color-selected objects such as LBGs. Nevertheless, we gather all available redshifts for our analysis.
In the GOODS-South, a total of 3166 spectroscopic redshifts are available. These are mostly taken with the FOcal Reducer and low dispersion Spectrograph (FORS2) and with the VIsible Multi-Object Spectrograph (VIMOS) on the VLT ( Vanzella et al. 2005 , 2006 , 2008 , 2009 ; Popesso et al. 2009 ; Balestra et al. 2010 ), by the ESO-GOODS team. The sample selection ranges from z -band selected, BzK selected, U -, B -, V -, i -band dropouts and X-ray sources. Many spectroscopic observations were carried out in the GOODS-North field too, with the number of spectroscopic redshifts being 2843. Barger et al. (2008) assembled a catalog containing various spectroscopic studies including redshift surveys of z -band or r -band limited objects (TKRS, e.g., Wirth et al. 2004; Cooper et al. 2011 ), and color-selected galaxies at high redshift (Reddy et al. 2006). To the redshifts available in the literature, they added redshifts from their own observations in an effort to increase the completeness of the spectroscopic coverage of the GOODS-North field. The redshift completeness is over 90% at z 850 = 23.3 mag, or at Ks = 21.6 mag.
Table 2 summarizes the spectroscopic redshift data in the GOODS fields.
Summary of the spectroscopic data
PPT Slide
Lager Image
a: Vanzella et al. (2005, 2006, 2008, 2009), Popesso et al. (2009), Balestra et al. (2010) b: Barger et al. (2008), Cooper et al. (2011), and references therein.
Figure 1 shows the Ks -band magnitudes versus spectroscopic redshifts of galaxies in the GOODS fields. Also plotted in the same figure are photometric redshifts, which are derived from the multi-wavelength imaging data as described in the next section. Using the photometric redshifts, we estimate the completeness of the spectroscopic redshift measurements in the redshift versus apparent magnitude parameter space, which is presented in Figure 2 for each GOODS field separately (see also Table 3 ). Later, this information will be used to add photometric redshift noise in the simulation data.
PPT Slide
Lager Image
The Ks-band magnitude versus spectroscopic red-shift (crosses enclosed by red diamonds) and photometric redshifts (crosses) for the GOODS-Norh (top) and -South (bottom).
PPT Slide
Lager Image
The completeness of the spectroscopic identification (the ratio of the number of spectroscopically confirmed galaxies to all galaxies) at different Ks magnitude bins as a function of redshift for the GOODS-Norh (top) and the GOODS-South (bottom).
Completeness of the spectroscopic identification
PPT Slide
Lager Image
a: The ratio of the number of spectroscopically confirmed galaxies to all galaxies at different Ks magnitude/redshift bins for the GOODS-North field. b: The ratio of the number of spectroscopically confirmed galaxies to all galaxies at different Ks magnitude/redshift bins for the GOODS-South field.
3. PHOTOMETRIC REDSHIFTS AND SED FITTING
We derive photometric redshifts of galaxies to aid the selection of MSGs. Galaxies are primarily drawn from a Ks -band limited sample, but are supplemented by z -band selected galaxies, the latter being basically U - and B -band dropout galaxies ( Giavalisco et al. 2004a ; Nonino et al. 2009 ). In the Ks -band limited sample, galaxies with Ks ≤ 23.8 mag (GOODS-South) and Ks ≤ 24 mag (GOODS-North) are used. Spectroscopic redshifts are available for 30% of the magnitude-limited sample, and for galaxies without spectroscopic redshifts, we derive photometric redshifts using U , B 435 , V 606 , i 775 , z 850 , J , H , Ks bands, and the Spitzer photometry at 3.6 through 8.0 µ m if available.
To derive photometric redshifts, we generated 742 SED templates using the stellar population synthesis code GALAXEV ( Bruzual & Charlot 2003 ). The assumed parameters for these SED templates are the following. For the star formation rates, we use a single burst, or a constant star formation rate (SFR) or an exponentially declining SFR in the form of SFR ∝ exp(-t/ τ ) for τ = 1, 2, 3, 5, 15, or 30 Gyr. For the metallicity, we assumed values of 0.4 Z , 1 Z or 2 Z , and we adopted the Salpeter initial mass function with a mass range of 0.1 to 100 M . The SED ages range from 100 Myr to 10 Gyr with an upper limit (the age of the Universe at the corresponding redshift). The reddening parameter AV is allowed to change from 0 to 1.8 with the optical depth, τ λ , following the two component power law model of Charlot & Fall (2000) . This model fits the extinction properties of star-forming galaxies well, and adopts two different optical depths depending on the age of stars (younger or older than 10 7 years). The Bayesian photometric redshift estimation (BPZ; Bentez 2000 ) code is used for deriving photometric redshifts. During the photometric redshift derivation, the best-fit templates are identified, which are used for estimating stellar masses of galaxies. We also obtain stellar masses of galaxies with spectroscopic redshifts by fitting their SEDs using the same code. The age of the template SEDs is restricted to the age of the Universe at the corresponding redshift when performing the SED fit.
To test the accuracy of the photometric redshifts, we compare our photometric redshifts ( z phot ) with the available spectroscopic redshifts ( z spec ). Here, we use only objects having spectroscopic redshifts with high quality flags (B or better: Vanzella et al. 2006 , 2008 , 2009 ; Poppesso et al. 2009 ; Balestra et al. 2010 , 3 or better: Barger et al. 2008 ; Cooper et al. 2011 ). The cross-match of the photometric redshift sample with the z spec data gives 6009 galaxies for this comparison. Figure 3 shows the comparison of the spectroscopic redshifts versus the photometric redshifts. In most cases, photometric redshifts are reasonably accurate with a small number of outliers, where the root-mean square (rms) scatter of z phot versus z spec , σz , is σz ≃ 0.12 at low redshift ( z < 1), and σz ≃ 0.2 at high redshift ( z > 2). At the intermediate redshift range, the photometric redshift accuracy is poor with a value of σz ≃ 0.4. Figure 4 shows the rms scatter of z phot against z spec as a function of z spec . We find that the average absolute scatter is 〈|∆z|〉 = 〈(| z spec z phot )/(1 + z spec )|〉 of 0.06 after excluding several clear outliers, which is slightly worse than, but comparable to the values (〈|∆ z |〉) from other photometric redshift studies ( Wuyts et al. 2008 ; Santini et al. 2009 ).
PPT Slide
Lager Image
Photometric redshift (zphot) versus spectroscopic redshift (zspec) over the combined areas of the two GOODS fields. The diamonds and crosses indicate the Ks-band and z-band limited samples, respectively.
PPT Slide
Lager Image
The rms scatter, σz, of the photometric redshift (σz) as a function of spectroscopic redshift. The σz values stay at around 0.1 to 0.2 for most redshift intervals, except for zspec = 1.5 − 2.0 where the σz value is quite large at around 0.4.
4. MSGS IN GOODS FIELDS
- 4.1. Search Method
We searched for MSGs by running circular top-hat filters with a 1 Mpc diameter in physical size on a two-dimensional spatial distribution of galaxies within a fixed redshift range. The procedure is described below, and graphically illustrated in Figure 5 . Note that the length scale of the analysis is set to physical scale rather than the comoving scale throughout the paper, unless mentioned otherwise.
PPT Slide
Lager Image
This figure illustrates how MSGs are identified. The left panel shows top-hat filters (circles) that are moved around over the two-dimensional distribution of galaxies in a specific redshift bin. We obtain the surface number density of galaxies in each circle by counting the number of objects within each redshift slice and dividing by the area covered by the circle. The surface number densities are computed over the entire image, and the distribution of the surface number density is plotted (right panel). The rms scatter, σrms of the distribution is determined from the distribution, and the significance of the overdensity, the overdensity factor, is determined as the number of σ’s from the mean surface number density. MSGs are defined as the location where the surface density exceeds 3.5 σrms (right panel).
The redshift bins are chosen to have a width of ±2000 km s −1 for galaxies with spectroscopic redshifts, or a width corresponding to the median photometric redshift accuracy at a certain redshift range ( σz ) as presented in Figure 4 .
The center of each successive redshift bin is shifted by ∆ v = 3000 km s −1 from 0.65 to 4.5. The photometric redshift uncertainty of 0.1 at z = 1 corresponds to a velocity uncertainty of 15000 km s −1 . So, a shift of ∆ v = 3000 km s −1 is fine enough not to miss any overdensities between two successive redshift bins. For the spectroscopic redshift sample, this shift allows an overlap between two successive redshift bins and the identification of MSGs that are situated between two successive redshift bins, since the velocity dispersion of a typical galaxy cluster is σv ~ 1000 km s −1 .
The density contour maps are generated by counting the number of objects around a point within a 1 Mpc diameter circle in the image. The top-hat filter is moved around over a 100 by 100 grid for each GOODS field. The grid dimension gives a separation of each contour map center of 0.2 arcmin or 0.1 Mpc at z = 2. The separation between the grid nodes is sufficiently smaller than the typical scales of groups/clusters. The distributions of the projected number densities are created, and the rms of the distribution ( σrms ) is determined after a σ clipping, with the clipping at 3- σ over 3 iterations. Note that we also obtain the standard deviation σ G of the same projected number density distribution by performing a Gaussian fit, and the σ G values derived in this way agree with the values of σrms . Finally we identify MSGs as the regions where the projected number density exceeds the mean projected number density by more than 3.5 σrms . Here, we define the overdensity factor as the number of σrms units the number density of an overdense region deviates from the mean number density.
The top-hat filter diameter may affect the search result by changing the overdensity factors. If the search radius is too large or too small, the overdensity factors may not reach the MSG identification threshold. Therefore, we performed the analysis by varying the top-hat filter diameter from 0.5 Mpc to 2 Mpc to check for any dependence of our result on the adopted diameter of the top-hat filter, but our test shows no signifi- cant difference in the overdensity factor of MSGs over this range in scale ( Figure 6 ). We also checked how the number of MSGs changes depending on the filter diameter. With a 0.5 Mpc diameter, the number of MSGs increases by ~10% at z > 2.9, while there is no change in the number of MSGs at lower redshifts. The increase in the MSG number is mostly at lower mass regions, i.e., we tend to find less massive MSGs when using a top-hat filter with a smaller diameter. On the other hand, if a large diameter of 2.0 Mpc is used, we miss ~ 30% of MSGs at z < 2, meaning that larger search diameters tend to miss concentrated structures. Hence, the 1 Mpc diameter appears to be optimal for finding MSGs.
PPT Slide
Lager Image
Examples showing how the overdensity factor changes as a function of a diameter of the top-hat circular filter. The overdensity factor on the y-axis is in units of σrms. From 0.5 Mpc to 2.0 Mpc, we find that the over-density factor is not very sensitive to the assumed diameter, although the adopted 1 Mpc seems to give an optimal result. The examples are given for typical MSGs at z = 0.95, 2.55 and 3.7 in GOODS-South.
Foreground or background halos or filamentary structures along the line of sight could produce strong clustering signal that can be identified as MSGs, especially when photometric redshifts are the main source of the redshift identification. We examined such a possibility using a simulation dataset (see Section 5 and Figure 12 ) that is similarly constructed to mimic the observation. We find that the identification of MSGs is mostly due to one strong massive halo (see the top and middle panels of Figure 12 ), rather than many halos with similar masses scattered along the line of sight. About 30% of the identified MSGs consist of several distinctive peaks along the line of sight, but even in such cases, the distinctive halos are closely located in velocity and real space (the bottom panel of Figure 12 ). However, the numerous minor halos along the line of sight can contribute significantly to the mass estimate of MSGs, but such effects are accounted for by estimating the interloper contribution (see Sections 4.2 and 5).
PPT Slide
Lager Image
Examples showing how MSG masses are determined in the simulation. Plotted are the redshifts and masses of halos of which at least one galaxy is identified to belong to an MSG during the MSG search. The diamonds are the points from the modified mock catalog, and the filled, inverted triangles are the points from the true redshift catalog. The vertical dashed lines indicate the redshifts at ∆v = ±2000 km/s from the redshift of the corresponding MSG and the vertical dotted lines indicate the median photometric redshift error (Figure 4) at the search redshift of each respective MSG. The horizontal lines indicate the halo mass of 1013M and 5 × 1013M, respectively, to help readers judge the mass of subhalos. The examples are for MSGs at z = 1.20, 2.49 and 3.67. We also note the true mass of the identified MSGs on top of each panel. To derive the true MSG mass, we add up the masses of the halos with true redshifts (blue points between the dashed lines). For the observed MSG mass before finter correction, halo masses represented by the diamonds between the red-dotted lines are added.
The redshift of each MSG is chosen to be the mean spectroscopic redshift of the galaxies in the MSG, within 1 Mpc of the projected distance from the peak of the distribution. If no z spec is available, we pick the central redshift of the redshift bin where the MSG is identified.
For this procedure, note that we used a 1 Mpc radius, rather than a 1 Mpc diameter, as used for the MSG search. While a top-hat filter with a 1 Mpc diameter is efficient for identifying MSGs, MSGs are in general more extended than 1 Mpc diameter. The larger circle includes more member galaxies while minimizing the number of interlopers.
- 4.2. MSGs in GOODS Fields
In Figures 7 and 8 , we present the projected surface number density contour maps of the MSGs (left panel) and the redshift distributions of galaxies that falls within 1 Mpc radius from each MSG center (right panel). The crosses indicate the spectroscopically con- firmed galaxies within a redshift slice of | v | < 2000 km s −1 at each redshift. The shaded histograms represent the spectroscopically identified galaxies within the same redshift window. When multiple MSGs are identified at the same redshift, MSGs are categorized as “A”, “B”, and “C”, which denote the significance level in alphabetical order (A being the most significant). The numerical digits distinguish MSGs of the same significance. Overall, we identify 34 and 25 MSGs in the GOODS-North and the GOODS-South, respectively. The identified MSGs are listed in Table 4 .
PPT Slide
Lager Image
Left: MSGs in the GOODS-North field. The contours indicate the surface number density of galaxies in each redshift slice in units of number arcmin−2. Both galaxies with zspec and zphot are included to make the surface density map. The crosses represent the spectroscopically confirmed galaxies in the rest-frame redshift interval of ±2000 km s−1 at each redshift. Member galaxies of MSGs are considered as those lying within 2 Mpc diameter circle from the number density peak (the dashed circles) at the MSG redshift. Right: The redshift distribution of MSG member galaxies. Galaxies with zspec are indicated using the shaded histograms, while galaxies with zphot are indicated with the open histograms.
PPT Slide
Lager Image
(Continued) The yellow open triangles represent radio galaxies (Morrison et al. 2010) within the same redshift interval.
PPT Slide
Lager Image
(Continued)
PPT Slide
Lager Image
(Continued)
PPT Slide
Lager Image
(Continued)
PPT Slide
Lager Image
(Continued)
PPT Slide
Lager Image
(Continued)
PPT Slide
Lager Image
(Continued) Note that the blue diamonds in the surface density map of MSG at z ~ 4.05 indicate the three submm galaxies at z = 4.05 discussed in Daddi et al. (2009).
PPT Slide
Lager Image
Same as Figure 7, but the GOODS-South field starts from the third panel in this figure. The filled red stars represent AGNs (Szokoly et al. 2004) and the yellow open triangles represent radio galaxies (Mainieri et al. 2008) within the same redshift interval.
PPT Slide
Lager Image
(Continued)
PPT Slide
Lager Image
(Continued)
PPT Slide
Lager Image
(Continued)
PPT Slide
Lager Image
(Continued)
PPT Slide
Lager Image
(Continued)
PPT Slide
Lager Image
(Continued)
Properties of MSGs in the GOODS fields
PPT Slide
Lager Image
Notes: (1) MSG name based on arcmin-precision coordinates; (2) redshift of the MSG, based on spectroscopic data; (3) redshift of MSG, based on the spectroscopic members, typical spectroscopic redshift accuracy is better than 0.01; (4) overdenisty factor; (5) number of spectroscopically confirmed member galaxies within rp × 1.5 Mpc radius circle with ∆v < 2000 km s−1 from the surface number density peak of MSGs; (6) number of galaxies with photometric redshift within rp × 1.5 Mpc radius circle from the surface number density peak of MSGs; (7) sum of stellar masses ×1/finter × fMF of each MSG; (8) mass of MSG; (9) ratio of stellar mass and total mass of MSG; (10) mass of MSGs from literature; (11) projected radius of MSG. References: a: Daddi et al. (2009); b: Salimbeni et al. (2009); c: Díaz-Sáanchez et al. (2007); d: Kurk et al. (2009); e: Castellano et al. (2007); f: Kang & Im (2009)
These 59 MSGs have overdensity factors of 3.5 − 8 σ , and some of them have already been reported in the literature as clusters/proto-clusters ( Dawson et al. 2001 ; Szokoly et al. 2004 ; Vanzella et al. 2006 ; Díaz-Sánchez et al. 2007 ; Kang & Im 2009 ; Salimbeni et al. 2009 ; Kurk et al. 2009 ; Daddi et al. 2009 ; Balestra et al. 2010 ). In Figures 7 an 8 , we also mark radio galaxies (yellow open triangles), submm galaxies (blue open diamonds), and AGNs (red filled stars) if there are such objects in the redshift range. The MSGs at z ~ 0.85, 1.02 and 1.14 in the GOODS-North and z ~ 0.67, 0.74, 1.09 and 1.61 in the GOODS-South include radio galaxies at the same redshift, and the MSGs at z = 0.95, 1.23, 1.61 and z = 3.7 are associated with AGNs in the GOODS-South. However, only 20% of the identified MSGs are found to be spatially associated with known AGN/radio sources. Therefore, MSG searches radio galaxies and AGNs as signposts can miss a large number of overdense areas.
- 4.3. Mass of MSGs
We estimate MSG masses (including dark matter) by adding up the stellar mass of galaxies within 1.5× the projected harmonic mean radius ( rp ) and converting it to the total mass assuming a M -to-M halo ratio of galaxies at z < 1. The projected harmonic mean radius rp may not be the best measure of cluster radius, but it is an effective way to measure the size of gravitationally bound systems with small number of member galaxies ( Carlberg et al. 1996 ), independently of the center definition. We use 1.5 × rp since this corresponds roughly to the distance from the MSG peaks where the number density approaches the 1- σ value of the surface number density distribution.
We determine rp using galaxies within 1 Mpc radius from the surface number density peak as follows:
PPT Slide
Lager Image
Here, N is the number of objects considered, and ri is the distance to the i th galaxy from the surface density peak measured in physical scale. This gives a range of projected harmonic radii rp = 0.3 − 0.7 Mpc for the MSGs.
The stellar mass of each member galaxy within 1.5 × rp is determined by performing an SED fit using the same templates as for the photometric redshift estimates. For sources with non-detections in the NIR or IRAC bands, we use upper limits for the SED fitting and the stellar mass estimation. Note that our templates assume the Salpeter IMF, therefore, the values can shift by 0.2 − 0.24 dex if different IMFs (Kroupa IMF, Chabrier IMF) are used. Next, these stellar masses are summed, and several corrections are made to the summed stellar mass of MSG in the following order: (1) the foreground or background galaxy contamination due to uncertainty in z phot ; (2) incompleteness of the galaxy stellar mass function in the magnitude-limited sample; and (3) conversion of the stellar mass to the total MSG mass.
First, the summed stellar mass is corrected for foreground or background galaxies in the photometric redshift sample. The level of contamination is determined from the simulation (See Section 5 for details), and de- fined as the mass of MSG member candidates within the 1.5 × rp radius versus the mass of the true members of MSGs. The interloper fraction, f inter , is about a factor of 2 - 3 for most redshift intervals, except for 1.5 < z < 2. At 1.5 < z < 2, the large uncertainty in z phot and the paucity of z spec leads to an increased level of contamination by interlopers (×5, Figure 9 ).
PPT Slide
Lager Image
The correction factor finter as a function of redshift. finter is derived from simulation data and corrects for the foreground or background galaxy contamination in the zphot sample.
Then, we make a correction for the incompleteness in the coverage of the faint (less massive) end of the luminosity (mass) function, since the MSG selection is based on a magnitude-limited sample. This correction factor, f MF , is determined from the simulation data (Section 5). We estimate this value using the ratio of the total stellar mass of all galaxies that belong to an MSG to the total stellar mass of the galaxies that are brighter than the limiting magnitude. Figure 10 shows f MF as a function of redshift.
PPT Slide
Lager Image
The correction factor fMF as a function of redshift. It is derived from simulation data and corrects for the incompleteness in the stellar mass function in a magnitudelimited sample.
With the corrections of f inter and f MF , we obtain the total stellar mass of an MSG. In the next step, M ∗,total is converted to the total MSG mass, including dark-matter, by adopting the stellar mass to the observed halo mass relation of galaxies at 0 < z < 1 ( Giodini et al. 2009 ; Andreon 2010 ). According to these and other relevant works, it is known that the ratio of the stellar mass to the cluster or halo mass changes with the halo mass. This factor, f ∗,MSG , is ( Giodini et al. 2009 ):
PPT Slide
Lager Image
If M ∗,total = 10 13.7 , then 1/ f ∗,MSG = 0.05 is adopted. Here, the derivation of the stellar mass assumes the Salpeter IMF. This relation is suggested to change little as a function of redshift out to z = 4 ( Behroozi et al. 2010 ; Moster et al. 2010 ), therefore we apply the same relation to galaxies at z > 1.
We note that some of the MSGs at higher redshifts are suspected to be a collection of small mass halos. In such a case, the procedure described above could underestimate the MSG mass since the smaller halos have smaller f ∗,MSG values. If an MSG is made of three halos with an equal mass, this could lead to an underestimation of M MSG by 25%. This amount is small compared to the uncertainties related to the correction procedures (a factor of a few), therefore it can be neglected. We also note from the simulation that most of the MSGs are the dominant density peaks within the search redshift range. Finally, M MSG is derived as
PPT Slide
Lager Image
where m ∗,i is the stellar mass of the i th galaxy identi- fied as a MSG member in the magnitude-limited sample. Note that the systematic error of the stellar mass estimate due to the assumed IMF gets cancels due to the multiplication of f ∗,MSG by m ∗,i , as long as an identical IMF is assumed for deriving both the f ∗,MSG and m ∗,i values. The derived M MSG and M ∗,total are listed in Table 4 .
Errors in the MSG mass estimates are dominated by dispersion in f inter , f MF , and f ∗,MSG , and errors in the stellar mass estimates. The dispersions in f inter and f MF are inferred from the simulation as described in the next section, and they amount to 0.16 dex and 0.1 dex respectively. The dispersion in f ∗,MSG is based on Eq. (2), and is about 0.1 dex. Therefore, the combined dispersion of the three factors is 0.21 dex. Typical random errors in galaxy stellar mass estimates are a factor of 2 - 3 ( Pérez-González et al. 2008 ), but the random errors cancel out with each other when stellar masses of many galaxies are added up. More worrisome is the systematic error in stellar mass estimates which could be as large as a factor of 2 (Kannappan & Gawiser 2007; Lee et al. 2009). As we shall see in the next section, the ratio between our M MSG estimates and M MSG estimates from other groups shows a dispersion of about 0.4 dex − 0.5 dex, with a systematic difference up to 0.2 dex. If two independent MSG estimates have equal uncertainties, this implies that the error in MMSG estimates is 0.25 dex - 0.35 dex. With these considerations, we conclude that the errors in M MSG are somewhere between 0.21 dex - 0.35 dex. As a fiducial error of M MSG , we take 0.28 dex which is the mid-point of 0.21 dex to 0.35 dex. The choice of the M MSG error becomes important when determining the expected number of MSGs, and we will investigate how the result changes for different values of the M MSG error.
- 4.4. Comparison with Previously Identified (Proto-) Clusters in GOODS Fields
Previous works have identified large scale structures and clusters/proto-clusters in the GOODS fields. Most of them are included as MSGs in our study, and they are identified in Table 4 by indicating their mass in the literature in column 9. Among 16 MSGs at z < 2 in GOODS-South, 9 were previously identified ( Vanzella et al. 2006 ; Díaz-Sánchez et al. 2007 ; Castellano et al. 2007 ; Kang & Im 2009 ; Salimbeni et al. 2009 ; Kurk et al. 2009 ; Balestra et al. 2010 ). For the GOODS-North field, only a few MSGs were previously identified as large scale structures ( Dawson et al. 2001 ; Daddi et al. 2009 ).
At high redshifts, we note that an MSG at z = 4.05 in the GOODS-North field was previously identified, and MSGs at z = 2.55 and 3.70 were previously identified in the GOODS-South field ( Gilli et al. 2003 ; Szokoly et al. 2004 ; Kang & Im 2009 ; Daddi et al. 2009 ). However, many of the MSGs found in this work are new identifications.
In Figure 11 , we compare the MSG masses in the literature and our MSG mass estimates. The comparison shows that our M MSG estimates correlates reasonably well with the literature values. The mean offset between our mass estimates and the literature values is between −0.03 dex and 0.21 dex (〈log 10 [ M MSG (ours)/ M MSG (literature)]〉), depending on the largest values and the smallest values taken from the literature respectively. The rms dispersion of the relation ranges between 0.4−0.5 dex, meaning that MSG mass estimates from different methods and groups are consistent with each other within a factor of a few.
PPT Slide
Lager Image
The comparison of MSG masses from this paper with those obtained by other groups and available in the literature.
5. ANALYSIS OF MILLENNIUM SIMULATION DATA
Since the identification and mass estimates of MSGs involve several steps that cannot be directly compared to theoretical predictions of massive halos, we make use of the simulation data. We compare the observed results with the theoretical predictions by mimicking the observation and our analysis method of the GOODS data as much as possible. To do so, we adopt the mock catalogs constructed by De Lucia & Blaizot (2007) and Kitzbichler & White (2007) which are based on the Millennium Simulation. The Millennium Simulation data ( Springel et al. 2005 ) follows trajectories 1 × 10 10 particles of mass 8.6 × 10 8 h −1 M in a periodic box with 500 h −1 Mpc on a side. The mock catalog identifies galaxies and assigns luminosity values in the filter bands that overlap with the GOODS data and therefore it is well suited for this purpose. The mock catalog includes information such as redshift (plus peculiar velocity), galaxy luminosity, galaxy stellar mass, halo mass where a galaxy belongs to, and the spatial information of galaxies. To perform a realistic analysis of the simulation data, we assign photometric redshifts for some galaxies in the simulation, and use the photometric redshift information rather than their true redshift in the mock catalog during the analysis. This procedure is necessary since the GOODS data include a fair amount of photometric redshift information, and the uncertainties in the photometric redshift values diffuse MSGs in the radial direction and make the identification of MSG members less unambiguous. The quantification of what fraction of galaxies should be given photometric redshifts is done following the spectroscopic redshift completeness in the observed data as a function of the apparent Ks -band magnitude and redshift ( Figure 2 ). Following this spectroscopic redshift completeness, we assign either spectroscopic redshifts ( z spec,sim ) or photometric redshifts ( z phot,sim ) to galaxies in the simulation within a certain apparent magnitude and redshift bins as
PPT Slide
Lager Image
where δz is a value randomly picked from a Gaussian distribution with a mean of zero and a standard deviation of σ zphot as in Figure 4 .
We name the mock catalog that implements the photometric redshifts as the “modified” mock catalog, to distinguish it from the “original” mock catalog ( De Lucia & Blaizot 2007 ; Kitzbichler & White 2007 ). MSGs are searched in the modified mock catalog using the same method as we used for the analysis of the GOODS data. For the MSGs identified in the simulation, we have the true redshift distribution of galaxies and the mass of halos associated with these galaxies and MSGs. These information is then used for deriving the correction factors that are used for the analysis of the observed data.
Below, we describe how the correction factor f inter that corrects for foreground or background interlopers is derived from the simulation. The factor f inter corresponds to the ratio between the MSG mass derived from the true redshift values given in the simulation and the MSG mass estimate based on the modified mock catalog before the interloper correction. Since the photometric redshifts have large errors compared to spectroscopic redshifts, interloper correction is important for a sample with many photometric redshifts. Figure 12 shows examples of the procedure for determining MSG masses and f inter at three different redshift bins. The figure shows the redshift versus the halo mass ( M halo ) of member galaxies that are selected as possible MSG members in the modified mock catalog by imposing the same observational selection constraints as the MSG search done in the GOODS fields. Two different sets of points are plotted in the figure. In one set, true redshifts in the original mock catalog are used (the inverted blue triangles). The other set uses a combination of true and photometric redshifts as in our modified mock catalog. When a halo includes multiple galaxies and the true redshifts are available for some of these galaxies, the true redshift is used to indicate the redshift of the halo. The two vertical dotted lines indicate the redshift interval over which the MSG search was performed and the two vertical dashed lines indicate the rest-frame redshift interval of ∆ v = ±2000 km / s centered on the redshift of the MSG found from the modi- fied mock catalog. Some of the blue triangles spill over the boundary defined by the two vertical dotted lines. They are galaxies that have z phot values in the modified mock catalog, and their z phot ’s fall within the boundary of the dotted lines although their true redshifts are outside the boundary. The “true” MSG mass without interlopers is determined by adding up the masses of halos whose true redshifts (the inverted triangles) lie within the redshift interval of ∆ v < 2000 km s −1 (the vertical dashed lines) and within 1.5× rp of the peak of the projected number density. By doing so, we exclude less significant peaks in the foreground and the background. Another MSG mass is estimated by summing up the masses of halos in the search redshift range (the open diamonds), which we call, “zphot-diluted” MSG mass. Then, finter is the ratio of the “true” MSG mass to the “ z phot -diluted” MSG mass. Figure 9 shows f inter values as a function of redshift. The figure indicates that MSG masses determined from the modified mock catalog are overestimated by a factor of 2 to 3 on average, if we do not take into account f inter . The situation becomes worse (> 4× mass overestimate) if both the accuracy of z phot and the fraction of z spec sample is low such as at the redshift interval of 1.5 < z < 2.
The factor f MF is estimated from the simulation too. We estimate this value as the ratio of the total stellar mass of all galaxies that belong to an MSG to the stellar mass of galaxies associated with the MSG that are brighter than the limiting magnitude used in the observations.
Finally, we derive the dispersion in f inter and f MF values, a measure of how large the errors in MSG mass estimates are in the GOODS data. The rms dispersions are found to be 0.16 dex and 0.1 dex, respectively.
For the MSGs found in the simulation, we add a random error of 0.28 dex to the derived MSG mass in order to mimic the observation. Adding the dispersion is important since we are studying the massive end of the halo mass function, and the number of very massive MSGs is sensitive to the measurement errors. The number of MSGs can increase significantly if the errors are large, since the errors are more likely to move less massive MSGs to the more massive regime than the other way around.
6. NUMBER DENSITY OF MSGS
In order to see if the number density of MSGs in the GOODS fields agrees with the simulation data, we derive the number density evolution of MSGs as a function of redshift, and compare it with the simulation result.
The number density of MSGs is derived by summing the number of MSGs within each redshift bin with a bin size of ∆ z ≃ 0.2×(1+ z ), and then dividing the number by the surveyed comoving volume covered by the redshift bin. We estimate the number densities only over the mass range where our MSG search scheme is effective at picking up MSGs (i.e., M > a few ×10 13 M ). A certain amount of incompleteness is expected in the number density at low mass end (by missing low mass MSGs), but we do not adopt any incompleteness correction since we expect the same level of incompleteness to be present for both the observations and the simulation data. Note that MSGs at 1.5 < z < 2 are excluded in this analysis due to their large mass uncertainties that arise from inaccuracy in photometric redshifts and a small fraction of galaxies with z spec .
Figure 13 shows the number density of MSGs at different redshifts with two mass thresholds (top: 2 × 10 13 M , bottom: 7 × 10 13 M ). The red thick squares indicate the observed number densities from the GOODS fields, and the blue triangles show the numbers derived from the analysis of the simulation data. When deriving the number densities from the simulation, a random error of 0.28 dex is added to MSG mass values as described in the previous section. At the lower mass cut of M > 2 × 10 13 M (top panel of Figure 13 ), the MSG number densities from z = 0.6 to z = 4.5 are marginally consistent with the simulation within a factor of two, with more MSGs in the GOODS fields than in the simulation. The discrepancy between the observations and the simulation becomes more prominent for a higher mass cut of M > 7 × 10 13 M (bottom panel of Figure 13 ). Here, we find that the number densities of MSGs from the GOODS data are significantly higher than those from the simulation at z > 2 by a factor of ~5. There are too many MSGs in the observed data at z > 2, compared to the simulation.
PPT Slide
Lager Image
The observed number densities of MSGs as a function of redshift (red squares), compared with the simulation result (blue triangles). The upper panel shows the number density of MSGs with M > 2 × 1013 M, while the bottom panels show the number density of MSGs with M > 7 × 1013 M. The error bars represent the Poisson errors. When constructing the simulation result, a random scatter of 0.28 dex is taken into account. For the case of the lower mass threshold, the observed MSG number densities are marginally consistent with the simulation. However, the number densities of massive MSGs are significantly greater than the simulation at z > 2 when the higher mass MSGs are considered.
Since the exact value of the error in the MSG mass estimates is uncertain, we examine in Figure 14 how the result changes if we vary the error added in the MSG mass (in the simulation). In Figure 14 , the random errors in MSG mass are changed to the values of 0, 0.21 dex, 0.35 dex, and 0.5 dex. As the random error increases, the discrepancy between the observations and the simulation narrows, but the difference never decreases enough to remove the discrepancy.
PPT Slide
Lager Image
The number densities of MSGs with assumed scatters of 0, 0.21 dex, 0.35 dex, and 0.50 dex, from top to bottom panel, respectively, for the simulation results. The red squares and blue triangles indicate MSGs having a mass exceeding 2 × 1013M (left) and 7 × 1013M (right) in the GOODS fields and in the simulation data, respectively.
We also investigate how systematic errors in the MSG mass estimates affect the result. Figure 15 shows the results, where the systematic offset in the MSG mass measurements are varied from 0.1 dex to 0.2 dex. These systematic offsets are applied to the masses of MSGs from the observed data (i.e., MSG masses are reduced by this amount). As in the case for Figure 13 , the random errors of 0.28 dex are assumed for the simulation. When the observed MSG masses are systematically reduced by 0.2 dex, we find that the discrepancy in the MSG number densities between the observation and the simulation disappears. However, we note that the systematic error we quoted earlier, when deriving the stellar masses through SED-fitting, does not affect the derived MSG mass. The conversion of the stellar mass to the halo mass also includes the same kind of systematics in the stellar mass, thus the systematic error cancels out when deriving the MSG mass.
PPT Slide
Lager Image
The number densities of MSGs with the systematic offsets with the reduction of the MSG masses from the observed data by 0.1 dex (upper panel) and 0.2 dex (bottom panel). A scatter of 0.28 dex is assumed for the simulation. The red squares and blue triangles indicate MSGs having a mass exceeding 2 × 1013M (left) and 7 × 1013M (right) in the GOODS fields and the simulation, respectively.
Overall, we conclude that there is an over-abundance of massive MSGs with M > 7 × 10 13 M at z > 2 by a factor of a few, compared with the Millennium Simulation. One might be able to reconcile the observation and the simulation if there is a systematic offset in the derived MSG mass by 0.2 dex, but it is not clear how such a systematic bias could arise.
7. DISCUSSION
- 7.1. Cosmic Variance
We consider here the possibility that our finding of too many MSGs is due to cosmic variance effect. If the GOODS fields have an anomalously large number of MSGs due to cosmic variance, then, the difference between the simulation and the observation can be alleviated. The moderate field of view of 320 arcmin 2 covers a large volume (~ 10 5 Mpc 3 ) when the redshift bins are of order ∆ z ~ 0.1, which is sufficient to negate the cosmic variation. Nevertheless, we can estimate from the simulation the chance coincidence that a survey area as small as the GOODS fields could have anomalously large number of MSGs. To do so, we divide the simulation area into 2500 sets of two 160 arcmin 2 , and investigate how many such sets can have the MSG number densities as large as our result. No subdivided regions reveal as many MSGs as what we found from the GOODS field, giving a probability of < 1/2500 for the chance coincidence that the GOODS fields happened be highly overdense areas.
- 7.2. Implications for Initial Density Fluctuations
Here, we discuss the implication of our results on the initial density fluctuations. We focus on σ 8 and ns , and the non-Gaussianity. The initial density fluctuation assumed in the Millennium Simulation is not up-to-date with the more recent observational constraints. For σ 8 at a scale of r = 8 h −1 Mpc and the primordial powerlaw index ns , they adopted the values of σ 8 = 0.9 and ns = 1 while the latest measurements from the WMAP 7-year data suggest the values of σ 8 = 0.816±0.024 and ns = 0.968±0.012, respectively ( Jarosik et al. 2011 ; Komatsu et al. 2011 ). The adoption of σ 8 = 0.9 and ns = 1 introduces a factor of ~ 2 − 3 increase of the halo number density in the Millennium Simulation compared to the halo number density from the Sheth-Tormen halo mass function with the WMAP values of σ 8 and ns ( Figure 16 ). This shows that the expected number density of MSGs in the simulation would become even less than what is found in the Millennium Simulation data, if we were to adopt the WMAP σ 8 and ns values.
PPT Slide
Lager Image
The ratio of the abundance of halos as a function of redshift for the cosmological parameters of the WMAP 7-year and of the Millennium Simulation with three different mass cuts. To create this figure, the Sheth-Tormen halo mass function is convolved with a lognormal distribution function in mass with a standard deviation of 0.28 dex.
Next, we consider non-Gaussianity in the initial density fluctuation as one possibility, since the abundance of massive halos at high redshifts is also sensitive to non-Gaussianity of initial density fluctuation spectrum. Recently, there is interest in the non-Gaussianity from a theoretical standpoint ( D’Amico et al. 2011 ; LoVerde & Smith 2011 ; Smith et al. 2012 ). Deviations from Gaussian initial conditions are parameterized in terms of the dimensionless fNL parameter ( Salopek & Bond 1990 ; Gangui et al. 1994 ; Verde et al. 2000 ; Komatsu & Spergel 2001 ). If the initial density fluctuation has a skewed distribution in which high density peaks are more abundant than what is expected in a Gaussian density fluctuation, one would expect that more massive halos can be more easily found early in the Universe. Such a case corresponds to a large positive value of fNL in the non-Gaussian framework. Adopting a large non-linearity parameter of fNL ~ +500 or ~ −500, one can boost or decrease the number density of massive halos by a factor of few at z > 2 ( Grossi et al. 2007 ). Therefore, our result may imply that the fNL value could be as large as +500. On the other hand, several recent works provide constraints on fNL to be about ~ 30 ( Smith et al. 2009 ; Curto et al. 2009 ; Komatsu 2010 ). If so, the number density of massive halos does not change more than a factor of 1.1, a negligible amount. Clearly, while further investigation is needed, but it seems that the non-Gaussianity cannot explain our result, assuming that the recent limits of fNL ~ 30 or less are valid.
- 7.3. Constraints on Hot Dark Matter
We discuss the implication of our result for hot dark matter. The streaming motion of hot dark matter such as massive neutrinos smooths out the density fluctuation, and suppresses the formation of massive structures leading to a reduction of the number of massive structures at high redshift. Therefore, massive hot dark matter particles are likely to reduce the number of MSGs in the simulation, rather than helping explain our result. Marulli et al. (2011) predict that the number density of halos with M = 10 14 M , in comparison to the standard Sheth-Tormen halo mass function, is reduced by a factor of ~ 10 and ~ 3 for Mν = 0.6 eV, and Mν = 0.3 eV respectively, where Mν represents the total neutrino mass of different species. If the neutrinos have masses as heavy as Mν ~ 0.3 eV, the discrepancy between the observed MSG number density and the number density in the simulation will become even more severe. The current experimental lower limit is Mν ~ 0.05 eV ( Lesgourgues & Pastor 2006 ), and our result suggests that Mν is not likely to be much higher than the lower limit.
- 7.4. Why Are There Too Many MSGs atz> 2?
The discussions above narrow down the possible causes for the over-abundance of MSGs in the GOODS fields to three remaining possibilities: (1) some baryonic physics that are not implemented well in the simulation; (2) systematic overestimate of MSG mass; and (3) the assumption of the ΛCDM cosmology in the simulation.
First, we consider a case that unknown baryonic physics in galaxy formation leads to play a role in the mismatch of the MSGs number density between the observation results and the simulation. Our search and mass estimation of MSGs rely on galaxies, therefore it is important that the simulation reproduces the observed galaxy properties well. The Millennium Simulation has been tested extensively against observational data, and in many cases, excellent agreement has been reported between the observation results and the simulation (e.g., Scoville et al. 2013 ). However, it is noteworthy that the Millennium Simulation has dif- ficulty explaining some of the galaxy properties such as the abundance of massive galaxies at high redshifts ( Collins et al. 2009 ), and a certain aspect of the SFRdensity relation ( Elbaz et al. 2007 ). A more comprehensive comparison of the observed data and simulated data, including spectroscopic data and newer simulations, should help us understand what kind of baryonic physics has been missed in the simulation.
Second, as we have shown earlier, systematic overestimation of MSG masses by a factor of 0.2 dex can explain why the MSG number density appears to be too high. This is perhaps the most likely explanation of our results, although quantification of this effect is difficult. Independent methods to estimate MSG masses would provide a useful consistency check of the MSG mass. We hope to produce such a study in the future, including more spectroscopic data for our MSG member candidates.
Finally, if all of the above possibilities fail to explain the result, the over-abundance of MSGs at high redshifts could present a problem for the ΛCDM cosmology. A different simulation may reveal that this is not the case, but at this moment, we suggest that the ΛCDM cosmology may be in trouble if the above solutions do not succeed to explain our results.
8. CONCLUSIONS
Using the multi-wavelength data, we identified 59 MSGs from the combined areas of the GOODS-South and the GOODS-North fields with significances of 3.5− 8 σ from z = 0.6 to z ~ 4.5. Among them, ~ 20% of MSGs show plausible associations with AGN/radio sources.
In comparison with a simulation data set, we find a discrepancy between the observed number densities of MSGs and those from the simulation at z > 1 ( M > 7 × 10 13 M ). The discrepancy becomes more significant at higher redshifts ( z > 2) by a factor of ~ 5 or more. Even after considering possible systematic effects, our result implies that there are too many massive structures at z > 2 compared to the ΛCDM prediction. By tweaking the conditions for the initial density fluctuation or baryonic physics in galaxy formation, one may be able to explain the result within the ΛCDM cosmology framework, but as of now the overabundance of MSGs at z > 2 stands as a challenge to the models based on the ΛCDM cosmology.
Acknowledgements
This work was supported by the Creative Initiative program, No. 2008-0060544, of the National Research Foundation of Korea (NRFK) funded by the Korean government (MSIP). We thank Ranga-Ram Chary for providing us with the IRAC catalog, and Lihwai Lin for the WIRCAM images of the GOODS-N field that were useful during the initial stage of the analysis. We thank Giulia Rodighiero, Ho Seong Hwang, Jae-Woo Kim, Juhan Kim, Marios Karouzos, and Rafael Gobat for useful communication.
References
Allen S. W. , Evrard A. E. , Mantz A. B. 2011 Cosmological Parameters from Observations of Galaxy Clusters ARAA 49 409 -    DOI : 10.1146/annurev-astro-081710-102514
Andreon S. 2010 The Stellar Mass Fraction and Baryon Content of Galaxy Clusters and Groups MNRAS 407 263 -    DOI : 10.1111/j.1365-2966.2010.16856.x
Bahcall N. A. , Ostriker J. P. , Perlmutter S. , Steinhardt P. J. 1999 The Cosmic Triangle: Revealing the State of the Universe Science 284 1481 -    DOI : 10.1126/science.284.5419.1481
Baldi M. 2012 Early Massive Clusters and the Bouncing Coupled Dark Energy MNRAS 420 430 -    DOI : 10.1111/j.1365-2966.2011.20048.x
Balestra I. 2010 The Great Observatories Origins Deep Survey. VLT/VIMOS Spectroscopy in the GOODS-South Field: Part II A&A 512 A12 - 12
Barger A. J. , Cowie L. L. , Wang W.-H. 2008 A Highly Complete Spectroscopic Survey of the GOODS-N Field1 ApJ 689 687 -    DOI : 10.1086/592735
Bayliss M. B. , Ashby M. L. N. , Ruel J. 2014 SPT-CLJ2040-4451: An SZ-Selected Galaxy Cluster at z = 478 With Significant Ongoing Star Formation ApJ 794 12 -    DOI : 10.1088/0004-637X/794/1/12
Behroozi P. S. , Conroy C. , Wechsler R. H. 2010 A Comprehensive Analysis of Uncertainties Affecting the Stellar Mass-Halo Mass Relation for 0    DOI : 10.1088/0004-637X/717/1/379
Benson B. A. , de Haan T. , Dudley J. P. 2013 Cosmological Constraints from Sunyaev-Zel’dovich-Selected Clusters with X-Ray Observations in the First 178 deg2 of the South Pole Telescope Survey ApJ 763 147 -    DOI : 10.1088/0004-637X/763/2/147
Benítez N. 2000 Bayesian Photometric Redshift Estimation ApJ 536 571 -    DOI : 10.1086/308947
Bertin E. , Arnouts S. 1996 SExtractor: Software for Source Extraction A&A Suppl. 117 393 -
Bielby R. M. , Finoguenov A. , Tanaka M. 2010 The WIRCAM Deep Infrared Cluster Survey. I. Groups and Clusters at z ≳ 1.1 A&A 523 A66 -
Biviano A. , Murante G. , Borgani S. , Diaferio A. , Dolag K. , Girardi M. 2006 On the Efficiency and Reliability of Cluster mass Estimates Based on Member Galaxies A&A 456 23 -
Bond J. R. , Efstathiou G. , Silk J. 1980 Massive Neutrinos and the Large-Scale Structure of the Universe Physical Review Letters 45 1980 -    DOI : 10.1103/PhysRevLett.45.1980
Bruzual G. , Charlot S. 2003 Stellar Population Synthesis at the Resolution of 2003 MNRAS 344 1000 -    DOI : 10.1046/j.1365-8711.2003.06897.x
Carlberg R. G. , Yee H. K. C. , Ellingson E. 1996 Galaxy Cluster Virial Masses and Ω ApJ 462 32 -    DOI : 10.1086/177125
Capak P. 2004 A Deep Wide-Field, Optical, and Near-Infrared Catalog of a Large Area around the Hubble Deep Field North AJ 127 180 -    DOI : 10.1086/380611
Capak P. L. , Riechers D. , Scoville N. Z. 2011 A Massive Protocluster of Galaxies at a Redshift of z ∼ 5.3 Nature 470 233 -    DOI : 10.1038/nature09681
Carlesi E. , Knebe A. , Yepes G. 2012 N-Body Simulations with a Cosmic Vector for Dark Energy MNRAS 424 699 -    DOI : 10.1111/j.1365-2966.2012.21258.x
Castellano M. 2007 A Photometrically Detected Forming Cluster of Galaxies at Redshift 1.6 in the GOODS Field ApJ 671 1497 -    DOI : 10.1086/521595
Cayón L. , Gordon C. , Silk J. 2011 Probability of the Most Massive Cluster under Non-Gaussian Initial Conditions MNRAS 415 849 -    DOI : 10.1111/j.1365-2966.2011.18770.x
Charlot S. , Fall S. M. 2000 A Simple Model for the Absorption of Starlight by Dust in Galaxies ApJ 539 718 -    DOI : 10.1086/309250
Chae K.-H. , Biggs A. D. , Blandford R. D. 2002 Constraints on Cosmological Parameters from the Analysis of the Cosmic Lens All Sky Survey Radio-Selected Gravitational Lens Statistics Physical Review Letters 89 151301 -    DOI : 10.1103/PhysRevLett.89.151301
Chiba M. , Yoshii Y. 1997 Do Lensing Statistics Rule Out a Cosmological Constant? ApJ 489 485 -    DOI : 10.1086/304831
Colless M. 2001 The 2dF Galaxy Redshift Survey: Spectra and Redshifts MNRAS 328 1039 -    DOI : 10.1046/j.1365-8711.2001.04902.x
Collins C. A. , Stott J. P. , Hilton M. 2009 Early Assembly of the Most Massive Galaxies Nature 458 7865 -
Cooper M. C. 2011 The DEEP3 Galaxy Redshift Survey: Keck/DEIMOS Spectroscopy in the GOODS-N Field ApJS 193 14 -    DOI : 10.1088/0067-0049/193/1/14
Coupon J. , Kilbinger M. , McCracken H. J. 2012 Galaxy Clustering in the CFHTLS-Wide: the Changing Relationship between Galaxies and Haloes Since z ∼ 1.2 A&A 542 A5 -
Cucciati O. , Zamorani G. , Lemaux B. C. 2014 Discovery of a Rich Proto-Cluster at z = 2.9 and Associated Diffuse Cold Gas in the VIMOS Ultra-Deep Survey(VUDS) A&A 570 16 -
Curto A. , Martínez-González E. , Mukherjee P. 2009 Wilkinson Microwave Anisotropy Probe 5-yr Constraints on fnl with Wavelets MNRAS 393 615 -    DOI : 10.1111/j.1365-2966.2008.14230.x
Daddi E. , Dannerbauer H. , Stern D. 2009 Two Bright Submillimeter Galaxies in a z = 4.05 Protocluster in Goods-North, and Accurate Radio-Infrared Photometric Redshifts ApJ 694 1517 -    DOI : 10.1088/0004-637X/694/2/1517
D’Amico G. , Musso M. , Norena J. , Paranjape A. 2011 An Improved Calculation of the Non-Gaussian Halo Mass Function JCAP 2 1 -
Dawson S. 2001 Serendipitously Detected Galaxies in the Hubble Deep Field AJ 122 598 -    DOI : 10.1086/321160
De Lucia G. , Blaizot J. 2007 The Hierarchical Formation of the Brightest Cluster Galaxies MNRAS 375 2 -    DOI : 10.1111/j.1365-2966.2006.11287.x
Díaz-Sánchez A. , Villo-Pérez I. , Pérez-Garrido A. , Rebolo R. 2007 A Low Mass Cluster of Extremely RedGalaxies at z = 1.10 in the GOODS Southern Field MNRAS 377 516 -    DOI : 10.1111/j.1365-2966.2007.11643.x
Doherty M. 2010 Optical and Near-IR Spectroscopy of Candidate Red Galaxies in Two z ∼ 2.5 Proto-Clusters A&A 509 A83 -
Eisenstein D. J. , Zehavi I. , Hogg D. W. 2005 Detection of the Baryon Acoustic Peak in the Large-Scale Correlation Function of SDSS Luminous Red Galaxies ApJ 633 560 -    DOI : 10.1086/466512
Elbaz D. , Daddi E. , Le Borgne D. 2007 The Reversal of the Star Formation-Density Relation in the Distance Universe A&A 468 33 -
Fassbender R. , Nastasi A. , Böhringer H. 2011 The X-Ray Luminous Galaxy Cluster XMMU J1007.4+1237 at z = 1.56. The Dawn of Starburst Activity in Cluster Cores A&A 527 L10 -
Gangui A. , Lucchin F. , Matarrese S. , Mollerach S. 1994 The Three-Point Correlation Function of the Cosmic Microwave Background in Inflationary Models ApJ 430 447 -    DOI : 10.1086/174421
Giavalisco M. 2004 The Great Observatories Origins Deep Survey: Initial Results from Optical and Near-Infrared Imaging ApJL 600 L93 -    DOI : 10.1086/379232
Gilbank D. G. , Gladders M. D. , Yee H. K. C. , Hsieh B. C. 2011 The Red-Sequence Cluster Survey-2 (RCS-2): Survey Details and Photometric Catalog Construction AJ 141 94 -    DOI : 10.1088/0004-6256/141/3/94
Gilli R. 2003 Tracing the Large-Scale Structure in the Chandra Deep Field South ApJ 592 721 -    DOI : 10.1086/375777
Giodini S. 2009 Stellar and Total Baryon Mass Fractions in Groups and Clusters Since Redshift 1 ApJ 703 982 -    DOI : 10.1088/0004-637X/703/1/982
Gladders M. D. , Yee H. K. C. 2005 The Red-Sequence Cluster Survey. I. The Survey and Cluster Catalogs for Patches RCS 0926+37 and RCS 1327+29 ApJS 157 1 -    DOI : 10.1086/427327
Gonzalez A. H. , Stanford S. A. , Browdin M. 2012 IDCS J1426.5+3508: Cosmological Implications of a Massive, Strong Lensing Cluster at z = 1.75 ApJ 753 163 -    DOI : 10.1088/0004-637X/753/2/163
Goto T. , Hanami H. , Im M. 2008 Galaxy Clusters at 0.9 < z < 1.7 in the AKARI NEP Deep Field PASJ 60 531 -
Grossi M. , Dolag K. , Branchini E. , Matarrese S. , Moscardini L. 2007 Evolution of Massive Haloes in Non-Gaussian Scenarios MNRAS 382 1261 -    DOI : 10.1111/j.1365-2966.2007.12458.x
Haiman Z. , Mohr J. J. , Holder G. P. 2001 Constraints on Cosmological Parameters from Future Galaxy Cluster Surveys ApJ 553 545 -    DOI : 10.1086/320939
Harrison I. , Coles P. 2012 Testing Cosmology with Extreme Galaxy Cluster MNRAS 421 L19 -    DOI : 10.1111/j.1745-3933.2011.01198.x
Harrison I. , Hotchkiss S. 2013 A Consistent Approach to Falsifying ΛCDM with Rare Galaxy Clusters JCAP 7 22 -
Hilton M. 2007 The XMM Cluster Survey: The Dynamical State of XMMXCS J2215.9-1738 at z = 1.457 ApJ 670 1000 -    DOI : 10.1086/522827
Holz D. E. , Perlmutter S. 2012 The Most Massive Objects in the Universe ApJL 755 L36 -    DOI : 10.1088/2041-8205/755/2/L36
Hotchkiss S. 2011 Quantifying the Rareness of Extreme Galaxy Clusters JCAP 7 4 -
Im M. , Simard L. , Faber S. M. 2002 The DEEP Groth Strip Survey. X. Number Density and LuminosityFunction of Field E/S0 Galaxies at z < 1 ApJ 571 136 -    DOI : 10.1086/339854
Im M. , Griffiths R. E. , Ratnatunga K. U. 1997 A Measurement of the Cosmological Constant Using Elliptical Galaxies as Strong Gravitational Lenses ApJ 475 457 -    DOI : 10.1086/303583
Jarosik N. 2011 Seven-year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Sky Maps, Systematic Errors, and Basic Results ApJS 192 14 -    DOI : 10.1088/0067-0049/192/2/14
Jee M. J. , Rosati P. , Ford H. C. 2009 Hubble Space Telescope Weak-Lensing Study of the Galaxy Cluster XMMU J2235.3 - 2557 at z ∼ 1.4: A Surprisingly Massive Galaxy Cluster When the Universe is One-third of its Current Age ApJ 704 672 -    DOI : 10.1088/0004-637X/704/1/672
Kajisawa M. 2006 Protoclusters with Evolved Populations around Radio Galaxies at z ∼ 2.5 MNRAS 371 577 -    DOI : 10.1111/j.1365-2966.2006.10704.x
Kang E. , Im M. 2009 Overdensities of Galaxies at z ∼3.7 in Chandra Deep Field-South ApJL 691 L33 -    DOI : 10.1088/0004-637X/691/1/L33
Kitzbichler M. G. , White S. D. M. 2007 The High Redshift Galaxy Population in Hierarchical Galaxy Formation Models MNRAS 376 2 -    DOI : 10.1111/j.1365-2966.2007.11458.x
Kravtsov A. V. , Borgani S. 2012 Formation of Galaxy Clusters ARAA 50 353 -    DOI : 10.1146/annurev-astro-081811-125502
Kodama T. , Tanaka I. , Kajisawa M. , Kurk J. , Venemans B. , De Breuck C. , Vernet J. , Lidman C. 2007 The First Appearance of the Red Sequence of Galaxies in Proto-Clusters at 2 ≲ z ≲ 3 MNRAS 377 1717 -    DOI : 10.1111/j.1365-2966.2007.11739.x
Komatsu E. , Spergel D. N. 2001 Acoustic Signatures in the Primary Microwave Background Bispectrum PRD 63 063002 -    DOI : 10.1103/PhysRevD.63.063002
Komatsu E. 2010 Hunting for Primordial Non-Gaussianity in the Cosmic Microwave Background Classical and Quantum Gravity 27 124010 -    DOI : 10.1088/0264-9381/27/12/124010
Komatsu E. 2011 Seven-Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Cosmological Interpretation ApJS 192 18 -    DOI : 10.1088/0067-0049/192/2/18
Kuiper E. , Hatch N. A. , Venemans B. P 2011 Discovery of a High-z Protocluster with Tunable Filters: the Case of 6C0140+326 at z = 4.4 MNRAS 417 1088 -    DOI : 10.1111/j.1365-2966.2011.19324.x
Kurk J. D. 2000 A Search for Clusters at High Redshift. I. Candidate Lyalpha Emitters near 1138-262 at z = 2.2 A&A 358 L1 -
Kurk J. 2009 GMASS Ultradeep Spectroscopy of Galaxies at z ∼ 2. V. Witnessing the Assembly at z = 1.6 of a Galaxy Cluster A&A 504 331 -
Le Fevre O. , Deltorn J. M. , Crampton D. , Dickinson M. 1996 Clustering around the Radio Galaxy MRC 0316-257 at Z = 3.14 ApJL 471 L11 -    DOI : 10.1086/310319
Lesgourgues J. , Pastor S. 2006 Massive Neutrinos and Cosmology PRep 429 307 -
Li I. H. , Yee H. K. C. , Hsieh B. C. , Gladders M. 2012 Evolution of Group Galaxies from the First Red-Sequence Cluster Survey ApJ 749 150 -    DOI : 10.1088/0004-637X/749/2/150
LoVerde M. , Smith K. M. 2011 The Non-Gaussian Halo Mass Function with fNL, gNL and τNL JCAP 8 3 -
Mainieri V. 2008 The VLA Survey of the Chandra Deep Field-South. II. Identification and Host Galaxy Properties of Submillijansky Sources ApJS 179 95 -    DOI : 10.1086/591053
Mantz A. , Allen S. W. , Rapetti D. , Ebeling H. 2010 The Observed Growth of Massive Galaxy clusters - I. Statistical Methods and Cosmological Constraints MNRAS 406 1759 -
Marulli F. , Carbone C. , Viel M. , Moscardini L. , Cimatti A. 2011 Effects of Massive Neutrinos on the Large-Scale Structure of the Universe MNRAS 418 346 -    DOI : 10.1111/j.1365-2966.2011.19488.x
McCarthy P. J. , Yan H. , Abraham R. G. 2007 A Compact Cluster of Massive Red Galaxies at a Redshift of 1.5 ApJL 664 L17 -    DOI : 10.1086/520575
Menanteau F. , Hughes J. P. , Sifon C. 2012 The Atacama Cosmology Telescope: ACT-CL J0102-4915 “El Gordo,” a Massive Merging Cluster at Redshift 0.87 ApJ 748 7 -    DOI : 10.1088/0004-637X/748/1/7
Mortonson M. J. , Hu W. , Huterer D. 2011 Simultaneous Falsification of CDM and Quintessence with Massive, Distant Clusters Physical Review D 83 02315 -
Moster B. P. , Somerville R. S. , Maulbetsch C. 2010 Constraints on the Relationship between Stellar Mass and Halo Mass at Low and High Redshift ApJ 710 903 -    DOI : 10.1088/0004-637X/710/2/903
Morrison G. E. , Owen F. N. , Dickinson M. , Ivison R. J. , Ibar E. 2010 Very Large Array 1.4 GHz Observations of the GOODS-North Field: Data Reduction and Analysis ApJS 188 178 -    DOI : 10.1088/0067-0049/188/1/178
Muzzin A. V. 2008 An Infrared Study of Distant Galaxy Clusters, Ph.D. Thesis
Muzzin A. , Wilson G. , Yee H. K. C. 2009 Two Massive Red-Sequence-Selected Galaxy Clusters at z ∼1.2 in the SpARCS-North Cluster Survey ApJ 698 1934 -    DOI : 10.1088/0004-637X/698/2/1934
Nastasi A. , Fassbender R. , Böhringer H. 2011 Discovery of the X-Ray Selected Galaxy Cluster XMMU J0338.8+0021 at z = 1.49. Indications of a Young System with a Brightest Galaxy in Formation A&A 532 L6 -
Nonino M. 2009 Deep U Band and R Imaging of Goods-South: Observations, Data Reduction and First Results ApJS 183 244 -    DOI : 10.1088/0067-0049/183/2/244
Papovich C. , Momcheva I. , Willmer C. N. A. 2010 A Spitzer-Selected Galaxy Cluster at z = 1.62 ApJ 716 1503 -    DOI : 10.1088/0004-637X/716/2/1503
Park C. , Choi Y.-Y. , Kim J. 2012 The Challenge of the Largest Structures in the Universe to Cosmology ApJL 759 L7 -    DOI : 10.1088/2041-8205/759/1/L7
Pentericci L. , Kurk J. D. , Röttgering H. J. A. 2000 A Search for Clusters at High Redshift. II. A Proto Cluster around a Radio Galaxy at z = 2.16 A&A 361 L25 -
Pérez-González P. G. , Rieke G. H. , Villar V. 2008 The Stellar Mass Assembly of Galaxies from z = 0 to z = 4: Analysis of a Sample Selected in the Rest-Frame Near-Infrared with Spitzer ApJ 675 234 -    DOI : 10.1086/523690
Perlmutter S. , Aldering G. , Goldhaber G. 1999 Measurements of Ω and Λ from 42 High-Redshift Supernovae ApJ 517 565 -    DOI : 10.1086/307221
Popesso P. 2009 The Great Observatories Origins Deep Survey. VLT/VIMOS Spectroscopy in the GOODS South Field A&A 494 443 -
Reed D. S. , Bower R. , Frenk C. S. 2007 The Halo Mass Function from the Dark Ages through the Present Day MNRAS 374 2 -    DOI : 10.1111/j.1365-2966.2006.11204.x
Retzlaff J. , Rosati P. , Dickinson M. , Vandame B. , Rité C. , Nonino M. , Cesarsky C. 2010 The Great Observatories Origins Deep Survey. VLT/ISAAC Near-Infrared Imaging of the GOODS-SouthField A&A 511 A50 -
Riess A. G. , Filippenko A. V. , Challis P. 1998 Observational Evidence from Supernovae for an Accelerating Universe and a Cosmological Constant AJ 116 1009 -    DOI : 10.1086/300499
Salimbeni S. 2009 A Comprehensive Study of Large-Scale Structures in the GOODS-SOUTH Field up to z ∼2.5 A&A 501 865 -
Salopek D. S. , Bond J. R. 1990 Nonlinear Evolution of Long-Wavelength Metric Fluctuations in Inflationary Models PRD 42 3936 -    DOI : 10.1103/PhysRevD.42.3936
Santini P. 2009 Star Formation and Mass Assembly in High Redshift Galaxies A&A 504 751 -
Scoville N. , Arnouts S. , Aussel H. 2013 Evolution of Galaxies and Their Environments at z = 0.1 − 3 in Cosmos ApJS 206 3 -    DOI : 10.1088/0067-0049/206/1/3
Seljak U. , Makarov A. , McDonald P. 2005 Cosmological Parameter Analysis Including SDSS Ly𝛼 Forest and Galaxy Bias: Constraints on the Primordial Spectrum of Fluctuations, Neutrino Mass, and Dark Energy PRD 71 103515 -    DOI : 10.1103/PhysRevD.71.103515
Smith K. M. , Ferraro S. , LoVerde M. 2012 Halo Clustering and gNL-type Primordial Non-Gaussianity JCAP 3 32 -
Smith K. M. , Senatore L. , Zaldarriaga M. 2009 Optimal Limits on fNL local from WMAP 5-Year Data JCAP 9 6 -
Spergel D. N. 2003 First-Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Determination of Cosmological Parameters ApJS 148 175 -    DOI : 10.1086/377226
Springel V. , White S. D. M. , Jenkins A. 2005 Simulations of the Formation, Evolution and Clustering of Galaxies and Quasars Nature 435 629 -    DOI : 10.1038/nature03597
Stalder B. , Ruel J. , Suhada R. 2013 SPT-CLJ0205-5829: A z = 1.32 Evolved Massive Galaxy Cluster in the South Pole Telescope Sunyaev-Zel’dovich Effect Survey ApJ 763 93 -    DOI : 10.1088/0004-637X/763/2/93
Stanford S. A. 2006 The XMM Cluster Survey: A Massive Galaxy Cluster at z = 1.45 ApJL 646 L13 -    DOI : 10.1086/506449
Szokoly G. P. 2004 The Chandra Deep Field-South: Optical Spectroscopy ApJS 155 271 -    DOI : 10.1086/424707
Toshikawa J. , Kashikawa N. , Ota K. 2012 Discovery of a Protocluster at z ∼ 6 ApJ 750 137 -    DOI : 10.1088/0004-637X/750/2/137
Vanzella E. 2005 The Great Observatories Origins Deep Survey. VLT/FORS2 Spectroscopy in the GOODS-SouthField A&A 434 53 -
Vanzella E. 2006 The Great Observatories Origins Deep Survey. VLT/FORS2 Spectroscopy in the GOODS-SouthField: Part II A&A 454 423 -
Vanzella E. 2008 The Great Observatories Origins Deep Survey. VLT/FORS2 Spectroscopy in the GOODS-SouthField: Part III A&A 478 83 -
Vanzella E. 2009 Spectroscopic Observations of Lyman Break Galaxies at Redshifts 4, 5, and 6 in the Goods-South Field A&A 695 1163 -
Venemans B. P. 2007 Protoclusters Associated with z > 2 Radio Galaxies . I. Characteristics of High Redshift Protoclusters A&A 461 823 -
Verde L. , Wang L. , Heavens A. F. , Kamionkowski M. 2000 Large-Scale Structure, the Cosmic Microwave Background and Primordial Non-Gaussianity MNRAS 313 141 -    DOI : 10.1046/j.1365-8711.2000.03191.x
Waizmann J.-C. , Ettori S. , Moscardini L. 2012 An Application of Extreme Value Statistics to the Most Massive Galaxy Clusters at Low and High Redshifts MNRAS 420 1754 -    DOI : 10.1111/j.1365-2966.2011.20171.x
Waizmann J.-C. , Ettori S. , Bartelmann M. 2013 Order Statistics Applied to the Most Massive and Most Distant Galaxy Clusters MNRAS 432 914 -    DOI : 10.1093/mnras/stt460
Wang W.-H. , Cowie L. L. , Barger A. J. , Keenan R. C. , Ting H.-C. 2010 Ultradeep KS Imaging in the GOODS-N ApJS 187 251 -    DOI : 10.1088/0067-0049/187/1/251
Williamson R. , Benson B. A. , High F. W. 2011 A Sunyaev-Zel’dovich-Selected Sample of the Most Massive Galaxy Clusters in the 2500 deg2 South Pole Telescope Survey ApJ 738 139 -    DOI : 10.1088/0004-637X/738/2/139
Wold M. 2003 Overdensities of Extremely Red Objects in the Fields of High-Redshift Radio-Loud Quasars AJ 126 1776 -    DOI : 10.1086/378362
Wuyts S. , Labbé I. , Schreiber N. M. F. 2008 FIRE-WORKS U38-to-24 μm Photometry of the GOODS Chandra Deep Field-South: Multiwavelength Catalog and Total Infrared Properties of Distant Ks-Selected Galaxies ApJ 682 985 -    DOI : 10.1086/588749