We have initiated a new survey for local extremely metal-poor galaxies (EMPGs) with Subaru/Hyper Suprime-Cam (HSC) large-area (~500 deg^2) optical images reaching a 5 sigma limit of ~26 magnitude, about 100 times deeper than the Sloan Digital Sky Survey (SDSS). To select Z/Z_sun<0.1 EMPGs from ~40 million sources detected in the Subaru images, we first develop a machine-learning (ML) classifier based on a deep neural network algorithm with a training data set consisting of optical photometry of galaxy, star, and QSO models. We test our ML classifier with SDSS objects having spectroscopic metallicity measurements, and confirm that our ML classifier accomplishes 86%-completeness and 46%-purity EMPG classifications with photometric data. Applying our ML classifier to the photometric data of the Subaru sources as well as faint SDSS objects with no spectroscopic data, we obtain 27 and 86 EMPG candidates from the Subaru and SDSS photometric data, respectively. We conduct optical follow-up spectroscopy for 10 out of our EMPG candidates with Magellan/LDSS-3+MagE, Keck/DEIMOS, and Subaru/FOCAS, and find that the 10 EMPG candidates are star-forming galaxies at z=0.007-0.03 with large H_beta equivalent widths of 104-265 A, stellar masses of log(M*/M_sun)=5.0-7.1, and high specific star-formation rates of ~300 Gyr^{-1}, which are similar to those of early galaxies at z>6 reported recently. We spectroscopically confirm that 3 out of 10 candidates are truly EMPGs with Z/Z_sun<0.1, one of which is HSC J1631+4426, the most metal-poor galaxy with Z/Z_sun=0.016 reported ever.

We present element abundance ratios and ionizing radiation of local young low-mass (~$10^{6}$ M_sun) extremely metal poor galaxies (EMPGs) with a 2% solar oxygen abundance (O/H)_sun and a high specific star-formation rate (sSFR~300 Gyr$^{-1}$), and other (extremely) metal poor galaxies, which are compiled from Extremely Metal-Poor Representatives Explored by the Subaru Survey (EMPRESS) and the literature. Weak emission lines such as [FeIII]4658 and HeII4686 are detected in very deep optical spectra of the EMPGs taken with 8m-class telescopes including Keck and Subaru (Kojima et al. 2019, Izotov et al. 2018), enabling us to derive element abundance ratios with photoionization models. We find that neon- and argon-to-oxygen ratios are comparable to those of known local dwarf galaxies, and that the nitrogen-to-oxygen abundance ratios (N/O) are lower than 20% (N/O)_sun consistent with the low oxygen abundance. However, the iron-to-oxygen abundance ratios (Fe/O) of the EMPGs are generally high; the EMPGs with the 2%-solar oxygen abundance show high Fe/O ratios of ~90-140% (Fe/O)_sun, which are unlikely explained by suggested scenarios of Type Ia supernova iron productions, iron's dust depletion, and metal-poor gas inflow onto previously metal-riched galaxies with solar abundances. Moreover, these EMPGs have very high HeII4686/H$\beta$ ratios of ~1/40, which are not reproduced by existing models of high-mass X-ray binaries whose progenitor stellar masses are less than 120 M_sun. Comparing stellar-nucleosynthesis and photoionization models with a comprehensive sample of EMPGs identified by this and previous EMPG studies, we propose that both the high Fe/O ratios and the high HeII4686/H$\beta$ ratios are explained by the past existence of super massive ($>$300 M_sun) stars, which may evolve into intermediate-mass black holes ($\gtrsim$100 M_sun).

We present the first stellar velocity dispersion measurement of a massive quenching galaxy at z=4.01. The galaxy is first identified as a massive z>~4 galaxy with suppressed star formation from photometric redshifts based on deep multi-band data in the UKIDSS Ultra Deep Survey field. A follow-up spectroscopic observation with MOSFIRE on Keck revealed strong multiple absorption features, which are identified as Balmer absorption lines, giving a secure redshift of z=4.01. Thanks to the high S/N of the spectrum, we are able to estimate the stellar velocity dispersion, sigma=268+/-59 km/s. This velocity dispersion is consistent with that of massive galaxies today, implying no significant evolution in stellar velocity dispersion over the last 12 Gyr. Based on an upper limit on its physical size from deep optical images (r_eff<1.3 kpc), we find that its dynamical mass is consistent with the stellar mass inferred from photometry. Furthermore, the galaxy is located on the mass fundamental plane extrapolated from lower redshift galaxies. Combining all these results, we find that the velocity dispersion does not significantly evolve with redshift, although the size and mass of massive quenched galaxies do. This suggests that the mass in the core of massive galaxies does not evolve significantly, while most of the mass growth occurs in the outskirts of the galaxies, which also increases the size. This picture is consistent with a two-phase formation scenario in which mass and size growth is due to accretion in the outskirts of galaxies via mergers.

We present the morphology and stellar population of 27 extremely metal poor galaxies (EMPGs) at $z\sim0$ with metallicities of $0.01-0.1$ Z$_{\odot}$. We conduct multi-component surface brightness (SB) profile fitting for the deep Subaru/HSC $i$-band images of the EMPGs with the Galfit software, carefully removing the SB contributions of potential associated galaxies (PAGs). We find that the EMPGs with a median stellar mass of $\log(M_{*}/{\rm M}_{\odot})=5.77$ have a median S\'ersic index of $n=1.08$ and a median effective radius of $r_{\rm e}=176$ pc, suggesting that typical EMPGs are very compact disk galaxies. We compare the EMPGs with $z\sim 6$ galaxies and local galaxies on the size-mass ($r_{\rm e}$--$M_*$) diagram, and identify that the majority of the EMPGs have a $r_{\rm e}$--$M_*$ relation similar to $z\sim0$ star-forming galaxies rather than $z\sim6$ galaxies. Not every EMPG is a local analog of high-$z$ young galaxies in the $r_{\rm e}$--$M_*$ relation. We also study the PAGs of our galaxies, and find that 23 out of the 27 EMPGs show detectable PAGs within the projected distance of 10 kpc. The PAGs have median values of $n=0.93$, $r_{\rm e}=1.41$ kpc, and $\log(M_{*}/{\rm M}_{\odot})=7.47$ that are similar to those of local dwarf irregulars and ultra-diffuse galaxies. A spectrum of one EMPG-PAG system, so far available, indicates that the PAG is dynamically related to the EMPG with a median velocity difference of $\Delta V=96.1$ km s$^{-1}$. This moderately-large $\Delta V$ cannot be explained by the dynamics of the PAG, but likely by the infall on the PAG. EMPGs may form in infalling gas, and become part of PAG disks undergoing further metal enrichment.

We construct a model of H$\alpha$ emitters (HAEs) based on a semi-analytic galaxy formation model, the New Numerical Galaxy Catalog ($\nu^2$GC). In this paper, we report our estimate for the field variance of the HAE distribution. By calculating the H$\alpha$ luminosity from the star-formation rate of galaxies, our model well reproduces the observed H$\alpha$ luminosity function (LF) at $z=0.4$. The large volume of the $\nu^2$GC makes it possible to examine the spatial distribution of HAEs over a region of (411.8 Mpc)$^3$ in the comoving scale. The surface number density of $z=0.4$ HAEs with $L_{\rm H\alpha} \geq 10^{40}$ erg s$^{-1}$ is 308.9 deg$^{-2}$. We have confirmed that the HAE is a useful tracer for the large-scale structure of the Universe because of their significant overdensity ($>$ 5$\sigma$) at clusters and the filamentary structures. The H$\alpha$ LFs within a survey area of $\sim$2 deg$^2$ (typical for previous observational studies) show a significant field variance up to $\sim$1 dex. Based on our model, one can estimate the variance on the H$\alpha$ LFs within given survey areas.

We have obtained spectroscopic confirmation with Hubble Space Telescope WFC3/G141 of a first sizeable sample of nine quiescent galaxies at 2.4<z<3.3. Their average near-UV/optical rest-frame spectrum is characterized by low attenuation (Av$\sim$0.6 mag) and a strong Balmer break, larger than the 4000 A break, corresponding to a fairly young age of $\sim$300 Myr. This formally classifies a substantial fraction of classically selected quiescent galaxies at z $\sim$ 3 as post-starbursts, marking their convergence to the quenching epoch. The rapid spectral evolution with respect to z $\sim$ 1.5 quiescent galaxies is not matched by an increase of residual star-formation, as judged from the weak detection of [O II]3727A emission, pointing to a flattening of the steep increase in gas fractions previously seen from z $\sim$ 0 to 1.8. However, radio 3GHz stacked emission implies either much stronger dust-obscured star formation or substantial further evolution in radio-mode AGN activity with respect to z $\sim$ 1.5.

We report two secure ($z=3.775, 4.012$) and one tentative ($z\approx3.767$) spectroscopic confirmations of massive and quiescent galaxies through $K$-band observations with Keck/MOSFIRE and VLT/X-Shooter. The stellar continuum emission, the absence of strong nebular emission lines and the lack of significant far-infrared detections confirm the passive nature of these objects, disfavoring the alternative solution of low-redshift dusty star-forming interlopers. We derive stellar masses of $\mathrm{log}(M_{\star}/M_\odot)\sim11$ and ongoing star formation rates placing these galaxies $\gtrsim 1-2$ dex below the main sequence at their redshifts. The adopted parametrization of the star formation history suggests that these sources experienced a strong ($\langle \rm SFR \rangle \sim 1200-3500\,M_\odot\,\mathrm{yr}^{-1}$) and short ($\sim 50$ Myr) burst of star formation, peaking $\sim 150-500$ Myr before the time of observation, all properties reminiscent of the characteristics of sub-millimeter galaxies (SMGs) at $z>4$. We investigate this connection by comparing the comoving number densities and the properties of these two populations. We find a fair agreement only with the deepest sub-mm surveys detecting not only the most extreme starbursts, but also more normal galaxies. We support these findings by further exploring the Illustris-TNG cosmological simulation, retrieving populations of both fully quenched massive galaxies at $z\sim3-4$ and SMGs at $z\sim4-5$, with number densities and properties in agreement with the observations at $z\sim3$, but in increasing tension at higher redshift. Nevertheless, as suggested by the observations, not all the progenitors of quiescent galaxies at these redshifts shine as bright SMGs in their past and, similarly, not all bright SMGs quench by $z\sim3$, both fractions depending on the threshold assumed to define the SMGs themselves.

We have obtained data for 41 star forming galaxies in the MUSE Atlas of Disks (MAD) survey with VLT/MUSE. These data allow us, at high resolution of a few 100 pc, to extract ionized gas kinematics ($V, \sigma$) of the centers of nearby star forming galaxies spanning 3 dex in stellar mass. This paper outlines the methodology for measuring the ionized gas kinematics, which we will use in subsequent papers of this survey. We also show how the maps can be used to study the kinematics of diffuse ionized gas for galaxies of various inclinations and masses. Using two different methods to identify the diffuse ionized gas, we measure rotation velocities of this gas for a subsample of 6 galaxies. We find that the diffuse ionized gas rotates on average slower than the star forming gas with lags of 0-10 km/s while also having higher velocity dispersion. The magnitude of these lags is on average 5 km/s lower than observed velocity lags between ionized and molecular gas. Using Jeans models to interpret the lags in rotation velocity and the increase in velocity dispersion we show that most of the diffuse ionized gas kinematics are consistent with its emission originating from a somewhat thicker layer than the star forming gas, with a scale height that is lower than that of the stellar disk.

This paper presents the second data release of the Hyper Suprime-Cam Subaru Strategic Program, a wide-field optical imaging survey on the 8.2 meter Subaru Telescope. The release includes data from 174 nights of observation through January 2018. The Wide layer data cover about 300 deg^2 in all five broadband filters (grizy) to the nominal survey exposure (10min in gr and 20min in izy). Partially observed areas are also included in the release; about 1100 deg^2 is observed in at least one filter and one exposure. The median seeing in the i-band is 0.6 arcsec, demonstrating the superb image quality of the survey. The Deep (26 deg^2) and UltraDeep (4 deg^2) data are jointly processed and the UltraDeep-COSMOS field reaches an unprecedented depth of i~28 at 5 sigma for point sources. In addition to the broad-bands, narrow-band data are also available in the Deep and UltraDeep fields. This release includes a major update to the processing pipeline, including improved sky subtraction, PSF modeling, object detection, and artifact rejection. The overall data quality has been improved, but this release is not without problems; there is a persistent deblender problem as well as new issues with masks around bright stars. The user is encouraged to review the issue list before utilizing the data for scientific explorations. All the image products as well as catalog products are available for download. The catalogs are also loaded to a database, which provides an easy interface for users to retrieve data for objects of interest. In addition to these main data products, detailed galaxy shape measurements withheld from the Public Data Release 1 (PDR1) are now available to the community. The shape catalog is drawn from the S16A internal release, which has a larger area than PDR1 (160 deg^2). All products are available at the data release site, https://hsc-release.mtk.nao.ac.jp/.

