Search

We perform a super-resolution analysis of the Subaru Hyper Suprime-Cam (HSC) images to estimate the major merger fractions of z~4-7 dropout galaxies at the bright end of galaxy UV luminosity functions (LFs). Our super-resolution technique improves the spatial resolution of the ground-based HSC images, from ~1" to ~90%. We apply the super-resolution technique to 6412, 16, 94, and 13 very bright dropout galaxies at z~4, 5, 6, and 7, respectively, in a UV luminosity range of L_UV~3-15 L_UV* corresponding to -24~10 L_UV*, even after the AGN contribution is subtracted. The results support scenarios in which other additional mechanisms such as insignificant mass quenching and low dust obscuration contribute to the DPL shape of galaxy UV LFs., Comment: 20 pages, 9 figures. Accepted for publication in PASJ

Hydrogen Lyman-$\alpha$ (Ly$\alpha$) emission has been one of the major observational probes for the high redshift universe, since the first discoveries of high-$z$ Ly$\alpha$ emitting galaxies in the late 1990s. Due to the strong Ly$\alpha$ emission originated by resonant scattering and recombination of the most-abundant element, Ly$\alpha$ observations witness not only HII regions of star formation and AGN but also diffuse HI gas in the circum-galactic medium (CGM) and the inter-galactic medium (IGM). Here we review Ly$\alpha$ sources, and present theoretical interpretations reached to date. We conclude that: 1) A typical Ly$\alpha$ emitter (LAE) at $z\gtrsim 2$ with a $L^*$ Ly$\alpha$ luminosity is a high-$z$ counterpart of a local dwarf galaxy, a compact metal-poor star-forming galaxy (SFG) with an approximate stellar (halo) mass and star-formation rate of $10^{8-9} M_\odot$ ($10^{10-11} M_\odot$) and $1-10 M_\odot$ yr$^{-1}$, respectively; 2) High-$z$ SFGs ubiquitously have a diffuse Ly$\alpha$ emitting halo in the CGM extending to the halo virial radius and beyond; 3) Remaining neutral hydrogen at the epoch of reionization makes a strong dimming of Ly$\alpha$ emission for galaxies at $z>6$ that suggest the late reionization history. The next generation large telescope projects will combine Ly$\alpha$ emission data with HI Ly$\alpha$ absorptions and 21cm radio data that map out the majority of hydrogen (HI+HII) gas, uncovering the exchanges of i) matter by outflow/inflow and ii) radiation, relevant to cosmic reionization, between galaxies and the CGM/IGM., Comment: 45 pages, 18 figures. Invited review published in Annual Reviews of Astronomy and Astrophysics. See https://www.annualreviews.org/doi/abs/10.1146/annurev-astro-032620-021859 for updated version

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

