Researcher profile

Yuichi Harikane

Yuichi Harikane contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
32works
0followers
4topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

32 published item(s)

preprint2026arXiv

A New Constraint on the Optical Depth from the Reionization History Independent of CMB Large-Scale E-Mode Polarization

Recent studies report a mild discrepancy between baryon acoustic oscillation (BAO) and cosmic microwave background (CMB) measurements within the $Λ$CDM framework. This discrepancy could be explained if the optical depth $τ$ inferred from the CMB large-scale E-mode polarization is underestimated, which may be biased by foreground-subtraction or instrumental systematics. In this work, we present a determination of $τ$ independent of the large-scale E-mode polarization, using the latest measurements of the redshift evolution of the neutral hydrogen fraction $x_\mathrm{HI}(z)$, which is constrained by Lyman-$α$ forest and damping-wing absorption measurements at $z\sim5$-$14$, based on ground-based optical and JWST observations. Combining $x_\mathrm{HI}(z)$ with the Planck CMB power spectra excluding the large-scale E-mode polarization, we obtain $τ=0.0552^{+0.0019}_{-0.0026}$, a stringent constraint consistent with the previous CMB results including the large-scale E-mode. We also evaluate a potential systematic error in our method associated with absorption modeling, obtaining $τ=0.0552^{+0.0075}_{-0.0049}$. Using this constraint on $τ$, we resolve the degeneracy in the $τ$-$Ω_m$ plane and find a $2.4σ$ tension with the DESI DR2 BAO results, thereby confirming the claimed mild discrepancy suggestive of physics beyond $Λ$CDM. Finally, we derive an upper limit on the sum of neutrino masses, $Σm_ν<0.0550\,(0.0717)$ eV at the 95% (99%) confidence level. This limit favors the normal mass ordering and, when combined with the lower limits from neutrino oscillation experiments, yields a further constraint, $Σm_ν=0.0594_{-0.0007}^{+0.0113}$ eV. However, the cosmological upper limit and the oscillation-based lower limit show a mild $2.2σ$ tension, providing an independent indication of possible physics beyond $Λ$CDM.

preprint2026arXiv

VENUS: Two Faint Little Red Dots Separated by $\sim70\,\mathrm{pc}$ Hidden in a Single Lensed Galaxy at $z\sim7$

We report the identification of a pair of faint little red dots (LRDs), dubbed Red Eyes, in a strongly-lensed galaxy at $z\sim7$ behind the PLCKG004.5-10.5 cluster, identified from the JWST Treasury program VENUS. Red Eyes are spatially resolved on the image plane with distinct colors, while the critical curve lies far north of Red Eyes, clearly requiring two different LRDs rather than a single LRD. Red Eyes is an extremely close pair of LRDs separated by $\sim70\,\mathrm{pc}$ in the source plane with a magnification of $μ\sim20$, which consistently explains another counter-image detected to the north-west. Red Eyes is hosted in a typical star-forming galaxy with $M_{\mathrm{UV,int}}\sim -19$, but its own UV emission is very faint ($M_{\mathrm{UV,int}} \gtrsim -16$). Moreover, Red Eyes does not reside at the galaxy center but lies at an offset position of approximately one effective radius $R_{\mathrm{e}}$ away from the galaxy center. If observed without lensing, Red Eyes would appear as a typical star-forming galaxy at $z\sim 7$ with $M_{\mathrm{UV}}\sim -19$, showing no apparent LRD signatures in either morphology or SED. These results suggest that multiple off-center LRDs, similar to Red Eyes, may be commonly hidden in a typical high-$z$ star-forming galaxy. In this case, various plausible scenarios may emerge, one of which is that intermediate-mass black holes (IMBHs) with $M_\mathrm{BH}\sim10^{4\text{--}6}\,M_\odot$ may form in star clusters on a stellar disk and contribute to the growth of the central supermassive black hole via mergers, with some IMBHs detectable as luminous LRDs in a sufficiently active and massive phase.

preprint2024arXiv

Census for the Rest-frame Optical and UV Morphologies of Galaxies at $z=4-10$: First Phase of Inside-Out Galaxy Formation

We present the rest-frame optical and UV surface brightness (SB) profiles for $149$ galaxies with $M_{\rm opt}< -19.4$ mag at $z=4$-$10$ ($29$ of which are spectroscopically confirmed with JWST NIRSpec), securing high signal-to-noise ratios of $10$-$135$ with deep JWST NIRCam $1$-$5μ$m images obtained by the CEERS survey. We derive morphologies of our high-$z$ galaxies, carefully evaluating the systematics of SB profile measurements with Monte Carlo simulations as well as the impacts of a) AGNs, b) multiple clumps including galaxy mergers, c) spatial resolution differences with previous HST studies, and d) strong emission lines, e.g., H$α$ and [OIII], on optical morphologies with medium-band F410M images. Conducting Sérsic profile fitting to our high-$z$ galaxy SBs with GALFIT, we obtain the effective radii of optical $r_{\rm e, opt}$ and UV $r_{\rm e, UV}$ wavelengths ranging $r_{\rm e, opt}=0.05$-$1.6$ kpc and $r_{\rm e, UV}=0.03$-$1.7$ kpc that are consistent with previous results within large scatters in the size luminosity relations. However, we find the effective radius ratio, $r_{\rm e, opt}/r_{\rm e, UV}$, is almost unity, $1.01^{+0.35}_{-0.22}$, over $z=4$-$10$ with no signatures of past inside-out star formation such found at $z\sim 0$-$2$. There are no spatial offsets exceeding $3σ$ between the optical and UV morphology centers in case of no mergers, indicative of major star-forming activity only found near a mass center of galaxies at $z\gtrsim 4$ probably experiencing the first phase of inside-out galaxy formation.