We investigate the contribution of clumps and satellites to the galaxy mass assembly. We analyzed spatially-resolved Hubble Space Telescope observations (imaging and slitless spectroscopy) of 53 star-forming galaxies at z ~ 1 - 3. We created continuum and emission line maps and pinpointed residual "blobs" detected after subtracting the galaxy disk. Those were separated into compact (unresolved) and extended (resolved) components. Extended components have sizes ~ 2 kpc and comparable stellar mass and age as the galaxy disks, whereas the compact components are 1.5 dex less massive and 0.4 dex younger than the disks. Furthermore the extended blobs are typically found at larger distances from the galaxy barycenter than the compact ones. Prompted by these observations and by the comparison with simulations, we suggest that compact blobs are in-situ formed clumps, whereas the extended ones are accreting satellites. Clumps and satellites enclose respectively ~ 20% and ~< 80% of the galaxy stellar mass, ~ 30% and ~ 20% of its star formation rate. Considering the compact blobs, we statistically estimated that massive clumps (Mstar >~ 10^9 Msun) have lifetimes of ~ 650 Myr, and the less massive ones (10^8 < Mstar < 10^9 Msun) of ~ 145 Myr. This supports simulations predicting long-lived clumps (lifetime > 100 Myr). Finally, ~< 30% (13%) of our sample galaxies are undergoing single (multiple) merger(s), they have a projected separation ~< 10 kpc, and the typical mass ratio of our satellites is 1:5 (but ranges between 1:10 and 1:1), in agreement with literature results for close pair galaxies.

We have obtained three-dimensional maps of the universe in $\sim200\times200\times80$ comoving Mpc$^3$ (cMpc$^3$) volumes each at $z=5.7$ and $6.6$ based on a spectroscopic sample of 179 galaxies that achieves $\gtrsim80$\% completeness down to the Ly$\alpha$ luminosity of $\log(L_{\rm Ly\alpha}/[\mathrm{erg\ s^{-1}}])=43.0$, based on our Keck and Gemini observations and the literature. The maps reveal filamentary large-scale structures and two remarkable overdensities made out of at least 44 and 12 galaxies at $z=5.692$ (z57OD) and $z=6.585$ (z66OD), respectively, making z66OD the most distant overdensity spectroscopically confirmed to date with $>10$ spectroscopically confirmed galaxies. We compare spatial distributions of submillimeter galaxies at $z\simeq 4-6$ with our $z=5.7$ galaxies forming the large-scale structures, and detect a $99.97\%$ signal of cross correlation, indicative of a clear coincidence of dusty star-forming galaxy and dust unobscured galaxy formation at this early epoch. The galaxies in z57OD and z66OD are actively forming stars with star formation rates (SFRs) $\gtrsim5$ times higher than the main sequence, and particularly the SFR density in z57OD is 10 times higher than the cosmic average at the redshift (a.k.a. the Madau-Lilly plot). Comparisons with numerical simulations suggest that z57OD and z66OD are protoclusters that are progenitors of the present-day clusters with halo masses of $\sim10^{14}\ \mathrm{M_\odot}$.

We present the results of our near-infrared (NIR) spectroscopic observations of a rich cluster candidate around a radio galaxy at $z=1.52$ (4C65.22) with Subaru/MOIRCS and LBT/LUCI. We observed 71 galaxies mostly on the star-forming main sequence selected by our previous broad-band (photo-$z$) and narrow-band H$\alpha$ imaging observation in this cluster environment. We successfully confirmed the redshifts of 39 galaxies, and conclude that this is a gravitationally bound, real cluster at $z=1.517$. Our spectroscopic data also suggest a hint of large-scale filaments or sheet-like three-dimensional structures crossing at the highest-density cluster core. By stacking the spectra to derive their average interstellar medium (ISM) gas-phase metallicity based on the [NII]/H$\alpha$ emission line flux ratio, we find that the mass-metallicity relation (MZR) in the 4C65.22 cluster environment is consistent with that of H$\alpha$-selected field galaxies at similar redshifts. Our results suggest that the environmental impacts on the MZR is small at high redshifts, but a larger sample of high-$z$ clusters and their member galaxies is still required to fully address the effect of environment as well as its cluster-cluster variation.

We investigate the multi-wavelength properties of host galaxies of 3701 X-ray-selected active galactic nuclei (AGNs) out to z~5 in the Chandra-COSMOS Legacy Survey. Thanks to the extensive multi-wavelength photometry available in the COSMOS field, we derive AGN luminosities, host stellar masses, and star formation rates (SFRs) via a multi-component SED fitting technique. Type 1 and Type 2 AGNs follow the same intrinsic Lx-L6um relation, suggesting that mid-infrared emission is a reasonably good measure of the AGN accretion power regardless of obscuration. We find that there is a strong increase in Type 1 AGN fraction toward higher AGN luminosity, possibly due to the fact that Type 1 AGNs tend to be hosted by more massive galaxies. The AGN luminosity and SFR are consistent with an increase toward high stellar mass, while both the Mstellar-dependence is weaker towards the high-mass end, which could be interpreted as a consequence of quenching both star formation and AGN activity in massive galaxies. AGN host galaxies tend to have SFRs that are consistent with normal star-forming galaxies, independent of AGN luminosities. We confirm that black hole accretion rate and SFR are correlated up to z~5, when forming stars. The majority (~73%) of our AGN sample are faint in the far-infrared, implying that the moderate-luminosity AGNs seem to be still active after the star formation is suppressed. It is not certain whether AGN activity plays a role in quenching the star formation. We conclude that both AGN activity and star formation might be more fundamentally related to host stellar mass.

