Search

Using the Subaru/FOCAS IFU capability, we examine the spatially resolved relationships between gas-phase metallicity, stellar mass, and star-formation rate surface densities (Sigma_* and Sigma_SFR, respectively) in extremely metal-poor galaxies (EMPGs) in the local universe. Our analysis includes 24 EMPGs, comprising 9,177 spaxels, which span a unique parameter space of local metallicity (12+log(O/H) = 6.9 to 7.9) and stellar mass surface density (Sigma_* ~ 10^5 to 10^7 Msun/kpc^2), extending beyond the range of existing large integral-field spectroscopic surveys. Through spatially resolved emission line diagnostics based on the [NII] BPT-diagram, we verify the absence of evolved active galactic nuclei in these EMPGs. Our findings reveal that, while the resolved mass-metallicity relation exhibits significant scatter in the low-mass regime, this scatter is closely correlated with local star-formation surface density. Specifically, metallicity decreases as Sigma_SFR increases for a given Sigma_*. Notably, half of the EMPGs show a distinct metal-poor horizontal branch on the resolved mass-metallicity relation. This feature typically appears at the peak clump with the highest Sigma_* and Sigma_SFR and is surrounded by a relatively metal-enriched ambient region. These findings support a scenario in which metal-poor gas infall fuels episodic star formation in EMPGs, consistent with the kinematic properties observed in these systems. In addition, we identify four EMPGs with exceptionally low central metallicities (12+log(O/H) <~ 7.2), which display only a metal-poor clump without a surrounding metal-rich region. This suggests that such ultra-low metallicity EMPGs, at less than a few percent of the solar metallicity, may serve as valuable analogs for galaxies in the early stages of galaxy evolution., Comment: 29 pages, 14 figures. Submitted to ApJ

We investigate UV continuum slopes $\beta$ of 974 galaxies at $z=4-14$ using archival JWST/NIRSpec PRISM spectra obtained from major JWST GTO, ERS, and GO programs, including JADES, CEERS, and UNCOVER. Among these galaxies, we identify a remarkable galaxy at $z=9.25$, dubbed EBG-1, with a significantly blue UV slope $\beta=-2.99\pm0.15$, unlike the rest of the galaxies that exhibit red continua or ambiguous blue continua hindered by large uncertainties. We confirm that the $\beta$ value negligibly changes by the data reduction and fitting wavelength ranges for UV emission/absorption line masking. The extreme blue slope, $\beta=-3.0$, rules out significant contributions from dust extinction or AGN activity. Comparing with stellar and nebular emission models, we find that such a blue UV slope cannot be reproduced solely by stellar models even with very young, metal-poor, or top-heavy contiguous star formation associated with strong nebular continua making the UV slopes red, but with a high ionizing photon escape fraction, $f_\mathrm{esc}^\mathrm{ion} \gtrsim 0.5$, for a weak nebular continuum. While the H$\beta$ emission line is not detected, likely due to the limited sensitivity of the spectrum, we find moderately weak [O III] $\lambda\lambda$4959,5007 emission lines for the given star-formation rate ($3\, \mathrm{M_\odot}$ yr$^{-1}$) and stellar mass ($10^{8.0} \, \mathrm{M_\odot}$) that are about three times weaker than the average emission lines, again suggestive of the high ionizing photon escape fraction, $f_\mathrm{esc}^\mathrm{ion} \sim 0.7$ or more. EBG-1 would provide crucial insights into stellar and nebular continuum emission in high-redshift galaxies, serving as an example of the ionizing photon escaping site at the epoch of reionization., Comment: 9 pages, 8 figures, submitted to ApJ

High-$z$ radio galaxies (HzRGs) are considered important objects for understanding the formation and evolution of massive galaxies in the early universe. However, till date, detailed studies of the stellar population of HzRGs such as the star-formation history have been scarce. Therefore, this study conducted a new survey to establish a less-biased sample of HzRGs and consequently investigate their properties. We utilized a sample of $g$-dropout Lyman break galaxies (LBGs) obtained from an optical wide and deep imaging survey made by Subaru Hyper Suprime-Cam (HSC). Based on the cross-matching of this LBG sample with the VLA FIRST radio survey data, we constructed a photometric sample of high-redshift radio galaxies (HzRGs) at $z \sim 4$ for $\sim$560 deg$^2$ survey field. Consequently, we identified 146 HzRG candidates. To analyze the characteristics of these candidates, we focus on objects exhibiting the near-infrared photometry of VIKING or UKIDSS and the mid-infrared photometry of unWISE (28 objects). The results indicate that 7 objects exhibit SEDs consistent with galaxies at $z \sim 4$. The HzRG candidates have very large stellar masses with $\sim 4.2 \times 10^{11} M_{\odot}$ on average. This stellar mass is similar to that of previously discovered USS HzRGs at $z \sim 4$, though our sample is affected by a sample selection bias that selects only HzRGs with $M_{\star} > 10^{11} M_{\odot}$. Further, the SEDs of those HzRG candidates suggest a past fast quenching with a rough timescale of $\sim$0.1 Gyr, as evidenced from the rest-frame UVJ diagram., Comment: 19 pages, 13 figures

We present the luminosity functions (LFs) and angular correlation functions (ACFs) derived from 18,960 Ly$\alpha$ emitters (LAEs) at $z=2.2-7.3$ over a wide survey area of $\lesssim24 {\rm deg^2}$ that are identified in the narrowband data of the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) and the Cosmic HydrOgen Reionization Unveiled with Subaru (CHORUS) surveys. Confirming the large sample with the 241 spectroscopically identified LAEs, we determine Ly$\alpha$ LFs and ACFs in the brighter luminosity range down to $0.5L_{\star}$, and confirm that our measurements are consistent with previous studies but offer significantly reduced statistical uncertainties. The improved precision of our ACFs allows us to clearly detect one-halo terms at some redshifts, and provides large-scale bias measurements that indicate hosting halo masses of $\sim 10^{11} M_\odot$ over $z\simeq 2-7$. By comparing our Ly$\alpha$ LF (ACF) measurements with reionization models, we estimate the neutral hydrogen fractions in the intergalactic medium to be $x_{\rm \HI}<0.05$ (=${0.06}^{+0.12}_{-0.03}$), $0.15^{+0.10}_{-0.08}$ (${0.21}^{+0.19}_{-0.14}$), $0.18^{+0.14}_{-0.12}$, and $0.75^{+0.09}_{-0.13}$ at $z=5.7$, $6.6$, $7.0$, and $7.3$, respectively. Our findings suggest that the neutral hydrogen fraction remains relatively low, $x_{\rm \HI} \lesssim 0.2$, at $z=5-7$, but increases sharply at $z > 7$, reaching $x_{\rm \HI} \sim 0.9$ by $z \simeq 8-9$, as indicated by recent JWST studies. The combination of our results from LAE observations with recent JWST observations suggests that the major epoch of reionization occurred around $z \sim 7-8$, likely driven by the emergence of massive sources emitting significant ionizing photons., Comment: Submitted to ApJS

We present deep Subaru/FOCAS spectra for two extreme emission line galaxies (EELGs) at $z\sim 1$ with strong {\sc[Oiii]}$\lambda$5007 emission lines, exhibiting equivalent widths (EWs) of $2905^{+946}_{-578}$ \AA\ and $2000^{+188}_{-159}$ \AA, comparable to those of EELGs at high redshifts that are now routinely identified with JWST spectroscopy. Adding a similarly large {\sc [Oiii]} EW ($2508^{+1487}_{-689}$ \AA) EELG found at $z\sim 2$ in the JWST CEERS survey to our sample, we explore for the physical origins of the large {\sc [Oiii]} EWs of these three galaxies with the Subaru spectra and various public data including JWST/NIRSpec, NIRCam, and MIRI data. While there are no clear signatures of AGN identified by the optical line diagnostics, we find that two out of two galaxies covered by the MIRI data show strong near-infrared excess in the spectral energy distributions (SEDs) indicating obscured AGN. Because none of the three galaxies show clear broad H$\beta$ lines, the upper limits on the flux ratios of broad-H$\beta$ to {\sc [Oiii]} lines are small, $\lesssim 0.15$ that are comparable with Seyfert $1.8-2.0$ galaxies. We conduct \texttt{Cloudy} modeling with the stellar and AGN incident spectra, allowing a wide range of parameters including metallicities and ionization parameters. We find that the large {\sc [Oiii]} EWs are not self-consistently reproduced by the spectra of stars or unobscured AGN, but obscured AGN that efficiently produces O$^{++}$ ionizing photons with weak nuclear and stellar continua that are consistent with the SED shapes., Comment: submitted to ApJ

