Possible Identification Of Massive and Evolved Galaxies At z\gtrsim 5

Possible Identification Of Massive and Evolved Galaxies At z≳5

Abstract

We report on the identification of the old stellar population galaxy candidates at . We developed a new infrared color selection scheme to isolate galaxies with the strong Balmer breaks at , and applied it to the ultra-deep and wide infrared survey data from the Extended Deep Survey (SEDS) and the UKIRT Infrared Deep Sky Survey. The eight objects satisfying and are selected in the 0.34 deg SEDS Ultra Deep Survey field. Rich multi-wavelength imaging data from optical to far-infrared are also used to reject blending sources and strong nebular line emitters, and we finally obtained the three most likely evolved galaxies at . Their stacked spectral energy distribution is fitted well with the old stellar population template with M M, star formation rate  M yr, dust , and age  Gyr at , where the dusty star-forming galaxies at is disfavored because of the faintness in the m. The stellar mass density of these evolved galaxy candidates,  M Mpc, is much lower than that of star-forming galaxies, but the non-zero fraction suggests that initial star-formation and quenching have been completed by .

reception date \Acceptedacception date \Publishedpublication date

\KeyWords

galaxies: evolution — galaxies: formation — galaxies: high-redshift

1 Introduction

When and how did massive quiescent or evolved galaxies like present-day elliptical galaxies form? According to the current cold dark matter paradigm, dark matter halos grow by assembling material from their neighborhood, or by merging with other halos. This view is supported by cosmological simulations (White & Rees, 1978; Davis et al., 1985; Springel et al., 2005; Mo, van den Bosch, & White, 2010). On the other hand, the growth of galaxies, the baryonic components of the halos, seems not to proceed in a simple bottom-up manner. Observations have shown that many of the most massive galaxies (M 10M) in the local universe assembled quickly and at early epochs () while less massive galaxies grew more gradually (Kodama et al., 1998; Marchesini et al., 2009, 2010; Ilbert et al., 2010; Brammer et al., 2011; Muzzin et al., 2013). It is clear that the baryonic physics of galaxy formation involves more than just cooling and accretion of gas at a rate dictated by gravity.

In order to understand the balance among several processes in galaxy formation [major/minor mergers, harassment, feedback by active galactic nuclei (AGNs) or supernovae, cold gas accretion, and so on], it is important to unravel the relationship among different galaxy populations such as actively star forming galaxies, quiescent galaxies, and AGNs. The first quiescent galaxies to emerge are especially interesting because their redshifts and masses will have much to say about the history of galaxy formation and star formation activity in the very early universe. However, our present knowledge of galaxies at these epochs are based primarily on star-forming galaxies and AGNs, due to the limited sensitivity and volume of surveys carried out to-date, as well as an intrinsic rarity of quiescent galaxies at high redshift (Kajisawa et al., 2011; Domínguez Sánchez et al., 2011; Muzzin et al., 2013). A new, larger sample of quiescent galaxies at higher redshift is needed to obtain a more complete picture of galaxy formation and evolution at early epochs.

In contrast with star-forming galaxies, which typically exhibit numerous emission lines or strong Lyman breaks, quiescent galaxies are relatively difficult to be identified at high redshift because of their nearly featureless spectral energy distributions (SEDs). The Balmer/4000-Å-break is the only strong feature in the SEDs of quiescent galaxies, and it has long been used to identify them at (Elston et al., 1988; Dickinson, 1995; Yamada et al., 1997, 2005; Franx et al., 2003; Hatch et al., 2011). For example, objects with or in the Vega magnitude system select elliptical galaxies with a strong 4000Å-break at (Extremely Red Objects: EROs; Thompson et al. (1999); MacCarthy et al. (2001); Cimatti et al. (2002); Caputi et al. (2004)). This ERO selection was extended to higher redshift, , by adopting the red color criterion (, IRAC-selected extremely red object (IERO); Yan et al. (2004)]. The near-infrared (NIR) color criterion (Vega) is often used to select quiescent galaxies at (Distant Red Galaxies: DRGs; van Dokkum et al. (2003); Kajisawa et al. (2006b); Kubo et al. (2013)). While dusty star-forming galaxies also satisfy these simple and color criteria, a substantial fraction of EROs and DRGs are spectroscopically confirmed to be passively evolved galaxies (Cimatti et al., 2002, 2004; Doherty et al., 2005; Kriek et al., 2006a, b). As spectroscopy is time-consuming, many authors have instead separated passive and dusty star-forming galaxies photometrically using two colors so that one color brackets the Balmer break and the other color represents the shape of spectrum at longer wavelengths. The color selection [ and in the Vega magnitude system] is similar to the DRG selection, but more sensitive to the evolved galaxy population at (Kajisawa et al., 2006a; Doherty et al., 2010). It is also known that passive and dust-obscured galaxies at separate well in the vs. two color diagram (Labbé et al., 2005; Papovich et al., 2006; Kubo et al., 2013). Recently, rest-frame and (rest-) color selection is favored, which allows us to distinguish passive galaxies from star-forming galaxies independently of their redshifts as long as SED fitting works well (Wuyts et al., 2007; Williams et al., 2009; Whitaker et al., 2011; Brammer et al., 2011; Morishita et al., 2014).