Based on optical/near-IR Magellan FIRE spectra of 25 starburst galaxies at 0.5 < z < 0.9, Calabr\`o et al.(2018) showed that their attenuation properties can be explained by a single-parameter sequence of total obscurations ranging from A(V)=2 to A(V)=30 towards the starburst core centers in a mixed stars and dust configuration. We investigate here the origin of this sequence for the same sample. We show that total attenuations anti-correlate with the starburst sizes in radio (3 GHz) with a significance larger than 5sigma and a scatter of 0.26 dex. More obscured and compact starbursts also show enhanced N2 (=[NII]/Halpha) ratios and larger line velocity widths that we attribute to an increasing shock contribution toward later merger phases, driven by deeper gravitational potential wells at the coalescence. Additionally, the attenuation is also linked to the equivalent width (EW) of hydrogen recombination lines, which is sensitive to the luminosity weighted age of the relatively unobscured stellar populations. Overall, the correlations among A(V), radio size, line width, N2 and EW of Balmer/Paschen lines converge towards suggesting an evolutionary sequence of merger stages: all of these quantities are likely to be good time-tracers of the merger phenomenon, and their large spanned range appears to be characteristic of the different merger phases. Half of our sample at higher obscurations have radio sizes approximately 3 times smaller than early type galaxies at the same redshift, suggesting that, in analogy with local Ultraluminous Infrared galaxies (ULIRGs), these cores cannot be directly forming elliptical galaxies. Finally, we detect mid-IR AGN torus for half of our sample and additional X-ray emission for 6 starbursts; intriguingly, the latter have systematically more compact sizes, suggestive of emerging AGNs towards later merger stages, possibly precursors of a later QSO phase.

We study the physical properties of the ionized gas in local disks using the sample of 38 nearby $\sim10^{8.5-11.2}$M$_\odot$ Star-Forming Main Sequence (SFMS) galaxies observed so far as part of the MUSE Atlas of Disks (MAD). Specifically, we use all strong emission lines in the MUSE wavelength range 4650-9300 \AA\ to investigate the resolved ionized gas properties on $\sim$100 pc scales. This spatial resolution enables us to disentangle HII regions from the Diffuse Ionized Gas (DIG) in the computation of gas metallicities and Star Formation Rates (SFRs) of star forming regions. The gas metallicities generally decrease with radius. The metallicity of the HII regions is on average $\sim$0.1 dex higher than that of the DIG, but the metallicity radial gradient in both components is similar. The mean metallicities within the inner galaxy cores correlate with the total stellar mass of the galaxies. On our <100 pc scales, we find two correlations previously reported at kpc scales: a spatially resolved Mass-Metallicity Relation (RMZR) and a spatially resolved SFMS (RSFMS). We find no secondary dependency of the RMZR with the SFR density. We find that both resolved relations have a local origin, as they do not depend on the total stellar mass. The observational results of this paper are consistent with the inside-out scenario for the growth of galactic disks.

We present an analysis of the molecular gas properties, based on CO(2 - 1) emission, of twelve starburst galaxies at z~1.6 selected by having a boost (>~4x) in their star formation rate (SFR) above the average star-forming galaxy at an equivalent stellar mass. ALMA observations are acquired of six additional galaxies than previously reported through our effort. As a result of the larger statistical sample, we significantly detect, for the first time at high-z, a systematically lower L'_CO/L_IR ratio in galaxies lying above the star-forming `main sequence' (MS). Based on an estimate of alpha_CO (i.e., the ratio of molecular gas mass to L'_CO(1-0)), we convert the observational quantities (e.g., L'_CO/L_IR) to physical units (M_gas/SFR) that represent the gas depletion time or its inverse, the star formation efficiency. We interpret the results as indicative of the star formation efficiency increasing in a continuous fashion from the MS to the starburst regime, whereas the gas fractions remain comparable to those of MS galaxies. Although, the balance between an increase in star-formation efficiency or gas fraction depends on the adopted value of alpha_CO as discussed.

We obtained optical/near-IR rest-frame Magellan FIRE spectra (including Pa$\beta$ and Pa$\gamma$) of 25 starburst galaxies at 0.5<z<0.9, with average star formation rates (SFR) x7 above the Main Sequence (MS). We find that Paschen-to-Balmer line ratios saturate around a constant value corresponding to $A_{\rm V}\sim$2-3 mag, while line to IR luminosity ratios suggest a large range of more extreme obscurations and appear to be uncorrelated to the former. This behavior is not consistent with standard attenuation laws derived for local and distant galaxies, while being remarkably consistent with observations of starburst cores in which young stars and dust are homogeneously mixed. This model implies $A_{\rm V}=$2-30 mag attenuation to the center of starburst cores, with a median of ~9 mag (a factor of 4000). X-ray hardness ratios for 6 AGNs in our sample and column densities derived from observed dust masses and radio sizes independently confirm this level of attenuation. In these conditions observed optical/near-IR emission comes from surface regions, while inner starburst cores are invisible. We thus attribute the high [NII]/H$\alpha$ ratios to widespread shocks from accretion, turbulence and dynamic disturbances rather than to AGNs. The large range of optical depths demonstrates that substantial diversity is present within the starburst population, possibly connected to different merger phases or progenitor properties. The majority of our targets are, in fact, morphologically classified as mergers. We argue that the extreme obscuration provides in itself smoking gun evidence of their merger origin, and a powerful tool for identifying mergers at even higher redshifts.

We present ALMA observations of the 870$\mu$m continuum and CO(4-3) line emission in the core of the galaxy cluster ClJ1449+0856 at z=2, a NIR-selected, X-ray detected system in the mass range of typical progenitors of today's massive clusters. The 870$\mu$m map reveals six F$_{870\mu m}$ > 0.5 mJy sources spread over an area of 0.07 arcmin$^2$, giving an overdensity of a factor ~10 (6) with respect to blank field counts down to F$_{870\mu m}$ > 1 (0.5) mJy. On the other hand, deep CO(4-3) follow-up confirms membership of three of these sources, but suggests that the remaining three, including the brightest 870$\mu$m sources in the field (F$_{870\mu m}\gtrsim$2 mJy), are likely interlopers. The measurement of 870$\mu$m continuum and CO(4-3) line fluxes at the positions of previously-known cluster members provides a deep probe of dusty star formation occurring in the core of this high-redshift structure, adding up to a total SFR~700$\pm$100 M$_{\odot}$/yr and yielding an integrated star formation rate density of ~10$^4$ M$_{\odot}$/yr/Mpc$^3$, five orders of magnitude larger than in the field at the same epoch, due to the concentration of star-forming galaxies in the small volume of the dense cluster core. The combination of these observations with previously available HST imaging highlights the presence in this same volume of a population of galaxies with already suppressed star formation. This diverse composition of galaxy populations in ClJ1449+0856 is especially highlighted at the very cluster center, where a complex assembly of quiescent and star-forming sources is likely forming the future Brightest Cluster Galaxy.

We use ALMA and JVLA observations of the galaxy cluster Cl J1449+0856 at z=1.99, in order to study how dust-obscured star-formation, ISM content and AGN activity are linked to environment and galaxy interactions during the crucial phase of high-z cluster assembly. We present detections of multiple transitions of $^{12}$CO, as well as dust continuum emission detections from 11 galaxies in the core of Cl J1449+0856. We measure the gas excitation properties, star-formation rates, gas consumption timescales and gas-to-stellar mass ratios for the galaxies. We find evidence for a large fraction of galaxies with highly-excited molecular gas, contributing $>$50% to the total SFR in the cluster core. We compare these results with expectations for field galaxies, and conclude that environmental influences have strongly enhanced the fraction of excited galaxies in this cluster. We find a dearth of molecular gas in the galaxies' gas reservoirs, implying a high star-formation efficiency (SFE) in the cluster core, and find short gas depletion timescales $\tau$<0.1-0.4 Gyrs for all galaxies. Interestingly, we do not see evidence for increased specific star-formation rates (sSFRs) in the cluster galaxies, despite their high SFEs and gas excitations. We find evidence for a large number of mergers in the cluster core, contributing a large fraction of the core's total star-formation compared with expectations in the field. We conclude that the environmental impact on the galaxy excitations is linked to the high rate of galaxy mergers, interactions and active galactic nuclei in the cluster core.

We investigate Lya, [OIII]5007, Ha, and [CII]158um emission from 1124 galaxies at z=4.9-7.0. Our sample is composed of 1092 Lya emitters (LAEs) at z=4.9, 5.7, 6.6, and 7.0 identified by Subaru/Hyper Suprime-Cam (HSC) narrowband surveys covered by Spitzer large area survey with Subaru/HSC (SPLASH) and 34 galaxies at z=5.148-7.508 with deep ALMA [CII]158um data in the literature. Fluxes of strong rest-frame optical lines of [OIII] and Ha (Hb) are constrained by significant excesses found in the SPLASH 3.6 and 4.5um photometry. At z=4.9, we find that the rest-frame Ha equivalent width and the Lya escape fraction f_Lya positively correlate with the rest-frame Lya equivalent width EW^0_Lya. The f_Lya-EW^0_Lya correlation is similarly found at z~0-2, suggesting no evolution of the correlation over z~0-5. The typical ionizing photon production efficiency of LAEs is logxi_ion/[Hz erg^-1]~25.5 significantly (60-100%) higher than those of LBGs at a given UV magnitude. At z=5.7-7.0, there exists an interesting turn-over trend that the [OIII]/Ha flux ratio increases in EW^0_Lya~0-30 A, and then decreases out to EW^0_Lya~130 A. We also identify an anti-correlation between a [CII] luminosity to star-formation rate ratio (L_[CII]/SFR) and EW^0_Lya at the >99% confidence level. We carefully investigate physical origins of the correlations with stellar-synthesis and photoionization models, and find that a simple anti-correlation between EW_Lya^0 and metallicity explains self-consistently all of the correlations of Lya, Ha, [OIII]/Ha, and [CII] identified in our study, indicating detections of metal-poor (~0.03 Zo) galaxies with EW^0_Lya~200 A.

We derive two-dimensional dust attenuation maps at $\sim1~\mathrm{kpc}$ resolution from the UV continuum for ten galaxies on the $z\sim2$ Star-Forming Main Sequence (SFMS). Comparison with IR data shows that 9 out of 10 galaxies do not require further obscuration in addition to the UV-based correction, though our sample does not include the most heavily obscured, massive galaxies. The individual rest-frame $V$-band dust attenuation (A$_{\rm V}$) radial profiles scatter around an average profile that gently decreases from $\sim1.8$ mag in the center down to $\sim0.6$ mag at $\sim3-4$ half-mass radii. We use these maps to correct UV- and H$\alpha$-based star-formation rates (SFRs), which agree with each other. At masses $<10^{11}~M_{\rm sun}$, the dust-corrected specific SFR (sSFR) profiles are on average radially constant at a mass-doubling timescale of $\sim300~\mathrm{Myr}$, pointing at a synchronous growth of bulge and disk components. At masses $>10^{11}~M_{\rm sun}$, the sSFR profiles are typically centrally-suppressed by a factor of $\sim10$ relative to the galaxy outskirts. With total central obscuration disfavored, this indicates that at least a fraction of massive $z\sim2$ SFMS galaxies have started their inside-out star-formation quenching that will move them to the quenched sequence. In combination with other observations, galaxies above and below the ridge of the SFMS relation have respectively centrally-enhanced and centrally-suppressed sSFRs relative to their outskirts, supporting a picture where bulges are built due to gas `compaction' that leads to a high central SFR as galaxies move towards the upper envelope of SFMS.

Hyper Suprime-Cam (HSC) is a wide-field imaging camera on the prime focus of the 8.2m Subaru telescope on the summit of Maunakea in Hawaii. A team of scientists from Japan, Taiwan and Princeton University is using HSC to carry out a 300-night multi-band imaging survey of the high-latitude sky. The survey includes three layers: the Wide layer will cover 1400 deg$^2$ in five broad bands ($grizy$), with a $5\,\sigma$ point-source depth of $r \approx 26$. The Deep layer covers a total of 26~deg$^2$ in four fields, going roughly a magnitude fainter, while the UltraDeep layer goes almost a magnitude fainter still in two pointings of HSC (a total of 3.5 deg$^2$). Here we describe the instrument, the science goals of the survey, and the survey strategy and data processing. This paper serves as an introduction to a special issue of the Publications of the Astronomical Society of Japan, which includes a large number of technical and scientific papers describing results from the early phases of this survey.

We present new results from near-infrared spectroscopy with Keck/MOSFIRE of [OIII]-selected galaxies at $z\sim3.2$. With our $H$ and $K$-band spectra, we investigate the interstellar medium (ISM) conditions, such as ionization states and gas metallicities. [OIII] emitters at $z\sim3.2$ show a typical gas metallicity of $\mathrm{12+log(O/H) = 8.07\pm0.07}$ at $\mathrm{log(M_*/M_\odot) \sim 9.0-9.2}$ and $\mathrm{12+log(O/H) = 8.31\pm0.04}$ at $\mathrm{log(M_*/M_\odot) \sim 9.7-10.2}$ when using the empirical calibration method. We compare the [OIII] emitters at $z\sim3.2$ with UV-selected galaxies and Ly$\alpha$ emitters at the same epoch and find that the [OIII]-based selection does not appear to show any systematic bias in the selection of star-forming galaxies. Moreover, comparing with star-forming galaxies at $z\sim2$ from literature, our samples show similar ionization parameters and gas metallicities as those obtained by the previous studies using the same calibration method. We find no strong redshift evolution in the ISM conditions between $z\sim3.2$ and $z\sim2$. Considering that the star formation rates at a fixed stellar mass also do not significantly change between the two epochs, our results support the idea that the stellar mass is the primary quantity to describe the evolutionary stages of individual galaxies at $z>2$.

The Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) is a three-layered imaging survey aimed at addressing some of the most outstanding questions in astronomy today, including the nature of dark matter and dark energy. The survey has been awarded 300 nights of observing time at the Subaru Telescope and it started in March 2014. This paper presents the first public data release of HSC-SSP. This release includes data taken in the first 1.7 years of observations (61.5 nights) and each of the Wide, Deep, and UltraDeep layers covers about 108, 26, and 4 square degrees down to depths of i~26.4, ~26.5, and ~27.0 mag, respectively (5sigma for point sources). All the layers are observed in five broad bands (grizy), and the Deep and UltraDeep layers are observed in narrow bands as well. We achieve an impressive image quality of 0.6 arcsec in the i-band in the Wide layer. We show that we achieve 1-2 per cent PSF photometry (rms) both internally and externally (against Pan-STARRS1), and ~10 mas and 40 mas internal and external astrometric accuracy, respectively. Both the calibrated images and catalogs are made available to the community through dedicated user interfaces and database servers. In addition to the pipeline products, we also provide value-added products such as photometric redshifts and a collection of public spectroscopic redshifts. Detailed descriptions of all the data can be found online. The data release website is https://hsc-release.mtk.nao.ac.jp/.

We study the history from $z\sim2$ to $z\sim0$ of the stellar mass assembly of quiescent and star-forming galaxies in a spatially resolved fashion. For this purpose we use multi-wavelength imaging data from the Hubble Space Telescope (HST) over the GOODS fields and the Sloan Digital Sky Survey (SDSS) for the local population. We present the radial stellar mass surface density profiles of galaxies with $M_{\ast}>10^{10} M_{\odot}$, corrected for mass-to-light ratio ($M_{\ast}/L$) variations, and derive the half-mass radius ($R_{m}$), central stellar mass surface density within 1 kpc ($\Sigma_{1}$) and surface density at $R_{m}$ ($\Sigma_{m}$) for star-forming and quiescent galaxies and study their evolution with redshift. At fixed stellar mass, the half-mass sizes of quiescent galaxies increase from $z\sim2$ to $z\sim0$ by a factor of $\sim3-5$, whereas the half-mass sizes of star-forming galaxies increase only slightly, by a factor of $\sim2$. The central densities $\Sigma_{1}$ of quiescent galaxies decline slightly (by a factor of $\lesssim1.7$) from $z\sim2$ to $z\sim0$, while for star-forming galaxies $\Sigma_{1}$ increases with time, at fixed mass. We show that the central density $\Sigma_{1}$ has a tighter correlation with specific star-formation rate (sSFR) than $\Sigma_{m}$ and for all masses and redshifts galaxies with higher central density are more prone to be quenched. Reaching a high central density ($\Sigma_{1} \gtrsim 10^{10} M_{\odot} \mathrm{kpc}^2$) seems to be a prerequisite for the cessation of star formation, though a causal link between high $\Sigma_{1}$ and quenching is difficult to prove and their correlation can have a different origin.

We present the deep and wide $V$ and $I_c$ photometry of the Sextans dwarf spheroidal galaxy (dSph) taken by Suprime-Cam imager on the Subaru Telescope, which extends out to the tidal radius. The colour-magnitude diagram (CMD) reaches two magnitudes below the main sequence (MS) turn-off, showing a steep red giant branch, blue and red horizontal branch (HB), sub-giant branch (SGB), MS, and blue stragglers (BS). We construct the radial profile of each evolutionary phase and demonstrate that blue HB stars are more spatially extended, while red HB stars are more centrally concentrated than the other components. The colour distribution of SGB stars also varies with the galactocentric distance; the inner SGB stars shift bluer than those in the outskirt. The radial differences in the CMD morphology indicate the existence of the age gradient. The relatively younger stars ($\sim10$ Gyr) are more centrally concentrated than the older ones ($\sim13$ Gyr). The spatial contour maps of stars in different age bins also show that the younger population has higher concentration and higher ellipticity than the older one. We also detect the centrally concentrated bright BS stars, the number of which is consistent with the idea that a part of these stars belongs to the remnant of a disrupted star cluster discovered in the previous spectroscopic studies.

We present a sample of 34 spectroscopically confirmed BzK-selected ~1e11 Msun quiescent galaxies (pBzK) in the COSMOS field. The targets were initially observed with VIMOS on the VLT to facilitate the calibration of the photometric redshifts of massive galaxies at z >~ 1.5. Here we describe the reduction and analysis of the data, and the spectrophotometric properties of these pBzK galaxies. In particular, using a spatially resolved median 2D spectrum, we find that the fraction of stellar populations with ages <1 Gyr is at least 3 times higher in the outer regions of the pBzK galaxies than in their cores. This results in a mild age gradient of ~<0.4 Gyr over ~6 kpc and suggests either the occurrence of widespread rejuvenation episodes or that inside-out quenching played a role in the passivization of this galaxy population. We also report on low-level star formation rates derived from the [OII]3727A emission line, with SFR_OII = 3.7-4.5 Msun/yr. This estimate is confirmed by an independent measurement on a separate sample of similarly-selected quiescent galaxies in the COSMOS field, using stacked far-infrared data (SFR_FIR = 2-4 Msun/yr). This second, photometric sample also displays significant excess at 1.4 GHz, suggestive of the presence of radio-mode AGN activity.

We investigate the physical conditions of ionized gas in high-z star-forming galaxies using diagnostic diagrams based on the rest-frame optical emission lines. The sample consists of 701 galaxies with an Ha detection at $1.4\lesssim z\lesssim1.7$, from the FMOS-COSMOS survey, that represent the normal star-forming population over the stellar mass range $10^{9.6} \lesssim M_\ast/M_\odot \lesssim 10^{11.6}$ with those at $M_\ast>10^{11}~M_\odot$ being well sampled. We confirm an offset of the average location of star-forming galaxies in the BPT diagram ([OIII]/Hb vs. [NII]/Ha), primarily towards higher [OIII]/Hb, compared with local galaxies. Based on the [SII] ratio, we measure an electron density ($n_e=220^{+170}_{-130}~\mathrm{cm^{-3}}$), that is higher than that of local galaxies. Based on comparisons to theoretical models, we argue that changes in emission-line ratios, including the offset in the BPT diagram, are caused by a higher ionization parameter both at fixed stellar mass and at fixed metallicity with additional contributions from a higher gas density and possibly a hardening of the ionizing radiation field. Ionization due to AGNs is ruled out as assessed with Chandra. As a consequence, we revisit the mass-metallicity relation using [NII]/Ha and a new calibration including [NII]/[SII] as recently introduced by Dopita et al. Consistent with our previous results, the most massive galaxies ($M_\ast\gtrsim10^{11}~M_\odot$) are fully enriched, while those at lower masses have metallicities lower than local galaxies. Finally, we demonstrate that the stellar masses, metallicities and star formation rates of the FMOS sample are well fit with a physically-motivated model for the chemical evolution of star-forming galaxies.

We use HST/WFC3 imaging to study the red population in the IR-selected, X-ray detected, low-mass cluster Cl J1449+0856 at z=2, one of the few bona-fide established clusters discovered at this redshift, and likely a typical progenitor of an average massive cluster today. This study explores the presence and significance of an early red sequence in the core of this structure, investigating the nature of red sequence galaxies, highlighting environmental effects on cluster galaxy populations at high redshift, and at the same time underlining similarities and differences with other distant dense environments. Our results suggest that the red population in the core of Cl J1449+0856 is made of a mixture of quiescent and dusty star-forming galaxies, with a seedling of the future red sequence already growing in the very central cluster region, and already characterising the inner cluster core with respect to lower density environments. On the other hand, the color-magnitude diagram of this cluster is definitely different from that of lower-redshift (z<1) clusters, as well as of some rare particularly evolved massive clusters at similar redshift, and it is suggestive of a transition phase between active star formation and passive evolution occurring in the proto-cluster and established lower-redshift cluster regimes.

We investigate the origin of the evolution of the population-averaged size of quenched galaxies (QGs) through a spectroscopic analysis of their stellar ages. The two most favoured scenarios for this evolution are either the size growth of individual galaxies through a sequence of dry minor merger events, or the addition of larger, newly quenched galaxies to the pre-existing population (i.e., a progenitor bias effect). We use the 20k zCOSMOS-bright spectroscopic survey to select bona fide quiescent galaxies at 0.2<z<0.8. We stack their spectra in bins of redshift, stellar mass and size to compute stellar population parameters in these bins through fits to the rest-frame optical spectra and through Lick spectral indices. We confirm a change of behaviour in the size-age relation below and above the ~10^11 MSun stellar mass scale: In our 10.5 < log M*/MSun < 11 mass bin, over the entire redshift window, the stellar populations of the largest galaxies are systematically younger than those of the smaller counterparts, pointing at progenitor bias as the main driver of the observed average size evolution at sub-10^11 MSun masses. In contrast, at higher masses, there is no clear trend in age as a function of galaxy size, supporting a substantial role of dry mergers in increasing the sizes of these most massive QGs with cosmic time. Within the errors, the [alpha/Fe] abundance ratios of QGs are (i) above-solar over the entire redshift range of our analysis, hinting at universally short timescales for the buildup of the stellar populations of QGs, and (ii) similar at all masses and sizes, suggesting similar (short) timescales for the whole QG population and strengthening the role of mergers in the buildup of the most massive QGs in the Universe.

We present the discovery of a giant $\gtrsim$100~kpc Ly$\alpha$ nebula detected in the core of the X-ray emitting cluster CL~J1449+0856 at $z=1.99$ through Keck/LRIS narrow-band imaging. This detection extends the known relation between Ly$\alpha$ nebulae and overdense regions of the Universe to the dense core of a $5-7\times10^{13}$ M$_{\odot}$ cluster. The most plausible candidates to power the nebula are two Chandra-detected AGN host cluster members, while cooling from the X-ray phase and cosmological cold flows are disfavored primarily because of the high Ly$\alpha$ to X-ray luminosity ratio ($L_{\mathrm{Ly\alpha}}/L_{\mathrm{X}} \approx0.3$, $\gtrsim10-1000\times$ higher than in local cool-core clusters) and by current modeling. Given the physical conditions of the Ly$\alpha$-emitting gas and the possible interplay with the X-ray phase, we argue that the Ly$\alpha$ nebula would be short-lived ($\lesssim10$ Myr) if not continuously replenished with cold gas at a rate of $\gtrsim1000$ M$_{\odot}$ yr$^{-1}$. We investigate the possibility that cluster galaxies supply the required gas through outflows and we show that their total mass outflow rate matches the replenishment necessary to sustain the nebula. This scenario directly implies the extraction of energy from galaxies and its deposition in the surrounding intracluster medium, as required to explain the thermodynamic properties of local clusters. We estimate an energy injection of the order of $\thickapprox2$ keV per particle in the intracluster medium over a $2$ Gyr interval. In our baseline calculation AGN provide up to $85$% of the injected energy and 2/3 of the mass, while the rest is supplied by supernovae-driven winds.

We present new Keck/MOSFIRE K-band spectroscopy for a sample of 14 faint, X-ray-selected active galactic nuclei (AGNs) in the COSMOS field. The data cover the spectral region surrounding the broad Balmer emission lines, which enables the estimation of black hole masses (M_BH) and accretion rates (in terms of L/L_Edd). We focus on 10 AGN at z~3.3, where we observe the Hbeta spectral region, while for the other four z~2.4 sources we use the Halpha broad emission line. Compared with previous detailed studies of unobscured AGNs at these high redshifts, our sources are fainter by an order of magnitude, corresponding to number densities of order ~10^-6--10^-5 Mpc^-3. The lower AGN luminosities also allow for a robust identification of the host galaxy emission, necessary to obtain reliable intrinsic AGN luminosities, BH masses and accretion rates. We find the AGNs in our sample to be powered by supermassive black holes (SMBHs) with a typical mass of M_BH~5*10^8 M_Sol -- significantly lower than the higher-luminosity, rarer quasars reported in previous studies. The accretion rates are in the range L/L_Edd~0.1-0.5, with an evident lack of sources with lower L/L_Edd (and higher M_BH), as found in several studies of faint AGNs at intermediate redshifts. Based on the early growth expected for the SMBHs in our sample, we argue that a significant population of faint z~6 AGNs, with M_BH~10^6 M_Sol, should be detectable in the deepest X-ray surveys available, but this is not observed. We discuss several possible explanations for the apparent absence of such a population, concluding that the most probable scenario involves an evolution in source obscuration and/or radiative efficiencies.

We study the relationship between stellar mass, star formation rate (SFR),ionization state, and gas-phase metallicity for a sample of 41 normal star-forming galaxies at $3 \lesssim z \lesssim 3.7$. The gas-phase oxygen abundance, ionization parameter, and electron density of ionized gas are derived from rest-frame optical strong emission lines measured on near-infrared spectra obtained with Keck/MOSFIRE. We remove the effect of these strong emission lines in the broad-band fluxes to compute stellar masses via spectral energy distribution fitting, while the SFR is derived from the dust-corrected ultraviolet luminosity. The ionization parameter is weakly correlated with the specific SFR, but otherwise the ionization parameter and electron density do not correlate with other global galaxy properties such as stellar mass, SFR, and metallicity. The mass-metallicity relation (MZR) at $z\simeq3.3$ shows lower metallicity by $\simeq 0.7$ dex than that at $z=0$ at the same stellar mass. Our sample shows an offset by $\simeq 0.3$ dex from the locally defined mass-metallicity-SFR relation, indicating that simply extrapolating such relation to higher redshift may predict an incorrect evolution of MZR. Furthermore, within the uncertainties we find no SFR-metallicity correlation, suggesting a less important role of SFR in controlling the metallicity at high redshift. We finally investigate the redshift evolution of the MZR by using the model by Lilly et al. (2013), finding that the observed evolution from $z=0$ to $z\simeq3.3$ can be accounted for by the model assuming a weak redshift evolution of the star formation efficiency.

We measure a relation between the depth of four prominent rest-UV absorption complexes and metallicity for local galaxies and verify it up to z~3. We then apply this relation to a sample of 224 galaxies at 3.5 < z < 6.0 (<z> = 4.8) in COSMOS, for which unique UV spectra from DEIMOS and accurate stellar masses from SPLASH are available. The average galaxy population at z~5 and log(M/Msun) > 9 is characterized by 0.3-0.4 dex (in units of 12+log(O/H)) lower metallicities than at z~2, but comparable to z~3.5. We find galaxies with weak/no Ly-alpha emission to have metallicities comparable to z~2 galaxies and therefore may represent an evolved sub-population of z~5 galaxies. We find a correlation between metallicity and dust in good agreement with local galaxies and an inverse trend between metallicity and star-formation rate (SFR) consistent with observations at z~2. The relation between stellar mass and metallicity (MZ relation) is similar to z~3.5, however, there are indications of it being slightly shallower, in particular for the young, Ly-alpha emitting galaxies. We show that, within a "bathtub" approach, a shallower MZ relation is expected in the case of a fast (exponential) build-up of stellar mass with an e-folding time of 100-200 Myr. Due to this fast evolution, the process of dust production and metal enrichment as a function of mass could be more stochastic in the first billion years of galaxy formation compared to later times.

Local starbursts have a higher efficiency of converting gas into stars, as compared to typical star-forming galaxies at a given stellar mass, possibly indicative of different modes of star formation. With the peak epoch of galaxy formation occurring at z > 1, it remains to be established whether such an efficient mode of star formation is occurring at high-redshift. To address this issue, we measure the molecular gas content of seven high-redshift (z ~ 1.6) starburst galaxies with the Atacama Large (sub-)Millimeter Array and IRAM/Plateau de Bure Interferometer. Our targets are selected from the sample of Herschel far-infrared detected galaxies having star formation rates (~300-800 Msolar/yr) elevated (>4x) above the star-forming `main sequence', and included in the FMOS-COSMOS near-infrared spectroscopic survey of star-forming galaxies at z ~ 1.6 with Subaru. We detect CO emission in all cases at high levels of significance, indicative of high gas fractions (~30-50%). Even more compelling, we firmly establish with a clean and systematic selection that starbursts, identified as main-sequence outliers, at high redshift generally have a lower ratio of CO to total infrared luminosity as compared to typical main-sequence star-forming galaxies, although with a smaller offset than expected based on past studies of local starbursts. We put forward a hypothesis that there exists a continuous increase in star formation efficiency with elevation from the main sequence with galaxy mergers as a possible physical driver. Along with a heightened star formation efficiency, our high-redshift sample is similar in other respects to local starbursts such as being metal rich and having a higher ionization state of the interstellar medium.

We present a spectroscopic survey of galaxies in the COSMOS field using the Fiber Multi-Object Spectrograph (FMOS), a near-infrared instrument on the Subaru Telescope. Our survey is specifically designed to detect the Halpha emission line that falls within the H-band (1.6-1.8 um) spectroscopic window from star-forming galaxies with 1.4 < z < 1.7 and M_stellar>~10^10 Msolar. With the high multiplex capability of FMOS, it is now feasible to construct samples of over one thousand galaxies having spectroscopic redshifts at epochs that were previously challenging. The high-resolution mode (R~2600) effectively separates Halpha and [NII]6585 thus enabling studies of the gas-phase metallicity and photoionization state of the interstellar medium. The primary aim of our program is to establish how star formation depends on stellar mass and environment, both recognized as drivers of galaxy evolution at lower redshifts. In addition to the main galaxy sample, our target selection places priority on those detected in the far-infrared by Herschel/PACS to assess the level of obscured star formation and investigate, in detail, outliers from the star formation rate - stellar mass relation. Galaxies with Halpha detections are followed up with FMOS observations at shorter wavelengths using the J-long (1.11-1.35 um) grating to detect Hbeta and [OIII]5008 that provides an assessment of extinction required to measure star formation rates not hampered by dust, and an indication of embedded Active Galactic Nuclei. With 460 redshifts measured from 1153 spectra, we assess the performance of the instrument with respect to achieving our goals, discuss inherent biases in the sample, and detail the emission-line properties. Our higher-level data products, including catalogs and spectra, are available to the community.

We investigate the stellar population properties of a sample of 24 massive quenched galaxies at $1.25<z_\mathrm{spec}<2.09$ identified in the COSMOS field with our Subaru/MOIRCS near-IR spectroscopic observations. Tracing the stellar population properties as close to their major formation epoch as possible, we try to put constraints on the star formation history, post-quenching evolution, and possible progenitor star-forming populations for such massive quenched galaxies. By using a set of Lick absorption line indices on a rest-frame optical composite spectrum, the average age, metallicity [Z/H], and $\alpha$-to-iron element abundance ratio [$\alpha$/Fe] are derived as $\log(\mathrm{age}/\mathrm{Gyr})=0.04_{-0.08}^{+0.10}$, $\mathrm{[Z/H]}=0.24_{-0.14}^{+0.20}$, and $[\alpha/\mathrm{Fe}]=0.31_{-0.12}^{+0.12}$, respectively. If our sample of quenched galaxies at $\langle z \rangle = 1.6$ is evolved passively to $z=0$, their stellar population properties will align in excellent agreement with local counterparts at similar stellar velocity dispersions, which qualifies them as progenitors of local massive early-type galaxies. Redshift evolution of stellar population ages in quenched galaxies combined with low redshift measurements from the literature suggests a formation redshift of $z_\mathrm{f} \sim 2.3$ around which the bulk of stars in these galaxies have been formed. The measured [$\alpha$/Fe] value indicates a star formation timescale of $\lesssim 1$ Gyr, which can be translated into a specific star formation rate of $\simeq 1\,\mathrm{Gyr}^{-1}$ prior to quenching. Based on these findings, we discuss identifying possible progenitor star-forming galaxies at $z \simeq 2.3$. We identify normal star-forming galaxies, i.e, those on the star-forming main sequence, followed by a rapid quenching event, as likely precursors of the quenched galaxies at $\langle z \rangle = 1.6$ presented here.

Gas and dust in star-forming galaxies at the peak epoch of galaxy assembly are presently the topic of intense study, but little is known about the interstellar medium (ISM) of distant, passively evolving galaxies. We report on a deep 3 mm-band search with IRAM/PdBI for molecular (H$_2$) gas in a massive ($M_{\star}{\sim}6{\times}10^{11}M_{\odot}$) elliptical galaxy at z=1.4277, the first observation of this kind ever attempted. We place a 3$\sigma$ upper limit of 0.32 Jy km/s on the flux of the CO($J$=2$\rightarrow$1) line or $L'_{\rm CO}$$<$8.8$\times$10$^{9}$ K km/s pc$^2$, assuming a disk-like CO-morphology and a circular velocity scaling with the stellar velocity dispersion as in local early-type galaxies (ETGs). This translates to an H$_2$ mass of $<$3.9$\times$10$^{10}$($\alpha_{\rm CO}$/4.4)$M_{\odot}$ or a gas fraction of $\lesssim$6% assuming a Salpeter initial mass function (IMF) and an ISM dominated by H$_2$, as observed in many local, high-mass ellipticals. This low value approaches that of local ETGs, suggesting that the low star formation activity in massive, high-z passive galaxies reflects a true dearth of gas and a lesser role for inhibitive mechanisms like morphological quenching.

We use deep panchromatic datasets in the GOODS-N field, from GALEX to the deepest Herschel far-infrared and VLA radio continuum imaging, to explore, using mass-complete samples, the evolution of the star formation activity and dust attenuation of star-forming galaxies to z~4. Our main results can be summarized as follows: i) the slope of the SFR-M correlation is consistent with being constant, and equal to ~0.8 at least up to z~1.5, while its normalization keeps increasing with redshift; ii) for the first time here we are able to explore the FIR-radio correlation for a mass-selected sample of star-forming galaxies: the correlation does not evolve up to z~4; iii) we confirm that galaxy stellar mass is a robust proxy for UV dust attenuation in star-forming galaxies, with more massive galaxies being more dust attenuated, strikingly we find that this attenuation relation evolves very weakly with redshift, the amount of dust attenuation increasing by less than 0.3 magnitudes over the redshift range [0.5-4] for a fixed stellar mass, as opposed to a tenfold increase of star formation rate; iv) the correlation between dust attenuation and the UV spectral slope evolves in redshift, with the median UV spectral slope of star-forming galaxies becoming bluer with redshift. By z~3, typical UV slopes are inconsistent, given the measured dust attenuation, with the predictions of commonly used empirical laws. Finally, building on existing results, we show that gas reddening is marginally larger (by a factor of around 1.3) than stellar reddening at all redshifts probed, and also that the amount of dust attenuation at a fixed ISM metallicity increases with redshift. We speculate that our results support evolving ISM conditions of typical star-forming galaxies such that at z~1.5 Main Sequence galaxies have ISM conditions getting closer to those of local starbursts.

