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

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

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

We present our new ALMA observations targeting [OIII]88um, [CII]158um, [NII]122um, and dust continuum emission for three Lyman break galaxies at z=6.0293-6.2037 identified in the Subaru/Hyper Suprime-Cam survey. We clearly detect [OIII] and [CII] lines from all of the galaxies at 4.3-11.8sigma levels, and identify multi-band dust continuum emission in two of the three galaxies, allowing us to estimate infrared luminosities and dust temperatures simultaneously. In conjunction with previous ALMA observations for six galaxies at z>6, we confirm that all the nine z=6-9 galaxies have high [OIII]/[CII] ratios of L[OIII]/L[CII]~3-20, ~10 times higher than z~0 galaxies. We also find a positive correlation between the [OIII]/[CII] ratio and the Lya equivalent width (EW) at the ~90% confidence level. We carefully investigate physical origins of the high [OIII]/[CII] ratios at z=6-9 using Cloudy, and find that high density of the interstellar medium, low C/O abundance ratio, and the cosmic microwave background attenuation are responsible to only a part of the z=6-9 galaxies. Instead, the observed high [OIII]/[CII] ratios are explained by 10-100 times higher ionization parameters or low photodissociation region (PDR) covering fractions of 0-10%, both of which are consistent with our [NII] observations. The latter scenario can be reproduced with a density bounded nebula with PDR deficit, which would enhance the Lya, Lyman continuum, and C+ ionizing photons escape from galaxies, consistent with the [OIII]/[CII]-Lya EW correlation we find.

We present an IGM HI tomography map in a survey volume of $16 \times 19 \times 131 \ h^{-3} {\rm comoving \ Mpc}^{3}$ (cMpc$^3$) centered at MAMMOTH-1 nebula and three neighbouring quasars at $z=2.3$. MAMMOTH-1 nebula is an enormous Ly$\alpha$ nebula (ELAN), hosted by a type-II quasar dubbed MAMMOTH1-QSO, that extends over $1\ h^{-1}$ cMpc with not fully clear physical origin. Here we investigate the HI-gas distribution around MAMMOTH1-QSO with the ELAN and three neighbouring type-I quasars, making the IGM HI tomography map with a spatial resolution of $2.6\ h^{-1}$ cMpc. Our HI tomography map is reconstructed with HI Ly$\alpha$ forest absorption of bright background objects at $z=2.4-2.9$: one eBOSS quasar and 16 Keck/LRIS galaxy spectra. We estimate the radial profile of HI flux overdensity for MAMMOTH1-QSO, and find that MAMMOTH1-QSO resides in a volume with significantly weak HI absorption. This suggests that MAMMOTH1-QSO has a proximity zone where quasar illuminates and photo-ionizes the surrounding HI gas and suppresses HI absorption, and that the ELAN is probably a photo-ionized cloud embedded in the cosmic web. The HI radial profile of MAMMOTH1-QSO is very similar to those of three neighbouring type-I quasars at $z=2.3$, which is compatible with the AGN unification model. We compare the distributions of the HI absorption and star-forming galaxies in our survey volume, and identify a spatial offset between density peaks of star-forming galaxies and HI gas. This segregation may suggest anisotropic UV background radiation created by star-forming galaxy density fluctuations.

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

We report on the detection of the [CII] 157.7 $\mu$m emission from the Lyman break galaxy (LBG) MACS0416_Y1 at z = 8.3113, by using the Atacama Large Millimeter/submillimeter Array (ALMA). The luminosity ratio of [OIII] 88 $\mu$m (from previous campaigns) to [CII] is 9.31 $\pm$ 2.6, indicative of hard interstellar radiation fields and/or a low covering fraction of photo-dissociation regions. The emission of [CII] is cospatial to the 850 $\mu$m dust emission (90 $\mu$m rest-frame, from previous campaigns), however the peak [CII] emission does not agree with the peak [OIII] emission, suggesting that the lines originate from different conditions in the interstellar medium. We fail to detect continuum emission at 1.5 mm (160 $\mu$m rest-frame) down to 18 $\mu$Jy (3$\sigma$). This nondetection places a strong limit on the dust spectrum, considering the 137 $\pm$ 26 $\mu$Jy continuum emission at 850 $\mu$m. This suggests an unusually warm dust component (T $>$ 80 K, 90% confidence limit), and/or a steep dust-emissivity index ($\beta_{\rm dust}$ $>$ 2), compared to galaxy-wide dust emission found at lower redshifts (typically T $\sim$ 30 - 50 K, $\beta_{\rm dust}$ $\sim$ 1 - 2). If such temperatures are common, this would reduce the required dust mass and relax the dust production problem at the highest redshifts. We therefore warn against the use of only single-wavelength information to derive physical properties, recommend a more thorough examination of dust temperatures in the early Universe, and stress the need for instrumentation that probes the peak of warm dust in the Epoch of Reionization.

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

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

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

