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

We 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 make use of the ASAGAO, deep 1.2 mm continuum observations of a 26 arcmin$^2$ region in the GOODS-South field obtained with ALMA, to probe dust-enshrouded star formation in $K$-band selected (i.e., stellar mass selected) galaxies, which are drawn from the ZFOURGE catalog. Based on the ASAGAO combined map, which was created by combining ASAGAO and ALMA archival data in the GOODS-South field, we find that 24 ZFOURGE sources have 1.2 mm counterparts with a signal-to-noise ratio $>$ 4.5 (1$\sigma\simeq$ 30 - 70 $\mu$Jy beam$^{-1}$ at 1.2 mm). Their median redshift is estimated to be $z_\mathrm{median}=$ 2.38 $\pm$ 0.14. They generally follow the tight relationship of the stellar mass versus star formation rate (i.e., the main sequence of star-forming galaxies). ALMA-detected ZFOURGE sources exhibit systematically larger infrared (IR) excess (IRX $\equiv L_\mathrm{IR}/L_\mathrm{UV}$) compared to ZFOURGE galaxies without ALMA detections even though they have similar redshifts, stellar masses, and star formation rates. This implies the consensus stellar-mass versus IRX relation, which is known to be tight among rest-frame-UV-selected galaxies, can not fully predict the ALMA detectability of stellar-mass-selected galaxies. We find that ALMA-detected ZFOURGE sources are the main contributors to the cosmic IR star formation rate density at $z$ = 2 - 3.

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.

Extensive surveys with the \textit{Hubble Space Telescope} (HST) over the past decade, targeting some of the most massive clusters in the sky, have uncovered dozens of galaxy-cluster strong lenses. The massive cluster strong-lens scale is typically $\theta_{E}\sim10\arcsec$ to $\sim30-35\arcsec$, with only a handful of clusters known with Einstein radii $\theta_{E}\sim40\arcsec$ or above (for $z_{source}=2$, nominally). Here we report another very large cluster lens, RXC J0032.1+1808 ($z=0.3956$), the second richest cluster in the redMapper cluster catalog and the 85th most massive cluster in the Planck Sunyaev-Zel'dovich catalog. With our Light-Traces-Mass and fully parametric (dPIEeNFW) approaches, we construct strong lensing models based on 18 multiple images of 5 background galaxies newly identified in the \textit{Hubble} data mainly from the \textit{Reionization Lensing Cluster Survey} (RELICS), in addition to a known sextuply imaged system in this cluster. Furthermore, we compare these models to Lenstool and GLAFIC models that were produced independently as part of the RELICS program. All models reveal a large effective Einstein radius of $\theta_{E}\simeq40\arcsec$ ($z_{source}=2$), owing to the obvious concentration of substructures near the cluster center. Although RXC J0032.1+1808 has a very large critical area and high lensing strength, only three magnified high-redshift candidates are found within the field targeted by RELICS. Nevertheless, we expect many more high-redshift candidates will be seen in wider and deeper observations with \textit{Hubble} or \emph{JWST}. Finally, the comparison between several algorithms demonstrates that the total error budget is largely dominated by systematic uncertainties.

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 compute the spherically-averaged power spectrum from four seasons of data obtained for the Epoch of Reionisation (EoR) project observed with the Murchison Widefield Array (MWA). We measure the EoR power spectrum over $k= 0.07-3.0~h$Mpc$^{-1}$ at redshifts $z=6.5-8.7$. The largest aggregation of 110 hours on EoR0 high-band (3,340 observations), yields a lowest measurement of (43~mK)$^2$ = 1.8$\times$10$^3$ mK$^2$ at $k$=0.14~$h$Mpc$^{-1}$ and $z=6.5$ (2$\sigma$ thermal noise plus sample variance). Using the Real-Time System to calibrate and the CHIPS pipeline to estimate power spectra, we select the best observations from the central five pointings within the 2013--2016 observing seasons, observing three independent fields and in two frequency bands. This yields 13,591 2-minute snapshots (453 hours), based on a quality assurance metric that measures ionospheric activity. We perform another cut to remove poorly-calibrated data, based on power in the foreground-dominated and EoR-dominated regions of the two-dimensional power spectrum, reducing the set to 12,569 observations (419 hours). These data are processed in groups of 20 observations, to retain the capacity to identify poor data, and used to analyse the evolution and structure of the data over field, frequency, and data quality. We subsequently choose the cleanest 8,935 observations (298 hours of data) to form integrated power spectra over the different fields, pointings and redshift ranges.

Gravitational lensing sometimes dominates the observed properties of apparently very bright objects. We present morphological properties in the high-resolution (FWHM $\sim0.''15$) Atacama Large Millimeter/submillimeter Array (ALMA) 1-mm map for an ultra-luminous quasar (QSO) at $z=6.30$, SDSS J010013.02+280225.8 (hereafter J0100+2802), whose black hole mass $M_{\rm BH}$ is the most massive ($\sim$ 1.2$\times10^{10}M_{\odot}$) at $z>6$ ever known. We find that the continuum emission of J0100+2802 is resolved into a quadruple system within a radius of $0.''2$, which can be interpreted as either multiple dusty star-forming regions in the host galaxy or multiple images due to strong gravitational lensing. The Mg$_{\rm II}$ absorption and the potential Ly$\alpha$ line features have been identified at $z=2.33$ in the near-infrared spectroscopy towards J0100+2802, and a simple mass model fitting well reproduces the positions and flux densities of the quadruple system, both of which are consistent with the latter interpretation. Although a high-resolution map taken in the Advanced Camera for Survey (ACS) on board Hubble Space Telescope (HST) shows a morphology with an apparently single component, in our fiducial lens mass model it can simply be explained by a $\sim50$ pc scale offset between the ALMA and HST emission regions. In this case, the magnification factor for the observed HST emission is obtained to $\sim450$, reducing the intrinsic $M_{\rm BH}$ estimate to even below $10^{9}\,M_{\odot}$. The confirmation or the rejection of the gravitational lensing scenario is important for our understanding of the super-massive black holes in the early Universe.

We present the results from the search for Lyman Alpha emitters (LAEs) in the proximity of 11 CIV absorption systems at $z > 4.7$ in the spectrum of the QSO J1030+0524, using data from MUSE. We have found multiple LAE candidates close to four CIV systems at $z_{CIV} = 4.94$-5.74 with log$_{10}$(N CIV[cm$^{-2}$])$ > 13.5.$ At z = 5-6, CIV systems with W0(CIV) $> 0.2$ {\AA} seem more likely to have galaxies with Ly${\alpha}$ emission within $\rho < 200$ proper kpc (4/5 cases), than the CIV systems with W0(CIV)$ < 0.2$ {\AA} (0/6 cases). The impact parameter of LAE-CIV systems with equivalent widths W0(CIV)$ > 0.5 $ {\AA} is in the range 11 $ < \rho <$ 200 proper kpc (pkpc). Furthermore, all candidates are in the luminosity range 0.18-1.15 L*$(Ly{ \alpha }, z = 5.7)$ indicating that the environment of CIV systems within 200 pkpc is populated by the faint end of the Ly${\alpha}$ luminosity function. We report a 0.28 L*$(Ly{\alpha})$ galaxy at a separation of $\rho = 10$ pkpc from a strong CIV absorption (log10(NCIV[cm$^{-2}$]) = 14.52) at $z_{CIV} = 5.72419$. The prevalence of sub-L*(Ly$\alpha$) galaxies in the proximity of z $> 4.9$ CIV systems suggest that the absorbing material is rather young, likely ejected in the recent past of the galaxies at these redshifts. The connection between faint LAEs and high-ionization absorption systems reported in this work, is potentially a consequence of the role of low mass galaxies in the early evolution of the circum-galactic and intergalactic media.

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

We present deep rest-frame UV spectroscopic observations using the Gran Telescopio Canarias of six gravitationally lensed Lya emitters (LAEs) at $2.36<z<2.82$ selected from the BELLS GALLERY survey. By taking the magnifications into account, we show that LAEs can be as luminous as L(Lya) = 30x10$^{42}$ erg s-1 and M(UV) = -23 (AB) without invoking an AGN component, in contrast with previous findings. We measure Lya rest-frame equivalent widths, EW(Lya), ranging from 16\AA to 50\AA and Lya escape fractions, fesc(Lya), from 10% to 40%. Large EW(Lya) and fesc(Lya) are found predominantly in LAEs showing weak low-ionization ISM absorption (EW < 1\AA) and narrow Lya profiles (< 300 km s-1 FWHM) with their peak close (< 80 km s-1) to their systemic redshifts, suggestive of less scatter from low HI column densities that favours the escape of Lya photons. We infer stellar metallicities of Z/Zsun ~ 0.2 in almost all LAEs by comparing the P-Cygni profiles of the wind lines NV1240\AA and CIV1549\AA with those from stellar synthesis models. We also find a trend between M(UV) and the velocity offset of ISM absorption lines, such as the most luminous LAEs experience stronger outflows. The most luminous LAEs show star formation rates up to 180 Msun yr-1, yet they appear relatively blue ($\beta$(UV) ~ -1.8 to -2.0) showing evidence of little dust attenuation (E(B-V) = 0.10-0.14). These luminous LAEs may be particular cases of young starburst galaxies that have had no time to form large amounts of dust. If so, they are ideal laboratories to study the early phase of massive star formation, stellar and dust mass growth, and chemical enrichment histories of starburst galaxies at high-z.

We present velocities of galactic outflows in seven star-forming galaxies at $z=$5-6 with stellar masses of $M_* \sim10^{10.1}\ \rm{M_\odot}$. Although it is challenging to observationally determine the outflow velocities, we overcome this by using ALMA [CII]158 $\mu$m emission lines for systemic velocities and deep Keck spectra with metal absorption lines for velocity profiles available to date. We construct a composite Keck spectrum of the galaxies at $z=$5-6 with the [CII]-systemic velocities, and fit outflow-line profiles to the SiII1260, CII1335, and SiIV1394,1403 absorption lines in the composite spectrum. We measure the maximum (90%) and central outflow velocities to be $v_\rm{max}=700^{+180}_{-110}\ \rm{km\ s^{-1}}$ and $v_\rm{out}= 400^{+100}_{-150}\ \rm{km\ s^{-1}}$ on average, respectively, showing no significant differences between the outflow velocities derived with the low to high-ionization absorption lines. For $M_* \sim10^{10.1}\ \rm{M_\odot}$, we find that the $v_\rm{max}$ value of our $z=$5-6 galaxies is 3 times higher than those of $z\sim0$ galaxies and comparable to $z\sim2$ galaxies. Estimating the halo circular velocity $v_\rm{cir}$ from the stellar masses and the abundance matching results, we investigate a $v_\rm{max}$-$v_\rm{cir}$ relation. Interestingly, $v_\rm{max}$ for galaxies with $M_*=10^{10.0-10.8}\ \rm{M_\odot}$ shows a clear positive correlation with $v_\rm{cir}$ and/or the galaxy star formation rate over $z=$0-6 with a small scatter of $\simeq \pm 0.1$ dex, which is in good agreement with theoretical predictions. This positive correlation suggests that the outflow velocity is physically related to the halo circular velocity, and that the redshift evolution of $v_\rm{max}$ at fixed $M_*$ is explained by the increase in $v_\rm{cir}$ toward high redshift.

We report the discovery of 10-kpc scale [CII] 158um halos surrounding star-forming galaxies in the early Universe. We choose deep ALMA data of 18 galaxies each with a star-formation rate of ~ 10-70 Msun with no signature of AGN whose [CII] lines are individually detected at z=5.153-7.142, and conduct stacking of the [CII] lines and dust-continuum in the uv-visibility plane. The radial profiles of the surface brightnesses show a 10-kpc scale [CII] halo at the 9.2sigma level significantly extended more than the HST stellar continuum data by a factor of ~5 on the exponential-profile basis, as well as the dust continuum. We also compare the radial profiles of [CII] and Lya halos universally found in star-forming galaxies at this epoch, and find that the scale lengths agree within the 1sigma level. While two independent hydrodynamical zoom-in simulations match the dust and stellar continuum properties, the simulations cannot reproduce the extended [CII] line emission. The existence of the extended [CII] halo is the evidence of outflow remnants in the early galaxies and suggest that the outflows may be dominated by cold-mode outflows expelling the neutral gas.

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

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

