quiescent galaxies in the high-redshift clusters

# Evolution of Star-formation Properties of High-redshift Cluster Galaxies since z=2

Seong-Kook Lee11affiliation: Center for the Exploration of the Origin of the Universe, Department of Physics and Astronomy, Seoul National University, Seoul, Korea 55affiliation: sklee@astro.snu.ac.kr, mim@astro.snu.ac.kr , Myungshin Im11affiliation: Center for the Exploration of the Origin of the Universe, Department of Physics and Astronomy, Seoul National University, Seoul, Korea 55affiliation: sklee@astro.snu.ac.kr, mim@astro.snu.ac.kr , Jae-Woo Kim11affiliation: Center for the Exploration of the Origin of the Universe, Department of Physics and Astronomy, Seoul National University, Seoul, Korea , Jennifer Lotz22affiliation: Space Telescope Science Institute, 3700 San Martin Drive, Baltimore, MD 21218, USA , Conor McPartland33affiliation: Institute for Astronomy, University of Hawaii, 2680 Woodlawn Drive, Honolulu, HI 96822, USA , Michael Peth44affiliation: Department of Physics and Astronomy, Johns Hopkins University, 3400 North Charles Street, Baltimore, MD 21218-2686, USA , Anton Koekemoer22affiliation: Space Telescope Science Institute, 3700 San Martin Drive, Baltimore, MD 21218, USA
###### Abstract

Using a stellar mass limited sample of galaxies () at , we show that the stellar mass, rather than the environment, is the main parameter controlling quenching of star formation in galaxies with out to . On the other hand, the environmental quenching becomes efficient at regardless of galaxy mass, and it serves as a main star formation quenching mechanism for lower mass galaxies. Our result is based on deep optical and near-infrared imaging data over 2800 arcmin, enabling us to negate cosmic variance and identify 46 galaxy cluster candidates with . From to , the fraction of quiescent galaxies increases by a factor of over the entire redshift range, but the difference between cluster and field environment is negligible. Rapid evolution in the quiescent fraction is seen from to for massive galaxies suggesting a build-up of massive quiescent galaxies at . For galaxies with at , the quiescent fraction is found to be as much as a factor of 2 larger in clusters than in field, showing the importance of environmental quenching in low mass galaxies at low redshift. Most high mass galaxies are already quenched at , therefore environmental quenching does not play a significant role for them, although the environmental quenching efficiency is nearly identical between high and low mass galaxies.

galaxies: clusters: general — galaxies: high redshift — galaxies: evolution — galaxies: stellar content — galaxies: star formation

## 1. Introduction

Galaxies define basic separate, independent entities in the universe and are the building blocks of her. Therefore, the study of their evolution throughout the history of the universe is the very core in our understanding of the universe.

The cold dark matter (CDM) cosmological models say that galaxies form through the gravitational collapse of baryonic matters — mostly in the form of gas — inside dark matter (DM) halos. As gas collects within DM halos, it begins to cool and form stars, which gives galaxies the very glittering look as we observe them. An interesting question regarding star-forming (SF) activity of galaxies is how this activity evolves and what governs this activity. Unlike the collapse of DM, which is mainly governed by gravity, SF activity is governed or affected by various baryonic physical processes — including cooling of gas as well as feedbacks from stars, supernovae (SNe), or active galacitc nuclei (AGN) — thus making a detailed understanding of the SF activity complicated.

Regarding the evolution of SF activity of galaxies, we know that galaxies in the local universe have much lower levels of SF activity than the ones in the past. First, the global SF history, revealed from various galaxy surveys (Hopkins & Beacom, 2006; Cucciati et al., 2012; Behroozi et al., 2013; Burgarella et al., 2013; Magnelli et al., 2013; Madau & Dickinson, 2014) shows significant drop from redshift, to — with a possible broad peak around . While this shows the behavior of the collective SF activity in the universe, we can still speculate from this that the SF activity in individual galaxies must also decrease significantly. Second, star-formation rates (SFRs) of normal SF galaxies at are significantly lower than SFRs of higher redshift () SF galaxies with similar mass (e.g., Daddi et al., 2007; Sobral et al., 2014). This suggests an average drop in SF activities of individual galaxies with time toward . Lastly, a large fraction of local galaxies — especially massive ones — are red, quiescent, and form a tight red sequence (e.g. Strateva et al., 2001; Blanton et al., 2003; Balogh et al., 2004). Therefore, we need to understand the physical mechanism(s) that can lower and eventually stop the SF activity of galaxies to explain this change of average SF properties of galaxies from high redshift to the local universe.

At low redshift, it has been well known that SFR, color and morphology of galaxies, show strong dependence on their mass (e.g. Kauffmann et al., 2003; Jimenez et al., 2005; Baldry et al., 2006), in such a way that less massive galaxies show more SF activities than massive ones. This clear mass-dependence seems to persist at higher redshift. By analysing the red-sequence luminosity function of a large sample () of galaxy clusters at , Gilbank et al. (2008) found a dearth of red galaxies at faint luminosity at higher redshift. Recently, using the UltraVISTA (McCracken et al., 2012) data, Ilbert et al. (2013) have shown the stellar-mass dependence in SF quenching at in a sense that stellar mass density of massive () quiescent galaxies remains unchanged while less massive galaxies keep being quenched — i.e., massive ones already quenched before , earlier than less massive ones. The origin of this mass-dependent SF cessation or quenching is still unclear. It could be the effect of feedback from AGN (Hopkins et al., 2006; Somerville et al., 2008), which is believed to selectively act for more massive galaxies. Another compelling explanation is the heating of accreted cold gas in massive halos (Birnboim & Dekel, 2003).

On the other hand, the SF activity — or related properties, such as color or morphology — shows environmental dependence as well. Locally, galaxies show distinct SF properties and colors in different environment (Lewis et al., 2002; Kauffmann et al., 2004; Blanton et al., 2005). Naturally arising questions are how and when these environmental trends have been developed. By analysing galaxy colors out to from VIMOS-VLT Deep Survey (VVDS), Cucciati et al. (2006) found that the color–density relation progressively weakens with increasing redshift and possibly reverses at . Grützbauch et al. (2011), based on the Palomar Observatory Wide-Field Infrared (POWIR; Conselice et al., 2008) survey data in the DEEP2 field, also argued that there is no or weak correlation between galaxies’ local number density and the color or the blue galaxy fraction for galaxies. In the case of the relation between SFR and density, Cooper et al. (2008) found the reversal of the relation between SFR and overdensity at using galaxy samples from the DEEP2 Galaxy Redshift Survey (GRS). With GOODS data, Elbaz et al. (2007) also argued that the relation between SFR and density was reversed at , showing a sharp contrast to the local trend.

On the contrary, Cooper et al. (2006) argued that the mean galaxy environment shows a strong dependence on galaxy color at , similarly with what is locally found. Patel et al. (2009), analysing galaxies, found that the total SFR–local-density relation still persists at this redshift. This environmental dependence of galaxy properties has been shown to persist even at higher redshift () by Quadri et al. (2012), through the analysis of quiescent galaxy fraction of UKIDSS/UDS galaxies, even though the environmental trend of quiescent fraction becomes weaker with redshift. Combining data from several surveys — UltraVISTA, 3D-HST (Brammer et al., 2012), Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS; Grogin et al., 2011; Koekemoer et al., 2011), and SDSS, Tal et al. (2014) have studied the evolution of quiescent galaxy fraction in groups over a wide range of redshift up to . Their finding is that quiescent fraction of group satellite galaxies is similar with background (i.e. field) galaxies at , while it increases faster than field galaxies at lower redshift down to . So, there seems to exist a certain degree of disagreement on how and how much the environmental trend changes as we go up to high redshift, even though there is a consensus on the fact that the environmental dependence weakens as redshift increases (e.g., Scoville et al., 2013). Therefore, we have not yet reached to a firm conclusion about the rate or amount of evolution in the effects of environment on the formation and evolution of the quiescent galaxy population.

In studying the role of environment in formation and evolution of quiescent galaxies, galaxy cluster offers a very unique laboratory. They define the densest galaxy environment, thus are expected to follow the highest density peaks in the universe. Besides the fact that denser structures collapse earlier in the CDM paradigm, galaxies in these dense environments would also be subject to several physical processes — such as ram-pressure stripping (Gunn & Gott, 1972), strangulation (Larson et al., 1980; Balogh et al., 2000), and harassment (Moore et al., 1998) — acting exclusively in group- or cluster-like environments. In this regard, studying the stellar population properties of galaxies in galaxy clusters at various redshifts is a useful and crucial test bed in the investigation of the environmental effects on the formation and evolution of quiescent galaxy population.

