The unique combination of superb spatial resolution, wide-area coverage, and deep depth of the optical imaging from the Hyper Suprime-Cam (HSC) Subaru Strategic Program is utilized to search for dual quasar candidates. Using an automated image analysis routine on 34,476 known SDSS quasars, we identify those with two (or more) distinct optical point sources in HSC images covering 796 deg^2. We find 421 candidates out to a redshift of 4.5 of which one hundred or so are more likely after filtering out contaminating stars. Angular separations of 0.6 - 4.0" correspond to projected separations of 3 - 30 kpc, a range relatively unexplored for population studies of luminous dual quasars. Using Keck-I/LRIS and Gemini-N/NIFS, we spectroscopically confirm three dual quasar systems at z < 1, two of which are previously unknown out of eight observed, based on the presence of characteristic broad emission lines in each component, while highlighting that the continuum of one object in one of the pairs is reddened. In all cases, the [OIII]5007 emission lines have mild velocity offsets, thus the joint [OIII] line profile is not double-peaked. We find a dual quasar fraction of 0.26+/-0.18% and no evidence for evolution. A comparison with the Horizon-AGN simulation seems to support the case of no evolution in the dual quasar fraction when broadly matching the quasar selection. These results may indicate a scenario in which the frequency of the simultaneous triggering of luminous quasars is not as sensitive as expected to the cosmic evolution of the merger rate or gas content of galaxies.

A proto-cluster core is the most massive dark matter halo (DMH) in a given proto-cluster. To reveal the galaxy formation in core regions, we search for proto-cluster cores at $z\sim 2$ in $\sim 1.5\, \mathrm{deg}^{2}$ of the COSMOS field. Using pairs of massive galaxies ($\log(M_{*}/M_{\odot})\geq11$) as tracers of cores, we find 75 candidate cores, among which 54\% are estimated to be real. A clustering analysis finds that these cores have an average DMH mass of $2.6_{-0.8}^{+0.9}\times 10^{13}\, M_{\odot}$, or $4.0_{-1.5}^{+1.8}\, \times 10^{13} M_{\odot}$ after contamination correction. The extended Press-Schechter model shows that their descendant mass at $z=0$ is consistent with Fornax-like or Virgo-like clusters. Moreover, using the IllustrisTNG simulation, we confirm that pairs of massive galaxies are good tracers of DMHs massive enough to be regarded as proto-cluster cores. We then derive the stellar mass function (SMF) and the quiescent fraction for member galaxies of the 75 candidate cores. We find that the core galaxies have a more top-heavy SMF than field galaxies at the same redshift, showing an excess at $\log(M_{*}/M_{\odot})\gtrsim 10.5$. The quiescent fraction, $0.17_{-0.04}^{+0.04}$ in the mass range $9.0\leq \log(M_{*}/M_{\odot})\leq 11.0$, is about three times higher than that of field counterparts, giving an environmental quenching efficiency of $0.13_{-0.04}^{+0.04}$. These results suggest that stellar mass assembly and quenching are accelerated as early as at $z\sim 2$ in proto-cluster cores.

Based mainly on X-ray observations, studies are made on interactions between the intra-cluster medium (ICM) in clusters of galaxies and their member galaxies. Through (magneto)hydrodynamic and gravitational channels, the moving galaxies are expected to drag the ICM around them, and transfer to the ICM some fraction of their dynamical energies on cosmological time scales. This hypothesis is in line with several observations, including the possible cosmological infall of galaxies towards the cluster center, found over redshifts of z~1 to z~0. Further assuming that the energy lost by the galaxies is first converted into ICM turbulence and then dissipated, this picture can explain the subsonic and uniform ICM turbulence, measured with Hitomi in the core region of the Perseus cluster. The scenario may also explain several other unanswered problems regarding clusters of galaxies, including what prevents the ICM from the expected radiative cooling, how the various mass components in nearby clusters have attained different radial distributions, and how a thermal stability is realized between hot and cool ICM components that co-exist around cD galaxies. This view is also considered to pertain to the general scenario of galaxy evolution, including their environmental effects.

We first present new Subaru narrow-band observations of the Ly$\alpha$ halo around the quasi-stellar object (QSO) CFHQ J232908$-$030158 at $z=6.42$, which appears the most luminous and extended halo at $z>5$ ($L_{Ly\alpha}=9.8\times10^{43}$ erg s$^{-1}$ within $37$ pkpc diameter). Then, combining these measurements with available data in the literature, we find two different evolutions of QSOs' Ly$\alpha$ halos. First is a possible short-term evolution with QSO age seen in four $z>6$ QSOs. We find the anti-correlation between the Ly$\alpha$ halo scales with QSOs' IR luminosity, with J2329-0301's halo being the brightest and largest. It indicates that ionizing photons escape more easily out to circum-galactic regions when host galaxies are less dusty. We also find a positive correlation between IR luminosity and black hole mass ($M_\text{BH}$). Given $M_\text{BH}$ as an indicator of QSO age, we propose a hypothesis that a large Ly$\alpha$ halo mainly exists around QSOs in the young phase of their activity due to a small amount of dust. The second is an evolution with cosmic time seen over $z\sim2-5$. We find the increase of surface brightness toward lower-redshift with a similar growth rate to that of dark matter halos (DHs) which evolve to $M_\text{DH}=10^{12}-10^{13}$ M$_\odot$ at $z=2$. The extent of Ly$\alpha$ halos is also found to increase at a rate scaling with the virial radius of growing DHs, $r_\text{vir} \propto M_\text{DH}^{1/3}$($1+z$)$^{-1}$. These increases are consistent with a scenario that the CGM around QSOs evolves in mass and size keeping pace with hosting DHs.

This is the second paper of a series to shed light on how galactic properties depend on the intergalactic medium (IGM) environment traced by the Ly$\alpha$ forest. In this paper, we observationally investigate the IGM--galaxy connection using the publicly available 3D IGM tomography data (CLAMATO: Lee et al. 2016, 2018) and several galaxy catalogs in the COSMOS field. We measure the cross-correlation function (CCF) for $570$ galaxies with spec-$z$ measurements, and detect a correlation with the IGM up to $50$ $h^{-1}$ comoving Mpc. We show that galaxies with stellar masses of $10^9-10^{10}$ M$_\odot$ are the dominant contributor to the total CCF signal. We also investigate CCFs for several galaxy populations: Ly$\alpha$ emitters (LAEs), H$\alpha$ emitters (HAEs), [OIII] emitters (O3Es), active galactic nuclei (AGNs), and submillimeter galaxies (SMGs), and detect the highest signal in AGNs and SMGs at large scales ($r\geq5$ $h^{-1}$ Mpc), but in LAEs at small scales ($r<5$ $h^{-1}$ Mpc). We find that they live in various IGM environments --HAEs trace the IGM in a similar manner to the continuum-selected galaxies, but LAEs and O3Es tend to reside in higher density regions. Additionally, LAEs' CCF is flat up to $r\sim3$ $h^{-1}$ Mpc, indicating that they tend to avoid the highest density regions. For AGNs and SMGs, the CCF peak at $r=5-6$ $h^{-1}$ Mpc implies that they tend to be in locally lower density regions. We suspect that it is due to the IGM HI photoionization by AGNs, i.e., the proximity effect.

This is the first paper in a series to systematically investigate the environmental dependence of galactic properties traced by Ly$\alpha$ forest absorption in the interagalactic medium (IGM). Using our cosmological hydrodynamic simulations, we investigate the IGM--galaxy connection at $z=2$ by two methods: (I) cross-correlation analysis between galaxies and the fluctuation of Ly$\alpha$ forest transmission ($\delta_\text{F}$); and (II) comparing the overdensity of neutral hydrogen (HI) and galaxies. Our simulations reproduce observed cross-correlation functions (CCF) between the Ly$\alpha$ forest and Lyman-break galaxies. We further investigate the variation of the CCF using subsamples divided by dark matter halo mass ($M_\text{DH}$), galaxy stellar mass ($M_*$), and star-formation rate (SFR), and find that the CCF signal becomes stronger with increasing $M_\text{DH}$, $M_*$, and SFR. The CCFs between galaxies and gas-density fluctuation are also found to have similar trends. Therefore, the variation of the $\delta_\text{F}$--CCF depending on $M_\text{DH}$, $M_*$, and SFR is due to varying gas density around galaxies. We conclude that $M_\text{DH}$ is the most sensitive parameter to characterize the correlation between galaxies and the IGM distribution. Method (II) also finds correlations between galaxies and HI overdensities. Overall, our analyses confirm the spatial correlation between galaxies and IGM HI, with more massive galaxies being clustered in higher-density regions as expected from the $\Lambda$CDM paradigm.

The CFHT Large Area U-band Deep Survey (CLAUDS) uses data taken with the MegaCam mosaic imager on CFHT to produce images of 18.60 deg2 with median seeing of FWHM=0.92 arcsec and to a median depth of U = 27.1 AB (5 sigma in 2 arcsec apertures), with selected areas that total 1.36 deg2 reaching a median depth of U=27.7 AB. These are the deepest U-band images assembled to date over this large an area. These data are located in four fields also imaged to comparably faint levels in grizy and several narrowband filters as part of the Hyper Suprime-Cam (HSC) Subaru Strategic Program (HSC-SSP). These CFHT and Subaru datasets will remain unmatched in their combination of area and depth until the advent of the Large Synoptic Survey Telescope (LSST). This paper provides an overview of the scientific motivation for CLAUDS and gives details of the observing strategy, observations, data reduction, and data merging with the HSC-SSP. Three early applications of these deep data are used to illustrate the potential of the dataset: deep U-band galaxy number counts, z~3 Lyman break galaxy (LBG) selection, and photometric redshifts improved by adding CLAUDS U to the Subaru HSC grizy photometry.

We report the Subaru Hyper Suprime-Cam (HSC) discovery of two Ly$\alpha$ blobs (LABs), dubbed z70-1 and z49-1 at $z=6.965$ and $z=4.888$ respectively, that are Ly$\alpha$ emitters with a bright ($\log L_{\rm Ly\alpha}/{\rm [erg\ s^{-1}]}>43.4$) and spatially-extended Ly$\alpha$ emission, and present the photometric and spectroscopic properties of a total of seven LABs; the two new LABs and five previously-known LABs at $z=5.7-6.6$. The z70-1 LAB shows the extended Ly$\alpha$ emission with a scale length of $1.4\pm 0.2$ kpc, about three times larger than the UV continuum emission, making z70-1 the most distant LAB identified to date. All of the 7 LABs, except z49-1, exhibit no AGN signatures such as X-ray emission, {\sc Nv}$\lambda$1240 emission, or Ly$\alpha$ line broadening, while z49-1 has a strong {\sc Civ}$\lambda$1548 emission line indicating an AGN on the basis of the UV-line ratio diagnostics. We carefully model the point-spread functions of the HSC images, and conduct two-component exponential profile fitting to the extended Ly$\alpha$ emission of the LABs. The Ly$\alpha$ scale lengths of the core (star-forming region) and the halo components are $r_{\rm c}=0.6-1.2$ kpc and $r_{\rm h}=2.0-13.8$ kpc, respectively. The average $r_{\rm h}$ of the LABs falls on the extrapolation of the $r_{\rm h}$-Ly$\alpha$ luminosity relation of the Ly$\alpha$ halos around VLT/MUSE star-forming galaxies at the similar redshifts, suggesting that typical LABs at $z\gtrsim5$ are not special objects, but star-forming galaxies at the bright end.

The physical origin of diffuse Ly$\alpha$ halos (LAHs) around star-forming galaxies is still a matter of debate. We present the dependence of LAH luminosity ($L({\rm Ly}\alpha)_H$) on the stellar mass ($M_\star$), $SFR$, color excess ($E(B-V)_\star$), and dark matter halo mass ($M_{\rm h}$) of the parent galaxy for $\sim 900$ Ly$\alpha$ emitters (LAEs) at $z\sim2$ divided into ten subsamples. We calculate $L({\rm Ly}\alpha)_H$ using the stacked observational relation between $L({\rm Ly}\alpha)_H$ and central Ly$\alpha$ luminosity by Momose et al. (2016), which we find agrees with the average trend of VLT/MUSE-detected individual LAEs. We find that our LAEs have relatively high $L({\rm Ly}\alpha)_H$ despite low $M_\star$ and $M_{\rm h}$, and that $L({\rm Ly}\alpha)_H$ remains almost unchanged with $M_\star$ and perhaps with $M_{\rm h}$. These results are incompatible with the cold stream (cooling radiation) scenario and the satellite-galaxy star-formation scenario, because the former predicts fainter $L({\rm Ly}\alpha)_H$ and both predict steeper $L({\rm Ly}\alpha)_H$ vs. $M_\star$ slopes. We argue that LAHs are mainly caused by Ly$\alpha$ photons escaping from the main body and then scattered in the circum-galactic medium. This argument is supported by LAH observations of H$\alpha$ emitters (HAEs). When LAHs are taken into account, the Ly$\alpha$ escape fractions of our LAEs are about ten times higher than those of HAEs with similar $M_\star$ or $E(B-V)_\star$, which may partly arise from lower HI gas masses implied from lower $M_{\rm h}$ at fixed $M_\star$, or from another Ly$\alpha$ source in the central part.

We report on the relation between the mass of supermassive black holes (SMBHs; M_BH) and that of hosting dark matter halos (M_h) for 49 z ~ 6 quasi-stellar objects (QSOs) with [CII]158um velocity-width measurements. Here, we estimate M_h assuming that the rotation velocity from FWHM_CII is equal to the circular velocity of the halo; we have tested this procedure using z ~ 3 QSOs that also have clustering-based M_h estimates. We find that a vast majority of the z ~ 6 SMBHs are more massive than expected from the local M_BH - M_h relation, with one-third of the sample by factors >~ 10^2. The median mass ratio of the sample, M_BH/M_h = 6 x 10^{-4}, means that 0.4% of the baryons in halos are locked up in SMBHs. The mass growth rates of our SMBHs amount to ~ 10% of the SFRs, or ~ 1% of the mean baryon accretion rates, of the hosting galaxies. A large fraction of the hosting galaxies are consistent with average galaxies in terms of SFR and perhaps of stellar mass and size. Our study indicates that the growth of SMBHs (M_BH ~ 10^{8-10} Msun) in luminous z ~ 6 QSOs greatly precedes that of hosting halos owing to efficient gas accretion even under normal star formation activities, although we cannot rule out the possibility that undetected SMBHs have local M_BH/M_h ratios. This preceding growth is in contrast to much milder evolution of the stellar-to-halo mass ratio.

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.

The survey of Lyman $\alpha$ emitters (LAEs) with Subaru Hyper Suprime-Cam, called SILVERRUSH (Ouchi et al.), is producing massive data of LAEs at $z\gtrsim6$. Here we present LAE simulations to compare the SILVERRUSH data. In 162$^3$ comoving Mpc$^3$ boxes, where numerical radiative transfer calculations of reionization were performed, LAEs have been modeled with physically motivated analytic recipes as a function of halo mass. We have examined $2^3$ models depending on the presence or absence of dispersion of halo Ly$\alpha$ emissivity, dispersion of the halo Ly$\alpha$ optical depth, $\tau_\alpha$, and halo mass dependence of $\tau_\alpha$. The unique free parameter in our model, a pivot value of $\tau_\alpha$, is calibrated so as to reproduce the $z=5.7$ Ly$\alpha$ luminosity function (LF). We compare our model predictions with Ly$\alpha$ LFs at $z=6.6$ and $7.3$, LAE angular auto-correlation functions (ACFs) at $z=5.7$ and $6.6$, and LAE fractions in Lyman break galaxies at $5<z<7$. The Ly$\alpha$ LFs and ACFs are reproduced by multiple models, but the LAE fraction turns out to be the most critical test. The dispersion of $\tau_\alpha$ and the halo mass dependence of $\tau_\alpha$ are essential to explain all observations reasonably. Therefore, a simple model of one-to-one correspondence between halo mass and Ly$\alpha$ luminosity with a constant Ly$\alpha$ escape fraction has been ruled out. Based on our best model, we present a formula to estimate the intergalactic neutral hydrogen fraction, $x_{\rm HI}$, from the observed Ly$\alpha$ luminosity density at $z\gtrsim6$. We finally obtain $x_{\rm HI}=0.5_{-0.3}^{+0.1}$ as a volume-average at $z=7.3$.

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

We construct $z\sim6-7$, 8, and 9 faint Lyman break galaxy samples (334, 61, and 37 galaxies, respectively) with accurate size measurements with the software $\texttt{glafic}$ from the complete Hubble Frontier Fields (FF) cluster and parallel fields data. These are the largest samples hitherto and reach down to the faint ends of recently obtained deep luminosity functions. At faint magnitudes, however, these samples are highly incomplete for galaxies with large sizes, implying that derivation of the luminosity function (LF) sensitively depends on the intrinsic size--luminosity (RL) relation. We thus conduct simultaneous maximum-likelihood estimation of LF and RL relation parameters from the observed distribution of galaxies on the RL plane with help of a completeness map as a function of size and luminosity. At $z\sim6-7$, we find that the intrinsic RL relation expressed as $r_\textrm{e} \propto L^\beta$ has a notably steeper slope of $\beta=0.46^{+0.08}_{-0.09}$ than those at lower redshifts, which in turn implies that the LF has a relatively shallow faint-end slope of $\alpha=-1.86^{+0.17}_{-0.18}$. This steep $\beta$ can be reproduced by a simple analytical model in which smaller galaxies have lower specific angular momenta. The $\beta$ and $\alpha$ values for the $z\sim8$ and 9 samples are consistent with those for $z\sim6-7$ but with larger errors. For all three samples there is a large, positive covariance between $\beta$ and $\alpha$, implying that the simultaneous determination of these two parameters is important. We also provide new strong lens mass models of Abell S1063 and Abell 370 as well as updated mass models of Abell 2744 and MACS J0416.1$-$2403.

We present UV luminosity functions of dropout galaxies at $z\sim6-10$ with the complete Hubble Frontier Fields data. We obtain a catalog of $\sim450$ dropout-galaxy candidates (350, 66, and 40 at $z\sim6-7$, 8, and 9, respectively), whose UV absolute magnitudes reach $\sim-14$ mag, $\sim2$ mag deeper than the Hubble Ultra Deep Field detection limits. We carefully evaluate number densities of the dropout galaxies by Monte-Carlo simulations, including all lensing effects such as magnification, distortion, and multiplication of images as well as detection completeness and contamination effects in a self-consistent manner. We find that UV luminosity functions at $z\sim6-8$ have steep faint-end slopes, $\alpha\sim-2$, and likely steeper slopes, $\alpha\lesssim-2$ at $z\sim9-10$. We also find that the evolution of UV luminosity densities shows a non-accelerated decline beyond $z\sim8$ in the case of $M_\mathrm{trunc}=-15$, while an accelerated in the case of $M_\mathrm{trunc}=-17$. We examine whether our results are consistent with the Thomson scattering optical depth from the Planck satellite and the ionized hydrogen fraction $Q_\mathrm{HII}$ at $z\lesssim7$ based on the standard analytic reionization model. We find that there exist reionizaiton scenarios that consistently explain all the observational measurements with the allowed parameters of $f_{\rm esc}=0.17^{+0.07}_{-0.03}$ and $M_\mathrm{trunc}>-14.0$ for $\log\xi_\mathrm{ion}/[\mathrm{erg}^{-1}\ \mathrm{Hz}]=25.34$, where $f_{\rm esc}$ is the escape fraction, $M_\mathrm{trunc}$ is the faint limit of the UV luminosity function, and $\xi_\mathrm{ion}$ is the conversion factor of the UV luminosity to the ionizing photon emission rate. The length of the reionization period is estimated to be $\Delta z=3.9^{+2.0}_{-1.6}$ (for $0.1<Q_\mathrm{HII}<0.99$), consistent with the recent estimate from Planck.

We present average stellar population properties and dark matter halo masses of $z \sim 2$ \lya emitters (LAEs) from SED fitting and clustering analysis, respectively, using $\simeq$ $1250$ objects ($NB387\le25.5$) in four separate fields of $\simeq 1$ deg$^2$ in total. With an average stellar mass of $10.2\, \pm\, 1.8\times 10^8\ {\mathrm M_\odot}$ and star formation rate of $3.4\, \pm\, 0.4\ {\mathrm M_\odot}\ {\rm yr^{-1}}$, the LAEs lie on an extrapolation of the star-formation main sequence (MS) to low stellar mass. Their effective dark matter halo mass is estimated to be $4.0_{-2.9}^{+5.1} \times 10^{10}\ {\mathrm M_\odot}$ with an effective bias of $1.22^{+0.16}_{-0.18}$ which is lower than that of $z \sim 2$ LAEs ($1.8\, \pm\, 0.3$), obtained by a previous study based on a three times smaller survey area, with a probability of $96\%$. However, the difference in the bias values can be explained if cosmic variance is taken into account. If such a low halo mass implies a low HI gas mass, this result appears to be consistent with the observations of a high \lya escape fraction. With the low halo masses and ongoing star formation, our LAEs have a relatively high stellar-to-halo mass ratio (SHMR) and a high efficiency of converting baryons into stars. The extended Press-Schechter formalism predicts that at $z=0$ our LAEs are typically embedded in halos with masses similar to that of the Large Magellanic Cloud (LMC); they will also have similar SHMRs to the LMC, if their SFRs are largely suppressed after $z \sim 2$ as some previous studies have reported for the LMC itself.

We combined deep U-band imaging from the KPNO-4m/MOSAIC camera with very deep multi-waveband data from the optical to infrared, to select Lyman Break Galaxies (LBGs) at z~3 using U-V and V-R colors in the Subaru Deep Field. With the resulting sample of 5161 LBGs, we construct the UV luminosity function down to $M_{UV} = -18$ and find a steep faint-end slope of $\alpha=-1.78 \pm 0.05$. We analyze rest-frame UV-to-IR spectral energy distributions generated from the median optical photometry and photometry on median-stacked IR images. In the stacks of faint LBGs, we find a background depression centered on the galaxy. This deficit results from the systematic difficulty of SExtractor in finding faint galaxies in regions with higher-than-average surface densities of foreground galaxies. We corrected our stacked magnitudes for this. Best-fit stellar population templates for the stacked LBG SEDs indicate stellar masses and star-formation rates of log M*/Msun = 10 and 50 M$_\odot$/yr at i' = 24, down to log M*/Msun = 8 and = 3 M$_\odot$/yr at i' = 27. For the faint stacked LBGs there is a 1-mag excess over the expected stellar continuum in the K-band, which we attribute to redshifted [OIII]4959+5007 and H$\beta$ lines. Their implied equivalent widths increase with decreasing mass, reaching $\rm{EW_0([O III]4959,5007+H\beta)}$ =1500A in the faintest bin. Such strong [OIII] emission is seen only in a miniscule fraction of the most extreme local emission-line galaxies, but it probably universal in the faint galaxies that reionized the universe. Finally, we analyze clustering by computing the angular correlation function and performing halo occupation distribution (HOD) analysis. We find a mean dark halo mass of log(Mhalo/h) Msun = 11.29$\pm 0.12$ for the full sample of LBGs, and log(Mhalo/h) Msun = 11.49$\pm 0.1$ for the brightest half.

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.

The stellar disk size of a galaxy depends on the ratio of the disk stellar mass to the halo mass, $m_\star \equiv M_\star/M_{\rm dh}$, and the fraction of the dark halo angular momentum transferred to the stellar disk, $j_\star \equiv J_\star/J_{\rm dh}$. Since $m_{\star}$ and $j_{\star}$ are determined by many star-formation related processes, measuring $j_\star$ and $m_\star$ at various redshifts is essential to understand the formation history of disk galaxies. We use the 3D-HST GOODS-S, COSMOS, and AEGIS imaging data and photo-$z$ catalog to examine $j_\star$ and $m_\star$ for star-forming galaxies at $z \sim$ 2, 3, and 4, when disks are actively forming. We find that the $j_\star/m_\star$ ratio is $\simeq 0.77\pm 0.06$ for all three redshifts over the entire mass range examined, $8\times 10^{10} < M_{\rm dh}/h^{-1} M_\odot < 2\times 10^{12}$, with a possible ($<30\%$) decrease with mass. This high ratio is close to those of local disk galaxies, descendants of our galaxies in terms of $M_{\rm dh}$ growth, implying a nearly constant $j_\star/m_\star$ over past 12 Gyr. These results are remarkable because mechanisms controlling angular momentum transfer to disks such as inflows and feedbacks depend on both cosmic time and halo mass and indeed theoretical studies tend to predict $j_\star/m_\star$ changing with redshift and mass. It is found that recent theoretical galaxy formation simulations predict smaller $j_{\star}/m_{\star}$ than our values. We also find that a significant fraction of our galaxies appears to be unstable against bar formation.