At , the Balmer/4000Å break is shifted to m, and detecting quiescent galaxies at those redshifts is a challenging task. However, some candidates at have been found in extremely deep NIR and IR surveys (Rodighiero et al., 2007; Wiklind et al., 2008; Mancini et al., 2009; Huang et al., 2011; Caputi et al., 2012; Nayyeri et al., 2014). Huang et al. (2011) analyzed very deep Hubble Space Telescope () and imaging in the Great Observatories Origins Deep Survey South (GOODS-South) field and identified four galaxies with . Huang et al. (2011) argued that the four extremely red objects are likely distant passive galaxies rather than nearby dusty galaxies, because their Multiband Imaging Photometer (MIPS) 24-m fluxes are weaker than expected if their rest-infrared SEDs are similar to those of the central region of M82 (a local dusty galaxy with ) or those of the dust obscured galaxies (DOGs) at . Employing a similar strategy, Caputi et al. (2012) identified 25 objects with over 180 arcmin of the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS; Grogin et al. (2011); Koekemoer et al. (2011)) Ultra Deep Survey (UDS) field. Their sample includes one massive galaxy candidate whose SED is fitted well with a 0.5-Gyr-old stellar population template with a stellar mass of  M. Many other red objects also prefer high-redshift photometric solutions (), but their red SEDs are consistent with significant dust extinction: (Caputi et al., 2012). Wiklind et al. (2008) proposed a color selection that is more sensitive to galaxies with strong Balmer breaks (Balmer Break Galaxies: BBGs) at , and which is based on the three colors, , , and . They adopted their selection to the multi-band imaging available in the GOODS-South field, and identified the 11 candidates of passive galaxies with ages  Gyr at . They mentioned that the number of candidates decrease to 4 if they allow only objects undetected in the MIPS 24m.

In this work, we applied a similar color selection of high-redshift evolved galaxies to the 3.6m and 4.5m imaging from the Extended Deep Survey (SEDS: Ashby et al. (2013)) and the -band imaging from the United Kingdom Infrared Telescope (UKIRT) Infrared Deep Sky Survey (UKIDSS: Lawrence et al. (2007)). The color is expected to be more sensitive to the Balmer break at than or because it spans a relatively short wavelength interval. The wide area of the UKIDSS UDS (Almaini et al. (in prep)) enables us to search for rare galaxies over the entire SEDS UDS field (0.34 deg), whereas the survey area of the previous studies using are relatively small. The aim of this work is to identify evolved or old stellar population galaxies at .

This paper is structured as follows. In Section 2 we summarize data used in this work. We describe the versus two-color selection scheme and further selection using multi-wavelength imaging in Section 3. Stacking analysis and results for the most likely evolved galaxy candidates are described in Section 4. In Section 5, we discuss the implications about the high- evolved galaxies. Our conclusions are summarized in Section 6. Unless specified otherwise, we use the AB magnitude system (Oke & Gunn, 1983) and adopt a cosmology with km s Mpc, , and (Komatsu et al., 2011).

2 Data

This research is primarily based on data from SEDS (Ashby et al., 2013), a very deep infrared survey using ’s Infrared Array Camera (IRAC: Fazio et al. (2004)) to cover five well-known extragalactic science fields. SEDS reaches a depth of 25.4 mag (0.25 Jy; 5) in both of the warm IRAC bands at 3.6 and 4.5 m. In this work we analyze the 0.34 deg of the SEDS UDS field, where rich multi-wavelength imaging is available from a diverse array of ground- and space-based facilities. The surveys used in this study are summarized in Table 1 and Figure 1.

The visible-wavelength photometry comes from the Subaru/XMM-Newton Deep Survey (SXDS) which covers the entire SEDS UDS area, reaching 5 limits (2 diameter aperture) of , , , (Furusawa et al., 2008). For -band imaging, we used the new data obtained in 2008 - 2009 (P.I.: H. Furusawa; Furusawa et al. (2016)) with Suprime-Cam (Miyazaki et al., 2002) on the Subaru Telescope. This new -band imaging reaches the 5 limits of (2 diameter aperture), which is 0.5 mag deeper than that of the SXDS -band imaging.

The depth and coverage of the near-infrared imaging are the most important factors in this study. Deep - or -band imaging is required to detect or constrain the flux at shorter wavelengths than the Balmer break of quiescent galaxies at . We obtained the deep -, -, and -band imaging from the UKIDSS UDS (DR10), which have 0.8 deg coverage containing nearly all of the SEDS UDS coverage. The UKIDSS project is defined in Lawrence et al. (2007). Further details on the UDS can be found in Almaini et al. (in prep). 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. (in prep) and Hambly et al. (2008). In principle, wide area coverage should increase the detection rate for quiescent galaxies at high redshift because their spatial distribution is expected to be highly biased, which is the confirmed trend at in the previous studies of EROs or DRGs (Daddi et al., 2000, 2003; Roche et al., 2002; Hamana et al., 2006; Tinker et al., 2010; Uchimoto et al., 2012). The UDS reaches depths of , , and (5). Compared with the CANDELS/ and -band imaging, the UDS -band imaging yields the advantage in the wide area and small wavelength interval between and 3.6m. We estimated the above 5 limiting magnitudes of the UDS by measuring the standard deviation of the background fluctuations with 2 diameter aperture. Our estimated limiting magnitudes are 0.3 - 0.4 mag shallower than the official values reported by the WFCAM Science Archive12. For the present work, we adopt our own measurements as the conservative limiting magnitudes.

While no survey achieved the same depth as SEDS in the mid-infrared range, shallower IRAC/5.8, 8.0, and MIPS/24m imaging from the UKIDSS Ultra Deep Survey (SpUDS; P.I.: J. Dunlop) may be useful for constraining the SEDs of galaxies. SpUDS is shallow, but covers a wide area ( deg), which contains the entire SEDS UDS coverage. We obtained the SpUDS imaging data through NASA/IPAC Infrared Science Archive (IRSA), which is operated by the Jet Propulsion Laboratory, California Institute of Technology, under contract with the National Aeronautics and Space Administration. We found that the limiting magnitudes are 21.7, 21.5, and 19.2 mag (7.6, 9.1, and 76 Jy; 5) in the 5.8, 8.0, and 24m bands by measuring the standard deviations of the background fluctuations, where we adopted 2 Full Width at Half Maximum (FWHM) of the Point Spread Function (PSF) diameter apertures (4.2, 4.6, and 12 arcsec).

