We report a discovery of $z = 4.72$ radio galaxy, HSC J083913.17+011308.1, by using the Lyman break technique with the Hyper Suprime-Cam Subaru Strategic Survey (HSC-SSP) catalog for VLA FIRST radio sources. The number of known high-$z$ radio galaxies (HzRGs) at $z > 3$ is quite small to constrain the evolution of HzRGs so far. The deep and wide-area optical survey by HSC-SSP enables us to apply the Lyman break technique to a large search for HzRGs. For an HzRG candidate among pre-selected $r$-band dropouts with a radio detection, a follow-up optical spectroscopy with GMOS/Gemini has been performed. The obtained spectrum presents a clear Ly$\alpha$ emission line redshifted to $z=4.72$. The SED fitting analysis with the rest-frame UV and optical photometries suggests the massive nature of this HzRG with $\log{M_*/M_{\odot}} = 11.4$. The small equivalent width of Ly$\alpha$ and the moderately red UV colors indicate its dusty host galaxy, implying a chemically evolved and dusty system. The radio spectral index does not meet a criterion for an ultra-steep spectrum: $\alpha^{325}_{1400}$ of $-1.1$ and $\alpha^{150}_{1400}$ of $-0.9$, demonstrating that the HSC-SSP survey compensates for a sub-population of HzRGs which are missed in surveys focusing on an ultra-steep spectral index.

We present our ALMA Band 8 observations of a damped Ly${\alpha}$ absorption (DLA) system at $z$=3.150 observed in the spectrum of the quasar Q2233+131 at $z$=3.295. The optical counterpart of this DLA has been identified and it shows a double-peaked Ly${\alpha}$ emission line. Since one possible origin of DLAs at high redshift is an outflowing gas from star-forming galaxies, DLA2233+131 provides a good laboratory to investigate the nature of high-$z$ DLAs. Motivated by this, we have carried out ALMA band 8 observations to study the [C II] line in this system. However, we do not detect any significant emission line in the observed pass bands. Instead, we have serendipitously found three submm continuum sources in the observed sky area. One appears to be the quasar Q2233+131 itself while the other two sources are newly identified submm galaxies (SMGs), called SMG1 and SMG2 in this paper. They are located at a separation of 400.7 and 800.1 from Q2233+131, respectively. Their 646 ${\mu}$m fluxes are 6.35 mJy and 6.43 mJy, respectively, being higher than that of Q2233+131, 3.62 mJy. Since these two SMGs are not detected in the optical images obtained with the Hubble Space Telescope and the Subaru Telescope, they have a very red spectral energy distribution. It is, therefore,suggested that they are high-redshift galaxies or very dusty galaxies at intermediate redshift although we cannot rule out the possibility that they are optically very faint SMG analogs at low redshift. Follow up observations will be necessary to explore the nature of this interesting region.

We search for galaxies with a strong Balmer break (Balmer Break Galaxies; BBGs) at $z \sim 6$ over a 0.41 deg$^2$ effective area in the COSMOS field. Based on rich imaging data, including data obtained with the Atacama Large Millimeter/submillimeter Array (ALMA), three candidates are identified by their extremely red $K - [3.6]$ colors as well as by non-detection in X-ray, optical, far-infrared (FIR), and radio bands. The non-detection in the deep ALMA observations suggests that they are not dusty galaxies but BBGs at $z \sim 6$, although contamination from Active Galactic Nuclei (AGNs) at $z \sim 0$ cannot be completely ruled out for the moment. Our spectral energy distribution (SED) analyses reveal that the BBG candidates at $z \sim 6$ have stellar masses of $\approx 5 \times 10^{10} M_{\odot}$ dominated by old stellar populations with ages of $\gtrsim 700$ Myr. Assuming that all the three candidates are real BBGs at $z \sim 6$, we estimate the stellar mass density (SMD) to be $2.4^{+2.3}_{-1.3} \times 10^{4} M_{\odot}$ Mpc$^{-3}$. This is consistent with an extrapolation from the lower redshift measurements. The onset of star formation in the three BBG candidates is expected to be several hundred million years before the observed epoch of $z \sim 6$. We estimate the star-formation rate density (SFRD) contributed by progenitors of the BBGs to be 2.4 -- 12 $\times 10^{-5} M_{\odot}$ yr$^{-1} $Mpc$^{-3}$ at $z > 14$ (99.7\% confidence range). Our result suggests a smooth evolution of the SFRD beyond $z = 8$.

Cosmological simulations predict the Universe contains a network of intergalactic gas filaments, within which galaxies form and evolve. However, the faintness of any emission from these filaments has limited tests of this prediction. We report the detection of rest-frame ultraviolet Lyman-alpha radiation from multiple filaments extending more than one megaparsec between galaxies within the SSA 22 proto-cluster at a redshift of 3.1. Intense star formation and supermassive black-hole activity is occurring within the galaxies embedded in these structures, which are the likely sources of the elevated ionizing radiation powering the observed Lyman-alpha emission. Our observations map the gas in filamentary structures of the type thought to fuel the growth of galaxies and black holes in massive proto-clusters.

Using the HST/ACS $I_{\rm F814W}$-band data, we investigated distribution of apparent axial ratios of $\sim21000$ galaxies with $M_{V}<-20$ at $0.2<z<1.0$ in the COSMOS field as a function of stellar mass, specific star formation rate (sSFR), and redshift. We statistically estimated intrinsic 3-dimensional shapes of these galaxies by fitting the axial-ratio distribution with triaxial ellipsoid models characterized by face-on (middle-to-long) and edge-on (short-to-long) axial ratios $B/A$ and $C/A$. We found that the transition from thin disk to thick spheroid occurs at $\Delta$MS $\sim-1$ dex, i.e., 10 times lower sSFR than that of the main sequence for galaxies with $M_{\rm star} = 10^{10}$--$10^{11} M_{\odot}$ at $0.2<z<1.0$. Furthermore, the intrinsic thickness ($C/A$) of passively evolving galaxies with $M_{\rm star}=10^{10}$--$10^{11}M_{\odot}$ significantly decreases with time from $C/A \sim 0.40$ -- $0.50$ at $z\sim 0.8$ to $C/A\sim0.33$ -- $0.37$ at $z\sim0.4$, while those galaxies with $M_{\rm star}>10^{11}M_{\odot}$ have $C/A\sim0.5$ irrespective of redshift. On the other hand, star-forming galaxies on the main sequence with $10^{9.5}$--$10^{11}M_{\odot}$ show no significant evolution in their shape at $0.2<z<1.0$, but their thickness depends on stellar mass;more massive star-forming galaxies tend to have lower $C/A$ (thinner shape) than low-mass ones. These results suggest that some fraction of star-forming galaxies with a thin disk, which started to appear around $z\sim1$, quench their star formation without violent morphological change, and these newly added quiescent galaxies with a relatively thin shape cause the significant evolution in the axial-ratio distribution of passively evolving galaxies with $M_{\rm star}<10^{11}M_{\odot}$ at $z<1$.

We present a study of the gas kinematics of star-forming galaxies associated with protocluster 4C 23.56 at $z=2.49$ using $0''.4$ resolution CO (4-3) data taken with ALMA. Eleven H$\alpha$ emitters (HAEs) are detected in CO (4-3), including six HAEs that were previously detected in CO (3-2) at a coarser angular resolution. The detections in both CO lines are broadly consistent in the line widths and the redshifts, confirming both detections. With an increase in the number of spectroscopic redshifts, we confirm that the protocluster is composed of two merging groups with a total halo mass of $\log{(M_{\rm cl}/M_{\odot})} =13.4-13.6$, suggesting that the protocluster would evolve into a Virgo-like cluster ($>10^{14} M_{\odot}$). We compare the CO line widths and the CO luminosities with galaxies in other (proto)clusters ($n_{\rm gal}=91$) and general fields ($n_{\rm gal}=80$) from other studies. The 4C23.56 protocluster galaxies have CO line widths and luminosities comparable to other protocluster galaxies on average. On the other hand, the CO line widths are on average broader by $\approx50\%$ compared to field galaxies, while the median CO luminosities are similar. The broader line widths can be attributed to both effects of unresolved gas-rich mergers and/or compact gas distribution, which is supported by our limited but decent angular resolution observations and the size estimate of three galaxies. Based on these results, we argue that gas-rich mergers may play a role in the retention of the specific angular momentum to a value similar to that of field populations during cluster assembly, though we need to verify this with a larger number of samples.

We present physical properties of radio galaxies (RGs) with $f_{\rm 1.4 GHz} >$ 1 mJy discovered by Subaru Hyper Supreme-Cam (HSC) and VLA Faint Images of the Radio Sky at Twenty-Centimeters (FIRST) survey. For 1056 FIRST RGs at $0 < z \leq 1.7$ with HSC counterparts in about 100 deg$^2$, we compiled multi-wavelength data of optical, near-infrared (IR), mid-IR, far-IR, and radio (150 MHz). We derived their color excess ($E (B-V)_{*}$), stellar mass, star formation rate (SFR), IR luminosity, the ratio of IR and radio luminosity ($q_{\rm IR}$), and radio spectral index ($\alpha_{\rm radio}$) that are derived from the SED fitting with CIGALE. We also estimated Eddington ratio based on stellar mass and integration of the best-fit SEDs of AGN component. We found that $E (B-V)_{*}$, SFR, and IR luminosity clearly depend on redshift while stellar mass, $q_{\rm IR}$, and $\alpha_{\rm radio}$ do not significantly depend on redshift. Since optically-faint ($i_{\rm AB} \geq 21.3$) RGs that are newly discovered by our RG survey tend to be high redshift, they tend to not only have a large dust extinction and low stellar mass but also have high SFR and AGN luminosity, high IR luminosity, and high Eddington ratio compared to optically-bright ones. The physical properties of a fraction of RGs in our sample seem to differ from a classical view of RGs with massive stellar mass, low SFR, and low Eddington ratio, demonstrating that our RG survey with HSC and FIRST provides us curious RGs among entire RG population.

We report on the optical properties of infrared (IR)-bright dust-obscured galaxies (DOGs) that are defined as (i-[22])_AB > 7.0. Because supermassive black holes (SMBHs) in IR-bright DOGs are expected to be rapidly growing in the major merger scenario, they provide useful clues for understanding the co-evolution of SMBHs and their host galaxies. However, the optical properties of IR-bright DOGs remain unclear because the optical emission of a DOG is very faint. By combining ~ 105 deg^2 images of the optical, near-IR, and mid-IR data obtained from the Subaru Hyper Suprime-Cam (HSC) survey, the VISTA VIKING survey, and the WISE all-sky survey, respectively, 571 IR-bright DOGs were selected. We found that IR-bright DOGs show a redder (g-z)_AB color than other populations of dusty galaxies, such as ultra-luminous IR galaxies (ULIRGs) at a similar redshift, with a significantly large dispersion. Among the selected DOGs, star-formation (SF) dominated DOGs show a relatively red color, while active galactic nucleus (AGN) dominated DOGs show a rather blue color in optical. This result is consistent with the idea that the relative AGN contribution in the optical emission becomes more significant at a later stage in the major merger scenario. We discovered eight IR-bright DOGs showing a significant blue excess in blue HSC bands (BluDOGs). This blue excess can be interpreted as a leaked AGN emission that is either a directly leaking or a scattered AGN emission, as proposed for some blue-excess Hot DOGs in earlier studies.

We report the result of optical identifications of FIRST radio sources with the Hyper Suprime-Cam Subaru Strategic Program survey (HSC-SSP). The positional cross-match within 1" between the FIRST and HSC-SSP catalogs (i ~< 26) produced more than 3600 optical counterparts in the 156 deg^2 of the HSC-SSP field. The matched counterparts account for more than 50% of the FIRST sources in the search field, which substantially exceed previously reported fractions of SDSS counterparts (i ~< 22) of ~30%. Among the matched sample, 9% are optically unresolved sources such as radio-loud quasars. The optically faint (i > 21) radio galaxies (RGs) show that the fitting linear function of the 1.4 GHz source counts has a slope that is flatter than that of the bright RGs, while optically faint radio quasars show a slope steeper than that of bright radio quasars. The optically faint RGs show a flat slope in the i-band number counts down to 24 mag, implying either less massive or distant radio-active galactic nuclei (AGNs) beyond 24 mag. The photometric redshift and the comparison of colors with the galaxy models show that most of the matched RGs are distributed at redshifts from 0 to 1.5. The optically faint sample includes the high radio-loudness sources that are not seen in the optically bright sample. Such sources are located at redshift z > 1. This study gives ~1500 radio AGNs lying at the optically faint end and high-redshift regime not probed by previous searches.