We present an [O/Fe] ratio of a luminous galaxy GN-z11 at $z=10.60$ derived with the deep public JWST/NIRSpec data. We fit the medium-resolution grating (G140M, G235M, and G395M) data with the model spectra consisting of BPASS-stellar and CLOUDY-nebular spectra in the rest-frame UV wavelength ranges with Fe absorption lines, carefully masking out the other emission and absorption lines in the same manner as previous studies conducted for lower redshift ($z\sim 2-6$) galaxies with oxygen abundance measurements. We obtain an Fe-rich abundance ratio $\mathrm{[O/Fe]}=-0.37^{+0.43}_{-0.22}$, which is confirmed with the independent deep prism data as well as by the classic 1978 index method. This [O/Fe] measurement is lower than measured for star-forming galaxies at $z\sim 2-3$. Because $z=10.60$ is an early epoch after the Big Bang ($\sim 430$ Myr) and the first star formation (likely $\sim 200$ Myr), it is difficult to produce Fe by Type Ia supernovae (SNeIa) requiring sufficient delay time for white-dwarf formation and gas accretion. The Fe-rich abundance ratio in GN-z11 suggests that the delay time is short, or that the major Fe enrichment is not accomplished by SNeIa but bright hypernovae (BrHNe) and/or pair-instability supernovae (PISNe), where the yield models of BrHNe and PISNe explain Fe, Ne, and O abundance ratios of GN-z11. The [O/Fe] measurement is not too low to rule out the connection between GN-z11 and globular clusters (GCs) previously suggested by the nitrogen abundance, but rather supports the connection with a GC population at high [N/O] if a metal dilution process exists., Comment: 13 pages, 11 figures, accepted for publication in ApJ

We present the number densities and physical properties of the bright galaxies spectroscopically confirmed at $z\sim7-14$. Our sample is composed of 60 galaxies at $z_\mathrm{spec}\sim7-14$, including recently-confirmed galaxies at $z_\mathrm{spec}=12.34-14.32$ with JWST, as well as new confirmations at $z_\mathrm{spec}=6.583-7.643$ with $-24< M_\mathrm{UV}< -21$ mag using ALMA and Keck. Our JWST/NIRSpec observations have also revealed that very bright galaxy candidates at $z\sim10-13$ identified from ground-based telescope images before JWST are passive galaxies at $z\sim3-4$, emphasizing the necessity of strict screening and spectroscopy in the selection of the brightest galaxies at $z>10$. The UV luminosity functions derived from these spectroscopic results are consistent with a double power-law function, showing tensions with theoretical models at the bright end. To understand the origin of the overabundance of bright galaxies, we investigate their morphologies using JWST/NIRCam high-resolution images obtained in various surveys including PRIMER and COSMOS-Web. We find that $\sim70\%$ of the bright galaxies at $z\sim7$ exhibit clumpy morphologies with multiple sub-components, suggesting merger-induced starburst activity, which is consistent with SED fitting results showing bursty star formation histories. At $z\gtrsim10$, bright galaxies are classified into two types of galaxies; extended ones with weak high-ionization emission lines, and compact ones with strong high-ionization lines including NIV]$\lambda$1486, indicating that at least two different processes (e.g., merger-induced starburst and compact star formation/AGN) are shaping the physical properties of the brightest galaxies at $z\gtrsim10$ and are responsible for their overabundance., Comment: 34 pages, 22 figures, Accepted for Publication in ApJ

We present HeI/H$\beta$-flux and He/H-abundance ratios in three JWST galaxies with significant constraints on N/O-abundance ratios, GS-NDG-9422, RXCJ2248-ID, and GLASS150008 at $z\sim 6$ mostly with the spectroscopic coverage from HeI$\lambda$4471 and HeII$\lambda$4686 to HeI$\lambda$7065, comparing with 68 local-dwarf galaxies. We find that these high-$z$ galaxies present strong HeI emission with HeI/H$\beta$ flux ratios generally larger than those of local-dwarf galaxies. We derive He/H with all of the detected HeI, HeII, and $2-3$ hydrogen Balmer lines in the same manner as the local He/H determination conducted for cosmology studies. These high-$z$ galaxies show He overabundance He/H$\gtrsim 0.10$ or high electron density $n_\mathrm{e}\sim 10^{3-4}$ cm$^{-3}$ much larger than local values at low O/H, $12+\log \mathrm{(O/H)}=7-8$. In contrast, we obtain low He/H and $n_\mathrm{e}$ values for our local-dwarf galaxies by the same technique with the same helium and hydrogen lines, and confirm that the difference between the high-$z$ and local-dwarf galaxies are not mimicked by systematics. While two scenarios of 1) He overabundance and 2) high electron density are not clearly concluded, we find that there is a positive correlation on the He/H-N/O or $n_\mathrm{e}$-N/O plane by the comparison of the high-$z$ and local-dwarf galaxies. The scenario 1) suggests that the overabundant helium and nitrogen are not explained by the standard chemical enrichment of core-collapse supernovae, but the CNO-cycle products and equilibrium ratios, respectively. The scenario 2) indicates that the strong helium lines are originated from the central dense clouds of the high-$z$ galaxies by excessive collisional excitation., Comment: 19 pages, 12 figures, accepted for publication in ApJ

We investigate the dynamics of GN-z11, a luminous galaxy at $z=10.60$, carefully analyzing the public deep integral field spectroscopy (IFS) data taken with JWST NIRSpec IFU. While the observations of the IFS data originally targeted a He II clump near GN-z11, we find that CIII]$\lambda\lambda$1907,1909 emission from ionized gas at GN-z11 is bright and spatially extended significantly beyond the point-spread function (PSF). The spatially extended CIII emission of GN-z11 shows a velocity gradient, red- and blue-shifted components in the north and south directions, respectively, which cannot be explained by the variation of [CIII]$\lambda$1907/CIII]$\lambda$1909 line ratios. Assuming the velocity gradient is produced by disk rotation, we perform forward modeling with GalPak$^{3D}$, including the effects of PSF smearing and line blending, and obtain a rotation velocity of $v_{rot}=257^{+138}_{-117}$ km s$^{-1}$, a velocity dispersion of $\sigma_v=91^{+18}_{-32}$ km s$^{-1}$, and a ratio of $v_{rot}/\sigma_v=2.83^{+1.82}_{-1.41}$. The $v_{rot}/\sigma_v$ value would suggest a rotation-dominated disk existing at $z>10$ albeit with the large uncertainties. The rotation velocity agrees with those of numerical simulations predicting a rotating disk formed in the early universe under the condition of mass compaction and weak feedback. While the velocity gradient is consistent with the rotating disk solution, we recognize that galactic outflows can also explain the velocity gradient as well as the extended morphology and the high velocity dispersion found in the outskirt. Higher S/N and resolution data are necessary to conclude the physical origin of the velocity gradient in GN-z11., Comment: Accepted for publication in ApJ

We present follow-up spectroscopy on a protocluster candidate selected from the wide-field imaging of the Hyper SuprimeCam Subaru Strategic Programme. The target protocluster candidate was identified as a $4.5\sigma$ overdense region of $g$-dropout galaxies, and the redshifts of $g$-dropout galaxies are determined by detecting their Ly$\alpha$ emission. Thirteen galaxies, at least, are found to be clustering in the narrow redshift range of $\Delta z<0.05$ at $z=3.699$. This is clear evidence of the presence of a protocluster in the target region. Following the discovery of the protocluster at $z=3.699$, the physical properties and three-dimensional distribution of its member galaxies are investigated. Based on spectroscopically-confirmed $g$-dropout galaxies, we find an overabundance of rest-frame ultraviolet (UV) bright galaxies in the protocluster. The UV brightest protocluster member turns out to be an active galactic nucleus, and the other UV brighter members tend to show smaller Ly$\alpha$ equivalent widths than field counterparts. The member galaxies tend to densely populate near the centre of the protocluster, but the separation from the nearest neighbour rather than the distance from the centre of the protocluster is more tightly correlated to galaxy properties, implying that the protocluster is still in an early phase of cluster formation and only close neighbours have a significant impact on the physical properties of protocluster members. The number density of massive galaxies, selected from an archival photometric-redshift catalogue, is higher near the centre of the protocluster, while dusty starburst galaxies are distributed on the outskirts. These observational results suggest that the protocluster consists of multiple galaxy populations, whose spatial distributions may hint at the developmental phase of the galaxy cluster., Comment: 11 pages, 8 figures, submitted to MNRAS