We also used SPIRE 250, 350, and 500m imaging from the Herschel Multi-tiered Extragalactic Survey (HerMES: Oliver et al. (2012)). In the second HerMES data release (DR2), SPIRE imaging in the L4-UDS field (Oliver et al., 2012) is available, which contains the entire SEDS UDS coverage. The HerMES data was accessed through the Herschel Database in Marseille (HeDaM)13 operated by CeSAM and hosted by the Laboratoire d’Astrophysique de Marseille. The depth of the SPIRE imaging is estimated by using the corresponding error maps (Levenson et al., 2010), finding that the limiting magnitudes are 14.7, 14.2, and 17.1 mJy () in the 250, 350, 500m bands. We used the HerMES band-merged source catalog (xID catalog: Wang et al. (2014)) for the 250, 350, and 500 m photometry.

3 Candidate selection

3.1 Multi-band photometry for SEDS-detected objects

We used the SEDS catalog (Ashby et al., 2013) for source extraction and photometry in the 3.6 and 4.5m. The SEDS catalog yields PSF-fitted total magnitudes, from which we estimated 2.4 arcsec aperture magnitudes using the aperture correction factors of Ashby et al. (2013). This small aperture size compared with the PSF (FWHM  arcsec) is recommended by the SEDS survey team. We also investigated signal to noise ratio () in the 3.6um photometry as a function of aperture size, the results of which were that decreases with increasing aperture size and the 2.4 arcsec aperture is better. In this study, we focus on relatively isolated sources among the cataloged objects, because source confusion makes it hard to compare or combine multi-band photometry. In order to select isolated sources, we performed source extraction again on the 3.6m image by using SExtractor (Bertin & Arnouts, 1996) version 2.5.0. It is suggested that SExtractor does not resolve heavy source blending at faint flux level in IRAC images (Sanders et al., 2007; Ashby et al., 2013). The software was configured carefully so that we could extract as faint isolated sources as possible. We show the SExtractor parameter settings in Table 2. The objects detected in both of the SEDS catalog and our own SExtractor catalog are selected for our analysis. The objects used for our analysis are shown in Figure 2, where all objects in the SEDS catalog are also shown for comparison. We failed to detect not only confused but also very faint sources, which does not change the result of this paper because such faint objects in 3.6m cannot be identified as red objects, which we focus in the following sections. The completeness of our source extraction in the 3.6m image relative to the SEDS catalog is shown in Figure 3. By multiplying this relative completeness and the completeness of the SEDS catalog (Ashby et al., 2013), we estimated the completeness of this work (Figure 3).

We measured the 2.4 aperture magnitudes using the photometry task “PHOT” of Image Reduction and Analysis Facility (IRAF) for the , , , , , , , , 5.8m, and 8.0m images, where each aperture was centered on the position of the object detected in the 3.6 m image. The -, -, -, -, -, -, -, and -band images were convolved to the same resolution as that of the 3.6 m and 4.5 m images (FWHM ). We corrected the , , , , , , , and magnitudes for the Galactic extinction with A, A, A, A, A, A, A, and A. These values were estimated for the position (RA, Dec: \timeform02h18m00s, \timeform-05D00’00” in J2000) in the UDS field based on the work by Schlegel et al. (1998), assuming RAE. When we measured the diameter aperture magnitudes in the 5.8 and 8.0m images, we corrected the difference of PSFs to be matched with those of the 3.6 and 4.5m images assuming the same surface brightness profile in all bands (“negative aperture correction”). For the purpose, we smoothed the 4.5m image to the same resolution as that of the 5.8 and 8.0m images (FWHM and ), respectively. The correction factors were then obtained from the flux ratio of the original and smoothed images for each individual object. The average correction factor of the 5.8 and 8.0m images are and mag, respectively.

We checked whether our objects are sufficiently bright to be detected in the SpUDS 24m and HerMES imaging. For the SpUDS 24m image, object detection and photometry with 12 diameter apertures was performed using SExtractor. For the HerMES images, the 250, 350, and 500m fluxes were obtained for the objects in the HerMES DR1 catalog (Oliver et al., 2012; Wang et al., 2014).

3.2 Color Selection of BBGs

Our sample selection of evolved galaxies with strong Balmer breaks (BBGs; Wiklind et al. (2008)) is based on the and colors. We examined the colors of various types of galaxies using the stellar synthesis models of GALAXEV (Bruzual & Charlot (2003), hereafter BC03). In Figure 4, we show the behavior of the colors of the BC03 model galaxies over a wide variety of redshift, star formation history (SFH), age, metallicity, and dust extinction, where initial mass function (IMF) and dust extinction law were fixed at a Chabrier IMF (Chabrier, 2003) and Calzetti law (Calzetti et al., 2000), respectively. It is notable that galaxies at with ages much older than their star-formation timescale () are red in but blue in while dusty galaxies undergoing active star formation are red in both colors. From these color behaviors, we determined the color criteria of the BBGs at as follows:

 K−[3.6]>1.3 K−[3.6]>2.4([3.6]−[4.5])+0.6

The diagonal boundary in the versus two-color diagram is set parallel to the reddening vector at (black arrows in Figure 4) so that this color classification is effective even if evolved galaxies have moderate amount of dust. Extremely metal-poor galaxies without dust do not satisfy this BBG criteria even if they are enough old compared with their star-formation timescale. Eight candidate BBGs are selected, which are shown in the two-color and color-magnitude diagrams of Figure 5.