We present results of the cross-correlation Lya intensity mapping with Subaru/Hyper Suprime-Cam (HSC) ultra-deep narrowband images and Lya emitters (LAEs) at $z=5.7$ and $6.6$ in a total area of $4$ deg$^2$. Although overwhelming amount of data quality controls have been performed for the narrowband images and the LAE samples, we further conduct extensive analysis evaluating systematics of large-scale point spread function wings, sky subtractions, and unknown errors on the basis of physically uncorrelated signals and sources found in real HSC images and object catalogs, respectively. Removing the systematics, we carefully calculate cross-correlations between Lya intensity of the narrowband images and the LAEs. We identify very diffuse Lya emission with the $3\sigma$ ($2\sigma$) significance at $>150$ comoving kpc (ckpc) far from the LAEs at $z = 5.7$ ($6.6$), beyond a virial radius of star-forming galaxies with $M_\mathrm{h} \sim 10^{11} M_\odot$. The diffuse Lya emission possibly extends up to $1,000$ ckpc with the surface brightness of $10^{-20}-10^{-19}$ erg s$^{-1}$ cm$^{-2}$ arcsec$^{-2}$. We confirm that the small-scale ($<150$ ckpc) Lya radial profiles of LAEs in our Lya intensity maps are consistent with those obtained by recent MUSE observations (Leclercq et al. 2017). Comparisons with numerical simulations suggest that the large-scale ($\sim 150-1,000$ ckpc) Lya emission are not explained by unresolved faint sources of neighboring galaxies including satellites, but by a combination of Lya photons emitted from the central LAE and other unknown sources, such as a cold-gas stream and galactic outflow. We find no evolution in the Lya radial profiles of our LAEs from $z=5.7$ to $6.6$, where theoretical models predict a flattening of the profile slope made by cosmic reionization, albeit with our moderately large observational errors.

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

We present the Atacama Large Millimeter/submillimeter Array (ALMA) detection of the [O III] 88 $\mu$m line and rest-frame 90 $\mu$m dust continuum emission in a Y-dropout Lyman break galaxy (LBG), MACS0416_Y1, lying behind the Frontier Field cluster MACS J0416.1-2403. This [O III] detection confirms the LBG with a spectroscopic redshift of $z = 8.3118 \pm 0.0003$, making this object one of the furthest galaxies ever identified spectroscopically. The observed 850 $\mu$m flux density of $137 \pm 26$ $\mu$Jy corresponds to a de-lensed total infrared (IR) luminosity of $L_{\rm IR} = (1.7 \pm 0.3) \times 10^{11} L_{\odot}$ if assuming a dust temperature of $T_{\rm dust} = 50$ K and an emissivity index of $\beta = 1.5$, yielding a large dust mass of $4 \times 10^6 M_{\odot}$. The ultraviolet-to-far IR spectral energy distribution modeling where the [O III] emissivity model is incorporated suggests the presence of a young ($\tau_{\rm age} \approx 4$ Myr), star-forming (SFR $\approx 60 M_{\odot}$ yr$^{-1}$), moderately metal-polluted ($Z \approx 0.2 Z_{\odot}$) stellar component with a mass of $M_{\rm star} = 3 \times 10^8 M_{\odot}$. An analytic dust mass evolution model with a single episode of star-formation does not reproduce the metallicity and dust mass in $\tau_{\rm age} \approx 4$ Myr, suggesting a pre-existing evolved stellar component with $M_{\rm star} \sim 3 \times 10^9 M_{\odot}$ and $\tau_{\rm age} \sim 0.3$ Gyr as the origin of the dust mass.