preprint2022arXiv

A Search for H-Dropout Lyman Break Galaxies at z~12-16

We present two bright galaxy candidates at z~12-13 identified in our H-dropout Lyman break selection with 2.3 deg2 near-infrared deep imaging data. These galaxy candidates, selected after careful screening of foreground interlopers, have spectral energy distributions showing a sharp discontinuity around 1.7 um, a flat continuum at 2-5 um, and non-detections at <1.2 um in the available photometric datasets, all of which are consistent with z>12 galaxy. An ALMA program targeting one of the candidates shows a tentative 4sigma [OIII]88um line at z=13.27, in agreement with its photometric redshift estimate. The number density of the z~12-13 candidates is comparable to that of bright z~10 galaxies, and is consistent with a recently proposed double power-law luminosity function rather than the Schechter function, indicating little evolution in the abundance of bright galaxies from z~4 to 13. Comparisons with theoretical models show that the models cannot reproduce the bright end of rest-frame ultraviolet luminosity functions at z~10-13. Combined with recent studies reporting similarly bright galaxies at z~9-11 and mature stellar populations at z~6-9, our results indicate the existence of a number of star-forming galaxies at z>10, which will be detected with upcoming space missions such as James Webb Space Telescope, Nancy Grace Roman Space Telescope, and GREX-PLUS.

preprint2022arXiv

Are the Newly-Discovered $z \sim 13$ Drop-out Sources Starburst Galaxies or Quasars?

The detection of two $z\sim 13$ galaxy candidates has opened a new window on galaxy formation at an era only $330$ Myr after the Big Bang. Here, we investigate the physical nature of these sources: are we witnessing star forming galaxies or quasars at such early epochs? If powered by star formation, the observed ultraviolet (UV) luminosities and number densities can be jointly explained if: (i) these galaxies are extreme star-formers with star formation rates $5-24\times$ higher than those expected from extrapolations of average lower-redshift relations; (ii) the star formation efficiency increases with halo mass and is countered by increasing dust attenuation from $z \sim 10-5$; (iii) they form stars with an extremely top-heavy initial mass function. The quasar hypothesis is also plausible, with the UV luminosity produced by black holes of $\sim 10^8 \, \mathrm{M_{\odot}}$ accreting at or slightly above the Eddington rate ($f_{\rm Edd}\sim 1.0$). This black hole mass at $z\sim13$ would require very challenging, but not implausible, growth parameters. If spectroscopically confirmed, these two sources will represent a remarkable laboratory to study the Universe at previously inaccessible redshifts.

preprint2022arXiv

Bridging Optical and Far-Infrared Emission-Line Diagrams of Galaxies from Local to the Epoch of Reionization: Characteristic High [O III] 88 $\mathrm{μm}$/SFR at $z > 6$

We present photoionization modeling of galaxy populations at $z\sim0$, $2$, and $> 6$ to bridge optical and far-infrared (FIR) emission-line diagrams. We collect galaxies with measurements of optical and/or FIR ([O III] 88 $\mathrm{μm}$ and [C II] 158 $\mathrm{μm}$) emission line fluxes and plot them on the [O III]$\lambda5007/\mathrm{Hβ}$--[N II]$\lambda6585/\mathrm{Hα}$ (BPT) and L([O III]88)/SFR--L([C II]158)/SFR diagrams, where SFR is the star-formation rate and L([O III]88) and L([C II]158) are the FIR line luminosities. We aim to explain the galaxy distributions on the two diagrams with photoionization models that employ three nebular parameters: the ionization parameter $U$, hydrogen density $n_\text{H}$, and gaseous metallicity $Z_\text{gas}$. Our models successfully reproduce the nebular parameters of local galaxies, and then predict the distributions of the $z\sim0$, $2$, and $> 6$ galaxies on the diagrams. The predicted distributions illustrate the redshift evolution on all the diagrams; e.g., [O III]$/\mathrm{Hβ}$ and [O III]88/[C II]158 ratios continuously decrease from $z > 6$ to $0$. Specifically, the $z > 6$ galaxies exhibit $\sim\!0.5$ dex higher $U$ than low-redshift galaxies at a given $Z_\text{gas}$ and show predicted flat distributions on the BPT diagram at $\log{\mathrm{[O III]/Hβ}} = 0.5$-$0.8$. We find that some of the $z > 6$ galaxies exhibit high L([O III]88)/SFR ratios. To explain these high ratios, our photoionization models require a low stellar-to-gaseous metallicity ratio or bursty/increasing star-formation history at $z > 6$. The James Webb Space Telescope will test the predictions and scenarios for the $z > 6$ galaxies proposed by our photoionization modeling.

preprint2022arXiv

CHORUS IV: Mapping the Spatially Inhomogeneous Cosmic Reionization with Subaru HSC

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