We study the UV luminosity functions (LFs) at $z\sim 4$, $5$, $6,$ and $7$ based on the deep large-area optical images taken by the Hyper Suprime-Cam (HSC) Subaru strategic program (SSP). On the 100 deg$^2$ sky of the HSC SSP data available to date, we make enormous samples consisting of a total of 579,565 dropout candidates at $z\sim 4-7$ by the standard color selection technique, 358 out of which are spectroscopically confirmed by our follow-up spectroscopy and other studies. We obtain UV LFs at $z \sim 4-7$ that span a very wide UV luminosity range of $\sim 0.002 - 100 \, L_{\rm UV}^\ast$ ($-26 < M_{\rm UV} < -14$ mag) by combining LFs from our program and the ultra-deep Hubble Space Telescope legacy surveys. We derive three parameters of the best-fit Schechter function, $\phi^\ast$, $M_{\rm UV}^\ast$, and $\alpha$, of the UV LFs in the magnitude range where the AGN contribution is negligible, and find that $\alpha$ and $\phi^\ast$ decrease from $z\sim 4$ to $7$ with no significant evolution of $M_{\rm UV}^\ast$. Because our HSC SSP data bridge the LFs of galaxies and AGNs with great statistical accuracy, we carefully investigate the bright end of the galaxy UV LFs that are estimated by the subtraction of the AGN contribution either aided with spectroscopy or the best-fit AGN UV LFs. We find that the bright end of the galaxy UV LFs cannot be explained by the Schechter function fits at $> 2 \sigma$ significance, and require either double power-law functions or modified Schechter functions that consider a magnification bias due to gravitational lensing.

We present clustering properties from 579,492 Lyman break galaxies (LBGs) at z~4-6 over the 100 deg^2 sky (corresponding to a 1.4 Gpc^3 volume) identified in early data of the Hyper Suprime-Cam (HSC) Subaru strategic program survey. We derive angular correlation functions (ACFs) of the HSC LBGs with unprecedentedly high statistical accuracies at z~4-6, and compare them with the halo occupation distribution (HOD) models. We clearly identify significant ACF excesses in 10"<$\theta$<90", the transition scale between 1- and 2-halo terms, suggestive of the existence of the non-linear halo bias effect. Combining the HOD models and previous clustering measurements of faint LBGs at z~4-7, we investigate dark-matter halo mass (Mh) of the z~4-7 LBGs and its correlation with various physical properties including the star-formation rate (SFR), the stellar-to-halo mass ratio (SHMR), and the dark matter accretion rate (dotMh) over a wide-mass range of Mh/M$_\odot$=4x10^10-4x10^12. We find that the SHMR increases from z~4 to 7 by a factor of ~4 at Mh~1x10^11 M$_\odot$, while the SHMR shows no strong evolution in the similar redshift range at Mh~1x10^12 M$_\odot$. Interestingly, we identify a tight relation of SFR/dotMh-Mh showing no significant evolution beyond 0.15 dex in this wide-mass range over z~4-7. This weak evolution suggests that the SFR/dotMh-Mh relation is a fundamental relation in high-redshift galaxy formation whose star formation activities are regulated by the dark matter mass assembly. Assuming this fundamental relation, we calculate the cosmic SFR densities (SFRDs) over z=0-10 (a.k.a. Madau-Lilly plot). The cosmic SFRD evolution based on the fundamental relation agrees with the one obtained by observations, suggesting that the cosmic SFRD increase from z~10 to 4-2 (decrease from z~4-2 to 0) is mainly driven by the increase of the halo abundance (the decrease of the accretion rate).

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

We 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.

Active galactic nuclei (AGNs) in the high-redshift Universe are thought to reside in overdense environments. However, recent works provide controversial results partly due to the use of different techniques and possible suppression of nearby galaxy formation by AGN feedback. We conducted deep and wide-field imaging observations with the Suprime-Cam on the Subaru telescope and searched for Lyman-alpha emitters (LAEs) around two QSOs (quasi-stellar objects) at $z\sim4.9$ and a radio galaxy at $z\sim4.5$ by using narrow-band filters to address these issues more robustly. In the QSO fields, we obtained additional broad-band images to select Lyman-break galaxies (LBGs) at $z\sim5$ for comparison. We constructed a photometric sample of 301 LAEs and 170 LBGs in total. A wide field of view (34arcmin$\times$27arcmin, corresponding to 80$\times$60 comoving Mpc$^2$) of the Suprime-Cam enabled us to probe galaxies in the immediate vicinities of the AGNs and in the blank fields simultaneously and compare various properties of them in a consistent manner. The two QSOs are located near local density peaks ($<2\sigma$) and one of the QSOs has a close companion LAE with projected separation of 80 physical kpc. The radio galaxy is found to be near a void of LAEs. The number densities of LAEs/LGBs in a larger spatial scale around the AGNs are not significantly different from those in blank fields. No sign of feedback is found down to $L_{Ly\alpha}\sim10^{41.8}\mathrm{~erg~s^{-1}}$. Our results suggest that high-redshift AGNs are not associated with extreme galaxy overdensity and that this cannot be attributed to the effect of AGN feedback.

We present physical properties of spectroscopically confirmed Ly$\alpha$ emitters (LAEs) with very large rest-frame Ly$\alpha$ equivalent widths EW$_{\rm 0}$(Ly$\alpha$). Although the definition of large EW$_{\rm 0}$(Ly$\alpha$) LAEs is usually difficult due to limited statistical and systematic uncertainties, we identify six LAEs selected from $\sim 3000$ LAEs at $z\sim 2$ with reliable measurements of EW$_{\rm 0}$ (Ly$\alpha$) $\simeq 200-400$ \AA\ given by careful continuum determinations with our deep photometric and spectroscopic data. These large EW$_{\rm 0}$(Ly$\alpha$) LAEs do not have signatures of AGN, but notably small stellar masses of $M_{\rm *} = 10^{7-8}$ $M_{\rm \odot}$ and high specific star-formation rates (star formation rate per unit galaxy stellar mass) of $\sim 100$ Gyr$^{-1}$. These LAEs are characterized by the median values of $L({\rm Ly\alpha})=3.7\times 10^{42}$ erg s$^{-1}$ and $M_{\rm UV}=-18.0$ as well as the blue UV continuum slope of $\beta = -2.5\pm0.2$ and the low dust extinction $E(B-V)_{\rm *} = 0.02^{+0.04}_{-0.02}$, which indicate a high median Ly$\alpha$ escape fraction of $f_{\rm esc}^{\rm Ly\alpha}=0.68\pm0.30$. This large $f_{\rm esc}^{\rm Ly\alpha}$ value is explained by the low {\sc Hi} column density in the ISM that is consistent with FWHM of the Ly$\alpha$ line, ${\rm FWHM (Ly\alpha)}=212\pm32$ km s$^{-1}$, significantly narrower than those of small EW$_{\rm 0}$(Ly$\alpha$) LAEs. Based on the stellar evolution models, our observational constraints of the large EW$_{\rm 0}$ (Ly$\alpha$), the small $\beta$, and the rest-frame He{\sc ii} equivalent width imply that at least a half of our large EW$_{\rm 0}$(Ly$\alpha$) LAEs would have young stellar ages of $\lesssim 20$ Myr and very low metallicities of $Z<0.02 Z_\odot$ regardless of the star-formation history.

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 present clustering analysis results from 10,381 Lyman break galaxies (LBGs) at z~ 4-7, identified in the Hubble legacy deep imaging and new complimentary large-area Subaru/Hyper Suprime-Cam data. We measure the angular correlation functions (ACFs) of these LBGs at z~4, 5, 6, and 7, and fit these measurements using halo occupation distribution (HOD) models that provide an estimate of halo masses, M_h~(1-20)x10^11 Msun. Our M_h estimates agree with those obtained by previous clustering studies in a UV-magnitude vs. M_h plane, and allow us to calculate stellar-to-halo mass ratios (SHMRs) of LBGs. By comparison with the z~0 SHMR, we identify evolution of the SHMR from z~0 to z~4, and z~4 to z~7 at the >98% confidence levels. The SHMR decreases by a factor of ~2 from z~0 to 4, and increases by a factor of ~4 from z~4 to 7. We compare our SHMRs with results of a hydrodynamic simulation and a semi-analytic model, and find that these theoretical studies do not predict the SHMR increase from z~4 to 7. We obtain the baryon conversion efficiency (BCE) of LBGs at z~4, and find that the BCE increases with increasing dark matter halo mass. Finally, we compare our clustering+HOD estimates with results from abundance matching techniques, and conclude that the M_h estimates of the clustering+HOD analyses agree with those of the simple abundance matching within a factor of 3, and that the agreement improves when using more sophisticated abundance matching techniques that include subhalos, incompleteness, and/or evolution in the star formation and stellar mass functions.

We present the Lya luminosity functions (LFs) derived by our deep Subaru narrowband survey that identifies a total of 3,137 Lya emitters (LAEs) at $z = 2.2$ in five independent blank fields. The sample of these LAEs is the largest, to date, and covers a very wide Lya luminosity range of $\log L_{Ly\alpha} = 41.7-44.4$ erg s$^{-1}$. We determine the Lya LF at $z = 2.2$ with unprecedented accuracies, and obtain the best-fit Schechter parameters of $L^{*}_{Ly\alpha} = 5.29^{+1.67}_{-1.13} \times 10^{42}$ erg s$^{-1}$, $\phi^{*}_{Ly\alpha} = 6.32^{+3.08}_{-2.31} \times 10^{-4}$ Mpc$^{-3}$, and $\alpha = -1.75^{+0.10}_{-0.09}$ showing a steep faint-end slope. We identify a significant hump at the LF bright end ($\log L_{Ly\alpha} > 43.4$ erg s$^{-1}$). Because all of the LAEs in the bright-end hump have (a) bright counterpart(s) either in the X-ray, UV, or radio data, this bright-end hump is not made by gravitational lensing magnification bias but AGNs. These AGNs allow us to derive the AGN UV LF at $z \sim 2$ down to the faint magnitude limit of $M_{UV} \simeq -22.5$, and to constrain the faint-end slope of AGN UV LF, $\alpha_{AGN}=-1.2 \pm 0.1$, that is flatter than those at $z > 4$. Based on the Lya and UV LFs from our and previous studies, we find the increase of Lya escape fraction $f^{Ly\alpha}_{esc}$ from $z \sim 0$ to $6$ by two orders of magnitude. This large $f^{Ly\alpha}_{esc}$ increase can be explained neither by the evolution of stellar population nor outflow alone, but the evolution of neutral hydrogen HI density in inter-stellar medium that enhances dust attenuation for Lya by resonance scattering. Our uniform expanding shell models suggest that the typical HI column density decreases from $N_{HI} \sim 7 \times 10^{19}$ ($z \sim 0$) to $\sim 1 \times 10^{18}$ cm$^{-2}$ ($z \sim 6$) to explain the large $f^{Ly\alpha}_{esc}$ increase.

We present the results of a Lya profile analysis of 12 Lya emitters (LAEs) at z = 2.2 with high-resolution Lya spectra. We find that all 12 objects have a Lya profile with the main peak redward of the systemic redshift defined by nebular lines, and five have a weak, secondary peak blueward of the systemic redshift (blue bump). The average velocity offset of the red main peak (the blue bump, if any) with respect to the systemic redshift is Delta_v_Lya,r = 174+/- 19 km s-1 (Delta_v_Lya,b = -316+/-45 km s-1), which is smaller than (comparable to) that of Lyman-break galaxies (LBGs). The outflow velocities inferred from metal absorption lines in three individual and one stacked spectra are comparable to those of LBGs. The expanding shell model constructed by Verhamme et al. (2006) reproduces not only the Lya profiles but also other observed quantities including the outflow velocity and the FWHM of nebular lines for the non-blue bump objects. On the other hand, the model predicts too high FWHMs of nebular lines for the blue bump objects, although this discrepancy may disappear if we introduce additional Lya photons produced by gravitational cooling. We show that the small Delta_v_Lya,r values of our sample can be explained by low neutral-hydrogen column densities of log(NHI) = 18.9 cm-2 on average. This value is more than one order of magnitude lower than those of LBGs but is consistent with recent findings that LAEs have high ionization parameters and low Hi gas masses. This result suggests that low NHI values, giving reduced numbers of resonant scattering of Lya photons, are the key to the strong Lya emission of LAEs.

We conduct precise strong lensing mass modeling of four ${\it Hubble}$ Frontier Fields (HFF) clusters, Abell$~$2744, MACS$~$J0416.1$-$2403, MACS$~$J0717.5$+$3745, and MACS$~$J1149.6$+$2223, for which HFF imaging observations are completed. We construct a refined sample of more than 100 multiple images for each cluster by taking advantage of the full depth HFF images, and conduct mass modeling using the ${\small \rm GLAFIC}$ software, which assumes simply parametrized mass distributions. Our mass modeling also exploits a magnification constraint from the lensed Type Ia supernova HFF14Tom for Abell$~$2744 and positional constraints from the multiple images S1$-$S4 of the lensed supernova SN Refsdal for MACS$~$J1149.6$+$2223. We find that our best-fitting mass models reproduce the observed image positions with RMS errors of $\sim 0.4$ arcsec, which are smaller than RMS errors in previous mass modeling that adopted similar numbers of multiple images. Our model predicts a new image of SN Refsdal with a relative time delay and magnification that are fully consistent with a recent detection of reappearance. We then construct catalogs of $z\sim 6-9$ dropout galaxies behind the four clusters and estimate magnification factors for these dropout galaxies with our best-fitting mass models. The dropout sample from the four cluster fields contains $\sim 120$ galaxies at $z\gtrsim 6$, about 20 of which are predicted to be magnified by a factor of more than 10. Some of the high-redshift galaxies detected in the HFF have lensing-corrected magnitudes of $M_{\rm UV}\sim -15$ to $-14$. Our analysis demonstrates that the HFF data indeed offer an ideal opportunity to study faint high-redshift galaxies. All lensing maps produced from our mass modeling will be made available on the STScI website.