We present the redshift evolution of the radial surface brightness (SB) profile of the rest-frame UV and optical stellar continua for 9119 Lya emitters (LAEs) at z~0-8 and 0-2, respectively. Using Hubble Space Telescope data and the LAE catalogs taken from the literature, we derive the structural quantities of the 9119 LAEs and ~180,000 comparison galaxies of photo-z star-forming galaxies (SFGs) and Lyman break galaxies (LBGs) by the well-tested profile fitting. From 936 well-fitted LAEs, we carefully define the homogeneous sample of LAEs falling in the same ranges of UV-continuum luminosity and Lya equivalent width over z~0-8, and evaluate the redshift evolution. We find that the effective radius r_e distribution is represented by a log-normal function, and that the median Sersic index is almost constant at n~1-1.5 for the LAEs over z~0-7, suggesting that typical LAEs have a stellar-disk morphology. The size-luminosity relation of the LAEs monotonically decreases towards high-z, following size-luminosity relations of SFGs and LBGs. The median r_e values of the LAEs significantly evolve as r_e~(1+z)^-1.37, similar to those of the SFGs and LBGs in the same luminosity range, in contrast with the claims of no evolution made by previous studies whose LAE samples are probably biased to faint sources at low-z. The r_e distribution, star-formation rate surface densities, and stellar-to-halo size ratios of the LAEs are comparable with those of the SFGs and LBGs, indicating that LAEs have stellar components similar to SFGs and LBGs with a Lya emissivity controlled by the non-stellar physics such as geometry, kinematics, and ionization states of the inter-stellar/circum-galactic medium.

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