We investigate the physical origins of the Balmer decrement anomalies in GS-NDG-9422 (Cameron et al. 2023) and RXCJ2248-ID (Topping et al. 2024) galaxies at $z\sim 6$ whose $\mathrm{H}\alpha/\mathrm{H}\beta$ values are significantly smaller than $2.7$, the latter of which also shows anomalous $\mathrm{H}\gamma/\mathrm{H}\beta$ and $\mathrm{H}\delta/\mathrm{H}\beta$ values beyond the errors. Because the anomalous Balmer decrements are not reproduced under the Case B recombination, we explore the nebulae with the optical depths smaller and larger than the Case B recombination by physical modeling. We find two cases quantitatively explaining the anomalies; 1) density-bounded nebulae that are opaque only up to around Ly$\gamma$-Ly8 transitions and 2) ionization-bounded nebulae partly/fully surrounded by optically-thick excited H{\sc i} clouds. The case of 1) produces more H$\beta$ photons via Ly$\gamma$ absorption in the nebulae, requiring fine tuning in optical depth values, while this case helps ionizing photon escape for cosmic reionization. The case of 2) needs the optically-thick excited H{\sc i} clouds with $N_2\simeq 10^{12}-10^{13}$ $\mathrm{cm^{-2}}$, where $N_2$ is the column density of the hydrogen atom with the principal quantum number of $n=2$. Interestingly, the high $N_2$ values qualitatively agree with the recent claims for GS-NDG-9422 with the strong nebular continuum requiring a number of $2s$-state electrons and for RXCJ2248-ID with the dense ionized regions likely coexisting with the optically-thick clouds. While the physical origin of the optically-thick excited H{\sc i} clouds is unclear, these results may suggest gas clouds with excessive collisional excitation caused by an amount of accretion and supernovae in the high-$z$ galaxies., Comment: 13 pages, 7 figures, accepted for publication in ApJ

We demonstrated the lattice structures and the superconducting phases of metallic double-chain based cuprate Pr2Ba4Cu7O15-delta exhibiting higher Tc. After the oxygen heat treatment on citrate pyrolysis precursors, their reduction treatment followed by a quench procedure caused higher Tc samples with 26.5-30 K. The crystal structural parameters for the superconducting sample (delta = 0.81) were analyzed from the powder synchrotron X-ray diffraction data using RIETAN-FP program. The effect of magnetic field on the superconducting phase of these samples with different oxygen defects (delta =0.73, 0.81 and 0.87) was examined, for our understanding of the superconducting magnetic field-temperature phase diagram. For delta = 0.87 sample with Tc = 30 K, the resistive critical field Hc* was estimated to be 13 T at 4.2 K. The oxygen deficiency dependence on Tc,on for our samples was compared with the data of several other groups., Comment: 8 pages, 6 figures

We present the evolution of Ly$\alpha$ emission derived from 53 galaxies at $z=6.6-13.2$ that are identified by multiple JWST/NIRSpec spectroscopy programs of ERS, GO, DDT, and GTO. These galaxies fall on the star-formation main sequence and are the typical star-forming galaxies with UV magnitudes of $-22.5\leq M_\mathrm{UV}\leq-17.0$. We find that 15 out of 53 galaxies show Ly$\alpha$ emission at the $>3\sigma$ levels, and obtain Ly$\alpha$ equivalent width (EW) measurements and stringent $3\sigma$ upper limits for the 15 and 38 galaxies, respectively. Confirming that Ly$\alpha$ velocity offsets and line widths of our galaxies are comparable with those of low-redshift Ly$\alpha$ emitters, we investigate the redshift evolution of the Ly$\alpha$ EW. We find that Ly$\alpha$ EWs statistically decrease towards high redshifts on the Ly$\alpha$ EW vs. $M_{\rm UV}$ plane for various probability distributions of the uncertainties. We then evaluate neutral hydrogen fractions $x_{\rm HI}$ with the Ly$\alpha$ EW redshift evolution and the cosmic reionization simulation results on the basis of a Bayesian inference framework, and obtain $x_{\rm HI}<0.79$, $=0.62^{+0.15}_{-0.36}$, and $0.93^{+0.04}_{-0.07}$ at $z\sim7$, $8$, and $9-13$, respectively. These moderately large $x_{\rm HI}$ values are consistent with the Planck CMB optical depth measurement and previous $x_{\rm HI}$ constraints from galaxy and QSO Ly$\alpha$ damping wing absorptions, and strongly indicate a late reionization history. Such a late reionization history suggests that major sources of reionization would emerge late and be hosted by moderately massive halos in contrast with the widely-accepted picture of abundant low-mass objects for the sources of reionization., Comment: Accepted for publication in ApJ

We present a protocluster search covering $z\sim3$ to $z\sim5$ based on the combination of the Hyper SuprimeCam Subaru Strategic Programme and the CFHT Large Area $U$-band Deep Survey. We identify about 30 protocluster candidates per unit redshift over the $\sim25\,\mathrm{deg^2}$ area of the Deep/Ultra-Deep layer. Protocluster candidates are selected as regions with a significantly enhanced surface density of dropout galaxies. With this large sample, we characterise the properties of their individual member galaxies. We compare the number counts of dropout galaxies in protocluster candidates with that of coeval field galaxies. Rest-frame UV bright galaxies are over-abundant in protocluster candidates, a trend seen across the full redshift range studied. We do not find evidence for their spatial distribution within protocluster candidates to be distinct from their fainter counterparts, nor for their UV colour to be different from that of field galaxies with the same brightness. Cosmological simulations predict this bright-end excess, with the main cause being a richer population of massive galaxies, with only a minor contribution from an enhancement in star formation activity (and therefore UV emission) at fixed mass. $U$-to-$K$ SED modelling of our observed samples supports this interpretation. This environmental differentiation in number counts is already in place at $z\sim5$, with no significant redshift dependence over the range in lookback times probed. These observational results and model predictions suggest that the cosmic clock is ahead in high-density environments., Comment: 17 pages, 13 figures, accepted for publication in MNRAS

We study outflows in 130 galaxies with -221000 km/s. With the velocity shift and line widths of the outflow broad lines, we obtain ~80-500 km/s for the outflow velocities. We find that the outflow velocities as a function of star-formation rate are comparable to or higher than those of galaxies at z~1, accounting for the selection bias, while the outflow velocities of AGN are large but not significantly different from the others. Interestingly, these outflow velocities are typically not high enough to escape from the galactic potentials, suggestive of fountain-type outflows, which are concluded on the basis of thorough comparisons with recent JWST results. We estimate mass loading factors ${\eta}$ to be 0.1-1 that are not particularly large, but comparable with those of z~1 outflows. The large fraction of galaxies with outflows (30% with high resolution data) provides constraints on outflow parameters, suggesting a wide opening angle of >45 deg and a large duty-cycle of >30%, which gives a picture of more frequent and spherical outflows in high-z galaxies., Comment: 17 pages, 12 figures, re-submitted to ApJ on May 2024

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\mu$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$\alpha$ and [OIII], on optical morphologies with medium-band F410M images. Conducting S\'ersic 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\sigma$ 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., Comment: 33 pages, 18 figures, 6 tables, accepted