While there is a sharp contrast in galaxy properties between galaxies in clusters and in field locally, as the cluster environment is dominated by the red, quiescent galaxies (e.g., Dressler, 1980), this picture seems to change as we go out to high redshift (Butcher & Oemler, 1984; Dressler et al., 1997; Poggianti et al., 2001). Combining their field sample with the cluster samples from Holden et al. (2007), van der Wel et al. (2007) showed that the difference in the red galaxy fraction between the cluster- and the field environments increases with decreasing redshift — i.e., the red galaxy fraction increases toward low redshift in the cluster environment while it remains nearly constant in the field environment. At higher redshift, , the story seems to become more complicated: Studies suggest that the galaxy SF properties (the SFR or the color) in the cluster environment at can be different from those of cluster galaxies at lower redshift and that there can be cluster-to-cluster variation in the properties of these high-redshift proto-clusters. Gobat et al. (2011) analysed colors of red galaxies in a proto-cluster at originally identified by their red /IRAC colors, and find that these red galaxies have old stellar populations. On the other hand, Zirm et al. (2008) found a well developed red-sequence in a proto-cluster around a radio galaxy and suggested that many of these red galaxies are porbably dusty star-forming galaxies. Fassbender et al. (2011) also reported the existence of the actively star-forming galaxies in the X-ray selected proto-cluster at . Tran et al. (2010) studied the color and SF properties of the member galaxies of -selected cluster (Papovich et al., 2010). They also found that the galaxies in this high redshift cluster, unlike the ones in lower redshift, are more dominated by blue galaxies. This, combined with the higher star-forming galaxy fraction within this cluster than in the lower redshift (), indicates that many cluster galaxies are still actively forming stars at . The results of Hilton et al. (2010), based on the X-ray and mid-IR (MIR) observation, also show the strong SF activity of the cluster member galaxies at . Strazzullo et al. (2013) have found that there are massive star-forming galaxies along with passive ones even in the core region of the IRAC-selected proto-cluster studied by (Gobat et al., 2011). All together, these observational results point to the redshift range of as the active formation era of the galaxy clusters as well as the epoch when many of cluster galaxies are still forming stars actively.

Summarizing, both the mass and environment of galaxies seem to affect the quenching of SF activity of galaxies and the formation of quiescent galaxies. Also, the environmental dependence of SF activity of galaxies seems change with redshift. Therefore, more investigation about the timing (i.e., the redshift dependence) as well as the relative contribution from these distinct mechanisms — mass versus environment — is still needed. Peng et al. (2010) have provided interesting results in this regard. Using SDSS (York et al., 2000) and zCOSMOS (Lilly et al., 2007) data, they analysed the effects of mass as well as environment on the evolution and star-formation activity of galaxies at redshifts up to . They found that “mass quenching” is more dominant for high stellar-mass ( ) galaxies, while lower mass galaxies are more affected by “environmental quenching”. They also suggested that the combination of these two quenching mechanism can explain the evolution of stellar mass function of star-forming as well as quiescent galaxies. Analysing H emitters from the High- Emission Line Survey (HiZELS), Sobral et al. (2011) have shown that the median SFR depends on galaxies’ environment for galaxies with their stellar mass lower than , while there is no environmental dependence for more massive galaxies.

In this work, we extend this kind of investigation to higher redshift (up to ), examining the effects of these two suggested drivers — i.e. mass and environment — in the formation and the evolution of quiescent galaxy population. It has been known that a tight red sequence already exists at (e.g. Im et al., 2002; Bell et al., 2004). Therefore, investigation at higher () redshift is crucial in catching the formation of these quiescent galaxies to place a meaningful constraints on their evolution. Our work is based on the deep near-infrared (NIR) data from the United Kingdom Infrared Telescope (UKIRT) Infrared Deep Sky Survey (UKIDSS) Ultra Deep Survey (UDS) as well as deep optical data from the Subaru/XMM-Newton Deep Survey (SXDS), combined with mid-IR (MIR) data from the Space Telescope. Deep UDS NIR data are crucial in robust estimation of SF or stellar population properties of galaxies, breaking the degeneracy between old and dusty star-forming populations, thus in a study of quiescent galaxy population evolution at redshift as high as . Our analysis is based on (1) an un-biased cluster detection (unlike, for example, the red-sequence method), (2) a relatively long redshift-baseline () and a large spatial area ( Mpc at ), and (3) robust estimation of SFR and stellar mass from SED-fitting. As demonstrated below, deep NIR data in the UDS field enables us the robust estimation of the stellar population properties of the galaxies up to .

We present the data set used in this work as well as the sample selection procedure in Section 2. In Section 3, we explain how the photometric redshifts as well as the stellar population properties of the UDS galaxies are estimated. We explain the identifying procedure of the high-redshift clusters and provide the properties of the selected cluster candidates in Section 4. We analyse the stellar population properties and the evolutionary trends of the quiescent galaxies, focusing on the drivers of this population in Section 5, and we summarize our results in Section 6. We adopt the standard flat CDM cosmology, with () = (0.3,0.7), and = 70 . All magnitudes are given in AB magnitude system (Oke, 1974) except when otherwise mentioned.

## 2. Data and Sample

The UDS (Almaini et al. in prep.), which is one of the five surveys in the UKIDSS (Lawrence et al., 2007), provides a very deep NIR imaging dataset over an area of degree located on the SXDS. The survey was carried out with the Wide Field Camera (WFCAM; Casali et al., 2007) on the UKIRT in three NIR broad-bands with the 5- limits of , , and (in Vega, DR10).

The SXDS provides the deep optical data from the SUPRIMECAM on the Subaru telescope from -band through -band. The 3- depths of the Subaru data are , , , , and . We use the released DR1 (-band detected) catalog (Furusawa et al., 2008), which is obtained from the SXDS DR1 release page, for this work. The UDS region has also been observed by the /IRAC as the SpUDS Spitzer Legacy Survey (PI:Dunlop), reaching magnitude (Channels 1 and 2). We obtain the catalog at the Spitzer Legacy Survey archive, and aperture-corrected total magnitudes are used for /IRAC data. We also use the publicly available spectroscopic redshift data for about 4000 objects (Simpson et al. in prep, Akiyama et al. in prep. and Smail et al. 2008)

We perform photometry on the UDS DR10 images () using SExtractor (Bertin & Arnouts, 1996) software. We use convolution mask of a Gaussian PSF with FWHM of 5 pixels, and set the detection minimum area as 6 pixels. The detection and photometry are done in each band separately, and then the match between different bands is done with a matching radius of and with the -band image as the reference image. About objects (both in - and -band catalogs) have multiple matches, and the nearest one is selected as the best match. To construct the panchromatic spectral energy distributions (SEDs) from the optical to the MIR in the observed frame, we match, first, the Subaru optical and UDS NIR catalogs, and then with the SpUDS IRAC catalog, with a matching radius of and the position in -band catalog is used as the reference position. For optical and NIR data, we use the auto magnitudes, and the aperture-corrected total magnitudes are used for IRAC data.

We check the validity of using auto magnitude () in crowded region, such as clusters, by comparing auto magnitude values with aperture (circular diameter = ) magnitude () in field and cluster regions. If values are biased in cluster regions, then we expect that to deviate from differently in cluster and field. At several magnitude bins (), we compare the distributions of between cluster and field region through the Kolmogorov-Smirnov test (K-S test). At each magnitude bin, the maximum difference, , was 0.073, 0.086, 0.071, 0.058, and the probability of null hypothesis, p, was 0.93, 0.62, 0.64, 0.94, showing that there is no systematic bias in auto-magnitude measurements that is specific in cluster region.

In our analysis, we use objects which are detected in both of optical (Subaru) and NIR (UKIRT) catalogs as our sample, but regardless of IRAC detection. There are about of -band detected objects not detected in band. The fraction is lower as if we only consider the sources with . These -band undetected sources have, on average, redder IR colors than the detected ones. The mean values of and are 1.0 and 1.8, respectively, while the corresponding values are 0.6 and 1.0 for -band detected objects. And, about of - and -bands detected objects are matched with IRAC sources.

From this ‘Subaru+UKIRT+Spitzer’ multi-band photometric catalog, we first cull out stars either based on the spectroscopic classification or broad-band color (BzK or -). Specifically, from total objects, we exclude: (1) spectroscopically classified stars, (2) stellar candidates based on their location in the (-)–(-) color–color plane, with the color criterion of , similarly with Daddi et al. (2004) or Hanami et al. (2012), and (3) further candidates with and , based on the fact that galaxies and stars form clearly separate branches in this versus plane.

We need to cull out active galactic nuclei (AGN) from our sample because the measurements of photometric redshift and stellar population properties, which are crucial in our study, would be affected by the AGN activity. To remove AGN candidates from this star-removed, multi-band catalog, we use three steps: First, we exclude 77 objects which are classified as quasars in the spectroscopic catalogs. Next, we match our multiband catalog with the XMM-Newton point source catalog (Ueda et al., 2008) to remove any matched object. At this step, about 220 objects are matched and excluded. Finally, for the objects with the matched IRAC photometry, we use the NIR MIR color criteria suggested by Messias et al. (2012) to find out any remaining AGN candidates. The applied criteria are and , where, and are the magnitudes in IRAC channel 2 () and channel 4 (), respectively. About 7,000 objects satisfy these color criteria, thus are excluded from the sample. These color criteria, similarly with other IRAC color criteria, cannot distinguish between AGN and high-redshift () galaxies. Therefore, by applying these criteria, we might lose some non-active galaxies at as well. However, this does not affect our results, because our target redshift range does not exceed .