We present the results of the deep Subaru/FOCAS and Keck/MOSFIRE spectroscopy for four spatially extended [O III]$\lambda\lambda$4959,5007 sources, dubbed [O III] blobs, at $z=0.6-0.8$ that are originally pinpointed by large-area Subaru imaging surveys. The line diagnostics of the rest-frame optical lines suggests that only one [O III] blob, OIIIB-3, presents an AGN signature, indicating that hot gas of the rest of the [O III] blobs is heated by star formation. One of such star-forming [O III] blobs, OIIIB-4, at $z=0.838$ has an [O III] equivalent width of $845\pm27$ \r{A} and an [O III] to [O II]$\lambda\lambda$3726,3729 ratio of [O III]/[O II]= $6.5\pm2.7$ that are as high as those of typical green peas (Cardamone et al. 2009). The spatially resolved spectrum of OIIIB-4 shows [O III]/[O II]= $5-10$ over $14$ kpc in the entire large [O III] extended regions of OIIIB-4, unlike the known green peas whose strong [O III] emission region is compact. Moreover, OIIIB-4 presents no high ionization emission lines, unlike green beans that have extended [O III] emission with a type-2 AGN. OIIIB-4 is thus a giant green pea, which is a low stellar mass ($7\times10^7$ $M_\odot$) galaxy with a very high specific star formation rate (sSFR = $2\times10^2\, {\rm Gyr}^{-1}$), a high ionization parameter ($q_{ion} \sim 3\times10^8\, {\rm cm\,s^{-1}}$), and a low metallicity similar to those of green peas. Neither an AGN-light echo nor a fast radiative shock likely takes place due to the line diagnostics for spatially-resolved components of OIIIB-4 and no detections of He II$\lambda$4686 or [Ne V]$\lambda\lambda$3346,3426 lines that are fast-radiative shock signatures. There is a possibility that the spatially-extended [O III] emission of OIIIB-4 is originated from outflowing gas produced by the intense star formation in a density-bounded ionization state.

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 discovery of the first spectroscopically resolved C II /C II* 1334, 1335A doublet in the Lyman-break galaxy J0215-0555 at z = 5.754. The separation of the resonant and fluorescent emission channels was possible thanks to the large redshift of the source and long integration time, as well as the small velocity width of the feature, 0.6 +- 0.2A. We model this emission and find that at least two components are required to reproduce the combination of morphologies of C II* emission, C II absorption and emission, and Lyman-alpha emission from the object. We suggest that the close alignment between the fluorescence and Lyman-alpha emission could indicate an ionisation escape channel within the object. While the faintness of such a C II /C II* doublet makes it prohibitively difficult to pursue for similar systems with current facilities, we suggest it can become a valuable porosity diagnostic in the era of JWST and the upcoming generations of ELTs.

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.

Measurements of stellar properties of galaxies when the universe was less than one billion years old yield some of the only observational constraints of the onset of star formation. We present here the inclusion of \textit{Spitzer}/IRAC imaging in the spectral energy distribution fitting of the seven highest-redshift galaxy candidates selected from the \emph{Hubble Space Telescope} imaging of the Reionization Lensing Cluster Survey (RELICS). We find that for 6/8 \textit{HST}-selected $z\gtrsim8$ sources, the $z\gtrsim8$ solutions are still strongly preferred over $z\sim$1-2 solutions after the inclusion of \textit{Spitzer} fluxes, and two prefer a $z\sim 7$ solution, which we defer to a later analysis. We find a wide range of intrinsic stellar masses ($5\times10^6 M_{\odot}$ -- $4\times10^9$ $M_{\odot}$), star formation rates (0.2-14 $M_{\odot}\rm yr^{-1}$), and ages (30-600 Myr) among our sample. Of particular interest is Abell1763-1434, which shows evidence of an evolved stellar population at $z\sim8$, implying its first generation of star formation occurred just $< 100$ Myr after the Big Bang. SPT0615-JD, a spatially resolved $z\sim10$ candidate, remains at its high redshift, supported by deep \textit{Spitzer}/IRAC data, and also shows some evidence for an evolved stellar population. Even with the lensed, bright apparent magnitudes of these $z \gtrsim 8$ candidates (H = 26.1-27.8 AB mag), only the \textit{James Webb Space Telescope} will be able further confirm the presence of evolved stellar populations early in the universe.

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

We present very faint dropout galaxies at z~6-9 with a stellar mass M* down to M*~10^6Mo that are found in deep optical/near-infrared (NIR) images of the full data sets of the Hubble Frontier Fields (HFF) program in conjunction with deep ground-based and Spitzer images and gravitational lensing magnification effects. We investigate stellar populations of the HFF dropout galaxies with the optical/NIR photometry and BEAGLE models made of self-consistent stellar population synthesis and photoionization models, carefully including strong nebular emission impacting on the photometry. We identify 357 galaxies with M*~10^6-10^9Mo, and find that a stellar mass to UV luminosity L_UV ratio M*/L_UV is nearly constant at M*~10^6-10^9Mo. Our best-estimate M*/L_UV function is comparable to a model of star-formation duration time of 100 Myr, but 2-7 times higher than the one of 10 Myr assumed in a previous study (at the 5sigma level) that would probably underestimate M* of faint galaxies. We derive the galaxy stellar mass functions (GSMFs) at z~6-9 that agree with those obtained by previous studies with no M*/L_UV assumptions at M*>~10^8Mo, and that extends to M*~10^6Mo. Estimating the stellar mass densities rho* with the GSMFs, we find that rho* smoothly increases from log(rho*/[Mo Mpc^(-3)])=5.91 +0.75/-0.65 at z~9 to 6.21 +0.39/-0.37 at z~6-7, which is consistent with the one estimated from star-formation rate density measurements. In conjunction with the estimates of the galaxy effective radii R_e on the source plane, we have pinpointed two objects with low stellar masses (M*<=10^7Mo) and very compact morphologies (R_e<=40 physical pc) that are comparable with those of globular clusters (GCs) in the Milky Way today. These objects are candidates of star clusters that should be a part or a dominant component of high-redshift low-mass galaxy, some of which may be related to GCs today.

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.

Strong gravitational lensing by clusters of galaxies probes the mass distribution at the core of each cluster and magnifies the universe behind it. MACS J0417.5-1154 at z=0.443 is one of the most massive clusters known based on weak lensing, X-ray, and Sunyaev-Zel'dovich analyses. Here we compute a strong lens model of MACS J0417 based on Hubble Space Telescope imaging observations collected, in part, by the Reionization Lensing Cluster Survey (RELICS), and recently reported spectroscopic redshifts from the MUSE instrument on the Very Large Telescope (VLT). We measure an Einstein radius of ThetaE=36'' at z = 9 and a mass projected within 200 kpc of M(200 kpc) = 1.78+0.01-0.03x10**14Msol. Using this model, we measure a ratio between the mass attributed to cluster-member galaxy halos and the main cluster halo of order 1:100. We assess the probability to detect magnified high-redshift galaxies in the field of this cluster, both for comparison with RELICS HST results and as a prediction for the James Webb Space Telescope (JWST) Guaranteed Time Observations upcoming for this cluster. Our lensing analysis indicates that this cluster has similar lensing strength to other clusters in the RELICS program. Our lensing analysis predicts a detection of at least a few z~6-8 galaxies behind this cluster, at odds with a recent analysis that yielded no such candidates in this field. Reliable strong lensing models are crucial for accurately predicting the intrinsic properties of lensed galaxies. As part of the RELICS program, our strong lensing model produced with the Lenstool parametric method is publicly available through the Mikulski Archive for Space Telescopes (MAST).

We report detections of two 1.2 mm continuum sources ($S_\mathrm{1.2mm}$ ~ 0.6 mJy) without any counterparts in the deep $H$- and/or $K$-band image (i.e., $K$-band magnitude $\gtrsim$ 26 mag). These near-infrared-dark faint millimeter sources are uncovered by ASAGAO, a deep and wide-field ($\simeq$ 26 arcmin$^2$) Atacama Large Millimeter/submillimeter Array (ALMA) 1.2 mm survey. One has a red IRAC (3.6 and 4.5 $\mu$m) counterpart, and the other has been independently detected at 850 and 870 $\mu$m using SCUBA2 and ALMA Band 7, respectively. Their optical to radio spectral energy distributions indicate that they can lie at $z \gtrsim$ 3-5 and can be in the early phase of massive galaxy formation. Their contribution to the cosmic star formation rate density is estimated to be ~ 1 $\times$ 10$^{-3}$ $M_\odot$ yr$^{-1}$ Mpc$^{-3}$ if they lie somewhere in the redshift range of $z$ ~ 3-5. This value can be consistent with, or greater than that of bright submillimeter galaxies ($S_\mathrm{870\mu m}>$ 4.2 mJy) at $z$ ~ 3-5. We also uncover 3 more candidates near-infrared-dark faint ALMA sources without any counterparts ($S_\mathrm{1.2mm}$ ~ 0.45-0.86 mJy). These results show that an unbiased ALMA survey can reveal the dust-obscured star formation activities, which were missed in previous deep optical/near-infrared surveys.

Large surveys of galaxy clusters with the Hubble and Spitzer Space Telescopes, including CLASH and the Frontier Fields, have demonstrated the power of strong gravitational lensing to efficiently deliver large samples of high-redshift galaxies. We extend this strategy through a wider, shallower survey named RELICS, the Reionization Lensing Cluster Survey. This survey, described here, was designed primarily to deliver the best and brightest high-redshift candidates from the first billion years after the Big Bang. RELICS observed 41 massive galaxy clusters with Hubble and Spitzer at 0.4-1.7um and 3.0-5.0um, respectively. We selected 21 clusters based on Planck PSZ2 mass estimates and the other 20 based on observed or inferred lensing strength. Our 188-orbit Hubble Treasury Program obtained the first high-resolution near-infrared images of these clusters to efficiently search for lensed high-redshift galaxies. We observed 46 WFC3/IR pointings (~200 arcmin^2) with two orbits divided among four filters (F105W, F125W, F140W, and F160W) and ACS imaging as needed to achieve single-orbit depth in each of three filters (F435W, F606W, and F814W). As previously reported by Salmon et al., we discovered 322 z ~ 6 - 10 candidates, including the brightest known at z ~ 6, and the most distant spatially-resolved lensed arc known at z ~ 10. Spitzer IRAC imaging (945 hours awarded, plus 100 archival) has crucially enabled us to distinguish z ~ 10 candidates from z ~ 2 interlopers. For each cluster, two HST observing epochs were staggered by about a month, enabling us to discover 11 supernovae, including 3 lensed supernovae, which we followed up with 20 orbits from our program. We delivered reduced HST images and catalogs of all clusters to the public via MAST and reduced Spitzer images via IRSA. We have also begun delivering lens models of all clusters, to be completed before the JWST GO call for proposals.

Strong gravitational lensing (SL) is a powerful means to map the distribution of dark matter. In this work, we perform a SL analysis of the prominent X-ray cluster RXJ0152.7-1357 (z=0.83, also known as CL 0152.7-1357) in \textit{Hubble Space Telescope} images, taken in the framework of the Reionization Lensing Cluster Survey (RELICS). On top of a previously known $z=3.93$ galaxy multiply imaged by RXJ0152.7-1357, for which we identify an additional multiple image, guided by a light-traces-mass approach we identify seven new sets of multiply imaged background sources lensed by this cluster, spanning the redshift range [1.79-3.93]. A total of 25 multiple images are seen over a small area of ~0.4 $arcmin^2$, allowing us to put relatively high-resolution constraints on the inner matter distribution. Although modestly massive, the high degree of substructure together with its very elongated shape make RXJ0152.7-1357 a very efficient lens for its size. This cluster also comprises the third-largest sample of z~6-7 candidates in the RELICS survey. Finally, we present a comparison of our resulting mass distribution and magnification estimates with those from a Lenstool model. These models are made publicly available through the MAST archive.

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

