The Environmental Dependence of Galaxy Properties at z=2
Key Words.:Galaxies : clusters : individual : PKS1138-26, Galaxies : formation, Galaxies : fundamental parameters
Context: We report on the environmental dependence of galaxy properties at . We construct multi-band photometric data sets in the (proto-)cluster PKS1138-26 field and in the GOODS field. We then fit spectral energy distributions of the galaxies with model templates generated with the latest stellar population synthesis code and derive physical properties of galaxies from the fits. To quantify the environmental dependence of galaxy properties, a special care is taken of systematic errors – we use data sets that have almost the same wavelength samplings, use the same code to fit SEDs with the same set of templates, and compare relative differences between the two samples. We find that the PKS1138 galaxies have similar ages, shorter star formation time scales, lower star formation rates, and weaker dust extinction compared to the GOODS galaxies at . This trend is similar to that observed locally, suggesting that the environmental dependence of galaxy properties is already partly in place as early as . We show that the PKS1138 galaxies assemble the bulk of their masses Gyr earlier than field galaxies, i.e., the galaxy formation depends on environment. Galaxy mergers should frequently occur during the first collapse of clusters and they might play an important role in driving the observed environmental dependence of galaxy properties at .
The formation and evolution of galaxies in the Universe are dependent on environment in which galaxies live. In the local Universe, red early-type galaxies are the dominant population in rich galaxy clusters, while blue late-type galaxies are the dominant population in the low-density field. Not only galaxy properties, but the formation epoch of galaxies also depends on environment in the sense that cluster galaxies form earlier than field galaxies (e.g., Kuntschner et al. 2002; Gebhardt et al. 2003; Thomas et al. 2005). Earlier studies concentrated on nearby galaxies, but with the recent advent of large telescopes, environment studies at became possible. Interestingly, the environmental dependence of galaxy properties observed at is already strong; clusters at are dominated by red early-type galaxies (e.g., Blakeslee et al. 2003; Nakata et al. 2005; Postman et al. 2005; Lidman et al. 2008; Mei et al. 2009, but see also Cucciati et al. 2006). Also, the formation epoch of cluster galaxies measured at is consistent with that observed locally (Gobat et al. 2008) . Although there is a clear sign of galaxy evolution between and 0 (e.g., Elbaz et al. 2007; Cooper et al. 2008), one has to observe galaxies at even higher redshifts to fully quantify the environmental dependence of galaxy formation and evolution.
Because of observational difficulties, only a few high redshift clusters are known so far, the highest redshift cluster being at (Stanford et al. 2006). Higher redshift galaxies appear fainter and their rest-frame optical light migrates to the near-IR, where the sky background is brighter and it is challenging to observe faint galaxies. Furthermore, high redshift clusters are poor clusters as they are still fast growing according to the dark matter halo growth models (Press & Schechter 1974; Springel et al. 2005). Such poor clusters are difficult to locate due to their weak density contrasts to the general field. There are a number of ways to find high redshift clusters, but one of the proven techniques is to look around high redshift radio galaxies (Miley & De Breuck 2008). While not all the radio galaxies are in over-density regions, many of them host clear over-densities of galaxies around them (Venemans et al. 2007) and they are called proto-clusters. Although over-densities of red massive galaxies are not necessarily confirmed around them, they likely virialize and evolve to clusters at lower redshifts. Among the several proto-clusters reported so far, PKS1138-26 at is one of the most promising proto-clusters for its clear over-density of spectroscopically confirmed galaxies by previous studies (Miley & De Breuck 2008).
Early studies of the the PKS1138-26 radio galaxy at were performed by Pentericci et al. (1997, 1998), who reported a clumpy morphology of the radio galaxy. Followed by these initial observations, Kurk et al. (2000) first reported an over-density of star forming galaxies around the radio galaxy. This region was then followed up by several authors. Pentericci et al. (2000) performed spectroscopic follow-up observations of Lyman emitters reported in Kurk et al. (2000) and confirmed 14 galaxies close to the radio galaxy redshift. Kurk et al. (2004b, a) carried out further imaging and spectroscopic observations of the field targeting H emitters and confirmed another 10 galaxies at the cluster redshift. An X-ray observation has also been performed (Pentericci et al. 2002) and follow-up spectroscopy confirmed at least 5 X-ray sources at the cluster redshift (Croft et al. 2005). By now, there are more than 20 objects confirmed at the cluster redshift. Detailed analyses of the radio galaxy and the surrounding region with the superb resolution imaging with HST have also been performed (Miley et al. 2006; Zirm et al. 2008; Hatch et al. 2008, 2009), which added further lines of evidence for the forming (proto-)cluster. In fact, Zirm et al. (2008) reported on the forming cluster red sequence. Recently, Doherty et al. (2009) carried out a near-IR spectroscopic follow-up observation and confirmed two massive red galaxies at the cluster redshift.
Given the the convincing over-density of galaxies and wealth of imaging and spectroscopic data available in the field, PKS1138 is an ideal sample to study the environmental dependence of galaxy properties at this high redshift. In this paper, we perform an extensive analysis of galaxies around PKS1138 to quantify the environmental dependence of galaxy evolution and formation at .
The layout of the paper is as follows. We summarize our data in Section 2. We then describe details of our method of fitting spectral energy distributions of galaxies in Section 3. Before presenting our results, we perform extensive sanity checks in Section 4. Section 5 presents physical parameters of galaxies obtained from the fits as a function of environment at and Section 6 discusses implications of our results for galaxy formation. Finally, we summarize the paper in Section 7.
Unless otherwise stated, we adopt H, , and . Magnitudes are on the AB system. We use the following abbreviations: IMF for initial mass function, SED for spectral energy distribution and SFR for star formation rate.
We use two data sets for the purpose of this paper. One is from the (proto-)cluster field PKS1138 to examine galaxies in an over-density region at . The other one is from GOODS for a field counterpart for comparison. We summarize the two sets of data in this section.
The field has been imaged by several instruments — -band with LRIS on Keck (Zirm et al. 2008), band with FORS on VLT (Kurk et al. 2000, 2004b), bands with MOIRCS on Subaru (Kodama et al. 2007), with IRAC and with MIPS on Spitzer (Seymour et al. 2007). In addition to these published data, we used -band data taken with FORS2 on VLT and -band data with SOFI on NTT, which were obtained as a filler target of other programs and have not been published elsewhere. The and -band data were reduced in a standard manner. We further supplemented the data set with superb resolution imaging by the ACS on-board the HST (Miley et al. 2006). We retrieved the pipeline reduced and band images from the Hubble Legacy Archive. Table 1 summarizes the data.
The photometric zero points for the ground-based optical-nearIR images were obtained from the standard star observations. As pointed out by Doherty et al. (2009), the and band photometric zero-points used in Kodama et al. (2007) were off by 0.3 mag. We re-measured the zero points from the standard stars observed in the same nights and adopted those revised zero points. For the ACS images, we used the zero points from Sirianni et al. (2005). We have checked all the zero points against stars from Gunn & Stryker (1983). We convolved the SEDs of stars from Gunn & Stryker (1983) with filter responses and atmosphere and derive synthesized magnitudes of stars. We then compared the observed sequence of stars with the synthesized sequence on color-color diagrams with various color combinations. We found small zero point offsets with respect to the Gunn & Stryker (1983) stars and corrected for them to better match with Gunn & Stryker (1983). The applied offsets were typically mag and the largest was mag to the -band.
Objects were detected using sextractor (Bertin & Arnouts 1996) in the individual optical and near-IR bands because the seeing sizes vary from band to band and we did not smooth them to a common seeing. We used MAG_AUTO with aperture corrections assuming point sources for the optical and near-infrared photometry. For the IRAC images, we performed aperture photometry using the -band image for object detections. We then applied aperture corrections to obtain the total fluxes.
The aperture corrections were estimated as follows. We first constructed an average PSF image using bright, unsaturated stars in each band. Then we randomly distributed PSF objects in the images and repeated the object detection and photometry. The differences between the input magnitudes and measured magnitudes were used for the aperture corrections. We also estimated detection limits for point sources and they are summarized in Table 1. The last column in the table shows the 50% detection limits. For images taken with ground-based facilities, apparent sizes of faint sources are small and magnitude limits for point sources are a reasonable proxy for extended sources. For the ACS data, the magnitude limits for point sources are extremely deep due to the superb resolution. We measured typical magnitudes of the detected extended sources where their errors become magnitude (i.e., limits) and adopted them as magnitude limits for extended sources.
All the photometric errors were estimated from the sky noise in the same aperture sizes as used for objects. The photometric catalogs were cross-correlated with the -band catalog within the seeing FWHM of each band, and a -band selected catalog was produced. The final catalog contains and IRAC photometry and the Galactic extinction was corrected for in each band using the dust map from Schlegel et al. (1998). Stars are removed from the catalog based on their compactness and colors (we use a catalog in which stars are not removed only in Section 4.1).
Note that we have an IRAC image, but we do not use it as strong dust emission (e.g., PAH emission at ) from low redshift galaxies fall in this band. Our models described below do not include dust emission. We could in principle run photo- without the 8 photometry, select galaxies, and then run the full SED fits to the selected galaxies including the 8um photometry, which does not probe PAH at . But, that complicates the error analysis as it is not very straightforward to quantify how the errors propagate if we use two different sets of templates.
We also have a MIPS image of the PKS1138 field, but we only use it for a sanity check in Section 4. We apply a magnitude cut of () to the catalog to avoid any significant incompleteness effects. Also, we use only galaxies detected in more than 5 bands to ensure reliable SED fitting.
We use data from the GOODS-MUSIC sample (Grazian et al. 2006; Santini et al. 2009) for a field counterpart of the PKS1138 data. There is a wealth of data in the GOODS field, but we restrict ourselves to and IRAC data for this work so that we have almost the same wavelength sampling as the PKS1138 data. This is a crucial point of this work – if we used more data in the GOODS field, then we would not be able to compare PKS1138 and GOODS due to the different levels of systematic biases in the analyses. Note that we will use the MIPS photometry just for a sanity check in Section 4. The data covers a 78 arcmin field, which is limited by the -band field coverage. We apply the magnitude cut of to the GOODS catalog and use galaxies with detections in more than 5 bands as done for PKS1138.
To make sure that GOODS is not a peculiar field, we derive galaxy density in the Subaru XMM-Newton Deep Field (SXDF) using the photometric data described in Finoguenov et al. (2009), which covered approximately . The density at turns out to be similar to the GOODS density; in SXDF and in GOODS. This suggests that GOODS samples a typical Universe at .
|Band||Instrument||PSF size||50% limit|
3 The SED Fitting Method
The idea of the work is to use the data from the two fields with almost the same wavelength coverage and sampling, feed the data to the same SED fitting code, fit them with the same set of templates, and discuss relative differences between the two samples. We presented our data sets in the last section. We here describe details of our SED fitting code. The SED fitting procedure follows the conventional minimizing statistics. We prepare model templates and compare those templates with the observed data. We then analyze the fits and derive physical parameters out of the fits such as SFR and dust extinction. We start with preparing templates and then describe details of the fitting procedure.
3.1 Model Templates
The model templates are generated using an updated version of
the Bruzual & Charlot (2003) population synthesis code,
which takes into account the effects of thermally pulsating AGB stars.
We adopt the Salpeter initial mass function
star formation time scale () assuming the exponentially declining SFR (i.e., )
optical depth of dust extinction in the -band ()
age, which is the time since the onset of star formation to the observed epoch.
We allow to vary between 0 (single burst) and (constant SFR), and between 0 and 10. For the single burst model (), we assume no dust to mimic the passive evolution. For the other models, we do not assume any correlation between and . We adopt the two component dust extinction model of Charlot & Fall (2000). The actual amount of dust applied in the templates is not simply but is dependent on star formation histories. We measure a -band magnitude of a template and compare it with the -band magnitude of dust-free template of the same age, , and metallicity. The difference between the two is . Note that we make fine model grids where galaxies populate. We implement effects of the intergalactic extinction following Furusawa et al. (2000), who used the recipe by Madau (1995) and relaxed the allowed range of extinctions, as we explore the Universe. Also, we use SWIRE AGN templates from Polletta et al. (2007) as we might expect an increased fraction of AGNs in the PKS1138 region (Pentericci et al. 2002; Galametz et al. 2009).
We convolve the template spectra redshifted to various redshifts with response functions of the detectors and filters with/without atmosphere (we do not convolve atmosphere for data from the space) and generate synthesized fluxes. We use a logical constraint on redshift that the age of a template spectrum must be younger than the age of the Universe at that redshift. We have confirmed that we obtain essentially the same results if we do not use the constraint as we discuss below. The library of the synthesized fluxes for all the template is used to fit the observed SEDs of galaxies as detailed below. We use 2 million templates in total for both PKS1138 and GOODS.
3.2 Fitting procedure
We use the conventional minimizing statistics to fit the models to the observed data. The fit itself is performed in fluxes, not in magnitudes. Non-detections are treated as zero flux with flux errors from the magnitude limits (here we take the 50% detection limits). In this way, we avoid using upper limits in the fits and keep the statistics simple. We impose a prior constraint that all galaxies have stellar masses lower than . This is a reasonable constraint – the most massive galaxies at have (Kriek et al. 2008).
The best fitting template gives photometric redshift and physical parameters of galaxies. To be specific, we measure , , , age, SFR, and stellar mass from the fits. If the best fitting template is an AGN template, then we can not extract these parameters except for because there is an non-negligible contribution of AGN to the overall SED. We only use their photometric redshifts.
An uncertainty on each parameter is estimated by marginalizing over the parameter of interest and taking an interval of . However, all the parameters are correlated and this is just a rough error estimate. It is not easy to handle all the co-variances correctly, but fortunately, accurate random error estimates are not crucial for our purpose. It is because we discuss relative, systematic differences between the two data sets. We should have the same level of co-variances between the parameters as we use the same code, templates and have almost the same wavelength sampling of SEDs. We therefore take this simple error estimate, but we should bear in mind that our errors are not accurate. Fig. 1 shows an example of our SED fits.
4 Sanity Checks
Before we present our results, we perform a few sanity checks to make sure that our results are robust. Here we discuss (1) photometric zero points, (2) accuracy of our photometric redshifts, (3) accuracy of the SFR estimates from the SED fits, and (4) secondary burst models.
4.1 Photometric Zero Points
We cannot avoid random errors in the photometric zero points, but it is important to make sure that we do not have any systematic zero point offsets as a function of wavelength, which can cause color stretches in the color space. For example, if the zero points are systematically fainter in the blue bands, we might get higher extinction and/or older ages from the SED fits.
The easiest way to check the color stretch is to look at the sequence
of stars in color-color diagrams.
Fig. 2 is one such plot.
The plot shows vs. of the detected objects brighter than
in the two samples.
One can easily identify the stellar sequence in the lower part of the plot.
There is no strong systematic offset between the sequence of
stars in the PKS1138 field and that in the GOODS field.
We have also checked various color combinations and found that the
offsets are mags.
To accommodate with any possible zero point errors,
we add 0.1 mag in the quadrature to all the magnitude uncertainties
in both PKS1138 and GOODS catalogs before performing the SED fits
described in the last section
4.2 Photometric Redshifts
Another good check will be to look at the accuracy of photometric redshifts. We want to make sure that we can cull the majority of the galaxies at with a reasonably small amount of fore-/background contamination.
Fig. 3 compares our photometric redshifts () with spectroscopic redshifts (). The spectroscopic redshifts in PKS1138 are from Pentericci et al. (2000), Kurk et al. (2004b), Croft et al. (2005), and Doherty et al. (2009) and those in GOODS are from the GOODS-MUSIC catalog (Grazian et al. 2006 and references therein). The correlation is reasonably good, although there are a number of galaxies with wrong photo-. We find that galaxies with poor fits (the reduced larger than 3) often have wrong photometric redshifts. We remove those galaxies from the main analysis. Our results remain unchanged if we do not apply this cut as we quantify later.
A fraction of photo- outliers () is 5% in GOODS over the entire redshift range. Our spectroscopic sample for PKS1138 is strongly skewed at particular redshifts, and we do not quote the numbers for comparison here. The photo-z selection criterion for the main analysis is a trade off between the completeness and contamination. The selection adopted for the main analysis, within , is primarily motivated by the rate of recovering the spectroscopically confirmed PKS1138 members (it recovers 6 out of 7 confirmed PKS1138 members brighter than ). We have confirmed that our results presented below do not change if we change the selection criterion within a reasonable range (e.g., ). The contamination of fore-/background galaxies in our sample used in the main analysis is 36% in GOODS, which is estimated from the spectroscopically confirmed galaxies at with photometric redshift consistent with within . Note that the numbers quoted here should be considered rough estimates because the spectroscopic redshifts are collected from various surveys and the sample is quite heterogeneous. As we will show later, contamination due to wrong photometric redshifts (e.g., foreground galaxies at but with ) do not strongly affect our conclusions. We note in passing that 6 out of 50 galaxies in PKS1138 and 8 out of 50 galaxies in GOODS that are used in the main analysis are spectroscopic members.
It is also important to verify that the parameters derived from the SED fits are not strongly affected by the uncertainties in the photometric redshifts. We take into account the random errors when we discuss differences between PKS1138 and GOODS in the main analysis, but one may worry that photometric redshifts may introduce systematics in the SED fits and in the derived parameters.
To quantify this, we fit SEDs at redshifts fixed at for the spectroscopic objects. Fig. 4 compares the parameters derived at and those with redshifts as a free parameter. We observe a good correlation between them. In particular, there is no strong systematic offsets as shown in the top panel of Table 2, where we show the median differences between the physical parameters derived at and those at for the spectroscopic members at . The values in the brackets show the rate at which the two parameters agree within . Spectroscopic redshifts have narrow errors on all the parameters, but photometric redshifts do not introduce any systematics. We will focus on systematic differences between the two samples in the next section and all the differences we observe are larger than the offsets listed in Table 2. Therefore, the uncertainties in the photometric redshifts do not affect our results.
We will not use parameters derived at in the main analysis because we would like to show what parameters the contaminant galaxies (i.e., objects with wrong photometric redshifts) typically have and to show that they do not affect our conclusions. Also, the spectroscopic samples are strongly biased towards particular types of galaxies and we do not want to put a weight on them. For example, many of the spectroscopically observed galaxies in PKS1138 are either emission line objects or X-ray sources (Pentericci et al. 2000; Kurk et al. 2004b) and they are not necessarily typical galaxies. We use photometric redshift as a free parameter for all the objects so as not to introduce any biases. We have confirmed that our results do not change if we use the parameters derived at where available, as expected from Fig. 4.
Finally, we show that the logical age constraint - the age of a model template must be younger than the universe at a given redshift - does not change our results. The middle panel of Table 2 quantifies the effect of the age constraint on the derived parameters. We observe no systematic differences and the derived parameters of most galaxies agree within . We can therefore safely apply the constraint in the main analysis. The reason why the observed offset is zero is because the model grids are discrete and many galaxies remain in the same grid. SFRs do not have grids as they are normalized by the observed fluxes of galaxies, but most galaxies have the same best-fitting templates and the median SFRs remain the same.
|age [Gyr]||(75%)||0 (81%)|
|SF time scale [Gyr]||0 (100%)||0 (70%)|
|[mag]||0 (100%)||0 (88%)|
|age [Gyr]||0 (92%)||0 (96%)|
|SF time scale [Gyr]||0 (90%)||0 (100%)|
|SFR ||0 (76%)||0 (92%)|
|[mag]||0 (80%)||0 (97%)|
|age [Gyr]||0 (70%)||0 (79%)|
|SF time scale [Gyr]||0 (72%)||0 (80%)|
|SFR ||(56%)||0 (73%)|
4.3 SFR from SED Fits
Next, we check the accuracy of the SFRs from the SED fits. There are many ways to derive SFRs of galaxies, but here we take those from the warm dust emission probed by the MIPS observations. We take the conversion from the MIPS flux to the total IR flux by Reddy et al. (2006) and derive SFRs using the formula given by Kennicutt (1998). We have a deep MIPS image in both PKS1138 (Seymour et al. 2007) and GOODS (Santini et al. 2009), reaching down to and , respectively, at . It should be possible to measure SFRs from various spectroscopic observations carried out in the GOODS field, but the spectra are not always public.
We plot in Fig. 5 SFRs from the SED fits against SFRs from MIPS. Note that we plot only galaxies that are detected in MIPS and are at within used in the main analyses presented in the following sections. For the GOODS sample, the agreement between and is relatively good, although there are some outliers. Roughly of the galaxies have consistent SFRs from the SED fits and MIPS photometry within 1 dex at . Note that itself must have a significant error given the large scatter in mid-far IR SEDs of galaxies (Elbaz et al. 2002).
The PKS1138 sample has a similar rate of catastrophic failures to the GOODS sample – 6 out of 24 galaxies () plotted fall below MIPS SFRs by more than 1 dex, although only one of them is inconsistent at . SFRs for the PKS1138 galaxies have larger errors compared to the GOODS SFRs. This is because some of our optical images in the PKS1138 are shallow compared to the GOODS images. The optical images probe rest-frame UV at and they are crucial to pin down accurate SFRs.
We have visually inspected the galaxies with (i.e., the upper error is lower than ) used in the main analysis on the MIPS images. 75% (9 out of 12) of galaxies with low SFRs in GOODS are detected in MIPS and have . The other three are also detected but have . This high fraction of the MIPS detection illustrates the fundamental difficulty in distinguishing dusty starbursts from passive galaxies. But, this can be due to wrong photometric redshifts, AGN contamination, or object blending/confusion in the MIPS image due to the large PSF. We have looked at X-ray point sources in the field (Lehmer et al. 2005) and found that none of the low SFR galaxies are detected. We note in passing that 3 out of 50 GOODS galaxies used in the main analysis are detected in X-rays, but they do not affect our results.
In contrast to GOODS, the PKS1138 galaxies used in the main analysis show a difference between low and high SFR galaxies – only of galaxies (4 galaxies out of 21) with are detected in the MIPS image. Our MIPS image reaches at and if all the low SFR galaxies in PKS1138 are forming stars at similar rates to the GOODS galaxies, they should have been detected. We will show that PKS1138 galaxies have lower SFRs than GOODS galaxies in the next section and the observed lower frequency of the MIPS detections in PKS1138 gives a strong support to this result. We note that Galametz et al. (2009) suggested an increased fraction of AGNs in proto-clusters and some of the PKS1138 galaxies with low SFRs may host an AGN, boosting the mid-IR flux. In fact, one of the low SFR galaxies that is used in the main analysis is detected in X-rays (Pentericci et al. 2002). This is the only object detected in X-rays and used in the main analysis and this object does not affect our result in any way.
We admit that the correlation between and those measured from other methods is not particularly good, and we do not attempt to discuss precise SFRs of individual galaxies. Instead, we classify galaxies into two coarse classes by their SFRs. We define high SFRs as , and low SFRs are lower than that. We do not try to further subdivide the SFRs. However, we will later introduce another class - starbursts - as just to illustrate starbursting populations at the time of cluster formation at very high redshifts.
4.4 Secondary burst models
As we discuss later, cluster galaxies may have experienced many mergers in the past and the models may not reproduce their star formation histories very well. We have implemented a secondary burst on top of the models and checked how the results change. We assume an instantaneous burst and add a secondary burst either at 0.2, 0.5, 1, 2, and 5 Gyr after the onset of star formation and the burst adds either 10% or 50% of stellar mass of the galaxy.
The bottom line of this exercise is that our conclusions remain unchanged. The bottom panel of Table 2 quantifies the difference between the secondary burst models and the models. We observe no strong differences in age and . The agreement of SFRs between the two models is not particularly good, but our results remain unchanged because the systematic differences are very small. The dust extinction shows a poor agreement and a systematic offset in . But, the amount of the offset is similar in PKS1138 and GOODS, and our conclusion will therefore remain the same as we discuss only relative differences between the two samples.
The secondary burst models are probably more realistic than the simple models. However, the addition of the secondary burst raises the number of free parameters in our SED fitting from 5 to 7 (or 8 if we allow the star formation time scale of the secondary burst to vary), while we only have 11 broad band photometry points. As this introduces too many degeneracies between the parameters and the secondary burst complicates the interpretation of the results, we will consider only the models for our main analysis in the following sections.
We now move on to present the results from the SED fits. We stress again that we use the two samples with almost the same wavelength samplings, feed them to the same code with the same templates, and discuss the relative differences between the two samples.
In what follows, we study galaxies with , detected in more than 5 bands, and being consistent with within (note the PKS1138 radio galaxy is at ). We exclude galaxies with poor fits () as they likely have wrong photometric redshifts as discussed in the last section. We take the same sample definition in both PKS1138 and GOODS, so that we can make a fair comparison. The physical parameters we focus on in this section are age, star formation time scale (), dust extinction, and SFR. We first present evidence that the PKS1138 is indeed an over-density region as suggested by the earlier studies. We then move on to discuss physical properties of galaxies from the SED fits as a function of environment.
5.1 Over-density of galaxies in PKS1138
Fig. 6 compares the distributions of galaxies around the radio galaxy redshift in the PKS1138 and GOODS fields. We could not obtain a good enough photometric redshift for the PKS1138 radio galaxy, (i.e., it is 2.3 away from the spectroscopic redshift) and it goes out of our main sample. But, the radio galaxy is a powerful AGN and our photometric redshift is relatively good for such an extreme object. We show the radio galaxy in Fig. 6 just to illustrate the location of the radio galaxy.
The galaxies around PKS1138 are arrayed in a filamentary structure in the east-west direction, and the radio galaxy is located at the center. There is an over-density of galaxies around the central radio galaxy. Interestingly, many of the galaxies around the radio galaxy show relatively low SFRs. On the other hand, the galaxy distribution in the GOODS field is quite sparse and there is no concentration of low SFR galaxies. The surface densities of galaxies () are and in PKS1138 and GOODS, respectively, suggesting an over-density of a factor of in PKS1138, which is in rough agreement with Kurk et al. (2004b) who found a factor of over-density of H galaxies. The over-density is more pronounced for bright galaxies ; and for (i.e., a factor of 9). All the above numbers for PKS1138 are averaged over the entire probed field and the over-density can be further pronounced if we take only the central part of the PKS1138 field – the density of PKS1138 galaxies within 1′ from the radio galaxy is for (i.e., a factor of 16) and for (i.e., a factor of 20).
One may wonder that nearby fore-/background structure may be contaminating the structure we observe in PKS1138. But, the contamination should be small because the observed galaxy distribution is similar to the distribution of H emitters at the radio galaxy redshift (Kurk et al. 2004b) – weak concentrations of galaxies on the West and South-East of the radio galaxy are both seen in the distributions of photo- selected and H selected galaxies. Also, structures at different redshifts are unlikely to make a coherent structure around the radio galaxies with a concentration of low SFR galaxies. Zirm et al. (2008) found a concentration of red galaxies selected in colors around the very vicinity of the radio galaxy (″), we do not observe such a strong concentration of red galaxies, but it might be a concentration of faint galaxies (we use , while they used , roughly corresponding to ).
Fig. 7 plots color magnitude diagrams of galaxies in the two fields. We model the location of the red sequence using the updated Bruzual & Charlot (2003) population synthesis models following the procedure described in Lidman et al. (2008). The most striking trend is that there is a hint of the red sequence, which is a ubiquitous feature of galaxy clusters at lower redshifts, at in PKS1138. The bright red galaxies are around the model red sequence formed at . This suggests that these red galaxies formed the bulk of their stars around that redshift. The possible red sequence in PKS1138 is populated half by galaxies with low SFRs () and half by galaxies with high SFRs, suggesting that the red sequence is being formed. We note that spectroscopically confirmed contaminant galaxies tend to be blue, and thus the possible red sequence in PKS1138 is not due to contamination from fore-/background galaxies. In contrast to PKS1138, there is no clear sign of a red sequence in the GOODS field, and there are very few galaxies brighter than .
To further illustrate the over-density of PKS1138, we plot stellar mass functions in Fig. 8. We do not derive volume densities of galaxies as we are probably looking at a coherent structure in redshift in PKS1138, whose size we do not yet know precisely. We simply derive a surface density. The over-density of PKS1138 is evident – PKS1138 hosts galaxies more densely than GOODS. We will focus on galaxies more massive than in the next section, for which our samples are nearly complete in both PKS1138 and GOODS.
To sum up, we observe the significant over-density of galaxies and the forming red sequence in the PKS1138. This suggests that there is a proto-cluster (or perhaps a real cluster) around the radio galaxy. But, we have to wait for intensive spectroscopic follow-up observations to confirm it. In contrast to PKS1138, we do not observe any possible (proto-)clusters nor red sequence in the GOODS field and the galaxy density is lower. Therefore, the two samples probe different environments at and they provide us with a unique opportunity to quantify the dependence of galaxy properties on density at this high redshift.
5.2 Physical parameters from the SED fits
We are now going to look into differences in galaxy properties between the two environments. We quantify four fundamental properties of galaxies – age (which is time since the onset of star formation), star formation time scale, SFR, and dust extinction as summarized in Fig. 9. We are aware of the extremely uncertain nature of various parameters gone into the SED fits. For example, initial mass function (IMF) is one of the most uncertain assumptions. A number of IMFs are suggested in the literature, each of which produce systematically different results, and to make matters worse, IMFs may depend on mode of star formation, e.g., IMFs could become top-heavy during starbursts. Also, the models we adopt may not be a good approximation of star formation histories of real galaxies, we assume all the stars have the same metallicities, and etc. Again, all the uncertainties are inherent in the both samples and we discuss only the relative differences between the two samples, which must be the most robust results from the SED fits.
We discuss the four parameters derived from the SED fits.
Since all the parameters are correlated with stellar masses of
galaxies and we plot them as a function of stellar mass.
Note that we miss a fraction of galaxies with masses lower than
in PKS1138, and
we focus on galaxies more massive than .
This is a conservative mass cut and our results below
do not suffer from any strong incompleteness effects.
But, we keep lower mass galaxies in the plots below to illustrate
the mass dependency of the four parameters.
AGE (Fig. 10a) :
The age distributions are not significantly different.
We find that the fractions of galaxies older than 2 Gyr are
and , respectively
We note that Steidel et al. (2005) studied a spectroscopically
confirmed redshift spike at and found that
galaxies in the spike are older than those outside the spike.
Their result may appear inconsistent with ours.
But, their sample is a collection of star forming galaxies,
while we probe the entire galaxy population.
A fair comparison cannot be made.
STAR FORMATION TIME SCALE, (Fig. 10b) :
PKS1138 galaxies tend to have shorter formation time scales.
The fractions of galaxies with Gyr are
and in PKS1138 and GOODS, respectively.
The Mann-Whitney test suggests that the distributions
are likely different (median probability of %).
We obtained probabilities in 52% of the Monte-Carlo realizations.
As seen in the GOODS data, the fore-/background contaminants do not
favor any particular , and the observed difference is not due
to the contamination in PKS1138.
Combining with the age distribution, most galaxies in both fields have
, suggesting that they already formed the bulk of
their stars by the time of the observation.
The shorter formation time scale of PKS1138 is one of the most
important results in this paper, and we will further discuss it
in the next section.
SFR (Fig. 10c) :
PKS1138 galaxies tend to show weaker star formation and lower specific SFRs.
The fractions of galaxies that have are
and in PKS1138 and GOODS, respectively.
The Mann-Whitney test suggests that they are not from the same parent
The fore-/background contamination tends to have high star formation.
They unlikely to contribute to the observed difference.
We recall that we obtained a lower frequency of MIPS detections
in PKS1138 than in GOODS, which gives a further support
to the lower SFRs in PKS1138.
DUST EXTINCTION (Fig. 10d) :
PKS1138 galaxies tend to have less dust, as expected from their low SFRs.
The fraction of PKS1138 galaxies with is
( is a typical amount of dust in local star forming galaxies).
Most of galaxies in GOODS are very dusty, and there seems a correlation
between stellar mass and dust extinction (e.g., Reddy et al. 2006).
The fraction of galaxies is .
The Mann-Whitney test shows that the distributions
are different (%).
To sum up, we find that galaxies in PKS1138 and GOODS have different properties. Interestingly, the trend we observe in clusters still qualitatively persists even at this high redshift — we observe galaxies with lower SFR and less dust in higher density regions, where we see clearer red sequence. The environmental dependence of galaxy properties is already in place at , at least partly. We will further extend the discussion and quantify the environmental dependence of galaxy formation in the next section.
From the extensive SED fits of galaxies in PKS1138 and GOODS, we have found that massive galaxies in PKS1138 tend to have
shorter star formation time scales
lower star formation rates
lower amounts of dust
compared to those in GOODS. The combination of first and second points is interesting. Galaxies in PKS1138 and GOODS start forming stars at a similar epoch in a statistical sense, but PKS1138 galaxies form more rapidly. We recall that our definition of age is the time since the onset of star formation. The third and fourth points are basically by-products of the first and second points. At , PKS1138 galaxies have already undergone intense star formation and their SFRs are rapidly declining, while GOODS galaxies are still actively forming stars due to their longer star formation time scales. This results in the lower SFRs in PKS1138 (third point). The fourth point can then be easily understood given the correlation between between SFR and dust amount (lower SFR galaxies have less dust, e.g., Hopkins et al. 2003).
Cluster galaxies have a shorter star formation time scale – this is the same trend as observed in clusters. For example, Gobat et al. (2008) showed the same trend in a cluster based on a photo-spectroscopic analysis. It is striking that the trend holds even at . The difference in the star formation time scale suggests that cluster and field galaxies may form in different ways. Let us discuss the formation and evolution of cluster and field galaxies in detail. We first introduce two parameters to further quantify the galaxy formation. We then extend the discussion and address the origin of the environmental dependence of galaxy properties observed at lower redshifts. We note a caveat here that our results are based on one (proto-)cluster only and the trends we observe may not represent global trends at . A larger sample of clusters should be investigated to draw a global picture.
6.1 The environmental dependence of galaxy formation
We have age, star formation time scale, and SFR for each galaxy. Assuming the exponentially decaying star formation histories, we can estimate the SFR at the onset of star formation, which we call initial SFR (; see Fig. 9). Fig. 11 plots as a function of stellar mass. The cluster galaxies tend to have higher – the fractions of are , and in PKS1138 and GOODS, respectively. The Mann-Whitney test supports this difference (the median probability is % and we obtain % probability in 74% of the realizations). PKS1138 galaxies have experienced more intense formation histories. Of course, all this is based on the assumption of the exponentially declining SFRs, and we have ignored effects of galaxy-galaxy mergers, which we will discuss later.
If galaxies continue to form stars following the exponential decay, we can derive the time when the galaxies form bulk of their stars. We plot in Fig. 12 the epoch when galaxies form 80% of the stars that they would have at (; see Fig. 9). PKS1138 galaxies typically formed around or higher, while GOODS galaxies typically formed below . We recall that we obtained from the location of the red sequence in PKS1138 in Section 5.1, which is in agreement with what we find here. We find that of galaxies in PKS1138 form the bulk of their stars at , while the fraction is in GOODS. The Mann-Whitney test suggests that they are likely different – the median probability is and we obtained % probabilities in 38% of the realization. Interestingly, the ages of the cluster and field galaxies are not very different (Fig. 10a). The difference in is therefore due to their shorter formation time scales. The formation of cluster galaxies is a more intense event, they form in a shorter time scale, and they assemble the bulk of their stars earlier by Gyr than the field galaxies.
To summarize, we measure the median and SFR for galaxies and illustrate the differences in star formation histories of galaxies in PKS1138 and GOODS in Fig. 13. This will the summary plot of the paper. As discussed above, PKS1138 galaxies have higher and shorter . PKS1138 experience much more intense galaxy formation at early times than GOODS galaxies. The galaxy assemblies are completed in a short time, and by the time we observe them ( Gyr after the onset of star formation as shown in Fig. 10a) their SFRs become low. The plot summarizes the points listed at the beginning of this section (except for the fourth one, which is a result of SFR-dust correlation).
The plot indicates that PKS1138 galaxies and GOODS galaxies did not form in the same way. The formation of cluster galaxies is a more intense event. But, how can massive galaxies form on such a short time scale? In general, star formation activities have negative feedback to themselves. Hot, young stars ionize surrounding gas, which prevents further star formation. Also, supernova explosions of massive stars give negative feedback. If galaxies undergo an AGN phase, it might also give negative feedback. The formation of massive galaxies on a very short time scale may be difficult to understand in this respect, and that makes us speculate that these massive galaxies might have taken another route – mergers.
Galaxy formation takes place in high peaks of density fluctuations in the Universe. Clusters that we observe today were regions of large-scale over-densities, where many density peaks were embedded. On the other hand, field galaxies form in more isolated peaks. The current galaxy formation theories predict that small galaxies form first and they form progressively massive galaxies via mergers. In the early Universe, cluster galaxies must have had a lot better chances to merge with other galaxies than field galaxies simply because there are more galaxies around them. This accelerated frequency of mergers might be the cause of the rapid formation of cluster galaxies.
A simple scenario of an equal mass merger can be helpful here to show that mergers tend to make higher and push to higher redshifts. Assuming that two equal mass galaxies merge into a single galaxy, the merger event will double the stellar mass of the galaxy. If the galaxy experiences starburts triggered by the merger, its stellar mass will be more than double and its SFR gets lower after the burst. We then observe this massive galaxy with low SFR, and depending on the time elapsed after the burst, we would fit a short to this galaxy in order to reproduce its low SFR and large stellar mass. In fact, we find that the secondary burst discussed in Section 4.4 reduces the median of the PKS1138 galaxis by , although it is still high, , which might indicate that galaxies experienced more than one merger. On the other hand, the of the GOODS galaxies remain the same within .
We are not yet sure if PKS1138 is a collapsed cluster or currently collapsing proto-cluster. But, even if it is a cluster, it must be a young system, and it is unlikely that galaxies had enough time to be affected by ’nurture’ effects such as ram-pressure stripping. We may suggest that mergers in the early times of the cluster formation may be one of the key processes to establish the environmental dependence of galaxy properties observed at lower redshifts. We will further pursue this point in the following section.
6.2 The Nature effects in the environmental dependence of galaxy properties
Galaxy properties such as colors and morphology are known to depend on environment in which galaxies reside. This environmental dependence is shaped by two effects – nature and nurture effects. That is, both how galaxies form and how they evolve are important.
Fig. 12 shows that only Gyr has passed since the PKS1138 galaxies formed bulk of their stars to the observed epoch. Assuming that PKS1138 is a virialized system, we apply a -clipped gapper method (Beers et al. 1990) and obtain a velocity dispersion of 400 using spectroscopic redshifts within 2 arcmin, corresponding to a physical scale of 1 Mpc, from the literature (Pentericci et al. 2000; Kurk et al. 2004b; Croft et al. 2005). The virial radius () of the cluster is 0.32 Mpc, giving a crossing time scale of Gyr. Therefore, the time elapsed since the bulk formation of the PKS1138 galaxies is comparable to the crossing time scale of this cluster. This suggests that nurture effects may not have had enough time to fully work. They may have affected a fraction of galaxies, but they probably could not change the average properties of galaxies. We may be witnessing the nature effects in shaping the environmental dependence in PKS1138.
We admit that it is not very straightforward to classify nature and nurture effects at this high redshift. Now, let us introduce two kinds of mergers; early-epoch mergers and late-epoch mergers. We refer to mergers occur during the first collapse of clusters as early-epoch mergers, and those occur afterwards as late-epoch mergers. This is to sort out the two effects, early-epoch mergers being nature effects (initial conditions) and late-epoch mergers being nurture effects (environmental effects). Of course, nurture effects also include ram-pressure stripping, harassment, etc. On the other hand, nature effects do not include these processes driven by the deep potential well or intracluster medium.
In this classification, the environmental dependence of PKS1138 is likely due to early-epoch mergers because Gyr time is probably not enough for nurture effects to fully work. Galaxy clusters at do not give useful information about nature effects because nurture effects have had enough time to fully operate ( Gyr since the formation epoch) and it is not straightforward to disentangle the two effects. Our results suggest that nature effects are strong and they form the basis of the environmental dependence of galaxy properties. A way to probe the significance of the nature effects is to quantify morphology of the PKS1138 galaxies. If early-epoch mergers are an important effect, then we expect to observe more signs of recent interactions in PKS1138 than in GOODS (see below). Also, we expect to observe post-starburst galaxies if galaxies undergone interaction triggered starbursts. A near-IR spectroscopic follow-up campaign of PKS1138 is currently underway. We will be able to study spectral properties of the PKS1138 galaxies, which allow us to look deeper into star formation histories.
To sum up, the strong nature effects may shape the environmental dependence of galaxy properties. An ultimate goal of environment studies will be to quantify the relative contribution of nature and nurture effects. But, that will require statistical work on (proto-)clusters. Although the possible significance of nature effects we suggest here is based only on one (proto-)cluster PKS1138, let us further discuss it. It has an interesting implication for the build-up of the cluster red sequence.
6.3 The massive end of the cluster red sequence
The sequence of red early-type galaxies is a ubiquitous feature of galaxy clusters. Over the last few years, there is an accumulating amount of evidence in the literature that the cluster red sequence grows from the massive end to the low-mass end (e.g., Tanaka et al. 2005, 2007; Koyama et al. 2007; Tanaka et al. 2008; Gilbank et al. 2008). The massive end of the red sequence cannot be formed via a simple fading of blue, star forming galaxies because such massive blue galaxies do not exist even at (see Fig. 7 and discussions in Faber et al. 2007). We need mergers to form it.
An interesting point here is that we observe the brightest tip of the red sequence in PKS1138, which was also noted by Zirm et al. (2008). The massive end of the red sequence in a young system – this might be due to early-epoch mergers. Early-epoch mergers might have formed very massive galaxies during the first gravitational collapse of clusters. Early-epoch mergers should occur more frequently in cluster environments than in the field, and that helps explain why we do not observe red sequence in GOODS. As suggested by Zirm et al. (2008), the red sequence in PKS1138 may be being formed or just formed at the time of observation. This view is further supported by our observation that roughly half of the red sequence galaxies have high SFRs (Fig. 10). This formation redshift of is in line with predictions from the build-up of the red sequence observed in lower redshift clusters (Tanaka et al. 2007).
We present in Fig. 14 ACS -band images of the bright red galaxies in PKS1138 ( and ). Here we only briefly discuss their morphologies and a detailed study will be presented elsewhere (Zirm et al. in prep). Interestingly, a half of the galaxies show disturbed morphology and/or have nearby companions, lending a support to the picture of the accelerated mergers in clusters. Two out of three apparently disturbed galaxies are detected in MIPS. Early-epoch mergers during the first collapse of clusters may form the brightest end of the red sequence, and at the same time, form the basis of the environmental dependence of galaxy properties.
The low-mass end of the red sequence is built up at later times. Nurture effects may come in there. Once clusters form, intracluster gas and deep potential field, in addition to late-epoch mergers, can affect galaxies and terminate their star formation activities. In this way, moderate-low mass galaxies could become red and form the low-mass end of the red sequence. The combination of nature and nurture effects may have conspired to produce the down-sizing behavior of the observed build-up of the red sequence (e.g., Tanaka et al. 2005).
6.4 Proto-clusters at higher redshifts
Finally, we finish the discussion with a forecast for future (proto-)cluster studies. We can go back in time and see how the PKS1138 field looked like at higher redshifts as we have star formation histories of individual galaxies. A caveat of course is that we ignore all the early-epoch mergers occurred before the time of the observation (the only way to recover that information is to resolve the galaxies into individual stars). Another caveat is that we cannot track spatial positions of galaxies back in time. But, PKS1138 is likely a collapsing/collapsed system and it may well have been an over-density region already at .
Fig. 15 shows distribution of galaxies and their star formation rates evolved back to . Compared to Fig. 6, galaxies are more actively forming stars on average. Interestingly, galaxies in PKS1138 are more actively forming stars and there are more starbursting galaxies with than in GOODS. This is in stark contrast to Fig. 6, where we saw that PKS1138 galaxies have lower SFRs than those in GOODS. We still observe a hint of a galaxy over density in PKS1138 at . The plots suggest that, as we approach the formation epoch of clusters (i.e., early phase of the gravitational collapse to a massive cluster halo), we expect to observe an over density of starbursting galaxies. Along with starbursting galaxies, low SFR galaxies already appear in PKS1138, while such galaxies are extremely rare in GOODS. At this redshift, low SFR galaxies and starbursting galaxies may co-exist in a forming cluster.
It is not easy to show how early-epoch mergers change the picture we see here. But, we still expect to observe starbursting galaxies triggered by early-epoch mergers in collapsing clusters with a higher over-densities of lower-mass, pre-merger galaxies. We deem that observations with existing/future sub-millimeter arrays would be able to discover many forming clusters at high redshifts. In fact, some submm observations of distant radio galaxies have found possible over-densities of dusty starburst populations around them (e.g., De Breuck et al. 2004; Greve et al. 2007). Full wavelength observations will be essential for future proto-cluster studies at very high redshifts.
We have studied the environmental dependence of galaxy properties at based on the multi-band data available in the (proto-)cluster field PKS1138 and in GOODS. We have performed the extensive SED fits with a special care for systematic biases between the two samples. The results from the careful SED fits suggest that the environmental dependence is at least partly in place at this high redshift.
We have first shown that PKS1138 is indeed an over-density region with an excess of red galaxies, forming the brightest tip of the red sequence, compared to GOODS. The red galaxies tend to cluster around the radio galaxy. These results support the claim that PKS1138 is a (proto-)cluster at . Interestingly, the red sequence is populated both by low SFR galaxies and high SFR galaxies, suggesting that the red sequence is being formed.
We then have looked into detailed properties of galaxies derived from the SED fits. PKS1138 galaxies have similar age (we define age as time since the onset of star formation), shorter star formation time scale, lower SFR, and less dust compared to those in GOODS at similar redshifts. The averaged star formation history of the PKS1138 suggests that the cluster galaxies form on a shorter time scale and they form the bulk of their stars Gyr earlier than the field galaxies, which is consistent with lower redshift observations (e.g., Thomas et al. 2005; Gobat et al. 2008).
The environmental dependence of galaxy properties should be shaped both by nature and nurture effects. The observed environmental dependence at suggests that nature effects may be a strong effect as the PKS1138 (proto-)cluster is likely a young system. Possibly an accelerated rate of mergers in collapsing clusters may have played a role there and they may also be the primary cause of the short formation time scale of cluster galaxies. However, further studies on (proto-)clusters will be needed to confirm this picture. We expect that star bursting galaxies populate in very high redshift proto-clusters and sub-millimeter observations may be useful to confirm very high- systems.
Acknowledgements.We thank Nick Seymour, Tadayuki Kodama and Andrew Zirm for providing us with their data on PKS1138, Chris Lidman for reducing the SOFI data obtained at the European Southern Observatory using the Very Large Telescope on Cerro Paranal through ESO program 66.A-0597 and 167.A-0409, and also using the New Technology Telescope on Cerro La Silla through ESO program 076.A-0670. We thank Ricardo Demarco for helping with the NTT observation, and Stephane Charlot and Gustavo Bruzual for providing us with their latest stellar population synthesis code. We are grateful to the anonymous referee for useful comments, which helped improve the paper. This work was supported by World Premier International Research Center Initiative (WPI Initiative), MEXT, Japan. JK thanks the DFG for support via German-Israeli Project Cooperation grant STE1869/1-1.GE625/15-1. This study is also based on observations made with the NASA/ESA Hubble Space Telescope, and obtained from the Hubble Legacy Archive, which is a collaboration between the Space Telescope Science Institute (STScI/NASA), the Space Telescope European Coordinating Facility (ST-ECF/ESA) and the Canadian Astronomy Data Centre (CADC/NRC/CSA).
- offprints: M. Tanaka
- The Chabrier initial mass function (Chabrier 2003) gives smaller SFRs and stellar masses by a factor of compared to the Salpeter IMF. Other than that, it gives the same results as those presented below.
- The addition of random errors is not an ideal way to remedy the systematic errors. But, we apply the additional errors to make sure that the systematic errors do not dominate the overall error budget. The price we have to pay is that the additional errors slightly degrade the accuracy of photometric redshifts.
- The top panel shows the median differences between the parameters derived at and those at . The numbers in the parenthesis shows the fraction at which the two parameters agree within . The middle panel shows the differences we obtain if we do not use the logical age constraint. The numbers in the bottom panel mean the median differences between parameters from the secondary burst models and those from the fiducial models.
- As discussed in Section 4.2, we do not use galaxies with poor fits (). The fractions do not change if we use them; and for PKS1138 and GOODS, respectively. The fractions quoted below for the other parameters also remain essentially the same: and for 0.5 Gyr, and for , and and for mag.
- Beers, T. C., Flynn, K., & Gebhardt, K. 1990, AJ, 100, 32
- Bertin, E. & Arnouts, S. 1996, A&AS, 117, 393
- Blakeslee, J. P., Franx, M., Postman, M., et al. 2003, ApJ, 596, L143
- Bruzual, G. & Charlot, S. 2003, MNRAS, 344, 1000
- Chabrier, G. 2003, PASP, 115, 763
- Charlot, S. & Fall, S. M. 2000, ApJ, 539, 718
- Cooper, M. C., Newman, J. A., Weiner, B. J., et al. 2008, MNRAS, 383, 1058
- Croft, S., Kurk, J., van Breugel, W., et al. 2005, AJ, 130, 867
- Cucciati, O., Iovino, A., Marinoni, C., et al. 2006, A&A, 458, 39
- De Breuck, C., Bertoldi, F., Carilli, C., et al. 2004, A&A, 424, 1
- Doherty et al. 2009, A&A, submitted
- Elbaz, D., Cesarsky, C. J., Chanial, P., et al. 2002, A&A, 384, 848
- Elbaz, D., Daddi, E., Le Borgne, D., et al. 2007, A&A, 468, 33
- Faber, S. M., Willmer, C. N. A., Wolf, C., et al. 2007, ApJ, 665, 265
- Furusawa, H., Shimasaku, K., Doi, M., & Okamura, S. 2000, ApJ, 534, 624
- Galametz, A., Stern, D., Eisenhardt, P. R. M., et al. 2009, ApJ, 694, 1309
- Gebhardt, K., Faber, S. M., Koo, D. C., et al. 2003, ApJ, 597, 239
- Gilbank, D. G., Yee, H. K. C., Ellingson, E., et al. 2008, ApJ, 673, 742
- Gobat, R., Rosati, P., Strazzullo, V., et al. 2008, A&A, 488, 853
- Grazian, A., Fontana, A., de Santis, C., et al. 2006, A&A, 449, 951
- Greve, T. R., Stern, D., Ivison, R. J., et al. 2007, MNRAS, 382, 48
- Gunn, J. E. & Stryker, L. L. 1983, ApJS, 52, 121
- Hatch, N. A., Overzier, R. A., Kurk, J. D., et al. 2009, MNRAS, 395, 114
- Hatch, N. A., Overzier, R. A., Röttgering, H. J. A., Kurk, J. D., & Miley, G. K. 2008, MNRAS, 383, 931
- Hopkins, A. M., Miller, C. J., Nichol, R. C., et al. 2003, ApJ, 599, 971
- Kennicutt, Jr., R. C. 1998, ARA&A, 36, 189
- Kodama, T., Tanaka, I., Kajisawa, M., et al. 2007, MNRAS, 377, 1717
- Koyama, Y., Kodama, T., Tanaka, M., Shimasaku, K., & Okamura, S. 2007, MNRAS, 382, 1719
- Kriek, M., van Dokkum, P. G., Franx, M., et al. 2008, ApJ, 677, 219
- Kuntschner, H., Smith, R. J., Colless, M., et al. 2002, MNRAS, 337, 172
- Kurk, J. D., Pentericci, L., Overzier, R. A., Röttgering, H. J. A., & Miley, G. K. 2004a, A&A, 428, 817
- Kurk, J. D., Pentericci, L., Röttgering, H. J. A., & Miley, G. K. 2004b, A&A, 428, 793
- Kurk, J. D., Röttgering, H. J. A., Pentericci, L., et al. 2000, A&A, 358, L1
- Lehmer, B. D., Brandt, W. N., Alexander, D. M., et al. 2005, ApJS, 161, 21
- Lidman, C., Rosati, P., Tanaka, M., et al. 2008, A&A, 489, 981
- Madau, P. 1995, ApJ, 441, 18
- Mei, S., Holden, B. P., Blakeslee, J. P., et al. 2009, ApJ, 690, 42
- Miley, G. & De Breuck, C. 2008, A&A Rev., 15, 67
- Miley, G. K., Overzier, R. A., Zirm, A. W., et al. 2006, ApJ, 650, L29
- Nakata, F., Kodama, T., Shimasaku, K., et al. 2005, MNRAS, 357, 1357
- Pentericci, L., Fan, X., Rix, H., et al. 2002, AJ, 123, 2151
- Pentericci, L., Kurk, J. D., Röttgering, H. J. A., et al. 2000, A&A, 361, L25
- Pentericci, L., Roettgering, H. J. A., Miley, G. K., Carilli, C. L., & McCarthy, P. 1997, A&A, 326, 580
- Pentericci, L., Roettgering, H. J. A., Miley, G. K., et al. 1998, ApJ, 504, 139
- Polletta, M., Tajer, M., Maraschi, L., et al. 2007, ApJ, 663, 81
- Postman, M., Franx, M., Cross, N. J. G., et al. 2005, ApJ, 623, 721
- Press, W. H. & Schechter, P. 1974, ApJ, 187, 425
- Reddy, N. A., Steidel, C. C., Fadda, D., et al. 2006, ApJ, 644, 792
- Salpeter, E. E. 1955, ApJ, 121, 161
- Santini, P., Fontana, A., Grazian, A., et al. 2009, ArXiv e-prints
- Schlegel, D. J., Finkbeiner, D. P., & Davis, M. 1998, ApJ, 500, 525
- Seymour, N., Stern, D., De Breuck, C., et al. 2007, ApJS, 171, 353
- Sirianni, M., Jee, M. J., Benítez, N., et al. 2005, PASP, 117, 1049
- Springel, V., White, S. D. M., Jenkins, A., et al. 2005, Nature, 435, 629
- Stanford, S. A., Romer, A. K., Sabirli, K., et al. 2006, ApJ, 646, L13
- Steidel, C. C., Adelberger, K. L., Shapley, A. E., et al. 2005, ApJ, 626, 44
- Tanaka, M., Finoguenov, A., Kodama, T., et al. 2008, A&A, 489, 571
- Tanaka, M., Kodama, T., Arimoto, N., et al. 2005, MNRAS, 362, 268
- Tanaka, M., Kodama, T., Kajisawa, M., et al. 2007, MNRAS, 377, 1206
- Thomas, D., Maraston, C., Bender, R., & Mendes de Oliveira, C. 2005, ApJ, 621, 673
- Venemans, B. P., Röttgering, H. J. A., Miley, G. K., et al. 2007, A&A, 461, 823
- Zirm, A. W., Stanford, S. A., Postman, M., et al. 2008, ApJ, 680, 224