We present a study of stellar populations in a sample of spectroscopically-confirmed Lyman-break galaxies (LBGs) and Ly$\alpha$ emitters (LAEs) at $5.7<z<7$. These galaxies have deep optical and infrared images from Subaru, $HST$, and $Spitzer$/IRAC. We focus on a subset of 27 galaxies with IRAC detections, and characterize their stellar populations utilizing galaxy synthesis models based on the multi-band data and secure redshifts. By incorporating nebular emission estimated from the observed Ly$\alpha$ flux, we are able to break the strong degeneracy of model spectra between young galaxies with prominent nebular emission and older galaxies with strong Balmer breaks. The results show that our galaxies cover a wide range of ages from several to a few hundred million years (Myr), and a wide range of stellar masses from $\sim10^8$ to $\sim10^{11}$ $M_{\odot}$. These galaxies can be roughly divided into an `old' subsample and a `young' subsample. The `old' subsample consists of galaxies older than 100 Myr, with stellar masses higher than $10^9$ $M_{\odot}$. The galaxies in the `young' subsample are younger than $\sim$30 Myr, with masses ranging between $\sim10^8$ and $\sim3\times10^9$ $M_{\odot}$. Both subsamples display a correlation between stellar mass and star-formation rate (SFR), but with very different normalizations. The average specific SFR (sSFR) of the `old' subsample is 3--4 Gyr$^{-1}$, consistent with previous studies of `normal' star-forming galaxies at $z\ge6$. The average sSFR of the `young' subsample is an order of magnitude higher, likely due to starburst activity. Our results also indicate little or no dust extinction in the majority of the galaxies, as already suggested by their steep rest-frame UV slopes. Finally, LAEs and LBGs with strong Ly$\alpha$ emission are indistinguishable in terms of age, stellar mass, and SFR.

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.

We present statistical properties of diffuse Lyman-alpha halos (LAHs) around high-$z$ star-forming galaxies with large Subaru samples of Lyman-alpha emitters (LAEs) at $z=2.2$. We make subsamples defined by the physical quantities of LAEs' central Lyman-alpha luminosities, UV magnitudes, Lyman-alpha equivalent widths, and UV slopes, and investigate LAHs' radial surface brightness (SB) profiles and scale lengths $r_n$ as a function of these physical quantities. We find that there exist prominent LAHs around LAEs with faint Lyman-alpha luminosities, bright UV luminosities, and small Lyman-alpha equivalent widths in cumulative radial Lyman-alpha SB profiles. We confirm this trend with the anti-correlation between $r_n$ and Lyman-alpha luminosities (equivalent widths) based on the Spearman's rank correlation coefficient that is $\rho=-0.9$ ($-0.7$) corresponding to the $96\%$ ($93\%$) confidence level, although the correlation between $r_n$ and UV magnitudes is not clearly found in the rank correlation coefficient. Our results suggest that LAEs with properties similar to typical Lyman-break galaxies (with faint Lyman-alpha luminosities and small equivalent widths) possess more prominent LAHs. We investigate scenarios for the major physical origins of LAHs with our results, and find that the cold stream scenario is not preferred, due to the relatively small equivalent widths up to $77$\AA\ in LAHs that include LAEs' central components. There remain two possible scenarios of Lyman-alpha scattering in circum-galactic medium and satellite galaxies that cannot be tested with our observational data.

We present results from Subaru/FMOS near-infrared (NIR) spectroscopy of 118 star-forming galaxies at $z\sim1.5$ in the Subaru Deep Field. These galaxies are selected as [OII]$\lambda$3727 emitters at $z\approx$ 1.47 and 1.62 from narrow-band imaging. We detect H$\alpha$ emission line in 115 galaxies, [OIII]$\lambda$5007 emission line in 45 galaxies, and H$\beta$, [NII]$\lambda$6584, and [SII]$\lambda\lambda$6716,6731 in 13, 16, and 6 galaxies, respectively. Including the [OII] emission line, we use the six strong nebular emission lines in the individual and composite rest-frame optical spectra to investigate physical conditions of the interstellar medium in star-forming galaxies at $z\sim$1.5. We find a tight correlation between H$\alpha$ and [OII], which suggests that [OII] can be a good star formation rate (SFR) indicator for galaxies at $z\sim1.5$. The line ratios of H$\alpha$/[OII] are consistent with those of local galaxies. We also find that [OII] emitters have strong [OIII] emission lines. The [OIII]/[OII] ratios are larger than normal star-forming galaxies in the local Universe, suggesting a higher ionization parameter. Less massive galaxies have larger [OIII]/[OII] ratios. With evidence that the electron density is consistent with local galaxies, the high ionization of galaxies at high redshifts may be attributed to a harder radiation field by a young stellar population and/or an increase in the number of ionizing photons from each massive star.

We investigate sizes of $z\sim6-8$ dropout galaxies using the complete data of the Abell 2744 cluster and parallel fields in the Hubble Frontier Fields program. By directly fitting light profiles of observed galaxies with lensing-distorted S\'ersic profiles on the image plane with the \texttt{glafic} software, we accurately measure intrinsic sizes of 31 $z\sim6-7$ and eight $z\sim8$ galaxies, including those as faint as $M_{\mathrm{UV}}\simeq-16.6$. We find that half-light radii $r_\mathrm{e}$ positively correlates with UV luminosity at each redshift, although the correlation is not very tight. Largest ($r_\mathrm{e}>0.8$ kpc) galaxies are mostly red in UV color while smallest ($r_\mathrm{e} < 0.08$ kpc) ones tend to be blue. We also find that galaxies with multiple cores tend to be brighter. Combined with previous results at $2.5\lesssim z\lesssim12$, our result confirms that the average $r_{\mathrm{e}}$ of bright ($(0.3-1)L^\star_{z=3}$) galaxies scales as $r_{\mathrm{e}}\propto(1+z)^{-m}$ with $m=1.24\pm0.1$. We find that the ratio of $r_\mathrm{e}$ to virial radius is virtually constant at $3.3\pm0.1\%$ over a wide redshift range, where the virial radii of hosting dark matter halos are derived based on the abundance matching. This constant ratio is consistent with the disk formation model by Mo et al. (1998) with $j_\mathrm{d}\sim m_\mathrm{d}$, where $j_\mathrm{d}$ and $m_\mathrm{d}$ are the fractions of the angular momentum and mass within halos confined in the disks. A comparison with various types of local galaxies indicates that our galaxies are most similar to circumnuclear star-forming regions of barred galaxies in the sense that a sizable amount of stars are forming in a very small area.

By stacking publicly available deep Spitzer/MIPS 24 $\mu$m and Herschel/PACS images for 213 $z \simeq 2.18$ Ly$\alpha$ Emitters (LAEs) in GOODS-South, we obtain a strong upper limit to the IR luminosity of typical LAEs and discuss their attenuation curve for the first time. The $3\sigma$ upper limit $L_{\rm TIR}^{3\sigma}= 1.1 \times 10^{10} L_\odot$, determined from the MIPS data providing the lowest limit, gives $IRX \equiv L_{\rm TIR}/L_{\rm UV} \leq 2.2$. Here we assume that the local calibration between the 8 $\mu$m emission and the dust SED shape and metallicity applies at high redshifts and that our LAEs have low metallicities as suggested by previous studies. The inferred escape fractions of Ly$\alpha$, $16$--$37%$, and UV continuum, $\ge 44%$, are higher than the cosmic averages at the same epoch. The SMC attenuation curve is consistent with the IRX and the UV slope $\beta = -1.4^{+0.2}_{-0.2}$ of our stacked LAE, while the Meurer's relation (Calzetti curve) predicts a 3.8 times higher IRX; we also discuss the validity of PACS-based $L_{\rm TIR}^{3\sigma}$ allowing the Meurer's relation. SED fitting using the Calzetti curve also gives a $\sim 10$ times higher SFR than from the $L_{\rm TIR}^{3\sigma}$ and $L_{\rm UV}$. With $M_{\star}=6.3^{+0.8}_{-2.0} \times10^8 \mathrm{M_{\odot}}$, our LAEs lie on a lower-mass extrapolation of the star formation main sequence at $z \sim 2$, suggesting that the majority of $z \sim 2$ LAEs are mildly star forming with relatively old ages of $\sim 200$ Myr. The faint $L_{\rm TIR}^{3\sigma}$ implies that LAEs contribute little to the faint ($\gtrsim 100 \mu$Jy) submm number counts by ALMA.

We present the discovery of one or two extremely faint z~6 quasars in 6.5 deg^2 utilizing a unique capability of the wide-field imaging of the Subaru/Suprime-Cam. The quasar selection was made in (i'-z_B) and (z_B-z_R) colors, where z_B and z_R are bandpasses with central wavelengths of 8842A and 9841A, respectively. The color selection can effectively isolate quasars at z~6 from M/L/T dwarfs without the J-band photometry down to z_R<24.0, which is 3.5 mag. deeper than SDSS. We have selected 17 promising quasar candidates. The follow-up spectroscopy for seven targets identified one apparent quasar at z=6.156 with M_1450=-23.10. We also identified one possible quasar at z=6.041 with a faint continuum of M_1450=-22.58 and a narrow Lyman-alpha emission with HWHM=427 km/s, which cannot be distinguished from Lyman-alpha emitters. We derive the quasar luminosity function at z~6 by combining our faint quasar sample with the bright quasar samples by SDSS and CFHQS. Including our data points invokes a higher number density in the faintest bin of the quasar luminosity function than the previous estimate employed. This suggests a steeper faint-end slope than lower-z, though it is yet uncertain based on a small number of spectroscopically identified faint quasars and several quasar candidates are still remain to be diagnosed. The steepening of the quasar luminosity function at the faint-end does increase the expected emission rate of the ionizing photon, however, it only changes by a factor of ~2-6. This was found to be still insufficient for the required photon budget of reionization at z~6.

We present the comprehensive analyses of faint dropout galaxies up to $z\sim10$ with the first full-depth data set of Abell 2744 lensing cluster and parallel fields observed by the Hubble Frontier Fields (HFF) program. We identify $54$ dropouts at $z\sim5-10$ in the HFF fields, and enlarge the size of $z\sim9$ galaxy sample obtained to date. Although the number of highly magnified ($\mu\sim10$) galaxies is small due to the tiny survey volume of strong lensing, our study reaches the galaxies' intrinsic luminosities comparable to the deepest-field HUDF studies. We derive UV luminosity functions with these faint dropouts, carefully evaluating the combination of observational incompleteness and lensing effects in the image plane by intensive simulations including magnification, distortion, and multiplication of images, with the evaluations of mass model dependences. Our results confirm that the faint-end slope, $\alpha$, is as steep as $-2$ at $z\sim6-8$, and strengthen the evidence of the rapid decrease of UV luminosity densities, $\rho_\mathrm{UV}$, at $z>8$ from the large $z\sim9$ sample. We examine whether the rapid $\rho_\mathrm{UV}$ decrease trend can reconcile with the large Thomson scattering optical depth, $\tau_\mathrm{e}$, measured by CMB experiments allowing a large space of free parameters such as average ionizing photon escape fraction and stellar-population dependent conversion factor. No parameter set can reproduce both the rapid $\rho_\mathrm{UV}$ decrease and the large $\tau_\mathrm{e}$. It is possible that the $\rho_\mathrm{UV}$ decrease moderates at $z\gtrsim11$, that the free parameters significantly evolve towards high-$z$, or that there exist additional sources of reionization such as X-ray binaries and faint AGNs.

We present the ultra-deep Subaru narrowband imaging survey for Lya emitters (LAEs) at $z=7.3$ in SXDS and COSMOS fields with a total integration time of 106 hours. Exploiting our new sharp bandwidth filter, NB101, installed on Suprime-Cam, we have reached $L(Lya)=2.4\times10^{42} \ erg \ s^{-1}$ ($5\sigma$) for $z=7.3$ LAEs, about 4 times deeper than previous Subaru $z \gtrsim 7$ studies, which allows us to reliably investigate evolution of Lya luminosity function (LF), for the first time, down to the luminosity limit same as those of Subaru $z=3.1-6.6$ LAE samples. Surprisingly, we only find three and four LAEs in SXDS and COSMOS fields, respectively, while one expects a total of $\sim 65$ LAEs by our survey in the case of no Lya LF evolution from $z=6.6$ to $7.3$.We identify a decrease of Lya LF from $z=6.6$ to $7.3$ at the $>90\%$ confidence level from our $z=7.3$ Lya LF.Moreover, the evolution of Lya LF is clearly accelerated at $z>6.6$ beyond the measurement uncertainties including cosmic variance. Because no such accelerated evolution of UV-continuum LF or cosmic star-formation rate (SFR) is found at $z\sim 7$, but suggested only at $z>8$ (Oesch et al. 2013, Bouwens et al. 2014), this accelerated Lya LF evolution is explained by physical mechanisms different from a pure SFR decrease but related to Lya production and escape in the process of cosmic reionization. Because a simple accelerating increase of IGM neutral hydrogen absorbing Lya would not reconcile with Thomson scattering optical depth measurements from WMAP and Planck, our findings may support new physical pictures suggested by recent theoretical studies, such as the existence of HI clumpy clouds within cosmic ionized bubbles selectively absorbing Lya and the large ionizing photon escape fraction of galaxies making weak Lya emission.

We performed a systematic spectroscopic observation of a protocluster at $z=6.01$ in the Subaru Deep Field. We took spectroscopy for all 53 $i'$-dropout galaxies down to $z'=27.09\,\mathrm{mag}$ in/around the protocluster region. From these observations, we confirmed that 28 galaxies are at $z\sim6$, of which ten are clustered in a narrow redshift range of $\Delta z<0.06$. To trace the evolution of this primordial structure, we applied the same $i'$-dropout selection and the same overdensity measurements used in the observations to a semi-analytic model built upon the Millennium Simulation. We obtain a relation between the significance of overdensities observed at $z\sim6$ and the predicted dark matter halo mass at $z=0$. This protocluster with $6\sigma$ overdensity is expected to grow into a galaxy cluster with a mass of $\sim5\times10^{14}\,\mathrm{M_\odot}$ at $z=0$. Ten galaxies within $10\,\mathrm{comoving\>Mpc}$ of the overdense region can, with more than an 80% probability, merge into a single dark matter halo by $z=0$. No significant differences appeared in UV and Ly$\alpha$ luminosities between the protocluster and field galaxies, suggesting that this protocluster is still in the early phase of cluster formation before the onset of any obvious environmental effects. However, further observations are required to study other properties, such as stellar mass, dust, and age. We do find that galaxies tend to be in close pairs in this protocluster. These pair-like subgroups will coalesce into a single halo and grow into a more massive structure. We may witness an onset of cluster formation at $z\sim6$ toward a cluster as seen in local universe.

Metal absorption systems are products of star formation. They are believed to be associated with massive star forming galaxies, which have significantly enriched their surroundings. To test this idea with high column density CIV absorption systems at z~5.7, we study the projected distribution of galaxies and characterise the environment of CIV systems in two independent quasar lines-of-sight: J103027.01+052455.0 and J113717.73+354956.9. Using wide field photometry (~80x60h$^{-1}$ comoving Mpc), we select bright (Muv(1350\AA)<-21.0 mag.) Lyman break galaxies (LBGs) at z~5.7 in a redshift slice \Delta z~0.2 and we compare their projected distribution with z~5.7 narrow-band selected Lyman alpha emitters (LAEs, \Delta z~0.08). We find that the CIV systems are located more than 10h$^{-1}$ projected comoving Mpc from the main concentrations of LBGs and no candidate is closer than ~5h$^{-1}$ projected comoving Mpc. In contrast, an excess of LAEs -lower mass galaxies- is found on scales of ~10h$^{-1}$ comoving Mpc, suggesting that LAEs are the primary candidates for the source of the CIV systems. Furthermore, the closest object to the system in the field J1030+0524 is a faint LAE at a projected distance of 212h$^{-1}$ physical kpc. However, this work cannot rule out undiscovered lower mass galaxies as the origin of these absorption systems. We conclude that, in contrast with lower redshift examples (z<3.5), strong CIV absorption systems at z~5.7 trace low-to-intermediate density environments dominated by low-mass galaxies. Moreover, the excess of LAEs associated with high levels of ionizing flux agrees with the idea that faint galaxies dominate the ionizing photon budget at this redshift.

We present diffuse Lyman-alpha halos (LAHs) identified in the composite Subaru narrowband images of 100-3600 Lyman-alpha emitters (LAEs) at z=2.2, 3.1, 3.7, 5.7, and 6.6. First, we carefully examine potential artifacts mimicking LAHs that include a large-scale point-spread function (PSF) made by instrumental and atmospheric effects. Based on our critical test with composite images of non-LAE samples whose narrowband-magnitude and source-size distributions are the same as our LAE samples, we confirm that no artifacts can produce a diffuse extended feature similar to our LAHs. After this test, we measure the scale lengths of exponential profile for the LAHs estimated from our z=2.2-6.6 LAE samples of L(Lyman-alpha) > 2 x 10^42 erg s^-1. We obtain the scale lengths of ~ 5-10 kpc at z=2.2-5.7, and find no evolution of scale lengths in this redshift range beyond our measurement uncertainties. Combining this result and the previously-known UV-continuum size evolution, we infer that the ratio of LAH to UV-continuum sizes is nearly constant at z=2.2-5.7. The scale length of our z=6.6 LAH is larger than 5-10 kpc just beyond the error bar, which is a hint that the scale lengths of LAHs would increase from z=5.7 to 6.6. If this increase is confirmed by future large surveys with significant improvements of statistical and systematical errors, this scale length change at z > 6 would be a signature of increasing fraction of neutral hydrogen scattering Lyman-alpha photons, due to cosmic reionization.

We present a statistical study of velocities of Lya, interstellar (IS) absorption, and nebular lines and gas covering fraction for Lya emitters (LAEs) at z~2. We make a sample of 22 LAEs with a large Lya equivalent width (EW) of > 50A based on our deep Keck/LRIS observations, in conjunction with spectroscopic data from the Subaru/FMOS program and the literature. We estimate the average velocity offset of Lya from a systemic redshift determined with nebular lines to be dv_Lya=234+-9 km s-1. Using a Kolmogorv-Smirnov test, we confirm the previous claim of Hashimoto et al. (2013) that the average dv_Lya of LAEs is smaller than that of LBGs. Our LRIS data successfully identify blue-shifted multiple IS absorption lines in the UV continua of four LAEs on an individual basis. The average velocity offset of IS absorption lines from a systemic redshift is dv_IS=204+-27 km s-1, indicating LAE's gas outflow with a velocity comparable to typical LBGs. Thus, the ratio, R^Lya_ IS = dv_Lya/dv_IS of LAEs, is around unity, suggestive of low impacts on Lya transmission by resonant scattering of neutral hydrogen in the IS medium. We find an anti-correlation between Lya EW and the covering fraction, f_c, estimated from the depth of absorption lines, where f_c is an indicator of average neutral hydrogen column density, N_HI. The results of our study support the idea that N_HI is a key quantity determining Lya emissivity.

We present the results of structure analyses for a large sample of 426 Lya emitters (LAEs) at z~2.2 that are observed with HST/ACS and WFC3-IR by deep extra-galactic legacy surveys. We confirm that the merger fraction and the average ellipticity of LAE's stellar component are 10-30 % and 0.4-0.6, respectively, that are comparable with previous study results. We successfully identify that some LAEs have a spatial offset between Lya and stellar-continuum emission peaks, d_Lya, by ~2.5-4 kpc beyond our statistical errors. To uncover the physical origin of strong Lya emission found in LAEs, we investigate Lya equivalent width (EW) dependences of these three structural parameters, merger fraction, d_Lya, and ellipticity of stellar distribution in the range of EW(Lya)=20-250A. Contrary to expectations, we find that merger fraction does not significantly increase with Lya EW. We reveal an anti-correlation between d_Lya and EW(Lya) by Kolmogorov-Smirnov (KS) test. There is a trend that the LAEs with a large Lya EW have a small ellipticity. This is consistent with the recent theoretical claims that Lya photons can more easily escape from face-on disks having a small ellipticity, due to less inter-stellar gas along the line of sight, although our KS test indicates that this trend is not statistically significant. Our results of Lya-EW dependence generally support the idea that an HI column density is a key quantity determining Lya emissivity.

We present the first systematic search for extended metal-line [OII]{\lambda}{\lambda}3726,3729 nebulae, or [OII] blobs (OIIBs), at z=1.2 using deep narrowband imaging with a survey volume of 1.9x10^5 Mpc^3 on the 0.62 deg^2 sky of Subaru-XMM Deep Survey (SXDS) field. We discover a giant OIIB, dubbed 'OIIB 1', with a spatial extent over ~75 kpc at a spectroscopic redshift of z=1.18, and also identify a total of twelve OIIBs with a size of >30 kpc. Our optical spectrum of OIIB 1 presents [NeV]{\lambda}3426 line at the 6{\sigma} level, indicating that this object harbors an obscured type-2 AGN. The presence of gas outflows in this object is suggested by two marginal detections of FeII{\lambda}2587 absorption and FeII*{\lambda}2613 emission lines both of which are blueshifted at as large as 500-600 km/s, indicating that the heating source of OIIB 1 is AGN or associated shock excitation rather than supernovae produced by starbursts. The number density of OIIB 1-type giant blobs is estimated to be ~5x10^{-6} Mpc^{-3} at z~1.2, which is comparable with that of AGNs driving outflow at a similar redshift, suggesting that giant OIIBs are produced only by AGN activity. On the other hand, the number density of small OIIBs, 6x10^{-5} Mpc^{-3}, compared to that of z~1 galaxies in the blue cloud in the same M_B range, may imply that 3% of star-forming galaxies at z~1 are quenching star formation through outflows involving extended [OII] emission.

Using deep narrow-band (NB) imaging and optical spectroscopy from the Keck telescope and MMT, we identify a sample of 20 emission-line galaxies (ELGs) at z=0.065-0.90 where the weak auroral emission line, [OIII]4363, is detected at >3\sigma. These detections allow us to determine the gas-phase metallicity using the "direct'' method. With electron temperature measurements and dust attenuation corrections from Balmer decrements, we find that 4 of these low-mass galaxies are extremely metal-poor with 12+log(O/H) <= 7.65 or one-tenth solar. Our most metal-deficient galaxy has 12+log(O/H) = 7.24^{+0.45}_{-0.30} (95% confidence), similar to some of the lowest metallicity galaxies identified in the local universe. We find that our galaxies are all undergoing significant star formation with average specific star formation rate (SFR) of (100 Myr)^{-1}, and that they have high central SFR surface densities (average of 0.5 Msun/yr/kpc^2. In addition, more than two-thirds of our galaxies have between one and four nearby companions within a projected radius of 100 kpc, which we find is an excess among star-forming galaxies at z=0.4-0.85. We also find that the gas-phase metallicities for a given stellar mass and SFR lie systematically below the local M-Z-(SFR) relation by \approx0.2 dex (2\sigma\ significance). These results are partly due to selection effects, since galaxies with strong star formation and low metallicity are more likely to yield [OIII]4363 detections. Finally, the observed higher ionization parameter and electron density suggest that they are lower redshift analogs to typical z>1 galaxies.

We report deep ALMA observations complemented with associated HST imaging for a luminous (m_uv=25) galaxy, `Himiko', at a redshift z=6.595. The galaxy is remarkable for its high star formation rate, 100 Mo/yr, securely estimated from our deep HST and Spitzer photometry, and the absence of any evidence for strong AGN activity or gravitational lensing magnification. Our ALMA observations probe an order of magnitude deeper than previous IRAM observations, yet fail to detect a 1.2mm dust continuum, indicating a flux <52uJy comparable with or weaker than that of local dwarf irregulars with much lower star formation rates. We likewise provide a strong upper limit for the flux of [CII] 158um, L([CII]) < 5.4x10^7 Lo, a diagnostic of the hot interstellar gas often described as a valuable probe for early galaxies. In fact, our observations indicate Himiko lies off the local L([CII]) - star formation rate scaling relation by a factor of more than 30. Both aspects of our ALMA observations suggest Himiko is an unique object with a very low dust content and perhaps nearly primordial interstellar gas. Our HST images provide unique insight into the morphology of this remarkable source, highlighting an extremely blue core of activity and two less extreme associated clumps. Himiko is undergoing a triple major merger event whose extensive ionized nebula of Lyman alpha emitting gas, discovered in our earlier work with Subaru, is powered by star formation and the dense circum-galactic gas. We are likely witnessing an early massive galaxy during a key period of its mass assembly close to the end of the reionization era.

We present a detailed structural and morphological study of a large sample of spectroscopically-confirmed galaxies at z >= 6, using deep HST near-IR broad-band images and Subaru optical narrow-band images. The galaxy sample consists of 51 Lyman-alpha emitters (LAEs) at z ~ 5.7, 6.5, and 7.0, and 16 Lyman-break galaxies (LBGs) at 5.9 < z < 6.5. These galaxies exhibit a wide range of rest-frame UV continuum morphology in the HST images, from compact features to multiple component systems. The fraction of merging/interacting galaxies reaches 40% ~ 50% at the brightest end of M_1500 <= -20.5 mag. The intrinsic half-light radii r_{hl,in}, after correction for PSF broadening, are roughly between r_{hl,in} ~ 0.05" (0.3 kpc) and 0.3" (1.7 kpc) at M_1500 <= -19.5 mag. The median r_{hl,in} value is 0.16" (~0.9 kpc). This is consistent with the sizes of bright LAEs and LBGs at z >= 6 in previous studies. In addition, more luminous galaxies tend to have larger sizes, exhibiting a weak size-luminosity relation r_{hl,in} \propto L^{0.14} at M_1500 <= -19.5 mag. The slope of 0.14 is significantly flatter than those in fainter LBG samples. We discuss the morphology of z >= 6 galaxies with nonparametric methods, including the CAS system and the Gini and M_20 parameters, and demonstrate their validity through simulations. We search for extended Lyman-alpha emission halos around LAEs at z ~ 5.7 and 6.5, by stacking a number of narrow-band images. We do not find evidence of extended halos predicted by cosmological simulations. Such Lyman-alpha halos, if they exist, could be weaker than predicted. Finally, we investigate any positional misalignment between UV continuum and Lyman-alpha emission in LAEs. While the two positions are generally consistent, several merging galaxies show significant positional differences. This is likely caused by a disturbed ISM distribution due to merging activity.

We present deep HST near-IR and Spitzer mid-IR observations of a large sample of spectroscopically-confirmed galaxies at z >= 6. The sample consists of 51 Lyman-alpha emitters (LAEs) at z ~ 5.7, 6.5, and 7.0, and 16 Lyman-break galaxies (LBGs) at 5.9 < z < 6.5. The near-IR images were mostly obtained with WFC3 in the F125W and F160W bands, and the mid-IR images were obtained with IRAC in the 3.6um and 4.5um bands. Our galaxies also have deep optical imaging data from Subaru Suprime-Cam. We utilize the multi-band data and secure redshifts to derive their rest-frame UV properties. These galaxies have steep UV continuum slopes roughly between beta ~ -1.5 and -3.5, with an average value of beta ~ -2.3, slightly steeper than the slopes of LBGs in previous studies. The slope shows little dependence on UV continuum luminosity except for a few of the brightest galaxies. We find a statistically significant excess of galaxies with slopes around beta ~ -3, suggesting the existence of very young stellar populations with extremely low metallicity and dust content. Our galaxies have moderately strong rest-frame Lyman-alpha equivalent width (EW) in a range of ~10 to ~200 \AA. The star-formation rates are also moderate, from a few to a few tens solar masses per year. The LAEs and LBGs in this sample share many common properties, implying that LAEs represent a subset of LBGs with strong Lyman-alpha emission. Finally, the comparison of the UV luminosity functions between LAEs and LBGs suggests that there exists a substantial population of faint galaxies with weak Lyman-alpha emission (EW < 20 \AA) that could be the dominant contribution to the total ionizing flux at z >= 6.

Based on optical and X-ray data for a sample of 34 relaxed rich clusters of galaxies with redshifts of 0.1-0.9, we studied relative spatial distributions of the two major baryon contents, the cluster galaxies and the hot plasmas. Using multi-band photometric data taken with the UH88 telescope, we determined the integrated (two dimensional) radial light profiles of member galaxies in each cluster using two independent approaches, i.e., the background subtraction and the color-magnitude filtering. The ICM mass profile of each cluster in our sample, also integrated in two dimensions, was derived from a spatially-resolved spectral analysis using XMM-Newton and Chandra data. Then, the radially-integrated light profile of each cluster was divided by its ICM mass profile, to obtain a profile of "galaxy light vs. ICM mass ratio". The ratio profiles over the central 0.65 R500 regions were found to steepen from the higher- to lower- redshift subsamples, meaning that the galaxies become more concentrated in the ICM sphere towards lower redshifts. The evolution is also seen in galaxy number vs. ICM mass ratio profiles. A range of systematic uncertainties in the galaxy light measurements, as well as many radius-/redshift- dependent biases to the galaxy vs. ICM profiles have been assessed, but none of them is significant against the observed evolution. Besides, the galaxy light vs. total mass ratio profiles also exhibit gradual concentration towards lower redshift. We interpret in the context that the galaxies, the ICM, and the dark matter components followed a similar spatial distribution in the early phase (z>0.5), while the galaxies have fallen towards the center relative to the others at a later phase.

We present results from Keck/NIRSPEC and Magellan/MMIRS follow-up spectroscopy of Lya emitters (LAEs) at z=2.2 identified in our Subaru narrowband survey. We successfully detect Ha emission from seven LAEs, and perform a detailed analysis of six LAEs free from AGN activity, two out of which, CDFS-3865 and COSMOS-30679, have [OII] and [OIII] line detections. They are the first [OII]-detected LAEs at high-z, and their [OIII]/[OII] ratios and R23-indices provide the first simultaneous determinations of ionization parameter and oxygen abundance for LAEs. CDFS-3865 has a very high ionization parameter (q_{ion}=2.5^{+1.7}_{-0.8}x10^8 cm s^{-1}) and a low oxygen abundance (12+log(O/H)=7.84^{+0.24}_{-0.25}) in contrast with moderate values of other high-z galaxies such as Lyman-break galaxies (LBGs). COSMOS-30679 also possesses a relatively high ionization parameter (q_{ion}=8^{+10}_{-4}x10^7 cm s^{-1}) and a low oxygen abundance (12+log(O/H)=8.18^{+0.28}_{-0.28}). Both LAEs appear to fall below the mass-metallicity relation of z~2 LBGs. Similarly, a low metallicity of 12+log(O/H)<8.4 is independently indicated for typical LAEs from a composite spectrum and the [NII]/Ha index. Such high ionization parameters and low oxygen abundances can be found in local star-forming galaxies, but this extreme local population occupies only ~0.06% of the SDSS spectroscopic galaxy sample with a number density ~100 times smaller than that of LAEs. With their high ionization parameters and low oxygen abundances, LAEs would represent an early stage of galaxy formation dominated by massive stars in compact star-forming regions. High-q_{ion} galaxies like LAEs would produce ionizing photons efficiently with a high escape fraction achieved by density-bounded HII regions, which would significantly contribute to cosmic reionization at z>6.

We present the results of Magellan/MMIRS and Keck/NIRSPEC spectroscopy for five Lya emitters (LAEs) at z=2.2 for which high-resolution FUV spectra from Magellan/MagE are available. We detect nebular emission lines including Ha on the individual basis and low-ionization interstellar (LIS) absorption lines in a stacked FUV spectrum, and measure average offset velocities of the Lya line, Delta_v_Lya, and LIS absorption lines, Delta_v_abs, with respect to the systemic velocity defined by the nebular lines. For a sample of eight z~2-3 LAEs without AGN from our study and the literature, we obtain Delta_v_Lya = 175+/-35 km s^{-1}, which is significantly smaller than that of Lyman-break Galaxies (LBGs), Delta_v_Lya=400 km s^{-1}. The stacked FUV spectrum gives Delta_v_abs = -179 +/- 73 km s^{-1}, comparable to that of LBGs. These positive Delta_v_Lya and negative Delta_v_abs suggest that LAEs also have outflows. In contrast to LBGs, however, the LAEs' Delta_v_Lya is as small as |Delta_v_abs|, suggesting low neutral hydrogen column densities. Such a low column density with a small number of resonant scattering may cause the observed strong Lya emission of LAEs. We find an anti-correlation between Lya equivalent width (EW) and Delta_v_Lya in a compilation of LAE and LBG samples. Although its physical origin is not clear, this anti-correlation result appears to challenge the hypothesis that a strong outflow, by means of a reduced number of resonant scattering, produces a large EW. If LAEs at z>6 have similarly small Delta_v_Lya values, constraints on the reionization history derived from the Lya transmissivity may need to be revised.

We analyze the redshift- and luminosity-dependent sizes of dropout galaxy candidates in the redshift range z~7-12 using deep images from the UDF12 campaign, data which offers two distinct advantages over that used in earlier work. Firstly, we utilize the increased S/N ratio offered by the UDF12 imaging to provide improved size measurements for known galaxies at z=6.5-8 in the HUDF. Specifically, we stack the new deep F140W image with the existing F125W data in order to provide improved measurements of the half-light radii of z-dropouts. Similarly we stack this image with the new deep UDF12 F160W image to obtain new size measurements for a sample of Y-dropouts. Secondly, because the UDF12 data have allowed the construction of the first robust galaxy sample in the HUDF at z>8, we have been able to extend the measurement of average galaxy size out to significantly higher redshifts. Restricting our size measurements to sources which are now detected at >15sigma, we confirm earlier indications that the average half-light radii of z~7-12 galaxies are extremely small, 0.3-0.4 kpc, comparable to the sizes of giant molecular associations in local star-forming galaxies. We also confirm that there is a clear trend of decreasing half-light radius with increasing redshift, and provide the first evidence that this trend continues beyond z~8. Modeling the evolution of the average half-light radius as a power-law (1+z)^s, we obtain a best-fit index of s=-1.28+/-0.13 over the redshift range z~4-12, mid-way between the physically expected evolution for baryons embedded in dark halos of constant mass (s=-1) and constant velocity (s=-1.5). A clear size-luminosity relation, such as that found at lower redshift, is also evident in both our z- and Y-dropout sample. This relation can be interpreted in terms of a constant surface density of star formation over a range in luminosity of 0.05-1.0L*_z=3.(abridged)

We have identified a very interesting Ly-alpha emitter, whose Ly-alpha emission line has an extremely large observed equivalent width of EW_0=436^{+422}_{-149}A, which corresponds to an extraordinarily large intrinsic rest-frame equivalent width of EW_0^{int}=872^{+844}_{-298}A after the average intergalactic absorption correction. The object was spectroscopically confirmed to be a real Ly-alpha emitter by its apparent asymmetric Ly-alpha line profile detected at z=6.538. The continuum emission of the object was definitely detected in our deep z'-band image; thus, its EW_0 was reliably determined. Follow-up deep near-infrared spectroscopy revealed emission lines of neither He II lambda1640 as an apparent signature of Population III, nor C IV lambda1549 as a proof of active nucleus. No detection of short-lived He II lambda1640 line is not necessarily inconsistent with the interpretation that the underlying stellar population of the object is dominated by Population III. We found that the observed extremely large EW_0 of the Ly-alpha emission and the upper limit on the EW_0 of the He II lambda1640 emission can be explained by population synthesis models favoring a very young age less than 2-4Myr and massive metal-poor (Z<10^{-5}) or even metal-free stars. The observed large EW_0 of Ly-alpha is hardly explained by Population I/II synthesis models with Z>10^{-3}. However, we cannot conclusively rule out the possibility that this object is composed of a normal stellar population with a clumpy dust distribution, which could enhance the Ly-alpha EW_0, though its significance is still unclear.

We present the first detailed study of the stellar populations of star-forming galaxies at z~1.5, which are selected by their [O II] emission line, detected in narrow-band surveys. We identified ~1,300 [O II] emitters at z=1.47 and z=1.62 in the Subaru Deep Field with rest-frame EWs above 13\AA. Optical and near-infrared spectroscopic observations for ~10% of our samples show that our separation of [O II] from [O III] emission-line galaxies in two-color space is 99% successful. We analyze the multi-wavelength properties of a subset of ~1,200 galaxies with the best photometry. They have average rest-frame EW of 45\AA, stellar mass of 3 x 10^9 M_sun, and stellar age of 100 Myr. In addition, our SED fitting and broad-band colors indicate that [O II] emitters span the full range of galaxy populations at z~1.5. We also find that 80% of [O II] emitters are also photometrically classified as "BX/BM" (UV) galaxies and/or the star-forming "BzK" (near-IR) galaxies. Our [O II] emission line survey produces a far more complete, and somewhat deeper sample of z~1.5 galaxies than either the BX/BM or sBzK selection alone. We constructed average SEDs and find that higher [O II] EW galaxies have somewhat bluer continua. SED model-fitting shows that they have on average half the stellar mass of galaxies with lower [O II] EW. The observed [O II] luminosity is well-correlated with the far-UV continuum with a logarithmic slope slightly 0f 0.89\pm0.22. The scatter of the [O II] luminosity against the far-UV continuum suggests that [O II] can be used as a SFR indicator with a reliability of 0.23 dex.

We have performed deep imaging surveys for LyA emitters (LAEs) at redshift ~7.3 in two blank fields, the Subaru Deep Field (SDF) and the Subaru/XMM-Newton Deep survey Field (SXDF), using the Subaru/Suprime-Cam equipped with new red-sensitive CCDs and a new narrow-band filter, NB1006 (lambda_c=10052 Ang, FWHM=214 Ang). We identified four objects as LAE candidates that exhibit luminosity excess in NB1006. By carrying out deep follow-up spectroscopy for three of them using Subaru/FOCAS and Keck/DEIMOS, a definitively asymmetric emission line is detected for one of them, SXDF-NB1006-2. Assuming this line is LyA, this object is a LAE at z=7.215 which has luminosity of 1.2^{+1.5}_{-0.6} x 10^43 [erg s-1] and a weighted skewness S_w=4.90+-0.86. Another object, SDF-NB1006-2, shows variable photometry and is thus probably a quasar (QSO) or an active galactic nucleus (AGN). It shows an asymmetric emission line at 10076 Ang, which may be due to either LyA at z=7.288 or [OII] at z=1.703. The third object, SDF-NB1006-1, is likely a galaxy with temporal luminosity enhancement associated with a supernova explosion, as the brightness of this object varies between the observed epochs. Its spectrum does not show any emission lines. The inferred decrease in the number density of LAEs toward higher redshift is n_LyA(z=7.3)/n_LyA(z=5.7) = 0.05^+0.11_-0.05 from z=5.7 to 7.3 down to L(LyA)=1.0 x 10^43 [erg s-1]. The present result is consistent with the interpretation in previous studies that the neutral hydrogen fraction is rapidly increasing from z=5.7 to 7.3.

We report the discovery of a protocluster at z~6 containing at least eight cluster member galaxies with spectroscopic confirmations in the wide-field image of the Subaru Deep Field (SDF). The overdensity of the protocluster is significant at the 6 sigma level, based on the surface number density of i'-dropout galaxies. The overdense region covers ~36 sq. arcmin, and includes 30 i'-dropout galaxies. Follow-up spectroscopy revealed that 15 of these are real z~6 galaxies (5.7 < z < 6.3). Eight of the 15 are clustering in a narrow redshift range centered at z=6.01, corresponding to a seven-fold increase in number density over the average in redshift space. We found no significant difference in the observed properties, such as Ly-alpha luminosities and UV continuum magnitudes, between the eight protocluster members and the seven non-members. The velocity dispersion of the eight protocluster members is 647 km/s, which is about three times higher than that predicted by the standard cold dark matter model. This discrepancy could be attributed to the distinguishing three-dimensional distribution of the eight protocluster members. We discuss two possible explanations for this discrepancy: either the protocluster is already mature, with old galaxies at the center, or it is still immature and composed of three subgroups merging to become a larger cluster. In either case, this concentration of z=6.01 galaxies in the SDF may be one of the first sites of formation of a galaxy cluster in the universe.

Using deep narrow-band and broad-band imaging, we identify 401 z~0.40 and 249 z~0.49 H-alpha line-emitting galaxies in the Subaru Deep Field. Compared to other H-alpha surveys at similar redshifts, our samples are unique since they probe lower H-alpha luminosities, are augmented with multi-wavelength (rest-frame 1000AA--1.5 microns) coverage, and a large fraction (20%) of our samples has already been spectroscopically confirmed. Our spectra allow us to measure the Balmer decrement for nearly 60 galaxies with H-beta detected above 5-sigma. The Balmer decrements indicate an average extinction of A(H-alpha)=0.7^{+1.4}_{-0.7} mag. We find that the Balmer decrement systematically increases with higher H-alpha luminosities and with larger stellar masses, in agreement with previous studies with sparser samples. We find that the SFRs estimated from modeling the spectral energy distribution (SED) is reliable---we derived an "intrinsic" H-alpha luminosity which is then reddened assuming the color excess from SED modeling. The SED-predicted H-alpha luminosity agrees with H-alpha narrow-band measurements over 3 dex (rms of 0.25 dex). We then use the SED SFRs to test different statistically-based dust corrections for H-alpha and find that adopting one magnitude of extinction is inappropriate: galaxies with lower luminosities are less reddened. We find that the luminosity-dependent dust correction of Hopkins et al. yields consistent results over 3 dex (rms of 0.3 dex). Our comparisons are only possible by assuming that stellar reddening is roughly half of nebular reddening. The strong correspondence argue that with SED modeling, we can derive reliable intrinsic SFRs even in the absence of H-alpha measurements at z~0.5.

We present the average metallicity and star-formation rate of Lya emitters (LAEs) measured from our large-area survey with three narrow-band (NB) filters covering the Lya, [OII]3727, and Ha+[NII] lines of LAEs at z=2.2. We select 919 z=2.2 LAEs from Subaru/Suprime-Cam NB data in conjunction with Magellan/IMACS spectroscopy. Of these LAEs, 561 and 105 are observed with KPNO/NEWFIRM near-infrared NB filters whose central wavelengths are matched to redshifted [OII] and Ha nebular lines, respectively. By stacking the near-infrared images of the LAEs, we successfully obtain average nebular-line fluxes of LAEs, the majority of which are too faint to be identified individually by narrow-band imaging or deep spectroscopy. The stacked object has an Ha luminosity of 1.7x10^{42} erg s^{-1} corresponding to a star formation rate (SFR) of 14 M_{sun} yr^{-1}. We place, for the first time, a firm lower limit to the average metallicity of LAEs of Z>~0.09 Z_{sun} (2sigma) based on the [OII]/(Ha+[NII]) index together with photo-ionization models and empirical relations. This lower limit of metallicity rules out the hypothesis that LAEs, so far observed at z~2, are extremely metal poor (Z<2x10^{-2} Z_{sun}) young galaxies at the 4sigma level. This limit is higher than a simple extrapolation of the observed mass-metallicity relation of z~2 UV-selected galaxies toward lower masses (5x10^{8} M_{sun}), but roughly consistent with a recently proposed fundamental mass-metallicity relation when the LAEs' relatively low SFR is taken into account. The Ha and Lya luminosities of our NB-selected LAEs indicate that the escape fraction of Lya photons is ~12-30 %, much higher than the values derived for other galaxy populations at z~2.

We present the results of our ultra-deep Keck/DEIMOS spectroscopy of z-dropout galaxies in the SDF and GOODS-N. For 3 out of 11 objects, we detect an emission line at ~ 1um with a signal-to-noise ratio of ~ 10. The lines show asymmetric profiles with high weighted skewness values, consistent with being Lya, yielding redshifts of z=7.213, 6.965, and 6.844. Specifically, we confirm the z=7.213 object in two independent DEIMOS runs with different spectroscopic configurations. The z=6.965 object is a known Lya emitter, IOK-1, for which our improved spectrum at a higher resolution yields a robust skewness measurement. The three z-dropouts have Lya fluxes of 3 x 10^-17 erg s^-1 cm^-2 and rest-frame equivalent widths EW_0^Lya = 33-43A. Based on the largest spectroscopic sample of 43 z-dropouts that is the combination of our and previous data, we find that the fraction of Lya-emitting galaxies (EW_0^Lya > 25A) is low at z ~ 7; 17 +- 10% and 24 +- 12% for bright (Muv ~= -21) and faint (Muv ~= -19.5) galaxies, respectively. The fractions of Lya-emitting galaxies drop from z ~ 6 to 7 and the amplitude of the drop is larger for faint galaxies than for bright galaxies. These two pieces of evidence would indicate that the neutral hydrogen fraction of the IGM increases from z ~ 6 to 7, and that the reionization proceeds from high- to low-density environments, as suggested by an inside-out reionization model.

We report the detection of an extremely bright ($\sim$37 mJy at 1100 $\mu$m and $\sim$91 mJy at 880 $\mu$m) submillimeter galaxy (SMG), AzTEC-ASTE-SXDF1100.001 (hereafter referred to as SXDF1100.001 or Orochi), discovered in 1100 $\mu$m observations of the Subaru/XMM-Newton Deep Field using AzTEC on ASTE. Subsequent CARMA 1300 $\mu$m and SMA 880 $\mu$m observations successfully pinpoint the location of Orochi and suggest that it has two components, one extended (FWHM of $\sim$ 4$^{\prime\prime}$) and one compact (unresolved). Z-Spec on CSO has also been used to obtain a wide band spectrum from 190 to 308 GHz, although no significant emission/absorption lines are found. The derived upper limit to the line-to-continuum flux ratio is 0.1--0.3 (2 $\sigma$) across the Z-Spec band. Based on the analysis of the derived spectral energy distribution from optical to radio wavelengths of possible counterparts near the SMA/CARMA peak position, we suggest that Orochi is a lensed, optically dark SMG lying at $z \sim 3.4$ behind a foreground, optically visible (but red) galaxy at $z \sim 1.4$. The deduced apparent (i.e., no correction for magnification) infrared luminosity ($L_{\rm IR}$) and star formation rate (SFR) are $6 \times 10^{13}$ $L_{\odot}$ and 11000 $M_{\odot}$ yr$^{-1}$, respectively, assuming that the $L_{\rm IR}$ is dominated by star formation. These values suggest that Orochi will consume its gas reservoir within a short time scale ($3 \times 10^{7}$ yr), which is indeed comparable to those in extreme starbursts like the centres of local ULIRGs.

We present Keck spectroscopic observations of z>6 Lyman-break galaxy (LBG) candidates in the Subaru Deep Field (SDF). The candidates were selected as i'-dropout objects down to z'=27 AB magnitudes from an ultra-deep SDF z'-band image. With the Keck spectroscopy we identified 19 LBGs with prominent Ly_alpha emission lines at 6< z < 6.4. The median value of the Ly_alpha rest-frame equivalent widths (EWs) is ~50 A, with four EWs >100 A. This well-defined spectroscopic sample spans a UV-continuum luminosity range of -21.8< M_{UV}<-19.5 (0.6~5 L*_{UV}) and a Ly_alpha luminosity range of (0.3~3) x 10^{43} erg s^{-1} (0.3~3 L*_ {Ly_alpha}). We derive the UV and Ly_alpha luminosity functions (LFs) from our sample at <z>~6.2 after we correct for sample incompleteness. We find that our measurement of the UV LF is consistent with the results of previous studies based on photometric LBG samples at 5<z<7. Our Ly_alpha LF is also generally in agreement with the results of Ly_alpha-emitter surveys at z~5.7 and 6.6. This study shows that deep spectroscopic observations of LBGs can provide unique constraints on both the UV and Ly_alpha LFs at z>6.

Several UV and near-infrared color selection methods have identified galaxies at z = 1-3. Since each method suffers from selection biases, we have applied three leading techniques (Lyman break, BX/BM, and BzK selection) simultaneously in the Subaru Deep Field. This field has reliable ({\Delta}z/(1 + z) = 0.02--0.09) photometric redshifts for ~53,000 galaxies from 20 bands (1500{\AA}--2.2{\mu}m). The BzK, LBG, and BX/BM samples suffer contamination from z<1 interlopers of 6%, 8%, and 20%, respectively. Around the redshifts where it is most sensitive (z~1.9 for star-forming BzK, z~1.8 for z~2 LBGs, z~1.6 for BM, and z~2.3 for BX), each technique finds 60-80% of the census of the three methods. In addition, each of the color techniques shares 75-96% of its galaxies with another method, which is consistent with previous studies that adopt identical criteria on magnitudes and colors. Combining the three samples gives a comprehensive census that includes ~90% of z-phot = 1-3 galaxies, using standard magnitude limits similar to previous studies. In fact, we find that among z = 1-2.5 galaxies in the color selection census, 81-90% of them can be selected by just combining the BzK selection with one of the UV techniques (z~2 LBG or BX and BM). The average galaxy stellar mass, reddening and SFRs all decrease systematically from the sBzK population to the LBGs, and to the BX/BMs. The combined color selections yield a total cosmic SFR density of 0.18 $\pm$ 0.03 M_sun yr^{-1} Mpc^{-3} for K_AB <= 24. We find that 65% of the star formation is in galaxies with E(B-V) > 0.25 mag, even though they are only one-fourth of the census by number.

We carried out extended spectroscopic confirmations of Ly-alpha emitters (LAEs) at z=6.5 and 5.7 in the Subaru Deep Field. Now, the total number of spectroscopically confirmed LAEs is 45 and 54 at z=6.5 and 5.7, respectively, and at least 81% (70%) of our photometric candidates at z=6.5 (5.7) have been spectroscopically identified as real LAEs. We made careful measurements of the Ly-alpha luminosity, both photometrically and spectroscopically, to accurately determine the Ly-alpha and rest-UV luminosity functions (LFs). The substantially improved evaluation of the Ly-alpha LF at z=6.5 shows an apparent deficit from z=5.7 at least at the bright end, and a possible decline even at the faint end, though small uncertainties remain. The rest-UV LFs at z=6.5 and 5.7 are in good agreement, at least at the bright end, in clear contrast to the differences seen in the Ly-alpha LF. These results imply an increase in the neutral fraction of the intergalactic medium from z=5.7 to 6.5. The rest-frame equivalent width (EW_0) distribution at z=6.5 seems to be systematically smaller than z=5.7, and it shows an extended tail toward larger EW_0. The bright end of the rest-UV LF can be reproduced from the observed Ly-alpha LF and a reasonable EW_0-UV luminosity relation. Integrating this rest-UV LF provides the first measurement of the contribution of LAEs to the photon budget required for reionization. The derived UV LF suggests that the fractional contribution of LAEs to the photon budget among Lyman break galaxies significantly increases towards faint magnitudes. Low-luminosity LAEs could dominate the ionizing photon budget, though this inference depends strongly on the uncertain faint-end slope of the Ly-alpha LF.

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

We present a wide-field H-alpha imaging survey of the rich cluster CL0939+4713 (Abell 851) at z=0.41 with Suprime-Cam on the Subaru Telescope, using the narrow-band filter NB921. The survey is sensitive to active galaxies with star formation rates down to ~0.3Msun/yr throughout the 27'x27' field, and we identified 445 H-alpha emitters along the large-scale structures around the cluster. Using this sample, we find that (1) the fraction of H-alpha emitters is a strong function of environment and shows a clear decline toward the cluster central region. We also find that (2) the color of H-alpha emitters is clearly dependent on environment. The majority of the H-alpha emitters have blue colors with B-I<2, but we find H-alpha emitters with red colors as well. Such red emitters are very rare in the cluster center or its immediate surrounding regions, while they are most frequently found in groups located far away from the cluster center. These groups coincide with the environment where a sharp transition in galaxy color distribution is seen. This may suggest that dusty star formation activity tends to be involved in galaxy truncation processes that are effective in groups, and it is probably related to the "pre-processing" that generates present-day cluster S0 galaxies. Finally, we confirm that (3) the mass-normalized integrated star formation rate in clusters (i.e. the total star formation rate within 0.5xR200 from the cluster center divided by the cluster dynamical mass) rapidly increases with look-back time following approximately ~(1+z)^6, and it is also correlated with the cluster mass.

We investigate mass-dependent galaxy evolution based on a large sample of (more than 50,000) K-band selected galaxies in a multi-wavelength catalog of the Subaru/XMM-Newton Deep Survey (SXDS) and the UKIRT Infrared Deep Sky Survey (UKIDSS)/Ultra Deep Survey (UDS). We employ the optical to near-infrared photometry to determine photometric redshifts of these galaxies. Then, we estimate the stellar mass of our sample galaxies using a standard fitting procedure. From the sample galaxies, we obtain the stellar mass function of galaxies and the cosmic stellar mass density up to z<4. Our results are consistent with previous studies and we find a considerable number of low-mass galaxies (M<10^{10.5}) at the redshift range 3<z<4. The stellar-mass dependent correlation functions of our sample galaxies show clear evolution and they connect to that in the local universe consistently. Also, the massive galaxies show strong clustering throughout our studied redshift range. The correlation length of massive galaxies rapidly decreases from z=4 to 2. We also find some high mass density regions of massive galaxies at 1.4<z<2.5 in our sample, which may be candidate progenitors of the present-day clusters of galaxies. At this redshift range, massive star-forming galaxies are the dominant population making up the structures and the passively evolving galaxies show stronger clustering and they may have formed earlier than those star-forming galaxies.

We investigate the stellar populations of Lyman alpha emitters (LAEs) at z=5.7 and 6.6 in a 0.65 deg^2 sky of the Subaru/XMM-Newton Deep Survey (SXDS) Field, using deep images taken with Subaru/Suprime-Cam, UKIRT/WFCAM, and Spitzer/IRAC. We produce stacked multiband images at each redshift from 165 (z=5.7) and 91 (z=6.6) IRAC-undetected objects, to derive typical spectral energy distributions (SEDs) of z~6-7 LAEs for the first time. The stacked LAEs have as blue UV continua as the HST/WFC3 z-dropout galaxies of similar Muv, with a spectral slope beta ~ -3, but at the same time they have red UV-to-optical colors with detection in the 3.6um band. Using SED fitting we find that the stacked LAEs have low stellar masses of ~(3-10)*10^7 Msun, very young ages of ~1-3 Myr, negligible dust extinction, and strong nebular emission from the ionized interstellar medium, although the z=6.6 object is fitted similarly well with high-mass models without nebular emission; inclusion of nebular emission reproduces the red UV-to-optical color while keeping the UV color sufficiently blue. We infer that typical LAEs at z~6-7 are building blocks of galaxies seen at lower redshifts. We find a tentative decrease in the Lyman alpha escape fraction from z=5.7 to 6.6, which may imply an increase in the intergalactic medium neutral fraction. From the minimum contribution of nebular emission required to fit the observed SEDs, we place an upper limit on the escape fraction of ionizing photons to be f_esc^ion~0.6 at z=5.7 and ~0.9 at z=6.6. We also compare the stellar populations of our LAEs with that of stacked HST/WFC3 z-dropout galaxies.

We present Lya luminosity function (LF), clustering measurements, and Lya line profiles based on the largest sample, to date, of 207 Lya emitters (LAEs) at z=6.6 on the 1-deg^2 sky of Subaru/XMM-Newton Deep Survey (SXDS) field. Our z=6.6 Lya LF including cosmic variance estimates yields the best-fit Schechter parameters of phi*=8.5 +3.0/-2.2 x10^(-4) Mpc^(-3) and L*(Lya)=4.4 +/-0.6 x10^42 erg s^(-1) with a fixed alpha=-1.5, and indicates a decrease from z=5.7 at the >~90% confidence level. However, this decrease is not large, only =~30% in Lya luminosity, which is too small to be identified in the previous studies. A clustering signal of z=6.6 LAEs is detected for the first time. We obtain the correlation length of r_0=2-5 h^(-1) Mpc and bias of b=3-6, and find no significant boost of clustering amplitude by reionization at z=6.6. The average hosting dark halo mass inferred from clustering is 10^10-10^11 Mo, and duty cycle of LAE population is roughly ~1% albeit with large uncertainties. The average of our high-quality Keck/DEIMOS spectra shows an FWHM velocity width of 251 +/-16 km s^(-1). We find no large evolution of Lya line profile from z=5.7 to 6.6, and no anti-correlation between Lya luminosity and line width at z=6.6. The combination of various reionization models and our observational results about the LF, clustering, and line profile indicates that there would exist a small decrease of IGM's Lya transmission owing to reionization, but that the hydrogen IGM is not highly neutral at z=6.6. Our neutral-hydrogen fraction constraint implies that the major reionization process took place at z>~7.

We conducted a deep narrowband NB973 (FWHM = 200 A centered at 9755 A) survey of z=7 Lyman alpha emitters (LAEs) in the Subaru/XMM-Newton Deep Survey Field, using the fully depleted CCDs newly installed on the Subaru Telescope Suprime-Cam, which is twice more sensitive to z=7 Lyman alpha at ~ 1 micron than the previous CCDs. Reaching the depth 0.5 magnitude deeper than our previous survey in the Subaru Deep Field that led to the discovery of a z=6.96 LAE, we detected three probable z=7 LAE candidates. Even if all the candidates are real, the Lyman alpha luminosity function (LF) at z=7 shows a significant deficit from the LF at z=5.7 determined by previous surveys. The LAE number and Lyman alpha luminosity densities at z=7 is ~ 7.7-54% and ~5.5-39% of those at z=5.7 to the Lyman alpha line luminosity limit of L(Ly-alpha) >~ 9.2 x 10^{42} erg s^{-1}. This could be due to evolution of the LAE population at these epochs as a recent galaxy evolution model predicts that the LAE modestly evolves from z=5.7 to 7. However, even after correcting for this effect of galaxy evolution on the decrease in LAE number density, the z=7 Lyman alpha LF still shows a deficit from z=5.7 LF. This might reflect the attenuation of Lyman alpha emission by neutral hydrogen remaining at the epoch of reionization and suggests that reionization of the universe might not be complete yet at z=7. If we attribute the density deficit to reionization, the intergalactic medium (IGM) transmission for Lyman alpha photons at z=7 would be 0.4 <= T_{Ly-alpha}^{IGM} <= 1, supporting the possible higher neutral fraction at the earlier epochs at z > 6 suggested by the previous surveys of z=5.7-7 LAEs, z ~ 6 quasars and z > 6 gamma-ray bursts.

We obtained mid-infrared 3.6 and 4.5 micron imaging of a z=6.96 Lyman alpha emitter (LAE) IOK-1 discovered in the Subaru Deep Field, using Spitzer Space Telescope Infrared Array Camera observations. After removal of a nearby bright source, we find that IOK-1 is not significantly detected in any of these infrared bands to m_3.6 ~ 24.00 and m_4.5 ~ 23.54 at 3 sigma. Fitting population synthesis models to the spectral energy distribution consisting of the upper limit fluxes of the optical to infrared non-detection images and fluxes in detection images, we constrain the stellar mass M* of IOK-1. This LAE could have either a mass as low as M* <~ 2-9 x 10^8 Msun for the young age (<~ 10 Myr) and the low dust reddening (A_V ~ 0) or a mass as large as M* <~ 1-4 x 10^{10} Msun for either the old age (> 100 Myr) or the high dust reddening (A_V ~ 1.5). This would be within the range of masses of z ~ 3-6.6 LAEs studied to date, ~ 10^6-10^{10} Msun. Hence, IOK-1 is not a particularly unique galaxy with extremely high mass or low mass but is similar to one of the LAEs seen at the later epochs.

We investigate the sky distribution of z ~ 6 Lyman break galaxies selected as i'-dropouts having i' - z' > 1.45 down to z' < 26.5 in the Subaru Deep Field (SDF). We discover 37 i'-dropouts clustered in a projected comoving 21.6 x 21.6 Mpc^2 region at z = 6, showing a local density excess. Carrying out follow-up spectroscopy, we identify four of them as Lyman-alpha emitters at z = 5.92, 6.01, 6.03 and 6.03 (spread over a distance of 46.6 Mpc). The number density of the cluster itself in SDF is ~ 2.2 x 10^{-7} Mpc^{-3}, smaller than those of protoclusters (i.e., forming galaxy clusters) at z ~ 2-5.7. Also, the structure shows ~4-21 times larger galaxy number density than those of z ~ 6 galaxies in a general field. It has a mass of M ~ 1.5^{+1.8}_{-0.5} x 10^{15}M_sun, comparable to those of z ~ 0-5 protoclusters. Since the contamination of our sample by interlopers is estimated to be quite low, 5.9%, most of the other unconfirmed i'-dropouts in the overdense region can be also z ~ 6 galaxies. Hence, it could be a candidate forming cluster at z ~ 6, representing a progenitor of galaxy clusters seen in the recent-day Universe.

(abridged) We present a study of stellar population of LAEs at z=4.86 in GOODS-N and its flanking field. With the publicly available IRAC data in GOODS-N and further IRAC observations in the flanking fields, we select five LAEs which are not contaminated by neighboring objects in IRAC images and construct their observed SEDs with I_c, z', IRAC 3.6micron, and 4.5micron band photometry. The SEDs cover the rest-frame UV to optical wavelengths. We derive stellar masses, ages, color excesses, and star formation rates of five LAEs using SED fitting method. Assuming the constant star formation history, we find that the stellar masses range from 10^8 to $10^{10} Msun with the median value of 2.5x10^9 Msun. The derived ages range from very young ages (7.4 Myr) to 437 Myr with a median age of 25 Myr. The color excess E(B-V) are between 0.1-0.4 mag. Star formation rates are 55-209 Msun/yr. A comparison of the stellar populations is made between three LAEs and 88 LBGs selected at the same redshift, in the same observed field, and down to the same limit of the rest-frame UV luminosity. These three LAEs are the brightest and reddest samples among the whole LAE samples at z=4.86. The LAEs distribute at the relatively faint part of UV-luminosity distribution of LBGs. Deriving the stellar properties of the LBGs by fitting their SEDs with the same model ensures that model difference does not affect the comparison. It is found that the stellar properties of the LAEs lie on distributions of those of LBGs. On average, the LAEs show less dust extinction, and lower star formation rates than LBGs, while the stellar mass of LAEs nearly lies in the middle part of the mass distribution of LBGs. However, the stellar properties of LAEs and LBGs are similar at the fixed UV or optical luminosity. We also examine the relations between the output properties from the SED fitting and the rest-frame Lya equivalent width.