We present deep $J$ and $H$-band images in the extended Great Observatories Origins Deep Survey-North (GOODS-N) field covering an area of 0.22 $\rm{deg}^{2}$. The observations were taken using WIRCam on the 3.6-m Canada France Hawaii Telescope (CFHT). Together with the reprocessed $K_{\rm s}$-band image, the $5\sigma$ limiting AB magnitudes (in 2" diameter apertures) are 24.7, 24.2, and 24.4 AB mag in the $J$, $H$, and $K_{\rm s}$ bands, respectively. We also release a multi-band photometry and photometric redshift catalog containing 93598 sources. For non-X-ray sources, we obtained a photometric redshift accuracy $\sigma_{\mathrm{NMAD}}=0.036$ with an outlier fraction $\eta = 7.3\%$. For X-ray sources, which are mainly active galactic nuclei (AGNs), we cross-matched our catalog with the updated 2M-CDFN X-ray catalog from Xue et al. (2016) and found that 658 out of 683 X-ray sources have counterparts. $GALEX$ UV data are included in the photometric redshift computation for the X-ray sources to give $\sigma_{\mathrm{NMAD}} = 0.040$ with $\eta=10.5\%$. Our approach yields more accurate photometric redshift estimates compared to previous works in this field. In particular, by adopting AGN-galaxy hybrid templates, our approach delivers photometric redshifts for the X-ray counterparts with fewer outliers compared to the 3D-HST catalog, which fit these sources with galaxy-only templates.

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

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

We present the survey design, data reduction, construction of images, and source catalog of the Atacama Large Millimeter/submillimeter Array (ALMA) twenty-six arcmin^2 survey of GOODS-S at one-millimeter (ASAGAO). ASAGAO is a deep (1sigma ~ 61 uJy/beam for a 250 klambda-tapered map with a synthesized beam size of 0.51" x 0.45") and wide area (26 arcmin^2) survey on a contiguous field at 1.2 mm. By combining with ALMA archival data in the GOODS-South field, we obtained a deeper map in the same region (1sigma ~ 30 uJy/beam for a deep region with a 250 klambda-taper, and a synthesized beam size of 0.59" x 0.53"), providing the largest sample of sources (25 sources at >=5.0sigma, 45 sources at >=4.5sigma) among ALMA blank-field surveys to date. The number counts shows that 52(+11 -8)% of the extragalactic background light at 1.2 mm is resolved into discrete sources at S1.2m > 135 uJy. We create infrared (IR) luminosity functions (LFs) in the redshift range of z = 1-3 from the ASAGAO sources with KS-band counterparts, and constrain the faintest luminosity of the LF at 2.0 < z < 3.0. The LFs are consistent with previous results based on other ALMA and SCUBA-2 observations, which suggest a positive luminosity evolution and negative density evolution with increasing redshift. We find that obscured star-formation of sources with IR luminosities of log(L(IR)/Lsun)} ~> 11.8 account for ~~60%-90% of the z ~ 2 cosmic star-formation rate density.

We present a strong-lensing analysis of four massive galaxy clusters imaged with the Hubble Space Telescope in the Reionization Lensing Cluster Survey. We use a Light-Traces-Mass technique to uncover sets of multiply images and constrain the mass distribution of the clusters. These mass models are the first published for Abell S295 and MACS J0159.8-0849, and are improvements over previous models for Abell 697 and MACS J0025.4-1222. Our analysis for MACS J0025.4-1222 and Abell S295 shows a bimodal mass distribution supporting the merger scenarios proposed for these clusters. The updated model for MACS J0025.4-1222 suggests a substantially smaller critical area than previously estimated. For MACS J0159.8-0849 and Abell 697 we find a single peak and relatively regular morphology, revealing fairly relaxed clusters. Despite being less prominent lenses, three of these clusters seem to have lensing strengths, i.e. cumulative area above certain magnification, similar to the Hubble Frontier Fields clusters (e.g., A($\mu>5$) $\sim 1-3$ arcmin$^2$, A($\mu>10$) $\sim 0.5-1.5$ arcmin$^2$), which in part can be attributed to their merging configurations. We make our lens models publicly available through the Mikulski Archive for Space Telescopes. Finally, using Gemini-N/GMOS spectroscopic observations we detect a single emission line from a high-redshift $J_{125}\simeq25.7$ galaxy candidate lensed by Abell 697. While we cannot rule out a lower-redshift solution, we interpret the line as Ly$\alpha$ at $z=5.800\pm 0.001$, in agreement with its photometric redshift and dropout nature. Within this scenario we measure a Ly$\alpha$ rest-frame equivalent width of $52\pm22$ \AA, and an observed Gaussian width of $117\pm 15$ km/s.

Detection of the redshifted 21cm-line signal from neutral hydrogen in the intergalactic medium (IGM) during the Epoch of Reionization (EoR) is complicated by intense foregrounds such as galactic synchrotron and extragalactic radio galaxies. The 21cm-Lyman-$\alpha$ emitter(LAE) cross-correlation is one of the tools available to reduce the foreground effects because the foreground emission from such radio sources is statistically independent of LAE distribution. LAE surveys during the EoR at redshifts $z=6.6$ and $7.3$ are ongoing by the Subaru Hyper Suprime-Cam (HSC). Additionally, Prime Focus Spectrograph (PFS) will provide precise redshift information of the LAEs discovered by the HSC survey. In this paper, we investigate the detectability of the 21cm signal with the 21cm-LAE cross-correlation by using our improved reionization simulations. We also focus on the error budget and evaluate it quantitatively in order to consider a strategy to improve the signal-to-noise ratio. In addition, we explore an expansion of the LAE survey to suggest optimal survey parameters and show a potential to measure a characteristic size of ionized bubbles via the turnover scale of the cross-power spectrum. As a result, we find that the Murchison Widefield Array (MWA) has ability to detect the cross-power spectrum signal on large scales by combining LAE Deep field survey of HSC. We also show that the sensitivity is improved dramatically at small scales by adding redshift information from the PFS measurements. The Square Kilometre Array (SKA) has a potential to measure the turnover scale with an accuracy of $6\times10^{-3}~{\rm Mpc^{-1}}$.

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

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.

A fundamental quest of modern astronomy is to locate the earliest galaxies and study how they influenced the intergalactic medium a few hundred million years after the Big Bang. The abundance of star-forming galaxies is known to decline from redshifts of about 6 to 10, but a key question is the extent of star formation at even earlier times, corresponding to the period when the first galaxies might have emerged. Here we present spectroscopic observations of MACS1149-JD1, a gravitationally lensed galaxy observed when the Universe was less than four per cent of its present age. We detect an emission line of doubly ionized oxygen at a redshift of $9.1096\pm0.0006$, with an uncertainty of one standard deviation. This precisely determined redshift indicates that the red rest-frame optical colour arises from a dominant stellar component that formed about 250 million years after the Big Bang, corresponding to a redshift of about 15. Our results indicate the it may be possible to detect such early episodes of star formation in similar galaxies with future telescopes.

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.