preprint2022arXiv

EMPRESS. V. Metallicity Diagnostics of Galaxies over 12+log(O/H)=~6.9-8.9 Established by a Local Galaxy Census: Preparing for JWST Spectroscopy

We present optical-line gas metallicity diagnostics established by the combination of local SDSS galaxies and the largest compilation of extremely metal-poor galaxies (EMPGs) including new EMPGs identified by the Subaru EMPRESS survey. A total of 103 EMPGs are included that cover a large parameter space of magnitude (Mi=-19 to -7) and H-beta equivalent width (10-600 Ang), i.e., wide ranges of stellar mass and star-formation rate. Using reliable metallicity measurements of the direct method for these galaxies, we derive the relationships between strong optical-line ratios and gas-phase metallicity over the range of 12+log(O/H)=~6.9-8.9 corresponding to 0.02-2 solar metallicity Zsun. We confirm that R23-index, ([OIII]+[OII])/H-beta, is the most accurate metallicity indicator with the metallicity uncertainty of 0.14 dex over the range among various popular metallicity indicators. The other metallicity indicators show large scatters in the metal-poor range (<0.1 Zsun). It is explained by our CLOUDY photoionization modeling that, unlike R23-index, the other metallicity indicators do not use a sum of singly and doubly ionized lines and cannot trace both low and high ionization gas. We find that the accuracy of the metallicity indicators is significantly improved, if one uses H-beta equivalent width measurements that tightly correlate with ionization states. In this work, we also present the relation of physical properties with UV-continuum slope beta and ionization production rate xi_ion derived with GALEX data for the EMPGs, and provide local anchors of galaxy properties together with the optical-line metallicity indicators that are available in the form of ASCII table and useful for forthcoming JWST spectroscopic studies.

preprint2022arXiv

EMPRESS. VI. Outflows Investigated in Low-Mass Galaxies with $M_*=10^4-10^7~M_\odot$: Weak Feedback in Low-Mass Galaxies?

We study emission line profiles of 21 nearby low-mass ($M_*=10^4-10^7~M_\odot$) galaxies in deep medium-high resolution spectra taken with Magellan/MagE. These low-mass galaxies are actively star-forming systems with high specific star-formation rates of $\mathrm{sSFR}\sim100-1000~\mathrm{Gyr}^{-1}$ that are well above the star-formation main sequence and its extrapolation. We identify broad-line components of H$α$ and [OIII]$λ5007$ emission in 14 out of the 21 galaxies that cannot be explained by the MagE instrumental profile or the natural broadening of line emission. We conduct double Gaussian profile fitting to the emission of the 14 galaxies, and find that the broad-line components have line widths significantly larger than those of the narrow-line components, indicative of galactic outflows. The board-line components have moderately large line widths of $\sim 100$ km s$^{-1}$. We estimate the maximum outflow velocities $v_\mathrm{max}$ and obtain values of $\simeq 60-200$ km s$^{-1}$, which are found to be comparable to or slightly larger than the escape velocities. Positive correlations of $v_\mathrm{max}$ with star-formation rates, stellar masses, and circular velocities, extend down into this low-mass regime. Broad- to narrow-line flux ratios BNRs are generally found to be smaller than those of massive galaxies. The small $v_\mathrm{max}$ and BNRs suggest that the mass loading factors $η$ can be as small as 0.1 - 1 or below, in contrast to the large $η$ of energy-driven outflows predicted by numerical simulations.

preprint2022arXiv

EMPRESS. VII. Ionizing Spectrum Shapes of Extremely Metal-Poor Galaxies: Uncovering the Origins of Strong HeII and the Impact on Cosmic Reionization

Strong high-ionization lines such as HeII of young galaxies are puzzling at high and low redshift. Although recent studies suggest the existence of non-thermal sources, whether their ionizing spectra can consistently explain multiple major emission lines remains a question. Here we derive the general shapes of the ionizing spectra for three local extremely metal-poor galaxies (EMPGs) that show strong HeII$λ$4686. We parameterize the ionizing spectra composed of a blackbody and power-law radiation mimicking various stellar and non-thermal sources. We use photoionization models for nebulae, and determine seven parameters of the ionizing spectra and nebulae by Markov Chain Monte Carlo methods, carefully avoiding systematics of abundance ratios. We obtain the general shapes of ionizing spectra explaining $\sim 10$ major emission lines within observational errors with smooth connections from observed X-ray and optical continua. We find that an ionizing spectrum of one EMPG has a blackbody-dominated shape, while the others have convex downward shapes at $>13.6$ eV, which indicate a diversity of the ionizing spectrum shapes. We confirm that the convex downward shapes are fundamentally different from ordinary stellar spectrum shapes, and that the spectrum shapes of these galaxies are generally explained by the combination of the stellar and ultra-luminous X-ray sources. Comparisons with stellar synthesis models suggest that the diversity of the spectrum shapes arises from differences in the stellar age. If galaxies at $z\gtrsim 6$ are similar to the EMPGs, high energy ($>54.4$ eV) photons of the non-stellar sources negligibly contribute to cosmic reionization due to relatively weak radiation.

preprint2022arXiv

First Identification of a CMB Lensing Signal Produced by $1.5$ Million Galaxies at $z\sim$4: Constraints on Matter Density Fluctuations at High Redshift