We present chemical abundance ratios of 70 star-forming galaxies at $z\sim4$-10 observed by the JWST/NIRSpec ERO, GLASS, and CEERS programs. Among the 70 galaxies, we have pinpointed 2 galaxies, CEERS_01019 at $z=8.68$ and GLASS_150008 at $z=6.23$, with extremely low C/N ([C/N]$\lesssim -1$), evidenced with CIII]$\lambda\lambda$1907,1909, NIII]$\lambda$1750, and NIV]$\lambda\lambda$1483,1486, which show high N/O ratios ([N/O]$\gtrsim 0.5$) comparable with the one of GN-z11 regardless of whether stellar or AGN radiation is assumed. Such low C/N and high N/O ratios found in CEERS_01019 and GLASS_150008 (additionally identified in GN-z11) are largely biased towards the equilibrium of the CNO cycle, suggesting that these 3 galaxies are enriched by metals processed by the CNO cycle. On the C/N vs. O/H plane, these 3 galaxies do not coincide with Galactic HII regions, normal star-forming galaxies, and nitrogen-loud quasars with asymptotic giant branch stars, but globular-cluster (GC) stars, indicating a connection with GC formation. We compare C/O and N/O of these 3 galaxies with those of theoretical models, and find that these 3 galaxies are explained by scenarios with dominant CNO-cycle materials, i.e. Wolf-Rayet stars, supermassive ($10^{3}-10^{5}\ M_{\odot}$) stars, and tidal disruption events, interestingly with a requirement of frequent direct collapses. For all the 70 galaxies, we present measurements of Ne/O, S/O, and Ar/O, together with C/O and N/O. We identify 4 galaxies with very low Ne/O, $\log(\rm Ne/O)<-1.0$, indicating abundant massive ($\gtrsim30\ M_\odot$) stars., Comment: Accepted for publication in ApJ

We present statistics of $z\sim 6-9$ galaxy outflows indicated by spatially-extended gas emission and broad lines. With a total of 61 spectroscopically confirmed galaxies at $z\sim 6-9$ in the JWST CEERS, GLASS, and ERO data, we find four galaxies with [O{\sc iii}]+H$\beta$ ionized gas emission significantly extended beyond the kpc-scale stellar components on the basis of the emission line images constructed by the subtraction of NIRCam broadband (line on/off-band) images. By comparison with low-$z$ galaxies, the fraction of galaxies with the spatially extended gas, 4/18, at $z\sim 6-9$ is an order of magnitude higher than those at $z\sim 0-1$, which can be explained by events triggered by frequent major mergers at high redshift. We also investigate medium- and high-resolution NIRSpec spectra of 30 galaxies at $z\sim 6-9$, and identify five galaxies with broad ($140-800$ km s$^{-1}$) lines in the [O{\sc iii}] forbidden line emission, suggestive of galaxy outflows. One galaxy at $z=6.38$ shows both the spatially-extended gas emission and the broad lines, while none of the galaxies with the spatially-extended gas emission or broad lines present a clear signature of AGN either in the line diagnostics or Type 1 AGN line broadening ($>1000$ km s$^{-1}$), which hint outflows mainly driven by stellar feedback. The existence of galaxies with/without spatially-extended gas emission or broad lines may be explained by different viewing angles towards outflows, or that these are galaxies in the early, late, post phases of galaxy outflows at high redshift, where the relatively large fractions of such galaxies indicate the longer-duration and/or more-frequent outflows at the early cosmic epoch., Comment: 16 pages, 11 figures, accepted for publication in ApJ

We present volume-averaged neutral hydrogen fractions $x_{\rm \HI}$ and ionized bubble radii $R_{\rm b}$ measured with Ly$\alpha$ damping wing absorption of galaxies at the epoch of reionization. We combine JWST/NIRSpec spectra taken by CEERS, GO-1433, DDT-2750, and JADES programs, and obtain a sample containing 27 bright UV-continuum ($M_{\rm UV}<-18.5~{\rm mag}$) galaxies at $7

We present new Spitzer Infrared Array Camera (IRAC) 3.6 and 4.5{\mu}m mosaics of three fields, E-COSMOS, DEEP2-F3, and ELAIS-N1. Our mosaics include both new IRAC observations as well as re-processed archival data in these fields. These fields are part of the HSC-Deep grizy survey and have a wealth of additional ancillary data. The addition of these new IRAC mosaics is critical in allowing for improved photometric redshifts and stellar population parameters at cosmic noon and earlier epochs. The total area mapped by this work is {\sim} 17 deg2 with a mean integration time of {\sim}1200s, providing a median 5{\sigma} depth of 23.7(23.3) at 3.6(4.5){\mu}m in AB. We perform SExtractor photometry both on the combined mosaics as well as the single-epoch mosaics taken {\sim}6 months apart. The resultant IRAC number counts show good agreement with previous studies. In combination with the wealth of existing and upcoming spectro-photometric data in these fields, our IRAC mosaics will enable a wide range of galactic evolution and AGN studies. With that goal in mind, we make the combined IRAC mosaics and coverage maps of these three fields publicly available. counts show good agreement with previous studies., Comment: accepted for publication in AJ

We present the quasar luminosity function (LF) at $z = 7$, measured with 35 spectroscopically confirmed quasars at $6.55 < z < 7.15$. The sample of 22 quasars from the Subaru High-$z$ Exploration of Low-Luminosity Quasars (SHELLQs) project, combined with 13 brighter quasars in the literature, covers an unprecedentedly wide range of rest-frame ultraviolet magnitudes over $-28 < M_{1450} < -23$. We found that the binned LF flattens significantly toward the faint end populated by the SHELLQs quasars. A maximum likelihood fit to a double power-law model has a break magnitude $M^*_{1450} = -25.60^{+0.40}_{-0.30}$, a characteristic density $\Phi^* = 1.35^{+0.47}_{-0.30}$ Gpc$^{-3}$ mag$^{-1}$, and a bright-end slope $\beta = -3.34^{+0.49}_{-0.57}$, when the faint-end slope is fixed to $\alpha = -1.2$ as observed at $z \le 6$. The overall LF shape remains remarkably similar from $z = 4$ to $7$, while the amplitude decreases substantially toward higher redshifts, with a clear indication of an accelerating decline at $z \ge 6$. The estimated ionizing photon density, $10^{48.2 \pm 0.1}$ s$^{-1}$ Mpc$^{-3}$, is less than 1 % of the critical rate to keep the intergalactic medium ionized at $z = 7$, and thus indicates that quasars are not a major contributor to cosmic reionization., Comment: The Astrophysical Journal Letters, in press

We present 20,567 Ly$\alpha$ emitters (LAEs) at $z=2.2-7.3$ that are photometrically identified by the SILVERRUSH program in a large survey area up to 25 deg$^2$ with deep images of five broadband filters (grizy) and seven narrowband filters targeting Ly$\alpha$ lines at $z=2.2$, $3.3$, $4.9$, $5.7$, $6.6$, $7.0$, and $7.3$ taken by the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) and the Cosmic HydrOgen Reionization Unveiled with Subaru (CHORUS) survey. We select secure $>5\sigma$ sources showing narrowband color excesses via Ly$\alpha$ break screening, taking into account the spatial inhomogeneity of limiting magnitudes. After removing spurious sources by careful masking and visual inspection of coadded and multi-epoch images obtained over the 7 yr of the surveys, we construct LAE samples consisting of 6995, 4641, 726, 6124, 2058, 18, and 5 LAEs at $z=2.2$, 3.3, 4.9, 5.7, 6.6, 7.0, and 7.3, respectively, although the $z=7.3$ candidates are tentative. Our LAE catalogs contain 241 spectroscopically confirmed LAEs at the expected redshifts from previous work. We demonstrate that the number counts of our LAEs are consistent with previous studies with similar LAE selection criteria. The LAE catalogs will be made public on our project webpage with detailed descriptions of the content and ancillary information about the masks and limiting magnitudes., Comment: 38 pages, 19 Figures, 5 Tables. Accepted for publication in ApJS

We present ionizing spectra estimated at 13.6--100 eV for ten dwarf galaxies with strong high ionization lines of He {\sc {ii}}$\lambda$4686 and [Ne {\sc{v}}]$\lambda$3426 ([Ne {\sc{iv}}]$\lambda$2424) at $z=0$ ($z=8$) that are identified in our Keck/LRIS spectroscopy and the literature (the JWST ERO program). With the flux ratios of these high ionization lines and $>10$ low-ionization lines of hydrogen, helium, oxygen, neon, and sulfur, we determine ionizing spectra consisting of stellar and non-thermal power-law radiation by photoionization modeling with free parameters of nebular properties including metallicity and ionization parameter, cancelling out abundance ratio differences. We find that all of the observed flux ratios are well reproduced by the photoinization models with the power law index $\alpha_{\rm EUV}$ of $\alpha_{\rm EUV}\sim (-1)-0$ and the luminosity $L_{\rm EUV}$ of $L_{\rm EUV}\sim 10^{40}-10^{42}$ erg s$^{-1}$ at $\sim 55-100$ eV for six galaxies, while four galaxies include large systematics in $\alpha_{\rm EUV}$ caused by stellar radiation contamination. We then compare $\alpha_{\rm EUV}$ and $L_{\rm EUV}$ of these six galaxies with those predicted by the black hole (BH) accretion disk models, and find that $\alpha_{\rm EUV}$ and $L_{\rm EUV}$ are similar to those of the intermediate mass black holes (IMBHs) in BH accretion disk models {albeit with possibilities of the other scenarios.} Confirming these results with a known IMBH having a mass $M_{\rm BH}$ of $M_{\rm BH}=10^{5.75} \ M_\odot$, we find that four local galaxies and one $z=7.665$ galaxy have ionizing spectra consistent with those of IMBHs with $M_{\rm BH} \sim 10^3-10^5 \ M_\odot$., Comment: Accepted for publication in ApJ and 25 pages, 12 figures