The BBG candidates have large uncertainties in the color measurements because the -band imaging is shallower than the SEDS 3.6m and 4.5m imaging. In the color-magnitude diagram of Figure 5 (right), we show the peaks of color distribution for all objects and the associated 4 uncertainties as a function of 3.6m magnitudes. The 4 color upper limits for the objects with the peak colors are below and the number of objects is less than 10000 mag down to , which implies that typical foreground or background objects are less likely to satisfy our color criteria.

We also checked stellar contamination in our sample. First, red giant stars in the Milky Way (MW) were ruled out because of their bright absolute magnitudes ranging from to in the 3.6m (Wainscoat et al., 1992; Johnson, 1966), which correspond to the apparent magnitudes () brighter than our BBGs out to 100 kpc distance. Second, we investigated the infrared colors of the brown dwarfs, which might have red colors due to the cool temperature and methane absorption (Leggett et al., 2010), using the catalog of Dupuy & Liu (2012). We found that the T7.5 or later type dwarf stars have colors redder than at the absolute magnitude range of , but their colors are extremely red (). We conclude that the BBG candidates are not MW stars but galaxies.

3.3 Multi-wavelength constraint

For the eight BBG candidates, we combined multi-band photometry from to 24m to make SEDs. We also checked the Herschel band-merged catalog (Wang et al., 2014), but none of the eight BBGs has counterpart down to 8.8 mJy at 250m ().

The eight BBG candidates may still include some contaminants. Several authors suggested that the strong emission lines from young star-forming galaxies potentially affect the broad-band photometry for high redshift galaxies (Schaerer & de Barros, 2009, 2010; Ono et al., 2010; Inoue, 2011; Stark et al., 2013). Especially, strong [OIII] lines from galaxies or H lines from galaxies can make their colors red. Also, some unforeseen effects caused by completely blending foreground objects might change the colors so that the objects satisfy the BBG color criteria. It is important to identify the most likely candidates of evolved galaxies using the multi-photometry data. For the eight color selected BBGs, we required them to be undetected in the , , , , , and 24m bands (). The observed , , and bands correspond to shorter wavelengths than the Lyman break in the rest frame for galaxies. The and non-detection requirement was adopted to reject the strongest nebular line emitters that have strong UV continuum or Ly line in those bands. The 24m non-detection requirement was adopted to reject local ultra-dusty galaxies () and Type-2 AGNs. After excluding objects detected in the , , , , , and 24m bands, we obtained only three BBGs as the most likely evolved galaxy candidates. The observed properties of these three most likely BBGs are shown in Table 3.3.

{longtable}

lccccc Observed properties of the most likely BBGs Name & RA & Dec & & &
\endhead\endlastfootSEDS_UDS_BBG-23 & \timeform02h16m38s.18 & \timeform-05D13’52”.1 & 25.240.43 & 23.750.13 & 23.730.13
SEDS_UDS_BBG-33 & \timeform02h17m20s.23 & \timeform-05D09’07”.5 & 25.450.52 & 23.980.14 & 23.660.13
SEDS_UDS_BBG-34 & \timeform02h17m15s.66 & \timeform-04D57’59”.8 & 25.5 () & 24.030.15 & 23.790.13

4 Stacking analysis

4.1 Image stacking

We made the median stacked images in the , , , , , , , , 3.6m, 4.5m, 5.8m, 8.0m, and m for the three most likely BBGs to reduce the background noise fluctuations and constrain their average SEDs. The stacking of only three objects may cause accidental superposition of very faint foreground objects in the photometric aperture. In order to check such superposition effect, we made the 800 different stacked images for each band by rotating the individual images randomly. For these 800 stacked images we performed photometry to evaluate the average magnitude and standard deviation, and finally selected one which yields the average magnitude. The standard deviations in magnitude, which are caused by image rotations, are combined with those of background noise fluctuations to estimate the photometric errors. The stacked images typically become mag deeper than the individual frames.

The stacked images are shown in Figure 6. We have no detection in the , , , , , , 5.8m, and 8.0m, but these non-detections yield strong constraint in the SED fitting below. The marginal detection in the stacked -band image is consistent with the picture that they are galaxies with the relatively faint UV continuum and weak Lyman breaks at , rather than actively star-forming galaxies with the strong UV continuum at the same redshift.

4.2 SED fitting

For the stacked images, we re-performed the photometry in the same manner as that for individual objects (subsection 3.1). After obtaining the multi-band aperture photometry, we performed SED fitting in the wavelength range between 0.4m and 8.0m to study the averaged properties of the BBGs at (redshift, SFH, age, stellar mass, and so on). We used the public SED fitting code Hyperz (Bolzonella et al., 2000). By setting and for the flux fainter than 2, we kept the number of the measurements larger than that of the fitting parameters and avoided the obvious overfitting.

Spectral templates consist of the BC03 models with exponentially declining star formation rate (SFR; = 0.01, 0.1, 1, and 10 Gyr) and the three different metallicity values (0.0001, 0.004, and 0.02), which are the same as the templates used in Figure 4. For them, redshift is changed from 0 to 7 and dust extinction is changed from 0 to 4. Model ages are forced to be less than the universe age at given redshift. We also used spectral templates with strong nebular emission lines (NL templates). In the NL templates, the BC03 model templates with constant SFR as the stellar continuum component. Our calculation of the nebular emission line flux is the same as that in Salmon et al. (2015). The number of ionizing photons was estimated from the SFH, age, and metallicity of each given model, from which we obtained H line flux assuming the case B recombination and an escape fraction of ionizing photons . The nebular emission line flux is calculated from their H line flux and the line ratio table of Inoue (2011) who used Cloudy 08.00 (Ferland et al., 1998) to derive nebular spectra as a function of metallicity. We attenuated both nebular line flux and stellar continuum with the same dust extinction, A, ranging from 0 to 2. Absorption of ionizing photons (LyC absorption) by dust in the ionized nebulae was neglected in this study. Let us take the cases where nebular emission lines contribute broad-band flux at the maximum ( and no LyC absorption) for the aim of this study. Our settings for the SED fitting are summarized in Table 3. We adopted the Calzetti dust extinction law and the Chabrier IMF for all templates.