We report the first detection of the dark matter distribution around Lyman break galaxies (LBGs) at high redshift through the Cosmic Microwave Background (CMB) lensing measurements with the public {\it Planck} PR3 $κ$ map. The LBG sample consists of 1,473,106 objects with the median redshift of $z \sim 4$ that are identified in a total area of 305 deg$^2$ observed by the Hyper Suprime-Cam (HSC) Strategic Survey Program (SSP) survey. After careful investigations of systematic uncertainties, such as contamination from foreground galaxies and Cosmic Infrared Background (CIB), we obtain the significant detection of the CMB lensing signal at $5.1σ$ that is dominated by 2-halo term signals of the LBGs. Fitting a simple model consisting of the Navarro-Frenk-White (NFW) profile and the linear-bias model, we obtain the typical halo mass of $2.9^{+9.5}_{-2.5} \times 10^{11} h^{-1} M_\odot$. Combining the CMB lensing and galaxy-galaxy clustering signals on the large scales, we demonstrate the first cosmological analysis at $z\sim4$ that constrains $(Ω_{{\rm m}0}$, $σ_8)$. We find that our constraint on $σ_8$ is roughly consistent with the {\it Planck} best-fit cosmology, while this $σ_8$ constraint is lower than the {\it Planck} cosmology over the $1σ$ level. This study opens up a new window for constraining cosmological parameters at high redshift by the combination of CMB and high-$z$ galaxies as well as studying the interplay between galaxy evolution and larges-scale structure at such high redshift, by upcoming CMB and optical and near-infrared imaging surveys.

preprint2021arXiv

A Wide and Deep Exploration of Radio Galaxies with Subaru HSC (WERGS). VI. Distant Filamentary Structures Pointed by High-z Radio Galaxies at z~4

We present the environmental properties around high-$z$ radio galaxies (HzRGs) at $z\sim4$, which have been poorly investigated because of their rarity. We use the largest samples of HzRGs and $g$-dropout galaxy overdense regions at $z\sim4$, which were constructed from Hyper Suprime-Cam Subaru Strategic Program, to characterize the HzRG environments statistically. We measure the $g$-dropout galaxy overdensities around 21 HzRGs whose rest-frame 1.4 GHz radio luminosities ($L_{1.4\mathrm{GHz}}$) are $10^{26-27}$ W Hz$^{-1}$. We find that the overdensities around the faint HzRGs with $L_{1.4\mathrm{GHz}}\sim10^{26.0-26.5}$ W Hz$^{-1}$ tend to be higher than that of the $g$-dropout galaxies. On the other hand, no significant difference of density environments is found between the luminous HzRGs with $ L_{1.4\mathrm{GHz}}\sim10^{26.5-27.0} $ W Hz$^{-1}$ and the $g$-dropout galaxies. The HzRGs are found to occupy more massive halos than $g$-dropout galaxies through a cross-correlation between the HzRGs and $g$-dropout galaxies. This trend is more pronounced in the faint HzRGs. These results are consistent with a scenario where HzRGs get older and more massive as the radio-luminosity decreases. The HzRGs are expected to trace the progenitors of local cluster halos from their calculated halo mass. In addition, we find that surrounding galaxies tend to distribute along the radio-jet major axis of the HzRGs at angular distances less than $\lesssim500$ physical kpc. Our findings imply the onset of the filamentary structures around the HzRGs at $z\sim4$.

preprint2021arXiv

EMPRESS. IV. Extremely Metal-Poor Galaxies (EMPGs) Including Very Low-Mass Primordial Systems with M*=10^4--10^5 M_sun and 2--3% (O/H)_sun: High (Fe/O) Suggestive of Metal Enrichment by Hypernovae/Pair-Instability Supernovae

We present Keck/LRIS follow-up spectroscopy for 13 photometric candidates of extremely metal poor galaxies (EMPGs) selected by a machine-learning technique applied to the deep (~26 AB mag) optical and wide-area (~500 deg^2) Subaru imaging data in the EMPRESS survey. Nine out of the 13 candidates are EMPGs with an oxygen abundance (O/H) less than ~10% solar value (O/H)_sun, and four sources are contaminants of moderately metal-rich galaxies or no emission-line objects. Notably, two out of the nine EMPGs have extremely-low stellar masses and oxygen abundances of 5*10^4--7*10^5 M_sun and 2--3% (O/H)_sun, respectively. With a sample of five EMPGs with (Fe/O) measurements, two (three) of which are taken from this study (the literature), we confirm that two EMPGs with the lowest (O/H) ratios of ~2% (O/H)_sun show high (Fe/O) ratios of ~0.1, close to the solar abundance ratio. Comparing galaxy chemical enrichment models, we find that the two EMPGs cannot be explained by a scenario of metal-poor gas accretion/episodic star-formation history due to their low (N/O) ratios. We conclude that the two EMPGs can be reproduced by an inclusion of bright hypernovae and/or hypothetical pair-instability supernovae (SNe) preferentially produced in a metal-poor environment. This conclusion implies that primordial galaxies at z~10 could have a high abundance of Fe that is not originated from Type Ia SNe with delays, and that Fe may not serve as a cosmic clock for primordial galaxies.