We present gas-phase elemental abundance ratios of 7 local extremely metal-poor galaxies (EMPGs) including our new Keck/LRIS spectroscopy determinations together with 33 JWST $z\sim 4-10$ star-forming galaxies in the literature, and compare chemical evolution models. We develop chemical evolution models with the yields of core-collapse supernovae (CCSNe), Type Ia supernovae, hypernovae (HNe), and pair-instability supernovae (PISNe), and compare the EMPGs and high-$z$ galaxies in conjunction with dust depletion contributions. We find that high Fe/O values of EMPGs can (cannot) be explained by PISN metal enrichments (CCSN/HN enrichments even with the mixing-and-fallback mechanism enhancing iron abundance), while that the observed Ar/O and S/O values are much smaller than the predictions of the PISN models. The abundance ratios of the EMPGs can be explained by the combination of Type Ia SNe and CCSNe/HNe whose inner layers of argon and sulfur mostly fallback, which are comparable with Sculptor stellar chemical abundance distribution, suggesting that early chemical enrichment is taken place in the EMPGs. Comparing our chemical evolution models with the star-forming galaxies at $z\sim 4-10$, we find that the Ar/O and S/O ratios of the high-$z$ galaxies are comparable with those of the CCSNe/HNe models, while majority of the high-$z$ galaxies do not have constraints good enough to rule out contributions from PISNe. The high N/O ratio recently reported in GN-z11 cannot be explained even by rotating PISNe, but could be reproduced by the winds of rotating Wolf Rayet stars that end up as a direct collapse.

We present pure spectroscopic constraints on the UV luminosity functions and cosmic star formation rate (SFR) densities from 25 galaxies at $z_\mathrm{spec}=8.61-13.20$. By reducing the JWST/NIRSpec spectra taken in multiple programs of ERO, ERS, GO, and DDT with our analysis technique, we independently confirm 16 galaxies at $z_\mathrm{spec}=8.61-11.40$ including new redshift determinations, and a bright interloper at $z_\mathrm{spec}=4.91$ that was claimed as a photometric candidate at z~16. In conjunction with nine galaxies at redshifts up to $z_\mathrm{spec}=13.20$ in the literature, we make a sample of 25 spectroscopically-confirmed galaxies in total and carefully derive the best estimates and lower limits of the UV luminosity functions. These UV luminosity function constraints are consistent with the previous photometric estimates within the uncertainties and indicate mild redshift evolution towards z~12 showing tensions with some theoretical models of rapid evolution. With these spectroscopic constraints, we obtain firm lower limits of the cosmic SFR densities and spectroscopically confirm a high SFR density at z~12 beyond the constant star-formation efficiency models, which supports earlier claims from the photometric studies. While there are no spectroscopically-confirmed galaxies with very large stellar masses violating the $\Lambda$CDM model due to the removal of the bright interloper, we confirm star-forming galaxies at $z_\mathrm{spec}=11-13$ with stellar masses much higher than model predictions. Our results indicate possibilities of high star-formation efficiency (>5%), hidden AGN, top-heavy initial mass function (possibly with Pop-III), and large scatter/variance. Having these successful and unsuccessful spectroscopy results, we suggest observational strategies for efficiently removing low redshift interlopers for future JWST programs., Comment: 27 pages, 19 figures, ApJ in press

Since the late 20th century, young star-forming dwarf galaxies with extremely low heavy-element abundances have been identified in the local universe\cite{1972ApJ...173...25S, 1990Natur.343..238I, 2022ApJS..262....3N}. Such a population of galaxies, represented by SBS 0335-052E, is intensively studied as a laboratory of galaxy formation, and interpreted as galaxies with hard emission and energetic outflows driven by massive stars\cite{1990Natur.343..238I, 2004ApJ...606..213T, 2009AA...503...61I}. Here we report the temporal flux variability of SBS 0335-052E in $3-4 \ \mu$m bands on timescale of months to years with dimming and brightening up to 50\% over 12 years. This is a clear signature of dust torus emission of an active massive black hole in SBS 0335-052E. The deep optical spectrum reveals a very broad component with $1.24\ (\pm 0.01) \times 10^{4}$ km s$^{-1}$ in H$\alpha$ emission, suggesting emission originated from broad line regions around the massive black hole. The black hole mass is estimated to be $\sim 100$ thousand solar masses from the spectral energy distribution. This is the first example of the massive black hole existence in a long-considered young star-forming dwarf galaxy., Comment: 32 pages, 6 figures, 3 tables

We present demography of the dynamics and gas-mass fraction of 33 extremely metal-poor galaxies (EMPGs) with metallicities of $0.015-0.195~Z_\odot$ and low stellar masses of $10^4-10^8~M_\odot$ in the local universe. We conduct deep optical integral-field spectroscopy (IFS) for the low-mass EMPGs with the medium high resolution ($R=7500$) grism of the 8m-Subaru FOCAS IFU instrument by the EMPRESS 3D survey, and investigate H$\alpha$ emission of the EMPGs. Exploiting the resolution high enough for the low-mass galaxies, we derive gas dynamics with the H$\alpha$ lines by the fitting of 3-dimensional disk models. We obtain an average maximum rotation velocity ($v_\mathrm{rot}$) of $15\pm3~\mathrm{km~s^{-1}}$ and an average intrinsic velocity dispersion ($\sigma_0$) of $27\pm10~\mathrm{km~s^{-1}}$ for 15 spatially resolved EMPGs out of the 33 EMPGs, and find that all of the 15 EMPGs have $v_\mathrm{rot}/\sigma_0<1$ suggesting dispersion dominated systems. There is a clear decreasing trend of $v_\mathrm{rot}/\sigma_0$ with the decreasing stellar mass and metallicity. We derive the gas mass fraction ($f_\mathrm{gas}$) for all of the 33 EMPGs, and find no clear dependence on stellar mass and metallicity. These $v_\mathrm{rot}/\sigma_0$ and $f_\mathrm{gas}$ trends should be compared with young high-$z$ galaxies observed by the forthcoming JWST IFS programs to understand the physical origins of the EMPGs in the local universe., Comment: 18 pages, 9 figures, accepted for publication in ApJ

We present a first statistical sample of faint type-1 AGNs at $z>4$ identified by JWST/NIRSpec deep spectroscopy. Among the 185 galaxies at $z_\mathrm{spec}=3.8-8.9$ confirmed with NIRSpec, our systematic search for broad-line emission reveals 10 type-1 AGNs at $z=4.015-6.936$ whose broad component is only seen in the permitted H$\alpha$ line and not in the forbidden [OIII]$\lambda$5007 line that is detected with greater significance than H$\alpha$. The broad H$\alpha$ line widths of $\mathrm{FWHM}\simeq1000-6000\ \mathrm{km\ s^{-1}}$ suggest that the AGNs have low-mass black holes with $M_\mathrm{BH}\sim10^6-10^8\ M_\odot$, remarkably lower than those of low-luminosity quasars previously identified at $z>4$ with ground-based telescopes. JWST and HST high-resolution images reveal that the majority of them show extended morphologies indicating significant contribution to the total lights from their host galaxies, except for three compact objects two of which show red SEDs, probably in a transition phase from faint AGNs to low luminosity quasars. Careful AGN-host decomposition analyses show that their host's stellar masses are systematically lower than the local relation between the black hole mass and the stellar mass, implying a fast black hole growth consistent with predictions from theoretical simulations. A high fraction of the broad-line AGNs ($\sim5\%$), higher than $z\sim0$, indicates that a number density of such faint AGNs is higher than an extrapolation of the quasar luminosity function, implying a large population of AGNs including type 1 and type 2 in the early universe. Such faint AGNs contribute to cosmic reionization, while the total contribution is not large, up to $\sim50\%$ at $z\sim6$, because of their faint nature., Comment: 24 pages, 16 figures, accepted for publication in ApJ