The measured SED with the acceptable-fitting templates and the fitting probability distribution which Hyperz software calculated are shown in Figure 7. The best-fitting template is the old stellar population model at with star-formation timescale  Gyr, age  Gyr, SFR  M yr, no dust, and no nebular line. However, the dusty model at with A is also acceptable. We further checked this dusty solution using the m flux which is not used in the SED fitting. Dusty galaxies with A should have the significant dust emissions at longer wavelength and may be detectable in the 24m image. We estimated the infrared luminosity (L) of the dusty template at assuming that the luminosity attenuated by dust at 0.01m 2.2m is re-radiated at 5m 1000m (LL), which results in L L. We redshifted the model SED of Luminous Infrared Galaxy (LIRG; Rieke et al. (2009)) with the same L to predict the 24m magnitude of the dusty template, and obtained  mag. This is inconsistent with the measured flux, (), which suggests that the dusty solution is unlikely.

We show the probability distribution for only the NL templates in Figure 7 (dashed line in the right-hand panel). The strong nebular line emitter templates at and are also acceptable. While the strong nebular line emitter model cannot be rejected, we do not consider that all of the BBGs are nebular line emitters. The nebular line emitter solutions for the stacked SED have high stellar mass (M M) and high dust extinction (A). Extremely strong nebular line emitters with M M and A are very rare at (van der Wel et al., 2011; Maseda et al., 2014) and negligible especially in the photometrically selected passive galaxy sample at (Straatman et al., 2014).

From these above, we conclude that our BBG sample is expected to contain evolved galaxies at . Their physical properties are expected from the best-fitting template for the stacked SED: age  Gyr, M M, SFR  M yr and A at , where the uncertainties correspond to the 1 confidence interval () for each quantity. They are massive, but are consistent with the stellar mass of the analogs identified at intermediate redshifts (: Muzzin et al. (2013); Spitler et al. (2014)).

5 Discussion

In the following discussion about the physical properties of the evolved galaxies in the early universe, we roughly assume that all of the three BBG candidates have a stellar masses of M M and ages older than 0.3 Gyr at the redshift . This assumption is well supported by the stacking analysis, while we cannot rule out the possibility that some of them are massive, dusty, but are extremely strong nebular line emitters (subsection 4.2).

5.1 Implication about the stellar mass density of evolved galaxies at z>5

We estimated the number and stellar mass density of evolved galaxies at . The survey volume was calculated from the SEDS area (0.34 deg) and the redshift range (). The survey incompleteness was also corrected by using the completeness of our source extraction at the 3.6m magnitude range of the BBGs (),  % (see subsection 3.1 and Figure 3). The estimated number and stellar mass density of evolved galaxies at are  Mpc and  M Mpc, respectively. These densities may be underestimated, because we might have missed some evolved galaxies as a result of photometric uncertainties. We found that about 30 % of objects with an intrinsic color of are not observed as BBGs after considering the random photometric errors which are extracted from Gaussian probability distributions. This effect, however, does not change the following discussions very much.

We compared this stellar mass density of evolved galaxies at with those at of Muzzin et al. (2013) and Straatman et al. (2014). Muzzin et al. (2013) used the COSMOS/UltraVISTA survey data to select both star-forming and quiescent galaxies up to and to construct their stellar mass functions. Their classification of quiescent and star-forming galaxies using rest and colors is qualitatively similar with our BBG selection, in the sense that both selections aim to isolate galaxies with red color by Balmer breaks and blue color in longer wavelength than Balmer breaks. We estimated the stellar mass densities at by integrating their stellar mass functions down to MM, which equals the stellar mass lower limit assumed for our BBG sample. Straatman et al. (2014) used the FourStar Galaxy Evolution Survey (ZFOURGE) data to construct quiescent galaxy sample down to MM at , where their selection method was the same as that in Muzzin et al. (2013) (rest- color selection). For the 15 quiescent galaxies at , we adopted our stellar mass limit to calculate the stellar mass density. We show these stellar mass densities of quiescent galaxies with MM in Figure 8. The stellar mass density at revealed in this paper is much lower than that at and consistent with the trend that the stellar mass density of quiescent galaxies decrease continuously with increasing redshift.

We also compared the stellar mass density of evolved galaxies at with those of Wiklind et al. (2008). The color selection of Wiklind et al. (2008) for BBGs at is similar to our BBG selection in which case they require red colors, but different from our because they remove the lower redshift contaminations by using the blueward colors ( and ) instead of . They identified four candidates for evolved galaxies at without the 24m detection in the 0.04 deg area of the GOODS-South field which have the physical properties similar to ours: and  Gyr. The stellar mass density after correcting the detection incompleteness is  M Mpc, which is much more than our measurement (Figure 8). We cannot give a conclusive explanation for this inconsistency, although it is partially due to the fact that they did not use the nebular line emitter templates.