After trimming out AGN candidates as well as stars as explained, we apply the -band magnitude cut (, 5 detection limit) and the number of remaining galaxies is about 115,000 among which 1,400 objects have spectroscopic redshift information. For these 115,000 objects, we derive photometric redshifts as well as stellar population properties through SED-fitting (Section 3). After selecting galaxies within the redshift range of and applying the stellar mass cut of log — which is the stellar mass limit at the highest redshift bin () of our study — our final sample size is 46,641.

## 3. Estimation of Photometric Redshifts and Stellar Population Properties

### 3.1. Photometric Redshift

We derive photometric redshifts using the EAZY code (Brammer et al., 2008). We take the -band magnitude prior, which slightly improves the photometric redshift accuracy. We fit ten broad-band photometric data points from the observed frame optical (Subaru) to MIR (/IRAC). During the fitting procedure, we require that at least five colors to be measured — i.e. it is detected at least at six bands. Therefore, any galaxy which is not detected at more than five bands is not included in the following analysis. Only of objects are excluded for this reason, and these are mostly faint objects — for among these, and for .

For a subset of galaxies for which spectroscopic redshifts are available and , we compare the spectroscopic and the photometric redshifts to estimate the accuracy of the photometric redshift (Figure 1, left panel). In UDS, there are not enough galaxies with for a fair estimation of photometric accuracy. To test the photometric redshift accuracy at , we use the GOODS-S data, which carries similar coverage of broad-bands data with the UDS while having more abundant spectroscopic redshift samples at . These GOODS-S spectroscopic data are from various spectroscopic surveys (Le Fèvre et al., 2004; Szokoly et al., 2004; Mignoli et al., 2005; Cimatti et al., 2008; Vanzella et al., 2008; Popesso et al., 2009; Balestra et al., 2010; Silverman et al., 2010). We measure signal-to-noise of UDS photometric data as a function of magnitude at each band, and add additional noises (randomly following a Gaussian distribution) to the GOODS photometry to simulate the depth of the UDS data. Then, photometric redshifts are derived using these noise-added photometry. Within the redshift range of , the redshift range of interest of this work, there are of objects whose redshift error, defined as , is greater than 0.15. Removing these objects with catastrophic redshift errors, the mean error is 0.028.

The right panels of Figure 1 show as well as the outlier fraction as functions of photometric redshift. We can see that our photometric redshift estimation is reliable within the redshift range of our study, , but its uncertainty is larger at . The outlier fraction increases at high redshift (), but it is kept at . The reason for the increase in the photometric redshift uncertainty at is that the spectral break at Å  — one of the clearest redshift indicators — moves out of the optical regime at . When that happens, it becomes more difficult to sample the exact location of the rest-frame 4000 Å  break in the observed frame due to the gaps between NIR filters and the reduced number of available filters above 4000 Å  break. Also, the fraction of star forming galaxies are higher at in the spectroscopic sample since it is easier to spot emission lines than absorption lines for faint galaxy spectra. Weak 4000 Å  breaks and strong emission lines of star forming galaxies make it difficult to trace the continuum shape, leading to a reduced accuracy in photometric redshifts (e.g., Damen et al., 2009; Ilbert et al., 2009; Quadri et al., 2012; Hartley et al., 2014; Yang et al., 2014; Kang & Im, 2015). Not surprisingly, photometric redshift outliers tend to be star-forming galaxies, and some of them can move into the red sequence due to the large errors in photometric redshift.

In Figure 2, we show the distribution of galaxies with several bins ( and 1.6) with . Figure 2 shows that systematic offsets between photometric and spectroscopic redshifts are small ( at some redshift bins). The effect of the small systematic shift in redshift is negligible in our analysis. For example, a systematic offset of causes a shift in the stellar mass estimate by .

### 3.2. Stellar Population Properties from SED-fitting

We perform the stellar population analysis through spectral energy distribution (SED) fitting methods for our sample, with the same procedure as explained in detail in Lee et al. (2009, 2010). For redshift, we adopt if available, but otherwise.

The stellar population synthesis models of Bruzual & Charlot (2003, BC03, hereafter) with Padova 1994 stellar evolutionary tracks are used for our SED-fitting analysis. We assume the Chabrier (2003) initial mass function (IMF) with the lower and upper mass cuts of 0.1 and 100 M.

As for the star-formation histories (SFHs), we use the delayed SFHs (Lee et al., 2010). This form of SFH is shown to provide better estimation of stellar population properties than other forms, for example, exponentially declining ones (e.g. Maraston et al., 2010; Papovich et al., 2011; Pforr et al., 2012; Lee et al., 2014). Studying galaxies from DEEP2 Redshift Survey (Newman et al., 2013), Pacifici et al. (2013) have also shown that these galaxies have SFHs similar with this form.

It has the following functional form.

 Ψ(t,τ)∝tτ2e−t/τ, (1)

where, is the instantaneous SFR for a given set of and . The parameter is the time since the onset of the star-formation and is the time-scale parameter which governs how fast (or slowly) the SFR reaches its peak value before starting to decline. We allow the value of to vary within a very broad range from Gyr to 10.0 Gyr, and within 0.1 Gyr , where is the age of the Universe at redshift .

After extracting spectral templates from BC03 with these allowed values of parameters, we apply the Calzetti et al. (2000) attenuation law to model the dust attenuation of SEDs due to the inter-stellar dust. We vary the amount of attenuation, parametrized via within 0.0 1.50 with the increment of 0.025. For the attenuation due to the line-of-sight neutral hydrogen in the inter-galactic medium (IGM), we apply the extinction law of Madau (1995). During the fitting, the redshift is fixed at the value of the spectroscopic redshift if available, otherwise at the photometric redshift.

We apply the ‘grid-scanning’ method which has been used in several previous works (Wiklind et al., 2008; Lee et al., 2009, 2010, 2014). As explained in these references, our method of skimming the entire parameter space is less subject to the effects of any prior. From this SED-fitting results, we derive for each galaxy, the stellar mass (), the SFR, the mean age of the stellar population, the amount of dust attenuation (), as well as the star-formation history — parametrized in terms of and . The SFR used in this work is defined as the SFR averaged over recent 100 Myr, based on the reasoning outlined in Lee et al. (2009). We also derive the rest-frame color — will be denoted as from now on — for each galaxy from its best-fit BC03 template.

The broad wavelength coverage out to IR regime as well as the deep optical data from Subaru in the UDS field enables us to estimate the stellar population properties of high-redshift galaxies with high confidence and distinguish dusty SF galaxies from old galaxies among red galaxies. Figure 3 shows the examples of the SED-fitting results for a representative set of galaxies at various redshifts. Also shown are the key SED-fit parameters. In this figure, we show the observed galaxy SEDs (red circles) as well as the best-fit BC03 spectra for typical blue SF galaxies (panels in left column), dusty SF ones (panels in middle column) and old quiescent ones (panels in right column) at redshifts , 1, 1.5.

Here, we define “quiescent galaxies” as galaxies with sSFR ] yr, where is the age of the universe at . This cut was used in several previous works (e.g., Damen et al., 2009), and approaches the local value for sSFR cut at (e.g., Gallazzi et al., 2009; Ko et al., 2014). Under a delayed SFH model (i.e., Equation (1)), 1/[] corresponds to of the peak SFR. While this sSFR cut select galaxies with low SF activity in comparison to the average SF galaxies at that redshift, the selected galaxies can have SFRs that are not quite “quiescent”. For example, at , galaxies at the sSFR cut have /yr — i.e., SFR of starburst galaxies in the local universe. If we fit the SEDs of a 1 Gyr-old, passively evolving galaxies after a spontaneous burst (a truly quiescent galaxy at high redshift), our procedure returns sSFR yr. We check the MIPS 24 m detection of our galaxies using SpUDS MIPS catalog. We find that small fraction () of red SF galaxies is detected, while no quiescent galaxy is detected at 24 m.

We test how the error in photometric redshift measurement may affect the estimation of stellar masses of galaxies using Monte Carlo simulation. Specifically, first, we scatter the redshift of each galaxy assuming a Gaussian distribution of redshift error. Then we repeat SED-fitting by fixing the redshift at this scattered value. We repeat this procedure 100 times each galaxy at and , within a redshift range . The results of this simulation are given in Figure 4, where we show the number counts of galaxies at stellar mass bins. The resulting number distributions with scattered redshift (red) show no significant difference from the ones without scatter (blue) at both redshift bins. For individual galaxies, the mean and the standard deviation of is and , at and , respectively. Here, is , and is the measured stellar mass without redshift scatter. From this, we can ensure that the stellar mass is not significantly affected by the error in photometric redshift measurement.

## 4. Massive Structures of Galaxies

### 4.1. Stellar Mass Limit

Based on the stellar mass derived from the SED-fitting, we select galaxies with their stellar masses greater than (which is the stellar mass limit at the highest redshift range, ) within the redshift range of .