When the cosmic star formation history peaks (z ~ 2), galaxies vigorously fed by cosmic reservoirs are gas dominated and contain massive star-forming clumps, thought to form by violent gravitational instabilities in highly turbulent gas-rich disks. However, a clump formation event has not been witnessed yet, and it is debated whether clumps survive energetic feedback from young stars, thus migrating inwards to form galaxy bulges. Here we report spatially resolved spectroscopy of a bright off-nuclear emission line region in a galaxy at z = 1.987. Although this region dominates the star formation in the galaxy disk, its stellar continuum remains undetected in deep imaging, revealing an extremely young (age < 10 Myr) massive clump, forming through the gravitational collapse of > 10$^9$ M$_{\odot}$ of gas. Gas consumption in this young clump is > 10 times faster than in the host galaxy, displaying high star formation efficiency during this phase, in agreement with our hydrodynamic simulations. The frequency of older clumps with similar masses coupled with our initial estimate of their formation rate (~ 2.5 Gyr$^{-1}$) supports long lifetimes (~ 500 Myr), favouring scenarios where clumps survive feedback and grow the bulges of present-day galaxies.

Most present-day galaxies with stellar masses $\geq10^{11}$ solar masses show no ongoing star formation and are dense spheroids. Ten billion years ago, similarly massive galaxies were typically forming stars at rates of hundreds solar masses per year. It is debated how star formation ceased, on which timescales, and how this "quenching" relates to the emergence of dense spheroids. We measured stellar mass and star-formation rate surface density distributions in star-forming galaxies at redshift 2.2 with $\sim1$ kiloparsec resolution. We find that, in the most massive galaxies, star formation is quenched from the inside out, on timescales less than 1 billion years in the inner regions, up to a few billion years in the outer disks. These galaxies sustain high star-formation activity at large radii, while hosting fully grown and already quenched bulges in their cores.