We present the Ly$\alpha$ luminosity function (LF) derived from 34 Ly$\alpha$ emitters (LAEs) at $z=7.0$ on the sky of $3.1$ deg$^2$, the largest sample compared to those in the literature obtained at a redshift $z\gtrsim7$. The LAE sample is made by deep large-area Subaru narrowband observations conducted by the Cosmic HydrOgen Reionization Unveiled with Subaru (CHORUS) project. The $z=7.0$ Ly$\alpha$ LF of our project is consistent with those of the previous DECam and Subaru studies at the bright and faint ends, respectively, while our $z=7.0$ Ly$\alpha$ LF has uncertainties significantly smaller than those of the previous study results. Exploiting the small errors of our measurements, we investigate the shape of the faint to bright-end Ly$\alpha$ LF. We find that the $z=7.0$ Ly$\alpha$ LF shape can be explained by the steep slope of $\alpha \simeq -2.5$ suggested at $z=6.6$, and that there is no clear signature of a bright-end excess at $z\simeq 7$ claimed by the previous work, which was thought to be made by the ionized bubbles around bright LAEs whose Ly$\alpha$ photons could easily escape from the partly neutral IGM at $z \simeq 7$. We estimate the Ly$\alpha$ luminosity densities (LDs) with Ly$\alpha$ LFs at $z\simeq 6-8$ given by our and the previous studies, and compare the evolution of the UV-continuum LD estimated with dropouts. The Ly$\alpha$ LD monotonically decreases from $z\sim 6$ to $8$, and evolves stronger than the UV-continuum LD, indicative of the Ly$\alpha$ damping wing absorption of the IGM towards the heart of the reionization epoch.

We present the Ly$\alpha$ luminosity functions (LFs) at $z=$5.7 and 6.6 derived from a new large sample of 1,266 Ly$\alpha$ emitters (LAEs) identified in total areas of 14 and 21 deg$^2$, respectively, based on the early narrowband data of the Subaru/Hyper Suprime-Cam (HSC) survey. Together with careful Monte-Carlo simulations that account for the incompleteness of the LAE selection and the flux estimate systematics in the narrowband imaging, we have determined the Ly$\alpha$ LFs with the unprecedentedly small statistical and systematic uncertainties in a wide Ly$\alpha$ luminosity range of $10^{42.8-43.8}$ erg s$^{-1}$. We obtain the best-fit Schechter parameters of $L^{*}_{{\rm Lya}}=1.6^{+2.2}_{-0.6} (1.7^{+0.3}_{-0.7}) \times10^{43}$ erg s$^{-1}$, $\phi^{*}_{{\rm Lya}}=0.85^{+1.87}_{-0.77}\ (0.47^{+1.44}_{-0.44})\times10^{-4}$ Mpc$^{-3}$, and $\alpha=-2.6^{+0.6}_{-0.4}\ (-2.5^{+0.5}_{-0.5})$ at $z=5.7$ ($6.6$). We confirm that our best-estimate Ly$\alpha$ LFs are consistent with the majority of the previous studies, but find that our Ly$\alpha$ LFs do not agree with the high number densities of LAEs recently claimed by Matthee/Santos et al.'s studies that may overcorrect the incompleteness and the flux systematics. Our Ly$\alpha$ LFs at $z=5.7$ and $6.6$ show an indication that the faint-end slope is very steep ($\alpha \simeq -2.5$), although it is also possible that the bright-end LF results are enhanced by systematic effects such as the contribution from AGNs, blended merging galaxies, and/or large ionized bubbles around bright LAEs. Comparing our Ly$\alpha$ LF measurements with four independent reionization models, we estimate the neutral hydrogen fraction of the IGM to be $x_{\rm HI}=0.3\pm0.2$ at $z=$6.6 that is consistent with the small Thomson scattering optical depth obtained by Planck 2016.

We present an unprecedentedly large catalog consisting of 2,230 > L* Lya emitters (LAEs) at z=5.7 and 6.6 on the 13.8 and 21.2 deg2 sky, respectively, that are identified by the SILVERRUSH program with the first narrowband imaging data of the Hyper Suprime-Cam (HSC) survey. We confirm that the LAE catalog is reliable on the basis of 96 LAEs whose spectroscopic redshifts are already determined by this program and the previous studies. This catalogue is also available on-line. Based on this catalogue, we derive the rest-frame Lya equivalent-width distributions of LAEs at z~5.7-6.6 that are reasonably explained by the exponential profiles with the scale lengths of ~120-170A, showing no significant evolution from z~5.7 to z~6.6. We find that 275 LAEs with a large equivalent width (LEW) of >240A are candidates of young-metal poor galaxies and AGNs. We also find that the fraction of LEW LAEs to all ones is 4% and 21% at z~5.7 and z~6.6, respectively. Our LAE catalog includes 11 Lya blobs (LABs) that are LAEs with spatially extended Lya emission whose profile is clearly distinguished from those of stellar objects at the >~ 3sigma level. The number density of the LABs at z=6-7 is ~ 10^-7-10^-6 Mpc^-3, being ~ 10-100 times lower than those claimed for LABs at z~ 2-3, suggestive of disappearing LABs at z>~6, albeit with the different selection methods and criteria for the low and high-z LABs.

We present Lya and UV-nebular emission line properties of bright Lya emitters (LAEs) at z=6-7 with a luminosity of log L_Lya/[erg s-1] = 43-44 identified in the 21-deg2 area of the SILVERRUSH early sample developed with the Subaru Hyper Suprime-Cam (HSC) survey data. Our optical spectroscopy newly confirm 21 bright LAEs with clear Lya emission, and contribute to make a spectroscopic sample of 96 LAEs at z=6-7 in SILVERRUSH. From the spectroscopic sample, we select 7 remarkable LAEs as bright as Himiko and CR7 objects, and perform deep Keck/MOSFIRE and Subaru/nuMOIRCS near-infrared spectroscopy reaching the 3sigma-flux limit of ~ 2x10^{-18} erg s-1 for the UV-nebular emission lines of He II1640, C IV1548,1550, and O III]1661,1666. Except for one tentative detection of C IV, we find no strong UV-nebular lines down to the flux limit, placing the upper limits of the rest-frame equivalent widths (EW_0) of ~2-4 A for He II, C IV, and O III] lines. Here we also investigate the VLT/X-SHOOTER spectrum of CR7 whose 6 sigma detection of He II is claimed by Sobral et al. Although two individuals and the ESO-archive service carefully re-analyze the X-SHOOTER data that are used in the study of Sobral et al., no He II signal of CR7 is detected, supportive of weak UV-nebular lines of the bright LAEs even for CR7. Spectral properties of these bright LAEs are thus clearly different from those of faint dropouts at z~7 that have strong UV-nebular lines shown in the various studies. Comparing these bright LAEs and the faint dropouts, we find anti-correlations between the UV-nebular line EW_0 and UV-continuum luminosity, which are similar to those found at z~2-3.

We present the result of our spectroscopic follow-up observation for faint quasar candidates at z~5 in a part of the Canada-France-Hawaii Telescope Legacy Survey wide field. We select nine photometric candidates and identify three z~5 faint quasars, one z~4 faint quasar, and a late-type star. Since two faint quasar spectra show Civ emission line without suffering from a heavy atmospheric absorption, we estimate the black hole mass (M$_{BH}$) and Eddington ratio (L/L$_{Edd}$) of them. The inferred log M$_{BH}$ are 9.04+/-0.14 and 8.53+/-0.20, respectively. In addition, the inferred log (L/L$_{Edd}$) are -1.00+/-0.15 and -0.42+/-0.22, respectively. If we adopt that L/L$_{Edd}$= constant or $\propto$ (1+z)^2, the seed black hole masses (M$_{seed}$) of our z~5 faint quasars are expected to be >10^5 M$_\odot$ in most cases. We also compare the observational results with a mass accretion model where angular momentum is lost due to supernova explosions (Kawakatu & Wada 2008). Accordingly, M$_{BH}$ of the z~5 faint quasars in our sample can be explained even if M$_{seed}$ is ~10^3M$_\odot$. Since z~6 luminous qusars and our z~5 faint quasars are not on the same evolutionary track, z~6 luminous quasars and our z~5 quasars are not the same populations but different populations, due to the difference of a period of the mass supply from host galaxies. Furthermore, we confirm that one can explain M$_{BH}$ of z~6 luminous quasars and our z~5 faint quasars even if their seed black holes of them are formed at z~7.

We present the SILVERRUSH program strategy and clustering properties investigated with $\sim 2,000$ Ly$\alpha$ emitters at $z=5.7$ and $6.6$ found in the early data of the Hyper Suprime-Cam (HSC) Subaru Strategic Program survey exploiting the carefully designed narrowband filters. We derive angular correlation functions with the unprecedentedly large samples of LAEs at $z=6-7$ over the large total area of $14-21$ deg$^2$ corresponding to $0.3-0.5$ comoving Gpc$^2$. We obtain the average large-scale bias values of $b_{\rm avg}=4.1\pm 0.2$ ($4.5\pm 0.6$) at $z=5.7$ ($z=6.6$) for $\gtrsim L^*$ LAEs, indicating the weak evolution of LAE clustering from $z=5.7$ to $6.6$. We compare the LAE clustering results with two independent theoretical models that suggest an increase of an LAE clustering signal by the patchy ionized bubbles at the epoch of reionization (EoR), and estimate the neutral hydrogen fraction to be $x_{\rm HI}=0.15^{+0.15}_{-0.15}$ at $z=6.6$. Based on the halo occupation distribution models, we find that the $\gtrsim L^*$ LAEs are hosted by the dark-matter halos with the average mass of $\log (\left < M_{\rm h} \right >/M_\odot) =11.1^{+0.2}_{-0.4}$ ($10.8^{+0.3}_{-0.5}$) at $z=5.7$ ($6.6$) with a Ly$\alpha$ duty cycle of 1 % or less, where the results of $z=6.6$ LAEs may be slightly biased, due to the increase of the clustering signal at the EoR. Our clustering analysis reveals the low-mass nature of $\gtrsim L^*$ LAEs at $z=6-7$, and that these LAEs probably evolve into massive super-$L^*$ galaxies in the present-day universe.

We investigate gas contents of star-forming galaxies associated with protocluster 4C23.56 at z = 2.49 by using the redshifted CO(3-2) and 1.1 mm dust continuum with the Atacama Large Millimeter/submillimeter Array. The observations unveil seven CO detections out of 22 targeted H$\alpha$ emitters (HAEs) and four out of 19 in 1.1 mm dust continuum. They have high stellar mass ($M_{\star}>4\times 10^{10}$ $M_{\odot}$) and exhibit a specific star-formation rate typical of main-sequence star forming galaxies at $z\sim2.5$. Different gas mass estimators from CO(3-2) and 1.1 mm yield consistent values for simultaneous detections. The gas mass ($M_{\rm gas}$) and gas fraction ($f_{\rm gas}$) are comparable to those of field galaxies, with $M_{\rm gas}=[0.3, 1.8]\times10^{11} \times (\alpha_{\rm CO}/(4.36\times A(Z)$)) M$_{\odot}$, where $\alpha_{\rm CO}$ is the CO-to-H$_2$ conversion factor and $A(Z)$ the additional correction factor for the metallicity dependence of $\alpha_{\rm CO}$, and $\langle f_{\rm gas}\rangle = 0.53 \pm 0.07$ from CO(3-2). Our measurements place a constraint on the cosmic gas density of high-$z$ protoclusters, indicating the protocluster is characterized by a gas density higher than that of the general fields by an order of magnitude. We found $\rho (H_2)\sim 5 \times 10^9 \,M_{\odot}\,{\rm Mpc^{-3}}$ with the CO(3-2) detections. The five ALMA CO detections occur in the region of highest galaxy surface density, where the density positively correlates with global star-forming efficiency (SFE) and stellar mass. Such correlations imply a potentially critical role of environment on early galaxy evolution at high-z protoclusters, although future observations are necessary for confirmation.

We present the near-infrared high resolution imaging of an extremely dense group of galaxies at the core of the protocluster at $z=3.09$ in the SSA22 field by using the adaptive optics AO188 and the Infrared Camera and Spectrograph (IRCS) on Subaru Telescope. Wide morphological variety of them suggests their on-going dramatic evolutions. One of the two quiescent galaxies (QGs), the most massive one in the group, is a compact elliptical with an effective radius $r_{e} = 1.37\pm0.75$ kpc. It supports the two-phase formation scenario of giant ellipticals today that a massive compact elliptical is formed at once and evolves in the size and stellar mass by series of mergers. Since this object is a plausible progenitor of a brightest cluster galaxy (BCG) of one of the most massive clusters today, it requires strong size ($\ga10$) and stellar mass ($\sim$ four times by $z=0$) growths. Another QG hosts an AGN(s) and is fitted with a model composed from an nuclear component and S\'ersic model. It shows spatially extended [O{\footnotesize III}]$\lambda$5007 emission line compared to the continuum emission, a plausible evidence of outflows. Massive star forming galaxies (SFGs) in the group are two to three times larger than the field SFGs at similar redshift. Although we obtained the $K$-band image deeper than the previous one, we found no candidate new members. This implies a physical deficiency of low mass galaxies with stellar mass $M_{\star}\la4\times10^{10}~M_{\odot}$ and/or poor detection completeness of them owing to their diffuse morphologies.

We study the properties of dark matter halos that contain star-forming galaxies at $1.43 \le z \le 1.74$ using the FMOS-COSMOS survey. The sample consists of 516 objects with a detection of the H$\alpha$ emission line, that represent the star-forming population at this epoch having a stellar mass range of $10^{9.57}\le M_\ast/M_\odot \lesssim 10^{11.4}$ and a star formation rate range of $15\lesssim \mathrm{SFR}/(M_\odot \mathrm{yr^{-1}}) \lesssim 600$. We measure the projected two-point correlation function while carefully taking into account observational biases, and find a significant clustering amplitude at scales of $0.04$-$10~h^{-1}~\mathrm{cMpc}$, with a correlation length $r_0 = 5.21^{+0.70}_{-0.67}~h^{-1}~\mathrm{cMpc}$ and a bias $b=2.59^{+0.41}_{-0.34}$. We interpret our clustering measurement using a halo occupation distribution model. The sample galaxies appear to reside in halos with mass $M_\mathrm{h} = 4.6^{+1.1}_{-1.6}\times10^{12}~h^{-1}M_\odot$ on average that will likely become present-day halos of mass $M_\mathrm{h} (z=0) \sim2\times10^{13}~h^{-1}M_\odot$, equivalent to the typical halo mass scale of galaxy groups. We then confirm the decline of the stellar-to-halo mass ratio at $M_\mathrm{h}<10^{12}~M_\odot$, finding $M_\ast/M_\mathrm{h} \approx 5\times10^{-3}$ at $M_\mathrm{h}=10^{11.86}~M_\odot$, which is lower by a factor of 2-4 than those measured at higher masses. Finally, we use our results to illustrate the future capabilities of Subaru's Prime-Focus Spectrograph, a next-generation instrument that will provide strong constraints on the galaxy-formation scenario by obtaining precise measurements of galaxy clustering at $z>1$.

In order to investigate the physical relationship between strong Ly-alpha absorbers (logN_HI>20.0 cm^-2) such as damped Ly-alpha absorption systems (DLAs) and young star-forming galaxies at high redshift, we have conducted narrow-band observations of Ly-alpha emitters (LAEs) in a concentrated region of strong Ly-alpha absorbers at z= 2.3, the J1230+34 field. Using a catalog of Ly-alpha absorbers with logN_HI>20.0 cm^-2 based on the baryon oscillation spectroscopic survey (BOSS), we found 6 fields where 3 or more absorbers are concentrated within a (50 Mpc)^3 cubic box in the comoving scale. Among them, we focus on the J1230+34 field, where 2 DLAs and 2 sub-DLAs present. Our narrow-band imaging observations with Subaru/Suprime-Cam using a custom-made filter, NB400 (lambda_c=4003 A and FWHM=92 A) yield a sample of 149 LAEs in this field. In the large scale (~50 Mpc), we have found no differences between the obtained Ly-alpha luminosity function and those in the blank fields at similar redshifts. We also compare the frequency distribution of the Ly-alpha rest-frame equivalent width (EW_0) in the target field and other fields including both overdensity region and blank field, but find no differences. On the other hand, in the small scale (~10 Mpc), we have found a possible overdensity of LAEs around a DLA with the highest HI column density (N_HI=21.08 cm^-2) in the target field while there are no density excess around the other absorbers with a lower N_HI.

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 present measurements of the clustering properties of a sample of infrared (IR) bright dust-obscured galaxies (DOGs). Combining 125 deg$^2$ of wide and deep optical images obtained with the Hyper Suprime-Cam on the Subaru Telescope and all-sky mid-IR (MIR) images taken with Wide-Field Infrared Survey Explorer, we have discovered 4,367 IR-bright DOGs with $(i - [22])_{\rm AB}$ $>$ 7.0 and flux density at 22 $\mu$m $>$ 1.0 mJy. We calculate the angular autocorrelation function (ACF) for a uniform subsample of 1411 DOGs with 3.0 mJy $<$ flux (22 $mu$m) $<$ 5.0 mJy and $i_{\rm AB}$ $<$ 24.0. The ACF of our DOG subsample is well-fit with a single power-law, $\omega (\theta)$ = (0.010 $\pm$ 0.003) $\theta^{-0.9}$, where $\theta$ in degrees. The correlation amplitude of IR-bright DOGs is larger than that of IR-faint DOGs, which reflects a flux-dependence of the DOG clustering, as suggested by Brodwin et al. (2008). We assume that the redshift distribution for our DOG sample is Gaussian, and consider 2 cases: (1) the redshift distribution is the same as IR-faint DOGs with flux at 22 $\mu$m $<$ 1.0 mJy, mean and sigma $z$ = 1.99 $\pm$ 0.45, and (2) $z$ = 1.19 $\pm$ 0.30, as inferred from their photometric redshifts. The inferred correlation length of IR-bright DOGs is $r_0$ = 12.0 $\pm$ 2.0 and 10.3 $\pm$ 1.7 $h^{-1}$ Mpc, respectively. IR-bright DOGs reside in massive dark matter halos with a mass of $\log [\langle M_{\mathrm{h}} \rangle / (h^{-1} M_{\odot})]$ = 13.57$_{-0.55}^{+0.50}$ and 13.65$_{-0.52}^{+0.45}$ in the two cases, respectively.

We study the extent to which very bright (-23.0 < MUV < -21.75) Lyman-break selected galaxies at redshifts z~7 display detectable Lya emission. To explore this issue, we have obtained follow-up optical spectroscopy of 9 z~7 galaxies from a parent sample of 24 z~7 galaxy candidates selected from the 1.65 sq.deg COSMOS-UltraVISTA and SXDS-UDS survey fields using the latest near-infrared public survey data, and new ultra-deep Subaru z'-band imaging (which we also present and describe in this paper). Our spectroscopy has yielded only one possible detection of Lya at z=7.168 with a rest-frame equivalent width EW_0 = 3.7 (+1.7/-1.1) Angstrom. The relative weakness of this line, combined with our failure to detect Lya emission from the other spectroscopic targets allows us to place a new upper limit on the prevalence of strong Lya emission at these redshifts. For conservative calculation and to facilitate comparison with previous studies at lower redshifts, we derive a 1-sigma upper limit on the fraction of UV bright galaxies at z~7 that display EW_0 > 50 Angstrom, which we estimate to be < 0.23. This result may indicate a weak trend where the fraction of strong Lya emitters ceases to rise, and possibly falls between z~6 and z~7. Our results also leave open the possibility that strong Lya may still be more prevalent in the brightest galaxies in the reionization era than their fainter counterparts. A larger spectroscopic sample of galaxies is required to derive a more reliable constraint on the neutral hydrogen fraction at z~7 based on the Lya fraction in the bright galaxies.

We investigate morphological properties of 61 Lyman-alpha emitters (LAEs) at z = 4.86 identified in the COSMOS field, based on Hubble Space Telescope Advanced Camera for Surveys (ACS) imaging data in the F814W-band. Out of the 61 LAEs, we find the ACS counterparts for the 54 LAEs. Eight LAEs show double-component structures with a mean projected separation of 0."63 (~ 4.0 kpc at z = 4.86). Considering the faintness of these ACS sources, we carefully evaluate their morphological properties, that is, size and ellipticity. While some of them are compact and indistinguishable from the PSF half-light radius of 0."07 (~ 0.45 kpc), the others are clearly larger than the PSF size and spatially extended up to 0."3 (~ 1.9 kpc). We find that the ACS sources show a positive correlation between ellipticity and size and that the ACS sources with large size and round shape are absent. Our Monte Carlo simulation suggests that the correlation can be explained by (1) the deformation effects via PSF broadening and shot noise or (2) the source blending in which two or more sources with small separation are blended in our ACS image and detected as a single elongated source. Therefore, the 46 single-component LAEs could contain the sources which consist of double (or multiple) components with small spatial separation (i.e., < 0."3 or 1.9 kpc). Further observation with high angular resolution at longer wavelengths (e.g., rest-frame wavelengths of > 4000 A) is inevitable to decipher which interpretation is adequate for our LAE sample.

We report the discovery of an extremely dense group of massive galaxies at the centre of the protocluster at $z=3.09$ in the SSA22 field from near-infrared spectroscopy conducted with the Multi-Object InfraRed Camera and Spectrograph (MOIRCS) equipped on the Subaru Telecope. The newly discovered group comprises seven galaxies confirmed at $z_{\rm spec}\approx3.09$ within 180 kpc including five massive objects with the stellar masses larger than $10^{10.5}~M_{\odot}$ and is associated with a bright sub-mm source SSA22-AzTEC14. The dynamical mass of the group estimated from the line-of-sight velocity dispersion of the members is $M_{\rm dyn}\sim1.6\pm0.3\times10^{13}~M_{\odot}$. Such a dense group is expected to be very rare at high redshift as we found only a few comparable systems in large-volume cosmological simulations. Such rare groups in the simulations are hosted in collapsed halos with $M_{\rm vir}=10^{13.4}-10^{14.0}~M_{\odot}$ and evolve into the brightest cluster galaxies (BCGs) of the most massive clusters at present. The observed AzTEC14 group at $z=3.09$ is therefore very likely to be a proto-BCG in the multiple merger phase. The observed total stellar mass of the group is $5.8^{+5.1}_{-2.0}\times10^{11}~M_{\odot}$. It suggests that over half the stellar mass of its descendant had been formed by $z=3$. Moreover, we identified over two members for each of the four Ly$\alpha$ blobs (LABs) using our new spectroscopic data. This verifies our previous argument that many of the LABs in the SSA22 protocluster associated with multiple developed stellar components.

We present both the observations and the data reduction procedures of the Subaru COSMOS 20 project that is an optical imaging survey of the HST COSMOS field, carried out by using Suprime-Cam on the Subaru Telescope with the following 20 optical filters: 6 broad-band (B, g', V, r', i', and z'), 2 narrow-band (NB711 and NB816), and 12 intermediate-band filters (IA427, IA464, IA484, IA505, IA527, IA574, IA624, IA679, IA709, IA738, IA767, and IA827). A part of this project is described in Taniguchi et al. (2007) and Capak et al. (2007) for the six broad-band and one narrow-band (NB816) filter data. In this paper, we present details of the observations and data reduction for remaining 13 filters (the 12 IA filters and NB711). In particular, we describe the accuracy of both photometry and astrometry in all the filter bands. We also present optical properties of the Suprime-Cam IA filter system in Appendix.

In order to investigate the origin of quasars, we estimate the bias factor for low-luminosity quasars at high redshift for the first time. In this study, we use the two-point angular cross-correlation function (CCF) for both low-luminosity quasars at $-24<M_{\rm 1450}<-22$ and Lyman-break galaxies (LBGs). Our sample consists of both 25 low-luminosity quasars (16 objects are spectroscopically confirmed low-luminosity quasars) in the redshift range $3.1<z<4.5$ and 835 color-selected LBGs with $z^{\prime}_{\rm LBG}<25.0$ at $z\sim4$ in the COSMOS field. We have made our analysis for the following two quasar samples; (1) the spectroscopic sample (the 16 quasars confirmed by spectroscopy), and (2) the total sample (the 25 quasars including 9 quasars with photometric redshifts). The bias factor for low-luminosity quasars at $z\sim4$ is derived by utilizing the quasar-LBG CCF and the LBG auto-correlation function. We then obtain the $86\%$ upper limits of the bias factors for low-luminosity quasars, that are 5.63 and 10.50 for the total and the spectroscopic samples, respectively. These bias factors correspond to the typical dark matter halo masses, log $(M_{\rm DM}/(h^{-1}M_{\odot}))=$$12.7$ and $13.5$, respectively. This result is not inconsistent with the predicted bias for quasars which is estimated by the major merger models.

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 report a discovery of 6 massive galaxies with both extremely large Lya equivalent width and evolved stellar population at z ~ 3. These MAssive Extremely STrong Lya emitting Objects (MAESTLOs) have been discovered in our large-volume systematic survey for strong Lya emitters (LAEs) with twelve optical intermediate-band data taken with Subaru/Suprime-Cam in the COSMOS field. Based on the SED fitting analysis for these LAEs, it is found that these MAESTLOs have (1) large rest-frame equivalent width of EW_0(Lya) ~ 100--300 A, (2) M_star ~ 10^10.5--10^11.1 M_sun, and (3) relatively low specific star formation rates of SFR/M_star ~ 0.03--1 Gyr^-1. Three of the 6 MAESTLOs have extended Ly$\alpha$ emission with a radius of several kpc although they show very compact morphology in the HST/ACS images, which correspond to the rest-frame UV continuum. Since the MAESTLOs do not show any evidence for AGNs, the observed extended Lya emission is likely to be caused by star formation process including the superwind activity. We suggest that this new class of LAEs, MAESTLOs, provides a missing link from star-forming to passively evolving galaxies at the peak era of the cosmic star-formation history.

The relation between the stellar mass and the star formation rate characterizes how the instantaneous star formation is determined by the galaxy past star formation history and by the growth of the dark matter structures. We deconstruct the M-SFR plane by measuring the specific SFR functions in several stellar mass bins from z=0.2 out to z=1.4. Our analysis is primary based on a MIPS 24$\mu m$ selected catalogue combining the COSMOS and GOODS surveys. We estimate the SFR by combining mid- and far-infrared data for 20500 galaxies. The sSFR functions are derived in four stellar mass bins within the range 9.5<log(M/Msun)<11.5. First, we demonstrate the importance of taking into account selection effects when studying the M-SFR relation. Secondly, we find a mass-dependent evolution of the median sSFR with redshift varying as $sSFR \propto (1+z)^{b}$, with $b$ increasing from $b=2.88$ to $b=3.78$ between $M=10^{9.75}Msun$ and $M=10^{11.1}Msun$, respectively. At low masses, this evolution is consistent with the cosmological accretion rate and predictions from semi-analytical models (SAM). This agreement breaks down for more massive galaxies showing the need for a more comprehensive description of the star-formation history in massive galaxies. Third, we obtain that the shape of the sSFR function is invariant with time at z<1.4 but depends on the mass. We observe a broadening of the sSFR function ranging from 0.28 dex at $M=10^{9.75}Msun$ to 0.46 dex at $M=10^{11.1}Msun$. Such increase in the scatter of the M-SFR relation suggests an increasing diversity of SFHs as the stellar mass increases. Finally, we find a gradual decline of the sSFR with mass as $log(sSFR) \propto -0.17M$. We discuss the numerous physical processes, as gas exhaustion in hot gas halos or secular evolution, which can gradually reduce the sSFR and increase the SFH diversity.

Using the HST/WFC3 and ACS multi-band imaging data taken in CANDELS and 3D-HST, we study the general properties and the diversity of the progenitors of the Milky Way (MWs) and local massive galaxy (MGs) at 0.5 < z < 3.0, based on a constant cumulative number density analysis. After careful data reduction and stacking analysis, we conduct a radially resolved pixel SED fitting to obtain the radial distributions of the stellar mass and rest-frame colors. The stellar mass of MWs increases in self-similar way, irrespective of the radial distance, while that of MGs grows in inside-out way where they obtain ~ 75% of the total mass at outer (> 2.5 kpc) radius since z ~ 2. Although the radial mass profiles evolve in distinct ways, the formation and quenching of the central dense region (or bulge) ahead of the outer disk formation are found to be common for both systems. The sudden reddening of bulge at z ~ 1.6 and z ~ 2.4 for MWs and MGs, respectively, suggests the formation of bulge and would give a clue to the different gas accretion histories and quenching. A new approach to evaluate the morphological diversity is conducted by using the average surface density profile and its dispersion. The variety of the radial mass profiles for MGs peaks at higher redshift (z > 2.8), and then rapidly converges to more uniform shape at z < 1.5, while that for MWs remains in the outer region over the redshift. Compared with the observed star formation rates and color profiles, the evolution of variety is consistently explained by the star formation activities.

Using the Hubble Space Telescope/Wide Field Camera 3 imaging data and multi-wavelength photometric catalog, we investigated the dust temperature of passively evolving and star-forming galaxies at 0.2<z<1.0 in the CANDELS fields. We estimated the stellar radiation field by low-mass stars from the stellar mass and surface brightness profile of these galaxies and then calculated their steady-state dust temperature. At first, we tested our method using nearby early-type galaxies with the deep FIR data by the Herschel Virgo cluster survey and confirmed that the estimated dust temperatures are consistent with the observed temperatures within the uncertainty. We then applied the method to galaxies at 0.2<z<1.0, and found that most of passively evolving galaxies with Mstar > 10^{10} Msun have a relatively high dust temperature of Tdust > 20 K, for which the formation efficiency of molecular hydrogen on the surface of dust grains in the diffuse ISM is expected to be very low from the laboratory experiments. The fraction of passively evolving galaxies strongly depends on the expected dust temperature at all redshifts and increases rapidly with increasing the temperature around Tdust ~ 20 K. These results suggest that the dust heating by low-mass stars in massive galaxies plays an important role for the continuation of their passive evolution, because the lack of the shielding effect of the molecular hydrogen on the UV radiation can prevent the gas cooling and formation of new stars.

We present the results of near-infrared spectroscopic observations of the $K$-band selected candidate galaxies in the protocluster at $z=3.09$ in the SSA22 field. We observed 67 candidates with $K_{\rm AB}<24$ and confirmed redshifts of the 39 galaxies at $2.0< z_{\rm spec}< 3.4$. Of the 67 candidates, 24 are certainly protocluster members with $3.04\leq z_{\rm spec}\leq 3.12$, which are massive red galaxies those have been unidentified in previous optical observations of the SSA22 protocluster. Many distant red galaxies (DRGs; $J-K_{\rm AB}>1.4$), hyper extremely red objects (HEROs; $J-K_{\rm AB}>2.1$), {\it Spitzer} MIPS 24 $\mu$m sources, active galactic nuclei (AGNs) as well as the counterparts of Ly$\alpha$ blobs and the AzTEC/ASTE 1.1-mm sources in the SSA22 field are also found to be the protocluster members. The mass of the SSA22 protocluster is estimated to be $\sim2-5\times10^{14}~M_{\odot}$ and this system is plausibly a progenitor of the most massive clusters of galaxies in the current Universe. The reddest ($J-K_{\rm AB}\geq 2.4$) protocluster galaxies are massive galaxies with $M_{\rm star}\sim10^{11}~M_{\odot}$ showing quiescent star formation activities and plausibly dominated by old stellar populations. Most of these massive quiescent galaxies host moderately luminous AGNs detected by X-ray. There are no significant differences in the [O{\footnotesize III}] $\lambda$5007/H$\beta$ emission line ratios, and [O{\footnotesize III}] $\lambda$5007 line widths and spatial extents of the protocluster galaxies from those of massive galaxies at $z\sim2-3$ in the general field.

We compared the number of faint stars detected in deep survey fields with the current stellar distribution model of the Galaxy and found that the detected number in the H band is significantly smaller than the predicted number. This indicates that M-dwarfs, the major component, are fewer in the halo and the thick disk. We used archived data of several surveys in both the north and south field of GOODS (Great Observatories Origins Deep Survey), MODS in GOODS-N, and ERS and CANDELS in GOODS-S. The number density of M-dwarfs in the halo has to be 20+/-13% relative to that in the solar vicinity, in order for the detected number of stars fainter than 20.5 mag in the H band to match with the predicted value from the model. In the thick disk, the number density of M-dwarfs must be reduced (52+/-13%) or the scale height must be decreased (~600 pc). Alternatively, overall fractions of the halo and thick disks can be significantly reduced to achieve the same effect, because our sample mainly consists of faint M-dwarfs. Our results imply that the M-dwarf population in regions distant from the Galactic plane is significantly smaller than previously thought. We then discussed the implications this has on the suitability of the model predictions for the prediction of non-companion faint stars in direct imaging extrasolar planet surveys by using the best-fit number densities.

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.

Using the Hubble Space Telescope/Advanced Camera for Surveys data in the COSMOS field, we systematically searched clumpy galaxies at 0.2<z<1.0 and investigated the fraction of clumpy galaxies and its evolution as a function of stellar mass, star formation rate (SFR), and specific SFR (SSFR). The fraction of clumpy galaxies in star-forming galaxies with Mstar > 10^9.5 Msun decreases with time from ~0.35 at 0.8<z<1.0 to ~0.05 at 0.2<z<0.4 irrespective of the stellar mass, although the fraction tends to be slightly lower for massive galaxies with Mstar > 10^10.5 Msun at each redshift. On the other hand, the fraction of clumpy galaxies increases with increasing both SFR and SSFR in all the redshift ranges we investigated. In particular, we found that the SSFR dependences of the fractions are similar among galaxies with different stellar masses, and the fraction at a given SSFR does not depend on the stellar mass in each redshift bin. The evolution of the fraction of clumpy galaxies from z~0.9 to z~0.3 seems to be explained by such SSFR dependence of the fraction and the evolution of SSFRs of star-forming galaxies. The fraction at a given SSFR also appears to decrease with time, but this can be due to the effect of the morphological K-correction. We suggest that these results are understood by the gravitational fragmentation model for the formation of giant clumps in disk galaxies, where the gas mass fraction is a crucial parameter.

We present a 4 deg^2 weak gravitational lensing survey of subhalos in the very nearby Coma cluster using the Subaru/Suprime-Cam. The large apparent size of cluster subhalos allows us to measure the mass of 32 subhalos detected in a model-independent manner, down to the order of 10^-3 of the virial mass of the cluster. Weak-lensing mass measurements of these shear-selected subhalos enable us to investigate subhalo properties and the correlation between subhalo masses and galaxy luminosities for the first time. The mean distortion profiles stacked over subhalos show a sharply truncated feature which is well-fitted by a Navarro-Frenk-White (NFW) mass model with the truncation radius, as expected due to tidal destruction by the main cluster. We also found that subhalo masses, truncation radii, and mass-to-light ratios decrease toward the cluster center. The subhalo mass function, dn/dln M_sub, in the range of 2 orders of magnitude in mass, is well described by a single power law or a Schechter function. Best-fit power indices of 1.09_-0.32^+0.42 for the former model and 0.99_-0.23^+0.34 for the latter, are in remarkable agreement with slopes of ~0.9-1.0 predicted by the cold dark matter paradigm. The tangential distortion signals in the radial range of 0.02-2Mpc/h from the cluster center show a complex structure which is well described by a composition of three mass components of subhalos, the NFW mass distribution as a smooth component of the main cluster, and a lensing model from a large scale structure behind the cluster. Although the lensing signals are 1 order of magnitude lower than those for clusters at z~0.2, the total signal-to-noise ratio, S/N=13.3, is comparable to, or higher, because the enormous number of background source galaxies compensates for the low lensing efficiency of the low lensing efficiency of the nearby cluster.

We analyze the recent released HST/WFC3 IR images in the GOODS-N region to study the formation and evolution of Quiescent galaxies (QGs). After examining the reliability with artificial galaxies, we obtain the morphological parameters with S'ersic profile of 299 QGs and 1,083 star-forming galaxies (SFGs) at z ~ 0.5-3.0, finding the evolution of re and n of massive (M* > 10^10.5 Msun) QGs while weaker evolution of SFGs and less massive (M* < 10^10.5 Msun) QGs. The regression of the size evolution of massive QGs follows re \propto (1 + z)-{\alpha}re with {\alpha}re = 1.06 \pm 0.19 (a factor of ~ 2.2 increase from z ~ 2.5 to ~ 0.5), which is consistent with the general picture of the significant size growth. For the further understanding of the evolution scenario, we study the evolution of S'ersic index, n, and find that of massive QGs to significantly evolve as n \propto (1 + z)-{\alpha}n with {\alpha}n = 0.74 \pm 0.23 (n ~ 1 at z ~ 2.5 to n ~ 4 at z ~ 0.5), while those of the other populations are unchanged (n ~ 1) over the redshift range. The results in the present study are consistent with both of observation and numerical simulations, where gas-poor minor merger is believed to be the main evolution scenario. By taking account of the connection with less massive QGs and SFGs, we discuss the formation and evolution of the massive QGs over "Cosmic High Noon", or the peak of star-formation in the universe.

In order to understand nature of building blocks of galaxies in the early universe, we investigate "genuine irregular galaxies (GIGs)" in the nearby universe. Here, GIGs are defined as isolated galaxies without regular structures (spheroid, bulge, disk, bar, spiral arm, and nucleus). Using the results of two excellent studies on galaxy morphology based on the Sloan Digital Sky Survey (SDSS), we obtain a sample of 66 irregular galaxies. We carry out new classification of them into GIGs and non-GIGs which have regular structure or show evidence for galaxy interaction, by using the SDSS Data Release 10 images. We then find that a half of these irregular galaxies (33/66) are GIGs and obtain an unambiguous sample of 33 GIGs for the first time. We discuss their observational properties by comparing them with those of elliptical, S0, spiral galaxies, and irregular galaxies without the GIGs. We find that our GIGs have smaller sizes, lower optical luminosities, bluer rest-frame optical colors, lower surface stellar mass densities, and lower gas metallicity than normal galaxies. All these properties suggest that they are in chemically and dynamically younger phases even in the nearby universe.

The properties of K-band selected galaxies (K_AB<24) in the z = 3.09 SSA22 protocluster field are studied. 430 galaxies at 2.6 < z_phot < 3.6 are selected as potential protocluster members in a 112 arcmin^2 area based on their photometric redshifts. We find that \approx 20% of the massive galaxies with stellar masses >10^11 M_sun at z_phot \sim 3.1 have colors consistent with those of quiescent galaxies with ages > 0.5 Gyr. This fraction increases to \approx 50% after correcting for unrelated foreground/background objects. We also find that 30% of the massive galaxies are heavily reddened dusty star-forming galaxies. Few such quiescent galaxies at similar redshifts are seen in typical survey fields. An excess surface density of 24\mu m sources at z_phot \sim 3.1 is also observed, implying the presence of dusty star-formation activity in the protocluster. Cross-correlation with the X-ray data indicates that the fraction of K-band selected protocluster galaxies hosting active galactic nuclei (AGN) is also high compared with the field. The sky distribution of the quiescent galaxies, the 24\mu m sources, and the X-ray AGNs show clustering around a density peak of z=3.1 Ly\alpha emitters (LAEs). A significant fraction of the massive galaxies have already become quiescent, while the dusty star-formation is still active in the SSA22 protocluster. These findings indicate that we are witnessing the formation epoch of massive early-type galaxies at the center of predecessors to present-day rich galaxy clusters.

Large-scale structures (LSS) out to z $< 3.0$ are measured in the Cosmic Evolution Survey (COSMOS) using extremely accurate photometric redshifts (photoz). The Ks-band selected sample (from Ultra-Vista) is comprised of 155,954 galaxies. Two techniques -- adaptive smoothing and Voronoi tessellation -- are used to estimate the environmental densities within 127 redshift slices. Approximately 250 statistically significant overdense structures are identified out to z $= 3.0$ with shapes varying from elongated filamentary structures to more circularly symmetric concentrations. We also compare the densities derived for COSMOS with those based on semi-analytic predictions for a $\Lambda$CDM simulation and find excellent overall agreement between the mean densities as a function of redshift and the range of densities. The galaxy properties (stellar mass, spectral energy distributions (SEDs) and star formation rates (SFRs)) are strongly correlated with environmental density and redshift, particularly at z $< 1.0 - 1.2$. Classifying the spectral type of each galaxy using the rest-frame b-i color (from the photoz SED fitting), we find a strong correlation of early type galaxies (E-Sa) with high density environments, while the degree of environmental segregation varies systematically with redshift out to z $\sim 1.3$. In the highest density regions, 80% of the galaxies are early types at z=0.2 compared to only 20% at z = 1.5. The SFRs and the star formation timescales exhibit clear environmental correlations. At z $> 0.8$, the star formation rate density (SFRD) is uniformly distributed over all environmental density percentiles, while at lower redshifts the dominant contribution is shifted to galaxies in lower density environments.

We investigated the fraction of [OII] emitters in galaxies at z~0.9 as a function of the local galaxy density in the Hubble Space Telescope (HST) COSMOS 2 square degree field. [OII] emitters are selected by the narrow-band excess technique with the NB711-band imaging data taken with Suprime-Cam on the Subaru telescope. We carefully selected 614 photo-z selected galaxies with M_U3500 < -19.31 at z=0.901-0.920, which includes 195 [OII] emitters, to directly compare results with our previous study at z~1.2. We found that the fraction is almost constant at 0.3 Mpc^-2 < Sigma_10th < 10 Mpc^-2. We also checked the fraction of galaxies with blue rest-frame colors of NUV-R < 2 in our photo-z selected sample, and found that the fraction of blue galaxies does not significantly depend on the local density. On the other hand, the semi-analytic model of galaxy formation predicted that the fraction of star-forming galaxies at z~0.9 decreases with increasing the projected galaxy density even if the effects of the projection and the photo-z error in our analysis were taken into account. The fraction of [OII] emitters decreases from ~60% at z~1.2 to ~30% at z~0.9 independent of the galaxy environment. The decrease of the [OII] emitter fraction could be explained mainly by the rapid decrease of the star formation activity in the universe from z~1.2 to z~0.9.