We present morphological properties of dusty star-forming galaxies at z=1-3 determined with high-resolution (FWHM~0"19) Atacama Large Milllimeter/submilimeter Array (ALMA) 1-mm band maps of our ASAGAO survey covering a 26-arcmin^2 area in GOODS-S. In conjunction with the ALMA archival data, the present sample consists of 42 ALMA sources with a wide rest-frame far-infrared (FIR) luminosity L_FIR range of ~10^11-10^13 Lo. To obtain an average rest-frame FIR profile, we perform individual measurements and careful stacking of the ALMA sources using the uv-visibility method that includes positional-uncertainty and smoothing-effect evaluations through Monte-Carlo simulations. We find that the dusty star-forming galaxies have the average FIR-wavelength Sersic index and effective radius of n_FIR=1.2+/-0.2 and R_e,FIR=1.0-1.3 kpc, respectively, additionally with a point source at the center, indicative of the existence of AGN. The average FIR profile agrees with a morphology of an exponential-disk clearly distinguished from a spheroidal profile (Sersic index of 4). We also examine the rest-frame optical Sersic index n_opt and effective radius R_e,opt with the deep Hubble Space Telescope (HST) images. Interestingly, we obtain n_opt=0.9+/-0.3 (~n_FIR) and R_e,opt=3.2+/-0.6 kpc (>R_e,FIR), suggesting that the FIR-emitting disk is embedded within a larger stellar disk. The rest-frame UV and FIR data of HST and ALMA provide us a radial surface density profile of the total star-formation rate (SFR), where the FIR SFR dominates over the UV SFR at the center. Under the simple assumption of a constant SFR, a compact stellar distribution found in z~1-2 compact quiescent galaxies (cQGs) is well reproduced, while a spheroidal stellar morphology of cQGs (n_opt=4) cannot, suggestive of other important mechanisms such as dynamical dissipation.

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

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

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

We present a multi-wavelength catalog in the Subaru-XMM Deep Field (SXDF) as part of the Spitzer Large Area Survey with Hyper-Suprime-Cam (SPLASH). We include the newly acquired optical data from the Hyper-Suprime Cam Subaru Strategic Program, accompanied by IRAC coverage from the SPLASH survey. All available optical and near-infrared data is homogenized and resampled on a common astrometric reference frame. Source detection is done using a multi-wavelength detection image including the $u$-band to recover the bluest objects. We measure multi-wavelength photometry and compute photometric redshifts as well as stellar masses for $\sim$1.17 million objects over $\sim$4.2 deg$^2$ with $\sim$800,000 objects in the 2.4 deg$^2$ HSC-UltraDeep coverage. Using the available spectroscopic redshifts from various surveys over the range of $0<z<6$, we verify the performance of the photometric redshifts and we find a normalized median absolute deviation of 0.023 and outlier fraction of 3.2%. The SPLASH-SXDF catalog is a valuable, publicly available resource that is perfectly suited for studying galaxies in the early universe and tracing their evolution through cosmic time. The catalog is available for download from https://z.umn.edu/SXDF .

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 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 conduct a systematic study of galactic outflows in star-forming galaxies at $z\sim0$-$2$ based on the absorption lines of optical spectra taken from SDSS DR7, DEEP2 DR4, and Keck Erb et al. We carefully make stacked spectra of homogeneous galaxy samples with similar stellar mass distributions at $z\sim0$-$2$, and perform the multi-component fitting of model absorption lines and stellar continua to the stacked spectra. We obtain the maximum ($v_\rm{max}$) and central ($v_\rm{out}$) outflow velocities, and estimate the mass loading factors ($\eta$), a ratio of the mass outflow rate to the star formation rate (SFR). Investigating the redshift evolution of the outflow velocities measured with the absorption lines whose depths and ionization energies are similar (Na I D and Mg I at $z\sim0$-$1$; Mg II and C II at $z\sim1$-$2$), we identify, for the first time, that the average value of $v_\rm{max}$ ($v_\rm{out}$) significantly increases by 0.05-0.3 dex from $z\sim0$ to $2$ at a given SFR. Moreover, we find that the value of $\eta$ increases from $z\sim0$ to $2$ by $\eta \propto (1 + z)^{1.2\pm0.3}$ at a given halo circular velocity $v_\rm{cir}$ , albeit with a potential systematics caused by model parameter choices. The redshift evolution of $v_\rm{max}$ ($v_\rm{out}$) and $\eta$ is consistent with the galaxy-size evolution and the local velocity-SFR surface density relation, and explained by high-gas fractions in high-redshift massive galaxies, which is supported by recent radio observations. We obtain a scaling relation of $\eta \propto v_\rm{cir}^a$ for $a = -0.2 \pm 1.1$ in our $z\sim0$ galaxies that agrees with the momentum-driven outflow model ($a = -1$) within the uncertainty.

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 the large statistics of the galaxy effective radius in the rest-frame far-infrared (FIR) wavelength Re(FIR) obtained from 1627 deep ALMA 1-mm band maps that become public by 2017 July. Our ALMA sample consists of 1034 sources at z=0-6 that typically have the star-formation rate of ~100-1000 M$_{\odot}$/yr and the stellar mass of ~10$^{10}$-10$^{11.5}$ M$_{\odot}$. We homogeneously derive Re(FIR) and FIR luminosity L_FIR of our ALMA sources via the same uv-visibility method with the exponential disk fitting, carefully evaluating the selection and measurement incompletenesses by realistic Monte-Carlo simulations. We find that there is a positive correlation between Re(FIR) and L_FIR at the >99% significance level. The best-fit power-law function, Re(FIR) $\propto$ L_FIR$^{\alpha}$, provides $\alpha=0.28+/-0.07$, and shows that the average Re(FIR) at a fixed L_FIR decreases toward high redshifts. The best-fit $\alpha$ and the redshift evolution of Re(FIR) are similar to those of the galaxy effective radius in the rest-frame UV (optical) wavelength Re(UV) (Re(Opt.)) revealed by Hubble Space Telescope (HST) studies. We identify that our ALMA sources have significant trends of Re(FIR) <~ Re(UV) and Re(Opt.), suggesting that dusty starbursts take place in a compact region. Moreover, Re(FIR) of our ALMA sources is comparable to Re(Opt.) of quiescent galaxies at z~1-3 as a function of stellar mass, supporting the evolutionary connection between these two galaxy populations. We also investigate rest-frame UV and optical morphologies of our ALMA sources with deep HST images, and find that ~30-40% of our ALMA sources are classified as major mergers. This indicates that dusty starbursts are triggered not only by the major mergers but also the other mechanism(s).

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

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

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

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

We present initial results from the Subaru Strategic Program (SSP) with Hyper Suprime-Cam (HSC) on a comprehensive survey of emission-line galaxies at z<1.5 based on narrowband (NB) imaging. The first Public Data Release (PDR1) provides us with data from two NB filters, specifically NB816 and NB921 over 5.7 deg$^2$ and 16.2 deg$^2$ respectively. The $5 \sigma$ limiting magnitudes are 25.2 (UDeep layer, 1.4 deg$^2$) and 24.8 (Deep layer, 4.3 deg$^2$) mag in NB816, and 25.1 (UDeep, 2.9 deg$^2$) and 24.6--24.8 (Deep, 13.3 deg$^2$) mag in NB921. The wide-field imaging allows us to construct unprecedentedly large samples of 8,054 H$\alpha$ emitters at z ~ 0.25 and 0.40, 8,656 [OIII] emitters at z ~ 0.63 and 0.84, and 16,877 [OII] emitters at z ~ 1.19 and 1.47. We map the cosmic web on scales out to about 50 comoving Mpc that includes galaxy clusters, identified by red sequence galaxies, located at the intersection of filamentary structures of star-forming galaxies. The luminosity functions of emission-line galaxies are measured with precision and consistent with published studies. The wide field coverage of the data enables us to measure the luminosity functions up to brighter luminosities than previous studies. The comparison of the luminosity functions between the different HSC-SSP fields suggests that a survey volume of $>5\times10^5$ Mpc$^3$ is essential to overcome cosmic variance. Since the current data have not reached the full depth expected for the HSC-SSP, the color cut in i-NB816 or z-NB921 induces a bias towards star-forming galaxies with large equivalent widths, primarily seen in the stellar mass functions for the H$\alpha$ emitters at z ~ 0.25--0.40. Even so, the emission-line galaxies clearly cover a wide range of luminosity, stellar mass, and environment, thus demonstrating the usefulness of the NB data from the HSC-SSP to investigate star-forming galaxies at z<1.5.

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

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

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 N/O abundance ratios, ionization parameters {$q_{\rm ion}$}, and oxygen abundance O/H for a total of 41 galaxies (11 individual galaxies and a 30-galaxy stack) including {Ly$\alpha$} emitters and Lyman break galaxies at $z\sim 2$, and investigate galaxy evolution from $z\sim 0$ to $2$ in conjunction with 208,529 local SDSS galaxies and 9 green pea galaxies (GPs). In contrast with most of the recent studies, we obtain the N/O ratio, {$q_{\rm ion}$}, and O/H measurements by direct $T_{\rm e}$ methods with {[O\,\textsc{iii}]} $\lambda$4363 and {O\,\textsc{iii}]} $\lambda$1665 lines. Based on these reliable measurements, we find that there exist $z\sim 2$ galaxies with an excess of N/O falling beyond the local average of N/O--O/H relation, while the majority of the $z\sim 2$ galaxies have the N/O ratios nearly comparable with $z\sim 0$ galaxies in the N/O--stellar mass relation. Our galaxies place the upper limit of N/O ratio $\log ({\rm N/O})\le -1.26$ on average, suggesting that the N/O ratio evolves, if at all, by $<0.17$ dex. Exploiting our reliable measurements free from the N/O--{$q_{\rm ion}$}--O/H degeneracies, we identify, for the first time, that $z\sim 2$ galaxies with offsets in the BPT diagram indicate either 1) only an N/O excess, 2) only a {$q_{\rm ion}$} excess, or 3) both N/O and {$q_{\rm ion}$} excesses. We argue that the BPT offsets at $z\sim 2$ are not made by one of 1)--3) galaxy populations alone, but the composite of 1)--3) populations. We confirm that these 1)--3) populations also exist at $z\sim 0$, such as GPs and SDSS low-mass and high-SFR galaxies (LMHSs).

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 present spatial correlations of galaxies and IGM HI in the COSMOS/UltraVISTA 1.62 deg$^2$ field. Our data consist of 13,415 photo-$z$ galaxies at $z\sim2-3$ with $K_s<23.4$ and the Ly$\alpha$ forest absorptions in the background quasar spectra selected from SDSS data with no signature of damped Ly$\alpha$ system contamination. We estimate a galaxy overdensity $\delta_{gal}$ in an impact parameter of 2.5 pMpc, and calculate the Ly$\alpha$ forest fluctuations $\delta_{\langle F\rangle}$ whose negative values correspond to the strong Ly$\alpha$ forest absorptions. We identify weak evidence of an anti-correlation between $\delta_{gal}$ and $\delta_{\langle F\rangle}$ with a Spearman's rank correlation coefficient of $-0.39$ suggesting that the galaxy overdensities and the Ly$\alpha$ forest absorptions positively correlate in space at the $\sim90\%$ confidence level. This positive correlation indicates that high-$z$ galaxies exist around an excess of HI gas in the Ly$\alpha$ forest. We find four cosmic volumes, dubbed $A_{obs}$-$D_{obs}$, that have extremely large (small) values of $\delta_{gal} \simeq0.8$ ($-1$) and $\delta_{\langle F\rangle}$ $\simeq0.1$ ($-0.4$), three out of which, $B_{obs}$-$D_{obs}$, significantly depart from the correlation, and weaken the correlation signal. We perform cosmological hydrodynamical simulations, and compare with our observational results. Our simulations reproduce the correlation, agreeing with the observational results. Moreover, our simulations have model counterparts of $A_{obs}$-$D_{obs}$, and suggest that the observations pinpoint, by chance, a galaxy overdensity like a proto-cluster, gas filaments lying on the sightline, a large void, and orthogonal low-density filaments. Our simulations indicate that the significant departures of $B_{obs}$-$D_{obs}$ are produced by the filamentary large-scale structures and the observation sightline effects.

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

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

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