Even 10 billion years ago, the cores of the first galaxy clusters are often found to host a characteristic population of massive galaxies with already suppressed star formation. Here we search for distant cluster candidates at z~2 using massive passive galaxies as tracers. With a sample of ~40 spectroscopically confirmed passive galaxies at 1.3<z<2.1, we tune photometric redshifts of several thousands passive sources in the full 2 sq.deg. COSMOS field. This allows us to map their density in redshift slices, probing the large scale structure in the COSMOS field as traced by passive sources. We report here on the three strongest passive galaxy overdensities that we identify in the redshift range 1.5<z<2.5. While the actual nature of these concentrations is still to be confirmed, we discuss their identification procedure, and the arguments supporting them as candidate galaxy clusters (likely mid-10^13 M_sun range). Although this search approach is likely biased towards more evolved structures, it has the potential to select still rare, cluster-like environments close to their epoch of first appearance, enabling new investigations of the evolution of galaxies in the context of structure growth.

We investigate the environmental effect on the metal enrichment of star-forming galaxies (SFGs) in the farthest spectroscopically confirmed and X-ray detected cluster, CL~J1449+0856 at $z=1.99$. We combined HST/WFC3~G141 slitless spectroscopic data, our 13-band photometry, and a recent Subaru/MOIRCS near infrared spectroscopic follow-up to constrain the physical properties of SFGs in CL~J1449+0856 and in a mass-matched field sample. After a conservative active galactic nuclei (AGN) removal, stacking individual MOIRCS spectra of $6$ ($31$) sources in the cluster (field) in the mass range $10 \leq \mathrm{log(M/M_{\odot}}) \leq 11$, we find a $\sim4\sigma$ significant lower [N II]/H$\alpha$ ratio in the cluster than in the field. Stacking a subsample of $16$ field galaxies with H$\beta$ and [O III] in the observed range, we measure a [O III]/H$\beta$ ratio fully compatible with the cluster value. Converting these ratios into metallicities, we find that the cluster SFGs are up to $0.25\,\mathrm{dex}$ poorer in metals than their field counterparts, depending on adopted calibrations. The low metallicity in cluster sources is confirmed using alternative indicators. Furthermore, we observe a significantly higher H$\alpha$ luminosity and equivalent width in the average cluster spectrum than in the field. This is likely due to enhanced specific star formation rate, even if lower dust reddening and/or an uncertain environmental dependence of the continuum-to-nebular emission differential reddening may play a role. Our findings might be explained by the accretion of pristine gas around galaxies at $z=2$ and from cluster-scale reservoirs, possibly connected with a phase of rapid halo mass assembly at $z>2$ and of high galaxy merging rate.

We investigate the relationships between stellar mass, gas-phase oxygen abundance (metallicity), star formation rate, and dust content of star-forming galaxies at z$\sim$1.6 using Subaru/FMOS spectroscopy in the COSMOS field. The mass-metallicity relation at $z\sim1.6$ is steeper than the relation observed in the local Universe. The steeper MZ relation at $z\sim1.6$ is mainly due to evolution in the stellar mass where the MZ relation begins to turnover and flatten. This turnover mass is 1.2 dex larger at $z\sim1.6$. The most massive galaxies at $z\sim1.6$ ($\sim 10^{11}M_\odot$) are enriched to the level observed in massive galaxies in the local Universe. The mass-metallicity relation we measure at $z\sim1.6$ supports the suggestion of an empirical upper metallicity limit that does not significantly evolve with redshift. We find an anti-correlation between metallicity and star formation rate for galaxies at a fixed stellar mass at $z\sim1.6$ which is similar to trends observed in the local Universe. We do not find a relation between stellar mass, metallicity and star formation rate that is independent of redshift; our data suggest that there is redshift evolution in this relation. We examine the relation between stellar mass, metallicity and dust extinction. We find that at a fixed stellar mass dustier galaxies tend to be more metal rich. From examination of the stellar masses, metallicities, SFRs and dust extinctions we conclude that stellar mass is most closely related to dust extinction.

We present the first results from a near-IR spectroscopic survey of the COSMOS field, using the Fiber Multi-Object Spectrograph on the Subaru telescope, designed to characterize the star-forming galaxy population at $1.4<z<1.7$. The high-resolution mode is implemented to detect H$\alpha$ in emission between $1.6{\rm -}1.8 \mathrm{\mu m}$ with $f_{\rm H\alpha}\gtrsim4\times10^{-17}$ erg cm$^{-2}$ s$^{-1}$. Here, we specifically focus on 271 sBzK-selected galaxies that yield a H$\alpha$ detection thus providing a redshift and emission line luminosity to establish the relation between star formation rate and stellar mass. With further $J$-band spectroscopy for 89 of these, the level of dust extinction is assessed by measuring the Balmer decrement using co-added spectra. We find that the extinction ($0.6\lesssim A_\mathrm{H\alpha} \lesssim 2.5$) rises with stellar mass and is elevated at high masses compared to low-redshift galaxies. Using this subset of the spectroscopic sample, we further find that the differential extinction between stellar and nebular emission \hbox{$E_\mathrm{star}(B-V)/E_\mathrm{neb}(B-V)$} is 0.7--0.8, dissimilar to that typically seen at low redshift. After correcting for extinction, we derive an H$\alpha$-based main sequence with a slope ($0.81\pm0.04$) and normalization similar to previous studies at these redshifts.

We present deep Hubble Space Telescope Wide Field Camera 3 slitless spectroscopic observations of the distant cluster Cl J1449+0856. These cover a single pointing with 18 orbits of G141 spectroscopy and F140W imaging, allowing us to derive secure redshifts down to m_140~25.5 AB and 3sigma line fluxes of 5*10^(-18) erg/s/cm^2. In particular, we were able to spectroscopically confirm 12 early-type galaxies in the field up to z~3, 6 of which in the cluster core, which represents the first direct spectroscopic confirmation of passive galaxies in a z=2 cluster environment. With 140 redshifts in a ~6 arcmin^2 field, we can trace the spatial and redshift galaxy distribution in the cluster core and background field. We find two strong peaks at z=2.00 and z=2.07, where only one was seen in our previously published ground-based data. Thanks to the spectroscopic confirmation of the cluster ETGs, we can now re-evaluate the redshift of Cl J1449+0856 at z=2.00, rather than z=2.07, with the background overdensity being revealed to be sparse and "sheet"-like. This presents an interesting case of chance alignment of two close yet unrelated structures, each one preferentially selected by different observing strategies. With 6 quiescent or early-type spectroscopic members and 20 star-forming ones, Cl J1449+0856 is now reliably confirmed to be at z=2.00. The identified members can now allow for a detailed study of galaxy properties in the densest environment at z=2.

Abridged. We use COSMOS to study in a self-consistent way the change in the number densities of quenched early-type galaxies (Q-ETGs) of a given size over the interval 0.2 < z < 1.0 to study the claimed size evolution of these galaxies. At 10^10.5<Mgalaxy<10^11 Msun, we see no change in the number density of compact Q-ETGs, while at >10^11 Msun we find a decrease by 30%. In both mass bins, the increase of the median sizes of Q-ETGs with time is primarily caused by the addition to the size function of larger and more diffuse Q-ETGs. At all masses, compact Q-ETGs become systematically redder towards later epochs, with a (U-V) difference consistent with passive evolution of their stellar populations, indicating that they are a population that does not appreciably evolve in size. At all epochs, the larger Q-ETGs (at least in the lower mass bin) have average rest-frame colors systematically bluer than those of the more compact Q-ETGs, suggesting that the former are younger than the latter. The idea that new, large, Q-ETGs are responsible for the observed growth in the median size of the population at a given mass is supported by the sizes and number of the star-forming galaxies that are expected to be progenitors of the new Q-ETGs over the same period. In the low mass bin, the new Q-ETG have 30% smaller sizes than their star-forming progenitors. This is likely due to the fading of their disks after they cease star-formation. Comparison with higher z shows that the median size of newly-quenched galaxies roughly scales, at constant mass, as (1+z)^-1. The dominant cause of the size evolution seen in the Q-ETG population is thus that the average sizes of individual Q-ETGs scale with the average density of the Universe at the time when they were quenched, with subsequent size changes in individual objects through eg merging of secondary importance, especially at masses <10^11 Msun.