The primordial He abundance $Y_\mathrm{P}$ is a powerful probe of cosmology. Currently, $Y_\mathrm{P}$ is best determined by observations of metal-poor galaxies, while there are only a few known local extremely metal-poor ($, Comment: 22 pages, 9 figures, accepted for publication in ApJ

We present the galaxy size-mass ($R_{e}-M_{\ast}$) distributions using a stellar-mass complete sample of $\sim1.5$ million galaxies, covering $\sim100$ deg$^2$, with $\log(M_{\ast}/M_{\odot})>10.2~(9.2)$ over the redshift range $0.2, Comment: Accepted for publication in ApJ. 45 pages, 17 figures

We report fourteen and twenty-eight protocluster candidates at z = 5.7 and 6.6 over 14 and 19 deg2 areas, respectively, selected from 2,230 Lyα emitters (LAEs) photometrically identified with Subaru/Hyper Suprime-Cam (HSC) deep images. Six out of the 42 protocluster candidates include at least 1 spectroscopically confirmed LAEs at redshifts up to z = 6.574. By the comparisons with the cosmological Lyα radiative transfer (RT) model reproducing LAEs with the reionization effects, we find that more than a half of these protocluster candidates might be progenitors of the present-day clusters with a mass of ≳ 1014M⊙. We also investigate the correlation between LAE overdensity and Lya rest-frame equivalent width (EW), because the cosmological Lyα RT model suggests that a slope of EW-overdensity 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 Lyα emission from the partly neutral intergalactic medium. The available HSC data suggest that the slope of the EW-overdensity correlation does not evolve from the post-reionization epoch z = 5.7 to the EoR z = 6.6 beyond the moderately large statistical errors.

We present ALMA detection of the [O iii] 88 μm line and 850 μm dust continuum emission in a Y-dropout Lyman break galaxy, MACS0416_Y1. The [O iii] detection confirms the object with a spectroscopic redshift to be z = 8.3118±0.0003. The 850 μm continuum intensity (0.14 mJy) implies a large dust mass on the order of 4×106M⊙. The ultraviolet-to-far infrared spectral energy distribution modeling, where the [O iii] emissivity model is incorporated, suggests the presence of a young (τage ≍ 4 Myr), star-forming (SFR ≍ 60M⊙yr−1), and moderately metal-polluted (Z ≍ 0.2Z⊙) stellar component with a stellar mass of 3 × 108M⊙. An analytic dust mass evolution model with a single episode of star formation does not reproduce the metallicity and dust mass in ≍ 4 Myr, suggesting an underlying evolved stellar component as the origin of the dust mass.

We conduct intensity mapping to probe for extended diffuse Lyα emission around Lyα emitters (LAEs) at z ∼2−7, exploiting very deep (∼26 mag at 5σ) and large-area (∼4.5 deg2) Subaru/Hyper Suprime-Cam narrowband (NB) images and large LAE catalogs consisting of a total of 1540 LAEs at z = 2.2, 3.3, 5.7, and 6.6 obtained by the HSC-SSP and CHORUS projects. We calculate the spatial correlations of these LAEs with ∼1–2 billion pixel flux values of the NB images, deriving the average Lyα surface brightness (SBLyα ) radial profiles around the LAEs. By carefully estimating systematics such as fluctuations of sky background and point-spread functions, we detect Lyα emission at 100–1000 comoving kpc around z = 3.3 and 5.7 LAEs at the 3.2σ and 3.7σ levels, respectively, and tentatively (=2.0σ) at z = 6.6. The emission is as diffuse as ∼10−20–10−19 erg s−1 cm−2 arcsec−2 and extended beyond the virial radius of a dark matter halo with a mass of 1011 M ⊙. While the observed SBLyα profiles have similar amplitudes at z = 2.2–6.6 within the uncertainties, the intrinsic SBLyα profiles (corrected for the cosmological dimming effect) increase toward high redshifts. This trend may be explained by increasing hydrogen gas density due to the evolution of the cosmic volume. Comparisons with theoretical models suggest that extended Lyα emission around an LAE is powered by resonantly scattered Lyα photons in the CGM and IGM that originate from the inner part of the LAE and/or neighboring galaxies around the LAE.

Onoue, Masafusa, et al., We present deep Keck/MOSFIRE near-infrared spectroscopy of a strong Lyα emitting source at z = 6.1292, HSC J142331.71-001809.1, which was discovered by the SHELLQS program from imaging data of the Subaru Hyper Suprime-Cam (HSC) survey. This source is one of five objects that show narrow (FWHM < 500 km s-1) and prominent (LLyα > 1044 erg s-1) Lyα emission lines at absolute 1450 Å continuum magnitudes of M1450∼-22 mag. Its rest-frame Lyα equivalent width (EW) is 370 ± 30 Å. In the 2 hr Keck/MOSFIRE spectrum in Y band, the high-ionization C IV λλ1548,1550 doublet emission line was clearly detected with = - FWHM 120+20 20 km s-1 and a total rest-frame EW of - 37+5 6 Å. We also report the detection of weak continuum emission, and the tentative detection of O III]λλ1661,1666 in the 4 hr J-band spectrum. Judging from the UV magnitude, line widths, luminosities, and EWs of Lyα and C IV, we suggest that this source is a reionization-era analog of classical type-II AGNs, although there is a possibility that it represents a new population of AGN/galaxy composite objects in the early universe. We compare the properties of J1423-0018 to intermediate-redshift type-II AGNs and C IV emitters seen in z = 6-7 galaxy samples. Further observations of other metal emission lines in the rest-frame UV or optical, along with X-ray follow-up observations of the z = 6-7 narrow-line quasars, are needed for more robust diagnostics and to determine their nature., Funding was contributed by the FIRST program from Japanese Cabinet Office, the Ministry of Education, Culture, Sports, Science and Technology (MEXT), the Japan Society for the Promotion of Science (JSPS), Japan Science and Technology Agency (JST), the Toray Science Foundation, NAOJ, Kavli IPMU, KEK, ASIAA, and Princeton University. This work was supported by JSPS KAKENHI grant Nos. 17H04830 (YM), 20H01949 (TN), 18J01050 (YT). M.O. was supported by the ERC Advanced Grant 740246 Cosmic gas. Y. M. was supported by the Mitsubishi Foundation grant No. 30140. T.H. was supported by Leading Initiative for Excellent Young Researchers, MEXT, Japan (HJH02007).

Adaptive optics (AO) is a powerful tool for correcting wavefront errors induced by complex structures of biological samples which significantly causes image degradation. A scene-based sensing technique is being popular in microscopic AO systems with Shack-Hartmann (SH) wavefront sensors. A problem in application of the technique is that the shapes of images observed on SAs vary dependently to their positions on the aperture, especially when using microscopic objectives with higher NAs. To mitigate this problem, a differential sensing technique is used that enables measuring image shifts with high correlations over the aperture. Experiments using an artificial testing target including fluorescent beads, which simulates the leaf of moss, were conducted to investigate imaging performances of the present AO system. Unbiased maximum ratios were measured from blurred and AO-corrected images, and then the Strehl ratios were derived from them. Resultant Strehl ratios were around 0.58.

We develop a machine learning (ML) software to estimate morphological parameters (e.g., the half-light radius re) of high redshift galaxies in the Subaru/Hyper Suprime-Cam data. To make the ML software capture simultaneously galaxy morphological features and point spread function (PSF) broadening effects, we implement a two-stream convolutional neural network (CNN) for inputs of galaxy and PSF images. Thanks to large training samples of galaxy and PSF images, the two-stream CNN estimates re more accurately than a single-stream CNN with only galaxy images. Our ML software would be a useful tool to investigate galaxy morphological properties with PSF-unstable images obtained in future large-area ground-based surveys.

We present the selection, spectroscopic identification, and physical properties of extreme emission line galaxies (EELGs) at $33$ with intense [OIII] emission, and H$\alpha$ emitters at $z\simeq 2$, respectively. These spectroscopically identified EELGs at $z\simeq 3.3$ show, on average, higher specific star formation rates (sSFR) than the star-forming main sequence, low dust attenuation of $E(B-V) \lesssim 0.1$ mag, and high [OIII]/[OII] ratios of $\gtrsim 3$. We also find that our EELGs at $z\simeq 3.3$ have higher hydrogen ionizing photon production efficiencies ($\xi_\mathrm{ion}$) than the canonical value ($\simeq 10^{25.2}$ Hz/erg), indicating that they are efficient in ionizing their surrounding interstellar medium. These physical properties suggest that they are low metallicity galaxies with higher ionization parameters and harder UV spectra than normal SFGs, which is similar to galaxies with Lyman continuum (LyC) leakage. Among our EELGs, those with the largest [OIII]/[OII] and EW([OIII]) values would be the most promising candidates to search for LyC leakage., Comment: 29 pages, 16 figures, 4 tables, accepted for publication in ApJ

Significance: A scene-based adaptive-optics (AO) system is developed and a method for investigating its imaging performance is proposed. The system enables derivation of Strehl ratios from observed images via collaboration with computer simulations. The resultant Strehl ratios are comparable with those of other current AO systems. Aim: For versatile and noninvasive AO microscopy, a scene-based wavefront-sensing technique working on a Shack–Hartmann wavefront sensor is developed in a modal control system. The purpose of the research is to clarify the imaging performance of the AO system via the derivation of Strehl ratios from observed images toward applications in microscopy of living cells and tissues. Approach: Two imaging metrics that can be directly measured from observed images (i.e., an energy concentration ratio and unbiased maximum ratio) are defined and related to the Strehl ratio via computer simulations. Experiments are conducted using artificial targets to measure the imaging metrics, which are then converted to Strehl ratios. Results: The resultant Strehl ratios are >0.7 and 0.5 in the cases of defocus and higher aberrations, respectively. The half-widths at half-maximum of the AO-corrected bead images are favorably comparable to those of on-focus images under simple defocus aberration, and the AO system works both under bright-field illumination and on fluorescent bead images. Conclusions: The proposed scene-based AO system is expected to work with a Strehl ratio of more than 0.5 when applied to high-resolution live imaging of cells and tissues under bright-field and fluorescence microscopies.

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., Comment: 20 pages, 18 figures, Accepted for publication in ApJ

The spatial inhomogeneity is one of the important features for understanding the reionization process; however, it has not yet been fully quantified. To map this inhomogeneous distribution, we simultaneously detect Ly$\alpha$ emitters (LAEs) and Lyman break galaxies (LBGs) at $z \sim 6.6$ from the Subaru/Hyper Suprime-Cam (HSC) large-area ($\sim1.5\,\mathrm{ deg}^2 = 34000\,\mathrm{cMpc}^2$) deep survey. We estimate the neutral fraction, $x_\mathrm{HI}$, from the observed number density ratio of LAEs to LBGs, $n(\mathrm{LAE})/n(\mathrm{LBG})$ based on numerical radiative transfer simulation, in which model galaxies are selected to satisfy the observed selection function. While the average $x_\mathrm{HI}$ within the field of view is found to be $x_\mathrm{HI} < 0.4$, which is consistent with previous studies, the variation of $n(\mathrm{LAE})/n(\mathrm{LBG})$ within the field of view for each $140\,\mathrm{pMpc}^2$ is found to be as large as a factor of three. This may suggest a spatially inhomogeneous topology of reionization, but it also leaves open the possibility that the variation is based on the inherent large-scale structure of the galaxy distribution. Based on the simulations, it may be difficult to distinguish between the two from the current survey. We also find that LAEs in the high LAE density region are more populate high $\mathrm{EW}_0$, supporting that the observed $n(\mathrm{LAE})/n(\mathrm{LBG})$ is more or less driven by the neutral fraction, though the statistical significance is not high., Comment: 20 pages, 22 figures, accepted for publication in ApJ

We present new measurements of rest-UV luminosity functions and angular correlation functions from 4,100,221 galaxies at z~2-7 identified in the Subaru/Hyper Suprime-Cam survey and CFHT Large-Area U-band Survey. The obtained luminosity functions at z~4-7 cover a very wide UV luminosity range of ~0.002-2000L*uv combined with previous studies, revealing that the dropout luminosity function is a superposition of the AGN luminosity function dominant at Muv-22 mag, consistent with galaxy fractions based on 1037 spectroscopically-identified sources. Galaxy luminosity functions estimated from the spectroscopic galaxy fractions show the bright end excess beyond the Schechter function at >2sigma levels, which is possibly made by inefficient mass quenching, low dust obscuration, and/or hidden AGN activity. By analyzing the correlation functions at z~2-6 with halo occupation distribution models, we find a weak redshift evolution (within 0.3 dex) of the ratio of the star formation rate (SFR) to the dark matter accretion rate, SFR/(dMh/dt), indicating the almost constant star formation efficiency at z~2-6, as suggested by our earlier work at z~4-7. Meanwhile, the ratio gradually increases with decreasing redshift at z10 with $\propto10^{-0.5(1+z)}$, which will be directly tested with JWST., Comment: 53 pages, 26 figures, 10 tables, accepted for publication in ApJS. Dropout galaxy catalogs are available on our project webpage (http://cos.icrr.u-tokyo.ac.jp/rush.html)

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 {\sc Galfit} software, carefully removing the SB contributions of tails. We find that the EMPGs with a median stellar mass of $\log(M_{*}/{\rm M}_{\odot})=6.0$ have a median S{\'e}rsic index of $n=1.1$ and a median effective radius of $r_{\rm e}=200$ pc, suggesting that typical EMPGs have very compact disk. We compare the EMPGs with $z\sim6$ 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. A spectrum of one pair of EMPG and tail, so far available, indicates that the tail is dynamically related to the EMPG with a median velocity difference of $\Delta V=101\pm32$ km s$^{-1}$. This moderately-large $\Delta V$ cannot be explained by the dynamics of the tail, but likely by the infall on the tail. For the first time, we may identify the metal-poor star-forming system just now infalling into the tail., Comment: ApJ in Press

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

To determine the dominant sources for cosmic reionization, the evolution history of the global ionizing fraction, and the topology of the ionized regions, we have conducted a deep imaging survey using four narrow-band (NB) and one intermediate-band (IB) filters on the Subaru/Hyper Suprime-Cam (HSC), called Cosmic HydrOgen Reionization Unveiled with Subaru (CHORUS). The central wavelengths and full-widths-at-half-maximum of the CHORUS filters are, respectively, 386.2 nm and 5.5 nm for NB387, 526.0 nm and 7.9 nm for NB527, 717.1 nm and 11.1 nm for NB718, 946.2 nm and 33.0 nm for IB945, and 971.2 nm and 11.2 nm for NB973. This combination, including NB921 (921.5 nm and 13.5 nm) from the Subaru Strategic Program with HSC (HSC SSP), are carefully designed, as if they were playing a chorus, to observe multiple spectral features simultaneously, such as Lyman continuum, Ly$\alpha$, C~{\sc iv}, and He~{\sc ii} for $z=2$--$7$. The observing field is the same as that of the deepest footprint of the HSC SSP in the COSMOS field and its effective area is about 1.6 deg$^2$. Here, we present an overview of the CHORUS project, which includes descriptions of the filter design philosophy, observations and data reduction, multiband photometric catalogs, assessments of the imaging quality, measurements of the number counts, and example use cases of the data. All the imaging data, photometric catalogs, masked pixel images, data of limiting magnitudes and point spread functions, results of completeness simulations, and source number counts are publicly available through the HSC SSP database., Comment: 17 pages, 6 figures, PASJ open access (https://academic.oup.com/pasj/advance-article-abstract/doi/10.1093/pasj/psaa100/5959796), all the data products are available at https://hsc-release.mtk.nao.ac.jp/doc/index.php/chorus

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

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}$., Comment: 17 pages, 11 figures, Accepted for Publication in ApJ

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$., Comment: 16 pages, 12 figures, published in ApJ

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., Comment: 24 pages, 13 figures, 9 tables. Accepted for publication in PASJ. In addition to the dust continuum detection in ALMA Band 6 reported in the original manuscript, we newly report the dust continuum detection in ALMA Band 8 based on our new ALMA Cycle 5 observations

We present results of the cross-correlation Ly$\alpha$ intensity mapping with Subaru/Hyper Suprime-Cam (HSC) ultra-deep narrowband images and Ly$\alpha$ 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, we further conduct extensive analysis evaluating systematics of large-scale point-spread-function wings, sky subtractions, and unknown errors based on physically uncorrelated signals and sources found in real HSC images and object catalogs, respectively. Removing the systematics, we carefully calculate cross-correlations between Ly$\alpha$ intensity of the narrowband images and the LAEs. We tentatively identify very diffuse Ly$\alpha$ emission with the $\simeq 3\sigma$ ($\simeq 2\sigma$) significance at $\gtrsim$ 100 comoving kpc (ckpc) far from the LAEs at $z=5.7$ ($6.6$), around and probably even beyond a virial radius of star-forming galaxies with $M_\mathrm{h}\sim10^{11}M_\odot$. The diffuse Ly$\alpha$ 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 ($, Comment: 11 pages, 6 figures, accepted for publication in ApJ

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., Comment: Submitted to ApJ

We present element abundance ratios and ionizing radiation of local young low-mass (∼106 M ⊙) extremely metal-poor galaxies (EMPGs) with a 2% solar oxygen abundance (O/H)⊙ 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 [Fe iii] λ4658 and He ii λ4686 are detected in very deep optical spectra of the EMPGs taken with 8 m class telescopes, including Keck and Subaru, enabling us to derive element abundance ratios with photoionization models. We find that neon-to-oxygen 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)⊙, 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)⊙, which are unlikely to be 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, the EMPGs with the 2%-solar oxygen abundance have very high He ii λ4686/Hβ ratios of ∼1/40, which are not reproduced by existing models of high-mass X-ray binaries with progenitor stellar masses M ⊙. 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 He ii λ4686/Hβ ratios are explained by the past existence of supermassive (>300 M ⊙) stars, which may evolve into intermediate-mass black holes (≳100 M ⊙).

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_sun6 reported recently. We spectroscopically confirm that 3 out of 10 candidates are truly EMPGs with Z/Z_sun, Comment: 30 pages, 26 figures, and 8 tables; (Revision: metallicities changed due to the use of a more pricise dust correction method, but conclusion does not change, Minor: a table, some figures, and sentences are added for a clear explanation.)

We present an IGM H I tomographic map in a survey volume of 16 x 19 x 131 h(-3) comoving Mpc(3) (cMpc(3)) centered at MAMMOTH-1 nebula and three neighboring quasars at z = 2.3. The 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 no clear physical origin. Here we investigate the H I-gas distribution around MAMMOTH1-QSO with the ELAN and three neighboring type-I quasars, making the IGM H I tomographic map with a spatial resolution of 2.6 h(-1) cMpc. Our H I tomographic map is reconstructed with H I Lya 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 H I overdensity for MAMMOTH1-QSO, and find that MAMMOTH1-QSO resides in a volume with fairly weak H I absorption. This suggests that MAMMOTH1-QSO may have a proximity zone where quasar illuminates and photoionizes the surrounding H I gas and suppresses H I absorption, and that the ELAN is probably a photoionized cloud embedded in the cosmic web. The H I radial profile of MAMMOTH1-QSO is very similar to those of three neighboring type-I quasars at z = 2.3, which is compatible with the AGN unification model. We compare the distributions of the H I absorption and star-forming galaxies in our survey volume, and identify a spatial offset between density peaks of star-forming galaxies and H I gas. This segregation may suggest anisotropic UV background radiation created by star-forming galaxy density fluctuations. (© 2020 The American Astronomical Society). ISSN:0004-637X ISSN:2041-8213

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

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., 24 pages, 21 figures, accepted in ApJ

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., Comment: 13 pages, 15 figures, accepted in ApJ

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, Comment: PASJ in press, 30 pages, 20 figures

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 are carrying out the study of the evolution of radial surface brightness profiles of galaxies from z = 0 to 2 by stacking analysis using data corrected by the Hyper Suprime-Cam (HSC) Subaru Strategic Program (SSP). This will allow us to constrain the large scale average profiles of various galaxy populations at high redshift. From the stacking analysis of galaxies selected based on their photometric redshifts, we successfully detected the outer components of galaxies at z > 1 extending to at least ~80 kpc, which imply an early formation for the galaxy outskirts.

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

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., Comment: 30 pages, 13 figures, accepted for publication in PASJ special issue

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., Comment: 14 pages, 13 figures, accepted for publication in ApJ

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., 19 pages, 11 figures. Accepted for publication in PASJ special issue

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., Comment: 26 pages, 12 figures. Accepted for publication in PASJ special issue

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., Comment: 16 pages, 10 figures. Accepted for publication in PASJ special issue