In Figure 8, the stellar mass densities of star-forming galaxies are also plotted; they were estimated by integrating the stellar mass functions of star-forming galaxies at (Muzzin et al., 2013) and Lyman Break Galaxies (LBGs) at (González et al., 2011) down to MM. From these, we estimated the quiescent fraction, which is defined as the fraction of stellar mass density of quiescent galaxies compared with the total stellar mass density. Here we assume the sum of stellar mass density of quiescent galaxies and star-forming galaxies (LBG at ) as the total stellar mass density. The quiescent fraction at is only , which is much lower than at . This suggests that evolved galaxies do not dominate the massive end of stellar mass functions at . Such rare but non-zero evolved galaxies are naturally considered to be embedded in prominent massive halos, because both theoretical works (e.g., De Lucia et al. (2006)) and observations (Tanaka et al., 2010; Kubo et al., 2013) suggested that mass assembly and quenching of galaxies preferentially proceeded in the proto-cluster regions. It is worthwhile to investigate the surrounding environments of the evolved galaxy candidates individually in future.

5.2 Possible progenitor of evolved galaxies at z>5

Finally, we discuss the progenitor of the evolved galaxies at . Our finding of evolved galaxies at strongly suggests intense starburst in these systems at much larger redshift. Submilimeter-selected galaxies (SMGs) are considered as possible progenitors, because their dust-obscured star formation is much more active than that of UV-selected star-forming galaxies which leads to rapid consumption of gas and subsequent quenching (Toft et al., 2014).

Although very few SMGs with spectroscopic redshifts are available above , the object HFLS3 (Riechers et al., 2013) may be a good example of a progenitor of the evolved galaxies at . HFLS3 was identified by ultra red colors in the Herschel SPIRE wavelength range (m), and its redshift was determined to be based on several emission lines (Riechers et al., 2013). They found that HFL3 has an extremely high SFR ( M yr), a large amount of gas ( M), and large stellar mass ( M). These quantities are not changed very much even if the lensing magnification by the two foreground galaxies is taken into account (Cooray et al., 2014). Although the gas depletion timescale of HFLS3 is 36 Myr, star formation may be suppressed on a shorter timescale ( Myr) because of supernova (and AGN) feedback. If HFLS3 continue to produce stars with SFR M yr for 10 Myr and experiences subsequent passive evolution, it should become BBG with M M and age  Gyr at , which is consistent with stellar mass and age of our evolved galaxy candidates within the uncertainties. This suggests that the most active starburst (or “maximum starburst”) SMGs can be progenitors of evolved galaxies at .

We also checked their abundances. The expected number density of progenitors depends on their duty cycle. We here define duty cycle as the ratio of the time duration for which galaxies can be detected to the time interval that corresponds to the searched redshift range. Once galaxies are quenched, they can be observed as evolved galaxies at all times, which means that the duty cycle of evolved galaxies is 100 %. If a progenitor has a duty cycle of 100 %, the expected number density should be same as the observed number density of evolved galaxies,  Mpc. But the duty cycle of a star-forming progenitor may be small, because the time duration for which it was in active star-forming phase is expected to be short. The number density of such a progenitor should be  Mpcduty-cycle. For “maximum starburst” SMGs such as HFLS3, we adopt the following simple assumptions. Intrinsic redshift distribution is flat at . Completeness of the Herschel observation of Riechers et al. (2013), which is a part of HerMES (Oliver et al., 2012) in the First Look Survey (FLS) field, Locman Hole, and GOODS-North fields, for galaxies analog to HFLS3 at is 0.6. This completeness is consistent with that estimated in Wang et al. (2014) for objects with flux of 30 mJy at 500m (flux limit of Riechers et al. (2013)). Observed time interval corresponding to is 0.42 Gyr, and star-formation timescale of HFLS3 is  Myr. These assumptions lead to a spatial number density of  Mpc and duty cycle of  % for “maximum starburst” SMGs at . The number densities of “maximum starburst” SMGs at and expected progenitor as a function of duty cycle are shown in Figure 9. The number density of “maximum starburst” SMGs is consistent with that of expected progenitors. Riechers et al. (2013) reported on the four analogs of HFLS3, which have similar red Herschel SPIRE colors. They may explain the small difference between the number density of “maximum starburst” SMGs and expected progenitors, while there are large uncertainties caused by the simple assumptions. The picture that star-formation in the progenitor of evolved galaxies at rapidly proceeds surrounded with dust is a natural extension of the picture for lower-redshift, quiescent galaxies (Toft et al., 2014; Straatman et al., 2014).

6 Summary

In this paper, we have investigated the nature of candidate evolved galaxies at . We identified eight galaxies with and down to  mag from the SEDS UDS 0.34 deg area. Among them, we further selected the three most likely BBGs at by requiring non-detection in the , , , , , and 24m bands. The best-fitting spectral template for the stacked SED is the old stellar population model with M M, SFR  M yr, and age  Gyr at . We estimated the number and stellar mass densities of evolved galaxies at down to MM assuming all of the three most likely BBGs are really such objects, which results in  Mpc and  M Mpc, respectively. This stellar mass density of evolved galaxies at is consistent with the declining trend previously observed within a range of . A very simple abundance-matching test implies that the evolved galaxies identified in this paper seem to be the descendants of “maximum starburst” SMGs, such as HFLS3 reported in Riechers et al. (2013). Our finding of evolved galaxy candidates suggests that the first starburst and subsequent quenching process have been completed by .

The BBGs are too red and too faint for us to further constrain their detailed properties against the photometric uncertainties. These galaxies will, however, be fully revealed with the next-generation infrared instruments on the future Wide-field Imaging Surveyor for High-redshift (WISH: Yamada et al. (2012)) or James Webb Space Telescope (JWST: Gardner et al. (2006)).