preprint2021arXiv

First HETDEX Spectroscopic Determinations of Ly$α$ and UV Luminosity Functions at $z=2-3$: Bridging a Gap Between Faint AGN and Bright Galaxies

We present Ly$α$ 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$α$ emission, $2126$ of which host type 1 AGN showing broad (FWHM$~>1000$ km s$^{-1}$) Ly$α$ emission lines. We derive the Ly$α$ (UV) LF over 2 orders of magnitude covering bright galaxies and AGN in $\log L_\mathrm{Lyα}/\mathrm{[erg~s^{-1}]}=43.3-45.5$ ($-27<M_\mathrm{UV}<-20$) by the $1/V_\mathrm{max}$ estimator. Our results reveal the bright-end hump of the Ly$α$ LF is composed of type 1 AGN. In conjunction with previous spectroscopic results at the faint end, we measure a slope of the best-fit Schechter function to be $α_\mathrm{Sch}=-1.70^{+0.13}_{-0.14}$, which indicates $α_\mathrm{Sch}$ steepens from $z=2-3$ towards high redshift. Our UV LF agrees well with previous AGN UV LFs, and extends to faint-AGN and bright-galaxy regimes. The number fraction of Ly$α$-emitting objects ($X_\mathrm{LAE}$) increases from $M_\mathrm{UV}^*\sim-21$ to bright magnitude due to the contribution of type 1 AGN, while previous studies claim that $X_\mathrm{Lyα}$ decreases from faint magnitude to $M_\mathrm{UV}^*$, suggesting a valley in the $X_\mathrm{Lyα}-$magnitude relation at $M_\mathrm{UV}^*$. Comparing our UV LF of type 1 AGN at $z=2-3$ with those at $z=0$, we find that the number density of faint ($M_\mathrm{UV}>-21$) type 1 AGN increases from $z\sim2$ to $z\sim0$ as opposed to the evolution of bright ($M_\mathrm{UV}<-21$) type 1 AGN, suggesting the AGN downsizing in the rest-frame UV luminosity.

preprint2021arXiv

Galaxy Morphologies Revealed with Subaru HSC and Super-Resolution Techniques I: Major Merger Fractions of L_UV~3-15 L_UV* Dropout Galaxies at z~4-7

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&#34; to <~0.&#34;1, which is comparable to that of the Hubble Space Telescope, allowing us to identify z~4-7 bright major mergers at a high completeness value of >~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<~M_UV<~-22. The major merger fractions are estimated to be f_merger~10-20% at z~4 and ~50-70% at z~5-7, which shows no f_merger difference compared to those of a control faint galaxy sample. Based on the f_merger estimates, we verify contributions of source blending effects and major mergers to the bright-end of double power-law (DPL) shape of z~4-7 galaxy UV LFs. While these two effects partly explain the DPL shape at L_UV~3-10 L_UV*, the DPL shape cannot be explained at the very bright end of L_UV>~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.

preprint2021arXiv

GOLDRUSH. IV. Luminosity Functions and Clustering Revealed with ~4,000,000 Galaxies at z~2-7: Galaxy-AGN Transition, Star Formation Efficiency, and Implication for Evolution at z>10

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<-24 mag and the galaxy 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 z<5 within 0.3 dex, which quantitatively reproduces the redshift evolution of the cosmic SFR density, suggesting that the evolution is primarily driven by the increase of the halo number density due to the structure formation, and the decrease of the accretion rate due to the cosmic expansion. Extrapolating this calculation to higher redshifts assuming the constant efficiency suggests a rapid decrease of the SFR density at z>10 with $\propto10^{-0.5(1+z)}$, which will be directly tested with JWST.

preprint2021arXiv

SILVERRUSH XI: Constraints on the Ly$α$ luminosity function and cosmic reionization at $z=7.3$ with Subaru/Hyper Suprime-Cam

The Ly$α$ luminosity function (LF) of Ly$α$ emitters (LAEs) has been used to constrain the neutral hydrogen fraction in the intergalactic medium (IGM) and thus the timeline of cosmic reionization. Here we present the results of a new narrow-band imaging survey for $z=7.3$ LAEs in a large area of $\sim 3\ \mathrm{deg}^2$ with Subaru/Hyper Suprime-Cam. No LAEs are detected down to $L_{\mathrm{Ly}α}\simeq 10^{43.2}\ \mathrm{erg\ s^{-1}}$ in an effective cosmic volume of $\sim 2\times 10^6$ Mpc$^3$, placing an upper limit to the bright part of the $z=7.3$ Ly$α$ LF for the first time and confirming a decrease in bright LAEs from $z=7.0$. By comparing this upper limit with the Ly$α$ LF in the case of the fully ionized IGM, which is predicted using an observed $z=5.7$ Ly$α$ LF on the assumption that the intrinsic Ly$α$ LF evolves in the same way as the UV LF, we obtain the relative IGM transmission $T^\mathrm{IGM}_{\mathrm{Ly}α}(7.3)/T^\mathrm{IGM}_{\mathrm{Ly}α}(5.7)<0.77$, and then the volume-averaged neutral fraction $x_\mathrm{HI}(7.3)>0.28$. Cosmic reionization is thus still ongoing at $z=7.3$, being consistent with results from other $x_\mathrm{HI}$ estimation methods. A similar analysis using literature Ly$α$ LFs finds that at $z=6.6$ and 7.0 the observed Ly$α$ LF agrees with the predicted one, consistent with full ionization.