We present the discovery of a large-scale structure of emission-line galaxies at redshift z=4.86 behind a massive cluster of galaxies, A1689. Previous spectroscopic observations of a galaxy, A1689-7.1 at z=4.87, near this structure, revealed a possible overdense region of inter-galactic medium (IGM) around the galaxy, which extends at least ~80 comoving Mpc along the line of sight. In order to investigate whether this z~5 IGM overdense region contains a galaxy overdensity, we undertook narrow- and broad-band imaging observations around A1689-7.1 with Subaru/Suprime-Cam. We detected 51 candidate Ly-alpha emitters at redshift z=4.86+-0.03 in the 32 x 24 arcmin^2 field of view. After correction for lensing by the foreground cluster, we found a large-scale (~20 x 60 comoving Mpc) overdense region of galaxies around A1689-7.1 in the source plane at z=4.86. The densest peak in this region has an overdensity of delta~4, suggesting that this structure is probably a good candidate for a protocluster. A1689-7.1 is located at the edge of this region, where the local galaxy density is ~1.6 times the mean density and is close to the density contrast in the IGM along the line of sight to A1689-7.1 estimated from the optical depth. The overdensities of galaxies we have found may suggest that at least some parts of the IGM overdense region have already started to form galaxies and moreover they relate to the formation of a protocluster. Although we lack information on the three dimensional distributions of both IGM and galaxy overdense regions, the similarity of the scales of both regions may suggest that the two are parts of a single large-scale structure, which would be an large edge-on sheet along the line of sight with a size of ~20 x 60 x 80 comoving Mpc.

We present the first wide-field H-alpha imaging survey around the distant cluster RXJ1716.4+6708 at z=0.81 with a narrow-band filter on MOIRCS/Subaru, down to SFR(Ha)~1Msun/yr. Combining with a wide-field MIR imaging survey with AKARI, we compare the unobscured and obscured star formation activities in the cluster. We find that both H-alpha emitters and MIR galaxies avoid the cluster centre and their spatial distribution is quite similar. Most of the H-alpha emitters show blue colours, but we find some H-alpha emitters on the red sequence. The MIR galaxies tend to be systematically redder than the H-alpha emitters probably due to heavy dust extinction. Interestingly, the red H-alpha emitters and the red MIR galaxies (i.e. dusty red galaxies) are most commonly seen in the medium-density environment such as cluster outskirts, groups and filaments, where optical colours of galaxies change. We also find that A(Ha) exceeds ~3 in extreme cases and that such very dusty galaxies are also located in the medium-density environment. These findings suggest that dusty star formation is triggered in the in-fall region of the cluster, implying a probable link between galaxy transition and dusty star formation. We finally calculate the cluster total SFR and find that the cluster total SFR based on H-alpha alone can be underestimated more than factor ~2 even after 1 mag extinction correction. We suggest that the mass-normalized cluster SFR rapidly declines since z~1 following ~(1+z)^6, although the uncertainty is still large.

We investigate the stellar populations of Lyman alpha emitters (LAEs) at z=3.1 and 3.7 in 0.65 deg^2 of the Subaru/XMM-Newton Deep Field, based on rest-frame UV-to-optical photometry obtained from the Subaru/XMM-Newton Deep Survey, the UKIDSS/Ultra Deep Survey, and the Spitzer legacy survey of the UKIDSS/UDS. Among a total of 302 LAEs (224 for z=3.1 and 78 for z=3.7), only 11 are detected in the K band, i.e., brighter than K(3sigma)=24.1 mag. Eight of the 11 K-detected LAEs are spectroscopically confirmed. We find that the K-undetected objects, which should closely represent the LAE population as a whole, have low stellar masses of ~ 10^8 - 10^8.5 Msun, modest SFRs of 1 - 100 Msun yr^-1, and modest dust extinction of E(B-V) < 0.2. The K-detected objects are massive, Mstar ~ 10^9 - 10^10.5 Msun, and have significant dust extinction with a median of E(B-V) ~= 0.3. Four K-detected objects with the reddest spectral energy distributions, two of which are spectroscopically confirmed, are heavily obscured with E(B-V) ~ 0.65, and their continua resemble those of some local ULIRGs. Interestingly, they have large Lyman alpha equivalent widths ~= 70 - 250 A. If these four are excluded, our sample has a weak anti-correlation between EW(Lya) and Mstar. We compare the stellar masses and the specific star formation rates (sSFR) of LAEs with those of Lyman-break galaxies, distant red galaxies, submillimetre galaxies, and I- or K-selected galaxies with z_phot ~ 3. We find that the LAE population is the least massive among all the galaxy populations in question, but with relatively high sSFRs, while NIR-detected LAEs have Mstar and sSFR similar to LBGs. Our reddest four LAEs have very high sSFRs in spite of large Mstar, thus occupying a unique region in the Mstar versus sSFR space. (abridged)

