Search

We present Ly$\alpha$ and ultraviolet-continuum (UV) luminosity functions (LFs) of galaxies and active galactic nuclei (AGN) at $z=2.0-3.5$ determined by the un-targetted optical spectroscopic survey of the Hobby-Eberly Telescope Dark Energy Experiment (HETDEX). We combine deep Subaru imaging with HETDEX spectra resulting in $11.4$ deg$^2$ of fiber-spectra sky coverage, obtaining $18320$ galaxies spectroscopically identified with Ly$\alpha$ emission, $2126$ of which host type 1 AGN showing broad (FWHM$~>1000$ km s$^{-1}$) Ly$\alpha$ emission lines. We derive the Ly$\alpha$ (UV) LF over 2 orders of magnitude covering bright galaxies and AGN in $\log L_\mathrm{Ly\alpha}/\mathrm{[erg~s^{-1}]}=43.3-45.5$ ($-27-21$) type 1 AGN increases from $z\sim2$ to $z\sim0$ as opposed to the evolution of bright ($M_\mathrm{UV}, Comment: 27 pages, 19 figures, accepted for publication in ApJ

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 present the first measurement of the metallicity gradient in extremely metal-poor galaxies (EMPGs). With Subaru/Faint Object Camera And Spectrograph (FOCAS) Integral Field Unit (IFU), we have observed a nearby, low-mass EMPG, HSC J1631+4426, whose oxygen abundance and stellar mass are known to be 12+log(O/H) $=6.9$ and $\log_{10}(M_*/{\rm M}_\odot)=5.8$, respectively. The measured metallicity gradient is $-0.36 \pm 0.04$ dex kpc$^{-1}$ corresponding to $-0.049 \pm 0.006$ dex R$_\mathrm{e}^{-1}$ for the continuum effective radius of $R_\mathrm{e} = 0.14$ kpc. Our observation has successfully demonstrated that three-dimensional spectroscopy with 8m-class telescopes is powerful enough to reveal the metallicity distribution in local EMPGs, providing precious information of the baryon cycle in local analogs of primordial galaxies in the early Universe., Comment: PASJ accepted, 7 pages, 4 figures

We present cosmological ($30-400$ Mpc) distributions of neutral hydrogen (HI) in the inter-galactic medium (IGM) traced by Ly$\alpha$ Emitters (LAEs) and QSOs at $z=2.1-2.5$, selected with the data of the on-going Hobby-Eberly Telescope Dark Energy Experiment (HETDEX) and the eBOSS survey. We investigate spatial correlations of LAEs and QSOs with HI tomography maps reconstructed from HI Ly$\alpha$ forest absorption in the spectra of background galaxies and QSOs obtained by the CLAMATO survey and this study, respectively. In the cosmological volume far from QSOs, we find that LAEs reside in regions of strong HI absorption, i.e. HI rich, which is consistent with results of previous galaxy-background QSO pair studies. Moreover, there is an anisotropy in the HI-distribution plot of transverse and line-of-sight distances; on average the HI absorption peak is blueshifted by $\sim 200$ km s$^{-1}$ from the LAE Ly$\alpha$ redshift, reproducing the known average velocity offset between the Ly$\alpha$ emission redshift and the galaxy systemic redshift. We have identified a $\sim$ 40-Mpc scale volume of HI underdensity that is a candidate for a giant HII bubble, where six QSOs and an LAE overdensity exist at $\left < z \right > =2.16$. The coincidence of the QSO and LAE overdensities with the HI underdensity indicates that the ionizing photon radiation of the QSOs has created a highly ionized volume of multiple proximity zones in a matter overdensity. Our results suggest an evolutionary picture where HI gas in an overdensity of galaxies becomes highly photoionized when QSOs emerge in the galaxies., Comment: 26 pages, 19 figures. Accepted for Publication in ApJ

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

We report a discovery of $z = 4.72$ radio galaxy, HSC J083913.17+011308.1, by using the Lyman break technique with the Hyper Suprime-Cam Subaru Strategic Survey (HSC-SSP) catalog for VLA FIRST radio sources. The number of known high-$z$ radio galaxies (HzRGs) at $z > 3$ is quite small to constrain the evolution of HzRGs so far. The deep and wide-area optical survey by HSC-SSP enables us to apply the Lyman break technique to a large search for HzRGs. For an HzRG candidate among pre-selected $r$-band dropouts with a radio detection, a follow-up optical spectroscopy with GMOS/Gemini has been performed. The obtained spectrum presents a clear Ly$\alpha$ emission line redshifted to $z=4.72$. The SED fitting analysis with the rest-frame UV and optical photometries suggests the massive nature of this HzRG with $\log{M_*/M_{\odot}} = 11.4$. The small equivalent width of Ly$\alpha$ and the moderately red UV colors indicate its dusty host galaxy, implying a chemically evolved and dusty system. The radio spectral index does not meet a criterion for an ultra-steep spectrum: $\alpha^{325}_{1400}$ of $-1.1$ and $\alpha^{150}_{1400}$ of $-0.9$, demonstrating that the HSC-SSP survey compensates for a sub-population of HzRGs which are missed in surveys focusing on an ultra-steep spectral index., Comment: 10 pages, 5 figures, accepted for publication in AJ

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

Tang, S., et al., We report on a spectroscopic program to search for dual quasars using Subaru Hyper Suprime-Cam (HSC) images of SDSS quasars, which represent an important stage during galaxy mergers. Using Subaru/FOCAS and Gemini-N/GMOS, we identify three new physically associated quasar pairs having projected separations less than 20 kpc, out of 26 observed candidates. These include the discovery of the highest-redshift (z = 3.1) quasar pair with a separation, Funding was contributed by the FIRST program from the 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

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

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

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

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

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

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

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