We thank Brett Salmon and Steven Finkelstein for the code they offered which calculates nebular emission line flux. We also appreciate Hisanori Furusawa for kindly offering the new Subaru -band imaging. This research was conducted as a part of K. M.’s phD thesis in Tohoku University. K. M. is grateful to Takashi Ichikawa, Masafumi Noguchi, and Masayuki Akiyama for valuable comments and suggestions. As this research started during K. M. stayed at the Harvard-Smithsonian Astrophysical Observatory, K. M. appreciates the SAO predoctoral program. This research is financially supported by Japan Society for the Promotion of Science through Brain Circulation Program (R2301) and Grant-in-aid for Scientific Research (13J02968 and 26400217).

Footnotes

1. affiliation: College of General Education, Osaka Sangyo University, 3-1-1, nakagaito, Daito, Osaka 574-8530, Japan
2. affiliation: Institute of Space Astronautical Science, Japan Aerospace Exploration Agency, Sagamihara, Kanagawa 252-5210, Japan
3. affiliation: Harvard-Smithsonian Center for Astrophysics, 60 Garden St., Cambridge, MA 02138, USA
4. affiliation: Harvard-Smithsonian Center for Astrophysics, 60 Garden St., Cambridge, MA 02138, USA
5. affiliation: Harvard-Smithsonian Center for Astrophysics, 60 Garden St., Cambridge, MA 02138, USA
6. footnotemark:
7. footnotemark:
8. footnotemark:
9. footnotemark:
10. footnotemark:
11. footnotemark:
12. http://surveys.roe.ac.uk/wsa/dr10plus_release.html
13. http://hedam.lam.fr
14. footnotemark:
15. footnotemark:
16. footnotemark:
17. footnotemark:
18. footnotemark:
19. footnotemark:
20. footnotemark:
21. footnotemark:
22. footnotemark:
23. footnotemark:
24. footnotemark:
25. footnotemark:
26. footnotemark:
27. footnotemark:

References