preprint2021arXiv

SILVERRUSH. XI. Intensity Mapping for Lya Emission Extending over $100-1000$ comoving kpc around $z\sim2-7$ LAEs with Subaru HSC-SSP and CHORUS Data

We conduct intensity mapping to probe for extended diffuse Ly$α$ emission around Ly$α$ emitters (LAEs) at $z\sim2-7$, exploiting very deep ($\sim26$ mag at $5σ$) and large-area ($\sim4.5$ deg$^2$) Subaru/Hyper Suprime-Cam narrow-band (NB) images and large LAE catalogs consisting of a total of 1781 LAEs at $z=2.2$, $3.3$, $5.7$, and $6.6$ obtained by the HSC-SSP SILVERRUSH and CHORUS projects. We calculate the spatial correlations of these LAEs with $\sim1-2$ billion pixel flux values of the NB images, deriving the average Ly$α$ surface brightness (${\rm SB_{Lyα}}$) radial profiles around the LAEs. By carefully estimating systematics such as fluctuations of sky background and point spread functions, we detect diffuse Ly$α$ emission ($\sim10^{-20}-10^{-19}$ erg s$^{-1}$ cm$^{-2}$ arcsec$^{-2}$) at $100-1000$ comoving kpc around $z=3.3$ LAEs at the $4.1σ$ level and tentatively ($\sim2σ$) at the other redshifts, beyond the virial radius of a dark-matter halo with a mass of $10^{11}\ M_\odot$. While the observed ${\rm SB_{Lyα}}$ profiles have similar amplitudes at $z=2.2-6.6$ within the uncertainties, the intrinsic ${\rm SB_{Lyα}}$ 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 a LAE is powered by resonantly scattered Ly$α$ photons in the CGM and IGM that originates from the inner part of the LAE, and/or neighboring galaxies around the LAE.

preprint2021arXiv

Subaru High-z Exploration of Low-Luminosity Quasars (SHELLQs) XII. Extended [C II] Structure (Merger or Outflow) in a z = 6.72 Red Quasar

We present ALMA [C II] 158 $μ$m line and far-infrared (FIR) continuum emission observations toward HSC J120505.09$-$000027.9 (J1205$-$0000) at $z = 6.72$ with the beam size of $\sim 0&#39;&#39;.8 \times 0&#39;&#39;.5$ (or 4.1 kpc $\times$ 2.6 kpc), the most distant red quasar known to date. Red quasars are modestly reddened by dust, and are thought to be in rapid transition from an obscured starburst to an unobscured normal quasar, driven by powerful active galactic nucleus (AGN) feedback which blows out a cocoon of interstellar medium (ISM). The FIR continuum of J1205$-$0000 is bright, with an estimated luminosity of $L_{\rm FIR} \sim 3 \times 10^{12}~L_\odot$. The [C II] line emission is extended on scales of $r \sim 5$ kpc, greater than the FIR continuum. The line profiles at the extended regions are complex and broad (FWHM $\sim 630-780$ km s$^{-1}$). Although it is not practical to identify the nature of this extended structure, possible explanations include (i) companion/merging galaxies and (ii) massive AGN-driven outflows. For the case of (i), the companions are modestly star-forming ($\sim 10~M_\odot$ yr$^{-1}$), but are not detected by our Subaru optical observations ($y_{\rm AB,5σ} = 24.4$ mag). For the case of (ii), our lower-limit to the cold neutral outflow rate is $\sim 100~M_\odot$ yr$^{-1}$. The outflow kinetic energy and momentum are both much smaller than what predicted in energy-conserving wind models, suggesting that the AGN feedback in this quasar is not capable of completely suppressing its star formation.

preprint2021arXiv

Subaru High-z Exploration of Low-Luminosity Quasars (SHELLQs). XVI. 69 New Quasars at 5.8 < z < 7.0

We present the spectroscopic discovery of 69 quasars at 5.8 < z < 7.0, drawn from the Hyper Suprime-Cam (HSC) Subaru Strategic Program (SSP) imaging survey data. This is the 16th publication from the Subaru High-z Exploration of Low-Luminosity Quasars (SHELLQs) project, and completes identification of all but the faintest candidates (i.e., i-band dropouts with zAB < 24 and y-band detections, and z-band dropouts with yAB < 24) with Bayesian quasar probability Pq > 0.1 in the HSC-SSP third public data release (PDR3). The sample reported here also includes three quasars with Pq < 0.1 at z ~ 6.6, which we selected in an effort to completely cover the reddest point sources with simple color cuts. The number of high-z quasars discovered in SHELLQs has now grown to 162, including 23 type-II quasar candidates. This paper also presents identification of seven galaxies at 5.6 < z < 6.7, an [O III] emitter at z = 0.954, and 31 Galactic cool stars and brown dwarfs. High-z quasars and galaxies comprise 75 % and 16 % respectively of all the spectroscopic SHELLQs objects that pass our latest selection algorithm with the PDR3 photometry. That is, a total of 91 % of the objects lie at z > 5.6. This demonstrates that the algorithm has very high efficiency, even though we are probing an unprecedentedly low-luminosity population down to M1450 ~ -21 mag.