We estimate the galaxy stellar mass function and stellar mass density for star-forming and quiescent galaxies with 0.2<z<4. We construct a deep K<24 sample of 220000 galaxies selected using the UltraVISTA DR1 data release. Our analysis is based on precise 30-band photometric redshifts. By comparing these photometric redshifts with 10800 spectroscopic redshifts from the zCOSMOS bright and faint surveys, we find a precision of sigma(dz/(1+z))=0.008 at i<22.5 and sigma(dz/(1+zs))=0.03 at 1.5<z<4. We derive the stellar mass function and correct for the Eddington bias. We find a mass-dependent evolution of the global and star-forming populations. This mass-dependent evolution is a direct consequence of the star formation being quenched in galaxies more massive than M>10^10.7Msun. For the mass function of the quiescent galaxies, we do not find any significant evolution of the high-mass end at z<1; however we observe a clear flattening of the faint-end slope. From z~3 to z~1, the density of quiescent galaxies increases over the entire mass range. Their comoving stellar mass density increases by 1.6 dex between z~3 and z~1 and by less than 0.2dex at z<1. We infer the star formation history from the mass density evolution and we find an excellent agreement with instantaneous star formation rate measurements at z<1.5, while we find differences of 0.2dex at z>1.5 consistent with the expected uncertainties. We also develop a new method to infer the specific star formation rate from the mass function of star-forming galaxies. We find that the specific star formation rate of 10^10Msun galaxies increases continuously in the redshift range 1<z<4. Finally, we compare our results with a semi-analytical model and find that these models overestimate the density of low mass quiescent galaxies by an order of magnitude, while the density of low-mass star-forming galaxies is successfully reproduced.

We present a study of galaxy populations in the central region of the IRAC-selected, X-ray detected galaxy cluster Cl J1449+0856 at z=2. Based on a sample of spectroscopic and photometric cluster members, we investigate stellar populations and morphological structure of cluster galaxies over an area of ~0.7Mpc^2 around the cluster core. The cluster stands out as a clear overdensity both in redshift space, and in the spatial distribution of galaxies close to the center of the extended X-ray emission. The cluster core region (r<200 kpc) shows a clearly enhanced passive fraction with respect to field levels. However, together with a population of massive passive galaxies mostly with early-type morphologies, it also hosts massive actively star-forming, often highly dust-reddened sources. Close to the cluster center, a multi-component system of passive and star-forming galaxies could be the future BCG still assembling. We observe a clear correlation between passive stellar populations and an early-type morphology, in agreement with field studies at similar redshift. Passive early-type galaxies in this clusters are typically a factor 2-3 smaller than similarly massive early-types at z~0, but also on average larger by a factor ~2 than their field analogs at z~2, lending support to recent claims of an accelerated structural evolution in high-redshift dense environments. These results point towards the early formation of a population of massive galaxies, already evolved both in their structure and stellar populations, coexisting with still-actively forming massive galaxies in the central regions of young clusters 10 billion years ago.

We present the discovery of a massive, quiescent galaxy at z=2.99. We have obtained a HST/WFC3 spectrum of this object and measured its redshift from the detection of a deep 4000A break consistent with an old population and a high metallicity. By stellar population modeling of both its grism spectrum and broad-band photometry, we derive an age of ~0.7 Gyr, implying a formation redshift of z>4, and a mass >10^11 Msun. Although this passive galaxy is the most distant confirmed so far, we find that it is slightly less compact than other z>2 early-types of similar mass, being overall more analogous to those z~1.6 field early-type galaxies. The discovery of this object shows that early-type galaxies are detectable to at least z=3 and suggests that the diversity of structural properties found in z=1.4-2 ellipticals to earlier epochs could have its origin in a variety of formation histories among their progenitors.

We present a spectroscopic analysis of the metallicities, ages, and alpha-elements of the globular clusters (GCs) in the giant elliptical galaxy (gE) NGC 4636 in the Virgo cluster. Line indices of the GCs are measured from the integrated spectra obtained with Faint Object Camera and Spectrograph (FOCAS) on the Subaru 8.2 m telescope. We derive [Fe/H] values of 59 GCs based on the Brodie & Huchra method, and [Z/H], age, and [a/Fe] values of 33 GCs from the comparison of the Lick line indices with single stellar population models. The metallicity distribution of NGC 4636 GCs shows a hint of a bimodality with two peaks at [Fe/H]=-1.23 (sigma=0.32) and -0.35 (sigma=0.19). The age spread is large from 2 Gyr to 15 Gyr and the fraction of young GCs with age < 5 Gyr is about 27%. The [a/Fe] of the GCs shows a broad distribution with a mean value [a/Fe]~0.14 dex. The dependence of these chemical properties on the galactocentric radius is weak. We also derive the metallicities, ages, and [a/Fe] values for the GCs in other nearby gEs (M87, M49, M60, NGC 5128, NGC 1399, and NGC 1407) from the line index data in the literature using the same methods as used for NGC 4636 GCs. The metallicity distribution of GCs in the combined sample of seven gEs including NGC 4636 is found to be bimodal, supported by the KMM test with a significance level of >99.9%. All these gEs harbor some young GCs with ages less than 5 Gyr. The mean age of the metal-rich GCs ([Fe/H] > -0.9) is about 3 Gyr younger than that of the metal-poor GCs. The mean value of [a/Fe] of the gE GCs is smaller than that of the Milky Way GCs. We discuss these results in the context of GC formation in gEs.

[Abridged] We present the results of new near-IR spectroscopic observations of passive galaxies at z>1.4 in a concentration of BzK-selected galaxies in the COSMOS field. The observations have been conducted with Subaru/MOIRCS, and have resulted in absorption lines and/or continuum detection for 18 out of 34 objects. This allows us to measure spectroscopic redshifts for a sample almost complete to K(AB)=21. COSMOS photometric redshifts are found in fair agreement overall with the spectroscopic redshifts, with a standard deviation of ~0.05; however, ~30% of objects have photometric redshifts systematically underestimated by up to ~25%. We show that these systematic offsets in photometric redshifts can be removed by using these objects as a training set. All galaxies fall in four distinct redshift spikes at z=1.43, 1.53, 1.67 and 1.82, with this latter one including 7 galaxies. SED fits to broad-band fluxes indicate stellar masses in the range of ~4-40x10^10Msun and that star formation was quenched ~1 Gyr before the cosmic epoch at which they are observed. The spectra of several individual galaxies have allowed us to measure their Hdelta_F and Dn4000 indices, which confirms their identification as passive galaxies, as does a composite spectrum resulting from the coaddition of 17 individual spectra. The effective radii of the galaxies have been measured on the HST/ACS F814W image, confirming the coexistence at these redshifts of passive galaxies which are substantially more compact than their local counterparts with others that follow the local size-stellar mass relation. For the galaxy with best S/N spectrum we were able to measure a velocity dispersion of 270+/-105 km/s, indicating that this galaxy lies closely on the virial relation given its stellar mass and effective radius.

We take deep images of four ultra faint dwarf (UFD) galaxies, Canes Venatici I (CVn I), Bootes I (Boo I), Canes Venatici II (CVn II), and Leo IV using the Suprime-Cam on the Subaru Telescope. The colour-magitude diagrams (CMDs) extend below the main sequence turn-offs (MSTOs) and yield measurements of the ages of stellar populations. The stellar populations of faint three galaxies, Boo I, CVn II and Leo IV dwarf spheroidal galaxies (dSphs) are estimated as old as a Galactic globular cluster M92. We confirm that Boo I dSph has no intrinsic colour spread in the MSTO, and no spatial difference in the CMD morphology, which indicates that Boo I dSph is composed of an old single stellar population. One of the brightest UFDs, CVn I dSph, shows a relatively younger age (-12.6 Gyr) with respect to Boo I, CVn II, and Leo IV dSphs, and the distribution of red horizontal branch (HB) stars is more concentrated toward the centre than that of blue HB stars, suggesting that the galaxy contains complex stellar populations. Boo I and CVn I dSphs show the elongated and distorted shapes. CVn II dSph has the smallest tidal radius as a Milky Way satellite and has distorted shape, while Leo IV dSph shows less concentrated spherical shape. The simple stellar population of faint UFDs indicates that the gases in their progenitors were removed more effectively than those of brighter dSphs at an occurrence of their initial star formation. This is reasonable if the progenitors of UFDs belong to the less massive halos than those of brighter dSphs.

We have developed high dispersion VPH (volume phase holographic) grisms with zinc selenide (ZnSe) prisms for the cryogenic optical system of MOIRCS (Multi-Object near InfraRed Camera and Spectrograph) for Y-, J-, H- and K- band observations. We fabricated the VPH gratings using a hologram resin. After several heat cycles at between room temperature and 120 K, the VPH gratings were assembled to grisms by gluing with two ZnSe prisms. Several heat cycles were also carried out for the grisms before being installed into MOIRCS. We measured the efficiencies of the VPH grisms in a laboratory, and found them to be 70% - 82%. The performances obtained by observations of MOIRCS with the 8.2 m Subaru Telescope have been found to be very consistent with the results in the laboratory test. This is the first astronomical application of cryogenic VPH grisms.

We study the evolution of quiescent galaxies at 0.5<z<2.5 as a function of stellar mass, using very deep NIR imaging data taken with the Multi-Object Infrared Camera and Spectrograph on the Subaru Telescope in the GOODS-North region. The deep NIR data allow us to construct a stellar mass-limited sample of quiescent galaxies down to ~10^{10} Msun even at z~2 for the first time. We selected quiescent galaxies with age/tau>6 by performing SED fitting of the multi broad-band photometry from the U to Spitzer 5.8um bands with the population synthesis model of Bruzual & Charlot (2003) where exponentially decaying star formation histories are assumed. The number density of quiescent galaxies increases by a factor of ~3 from 1.0<z<1.5 to 0.5<z<1.0, and by a factor of ~10 from 1.5<z<2.5 to 0.5<z<1.0, while that of star-forming galaxies with age/tau<4 increases only by factors of ~2 and ~3 in the same redshift ranges. At 0.5<z<2.5, the low-mass slope of the stellar mass function of quiescent galaxies is alpha ~ 0 -- 0.6, which is significantly flatter than those of star-forming galaxies (alpha ~ -1.3 -- -1.5). As a result, the fraction of quiescent galaxies in the overall galaxy population increases with stellar mass in the redshift range. The fraction of quiescent galaxies at 10^{11}-10^{11.5} Msun increases from ~20-30% at z~2 to ~40-60% at z~0.75, while that at 10^{10}-10^{10.5} Msun increases from <~ 5% to ~15% in the same redshift range. These results could suggest that the quenching of star formation had been more effective in more massive galaxies at 1<~z<~2. Such a mass-dependent quenching could explain the rapid increase of the number density of ~M* galaxies relative to lower-mass galaxies at z >~ 1-1.5.

We present deep J-, H-, and Ks-band imaging data of the MOIRCS Deep Survey (MODS), which was carried out with Multi-Object Infrared Camera and Spectrograph (MOIRCS) mounted on the Subaru telescope in the GOODS-North region. The data reach 5sigma total limiting magnitudes for point sources of J=23.9, H=22.8, and Ks=22.8 (Vega magnitude) over 103 arcmin^2 (wide field). In 28 arcmin^2 of the survey area, which is ultra deep field of the MODS (deep field), the data reach the 5sigma depths of J=24.8, H=23.4, and Ks=23.8. The spatial resolutions of the combined images are FWHM ~ 0.6 arcsec and ~ 0.5 arcsec for the wide and deep fields in all bands, respectively. Combining the MODS data with the multi-wavelength public data taken with the HST, Spitzer, and other ground-based telescopes in the GOODS field, we construct a multi-wavelength photometric catalog of Ks-selected sources. Using the catalog, we present Ks-band number counts and near-infrared color distribution of the detected objects, and demonstrate some selection techniques with the NIR colors for high redshift galaxies. These data and catalog are publicly available via internet.

We report evidence of a fully established galaxy cluster at z=2.07, consisting of a ~20sigma overdensity of red, compact spheroidal galaxies spatially coinciding with extended X-ray emission detected with XMM-Newton. We use VLT VIMOS and FORS2 spectra and deep Subaru, VLT and Spitzer imaging to estimate the redshift of the structure from a prominent z=2.07 spectroscopic redshift spike of emission-line galaxies, concordant with the accurate 12-band photometric redshifts of the red galaxies. Using NICMOS and Keck AO observations, we find that the red galaxies have elliptical morphologies and compact cores. While they do not form a tight red sequence, their colours are consistent with that of a >1.3$~Gyr population observed at z~2.1. From an X-ray luminosity of .2*10^43 erg s^-1 and the stellar mass content of the red galaxy population, we estimate a halo mass of 5.3-8*10^13 Msun, comparable to the nearby Virgo cluster. These properties imply that this structure could be the most distant, mature cluster known to date and that X-ray luminous, elliptical-dominated clusters are already forming at substantially earlier epochs than previously known.