We investigate the stellar mass - black hole mass ($\mathcal{M}_*-\mathcal{M}_\mathrm{BH}$) relation with type 1 AGN down to $\mathcal{M}_\mathrm{BH}=10^7 M_\odot$, corresponding to a $\simeq -21$ absolute magnitude in rest-frame ultraviolet (UV), at $z = 2-2.5$. Exploiting the deep and large-area spectroscopic survey of the Hobby-Eberly Telescope Dark Energy Experiment (HETDEX), we identify 66 type 1 AGN with $\mathcal{M}_\mathrm{BH}$ ranging from $10^7$ to $10^{10} M_\odot$ that are measured with single-epoch virial method using C{\sc iv} emission lines detected in the HETDEX spectra. $\mathcal{M}_*$ of the host galaxies are estimated from optical to near-infrared photometric data taken with Spitzer, WISE, and ground-based 4-8m class telescopes by CIGALE SED fitting. We further assess the validity of SED fitting in two cases by host-nuclear decomposition performed through surface brightness profile fitting on spatially-resolved host galaxies with JWST/NIRCam CEERS data. We obtain the $\mathcal{M}_*-\mathcal{M}_\mathrm{BH}$ relation covering the unexplored low-mass ranges of $\mathcal{M}_\mathrm{BH}~\sim~10^7-10^8~M_\odot$, and conduct forward modelling to fully account for the selection biases and observational uncertainties. The intrinsic $\mathcal{M}_*-\mathcal{M}_\mathrm{BH}$ relation at $z\sim 2$ has a moderate positive offset of $0.52\pm0.14$~dex from the local relation, suggestive of more efficient black hole growth at higher redshift even in the low-mass regime of $\mathcal{M}_\mathrm{BH}~\sim~10^7-10^8~M_\odot$. Our $\mathcal{M}_*-\mathcal{M}_\mathrm{BH}$ relation is inconsistent with the $\mathcal{M}_\mathrm{BH}$ suppression at the low-$\mathcal{M}_*$ regime predicted by recent hydrodynamic simulations at a $98\%$ confidence level, suggesting that feedback in the low-mass systems may be weaker than those produced in hydrodynamic simulations., Comment: 16 pages, 8 figures, accepted for publication in ApJ

We search for local extremely metal-poor galaxies (EMPGs), selecting photometric candidates by broadband color excess and machine-learning techniques with the SDSS photometric data. After removing stellar contaminants by shallow spectroscopy with Seimei and Nayuta telescopes, we confirm that three candidates are EMPGs with 0.05--0.1 $Z_\odot$ by deep Magellan/MagE spectroscopy for faint {\sc[Oiii]}$\lambda$4363 lines. Using a statistical sample consisting of 105 spectroscopically-confirmed EMPGs taken from our study and the literature, we calculate cross-correlation function (CCF) of the EMPGs and all SDSS galaxies to quantify environments of EMPGs. Comparing another CCF of all SDSS galaxies and comparison SDSS galaxies in the same stellar mass range ($10^{7.0}-10^{8.4} M_\odot$), we find no significant ($>1\sigma$) difference between these two CCFs. We also compare mass-metallicity relations (MZRs) of the EMPGs and those of galaxies at $z\sim$ 0--4 with a steady chemical evolution model and find that the EMPG MZR is comparable with the model prediction on average. These clustering and chemical properties of EMPGs are explained by a scenario of stochastic metal-poor gas accretion on metal-rich galaxies showing metal-poor star formation. Extending the broadband color-excess technique to a high-$z$ EMPG search, we select 17 candidates of $z\sim$ 4--5 EMPGs with the deep ($\simeq30$ mag) near-infrared JWST/NIRCam images obtained by ERO and ERS programs. We find galaxy candidates with negligible {\sc[Oiii]}$\lambda\lambda$4959,5007 emission weaker than the local EMPGs and known high-$z$ galaxies, suggesting that some of these candidates may fall in 0--0.01 $Z_\odot$, which potentially break the lowest metallicity limit known to date.

We present the evolution of the mass-metallicity (MZ) relations at z=4-10 derived with 135 galaxies identified in the JWST/NIRSpec data taken from the three major public spectroscopy programs of ERO, GLASS, and CEERS. Because there are many discrepancies between flux measurements reported by early ERO studies, we first establish our NIRSpec data reduction procedure for reliable emission-line flux measurements and errors successfully explaining Balmer decrements with no statistical tensions via thorough comparisons of the early ERO studies. Applying the reduction procedure to the 135 galaxies, we obtain emission-line fluxes for physical property measurements. We confirm that 10 out of the 135 galaxies with [OIII]4363-lines have electron temperatures of (1.1-2.3)*10^4K, similar to lower-z star-forming galaxies, that can be explained by heating of young massive stars. We derive metallicities of the 10 galaxies by the direct method and the rest of the galaxies with strong lines by the metallicity calibrations of Nakajima et al. (2022) applicable for these low-mass metal-poor galaxies, anchoring the metallicities with the direct-method measurements. We thus obtain MZ relations and star-formation rate (SFR)-MZ relations over z=4-10. We find that there is a small evolution of the MZ relation from z~2-3 to z=4-10, while interestingly that the SFR-MZ relation shows no evolution up to z~8 but a significant decrease at z>8 beyond the error. This SFR-MZ relation decrease at z>8 may suggest a break of the metallicity equilibrium state via star-formation, inflow, and outflow, while further statistical and local-baseline studies are needed for a conclusion., Comment: 39 pages, 16 figures in the main text and 6 figures in Appendices. Accepted for publication in ApJ Supplement Series. Visit https://sites.google.com/view/kimihiko-nakajima/research for the summary table of the JWST objects (Table D1)

We present electron densities $n_{\rm e}$ in the inter-stellar medium (ISM) of star-forming galaxies at $z=4-9$ observed by the JWST/NIRSpec GLASS, ERO, and CEERS programs. We carefully evaluate line-spread functions of the NIRSpec instrument as a function of wavelength with the calibration data of a planetary nebula taken onboard, and obtain secure [OII]$\lambda\lambda$3726,3729 doublet fluxes for 14 galaxies at $z=4.02-8.68$ falling on the star-formation main sequence with the NIRSpec high and medium resolution spectra. We thus derive the electron densities of singly-ionized oxygen nebulae with the standard $n_{\rm e}$ indicator of [OII] doublet, and find that the electron densities of the $z=4-9$ galaxies are $n_{\rm e}\gtrsim 300$ cm$^{-3}$ significantly higher than those of low-$z$ galaxies at a given stellar mass, star-formation rate (SFR), and specific SFR. Interestingly, typical electron densities of singly ionized nebulae increase from $z=0$ to $z=1-3$ and $z=4-9$, which is approximated by the evolutionary relation of $n_{\rm e}\propto(1+z)^{p}$ with $p\sim 1-2$. Although it is not obvious that the ISM property of $n_{\rm e}$ is influenced by global galaxy properties, these results may suggest that nebula densities of high-$z$ galaxies are generally high due to the compact morphologies of high-$z$ galaxies evolving by $r_{\rm e}$ approximately proportional to $(1+z)^{-1}$ ($r_{\rm vir} \propto (1+z)^{-1}$) for a given stellar (halo) mass whose inverse square corresponds to the $p\sim 2$ evolutionary relation. The $p\sim 1-2$ evolutionary relation can be explained by a combination of the compact morphology and the reduction of $n_{\rm e}$ due to the high electron temperature of the high-$z$ metal poor nebulae., Comment: Accepted for publication in ApJ

We present cosmological-scale 3-dimensional (3D) neutral hydrogen ({\sc Hi}) tomographic maps at $z=2-3$ over a total of 837 deg$^2$ in two blank fields that are developed with Ly$\alpha$ forest absorptions of 14,736 background Sloan Digital Sky Survey (SDSS) quasars at $z$=2.08-3.67. Using the tomographic maps, we investigate the large-scale ($\gtrsim 10$ $h^{-1}$cMpc) average {\sc Hi} radial profiles and two-direction profiles of the line-of-sight (LoS) and transverse (Trans) directions around galaxies and AGN at $z=2-3$ identified by the Hobby-Eberly Telescope Dark Energy eXperiment (HETDEX) and SDSS surveys, respectively. The peak of the {\sc Hi} radial profile around galaxies is lower than the one around AGN, suggesting that the dark-matter halos of galaxies are less massive on average than those of AGN. The LoS profile of AGN is narrower than the Trans profile, indicating the Kaiser effect. There exist weak absorption outskirts at $\gtrsim 30$ $h^{-1}$cMpc beyond {\sc Hi} structures of galaxies and AGN found in the LoS profiles that can be explained by the {\sc Hi} gas at $\gtrsim 30$ $h^{-1}$cMpc falls toward the source positions. Our findings indicate that the {\sc Hi} radial profile of AGN has transitions from proximity zones ($\lesssim$ a few $h^{-1}$cMpc) to the {\sc Hi} structures ($\sim 1-30$ $h^{-1}$cMpc) and the weak absorption outskirts ($\gtrsim 30$ $h^{-1}$cMpc). Although there is no significant dependence of AGN types (type-1 vs. type-2) on the {\sc Hi} profiles, the peaks of the radial profiles anti-correlate with AGN luminosities, suggesting that AGN's ionization effects are stronger than the gas mass differences., Comment: 26 pages, 32 figures, accepted to ApJ

We present a novel technique called "photometric IGM tomography" to map the intergalactic medium (IGM) at $z\simeq4.9$ in the COSMOS field. It utilizes deep narrow-band (NB) imaging to photometrically detect faint Ly$\alpha$ forest transmission in background galaxies across the Subaru/Hyper-Suprime Cam (HSC)'s $1.8\rm\,sq.\,deg$ field of view and locate Ly$\alpha$ emitters (LAEs) in the same cosmic volume. Using ultra-deep HSC images and Bayesian spectral energy distribution fitting, we measure the Ly$\alpha$ forest transmission at $z\simeq4.9$ along a large number ($140$) of background galaxies selected from the DEIMOS10k spectroscopic catalogue at $4.98

We study the formation of the first galaxies in overdense regions modelled by the FORmation and EVolution of galaxies in Extremely overdense Regions motivated by SSA22 (FOREVER22) simulation project. Our simulations successfully reproduce the star formation rates and the $M_{\rm UV}-M_{\rm star}$ relations of candidate galaxies at $z \sim 10-14$ observed by the James Webb Space Telescope (JWST). We suggest that the observed galaxies are hosted by dark-matter haloes with $M_{\rm h} \gtrsim 10^{10}~{\rm M_{\odot}}$ and are in short-period starburst phases. On the other hand, even simulated massive galaxies in overdense regions cannot reproduce the intense star formation rates and the large stellar masses of observed candidates at $z \sim 16$. Also, we show that the contribution of population III stars to the UV flux decreases as the stellar mass increases and it is a few percent for galaxies with $M_{\rm star} \sim 10^{7}~{\rm M_{\odot}}$. Therefore, a part of the observed flux by JWST could be the light from population III stars. Our simulations suggest that the UV flux can be dominated by population III stars and the UV-slope shows $\beta \lesssim -3$ if future observations would reach galaxies with $M_{\rm stars} \sim 10^{5}~{\rm M_{\odot}}$ at $z \sim 20$ of which the mass fraction of population III stars can be greater than 10 percent., Comment: 9 pages, 6 figures, accepted for publication in MNRAS

We present our new Atacama Large Millimeter/submillimeter Array (ALMA) observations targeting CO(6-5) emission from three luminous Lyman break galaxies (LBGs) at $z_{\rm spec} = 6.0293$-$6.2037$ found in the Subaru/Hyper Suprime-Cam survey, whose [OIII]$88\mu$m and [CII]$158\mu$m emission have been detected with ALMA. We find a marginal detection of the CO(6-5) line from one of our LBGs, J0235-0532, at the $\simeq 4 \sigma$ significance level and obtain upper limits for the other two LBGs, J1211-0118 and J0217-0208. Our $z=6$ luminous LBGs are consistent with the previously found correlation between the CO luminosity and the infrared luminosity. The unique ensemble of the multiple far-infrared emission lines and underlying continuum fed to a photodissociation region model reveal that J0235-0532 has a relatively high hydrogen nucleus density that is comparable to those of low-$z$ (U)LIRGs, quasars, and Galactic star-forming regions with high $n_{\rm H}$ values, while the other two LBGs have lower $n_{\rm H}$ consistent with local star-forming galaxies. By carefully taking account of various uncertainties, we obtain total gas mass and gas surface density constraints from their CO luminosity measurements. We find that J0235-0532 locates below the Kennicutt-Schmidt (KS) relation, comparable to the previously CO(2-1) detected $z=5.7$ LBG, HZ10. Combined with previous results for dusty starbursts at similar redshifts, the KS relation at $z=5$-$6$ is on average consistent with the local one., Comment: 33 pages, 12 figures, accepted for publication in ApJ

We report on the galaxy density environment around a high-z radio galaxy (HzRG) at z=4.72, HSC J083913.17+011308.1 (HSC J0839+0113), probed using an r-dropout Lyman break galaxy (LBG) sample from the Hyper Suprime-Cam Subaru Strategic Program data. We find that HSC J0839+0113 resides in the outskirt of an overdense region identified by the r-dropout galaxies at a 4.7 sigma significance level. The projected distance between HSC J0839+0113 and the peak position of the overdense region is 0.4 physical Mpc which is shorter than the typical protocluster radius in this epoch. According to the extended Press Schechter and the light cone models, the HSC J0839+0113-hosted overdense region is expected to evolve into a halo > 10^14 Msun at z=0 with a high probability of >80 %. These findings suggest that HSC J0839+0113 is associated with a protocluster. The HSC J0839+0113 rich-system is the most overdense region of LBGs among the known protoclusters with LBGs in the same cosmic epoch., Comment: 7 pages, 3 figures, accepted in PASJ

We present morphologies of galaxies at $z \gtrsim 9$ resolved by JWST/NIRCam $2$-$5\mu$m imaging. Our sample consists of $22$ galaxy candidates identified by stringent dropout and photo-$z$ criteria in GLASS, CEERS, SMACS J0723, and Stephan's Quintet flanking fields, one of which has been spectroscopically identified at $z=11.44$. We perform surface brightness (SB) profile fitting with GALFIT for $6$ bright galaxies with S/N $=10$-$40$ on an individual basis and for stacked faint galaxies with secure point-spread functions (PSFs) of the NIRCam real data, carefully evaluating systematics by Monte-Carlo simulations. We compare our results with those of previous JWST studies, and confirm that effective radii $r_{\rm e}$ of our measurements are consistent with those of previous measurements at $z\sim 9$. We obtain $r_{\rm e}\simeq 200$-$300$ pc with the exponential-like profiles, S\'ersic indexes of $n\simeq 1$-$1.5$, for galaxies at $z\sim 12$-$16$, indicating that the relation of $r_{\rm e}\propto (1+z)^s$ for $s=-1.22^{+0.17}_{-0.16}$ explains cosmic evolution over $z\sim 0$-$16$ for $\sim L^*_{z=3}$ galaxies. One bright ($M_{\rm UV}=-21$ mag) galaxy at $z\sim 12$, GL-z12-1, has an extremely compact profile with $r_{\rm e}=39 \pm 11$ pc that is surely extended over the PSF. Even in the case that the GL-z12-1 SB is fit by AGN$+$galaxy composite profiles, the best-fit galaxy component is again compact, $r_{\rm e}=48^{+38}_{-15}$ pc that is significantly ($>5\sigma$) smaller than the typical $r_{\rm e}$ value at $z\sim 12$. Comparing with numerical simulations, we find that such a compact galaxy naturally forms at $z\gtrsim 10$, and that frequent mergers at the early epoch produce more extended galaxies following the $r_{\rm e}\propto (1+z)^s$ relation., Comment: 29 pages, 25 figures, accepted for publication in ApJ