During the course of our Lyman continuum imaging survey, we found that the spectroscopically confirmed Lya emitter LAE 221724+001716 at z = 3.10 in the SSA 22 field shows strong Lyman continuum emission (lambda_rest ~ 900 A) that escapes from this galaxy. However, another recent spectroscopic survey revealed that the supposed Lyman continuum emission could arise from a foreground galaxy at z = 1.76 if the emission line newly detected from the galaxy at lambda_obs ~ 3360 A is Lya. If this is the case, as the angular separation between these two galaxies is very small (~ 0.6"), LAE 221724+001716 at z = 3.10 could be amplified by the gravitational lensing caused by this intervening galaxy. Here we present a possible gravitational lens model for the system of LAE 221724+001716. First, we estimate the stellar mass of the intervening galaxy as Mstar ~ 3.5x10^9 Msun from its UV luminosity and ~ 3.0x10^7--2.4x10^9 Msun through the SED fitting. Then, we find that the gravitational magnification factor ranges from 1.01 to 1.16 using the so-called singular isothermal sphere model for strong lensing. While LAE 221724+001716 is the first system of an LAE-LAE lensing reported so far, the estimated magnification factor is not so significant because the stellar mass of the intervening galaxy is small.

The M82 `cap' is a gas cloud at a projected radius of 11.6 kpc along the minor axis of this well known superwind source. The cap has been detected in optical line emission and X-ray emission and therefore provides an important probe of the wind energetics. In order to investigate the ionization source of the cap, we observed it with the Kyoto3DII Fabry-Perot instrument mounted on the Subaru Telescope. Deep continuum, Ha, [NII]6583/Ha, and [SII]6716,6731/Ha maps were obtained with sub-arcsecond resolution. The superior spatial resolution compared to earlier studies reveals a number of bright Ha emitting clouds within the cap. The emission line widths (< 100 km s^-1 FWHM) and line ratios in the newly identified knots are most reasonably explained by slow to moderate shocks velocities (v_shock = 40--80 km s^-1) driven by a fast wind into dense clouds. The momentum input from the M82 nuclear starburst region is enough to produce the observed shock. Consequently, earlier claims of photoionization by the central starburst are ruled out because they cannot explain the observed fluxes of the densest knots unless the UV escape fraction is very high (f_esc > 60%), i.e., an order of magnitude higher than observed in dwarf galaxies to date. Using these results, we discuss the evolutionary history of the M82 superwind. Future UV/X-ray surveys are expected to confirm that the temperature of the gas is consistent with our moderate shock model.

We present the result of our low-luminosity quasar survey in the redshift range of 4.5 < z < 5.5 in the COSMOS field. Using the COSMOS photometric catalog, we selected 15 quasar candidates with 22 < i' < 24 at z~5, that are ~ 3 mag fainter than the SDSS quasars in the same redshift range. We obtained optical spectra for 14 of the 15 candidates using FOCAS on the Subaru Telescope and did not identify any low-luminosity type-1 quasars at z~5 while a low-luminosity type-2 quasar at z~5.07 was discovered. In order to constrain the faint end of the quasar luminosity function at z~5, we calculated the 1sigma confidence upper limits of the space density of type-1 quasars. As a result, the 1sigma confidence upper limits on the quasar space density are Phi< 1.33*10^{-7} Mpc^{-3} mag^{-1} for -24.52 < M_{1450} < -23.52 and Phi< 2.88*10^{-7} Mpc^{-3} mag^{-1} for -23.52 < M_{1450} < -22.52. The inferred 1sigma confidence upper limits of the space density are then used to provide constrains on the faint-end slope and the break absolute magnitude of the quasar luminosity function at z~5. We find that the quasar space density decreases gradually as a function of redshift at low luminosity (M_{1450} ~ -23), being similar to the trend found for quasars with high luminosity (M_{1450}<-26). This result is consistent with the so-called downsizing evolution of quasars seen at lower redshifts.

We present our new deep optical imaging and long-slit spectroscopy for Arp 220 that is the archetypical ULIRG in the local universe. Our sensitive Ha imaging has newly revealed large-scale, Ha absorption, i.e., post-starburst regions in this merger; one is found in the eastern superbubble and the other is in the two tidal tails that are clearly reveled in our deep optical imaging. The size of Ha absorption region in the eastern bubble is 5 kpc x 7.5 kpc and the observed Ha equivalent widths are ~2 A +- 0.2 A. The sizes of the northern and southern Ha-absorption tidal tails are ~5 kpc x 10 kpc and ~6 kpc x 20 kpc, respectively. The observed Ha equivalent widths range from 4 A to 7 A. In order to explain the presence of the two post-starburst tails, we suggest a possible multiple-merger scenario for Arp 220 in which two post-starburst disk-like structures merged into one, and then caused the two tails. This favors that Arp 220 is a multiple merging system composed of four or more galaxies, arising from a compact group of galaxies. Taking our new results into account, we discuss a star formation history in the last 1 Gyr in Arp 220.

We present the results of wide-field deep JHK imaging of the SSA22 field using MOIRCS instrument equipped with Subaru telescope. The observed field is 112 arcmin^2 in area, which covers the z=3.1 protocluster characterized by the overdensities of Ly Alpha emitters (LAEs) and Ly Alpha Blobs (LABs). The 5 sigma limiting magnitude is K_{AB} = 24.3. We extract the potential protocluster members from the K-selected sample by using the multi-band photometric-redshift selection as well as the simple color cut for distant red galaxies (DRGs; J-K_{AB}>1.4). The surface number density of DRGs in our observed fields shows clear excess compared with those in the blank fields, and the location of the densest area whose projected overdensity is twice the average coincides with the large-scale density peak of LAEs. We also found that K-band counterparts with z_{phot} = 3.1 are detected for 75% (15/20) of the LABs within their Ly Alpha halo, and the 40 % (8/20) of LABs have multiple components, which gives a direct evidence of the hierarchical multiple merging in galaxy formation. The stellar mass ofLABs correlates with their luminosity, isophotal area, and the Ly Alpha velocity widths, implying that the physical scale and the dynamical motion of Ly Alpha emission are closely related to their previous star-formation activities. Highly dust-obscured galaxies such as hyper extremely red objects (HEROs; J-K_{AB}>2.1) and plausible K-band counterparts of submillimeter sources are also populated in the high density region.

Using the very deep Subaru images of the GOODS-N region, from the MOIRCS Deep Survey and images from the HST/ACS, we have measured the Luminosity Ratio (LR) of the outer to the central regions of massive (M>10^{10.5}M_{Sun}) galaxies at fixed radii in a single rest-frame for z<3.5 as a new approach to the problem of size evolution. We didn't observe any evolution in the median LR. Had a significant size growth occurred, the outer to central luminosity ratios would have demonstrated a corresponding increase with a decrease in redshift.

We present scaling relations between stellar-mass (Mstar) and the size of galaxies at 0.3 < z < 3 for half- (R_50) and 90 percent-light (R_90) radii, using a deep K-band selected catalogue taken with the Subaru Telescope and MOIRCS in the GOODS-North region. The logarithmic slope R \propto Mstar^{0.1-0.2} is independent of redshift in a wide mass range of Mstar ~ 10^8-10^11 Msun, irrespective of galaxy populations (star-forming, quiescent). The offset change is < 50 percent. Provided that optical light in the rest frame traces the stellar mass of galaxies, the universal relation demonstrates that the stellar mass was built up in galaxies over their cosmic histories in a similar manner on average irrelevant to galaxy mass. The small offset in each stellar mass bin from the universal relation shows weak size evolution at a given mass. There is a moderate increase of 30-50 percent for R_50 and R_90 for less massive galaxies (Mstar < 10^10 Msun) from z~3 to z~1, while the sizes remains unchanged or slightly decrease towards z~0.3. For massive galaxies (Mstar > 10^11 Msun), the evolution is ~70-80 % increase in R_90 from z~3 to z~0.3, though that in R_50 is weaker. The evolution of compactness factor, R_50/R_90, which becomes smaller at lower redshift, is suggestive of minor merging effect in the outer envelope of massive galaxies.

In order to understand environmental effects on star formation in high-redshift galaxies, we investigate the physical relationships between the star formation activity, stellar mass, and environment for z ~1.2 galaxies in the 2 deg^2 COSMOS field. We estimate star formation using the [OII] emission line and environment from the local galaxy density. Our analysis shows that for massive galaxies M_*>10^10 M_sun, the fraction of [OII] emitters in high-density environments is 1.7 times higher than in low-density environments, while the [OII] emitter fraction does not depend on environment for low-mass M_* < 10^10 M_sun galaxies. In order to understand what drives these trends, we investigate the role of companion galaxies in our sample. We find that the fraction of [OII] emitters in galaxies with companions is 2.4 times as high as that in galaxies without companions at M_* > 10^10 M_sun. In addition, massive galaxies are more likely to have companions in high-density environments. However, although the "number" of star forming galaxies increases for massive galaxies with close companions and in dense environments, the "average" star formation rate of star forming galaxies at a given mass is independent of environment and the presence/absence of a close companion. These results suggest that interactions and/or mergers in high-density environment could induce star formation in massive galaxies at z~1.2, increasing the fraction of star-forming galaxies with M_* > 10^10 M_sun.

We study structure of star-forming galaxies at z~2 in GOODS-N field selected as sBzK galaxies down to K_{AB} <24.0 mag. Among 1029 sBzK galaxies, 551 galaxies (54%) show a single component in ACS/F850LP image obtained with the Hubble Space Telescope; the rest show multiple components. We fit the single-component sBzK galaxies with the single S\'ersic profile using the ACS/F850LP image and find that a majority of them (64%) show S\'ersic index of n=0.5-2.5, indicating that they have a disk-like structure. The resulting effective radii typically range from 1.0 to 3.0 kpc in the rest-frame UV wavelength. After correcting the effective radii to those in the rest-frame optical wavelength, we find that the single-component sBzK galaxies locate in the region where the local and z~1 disk galaxies distribute in the stellar mass-size diagram, suggesting comparable surface stellar mass density between the sBzK and z~0-1 disk galaxies. All these properties suggest that the single-component sBzK galaxies are progenitors of the present-day disk galaxies. However, by studying their intrinsic shape through comparison between the observed distribution of apparent axial ratios and the distribution for triaxial models with axes (A>B>C), we find that the mean B/A ratio is 0.61^{+0.05}_{-0.08} and disk thickness C/A is 0.28^{+0.03}_{-0.04}. This indicates that the single-component sBzK galaxies at z~2 have a bar-like or oval shape rather than a round disk shape. The shape seems to resemble to a bar/oval structure that form through bar instability; if it is the case, the intrinsic shape may give us a clue to understand dynamical evolution of baryonic matter in a dark matter halo.

Based on the ultraviolet to far-infrared photometry already compiled and presented in a companion paper (Barro et al. 2011a, Paper I), we present a detailed SED analysis of nearly 80,000 IRAC 3.6+4.5 micron selected galaxies in the Extended Groth Strip. We estimate photometric redshifts, stellar masses, and star formation rates separately for each galaxy in this large sample. The catalog includes 76,936 sources with [3.6] < 23.75 (85% completeness level of the IRAC survey) over 0.48 square degrees. The typical photometric redshift accuracy is Delta z/(1+z)=0.034, with a catastrophic outlier fraction of just 2%. We quantify the systematics introduced by the use of different stellar population synthesis libraries and IMFs in the calculation of stellar masses. We find systematic offsets ranging from 0.1 to 0.4 dex, with a typical scatter of 0.3 dex. We also provide UV- and IR-based SFRs for all sample galaxies, based on several sets of dust emission templates and SFR indicators. We evaluate the systematic differences and goodness of the different SFR estimations using the deep FIDEL 70 micron data available in the EGS. Typical random uncertainties of the IR-bases SFRs are a factor of two, with non-negligible systematic effects at z$\gtrsim$1.5 observed when only MIPS 24 micron data is available. All data products (SEDs, postage stamps from imaging data, and different estimations of the photometric redshifts, stellar masses, and SFRs of each galaxy) described in this and the companion paper are publicly available, and they can be accessed through our the web-interface utility Rainbow-navigator

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 an IRAC 3.6+4.5 microns selected catalog in the Extended Groth Strip (EGS) containing photometry from the ultraviolet to the far-infrared and stellar parameters derived from the analysis of the multi-wavelength data. In this paper, we describe the method used to build coherent spectral energy distributions (SEDs) for all the sources. In a companion paper, we analyze those SEDs to obtain robust estimations of stellar parameters such as photometric redshifts, stellar masses, and star formation rates. The catalog comprises 76,936 sources with [3.6]<23.75 mag (85% completeness level of the IRAC survey in the EGS) over 0.48 square degrees. For approximately 16% of this sample, we are able to deconvolve the IRAC data to obtain robust fluxes for the multiple counterparts found in ground-based optical images. Typically, the SEDs of the IRAC sources in our catalog count with more than 15 photometric data points, spanning from the UV to the FIR. Approximately 95% and 90% of all IRAC sources are detected in the deepest optical and near-infrared bands. Only 10% of the sources have optical spectroscopy and redshift estimations. Almost 20% and 2% of the sources are detected by MIPS at 24 and 70 microns, respectively. We also cross-correlate our catalog with public X-ray and radio catalogs. Finally, we present the Rainbow Navigator public web-interface utility designed to browse all the data products resulting from this work, including images, spectra, photometry, and stellar parameters.

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 the discovery of a significant excess of candidate Halpha emitters (HAEs) in the field of the radio galaxy 4C 23.56 at z=2.483. Using the MOIRCS near-infrared imager on the Subaru Telescope we found 11 candidate emission-line galaxies to a flux limit of ~7.5 10^-17 erg s-1 cm-2, which is about 5 times excess from the expected field counts with ~3-sigma significance. Three of these are spectroscopically confirmed as redshifted Halpha at z=2.49. The distribution of candidate emitters on the sky is tightly confined to a 1.2-Mpc-radius area at z=2.49, locating 4C 23.56 at the western edge of the distribution. Analysis of the deep Spitzer MIPS 24 mu m imaging shows that there is also an excess of faint MIPS sources. All but two of the 11 HAEs are also found in the MIPS data. The inferred star-formation rate (SFR) of the HAEs based on the extinction-corrected Halpha luminosity (median SFR >~100 M_solar yr-1) is similar to those of HAEs in random fields at z~2. On the other hand, the MIPS-based SFR for the HAEs is on average 3.6 times larger, suggesting the existence of the star-formation significanly obscured by dust. The comparison of the Halpha-based star-formation activities of the HAEs in the 4C 23.56 field to those in another proto-cluster around PKS 1138-262 at z=2.16 reveals that the latter tend to have fainter Halpha emission despite similar K-band magnitudes. This suggests that star-formation may be suppressed in the PKS 1138-262 protocluster relative to the 4C 23.56 protocluster. This difference among the HAEs in the two proto-clusters at z > 2 may imply that some massive cluster galaxies are just forming at these epochs with some variation among clusters.

We present results from Submillimeter Array (SMA) 860-micron sub-arcsec astrometry and multiwavelength observations of the brightest millimeter (S_1.1mm = 8.4 mJy) source, SSA22-AzTEC1, found near the core of the SSA22 protocluster that is traced by Ly\alpha emitting galaxies at z = 3.09. We identify a 860-micron counterpart with a flux density of S_860um = 12.2 +/- 2.3 mJy and absolute positional accuracy that is better than 0.3". At the SMA position, we find radio to mid-infrared counterparts, whilst no object is found in Subaru optical and near-infrared deep images at wavelengths \le 1 micron (J > 25.4 in AB, 2\sigma). The photometric redshift estimate, using flux densities at \ge 24 microns, indicates z_phot = 3.19^{+0.26}_{-0.35}, consistent with the protocluster redshift. We then model the near-to-mid-infrared spectral energy distribution (SED) of SSA22-AzTEC1, and find that the SED modeling requires a large extinction (A_V \approx 3.4 mag) of starlight from a stellar component with M_star ~ 10^{10.9} M_sun, assuming z = 3.1. Additionally, we find a significant X-ray counterpart with a very hard spectrum (Gamma_eff = -0.34 ^{+0.57}_{-0.61}), strongly suggesting that SSA22-AzTEC1 harbors a luminous AGN (L_X ~ 3*10^{44} ergs s^{-1}) behind a large hydrogen column (N_H ~ 10^{24} cm^{-2}). The AGN, however, is responsible for only ~10% of the bolometric luminosity of the host galaxy, and therefore the star-formation activity likely dominates the submillimeter emission. It is possible that SSA22-AzTEC1 is the first example of a protoquasar growing at the bottom of the gravitational potential underlying the SSA22 protocluster.

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 near-infrared (NIR) spectroscopy for 37 BzK-color-selected star-forming galaxies with MOIRCS on Subaru. The sample is drawn from the Ks-band selected catalog of the MOIRCS Deep Survey (MODS) in the GOODS-N region. About half of our samples are selected from the MIPS 24um-source catalog. Ha emission lines are detected from 23 galaxies, of which the median redshift is 2.12. We derived the star formation rates (SFRs) from extinction-corrected Ha luminosities. The extinction correction is estimated from the SED fitting of multi-band photometric data covering UV to NIR wavelengths. The Balmer decrement of the stacked emission lines shows that the amount of extinction for the ionized gas is larger than that for the stellar continuum. From a comparison of the extinction corrected Ha luminosity and other SFR indicators we found that the relation between the dust properties of stellar continuum and ionized gas is different depending on the intrinsic SFR. The comparison between the Ha SFR and stellar mass estimated from SED fitting shows no correlation between them. Some galaxies with stellar mass smaller than ~10^10 Msun show SFRs higher than ~100 Msun/yr. The specific SFRs (SSFRs) of these galaxies are remarkably high. The average stellar-population age of these high-SSFR galaxies is younger than 100 Myr, which is consistent with the implied high SSFR. The large SFR implies the possibility that the high SSFR galaxies significantly contribute to the cosmic SFR density of the universe at z~2. The total SFR density estimated from the Ha emission line galaxies is 0.089-0.136 Msun/yr/Mpc^3, which is consistent with the total SFR densities in the literature. The metallicity of the high-SSFR galaxies is larger than that expected from the mass-metallicity relation of UV-selected galaxies at z~2 by Erb et al. (2006a). [Abridged]

We present a population study of several types of galaxies within the protocluster surrounding the radio galaxy MRC0316-257 at z~3.1. In addition to the known population of Ly_alpha emitters (LAEs) and [OIII] emitters, we use colour selection techniques to identify protocluster candidates that are Lyman break galaxies (LBG) and Balmer break galaxies (BBGs). The radio galaxy field contains an excess of LBG candidates, with a surface density 1.6\pm0.3 times larger than found for comparable blank fields. This surface overdensity corresponds to an LBG volume overdensity of ~8\pm4. The BBG photometric redshift distribution peaks at the protocluster's redshift, but we detect no significant surface overdensity of BBG. This is not surprising because a volume overdensity similar to the LBGs would have resulted in a surface density of ~1.2 that found in the blank field. This could not have been detected in our sample. Masses and star formation rates of the candidate protocluster galaxies are determined using SED fitting. These properties are not significantly different from those of field galaxies. The galaxies with the highest masses and star formation rates are located near the radio galaxy, indicating that the protocluster environment influences galaxy evolution at z~3. We conclude that the protocluster around MRC0316-257 is still in the early stages of formation.

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 spectroscopic campaign to follow-up red colour-selected candidate massive galaxies in two high redshift proto-clusters surrounding radio galaxies. We observed a total of 57 galaxies in the field of MRC0943-242 (z=2.93) and 33 in the field of PKS1138-262 (z=2.16) with a mix of optical and near-infrared multi-object spectroscopy. We confirm two red galaxies in the field of PKS1138-262 at the redshift of the radio galaxy. Based on an analysis of their spectral energy distributions, and their derived star formation rates from the H-alpha and 24um flux, one object belongs to the class of dust-obscured star-forming red galaxies, while the other is evolved with little ongoing star formation. This result represents the first red and mainly passively evolving galaxy to be confirmed as companion galaxies in a z>2 proto-cluster. Both red galaxies in PKS1138-262 are massive, of the order of 4-6x10^11 M_Sol. They lie along a Colour-Magnitude relation which implies that they formed the bulk of their stellar population around z=4. In the MRC0943-242 field we find no red galaxies at the redshift of the radio galaxy but we do confirm the effectiveness of our JHK_s selection of galaxies at 2.3<z<3.1, finding that 10 out of 18 (56%) of JHK_s-selected galaxies whose redshifts could be measured fall within this redshift range. We also serendipitously identify an interesting foreground structure of 6 galaxies at z=2.6 in the field of MRC0943-242. This may be a proto-cluster itself, but complicates any interpretation of the red sequence build-up in MRC0943-242 until more redshifts can be measured.

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 investigate the X-ray properties of the K-band-selected galaxies at redshift 2 < z < 4 by using our deep near-infrared images obtained in the MOIRCS Deep Survey project and the published Chandra X-ray source catalog. 61 X-ray sources with the 2-10 keV luminosity L_X = 10^{42}-10^{44} erg/s are identified with the K-selected galaxies and we found that they are exclusively (90%) associated with the massive objects with stellar mass larger than 10^{10.5} Msun. Our results are consistent with the idea that the M_BH/M_str ratio of the galaxies at z=2-4 is similar to the present-day value. On the other hand, the AGN detection rate among the very massive galaxies with the stellar mass larger than 10^{11} Msun is high, 33% (26/78). They are active objects in the sense that the black-hole mass accretion rate is ~ 1-50% of the Eddington limit if they indeed have similar M_BH/M_str ratio with those observed in the local universe. The active duration in the AGN duty cycle of the high-redshift massive galaxies seems large.

We present the results of deep near-infrared imaging observations of the z=3.1 proto-cluster region in the SSA22a field taken by MOIRCS mounted on the Subaru Telescope. We observed a 21.7 arcmin^2 field to the depths of J=24.5, H=24.3, and K=23.9 (5 sigma). We examine the distribution of the K-selected galaxies at z~3 by using the simple color cut for distant red galaxies (DRGs) as well as the photometric-redshift selection technique. The marginal density excess of DRGs and the photo-z selected objects are found around the two most luminous Ly alpha blobs (LABs). We investigate the correlation between the K-selected objects and the LABs, and find that several galaxies with stellar mass M_* = 10^9-10^11 M_solar exist in vicinity of LABs, especially around the two most luminous ones. We also find that 7 of the 8 LABs in the field have plausible K_s-band counterparts and the sum of the stellar mass possibly associated with LABs correlates with the luminosity and surface brightness of them, which implies that the origin of Ly alpha emission may be closely correlated with their stellar mass or their previous star formation phenomena.

We present multi-waveband optical imaging data obtained from observations of the Subaru/XMM-Newton Deep Survey (SXDS). The survey field, centered at R.A.=02:18:00, decl.=-05:00:00, has been the focus of a wide range of multi-wavelength observing programs spanning from X-ray to radio wavelengths. A large part of the optical imaging observations are carried out with Suprime-Cam on Subaru Telescope at Mauna Kea in the course of Subaru Telescope Observatory Projects. This paper describes our optical observations, data reduction and analysis procedures employed, and the characteristics of the data products. A total area of 1.22 sqdeg is covered in five contiguous sub-fields, each of which corresponds to a single Suprime-Cam field of view (34'x27'), in five broad-band filters B, V, Rc, i', z' to the depths of B=28.4, V=27.8, Rc=27.7, i'=27.7 and z'=26.6 (AB, 3-sigma, 2-arcsec aperture). The data are reduced and compiled into five multi-waveband photometric catalogs, separately for each Suprime-Cam pointing. The i'-band catalogs contain about 900,000 objects, making the SXDS catalogs one of the largest multi-waveband catalogs in corresponding depth and area coverage. The SXDS catalogs can be used for an extensive range of astronomical applications such as the number density of the Galactic halo stars to the large scale structures at the distant universe. The number counts of galaxies are derived and compared with those of existing deep extragalactic surveys. The optical data, the source catalogs, and configuration files used to create the catalogs are publicly available via the SXDS web page (http://www.naoj.org/Science/SubaruProject/SXDS/index.html)

We explore the evolved galaxy population in the proto-clusters around four high-z radio galaxies at 2<~z<~3 based on wide-field near-infrared imaging. Three of the four fields are known proto-clusters as demonstrated by overdensities of line emitting galaxies at the same redshifts as the radio galaxies. We imaged the fields of three targets (PKS1138-262, USS0943-242 and MRC0316-257) to a depth of Ks~22 (5sigma) over a 4'x7' area centered on the radio galaxies with a new wide-field NIR camera, MOIRCS, on the Subaru Telescope. Another target (USS1558-003) was observed with SOFI on the NTT to a depth of Ks=20.5 over a 5'x5' area. We apply colour cuts in J-Ks and/or JHKs in order to exclusively search for galaxies located at high redshifts: z>2. To the 5sigma limiting magnitudes, we see a significant excess of NIR selected galaxies by a factor of two to three compared to those found in the field of GOODS-South. The spatial distribution of these NIR selected galaxies is not uniform and traces structures similar to those of emission line galaxies, although the samples of NIR selected galaxies and emitters show little overlap. We focus on the NIR colour-magnitude sequence of the evolved population and find that the bright-end (M_{stars}>10^{11}Msun) of the red sequence is well populated by z~2 but much less so in the z~3 proto-clusters. This may imply that the bright-end of the colour-magnitude sequence first appeared between z=3 and 2, an era coinciding with the appearance of submm galaxies and the peak of the cosmic star formation rate. Our observations show that during the same epoch, massive galaxies are forming in high density environments by vigorous star formation and assembly.

(Abridged) We report a discovery of possible large-scale structures around the RDCS J1252.9-2927 cluster at z=1.24 based on photometric redshifts. We carried out multi-band wide-field imaging with Suprime-Cam on the Subaru Telescope and WFCAM on the United Kingdom Infra-Red Telescope (UKIRT). The distribution of photo-z selected galaxies reveals clumpy structures surrounding the central cluster. We compare the observed structure with an X-ray map and find that two of the four plausible clumps show significant X-ray emissions and one with a marginal detection, which strongly suggest that they are dynamically bound systems. Following the discovery of the possible large-scale structure, we carried out deeper SOFI K_s-band imaging with New Technology Telescope on the four plausible clumps. We construct the optical-to-near-infrared colour-magnitude diagrams of the galaxies in the clumps, and find that the colour-magnitude relation (CMR) of the red galaxies in the clumps is sharply truncated below K_s=22. Interestingly, the main cluster shows a clear relation down to K_s=23 (Lidman et al. 2004). We suggest that galaxies follow the 'environment-dependent down-sizing' evolution. Massive galaxies in high density environments first stop forming stars and become red. Less massive galaxies in less dense environments become red at later times. Based on a few assumptions, we predict that the brightest tip of the CMR appears at z~2.5.

We present the first measurement of clustering properties of low mass galaxies with a stellar mass down to M_*~10^9 Msun at 1<z<4 in 24.4 arcmin^2 of the GOODS-North region with a depth of K_{AB}~25, based on the near infrared observations performed with MOIRCS at the Subaru Telescope. The correlation amplitude strongly depends on the K-band flux, color, and stellar mass of the galaxies. We find that K-band luminous galaxies have a larger correlation length than K-band faint galaxies. For color selected samples at 2<z<4, distant red galaxies with J-K>1.3 show a large bias of b~7.2+-1.3 on scales up to \theta~100" or 3.1 comoving Mpc, while blue galaxies with 0.5<J-K<1.3, in which most Lyman break galaxies are populated, have a weak clustering signal on large scales, but a possible strong small scale excess at \theta<10". For massive galaxies with M_*>~10^{10} Msun, we estimate a correlation length and bias to be r_0~4.5 h^{-1} Mpc and b=1.9-3.5, which are much larger than those of low mass (M_*~10^9-10^{10} Msun) galaxies. The comparison of our measurements with analytic CDM models constrains the properties of hosting dark halos, and indicates that the low mass galaxies would be progenitors of galaxies with a typical luminosity of L<~L_* in the local Universe. The blue galaxies in low mass samples are more strongly clustered in more massive halos with higher occupation numbers than low mass red galaxies. This fact suggests an environment effect due to the halo mass on star formation activity at high-z.

We use very deep near-infrared imaging data taken with Multi-Object InfraRed Camera and Spectrograph (MOIRCS) on the Subaru Telescope to investigate the number counts of Distant Red Galaxies (DRGs). We have observed a 4x7 arcmin^2 field in the Great Observatories Origins Deep Survey North (GOODS-N), and our data reach J=24.6 and K=23.2 (5sigma, Vega magnitude). The surface density of DRGs selected by J-K>2.3 is 2.35+-0.31 arcmin^-2 at K<22 and 3.54+-0.38 arcmin^-2 at K<23, respectively. These values are consistent with those in the GOODS-South and FIRES. Our deep and wide data suggest that the number counts of DRGs turn over at K~22, and the surface density of the faint DRGs with K>22 is smaller than that expected from the number counts at the brighter magnitude. The result indicates that while there are many bright galaxies at 2<z<4 with the relatively old stellar population and/or heavy dust extinction, the number of the faint galaxies with the similar red color is relatively small. Different behaviors of the number counts of the DRGs and bluer galaxies with 2<z_phot<4 at K>22 suggest that the mass-dependent color distribution, where most of low-mass galaxies are blue while more massive galaxies tend to have redder colors, had already been established at that epoch.

We use deep multi-band optical and near-infrared data for four general fields, GOODS-South, HDF North/South, and IRAC UDF in GOODS-North to investigate the evolution of the observed rest-frame U-V color of field galaxies as a function of the stellar mass. In these four fields, the U-V color distributions of the galaxies at each mass and redshift interval are very similar. At 0.3<z<2.7, we found that more massive galaxies always tend to have a redder U-V color. High- and low-mass galaxies exhibit quite different color evolutions. As seen in our previous study in HDF-N, the color distribution of low-mass (M_stellar < 3x10^9M_solar) galaxies becomes significantly bluer with an increase in the redshift. This evolution of the average color can be explained by a constant star formation rate model with z_form~4. On the other hand, the average color of high-mass galaxies (M_stellar > 3x10^10M_solar) evolves more strongly at a high redshift. Such mass-dependent color distribution and its evolution indicate that galaxies with a larger stellar mass appear to have shorter star-formation timescales, and on an average, they form the larger fraction of their stars in the earlier epoch.

We report a discovery of proto-cluster candidates around high redshift radio galaxies at z~2.5 on the basis of clear statistical excess of colour-selected galaxies around them seen in the deep near-infrared imaging data obtained with CISCO on Subaru Telescope. We have observed six targets, all at similar redshifts at z~2.5, and our data reach to J=23.5, H=22.6 and K=21.8 (5 sigma) and cover a 1.6' x 1.6' field centered on each radio galaxy. We apply colour cuts in JHKs in order to exclusively search for galaxies located at high redshifts, z>2. Over the magnitude range of 19.5<K<21.5 we see a significant excess of red galaxies with J-K>2.3 by a factor of two around the combined radio galaxies fields compared to those found in the general field of GOODS South. The excess of galaxies around the radio galaxies fields becomes more than factor of three around 19.5<K<20.5 when the two-colours cuts are applied with JHKs. Such overdensity of the colour-selected galaxies suggest that those fields tend to host high density regions at high redshifts, although there seems to be the variety of the density of the colour-selected galaxies in each field. In particular, two radio galaxies fields out of the six observed fields show very strong density excess and these are likely to be proto-clusters associated to the radio galaxies which would evolve into rich clusters of galaxies dominated by old passively evolving galaxies.

We obtained the number counts and the rest-frame B-band luminosity function of the color-selected old passively-evolving galaxies (OPEGs) at z=1 with very high statistical accuracy using a large and homogeneous sample of about 4000 such objects with z' <25 detected in the area of 1.03 deg^2 in the Subaru/XMM-Newton Deep Survey (SXDS) field. Our selection criteria are defined on the i'-z' and R-z' color-magnitude plane so that OPEGs at z=0.9-1.1 with formation redshift z_f=2-10 are properly sampled. The limiting magnitude corresponds to the luminosity of galaxies with M_*+3 at z=0. We made a pilot redshift observations for 99 OPEG candidates with 19 < z' < 22 and found that at least 78% (73/93) of the entire sample, or 95% (73/77) of these whose redshifts were obtained are indeed lie between z=0.87 and 1.12 and the most of their spectra show the continuum break and strong Ca H and K lines, indicating that these objects are indeed dominated by the old stellar populations. We then compare our results with the luminosity functions of the color- or the morphologically-selected early type galaxies at z=0 taking the evolutionary factor into account and found that the number density of old passive galaxies with sim M_* magnitude at z~1 averaged over the SXDS area is 40-60% of the equivalently red galaxies and 60-85% of the morphologically-selected E/S0 galaxies at z=0 depending on their luminosity evolution. It is revealed that more than half, but not all, of the present-day early-type galaxies had already been formed into quiescent passive galaxies at z=1.

Using the Subaru very deep K'-band imaging and HST WFPC2/NICMOS archival data of the Hubble Deep Field North, we investigate the evolution of the stellar mass, color, morphology of galaxies to z ~ 3. We mainly examine the rest-frame U-V color distribution of galaxies as a function of stellar mass. At 0.3<z<2, galaxies seem to be divided into the two populations at around the stellar mass of ~5x10^9 M_solar. The low-mass galaxies have relatively bluer rest U-V color and their color does not show clear correlation with stellar mass over the range of 10^8-5x10^9 M_solar. On the other hand, at higher mass, the more massive galaxies tend to have the redder U-V color. The average U-V color of the low-mass galaxies becomes bluer gradually with redshift, from U-V ~ 0.2 at z ~ 0.5 to U-V ~ -0.2 at z ~ 2. On the contrary, the correlation between the stellar mass and rest U-V color of the high-mass population does not seem to change significantly between z ~ 0.3 and z ~ 2. At z>2, it is seen that more massive galaxies tend to have redder U-V color over the range of 10^9-10^10 M_solar. These results suggests that the star formation history of galaxies depends on their stellar mass very much. The low-mass population is likely to have relatively long star formation timescale. At the stellar mass larger than ~5x10^9 M_solar, there must be some mechanisms which suppress the star formation in galaxies at 0<z<2.

We have obtained a new high-resolution K'-band image of the central region of the rich X-ray cluster RX J0848.9+4452 at z=1.26. We found that the brightest cluster galaxy (BCG) in the cluster is clearly separated into two distinct objects. Whereas the optical to near-infrared colors of the objects are consistent with the predictions of passive evolution models for galaxies formed at high redshift, the luminosities of the two galaxies are both considerably fainter than predicted by passive evolution of BCG's in low and intermediate redshift clusters. We argue that this is evidence of an on-going merger of normal cluster ellipticals to form the dominant galaxy in the core of RX J0848.9+4452. The two galaxies appear to point towards the nearby cluster ClG J0848+4453 and are aligned with the outer X-ray contour of their parent cluster, supporting a model of BCG formation by collimated infall along the surrounding large-scale structure.

Despite extensive observational efforts, the brightest sub--mm source in the Hubble Deep Field, HDF850.1, has failed to yield a convincing optical/infrared identification almost 4 years after its discovery. This failure is all the more notable given the availability of supporting multi-frequency data of unparalleled depth, and sub-arcsec positional accuracy for the sub-mm/mm source. Consequently, HDF850.1 has become a test case of the possibility that the most violently star-forming objects in the universe are too red and/or distant to be seen in the deepest optical images. Here we report the discovery of the host galaxy of HDF850.1. This object has been revealed by careful analysis of a new, deep K-prime image of the HDF obtained with the Subaru 8.2-m telescope. Its reality is confirmed by a similar analysis of the HST NICMOS F160W image of the same region. This object is extremely faint (K=23.5), clumpy (on sub-arcsec scales) and very red (I-K > 5.2; H-K = 1.4 +- 0.35). The likelihood that it is the correct identification is strongly reinforced by a reanalysis of the combined MERLIN+VLA 1.4-GHz map of the field which yields a new radio detection of HDF850.1 only 0.1 arcsec from the new near-ir counterpart, and with sufficient positional accuracy to exclude all previously considered alternative optical candidates. We have calculated new confidence limits on the estimated redshift of HDF850.1 and find z = 4.1 +- 0.5. We also calculate that the flux density of HDF850.1 has been boosted by a factor of ~3 through lensing by the intervening elliptical 3-586.0, consistent with predictions that a small but significant fraction of blank-field sub-mm sources are lensed by foreground galaxies. We discuss the wider implications of these results for the sub-mm population and cosmic star-formation history.

We present the results of the deep NIR imaging observations using the Subaru telescope. At the Hubble Deep Field North, we obtained the new and deepest K'-band image (10 hours) to study the rest-frame optical morphology of galaxies at z > 3 and the stellar mass distribution of galaxies at 2 < z < 4.5. We also study the rest-frame optical properties of the 'building blocks' in the field of 53W002 and overall color distribution there. We then combined these data to compare the color and magnitude distribution in these two fields. The opt-NIR color-magnitude distributions in the two fields looks very similar and we report the conspicuous color change below K_AB =22 and discuss the cause of the feature.

We discuss the properties of galaxies around the radio galaxy 3C 324 at z=1.2 based on BVRIK' multi-band imaging data. We have applied a photometric-redshift technique to objects in the 3C 324 field, and identified 35 objects as plausible cluster members. We have found that red and luminous members are concentrated in a small region enclosed by a circle of 40" radius (0.33 Mpc at z=1.2 for \Omega_0=0.3, \lambda_0=0.7, H_0=70km/s/Mpc cosmology) from the 3C 324 galaxy. The 3C 324 cluster is probably much more compact in size compared with the local clusters. We constructed a K'-band luminosity function of the cluster members and fit a Schechter function, and found the characteristic magnitude to be K*_{AB}=20.2+-0.6. This value is consistent with the extrapolation of the pure passive evolution seen for z<1 clusters. We have identified eight bright galaxies which form a red color-magnitude sequence. The slope of the sequence is consistent with the passive evolution model down to K'_{AB}<22; we also found that there is no clear age variation in these bright red galaxies. However, seven out of these eight galaxies exhibit a significant excess in the rest UV light with respect to the passive evolution model. This may suggest that the massive early-type galaxies in this high-redshift cluster are still forming stars to some extent. We have confirmed a truncation of the color-magnitude sequence at K'_{AB}\sim22; faint passively-evolving galaxies may not yet be present in this cluster at z\sim1.2. The overall color distribution of the cluster members, selected by the photometric redshift technique, is found to be very broad. We derived the fraction of blue galaxies in this cluster following a

Using the HST WFPC2/NICMOS archival data of the Hubble Deep Field North, we construct the nearly complete sample of the Mv <-20 (~L*+1) galaxies to z=2, and investigate when the Hubble sequence appeared, namely, the evolution of the morphology, colors, and the comoving number density of the sample. Even if taking into account of the uncertainty of the photometric redshift technique, the number density of relatively bright bulge-dominated galaxies in the HDF-N decrease significantly at z>1, and their rest-frame U-V color distribution is wide-spread over 0.5<z<2. On the other hand, while the number density of both disk-dominated and irregular galaxies does not show significant change at z<2, their distribution of the rest-frame U-V color alters at z~1.5: there is no relatively red (rest U-V>0.3) galaxies at z>1.5, while the significant fraction of these red disk-dominated or irregular galaxies exists at z<1.5. These results suggest that the significant evolution of the Hubble sequence which is seen in the present Universe occurs at 1<z<2.

We have obtained high-resolution K'-band images of the powerful z=1.206 radio galaxy 3C 324 with the Subaru telescope under seeing conditions of 0.3--0.4 arcsec. We clearly resolved the galaxy and directly compared it to the optical images obtained with the Hubble Space Telescope. The host galaxy of 3C 324 is revealed to be a moderately luminous elliptical galaxy with a smooth light profile. The effective radius of the galaxy, as determined by profile fitting, is 1.3+-0.1 arcsec (1.2 kpc), which is significantly smaller than the value of 2.2 arcsec, published in Best et al. (1998, MNRAS, 292, 758). The peak of the K'-band light coincides with the position of the radio core, which implies that the powerful AGN lies at the nucleus of the host galaxy. The peak also coincides with the gap in the optical knotty structures which may be a dust lane hiding the UV-optical emission of the AGN from our line of sight; it is very likely that we are seeing the obscuring structure almost edge-on. We clearly detected the `aligned component' in the K'-band image by subtracting a model elliptical galaxy from the observed image. The red R_F702W-K color of the outer region of the galaxy avoiding the aligned component indicates that the near infrared light of the host galaxy is dominated by an old stellar population.

We have investigated the optical and near-infrared colors of K'-selected galaxies in clusters at z ~ 1.2 near to the radio galaxy 3C 324 using images obtained with the Subaru telescope and archival HST data. The distribution of colors of the galaxies in the cluster region is found to be fairly broad, and it may imply significant scatter in their star-formation histories, although the effect of contamination of field galaxies is uncertain. The red sequence of galaxies whose R-K colors are consistent with passive evolution models for old galaxies is found to be truncated at K' ~ 20 mag, and there are few fainter galaxies with similar red colors in the cluster region. We find that the bulge-dominated galaxies selected by quantitative morphological classification form a broad sequence in the color-magnitude diagram, whose slope is much steeper than that expected from metallicity variations within a passively evolving coeval galaxy population. We argue that the observed color-magnitude sequence can be explained by metallicity and age variations, and the fainter galaxies with K' > 20 mag may be 1-2 Gyr younger than the brighter galaxies. Some spatial segregation of the color and K'-band luminosity is seen in the sky distribution; the redder and the brighter objects tend to be located near 3C 324.

We investigate the $K$-band luminosity distribution of galaxies in the region of clusters at $z\sim1.2$ near to the radio galaxy 3C 324. The imaging data were obtained during the commissioning period of the Subaru telescope. There is a significant excess of the surface number density of the galaxies with $K =$ 17--20 mag in the region within $\sim$ 40'' from 3C 324. At this bright end, the measured luminosity distribution shows a drop, which can be represented by the exponential cut off of the Schechter-function formula; the best-fitted value of the characteristic magnitude, $K^{*}$, is $\sim 18.4\pm0.8$. This measurement follows the evolutionary trend of the $K^*$ of the rich clusters observed at an intermediate redshift, which is consistent with passive evolution models with a formation redshift $z_f \gtsim 2$. At $K \gtsim 20$ mag, however, the excess of the galaxy surface density in the region of the clusters decreases abruptly, which may imply that the luminosity function of the cluster galaxies has a negative slope at the faint end. This may imply strong luminosity segregation between the inner and outer parts of the clusters, or some deficit of faint galaxies in the cluster central region of the cluster.

We investigated the color and morphology of optically selected galaxies in the region of clusters at z $\sim$ 1.2 near to the radio galaxy 3C 324 using archived data taken with the Hubble Space Telescope. The faint galaxies selected at the HST F702W band that contribute to the surface-density excess of the region have wide ranges of color, size, and morphology, which are not likely to be due to contamination by foreground galaxies. Namely, the rest-frame ultraviolet emission properties of the galaxies in the clusters are not very homogeneous; various amounts of star-formation activity may occur in a significant fraction of them. Although our analysis is purely statistical, we find that typical star-forming galaxies with blue colors have a relatively late-type morphology compared to the red quiescent population in the systems.