We present an unbiased deep [OII] emission survey of a cluster XMMXCS J2215.9-1738 at z=1.46, the most distant cluster to date with a detection of extended X-ray emission. With wide-field optical and near-infrared cameras (Suprime-Cam and MOIRCS, respectively) on Subaru telescope, we performed deep imaging with a narrow-band filter NB912 (lambda_c=9139A, Delta_lambda=134A) as well as broad-band filters (B, z', J and Ks). From the photometric catalogues, we have identified 44 [OII] emitters in the cluster central region of 6'x6' down to a dust-free star formation rate of 2.6 Msun/yr (3 sigma). Interestingly, it is found that there are many [OII] emitters even in the central high density region. In fact, the fraction of [OII] emitters to the cluster members as well as their star formation rates and equivalent widths stay almost constant with decreasing cluster-centric distance up to the cluster core. Unlike clusters at lower redshifts (z<1) where star formation activity is mostly quenched in their central regions, this higher redshift 2215 cluster shows its high star formation activity even at its centre, suggesting that we are beginning to enter the formation epoch of some galaxies in the cluster core eventually. Moreover, we find a deficit of galaxies on the red sequence at magnitudes fainter than ~M*+0.5 on the colour-magnitude diagram. This break magnitude is brighter than that of lower redshift clusters, and it is likely that we are seeing the formation phase of more massive red galaxies in the cluster core at z~1. These results may indicate inside-out and down-sizing propagation of star formation activity in the course of cluster evolution.

We present results of our large-area survey for z'-band dropout galaxies at z=7 in a 1568 arcmin^2 sky area covering the SDF and GOODS-N fields. Combining our ultra-deep Subaru/Suprime-Cam z'- and y-band (lambda_eff=1um) images with legacy data of Subaru and HST, we have identified 22 bright z-dropout galaxies down to y=26, one of which has a spectroscopic redshift of z=6.96 determined from Lya emission. The z=7 luminosity function (LF) yields the best-fit Schechter parameters of phi*=0.69 +2.62/-0.55 x10^(-3) Mpc^(-3), Muv*=-20.10 +/-0.76 mag, and alpha=-1.72 +/-0.65, and indicates a decrease from z=6 at a >95% confidence level. This decrease is beyond the cosmic variance in our two fields, which is estimated to be a factor of <~2. We have found that the cosmic star formation rate density drops from the peak at z=2-3 to z=7 roughly by a factor of ~10 but not larger than ~100. A comparison with the reionization models suggests either that the Universe could not be totally ionized by only galaxies at z=7, or more likely that properties of galaxies at z=7 are different from those at low redshifts having, e.g., a larger escape fraction (>~0.2), and/or a flatter IMF. Our SDF z-dropout galaxies appear to form 60-Mpc long filamentary structures, and the z=6.96 galaxy with Lya emission is located at the center of an overdense region consisting of four UV bright dropout candidates, which might suggest an existence of a well-developed ionized bubble at z=7.

Abridged: A photometric sample of ~7100 V<25.3 Lyman break galaxies (LBGs) has been selected by combining Subaru/Suprime-Cam BVRci'z' data with deep GALEX/NUV imaging of the Subaru Deep Field. Follow-up spectroscopy confirmed 24 LBGs at 1.5<z<2.7. Among the optical spectra, 12 have Ly-alpha emission with rest-frame equivalent widths of ~5-60AA. The success rate for identifying LBGs as NUV-dropouts at 1.5<z<2.7 is 86%. The rest-frame UV (1700AA) luminosity function (LF) is constructed from the photometric sample with corrections for stellar contamination and z<1.5 interlopers. The LF is 1.7+/-0.1 times higher than those of z~2 BXs and z~3 LBGs. Three explanations were considered, and it is argued that significantly underestimating low-z contamination or effective comoving volume is unlikely: the former would be inconsistent with the spectroscopic sample at 93% confidence, and the second explanation would not resolve the discrepancy. The third scenario is that different photometric selection of the samples yields non-identical galaxy populations, such that some BX galaxies are LBGs and vice versa. This argument is supported by a higher surface density of LBGs at all magnitudes while the redshift distribution of the two populations is nearly identical. This study, when combined with other star-formation rate (SFR) density UV measurements from LBG surveys, indicates that there is a rise in the SFR density: a factor of 3-6 (3-10) increase from z~5 (z~6) to z~2, followed by a decrease to z~0. This result, along with past sub-mm studies that find a peak at z~2 in their redshift distribution, suggest that z~2 is the epoch of peak star-formation. Additional spectroscopy is required to characterize the complete shape of the z~2 LBG UV LF via measurements of contamination and accurate distances.

We report the discovery of a giant Lya emitter (LAE) with a Spitzer/IRAC counterpart near the reionization epoch at z=6.595. The giant LAE is found from the extensive 1 deg^2 Subaru narrow-band survey for z=6.6 LAEs in the Subaru/XMM-Newton Deep Survey (SXDS) field, and subsequently identified by deep spectroscopy of Keck/DEIMOS and Magellan/IMACS. Among our 207 LAE candidates, this LAE is not only the brightest narrow-band object with L(Lya) = 3.9+/-0.2 x 10^43 erg/s in our survey volume of 10^6 Mpc^3, but also a spatially extended Lya nebula with the largest isophotal area whose major axis is at least ~3''. This object is more likely to be a large Lya nebula with a size of >~17 kpc than to be a strongly-lensed galaxy by a foreground object. Our Keck spectrum with medium-high spectral and spatial resolutions suggests that the velocity width is v(FWHM)=251+/-21 km/s, and that the line-center velocity changes by ~60 km/s in a 10-kpc range. The stellar mass and star-formation rate are estimated to be 0.9-5.0 x 10^10 Mo and >34 Mo/yr, respectively, from the combination of deep optical to infrared images of Subaru, UKIDSS-Ultra Deep Survey, and Spitzer/IRAC. Although the nature of this object is not yet clearly understood, this could be an important object for studying cooling clouds accreting onto a massive halo, or forming-massive galaxies with significant outflows contributing to cosmic reionization and metal enrichment of inter-galactic medium.

Efficient selection of emission line galaxies at z > 1 by photometric information in wide field surveys is one of the keys for future spectroscopic surveys to constrain dark energy using the baryon acoustic oscillation (BAO) signature. Here we estimate the H alpha and [O II] line luminosity functions of galaxies at z = 0.5-1.7 using a novel approach where multi-wavelength imaging data is used to jointly estimate both photometric redshifts and star-formation rates. These photometric estimates of line luminosities at high-redshift use the large data sets of the Subaru Deep Field and Subaru XMM-Newton Deep Field (covering \sim 1 deg^2) and are calibrated with the spectroscopic data of the local Sloan Digital Sky Survey galaxies. The derived luminosity functions (especially H alpha) are in reasonable agreement with the past estimates based on spectroscopic or narrow-band-filter surveys. This dataset is useful for examining the photometric selection of target galaxies for BAO surveys because of the large cosmological volume covered and the large number of galaxies with detailed photometric information. We use the sample to derive the photometric and physical properties of emission line galaxies to assist planning for future spectroscopic BAO surveys. We also show some examples of photometric selection procedures which can efficiently select these emission line galaxies.

We present a search for CO(3-2) emission in SDF-26821, a BzK-selected star-forming galaxy (sBzK) at z = 2.044, using the 45-m telescope of the Nobeyama Radio Observatory and the Nobeyama Millimeter Array. We do not detect significant emission and derive 2 \sigma limits: the CO luminosity of L'CO < 3.1 x 10^10 K km s^{-1} pc^{-2}, the ratio of far-infrared luminosity to CO luminosity of L_FIR/L'CO > 57 Lsun (K km s^{-1} pc^{-2})^{-1}, and the molecular gas mass of M_H2 < 2.5 x 10^10 Msun, assuming a velocity width of 200 km s^{-1} and a CO-to-H2 conversion factor of alpha_CO=0.8 Msun (K km s^{-1} pc^{-2})^{-1}. The ratio of L_FIR/L'CO, a measure of star formation efficiency (SFE), is comparable to or higher than the two z ~ 1.5 sBzKs detected in CO(2-1) previously, suggesting that sBzKs can have a wide range of SFEs. Comparisons of far-infrared luminosity, gas mass, and stellar mass among the sBzKs suggest that SDF-26821 is at an earlier stage of forming stars with a similar SFE and/or more efficiently forming stars than the two z ~ 1.5 sBzKs. The higher SFEs and specific star formation rates of the sBzKs compared to local spirals are indicative of the difference in star formation modes between these systems, suggesting that sBzKs are not just scaled-up versions of local spirals.

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

We present environmental dependence of dusty star forming activity in and around the cluster RXJ1716.4+6708 at z=0.81 based on wide-field and multi-wavelength observations with Suprime-Cam on the Subaru telescope and IRC onboard the AKARI satellite. Our optical data shows that the optical colour distribution of galaxies starts to dramatically change from blue to red at the medium-density environment such as cluster outskirts, groups and filaments. By combining with infrared data, we find that 15 micron galaxies tend to have optical colours between the red sequence and the blue cloud with a tail into the red sequence. The spatial distribution of the 15 micron galaxies over ~200 arcmin^2 around the cluster reveals that few 15 micron galaxies are detected in the cluster central region. This is probably due to the low star forming activity in the cluster core. However, interestingly, the fraction of 15 micron galaxies in the medium-density environments is as high as in the low-density field, despite the fact that the optical colours start to change in the medium-density environments. Furthermore, we find that 15 micron galaxies which have optically red colours (candidates for dusty red galaxies) and galaxies with high specific star formation rates are also concentrated in the medium-density environment. These results imply that the star forming activity in galaxies in groups and filaments is enhanced due to some environmental effects specific to the medium-density environment, and such a phenomenon is probably directly connected to the truncation of star forming activity in galaxies seen as the dramatic change in optical colours in such environments.

We investigate the relationship between two massive star-forming galaxy populations at redshift z~ 2; i.e. submillimetre galaxies (SMGs) and BzK-selected galaxies (BzKs). Out of 60 SMGs found in the Subaru/XMM-Newton deep field, we collect optical--NIR photometry of 28 radio counterparts for 24 SMGs, based on refined sky positions with a radio map for 35 SMGs (Ivison et al. 2007). We find a correlation between their K-band magnitudes and BzK [= (z-K)-(B-z)] colours: almost all of the K-faint (K_AB > 21.3) radio-detected SMGs have BzK>-0.2, and therefore BzKs. This result gives strong support to perform direct optical identification of SMGs by searching for BzKs around SMGs. We calculate the formal significance (P' value) for each of the BzK associations around radio-undetected SMGs, and find 6 new robust identifications, including one double identification. From this analysis, we obtain the current best estimate on the surface density of BzK-selected SMGs, which indicates that only ~1 per cent of BzKs are SMGs. If BzKs are normal disk-like galaxies at z~ 2 as indicated by the correlation between their star formation rate (SFR) and stellar mass and also by dynamical properties, SMGs are likely to be merging BzKs. In this case, a typical enhancement of SFR due to merging is only a factor of ~ 3, which is an order of magnitude lower than that of local ULIRGs. This may indicate that most of the merging BzKs could be observed as SMGs. Considering a possible high fraction of mergers at z~ 2 (at least it would be higher than the fraction at z> 1 of ~ 10 per cent), it is rather puzzling to find such a low fraction of SMGs in the progenitor population, i.e. BzKs.

We present luminosity functions (LFs) and various properties of Lya emitters (LAEs) at z=3.1, 3.7, and 5.7, in a 1 deg^2 sky of the Subaru/XMM-Newton Deep Survey (SXDS) Field. We obtain a photometric sample of 858 LAE candidates based on deep Subaru/Suprime-Cam imaging data, and a spectroscopic sample of 84 confirmed LAEs from Subaru/FOCAS and VLT/VIMOS spectroscopy in a survey volume of ~10^6 Mpc^3 with a limiting Lya luminosity of ~3x10^42 erg/s. We derive the LFs of Lya and UV-continuum (~1500 \AA) for each redshift, taking into account the statistical error and the field-to-field variation. We find that the apparent Lya LF shows no significant evolution between z=3.1 and 5.7 within factors of 1.8 and 2.7 in L* and phi*, respectively. On the other hand, the UV LF of LAEs increases from z=3.1 to 5.7, indicating that galaxies with Lya emission are more common at earlier epochs. We identify six LAEs with AGN activities from our spectra combined with VLA, Spitzer, and XMM-Newton data. Among the photometrically selected LAEs at z=3.1 and 3.7, only ~1 % show AGN activities, while the brightest LAEs with logL(Lya) >~ 43.4-43.6 erg/s appear to always host AGNs. Our LAEs are bluer in UV-continuum color than dropout galaxies, suggesting lower extinction and/or younger stellar populations. Our stacking analyses provide upper limits to the radio luminosity and the f(HeII)/f(Lya) line fraction, and constrain the hidden star formation (+low-luminosity AGN) and the primordial population in LAEs.

We made a narrowband NB973 (bandwidth of 200A at 9755A) imaging of the Subaru Deep Field (SDF) and found two z=7 Lyman alpha emitter (LAE) candidates down to NB973=24.9. Carrying out deep follow-up spectroscopy, we identified one of them as a real z=6.96 LAE. This has shown that galaxy formation was in progress just 750 Myr after the Big Bang. Meanwhile, the Lyman alpha line luminosity function of LAE is known to decline from z=5.7 to 6.6 in the SDF. L* at z=6.6 is 40-60% of that at z=5.7. We also confirm that the number density of z=7 LAE is only 17% of the density at z=6.6 comparing the latest SDF LAE samples. This series of significant decreases in LAE density with increasing redshift can be the result of galaxy evolution during these epochs. However, using the UV continuum luminosity functions of LAEs, those of Lyman break galaxies, and a LAE evolution model based on the hierarchical clustering, we find that galaxy evolution alone cannot explain all the decrease in density. This extra density deficit can be interpreted as the attenuation of the Lyman alpha photons from LAEs due to a rapid evolution of neutral hydrogen fraction during the ongoing cosmic reionization at z~6.6-7.

We investigate clustering properties of Lyman-break galaxies (LBGs) at z~3 based on deep multi-waveband imaging data from optical to near-infrared wavelengths in the Subaru/XMM-Newton Deep Field. The LBGs are selected by U-V and V-z' colors in one contiguous area of 561 arcmin^2 down to z'=25.5. We study the dependence of the clustering strength on rest-frame UV and optical magnitudes, which can be indicators of star formation rate and stellar mass, respectively. The correlation length is found to be a strong function of both UV and optical magnitudes with brighter galaxies being more clustered than faint ones in both cases. Furthermore, the correlation length is dependent on a combination of UV and optical magnitudes in the sense that galaxies bright in optical magnitude have large correlation lengths irrespective of UV magnitude, while galaxies faint in optical magnitude have correlation lengths decreasing with decreasing UV brightness. These results suggest that galaxies with large stellar masses always belong to massive halos in which they can have various star formation rates, while galaxies with small stellar masses reside in less massive halos only if they have low star formation rates. There appears to be an upper limit to the stellar mass and the star formation rate which is determined by the mass of hosting dark halos.

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 present environmental dependence of the build-up of the colour-magnitude relation (CMR) at z ~ 0.8. It is well established that massive early-type galaxies exhibit a tight CMR in clusters up to at least z ~ 1. The faint end of the relation, however, has been much less explored especially at high redshifts primarily due to limited depths of the data. Some recent papers have reported a deficit of the faint red galaxies on the CMR at 0.8 < z < 1, but this has not been well confirmed yet and is still controversial. Using a deep, multi-colour, panoramic imaging data set of the distant cluster RXJ1716.4+6708 at z=0.81, newly taken with the Prime Focus Camera (Suprime-Cam) on the Subaru Telescope, we carry out an analysis of faint red galaxies with a care for incompleteness. We find that there is a sharp decline in the number of red galaxies toward the faint end of the CMR below M*+2. We compare our result with those for other clusters at z ~ 0.8 taken from the literature, which show or do not show the deficit. We suggest that the "deficit" of faint red galaxies is dependent on the richness or mass of the clusters, in the sense that poorer systems show stronger deficits. This indicates that the evolutionary stage of less massive galaxies depends critically on environment.

We have studied the submillimetre (submm) properties of the following classes of near-infrared (NIR)-selected massive galaxies at high redshifts: BzK-selected star-forming galaxies (BzKs); distant red galaxies (DRGs); and extremely red objects (EROs). We used the SCUBA HAlf Degree Extragalactic Survey (SHADES), the largest uniform submm survey to date. Partial overlap of SIRIUS/NIR images and SHADES in SXDF has allowed us to identify 4 submm-bright NIR-selected galaxies, which are detected in the mid-infrared, 24 micron, and the radio, 1.4 GHz. We find that all of our submm-bright NIR-selected galaxies satisfy the BzK selection criteria, except for one galaxy whose B-z and z-K colours are however close to the BzK colour boundary. Two of the submm-bright NIR-selected galaxies satisfy all of the selection criteria we considered, i.e. they belong to the BzK-DRG-ERO overlapping population, or `extremely red' BzKs. Although these extremely red BzKs are rare (0.25 arcmin^{-2}), up to 20 % of this population could be submm galaxies. This fraction is significantly higher than that found for other galaxy populations studied here. Via a stacking analysis, we have detected the 850 micron flux of submm-faint BzKs and EROs in our SCUBA maps. While the contribution of z~2 BzKs to the submm background is about 10--15 % and similar to that from EROs typically at z~1, BzKs have a higher fraction (~30 %) of submm flux in resolved sources compared with EROs and submm sources as a whole. From the SED fitting analysis for both submm-bright and submm-faint BzKs, we found no clear signature that submm-bright BzKs are experiencing a specifically luminous evolutionary phase, compared with submm-faint BzKs. An alternative explanation might be that submm-bright BzKs are more massive than submm-faint ones.

We carried out a wide-field V, I imaging survey of the Local Group dwarf spheroidal galaxy Leo II using the Subaru Prime Focus Camera on the 8.2-m Subaru Telescope. The survey covered an area of 26.67 x 26.67 arcmin^2, far beyond the tidal radius of Leo II (8.63 arcmin), down to the limiting magnitude of V ~26, which is roughly 1 mag deeper than the turn-off point of the main sequence stars of Leo II. Radial number density profiles of bright and faint red giant branch (RGB) stars were found to change their slopes at around the tidal radius, and extend beyond the tidal radius with shallower slopes. A smoothed surface brightness map of Leo II suggests the existence of a small substructure of globular cluster luminosity beyond the tidal radius. We investigated the properties of the stellar population by means of the color-magnitude diagram. The horizontal branch (HB) morphology index shows a radial gradient in which red HB stars are more concentrated than blue HB stars, which is common to many Local Group dwarf spheroidal galaxies. The color distribution of RGB stars around the mean RGB sequence shows a larger dispersion at the center than in the outskirts, indicating a mixture of stellar populations at the center and a more homogeneous population in the outskirts. Based on the age estimation using subgiant branch (SGB) stars, we found that although the major star formation took place ~8 Gyr ago, a considerable stellar population younger than 8 Gyr is found at the center; such a younger population is insignificant in the outskirts.

Spatially extended Ly-alpha sources that are faint and/or compact in coninuum are candidates for extremely young (~< 10^7 yrs) galaxies at high redshifts. We present medium-resolution (R~2000) spectroscopy of such extended Ly-alpha sources found in our previous study at z~3-5, using VLT/VIMOS. The deep spectroscopy showed that all 18 objects we observed have large equivalent widths (EWs) exceeding 100 A. For about 30% of our sample (five objects), we identified conspicuous asymmetry on the profiles of the Ly-alpha line. They show broad wing emission components on the red side, and sharp cut-off on the blue side of the Ly-alpha line. Such asymmetry is often seen in superwind galaxies known to date, and also consistent with a theoretical prediction of superwind activity. There are eight objects (8/18 ~ 40%) that have large EWs exceeding 200 A, and no clear signature of superwind activities. Such large EWs cannot be explained in terms of photo-ionization by a moderately old (>10^7 yrs) stellar population, even with a top-heavy IMF or an extremely low metallicity. These eight objects clearly show a positive correlation between the Ly-alpha luminosity and the velocity width. This suggests that these eight objects are good candidates for forming-galaxies in a gas-cooling phase.

We use the BzK color selection proposed by Daddi et al. (2004) to obtain a sample of 1092 faint star-forming galaxies (hereafter sBzKs) from 180 arcmin^2 in the Subaru Deep Field. This sample represents star-forming galaxies at 1.4 < z < 2.5 down to K(AB)=23.2, which roughly corresponds to a stellar-mass limit of ~ 1 x 10^{10} Msun. We measure the angular correlation function (ACF) of these sBzKs to be w(theta) = (0.58 +- 0.13) x theta["]^{-0.8} and translate the amplitude into the correlation length assuming a reasonable redshift distribution. The resulting value, r0 = 3.2^{+0.6}_{-0.7} h^{-1} Mpc, suggests that our sBzKs reside in haloes with a typical mass of 2.8 x 10^{11} Msun. Combining this halo mass estimate with those for brighter samples of Kong et al. (2006), we find that the mass of dark haloes largely increases with K brightness, a measure of the stellar mass. Comparison with other galaxy populations suggests that faint sBzKs (K(AB)<23.2) and Lyman Break Galaxies at z ~ 2 are similar populations hosted by relatively low-mass haloes, while bright sBzKs (K(AB)<21) reside in haloes comparable to or more massive than those of Distant Red Galaxies and Extremely Red Objects. Using the extended Press-Schechter formalism, we predict that present-day descendants of haloes hosting sBzKs span a wide mass range depending on K brightness, from lower than that of the Milky Way up to those of richest clusters.

We study the galaxy population in the central region and a region about 0.6 Mpc away from the center of the Hydra I cluster in B- and Rc-bands down to M~-10 using the Subaru Suprime-Cam photometry. We find that the luminosity function of the entire population has a slightly steeper slope (alpha~-1.6) in the range of -20<M<-10 than those reported for other clusters in slightly brighter ranges. The slope appears to be steeper in poorer clusters. The number of faint galaxies (M>-14) increases in the cluster center as well as bright galaxies. The Hydra I cluster is dominated by red dwarfs and the luminosity function shows a significant upturn at M~-16 as is seen in several other nearby clusters, but not in the field. This upturn and the variation in the faint-end slope of the luminosity function may be caused by the cluster environment or the evolution history of individual clusters.

Abridged: Subaru Deep Field line-emitting galaxies in four narrow-band filters at low and intermediate redshifts are presented. Broad-band colors, follow-up optical spectroscopy, and multiple NB filters are used to distinguish Ha, [O II], and [O III] emitters between redshifts of 0.07 and 1.47 to construct their luminosity functions (LFs). These LFs are derived down to faint magnitudes, which allows for a more accurate determination of the faint end slope. With a large (N~200-900) sample for each redshift interval, a Schechter profile is fitted to each LF. Prior to dust extinction corrections, the [O III] and [O II] LFs reported in this paper agree reasonably well with those of Hippelein et al. The z=0.08 Ha LF, which reaches two orders of magnitude fainter than Gallego et al., is steeper by 25%. This indicates that there are more low luminosity star-forming galaxies for z<0.1. The faint end slope \alpha and \phi_{\star} show a strong evolution with redshift while L_{\star} show little evolution. The evolution in \alpha indicates that low-luminosity galaxies have a stronger evolution compared to brighter ones. Integrated star formation rate densities are derived via Ha, [O III], and [O II] for 0.07<z<1.47. A steep increase in the star-formation rate density, as a function of redshift, is seen for 0.4<z<0.9. For z>1, the star-formation rate densities are more or less constant. The latter is consistent with previous UV and [O II] measurements. Below z<0.4, the SFR densities are consistent with several Ha, [O II], and UV measurements, but other measurements are a factor of two higher. For example, the z=0.066-0.092 Ha LF agrees with those of Jones & Bland-Hawthorn, but at z=0.24 and 0.40, their number density is higher by a factor of two. This discrepancy can be explained by cosmic variance.

When galaxy formation started in the history of the Universe remains unclear. Studies of the cosmic microwave background indicate that the Universe, after initial cooling (following the Big Bang), was reheated and reionized by hot stars in newborn galaxies at a redshift in the range 6 < z < 14 (ref. 1). Though several candidate galaxies at redshift z > 7 have been identified photometrically (refs 2,3), galaxies with spectroscopically confirmed redshifts have been confined to z < 6.6 (refs. 4-8). Here we report a spectroscopic redshift of z = 6.96 (corresponding to just 750 Myr after the Big Bang) for a galaxy whose spectrum clearly shows Lyman-alpha emission at 9,682 A, indicating active star formation at a rate of about 10 M_sun/yr, where M_sun us the mass of the Sun. This demonstrates that galaxy formation was under way when the Universe was only about 6 per cent of its present age. The number density of galaxies at z = 7 seems to be only 18-36 per cent of the density at z = 6.6.

We investigate the luminosity functions of Lyman-break galaxies (LBG) at z~4 and 5 based on the optical imaging data obtained in the Subaru Deep Field Project. Three samples of LBGs in a contiguous 875 arcmin^2 area are constructed. One consists of 3,808 LBGs at z~4 down to i'=26.85 selected with the B-R vs R-i' diagram. The other two consist of 539 and 240 LBGs at z~5 down to z'=26.05 selected with two kinds of two-color diagrams: V-i' vs i'-z' and R-i' vs i'-z'. The adopted selection criteria are proved to be fairly reliable by spectroscopic observation. We derive the luminosity functions of the LBGs at rest-frame ultraviolet wavelengths down to M_{UV}=-19.2 at z~4 and M_{UV}=-20.3 at z~5. We find clear evolution of the luminosity function over the redshift range of 0<z<6, which is accounted for by a sole change in the characteristic magnitude, M^*. The cosmic star formation rate (SFR) density at z~4 and z~5 is measured from the luminosity functions. We examine the evolution of the cosmic SFR density and its luminosity dependence over 0<z<6. The SFR density contributed from brighter galaxies is found to change more drastically with cosmic time. The contribution from brighter galaxies has a sharp peak around z=3-4, while that from fainter galaxies evolves relatively mildly with a broad peak at earlier epoch. Combining the observed SFR density with the standard Cold Dark Matter model, we compute the cosmic SFR per unit baryon mass in dark haloes, i.e., the specific SFR. The specific SFR is found to scale with redshift as (1+z)^3 up to z~4, implying that the efficiency of star formation is on average higher at higher redshift in proportion to the cooling rate within dark haloes, while this is not simply the case at z>4.

We report an extensive search for Lyman-alpha emitters (LAEs) at z=6.5 in the Subaru Deep Field. Subsequent spectroscopy with Subaru and Keck identified eight more LAEs, giving a total of 17 spectroscopically confirmed LAEs at z=6.5. Based on this spectroscopic sample of 17, complemented by a photometric sample of 58 LAEs, we have derived a more accurate Lyman-alpha luminosity function of LAEs at z=6.5, which reveals an apparent deficit at the bright end of ~0.75 mag fainter L*, compared with that observed at z=5.7. The difference in the LAE luminosity functions between z=5.7 and 6.5 is significant at the 3-sigma level, which is reduced to 2-sigma when cosmic variance is taken into account. This result may imply that the reionization of the universe has not been completed at z=6.5. We found that the spatial distribution of LAEs at z=6.5 was homogeneous over the field. We discuss the implications of these results for the reionization of the universe.

Spatially extended Ly-alpha sources which are faint and compact in coninuum are candidates for extremely young galaxies (age of ~<10e+7 yrs) at high redshifts. We carried out a systematic survey for extended Ly-alpha sources, using deep intermediate-band imaging data taken with the Suprime-Cam on the Subaru Telescope. Our survey covers a field of view of 33' * 25' and a redshift range of 3.24<z<4.95, down to a Ly-alpha flux of ~1e-17 ergs/cm^2/s. We identified 41 extended Ly-alpha sources. The redshift distribution of these sources shows that this kind of objects are common in the early universe over the surveyed redshift range. The identified objects have typical sizes of \~10-15 kpc and luminosities of ~10^42 ergs/s. Follow-up spectroscopy made for 7 of the 41 objects showed that our sample suffers from little contamination. All 7 objects have large equivalent widths of Ly-alpha emission line, all but one exceeding 240A in the rest frame. The large equivalent widths suggest that their extended Ly-alpha emissions are unlikely to be due to normal starbursts, but possibly originated from very young galaxies. All 41 objects in our sample have much smaller Ly-alpha luminosities than the two Ly-alpha Blobs (LABs) found by Steidel et al. (2000) in spite of our much larger survey volume. This suggests that large and luminous extended Ly-alpha objects like the two LABs are very rare and are clustered in overdense regions.

We present the properties of Ly Alpha emitters (LAEs) at z = 5.7 in the Subaru Deep Field. A photometric sample of 89 LAE candidates is constructed from narrow-band (NB816) data down to NB816 = 26.0 (AB) in a continuous 725 arcmin^2 area. Spectra of 39 objects satisfying the photometric selection criteria for LAEs were obtained with Subaru and Keck II Telescopes, among which 28 were confirmed LAEs, one was a nearby galaxy, and eight were unclassified. We also obtained spectra of another 24 NB816-excess objects in the field, identifying six additional LAEs. We find that the Ly Alpha luminosity function derived from the photometric sample is reproduced well by a Schechter function with L* = (7.9+3.0-2.2) x 10^42 erg/s and phi* = (6.3+3.0-2.0) x 10^-4 Mpc^-3 for alpha = -1.5 (fixed) over the whole luminosity range of L ~= 3x10^42 - 3x10^43 erg/s. We then measure rest-frame Ly Alpha equivalent widths (EWs) for the confirmed LAEs, to find that the median among the 28 objects satisfying the photometric selection criteria is W_0^i = 233 A. We infer that 30% - 40% of LAEs at z=5.7 exceed W_0^i = 240 A. These large-EW objects probably cannot be accounted for by ordinary star-forming populations with a Salpeter IMF. We also find that LAEs with fainter far-UV luminosities have larger EWs. Finally, we derive the far-UV luminosity function of LAEs down to M_UV ~= -19.6 using the photometric sample, and compare it with that of Lyman-break galaxies (LBGs). We find that as high as about 80% of LBGs at z ~ 6 have W_0^i >= 100 A, in sharp contrast to lower-z counterparts.

We explored the clustering properties of Lyman Break Galaxies (LBGs) at z=4 and 5 with an angular two-point correlation function on the basis of the very deep and wide Subaru Deep Field data. We found an apparent dependence of the correlation function slope on UV luminosity for LBGs at both z=4 and 5. More luminous LBGs have a steeper correlation function. To compare these observational results, we constructed numerical mock LBG catalogs based on a semianalytic model of hierarchical clustering combined with high-resolution N-body simulation, carefully mimicking the observational selection effects. The luminosity functions for LBGs predicted by this mock catalog were found to be almost consistent with the observation. Moreover, the overall correlation functions of LBGs were reproduced reasonably well. The observed dependence of the clustering on UV luminosity was not reproduced by the model, unless subsamples of distinct halo mass were considered. That is, LBGs belonging to more massive dark haloes had steeper and larger-amplitude correlation functions. With this model, we found that LBG multiplicity in massive dark halos amplifies the clustering strength at small scales, which steepens the slope of the correlation function. The hierarchical clustering model could therefore be reconciled with the observed luminosity-dependence of the angular correlation function, if there is a tight correlation between UV luminosity and halo mass. Our finding that the slope of the correlation function depends on luminosity could be an indication that massive dark halos hosted multiple bright LBGs (abridged).

We report angular correlation function (ACF) of Lyman Break Galaxies (LBGs) with unprecedented statistical quality on the basis of 16,920 LBGs at z=4 detected in the 1 deg^2 sky of the Subaru/XMM-Newton Deep Field. The ACF significantly departs from a power law, and shows an excess on small scale. Particularly, the ACF of LBGs with i'<27.5 have a clear break between the small and large-scale regimes at the angular separation of ~7'' whose projected length corresponds to the virial radius of dark halos with a mass of 10^11-12 Mo, indicating multiple LBGs residing in a single dark halo. Both on small (2''<theta<3'') and large (40''<theta<400'') scales, clustering amplitudes monotonically increase with luminosity for the magnitude range of i'=24.5-27.5, and the small-scale clustering shows a stronger luminosity dependence than the large-scale clustering. The small-scale bias reaches b~10-50, and the outskirts of small-scale excess extend to a larger angular separation for brighter LBGs. The ACF and number density of LBGs can be explained by the cold dark matter model.

We report the discovery of a large-scale coherent filamentary structure of Lyman alpha emitters in a redshift space at z=3.1. We carried out spectroscopic observations to map the three dimensional structure of the belt-like feature of the Lyman alpha emitters discovered by our previous narrow-band imaging observations centered on the protocluster at z=3.1. The feature was found to consist of at least three physical filaments connecting with each other. The result is in qualitative agreement with the prediction of the 'biased' galaxy-formation theories that galaxies preferentially formed in large-scale filamentary or sheet-like mass overdensities in the early Universe. We also found that the two known giant Lyman alpha emission-line nebulae showing high star-formation activities are located near the intersection of these filaments, which presumably evolves into a massive cluster of galaxies in the local Universe. This may suggest that massive galaxy formation occurs at the characteristic place in the surrounding large-scale structure at high redshift.

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.

We report new spectroscopy of two star-forming galaxies with strong Ly_alpha emission at z=6.03 and z=6.04 in the Subaru Deep Field. These two objects are originally selected as i'-dropouts (i'-z' > 1.5) showing an interesting photometric property, the ``NB921 depression''. The NB921-band (centered at 9196A) magnitude is significantly depressed with respect to the z'-band magnitude. The optical spectra of these two objects exhibit asymmetric emission-lines at lambda_obs ~ 8540A and ~ 8560A, suggesting that these objects are Ly_alpha emitters at z~6. The rest-frame equivalent widths of the Ly_alpha emission of the two objects are 94A and 236A; the latter one is the Ly_alpha emitter with the largest Ly_alpha equivalent width at z > 6 ever spectroscopically confirmed. The spectroscopically measured Ly_alpha fluxes of these two objects are consistent with the interpretation that the NB921 depression is caused by the contribution of the strong Ly_alpha emission to the z'-band flux. Most of the NB921-depressed i'-dropout objects are thought to be strong Ly_alpha emitters at 6.0 < z < 6.5; Galactic L and T dwarfs and NB921-dropout galaxies at z > 6.6 do not dominate the NB921-depressed i'-dropout sample. Thus the NB921-depression method is very useful for finding high-z Ly_alpha emitters with a large Ly_alpha equivalent width over a large redshift range, 6.0 < z < 6.5. Although the broadband-selected sample at z ~ 3 contains only a small fraction of objects with a Ly_alpha equivalent width larger than 100A, the i'-dropout sample of the Subaru Deep Field contains a much larger fraction of such strong Ly_alpha emitters. This may imply a strong evolution of the Ly_alpha equivalent width from z > 6 to z ~ 3.

We discuss the environmental dependence of galaxy evolution based on deep panoramic imaging of two distant clusters taken with Suprime-Cam as part of the PISCES project. By combining with the SDSS data as a local counterpart for comparison, we construct a large sample of galaxies that spans wide ranges in environment, time, and stellar mass (or luminosity). We find that colours of galaxies, especially those of faint galaxies ($M_V>M_V^*+1$), change from blue to red at a break density as we go to denser regions. Based on local and global densities of galaxies, we classify three environments: field, groups, and clusters. We show that the cluster colour-magnitude relation is already built at $z=0.83$. In contrast to this, the bright-end of the field colour-magnitude relation has been vigorously built all the way down to the present-day and the build-up at the faint-end has not started yet. A possible interpretation of these results is that galaxies evolve in the 'down-sizing' fashion. That is, massive galaxies complete their star formation first and the truncation of star formation is propagated to smaller objects as time progresses. This trend is likely to depend on environment since the build-up of the colour-magnitude relation is delayed in lower-density environments. Therefore, we may suggest that the evolution of galaxies took place earliest in massive galaxies and in high density regions, and it is delayed in less massive galaxies and in lower density regions.

We report on the bright Lyman-break galaxies (LBGs) selected in a 767 arcmin^2 area of the Subaru Deep Field. The selection is made in the i-zR vs zB-zR plane, where zB and zR are new bandpasses with a central wavelength of 8842A and 9841A, respectively. This set of bandpasses enables us to separate well z~6 LBGs from foreground galaxies and Galactic cool stars. We detect 12 LBG candidates down to zR=25.4, and calculate the normalization of the rest-frame far-ultraviolet (FUV: 1400A) luminosity function at MFUV = -21.6 to be \phi(-21.6) = (2.6+/-0.7) x 10^{-5} mag^{-1} Mpc^{-3}. This must be the most reliable measurement ever obtained of the number density of bright z~6 LBGs, because it is more robust against both contamination and cosmic variance than previous values. The FUV luminosity density contributed from LBGs brighter than MFUV = -21.3 is (2.8+/-0.8) x 10^{24} ergs/s/Hz/Mpc^3, which is equivalent to a star formation rate density of (3.5+/-1.0) x 10^{-4} Msun/yr/Mpc^3. Combining our measurement with those at z<6 in the literature, we find that the FUV luminosity density of bright galaxies increases by an order of magnitude from z~6 to z~3 and then drops by 10^3 from z~3 to the present epoch, while the evolution of the total luminosity density is much milder. The evolutionary behavior of bright LBGs resembles that of luminous dusty star-forming galaxies and bright QSOs. The redshift of z~3 appears to be a remarkable era in the cosmic history when massive galaxies were being intensively formed.

We present the large-scale correlation function measured from a spectroscopic sample of 46,748 luminous red galaxies from the Sloan Digital Sky Survey. The survey region covers 0.72 h^{-3} Gpc^3 over 3816 square degrees and 0.16<z<0.47, making it the best sample yet for the study of large-scale structure. We find a well-detected peak in the correlation function at 100h^{-1} Mpc separation that is an excellent match to the predicted shape and location of the imprint of the recombination-epoch acoustic oscillations on the low-redshift clustering of matter. This detection demonstrates the linear growth of structure by gravitational instability between z=1000 and the present and confirms a firm prediction of the standard cosmological theory. The acoustic peak provides a standard ruler by which we can measure the ratio of the distances to z=0.35 and z=1089 to 4% fractional accuracy and the absolute distance to z=0.35 to 5% accuracy. From the overall shape of the correlation function, we measure the matter density Omega_mh^2 to 8% and find agreement with the value from cosmic microwave background (CMB) anisotropies. Independent of the constraints provided by the CMB acoustic scale, we find Omega_m = 0.273 +- 0.025 + 0.123 (1+w_0) + 0.137 Omega_K. Including the CMB acoustic scale, we find that the spatial curvature is Omega_K=-0.010+-0.009 if the dark energy is a cosmological constant. More generally, our results provide a measurement of cosmological distance, and hence an argument for dark energy, based on a geometric method with the same simple physics as the microwave background anisotropies. The standard cosmological model convincingly passes these new and robust tests of its fundamental properties.

We report the discovery of primeval large-scale structures (LSSs) including two proto-clusters in a forming phase at z=5.7. We carried out extensive deep narrow-band imaging in the 1 deg^2 sky of the Subaru/XMM-Newton Deep Field, and obtained a cosmic map of 515 Ly alpha emitters (LAEs) in a volume with a transverse dimension of 180 Mpc x 180 Mpc and a depth of ~40 Mpc in comoving units. This cosmic map shows filamentary LSSs, including clusters and surrounding 10-40 Mpc scale voids, similar to the present-day LSSs. Our spectroscopic follow-up observations identify overdense regions in which two dense clumps of LAEs with a sphere of 1-Mpc diameter in physical units are included. These clumps show about 130 times higher star-formation rate density, mainly due to a large overdensity, ~80, of LAEs. These clumps would be clusters in a formation phase involving a burst of galaxy formation.

We report the discovery of a probable large scale structure composed of many galaxy clumps around the known twin clusters at z=1.26 and z=1.27 in the Lynx region. Our analysis is based on deep, panoramic, and multi-colour imaging (26.4'x24.1') in VRi'z'-bands with the Suprime-Cam on the 8.2m Subaru telescope. This unique, deep and wide-field imaging data-set allows us for the first time to map out the galaxy distribution in the highest redshift supercluster known. We apply a photometric redshift technique to extract plausible cluster members at z~1.27 down to i'=26.15 (5sigma) corresponding to \~M*+2.5 at this redshift. From the 2-D distribution of these photometrically selected galaxies, we newly identify seven candidates of galaxy groups or clusters where the surface density of red galaxies is significantly high (>5sigma), in addition to the two known clusters. These candidates show clear red colour-magnitude sequences consistent with a passive evolution model, which suggests the existence of additional high density regions around the Lynx superclusters.

We investigate the environmental dependence of star formation and the morphology of galaxies in the local universe based on a volume-limited sample (0.030<z<0.065 and M_r<M_r^*+2) constructed from the data of the Sloan Digital Sky Survey. First, we focus on how galaxy properties change with local density. The star formation and morphology show a 'break' at the critical local density of log Sigma_crit~0.4 galaxies h_75^2 Mpc^-2, which is in agreement with previous studies. However, the break can be seen only for faint galaxies (M^*_r+1 < M_r < M^*_r+2), and bright galaxies (M_r < M^*_r+1) show no break. Thus, galaxies of different luminosities are found to show different environmental dependencies. Next, we examine dependencies on richness of galaxy systems. Star formation activities of galaxies are different from those of field galaxies even in systems as poor as sigma~100 km/s. This result suggests that environmental mechanisms that are effective only in rich systems, such as ram-pressure stripping of cold gas and harassment, have not played a major role in transforming galaxies into red early-type galaxies. Strangulation and interactions between galaxies, however, remain candidates of the driver of the environmental dependence. In the dense environments in the local universe, the slow transformation of faint galaxies occurs to some extent, but the transformation of bright galaxies is not clearly visible. We suggest that the evolution of bright galaxies is not strongly related to galaxy system, such as groups and clusters, while the evolution of faint galaxies is likely to be closely connected to galaxy system.

The Subaru Deep Field (SDF) project is a program of Subaru Observatory to carry out a deep galaxy survey over a blank field as large as 34'x27'. The program consists of very deep multi-band optical imaging, near infrared imaging for smaller portions of the field and follow-up optical spectroscopy. Major scientific goals of the project are to construct large samples of Lyman-break galaxies at z~4-5 and Lyman alpha emitters at z~5.7 and 6.6, and to make detailed studies these very high-redshift galaxy populations. In this paper, we describe the optical imaging observations and data reduction, presenting mosaicked images and object catalogs in seven bandpasses.The optical imaging was made through five broad-band filters, B, V, R, i', z', and two narrow-band filters, NB816 (lambda_c=8150A) and NB921 (lambda_c=9196A) with almost 10 hours long integrations for each band. The limiting magnitudes measured at 3-sigma on a 2" aperture are B=28.45, V=27.74, R=27.80, i'=27.43, z'=26.62, NB816=26.63, and NB921=26.54 in the AB system. The object catalog constructed for each of the seven bands contains more than 10^5 objects. The galaxy number counts corrected for detection incompleteness and star count contribution are found to be consistent with previous results in the literature. The mosaicked images and catalogs of all the bands have been made open to the public on Oct. 1, 2004 on the SDF project website at http://soaps.naoj.org/sdf/.

We report on the discovery of a star-forming galaxy at z=6.33 in the Subaru Deep Field. This object is selected as a candidate of an i'-dropout, high-redshift galaxy around z=6 because of its red i'-z' color in our deep optical imaging survey in the Subaru Deep Field. Our follow up optical spectroscopy reveals that this object is a strong Ly-alpha emitter with only very faint ultraviolet continuum. The rest-frame equivalent width of the detected Ly-alpha emission is as much as 130 A. Thus the light detected in our z' image is largely attributed to the Ly-alpha emission, i.e., ~40% of the z'-band flux is the strong Ly-alpha emission, giving a very red i'-z' color. This is consistent with the photometric property of this object because the narrow-band data obtained with the NB921 filter shows a significant depression, z'-NB921 = -0.54 mag. By using the photometric data, we show that some other objects among the 48 i'-dropout high-redshift galaxy candidates found in the Subaru Deep Field also show a significant NB921 depression. We briefly discuss the nature of these NB921-depressed objects.

We report the properties of the 35 robust candidates of Ly-alpha blobs (LABs), which are larger than 16 arcsec^2 in isophotal area and brighter than 0.7 x 10^-16 ergs s^-1 cm^-2, searched in and around the proto-cluster region at redshift z=3.1 discovered by Steidel et al. in the SSA22 field, based on wide-field (31'x23') and deep narrow-band (NB497; 4977/77) and broad-band (B,V, and R) images taken with the prime-focus camera on the Subaru telescope. The two previously known giant LABs are the most luminous and the largest ones in our survey volume of 1.3 x 10^5 Mpc^3. We revealed the internal structures of the two giant LABs and discovered some bubble-like features, which suggest that intensive starburst and galactic superwind phenomena occurred in these objects in the past. The rest 33 LABs have isophotal area of about 16-78 arcsec^2 and flux of 0.7-7 x 10^-16 ergs s^-1 cm^-2. These 35 LABs show a continuous distribution of isophotal area and emission line flux. The distributions of average surface brightness and morphology are widespread from relatively compact high surface brightness objects to very diffuse low surface brightness ones. The physical origins of these LABs may be (i) photo-ionization by massive stars, or active galactic nuclei, or (ii) cooling radiation from gravitationally heated gas, or (iii) shock heating by starburst driven galactic superwind. One third of them are apparently not associated with ultra-violet continuum sources that are bright enough to produce Ly-alpha emission, assuming a Salpeter initial mass function. The 90% of these LABs are located inside the high surface density region of the 283 relatively compact and strong Ly-alpha emitters selected in our previous study. This suggests that these LABs may be the phenomena related to dense environment at high redshift.

We investigate the clustering properties of 2,600 Lyman Break Galaxies (LBGs) at z=3.5-5.2 in two large blank fields, the Subaru Deep Field and the Subaru/XMM Deep Field (600arcmin^2 each). The angular correlation functions of these LBGs show a clear clustering at both z~4 and 5. The correlation lengths are r_0= 4.1^{+0.2}_{-0.2} and 5.9^{+1.3}_{-1.7} h_{100}^{-1} Mpc (r_0= 5.1^{+1.0}_{-1.1} and 5.9^{+1.3}_{-1.7} h_{100}^{-1} Mpc) for all the detected LBGs (for L>L* LBGs) at z~4 and 5, respectively. These correlation lengths correspond to galaxy-dark matter biases of b_g= 2.9^{+0.1}_{-0.1} and 4.6^{+0.9}_{-1.2} (b_g=3.5^{+0.6}_{-0.7} and 4.6^{+0.9}_{-1.2}), for all the detected LBGs (for L>L^* LBGs) at z~4 and 5, respectively. These results, combined with estimates for z~3 LBGs in the literature, show that the correlation length of L>L^* LBGs is almost constant, ~5 h_{100}^{-1} Mpc, over z~3-5, while the bias monotonically increases with redshift at z>3. We also find that for LBGs at z~4 the clustering amplitude increases with UV-continuum luminosity and with the amount of dust extinction. We estimate the mass of dark halos hosting various kinds of high-z galaxies including LBGs with the analytic model given by Sheth & Tormen (1999). We find that the typical mass of dark halos hosting L>L^* LBGs is about 1x10^{12} h_{70}^{-1}Msol over z~3-5, which is comparable to that of the Milky Way Galaxy. A single dark halo with ~10^{12} h_{70}^{-1} Msol is found to host 0.1-0.3 LBG on average but host about four K-band selected galaxies.

(abridged) We investigate photometric properties of Lyman Break Galaxies (LBGs) at z=3.5-5.2 based on large samples of 2,600 LBGs detected in deep (i'~27) and wide-field (1,200 arcmin^2) images taken in the Subaru Deep Field (SDF) and the Subaru/XMM Deep Field (SXDF). The selection criteria for the LBG samples are examined with 85 spectroscopically identified objects and by Monte Carlo simulations. We find in the luminosity functions of LBGs (i) that the number density of bright galaxies (M_{1700}<-22; corresponding to SFR_{corr}>100 Msolar yr^{-1}) decreases significantly from z=4 to 5 and (ii) that the faint-end slope of the luminosity function may become steeper towards higher redshifts. We estimate dust extinction of z=4 LBGs with M<M^* from UV slopes, and obtain E(B-V)=0.15+/-0.03 as the mean value. The dust extinction remains constant with apparent luminosity, but increases with intrinsic luminosity. We find no evolution in dust extinction between LBGs at z=3 and 4. We investigate the evolution of UV-luminosity density at 1700A, rho, and find that rho does not significantly change from z=3 to z=5, i.e., rho(z=4)/rho(z=3)=1.0+/-0.2 and rho(z=5)/rho(z=3)=0.8+/-0.4, thus the cosmic star-formation rate (SFR) density remains constant. We find that the stellar mass density estimated from the cosmic SFR is consistent with those derived directly from the stellar mass function at z=0-1, but exceeds those at z~3 by a factor of 3. We find that the ratio of the UV-luminosity density of Ly-a emitters (LAEs) to that of LBGs is ~60% at z=5, and thus about a half of the star formation at z=5 probably occurs in LAEs. We obtain a constraint on the escape fraction of UV-ionizing photons produced by LBGs, f_{esc}>~0.13.

We reported in a previous paper the discovery of large-scale structure of Lyman Alpha emitters (LAEs) at z=4.86+-0.03 with a projected size of 20 Mpc x 50 Mpc in narrow-band data of a 25' x 45' area of the Subaru Deep Field (Omega_0=0.3, lambda_0=0.7, H0=70 km/s/Mpc). However, the surveyed area, which corresponds to 55 Mpc x 100 Mpc, was not large enough that we can conclude that we are seeing a typical distribution of z~5 LAEs. In this Letter, we report the results of follow-up imaging of the same sky area using a new narrow-band filter (NB704, lambda_c=7046 A and FWHM=100 A) to detect LAEs at z=4.79, i.e., LAEs lying closer to us by 39 Mpc on average than the z=4.86 objects. We detect 51 LAEs at z=4.79+-0.04 down to NB704=25.7, and find that their sky distribution is quite different from the z=4.86 LAEs'. The clustering of z=4.79 LAEs is very weak on any scales and there is no large-scale high- contrast structure. The shape and the amplitude of the angular correlation function are thus largely different between the two samples. These results demonstrate a large cosmic variance in the clustering properties of LAEs on scales of ~ 50 Mpc.

We present an optical multicolor-imaging study of the galaxy cluster Abell 521 at $z = 0.25$, using Suprime-Cam on the Subaru Telescope, covering an area of $32 \times 20$ arcmin$^2$ ($9.4 \times 5.8 h_{50}^{-2}$ Mpc$^2$ at $z = 0.25$). Our imaging data taken with both a narrow-band filter, $NB816$ ($\lambda_0 = 8150$\AA and $\Delta \lambda = 120$\AA), and broad-band filters, $B,V,R_{\rm C}, i^\prime$, and $z^\prime$ allow us to find 165 H$\alpha$ emitters. We obtain the H$\alpha$ luminosity function (LF) for the cluster galaxies within 2 Mpc; the Schechter parameters are $\alpha = -0.75 \pm 0.23$, $\phi^\star = 10^{-0.25 \pm 0.20}$ Mpc$^{-3}$, and $L^\star = 10^{42.03 \pm 0.17}$ erg s$^{-1}$. Although the faint end slope, $\alpha$, is consistent with that of the local cluster H$\alpha$ LFs, the characteristic luminosity, $L^\star$, is about 6 times (or $\approx 2$ mag) brighter. This strong evolution implies that Abell 521 contains more active star-forming galaxies than the local clusters, being consistent with the observed Butcher-Oemler effect. However, the bright $L^\star$ of Abell 521 may be, at least in part, due to the dynamical condition of this cluster.

We report here the results of deep optical spectroscopy of the very extended emission-line region (VEELR) found serendipitously around the Seyfert 2 galaxy NGC 4388 in the Virgo cluster using the Subaru Telescope. The H-alpha recession velocities of most of the filaments of the region observed are highly blue-shifted with respect to the systemic velocity of the galaxy. The velocity field is complicated, and from the kinematic and morphological points of view, there seem to be several streams of filaments: low velocity (v = -100 km/s) filaments, high velocity (v = -300 km/s) filaments, and a very high velocity (v > -500 km/s) cloud complex. The emission-line ratios of the VEELR filaments are well explained by power-law photoionization models with solar abundances, suggesting that the Seyfert nucleus of NGC 4388 is the dominant ionization source of the VEELR and that the VEELR gas has moderate metallicity. In addition to photoionization, shock heating probably contributes to the ionization of the gas. In particular, the filaments outside the ionization cone of the Seyfert nucleus are mainly excited by shocks. We conclude that the VEELR was formerly the disk gas of NGC 4388, which has been stripped by ram pressure due to the interaction between the hot intra-cluster medium (ICM) and the galaxy. The velocity field and the morphology of the VEELR closely resemble snapshots from some numerical simulations of this process. In the case of NGC 4388, the ram pressure-stripped gas, which is normally seen as extended HI filaments, happens to be exposed and ionized by the radiation from the AGN, and so can be seen as optical emission-line gas.

We explore empirical relations between three different populations of high-redshift galaxies and their hosting dark halos employing the halo model approach. Specifically we consider Lyman-break galaxies at z\sim4 and at z\sim5, and Lyman-Alpha emitters at z\simeq 4.86. We adopt a halo occupation function prescription to parameterize the properties of their hosting halos and the efficiency of halo-dependent star formation. We find that the two LBG samples are well described by the halo model with an appropriate HOF. We obtain constraints on properties of their hosting halos. A typical mass of hosting halos for LBGs is 5\times10^{11}h^{-1}M_\odot and the expected number of LBGs per halo is \sim0.5, therefore there is an approximate one-to-one correspondence between halos and LBGs. We also find a sign of the minimum mass of LBG hosting halos decreasing with time. We discuss implications of these findings on the star formation history of LBGs. On the other hand, for LAEs, our simple HOF prescription fails to reproduce simultaneously the observed angular correlation function and the number density. This might imply either that the distribution of LAEs within hosting halos differs from that of dark matter, or that the strong large-scale correlation is due to the existence of an unusual, large overdense region, and so the survey region is not a representative of the z\sim5 universe, the definite answer should wait for a much wider survey of LAEs at high redshifts.

The distribution of early-type galaxy velocity dispersions, phi(sigma), is measured using a sample drawn from the SDSS database. Its shape differs significantly from that which one obtains by simply using the mean correlation between luminosity, L, and velocity dispersion, sigma, to transform the luminosity function into a velocity function: ignoring the scatter around the mean sigma-L relation is a bad approximation. An estimate of the contribution from late-type galaxies is also made, which suggests that phi(sigma) is dominated by early-type galaxies at velocities larger than ~ 200 km/s.

We present measurements of the cosmic shear correlation in the shapes of galaxies in the Suprime-Cam 2.1 deg^2 R_c-band imaging data. As an estimator of the shear correlation originated from the gravitational lensing, we adopt the aperture mass variance. We detect a non-zero E mode variance on scales between 2 and 40arcmin. We also detect a small but non-zero B mode variance on scales larger than 5arcmin. We compare the measured E mode variance to the model predictions in CDM cosmologies using maximum likelihood analysis. A four-dimensional space is explored, which examines sigma_8, Omega_m, Gamma and zs (a mean redshift of galaxies). We include three possible sources of error: statistical noise, the cosmic variance estimated using numerical experiments, and a residual systematic effect estimated from the B mode variance. We derive joint constraints on two parameters by marginalizing over the two remaining parameters. We obtain an upper limit of Gamma<0.5 for zs>0.9 (68% confidence). For a prior Gamma\in[0.1,0.4] and zs\in[0.6,1.4], we find sigma_8=(0.50_{-0.16}^{+0.35})Omega_m^{-0.37} for flat cosmologies and sigma_8=(0.51_{-0.16}^{+0.29})Omega_m^{-0.34}$ for open cosmologies (95% confidence). If we take the currently popular LCDM model, we obtain a one-dimensional confidence interval on sigma_8 for the 95.4% level, 0.62<\sigma_8<1.32 for zs\in[0.6,1.4]. Information on the redshift distribution of galaxies is key to obtaining a correct cosmological constraint. An independent constraint on Gamma from other observations is useful to tighten the constraint.

We present further observations of the extremely red object ERO J094258+4659.2, identified by \citet{iye00} as ERO R1 in their deep images of the cluster A851. We estimate its redshift independently by eight-band photometric redshift determination and cross-correlation of a new H-band spectrum with the optical spectra of local E/S0 galaxies, and conclude that it lies at $z \sim 1.5$. Although its colors are consistent both with an elliptical galaxy and an S0 galaxy at that redshift, its elongated shape and exponential luminosity profile suggest the presence of an evolved stellar disk component. We rule out the possibility that these properties are strongly influenced by gravitational lensing by the foreground cluster, and therefore conclude that this object is more likely to be an S0-like galaxy, rather than a lensed elliptical. The H-band spectrum does not show strong H$\alpha$ emission and the star formation rate therefore appears to be very modest. The presence of such a galaxy with an apparently relaxed disk of stars at this high redshift provides a new and strong constraint on theoretical models which aim to explain the formation and evolution of galaxies.

We report the discovery of a large-scale structure of Lyman alpha emitters (LAEs) at z=4.86 based on wide-field imaging with the prime-focus camera (Suprime-Cam) on the Subaru telescope. We observed a 25' x 45' area of the Subaru Deep Field in a narrow band (NB, lambdaC=7126 A and FWHM=73 A) together with R and i'. We isolate from these data 43 LAE candidates down to NB=25.5 mag using color criteria. Follow-up spectroscopy of five candidates suggests the contamination by low-z objects to be ~ 20%. We find that the LAE candidates are clustered in an elongated region on the sky of 20 Mpc in width and 50 Mpc in length at z=4.86, which is comparable in size to present-day large-scale structures (we adopt H0=70 km/s/Mpc, Omega0=0.3, lambda0=0.7). This elongated region includes a circular region of 12 Mpc radius of higher surface overdensity (delta=2), which may be the progenitor of a cluster of galaxies. Assuming this circular region to be a sphere with a spatial overdensity of 2, we compare our observation with predictions by Cold Dark Matter models. We find that an Omega0=0.3 flat model with sigma8=0.9 predicts the number of such spheres consistent with the observed number (one sphere in our survey volume) if the bias parameter of LAEs is b \simeq 6. This value suggests that the typical mass of dark haloes hosting LAEs at z \simeq 5 is of the order of 10^12 Msolar. Such a large mass poses an interesting question about the nature of LAEs.

The morphological dependence of the luminosity function is studied using a sample containing approximately 1500 bright galaxies classified into Hubble types by visual inspections for a homogeneous sample obtained from the Sloan Digital Sky Survey (SDSS) northern equatorial stripes. Early-type galaxies are shown to have a characteristic magnitude by 0.45 mag brighter than spiral galaxies in the $r^{\ast}$ band, consistent with the `universal characteristic luminosity' in the $B$ band. The shape of the luminosity function differs rather little among different morphological types: we do not see any symptoms of the sharp decline in the faint end for the luminosity function for early-type galaxies at least 2 mag fainter than the characteristic magnitude, although the faint end behaviour shows a slight decline ($\alpha\lsim -1$) compared with the total sample. We also show that a rather flat faint end slope for early-type galaxies is not due to an increasing mixture of the dwarf galaxies which have softer cores. This means that there are numerous faint early-type galaxies with highly concentrated cores.

We have built an 80-mega pixels (10240 X 8192) mosaic CCD camera, called Suprime-Cam, for the wide-field prime focus of the 8.2 m Subaru telescope. Suprime-Cam covers a field of view 34' X 27', a unique facility among the the 8-10 m class telescopes, with a resolution of 0.''202 per pixel. The focal plane consists of ten high-resistivity 2k X 4k CCDs developed by MIT Lincoln Laboratory, which are cooled by a large stirling-cycle cooler. The CCD readout electronics was designed to be scalable, which allows the multiple read-out of tens of CCDs. It takes 50 seconds to readout entire arrays. We designed a filter-exchange mechanism of the jukebox type that can hold up to ten large filters (205 X 170 X 15 mm^3). The wide-field corrector is basically a three-lens Wynne-type, but has a new type of atmospheric dispersion corrector. The corrector provides a flat focal plane and an un-vignetted field of view of 30' in diameter. The achieved co-planarity of the focal array mosaic is smaller than 30 micron peak-to-peak, which realizes mostly the seeing limited image over the entire field. The median seeing in the Ic-band, measured over one year and a half, is 0.''61. The PSF anisotropy in Suprime-Cam images, estimated by stellar ellipticities, is about 2% under this median seeing condition. At the time of commissioning, Suprime-Cam had the largest survey speed, which is defined as the field of view multiplied by the primary mirror area of the telescope, among those cameras built for sub-arcsecond imaging.

We have identified 38 candidates of intracluster planetary nebulae (ICPNe) in a $34'\times27'$ field in the core of the Virgo cluster based on the Suprime-Cam imaging through two narrow-band filters centered at the redshifted wavelengths of the [OIII]$\lambda=5007\AA$ and the Halpha \lambda=6563\AA lines. Broad-band images in V and R bands are used to check for any emissions in the adjacent continuum. We describe the method briefly and present the list of intracluster planetary nebulae candidates, together with their finding charts. The ICPN candidates show a highly inhomogeneous distribution, which may suggest an association with the M86-M84 subcluster. Fraction of diffuse intracluster light with respect to total light in galaxies is estimated to be about 10%, leading to an estimate of about 20% for the baryon fraction. Spectroscopic follow up and a wider survey are critical to reveal the nature of intracluster stellar population.

We have identified intracluster planetary nebulae in a Virgo cluster core field by imaging with the Subaru Suprime-Cam through two narrow band filters centered at the redshifted wavelengths of the [OIII] $\lambda=5007$ \AA and the H$\alpha$ $\lambda = 6563$ \AA lines; broad-band images in V and R were acquired to check for emission in the adjacent continuum. Emission line objects in Virgo are then selected from the two-color diagram ([OIII] -- Halpha) vs. ([OIII] -- (V+R)), which is calibrated using PNe in M84 (Jacoby et al. 1990). Using both [OIII] and Halpha allows us to distinguish bona-fide planetary nebulae from high redshift emission-line galaxies at the bright end of the [OIII] luminosity function. Spectroscopic observations of a subsample of these objects were made at the TNG and at the VLT, in a region around M84 and in an intracluster field respectively. The observations confirm the efficiency of the combined [OIII]+H$\alpha$ imaging to identify true PNe. We also obtained the first spectrum of an intracluster PN which shows the [OIII] doublet with S/N > 10 and its Halpha emission. From the results based on the spectroscopic follow-up, we derive a lower limit to the fraction of the Virgo cluster light contributed by the intracluster stars at the surveyed position in the cluster core: it amounts to 10%.

We study the evolution of elliptical galaxies at z >~ 1 on the basis of a sample of 247 Extremely Red Objects (EROs) with R-Ks >= 3.35 (AB) and Ks <= 22.1 (AB) constructed from BVRi'z'JHKs multicolor data of a 114 arcmin2 area in the Subaru/XMM Deep Survey Field. By fitting template spectra of old ellipticals (OEs) and young, dusty starbursts (DSs) to the multicolor data, we classify EROs into these two classes and estimate their redshifts. We find that 58% of the EROs in our sample belong to the OE class and that 24% of the OEs are fit by a spectrum having a disk component with the B-band bulge-to-total luminosity ratio of <= 0.9. We compare the observed surface density of OEs with predictions from passive evolution models. We then derive rest-frame B-band luminosity functions (LFs) of OEs in our sample at z=1-1.5 and 1.5-2.5. We find that the LF at z=1-1.5 roughly agrees with the LF of local ellipticals if a dimming of 1.3 mag from z=1.25 to the present epoch is assumed. On the other hand, the amplitude of the LF at z=1.5-2.5 is found to be lower than that of the local LF by a factor of ~ 3 over the whole range of magnitude observed. Taking account of a strong decrease in the number density of morphologically classified early-type galaxies at z >~ 1.5 found by several authors, we conclude that the majority of ellipticals seen at present have not established either a red color or a smooth 1/4-law profile before z ~ 1.5$. (abridged)

We report the first result of weak gravitational lensing survey on a 2.1 sq deg Rc-band image taken with a wide field camera (Suprime-Cam) on the prime focus of 8.2 m Subaru Telescope. The weak lensing mass reconstruction is applied to the data to search for dark matter halos of cluster scale; M >= 10^14 solar mass. The reconstructed convergence field is divided by 1-sigma noise to obtain the signal-to-noise ratio map (S/N-map) of the detection. Local maxima and minima are searched on the S/N-map and the probability distribution function (PDF) of the peaks are created to compare with model predictions. We found excess over noise PDF created from the randomized realization on both positive and negative sides. Negative peaks imply the presence of voids in the dark matter distribution and this is the first report of the detection. Positive peaks, on the other hand, represent the dark matter halos and the number count of the halos on the 2.1 sq deg image is 4.9 +- 2.3 for S/N > 5 where the Gaussian smoothing radius of the convergence map is 1'. The result is consistent with the prediction assuming the Press-Schechter mass function and the NFW halo profile under the cluster normalized CDM cosmology. Although the present statistics is not enough due to the limited field of view, this work demonstrates that dark matter halo count via weak lensing is a promising way to test the paradigm of structure formation and cosmological model.

We have constructed a very deep K'-selected multicolor BVRIz'JK' sample of 439 field galaxies. Photometric redshifts for sample galaxies were estimated. The overall redshift distribution N(z) for the K'<21.0 sample is consistent with previous observations, and for the first time we derive the N(z) down to K'=24.0. After taking account of the dust extinction and selection effects of the sample, the observed N(z) distribution is well described with the PLE model, while the hierarchical galaxy formation model shows an apparent deficiency of galaxies especially at z>2. The photometric redshift and the best-fit SED model evaluations allow the derivation of the rest-frame K', B, and UV(2000A) luminosity functions (LFs) and their evolutions. The rest-K' LF shows almost no evolution up to z=3, while the rest B LF shows mild luminosity evolution, and the rest UV LF shows strong luminosity evolution. This trend corresponds to the evolution of the rest UV-K' color, which gets bluer with increasing redshift. We also found that more massive galaxies are redder in this rest-frame color in every epoch. The rest-UV LF of our K'-selected galaxies shows a much shallower faint end slope at z=3 than that of previous estimations for rest-UV selected Lyman break galaxies. As a consequence, the contribution to the global star formation rate of our K'-selected galaxies is about 42% of that derived from the integration of LF of Lyman break galaxies at z=3. This result suggests that a large fraction of the star formation rate density at z>1.5 comes from the contribution from the faint (M_2000A>-20) blue galaxy population at high redshift universe that have not yet obviously been identified. (abridged)

Using a catalog of 147,986 galaxy redshifts and fluxes from the Sloan Digital Sky Survey (SDSS) we measure the galaxy luminosity density at z=0.1 in five optical bandpasses corresponding to the SDSS bandpasses shifted to match their restframe shape at z=0.1. We denote the bands {0.1}{u}, {0.1}{g}, {0.1}{r}, {0.1}{i}, {0.1}{z}, with \lambda_{eff} = [3216, 4240, 5595, 6792, 8111] Angstroms respectively. We use a maximum likelihood method which allows for a general form for the shape of the luminosity function, simple luminosity and number evolution, incorporates flux uncertainties, and accounts for the flux limits of the survey. We find luminosity densities at z=0.1 in absolute AB magnitudes in a Mpc^3 of [-14.10 \pm 0.15, -15.18 \pm 0.03, -15.90 \pm 0.03, -16.24 \pm 0.03, -16.56 \pm 0.02] in [{0.1}{u}, {0.1}{g}, {0.1}{r}, {0.1}{i}, {0.1}{z}], respectively, using \Omega_0 =0.3, \Omega_\Lambda=0.7, and h=1, and using Petrosian magnitudes. Similar results are obtained using Sersic model magnitudes, suggesting that flux from outside the Petrosian apertures is not a major correction. In the {0.1}{r} band, the best fit Schechter function to our results has \phi_\ast = (1.49 \pm 0.04) \times 10^{-2} h^3 Mpc^{-3}, M_\ast - 5\log_{10} h = -20.44 \pm 0.01, and \alpha = -1.05\pm 0.01. In solar luminosities, the luminosity density in {0.1}{r} is (1.84 \pm 0.04) h 10^8 L_{{0.1}{r},\odot} Mpc^{-3}. Our results are consistent with other estimates of the luminosity density, from the Two-degree Field Galaxy Redshift Survey and the Millenium Galaxy Catalog. They represent a substantial change (\sim 0.5 mag) from earlier SDSS luminosity density results based on commissioning data, almost entirely because of the inclusion of evolution in the luminosity function model.

We present the results of a survey for emission-line objects based on optical intermediate-band ($\lambda_{\rm c}$ = 5736 \AA ~ and $\Delta\lambda$ = 280 \AA) and broad-band ($B$, $V$, $R$, and $i^\prime$) observations of the Subaru/XMM-Newton Deep Field on the 8.2 m Subaru telescope with the Subaru Prime Focus Camera, Suprime-Cam. All the data were obtained during the guaranteed time observations of the Suprime-Cam instrument. The intermediate-band image covered a sky area with 10\minpoint62 $\times 12\minpoint40 \approx 132$ arcmin$^2$ in the Subaru/XMM-Newton Deep Field (Ouchi et al.). Using this image, we have found 23 emission-line sources whose observed emission-line equivalent widths are greater than 250 \AA. Their optical multicolor properties indicate that six emission-line sources are Ly$\alpha$ emitters at $z \approx$ 3.7 ($\Delta z \approx 0.22$). They are either intense starburst galaxies or active galactic nuclei like quasars at $z \approx$ 3.7. Two more emission-line sources may also be Ly$\alpha$ emitters at $z \approx$ 3.7 although their multicolor properties are marginal. Among the remaining 15 emission-line objects, eight objects appear strong emission-line galaxies at lower redshift; e.g., [O {\sc ii}] $\lambda$3727 emitters at $z \approx 0.54$, H$\beta$ at $z \approx 0.18$, or [O {\sc iii}]$\lambda$5007 emitters at $z \approx 0.15$. The remaining seven objects are unclassified because they are too faint to be detected in broad-band images. We discuss observational properties of these strong emission-line sources. In particular, our data allow us to estimate the star formation density at $z \approx 3.7$ for the first time.

Using photometry and spectroscopy of 144,609 galaxies from the Sloan Digital Sky Survey, we present bivariate distributions of pairs of seven galaxy properties: four optical colors, surface brightness, radial profile shape as measured by the Sersic index, and absolute magnitude. In addition, we present the dependence of local galaxy density (smoothed on 8 h^{-1} Mpc scales) on all of these properties. Several classic, well-known relations among galaxy properties are evident at extremely high signal-to-noise ratio: the color-color relations of galaxies, the color-magnitude relations, the magnitude-surface brightness relation, and the dependence of density on color and absolute magnitude. We show that most of the i-band luminosity density in the universe is in the absolute magnitude and surface brightness ranges used. Some of the relationships between parameters, in particular the color--magnitude relations, show stronger correlations for exponential galaxies and concentrated galaxies taken separately than for all galaxies taken together. We provide a simple set of fits of the dependence of galaxy properties on luminosity for these two sets of galaxies.

We report on early results from deep and wide-field narrow-band imaging on a 543 arcmin^2 area of the Subaru Deep Field. We find 87 Lya emitters (LAEs) at z=4.86(+/-0.03) which are photometrically selected by a combination of two broad bands (R and i') and one narrow band (NB711; 7126A FWHM=73A). We derive the luminosity functions (LFs) of the LAEs at Lya luminosity and at UV-continuum (1700A) luminosity. The LFs show little evolution between z=3.4 and z=4.86 either in Lya or UV-continuum emission. The amplitude of the LAE LF tends to decline at the bright magnitudes more rapidly than that of the LBG LF at similar redshifts. We calculate the angular correlation function of our LAEs up to ~15 arcmin separations. The angular correlation w(theta) is found to increase with decreasing angular separations, showing a clear signal of clustering. It is also found that the distribution of LAEs shows a large density gradient with a scale of >~15 arcmin, which would indicate the existence of a large-scale structure of LAEs on >~20h^(-1) Mpc scales.We fit the observed correlation function by Aw theta^(-0.8) to find Aw=29arcsec^(0.8). The estimated correlation length is r0=3.5(+/-0.3)h^(-1)Mpc in comoving units (Omega_m=0.3 and Omega_Lambda=0.7), which is slightly larger than the value for z~4 LBGs with i'<26. We calculate the angular correlation function for two sub-samples of the 87 LAEs divided by Lya luminosity, UV-continuum luminosity, and Lya equivalent width (EW). The Lya bright sub-sample shows a larger correlation amplitude than the Lya faint sub-sample, while no significant difference is found for the sub-samples divided by UV-continuum luminosity or EW.

We present the 158 standard stars that define the u'g'r'i'z' photometric system. These stars form the basis for the photometric calibration of the Sloan Digital Sky Survey (SDSS). The defining instrument system and filters, the observing process, the reduction techniques, and the software used to create the stellar network are all described. We briefly discuss the history of the star selection process, the derivation of a set of transformation equations for the UBVRcIc system, and plans for future work.

We found a very large, 35 kpc, emission-line region around the Seyfert type 2 galaxy NGC 4388, using deep narrow-band imaging with the prime focus camera (Suprime-Cam) of the Subaru telescope. This region consists of many faint gas clouds or filaments, and extends northeastwards from the galaxy. The typical H-alpha luminosity L(H-alpha) of the clouds is 10^37 erg s^-1, and the total L(H-alpha) of the region within 10 kpc from the nucleus is 2 x 10^38 erg s^-1, which corresponds to an ionized gas mass of 10^5 M_solar. The map of the emission-line intensity ratio I([OIII])/I(H-alpha) indicates that the inner (r < 12 kpc) region of the VEELR may be excited by nuclear ionizing radiation. The excitation mechanism of the outer (r > 12 kpc) region is unclear, but it is likely that the nuclear radiation is also a dominant source of its ionization. We discuss the origin of the ionized gas. Two plausible origins of the gas in the VEELR are (i) the tidal debris resulting from a past interaction with a gas-rich dwarf galaxy, i.e., a minor merger, or (ii) the interstellar medium of NGC 4388, stripped by the ram pressure of the hot intracluster medium of the Virgo cluster.

We investigate the radial dependence of the spectroscopic properties, in particular the Mg2, <Fe> and H beta spectroscopic indices, in a sample of galaxies spanning a wide range of absolute luminosity in the Coma cluster. After allowing for the magnitude dependence of these indices, we find a significant gradient in Mg2, in the sense that galaxies in the core of the cluster have stronger Mg2. We find only weak gradients in <Fe> and H beta. Using the model grids presented in an earlier paper in this series, we attribute the Mg2 gradient to changes in metal abundance. One possible mechanism to create this abundance gradient is pressure confinement by the intracluster medium of material from Supernova driven winds early in the history of the galaxies.

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

(abr.) We present the analysis of the spectroscopic catalog of galaxies in the Coma cluster from Mobasher et al. This catalog comprises ~300 spectra of cluster members with magnitudes M_B=-20.5 to -14 in two areas of ~1 X 1.5 Mpc towards the center and the SW region of the cluster. In this paper the analysis is restricted to the 257 galaxies with no emission lines in their spectra. The strength of the age-sensitive indices is found to correlate with galaxy magnitude over the whole magnitude range. Similarly, the metallicity-sensitive indices anticorrelate with magnitude. By comparing the observed indices with model grids based on the Padova isochrones, we derive luminosity-weighted ages and metallicities. We present the distributions of ages and metallicities for galaxies in various magnitude bins. The mean metallicity decreases with galaxy magnitude and, at a given luminosity, appears to be generally lower for galaxies in the SW region of Coma as compared to the center of the cluster. A broad range of ages, from younger than 3 Gyr to older than 9 Gyr, is found in galaxies of any magnitude. However, systematic trends of age with luminosity are present among galaxies in the central field. In the central Mpc of Coma, a large fraction of galaxies at any luminosity (50-60% of the giants, >30 % of the dwarfs) show no evidence in their central regions of star formation occurred at redshift z<2, while the proportion of galaxies with significant star formation occurring at intermediate (0.35<z<2) and low (z<0.35) redshifts is found to depend on galaxy luminosity. Faint galaxies with young luminosity-weighted ages appear to have a bimodal metallicity distribution that, if confirmed, would point to a composite formation scenario. Finally, we present an interpretation of the index-magnitude relations observed.

We study the luminosity function and the correlation function of about 1200 z~4 Lyman break galaxies (LBGs) with i'<26 that are photometrically selected from deep BRi' imaging data of a 618 arcmin^2 area in the Subaru/XMM-Newton Deep Field taken with Subaru Prime Focus Camera. The contamination and completeness of our LBG sample are evaluated, on the basis of the Hubble Deep Field-North (HDF-N) objects, to be 17% and 45%, respectively. We derive the UV (rest 1700A) luminosity functions (LFs) and find a large population of UV-luminous galaxies at z~4. The LFs of the red and blue subsamples imply that the bright LBGs are redder in the UV continuum than the average color of the LBGs. Then we calculate the correlation function over theta = 2''-1000'' and find that it is fitted fairly well by a power law, omega(theta)=A_omega theta^(-0.8), with A_omega=0.71 +/- 0.26. We estimate the correlation length r_0 (in comoving units) of the two-point spatial correlation function xi(r) = (r/r_0)^(-1.8) to be r_0=2.7 +0.5/-0.6 h^(-1) Mpc (Omega_m=0.3 and Omega_Lambda=0.7). The correlation function shows an excess of omega (theta) on small scales (theta < 5''), departing from the power-law fit at > 3 sigma significance level. Interpreting this as being due to galaxy mergers, we evaluate the fraction of galaxies undergoing mergers to be 3.0 +/- 0.9%, which is significantly smaller than those of galaxies at intermediate redshifts.

The ages of stellar populations in 52 elliptical and S0 galaxies in the Coma cluster are investigated using a new spectroscopic sample. More than 40% of the S0s are found to have experienced recent star formation in their central regions during the last ~5 Gyrs, while such activity is absent in the ellipticals. Galaxies in this sample have absolute magnitudes in the range -20.5 < M_B < -17.5, and the fraction of S0 galaxies with recent star formation is higher at fainter luminosities. The observed luminosity range of S0 galaxies with signs of recent star formation activity is consistent with them being the descendants of typical star-forming spirals at intermediate redshift whose star formation has been halted as a consequence of the dense environment.

We study the clustering properties of about 1200 z~4 Lyman Break Galaxy (LBG) candidates with i'<26 which are selected by color from deep BRi' imaging data of a 618 arcmin^2 area in the Subaru/XMM-Newton Deep Field taken with Subaru Prime Focus Camera. The contamination and completeness of our LBG sample are evaluated, on the basis of the Hubble Deep Field North (HDFN) objects, to be 17% and 45%, respectively. We derive the angular correlation function over theta = 2''-1000'', and find that it is fitted fairly well by a power law, omega(theta)=A_omega theta^{-0.8}, with A_omega = 0.71 +/- 0.26. We then calculate the correlation length r0 (in comoving units) of the two-point spatial correlation function xi(r) = (r/r0)^{-1.8} from A_omega using the redshift distribution of LBGs derived from the HDFN, and find r0=2.7 (+0.5/-0.6) h^{-1} Mpc in a Lambda-dominated universe (Omega_m=0.3 and Omega_Lambda=0.7). This is twice larger than the correlation length of the dark matter at z~4 predicted from an analytic model by Peacock & Dodds but about twice smaller than that of bright galaxies predicted by a semi-analytic model of Baugh et al. We find an excess of omega(theta) on small scales (theta < 5'') departing from the power law fit over 3 sigma significance levels. Interpreting this as due to galaxy mergers, we estimate the fraction of galaxies undergoing mergers in our LBG sample to be 3.0 +/- 0.9%, which is significantly smaller than those of galaxies at intermediate redshifts.

This is the second paper in a series studying the photometric and spectroscopic properties of galaxies of different luminosities in the Coma cluster. The sample selection, spectroscopic observations and completeness functions are presented here. To study the spectral properties of galaxies as a function of their local environment, two fields were selected for spectroscopic observations to cover both the core (Coma1) and outskirts (ie. south-west of the core and centered on NGC4839)- (Coma3) of the cluster. Medium resolution spectroscopy (6-9 \AA) was carried out for a total of 490 galaxies in both fields (302 in Coma1 and 188 in Coma3), using the WYFFOS multi-fiber spectrograph on the William Herschel Telescope. The galaxies cover a range of $12 < R < 20$, corresponding to -23 < M_R < -15 (H0=65 km/sec/Mpc). The redshifts are measured with an accuracy of 100 km/sec. The spectral line strengths and equivalent widths are also measured for the same galaxies and analysed in Poggianti et al (2001- paper III). A total of 189 (Coma1) and 90 (Coma3) galaxies are identified as members of the Coma cluster. An analysis of the colors show that only two members of the Coma cluster in our sample have B-R > 2. The completeness functions for the spectroscopic sample is presented.

We present bright galaxy number counts in five broad bands ($u', g', r', i', z'$) from imaging data taken during the commissioning phase of the Sloan Digital Sky Survey (SDSS). The counts are derived from two independent stripes of imaging scans along the Celestial Equator, one each toward the North and the South Galactic cap, covering about 230 and 210 square degrees, respectively. A careful study is made to verify the reliability of the photometric catalog. For galaxies brighter than $r^* = 16$, the catalog produced by automated software is examined against eye inspection of all objects. Statistically meaningful results on the galaxy counts are obtained in the magnitude range $12 \le r^* \le 21$, using a sample of 900,000 galaxies. The counts from the two stripes differ by about 30% at magnitudes brighter than $r^*= 15.5$, consistent with a local $2\sigma$ fluctuation due to large scale structure in the galaxy distribution. The shape of the number counts-magnitude relation brighter than $r^* = 16$ is well characterized by $N \propto 10^{0.6m}$, the relation expected for a homogeneous galaxy distribution in a ``Euclidean'' universe. In the magnitude range $ 16 < r^* < 21$, the galaxy counts from both stripes agree very well, and follow the prediction of the no-evolution model, although the data do not exclude a small amount of evolution. We use empirically determined color transformations to derive the galaxy number counts in the $B$ and $I_{814}$ bands. We compute the luminosity density of the universe at zero redshift in the five SDSS bands and in the $B$ band. We find ${\cal L}_{B} = 2.4 \pm 0.4 \times 10^8L_\odot h $Mpc$^{-3}$, for a reasonably wide range of parameters of the Schechter luminosity function in the $B$ band.

We investigate the photometric properties of 456 bright galaxies using imaging data recorded during the commissioning phase of the Sloan Digital Sky Survey (SDSS). Morphological classification is carried out by correlating results of several human classifiers. Our purpose is to examine the statistical properties of color indices, scale lengths, and concentration indices as functions of morphology for the SDSS photometric system. We find that $u'-g'$, $g'-r'$, and $r'-i'$ colors of SDSS galaxies match well with those expected from the synthetic calculation of spectroscopic energy distribution of template galaxies and with those transformed from $UBVR_CI_C$ color data of nearby galaxies. The agreement is somewhat poor, however, for $i'-z'$ color band with a discrepancy of $0.1-0.2$ mag. With the aid of the relation between surface brightness and radius obtained by Kent (1985), we estimate the averages of the effective radius of early type galaxies and the scale length of exponential disks both to be 2.6 kpc for $L^*$ galaxies. We find that the half light radius of galaxies depends slightly on the color bands, consistent with the expected distribution of star-forming regions for late type galaxies and with the known color gradient for early type galaxies. We also show that the (inverse) concentration index, defined by the ratio of the half light Petrosian radius to the 90% light Petrosian radius, correlates tightly with the morphological type; this index allows us to classify galaxies into early (E/S0) and late (spiral and irregular) types, allowing for a 15-20% contamination from the opposite class compared with eye-classified morphology.

The results of deep-imaging observations of a candidate of an absorbed QSO at z=0.653, AX J131831+3341, are presented. AX J131831+3341 was found during the course of optical follow-up observations of the ASCA Large Sky Survey, and has an X-ray luminosity of 10^{45} erg s^{-1} (2-10keV), which corresponds to those of QSOs. Its optical spectrum shows no significant broad H_beta emission line, suggesting that the object is an absorbed QSO. Deep R and V band images reveal the presence of a point-like nucleus and an asymmetric extended component. The nuclear component has a blue color, and the optical magnitude is much fainter than that expected from the observed X-ray flux for typical type-1 AGNs. These photometric properties and the presence of broad MgII 2800A emission can be explained simultaneously if the observed nuclear light is dominated by scattered nuclear light, though there is a possibility that the nuclear component is a slightly absorbed nucleus if its intrinsic X-ray to optical flux ratio is the largest among X-ray selected AGNs. The size of the extended component, which is thought to be the host galaxy of the QSO, is larger than those of normal disk galaxies at z=0-0.75, and the absolute magnitude is similar to those of the brightest host galaxies of QSOs at redshifts smaller than 0.5. The V-R and R-I colors of the component are consistent with a 1 Gyr-old stellar population model without absorption.

Deep near-infrared images of a blank 2'x2' section of sky near the Galactic north pole taken by Subaru Telescope are presented. The total integration times of the J and K' bands are 12.1 hours and 9.7 hours, resulting in 5-sigma limiting magnitudes of 25.1 and 23.5 mag, respectively. The numbers of sources within these limiting magnitudes found with an automated detection procedure are 385 in the J band and 350 in K'. Based on photometric measurements of these sources, we present number count vs. magnitude relations, color vs. magnitude diagrams, size vs. color relationships, etc. The slope of the galaxy number count plotted against the AB magnitude scale is about 0.23 in the 22 to 26 AB magnitude range of both bands. The spatial number density of galaxies as well as the slopes in the faint-end region given by the Subaru Deep Field (SDF) survey is consistent with those given by HST-NICMOS surveys as expressed on the AB magnitude diagram. Several sources having very large J-K' color are found including a few K' objects without detection at J. In addition, a number of faint Galactic stars are also detected, most of which are assigned to M-subdwarfs, together with a few brown dwarf candidates.

We derive the ratio of total blue luminosity to total baryon mass, LB/Mb, for massive (Mgas at the Abell radius is \ge 1 \times 10^{13} h^{-2.5} \Msolar) clusters of galaxies up to z \simeq 1 from the literature. Twenty-two clusters in our sample are at z > 0.1. Assuming that the relative mix of hot gas and galaxies in clusters does not change during cluster evolution, we use LB/Mb to probe the star formation history of the galaxy population as a whole in clusters. We find that LB/Mb of clusters increases with redshift from LB/Mb=0.024 (solar units) at z = 0 to \simeq 0.06 at z=1, indicating a factor of 2-3 brightening (we assume H0=70 km/s/Mpc). This amount of brightening is almost identical to the brightening of the M/LB ratio of early-type galaxies in clusters at 0.02 \le z \le 0.83 reported by van Dokkum et al. (1998). We compare the observed brightening of LB/Mb with luminosity evolution models for the galaxy population as a whole, changing the e-folding time of star formation \tau by 0.1 \le \tau \le 5 Gyr and the formation redshift \zF by 2 \le \zF < \infty. We find that \tau=0.1 Gyr 'single burst' models with \zF \ge 3 and \tau=5 Gyr 'disk' models with arbitrary \zF are consistent with the observed brightening, while models with \tau=1-2 Gyr tend to predict too steep brightening. We also derive the ratio of blue luminosity density to baryon density for field galaxies, adopting \Omega_b h^2 = 0.02, and find that blue luminosity per unit baryon is similar in clusters and in fields up to z \simeq 1 within the observational uncertainties.

We report new improved photometric redshifts of 1048 galaxies in the Hubble Deep Field (HDF). A standard chi^2 minimizing method is applied to seven-color UBVIJHK photometry by Fernandez-Soto, Lanzetta, & Yahil (1999). We use 187 template SEDs representing a wide variety of morphology and age of observed galaxies based on a population synthesis model by Kodama & Arimoto (1997). We introduce two new recipes. First, the amount of the internal absorption is changed as a free parameter in the range of E(B-V)=0.0 to 0.5 with an interval of 0.1. Second, the absorption due to intergalactic HI clouds is also changed by a factor of 0.5, 1.0, and 1.5 around the opacity given by Madau (1995). The total number of template SEDs is thus 187x6x3=3,366, except for the redshift grid. The dispersion sigma_z of our photometric redshifts with respect to spectroscopic redshifts is sigma_z=0.08 and 0.24 for z<2 and z>2, respectively, which are smaller than the corresponding values (sigma_z=0.09 and 0.45) by Fernandez-Soto et al. Improvement is significant, especially in z>2. This is due to smaller systematic errors which are largely reduced mainly by including three opacities due to intergalactic HI clouds. We discuss redshift distribution N(z) and cosmic star formation rate based on our new photometric redshifts.

We investigate pure luminosity evolution models for early-type (elliptical and S0) galaxies (i.e., no number density change or morphology transition), and examine whether these models are consistent with observed number counts in the B, I and K bands and redshift distributions of two samples of faint galaxies selected in the I and K bands. The models are characterized by the star formation time scale $\tau_{SF}$ and the time $t_{gw}$ when galactic wind blows in addition to several other conventional parameters. We find the single-burst model ($\tau_{SF}$=0.1 Gyr and $t_{gw}$=0.353 Gyr), which is known to reproduce the photometric properties of early-type galaxies in clusters, is inconsistent with redshift distributions of early-type galaxies in the field environment due to overpredictions of galaxies at $z\gsim1.4$ even with strong extinction which is at work until $t_{gw}$. In order for dust extinction to be more effective, we change $\tau_{SF}$ and $t_{gw}$ as free parameters, and find that models with $\tau_{SF}\gsim0.5$ Gyr and $t_{gw}>1.0$ Gyr can be made consistent with both the observed redshift distributions and number counts, if we introduce strong extinction $(E(B-V)\geq1$ as a peak value). These results suggest that early-type galaxies in the field environment do not have the same evolutionary history as described by the single-burst model.

We show that the growth curves of galaxies can be used to determine their bulge and disk parameters and bulge-to-total luminosity ratios, in addition to their conventional asymptotic magnitudes, provided that the point spread function is accurately known and signal-to-noise ratio is modest (S/N$\gtrsim30$). The growth curve is a fundamental quantity that most future large galaxy imaging surveys will measure. Bulge and disk parameters retrieved from the growth curve will enable us to perform statistical studies of luminosity structure for a large number of galaxies.

We carry out CCD photometry of galaxies in the 5.25 square region centered on Coma cluster down to $M_R=-16.0$, beyond the limit of conventional morphological classification. We use the angular two-point correlation function as well as radial profiles in order to characterize the luminosity segregation. We find strong luminosity segregation for our total sample over the magnitude range of $-20 \leq M_R \leq -16$, which is not entirely accounted for in terms of the morphology-density relation that is known to exist only for bright galaxies. We use a single consistent parameter, the degree of luminosity concentration, to parameterize the morphology of galaxies over the wide magnitude range, where both giant and dwarf galaxies are included. Galaxies with high central concentration (HCC) show strong luminosity segregation, i.e. their clustering strength depends strongly on luminosity while those with low central concentration (LCC) show almost no luminosity segregation. Radial density profile shows that brighter HCC-type galaxies tend to more strongly concentrate near the cluster center while LCC-type galaxies do not show such a dependence on luminosity. We show that these results are tenable against the contamination by field galaxies and uncertainties in our method of classification.

We describe an objective and automated method for detecting clusters of galaxies from optical imaging data. This method is a variant of the so-called `matched-filter' technique pioneered by Postman et al. (1996). With simultaneous use of positions and apparent magnitudes of galaxies, this method can, not only find cluster candidates, but also estimate their redshifts and richnesses as byproducts of detection. We examine errors in the estimation of cluster's position, redshift, and richness with a number of Monte Carlo simulations. No systematic discrepancies between the true and estimated values are seen for either redshift or richness. Spurious detection rate of the method is about less than 10% of those of conventional ones which use only surface density of galaxies. A cluster survey in the North Galactic Pole is executed to verify the performance characteristics of the method with real data. Two known real clusters are successfully detected. We expect these methods based on `matched-filter' technique to be essential tools for compiling large and homogeneous optically-selected cluster catalogs.

(Abridged) A simple quantitative model is presented for the history of galaxies to explain galaxy number counts, redshift distributions and some other related observations. We first infer that irregular galaxies and the disks of spiral galaxies are young, probably formed at $z\approx 0.5-2$ from a simultaneous consideration of colours and gas content under a moderate assumption on the star formation history. Assuming that elliptical galaxies and bulges of spiral galaxies, both called spheroids in the discussion, had formed early in the universe, the resulting scenario is that spiral galaxies formed as intergalactic gas accreting onto pre-existing bulges mostly at $z\approx 1-2$; irregular galaxies as seen today formed by aggregation of clouds at $z\approx 0.5-1.5$. Taking the formation epochs thus estimated into account, we construct a model for the history of galaxies employing a stellar population synthesis model. We assume that the number of galaxies does not change except that some of them (irregulars) were newly born, and use a morphology-dependent local luminosity function to constrain the number of galaxies. The predictions of the model are compared with the observation of galaxy number counts and redshift distributions for the $B$, $I$ and $K$ colour bands. It is shown that young irregular galaxies cause the steep slope of the $B$-band counts. The fraction of irregular galaxies increases with decreasing brightness: at $B=24$ mag, they contribute as much as spiral galaxies. Thus, ``the faint blue galaxy problem'' is solved by invoking young galaxies. This interpretation is corroborated by a comparison of our prediction with the morphologically-classified galaxy counts in the $I$ band.

We investigate the gas mass function of clusters of galaxies to measure the density fluctuation spectrum on cluster scales. The baryon abundance confined in rich clusters is computed from the gas mass function and compared with the mean baryon density in the universe which is predicted by the Big Bang Nucleosynthesis. This baryon fraction and the slope of the gas mass function put constraints on $\sigma_8$, the rms linear fluctuation on scales of $8h^{-1}\Mpc$, and the slope of the fluctuation spectrum, where $h$ is the Hubble constant in units of 100 $\kms \oMpc$. We find $\sigma_8 = 0.80 \pm 0.15$ and $n \sim -1.5$ for $0.5 \le h \le 0.8$, where we assume that the density spectrum is approximated by a power law on cluster scales: $\sigma(r) \propto r^{-{3+n\over{2}}}$. Our value of $\sigma_8$ is independent of the density parameter, $\Omega_0$, and thus we can estimate $\Omega_0$ by combining $\sigma_8$ obtained in this study with those from $\Omega_0$-dependent analyses to date. We find that $\sigma_8(\Omega_0)$ derived from the cluster abundance such as the temperature function gives $\Omega_0 \sim 0.5$ while $\sigma_8(\Omega_0)$ measured from the peculiar velocity field of galaxies gives $\Omega_0 \sim 0.2-1$, depending on the technique used to analyze peculiar velocity data. Constraints are also derived for open, spatially flat, and tilted Cold Dark Matter models and for Cold + Hot Dark Matter models.