preprint2021arXiv

Third Data Release of the Hyper Suprime-Cam Subaru Strategic Program

The paper presents the third data release of Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP), a wide-field multi-band imaging survey with the Subaru 8.2m telescope. HSC-SSP has three survey layers (Wide, Deep, and UltraDeep) with different area coverages and depths, designed to address a wide array of astrophysical questions. This third release from HSC-SSP includes data from 278 nights of observing time and covers about 670 square degrees in all five broad-band filters at the full depth ($\sim26$~mag at $5σ$) in the Wide layer. If we include partially observed area, the release covers 1,470 square degrees. The Deep and UltraDeep layers have $\sim80\%$ of the originally planned integration times, and are considered done, as we have slightly changed the observing strategy in order to compensate for various time losses. There are a number of updates in the image processing pipeline. Of particular importance is the change in the sky subtraction algorithm; we subtract the sky on small scales before the detection and measurement stages, which has significantly reduced false detections. Thanks to this and other updates, the overall quality of the processed data has improved since the previous release. However, there are limitations in the data (for example, the pipeline is not optimized for crowded fields), and we encourage the user to check the quality assurance plots as well as a list of known issues before exploiting the data. The data release website is https://hsc-release.mtk.nao.ac.jp/.

preprint2020arXiv

3D Distribution Map of HI Gas and Galaxies Around an Enormous Ly$α$ Nebula and Three QSOs at $z=2.3$ Revealed by the HI Tomographic Mapping Technique

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

preprint2020arXiv

Faint Quasars Live in the Same Number Density Environments as Lyman Break Galaxies at z ~ 4

Characterizing high-z quasar environments is key to understanding the co-evolution of quasars and the surrounding galaxies. To restrict their global picture, we statistically examine the g-dropout galaxy overdensity distribution around 570 faint quasar candidates at z ~ 4, based on the Hyper Suprime-Cam Subaru Strategic Program survey. We compare the overdensity significances of g-dropout galaxies around the quasars with those around g-dropout galaxies, and find no significant difference between their distributions. A total of 4 (22) out of the 570 faint quasars, 0.7_{-0.4}^{+0.4} (3.9_{-0.8}^{+0.8}) %, are found to be associated with the > 4 sigma overdense regions within an angular separation of 1.8 (3.0) arcmin, which is the typical size of protoclusters at this epoch. This is similar to the fraction of g-dropout galaxies associated with the > 4 sigma overdense regions. This result is consistent with our previous work that 1.3_{-0.9}^{+0.9} % and 2.0_{-1.1}^{+1.1} % of luminous quasars detected in the Sloan Digital Sky Survey exist in the > 4 sigma overdense regions within 1.8 and 3.0 arcmin separations, respectively. Therefore, we suggest that the galaxy number densities around quasars are independent of their luminosity, and most quasars do not preferentially appear in the richest protocluster regions at z ~ 4. The lack of an apparent positive correlation between the quasars and the protoclusters implies that: i) the gas-rich major merger rate is relatively low in the protocluster regions, ii) most high-z quasars may appear through secular processes, or iii) some dust-obscured quasars exist in the protocluster regions.

preprint2020arXiv

FIR-luminous [CII] emitters in the ALMA-SCUBA-2 COSMOS survey (AS2COSMOS): The nature of submillimeter galaxies in a 10 comoving Mpc-scale structure at z~4.6

We report the discovery of a 10 comoving Mpc-scale structure traced by massive submillimeter galaxies (SMGs) at z~4.6. These galaxies are selected from an emission line search of ALMA Band 7 observations targeting 184 luminous submillimeter sources ($S_{850μ{\rm m}}\geq$ 6.2 mJy) across 1.6 degrees$^2$ in the COSMOS field. We identify four [CII] emitting SMGs and two probable [CII] emitting SMG candidates at z=4.60-4.64 with velocity-integrated signal-to-noise ratio of SNR>8. Four of the six emitters are near-infrared blank SMGs. After excluding one SMG whose emission line is falling at the edge of the spectral window, all galaxies show clear velocity gradients along the major axes that are consistent with rotating gas disks. The estimated rotation velocities of the disks are 330-550 km s$^{-1}$ and the inferred host dark-matter halo masses are ~2-8 $\times$ 10$^{12}$M$_{\odot}$. From their estimated halo masses and [CII] luminosity function, we suggest that these galaxies have a high (50-100%) duty cycle and high (~0.1) baryon conversion efficiency (SFR relative to baryon accretion rate), and that they contribute $\simeq$2% to the total star-formation rate density at z=4.6. These SMGs are concentrated within just 0.3% of the full survey volume, suggesting they are strongly clustered. The extent of this structure and the individual halo masses suggest that these SMGs will likely evolve into members of a ~10$^{15}$M$_{\odot}$ cluster at z=0. This survey reveals synchronized dusty starburst in massive halos at z>4, which could be driven by mergers or fed by smooth gas accretion.

preprint2020arXiv

Large Population of ALMA Galaxies at z>6 with Very High [OIII]88um to [CII]158um Flux Ratios: Evidence of Extremely High Ionization Parameter or PDR Deficit?

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.

preprint2020arXiv