1. Almaini, O., et al., in prep.
2. Ashby, M. L. N., et al., 2013, \apj, 769, 80
3. Bertin, E., Arnouts, S., 1996, A&AS, 117, 393
4. Bolzonella, M., Miralles, J. -M., Pelló, R., 2000, A&A, 363, 476
5. Brammer, G. B., Whitaker, K. E., van Dokkum, P. G., et al. 2011, \apj, 739, 24
6. Bruzual, G. & Charlot, S., 2003, \apj, 405, 538 (BC03)
7. Calzetti, D., et al., 2000, \apj, 533, 682
8. Caputi, K. I., et al., 2012, \apj, 750, 20
9. Caputi, K. I., et al., 2004, \mnras, 353, 30
10. Casali, M., Adamson, A., Alves de Oliveira, C., et al. 2007, \aap, 467, 777
11. Chabrier, G., 2003, PASP, 115, 763
12. Cimatti, A., Daddi, E., Renzini, A., et al. 2004, \nat, 430, 184
13. Cimatti, A., et al., 2002, A&A, 381, 68
14. Cooray, A., Calanog, J., Wardlow, J. L., et al. 2014, \apj, 790, 40
15. Daddi, E., et al, \apj, 588, 50
16. Daddi, E., Cimatti, A., Pozzetti, L., Hoekstra, H., Röttgering, H. J. A., Renzini, A., Zamorani, G., & Mannucci, F., \apj, 292, 371
17. Davis, M., et al., 1985, \apj, 292, 371
18. De Lucia, G., Springel, V., White, S. D. M., Croton, D., & Kauffmann, G. 2006, \mnras, 366, 499
19. Dickinson, M., 1995, in ASP Conf. Ser. 86, Fresh View of Elliptical Galaxies, ed. A. Buzzoni, A. Renzini, & A. Serrano (San Francisco: ASP), 283
20. Doherty, M., et al., 2010, A&A, 509, 83
21. Doherty, M., et al., 2005, \mnras, 361, 525
22. Domínguez Sánchez, H., et al., 2011, \mnras, 417, 900
23. Dupuy, T. J., & Liu, M. C. 2012, \apjs, 201, 19
24. Elston, R., Rieke, G. H., & Rieke, M. J., 1988, \apj, 331, 77
25. Fazio, G. G., et al., 2004, \apjs, 154, 10
26. Ferland, G. J., Korista, K. T., Verner, D. A., Ferguson, J. B., Verner, E. M., 1998, PASP, 110, 761
27. Franx, M., et al., 2003, \apj, 587, 79
28. Furusawa, H., et al., 2008, \apjs, 176, 1
29. Furusawa, H., Kashikawa, N., Kobayashi, M. A. R., et al. 2016, arXiv:1604.02214
30. Gardner, J. P., Mather, J. C., Clampin, M., et al. 2006, \ssr, 123, 485
31. González, V., Labbé, I., Bouwens, R. J., et al. 2011, \apjl, 735, L34
32. Grogin, N. A., et al., 2011, \apjs, 197, 35
33. Hamana, T., Yamada, T., Ouchi, M., Iwata, I., & Kodama, T. 2006, \mnras, 369, 1929
34. Hambly, N. C., Collins, R. S., Cross, N. J. G., et al. 2008, \mnras, 384, 637
35. Hatch, N. A., et al., 2011, \mnras, 410, 1537
36. Hewett, P. C., Warren, S. J., Leggett, S. K., & Hodgkin, S. T. 2006, \mnras, 367, 454
37. Hodgkin, S. T., Irwin, M. J., Hewett, P. C., & Warren, S. J. 2009, \mnras, 394, 675
38. Huang, J. -S., Zheng, X. Z., Rigopoulou, D., Magdis, G., Fazio, G. G., Wang, T., 2011, \apj, 742, 13
39. Ilbert, O., et al., 2010, \apj, 709, 644
40. Inoue, A. K., 2011, \mnras, 415, 2920
41. Johnson, H. L. 1966, \araa, 4, 193
42. Kajisawa, M., Kodama, T., Tanaka, I., Yamada, T., & Bower, R., 2006a, \mnras, 371, 577
43. Kajisawa, M., et al., 2006b, PASJ, 58, 951
44. Kajisawa, M., et al., 2011, PASJ, 63, 403
45. Kodama, T., et al., 1998, A&A, 334, 99
46. Koekemoer, A. M., et al., 2011, \apjs, 197, 36
47. Komatsu, E., et al., 2011, \apjs, 192, 18
48. Kriek, M., et al., 2006a, \apj, 645, 44
49. Kriek, M., et al., 2006b, \apj, 649, 71
50. Kubo, M., et al., 2013, \apj, 778, 170
51. Labbé, I., Huang, J., Franx, M., et al. 2005, \apjl, 624, L81
52. Lawrence, A., et al., 2007, \mnras, 379, 1599
53. Leggett, S. K., Burningham, B., Saumon, D., et al. 2010, \apj, 710, 1627
54. Levenson, L., Marsden, G., Zemcov, M., et al. 2010, \mnras, 409, 83
55. Mancini, C., Matute, I., Cimatti, A., et al. 2009, \aap, 500, 705
56. Marchesini, D., et al., 2009, \apj, 701, 1765
57. Marchesini, D., et al., 2010, \apj, 725, 1277
58. Maseda, M. V., van der Wel, A., Rix, H.-W., et al. 2014, \apj, 791, 17
59. MacCarthy, P. J., et al., 2001, \apj, 560, 131
60. Miyazaki, S., Komiyama, Y., Sekiguchi, M., et al. 2002, \pasj, 54, 833
61. Mo, H., van den Bosch, F. C., & White, S., 2010, Galaxy Formation and Evolution, Cambridge: Cambridge Univ. Press
62. Morishita, T., Ichikawa, T., & Kajisawa, M. 2014, \apj, 785, 18
63. Muzzin, A., Marchesini, D., Stefanon, M., et al. 2013, \apj, 777, 18
64. Nayyeri, H., Mobasher, B., Hemmati, S., et al. 2014, \apj, 794, 68
65. Oke, J. B. & Gunn, J. E., 1983, \apj, 266, 713
66. Oliver, S. J., et al., 2012, \mnras, 424, 1614
67. Ono, Y., Ouchi, M., Shimasaku, K., Dunlop, J., Farrah, D., McLure, R., Okamura, S., 2010, \apj, 724, 1524
68. Papovich, C., Moustakas, L. A., Dickinson, M., et al. 2006, \apj, 640, 92
69. Riechers, D. A., Bradford, C. M., Clements, D. L., et al. 2013, \nat, 496, 329
70. Rieke, G. H., Alonso-Herrero, A., Weiner, B. J., et al. 2009, \apj, 692, 556
71. Roche, N. D., Almani, O., Dunlop, J., Ivison, R. J., Willott, C. J., 2002, \mnras, 337, 1282
72. Rodighiero, G., Cimatti, A., Franceschini, A., et al. 2007, \aap, 470, 21
73. Salmon, B., Papovich, C., Finkelstein, S. L., et al. 2015, \apj, 799, 183
74. Sanders, D. B., et al., 2007, \apjs, 172, 86
75. Schaerer, D. & de Barros, S., 2009, A&A, 502, 423
76. Schaerer, D. & de Barros, S., 2010, A&A, 515, 73
77. Schlegel, D. J., Finkbeiner, D. P., & Davis, M. 1998, \apj, 500, 525
78. Spitler, L. R., Straatman, C. M. S., Labbé, I., et al. 2014, \apjl, 787, L36
79. Springel, V., et al., 2005, Nature, 435, 629
80. Stark, D. P., Schenker, M. A., Ellis, R., Robertson, B., McLure, R., & Dunlop, J., 2013, \apj, 763, 129
81. Straatman, C. M. S., Labbé, I., Spitler, L. R., et al. 2014, \apjl, 783, L14
82. Tanaka, M., De Breuck, C., Venemans, B., & Kurk, J. 2010, \aap, 518, AA18
83. Thompson, D., et al., 1999, \apj, 523, 100
84. Tinker, J. L., Wechsler, R. H., Zheng, Z., 2010, \apj, 709, 67
85. Toft, S., Smolčić, V., Magnelli, B., et al. 2014, \apj, 782, 68
86. Uchimoto, Y. K., et al., 2012, \apj, 750, 116
87. van Dokkum, P. G., et al., 2003, \apj, 587, 83
88. van der Wel, A., Straughn, A. N., Rix, H.-W., et al. 2011, \apj, 742, 111
89. Wainscoat, R. J., Cohen, M., Volk, K., Walker, H. J., & Schwartz, D. E. 1992, \apjs, 83, 111
90. Wang, L., Viero, M., Clarke, C., et al. 2014, \mnras, 444, 2870
91. Whitaker, K. E., Labbé, I., van Dokkum, P. G., et al. 2011, \apj, 735, 86
92. White, S. D. M. & Rees, M. J., 1978, \mnras, 183, 341
93. Wiklind, T., Dickinson, M., Ferguson, H. C., et al. 2008, \apj, 676, 781
94. Williams, R. J., Quadri, R. F., Franx, M., van Dokkum, P., & Labbé, I. 2009, \apj, 691, 1879
95. Wuyts, S., Labbé, I., Franx, M., et al. 2007, \apj, 655, 51
96. Yamada, T., et al., 1997, \apj, 487, 125
97. Yamada, T., Kodama, T., Akiyama, M., et al. 2005, \apj, 634, 861
98. Yamada, T., Iwata, I., Ando, M., et al. 2012, \procspie, 8442,
99. Yan, H., Dickinson, M., Eisenhardt, P. R. M., et al. 2004, \apj, 616, 63
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