We have checked if the magnitude limit is deep enough to detect galaxies with stellar mass of out to , using galaxy template spectra of various galaxy types extracted from the BC03 library. In Figure 5, we show these spectra — young (solid blue curve), dusty young (solid green curve), and old (solid red curve) — with stellar mass of and redshifted to . We also show the limiting flux of the UDS data set in 5 Subaru optical filters (, , , , and ), in 3 UKIRT filters (, , and ) and in 2 IRAC channels ( and ) shown as black horizontal bars. We find that SEDs of all the galaxy types are above the detection limits in multiple bands, demonstrating that our dataset is nearly suitable for selecting galaxies with log () .

Figure 6 shows the redshift versus stellar mass of galaxies in the -band limited sample. In this figure, the green vertical lines show the sample redshift range of this study and the green horizontal lines show the stellar mass limits. Blue curve shows the stellar masses at different redshifts of star-forming galaxy templates from BC03 with formation redshift of . Yellow curve corresponds to quiescent galaxy templates with , which is a reasonable formation redshift for faint/low-mass galaxies, while red curve is for galaxy templates with extremely high formation redshift (). As shown in this figure, within the given redshift range, we can detect log () galaxies (corresponding lower green horizontal line) either quiescent (yellow curve) or star-forming (blue curve). It is still possible that we may miss some low mass () galaxies with extremely old () population, even though it is reasonable to assume that there are only few, if any, galaxies with these extreme properties in reality. For this reason, we apply the stellar mass cut of (upper green horizontal line) when we deal with quantities which might be affected by this missing of low-mass old galaxies. Within this redshift range and above the stellar mass limit (log () = 9.1), we have 46,641 galaxies.

### 4.2. Method and Procedure

