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.

We report the rest-frame ultraviolet luminosity function of $g$-dropout galaxies in 177 protocluster candidates (PC UVLF) at $z\sim4$ selected in the Hyper Suprime-Cam Subaru Strategic Program. Comparing with the UVLF of field galaxies at the same redshift, we find that the PC UVLF shows a significant excess towards the bright-end. This excess can not be explained by the contribution of only active galactic nuclei, and we also find that this is more significant in higher dense regions. Assuming that all protocluster members are located on the star formation main sequence, the PC UVLF can be converted into a stellar mass function. Consequently, our protocluster members are inferred to have a 2.8 times more massive characteristic stellar mass than that of the field Lyman break galaxies at the same redshift. This study, for the first time, clearly shows that the enhancement in star formation or stellar mass in overdense regions can generally be seen as early as at $z\sim4$. We also estimate the star formation rate density (SFRD) in protocluster regions as $\simeq 6-20\%$ of the cosmic SFRD, based on the measured PC UVLF after correcting for the selection incompleteness in our protocluster sample. This high value suggests that protoclusters make a non-negligible contribution to the cosmic SFRD at $z\sim4$, as previously suggested by simulations. Our results suggest that protoclusters are essential components for the galaxy evolution at $z\sim4$.

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.

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.

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

We present evidence that a region of high effective Ly$\alpha$ optical depth at $z\sim5.7$ is associated with an underdense region at the tail end of cosmic reionization. We carried out a survey of Lyman-break Galaxies (LBGs) using Subaru Hyper Suprime-Cam in the field of the $z=5.98$ quasar J0148+0600, whose spectrum presents an unusually long ($\sim160 \mathrm{cMpc}$) and opaque ($\tau \gtrsim7$) Ly$\alpha$ trough at $5.5\le z\le 5.9$. LBG candidates were selected to lie within the redshift range of the trough, and the projected number densities were measured within 90~cMpc of the quasar sightline. The region within $8'$ (or $\approx 19~\mathrm{cMpc}$) of the quasar position is the most underdense of the whole field. The significance of the presence of the void is estimated to be $99\%$. This is consistent with the significant deficit of Ly$\alpha$ emitters (LAEs) at $z=5.72$ reported by Becker et al. and suggests that the paucity of LAEs is not purely due to the removal of the Ly$\alpha$ emission by the high opacity but reflects a real coherent underdensity of galaxies across the entire redshift range of the trough. These observations are consistent with scenarios in which large optical depth fluctuations arise due to fluctuations in the galaxy-dominant UV background or due to residual neutral islands that are expected from reionization that is completed at redshifts as low as $z\lesssim5.5$.

We present the clustering properties of low-$z$ $(z\leq1.4)$ galaxies selected by the Hyper Suprime-Cam Subaru Strategic Program Wide layer over $145$ deg$^{2}$. The wide-field and multi-wavelength observation yields $5,064,770$ galaxies at $0.3\leq z\leq1.4$ with photometric redshifts and physical properties. This enables the accurate measurement of angular correlation functions and subsequent halo occupation distribution (HOD) analysis allows the connection between baryonic properties and dark halo properties. The fraction of less-massive satellite galaxies at $z\lesssim1$ is found to be almost constant at $\sim20\%$, but it gradually decreases beyond $M_{\star} \sim 10^{10.4}h^{-2}M_{\odot}$. However, the abundance of satellite galaxies at $z>1$ is quite small even for less-massive galaxies due to the rarity of massive centrals at high-$z$. This decreasing trend is connected to the small satellite fraction of Lyman break galaxies at $z>3$. The stellar-to-halo mass ratios at $0.3\leq z\leq1.4$ are almost consistent with the predictions obtained using the latest empirical model; however, we identify small excesses from the theoretical model at the massive end. The pivot halo mass is found to be unchanged at $10^{11.9-12.1}h^{-1}M_{\odot}$ at $0.3\leq z\leq1.4$, and we systematically show that $10^{12}h^{-1}M_{\odot}$ is a universal pivot halo mass up to $z\sim5$ that is derived using only the clustering/HOD analyses. Nevertheless, halo masses with peaked instantaneous baryon conversion efficiencies are much smaller than the pivot halo mass regardless of a redshift, and the most efficient stellar-mass assembly is thought to be in progress in $10^{11.0-11.5}h^{-1}M_{\odot}$ dark haloes.

We have carried out follow-up spectroscopy on three overdense regions of $g$- and $r$-dropout galaxies in the Canada-France-Hawaii Telescope Legacy Survey Deep Fields, finding two new protoclusters at $z=4.898$, 3.721 and a possible protocluster at $z=3.834$. The $z=3.721$ protocluster overlaps with a previously identified protocluster at $z=3.675$. The redshift separation between these two protoclusters is $\Delta z=0.05$, which is slightly larger than the size of typical protoclusters. Therefore, if they are not the progenitors of a $>10^{15}\,\mathrm{M_\odot}$ halo, they would grow into closely-located independent halos like a supercluster. The other protocluster at $z=4.898$ is also surrounded by smaller galaxy groups. These systems including protoclusters and neighboring groups are regarded as the early phase of superclusters. We quantify the spatial distribution of member galaxies of the protoclusters at $z=3.675$ and 3.721 by fitting triaxial ellipsoids, finding a tentative difference: one has a pancake-like shape while the other is filamentary. This could indicate that these two protoclusters are in different stages of formation. We investigate the relation between redshift and the velocity dispersion of protoclusters, including other protoclusters from the literature, in order to compare their dynamical states. Although there is no significant systematic trend in the velocity dispersions of protoclusters with redshift, the distribution is skewed to higher velocity dispersion over the redshift range of $z=2\mathrm{-}6$. This could be interpreted as two phases of cluster formation, one dominated by the steady accretion of galaxies, and the other by the merging between group-size halos, perhaps depending on the surrounding large-scale environments.

We perform a stacking analysis of {\it Planck}, {\it AKARI}, Infrared Astronomical Satellite ($IRAS$), Wide-field Infrared Survey Eplorer ($WISE$), and {\it Herschel} images of the largest number of (candidate) protoclusters at $z\sim3.8$ selected from the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP). Stacking the images of the $179$ candidate protoclusters, the combined infrared (IR) emission of the protocluster galaxies in the observed $12-850~\mu$m wavelength range is successfully detected with $>5\sigma$ significance (at $Planck$). This is the first time that the average IR spectral energy distribution (SED) of a protocluster has been constrained at $z\sim4$. The observed IR SEDs of the protoclusters exhibit significant excess emission in the mid-IR compared to that expected from typical star-forming galaxies (SFGs). They are reproduced well using SED models of intense starburst galaxies with warm/hot dust heated by young stars, or by a population of active galactic nuclei (AGN)/SFG composites. For the pure star-forming model, a total IR (from 8 to 1000 $\mu$m) luminosity of $19.3_{-4.2}^{+0.6}\times10^{13}~L_{\odot}$ and a star formation rate (SFR) of $16.3_{-7.8}^{+1.0}\times10^3~M_{\odot}$ yr$^{-1}$ are found whereas for the AGN/SFG composite model, $5.1_{-2.5}^{+2.5}\times10^{13}~L_{\odot}$ and $2.1^{+6.3}_{-1.7}\times10^3~M_{\odot}$ yr$^{-1}$ are found. Uncertainty remaining in the total SFRs; however, the IR luminosities of the most massive protoclusters are likely to continue increasing up to $z\sim4$. Meanwhile, no significant IR flux excess is observed around optically selected QSOs at similar redshifts, which confirms previous results. Our results suggest that the $z\sim4$ protoclusters trace dense, intensely star-forming environments that may also host obscured AGNs missed by the selection in the optical.

We present VLT/SINFONI observations of 35 quasars at 2.1 < z < 3.2, the majority of which were selected from the Clusters Around Radio-Loud AGN (CARLA) survey. CARLA quasars have large CIV-based black hole (BH) masses (M(BH) > 10^9 Msun) and powerful radio emission (P(500MHz) > 27.5 W/Hz). We estimate Ha-based M(BH), finding a scatter of 0.35 dex compared to CIV. We evaluate several recipes for correcting CIV-based masses, which reduce the scatter to 0.24 dex. The radio power of the radio-loud quasars is at most weakly correlated with the interconnected quantities Ha-width, L(5100A) and M(BH), suggesting that it is governed by different physical processes. However, we do find a strong inverse correlation between CIV blueshift and radio power linked to higher Eddington ratios and L(5100A). Under standard assumptions, the BH growth time is longer than the cosmic age for many CARLA quasars, suggesting that they must have experienced more efficient growth in the past. If these BHs were growing from seeds since the epoch of reionization, it is possible that they grew at the Eddington limit like the quasars at z ~ 6-7, and then continued to grow at the reduced rates observed until z ~ 2-3. Finally, we study the relation between M(BH) and environment, finding a weak positive correlation between M(BH) and galaxy density measured by CARLA.

We present ALMA [CII] line and far-infrared (FIR) continuum observations of three $z > 6$ low-luminosity quasars ($M_{\rm 1450} > -25$) discovered by our Subaru Hyper Suprime-Cam (HSC) survey. The [CII] line was detected in all three targets with luminosities of $(2.4 - 9.5) \times 10^8~L_\odot$, about one order of magnitude smaller than optically luminous ($M_{\rm 1450} \lesssim -25$) quasars. The FIR continuum luminosities range from $< 9 \times 10^{10}~L_\odot$ (3$\sigma$ limit) to $\sim 2 \times 10^{12}~L_\odot$, indicating a wide range in star formation rates in these galaxies. Most of the HSC quasars studied thus far show [CII]/FIR luminosity ratios similar to local star-forming galaxies. Using the [CII]-based dynamical mass ($M_{\rm dyn}$) as a surrogate for bulge stellar mass ($M_{\rm bulge}$), we find that a significant fraction of low-luminosity quasars are located on or even below the local $M_{\rm BH} - M_{\rm bulge}$ relation, particularly at the massive end of the galaxy mass distribution. In contrast, previous studies of optically luminous quasars have found that black holes are overmassive relative to the local relation. Given the low luminosities of our targets, we are exploring the nature of the early co-evolution of supermassive black holes and their hosts in a less biased way. Almost all of the quasars presented in this work are growing their black hole mass at much higher pace at $z \sim 6$ than the parallel growth model, in which supermassive black holes and their hosts grow simultaneously to match the local $M_{\rm BH} - M_{\rm bulge}$ relation at all redshifts. As the low-luminosity quasars appear to realize the local co-evolutionary relation even at $z \sim 6$, they should have experienced vigorous starbursts prior to the currently observed quasar phase to catch up with the relation.

We report the discovery of 28 quasars and 7 luminous galaxies at 5.7 $\le$ z $\le$ 7.0. This is the tenth in a series of papers from the Subaru High-z Exploration of Low-Luminosity Quasars (SHELLQs) project, which exploits the deep multi-band imaging data produced by the Hyper Suprime-Cam (HSC) Subaru Strategic Program survey. The total number of spectroscopically identified objects in SHELLQs has now grown to 93 high-z quasars, 31 high-z luminous galaxies, 16 [O III] emitters at z ~ 0.8, and 65 Galactic cool dwarfs (low-mass stars and brown dwarfs). These objects were found over 900 deg2, surveyed by HSC between 2014 March and 2018 January. The full quasar sample includes 18 objects with very strong and narrow Ly alpha emission, whose stacked spectrum is clearly different from that of other quasars or galaxies. While the stacked spectrum shows N V 1240 emission and resembles that of lower-z narrow-line quasars, the small Ly alpha width may suggest a significant contribution from the host galaxies. Thus these objects may be composites of quasars and star-forming galaxies.

We present deep near-infrared spectroscopy of six quasars at 6.1<z<6.7 with VLT/X-Shooter and Gemini-N/GNIRS. Our objects, originally discovered through a wide-field optical survey with the Hyper Suprime-Cam (HSC) Subaru Strategic Program (HSC-SSP), have the lowest luminosities (-25.5< M1450<-23.1 mag) of the z>5.8 quasars with measured black hole masses. From single-epoch mass measurements based on MgII2798, we find a wide range in black hole masses, from M_BH=10^7.6 to 10^9.3 Msun. The Eddington ratios L_bol/L_Edd range from 0.16 to 1.1, but the majority of the HSC quasars are powered by M_BH=10^9 Msun supermassive black holes (SMBHs) accreting at sub-Eddington rates. The Eddington ratio distribution of the HSC quasars is inclined to lower accretion rates than those of Willott et al. (2010a), who measured the black hole masses for similarly faint z=6 quasars. This suggests that the global Eddington ratio distribution is wider than has previously been thought. The presence of M_BH=10^9 Msun SMBHs at z=6 cannot be explained with constant sub-Eddington accretion from stellar remnant seed black holes. Therefore, we may be witnessing the first buildup of the most massive black holes in the first billion years of the universe, the accretion activity of which is transforming from active growth to a quiescent phase. Measurements of a larger complete sample of z>6 low-luminosity quasars, as well as deeper observations with future facilities will enable us to better understand the early SMBH growth in the reionization epoch.

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

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

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

We present new ALMA observations and physical properties of a Lyman Break Galaxy at z=7.15. Our target, B14-65666, has a bright ultra-violet (UV) absolute magnitude, $M_{\rm UV}\approx-22.4$, and has been spectroscopically identified in Ly$\alpha$ with a small rest-frame equivalent width of $\approx4$ \AA. Previous HST image has shown that the target is comprised of two spatially separated clumps in the rest-frame UV. With ALMA, we have newly detected spatially resolved [OIII] 88 $\mu$m, [CII] 158 $\mu$m, and their underlying dust continuum emission. In the whole system of B14-65666, the [OIII] and [CII] lines have consistent redshifts of $7.1520\pm0.0003$, and the [OIII] luminosity, $(34.4\pm4.1)\times10^{8}L_{\rm \odot}$, is about three times higher than the [CII] luminosity, $(11.0\pm1.4)\times10^{8}L_{\rm \odot}$. With our two continuum flux densities, the dust temperature is constrained to be $T_{\rm d}\approx50-60$ K under the assumption of the dust emissivity index of $\beta_{\rm d}=2.0-1.5$, leading to a large total infrared luminosity of $L_{\rm TIR}\approx1\times10^{12}L_{\rm \odot}$. Owing to our high spatial resolution data, we show that the [OIII] and [CII] emission can be spatially decomposed into two clumps associated with the two rest-frame UV clumps whose spectra are kinematically separated by $\approx200$ km s$^{-1}$. We also find these two clumps have comparable UV, infrared, [OIII], and [CII] luminosities. Based on these results, we argue that B14-65666 is a starburst galaxy induced by a major-merger. The merger interpretation is also supported by the large specific star-formation rate (defined as the star-formation rate per unit stellar mass), sSFR$=260^{+119}_{-57}$ Gyr$^{-1}$, inferred from our SED fitting. Probably, a strong UV radiation field caused by intense star formation contributes to its high dust temperature and the [OIII]-to-[CII] luminosity ratio.

We present the results of a survey of the brightest UV-selected galaxies in protoclusters. These proto-brightest cluster galaxy (proto-BCG) candidates are drawn from 179 overdense regions of $g$-dropout galaxies at $z\sim4$ from the Hyper Suprime-Cam Subaru Strategic Program identified previously as good protocluster candidates. This study is the first to extend the systematic study of the progenitors of BCGs from $z\sim2$ to $z\sim4$. We carefully remove possible contaminants from foreground galaxies and, for each structure, we select the brightest galaxy that is at least 1 mag brighter than the fifth brightest galaxy. We select 63 proto-BCG candidates and compare their properties with those of galaxies in the field and those of other galaxies in overdense structures. The proto-BCG candidates and their surrounding galaxies have different rest-UV color $(i - z)$ distributions to field galaxies and other galaxies in protoclusters that do not host proto-BCGs. In addition, galaxies surrounding proto-BCGs are brighter than those in protoclusters without proto-BCGs. The image stacking analysis reveals that the average effective radius of proto-BCGs is $\sim28\%$ larger than that of field galaxies. The $i-z$ color differences suggest that proto-BCGs and their surrounding galaxies are dustier than other galaxies at $z\sim4$. These results suggest that specific environmental effects or assembly biasses have already emerged in some protoclusters as early as $z \sim 4$, and we suggest that proto-BCGs have different star formation histories than other galaxies in the same epoch.

The large-scale overdensities of galaxies at z~2-7 known as protoclusters are believed to be the sites of cluster formation, and deep, wide survey projects such as the Large Synoptic Survey Telescope (LSST) and the Wide Field Infrared Survey Telescope (WFIRST) will deliver significant numbers of these interesting structures. Spectroscopic confirmation and interpretation of these targets, however, is still challenging, and will require wide-field multi-plexed spectroscopy on >20 m-class telescopes in the optical and near-infrared. In the coming decade, detailed studies of protoclusters will enable us, for the first time, to systematically connect these cluster progenitors in the early universe to their virialized counterparts at lower redshifts. This will allow us to address observationally the formation of brightest cluster galaxies and other cluster galaxy populations, the buildup of the intra-cluster light, the chemical enrichment history of the intra-cluster medium, and the formation and triggering of supermassive black holes in dense environments, all of which are currently almost exclusively approached either through the fossil record in clusters or through numerical simulations. Furthermore, at the highest redshifts (z~5-10), these large extended overdensities of star-forming galaxies are believed to have played an important role in the reionization of the universe, which needs to be tested by upcoming experiments. Theory and recent simulations also suggest important links between these overdensities and the formation of supermassive black holes, but observational evidence is still lacking. In this white paper we review our current understanding of this important phase of galaxy cluster history that will be explored by the next generation of large aperture ground-based telescopes GMT and TMT.

In the SSA22 field which exhibits a large-scale proto-cluster at $z=3.1$, we carried out a spectroscopic survey for Lyman Break Galaxies (LBGs) with the VLT/VIMOS and identified 78 confident LBGs at $z=2.5$--4. We stacked their spectra in the observer's frame by using a sophisticated method. Analyzing the composite spectrum, we have revealed that the large-scale proto-cluster at $z=3.1$ has a strong HI absorption dip of rest-frame equivalent width of $-1.7$ A. Another strong absorption dip found at $z=3.28$ is associated with a modestly high-density LBG peak, similar to that at $z=3.1$. We have also detected an HI transparency peak at $z=2.98$ in the composite spectrum, coincident with a void in the LBG distribution. In this paper, we also investigated the relation between LBGs, HI gas and AGNs at $z=3$--4 in the SSA22 field. Two AGNs at $z=3.353$ and 3.801 are, respectively, associated with the LBG concentration of an overdensity factor $\delta_{\rm LBG}\simeq2$ in the present statistics. Another structure at $z=3.453$ is remarkable: 20 comoving Mpc scale dense HI gas which is not associated with any apparent LBG overdensity but involving a pair of AGNs. Such structure may be a new type of the AGN-matter correlation. If the inhomogeneous structures over a comoving Gpc scale found in this paper are confirmed with sufficient statistics in the future, the SSA22 field will become a key region to test the standard cold dark matter structure formation scenario.

We report the discovery of a quasar at z = 7.07, which was selected from the deep multi-band imaging data collected by the Hyper Suprime-Cam (HSC) Subaru Strategic Program survey. This quasar, HSC J124353.93+010038.5, has an order of magnitude lower luminosity than do the other known quasars at z > 7. The rest-frame ultraviolet absolute magnitude is M1450 = -24.13 +/- 0.08 mag and the bolometric luminosity is Lbol = (1.4 +/- 0.1) x 10^{46} erg/s. Its spectrum in the optical to near-infrared shows strong emission lines, and shows evidence for a fast gas outflow, as the C IV line is blueshifted and there is indication of broad absorption lines. The Mg II-based black hole mass is Mbh = (3.3 +/- 2.0) x 10^8 Msun, thus indicating a moderate mass accretion rate with an Eddington ratio 0.34 +/- 0.20. It is the first z > 7 quasar with sub-Eddington accretion, besides being the third most distant quasar, known to date. The luminosity and black hole mass are comparable to, or even lower than, those measured for the majority of low-z quasars discovered by the Sloan Digital Sky Survey, and thus this quasar likely represents a z > 7 counterpart to quasars commonly observed in the low-z universe.

Enormous Ly$\alpha$ Nebulae (ELANe), unique tracers of galaxy density peaks, are predicted to lie at the nodes and intersections of cosmic filamentary structures. Previous successful searches for ELANe have focused on wide-field narrowband surveys, or have targeted known sources such as ultraluminous quasi-stellar-objects (QSOs) or radio galaxies. Utilizing groups of coherently strong Ly$\alpha$ absorptions (CoSLAs), we have developed a new method to identify high-redshift galaxy overdensities and have identified an extremely massive overdensity, BOSS1441, at $z=2-3$ (Cai et al. 2016a). In its density peak, we discover an ELAN that is associated with a relatively faint continuum. To date, this object has the highest diffuse Ly$\alpha$ nebular luminosity of $L_{\rm{nebula}}=5.1\pm0.1\times10^{44}$ erg s$^{-1}$. Above the 2$\sigma$ surface brightness limit of SB$_{\rm{Ly\alpha}}= 4.8\times10^{-18}$ erg s$^{-1}$ cm$^{-2}$ arcsec$^{-2}$, this nebula has an end-to-end spatial extent of 442 kpc. This radio-quiet source also has extended \civ\ $\lambda1549$ and \heii\ $\lambda1640$ emission on $\gtrsim30$ kpc scales. Note that the Ly$\alpha$, \heii\ and \civ\ emission all have double-peaked line profiles. Each velocity component has a full-width-half-maximum (FWHM) of $\approx700 - 1000$ km s$^{-1}$. We argue that this Ly$\alpha$ nebula could be powered by shocks due to an AGN-driven outflow or/and photoionization by a strongly obscured source.

We have carried out deep and wide field imaging observations with narrow bands targeting 11 quasar fields to systematically study the possible photoevaporation effect of quasar radiation on surrounding low mass galaxies at $z\sim2-3$. We focused on Lyman alpha emitters (LAEs) at the same redshifts as the quasars that lie within the quasar proximity zones, where the ultra-violet radiation from the quasars is higher than the average background at that epoch. We found that LAEs with high rest-frame equivalent width of Ly$\alpha$ emission ($EW_0$) of $\gtrsim 150$\AA$~$ with low stellar mass ($\lesssim 10^8 M_{\odot}$), are predominantly scarce in the quasar proximity zones, suggesting that quasar photoevaporation effects may be taking place. The halo mass of LAEs with $EW_0>150$\AA$~$ is estimated to be $3.6^{+12.7}_{-2.3}\times10^9 M_{\odot}$ either from the Spectral Energy Distribution (SED) fitting or the main sequence. Based on a hydrodynamical simulation, the predicted delay in star formation under a local UV background intensity with $J (\nu_L)\gtrsim10^{-21}$ erg s$^{-1}$ cm$^{-2}$ Hz$^{-1}$ sr$^{-1}$ for galaxies having less than this halo mass is about $>20$ Myr, which is longer than the expected age of LAEs with $EW_0>150$\AA. On the other hand, the photoevaporation seems to be less effective around very luminous quasars, which is consistent with the idea that these quasars are still in an early stage of quasar activity.

We present new measurements of the quasar luminosity function (LF) at $z \sim 6$, over an unprecedentedly wide range of the rest-frame ultraviolet luminosity $M_{1450}$ from $-30$ to $-22$ mag. This is the fifth in a series of publications from the Subaru High-$z$ Exploration of Low-Luminosity Quasars (SHELLQs) project, which exploits the deep multi-band imaging data produced by the Hyper Suprime-Cam (HSC) Subaru Strategic Program survey. The LF was calculated with a complete sample of 110 quasars at $5.7 \le z \le 6.5$, which includes 48 SHELLQs quasars discovered over 650 deg$^2$, and 63 brighter quasars discovered by the Sloan Digital Sky Survey and the Canada-France-Hawaii Quasar Survey (including one overlapping object). This is the largest sample of $z \sim 6$ quasars with a well-defined selection function constructed to date, and has allowed us to detect significant flattening of the LF at its faint end. A double power-law function fit to the sample yields a faint-end slope $\alpha = -1.23^{+0.44}_{-0.34}$, a bright-end slope $\beta = -2.73^{+0.23}_{-0.31}$, a break magnitude $M_{1450}^* = -24.90^{+0.75}_{-0.90}$, and a characteristic space density $\Phi^* = 10.9^{+10.0}_{-6.8}$ Gpc$^{-3}$ mag$^{-1}$. Integrating this best-fit model over the range $-18 < M_{1450} < -30$ mag, quasars emit ionizing photons at the rate of $\dot{n}_{\rm ion} = 10^{48.8 \pm 0.1}$ s$^{-1}$ Mpc$^{-3}$ at $z = 6.0$. This is less than 10 % of the critical rate necessary to keep the intergalactic medium ionized, which indicates that quasars are not a major contributor to cosmic reionization.

Sample return from the near-Earth asteroid known as 25143 Itokawa was conducted as part of the Hayabusa mission, with a large number of scientific findings being derived from the returned samples. Following the Hayabusa mission, Hayabusa2 was planned, targeting sample return from a primitive asteroid. The primary target body of Hayabusa2 was asteroid 162173 Ryugu; however, it was also necessary to gather physical information for backup target selection. Therefore, we examined five asteroids spectroscopically, 43 asteroids spectrophotometrically, and 41 asteroids through periodic analysis. Hence, the physical properties of 74 near-Earth asteroids were obtained, which helped the Hayabusa2 backup target search and, also, furthered understanding of the physical properties of individual asteroids and their origins.

We search for velocity changes (i.e., acceleration/deceleration) of narrow absorption lines (NALs) that are intrinsic to the quasars, using spectra of 6 bright quasars that have been observed more than once with 8-10m class telescopes. While variations in line strength and profile are frequently reported (especially in broader absorption lines), definitive evidence for velocity shifts has not been found with only a few exceptions. Direct velocity shift measurements are valuable constraints on the acceleration mechanisms. In this study, we determine velocity shifts by comparing the absorption profiles of NALs at two epochs separated by more than 10 years in the observed frame, using the cross-correlation function method and we estimate the uncertainties using Monte Carlo simulations. We do not detect any significant shifts but we obtain 3$\sigma$ upper limits on the acceleration of intrinsic NALs (compared to intervening NALs in same quasars) of $\sim$0.7 km s$^{-1}$ yr$^{-1}$ ($\sim$0.002 cm s$^{-2}$). We discuss possible scenarios for non-detection of NAL acceleration/deceleration and examine resulting constraints on the physical conditions in accretion disk wind.

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

We report discovery of 41 new high-z quasars and luminous galaxies, which were spectroscopically identified at 5.7 < z < 6.9. This is the fourth in a series of papers from the Subaru High-z Exploration of Low-Luminosity Quasars (SHELLQs) project, based on the deep multi-band imaging data collected by the Hyper Suprime-Cam (HSC) Subaru Strategic Program survey. We selected the photometric candidates by a Bayesian probabilistic algorithm, and then carried out follow-up spectroscopy with the Gran Telescopio Canarias and the Subaru Telescope. Combined with the sample presented in the previous papers, we have now spectroscopically identified 137 extremely-red HSC sources over about 650 deg2, which include 64 high-z quasars, 24 high-z luminous galaxies, 6 [O III] emitters at z ~ 0.8, and 43 Galactic cool dwarfs (low-mass stars and brown dwarfs). The new quasars span the luminosity range from M1450 ~ -26 to -22 mag, and continue to populate a few magnitude lower luminosities than have been probed by previous wide-field surveys. In a companion paper, we derive the quasar luminosity function at z ~ 6 over an unprecedentedly wide range of M1450 ~ -28 to -21 mag, exploiting the SHELLQs and other survey outcomes.

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.

Quasars (QSOs) hosting supermassive black holes are believed to reside in massive halos harboring galaxy overdensities. However, many observations revealed average or low galaxy densities around $z\gtrsim6$ QSOs. This could be partly because they measured galaxy densities in only tens of arcmin$^2$ around QSOs and might have overlooked potential larger scale galaxy overdensities. Some previous studies also observed only Lyman break galaxies (LBGs, massive older galaxies) and missed low mass young galaxies like Ly$\alpha$ emitters (LAEs) around QSOs. Here we present observations of LAE and LBG candidates in $\sim700$ arcmin$^2$ around a $z=6.61$ luminous QSO using Subaru Telescope Suprime-Cam with narrow/broadbands. We compare their sky distributions, number densities and angular correlation functions with those of LAEs/LBGs detected in the same manner and comparable data quality in our control blank field. In the QSO field, LAEs and LBGs are clustering in 4-20 comoving Mpc angular scales, but LAEs show mostly underdensity over the field while LBGs are forming $30\times60$ comoving Mpc$^2$ large scale structure containing 3-$7\sigma$ high density clumps. The highest density clump includes a bright (23.78 mag in the narrowband) extended ($\gtrsim 16$ kpc) Ly$\alpha$ blob candidate, indicative of a dense environment. The QSO could be part of the structure but is not located exactly at any of the high density peaks. Near the QSO, LAEs show underdensity while LBGs average to $4\sigma$ excess densities compared to the control field. If these environments reflect halo mass, the QSO may not be in the most massive halo, but still in a moderately massive one.

We present our ALMA Cycle 4 measurements of the [CII] emission line and the underlying far-infrared (FIR) continuum emission from four optically low-luminosity ($M_{\rm 1450} > -25$) quasars at $z \gtrsim 6$ discovered by the Subaru Hyper Suprime Cam (HSC) survey. The [CII] line and FIR continuum luminosities lie in the ranges $L_{\rm [CII]} = (3.8-10.2) \times 10^8~L_\odot$ and $L_{\rm FIR} = (1.2-2.0) \times 10^{11}~L_\odot$, which are at least one order of magnitude smaller than those of optically-luminous quasars at $z \gtrsim 6$. We estimate the star formation rates (SFR) of our targets as $\simeq 23-40~M_\odot ~{\rm yr}^{-1}$. Their line and continuum-emitting regions are marginally resolved, and found to be comparable in size to those of optically luminous quasars, indicating that their SFR or likely gas mass surface densities (key controlling parameter of mass accretion) are accordingly different. The $L_{\rm [CII]}/L_{\rm FIR}$ ratios of the hosts, $\simeq (2.2-8.7) \times 10^{-3}$, are fully consistent with local star-forming galaxies. Using the [CII] dynamics, we derived their dynamical masses within a radius of 1.5-2.5 kpc as $\simeq (1.4-8.2) \times 10^{10}~M_\odot$. By interpreting these masses as stellar ones, we suggest that these faint quasar hosts are on or even below the star-forming main sequence at $z \sim 6$, i.e., they appear to be transforming into quiescent galaxies. This is in contrast to the optically luminous quasars at those redshifts, which show starburst-like properties. Finally, we find that the ratios of black hole mass to host galaxy dynamical mass of the most of low-luminosity quasars including the HSC ones are consistent with the local value. The mass ratios of the HSC quasars can be reproduced by a semi-analytical model that assumes merger-induced black hole-host galaxy evolution.

We report fourteen and twenty-eight protocluster candidates at z=5.7 and 6.6 over 14 and 19 deg^2 areas, respectively, selected from 2,230 (259) Lya emitters (LAEs) photometrically (spectroscopically) identified with Subaru/Hyper Suprime-Cam (HSC) deep images (Keck, Subaru, and Magellan spectra and the literature data). Six out of the 42 protocluster candidates include 1-12 spectroscopically confirmed LAEs at redshifts up to z=6.574. By the comparisons with the cosmological Lya radiative transfer (RT) model reproducing LAEs with the reionization effects, we find that more than a half of these protocluster candidates are progenitors of the present-day clusters with a mass of > 10^14 M_sun. We then investigate the correlation between LAE overdensity delta and Lya rest-frame equivalent width EW_Lya^rest, because the cosmological Lya RT model suggests that a slope of EW_Lya^rest-delta relation is steepened towards the epoch of cosmic reionization (EoR), due to the existence of the ionized bubbles around galaxy overdensities easing the escape of Lya emission from the partly neutral intergalactic medium (IGM). The available HSC data suggest that the slope of the EW_Lya^rest-delta correlation does not evolve from the post-reionization epoch z=5.7 to the EoR z=6.6 beyond the moderately large statistical errors. There is a possibility that we would detect the evolution of the EW_Lya^rest - delta relation from z=5.7 to 7.3 by the upcoming HSC observations providing large samples of LAEs at z=6.6-7.3.

We study the internal structure of the Circum-Galactic Medium (CGM), using 29 spectra of 13 gravitationally lensed quasars with image separation angles of a few arcseconds, which correspond to 100 pc to 10 kpc in physical distances. After separating metal absorption lines detected in the spectra into high-ions with ionization parameter (IP) $>$ 40 eV and low-ions with IP $<$ 20 eV, we find that i) the fraction of absorption lines that are detected in only one of the lensed images is larger for low-ions ($\sim$16%) than high-ions ($\sim$2%), ii) the fractional difference of equivalent widths ($EW$s) between the lensed images is almost same (${\rm d}EW$ $\sim$ 0.2) for both groups although the low-ions have a slightly larger variation, and iii) weak low-ion absorbers tend to have larger ${\rm d}EW$ compared to weak high-ion absorbers. We construct simple models to reproduce these observed properties and investigate the distribution of physical quantities such as size and location of absorbers, using some free parameters. Our best models for absorbers with high-ions and low-ions suggest that i) an overall size of the CGM is at least $\sim$ 500 kpc, ii) a size of spherical clumpy cloud is $\sim$ 1 kpc or smaller, and iii) only high-ion absorbers can have diffusely distributed homogeneous component throughout the CGM. We infer that a high ionization absorber distributes almost homogeneously with a small-scale internal fluctuation, while a low ionization absorber consists of a large number of small-scale clouds in the diffusely distributed higher ionized region. This is the first result to investigate the internal small-scale structure of the CGM, based on the large number of gravitationally lensed quasar spectra.

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 report the study of far-IR sizes of submillimeter galaxies (SMGs) in relation to their dust-obscured star formation rate (SFR) and active galactic nuclei (AGN) presence, determined using mid-IR photometry. We determined the millimeter-wave ($\lambda_{\rm obs}=1100 \mu$m) sizes of 69 ALMA-identified SMGs, selected with $\geq10$$\sigma$ confidence on ALMA images ($F_{\rm 1100 \mu m}=1.7$--7.4 mJy). We found that all the SMGs are located above an avoidance region in the millimeter size-flux plane, as expected by the Eddington limit for star formation. In order to understand what drives the different millimeter-wave sizes in SMGs, we investigated the relation between millimeter-wave size and AGN fraction for 25 of our SMGs at $z=1$--3. We found that the SMGs for which the mid-IR emission is dominated by star formation or AGN have extended millimeter-sizes, with respective median $R_{\rm c,e} = 1.6^{+0.34}_{-0.21}$ and 1.5$^{+0.93}_{-0.24}$ kpc. Instead, the SMGs for which the mid-IR emission corresponds to star-forming/AGN composites have more compact millimeter-wave sizes, with median $R_{\rm c,e}=1.0^{+0.20}_{-0.20}$ kpc. The relation between millimeter-wave size and AGN fraction suggests that this size may be related to the evolutionary stage of the SMG. The very compact sizes for composite star-forming/AGN systems could be explained by supermassive black holes growing rapidly during the SMG coalescing, star-formation phase.

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 the cross-correlation between 151 luminous quasars ($M_{ \mathrm{UV}} < -26$) and 179 protocluster candidates at $z \sim 3.8$, extracted from the Wide imaging survey ($ \sim 121~ $deg$^2$) performed with a part of the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP). We find that only two out of 151 quasars reside in regions that are more overdense compared to the average field at $ > 4 \sigma $. The distributions of the distance between quasars and the nearest protoclusters and the significance of the overdensity at the position of quasars are statistically identical to those found for $g$-dropout galaxies, suggesting that quasars tend to reside in almost the same environment as star-forming galaxies at this redshift. Using stacking analysis, we find that the average density of $g$-dropout galaxies around quasars is slightly higher than that around $g$-dropout galaxies on $1.0 - 2.5$ pMpc scales, while at $ < 0.5$ pMpc that around quasars tends to be lower. We also find that quasars with higher UV-luminosity or with more massive black holes tend to avoid the most overdense regions, and that the quasar near zone sizes are anti-correlated with overdensity. These findings are consistent with a scenario in which the luminous quasar at $z \sim4 $ resides in structures that are less massive than those expected for the progenitors of today's rich clusters of galaxies, and possibly that luminous quasars may be suppressing star formation in their close vicinity.

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

We constrain the quasar contribution to cosmic reionization based on our deep optical survey of z~6 quasars down to z_R=24.15 using Subaru/Suprime-Cam in three UKIDSS-DXS fields covering 6.5 deg^2. In Kashikawa et al. (2015), we select 17 quasar candidates and report our initial discovery of two low-luminosity quasars (M_1450~ -23) from seven targets, one of which might be a Lyman alpha emitting galaxy. From an additional optical spectroscopy, none of the four candidates out of the remaining ten turn out to be genuine quasars. Moreover, the deeper optical photometry provided by the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) shows that, unlike the two already-known quasars, the i-z and z-y colors of the last six candidates are consistent with M- or L-type brown dwarfs. Therefore, the quasar luminosity function (QLF) in the previous paper is confirmed. Compiling QLF measurements from the literature over a wide magnitude range, including an extremely faint AGN candidate from Parsa et al. (2017}, to fit them with a double power-law, we find that the best-fit faint-end slope is alpha=-2.04^+0.33_-0.18 (-1.98^+0.48_-0.21) and characteristic magnitude is M_1450^*=-25.8^+1.1_-1.9 (-25.7^+1.0_-1.8) in the case of two (one) quasar detection. Our result suggests that, if the QLF is integrated down to M_1450=-18, quasars produce ~1-12% of the ionizing photons required to ionize the whole universe at z~6 with 2sigma confidence level, assuming that the escape fraction is f_esc=1 and the IGM clumpy factor is C=3. Even when the systematic uncertainties are taken into account, our result supports the scenario that quasars are the minor contributors of reionization.

We conduct a systematic search for galaxy protoclusters at $z\sim3.8$ based on the latest internal data release (S16A) of the Hyper SuprimeCam Subaru strategic program (HSC-SSP). In the Wide layer of the HSC-SSP, we investigate the large-scale projected sky distribution of $g$-dropout galaxies over an area of $121\,\mathrm{deg^2}$, and identify 216 large-scale overdense regions ($>4\sigma$ overdensity significance) that are good protocluster candidates. Of these, 37 are located within $8\,\mathrm{arcmin}$ ($3.4\,\mathrm{physicalMpc}$) from other protocluster candidates of higher overdensity, and are expected to merge into a single massive structure by $z=0$. Therefore, we find 179 unique protocluster candidates in our survey. A cosmological simulation that includes projection effects predicts that more than 76\% of these candidates will evolve into galaxy clusters with halo masses of at least $10^{14}\,M_{\odot}$ by $z=0$. The unprecedented size of our protocluster candidate catalog allowed us to perform, for the first time, an angular clustering analysis of the systematic sample of protocluster candidates. We find a correlation length of $35.0\,h^{-1}\,\mathrm{Mpc}$. The relation between correlation length and number density of $z\sim3.8$ protocluster candidates is consistent with the prediction of the $\Lambda$CDM model, and the correlation length is similar to that of rich clusters in the local universe. This result suggests that our protocluster candidates are tracing similar spatial structures as those expected of the progenitors of rich clusters and enhances the confidence that our method to identify protoclusters at high redshifts is robust. In the coming years, our protocluster search will be extended to the entire HSC-SSP Wide sky coverage of $\sim1400\,\mathrm{deg^2}$ to probe cluster formation over a wide redshift range of $z\sim2\mathrm{-}6$.

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 investigate the galaxy overdensity around proto-cluster scale quasar pairs at high (z>3) and low (z~1) redshift based on the unprecedentedly wide and deep optical survey of the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP). Using the first-year survey data covering effectively ~121 deg^2 with the 5sigma depth of i~26.4 and the SDSS DR12Q catalog, we find two luminous pairs at z~3.3 and 3.6 which reside in >5sigma overdense regions of g-dropout galaxies at i<25. The projected separations of the two pairs are R_perp=1.75 and 1.04 proper Mpc, and their velocity offsets are Delta V=692 and 1448 km s^{-1}, respectively. This result is in clear contrast to the average z~4 quasar environments as discussed in Uchiyama et al. (2017) and implies that the quasar activities of the pair members are triggered via major mergers in proto-clusters, unlike the vast majority of isolated quasars in general fields that may turn on via non-merger events such as bar and disk instabilities. At z~1, we find 37 pairs with R_perp<2 pMpc and Delta V<2300 km s^{-1} in the current HSC-Wide coverage, including four from Hennawi et al. (2006). The distribution of the peak overdensity significance within two arcminutes around the pairs has a long tail toward high density (>4sigma) regions. Thanks to the large sample size, we find a statistical evidence that this excess is unique to the pair environments when compared to single quasar and randomly selected galaxy environments at the same redshift range. Moreover, there are nine small-scale (R_perp<1 pMpc) pairs, two of which are found to reside in cluster fields. Our results demonstrate that <2 pMpc-scale quasar pairs at both redshift range tend to occur in massive haloes, although perhaps not the most massive ones, and that they are useful to search for rare density peaks.

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 detect 20 $z=7.0$ Ly$\alpha$ emitter (LAE) candidates to $L({\rm Ly}\alpha) \geq 2 \times 10^{42}$ erg s$^{-1}$ or $0.3$ $L^*_{z=7}$ and in $6.1\times 10^5$ Mpc$^3$ volume in the Subaru Deep Field and the Subaru/XMM-Newton Deep Survey field by 82 and 37 hours of Subaru Suprime-Cam narrowband NB973 and reddest optical $y$-band imaging. We compare their Ly$\alpha$ and UV luminosity functions (LFs) and densities and Ly$\alpha$ equivalent widths (EWs) to those of $z=5.7$, 6.6 and 7.3 LAEs from previous Suprime-Cam surveys. The Ly$\alpha$ LF (density) rapidly declines by a factor of $\times$1.5 (1.9) in $L({\rm Ly}\alpha)$ at $z=5.7-6.6$ (160 Myr), $\times$1.5 (1.6) at $z=6.6-7.0$ (60 Myr) at the faint end and $\times$2.0 (3.8) at $z=7.0-7.3$ (40 Myr). Also, in addition to the systematic decrease in EW at $z=5.7-6.6$ previously found, 2/3 of the $z=7.0$ LAEs detected in the UV continuum exhibit lower EWs than the $z=6.6$ ones. Moreover, while the UV LF and density do not evolve at $z=5.7-6.6$, they modestly decline at $z=6.6-7.0$, implying galaxy evolution contributing to the decline of the Ly$\alpha$ LF. Comparison of the $z=7.0$ Ly$\alpha$ LF to the one predicted by an LAE evolution model further reveals that galaxy evolution alone cannot explain all the decline of Ly$\alpha$ LF. If we attribute the discrepancy to Ly$\alpha$ attenuation by neutral hydrogen, the intergalactic medium transmission of Ly$\alpha$ photons at $z=7.0$ would be $T_{{\rm Ly}\alpha}^{\rm IGM} \leq 0.6-0.7$. It is lower (higher) than $T_{{\rm Ly}\alpha}^{\rm IGM}$ at $z=6.6$ (7.3) derived by previous studies, suggesting rapid increase in neutral fraction at $z > 6$.

We present spectroscopic identification of 32 new quasars and luminous galaxies discovered at 5.7 < z < 6.8. This is the second in a series of papers presenting the results of the Subaru High-z Exploration of Low-Luminosity Quasars (SHELLQs) project, which exploits the deep multi-band imaging data produced by the Hyper Suprime-Cam (HSC) Subaru Strategic Program survey. The photometric candidates were selected by a Bayesian probabilistic algorithm, and then observed with spectrographs on the Gran Telescopio Canarias and the Subaru Telescope. Combined with the sample presented in the previous paper, we have now identified 64 HSC sources over about 430 deg2, which include 33 high-z quasars, 14 high-z luminous galaxies, 2 [O III] emitters at z ~ 0.8, and 15 Galactic brown dwarfs. The new quasars have considerably lower luminosity (M1450 ~ -25 to -22 mag) than most of the previously known high-z quasars. Several of these quasars have luminous (> 10^(43) erg/s) and narrow (< 500 km/s) Ly alpha lines, and also a possible mini broad absorption line system of N V 1240 in the composite spectrum, which clearly separate them from typical quasars. On the other hand, the high-z galaxies have extremely high luminosity (M1450 ~ -24 to -22 mag) compared to other galaxies found at similar redshift. With the discovery of these new classes of objects, we are opening up new parameter spaces in the high-z Universe. Further survey observations and follow-up studies of the identified objects, including the construction of the quasar luminosity function at z ~ 6, are ongoing.