Strong gravitational lensing by galaxy clusters has become a powerful tool for probing the high-redshift Universe, magnifying distant and faint background galaxies. Reliable strong lensing (SL) models are crucial for determining the intrinsic properties of distant, magnified sources and for constructing their luminosity function. We present here the first SL analysis of MACS J0308.9+2645 and PLCK G171.9-40.7, two massive galaxy clusters imaged with the Hubble Space Telescope in the framework of the Reionization Lensing Cluster Survey (RELICS). We use the Light-Traces-Mass modeling technique to uncover sets of multiply imaged galaxies and constrain the mass distribution of the clusters. Our SL analysis reveals that both clusters have particularly large Einstein radii ($\theta_E>30"$ for a source redshift of $z_s=2$), providing fairly large areas with high magnifications, useful for high-redshift galaxy searches ($\sim2$ arcmin$^{2}$ with $\mu>5$ to $\sim1$ arcmin$^{2}$ with $\mu>10$, similar to a typical \textit{Hubble Frontier Fields} cluster). We also find that MACS J0308.9+2645 hosts a promising, apparently bright (J$\sim23.2-24.6$ AB), multiply imaged high-redshift candidate at $z\sim6.4$. These images are amongst the brightest high-redshift candidates found in RELICS. Our mass models, including magnification maps, are made publicly available for the community through the Mikulski Archive for Space Telescopes.

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 .

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

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

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

The most distant galaxies known are at z~10-11, observed 400-500 Myr after the Big Bang. The few z~10-11 candidates discovered to date have been exceptionally small- barely resolved, if at all, by the Hubble Space Telescope. Here we present the discovery of SPT0615-JD, a fortuitous z~10 (z_phot=9.9+/-0.6) galaxy candidate stretched into an arc over ~2.5" by the effects of strong gravitational lensing. Discovered in the Reionization Lensing Cluster Survey (RELICS) Hubble Treasury program and companion S-RELICS Spitzer program, this candidate has a lensed H-band magnitude of 25.7+/-0.1 AB mag. With a magnification of \mu~4-7 estimated from our lens models, the de-lensed intrinsic magnitude is 27.6+/-0.3 AB mag, and the half-light radius is r_e<0.8 kpc, both consistent with other z>9 candidates. The inferred stellar mass (log [M* /M_Sun]=9.7^{+0.7}_{-0.5}) and star formation rate (\log [SFR/M_Sun yr^{-1}]=1.3^{+0.2}_{-0.3}) indicate that this candidate is a typical star-forming galaxy on the z>6 SFR-M* relation. We note that three independent lens models predict two counterimages, at least one of which should be of a similar magnitude to the arc, but these counterimages are not yet detected. Counterimages would not be expected if the arc were at lower redshift. However, the only spectral energy distributions capable of fitting the Hubble and Spitzer photometry well at lower redshifts require unphysical combinations of z~2 galaxy properties. The unprecedented lensed size of this z~10 candidate offers the potential for the James Webb Space Telescope to study the geometric and kinematic properties of a galaxy observed 500 Myr after the Big Bang.

We investigate the X-ray active galactic nucleus (AGN) properties of millimeter galaxies in the Great Observatories Origins Deep Survey South (GOODS-S) field detected with the Atacama Large Millimeter/submillimeter Array (ALMA), by utilizing the Chandra 7-Ms data, the deepest X-ray survey to date. Our millimeter galaxy sample comes from the ASAGAO survey covering 26 arcmin$^2$ (12 sources at a 1.2-mm flux-density limit of $\approx$0.6 mJy), supplemented by the deeper but narrower 1.3-mm survey of a part of the ASAGAO field by Dunlop et al.\ (2017). Fourteen out of the total 25 millimeter galaxies have Chandra counterparts. The observed AGN fractions at $z=1.5-3$ is found to be 90$^{+8}_{-19}$\% and $57^{+23}_{-25}$\% for the ultra/luminous infrared galaxies with $\log L_{\rm IR}/L_{\odot} = 12-12.8$ and $\log L_{\rm IR}/L_{\odot} = 11.5-12$, respectively. The majority ($\sim$2/3) of the ALMA and/or Herschel detected X-ray AGNs at $z=1.5-3$ appear to be star-formation dominant populations, having $L_{\rm X}$/ $L_{\rm IR}$ ratios smaller than the "simultaneous evolution" value expected from the local black-hole mass to stellar mass ($M_{\rm BH}$-$M_*$) relation. On the basis of the $L_{\rm X}$ and stellar mass relation, we infer that a large fraction of star-forming galaxies at $z=1.5-3$ have black hole masses smaller than those expected from the local $M_{\rm BH}$-$M_*$ relation. This is opposite to previous reports on luminous AGN at same redshifts detected in wider and shallower surveys, which are subject to selection biases against lower luminosity AGN. Our results are consistent with an evolutionary scenario that star formation occurs first, and an AGN-dominant phase follows later, in objects finally evolving into galaxies with classical bulges.

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

We study the heating mechanisms and Ly{\alpha} escape fractions of 35 Ly{\alpha} blobs (LABs) at z = 3.1 in the SSA22 field. Dust continuum sources have been identified in 11 of the 35 LABs, all with star formation rates (SFRs) above 100 Msun/yr. Likely radio counterparts are detected in 9 out of 29 investigated LABs. The detection of submm dust emission is more linked to the physical size of the Ly{\alpha} emission than to the Ly{\alpha} luminosities of the LABs. A radio excess in the submm/radio detected LABs is common, hinting at the presence of active galactic nuclei. Most radio sources without X-ray counterparts are located at the centers of the LABs. However, all X-ray counterparts avoid the central regions. This may be explained by absorption due to exceptionally large column densities along the line-of-sight or by LAB morphologies, which are highly orientation dependent. The median Ly{\alpha} escape fraction is about 3\% among the submm-detected LABs, which is lower than a lower limit of 11\% for the submm-undetected LABs. We suspect that the large difference is due to the high dust attenuation supported by the large SFRs, the dense large-scale environment as well as large uncertainties in the extinction corrections required to apply when interpreting optical data.

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

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

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

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

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

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

Lyman alpha halos are observed ubiquitously around star-forming galaxies at high redshift, but their origin is still a matter of debate. We demonstrate that the emission from faint unresolved satellite sources, $M_{\rm UV} \gtrsim -17$, clustered around the central galaxies may play a major role in generating spatially extended Ly$\alpha$, continuum (${\rm UV + VIS}$) and H$\alpha$ halos. We apply the analytic formalism developed in Mas-Ribas & Dijkstra (2016) to model the halos around Lyman Alpha Emitters (LAEs) at $z=3.1$, for several different satellite clustering prescriptions. In general, our UV and Ly$\alpha$ surface brightness profiles match the observations well at $20\lesssim r \lesssim 40$ physical kpc from the centers of LAEs. We discuss how our profiles depend on various model assumptions and how these can be tested and constrained with future H$\alpha$ observations by the James Webb Space Telescope (JWST). Our analysis shows how spatially extended halos constrain (i) the presence of otherwise undetectable satellite sources, (ii) the integrated, volumetric production rates of Ly$\alpha$ and LyC photons, and (iii) their population-averaged escape fractions. These quantities are all directly relevant for understanding galaxy formation and evolution and, for high enough redshifts, cosmic reionization.

We conduct a systematic search for galaxies at $z=0.1-1.5$ with [OII]$\lambda3727$, [OIII]$\lambda5007$, or H$\alpha\lambda6563$ emission lines extended over at least 30 kpc by using deep narrowband and broadband imaging in Subaru-XMM Deep Survey (SXDS) field. These extended emission-line galaxies are dubbed [OII], [OIII], or H$\alpha$ blobs. Based on a new selection method that securely selects extended emission-line galaxies, we find 77 blobs at $z=0.40-1.46$ with the isophotal area of emission lines down to $1.2\times10^{-18}$ erg s$^{-1}$ cm$^{-2}$ kpc$^{-2}$. Four of them are spectroscopically confirmed to be [OIII] blobs at $z=0.83$. We identify AGN activities in 8 blobs with X-ray and radio data and find that the fraction of AGN contribution increases with increasing isophotal area of the extended emission. With the Kolmogorov-Smirnov (KS) and Anderson-Darling tests, we confirm that the stellar-mass distributions of H$\alpha$ and [OII] blobs are not drawn from those of the emitters at the $>90$% confidence level in that H$\alpha$ and [OII] blobs are located at the massive end of the distributions, but cannot reject null hypothesis of being the same distributions in terms of the specific star formation rates. It is suggested that galactic-scale outflows tend to be more prominent in more massive star-forming galaxies. Exploiting our sample homogeneously selected over the large area, we derive the number densities of blobs at each epoch. The number densities of blobs decrease drastically with redshifts at the rate that is larger than that of the decrease of cosmic star formation densities.

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

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

We present 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 report the serendipitous discovery of a quadruply lensed source at $z_{\rm s}=3.76$, HSC~J115252+004733, from the Hyper Suprime-Cam (HSC) Survey. The source is lensed by an early-type galaxy at $z_{\rm l}=0.466$ and a satellite galaxy. Here, we investigate the properties of the source by studying its size and luminosity from the imaging and the luminosity and velocity width of the Ly-$\alpha$ line from the spectrum. Our analyses suggest that the source is most probably a low-luminosity active galactic nucleus (LLAGN) but the possibility of it being a compact bright galaxy (e.g., a Lyman-$\alpha$ emitter or Lyman Break Galaxy) cannot be excluded. The brighter pair of lensed images appears point-like except in the HSC $i$-band (with a seeing $\sim0.5"$). The extended emission in the $i$-band image could be due to the host galaxy underneath the AGN, or alternatively, due to a highly compact lensed galaxy (without AGN) which appears point-like in all bands except in $i$-band. We also find that the flux ratio of the brighter pair of images is different in the Ks-band compared to optical wavelengths. Phenomena such as differential extinction and intrinsic variability cannot explain this chromatic variation. While microlensing from stars in the foreground galaxy is less likely to be the cause, it cannot be ruled out completely. If the galaxy hosts an AGN, then this represents the highest redshift quadruply imaged AGN known to date, enabling study of a distant LLAGN. Discovery of this unusually compact and faint source demonstrates the potential of the HSC survey.

Using galaxies as background light sources to map intervening Lya absorption is a novel approach to study the interplay among galaxies, the circum-galactic medium (CGM), and the intergalactic medium (IGM). Introducing a new measure of z = 3.1$ HI Lya absorption relative to the cosmic mean, Delta_NB497, estimated from photometric data of star-forming galaxies at 3.3 < z < 3.5, we have made two-dimensional Delta_NB497 maps in the z = 3.1 SSA22 proto-cluster region and two control fields (SXDS and GOODS-N fields) with a spatial resolution of ~ 5 comoving Mpc. The Delta_NB497 measurements in the SSA22 field are systematically larger than those in the control fields, and this HI absorption enhancement extends more than 50 comoving Mpc. The field-averaged (i.e., ~50 comoving Mpc scale) Delta_NB497 and the overdensity of Lya emitters (LAEs) seem to be correlated, while there is no clear dependency of the Delta_NB497 on the local LAE overdensity in a few comoving Mpc scale. These results suggest that diffuse HI gas spreads out in/around the SSA22 proto-cluster. We have also found an enhancement of Delta_NB497 at a projected distance < 100 physical kpc from the nearest z = 3.1 galaxies at least in the SSA22 field, which is probably due to HI gas associated with the CGM of individual galaxies. The HI absorption enhancement in the CGM-scale tends to be weaker around galaxies with stronger Lya emission, which suggests that the Lya escape fraction from galaxies depends on hydrogen neutrality in the CGM.

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 report the result from observations conducted with the Atacama Large Millimeter/submillimeter Array (ALMA) to detect [CII] 158 um fine structure line emission from galaxies embedded in one of the most spectacular Lyman-alpha blobs (LABs) at z=3.1, SSA22-LAB1. Of three dusty star-forming galaxies previously discovered by ALMA 860 um dust continuum survey toward SSA22-LAB1, we detected the [CII] line from one, LAB1-ALMA3 at z=3.0993+/-0.0004. No line emission was detected, associated with the other ALMA continuum sources or from three rest-frame UV/optical selected z_spec~3.1 galaxies within the field of view. For LAB1-ALMA3, we find relatively bright [CII] emission compared to the infrared luminosity (L_[CII]/L_[CII]) and an extremely high [CII] 158 um and [NII] 205 um emission line ratio (L_[CII]/L_[NII]>55). The relatively strong [CII] emission may be caused by abundant photodissociation regions and sub-solar metallicity, or by shock heating. The origin of the unusually strong [CII] emission could be causally related to the location within the giant LAB, although the relationship between extended Lyman-alpha emission and ISM conditions of associated galaxies is yet to be understand.

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

We present Keck/MOSFIRE spectra of the diagnostic nebular emission lines [OIII]5007,4959, [OII]3727, and H-beta for a sample of 15 redshift z=3.1-3.7 Ly-alpha emitters (LAEs) and Lyman break galaxies (LBGs). In conjunction with spectra from other surveys, we confirm earlier indications that LAEs have a much higher [OIII]/[OII] line ratio than is seen in similar redshift LBGs. By comparing their distributions on a [OIII]/[OII] versus R23 diagram, we demonstrate that this difference cannot arise solely because of their lower metallicities but most likely is due to a harder ionizing spectrum. Using measures of H-beta and recombination theory, we demonstrate, for a subset of our LAEs, that xi_ion - the number of Lyman continuum photons per UV luminosity - is indeed 0.2-0.5 dex larger than for typical LBGs at similar redshifts. Using photoionization models we estimate the effect this would have on both [OIII]/[OII] and R23 and conclude such a hard spectrum can only partially explain such intense line emission. The additional possibility is that such a large [OIII]/[OII] ratio is in part due to density rather than ionization bound nebular regions, which would imply a high escape fraction of ionizing photons. We discuss how further observations could confirm this possibility. Clearly LAEs with intense [OIII] emission represent a promising analog of those z>7 sources with similarly strong lines which are thought to be an important contributor to cosmic reionization.

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

We present a theoretical model aimed at explaining the IRX-$\beta$ relation for high redshift (z >5) galaxies. Recent observations (Capak+2015; Bouwens+2016) have shown that early Lyman Break Galaxies, although characterized by a large UV attenuation (e.g. flat UV beta slopes), show a striking FIR deficit, i.e. they are "infrared-dark". This marked deviation from the local IRX-beta relation can be explained by the larger molecular gas content of these systems. While dust in the diffuse ISM attains relatively high temperatures (Td = 45 K for typical size a=0.1 um; smaller grains can reach Td = 60 K), a sizable fraction of the dust mass is embedded in dense gas, and therefore remains cold. If confirmed, the FIR deficit might represent a novel, powerful indicator of the molecular content of high-z galaxies which can be used to pre-select candidates for follow-up deep CO observations. Thus, high-z CO line searches with ALMA might be much more promising than currently thought.

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

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

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

We report the first identification of a compact dense galaxy overdensity at $z\sim8$ dubbed A2744z8OD. A2744z8OD consists of eight $Y$-dropout galaxies behind Abell 2744 that is originally pinpointed by Hubble Frontier Fields studies. However, no studies have, so far, derived basic physical quantities of structure formation or made comparisons with theoretical models. We obtain a homogeneous sample of dropout galaxies at $z\sim8$ from eight field data of Hubble legacy images that are as deep as the A2744z8OD data. Using the sample, we find that a galaxy surface overdensity value of A2744z8OD is very high $\delta\simeq130$, where $\delta$ is defined by an overdensity in a small circle of $6"$ ($\simeq30$ physical kpc) radius. Because there is no such a large $\delta$ value reported for high-$z$ overdensities to date, A2744z8OD is a system clearly different from those found in previous high-$z$ overdensity studies. In the galaxy+structure formation models of Henriques et al. (2015), there exist a very similar overdensity, Modelz8OD, that is made of eight model dropout galaxies at $z\sim8$ in a $6"$-radius circle. Modelz8OD is a progenitor of a today's $10^{14}M_\odot$ cluster, and more than a half of the seven Modelz8OD galaxies are merged into the brightest cluster galaxy (BCG) of the cluster. If Modelz8OD is a counterpart of A2744z8OD, the models suggest that A2744z8OD would be a part of a forming cluster core of a today's $10^{14}M_\odot$ cluster that started star formation at $z>12$.

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

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

We present 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 new results on the evolution of the cosmic star formation rate as a function of stellar mass in the SXDS-UDS field. We make use of narrow-band selected emission line galaxies in four redshift slices between z = 1.46 and z = 0.63, and compute stellar masses by fitting a series of templates to recreate each galaxy's star formation history. We determine mass-binned luminosity functions in each redshift slice, and derive the star formation rate density (rhoSFR) as a function of mass using the [OIII] or [OII] emission lines. We calculate dust extinction and metallicity as a function of stellar mass, and investigate the effect of these corrections on the shape of the overall rhoSFR(M). We find that both these corrections are crucial for determining the shape of the rhoSFR(M), and its evolution with redshift. The fully corrected rhoSFR(M) is a relatively flat distribution, with the normalisation moving towards lower values of rhoSFR with increasing cosmic time/decreasing redshift, and requiring star formation to be truncated across all masses studied here. The peak of rhoSFR(M) is found in the 10^10.5<Msun<10^11.0 mass bin at z = 1.46. In the lower redshift slices the location of the peak is less certain, however low mass galaxies in the range 10^7.0<Msun<10^8.0 play an important part in the overall rhoSFR(M) out to at least z ~ 1.2.

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

We report the multi-wavelength identification of the X-ray sources found in the Subaru-XMM-Newton Deep Survey (SXDS) using deep imaging data covering the wavelength range between the far-UV to the mid-IR. We select a primary counterpart of each X-ray source by applying the likelihood ratio method to R-band, 3.6micron, near-UV, and 24micron source catalogs as well as matching catalogs of AGN candidates selected in 1.4GHz radio and i'-band variability surveys. Once candidates of Galactic stars, ultra-luminous X-ray sources in a nearby galaxy, and clusters of galaxies are removed there are 896 AGN candidates in the sample. We conduct spectroscopic observations of the primary counterparts with multi-object spectrographs in the optical and NIR; 65\% of the X-ray AGN candidates are spectroscopically-identified. For the remaining X-ray AGN candidates, we evaluate their photometric redshift with photometric data in 15 bands. Utilising the multi-wavelength photometric data of the large sample of X-ray selected AGNs, we evaluate the stellar masses, M*, of the host galaxies of the narrow-line AGNs. The distribution of the stellar mass is remarkably constant from z=0.1 to 4.0. The relation between M* and 2--10 keV luminosity can be explained with strong cosmological evolution of the relationship between the black hole mass and M*. We also evaluate the scatter of the UV-MIR spectral energy distribution (SED) of the X-ray AGNs as a function of X-ray luminosity and absorption to the nucleus. The scatter is compared with galaxies which have redshift and stellar mass distribution matched with the X-ray AGN. The UV-NIR SEDs of obscured X-ray AGNs are similar to those of the galaxies in the matched sample. In the NIR-MIR range, the median SEDs of X-ray AGNs are redder, but the scatter of the SEDs of the X-ray AGN broadly overlaps that of the galaxies in the matched sample.

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 deep 10h VLT/XSHOOTER spectroscopy for an extraordinarily luminous and extended Lya emitter at z=6.595 referred to as Himiko and first discussed by Ouchi et al. (2009), with the purpose of constraining the mechanisms powering its strong emission. Complementary to the spectrum, we discuss NIR imaging data from the CANDELS survey. We find neither for HeII nor any metal line a significant excess, with 3 sigma upper limits of 6.8, 3.1, and 5.8x10^{-18} erg/s/cm^2 for CIV $\lambda$1549, HeII $\lambda$1640, CIII] $\lambda$1909, respectively, assuming apertures with 200 km/s widths and offset by -250 km/s w.r.t to the peak Lya redshift. These limits provide strong evidence that an AGN is not a major contribution to Himiko's Lya flux. Strong conclusions about the presence of PopIII star-formation or gravitational cooling radiation are not possible based on the obtained HeII upper limit. Our Lya spectrum confirms both spatial extent and flux (8.8+/-0.5x10^{-17} erg/s/cm^2) of previous measurements. In addition, we can unambiguously exclude any remaining chance of it being a lower redshift interloper by significantly detecting a continuum redwards of Lya, while being undetected bluewards.

Ly$\alpha$ photons scattered by neutral hydrogen atoms in the circumgalactic media or produced in the halos of star-forming galaxies are expected to lead to extended Ly$\alpha$ emission around galaxies. Such low surface brightness Ly$\alpha$ halos (LAHs) have been detected by stacking Ly$\alpha$ images of high-redshift star-forming galaxies. We study the origin of LAHs by performing radiative transfer modeling of nine $z=3.1$ Lyman-Alpha Emitters (LAEs) in a high resolution hydrodynamic cosmological galaxy formation simulation. We develop a method of computing the mean Ly$\alpha$ surface brightness profile of each LAE by effectively integrating over many different observing directions. Without adjusting any parameters, our model yields an average Ly$\alpha$ surface brightness profile in remarkable agreement with observations. We find that observed LAHs cannot be accounted for solely by photons originating from the central LAE and scattered to large radii by hydrogen atoms in the circumgalactic gas. Instead, Ly$\alpha$ emission from regions in the outer halo is primarily responsible for producing the extended LAHs seen in observations, which potentially includes both star-forming and cooling radiation. With the limit on the star formation contribution set by the ultra-violet (UV) halo measurement, we find that cooling radiation can play an important role in forming the extended LAHs. We discuss the implications and caveats of such a picture.

The flow of baryons to and from a galaxy, which is fundamental for galaxy formation and evolution, can be studied with galaxy-metal absorption system pairs. Our search for galaxies around CIV absorption systems at $z\sim5.7$ showed an excess of photometric Lyman-$\alpha$ emitter (LAE) candidates in the fields J1030+0524 and J1137+3549. Here we present spectroscopic follow-up of 33 LAEs in both fields. In the first field, three out of the five LAEs within 10$h^{{-}1}$ projected comoving Mpc from the CIV system are within $\pm500$ km s$^{{-}1}$ from the absorption at $z_{\text{CIV}}=5.7242\pm0.0001$. The closest candidate (LAE 103027+052419) is robustly confirmed at $212.8^{+14}_{-0.4}h^{-1}$ physical kpc from the CIV system. In the second field, the LAE sample is selected at a lower redshift ($\Delta z\sim0.04$) than the CIV absorption system as a result of the filter transmission and, thus, do not trace its environment. The observed properties of LAE 103027+052419 indicate that it is near the high-mass end of the LAE distribution, probably having a large HI column density and large-scale outflows. Therefore, our results suggest that the CIV system is likely produced by a star-forming galaxy which has been injecting metals into the intergalactic medium since $z>6$. Thus, the CIV system is either produced by LAE 103027+052419, implying that outflows can enrich larger volumes at $z>6$ than at $z\sim3.5$, or an undetected dwarf galaxy. In either case, CIV systems like this one trace the ionized intergalactic medium at the end of cosmic hydrogen reionization and may trace the sources of the ionizing flux density.

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

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

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

Using a sample of 299 Ha-selected galaxies at z~0.8, we study the relationship between galaxy stellar mass, gas-phase metallicity, and star formation rate (SFR), and compare to previous results. We use deep optical spectra obtained with the IMACS spectrograph at the Magellan telescope to measure strong oxygen lines. We combine these spectra and metallicities with (1) rest-frame UV-to-optical imaging, which allows us to determine stellar masses and dust attenuation corrections, and (2) Ha narrowband imaging, which provides a robust measure of the instantaneous SFR. Our sample spans stellar masses of 10^9 to 6*10^11 solar masses, SFRs of 0.4 to 270 solar masses per year, and metal abundances of 12+log(O/H)~8.3-9.1 (~0.4-2.6 solar metallicity). The correlations that we find between the Ha-based SFR and stellar mass (i.e., the star-forming "main sequence"), and between the stellar mass and metallicity, are both consistent with previous z~1 studies of star-forming galaxies. We then study the relationship between the three properties using various plane-fitting techniques (Lara-Lopez et al.) and a curve-fitting projection (Mannucci et al.). In all cases, we exclude strong dependence of the M-Z relation on SFR, but are unable to distinguish between moderate and no dependence. Our results are consistent with previous mass-metallicity-SFR studies. We check whether dataset limitations may obscure a strong dependence on the SFR by using mock samples drawn from the SDSS. These experiments reveal that the adopted signal-to-noise cuts may have a significant effect on the measured dependence. Further work is needed to investigate these results, and to test whether a "fundamental metallicity relation" or a "fundamental plane" describes star-forming galaxies across cosmic time.

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

We present the statistics of faint submillimeter/millimeter galaxies (SMGs) and serendipitous detections of a submillimeter/millimeter line emitter (SLE) with no multi-wavelength continuum counterpart revealed by the deep ALMA observations. We identify faint SMGs with flux densities of 0.1-1.0 mJy in the deep Band 6 and Band 7 maps of 10 independent fields that reduce cosmic variance effects. The differential number counts at 1.2 mm are found to increase with decreasing flux density down to 0.1 mJy. Our number counts indicate that the faint (0.1-1.0 mJy, or SFR_IR ~ 30-300 Msun/yr) SMGs contribute nearly a half of the extragalactic background light (EBL), while the remaining half of the EBL is mostly contributed by very faint sources with flux densities of <0.1 mJy (SFR_IR <~ 30 Msun/yr). We conduct counts-in-cells analysis with the multifield ALMA data for the faint SMGs, and obtain a coarse estimate of galaxy bias, b_g <4. The galaxy bias suggests that the dark halo masses of the faint SMGs are <~ 7x10^12 Msun, which is smaller than those of bright (>1 mJy) SMGs, but consistent with abundant high-z star-forming populations such as sBzKs, LBGs, and LAEs. Finally, we report the serendipitous detection of SLE-1 with continuum counterparts neither in our 1.2 mm-band nor multi-wavelength images including ultra deep HST/WFC3 and Spitzer data. The SLE has a significant line at 249.9 GHz with a signal-to-noise ratio of 7.1. If the SLE is not a spurious source made by unknown systematic noise of ALMA, the strong upper limits of our multi-wavelength data suggest that the SLE would be a faint galaxy at z >~ 6.

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 report our Keck/MOSFIRE and Magellan/LDSS3 spectroscopy for an [OII] Blob, OIIB10, that is a high-$z$ galaxy with spatially extended [OII]$\lambda\lambda3726,3729$ emission over 30 kpc recently identified by a Subaru large-area narrowband survey. The systemic redshift of OIIB10 is $z=1.18$ securely determined with [OIII]$\lambda\lambda4959,5007$ and H$\beta$ emission lines. We identify FeII$\lambda$2587 and MgII$\lambda\lambda$2796,2804 absorption lines blueshifted from the systemic redshift by $80\pm50$ and $260\pm40$ km s$^{-1}$, respectively, which indicate gas outflow from OIIB10 with the velocity of $\sim 80-260$ km s$^{-1}$. This outflow velocity is comparable with the escape velocity, $250\pm140$ km s$^{-1}$, estimated under the assumption of a singular isothermal halo potential profile. Some fraction of the outflowing gas could escape from the halo of OIIB10, suppressing OIIB10's star-formation activity. We estimate a mass loading factor, $\eta$, that is a ratio of mass outflow rate to star-formation rate, and obtain $\eta>0.8\pm 0.1$ which is relatively high compared with low-$z$ starbursts including U/LIRGs and AGNs. The major energy source of the outflow is unclear with the available data. Although no signature of AGN is found in the X-ray data, OIIB10 falls in the AGN/star-forming composite region in the line diagnostic diagrams. It is possible that the outflow is powered by star formation and a type-2 AGN with narrow FWHM emission line widths of $70-130$ km s$^{-1}$. This is the first detailed spectroscopic study of oxygen-line blobs, which includes the analyses of the escape velocity, the mass loading factor, and the presence of an AGN, and a significant step to understanding the nature of oxygen-line blobs and the relation with gas outflow and star-formation quenching at high redshift.

Core-collapse supernovae (SNe) are believed to be the first significant source of dust in the Universe. Such SNe are expected to be the main dust producers in young high-redshift Lyman $\alpha$ emitters (LAEs) given their young ages, providing an excellent testbed of SN dust formation models during the early stages of galaxy evolution. We focus on the dust enrichment of a specific, luminous LAE (Himiko, $z\simeq 6.6$) for which a stringent upper limit of $52.1~\mu$Jy ($3\sigma$) has recently been obtained from ALMA continuum observations at 1.2 mm. We predict its submillimetre dust emission using detailed models that follow SN dust enrichment and destruction and the equilibrium dust temperature, and obtain a plausible upper limit to the dust mass produced by a single SN: $m_\mathrm{d,SN} < 0.15$--0.45 M$_\odot$, depending on the adopted dust optical properties. These upper limits are smaller than the dust mass deduced for SN 1987A and that predicted by dust condensation theories, implying that dust produced in SNe are likely to be subject to reverse shock destruction before being injected into the interstellar medium. Finally, we provide a recipe for deriving $m_\mathrm{d,SN}$ from submillimetre observations of young, metal poor objects wherein condensation in SN ejecta is the dominant dust formation channel.

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

We present a systematic study for ionization state of inter-stellar medium in galaxies at z=0-3 with ~140,000 SDSS galaxies and 108 intermediate to high redshift galaxies from the literature, using an ionization-parameter sensitive line ratio of [OIII]5007/[OII]3727 and photoionization models. We confirm that z=2-3 galaxies show an [OIII]/[OII] ratio significantly higher than a typical star-forming galaxy of SDSS by a factor of >~10, and the photoionization models reveal that these high-z galaxies have an ionization parameter of log(qion/cm s^{-1})~7.6-9.0, a factor of ~4-10 higher than local galaxies. For galaxies at any redshift, we identify a correlation between the [OIII]/[OII] ratio and galaxy global properties of star-formation rate (SFR), stellar mass (M_star), and metallicity (Z). We extend the fundamental metallicity relation (FMR) and develop the fundamental ionization relation (FIR), a four-dimensional relation of ionization parameter, SFR, M_star, and Z. The intermediate and high-z galaxies up to z~3 follow the FIR defined with the local galaxies, in contrast with the FMR whose possible evolution from z~2 to 3 is reported. We find that the FMR evolution of z~2-3 appears, if one omits ionization parameter differences, and that the FMR evolution does not exist for an average metallicity solution of z~3 galaxies with a high-ionization parameter. Interestingly, all of two local Lyman-continuum emitting galaxies (LyC leakers) have a high [OIII]/[OII] ratio, indicating a positive correlation between [OIII]/[OII] and ionizing photon escape fraction (fesc), which is successfully explained by our photoionization models. Because [OIII]/[OII] ratios of z~2-3 galaxies, especially Ly-alpha emitters (LAEs), are comparable to, or higher than, those of the local LyC leakers, these high-z galaxies are candidates of Lyman-continuum emitting objects. (abridged)

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

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

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

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

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

High-redshift quasars are believed to reside in massive halos in the early universe and should therefore be located in fields with overdensities of galaxies, which are thought to evolve into galaxy clusters seen in the local universe. However, despite many efforts, the relationship between galaxy overdensities and z~6 quasars is ambiguous. This can possibly be attributed to the difficulty of finding galaxies with accurate redshifts in the vicinity of z~6 quasars. So far, overdensity searches around z~6 quasars have been based on studies of Lyman break galaxies (LBGs), which probe a redshift range of Delta z ~ 1. This range is large enough to select galaxies that may not be physically related to the quasar. We use deep narrow- and broadband imaging to study the environment of the z=5.72 quasar ULAS J0203+0012. The redshift range probed by our narrow-band selection of Lyman alpha emitters (LAEs) is Delta z ~ 0.1, which is significantly narrower than the LBG searches. This is the first time that LAEs were searched for near a z~6 quasar, in an effort to provide clues about the environments of quasars at the end of the epoch of reionization. We find no enhancement of LAEs in the surroundings of ULAS J0203+0012 in comparison with blank fields. We explore different explanations and interpretations for this non-detection of a galaxy overdensity, including that (1) the strong ionization from the quasar may prevent galaxy formation in its immediate vicinity and (2) high-redshift quasars may not reside in the center of the most massive dark matter halos.

The Subaru Prime Focus Spectrograph (PFS) is a massively-multiplexed fiber-fed optical and near-infrared 3-arm spectrograph (N_fiber=2400, 380<lambda<1260nm, 1.3 degree diameter FoV), offering unique opportunities in survey astronomy. Here we summarize the science case feasible for a survey of Subaru 300 nights. We describe plans to constrain the nature of dark energy via a survey of emission line galaxies spanning a comoving volume of 9.3 (Gpc/h)^3 in the redshift range 0.8<z<2.4. In each of 6 redshift bins, the cosmological distances will be measured to 3% precision via BAO, and redshift-space distortions will be used to constrain structure growth to 6% precision. In the GA program, radial velocities and chemical abundances of stars in the Milky Way and M31 will be used to infer the past assembly histories of spiral galaxies and the structure of their dark matter halos. Data will be secured for 10^6 stars in the Galactic thick-disk, halo and tidal streams as faint as V~22, including stars with V < 20 to complement the goals of the Gaia mission. A medium-resolution mode with R = 5000 to be implemented in the red arm will allow the measurement of multiple alpha-element abundances and more precise velocities for Galactic stars, elucidating the detailed chemo-dynamical structure and evolution of each of the main stellar components of the Milky Way Galaxy and of its dwarf spheroidal galaxies. For the extragalactic program, our simulations suggest the wide avelength range will be powerful in probing the galaxy population and its clustering over a wide redshift range. We propose to conduct a color-selected survey of 1<z<2 galaxies and AGN over 16 deg^2 to J~23.4, yielding a fair sample of galaxies with stellar masses above ~10^{10}Ms at z~2. A two-tiered survey of higher redshift LBGs and LAEs will quantify the properties of early systems close to the reionization epoch.

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

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

We present new results on the cosmic star formation history in the SXDS-UDS field out to z=1.6. We compile narrow-band data from the Subaru Telescope and the Visible and Infrared Survey Telescope for Astronomy (VISTA) in conjunction with broad-band data from the SXDS and UDS, to make a selection of 5725 emission-line galaxies in 12 redshift slices, spanning 10 Gyr of cosmic time. We determine photometric redshifts for the sample using 11-band photometry, and use a spectroscopically confirmed subset to fine tune the resultant redshift distribution. We use the maximum-likelihood technique to determine luminosity functions in each redshift slice and model the selection effects inherent in any narrow-band selection statistically, to obviate the retrospective corrections ordinarily required. The deep narrow-band data are sensitive to very low star formation rates (SFRs), and allow an accurate evaluation of the faint end slope of the Schechter function, alpha. We find that alpha is particularly sensitive to the assumed faintest broadband magnitude of a galaxy capable of hosting an emission line, and propose that this limit should be empirically motivated. For this analysis we base our threshold on the limiting observed equivalent widths of emission lines in the local Universe. We compute the characteristic SFR of galaxies in each redshift slice, and the integrated SFR density, rho_{SFR}. We find our results to be in good agreement with the literature and parametrize the evolution of the SFR density as rho_{SFR} proportional to (1 + z)^{4.58} confirming a steep decline in star formation activity since z = 1.6.

We use the new ultra-deep, near-infrared imaging of the Hubble Ultra-Deep Field (HUDF) provided by our UDF12 HST WFC3/IR campaign to explore the rest-frame UV properties of galaxies at redshifts z > 6.5. We present the first unbiased measurement of the average UV power-law index, beta, for faint galaxies at z ~ 7, the first meaningful measurements of beta at z ~ 8, and tentative estimates for a new sample of galaxies at z ~ 9. Utilising galaxy selection in the new F140W imaging to minimize colour bias, and applying both colour and power-law estimators of beta, we find beta = -2.1 (+/-0.2) at z ~ 7 for galaxies with M_UV ~ -18. This means that the faintest galaxies uncovered at this epoch have, on average, UV colours no more extreme than those displayed by the bluest star-forming galaxies at low redshift. At z ~ 8 we find a similar value, beta = -1.9 (+/-0.3). At z ~ 9, we find beta = -1.8 (+/-0.6), essentially unchanged from z ~ 6 - 7 (albeit highly uncertain). Finally, we show that there is as yet no evidence for a significant intrinsic scatter in beta within our new, robust z ~ 7 galaxy sample. Our results are most easily explained by a population of steadily star-forming galaxies with either ~ solar metallicity and zero dust, or moderately sub-solar (~ 10-20%) metallicity with modest dust obscuration (A_V ~ 0.1-0.2). This latter interpretation is consistent with the predictions of a state-of-the-art galaxy-formation simulation, which also suggests that a significant population of very-low metallicity, dust-free galaxies with beta ~ -2.5 may not emerge until M_UV > -16, a regime likely to remain inaccessible until the James Webb Space Telescope.

We present a new determination of the UV galaxy luminosity function (LF) at redshift z ~ 7 and z ~ 8, and a first estimate at z ~ 9. An accurate determination of the form and evolution of the LF at high z is crucial for improving our knowledge of early galaxy evolution and cosmic reionization. Our analysis exploits fully the new, deepest WFC3/IR imaging from our HST UDF12 campaign, and includes a new, consistent analysis of all appropriate, shallower/wider-area HST data. Our new measurement of the evolving LF at z ~ 7-8 is based on a final catalogue of ~600 galaxies, and involves a step-wise maximum likelihood determination based on the redshift probability distribution for each object; this makes full use of the 11-band imaging now available in the HUDF, including the new UDF12 F140W data, and the deep Spitzer IRAC imaging. The final result is a determination of the z ~ 7 LF extending down to M_UV = -16.75, and the z ~ 8 LF down to M_UV = -17.00. Fitting a Schechter function, we find M* = -19.90 (+0.23/-0.28), log phi* = -2.96 (+0.18/-0.23), and a faint-end slope alpha=-1.90 (+0.14/-0.15) at z~7, and M* = -20.12 (+0.37/-0.48), log phi* = -3.35 (+0.28/-0.47), alpha=-2.02 (+0.22/-0.23) at z~8. These results strengthen suggestions that the evolution at z > 7 is more akin to `density evolution' than the apparent `luminosity evolution' seen at z ~ 5-7. We also provide the first meaningful information on the LF at z ~ 9, explore alternative extrapolations to higher z, and consider the implications for the evolution of UV luminosity density. Finally, we provide catalogues (including z_phot, M_UV and all photometry) for the 100 most robust z~6.5-11.9 galaxies in the HUDF used in this analysis. We discuss our results in the context of earlier work and the results of an independent analysis of the UDF12 data based on colour-colour selection (Schenker et al. 2013).

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

Understanding cosmic reionization requires the identification and characterization of early sources of hydrogen-ionizing photons. The 2012 Hubble Ultra Deep Field (UDF12) campaign has acquired the deepest infrared images with the Wide Field Camera 3 aboard Hubble Space Telescope and, for the first time, systematically explored the galaxy population deep into the era when cosmic microwave background (CMB) data indicates reionization was underway. The UDF12 campaign thus provides the best constraints to date on the abundance, luminosity distribution, and spectral properties of early star-forming galaxies. We synthesize the new UDF12 results with the most recent constraints from CMB observations to infer redshift-dependent ultraviolet (UV) luminosity densities, reionization histories, and electron scattering optical depth evolution consistent with the available data. Under reasonable assumptions about the escape fraction of hydrogen ionizing photons and the intergalactic medium clumping factor, we find that to fully reionize the universe by redshift z~6 the population of star-forming galaxies at redshifts z~7-9 likely must extend in luminosity below the UDF12 limits to absolute UV magnitudes of M_UV\sim -13 or fainter. Moreover, low levels of star formation extending to redshifts z~15-25, as suggested by the normal UV colors of z\simeq7-8 galaxies and the smooth decline in abundance with redshift observed by UDF12 to z\simeq10, are additionally likely required to reproduce the optical depth to electron scattering inferred from CMB observations.

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

We present a catalog of high redshift star-forming galaxies selected to lie within the redshift range z ~ 7-8 using the Ultra Deep Field 2012 (UDF12), the deepest near-infrared (near-IR) exposures yet taken with the Hubble Space Telescope. As a result of the increased near-infrared exposure time compared to previous HST imaging in this field, we probe 0.65 (0.25) mag fainter in absolute UV magnitude, at z ~ 7 (8), which increases confidence in a measurement of the faint end slope of the galaxy luminosity function. Through a 0.7 mag deeper limit in the key F105W filter that encompasses or lies just longward of the Lyman break, we also achieve a much-refined color-color selection that balances high redshift completeness and a low expected contamination fraction. We improve the number of drop-out selected UDF sources to 47 at z ~ 7 and 27 at z ~ 8. Incorporating brighter archival and ground-based samples, we measure the z ~ 7 UV luminosity function to an absolute magnitude limit of M_UV = -17 and find a faint end Schechter slope of \alpha = -1.87+/- 0.18. Using a similar color-color selection at z ~ 8 that takes account of our newly-added imaging in the F140W filter, and incorporating archival data from the HIPPIES and BoRG campaigns, we provide a robust estimate of the faint end slope at z ~ 8, \alpha = -1.94 +/- 0.23. We briefly discuss our results in the context of earlier work and that derived using the same UDF12 data but with an independent photometric redshift technique (McLure et al 2012).

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

We present the 2012 Hubble Ultra Deep Field campaign (UDF12), a large 128-orbit Cycle 19 \HST\ program aimed at extending previous WFC3/IR observations of the UDF by quadrupling the exposure time in the F105W filter, imaging in an additional F140W filter, and extending the F160W exposure time by 50%. The principal scientific goal of this project is to determine whether galaxies reionized the universe; our observations are designed to provide a robust determination of the star formation density at $z$$\,\gtrsim\,$8, improve measurements of the ultraviolet continuum slope at $z$$\,\sim\,7\,-\,$8, facilitate the construction of new samples of $z$$\,\sim\,9\,-\,$10 candidates, and enable the detection of sources up to $z$$\,\sim\,$12. For this project we committed to combining these and other WFC3/IR imaging observations of the UDF area into a single homogeneous dataset, to provide the deepest near-infrared observations of the sky currently achievable. In this paper we present the observational overview of the project, motivated by its scientific goals, and describe the procedures used in reducing the data as well as the final products that are produced. We have used the most up up-to-date methods for calibrating and combining the images, in particular paying attention to correcting several instrumental effects. We release the full combined mosaics, comprising a single, unified set of mosaics of the UDF, providing the deepest near-infrared blank-field view of the universe obtained to date, reaching magnitudes as deep as AB$\,\sim\,$30 in the near-infrared, and yielding a legacy dataset on this field of lasting scientific value to the community.

We present the results of the deepest search to date for star-forming galaxies beyond a redshift z~8.5 utilizing a new sequence of near-infrared Wide Field Camera 3 images of the Hubble Ultra Deep Field. This `UDF12' campaign completed in September 2012 doubles the earlier exposures with WFC3/IR in this field and quadruples the exposure in the key F105W filter used to locate such distant galaxies. Combined with additional imaging in the F140W filter, the fidelity of high redshift candidates is greatly improved. Using spectral energy distribution fitting techniques on objects selected from a deep multi-band near-infrared stack we find 7 promising z>8.5 candidates. As none of the previously claimed UDF candidates with 8.5<z<10 is confirmed by our deeper multi-band imaging, our campaign has transformed the measured abundance of galaxies in this redshift range. Although we recover the candidate UDFj-39546284 (previously proposed at z=10.3), it is undetected in the newly added F140W image, implying it lies at z=11.9 or is an intense emission line galaxy at z~2.4. Although no physically-plausible model can explain the required line intensity given the lack of Lyman alpha or broad-band UV signal, without an infrared spectrum we cannot rule out an exotic interloper. Regardless, our robust z ~ 8.5 - 10 sample demonstrates a luminosity density that continues the smooth decline observed over 6 < z < 8. Such continuity has important implications for models of cosmic reionization and future searches for z>10 galaxies with JWST.

In order to investigate the growth of super-massive black holes (SMBHs), we construct the black hole mass function (BHMF) and Eddington ratio distribution function (ERDF) of X-ray-selected broad-line AGNs at z~1.4 in the Subaru XMM-Newton Deep Survey field. In this redshift range, a significant part of the accretion growth of SMBHs is thought to be taking place. Black hole masses of X-ray-selected broad-line AGNs are estimated using the width of the broad MgII line and the 3000A monochromatic luminosity. We supplement the MgII FWHM values with the Ha FWHM obtained from our NIR spectroscopic survey. Using the black hole masses of broad-line AGNs at redshifts between 1.18 and 1.68, the binned broad-line AGN BHMF and ERDF are calculated using the Vmax method. To properly account for selection effects that impact the binned estimates, we derive the corrected broad-line AGN BHMF and ERDF by applying the Maximum Likelihood method, assuming that the ERDF is constant regardless of the black hole mass. We do not correct for the non-negligible uncertainties in virial BH mass estimates. If we compare the corrected broad-line AGN BHMF with that in the local Universe, the corrected BHMF at z~1.4 has a higher number density above 10^8 Msolar but a lower number density below that mass range. The evolution may be indicative of a down-sizing trend of accretion activity among the SMBH population. The evolution of broad-line AGN ERDF from z=1.4 to 0 indicates that the fraction of broad-line AGNs with accretion rate close to the Eddington-limit is higher at higher redshifts.

We present the Spitzer Extragalactic Representative Volume Survey (SERVS), an 18 square degrees medium-deep survey at 3.6 and 4.5 microns with the post-cryogenic Spitzer Space Telescope to ~2 microJy (AB=23.1) depth of five highly observed astronomical fields (ELAIS-N1, ELAIS-S1, Lockman Hole, Chandra Deep Field South and XMM-LSS). SERVS is designed to enable the study of galaxy evolution as a function of environment from z~5 to the present day, and is the first extragalactic survey both large enough and deep enough to put rare objects such as luminous quasars and galaxy clusters at z>1 into their cosmological context. SERVS is designed to overlap with several key surveys at optical, near- through far-infrared, submillimeter and radio wavelengths to provide an unprecedented view of the formation and evolution of massive galaxies. In this paper, we discuss the SERVS survey design, the data processing flow from image reduction and mosaicing to catalogs, as well as coverage of ancillary data from other surveys in the SERVS fields. We also highlight a variety of early science results from the survey.

We present measurements of the dust attenuation of H\alpha-selected emission-line galaxies at z=0.8 from the NewH\alpha\ narrowband survey. The analysis is based on deep follow-up spectroscopy with Magellan/IMACS, which captures the strong rest-frame optical emission lines from [OII] \lambda 3727 to [OIII] \lambda 5007. The spectroscopic sample used in this analysis consists of 341 confirmed H\alpha\ emitters. We place constraints on the AGN fraction using diagnostics which can be applied at intermediate redshift. We find that at least 5% of the objects in our spectroscopic sample can be classified as AGN and 2% are composite, i.e. powered by a combination of star-formation and AGN activity. We measure the dust attenuation for individual objects from the ratios of the higher order Balmer lines. The H\beta\ and H\gamma\ pair of lines is detected with S/N>5 in 55 individual objects and the H\beta\ and H\delta\ pair is detected in 50 individual objects. We also create stacked spectra to probe the attenuation in objects without individual detections. The median attenuation at H\alpha\ based on the objects with individually detected lines is A(H\alpha)=0.9+-1.0 magnitudes, in good agreement with the attenuation found in local samples of star-forming galaxies. We find that the z=0.8 galaxies occupy a similar locus of attenuation as a function of magnitude, mass and SFR as a comparison sample drawn from the SDSS DR4. Both the results from the individual z=0.8 galaxies and from the stacked spectra show consistency with the mass -- attenuation and SFR -- attenuation relations found in the local Universe, indicating that these relations are also applicable at intermediate redshift.

We present first results from a narrowband imaging program for intermediate redshift emission-line galaxies using the newly commissioned FourStar infrared camera at the 6.5m Magellan telescope. To enable prompt identification of H\alpha\ emitters, a pair of custom 1% filters, which sample low-airglow atmospheric windows at 1.19 \mu m and 2.10 \mu m, is used to detect both H\alpha\ and [OII]\lambda 3727 emission from the same redshift volume at z=2.2. Initial observations are taken over a 130 arcmin^2 area in the CANDELS-COSMOS field. The exquisite image quality resulting from the combination of the instrument, telescope, and standard site conditions (~0.55" FWHM) allows the 1.19 \mu m and 2.10 \mu m data to probe 3\sigma\ emission-line depths down to 1.0e-17 erg/s/cm^2 and 1.2e-17 erg/s/cm^2 respectively, in less than 10 hours of integration time in each narrowband. For H\alpha\ at z=0.8 and z=2.2, these fluxes correspond to observed star formation rates of ~0.3 and ~4 Msun/yr respectively. We find 122 sources with a 1.19 \mu m excess, and 136 with a 2.10 \mu m excess, 41 of which show an excess in both bands. The dual narrowband technique, as implemented here, is estimated to identify about >80% of z=2.2 H\alpha\ emitters in the narrowband excess population. With the most secure such sample obtained to-date, we compute constraints on the faint-end slope of the z=2.2 H\alpha\ luminosity function. These "narrow-deep" FourStar observations have been obtained as part of the larger NewH\alpha\ Survey, which will combine the data with "wide-shallow" imaging through a similar narrowband filter pair with NEWFIRM at the KPNO/CTIO 4m telescopes, to enable study of both luminous (but rare) and faint emission-line galaxies in the intermediate redshift universe. [Abridged]

Using stacks of Ly-a images of 2128 Ly-a emitters (LAEs) and 24 protocluster UV-selected galaxies (LBGs) at z=3.1, we examine the surface brightness profiles of Ly-a haloes around high-z galaxies as a function of environment and UV luminosity. We find that the slopes of the Ly-a radial profiles become flatter as the Mpc-scale LAE surface densities increase, but they are almost independent of the central UV luminosities. The characteristic exponential scale lengths of the Ly-a haloes appear to be proportional to the square of the LAE surface densities (r(Lya) \propto Sigma(LAE)^2). Including the diffuse, extended Ly-a haloes, the rest-frame Ly-a equivalent width of the LAEs in the densest regions approaches EW_0(Lya) ~ 200 A, the maximum value expected for young (< 10^7 yr) galaxies. This suggests that Ly-a photons formed via shock compression by gas outflows or cooling radiation by gravitational gas inflows may partly contribute to illuminate the Ly-a haloes; however, most of their Ly-a luminosity can be explained by photo-ionisation by ionising photons or scattering of Ly-a photons produced in HII regions in and around the central galaxies. Regardless of the source of Ly-a photons, if the Ly-a haloes trace the overall gaseous structure following the dark matter distributions, it is not surprising that the Ly-a spatial extents depend more strongly on the surrounding Mpc-scale environment than on the activities of the central galaxies.

We present a search for the [CII] 158micron fine structure line (a main cooling line of the interstellar medium) and the underlying far-infrared (FIR) continuum in three high-redshift (6.6<z<8.2) star-forming galaxies using the IRAM Plateau de Bure interferometer. We targeted two Lyman-Alpha-selected galaxies (Lyman-Alpha-Emitters, LAEs) with moderate UV-based star formation rates (SFR~20 M_sun/yr; Himiko at z=6.6 and IOK-1 at z=7.0) and a Gamma Ray Burst (GRB) host galaxy (GRB 090423 at z~8.2). Based on our 3 sigma rest-frame FIR continuum limits, previous (rest-frame) UV continuum measurements and spectral energy distribution (SED) fitting, we rule out SED shapes similar to highly obscured galaxies (e.g. Arp220, M82) and less extreme dust-rich nearby spiral galaxies (e.g. M51) for the LAEs. Conservatively assuming a SED shape typical of local spiral galaxies we derive upper limits for the FIR-based star formation rates (SFRs) of ~70 M_sun/yr, ~50 M_sun/yr and ~40 M_sun/yr for Himiko, IOK-1 and GRB 090423, respectively. For the LAEs these limits are only a factor ~3 higher than the published UV-based SFRs (uncorrected for extinction). This indicates that the dust obscuration in the z>6 LAEs studied here is lower by a factor of a few than what has recently been found in some LAEs at lower redshift (2<z<3.5) with similar UV-based SFRs. A low obscuration in our z>6 LAE sample is consistent with recent rest-frame UV studies of z~7 Lyman-Break-Galaxies (LBGs).

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

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

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

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

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

We present results of a survey for giant Ly-alpha nebulae (LABs) at z=3 with Subaru/Suprime-Cam. We obtained Ly-alpha imaging at z=3.09+-0.03 around the SSA22 protocluster and in several blank fields. The total survey area is 2.1 square degrees, corresponding to a comoving volume of 1.6 x 10^6 Mpc^3. Using a uniform detection threshold of 1.4 x 10^{-18} erg s^{-1} cm^{-2} arcsec^{-2} for the Ly-alpha images, we construct a sample of 14 LAB candidates with major-axis diameters larger than 100 kpc, including five previously known blobs and two known quasars. This survey triples the number of known LABs over 100 kpc. The giant LAB sample shows a possible "morphology-density relation": filamentary LABs reside in average density environments as derived from compact Ly-alpha emitters, while circular LABs reside in both average density and overdense environments. Although it is hard to examine the formation mechanisms of LABs only from the Ly-alpha morphologies, more filamentary LABs may relate to cold gas accretion from the surrounding inter-galactic medium (IGM) and more circular LABs may relate to large-scale gas outflows, which are driven by intense starbursts and/or by AGN activities. Our survey highlights the potential usefulness of giant LABs to investigate the interactions between galaxies and the surrounding IGM from the field to overdense environments at high-redshift.

As Lyman-alpha photons are scattered by neutral hydrogen, a change with redshift in the Lyman-alpha equivalent width distribution of distant galaxies offers a promising probe of the degree of ionization in the intergalactic medium and hence when cosmic reionization ended. This simple test is complicated by the fact that Lyman-alpha emission can also be affected by the evolving astrophysical details of the host galaxies. In the first paper in this series, we demonstrated both a luminosity and redshift dependent trend in the fraction of Lyman-alpha emitters seen within color-selected Lyman-break galaxies (LBGs) over the range 3<z<6; lower luminosity galaxies and those at higher redshift show an increased likelihood of strong emission. Here we present the results from much deeper 12.5 hour exposures with the Keck DEIMOS spectrograph focused primarily on LBGs at z~6 which enable us to confirm the redshift dependence of line emission more robustly and to higher redshift than was hitherto possible. We find 54+/-11% of faint z~6 LBGs show strong (W_0>25 A) emission, an increase of 1.6x from a similar sample observed at z~4. With a total sample of 74 z~6 LBGs, we determine the luminosity-dependent Lyman-alpha equivalent width distribution. Assuming continuity in these trends to the new population of z~7 sources located with the Hubble WFC3/IR camera, we predict that unless the neutral fraction rises in the intervening 200 Myr, the success rate for spectroscopic confirmation using Lyman-alpha emission should be high.

We study the properties of Lyman-alpha emitters (LAEs) and Lyman-break galaxies (LBGs) at z=3-6 using cosmological SPH simulations. We investigate two simple scenarios for explaining the observed Ly-a and rest-frame UV luminosity functions (LFs) of LAEs: (i) the "escape fraction" scenario, in which the "effective" escape fraction (including the IGM attenuation) of Ly-a photons is f_Lya ~0.1 (0.15) at z=3 (6), and (ii) the "stochastic" scenario, in which the fraction of LAEs that are turned on at z=3 (6) is \Cstoc ~0.07 (0.2) after correcting for the IGM attenuation. Our comparisons with a number of different observations suggest that the stochastic scenario is preferred over the escape fraction scenario. We find that the mean values of stellar mass, metallicity and black hole mass hosted by LAEs are all smaller in the stochastic scenario than in the escape fraction scenario. In our simulations, the galaxy stellar mass function evolves rapidly, as expected in hierarchical structure formation. However, its evolution is largely compensated by a beginning decline in the specific star formation rate, resulting in little evolution of the rest-frame UV LF from z=6 to 3. The rest-frame UV LF of both LAEs and LBGs at z=3 & 6 can be described well by the stochastic scenario provided the extinction is moderate, E(B-V) ~0.15, for both populations, although our simulation might be overpredicting the number of bright LBGs at z=6. We also discuss the correlation function and bias of LAEs. The Ly-a LFs at z=6 in a field-of-view of 0.2 deg^2 show a significantly larger scatter owing to cosmic variance relative to that in a 1 deg^2 field, and the scatter seen in the current observational estimates of the Ly-a LF can be accounted for by cosmic variance.

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

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

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

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

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

We present results from a new Keck spectroscopic survey of UV-faint LBGs in the redshift range 3<z<7. Combined with earlier Keck and published ESO VLT data, our sample contains more than 600 dropouts, offering new insight into the nature of sub-L* sources typical of those likely to dominate the cosmic reionisation process. Here we use this sample to characterise the fraction of strong Lya emitters within the continuum-selected dropouts. By quantifying how the "Lya fraction" varies with redshift, we seek to constrain changes in Lya transmission associated with reionisation. In order to distinguish the effects of reionisation from other factors which affect the Lya fraction (e.g. dust, ISM kinematics), we study the luminosity and redshift-dependence of the Lya fraction over 3<z<6, when the IGM is known to be ionised. These results reveal that low luminosity galaxies show strong Lya emission much more frequently than luminous systems, and that at fixed luminosity, the prevalence of strong Lya emission increases moderately with redshift over 3 < z < 6. Based on the correlation between blue UV slopes and strong Lya emitting galaxies in our dataset, we argue that the Lya fraction trends are governed by redshift and luminosity-dependent variations in the dust obscuration, with likely additional contributions from trends in the kinematics and covering fraction of neutral hydrogen. We find a tentative decrease in the Lya fraction at z~7 based on the limited IR spectra for candidate z~7 lensed LBGs, a result which, if confirmed with future surveys, would suggest an increase in the neutral fraction by this epoch. Given the supply of z and Y-drops now available from Hubble WFC3/IR surveys, we show it will soon be possible to significantly improve estimates of the Lya fraction using optical and near-IR spectrographs, thereby extending the study conducted in this paper to 7<z<8.

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

We present broad-band imaging with the Subaru Telescope of a 25'x25' field surrounding the radio galaxy TN J1338-1942 at redshift z=4.1. The field contains excesses of Lyman-alpha emitters (LAEs) and Lyman break galaxies (LBGs) identified with a protocluster surrounding the radio galaxy. Our new wide-field images provide information about the boundary of the protocluster and its surroundings. There are 874 candidate LBGs within our field, having redshifts in the range z=3.5-4.5. An examination of the brightest of these (with i'< 25.0) shows that the most prominent concentration coincides with the previously discovered protocluster. The diameter of this galaxy overdensity corresponds to ~2 Mpc at z=4, consistent with the previous estimation using LAEs. Several other concentrations of LBGs are observed in the field, some of which may well be physically connected with the z=4.1 protocluster. The observed structure in the smoothed LBG distribution can be explained as the projection of large-scale structure, within the redshift range z=3.5-4.5, comprising compact overdensities and prominent larger voids. If the 5-8 observed compact overdensities are associated with protoclusters, the observed protocluster volume density is ~5x10^-6 Mpc^-3, similar to the volume density of rich clusters in the local Universe.

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

Using the conditional luminosity function (CLF) -- the luminosity distribution of galaxies in a dark matter halo -- as a way to model galaxy statistics, we study how z=4 Lyman Break Galaxies (LBGs) are distributed in dark matter halos. For this purpose, we measure luminosity-dependent clustering of LBGs in the Subaru/XMM-Newton Deep Field by separating a sample of 16,920 galaxies to three magnitude bins in i'-band between 24.5 and 27.5. Our models fits to data show a possible trend for more luminous galaxies to appear as satellites in more massive halos. The satellite fraction of galaxies at z=4 in these magnitude bins is 0.13 to 0.3, 0.09 to 0.22, and 0.03 to 0.14, respectively, where the 1 sigma ranges account for differences coming from two different estimates of the z=4 LF from the literature. To jointly explain the LF and the large-scale linear bias factor of z=4 LBGs as a function of rest-UV luminosity requires central galaxies to be brighter in UV at z =4 than present-day galaxies in same dark matter mass halos. Moreover, UV luminosity of central galaxies in halos with total mass greater than roughly 10^{12} M_sun must decrease from z=4 to today by an amount more than the luminosity change for galaxies in halos below this mass. This mass-dependent luminosity evolution is preferred at more than 3 sigma confidence level compared to a pure-luminosity evolution scenario where all galaxies decrease in luminosity by the same amount from z=4 to today. The scenario preferred by the data is consistent with the ``down-sizing'' picture of galaxy evolution.

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

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

We report on the angular correlation function of Lyman-break galaxies (LBGs) at z~4 and 5 from deep samples obtained from the Great Observatories Deep Origins Survey (GOODS). Similar to LBGs at z~3, the shape of w(theta) of the GOODS LBGs is well approximated by a power-law with slope beta~0.6 at angular separation theta > 10 arcsec. The clustering strength of z~4, 5 LBGs also depends on the rest-frame UV luminosity, with brighter galaxies more strongly clustered than fainter ones, implying a general correlation between halos' mass and LBGs' star-formation rate. At smaller separations, w(theta) of deep samples significantly exceeds the extrapolation of the large-scale power-law fit, implying enhanced spatial clustering at scales r < 1 Mpc. We also find that bright LBGs statistically have more faint companions on scales theta < 20 arcsec than fainter ones, showing that the enhanced small-scale clustering is very likely due to sub-structure, namely the fact that massive halos can host multiple galaxies. A simple model for the halo occupation distribution and the CDM halo mass function reproduce well the observed w(theta). The scaling relationship of the clustering strength with volume density and with redshift is quantitatively consistent with that of CDM halos. A comparison of the clustering strength of three samples of equal luminosity limit at z ~ 3, 4 and 5 shows that the LBGs at z~5 are hosted in halos about one order of magnitude less massive than those in the lower redshift bins, suggesting that star-formation was more efficient at higher-redshift.

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

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

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

We examine masses of hosting haloes of two photometrically-selected high-z galaxy samples: the old passively-evolving galaxies (OPEGs) and Lyman Break Galaxies (LBGs) both taken from the Subaru/XMM-Newton Deep Survey (SXDS). The large survey area of the SXDS (1sq deg) allows us to measure the angular two-point correlation functions to a wide separation of >10 arcmin with a good statistical quality. We utilize the halo model prescription for estimating characteristic masses of hosting haloes from the measured large-scale clustering amplitudes. It is found that the hosting halo mass positively correlates with the luminosity of galaxies. Then, adopting the extended Press-Schechter model (EPS), we compute the predictions for the mass evolution of the hosting haloes in the framework of the cold dark matter (CDM) cosmology in order to make an evolutionary link between the two galaxy samples at different redshifts and to identify their present-day descendants by letting their haloes evolve forward in time. It is found that, in the view of the mass evolution of hosting haloes in the CDM model, bright LBGs are consistent with being the progenitor of the OPEGs, whereas it is less likely that the LBG population, as a whole, have evolved into the OPEG population. It is also found that the present-day descendants of both the bright LBGs and OPEGs are likely to be located in massive systems such as groups of galaxies or clusters of galaxies. Finally, we estimate the hosting halo mass of local early-type galaxy samples from the 2dF and SDSS based on the halo model and it turns out that their expected characteristic mass of hosting haloes is in good agreement with the EPS predictions for the descendant's mass of both the bright LBGs and OPEGs.

We use data from the first epoch of observations with the IRAC/Spitzer for the GOODS to detect and study a collection of LBGs at z ~ 6 to 5 in the HUDF, six of which have spectroscopic confirmation. At these redshifts, IRAC samples rest-frame optical light in the range 0.5 to 0.8 um, where the effects of dust extinction are smaller and the sensitivity to light from evolved stars is greater than at shorter, rest-frame UV observable from the ground or with the HST. As such, it provides useful constraints on the ages and masses of these galaxies' stellar populations. We find that the SEDs for many of these galaxies are best fitted by models of stellar populations with masses of a few x 1e10 M_sun, and with ages of a few hundred Myrs, values quite similar to those derived for typical LBGs at z ~ 3. When the universe was only 1 Gyr old, some galaxies had already formed a mass of stars approaching that of the present-day Milky Way, and that they started forming those stars at z > 7, and in some cases much earlier. We find that the lower limits to the space density for galaxies in this mass range are consistent with predictions from recent hydrodynamic simulations of structure formation in a LCDM universe. All objects in our samples are consistent with having solar metallicity, suggesting that they might have already been significantly polluted by metals. The values for dust reddening derived from the model fitting are low or zero, and we find that some of the galaxies have rest-frame UV colors that are even bluer than those predicted by the stellar population models to which we compare them. These colors might be attributed to the presence of very massive stars (> 100 M_sun), or by weaker intergalactic HI absorption than what is commonly assumed.

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

We investigate expected sub-mm emission and dust properties of the Lyman Break Galaxies (LBGs) in the Hubble Deep Field (HDF). The SCUBA Deep Survey (Hughes et al., 1998) provides an upper limit of the 850 micron flux densities of the LBGs. With this constraint, we argue that a typical ultra-violet to far-infrared spectral shape of the high-redshift LBGs is rather close to a template spectrum of low-reddening starburst galaxies in the local universe but different from that of heavily dust-enshrouded ultra-luminous far-infrared galaxies like Arp220. We also evaluate the lower-limit temperature of dust in LBGs assuming single- and two-component modified black-body spectrum. To estimate the total amount of energy re-emitted in FIR wavelength, we take two different approaches; model fitting of the UV spectra of LBGs and an empirical method using the relationship between UV spectral shape and UV/FIR flux ratio observed for local starburst galaxies. Both methods give lower-limit temperature of the LBGs as around 40 K, which is higher than the typical dust temperature of local optical- and infrared-selected galaxies. This result is also supported by the comparison of the expected sub-mm flux of the LBGs with the cosmic FIR background radiation. The high dust temperature may indicate the effective massive-star formation or different dust properties in the high-redshift LBGs.