We investigate evolution of clumpy galaxies with the Hubble Space Telescope (HST) samples of ~17,000 photo-z and Lyman break galaxies at z~0-8. We detect clumpy galaxies with off-center clumps in a self-consistent algorithm that is well tested with previous study results, and measure the number fraction of clumpy galaxies at the rest-frame UV, f_clumpy^UV. We identify an evolutionary trend of f_clumpy^UV over z~0-8 for the first time: f_clumpy^UV increases from z~8 to z~1-3 and subsequently decreases from z~1 to z~0, which follows the trend of Madau-Lilly plot. A low average Sersic index of n~1 is found in the underlining components of our clumpy galaxies at z~0-2, indicating that typical clumpy galaxies have disk-like surface brightness profiles. Our f_clumpy^UV values correlate with physical quantities related to star formation activities for star-forming galaxies at z~0-7. We find that clump colors tend to be red at a small galactocentric distance for massive galaxies with log(M_*/M_sun)>~11. All of these results are consistent with a picture that a majority of clumps form in the violent disk instability and migrate into the galactic centers.

We present statistics of 133 faint 1.2-mm continuum sources detected in about 120 deep ALMA pointing data that include all the archival deep data available by 2015 June. We derive number counts of 1.2 mm continuum sources down to 0.02 mJy partly with the assistance of gravitational lensing, and find that the total integrated 1.2 mm flux of the securely identified sources is 22.9^(+6.7)_(-5.6) Jy deg^(-2) that corresponds to 104^(+31)_(-25) % of the extragalactic background light (EBL) measured by COBE observations. These results suggest that the major 1.2 mm EBL contributors are sources with 0.02 mJy, and that very faint 1.2 mm sources with ~< 0.02 mJy contribute negligibly to the EBL with the possible flattening and/or truncation of number counts in this very faint flux regime. To understand the physical origin of our faint ALMA sources, we measure the galaxy bias bg by the counts-in-cells technique, and place a stringent upper limit of bg < 3.5 that is not similar to bg values of massive DRGs and SMGs but comparable to those of UV-bright sBzKs and LBGs. Moreover, in optical and near-infrared (NIR) deep fields, we identify optical-NIR counterparts for 59% of our faint ALMA sources, majority of which have luminosities, colors, and the IRX-beta relation same as sBzKs and LBGs. We thus conclude that about a half of our faint ALMA sources are dust-poor high-z galaxies as known as sBzKs and LBGs in optical studies, and that these faint ALMA sources are not miniature (U)LIRGs simply scaled down with the infrared brightness.

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