Wide and Deep Exploration of Radio Galaxies with Subaru HSC (WERGS). III. Discovery of a z = 4.72 Radio Galaxy with Lyman Break Technique

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$α$ 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$α$ 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: $α^{325}_{1400}$ of $-1.1$ and $α^{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.

preprint2019arXiv

Balmer Break Galaxy Candidates at $z \sim 6$: a Potential View on the Star-Formation Activity at $z \gtrsim 14$

We search for galaxies with a strong Balmer break (Balmer Break Galaxies; BBGs) at $z \sim 6$ over a 0.41 deg$^2$ effective area in the COSMOS field. Based on rich imaging data, including data obtained with the Atacama Large Millimeter/submillimeter Array (ALMA), three candidates are identified by their extremely red $K - [3.6]$ colors as well as by non-detection in X-ray, optical, far-infrared (FIR), and radio bands. The non-detection in the deep ALMA observations suggests that they are not dusty galaxies but BBGs at $z \sim 6$, although contamination from Active Galactic Nuclei (AGNs) at $z \sim 0$ cannot be completely ruled out for the moment. Our spectral energy distribution (SED) analyses reveal that the BBG candidates at $z \sim 6$ have stellar masses of $\approx 5 \times 10^{10} M_{\odot}$ dominated by old stellar populations with ages of $\gtrsim 700$ Myr. Assuming that all the three candidates are real BBGs at $z \sim 6$, we estimate the stellar mass density (SMD) to be $2.4^{+2.3}_{-1.3} \times 10^{4} M_{\odot}$ Mpc$^{-3}$. This is consistent with an extrapolation from the lower redshift measurements. The onset of star formation in the three BBG candidates is expected to be several hundred million years before the observed epoch of $z \sim 6$. We estimate the star-formation rate density (SFRD) contributed by progenitors of the BBGs to be 2.4 -- 12 $\times 10^{-5} M_{\odot}$ yr$^{-1} $Mpc$^{-3}$ at $z > 14$ (99.7\% confidence range). Our result suggests a smooth evolution of the SFRD beyond $z = 8$.

preprint2019arXiv

CHORUS. III. Photometric and Spectroscopic Properties of Ly$α$ Blobs at $z=4.9-7.0$

We report the Subaru Hyper Suprime-Cam (HSC) discovery of two Ly$α$ blobs (LABs), dubbed z70-1 and z49-1 at $z=6.965$ and $z=4.888$ respectively, that are Ly$α$ emitters with a bright ($\log L_{\rm Lyα}/{\rm [erg\ s^{-1}]}>43.4$) and spatially-extended Ly$α$ 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$α$ 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}$λ$1240 emission, or Ly$α$ line broadening, while z49-1 has a strong {\sc Civ}$λ$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$α$ emission of the LABs. The Ly$α$ 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$α$ luminosity relation of the Ly$α$ 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.

preprint2019arXiv

Early Low-Mass Galaxies and Star-Cluster Candidates at z~6-9 Identified by the Gravitational Lensing Technique and Deep Optical/Near-Infrared Imaging

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

preprint2019arXiv

Subaru High-z Exploration of Low-Luminosity Quasars (SHELLQs) VIII. A less biased view of the early co-evolution of black holes and host galaxies

We present ALMA [CII] line and far-infrared (FIR) continuum observations of three $z > 6$ low-luminosity quasars ($M_{\rm 1450} > -25$) discovered by our Subaru Hyper Suprime-Cam (HSC) survey. The [CII] line was detected in all three targets with luminosities of $(2.4 - 9.5) \times 10^8~L_\odot$, about one order of magnitude smaller than optically luminous ($M_{\rm 1450} \lesssim -25$) quasars. The FIR continuum luminosities range from $< 9 \times 10^{10}~L_\odot$ (3$σ$ limit) to $\sim 2 \times 10^{12}~L_\odot$, indicating a wide range in star formation rates in these galaxies. Most of the HSC quasars studied thus far show [CII]/FIR luminosity ratios similar to local star-forming galaxies. Using the [CII]-based dynamical mass ($M_{\rm dyn}$) as a surrogate for bulge stellar mass ($M_{\rm bulge}$), we find that a significant fraction of low-luminosity quasars are located on or even below the local $M_{\rm BH} - M_{\rm bulge}$ relation, particularly at the massive end of the galaxy mass distribution. In contrast, previous studies of optically luminous quasars have found that black holes are overmassive relative to the local relation. Given the low luminosities of our targets, we are exploring the nature of the early co-evolution of supermassive black holes and their hosts in a less biased way. Almost all of the quasars presented in this work are growing their black hole mass at much higher pace at $z \sim 6$ than the parallel growth model, in which supermassive black holes and their hosts grow simultaneously to match the local $M_{\rm BH} - M_{\rm bulge}$ relation at all redshifts. As the low-luminosity quasars appear to realize the local co-evolutionary relation even at $z \sim 6$, they should have experienced vigorous starbursts prior to the currently observed quasar phase to catch up with the relation.

preprint2018arXiv

Morphologies of ~190,000 Galaxies at z=0-10 Revealed with HST Legacy Data. III. Continuum Profile and Size Evolution of Lya Emitters

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.

preprint2018arXiv

The onset of star formation 250 million years after the Big Bang

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