We investigate rest-frame near-infrared (NIR) morphologies of a sample of 139 galaxies with M_{s} >= 1 x 10^{10} M_{sun} at z=0.8-1.2 in the GOODS-North field using our deep NIR imaging data (MOIRCS Deep Survey, MODS). We focus on Luminous Infrared Galaxies (LIRGs), which dominate high star formation rate (SFR) density at z~1, in the sample identified by cross-correlating with the Spitzer/MIPS 24um source catalog. We perform two-dimensional light profile fitting of the z~1 galaxies in the Ks-band (rest-frame J-band) with a single component Sersic model. We find that at z~1, ~90% of LIRGs have low Sersic indices (n<2.5, similar to disk-like galaxies) in the Ks-band, and those disk-like LIRGs consist of ~60% of the whole disk-like sample above M_{s} >= 3 x 10^{10} M_{sun}. The z~1 disk-like LIRGs are comparable or ~20% small at a maximum in size compared to local disk-like galaxies in the same stellar mass range. If we examine rest-frame UV-optical morphologies using the HST/ACS images, the rest-frame B-band sizes of the z~1 disk-like galaxies are comparable to those of the local disk-like galaxies as reported by previous studies on size evolution of disk-like galaxies in the rest-frame optical band. Measuring color gradients (galaxy sizes as a function of wavelength) of the z~1 and local disk-like galaxies, we find that the z~1 disk-like galaxies have 3-5 times steeper color gradient than the local ones. Our results indicate that (i) more than a half of relatively massive disk-like galaxies at z~1 are in violent star formation epochs observed as LIRGs, and also (ii) most of those LIRGs are constructing their fundamental disk structure vigorously. The high SFR density in the universe at z~1 may be dominated by such star formation in disk region in massive galaxies.

We study the evolution of star formation activity of galaxies at 0.5<z<3.5 as a function of stellar mass, using very deep NIR data taken with Multi-Object Infrared Camera and Spectrograph (MOIRCS) on the Subaru telescope in the GOODS-North region. The NIR imaging data reach K ~ 23-24 Vega magnitude and they allow us to construct a nearly stellar mass-limited sample down to ~ 10^{9.5-10} Msun even at z~3. We estimated star formation rates (SFRs) of the sample with two indicators, namely, the Spitzer/MIPS 24um flux and the rest-frame 2800A luminosity. The SFR distribution at a fixed Mstar shifts to higher values with increasing redshift at 0.5<z<3.5. More massive galaxies show stronger evolution of SFR at z>~1. We found galaxies at 2.5<z<3.5 show a bimodality in their SSFR distribution, which can be divided into two populations by a constant SSFR of ~2 Gyr^{-1}. Galaxies in the low-SSFR group have SSFRs of ~ 0.5-1.0 Gyr^{-1}, while the high-SSFR population shows ~10 Gyr^{-1}. The cosmic SFRD is dominated by galaxies with Mstar = 10^{10-11} Msun at 0.5<z<3.5, while the contribution of massive galaxies with Mstar = 10^{11-11.5} Msun shows a strong evolution at z>1 and becomes significant at z~3, especially in the case with the SFR based on MIPS 24um. In galaxies with Mstar = 10^{10-11.5} Msun, those with a relatively narrow range of SSFR (<~1 dex) dominates the cosmic SFRD at 0.5<z<3.5. The SSFR of galaxies which dominate the SFRD systematically increases with redshift. At 2.5<z<3.5, the high-SSFR population, which is relatively small in number, dominates the SFRD. Major star formation in the universe at higher redshift seems to be associated with a more rapid growth of stellar mass of galaxies.

We present observations of a very massive galaxy at z=1.82 which show that its morphology, size, velocity dispersion and stellar population properties that are fully consistent with those expected for passively evolving progenitors of today's giant ellipticals. These findings are based on a deep optical rest-frame spectrum obtained with the Multi-Object InfraRed Camera and Spectrograph (MOIRCS) on the Subaru telescope of a high-z passive galaxy candidate (pBzK) from the COSMOS field, for which we accurately measure its redshift of z=1.8230 and obtain an upper limit on its velocity dispersion sigma_star<326 km/s. By detailed stellar population modeling of both the galaxy broad-band SED and the rest-frame optical spectrum we derive a star-formation-weighted age and formation redshift of t_sf~1-2 Gyr and z_form~2.5-4, and a stellar mass of M_star~(3-4)x10^{11} M_sun. This is in agreement with a virial mass limit of M_vir<7x10^{11}M_sun, derived from the measured sigma_star value and stellar half-light radius, as well as with the dynamical mass limit based on the Jeans equations. In contrast with previously reported super-dense passive galaxies at z~2, the present galaxy at z=1.82 appears to have both size and velocity dispersion similar to early-type galaxies in the local Universe with similar stellar mass. This suggests that z~2 massive and passive galaxies may exhibit a wide range of properties, then possibly following quite different evolutionary histories from z~2 to z=0.

Results are presented from NIR spectroscopy of a sample of BzK-selected, massive star-forming galaxies (sBzKs) at 1.5<z<2.3 that were obtained with OHS/CISCO at Subaru and with SINFONI at VLT. Among the 28 sBzKs observed, Ha emission was detected in 14 objects, and for 11 of them the [NII]6583 was also measured. Multiwavelength photometry was also used to derive stellar masses and extinction parameters, whereas Ha and [NII] have allowed us to estimate SFR, metallicities, ionization mechanisms, and dynamical masses. In order to enforce agreement between SFRs from Ha with those derived from rest-frame UV and MIR, additional obscuration for the emission lines (that originate in HII regions) was required compared to the extinction derived from the slope of the UV continuum. We have also derived the stellar mass-metallicity relation, as well as the relation between stellar mass and specific SFR, and compared them to the results in other studies. At a given stellar mass, the sBzKs appear to have been already enriched to metallicities close to those of local star-forming galaxies of similar mass. The sBzKs presented here tend to have higher metallicities compared to those of UV-selected galaxies, indicating that NIR selected galaxies tend to be a chemically more evolved population. The sBzKs show specific SFRs that are systematically higher, by up to ~2 orders of magnitude, compared to those of local galaxies of the same mass. The empirical correlations between stellar mass and metallicity, and stellar mass and specific SFR are then compared with those of population synthesis models constructed either with the simple closed-box assumption, or within an infall scenario. Within the assumptions that are built-in such models, it appears that a short timescale for the star-formation (~100 Myr) and large initial gas mass appear to be required if one wants to reproduce both relations simultaneously.

We present evidence for very high gas fractions and extended molecular gas reservoirs in normal, near-infrared selected (BzK) galaxies at z~1.5, based on multi-configuration CO[2-1] observations obtained at the IRAM PdBI. Six of the six galaxies observed were securely detected. High resolution observations resolve the CO emission in four of them, implying sizes of order of 6-11 kpc and suggesting the presence of rotation. The UV morphologies are consistent with clumpy, unstable disks, and the UV sizes are consistent with the CO sizes. The star formation efficiencies are homogeneously low and similar to local spirals - the resulting gas depletion times are ~0.5 Gyr, much higher than what is seen in high-z submm galaxies and quasars. The CO luminosities can be predicted to within 0.15 dex from the star formation rates and stellar masses, implying a tight correlation of the gas mass with these quantities. We use dynamical models of clumpy disk galaxies to derive dynamical masses. These models are able to reproduce the peculiar spectral line shapes of the CO emission. After accounting for the stellar and dark matter masses we derive gas masses of 0.4-1.2x10^11 Msun. The conversion factor is very high: alpha_CO=3.6+-0.8, consistent with the Galaxy but four times higher than that of local ultra-luminous IR galaxies. The gas accounts for an impressive 50-65% of the baryons within the galaxies' half light radii. We are witnessing truly gas-dominated galaxies at z~1.5, a finding that explains the high specific SFRs observed for z>1 galaxies. The BzK galaxies can be viewed as scaled-up versions of local disk galaxies, with low efficiency star formation taking place inside extended, low excitation gas disks. They are markedly different than local ULIRGs and high-z submm galaxies, which have more excited and compact gas.

We present a universal linear correlation between the stellar mass and surface brightness (SB) of galaxies at 0.3<z<3, using a deep K-band selected catalog in the GOODS-North region. The correlation has a nearly constant slope, independent of redshift and color of galaxies in the rest-z frame. Considering unresolved compact galaxies, the tight correlation gives a lower boundary of SB for a given stellar mass; lower SB galaxies are prohibited over the boundary. The universal slope suggests that the stellar mass in galaxies was build up over their cosmic histories in a similar manner irrelevant to galaxy mass, as oppose to the scenario that massive galaxies mainly accumulated their stellar mass by major merging. In contrast, SB shows a strong dependence on redshift for a given stellar mass. It evolves as (1+z)^(-2.0~-0.8), in addition to dimming as (1+z)^4 by the cosmological expansion effect. The brightening depends on galaxy color and stellar mass. The blue population (rest-frame U-V<0), which is dominated by young and star-forming galaxies, evolves as ~(1+z)^(-0.8 +-0.3) in the rest-V band. On the other hand, the red population (U-V>0) and the massive galaxies (M_*>10^(10)M_sun) shows stronger brightening, (1+z)^(-1.5+-0.1). Based on the comparison with galaxy evolution models, we find that the phenomena are well explained by the pure luminosity evolution of galaxies out to z~3.

We present a kinematic analysis of the globular cluster(GC) system in the giant elliptical galaxy (gE) NGC 4636 in the Virgo cluster. Using the photometric and spectroscopic database of 238 GCs, we have investigated the kinematics of the GC system. The NGC 4636 GC system shows weak overall rotation, which is dominated by the red GCs. However, both the blue GCs and red GCs show some rotation in the inner region at R<4.3'. The velocity dispersion for all the GCs is derived to be sigma_p = 225{+12-9} km/s. The velocity dispersion for the blue GCs (sig=251 km/s) is slightly larger than that for the red GCs (sig=205 km/s). The velocity dispersions for the blue GCs about the mean velocity and about the best fit rotation curve have a significant variation depending on the galactocentric radius. Comparison of observed stellar and GC velocity dispersion profiles with the velocity dispersion profiles calculated from the stellar mass profile shows that the mass-to-light ratio should increase as the galactocentric distance increases, indicating the existence of an extended dark matter halo. From the comparison of the observed GC velocity dispersion profiles and the velocity dispersion profiles calculated for the X-ray mass profiles in the literature, we find that the orbit of the GC system is tangential, and that the orbit of the red GCs is slightly more tangential than that of the blue GCs. We compare the GC kinematics of NGC 4636 with those of other six gEs, finding that the kinematic properties of the GCs are diverse among gEs. We find several correlations between the kinematics of the GCs and the global parameters of their host galaxies. We discuss the implication of the results for the formation models of the GC system in gEs, and suggest a mixture scenario for the origin of the GCs in gEs.

We present a spectroscopic study of the globular clusters (GCs) in the giant elliptical galaxy NGC 4636 in the Virgo cluster. We selected target GC candidates using the Washington photometry derived from the deep CCD images taken at the KPNO 4m. Then we obtained the spectra of 164 target objects in the field of NGC 4636 using the Multi-Object Spectroscopy (MOS) mode of Faint Object Camera and Spectrograph (FOCAS) on the SUBARU 8.2m Telescope. We have measured the velocities for 122 objects: 105 GCs in NGC 4636, the nucleus of NGC 4636, 11 foreground stars, 2 background galaxies, and 3 probable intracluster GCs in the Virgo cluster. The GCs in NGC 4636 are located in the projected galactocentric radius within 10arcmin (corresponding to 43 kpc). The measured velocities for the GCs range from 300km/s to 1600km/s, with a mean value of 932_{-22}^{+25} km/s, which is in good agreement with the velocity for the nucleus of NGC 4636, 928\pm 45 km/s. The velocity dispersion of the GCs in NGC 4636 is derived to be 231_{-17}^{+15} km/s and the velocity dispersion of the blue GCs is slightly larger than that of the red GCs. Combining our results with data in the literature, we produce a master catalog of radial velocities for 238 GCs in NGC 4636. The velocity dispersion of the GCs in the master catalog is found to be 225_{-9}^{+12} km/s for the entire sample, 251_{-12}^{+18} km/s for 108 blue GCs, and 205_{-13}^{+11} km/s for 130 red GCs.

We present the SINS survey with SINFONI of high redshift galaxies. With 80 objects observed and 63 detected, SINS is the largest survey of spatially resolved gas kinematics, morphologies, and physical properties of star-forming galaxies at z~1-3. We describe the selection of the targets, the observations, and the data reduction. We then focus on the "SINS Halpha sample" of 62 rest-UV/optically-selected sources at 1.3<z<2.6 for which we targeted primarily the Halpha and [NII] emission lines. Only 30% of this sample had previous near-IR spectroscopic observations. As a whole, the SINS Halpha sample covers a reasonable representation of massive log(M*/Msun)>~10 star-forming galaxies at z~1.5-2.5, with some bias towards bluer systems compared to pure K-selected samples due to the requirement of secure optical redshift. The sample spans two orders of magnitude in stellar mass and in absolute and specific star formation rates, with median values of approximately log(M*/Msun) = 10.5, 70 Msun/yr, and 3/Gyr. The ionized gas distribution and kinematics are spatially resolved on scales ranging from 1.5 kpc for adaptive optics assisted observations to typically 4-5 kpc for seeing-limited data. The Halpha morphologies tend to be irregular and/or clumpy. About one-third are rotation-dominated yet turbulent disks, another third comprises compact and velocity dispersion-dominated objects, and the remaining galaxies are clear interacting/merging systems; the fraction of rotation-dominated systems increases among the more massive part of the sample. The Halpha luminosities and equivalent widths suggest on average roughly twice higher dust attenuation towards the HII regions relative to the bulk of the stars, and comparable current and past-averaged star formation rates. [Abridged]