We present redshift evolution of galaxy effective radius r_e obtained from the HST samples of ~190,000 galaxies at z=0-10. Our HST samples consist of 176,152 photo-z galaxies at z=0-6 from the 3D-HST+CANDELS catalogue and 10,454 LBGs at z=4-10 identified in CANDELS, HUDF09/12, and HFF parallel fields, providing the largest data set to date for galaxy size evolution studies. We derive r_e with the same technique over the wide-redshift range of z=0-10, evaluating the optical-to-UV morphological K-correction and the selection bias of photo-z galaxies+LBGs as well as the cosmological surface brightness dimming effect. We find that r_e values at a given luminosity significantly decrease towards high-z, regardless of statistics choices. For star-forming galaxies, there is no evolution of the power-law slope of the size-luminosity relation and the median Sersic index (n~1.5). Moreover, the r_e-distribution is well represented by log-normal functions whose standard deviation \sigma_{\ln{r_e}} does not show significant evolution within the range of \sigma_{\ln{r_e}}~0.45-0.75. We calculate the stellar-to-halo size ratio from our r_e measurements and the dark-matter halo masses estimated from the abundance matching study, and obtain a nearly constant value of r_e/r_vir=1.0-3.5% at z=0-8. The combination of the r_e-distribution shape+standard deviation, the constant r_e/r_vir, and n~1.5 suggests a picture that typical high-z star-forming galaxies have disk-like stellar components in a sense of dynamics and morphology over cosmic time of z~0-6. If high-z star-forming galaxies are truly dominated by disks, the r_e/r_vir value and the disk formation model indicate that the specific angular momentum of the disk normalized by the host halo is j_d/m_d=0.5-1. These are statistical results for galaxies' major stellar components, and the detailed study of clumpy sub-components is presented in the paper II.

We present new observations of the highest-redshift quasar known, ULAS J1120+0641, redshift $z=7.084$, obtained in the optical, at near-, mid-, and far-infrared wavelengths, and in the sub-mm. We combine these results with published X-ray and radio observations to create the multiwavelength spectral energy distribution (SED), with the goals of measuring the bolometric luminosity $L_{\rm bol}$, and quantifying the respective contributions from the AGN and star formation. We find three components are needed to fit the data over the wavelength range $0.12-1000\,\mu$m: the unobscured quasar accretion disk and broad-line region, a dusty clumpy AGN torus, and a cool 47K modified black body to characterise star formation. Despite the low signal-to-noise ratio of the new long-wavelength data, the normalisation of any dusty torus model is constrained within $\pm40\%$. We measure a bolometric luminosity $L_{\rm bol}=2.6\pm0.6\times10^{47}\,$erg$\,$s$^{-1}=6.7 \pm 1.6\times10^{13}L_{\odot}$, to which the three components contribute $31\%,32\%,3\%$, respectively, with the remainder provided by the extreme UV $<0.12\,\mu$m. We tabulate the best-fit model SED. We use local scaling relations to estimate a star formation rate (SFR) in the range $60-270\,{\rm M}_\odot$/yr from the [C$\,{\scriptsize \rm II}$] line luminosity and the $158\,\mu$m continuum luminosity. An analysis of the equivalent widths of the [C$\,{\scriptsize \rm II}$] line in a sample of $z>5.7$ quasars suggests that these indicators are promising tools for estimating the SFR in high-redshift quasars in general. At the time observed the black hole was growing in mass more than 100 times faster than the stellar bulge, relative to the mass ratio measured in the local universe, i.e. compared to ${M_{\rm BH}}/{M_{\rm bulge}} \simeq 1.4\times10^{-3}$, for ULAS J1120+0641 we measure ${\dot{M}_{\rm BH}}/{\dot{M}_{\rm bulge}} \simeq 0.2$.

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

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

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

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

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

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

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

A non-LTE analysis of K I resonance lines at 7664.91 and 7698.97 A was carried out for 15 red giants belonging to three globular clusters of different metallicity (M 4, M 13, and M 15) along with two reference early-K giants (rho Boo and alpha Boo), in order to check whether the K abundances are uniform within a cluster and to investigate the behavior of [K/Fe] ratio at the relevant metallicity range of -2.5 <[Fe/H] < -1. We confirmed that [K/H] (as well as [Fe/H]) is almost homogeneous within each cluster to a precision of < ~0.1 dex, though dubiously large deviations are exceptionally seen for two peculiar stars showing signs of considerably increased turbulence in the upper atmosphere. The resulting [K/Fe] ratios are mildly supersolar by a few tenths of dex for three clusters, tending to gradually increase from ~+0.1-0.2 at [Fe/H] ~-1 to ~+0.3 at [Fe/H] ~ -2.5. This result connects reasonably well with the [K/Fe] trend of disk stars (-1 < [Fe/H]) and that of extremely metal-poor stars (-4 <[Fe/H] < -2.5). That is, [K/Fe] appears to continue a gradual increase from [Fe/H]~0 toward a lower metallicity regime down to [Fe/H]~-3, where a broad maximum of [K/Fe]~+0.3-0.4 is attained, possibly followed by a slight downturn at [Fe/H]<~-3.