Our method of finding galaxy cluster candidates is based on the photometric redshift and similar with the method used in Kang & Im (2009, 2015). Following is the detailed explanation of the procedure we use to identify the cluster candidates. First, we divide the projected sky area into retangular grids with a width of , which corresponds to , 102, and 100 kpc at , 1.5, and 2. Then, at a given redshift bin and at each grid point, we count the galaxies with log ( within a radius of kpc ( at ) from the grid point within the redshift range of the typical redshift error in the UDS (i.e., ). We repeat this procedure at redshift bins with an increment of = 0.02.

Next, we find the mean () and the standard deviation () of the galaxy number counts () through the Gaussian fitting at each redshift bin. Then, we identify the - grid points with their galaxy number counts that exceed the 4- level from the mean, i.e.,

 N ≥ ¯N+4×σN. (2)

If we see the spatial distribution of these - grid points with extreme galaxy number counts (), we find that some points are connected to other points making large structure while some points are relatively isolated. In our case, the spacing between the grid points ( kpc) is much smaller than the typical size of galaxy clusters ( Mpc). Therefore, a cluster of galaxies would appear as a structure with connected - grid points. For example, if 9 grid points are connected with each other with a square shape (i.e., ) or if 10 grid points are connected in a rectangular shape (i.e., ), this connected structure would represent an overdense structure with a radius of kpc. Based on this reason, we exclude any - grid points when the number of connected grid points is smaller than 10 at each redshift bin. Then, among the remaining connected structures, we only select the structures as our final (proto-)cluster candidates only when these structures are identified at the same sky location (within the displacement, Mpc) in at least three successive redshift bins. In other words, if the connected structure is found only in one or two redshift bin(s), we discard this structure. This criterion has been applied to minimize any false detection of cluster candidate. By applying these two, rather conservative criteria in selecting cluster candidates, we try to make our selection of cluster sample as robust as possible. By comparing with a mock catalog from the GALFORM simulation (Merson et al., 2013), we find that our cluster finding method is nearly complete at selecting clusters at and . A detailed description of this method applied on simulation data can be found in Kang & Im (2015).

### 4.3. Galaxy Cluster Candidates

Following the procedure as explained in the previous section, we find 46 cluster candidates within the redshift range of . Among these, 19 clusters are the ones already identified by other authors (van Breukelen et al., 2006, 2007; Finoguenov et al., 2010; Papovich et al., 2010) in the same UDS-field, and 27 clusters are newly identified in this work. Figure 7 shows the number of the cluster candidates at each redshift bin.

In Table 1, we present the basic properties of these cluster candidates, including the number of member galaxies (Column 4), total stellar mass of member galaxies (Column 5), and the overdensity measure, defined as (Column 6).

After finding the cluster candidates, we define the spatial center and the redshift of the candidate clusters through the following procedure. At each redshift bin, the mean position, weighted by the galaxy count (), is found for each connected over-dense structure. Next, the center and the redshift of each cluster candidate is assigned in a similar manner, i.e, as the -weighted mean center and redshift, for each of the connected redshift bins. Then, we find the galaxies within the radius, Mpc — which corresponds to, or is slightly smaller than, a typical radius of a galaxy cluster — from the center and within the redshift interval, — which is the typical redshift error of our data. Finally, we re-calculate the center and the redshift of each cluster candidate as the stellar-mass weighted mean center and redshift of the member galaxies belonging to the overdense area. These center and redshift are the values listed in Columns (1), (2) and (3) in Table 1. The number of member galaxies varies within the range between 30 and 90 (the left panel of Figure 8). The right panel of the same figure shows the over-density measure, defined as ()/, and shown as (-axis) in the figure. Most clusters have values between 4 and 8, and 4 clusters at redshift have higher .

Out of 46 clusters, 13 clusters are also found by Finoguenov et al. (2010) based on the -ray observation, and they provide the total mass () measured from the -ray data, and the corresponding . We compare the sum of the stellar masses of the cluster member galaxies and their total (halo) mass for these 13 clusters, and these two mass measures — halo mass () and total stellar mass () — show a positive correlation as shown in the left pane of Figure 9. In the middle panel of this figure, we show the redshift versus the mass ratio () of clusters. Here, we cannot see any clear redshift-dependent trend, and the mean mass ratio is 0.013, shown as a blue horizontal line in this panel. The right panel shows the redshift versus of our 46 clusters. Here, red large circles are the 13 clusters with their taken from Finoguenov et al. (2010). Blue small circles are the remaining 33 clusters, for which we assign a halo mass assuming the stellar-to-halo mass ratio of 0.013. These clusters have halo masses of to , and total stellar masses of to .

## 5. Color and Star Formation Properties

In this section, we take a close look at the color and SF properties of galaxies belonging to the cluster candidates, compare these properties with those of field galaxies in the same UDS field, and investigate their evolutionary trend.

### 5.1. Optical Color and its Evolution

#### 5.1.1 Color of Cluster and Field Galaxies

Figure 10 shows color versus stellar mass diagrams of member galaxies of each cluster candidate. As shown in this figure, the red sequence is formed throughout the whole redshift range, while the fraction of the galaxies with red color decreases with increasing redshift. Hereafter, we refer the galaxies with as red galaxies, and the ones with as blue galaxies. In clusters, a significant of red galaxies starts to appear at redshift as high as .

At the very highest redshift (), only several massive galaxies have red () colors while progressively lower mass galaxies join the red sequence with decreasing redshift. At , we can see that the red sequence is well extended down to very low stellar mass (log ) in most of the clusters. This is a version of downsizing phenomenon, which is also shown in the study of two galaxy clusters (at and ) by Nantais et al. (2013). We confirm this trend of mass-dependent timing of red sequence formation with a significantly larger sample. A similar, luminosity-dependent trend in a deficit of red galaxies was suggested earlier by Gilbank et al. (2008). Please note that we can detect red galaxies down to a low-mass limit () throughout the redshift range (), thanks to the very deep optical data from SXDS. Therefore, the dearth of low-mass red galaxies at highest redshift bins is real — i.e., not affected by the magnitude limit of the data.

In Figure 12, we present the distributions of color at five redshift bins (, , , , and ) for cluster (top row) and field galaxies (bottom row). As can be seen in the color- diagrams, color shows a clear bimodal distribution, which suggests that the color transition of galaxies occurs in a short time scale. At redshifts greater than , the clusters are dominated by blue galaxies, while at , clusters are - — i.e., majority of the cluster galaxies are red at . But, red galaxies exist even at the highest redshift bin (). We can also observe the reddening of the peak color of blue galaxies with decreasing redshift, which reflects the evolution of the average stellar population of blue galaxies with time.

When we compare the color distributions of the cluster galaxies (top row) and of the field galaxies (bottom row), we can see that the color distribution is similar in clusters and in field at the highest redshift bin (). At , the color distribution in clusters becomes distinguished from that of the field galaxies. Clusters become to be more - compared to field at the same redshift. At the lowest redshift bin, blue galaxies are still the major population in field, while red galaxies clearly dominate in clusters. From this, we can conclude that the color transition of galaxies are accelerated in clusters at redshift lower than .

#### 5.1.2 Color Evolution in Clusters and Field

In Figure 13, we show the redshift evolution of the median color of blue galaxies, both for clusters and field. In this figure, the small green circles show the median color of blue galaxies and the redshift for each cluster candidate, and the red, large circles show the median colors and the median redshifts of blue cluster galaxies summed in each redshift bin. The blue diamonds represent the median values of field blue galaxies at each redshift bin. The corresponding SDSS value, which is derived from the MPA-JHU catalog, is shown as the magenta circle. The expected color evolution of galaxies from the BC03 model with various (delayed) SFHs are also shown as the green dotted lines ( Gyr, ), blue dashed lines ( Gyr, ), and red solid lines ( Gyr, ). Here, is the formation redshift. For each color, the lines from the bottom to the top are with the increasing amount of dust attenuation from with the increment with 0.1.

As already mentioned, the blue peak evolves to become redder with decreasing redshift, reflecting the average aging of the stellar population in the blue galaxies. Apparently, the BC03 models with Gyr, (with and 0.1 — i.e., upper two red solid lines) or the models with Gyr, (blue dashed curves) seem to well bracket the observational trend of blue galaxies. This does not necessarily mean that the individual galaxies would evolve with this SFH, even though it can be a representation of global SFH within this redshift range. Interestingly, these SFHs have a peak of SFR at or 2, showing a coincidence with the peak of the global SFR density evolution. It should be noted that the difference between cluster- and field-galaxies is not significant, considering the associated ranges of error bars and the fact that median color can be affected by non-negligible fraction of green galaxies.

In interpreting this color evolution of blue galaxies, we should consider the fact that the blue galaxy population is not a closed one. As galaxies evolve to become redder, these galaxies would move out to the red galaxy population. Also, low mass blue galaxies that were originally out of our mass-cut would be progressively included into the blue galaxy population as their stellar mass grows. The effect of these is to prevent the median color of the blue galaxy population from being reddened quickly. Therefore, it is very likely that the SFRs of the blue galaxies would evolve more rapidly (i.e., having smaller value of ) than the one represented by the red curves in Figure 13. Also, we cannot rule out the possibility of the change in the median dust attenuation values of the blue galaxies with redshift. The distributions of show little, if any, evolution within the redshift range, , with a hint of slight decrease of the mean at the highest redshift bins (). If we accept this small amount of evolution of dust attenuation, it means that the higher dust attenuation at lower redshift contributes, to some extent, to the redder color (of the blue galaxies) at lower redshift, requiring larger value of . This will compensate, to some extent, the effects of migration of fading galaxies out of the blue population, making our estimation of the (global representative) SFH more robust.

The evolution of the color of quiescent galaxies is shown in Figure 14. Unlike blue galaxies, the quiescent galaxy population show little color evolution within the redshift range of either in clusters or in field. This no-evolution of color reflects the fact that quiescent galaxies (with the colors bluer than the quiescent ones) keep joining the quiescent population while the color of the quiescent galaxies becomes redder with time. The dashed curves with various colors in Figure 14 are the evolutionary path of BC03 quiescent galaxies with various formation redshift. These curves are shown to guide how the colors would evolve if no additional quiescent galaxies with bluer color join the existing quiescent population. Purple circles are the corresponding SDSS colors in various stellar mass bins — log for the bottom one to log for the top one with an increment of 0.2 in logarithmic scale. Analysing red galaxy spectra at from the Deep Extragalactic Evolutionary Probe 2 (DEEP2), Schiavon et al. (2006) found relatively young ages ( Gyr) for these red galaxies with little or no emission lines. One of their interpretation is that blue SF galaxies keep joining this red population continuously (see also Harker et al., 2006; Faber et al., 2007) maintaining their mean stellar population ages young.

#### 5.1.3 Red Star-forming Galaxies

The color and sSFR of galaxies show a broad correlation in a sense that redder galaxies have in general lower sSFR. However, there are a non-negligible fraction of galaxies whose color is red () but are still forming stars actively. As shown in Figure 15, these red star-forming (SF) galaxies are, on average, dustier, older, and have lower sSFR than blue SF galaxies. No red SF galaxy has sSFR yr to be classified as ‘starbursts’. These indicate that red SF galaxies are in the stage, and migrating into the red quiescent population.

We investigate the morphological properties of red SF galaxies and compare them with the other galaxy populations to see how many red SF galaxies could be starburst galaxies resulting from gas-rich (wet) merger. A portion of the UDS region () is covered by the CANDELS (Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey; Grogin et al., 2011; Koekemoer et al., 2011) program, and observed with the WFC3 (Wide Field Camera 3) on board the . From the deep WFC3 H-band (F160W) image, which reaches to the 5- magnitude limit of 27.45 with FWHM of (Galametz et al., 2013), the Gini coefficients, , as well as the parameters of the galaxies are measured. Figure 16 shows the distributions of the Gini coefficients and the values of the blue, the red SF, and the red quiescent galaxies. The dividing lines between the early-, late-type and major mergers are from Lotz et al. (2008). There are only a few merger candidates among red SF galaxies (in the middle panels, represented as green symbols) and most of the red SF galaxies lie around the boundary between blue (left panels) and red quiescent galaxies (right panels), supporting the idea that these red SF galaxies are in transition phase from blue galaxies to red quiescent ones.

We also visually inspect the /WFC3 images of our galaxies. The blue SF population mostly consists of extended objects with disk-like morphology with small fraction of disturbed ones hinting the recent experience of merger or interaction. On the other hand, red SF galaxies are mixture of disky, spherical, and compact objects, reinforcing that these objects are in transition from blue galaxies into red quiescent ones. In Figure 17, we show example /WFC3 -band images of cluster galaxies at .

### 5.2. SF properties and Evolution of Cluster and Field Galaxies

In previous section, we investigated the rest-frame optical color properties of galaxies in clusters and field. Now, in this section, we concentrate on the SF properties of galaxies and analyse how the SF properties of galaxies within the high-redshift clusters evolve with time and also compare this with that of the field galaxies in the same redshift range.

In Figure 18, we compare the distributions of sSFR of the cluster- and the field galaxies at five redshift bins with the same redshift binning as in the Figure 12. In this figure, we assign sSFR = yr to all quiescent galaxies with sSFR yr. Similarly with color distribution, the sSFR distribution of SF galaxies in clusters and in the field are similar at high redshift bins (), while it shows a clear difference from that of the field galaxies at the lowest redshift bin (). For the quiescent galaxy fraction (over SF ones), the difference between clusters and field starts to show up from . While the evolutionary trends of sSFR distributions is similar with that of the color distributions shown in Figure 12, the sSFR distribution changes more slowly (especially for the cluster galaxies) than the color distribution. For example, at , the color distributions are already different between cluster and field, while the sSFR distributions of SF galaxies still remain similar in field- and cluster-environments.

Next, we compare the evolution of the quiescent galaxy fraction in cluster and in field in Figure 19, dividing the sample into three stellar-mass bins : log (left panel), log (middle), and log (right). In this figure, the filled blue diamonds represent the quiescent fraction of field galaxies at each redshift bin, while the quiescent fraction of individual clusters is shown as magenta dots. At several redshift bins, we sum the cluster galaxies in each -bin, and calculate the mean and the standard deviation. This is shown as the red diamonds with error bar in the figure. The red and blue circles are for the SDSS (Sloan Digital Sky Survey) galaxies from Baldry et al. (2006). We derive these values from their Equation (9), and choose the values with as field values and as cluster values in that equation. And, we apply the correction for red SF galaxies based on Haines et al. (2008).

From Figure 19, we find several interesting aspects in the evolution of quiescent galaxy fraction. First, the increase of the quiescent fraction slows down from the redshift for galaxies with their stellar mass, log , more significantly for the field galaxies, while it evolves fast from down to for massive galaxies (). At the two high stellar-mass bins, the quiescent fraction remains nearly unchanged from both in field and in cluster. Comparison with the local values (the red and the blue circles) shows that the quiescent fraction has already reached to the local value for the most massive galaxies (log ) at both in clusters and in field (the right panel), while the growth of quiescent fraction must be accelerated at to match the SDSS values in the case of low mass galaxies with log . This near-constant quiescent fraction of massive galaxies is in good agreement with the previous studies (e.g., Im et al., 2002; Ilbert et al., 2013). This change in the increase of quiescent fraction indicates that the redshift range defines the era of the rapid build-up of massive quiescent galaxies, which is in broad agreement with the finding of actively star-forming galaxy clusters at by several authors (e.g., Tran et al. 2010; Zeimann et al. 2012; Santos et al. 2014).

Second, the fraction of quiescent galaxies is higher in clusters than in field at for galaxies with , even though the scatter among the individual cluster candidates is quite large. From this, we can speculate that the environmental quenching becomes to work more strongly at redshift, , but with a certain amount of cluster-to-cluster variation. Also, we can see that the difference in the quiescent galaxy fraction between the cluster and the field environments is more significant for low-mass (log ) galaxies (the left panel) at redshift .

Last, as can be seen clearly in this figure, the quiescent fraction is a strong function of the stellar mass rather than their environment. At the highest mass bin, the -, which we define as the redshift or the epoch when the quiescent fraction starts to exceed , occurs already at . In sharp contrast, in the lowest mass bin, the quiescent fraction never reaches the - down to both in field and in cluster. Actually, the quiescent fraction is throughout the redshift range in field. This strong mass-dependence of quiescent fraction evolution well agrees with the recent results of Moustakas et al. (2013) — who found a nearly constant number density of massive quiescent galaxies at while the corresponding value rises rapidly with decreasing redshift for less massive galaxies, and of Huang et al. (2013). This strong mass dependence and the weaker environmental dependence dictates that the stellar mass plays a more dominant role in shaping the SFH of galaxies than the cluster-specific processes.

To gain an insight how rapidly (or how slowly) star formation is quenched for galaxies with different stellar masses, we devise a simple model explaining the evolution of star-forming galaxy fraction, , as follows.

 fsf(t)=fsf(0)×(@frac12)t/γ, (3)

where, is the fraction of SF galaxies at an initial time . In this model, is a - of SF galaxy population, which means that the fraction of SF galaxies becomes half of its initial values after the time . This simplified model assumes that the total number of galaxies does not change with time — i.e., no galaxy is added to or removed from the sample, and no merger occurs between galaxies.

Then, the quiescent fraction, would increase as,

 fquies.(t)=1−fsf(0)×(@frac12)t/γ. (4)

We fit our data points to this simple model, and we show this for the best-fit values of ’s as solid curves in each panel in the upper row of Figure 19. In each panel, the red and the blue curves are for cluster and field galaxies, respectively. For cluster galaxies at (solid red curves), the best-fit values for are , and Gyr for the lowest, middle and highest stellar-mass bins. In the case of field galaxies (solid blue curves), the corresponding values are , , and Gyr, each. This again demonstrates that: (1) massive galaxies become quiescent more rapidly both in clusters and in field (i.e., having smaller value of ) than less massive galaxies, and (2) star formation is more rapidly quenched in clusters than in field.

The lower panels of Figure 19 show the ratio of the quiescent fraction in clusters () to that in the field (). Here, we take an average value of the highest three redshift bins (shown as the black circles) and show it as a red diamond in each panel. In the two high mass bins (middle and right panel), the fraction is always with little evolution. However, this fraction increases rapidly with decreasing redshift for low mass galaxies (log ), and is higher than that in the higher mass bins at . This also indicates that the environmental quenching becomes more significant for low mass galaxies at .

To see the effect of stellar mass in determining the quiescent fraction more clearly, we show the stellar-mass dependent quiescent fraction at three redshift bins in Figure 20. In this figure, we can clearly see the strong stellar-mass dependence of the quiescent galaxy fraction at all three redshift bins as well as a clear excess in the quiescent galaxy fraction in clusters compared to the field for low-mass galaxies (log ) at the lowest redshift bin (left panel).

Summarizing, our investigation of the quiescent galaxy fraction and its evolution reveals several important aspects about quiescent galaxy formation: (1) At , quiescent, massive galaxies were built up rapidly, with a similar rate in both clusters and field. By z   1.3, most of the massive,quiescent galaxies were built up. (2) At , the environmental dependence of quiescent fraction become visible. The difference between cluster- and field-environment is clearer for low-mass (log ) galaxies. (3) Stellar mass plays more dominant role in determining the quiescent fraction than the environment throughout the entire redshift range.

### 5.3. Quenching Efficiency

In the previous section, we have shown that the quiescent galaxy fraction shows a clear stellar-mass–dependent trend, in a sense that the fraction of quiescent galaxies is higher for more massive galaxies throughout the redshift range, (upper panels of Figure 19), and also that the excess of the quiescent galaxy fraction in clusters over field environment is higher for low-mass (M) galaxies, especially at redshifts lower than (lower panels of Figure 19). Now, we investigate the environmental quenching efficiency (e.g. van den Bosch et al., 2008; Peng et al., 2010; Quadri et al., 2012). We define this environmental quenching efficiency as,

 (fq,cluster−fq,field)/fsf,field, (5)

where and are the fraction of quiescent galaxies in clusters and in field, respectively, and is the fraction of star-forming galaxies in field (). This quantity measures the fraction of SF galaxies in field that would have become quiescent if they were in cluster. Figure 21 shows the evolution of the environmental quenching efficiency in different stellar mass bins. Our results shows no clear stellar-mass dependent trend — even though the scatter is larger in the highest mass bin due to the small number of massive galaxies and the small values in this mass bin — confirming previous results (Peng et al., 2010; Quadri et al., 2012) but extending the probed redshift range to .

How can this result — no significant mass dependence in environmental quenching efficiency — be reconciled with the results in the previous section (e.g., Figure 19) — the high excess of the quiescent fraction in clusters in the lowest mass bin? This difference arises because there is another quenching mechanism — i.e., mass quenching. Massive galaxies are affected strongly by mass quenching, and they are mostly quiescent already at , leaving little room for the environmental quenching to play a significant role. On the other hand, the lower mass galaxies are quenched only mainly through the environment effects that becomes significant at .

## 6. Summary and Conclusion

In this paper, using deep optical to MIR data in the UKIDSS/UDS field, we have found 46 high-redshift galaxy clusters up to , among which 27 are newly found. We analyse the stellar population properties, such as color and SFR, of galaxies both in cluster and field environments over a wide redshift range to understand the effects of cluster environment on the galaxy evolution.

Through this analysis, we have found that the quiescent galaxy fraction increases rapidly with decreasing redshift at , and the increase is slowed down at lower redshift for massive (log ) galaxies. This trend points to the redshift range as the era of the rapid build-up of massive quiescent galaxies — the epoch when many massive galaxies stop their star formation and become quiescent. The difference in the quiescent galaxy fraction between clusters and the field increases at -1.4, which coincides with the epoch when the increase in the quiescent galaxy fraction begins to slow down. In this redshift range (), the quiescent galaxy fraction remains nearly unchanged in field, while it keeps increasing in clusters (but more slowly than at ). This difference is only significant for low mass galaxies with log () .

The environmental quenching efficiency shows no clear stellar-mass dependence, which is in agreement with previous results (e.g., Peng et al., 2010; Quadri et al., 2012), while our results extend this up to higher redshift (). However, the effects of environmental quenching appears more significantly for low-mass galaxies, because massive galaxies are affected by another quenching mechanism — i.e., mass-quenching — and most of these massive galaxies are already quenched at .

At , the quiescent fraction of galaxies shows no clear dependence on their environment, in contrast to its strong dependence on stellar mass. This infers that the SFH of galaxies is mainly shaped by their stellar mass during the early phase of evolution. The effects of the cluster environment on determining the quiescent fraction of galaxies seem to be more significant for less massive galaxies with , for which mass quenching does not yet affect their SF activity, at . The delayed appearance of the cluster-environment effects on the quiescent galaxy formation for low-mass galaxies implies that the cluster-specific processes which are mainly responsible for the excess of the quiescent fraction in clusters compared to the field region may be gradual, time-taking processes — like the strangulation (Larson et al., 1980; Balogh et al., 2000) — rather than abrupt or violent processes — like ram-pressure stripping, for low-mass galaxies. The investigation of the morphology of red SF galaxies — which shows that the most transition from the blue SF galaxies to red quiescent ones in clusters does not affecting their morphological appearance significantly — supports this idea.

One interesting question is why the SFHs of galaxies are affected by their stellar mass. More specifically, why is the quiescent galaxy fraction higher for more massive galaxies at all environment over a wide range of redshift () and without strong environmental dependence? One possible answer is that the dynamical time-scale would be shorter for the galaxies with greater stellar masses than less massive galaxies. If the SFR of galaxies depends on this dynamical time-scale, in a sense that shorter dynamical time-scale leads to higher SFR — as assumed in many galaxy formation models (e.g., Somerville et al., 2012), more massive galaxies will consume their gas more rapidly to arrive the red quiescent galaxy sequence earlier than low mass galaxies. Major mergers between massive SF galaxies can also accelerate this fast gas consumption. Another possibility is that any negative feedback which works preferentially for massive galaxies — such as AGN feedback (e.g. Hopkins et al., 2006; Somerville et al., 2008; Santos et al., 2014) or halo mass quenching (Birnboim & Dekel, 2003) — is the driver of the mass-quenching phenomena.

One of the interesting results of this work is that the quiescent galaxy fraction drops rapidly as we approach the redshift, , independent of the environment in which they reside. This result indicates that the majority of galaxies are actively forming stars at . Combined with the fact that the quiescent fraction within clusters and in the field are similar at this epoch, this indicates that studies of clusters or proto-clusters at would be crucial in revealing the very initial properties of the forming galaxy clusters as well as the evolution of galaxies in these massive structures. Also, we can expect that any cluster finding methods using the presence of old and quiescent galaxies within clusters (for example, like red-sequence technique) would miss many clusters or proto-clusters at . In our work, both of the identification of the galaxy clusters at , as well as the reliable estimation of the SFR of galaxies through the SED-fitting are possible thanks to the deep NIR data from the UKIDSS. This implies that reliable NIR (photometric or spectroscopic) data are essential for the future search and the study of the high-redshift galaxy clusters during this important epoch (). Therefore, we can expect a big leap in our understanding of the properties and the evolution of the high-redshift galaxy clusters or proto-clusters as well as of the galaxy evolution with the near-future NIR facilities either in space — like JWST (James Webb Space Telescope) — or on the ground — like GMT (Giant Magellan Telescope) with the GMACS (with the NIR-arm) and later with the NIRMOS.

This work was supported by the National Research Foundation of Korea (NRF) grant, No. 2008-0060544, funded by the Korea government (MSIP). The UKIDSS project is defined in Lawrence et al. (2007). The UKIDSS uses the UKIRT Wide Field Camera (WFCAM; Casali et al., 2007). The photometric system is described in Hewett et al. (2006), and the calibration is described in Hodgkin et al. (2009). The pipeline processing and science archive are described in Irwin et al (2009, in prep) and Hambly et al. (2008). Part of this work is based on observations taken by the CANDELS Multi-Cycle Treasury Program with the NASA/ESA HST, which is operated by the Association of Universities for Research in Astronomy, Inc., under NASA contract NAS5-26555. Facilities: UKIRT (WFCAM), Subaru (SUPRIMECAM), HST (WFC3).

## Appendix A Effects of different sSFR cut

Since galaxies selected with our redshift-dependent cut of are not quite “quiescent” (i.e., no SF) at high redshift, we examine here how a constant sSFR cut changes our main conclusion regarding SF quenching. Here, we choose a cut of sSFR = yr, which corresponds to the cut for local quiescent galaxies (e.g. Gallazzi et al., 2009; Ko et al., 2014).

Applying this constant sSFR cut, we first analyse the evolution of quiescent galaxies fraction. In Figure A1 — which corresponds to upper panels of Figure 18, we show the evolution of quiescent galaxy fraction in different mass bins. Not surprisingly, the quiescent galaxy fraction is overall lower than when we apply the 1/[]-cut. The most significant difference is that the quiescent galaxy fraction shows steeper increase from to , for galaxies with (middle and right panels), especially for , than in Figure 18. This indicates that SF quenching (to local quiescent galaxy level) is faster for more massive galaxies. This fast increase of quiescent fraction at high- agrees with Domínguez Sánchez et al. (2011) result.

While there are interesting differences in quiescent fraction evolution with different sSFR cuts, one important trend seems to hold: the quiescent galaxy fraction shows clearer dependence on their stellar mass than the environment. The environmental quenching efficiency (Figure A2) is not mass-dependent, but rises toward lower redshift (at ), similarly with Figure 20.

## References

• Baldry et al. (2006) Baldry, I. K. et al. 2006, MNRAS, 373, 469
• Balestra et al. (2010) Balestra, I. et al. 2010, A&A, 512, A12
• Balogh et al. (2000) Balogh, M. L., Navarro, J. F., & Morris, S. L. 2000, ApJ, 540, 113
• Balogh et al. (2004) Balogh, M. L. et al. 2004, MNRAS, 348, 1355
• Behroozi et al. (2013) Behroozi, P. S., Wechsler, R. H., & Conroy, C. 2013, ApJ, 770, 57
• Bell et al. (2004) Bell, E. F., Wolf, C., Meisenheimer, K., Rix, H.-W., Borch, A, Dye, S., Kleinheinrich, M., Wisotzki, L., & McIntosh, D. H. 2004, ApJ, 608, 752
• Bertin & Arnouts (1996) Bertin, E. & Arnouts, S. 1996, A&AS, 117, 393
• Birnboim & Dekel (2003) Birnboim, Y. & Dekel, A. 2003, MNRAS, 345, 349
• Blanton et al. (2005) Blanton, M. R., Eisenstein, D., Hogg, D. W., Schlegel, D. J., & Brinkmann, J. 2005, ApJ, 629, 143
• Blanton et al. (2003) Blanton, M. R. et al. 2003, ApJ, 594, 186
• Brammer et al. (2008) Brammer, G. B., van Dokkum, P. G., & Coppi, P. 2008, ApJ, 686, 1503
• Brammer et al. (2012) Brammer, G. B., van Dokkum, P. G., Franx, M., et al. 2012, ApJS, 200, 13
• Bruzual & Charlot (2003) Bruzual,G & Charlot, S. 2003, MNRAS, 344, 1000 (BC03)
• Burgarella et al. (2013) Burgarella, D. et al. 2013, A&A, 554, A70
• Butcher & Oemler (1984) Butcher, H. & Oemler, A. Jr. 1984, ApJ, 285, 426
• Calzetti et al. (2000) Calzetti, D. et al. 2000, ApJ, 533, 682
• Casali et al. (2007) Casali, M. et al. 2007, A&A, 467, 777
• Chabrier (2003) Chabrier, G. 2003, PASP, 115, 763
• Cimatti et al. (2008) Cimatti, A. et al. 2008, A&A, 482, 21
• Conselice et al. (2008) Conselice, C. J., Bundy, K., U, V., Eisenhardt, P., Lotz, J., & Newman J., 2008, MNRAS, 383, 1366
• Cooper et al. (2006) Cooper, M. C., et al. 2006, MNRAS, 370, 198
• Cooper et al. (2008) Cooper, M. C., et al. 2008, MNRAS, 383, 1058
• Cucciati et al. (2006) Cucciati, O., et al. 2006, A&A, 458, 39
• Cucciati et al. (2012) Cucciati, O., et al. 2012, A&A, 539, 31
• Daddi et al. (2004) Daddi, E. et al. 2004, ApJ, 617, 746
• Daddi et al. (2007) Daddi, E. et al. 2007, ApJ, 670, 156
• Damen et al. (2009) Damen, M. et al. 2009, ApJ, 690, 937
• Domínguez Sánchez et al. (2011) Domínguez Sánchez, H., Pozzi, F., Gruppioni, C. et al. 2011, MNRAS, 417, 900
• Dressler (1980) Dressler, A. 1980, ApJ, 236, 251
• Dressler et al. (1997) Dressler, A. et al. 1997, ApJ, 490, 577
• Elbaz et al. (2007) Elbaz, D., et al. 2007, A&A, 468, 33
• Faber et al. (2007) Faber, S. M. et al. 2007, ApJ, 665, 265
• Fassbender et al. (2011) Fassbender, R. et al. 2011, A&A, 527, L10
• Finoguenov et al. (2010) Finoguenov, A. et al. 2010, MNRAS, 403, 2063
• Furusawa et al. (2008) Furusawa, H., et al. 2008, ApJS, 176, 1
• Galametz et al. (2013) Galametz, A. et al. 2013, ApJS, 206, 10
• Gallazzi et al. (2009) Gallazzi, A. et al. 2009, ApJ, 690, 1883
• Gilbank et al. (2008) Gilbank, D. G., Yee, H. K. C., Ellingson, E., Gladders, M. D., Loh, Y.-S., Barrientos, L. F., & Barkhouse, W. A. 2008, ApJ, 673, 742
• Gobat et al. (2011) Gobat, R. et al. 2011, A&A, 526, 133
• Grogin et al. (2011) Grogin, N. et al. 2011, ApJS, 197, 35
• Grützbauch et al. (2011) Grützbauch, R., Conselice, C. J., Varela, J., Bundy, K., Cooper, M. C., Skibba, R., & Willmer, C. N. A. 2011, MNRAS, 411, 929
• Gunn & Gott (1972) Gunn, J. E. & Gott, J. R. III, 1972, ApJ, 176, 1
• Haines et al. (2008) Haines, C. P., Gargiulo, A. & Merluzzi, P. 2008, MNRAS, 385, 1201
• Hambly et al. (2008) Hambly, N. C. et al. 2008, MNRAS, 384, 637
• Hanami et al. (2012) Hanami, H. et al. 2012, PASJ, 64, 70
• Harker et al. (2006) Harker, J. J., Schiavon, R. P., Weiner, B. J., & Faber, S. M. 2006, ApJ, 647, L103
• Hartley et al. (2014) Hartley, W. G., Conselice, C. J., Mortlock, A. et al. 2014, arXiv:1406.6058
• Hartley et al. (2013) Hartley, W. G. et al. 2013, MNRAS, 431, 3045
• Hewett et al. (2006) Hewett, P. C., Warren, J. J., Leggett, S. K. & Hodgkins, S. T. 2006, MNRAS, 367, 454
• Hilton et al. (2010) Hilton, M. et al. 2010, ApJ, 718, 133
• Hodgkin et al. (2009) Hodgkin, S. T., Irwin, M. J., Hewett, P. C. & Warren, J. J. 2009, MNRAS, 394, 675
• Holden et al. (2007) Holden, B. P. et al. 2007, ApJ, 670, 190
• Hopkins & Beacom (2006) Hopkins, A. M. & Beacom, J. F. 2006, ApJ, 651, 142
• Hopkins et al. (2006) Hopkins, P. F., Hernquist, L., Cox, T. J., Di Matteo, T., Robertson, B., & Springel, V. 2006, ApJS, 163, 1
• Huang et al. (2013) Huang, J.-S. et al. 2013, ApJ, 766, 21
• Ilbert et al. (2009) Ilbert, O., Capak, P., Salvato, M. et al. 2009, ApJ, 690, 1236
• Ilbert et al. (2013) Ilbert, O. et al. 2013, arXiv:1301.3157
• Im et al. (2002) Im, M. et al. 2002, ApJ, 571, 136
• Jimenez et al. (2005) Jimenez, R., Panter, B., Heavens, A. F., & Verde, L. 2005, MNRAS, 356, 495
• Kang & Im (2009) Kang, E. & Im, M. 2009, ApJ, 691, L33
• Kang & Im (2015) Kang, E. & Im, M. 2015, JKAS, 48, 21
• Kauffmann et al. (1993) Kauffmann, G., White, S. D. M., & Guiderdoni, B. 1993, MNRAS, 264, 201
• Kauffmann et al. (2003) Kauffmann, G. et al. 2003, MNRAS, 341, 54
• Kauffmann et al. (2004) Kauffmann, G., et al. 2004, MNRAS, 353, 713
• Ko et al. (2014) Ko, J. et al. 2014, ApJ, 791, 134
• Kodama & Arimoto (1997) Kodama, T. & Arimoto, N. 1997, A&A, 320, 41
• Koekemoer et al. (2011) Koekemoer, A. M. et al. 2011, ApJS, 197, 36
• Koo et al. (2005) Koo, D. C. et al. 2005, ApJS, 157, 175
• Larson (1974) Larson, R. B. 1974, MNRAS, 166, 385
• Larson et al. (1980) Larson, R. B., Tinsley, B. M., & Caldwell, C. N. 1980, ApJ, 237, 692
• Lawrence et al. (2007) Lawrence, A. et al. 2007, MNRAS, 379, 1599
• Le Fèvre et al. (2004) Le Fèvre, O. et al. 2004, A&A, 428, 1043
• Lee et al. (2014) Lee, S.-K., Ferguson, H. C., Somerville, R. S., Giavalisco, M., Wiklind, T., & Dahlen, T. 2014, ApJ, 783, 81
• Lee et al. (2010) Lee, S.-K., Ferguson, H. C., Somerville, R. S., Wiklind, T., & Giavalisco, M. 2010, ApJ, 725, 1644
• Lee et al. (2009) Lee, S.-K., Idzi, R., Ferguson, H. C., Somerville, R. S., Wiklind, T., & Giavalisco, M. 2009, ApJS, 184, 100
• Lewis et al. (2002) Lewis, I. et al. 2002, MNRAS, 334, 673
• Lilly et al. (2007) Lilly, S. J. et al. 2007, ApJS, 172, 70
• Lotz et al. (2008) Lotz, J. M. et al. 2008, ApJ, 672,177
• Madau & Dickinson (2014) Madau, P. & Dickinson, M. 2014, ARA&A, 52, 415
• Magnelli et al. (2013) Magnelli, B. et al. 2013, A&A, 553, A132
• Maraston et al. (2010) Maraston, C., Pforr, J., Renzini, A., Daddi, E., Dickinson, M., Cimatti, A., & Tonini, C. MNRAS, 2010, 407, 830
• McCracken et al. (2012) McCracken, H. J., Milvang-Jensen, B., Dunlop, J., et al. 2012, A&A, 544,1 56
• Merson et al. (2013) Merson, A. I. et al. 2013, MNRAS, 429, 556
• Messias et al. (2012) Messias, H., Afonso, J., Salvato, M., Mobasher, B., & Hopkins, A. M. 2012, ApJ, 754, 120
• Mignoli et al. (2005) Mignoli, M. et al. 2005, A&A, 437, 883
• Moore et al. (1998) Moore, B., Lake, G., & Katz, N. 1998, ApJ, 495, 139
• Moustakas et al. (2013) Moustakas, J. et al. 2013, ApJ, 767, 50
• Muzzin et al. (2012) Muzzin, A. et al. 2012, ApJ, 746, 188
• Nantais et al. (2013) Nantais, J. B. et al. A&A, 556, 112
• Newman et al. (2013) Newman, J. A., Cooper, M. C., Davis, M., et al. 2013, ApJS, 208, 5
• Oke (1974) Oke, J. B. 1974, ApJS, 27, 21
• Pacifici et al. (2013) Pacifici, C., Kassin, S. A., Weiner, B., Charlot, S., & Gardner, J. P. 2013, ApJ, 762, L15
• Papovich et al. (2011) Papovich, C., Finkelstein, S. L., Ferguson, H. C., Lotz, J. M., & Giavalisco, M. 2011, MNRAS, 412, 1123
• Papovich et al. (2010) Papovich, C. et al. 2010, ApJ, 716, 1503
• Patel et al. (2009) Patel, S. G., Holden, B. P., Kelson, D. D., Illingworth, G. D., & Franx, M. 2009, ApJ, 705, 67
• Peng et al. (2010) Peng, Y.-J., Lilly, S. J., Kovač, K., et al. 2010, ApJ, 721, 193
• Pforr et al. (2012) Pforr, J., Maraston, C., & Tonini, C. 2012, MNRAS, 422, 3285
• Poggianti et al. (2001) Poggianti, B. M. et al. 2001, ApJ, 563, 118
• Popesso et al. (2009) Popesso, P., Dickinson, M., Nonino, M., et al. 2009, A&A, 494, 443
• Pozzetti et al. (2003) Pozzetti, L. et al. 2003, A&A, 402, 837
• Quadri et al. (2012) Quadri, R. F., Williams, R. J., Franx, M., & Hildebrandt, H. 2012, ApJ, 744, 88
• Santos et al. (2014) Santos, J. S., Altieri, B., Tanaka, M. et al. 2014, MNRAS, 438, 2565
• Schiavon et al. (2006) Schiavon, R. P. et al. 2006, ApJ, 651, L93
• Scoville et al. (2013) Scoville, N. et al. 2013, ApJS, 206, 3
• Silverman et al. (2010) Silverman, J. D. et al. 2010, ApJS, 191, 124
• Smail et al. (2008) Smail, I., Sharp, R., Swinbank, A. M., Akiyama, M., Ueda, Y., Foucaud, S., Almaini, O. & Croom, S. 2008, MNRAS, 389, 407
• Smith et al. (2005) Smith, G. P., Treu, T., Ellis, R. S., Moran, S. M., & Dressler, A. 2005, ApJ, 620, 78
• Sobral et al. (2011) Sobral, D., Best, P. N., Smail, I., Geach, J. E., Cirasuolo, M., Garn, T., & Dalton, G. B. 2011, MNRAS, 411, 675
• Sobral et al. (2014) Sobral, D. et al. 2014, MNRAS, 437, 3516
• Somerville & Primack (1999) Somerville, R. S. & Primack, J. R. 1999, MNRAS, 310, 1087
• Somerville et al. (2012) Somerville, R. S., Gilmore, R. C., Primack, J. R., & Domínguez, A. 2012, MNRAS, 423, 1992
• Somerville et al. (2008) Somerville, R. S., Hopkins, P. F., Cox, T. J., Robertson, B. E., & Hernquist, L. 2008, MNRAS, 391, 481
• Strateva et al. (2001) Strateva, I. et al. 2001, AJ, 122, 1861
• Strazzullo et al. (2013) Strazzullo, V. et al. 2013, ApJ, 772, 118
• Szokoly et al. (2004) Szokoly, G. P., Bergeron, J., Hasinger, G., et al. 2004, ApJS, 155, 271
• Tal et al. (2014) Tal, T. et al. 2014, ApJ, 789, 164
• Tran et al. (2010) Tran, K.-V. et al. 2010, ApJ, 719, L126
• Ueda et al. (2008) Ueda, Y., et al. 2008, ApJS, 179, 124
• van Breukelen et al. (2006) van Breukelen, C. et al. 2006, MNRAS, 373, L26
• van Breukelen et al. (2007) van Breukelen, C. et al. 2007, MNRAS, 382, 971
• van den Bosch et al. (2008) van den Bosch, F. C., Aquino, D., Yang, X. et al. 2008, MNRAS, 387, 79
• van der Wel et al. (2007) van der Wel, A. et al. 2007, ApJ, 670, 206
• Vanzella et al. (2008) Vanzella, E., Cristiani, S., Dickinson, M., et al. 2008, A&A, 478, 83
• Wiklind et al. (2008) Wiklind, T., Dickinson, M., Ferguson, H. C., Giavalisco, M., Mobasher, B., Grogin, N. A., & Panagia, N. 2008 ApJ, 676,781
• Wolf et al. (2003) Wolf, C., Meisenheimer, K., Rix, H.-W., Borch, A., Dye, S., & Kleinheinrich, M. 2003, A&A, 401, 73
• Yang et al. (2014) Yang, G. et al. 2014, ApJS, 215, 27
• York et al. (2000) York, D. G., Adelman, J., Anderson, J. E., et al. 2000, AJ, 120, 1579
• Zeimann et al. (2012) Zeimann, G. R. et al. 2012, ApJ, 756, 115
• Zirm et al. (2008) Zirm, A. W. et al. 2008, ApJ, 680, 224
You are adding the first comment!
How to quickly get a good reply:
• Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
• Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
• Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters