• 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 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 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 agree 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 streams (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 a strong-lensing analysis of four massive galaxy clusters imaged with the Hubble Space Telescope in the framework of the Reionization Lensing Cluster Survey (RELICS). We use a Light-Traces-Mass modeling technique to uncover sets of multiply imaged galaxies, and constrain the mass distribution and strong-lensing properties of the clusters. The mass models we present here are the first published for Abell S295 and MACS J0159.8-0849. For Abell 697 (the tenth-highest Sunyaev-Zel'dovich mass cluster in the Planck catalog) and MACS J0025.4-1222 (the "baby bullet" cluster), thanks to RELICS data we are able to improve upon previous models. Our analysis for MACS J0025.4-1222 and Abell S295 shows a bimodal mass distribution following the cluster galaxy concentrations, in support of the merger scenarios proposed in previous studies for these clusters. In addition, 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, suggesting these are fairly relaxed clusters. Despite being smaller and less prominent lenses on average, three of the four clusters we analyze here seem to have lensing strengths similar to the typical Hubble Frontier Fields cluster in terms of the cumulative area above a certain magnification value (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, including mass-density, deflection, shear and magnifications maps.
  • 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 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 and the ALMA archive. 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.5 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), and supports a positive correlation of the FIR size-luminosity relation. We also examine the rest-frame optical Sersic index $n$_opt and effective radius $R$_e,opt with the deep Hubble Space Telescope (HST) images. Interestingly, we obtain $n$_opt=0.9+/-0.3 (~ $n$_FIR) and $R$_e,opt=3.2+/-0.6 kpc ($>R$_e,FIR), suggesting that the FIR-emitting disk is embedded within a larger stellar disk. The rest-frame UV and FIR data of HST and ALMA provide us a radial surface density profile of the total star-formation rate (SFR), where the FIR SFR dominates over the UV SFR at the center. Under the simple assumption of a constant SFR, a compact stellar distribution found in z~1-2 compact quiescent galaxies (cQGs) is well reproduced, while a spheroidal stellar morphology of cQGs ($n$_opt=4) cannot, suggestive of other important mechanisms such as dynamical dissipation.
  • We 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 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 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 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.
  • Hiroaki Aihara, Robert Armstrong, Steven Bickerton, James Bosch, Jean Coupon, Hisanori Furusawa, Yusuke Hayashi, Hiroyuki Ikeda, Yukiko Kamata, Hiroshi Karoji, Satoshi Kawanomoto, Michitaro Koike, Yutaka Komiyama, Robert H. Lupton, Sogo Mineo, Hironao Miyatake, Satoshi Miyazaki, Tomoki Morokuma, Yoshiyuki Obuchi, Yukie Oishi, Yuki Okura, Paul A. Price, Tadafumi Takata, Manobu M. Tanaka, Masayuki Tanaka, Yoko Tanaka, Tomohisa Uchida, Fumihiro Uraguchi, Yousuke Utsumi, Shiang-Yu Wang, Yoshihiko Yamada, Hitomi Yamanoi, Naoki Yasuda, Nobuo Arimoto, Masashi Chiba, Francois Finet, Hiroki Fujimori, Seiji Fujimoto, Junko Furusawa, Tomotsugu Goto, Andy Goulding, James E. Gunn, Yuichi Harikane, Takashi Hattori, Masao Hayashi, Krzysztof G. Helminiak, Ryo Higuchi, Chiaki Hikage, Paul T.P. Ho, Bau-Ching Hsieh, Kuiyun Huang, Song Huang, Masatoshi Imanishi, Ikuru Iwata, Anton T. Jaelani, Hung-Yu Jian, Nobunari Kashikawa, Nobuhiko Katayama, Takashi Kojima, Akira Konno, Shintaro Koshida, Haruka Kusakabe, Alexie Leauthaud, C.-H. Lee, Lihwai Lin, Yen-Ting Lin, Rachel Mandelbaum, Yoshiki Matsuoka, Elinor Medezinski, Shoken Miyama, Rieko Momose, Anupreeta More, Surhud More, Shiro Mukae, Ryoma Murata, Hitoshi Murayama, Tohru Nagao, Fumiaki Nakata, Hiroko Niikura, Atsushi J. Nishizawa, Masamune Oguri, Nobuhiro Okabe, Yoshiaki Ono, Masato Onodera, Masafusa Onoue, Masami Ouchi, Tae-Soo Pyo, Takatoshi Shibuya, Kazuhiro Shimasaku, Melanie Simet, Joshua Speagle, David N. Spergel, Michael A. Strauss, Yuma Sugahara, Naoshi Sugiyama, Yasushi Suto, Nao Suzuki, Philip J. Tait, Masahiro Takada, Tsuyoshi Terai, Yoshiki Toba, Edwin L. Turner, Hisakazu Uchiyama, Keiichi Umetsu, Yuji Urata, Tomonori Usuda, Sherry Yeh, Suraphong Yuma
    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.