We conduct a comprehensive study on dropout galaxy candidates at $z\sim9-16$ using the first 90 arcmin$^2$ JWST/NIRCam images taken by the early release observations (ERO) and early release science (ERS) programs. With the JWST simulation images, we find that a number of foreground interlopers are selected with a weak photo-$z$ determination ($\Delta\chi^2>4$). We thus carefully apply a secure photo-$z$ selection criterion ($\Delta\chi^2>9$) and conventional color criteria with confirmations of the ERO NIRSpec spectroscopic redshifts, and obtain a total of 23 dropout galaxies at $z\sim9-16$, including two candidates at $z_\mathrm{phot}=16.25_{-0.46}^{+0.24}$ and $16.41_{-0.55}^{+0.66}$. We perform thorough comparisons of dropout galaxies found in our work with recent JWST studies, and conclude that our galaxy sample is reliable enough for statistical analyses. We derive the UV luminosity functions at $z\sim9-16$, and confirm that our UV luminosity functions at $z\sim 9$ and $12$ agree with those determined by other HST and JWST studies. The cosmic star-formation rate density decreases from $z\sim9$ to $12$, and perhaps to $16$, but the densities at $z\sim12-16$ are higher than the constant star formation efficiency model. Interestingly, there are six bright galaxy candidates at $z\sim10-16$ with $M_\mathrm{UV}<-19.5$ mag and $M_*\sim10^{8-9} M_\odot$. Because a majority ($\sim80$\%) of these galaxies show no signatures of AGNs in their morphologies, the high cosmic star-formation rate densities and the existence of these UV-luminous galaxies are explained by no suppression of star-formation by the UV background radiation at the pre-reionization epoch and/or an efficient UV radiation production by a top-heavy IMF with Population III-like star formation., Comment: 38 pages, 21 figures. ApJS in press. All of the results are updated based on the latest NIRCam flux calibration (jwst_0995.pmap)

We present kinematics of 6 local extremely metal-poor galaxies (EMPGs) with low metallicities ($0.016-0.098\ Z_{\odot}$) and low stellar masses ($10^{4.7}-10^{7.6} M_{\odot}$). Taking deep medium-high resolution ($R\sim7500$) integral-field spectra with 8.2-m Subaru, we resolve the small inner velocity gradients and dispersions of the EMPGs with H$\alpha$ emission. Carefully masking out sub-structures originated by inflow and/or outflow, we fit 3-dimensional disk models to the observed H$\alpha$ flux, velocity, and velocity-dispersion maps. All the EMPGs show rotational velocities ($v_{\rm rot}$) of 5--23 km s$^{-1}$ smaller than the velocity dispersions ($\sigma_{0}$) of 17--31 km s$^{-1}$, indicating dispersion-dominated ($v_{\rm rot}/\sigma_{0}=0.29-0.80<1$) systems affected by inflow and/or outflow. Except for two EMPGs with large uncertainties, we find that the EMPGs have very large gas-mass fractions of $f_{\rm gas}\simeq 0.9-1.0$. Comparing our results with other H$\alpha$ kinematics studies, we find that $v_{\rm rot}/\sigma_{0}$ decreases and $f_{\rm gas}$ increases with decreasing metallicity, decreasing stellar mass, and increasing specific star-formation rate. We also find that simulated high-$z$ ($z\sim 7$) forming galaxies have gas fractions and dynamics similar to the observed EMPGs. Our EMPG observations and the simulations suggest that primordial galaxies are gas-rich dispersion-dominated systems, which would be identified by the forthcoming James Webb Space Telescope (JWST) observations at $z\sim 7$., Comment: ApJ in Press

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., Comment: 33 pages, 14 figures. Accepted for publication in ApJ Supplement Series. Visit https://sci.nao.ac.jp/MEMBER/ouchi/members/nakajima.html for the publicly available table of EMPGs (Table 6)

The primordial He abundance $Y_\mathrm{P}$ is a powerful probe of cosmology. Currently, $Y_\mathrm{P}$ is best determined by observations of metal-poor galaxies, while there are only a few known local extremely metal-poor ($<0.1 Z_\odot$) galaxies (EMPGs) having reliable He/H measurements with HeI$\lambda$10830 near-infrared (NIR) emission. Here we present deep Subaru NIR spectroscopy for 10 EMPGs. Combining the existing optical data, He/H values of 5 out of the 10 EMPGs are reliably derived by the Markov chain Monte Carlo algorithm. Adding the existing 3 EMPGs and 51 moderately metal-poor ($0.1-0.4 Z_\odot$) galaxies with reliable He/H estimates, we obtain $Y_\mathrm{P}=0.2370^{+0.0034}_{-0.0033}$ by linear regression in the $\mathrm{(He/H)}-\mathrm{(O/H)}$ plane, where we increase the number of EMPGs from 3 to 8 anchoring He/H of the most metal-poor gas in galaxies. Although our $Y_\mathrm{P}$ measurement and previous measurements are consistent, our result is slightly ($\sim 1\sigma$) smaller due to our EMPGs. With our $Y_\mathrm{P}$ and the existing primordial deuterium $D_\mathrm{P}$ measurement, we constrain the effective number of neutrino species $N_\mathrm{eff}$ and the baryon-to-photon ratio $\eta$ showing $\gtrsim 1-2\sigma$ tensions with the Standard Model and Planck Collaboration et al. (2020). Motivated by the tensions, we allow the degeneracy parameter of electron-neutrino $\xi_e$ to vary as well as $N_\mathrm{eff}$ and $\eta$. We obtain $\xi_e = 0.05^{+0.03}_{-0.02}$, $N_\mathrm{eff}=3.11^{+0.34}_{-0.31}$, and $\eta\times10^{10}=6.08^{+0.06}_{-0.06}$ from the $Y_\mathrm{P}$ and $D_\mathrm{P}$ measurements with a prior of $\eta$ taken from Planck Collaboration et al. (2020). Our constraints suggest a lepton asymmetry and allow for a high value of $N_\mathrm{eff}$ within the $1\sigma$ level, which could mitigate the Hubble tension., Comment: 22 pages, 9 figures, accepted for publication in ApJ

We report ALMA Band 3 observations of CO(6-5), CO(7-6), and [CI](2-1) in B14-65666 (``Big Three Dragons''), one of the brightest Lyman-Break Galaxies at $z>7$ in the rest-frame ultraviolet continuum, far-infrared continuum, and emission lines of [OIII] 88 $\mu$m and [CII] 158 $\mu$m. CO(6-5), CO(7-6), and [CI](2-1), whose $3\sigma$ upper limits on the luminosities are approximately 40 times fainter than the [CII] luminosity, are all not detected. The $L_{\rm [CII]}$/$L_{\rm CO(6-5)}$ and $L_{\rm [CII]}$/$L_{\rm CO(7-6)}$ ratios are higher than the typical ratios obtained in dusty star-forming galaxies or quasar host galaxies at similar redshifts, and they may suggest a lower gas density in the photodissociated region in B14-65666. By using the (1) [CII] luminosity, (2) dust mass-to-gas mass ratio, and (3) a dynamical mass estimate, we find that the molecular gas mass ($M_{\rm{mol}}$) is $(0.05-11)\times10^{10}$ $M_{\rm \odot}$. This value is consistent with the upper limit inferred from the nondetection of mid-$J$ CO and [CI](2-1). Despite the large uncertauinty in $M_{\rm mol}$, we estimate a molecular gas-to-stellar mass ratio ($\mu_{\rm{gas}}$) of $0.65-140$ and a gas depletion time ($\tau_{\rm dep}$) of $2.5-550$ Myr; these values are broadly consistent with those of other high-redshift galaxies. B14-65666 could be an ancestor of a passive galaxy at $z\gtrsim4$ if no gas is fueled from outside the galaxy., Comment: 15 pages, 6 figures, 4 tables. Accepted for publication in ApJ

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$\lambda$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., Comment: Accepted for publication by ApJ

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$\alpha$ and [OIII]$\lambda 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 $\eta$ can be as small as 0.1 - 1 or below, in contrast to the large $\eta$ of energy-driven outflows predicted by numerical simulations., Comment: 22 pages, 11 figures, Accepted for publication by ApJ

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$., Comment: 22 pages, 8 figures, accepted in ApJ

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., Comment: Accepted for publication in The Astrophysical Journal Supplement Series

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 describe the survey design, calibration, commissioning, and emission-line detection algorithms for the Hobby-Eberly Telescope Dark Energy Experiment (HETDEX). The goal of HETDEX is to measure the redshifts of over a million Ly$\alpha$ emitting galaxies between 1.88

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\sigma$) 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/., Comment: 25 pages, 19 figures, submitted to PASJ. Data available at https://hsc-release.mtk.nao.ac.jp/