We examine the clustering of quasars over a wide luminosity range, by utilizing 901 quasars at $\overline{z}_{\rm phot}\sim3.8$ with $-24.73<M_{\rm 1450}<-22.23$ photometrically selected from the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) S16A Wide2 date release and 342 more luminous quasars at $3.4<z_{\rm spec}<4.6$ having $-28.0<M_{\rm 1450}<-23.95$ from the Sloan Digital Sky Survey (SDSS) that fall in the HSC survey fields. We measure the bias factors of two quasar samples by evaluating the cross-correlation functions (CCFs) between the quasar samples and 25790 bright $z\sim4$ Lyman Break Galaxies (LBGs) in $M_{\rm 1450}<-21.25$ photometrically selected from the HSC dataset. Over an angular scale of \timeform{10.0"} to \timeform{1000.0"}, the bias factors are $5.93^{+1.34}_{-1.43}$ and $2.73^{+2.44}_{-2.55}$ for the low and high luminosity quasars, respectively, indicating no luminosity dependence of quasar clustering at $z\sim4$. It is noted that the bias factor of the luminous quasars estimated by the CCF is smaller than that estimated by the auto-correlation function (ACF) over a similar redshift range, especially on scales below \timeform{40.0"}. Moreover, the bias factor of the less-luminous quasars implies the minimal mass of their host dark matter halos (DMHs) is $0.3$-$2\times10^{12}h^{-1}M_{\odot}$, corresponding to a quasar duty cycle of $0.001$-$0.06$.

We present the luminosity function of z=4 quasars based on the Hyper Suprime-Cam Subaru Strategic Program Wide layer imaging data in the g, r, i, z, and y bands covering 339.8 deg^2. From stellar objects, 1666 z~4 quasar candidates are selected by the g-dropout selection down to i=24.0 mag. Their photometric redshifts cover the redshift range between 3.6 and 4.3 with an average of 3.9. In combination with the quasar sample from the Sloan Digital Sky Survey in the same redshift range, the quasar luminosity function covering the wide luminosity range of M1450=-22 to -29 mag is constructed. It is well described by a double power-law model with a knee at M1450=-25.36+-0.13 mag and a flat faint-end slope with a power-law index of -1.30+-0.05. The knee and faint-end slope show no clear evidence of redshift evolution from those at z~2. The flat slope implies that the UV luminosity density of the quasar population is dominated by the quasars around the knee, and does not support the steeper faint-end slope at higher redshifts reported at z>5. If we convert the M1450 luminosity function to the hard X-ray 2-10keV luminosity function using the relation between UV and X-ray luminosity of quasars and its scatter, the number density of UV-selected quasars matches well with that of the X-ray-selected AGNs above the knee of the luminosity function. Below the knee, the UV-selected quasars show a deficiency compared to the hard X-ray luminosity function. The deficiency can be explained by the lack of obscured AGNs among the UV-selected quasars.

We present the results of clustering analyses of Lyman break galaxies (LBGs) at $z\sim3$, $4$, and $5$ using the final data release of the Canada--France--Hawaii Telescope Legacy Survey (CFHTLS). Deep- and wide-field images of the CFHTLS Deep Survey enable us to obtain sufficiently accurate two-point angular correlation functions to apply a halo occupation distribution analysis. Mean halo masses, calculated as $\langle M_{h} \rangle=10^{11.7}-10^{12.8}h^{-1}M_{\odot}$, increase with stellar-mass limit of LBGs. The threshold halo mass to have a central galaxy follows the same increasing trend with the low-$z$ results, whereas the threshold halo mass to have a satellite galaxy shows higher values at $z=3-5$ than $z=0.5-1.5$ over the entire stellar mass range. Satellite fractions of dropout galaxies, even at less massive haloes, are found to drop sharply from $z=2$ down to less than $0.04$ at $z=3-5$. These results suggest that satellite galaxies form inefficiently within dark haloes at $z=3-5$ even for less massive satellites with $M_{\star}<10^{10}M_{\odot}$. We compute stellar-to-halo mass ratios (SHMRs) assuming a main sequence of galaxies, which is found to provide consistent SHMRs with those derived from a spectral energy distribution fitting method. The observed SHMRs are in good agreement with the model predictions based on the abundance-matching method within $1\sigma$ confidence intervals. We derive observationally, for the first time, $M_{{\rm h}}^{{\rm pivot}}$, which is the halo mass at a peak in the star-formation efficiency, at $3<z<5$, and it shows a little increasing trend with cosmic time at $z>3$. In addition, $M_{{\rm h}}^{{\rm pivot}}$ and its normalization are found to be almost unchanged during $0<z<5$. Our study shows an observational evidence that galaxy formation is ubiquitously most efficient near a halo mass of $M_{{\rm h}}\sim10^{12}M_{\odot}$ over cosmic time.

We surveyed a central 0.6 deg^2 region of Abell 1367 cluster for extended ionized gas clouds (EIGs) using the Subaru prime-focus camera (Suprime-Cam) with a narrow-band filter that covers Halpha. We discovered six new EIGs in addition to five known EIGs. We also found that the Halpha tail from the blue infalling group (BIG) is extended to about 330 kpc in projected distance, which is about twice longer than previously reported. Candidates of star-forming blobs in the tail are detected. The properties of the EIG parent galaxies in Abell 1367 basically resemble those in the Coma cluster. A noticeable difference is that the number of detached EIGs is significantly fewer in Abell 1367, while the fraction of blue member galaxies is higher. The results suggest a difference in the evolutionary stage of the clusters; Abell 1367 is at an earlier stage than the Coma cluster.

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

We have carried out a panoramic Ly$\alpha$ narrowband imaging with Suprime-Cam on Subaru towards the known protocluster USS1558--003 at $z=2.53$. Our previous narrowband imaging at near-infrared has identified multiple dense groups of H$\alpha$ emitters (HAEs) within the protocluster. We have now identified the large-scale structures across a $\sim$50 comoving Mpc scale traced by Ly$\alpha$ emitters (LAEs) in which the protocluster traced by the HAEs is embedded. On a smaller scale, however, there are remarkably few LAEs in the regions of HAE overdensities. Moreover, the stacking analyses of the images show that HAEs in higher-density regions show systematically lower escape fractions of Ly$\alpha$ photons than those of HAEs in lower-density regions. These phenomena may be driven by the extra depletion of Ly$\alpha$ emission lines along our line of sight by more intervening cold circumgalactic/intergalactic medium and/or dust existing in the dense core. We also caution that all the past high-$z$ protocluster surveys using LAEs as the tracers would have largely missed galaxies in the very dense cores of the protoclusters where we would expect to see any early environmental effects.

We conducted a deep narrow-band imaging survey with the Subaru Prime Focus Camera on the Subaru Telescope and constructed a sample of Ly$\alpha$ emitters (LAEs) at z=2.53 in the UDS-CANDELS field where a sample of H$\alpha$ emitters (HAEs) at the same redshift is already obtained from our previous narrow-band observation at NIR. The deep narrow-band and multi broadband data allow us to find LAEs of stellar masses and star-formation rates (SFRs) down to $\gtrsim$$10^8$ M$_\odot$ and $\gtrsim$0.2 M$_\odot$/yr, respectively. We show that the LAEs are located along the same mass-SFR sequence traced by normal star-forming galaxies such as HAEs, but towards a significantly lower mass regime. Likewise, LAEs seem to share the same mass--size relation with typical star-forming galaxies, except for the massive LAEs, which tend to show significantly compact sizes. We identify a vigorous mass growth in the central part of LAEs: the stellar mass density in the central region of LAEs increases as their total galaxy mass grows. On the other hand, we see no Ly$\alpha$ line in emission for most of the HAEs. Rather, we find that the Ly$\alpha$ feature is either absent or in absorption (Ly$\alpha$ absorbers; LAAs), and its absorption strength may increase with reddening of the UV continuum slope. We demonstrate that a deep Ly$\alpha$ narrow-band imaging like this study is able to search for not only LAEs but also LAAs in a certain redshift slice. This work suggests that LAEs trace normal star-forming galaxies in the low-mass regime, while they remain as a unique population because the majority of HAEs are not LAEs.

We discuss the host galaxy metallicity distribution of all long gamma-ray bursts (GRBs) whose redshifts are known to be $< 0.4$, including newly obtained spectroscopic datasets of the host galaxies of GRB 060614, 090417B, and 130427A. We compare the metallicity distribution of the low-redshift sample to the model predictions, and constrain the relation between metallicity and GRB occurrence. We take account of spatial variation of metallicities among star forming regions within a galaxy. We found that the models, in which only low-metallicity stars produce GRBs with a sharp cutoff of GRB production efficiency around 12+log(O/H) $\sim$ 8.3, can well reproduce the observed distribution, while the models with no metallicity dependence are not consistent with the observations. We also discuss possible sampling biases we may suffer by collecting long GRBs whose redshifts are known, presenting the photometric observations of the host galaxy of GRB 111225A at $z = 0.297$ whose redshift has been undetermined until $\sim$ 2.3 years after the burst.

Deep rest-frame optical spectroscopy is critical for characterizing and understanding the physical conditions and properties of the ionized gas in galaxies. Here, we present a new spectroscopic survey called "Metal Abundances across Cosmic Time" or $\mathcal{MACT}$, which will obtain rest-frame optical spectra for $\sim$3000 emission-line galaxies. This paper describes the optical spectroscopy that has been conducted with MMT/Hectospec and Keck/DEIMOS for $\approx$1900 $z=0.1-1$ emission-line galaxies selected from our narrowband and intermediate-band imaging in the Subaru Deep Field. In addition, we present a sample of 164 galaxies for which we have measured the weak [OIII]$\lambda$4363 line (66 with at least 3$\sigma$ detections and 98 with significant upper limits). This nebular emission line determines the gas-phase metallicity by measuring the electron temperature of the ionized gas. This paper presents the optical spectra, emission-line measurements, interstellar properties (e.g., metallicity, gas density), and stellar properties (e.g., star formation rates, stellar mass). Paper II of the $\mathcal{MACT}$ survey (Ly et al.) presents the first results on the stellar mass--gas metallicity relation at $z\lesssim1$ using the sample with [OIII]$\lambda$4363 measurements.

Cosmological simulations suggest a strong correlation between high optical-depth Ly$\alpha$ absorbers, which arise from the intergalactic medium (IGM), and 3-D mass overdensities on scales of $10-30$ $h^{-1}$ comoving Mpc. By examining the absorption spectra of $\sim$ 80,000 QSO sight-lines over a volume of 0.1 Gpc$^3$ in the Sloan Digital Sky Survey III (SDSS-III), we have identified an extreme overdensity, BOSS1441, which contains a rare group of strong Ly$\alpha$ absorbers at $z=2.32\pm 0.02$. This absorber group is associated with six QSOs at the same redshift on a 30 comoving Mpc scale. Using Mayall/MOSAIC narrowband and broadband imaging, we detect Ly$\alpha$ emitters (LAEs) down to $0.7\times L_{\rm{Ly\alpha}}^*$, and reveal a large-scale structure of Ly$\alpha$ emitters (LAEs) in this field. Our follow-up Large Binocular Telescope (LBT) observations have spectroscopically confirmed 19 galaxies in the density peak. We show that BOSS1441 has an LAE overdensity of $10.8\pm 2.6$ on a 15 comoving Mpc scale which could collapse to a massive cluster with $M\gtrsim10^{15}$ M$_\odot$ at $z\sim0$. This overdensity is among the most massive large-scale structures at $z\sim2$ discovered to date.

We report the discovery of 15 quasars and bright galaxies at 5.7 < z < 6.9. This is the initial result from the Subaru High-z Exploration of Low-Luminosity Quasars (SHELLQs) project, which exploits the exquisite multiband imaging data produced by the Subaru Hyper Suprime-Cam (HSC) Strategic Program survey. The candidate selection is performed by combining several photometric approaches including a Bayesian probabilistic algorithm to reject stars and dwarfs. The spectroscopic identification was carried out with the Gran Telescopio Canarias and the Subaru Telescope for the first 80 deg2 of the survey footprint. The success rate of our photometric selection is quite high, approaching 100 % at the brighter magnitudes (zAB < 23.5 mag). Our selection also recovered all the known high-z quasars on the HSC images. Among the 15 discovered objects, six are likely quasars, while the other six with interstellar absorption lines and in some cases narrow emission lines are likely bright Lyman-break galaxies. The remaining three objects have weak continua and very strong and narrow Ly alpha lines, which may be excited by ultraviolet light from both young stars and quasars. These results indicate that we are starting to see the steep rise of the luminosity function of z > 6 galaxies, compared with that of quasars, at magnitudes fainter than M1450 ~ -22 mag or zAB ~24 mag. Follow-up studies of the discovered objects as well as further survey observations are ongoing.

The physical properties and elemental abundances of the interstellar medium in galaxies during cosmic reionization are important for understanding the role of galaxies in this process. We report the Atacama Large Millimeter/submillimeter Array detection of an oxygen emission line at a wavelength of 88 micrometers from a galaxy at an epoch about 700 million years after the Big Bang. The oxygen abundance of this galaxy is estimated at about one-tenth that of the Sun. The non-detection of far-infrared continuum emission indicates a deficiency of interstellar dust in the galaxy. A carbon emission line at a wavelength of 158 micrometers is also not detected, implying an unusually small amount of neutral gas. These properties might allow ionizing photons to escape into the intergalactic medium.

We present the discovery of three protoclusters at $z\sim3\mathrm{-}4$ with spectroscopic confirmation in the Canada-France-Hawaii Telescope (CFHT) Legacy Survey Deep Fields. In these fields, we investigate the large-scale projected sky distribution of $z\sim3\mathrm{-}6$ Lyman break galaxies and identify 21 protocluster candidates from regions that are overdense at more than $4\sigma$ overdensity significance. Based on cosmological simulations, it is expected that more than $76\%$ of these candidates will evolve into a galaxy cluster of at least a halo mass of $10^{14}\,\mathrm{M_\odot}$ at $z=0$. We perform follow-up spectroscopy for eight of the candidates using Subaru/FOCAS, KeckII/DEIMOS, and Gemini-N/GMOS. In total we target 462 dropout candidates and obtain 138 spectroscopic redshifts. We confirm three real protoclusters at $z=3\mathrm{-}4$ with more than five members spectroscopically identified, and find one to be an incidental overdense region by mere chance alignment. The other four candidate regions at $z\sim5\mathrm{-}6$ require more spectroscopic follow-up in order to be conclusive. A $z=3.67$ protocluster, which has eleven spectroscopically confirmed members, shows a remarkable core-like structure composed of a central small region ($<0.5\,\mathrm{physical\>Mpc}$) and an outskirts region ($\sim1.0\,\mathrm{physical\>Mpc}$). The Ly$\alpha$ equivalent widths of members of the protocluster are significantly smaller than those of field galaxies at the same redshift while there is no difference in the UV luminosity distributions. These results imply that some environmental effects start operating as early as at $z\sim4$ along with the growth of the protocluster structure.

Using VLT/X-shooter we searched for emission line galaxies associated to four damped Lyman-$\alpha$ systems (DLAs) and one sub-DLA at 2.73<=z<=3.25 towards QSO J2358+0149. We detect [O III] emission from a "low-cool" DLA at z_abs = 2.9791 (having log N(HI)=21.69+\-0.10, [Zn/H] = -1.83+\-0.18) at an impact parameter of, $\rho$ ~12 kpc. The associated galaxy is compact with a dynamical mass of (1-6)x10^9 M_solar, very high excitation ([O III]/[O II] and [O III]/[H$\beta$] both greater than 10), 12+[O/H]<=8.5 and moderate star formation rate (SFR <=2 M_solar yr^{-1}). Such properties are typically seen in the low-z extreme blue compact dwarf galaxies. The kinematics of the gas is inconsistent with that of an extended disk and the gas is part of either a large scale wind or cold accretion. We detect Ly$\alpha$ emission from the z_abs = 3.2477 DLA (having log N(HI)=21.12+\-0.10 and [Zn/H]=-0.97+\-0.13).The Ly$\alpha$ emission is redshifted with respect to the metal absorption lines by 320 km s^{-1}, consistent with the location of the red hump expected in radiative transport models. We derive SFR ~0.2-1.7 M_solar yr^{-1} and Ly$\alpha$ escape fraction of >=10 per cent. No other emission line is detected from this system. Because the DLA has a small velocity separation from the quasar (~500 km s^{-1}) and the DLA emission is located within a small projected distance ($\rho<5$ kpc), we also explore the possibility that the Ly$\alpha$ emission is being induced by the QSO itself. QSO induced Ly$\alpha$ fluorescence is possible if the DLA is within a physical separation of 340 kpc to the QSO. Detection of stellar continuum light and/or the oxygen emission lines would disfavor this possibility. We do not detect any emission line from the remaining three systems.

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 the results of clustering analysis on the $z \sim 2$ star-forming galaxies. By combining our data with data from publicly available archives, we collect $g$-, $\zb / z$-, and $K$-band imaging data over 5.2 deg$^{2}$, which represents the largest area BzK/gzK survey. We apply colour corrections to translate our filter-set to those used in the original BzK selection for the gzK selection. Because of the wide survey area, we obtain a sample of 41,112 star-forming gzK galaxies at $z \sim 2$ (sgzKs) down to $\KAB < 23.0$, and determine high-quality two-point angular correlation functions (ACFs). Our ACFs show an apparent excess from power-law behaviour at small angular scale $(\theta \la 0.01^{\circ})$, which corresponds the virial radius of a dark halo at $z \sim 2$ with a mass of $\sim 10^{13} \Msun$. We find that the correlation lengths are consistent with the previous estimates over all magnitude range; however, our results are evaluated with a smaller margin of error than that in previous studies. The large amount of data enables us to determine ACFs differentially depending on the luminosity of the subset of the data. The mean halo mass of faint sgzKs $(22.0 < K \leq 23.0)$ was found to be $\Mh = (1.32^{+0.09}_{-0.12}) \times 10^{12} h^{-1} \Msun$, whereas bright sgzKs ($18.0 \leq K \leq 21.0)$ were found to reside in dark haloes with a mass of $\Mh = (3.26^{+1.23}_{-1.02}) \times 10^{13} h^{-1} \Msun$.

We present the results of an halo occupation distribution (HOD) analysis of star-forming galaxies at $z \sim 2$. We obtained high-quality angular correlation functions based on a large sgzK sample, which enabled us to carry out the HOD analysis. The mean halo mass and the HOD mass parameters are found to increase monotonically with increasing $K$-band magnitude, suggesting that more luminous galaxies reside in more massive dark haloes. The luminosity dependence of the HOD mass parameters was found to be the same as in the local Universe; however, the masses were larger than in the local Universe over all ranges of magnitude. This implies that galaxies at $z \sim 2$ tend to form in more massive dark haloes than in the local Universe, a process known as downsizing. By analysing the dark halo mass evolution using the extended Press--Schechter formalism and the number evolution of satellite galaxies in a dark halo, we find that faint Lyman break galaxies at $z \sim 4$ could evolve into the faintest sgzKs $(22.0 < K \leq 23.0)$ at $z \sim 2$ and into the Milky-Way-like galaxies or elliptical galaxies in the local Universe, whereas the most luminous sgzKs $(18.0 \leq K \leq 21.0)$ could evolve into the most massive systems in the local Universe. The stellar-to-halo mass ratio (SHMR) of the sgzKs was found to be consistent with the prediction of the model, except that the SHMR of the faintest sgzKs was smaller than the prediction at $z \sim 2$. This discrepancy may be attributed that our samples are confined to star-forming galaxies.

We revealed the detailed structure of a vastly extended H$\alpha$-emitting nebula ("H$\alpha$ nebula") surrounding the starburst/merging galaxy NGC 6240 by deep narrow-band imaging observations with the Subaru Suprime-Cam. The extent of the nebula is $\sim$90 kpc in diameter and the total H$\alpha$ luminosity amounts to $\approx 1.6 \times 10^{42}$ erg s$^{-1}$. The volume filling factor and the mass of the warm ionized gas are $\sim$10$^{-4}$--10$^{-5}$ and $\sim$$5 \times 10^8$ $M_\odot$, respectively. The nebula has a complicated structure, which includes numerous filaments, loops, bubbles, and knots. We found that there is a tight spatial correlation between the H$\alpha$ nebula and the extended soft X-ray-emitting gas, both in large and small scales. The overall morphology of the nebula is dominated by filamentary structures radially extending from the center of the galaxy. A large-scale bi-polar bubble extends along the minor axis of the main stellar disk. The morphology strongly suggests that the nebula was formed by intense outflows - superwinds - driven by starbursts. We also found three bright knots embedded in a looped filament of ionized gas that show head-tail morphologies in both emission-line and continuum, suggesting close interactions between the outflows and star forming regions. Based on the morphology and surface brightness distribution of the H$\alpha$ nebula, we propose the scenario that three major episodes of starburst/superwind activities which were initiated $\sim$10$^2$ Myr ago formed the extended ionized gas nebula of NGC 6240.

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.

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

The TMT Detailed Science Case describes the transformational science that the Thirty Meter Telescope will enable. Planned to begin science operations in 2024, TMT will open up opportunities for revolutionary discoveries in essentially every field of astronomy, astrophysics and cosmology, seeing much fainter objects much more clearly than existing telescopes. Per this capability, TMT's science agenda fills all of space and time, from nearby comets and asteroids, to exoplanets, to the most distant galaxies, and all the way back to the very first sources of light in the Universe. More than 150 astronomers from within the TMT partnership and beyond offered input in compiling the new 2015 Detailed Science Case. The contributing astronomers represent the entire TMT partnership, including the California Institute of Technology (Caltech), the Indian Institute of Astrophysics (IIA), the National Astronomical Observatories of the Chinese Academy of Sciences (NAOC), the National Astronomical Observatory of Japan (NAOJ), the University of California, the Association of Canadian Universities for Research in Astronomy (ACURA) and US associate partner, the Association of Universities for Research in Astronomy (AURA).

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 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 ALMA observations of the [CII] line and far-infrared (FIR) continuum of a normally star-forming galaxy in the reionization epoch, the z=6.96 Ly-alpha emitter (LAE) IOK-1. Probing to sensitivities of sigma_line = 240 micro-Jy/beam (40 km/s channel) and sigma_cont = 21 micro-Jy/beam, we found the galaxy undetected in both [CII] and continuum. Comparison of UV - FIR spectral energy distribution (SED) of IOK-1, including our ALMA limit, with those of several types of local galaxies (including the effects of the cosmic microwave background, CMB, on the FIR continuum) suggests that IOK-1 is similar to local dwarf/irregular galaxies in SED shape rather than highly dusty/obscured galaxies. Moreover, our 3 sigma FIR continuum limit, corrected for CMB effects, implies intrinsic dust mass M_dust < 6.4 x 10^7 M_sun, FIR luminosity L_FIR < 3.7 x 10^{10} L_sun (42.5 - 122.5 micron), total IR luminosity L_IR < 5.7 x 10^{10} L_sun (8 - 1000 micron) and dust-obscured star formation rate (SFR) < 10 M_sun/yr, if we assume that IOK-1 has a dust temperature and emissivity index typical of local dwarf galaxies. This SFR is 2.4 times lower than one estimated from the UV continuum, suggesting that < 29% of the star formation is obscured by dust. Meanwhile, our 3 sigma [CII] flux limit translates into [CII] luminosity, L_[CII] < 3.4 x 10^7 L_sun. Locations of IOK-1 and previously observed LAEs on the L_[CII] vs. SFR and L_[CII]/L_FIR vs. L_FIR diagrams imply that LAEs in the reionization epoch have significantly lower gas and dust enrichment than AGN-powered systems and starbursts at similar/lower redshifts, as well as local star-forming galaxies.

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.

We searched for star formation activity associated with high-z Damped Lyman-alpha systems (DLAs) with Subaru telescope. We used a set of narrow-band (NB) filters whose central wavelengths correspond to the redshifted Lyman-alpha emission lines of targeted DLA absorbers at 3<z<4.5. We detected one apparent NB-excess object located 3.80 arcsec (~28kpc) away from the quasar SDSS J031036.84+005521.7. Follow-up spectroscopy revealed an asymmetric Lyman-alpha emission at z_em=3.115+/-0.003, which perfectly matches the sub-DLA trough at z_abs=3.1150 with logN(HI)/cm^-2=20.05. The Lyman-alpha luminosity is estimated to be L(LyA)=1.07x10^42 erg s^-1, which corresponds to a star formation rate of 0.97 M_\odot yr^-1. Interestingly, the detected Lyman-alpha emission is spatially extended with a sharp peak. The large extent of the Lyman-alpha emission is remarkably one-sided toward the quasar line-of-sight, and is redshifted. The observed spatially asymmetric surface brightness profile can be qualitatively explained by a model of a DLA host galaxy, assuming a galactic outflow and a clumpy distribution of HI clouds in the circumgalactic medium. This large Lyman-alpha extension, which is similar to those found in Rauch et al. (2008), could be the result of complicated anisotropic radiative transfer through the surrounding neutral gas embedded in the DLA.

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

We present observations of SDF-05M05, an unusual optical transient discovered in the Subaru Deep Field (SDF). The duration of the transient is > ~800 d in the observer frame, and the maximum brightness during observation reached approximately 23 mag in the i' and z' bands. The faint host galaxy is clearly identified in all 5 optical bands of the deep SDF images. The photometric redshift of the host yields z~0.6 and the corresponding absolute magnitude at maximum is ~-20. This implies that this event shone with an absolute magnitude brighter than -19 mag for approximately 300 d in the rest frame, which is significantly longer than a typical supernova and ultra-luminous supernova. The total radiated energy during our observation was 1x10^51 erg. The light curves and color evolution are marginally consistent with some of luminous IIn supernova. We suggest that the transient may be a unique and peculiar supernova at intermediate redshift.

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 derived the luminosity function (LF) of dwarf galaxies in the Coma Cluster down to M_R=-10 at three fields located at the center, intermediate, and outskirt. The LF (-19<M_R<-10) shows no significant differences among the three fields. It shows a clear dip at M_R\sim-13, and is composed of two distinct components of different slopes; the bright component with -19<M_R<-13 has a flatter slope than the faint component with -13<M_R<-10 which has a steep slope. The bright component (-19<M_R<-13) consists of mostly red extended galaxies including few blue galaxies whose colors are typical of late-type galaxies. On the other hand, the faint component (-13<M_R<-10) consists of largely PSF-like compact galaxies. We found that both these compact galaxies and some extended galaxies are present in the center while only compact galaxies are seen in the outskirt. In the faint component, the fraction of blue galaxies is larger in the outskirt than in the center. We suggest that the dwarf galaxies in the Coma Cluster, which make up the two components in the LF, are heterogeneous with some different origins.

We searched for z=7.3 Lya emitters (LAEs) behind two lensing clusters, Abell 2390 and CL 0024, with the Subaru Telescope Suprime-Cam and a narrowband NB1006 (FWHM ~ 21 nm centered at 1005 nm). We investigated if there exist objects consistent with the color of z=7.3 LAEs behind the clusters but could not detect any LAEs to the unlensed line limit F(Lya) ~ 6.9 x 10^{-18} erg/s/cm^2. Using several z=7 Lya luminosity functions (LFs) from the literature, we estimated and compared the expected detection numbers of z ~ 7 LAEs in lensing and blank field surveys in the case of using an 8m class ground based telescope. Given the steep bright-end slope of the LFs, when the detector field-of view (FOV) is comparable to the angular extent of a massive lensing cluster, imaging cluster(s) is more efficient in detecting z ~ 7 LAEs than imaging a blank field. However, the gain is expected to be modest, a factor of two at most and likely much less depending on the adopted LFs. The main advantage of lensing-cluster survey, therefore, remains to be the gain in depth and not necessarily in detection efficiency. For much larger detectors, the lensing effect becomes negligible and the efficiency of LAE detection is proportional to the instrumental FOV. We also inspected NB1006 images of three z ~ 7 z-dropouts previously detected in Abell 2390 and found that none of them are detected in NB1006. Two of them are consistent with predictions from the previous studies that they would be at lower redshifts. The other one has a photometric redshift of z ~ 7.3, and if it is at z=7.3, its unlensed Lya line flux would be very faint: F(Lya) < 4.4 x 10^{-18} erg/s/cm^2 (1 sigma upper limit) or rest frame equivalent width of W(Lya) < 26A. Its Lya emission might be attenuated by neutral hydrogen, as recent studies show that the fraction of Lyman break galaxies displaying strong Lya emission is lower at z ~ 7 than at z <~ 6.

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.

Detailed analysis of the substructure of Lya nebulae can put important constraints on the physical mechanisms at work and the properties of galaxies forming within them. Using high resolution HST imaging of a Lya nebula at z~2.656, we have taken a census of the compact galaxies in the vicinity, used optical/near-infrared colors to select system members, and put constraints on the morphology of the spatially-extended emission. The system is characterized by (a) a population of compact, low luminosity (~0.1 L*) sources --- 17 primarily young, small (Re~1-2 kpc), disky galaxies including an obscured AGN --- that are all substantially offset (>20 kpc) from the line-emitting nebula; (b) the lack of a central galaxy at or near the peak of the Lya emission; and (c) several nearly coincident, spatially extended emission components --- Lya, HeII, and UV continuum --- that are extremely smooth. These morphological findings are difficult to reconcile with theoretical models that invoke outflows, cold flows, or resonant scattering, suggesting that while all of these physical phenomena may be occurring, they are not sufficient to explain the powering and large extent of Lya nebulae. In addition, although the compact galaxies within the system are irrelevant as power sources, the region is significantly overdense relative to the field galaxy population (by at least a factor of 4). These observations provide the first estimate of the luminosity function of galaxies within an individual Lya nebula system, and suggest that large Lya nebulae may be the seeds of galaxy groups or low-mass clusters.

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.

We present the results of deep imaging and spectroscopic observations of very extended ionized gas (EIG) around four member galaxies of the Coma cluster of galaxies: RB199, IC4040, GMP2923 and GMP3071. The EIGs were serendipitously found in an H-alpha narrow band imaging survey of the central region of the Coma cluster. The relative radial velocities of the EIGs with respect to the systemic velocities of the parent galaxies from which they emanate increase almost monotonically with the distance from the nucleus of the respective galaxies, reaching -400 - -800 km/s at around 40 - 80 kpc from the galaxies. The one-sided morphologies and the velocity fields of the EIGs are consistent with the predictions of numerical simulations of ram pressure stripping. We found a very low-velocity filament (v_rel = -1300 km/s) at the southeastern edge of the disk of IC4040. Some bright compact knots in the EIGs of RB199 and IC4040 exhibit blue continuum and strong H-alpha emission. The equivalent widths of the H-alpha emission exceed 200 A, and are greater than 1000 A for some knots. The emission line intensity ratios of the knots are basically consistent with those of sub-solar abundance HII regions. These facts indicate that intensive star formation occurs in the knots. Some filaments, including the low velocity filament of the IC4040 EIG, exhibit shock-like emission line spectra, suggesting that shock heating plays an important role in ionization and excitation of the EIGs.

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.

Using the sample of Type Ia supernovae (SNe Ia) discovered by the Hubble Space Telescope (HST) Cluster Supernova Survey and augmented with HST-observed SNe Ia in the GOODS fields, we search for correlations between the properties of SNe and their host galaxies at high redshift. We use galaxy color and quantitative morphology to determine the red sequence in 25 clusters and develop a model to distinguish passively evolving early-type galaxies from star-forming galaxies in both clusters and the field. With this approach, we identify six SN Ia hosts that are early-type cluster members and eleven SN Ia hosts that are early-type field galaxies. We confirm for the first time at z>0.9 that SNe Ia hosted by early-type galaxies brighten and fade more quickly than SNe Ia hosted by late-type galaxies. We also show that the two samples of hosts produce SNe Ia with similar color distributions. The relatively simple spectral energy distributions (SEDs) expected for passive galaxies enable us to measure stellar masses of early-type SN hosts. In combination with stellar mass estimates of late-type GOODS SN hosts from Thomson & Chary (2011), we investigate the correlation of host mass with Hubble residual observed at lower redshifts. Although the sample is small and the uncertainties are large, a hint of this relation is found at z>0.9. By simultaneously fitting the average cluster galaxy formation history and dust content to the red-sequence scatters, we show that the reddening of early-type cluster SN hosts is likely E(B-V) <~ 0.06. The similarity of the field and cluster early-type host samples suggests that field early-type galaxies that lie on the red sequence may also be minimally affected by dust. Hence, the early-type hosted SNe Ia studied here occupy a more favorable environment to use as well-characterized high-redshift standard candles than other SNe Ia.

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 a measurement of the Type Ia supernova (SN Ia) rate in galaxy clusters at 0.9 < z < 1.45 from the Hubble Space Telescope (HST) Cluster Supernova Survey. This is the first cluster SN Ia rate measurement with detected z > 0.9 SNe. Finding 8 +/- 1 cluster SNe Ia, we determine a SN Ia rate of 0.50 +0.23-0.19 (stat) +0.10-0.09 (sys) SNuB (SNuB = 10^-12 SNe L_{sun,B}^-1 yr^-1). In units of stellar mass, this translates to 0.36 +0.16-0.13 (stat) +0.07-0.06 (sys) SNuM (SNuM = 10^-12 SNe M_sun^-1 yr^-1). This represents a factor of approximately 5 +/- 2 increase over measurements of the cluster rate at z < 0.2. We parameterize the late-time SN Ia delay time distribution with a power law (proportional to t^s). Under the assumption of a cluster formation redshift of z_f = 3, our rate measurement in combination with lower-redshift cluster SN Ia rates constrains s = -1.41 +0.47/-0.40, consistent with measurements of the delay time distribution in the field. This measurement is generally consistent with expectations for the "double degenerate" scenario and inconsistent with some models for the "single degenerate" scenario predicting a steeper delay time distribution at large delay times. We check for environmental dependence and the influence of younger stellar populations by calculating the rate specifically in cluster red-sequence galaxies and in morphologically early-type galaxies, finding results similar to the full cluster rate. Finally, the upper limit of one host-less cluster SN Ia detected in the survey implies that the fraction of stars in the intra-cluster medium is less than 0.47 (95% confidence), consistent with measurements at lower redshifts.

We present a measurement of the volumetric Type Ia supernova (SN Ia) rate out to z ~ 1.6 from the Hubble Space Telescope Cluster Supernova Survey. In observations spanning 189 orbits with the Advanced Camera for Surveys we discovered 29 SNe, of which approximately 20 are SNe Ia. Twelve of these SNe Ia are located in the foregrounds and backgrounds of the clusters targeted in the survey. Using these new data, we derive the volumetric SN Ia rate in four broad redshift bins, finding results consistent with previous measurements at z > 1 and strengthening the case for a SN Ia rate that is equal to or greater than ~0.6 x 10^-4/yr/Mpc^3 at z ~ 1 and flattening out at higher redshift. We provide SN candidates and efficiency calculations in a form that makes it easy to rebin and combine these results with other measurements for increased statistics. Finally, we compare the assumptions about host-galaxy dust extinction used in different high-redshift rate measurements, finding that different assumptions may induce significant systematic differences between measurements.

We report a search for z=7 Ly-alpha emitters (LAEs) using a custom-made Narrow-Band filter, centered at 9755 Angstroms, with the instrument Suprime-Cam installed at the Subaru telescope. We observed two different fields and obtained two sample of 7 Ly-alpha emitters of which 4 are robust in each field. We are covering the luminosity range of 9.10^{42} - 2.10^{43} erg/s in comoving volumes of ~ 4.10^{5} and 4.3.10^{5} Mpc^{3}. From this result, we derived possible z~7 Ly-alpha luminosity functions for the full samples and for a subsample of 4 objects in each field. We do not observe, in each case, any strong evolution between the z=6.5 and z~7 Ly-alpha luminosity functions. Spectroscopic confirmation for these candidate samples is required to establish a definitive measure of the luminosity function at z~7.

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 present ACS, NICMOS, and Keck AO-assisted photometry of 20 Type Ia supernovae SNe Ia from the HST Cluster Supernova Survey. The SNe Ia were discovered over the redshift interval 0.623 < z < 1.415. Fourteen of these SNe Ia pass our strict selection cuts and are used in combination with the world's sample of SNe Ia to derive the best current constraints on dark energy. Ten of our new SNe Ia are beyond redshift $z=1$, thereby nearly doubling the statistical weight of HST-discovered SNe Ia beyond this redshift. Our detailed analysis corrects for the recently identified correlation between SN Ia luminosity and host galaxy mass and corrects the NICMOS zeropoint at the count rates appropriate for very distant SNe Ia. Adding these supernovae improves the best combined constraint on the dark energy density \rho_{DE}(z) at redshifts 1.0 < z < 1.6 by 18% (including systematic errors). For a LambdaCDM universe, we find \Omega_\Lambda = 0.724 +0.015/-0.016 (68% CL including systematic errors). For a flat wCDM model, we measure a constant dark energy equation-of-state parameter w = -0.985 +0.071/-0.077 (68% CL). Curvature is constrained to ~0.7% in the owCDM model and to ~2% in a model in which dark energy is allowed to vary with parameters w_0 and w_a. Tightening further the constraints on the time evolution of dark energy will require several improvements, including high-quality multi-passband photometry of a sample of several dozen z>1 SNe Ia. We describe how such a sample could be efficiently obtained by targeting cluster fields with WFC3 on HST.

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.

Faint Object Camera and Spectrograph (FOCAS) is a versatile common-use optical instrument for the 8.2m Subaru Telescope, offering imaging and spectroscopic observations. FOCAS employs grisms with resolving powers ranging from 280 to 8200 as dispersive optical elements. A grism is a direct-vision grating composed of a transmission grating and prism(s). FOCAS has five grisms with replica surface-relief gratings including an echelle-type grism, and eight grisms with volume-phase holographic (VPH) gratings. The size of these grisms is 110 mm X 106 mm in aperture with a maximum thickness of 110 mm. We employ not only the dichromated gelatin, but also the hologram resin as a recording material for VPH gratings. We discuss the performance of these FOCAS grisms measured in the laboratory, and verify it by test observations, and show examples of astronomical spectroscopic observations.

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 images of extended H-alpha clouds associated with 14 member galaxies in the Coma cluster obtained from deep narrow band imaging observations with Suprime-Cam at the Subaru Telescope. The parent galaxies of the extended H-alpha clouds are distributed farther than 0.2 Mpc from the peak of X-ray emission of the cluster. Most of the galaxies have colors bluer than g-r approx 0.5 and they account for 57% of the blue (g-r<0.5) bright (r<17.8 mag) galaxies in the central region of the Coma cluster. They reside near the red- and blue-shifted edges of the Coma cluster's radial velocity distribution. These findings suggest that the most of the parent galaxies were recently captured by the Coma cluster potential and are now infalling toward the cluster center with their disk gas being stripped off and producing the observed H-alpha clouds.

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

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.

QSOs have been thought to be important for tracing highly biased regions in the early universe, from which the present-day massive galaxies and galaxy clusters formed. While overdensities of star-forming galaxies have been found around QSOs at 2<z<5, the case for excess galaxy clustering around QSOs at z>6 is less clear. Previous studies with HST have reported the detection of small excesses of faint dropout galaxies in some QSO fields, but these surveys probed a relatively small region surrounding the QSOs. To overcome this problem, we have observed the most distant QSO at z=6.4 using the large field of view of the Suprime-Cam (34' x 27'). Newly-installed CCDs allowed us to select Lyman break galaxies (LBG) at z~6.4 more efficiently. We found seven LBGs in the QSO field, whereas only one exists in a comparison field. The significance of this apparent excess is difficult to quantify without spectroscopic confirmation and additional control fields. The Poisson probability to find seven objects when one expects four is ~10%, while the probability to find seven objects in one field and only one in the other is less than 0.4%, suggesting that the QSO field is significantly overdense relative to the control field. We find some evidence that the LBGs are distributed in a ring-like shape centered on the QSO with a radius of ~3 Mpc. There are no candidate LBGs within 2 Mpc from the QSO, i.e., galaxies are clustered around the QSO but appear to avoid the very center. These results suggest that the QSO is embedded in an overdense region when defined on a sufficiently large scale. This suggests that the QSO was indeed born in a massive halo. The central deficit of galaxies may indicate that (1) the strong UV radiation from the QSO suppressed galaxy formation in its vicinity, or (2) that star-formation closest to the QSO occurs mostly in an obscured mode that is missed by our UV selection.

We report the results of a spectropolarimetric observation of the C IV mini-BAL in the quasar HS1603+3820. The observations were carried out with the FOCAS instrument on the Subaru telescope and yielded an extremely high polarization sensitivity of 0.1%, at a resolving power of 1500. HS1603+3820 has been the target of a high-resolution spectroscopic monitoring campaign for more than four years, aimed at studying its highly variable C IV mini-BAL profile. Using the monitoring observations, in an earlier paper we were able to narrow down the causes of the variability to the following two scenarios: (1) scattering material of variable optical depth redirecting photons around the absorber, and (2) a variable, highly-ionized screen between the continuum source and the absorber which modulates the UV continuum incident on the absorber. The observations presented here provide a crucial test of the scattering scenario and lead us to disfavor it because (a) the polarization level is very small (p~0.6%) throughout the spectrum, and (b) the polarization level does not increase across the mini-BAL trough. Thus, the variable screen scenario emerges as our favored explanation of the C IV mini-BAL variability. Our conclusion is bolstered by recent X-ray observations of nearby mini-BAL quasars, which show a rapidly variable soft X-ray continuum that appears to be the result of transmission through an ionized absorber of variable ionization parameter and optical depth.

We report on work to increase the number of well-measured Type Ia supernovae (SNe Ia) at high redshifts. Light curves, including high signal-to-noise HST data, and spectra of six SNe Ia that were discovered during 2001 are presented. Additionally, for the two SNe with z>1, we present ground-based J-band photometry from Gemini and the VLT. These are among the most distant SNe Ia for which ground based near-IR observations have been obtained. We add these six SNe Ia together with other data sets that have recently become available in the literature to the Union compilation (Kowalski et al. 2008). We have made a number of refinements to the Union analysis chain, the most important ones being the refitting of all light curves with the SALT2 fitter and an improved handling of systematic errors. We call this new compilation, consisting of 557 supernovae, the Union2 compilation. The flat concordance LambdaCDM model remains an excellent fit to the Union2 data with the best fit constant equation of state parameter w=-0.997^{+0.050}_{-0.054} (stat) ^{+0.077}_{-0.082} (stat+sys\ together) for a flat universe, or w=-1.035^{+0.055}_{-0.059} (stat)^{+0.093}_{-0.097} (stat+sys together) with curvature. We also present improved constraints on w(z). While no significant change in w with redshift is detected, there is still considerable room for evolution in w. The strength of the constraints depend strongly on redshift. In particular, at z > 1, the existence and nature of dark energy are only weakly constrained by the data.

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 spectra of high-redshift supernovae (SNe) that were taken with the Subaru low resolution optical spectrograph, FOCAS. These SNe were found in SN surveys with Suprime-Cam on Subaru, the CFH12k camera on the Canada-France-Hawaii Telescope (CFHT), and the Advanced Camera for Surveys (ACS) on the Hubble Space Telescope (HST). These SN surveys specifically targeted z>1 Type Ia supernovae (SNe Ia). From the spectra of 39 candidates, we obtain redshifts for 32 candidates and spectroscopically identify 7 active candidates as probable SNe Ia, including one at z=1.35, which is the most distant SN Ia to be spectroscopically confirmed with a ground-based telescope. An additional 4 candidates are identified as likely SNe Ia from the spectrophotometric properties of their host galaxies. Seven candidates are not SNe Ia, either being SNe of another type or active galactic nuclei. When SNe Ia are observed within a week of maximum light, we find that we can spectroscopically identify most of them up to z=1.1. Beyond this redshift, very few candidates were spectroscopically identified as SNe Ia. The current generation of super red-sensitive, fringe-free CCDs will push this redshift limit higher.

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 present mid-infrared imaging at 11.7 $\mu$m for the quadruple lens systems, MG0414+0534 and Q2237+030, using the cooled mid-infrared camera and spectrometer (COMICS) attached on the Subaru telescope. MG0414+0534 is characterized by a bright pair of lensed images (A1, A2) and their optical flux ratio A2/A1 deviates significantly from the prediction of a smooth lens model. Q2237+030 is comprised of four lensed images, which are significantly affected by microlensing in a foreground lensing galaxy. Our mid-infrared observations of these lensed images have revealed that the mid-infrared flux ratio for A2/A1 of MG0414+0534 is nearly unity (0.90 +- 0.04). We find that this flux ratio is systematically small, at 4 to 5 sigma level, compared with the prediction of a best smooth lens model (1.09) represented by a singular isothermal ellipsoid and external shear, thereby suggesting the presence of substructures to explain our observational result. In contrast, for Q2237+030, our high signal-to-noise observation indicates that the mid-infrared flux ratios between all the four images of Q2237+030 are virtually consistent with the prediction of a smooth lens model. Based on the size estimate of the dust torus surrounding the nuclei of these QSOs, we set limits on the mass of a substructure in these lens systems, which can cause anomalies in the flux ratios. For MG0414+0534, since the required mass of a substructure inside its Einstein radius is > 360 M_sun, millilensing by a CDM substructure is most likely. If it is modeled as a singular isothermal sphere, the mass inside radius of 100 pc is given as > 1.0x10^5 M_sun. For Q2237+030, there is no significant evidence of millilensing, so the reported anomalous flux ratios in shorter wavelengths are entirely caused by microlensing by stars (abridged).

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 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 the first results of a MOS campaign to follow up cluster candidates located via weak lensing. Our main goals are to search for spatial concentrations of galaxies that are plausible optical counterparts of the weak lensing signals, and to determine the cluster redshifts from those of member galaxies. Around each of 36 targeted cluster candidates, we obtain 15-32 galaxy redshifts. For 28 of these targets, we confirm a secure cluster identification. This includes three cases where two clusters at different redshifts are projected along the same line-of-sight. In 6 of the 8 unconfirmed targets, we find multiple small galaxy concentrations at different redshifts. In both the remaining two targets, a single small galaxy concentration is found. We evaluate the weak lensing mass of confirmed clusters. For a subsample of our most cleanly measured clusters, we investigate the statistical relation between their weak lensing mass and the velocity dispersion of their member galaxies, comparing our sample with optically and X-ray selected samples from the literature. Our lensing-selected clusters are consistent with sigma_v=sigma_sis, with a similar scatter to the optically and X-ray selected clusters. We thus find no evidence of selection bias compared to these other techniques. We also derive an empirical relation between the cluster mass and the galaxy velocity dispersion, which is in reasonable agreement with the prediction of N-body simulations in the LCDM cosmology.

We found an unusual complex of narrow blue filaments, bright blue knots, and H-alpha emitting filaments and clouds, which morphologically resembled a complex of ``fireballs,'' extending up to 80 kpc south from an E+A galaxy RB199 in the Coma cluster. The galaxy has a highly disturbed morphology indicative of a galaxy--galaxy merger remnant. The narrow blue filaments extend in straight shapes toward the south from the galaxy, and several bright blue knots are located at the southern ends of the filaments. The Rc band absolute magnitudes, half light radii and estimated masses of the bright knots are -12 - -13 mag, 200 - 300 pc and 10^6-7 Msolar, respectively. Long, narrow H-alpha emitting filaments are connected at the south edge of the knots. The average color of the fireballs is B - Rc = 0.5, which is bluer than RB199 (B - R = 0.99), suggesting that most of the stars in the fireballs were formed within several times 10^8 yr. The narrow blue filaments exhibit almost no H-alpha emission. Strong H-alpha and UV emission appear in the bright knots. These characteristics indicate that star formation recently ceased in the blue filaments and now continues in the bright knots. The gas stripped by some mechanism from the disk of RB199 may be traveling in the intergalactic space, forming stars left along its trajectory. The most plausible fireball formation mechanism is ram pressure stripping by high-speed collision between the galaxy and the hot intra-cluster medium. The fireballs may be a snapshot of diffuse intra-cluster population formation, or halo star population formation in a cluster galaxy.

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.

To set useful limits on the abundance of small-scale dark matter halos (subhalos) in a galaxy scale, we have carried out mid-infrared imaging and integral-field spectroscopy for a sample of quadruple lens systems showing anomalous flux ratios. These observations using Subaru have been successful for distinguishing millilensing by subhalos from microlensing by stars. Current status for our lensing analysis of dark matter substructure is reported.

Large nebulae (>50 kpc) emitting strongly in Lya (also known as Lya ``blobs'') are likely signposts of ongoing massive galaxy formation. The relative rarity of these sources and their discovery in well-studied galaxy overdensities suggest that they may be associated with regions of high galaxy density. One of the largest Lya nebulae, discovered at a redshift of z~2.7 via its strong mid-infrared emission, provides an unbiased test of this association. We have carried out a deep intermediate-band imaging survey for Lya-emitting galaxies (LAEs) within a 30'x26' field of view around this Lya nebula. This is the first study of the environment of a Lya nebula found without a priori knowledge of its surroundings. We find that the nebula is located in an overdense region, at least 20x50 1/h_70 comoving Mpc in size, showing a factor of ~3 LAE number density enhancement relative to the edge of the field. Given the predicted number of such overdensities, we rule out the possibility of a chance coincidence at the <1% level. This study, in conjunction with previous work, provides strong confirmation of the association between the largest Lya nebulae and overdense regions of the Universe.

We present a new photometric search for high-z galaxies hosting Population III (PopIII) stars based on deep intermediate-band imaging observations obtained in the Subaru Deep Field (SDF), by using Suprime-Cam on the Subaru Telescope. By combining our new data with the existing broad-band and narrow-band data, we searched for galaxies which emit strongly both in Ly_alpha and in HeII 1640 (``dual emitters'') that are promising candidates for PopIII-hosting galaxies, at 3.93<z<4.01 and 4.57<z<4.65. Although we found 10 ``dual emitters'', most of them turn out to be [OII]-[OIII] dual emitters or H_beta-(H_alpha+[NII]) dual emitters at z<1, as inferred from their broad-band colors and from the ratio of the equivalent widths. No convincing candidate of Ly_alpha-HeII dual emitter of SFR_PopIII > 2 Msun/yr was found by our photometric search in 4.03 x 10^5 Mpc^3 in the SDF. This result disfavors low feedback models for PopIII star clusters, and implies an upper-limit of the PopIII SFR density of SFRD_PopIII < 5 x 10^-6 Msun/yr/Mpc^3. This new selection method to search for PopIII-hosting galaxies should be useful in future narrow-band surveys to achieve the first observational detection of PopIII-hosting galaxies at high redshifts.

We present our survey for optically faint variable objects using multi-epoch (8-10 epochs over 2-4 years) $i'$-band imaging data obtained with Subaru Suprime-Cam over 0.918 deg$^2$ in the Subaru/XMM-Newton Deep Field (SXDF). We found 1040 optically variable objects by image subtraction for all the combinations of images at different epochs. This is the first statistical sample of variable objects at depths achieved with 8-10m class telescopes or HST. The detection limit for variable components is $i'_{\rm{vari}}\sim25.5$ mag. These variable objects were classified into variable stars, supernovae (SNe), and active galactic nuclei (AGN), based on the optical morphologies, magnitudes, colors, and optical-mid-infrared colors of the host objects, spatial offsets of variable components from the host objects, and light curves. Detection completeness was examined by simulating light curves for periodic and irregular variability. We detected optical variability for $36\pm2%$ ($51\pm3%$ for a bright sample with $i'<24.4$ mag) of X-ray sources in the field. Number densities of variable obejcts as functions of time intervals $\Delta{t}$ and variable component magnitudes $i'_{\rm{vari}}$ are obtained. Number densities of variable stars, SNe, and AGN are 120, 489, and 579 objects deg$^{-2}$, respectively. Bimodal distributions of variable stars in the color-magnitude diagrams indicate that the variable star sample consists of bright ($V\sim22$ mag) blue variable stars of the halo population and faint ($V\sim23.5$ mag) red variable stars of the disk population. There are a few candidates of RR Lyrae providing a possible number density of $\sim10^{-2}$ kpc$^{-3}$ at a distance of $>150$ kpc from the Galactic center.

We present the properties of active galactic nuclei (AGN) selected by optical variability in the Subaru/XMM-Newton Deep Field (SXDF). Based on the locations of variable components and light curves, 211 optically variable AGN were reliably selected. We made three AGN samples; X-ray detected optically non-variable AGN (XA), X-ray detected optically variable AGN (XVA), and X-ray undetected optically variable AGN (VA). In the VA sample, we found a bimodal distribution of the ratio between the variable component flux and the host flux. One of these two components in the distribution, a class of AGN with a faint variable component $i'_{\rm{vari}}\sim25$ mag in bright host galaxies $i'\sim21$ mag, is not seen in the XVA sample. These AGN are expected to have low Eddington ratios if we naively consider a correlation between bulge luminosity and black hole mass. These galaxies have photometric redshifts $z_{\rm{photo}}\sim0.5$ and we infer that they are low-luminosity AGN with radiatively inefficient accretion flows (RIAFs). The properties of the XVA and VA objects and the differences from those of the XA objects can be explained within the unified scheme for AGN. Optical variability selection for AGN is an independent method and could provide a complementary AGN sample which even deep X-ray surveys have not found.

The nebular spectra of the broad-lined, SN 1998bw-like Type Ic SN 2002ap are studied by means of synthetic spectra. Two different modelling techniques are employed. In one technique, the SN ejecta are treated as a single zone, while in the other a density and abundance distribution in velocity is used from an explosion model. In both cases, heating caused by gamma-ray and positron deposition is computed (in the latter case using a Monte Carlo technique to describe the propagation of gamma-rays and positrons), as is cooling via forbidden-line emission. The results are compared, and although general agreement is found, the stratified models are shown to reproduce the observed line profiles much more accurately than the single-zone model. The explosion produced ~ 0.1 Msun of 56Ni. The distribution in velocity of the various elements is in agreement with that obtained from the early-time models, which indicated an ejected mass of ~ 2.5 Msun with a kinetic energy of 4 x 10^{51} erg. Nebular spectroscopy confirms that most of the ejected mass (~ 1.2 Msun) was oxygen. The presence of an oxygen-rich inner core, combined with that of 56Ni at high velocities as deduced from early-time models, suggests that the explosion was asymmetric, especially in the inner part.

We present the results of an ongoing weak lensing survey conducted with the Subaru telescope whose initial goal is to locate and study the distribution of shear-selected structures or halos. Using a Suprime-cam imaging survey spanning 21.82 square degree, we present a catalog of 100 candidate halos located from lensing convergence maps. Our sample is reliably drawn from that subset of our survey area, (totaling 16.72 square degree) uncontaminated by bright stars and edge effects and limited at a convergence signal to noise ratio of 3.69. To validate the sample detailed spectroscopic measures have been made for 26 candidates using the Subaru multi-object spectrograph, FOCAS. All are confirmed as clusters of galaxies but two arise as the superposition of multiple clusters viewed along the line of sight. Including data available in the literature and an ongoing Keck spectroscopic campaign, a total of 41 halos now have reliable redshifts. For one of our survey fields, the XMM LSS (Pierre et al. 2004) field, we compare our lensing-selected halo catalog with its X-ray equivalent. Of 15 halos detected in the XMM-LSS field, 10 match with published X-ray selected clusters and a further 2 are newly-detected and spectroscopically confirmed in this work. Although three halos have not yet been confirmed, the high success rate within the XMM-LSS field (12/15) confirms that weak lensing provides a reliable method for constructing cluster catalogs, irrespective of the nature of the constituent galaxies or the intracluster medium.

We list and analyze H I absorption lines at redshifts 2 < z < 4 with column density (12 < log(N_HI) < 19) in 40 high-resolutional (FWHM = 8.0 km/s) quasar spectra obtained with the Keck+HIRES. We de-blend and fit all H I lines within 1,000 km/s of 86 strong H I lines whose column densities are log(N_HI/[cm^-2]) > 15. Unlike most prior studies, we use not only Lya but also all visible higher Lyman series lines to improve the fitting accuracy. This reveals components near to higher column density systems that can not be seen in Lya. We list the Voigt profile fits to the 1339 H I components that we found. We examined physical properties of H I lines after separating them into several sub-samples according to their velocity separation from the quasars, their redshift, column density and the S/N ratio of the spectrum. We found two interesting trends for lines with 12 < log(N_HI) < 15 which are within 200-1000 km/s of systems with log(N_HI) > 15. First, their column density distribution becomes steeper, meaning relatively fewer high column density lines, at z < 2.9. Second, their column density distribution also becomes steeper and their line width becomes broader by about 2-3 km/s when they are within 5,000 km/s of their quasar.

We report our first results of photometric and spectroscopic observations for Tycho's supernova remnant (SNR Tycho) to search for the companion star of a type Ia supernova (SN Ia). From photometric observations using Suprime-Cam on the Subaru Telescope, we have picked up stars brighter than 22 mag (in $V$-band) for spectroscopy, which are located within a circular region with the radius of 30" around the center of SNR Tycho. If the ejecta of young supernova remnants, such as SNR Tycho, have a sufficient amount of Fe I, we should be able to detect absorption lines at 3720 \AA and 3860 \AA associated with transitions from the ground state of Fe I in the spectrum of the companion star. To identify the companion star of a SN Ia using these characteristic absorption lines of Fe I, we made optical low-resolution spectroscopy of their targets using FOCAS on the Subaru Telescope. In our spectroscopic observations, we obtained spectra of 17 stars in the SNR Tycho region and compare them with template stellar spectra. We detect significant absorption lines from two stars at 3720 \AA. Since widths of their absorption lines are broad, it is likely that the detected absorptions are due to Fe I in the expanding ejecta of SNR Tycho. However, none of stars exhibits a clear red wing in the observed profiles of the absorption, though a star in the background of the SNR should show it. Hence, we suggest another interpretation that the detected absorption lines might be caused by the peculiarity of stars. A star named Tycho(E) has the absorption line at 3720 \AA and its projected position is close to the center of SNR Tycho. Based on our observations, Tycho(E) is a new candidate as the companion star of Tycho's supernova.

We report a definitive confirmation of a large-scale structure around the super rich cluster CL0016+1609 at z=0.55. We made spectroscopic follow-up observations with FOCAS on Subaru along the large filamentary structure identified in our previous photometric studies. We have confirmed the physical connection of the huge filament extending over 20 Mpc in the N-S direction, and another filament extending from the main cluster to the East. Based on a simple energy argument, we show that it is likely that most of the clumps are bound to the main CL0016 cluster. This structure is surely one of the most prominent confirmed structures ever identified in the distant Universe, which then serves as an ideal laboratory to examine the environmental variation of galaxy properties. We draw star formation histories of galaxies from the composite spectra of red galaxies in field, group, and cluster environments. Combining the results from our previous studies, we find that red galaxies in groups at z~0.8 and red field galaxies at z~0.5 show strong Hd absorption lines for their D4000 indices. These are the environments in which we observed the on-going build-up of the colour-magnitude relation in our previous photometric analyses. The strong Hd absorption lines imply that their star formation is truncated on a relatively short time scale. We suggest that a galaxy-galaxy interaction is the most likely physical driver of the truncation of star formation and thus responsible for the build-up of the colour-magnitude relation since z~1. (Abridged)

We carried out a target survey for Lyman break galaxies (LBGs) and Lyman alpha emitters (LAEs) around QSO SDSS J0211-0009 at z=4.87. The deep and wide broadband and narrowband imaging simultaneously revealed the perspective structure of these two high-z populations. The LBGs without Ly-alpha emission form a filamentary structure including the QSO, while the LAEs are distributed around the QSO but avoid it within a distance of ~4.5Mpc. On the other hand, we serendipitously discovered a protocluster with a significant concentration of LBGs and LAEs where no strongly UV ionizing source such as a QSO or radio galaxy is known to exist. In this cluster field, two populations are spatially cross-correlated with each other. The relative spatial distribution of LAEs to LBGs is in stark contrast between the QSO and the cluster fields. We also found a weak trend showing that the number counts based on Ly-alpha and UV continuum fluxes of LAEs in the QSO field are slightly lower than in the cluster field, whereas the number counts of LBGs are almost consistent with each other. The LAEs avoid the nearby region around the QSO where the local UV background radiation could be ~100 times stronger than the average for the epoch. The clustering segregation between LBGs and LAEs seen in the QSO field could be due to either enhanced early galaxy formation in an overdense environment having caused all the LAEs to evolve into LBGs, or local photoionization due to the strong UV radiation from the QSO effectively causing a deficit in low-mass galaxies like LAEs.

We present Keck diffraction limited H-band photometry of a z~1.3 Type Ia supernova (SN) candidate, first identified in a Hubble Space Telescope (HST) search for SNe in massive high redshift galaxy clusters. The adaptive optics (AO) data were obtained with the Laser Guide Star facility during four observing runs from September to November 2005. In the analysis of data from the observing run nearest to maximum SN brightness, the SN was found to have a magnitude H=23.9 +/- 0.14 (Vega). We present the H-band (approximately rest-frame R) light curve and provide a detailed analysis of the AO photometric uncertainties. By constraining the aperture correction with a nearby (4" separation) star we achieve 0.14 magnitude photometric precision, despite the spatially varying AO PSF.

The Coma cluster is the richest and most compact of the nearby clusters, yet there is growing evidence that its formation is still on-going. With a new multi-slit imaging spectroscopy technique pioneered at the 8.2 m Subaru telescope and FOCAS, we have detected and measured the line-of-sight velocities of 37 intracluster planetary nebulae associated with the diffuse stellar population of stars in the Coma cluster core, at 100 Mpc distance. We detect clear velocity substructures within a 6 arcmin diameter field. A substructure is present at ~5000 km/s, probably from in-fall of a galaxy group, while the main intracluster stellar component is centered around ~6500 km/s, ~700 km/s offset from the nearby cD galaxy NGC 4874. The kinematics and morphology of the intracluster stars show that the cluster core is in a highly dynamically evolving state. In combination with galaxy redshift and X-ray data this argues strongly that the cluster is currently in the midst of a subcluster merger, where the NGC 4874 subcluster core may still be self-bound, while the NGC 4889 subcluster core has probably dissolved. The NGC 4889 subcluster is likely to have fallen into Coma from the eastern A2199 filament, in a direction nearly in the plane of the sky, meeting the NGC 4874 subcluster arriving from the west. The two inner subcluster cores are presently beyond their first and second close passage, during which the elongated distribution of diffuse light has been created. We predict the kinematic signature expected in this scenario, and argue that the extended western X-ray arc recently discovered traces the arc shock generated by the collision between the two subcluster gas halos. Any preexisting cooling core region would have been heated by the subcluster collision.

We develop a method to measure the strength of the absorption features in Type Ia supernova (SN Ia) spectra and use it to make a quantitative comparison between the spectra of Type Ia supernovae at low and high redshifts. In this case study, we apply the method to 12 high-redshift (0.212 < z < 0.912) SNe Ia observed by the Supernova Cosmology Project . Through measurements of the strengths of these features and of the blueshift of the absorption minimum in Ca II H&K, we show that the spectra of the high-redshift SNe Ia are quantitatively similar to spectra of nearby SNe Ia (z < 0.15). One supernova in our high redshift sample, SN 2002fd at z=0.279, is found to have spectral characteristics that are associated with peculiar SN 1991T/SN 1999aa-like supernovae.

The Coma cluster is the richest and most compact of the nearby clusters, yet there is growing evidence that its formation is still on-going. A sensitive probe of this evolution is the dynamics of intracluster stars, which are unbound from galaxies while the cluster forms, according to cosmological simulations. With a new multi-slit imaging spectroscopy technique pioneered at the 8.2 m Subaru telescope and FOCAS, we can now detect and measure the line-of-sight velocities of the intracluster planetary nebulae which are associated with the diffuse stellar population of stars, at 100 Mpc distance. We detect significant velocity substructures within a 6 arcmin diameter field, centred on the Coma X-ray cluster emission. One substructure ispresent at ~5000 km/s, probably from infall of a galaxy group, while the main intracluster stellar component moves at ~ 6500 km/s. Hence the ICPNs associated with the diffuse light at the position of the MSIS field are not bound to the nearby cD galaxy NGC 4874, whose radial velocity is ~ 700 km/s greater.

We have observed the quadruply lensed quasar 1RXS J1131-1231 with the integral field spectrograph mode of the Kyoto Tridimensional Spectrograph II mounted on the Subaru telescope. Its field of view has covered simultaneously the three brighter lensed images A, B, and C, which are known to exhibit anomalous flux ratios in their continuum emission. We have found that the [OIII] line flux ratios among these lensed images are consistent with those predicted by smooth-lens models. The absence of both microlensing and millilensing effects on this [OIII] narrow line region sets important limits on the mass of any substructures along the line of sight, which is expressed as M_E < 10^5 M_solar for the mass inside an Einstein radius. In contrast, the H_beta line emission, which originates from the broad line region, shows an anomaly in the flux ratio between images B and C, i.e., a factor two smaller C/B ratio than predicted by smooth-lens models. The ratio of A/B in the H_beta line is well reproduced. We show that the anomalous C/B ratio for the H_beta line is caused most likely by micro/milli-lensing of image C. This is because other effects, such as the differential dust extinction and/or arrival time difference between images B and C, or the simultaneous lensing of another pair of images A and B, are all unlikely. In addition, we have found that the broad H_beta line of image A shows a slight asymmetry in its profile compared with those in the other images, which suggests the presence of a small microlensing effect on this line emitting region of image A.

We report new follow-up spectroscopy of i-dropout galaxies with an NB921-band depression found in the Subaru Deep Field. The NB921-depressed i-dropout selection method is expected to select galaxies with large equivalent width Ly alpha emission over a wide redshift range, 6.0<z<6.5. Two of four observed targets show a strong emission line with a clear asymmetric profile, identified as Ly alpha emitters at z=6.11 and 6.00. Their rest-frame equivalent widths are 153A and 114A, which are lower limits on the intrinsic equivalent widths. Through our spectroscopic observations (including previous ones) of NB921-depressed i-dropout galaxies, we identified 5 galaxies in total with a rest-frame equivalent width larger than 100A at 6.0<z<6.5 out of 8 photometric candidates, which suggests that the NB921-depressed i-dropout selection method is possibly an efficient way to search for Ly alpha emitters with a large Ly alpha equivalent width, in a wider redshift range than usual narrow-band excess techniques. By combining these findings with our previous observational results, we infer that the fraction of broad-band selected galaxies having a rest-frame equivalent width larger than 100A is significantly higher at z~6 (the cosmic age of ~1 Gyr) than that at z~3 (~2 Gyr), being consistent with the idea that the typical stellar population of galaxies is significantly younger at z~6 than that at z~3. The NB921-depressed i-dropout galaxies may be interesting candidates for hosts of massive, zero-metallicity Population III stars.

We present six new and two previously published high-resolution spectra of the quasar HS 1603+3820 (z_em = 2.542) taken over an interval of 4.2 years (1.2 years in the quasar rest frame). The observations were made with the High-Dispersion Spectrograph on the Subaru telescope and Medium-Resolution Spectrograph on the Hobby-Eberly Telescope. The purpose was to study the narrow absorption lines (NALs). We use time variability as well as coverage fraction analysis to separate intrinsic absorption lines, which are physically related to the quasar, from intervening absorption lines. By fitting models to the line profiles, we derive the parameters of the respective absorbers as a function of time. Only the mini-BAL system at z_abs ~ 2.43 (v_off ~ 9,500 km/s) shows both partial coverage and time variability, although two NAL systems possibly show evidence of partial coverage. We find that all the troughs of the mini-BAL system vary in concert and its total equivalent width variations resemble those of the coverage fraction. However, no other correlations are seen between the variations of different model parameters. Thus, the observed variations cannot be reproduced by a simple change of ionization state nor by motion of a homogeneous parcel of gas across the cylinder of sight. We propose that the observed variations are a result of rapid continuum fluctuations, coupled with coverage fraction fluctuations caused by a clumpy screen of variable optical depth located between the continuum source and the mini-BAL gas. An alternative explanation is that the observed partial coverage signature is the result of scattering of continuum photons around the absorber, thus the equivalent width of the mini-BAL can vary as the intensity of the scattered continuum changes.

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.

In the deep narrow band image of the Coma Cluster taken with Suprime-Cam of the Subaru telescope, we found an extremely long and narrow (~ 60 kpc x 2 kpc) Halpha emitting region associated with a poststarburst galaxy (D100). Follow up spectroscopy shows that the region has the same redshift as D100. The surface brightness of the region is typically 25 mag(AB) arcsec^-2 in Halpha, which corresponds to 0.5--4x10^-17 erg s^-1 cm^-2 arcsec^-2. We set two possible explanations for the origin of the region; gas stripped off from a merged dwarf, or gas stripped off from D100 by ram pressure. Either scenario has difficulty to fully explain all the observed characteristics of the region.

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.

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 and spectroscopic confirmation of three gravitationally-lensed images of a galaxy at z=3.9 in the background of a distant, rich cluster of galaxies at z=0.83, on the basis of observations with Faint Object Camera And Spectrograph (FOCAS) on the Subaru telescope. We construct a simple lens model of the cluster mass distribution based on Jee et al.'s weak lensing mass estimates from deep, high-resolution images by Advanced Camera for Surveys (ACS) on the Hubble Space Telescope. The lens model can account simultaneously for the observed image configuration and the flux ratio of the closer pair located close to the critical curve. The parities of the three images are also consistent with the lensing hypothesis. Since this galaxy is apparently bright (i'_{AB} \sim 23.7) for its redshift due to the magnification, it serves as a good high redshift target on which we can make extensive and detailed studies based on multi-wavelength imaging and spectroscopy.

We carried out a Subaru and UKIRT near infrared imaging survey for H-alpha emitting galaxies around two pair quasar systems (Q0301-005/Q0302-003 and Q2343+125/Q2344+125), and a triple quasar system (KP76/KP77/KP78). Narrow band near infrared filters covering the H-alpha emission expected for galaxies at the confirmed C IV absorption redshift toward the quasar systems were used for this survey. These quasar pairs or triplet are separated at most by 17 arcmins (~5 h^-1 Mpc in proper distance) from each other on the sky, and have common C IV absorption lines at almost identical redshifts at z=2.24-2.43, which suggests there could be a Mpc-scale absorbing systems such as a cluster, or a group, of galaxies that cover all the line-of-sights to the pair/triple quasars. Using narrow-band deep images, we detected five candidates for H-alpha emitting galaxies around two of the six fields, Q2343+125 and Q2344+125, whose apparent star formation rates are, extremely high, 20-466 M_solar/year. However, all or most of them are not likely to be galaxies at the absorption redshift but galaxies at lower redshift, because of their extreme brightness. In the fields of the other quasars, we detected no star-forming galaxies, nor did we find any number excess of galaxy counts around them. This no-detection results could be because the luminosities and star formation rates of galaxies are lower than the detection limits of our observations (K' > 21 and SFR < 1.8-240 h^-2 M_solar/year). They could be located outside of the observed fields by chance. Otherwise, most C IV absorption lines could be ascribed not to cluster of galaxies, but to isolated star forming pockets far from bright galaxies and could be analogous objects to weak Mg II absorbers.

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.

In order to search for high-redshift galaxies beyond $z = 6.6$ in the Subaru Deep Field, we have investigated NB921-dropout galaxies where NB921 is the narrowband filter centered at 919.6 nm with FWHM of 13.2 nm for the Suprime-Cam on the Subaru Telescope. There are no secure NB921-dropout candidates brighter than $z^\prime = 25.5$. Based on this result, we discuss the UV luminosity function of star-forming galaxies at $z > 6.6$.

We present Subaru spectroscopy of the early-type lensing galaxies in the lens systems HST14113+5211 (z=0.464) and B2045+265 (z=0.868), being aimed at measuring the velocity dispersion of the lensing galaxies as an important component to their mass distributions and internal dynamics. For HST14113+5211 we have obtained 174+/-20km/s (1-sigma) inside an aperture of 0.9", and for B2045+265 we have obtained 213+/-23km/s inside an aperture of 1.2". To extract the significance of these information on the mass distributions and stellar dynamics of the lensing galaxies, we construct two alternative models for both the mass and velocity distributions. It is found that the mass-to-light ratios as derived from a joint lensing and stellar dynamics analysis are virtually consistent with those obtained from the use of the Fundamental Plane. Also, for HST14113+5211, we find that the total mass distribution is well reproduced by a power-law density profile with an index of 2, thereby suggesting that a singular isothermal model is a good fit to the lensing galaxy. In contrast, the corresponding slope for B2045+265 is flatter than isothermal, suggesting that additional contributions to lensing mass from surrounding structures (possibly the group of galaxies) provide the observed angular separations of multiple images.

In order to search for high-redshift galaxies around $z \sim 6$ in the Subaru Deep Field, we have investigated NB816-dropout galaxies where NB816 is the narrowband filter centered at 815nm with FWHM of 12.5 nm for the Suprime-Cam on the Subaru Telescope. Since the NB816 imaging is so deep we can detect 10 well-defined NB816-dropout galaxies that we identify to lie at $z$ = 5.8 -- 6.5. We discuss their observational properties.

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 simulations of quasar colors, magnitudes, and numbers at redshifts 5<z<10 based on our discovery of ten new high-redshift quasars and the cloning of lower redshift Sloan Digital Sky Survey (SDSS) quasars. The ten quasars have redshifts ranging from z=4.7 to z=5.3 and i-magnitudes of 20.21 to 20.94. The natural diversity of spectral features in the cloned sample allows more realistic simulation of the quasar locus width than previously possible with synthetic template spectra. Colors are generated for the z>6 epoch taking advantage of the new UKIDSS near-infrared filter set, and we examine the redshift intervals of maximum productivity, discussing color selection and survey depth issues. On the basis of the SDSS sample, we find that the surface density of z>4.7 quasars increases by a factor of 3X by extending 0.7 i-magnitudes deeper than the SDSS spectroscopic survey limit of i=20.2 -- correspondingly we predict a total of ~400 faint quasars in the SDSS main area that have redshift z>4.7 and magnitudes $<20.9.

We present mid-infrared imaging at 11.7 mu m for the quadruple lens systems, PG1115+080 and B1422+231, using the cooled mid-infrared camera and spectrometer (COMICS) attached on the Subaru telescope. These lensed QSOs are characterized by their anomalous optical and radio flux ratios as obtained for (A1, A2) images of PG1115+080 and (A, B, C) images of B1422+231, respectively, i.e., such flux ratios are hardly reproduced by lens models with smooth mass distribution. Our mid-infrared observations for these images have revealed that the mid-infrared flux ratio A2/A1 of PG1115+080 is virtually consistent with smooth lens models (but inconsistent with the optical flux ratio), whereas for B1422+231, the mid-infrared flux ratios among (A, B, C) are in good agreement with the radio flux ratios. We also identify a clear infrared bump in the spectral energy distributions of these QSOs, thereby indicating that the observed mid-infrared fluxes originate from a hot dust torus around a QSO nucleus. Based on the size estimate of the dust torus, we place limits on the mass of a substructure in these lens systems, causing the anomalous optical or radio flux ratios. For PG1115+080, the mass of a substructure inside an Einstein radius, M_E, is < 16 Msun, corresponding to either a star or a low-mass CDM subhalo having the mass of M_{100}^{SIS} < 2.2 * 10^4 Msun inside radius of 100 pc if modeled as a singular isothermal sphere (SIS). For B1422+231, we obtain M_E > 209 Msun, indicating that a CDM subhalo is more likely, having the mass of M_{100}^{SIS} > 7.4 * 10^4 Msun

[OIII] lambda 5007AA emission lines of 16 intracluster planetary nebulae candidates in the Coma cluster were detected with a Multi-Slit Imaging Spectroscopy (MSIS) technique using FOCAS on the Subaru telescope. The identification of these faint emission sources as PNe is supported by (i) their point-like flux distribution in both space and wavelength, with tight limits on the continuum flux; (ii) the identification of the second [OIII] lambda 4959 line in the only object at high enough velocity that this line too falls into the filter bandpass; (iii) emission line fluxes consistent with PNe at 100 Mpc distance, in the range 2.8 x 10^{-19} - 1.7 x 10^{-18} erg/s/cm^2; and (iv) a narrow velocity distribution approximately centered on the systemic velocity of the Coma cluster. Comparing with the velocities of galaxies in our field, we conclude that the great majority of these candidates would be intracluster PNe, free floating in the Coma cluster core. Their velocity dispersion is ~760 km/s, and their mean velocity is lower than that of the galaxies. The velocity distribution suggests that the intracluster stellar population has different dynamics from the galaxies in the Coma cluster core.

We present new results of a deep optical imaging survey using a narrowband filter ($NB921$) centered at $\lambda =$ 9196 \AA ~ together with $B$, $V$, $R$, $i^\prime$, and $z^\prime$ broadband filters in the sky area of the Subaru Deep Field which has been promoted as one of legacy programs of the 8.2m Subaru Telescope. We obtained a photometric sample of 58 Ly$\alpha$ emitter candidates at $z \approx$ 6.5 -- 6.6 among $\sim 180$ strong $NB921$-excess ($z^\prime - NB921 > 1.0$) objects together with a color criterion of $i^\prime - z^\prime > 1.3$. We then obtained optical spectra of 20 objects in our $NB921$-excess sample and identified at least nine Ly$\alpha$ emitters at $z \sim 6.5$ -- 6.6 including the two emitters reported by Kodaira et al. (2003). Since our Ly$\alpha$ emitter candidates are free from strong amplification of gravitational lensing, we are able to discuss their observational properties from a statistical point of view. Based on these new results, we obtain a lower limit of the star formation rate density of $\rho_{\rm SFR} \simeq 5.5 \times 10^{-4}$ $h_{0.7}$ $M_\odot$ yr$^{-1}$ Mpc$^{-3}$ at $z \approx 6.6$, being consistent with our previous estimate. We discuss the nature of star-formation activity in galaxies beyond $z=6$.

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 have carried out near-infrared JHK spectroscopy of a gravitationally lensed submillimeter galaxy SMM J14011+0252 at z=2.565, using OHS and CISCO on the Subaru telescope. This object consists of two optical components, J1 and J2, which are lensed by the cluster Abell 1835. J1 suffers additional strong lensing by a foreground galaxy at z=0.25 in the cluster. The rest-optical H-alpha, H-beta, and [O II]3727 lines are detected in both J1 and J2, and [N II]6548,6583 lines are also detected in J1. A diagnosis of emission-line ratios shows that the excitation source of J1 is stellar origin, consistent with previous X-ray observations. The continua of J1 and J2 show breaks at rest 4000A, indicating relatively young age. Combined with optical photometry, we have carried out model spectrum fitting of J2 and find that it is a very young (~50 Myr) galaxy of rather small mass (~10e8 M_sol) which suffers some amount of dust extinction. A new gravitational lensing model is constructed to assess both magnification factor and contamination from the lensing galaxy of the component J1, using HST-F702W image. We have found that J1 suffers strong lensing with magnification of ~30, and its stellar mass is estimated to be < 10e9 M_sol. These results suggest that SMM J14011+0252 is a major merger system at high redshift that undergoes intense star formation, but not a formation site of a giant elliptical. Still having plenty of gas, it will transform most of the gas into stars and will evolve into a galaxy of < 10e10 M_sol. Therefore, this system is possibly an ancestor of a less massive galaxy such as a mid-sized elliptical or a spiral at the present.

We present VLT FORS1 and FORS2 spectra of 39 candidate high-redshift supernovae that were discovered as part of a cosmological study using Type Ia supernovae (SNe Ia) over a wide range of redshifts. From the spectra alone, 20 candidates are spectrally classified as SNe Ia with redshifts ranging from z=0.212 to z=1.181. Of the remaining 19 candidates, 1 might be a Type II supernova and 11 exhibit broad supernova-like spectral features and/or have supernova-like light curves. The candidates were discovered in 8 separate ground-based searches. In those searches in which SNe Ia at z ~ 0.5 were targeted, over 80% of the observed candidates were spectrally classified as SNe Ia. In those searches in which SNe Ia with z > 1 were targeted, 4 candidates with z > 1 were spectrally classified as SNe Ia and later followed with ground and space based observatories. We present the spectra of all candidates, including those that could not be spectrally classified as supernova.

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

Using H I absorption alone, we attempt to separate H I absorption lines in quasar spectra into two categories; HDLs (Higher Density Lines) and LDLs (Lower Density Lines), and we discuss the difference in their physical properties. We deblend and fit all H I lines with Voigt profiles, and make an unbiased sample of H I lines covering a wide column density range (12 < log N_HI < 19 cm^-2). To reduce the influence of line blending, we simultaneously fit several Lyman series lines. As a result of a two-point correlation analysis, we found that higher column density H I lines are clustering at $\Delta$ v < 200 km/s, while lower ones at $\Delta$ v < 100 km/s. We define HDLs as H I lines with 15 < log N_HI < 19 cm^-2 and all H I lines within +/- 200 km/s of a line with log N_HI > 15 cm^-2, and LDLs as others with 12 < log N_HI < 15 cm^-2. We found that the HDLs have smaller minimum b-values for a given column density than the LDLs. This difference is successfully reproduced by our Hydrodynamic simulation. The LDLs seem to be cool or shock-heated diffuse IGM gas, while the HDLs are likely to be cooler dense gas near to galaxies.

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.

A main topic at this meeting is how galaxies are affected when they enter for the first time the cluster environment from the outskirts. Most of the times we are forced to infer the environmental effects indirectly, relying on systematic variations of galaxy properties with environment, but there aren't many examples of direct observations able to unveil ongoing transformations taking place, and the corresponding mechanism producing it. We present a case in which it is possible to identify the cluster environment, and in particular the intracluster medium and the recent infall history of galaxies onto the cluster, as the cause for a recent, abrupt change in the evolutionary history of galaxies.

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 present an investigation of the environment (\approx 600 kpc radius) of a pair of luminous z=4.25 quasars, SDSS J1439-0034 A,B, separated by 33''. An analysis of high-quality Subaru spectra of the quasars suggests that this configuration is indeed a physical pair and not a gravitational lens; the redshifts are slightly different (although marginally consistent with being the same), and the two spectra have strikingly different features. We search for bright galaxies (L\gsim 0.4 L^*) having similar redshifts using the V dropout technique and semi-narrow band imaging looking for strong Ly alpha emission. We find no enhancement in the galaxy density around the quasar pair; its environment differs very little from a general field, with the upper limit of the density enhancement being about 3 at a 90% confidence. We infer that bright quasars happened to appear in two normal galaxies in a general field.

A satisfactory understanding of the origin of the dependence of galaxy properties on their environment has remained, so far, out of reach. In the light of numerous observational results and substantial theoretical progress obtained for clusters of galaxies in the last years, a primary goal is to understand how the star formation activity depends on cluster substructure, i.e. on the merging/accretion history of a cluster. In this contribution we present a case in which it is possible to identify the cluster environment, and in particular the intracluster medium and the recent infall history of galaxies onto the cluster, as the cause for an abrupt change in the star formation histories of a subset of galaxies in the Coma cluster.

A surprising result of our recent spectroscopic survey of galaxies in the Coma cluster has been the discovery of a possible bimodal distribution in the metallicities of faint galaxies at $M_B>-17$. We identified a group of dwarfs with luminosity-weighted metallicities around solar and a group with [M/H] around -1.5. A metallicity bimodality among galaxies of similar luminosities is unexpected and suggests that faint cluster galaxies could be an heterogeneous population that formed through more than one evolutionary path, possibly as a consequence of the cluster environment.

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.

The spectroscopic properties of galaxies in the Coma cluster are compared with those of galaxies in rich clusters at $z \sim 0.5$, to investigate the evolution of the star formation history in clusters. Luminous galaxies with $M_V \leq -20$ and post-starburst/post-starforming (k+a) spectra which constitute a significant fraction of galaxies in distant cluster samples are absent in Coma, where spectacular cases of k+a spectra are found instead at $M_V>-18.5$ and represent a significant proportion of the cluster dwarf galaxy population. A simple inspection of their positions on the sky indicates that this type of galaxy does not show a preferential location within the cluster, but the bluest and strongest-lined group of k+a's lies in projection towards the central 1.4 Mpc of Coma and have radial velocities significantly higher than the cluster mean. We find a striking correlation between the positions of these young and strong post-starburst galaxies and substructure in the hot intracluster medium (ICM) identified from {\it XMM-Newton} data, with these galaxies lying close to the edges of two infalling substructures. This result strongly suggests that the interaction with the dense ICM could be responsible for the quenching of the star formation (thus creating the k+a spectrum), and possibly, for any previous starburst. The evolution with redshift of the luminosity distribution of k+a galaxies can be explained by a ``downsizing effect'', with the maximum luminosity/mass of actively star-forming galaxies infalling onto clusters decreasing at lower redshift. We discuss the possible physical origin of this downsizing effect and the implications of our results for current scenarios of environmental effects on the star formation in galaxies.

We performed deep K'-band imaging observations of 2 massive clusters, MS 0451.6-0305 at z = 0.55 and MS 0440.5+0204 at z = 0.19, for searching counterparts of the faint sub-mm sources behind these clusters, which would provide one of the deepest extremely red object(ERO) samples. Comparing our near-infrared images with optical images taken by the Hubble Space Telescope and by the Subaru Telescope, we identified 13 EROs in these fields. The sky distributions of EROs are consistent with the previous results, that there is a sign of strong clustering among detected EROs. Also, the surface density with corrected lensing amplification factors in both clusters are in good agreement with that derived from previous surveys. We found 7 EROs and 3 additional very red objects in a small area (\sim 0.6 arcmin^2) of the MS 0451.6-0305 field around an extended SCUBA source. Many of their optical and near-infrared colors are consistent with dusty star-forming galaxies at high redshifts(z \sim 1.0-4.0), and they may be constituting a cluster of dusty starburst galaxies and/or lensed star-forming galaxies at high redshift. Their red J-K' colors and faint optical magnitudes suggest they are relatively old massive stellar systems with ages(>300 Mega years) suffering from dust obscuration. We also found a surface-density enhancement of EROs around the SCUBA source in the MS 0440.5+0204 field.

Spectroscopic and spectropolarimetric observations of SN 2003dh/GRB 030329 obtained in 2003 May using the Subaru 8.2 m telescope are presented. The properties of the SN are investigated through a comparison with spectra of the Type Ic hypernovae SNe 1997ef and 1998bw. (Hypernovae being a tentatively defined class of SNe with very broad absorption features: these features suggest a large velocity of the ejected material and possibly a large explosion kinetic energy.) Comparison with spectra of other hypernovae shows that the spectrum of SN 2003dh obtained on 2003 May 8 and 9, i.e., 34-35 rest-frame days after the GRB (for z=0.1685), are similar to those of SN 1997ef obtained ~34-42 days after the fiducial time of explosion of that SN. The match with SN 1998bw spectra is not as good (at rest 7300-8000 A, but again spectra obtained ~33-43 days after GRB 980425 are preferred. This indicates that the SN may have intermediate properties between SNe 1997ef and 1998bw. Based on the analogy with the other hypernovae, the time of explosion of SN 2003dh is then constrained to be between -8 and +2 days of the GRB. The Si and O P-Cygni lines of SN 2003dh seem comparable to those of SN 1997ef, which suggests that the ejected mass in SN 2003dh may match that in SN 1997ef. Polarization was marginally detected at optical wavelengths. This is consistent with measurements of the late afterglow, implying that it mostly originated in the interstellar medium of the host galaxy.

Several high-z (z > 5.7) quasars have been found in the course of Sloan Digital Sky Survey. The presence of such very high-z quasars is expected to give constraints on early structure formation. On one hand, it is suggested that these most luminous objects at high redshift are biased toward the highly magnified objects by gravitational lensing. To clarify the effect of gravitational lensing on the high-z quasars, we began the imaging survey of intervening lensing galaxies. Indeed our previous optical image showed that SDSSp J104433.04+012502.2 at z=5.74 is gravitationally magnified by a factor 2. In this paper, we report our new optical imaging of other two high-z quasars, SDSSp J103027.10+052455.0 at z=6.28 and SDSSp J130608.26+035626.3 at z=5.99. Since we find neither intervening galaxy nor counter image with i^{\prime} < 25.4-25.8 around each quasar, we conclude that they are not strongly magnified regardless that a lens galaxy is dusty.

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 have performed a deep optical imaging survey using a narrowband filter ($NB921$) centered at $\lambda =$ 9196 \AA ~ together with $i^\prime$ and $z^\prime$ broadband filters covering an 814 arcmin$^2$ area of the Subaru Deep Field. We obtained a sample of 73 strong $NB921$-excess objects based on the following two color criteria; $z^\prime - NB921 > 1$ and $i^\prime - z^\prime > 1.3$. We then obtained optical spectroscopy of nine objects in our $NB921$-excess sample, and identified at least two Ly$\alpha$ emitters at$z=6.541 \pm 0.002$ and $z=6.578 \pm 0.002$, each of which shows the characteristic sharp cutoff together with the continuum depression at wavelengths shortward of the line peak. The latter object is more distant than HCM-6A at $z=6.56$ and thus this is the most distant known object found so far. These new data allow us to estimate the first meaningful lower limit of the star formation rate density beyond redshift 6; $\rho_{\rm SFR} \sim 5.2 \times 10^{-4} M_\odot$ yr$^{-1}$ Mpc$^{-3}$. Since it is expected that the actual density is higher by a factor of several than this value, our new observation reveals that a moderately high level of star formation activity already occurred at $z \sim$ 6.6.

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.

A large spectroscopic survey is constructed of galaxies in the Coma cluster. The survey covers a wide area (1 deg$^2$) to deep magnitudes (R ~ 19.5), covering both the core (high density) and outskirts (intermediate to low density) of the cluster. The spectroscopic sample consists of a total of 1191 galaxies, of which, 760 galaxies are confirmed members of the Coma cluster. A statistical technique is developed to correct the spectroscopic sample for incompletness. The corrected sample is then used to construct R-band luminosity function (LF) spanning a range of 7 magnitudes (-23 < M_R < -16) both at the core and outskirts of the cluster. Dependence of the LF on local environment in Coma is explored. The LFs are found to be the same, within the errors, between the inner and outer regions and close to those from recent measurements for field galaxies. The total B-band LF for the Coma cluster, fitted to a Schecter form is also measured and shows a dip at M(B) = -18 mag., in agreement with previous studies. The implications of this feature are discussed. The LF is studied in B-R color intervals and shows a steep faint-end slope for red (B-R > 1.35) galaxies, both at the core and outskirts of the cluster. This population of low luminosity red galaxies has a higher surface density than the blue (B-R < 1.35) star-forming population, dominating the faint-end of the Coma cluster LF. It is found that relative number of high surface brightness galaxies is larger at the cluster core, implying destruction of low surface brightness galaxies in dense core environment.

We present a high resolution spectrum of the quasar, HS1603+3820 (z_em=2.542), observed with the High Dispersion Spectrograph (HDS) on Subaru Telescope. This quasar, first discovered in the Hamburg/CfA Quasar Survey, has 11 C IV lines at 1.96 < z_abs < 2.55. Our spectrum covers 8 of the 11 C IV lines at z_abs > 2.29 and resolves some of them into multiple narrow components with b < 25 km/s because of the high spectral resolution R=45000, while other lines show broad profiles (b > 65 km/s). We use three properties of C IV lines, specifically, time variability, covering factor, and absorption line profile, to classify them into quasar intrinsic absorption lines (QIALs) and spatially intervening absorption lines (SIALs). The C IV lines at 2.42 < z_abs < 2.45 are classified as QIALs in spite of their large velocity shifts from the quasar. Perhaps they are produced by gas clouds ejected from the quasar with the velocity of v_ej = 8000 km/s -- 10000 km/s. On the other hand, three C IV lines at 2.48 < z_abs < 2.55 are classified as SIALs, which suggests there exist intervening absorbers near the quasar. We, however, cannot rule out QIALs for the two lines at z_abs ~ 2.54 and 2.55, because their velocity shifts, 430 km/s blueward and 950 km/s redward of the quasar, are very small. The C IV line at z_abs ~ 2.48 consists of many narrow components, and has also corresponding low-ionization metal lines (Al II, Si II, and Fe II). The velocity distribution of these low-ionization ions is concentrated at the center of the system compared to that of high-ionization C IV ion. Therefore we ascribe this system of absorption lines to an intervening galaxy.

We present spectropolarimetry of the Type Ic supernova SN 2002ap and give a preliminary analysis: the data were taken at two epochs, close to and one month later than the visual maximum (2002 February 8). In addition we present June 9 spectropolarimetry without analysis. The data show the development of linear polarization. Distinct polarization profiles were seen only in the O I \lambda 7773 multiplet/Ca II IR triplet absorption trough at maximum light and in the Ca II IR triplet absorption trough a month later, with the latter showing a peak polarization as high as ~2 %. The intrinsic polarization shows three clear position angles: 80 degs for the February continuum, 120 degs for the February line feature, and 150 degs for the March data. We conclude that there are multiple asymmetric components in the ejecta. We suggest that the supernova has a bulk asymmetry with an axial ratio projected on the sky that is different from 1 by of order 10 %. Furthermore, we suggest very speculatively that a high velocity ejecta component moving faster than ~0.115c (e.g., a jet) contributes to polarization in the February epoch.

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)

We present new optical images ($B$, $V$, and H$\alpha$) of the archetypical starburst/superwind galaxy M82 obtained with the 8.2 m Subaru Telescope to reveal new detailed structures of the superwind-driven nebula and the high-latitude dark lanes. The emission-line nebula is decomposed into (1) a ridge-dominated component comprising numerous filament/loop sub-structures whose overall morphology appears as a pair of narrow cylinders, and (2) a diffuse component extended over much wider opening angle from the nucleus. We suggest that these two components have different origins. The ridge-dominated component appears as a pair of cylinders rather than a pair of cones. Since this morphological property is similar to that of hot plasma probed by soft X-ray, this component seems to surround the hot plasma. On the other hand, the diffuse component may arise from dust grains which scatter stellar light from the galaxy. Since inner region of this component is seen over the prominent ^^ ^^ X"-shaped dark lanes streaming out from the nuclear region and they can be reproduced as a conical distribution of dust grains, there seems to be a dusty cold outflow as well as the hot one probed by soft X-ray and shock-excited optical emission lines. If this is the case, the presence of such high-latitude dust grains implies that neutral gaseous matter is also blown out during the course of the superwind activity.

During the course of our deep optical imaging survey for Ly alpha emitters at z approximately 5.7 in the field around the z=5.74 quasar SDSSp J104433.04-012502.2, we have found a candidate strong emission-line source. Follow-up optical spectroscopy shows that the emission line profile of this object is asymmetric, showing excess red-wing emission. These properties are consistent with an identification of Ly alpha emission at a redshift of z=5.687 +/- 0.002. The observed broad line width, Delta V_{FWHM} ~= 340 km s^{-1} and excess red-wing emission also suggest that this object hosts a galactic superwind.

Photometric and spectroscopic data of the energetic Type Ic supernova (SN) 2002ap are presented, and the properties of the SN are investigated through models of its spectral evolution and its light curve. The SN is spectroscopically similar to the "hypernova" SN 1997ef. However, its kinetic energy [$\sim (4-10) \times 10^{51}$ erg] and the mass ejected (2.5-5 $M_{\odot}$) are smaller, resulting in a faster-evolving light curve. The SN synthesized $\sim 0.07 M_{\odot}$ of $^{56}$Ni, and its peak luminosity was similar to that of normal SNe. Brightness alone should not be used to define a hypernova, whose defining character, namely very broad spectral features, is the result of a high kinetic energy. The likely main-sequence mass of the progenitor star was 20-25 $M_{\odot}$, which is also lower than that of both hypernovae SNe 1997ef and 1998bw. SN 2002ap appears to lie at the low-energy and low-mass end of the hypernova sequence as it is known so far. Observations of the nebular spectrum, which is expected to dominate by summer 2002, are necessary to confirm these values.

We present new optical spectroscopy of the lens elliptical galax in the ``Einstein Cross'' lens system HST 14176+5226, using the Faint Object Camera and Spectrograph (FOCAS) of the Subaru t Our spectroscopic observations are aimed at measuring the stella dispersion of the lens galaxy, located at high redshift of z_L= as an important component to lens models. We have measured this 230 +- 14 km s^{-1} (1 sigma) inside 0.35 effective radi based on the comparison between the observed galaxy spectrum and templates of three G-K giants by means of the Fourier cross-corr To extract the significance of this information on the geometry universe which also affects the lensing of the background image, to fit three different lens models to the available data of the Provided that the lens galaxy has the structural and dynamical p (i.e., its radial density profile, core radius, and velocity ani similar to those of local elliptical galaxies, we calculate the function for the simultaneous reproduction of both the observed and newly measured velocity dispersion of the lens. Although the interval depends rather sensitively on the adopted lens models o parameters, our experiments suggest the larger likelihood for a cosmological constant, Omega_Lambda: formal 1 sigma lower Omega_Lambda in the flat universe ranges 0.73 to 0.97, where lower limit is basically unavailable. This method for determinin model is thus dependent on lens models but is insensitive to oth ambiguities, such as the dust absorption or the evolutionary eff galaxies. Exploring spectroscopic observations of more lens gala redshift may minimize the model uncertainties and thus place a m constraint on Omega_Lambda.

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.

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

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 present new optical spectroscopy of the high-redshift powerful radio galaxy MRC 0406$-$244 at redshift of 2.429. We find that the two extensions toward NW and SE probed in the rest-frame ultraviolet image are heated mainly by the nonthermal continuum of the active galactic nucleus. However, each extension shows a shell-like morphology, suggesting that they are a pair of superbubbles induced by the superwind activity rather than by the interaction between the radio jet and the ambient gas clouds. If this is the case, the intense starburst responsible for the formation of superbubbles could occur $\sim 1 \times 10^9$ yr ago. On the other hand, the age of the radio jets may be of the order of $\sim 10^6$ yr, being much shorter than the starburst age. Therefore, the two events, i.e., the starburst and the radio-jet activities, are independent phenomena. However, their directions of the expanding motions could be governed by the rotational motion of the gaseous component in the host galaxy. This idea appears to explain the alignment effect of MRC 0406$-$244.

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.

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 construct infrared star count models of the Galaxy applicable at faint magnitudes. In addition to dwarfs and subdwarfs, we also take into account L dwarfs, T dwarfs, and white dwarfs of both disk and halo. As a test case of the models, we analyze the infrared imaging data at J and K' obtained during the Subaru Deep Field Survey to study stellar objects. Infrared star count models we have obtained predict that the Next Generation Space Telescope will primarily see T dwarfs, M subdwarfs, and old halo white dwarfs at faint magnitudes.

We present results of optical identification of the X-ray sources detected in the ASCA Large Sky Survey. Optical spectroscopic observations were done for 34 X-ray sources which were detected with the SIS in the 2-7 keV band above 3.5 sigma. The sources are identified with 30 AGNs, 2 clusters of galaxies, and 1 galactic star. Only 1 source is still unidentified. The flux limit of the sample corresponds to 1 x 10^{-13} erg s^{-1} cm^{-2} in the 2-10 keV band. Based on the sample, the paper discusses optical and X-ray spectral properties of the AGNs, contribution of the sources to the Cosmic X-ray Background, and redshift and luminosity distributions of the AGNs. An interesting result is that the redshift distribution of the AGNs suggests a deficiency of high-redshift (0.5 < z < 2) and luminous (Lx(2-10 keV) > 10^{44} erg s^{-1}) absorbed narrow-line AGNs (so called type 2 QSOs).

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 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 report the optical identification of the hardest X-ray source (AX J131501+3141) detected in an unbiased wide-area survey in the 0.5--10 keV band, the ASCA Large Sky Survey. The X-ray spectrum of the source is very hard and is well reproduced by a power-law component (Gamma = 1.5^+0.7_-0.6) with N_H = 6^+4_-2 *10^22 cm^-2 (Sakano et al. 1998). We have found a galaxy with R=15.62 mag near the center of the error circle for the X-ray source. The optical spectrum of the galaxy shows only narrow emission lines whose ratios correspond to those of a type 2 Seyfert galaxy at z = 0.072, implying an absorption-corrected X-ray luminosity of 2*10^43 erg sec^-1 (2--10 keV) and M_B = -20.93 mag. A radio point source is also associated with the center of the galaxy. We thus identify the X-ray source with this galaxy as an obscured AGN. The hidden nature of the nucleus of the galaxy in the optical band is consistent with the X-ray spectrum. These results support the idea that the obscured AGNs/QSOs contribute significantly to the cosmic X-ray background in the hard band at the faint flux level.