We investigate the properties of 12 ultra-massive passively evolving early type galaxies (ETGs) at z_phot>1.4 in the COSMOS 2 deg^2 field. These 12 ETGs were selected as pBzKs, have accurate 1.4<= z_phot <=1.7,high Sersic index profiles typical of ellipticals, no detection at 24 micron, resulting in a complete ETG sample at M*>2.5x10^11 M_sun (Chabrier IMF). Contrary to previous claims, the half light radii estimated in very high S/N imaging data from HST+ACS are found to be large for most of the sample, consistent with local ellipticals. If the high redshift ETGs with M*<2.5x10^11 M_sun are really small in size and compact as reported in previous studies, our result may suggest a "downsizing" scenario, whereby the most massive ETGs reach their final structure earlier and faster than lower mass ones. However, simulating galaxies with morphological properties fixed to those of local ETGs with the same stellar mass show that the few compact galaxies that we still recover in our sample can be understood in term of fluctuations due to noise preventing the recovery of the extended low surface brightness halos in the light profile. Such halos, typical of Sersic profiles, extending even up to 40 kpc, are indeed seen in our sample.

We use very deep near-infrared (NIR) imaging data obtained in MOIRCS Deep Survey (MODS) to investigate the evolution of the galaxy stellar mass function back to z~3. The MODS data reach J=24.2, H=23.1, K=23.1 (5sigma, Vega magnitude) over 103 arcmin^2 (wide) and J=25.1, H=23.7, K=24.1 over 28 arcmin^2 (deep) in the GOODS-North region. The wide and very deep NIR data allow us to measure the number density of galaxies down to low stellar mass (10^9-10^10 Msun) even at high redshift with high statistical accuracy. The normalization of the mass function decreases with redshift and the integrated stellar mass density becomes ~ 8-18% of the local value at z~2 and ~ 4-9% at z~3, which are consistent with results of previous studies in general fields. Furthermore, we found that the low-mass slope becomes steeper with redshift from alpha ~- 1.3 at z~1 to alpha ~- 1.6 at z~3, and that the evolution of the number density of low-mass (10^9-10^10 Msun) galaxies is weaker than that of M* (~10^11 Msun) galaxies. This indicates that the contribution of low-mass galaxies to the total stellar mass density has been significant at high redshift. The steepening of the low-mass slope with redshift is opposite trend expected from the stellar mass dependence of the specific star formation rate reported in previous studies. The present result suggests that the hierarchical merging process overwhelmed the effect of the stellar mass growth by star formation and was very important for the stellar mass assembly of these galaxies at 1<~z<~3.

We present spectroscopy of 15 star-forming BzK galaxies (sBzKs) with K(AB)<23 in the Subaru Deep Field, for which Halpha and some other emission lines are detected in 0.9 to 2.3 micron spectra with a resolution of R=500. Using Halpha luminosities, we obtain star formation rates (SFRs), and then specific SFRs (SSFRs) dividing SFRs by stellar masses, which are derived from SED fitting to BVRi'z'K photometry. It is found that sBzKs with higher stellar masses have larger SFRs. A negative correlation is seen between stellar mass and SSFR, which is consistent with the previous results for z~2 galaxies. This implies that a larger growth of stellar mass occurs in less massive galaxies. In addition, gas-phase oxygen abundances, 12+log(O/H), are derived from the ratio of NII(lambda 6584) to Halpha using the N2 index method. We have found a correlation between stellar mass and oxygen abundance in the sense that more massive sBzKs tend to be more metal rich, which is qualitatively consistent with the relation for UV-selected z~2 galaxies. However, the metallicity of the sBzKs is ~0.2 dex higher than that of UV-selected galaxies with similar stellar masses, which is significant considering the small uncertainties. The sBzKs in our sample have redder R-K colors than the UV-selected galaxies. This galaxy color-dependence in the oxygen abundance may be caused by older or dustier galaxies having higher metallicities at z~2.

We present deep and wide V, I CCD photometry of Ursa Major I (UMa I) dwarf spheroidal galaxy (dSph) in Local Group. The images of the galaxy were taken by Subaru/Suprime-Cam wide field camera, covering a field of 34\arcmin $\times$ 27\arcmin located at the centre of the galaxy. Colour-magnitude diagram (CMD) of the UMa I dSph shows a steep and narrow red giant branch (RGB), blue and red horizontal branch (HB), and main sequence (MS) stars. A well-defined main sequence turn-off (MSTO) is found to be located at V$_{0,MSTO}\sim$23.5 mag. The distance modulus is derived as $(m-M)_0=19.93\pm0.1$ (corresponding to a distance D$=96.8\pm4$ kpc) from the V-band magnitude of the horizontal branch (V$_{0,HB}=20.45\pm0.02$). The mean metallicity of the RGB stars is estimated by the V$-$I colour to [Fe/H]$\sim-2.0$. The turn-off age estimated by overlaying the theoretical isochrones reveals that most of stars in the UMa I dSph are formed at very early epoch ($\sim12$Gyrs ago). The isopleth map of stellar number density of the UMa I dSph, based upon the resolved star counts of MS, RGB, HB stars as well as blue stragglers (BS), shows that the morphology of the UMa I dSph is quite irregular and distorted, suggesting that the galaxy is in a process of disruption. The very old and metal-poor nature of the stellar population implies that the star formation history of this newly discoverd faint dSph may have been different from other well-known `classical' dSphs, which show significant stellar population of intermediate age. The stellar population of the UMa I dSph closely resembles that of Galactic old metal-poor globular cluster, but its size is typical of Galactic dSphs (r$_{e}$=188 [pc], r$_{1/2}$=300 [pc]), and the shape of its spatial density contours suggests that it is undergoing tidal disruption.

We present a simple set of kinematic criteria that can distinguish between galaxies dominated by ordered rotational motion and those involved in major merger events. Our criteria are based on the dynamics of the warm ionized gas (as traced by H-alpha) within galaxies, making this analysis accessible to high-redshift systems, whose kinematics are primarily traceable through emission features. Using the method of kinemetry (developed by Krajnovic and co-workers), we quantify asymmetries in both the velocity and velocity dispersion maps of the warm gas, and the resulting criteria enable us to empirically differentiate between non-merging and merging systems at high redshift. We apply these criteria to 11 of our best-studied rest-frame UV/optical-selected z~2 galaxies for which we have near infrared integral field spectroscopic data from SINFONI on the VLT. Of these 11 systems, we find that >50% have kinematics consistent with a single rotating disk interpretation, while the remaining systems are more likely undergoing major mergers. This result, combined with the short formation timescales of these systems, provides evidence that rapid, smooth accretion of gas plays a significant role in galaxy formation at high redshift.

In this paper and a companion paper, we report on a wide-field imaging survey of the globular cluster (GC) populations around M87 carried out with Suprime-Cam on the 8.2m Subaru telescope. Here we describe the observations, data reduction, and data analysis and present luminosity functions of GC populations around M87 and NGC 4552, another luminous Virgo elliptical in our survey field. The imaging data were taken in the B, V, and I bands with a sky coverage out to 0.5 Mpc from the M87 centre. GC candidates were selected by applying a colour criterion on the B-V and V-I diagram to unresolved objects, which greatly reduces contamination. The data from control fields taken with Subaru/Suprime-Cam were also analyzed for subtraction of contamination in the GC sample. We investigate GC luminosity functions (GCLFs) within 10' (~ 45 kpc) from the host galaxy centre in detail. The V-band turnover magnitude (V_TO) is estimated to be 23.62 +-0.06 mag and 23.56 +- 0.20 mag for the GC population in M87 and NGC 4552, respectively. The GCLF is found to be a function of GC colour; V_TO of the red GC subpopulation (V-I > 1.1) is fainter than that of the blue GC subpopulation (V-I < 1.1) in both M87 and NGC 4552, as expected if the colour differences are primarily due to a metallicity effect and the mass functions of the two subpopulations are similar. The radial dependence of the GCLF is also investigated for the GC population in M87. The GCLF of each subpopulation at 1' < R < 5' is compared with that at 5' < R < 10' but no significant trend with distance is found in the shape of the GCLF. We also estimate GC specific frequencies (S_N) for M87 and NGC 4552. The S_N of the M87 GC population is estimated to be 12.5 +- 0.8 within 25'. The S_N value of the NGC 4552 GC population is estimated to be 5.0 +- 0.6 within 10'.

We have performed a wide-field imaging survey of the globular cluster (GC) populations around M87 out to a radius of ~ 0.5 Mpc with Suprime-Cam on the 8.2m Subaru Telescope. By investigating the GC colour distribution as a function of distance from M87 and NGC 4552, another luminous Virgo elliptical in our survey field, it is found that clear bimodality ((V-I)_peak ~ 1.0 and 1.2) is seen only in the inner regions (< 10 kpc) of the host galaxies and that it becomes weaker with radius due to the decreasing contribution of the red GC (V-I > 1.1) subpopulation. It is also found (both around M87 and NGC 4552) that while the spatial distribution of the red GCs is as centrally concentrated as the host galaxy halo light distribution, the distribution of the blue GCs (V-I < 1.1) tends to be more extended. However, the blue GC distribution around M87 is not as extended as the Virgo cluster mass density profile. Based on these facts, we argue that the contribution of GCs associated with the Virgo cluster (e.g., intergalactic GCs) is not significant around M87 and is not the primary reason for the high S_N value of M87. Instead, most of the blue GCs around luminous ellipticals, as well as the red GCs, are presumed to be associated with the host galaxy. We model the radial profile of GC surface densities out to ~ 0.5 Mpc from M87 by a superposition of the GC populations associated with M87 and with NGC 4552, suggesting the existence of an additional intergalactic GC (i-GC) population independent of the luminous ellipticals. We obtain a crude estimate of i-GC specific frequency S_N = 2.9 (+4.2, -1.5), which may indicate the contribution of tidally captured GCs in the GC population of M87 is low to be consistent with the high S_N value of M87.

We have combined deep BRIz' imaging over 2x940 arcmin^2 fields obtained with the Suprime-Cam on the Subaru telescope with JKs imaging with the SOFI camera at the New Technology Telescope to search for high-redshift massive galaxies. K-band selected galaxies have been identified over an area of ~920 arcmin^2 to K_Vega=19.2, of which 320 arcmin^2 are complete to K_Vega=20. The BzK selection technique was used to obtain complete samples of ~500 candidate massive star-forming galaxies (sBzKs) and ~160 candidate massive, passively-evolving galaxies (pBzKs), both at 1.4<z<2.5. With the (R-K)_Vega > 5 criterion we also identified ~850 extremely red objects (EROs). The surface density of sBzKs and pBzKs is found to 1.20+/-0.05 arcmin^{-2} and 0.38+/-0.03 arcmin^{-2}, respectively. Both sBzKs and pBzKs are strongly clustered, at a level at least comparable to that of EROs, with pBzKs appearing more clustered than sBzKs. We estimate the reddening, star formation rates (SFRs) and stellar masses (M_*) of the sBzKs, confirming that to K_Vega~20 median values are M_*~10^{11}M_sun, SFR 190M_sun yr^{-1}, and E(B-V)~0.44. The most massive sBzKs are also the most actively star-forming, an effect which can be seen as a manifestation of downsizing at early epochs. The space density of massive pBzKs at z~1.4-2 is 20%+/-7% that of similarly massive early-type galaxies at z~0, and similar to that of sBzKs of the same mass. We argue that star formation quenching in these sBzKs will result in nearly doubling the space density of massive early-type galaxies, thus matching their local density.

We present MAMBO 1.2 mm observations of five BzK-pre-selected vigorous starburst galaxies at z~2. Two of these were detected at more than 99.5% confidence levels, with 1.2 mm fluxes around 1.5 mJy. These millimeter fluxes imply vigorous activity with star-formation rates (SFRs) approx. 500-1500 Msun/yr, confirmed also by detections at 24 microns with the MIPS camera on board of the Spitzer satellite. The two detected galaxies are the ones in the sample with the highest SFRs estimated from the rest-frame UV, and their far-IR- and UV-derived SFRs agree reasonably well. This is different from local ULIRGs and high-z submm/mm selected galaxies for which the UV is reported to underestimate SFRs by factors of 10-100, but similar to the average BzK-ULIRG galaxy at z~2. The two galaxies detected at 1.2 mm are brighter in K than the typical NIR-counterparts of MAMBO and SCUBA sources, implying also a significantly different K-band to submm/mm flux ratio. This suggests a scenario in which z~2 galaxies, after their rapid (sub)mm brightest phase opaque to optical/UV light, evolve into a longer lasting phase of K-band bright and massive objects. Targeting the most UV active BzKs could yield substantial detection rates at submm/mm wavelengths.