Source author record

Akio K. Inoue

Akio K. Inoue appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

47works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

47 published item(s)

preprint2026arXiv

A z ~ 0.4 Galaxy Reflecting the high-redshift Little Red Dots: An Extended Starburst with an Overmassive Black Hole

One of the most remarkable discoveries of JWST is a population of compact, red sources at z > 4, commonly referred to as Little Red Dots (LRDs). Spectroscopic identifications reported that most LRDs are active galactic nuclei (AGNs), which are preferentially found around z ~ 6 and could imply a key phase in the formation and growth of black holes (BHs) in the early universe. Photometric surveys at lower redshift have recently been carried out to trace their evolution across cosmic time, and a small number of LRD-like galaxies have been spectroscopically identified at both Cosmic Noon and in the local universe. Here we report the discovery of one of the lowest-z LRD-like galaxies, J204837.26-002437.2 (hereafter J2048) at z = 0.4330, using new Gemini-N/GMOS IFU observations combined with archival multi-band photometric SED data. The GMOS data reveal extended blue emission from starburst with a star formation rate of 400 solar mass per year, together with an extended, highly fast ionized outflow. This is the first spectroscopic confirmation of extended host emission and outflow in an LRD-like galaxy, providing a unique laboratory for understanding the nature of their high-redshift counterparts. Moreover, J2048 would host an extremely overmassive BH with a BH-to-stellar mass ratio of ~ 60%, with the BH mass and host stellar mass estimated to be 10^10.2 and 10^10.4 solar masses, respectively. We discuss the origin and evolutionary fate of J2048, and the implications that such low-z analogs have for interpreting the properties of high-z LRDs.

preprint2022arXiv

A Candidate of a Least-Massive Black Hole at the First 1.1 Billion Years of the Universe

We report a candidate of a low-luminosity active galactic nucleus (AGN) at z=5 that was selected from the first near-infrared images of the JWST CEERS project. This source, named CEERS-AGN-z5-1 at absolute 1450 A magnitude M1450=-19.5 +/- 0.3, was found via a visual selection of compact sources from a catalog of Lyman break galaxies at z>4, taking advantage of the superb spatial resolution of the JWST/NIRCam images. The 20 photometric data available from CFHT, HST, Spitzer, and JWST suggest that the continuum shape of this source is reminiscent of that for an unobscured AGN, and there is a clear color excess in the filters where the redshifted Hbeta+[OIII] and Halpha are covered. The estimated line luminosity is L_Hbeta+[OIII] =10^43.0 erg s-1 and L_Halpha =10^42.9 erg s-1 with the corresponding rest-frame equivalent width EW_{Hbeta+[OIII]} =1100 A and EW_Halpha =1600 A, respectively. Our SED fitting analysis favors the scenario that this object is either a strong broad-line emitter or even a super-Eddington accreting black hole (BH), although a possibility of an extremely young galaxy with moderate dust attenuation is not completely ruled out. The bolometric luminosity, L_bol=2.5 +/- 0.3 \times 10^44 erg s-1, is consistent with those of z<0.35 broad-line AGNs with M_BH = 10^6 M_sun accreting at the Eddington limit. This new AGN population at the first 1.1 billion years of the universe may close the gap between the observed BH mass range at high redshift and that of BH seeds. Spectroscopic confirmation is awaited to secure the redshift and its AGN nature.

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

FOREVER22: galaxy formation in protocluster regions

We present results from a new cosmological hydrodynamics simulation campaign of protocluster (PC) regions, FOREVER22: FORmation and EVolution of galaxies in Extremely-overdense Regions motivated by SSA22. The simulations cover a wide range of cosmological scales using three different zoom set-ups in a parent volume of $(714.2~\rm cMpc)^{3}$: PCR (Proto-Cluster Region; $V= (28.6~{\rm cMpc})^{3} $, SPH particle mass, $m_{\rm{SPH}} = 4.1 \times 10^{6}~\rm M_{\odot}$ and final redshift, $z_{\rm end}=2.0$), BCG (Brightest proto-Cluster Galaxy; $V \sim (10~{\rm cMpc})^{3} $, $m_{\rm SPH} = 5.0\times10^{5}~\rm M_{\odot}$ and $z_{\rm end}=4.0$ ), and First ( $V \sim (3~{\rm cMpc})^{3} $, $m_{\rm SPH} = 7.9 \times 10^{3}~\rm M_{\odot}$ and $z_{\rm end}=9.5$) runs, that allow to focus on different aspects of galaxy formation. In the PCR runs, we follow 10 PCs, each harbouring 1 - 4 SMBHs with $M_{\rm BH} \ge 10^{9}~\rm M_{\odot}$. One of the PC cores shows a spatially close arrangement of seven starburst galaxies with ${\rm SFR} \gtrsim 100~\rm M_{\odot}~yr^{-1}$ each, that are dust-obscured and would appear as submillimeter galaxies with flux $\gtrsim 1~$ mJy at $1.1~ \rm mm$ in observations. The BCG runs show that the total SFRs of haloes hosting BCGs are affected by AGN feedback, but exceed $1000~\rm M_{\odot}~yr^{-1}$ at $z \lesssim 6$. The First runs resolve mini-haloes hosting population (Pop) III stars and we show that, in PC regions, the dominant stellar population changes from Pop III to Pop II at $z \gtrsim 20$, and the first galaxies with ${\rm SFR} \gtrsim 18~\rm M_{\odot}~yr^{-1}$ form at $z \sim 10$. These can be prime targets for future observations with the James Webb Space Telescope. Our simulations successfully reproduce the global star formation activities in observed PCs and suggest that PCs can kickstart cosmic reionization.

preprint2022arXiv

Possible Systematic Rotation in the Mature Stellar Population of a $z=9.1$ Galaxy

We present new observations with the Atacama Large Millimeter/submillimeter Array for a gravitationally-lensed galaxy at $z=9.1$, MACS1149-JD1. [O III] 88-$μ$m emission is detected at 10$σ$ with a spatial resolution of $\sim0.3$ kpc in the source plane, enabling the most distant morpho-kinematic study of a galaxy. The [O III] emission is distributed smoothly without any resolved clumps and shows a clear velocity gradient with $ΔV_{\rm obs}/2σ_{\rm tot}=0.84\pm0.23$, where $ΔV_{\rm obs}$ is the observed maximum velocity difference and $σ_{\rm tot}$ is the velocity dispersion measured in the spatially-integrated line profile, suggesting a rotating system. Assuming a geometrically thin self-gravitating rotation disk model, we obtain $V_{\rm rot}/σ_V=0.67_{-0.26}^{+0.73}$, where $V_{\rm rot}$ and $σ_V$ are the rotation velocity and velocity dispersion, respectively, still consistent with rotation. The resulting disk mass of $0.65_{-0.40}^{+1.37}\times10^{9}$ M$_\odot$ is consistent with being associated with the stellar mass identified with a 300 Myr-old stellar population independently indicated by a Balmer break in the spectral energy distribution. We conclude that the most of the dynamical mass is associated with the previously-identified mature stellar population that formed at $z\sim15$.

preprint2022arXiv

The Age of Discovery with the James Webb: Excavating the Spectral Signatures of the First Massive Black Holes

The James Webb Space Telescope (JWST) will open a new window of the most distant universe and unveil the early growth of supermassive black holes (BHs) in the first galaxies. In preparation for deep JWST imaging surveys, it is crucial to understand the color selection of high-redshift accreting seed BHs. We model the spectral energy distribution of super-Eddington accreting BHs with millions of solar masses in metal-poor galaxies at $z\gtrsim 8$, applying post-process line transfer calculations to radiation hydrodynamical simulation results. Ten kilosecond exposures with the NIRCam and MIRI broad-band filters are sufficient to detect the radiation flux from the seed BHs with bolometric luminosities of $L_{\rm bol}\simeq 10^{45}~{\rm erg~s}^{-1}$. While the continuum colors are similar to those of typical low-$z$ quasars, strong H$α$ line emission with a rest-frame equivalent width ${\rm EW}_{\rm rest}\simeq 1300~Å$ is so prominent that the line flux affects the broad-band colors significantly. The unique colors, for instance F356W$-$F560W $\gtrsim 1$ at $7<z<8$ and F444W$-$F770W $\gtrsim 1$ at $9<z<12$, provide robust criteria for photometric selection of the rapidly growing seed BHs. Moreover, NIRSpec observations of low-ionization emission lines can test whether the BH is fed via a dense accretion disk at super-Eddington rates.

preprint2022arXiv

Where's Swimmy?: Mining unique color features buried in galaxies by deep anomaly detection using Subaru Hyper Suprime-Cam data

We present the Swimmy (Subaru WIde-field Machine-learning anoMalY) survey program, a deep-learning-based search for unique sources using multicolored ($grizy$) imaging data from the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP). This program aims to detect unexpected, novel, and rare populations and phenomena, by utilizing the deep imaging data acquired from the wide-field coverage of the HSC-SSP. This article, as the first paper in the Swimmy series, describes an anomaly detection technique to select unique populations as "outliers" from the data-set. The model was tested with known extreme emission-line galaxies (XELGs) and quasars, which consequently confirmed that the proposed method successfully selected 60-70% of the quasars and 60% of the XELGs without labeled training data. In reference to the spectral information of local galaxies at $z=$0.05-0.2 obtained from the Sloan Digital Sky Survey, we investigated the physical properties of the selected anomalies and compared them based on the significance of their outlier values. The results revealed that XELGs constitute notable fractions of the most anomalous galaxies, and certain galaxies manifest unique morphological features. In summary, a deep anomaly detection is an effective tool that can search rare objects, and ultimately, unknown unknowns with large data-sets. Further development of the proposed model and selection process can promote the practical applications required to achieve specific scientific goals.

preprint2021arXiv

A Morphological Study on Galaxies Hosting Optical Variability-Selected AGNs in the COSMOS Field

The morphological study is crucial to investigate the connections between active galactic nuclei (AGN) activities and the evolution of galaxies. Substantial studies have found that radiative-mode AGNs primarily reside in disk galaxies, questioning the merger-driven mechanism of AGN activities. In this study, through S{é}rsic profile fitting and non-parametric morphological parameter measurements, we investigated the morphology of host galaxies of 485 optical variability-selected low luminosity AGNs at $z\lesssim4.26$ in the COSMOS field. We analyzed high-resolution images of the Hubble Space Telescope to measure these morphological parameters. We only successfully measured the morphological parameters for 76 objects and most AGN hosts ($\sim70\%$) were visually compact point-like sources. We examined the obtained morphological information as a function of redshift and compared them with literature data. We found that these AGN host galaxies showed no clear morphological preference. However, the merger rate increased with the higher hosts' SFRs and AGN luminosity. Interestingly, we found ongoing star formation consistent with the typical star forming populations in both elliptical and spiral galaxies while these two types of galaxies were more symmetric than normal star forming galaxies. These results suggested that optical variability-selected AGNs have} higher probabilities to reside in elliptical galaxies than infrared-selected AGNs (IR-AGNs), whose host galaxies had a strong disk-dominance, and supported recent studies that the AGN feedback could enhance star forming activities in host galaxies.

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

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. 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''.8 \times 0''.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.

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

A 16 deg$^2$ survey of emission-line galaxies at z<1.6 from HSC-SSP PDR2 and CHORUS

We have conducted a comprehensive survey of emission-line galaxies at $z\lesssim1.6$ based on narrowband (NB) imaging data taken with Hyper Suprime-Cam (HSC) on the Subaru telescope. In this paper, we update the catalogs of H$α$, [OIII], and [OII] emission-line galaxies using the data from the second Public Data Release (PDR2) of Subaru Strategic Program (SSP) of HSC and Cosmic HydrOgen Reionization Unveiled with Subaru (CHORUS) survey along with the spectroscopic redshifts for 2,019 emission-line galaxies selected with the PDR1 data. The wider effective coverage of NB816 and NB921, 16.3 deg$^2$ and 16.9 deg$^2$ respectively, are available in the Deep and UltraDeep layers of HSC-SSP from the PDR2. The CHORUS survey provides us with data with additional three NBs (NB527, NB718, and NB973) in the COSMOS field in the UltraDeep layer (1.37 deg$^2$). The five NB datasets allow us to investigate the star-forming galaxies presenting emission-lines at 14 specific redshifts ranging from $z\sim1.6$ down to $z\sim0.05$. We revisit the distribution of large-scale structures and luminosity functions (LFs) for the emission-line galaxies with the large samples of 75,377 emission-line galaxies selected. The redshift revolution of LFs shows that the star formation rate densities (SFRDs) decreases monotonically from $z\sim1.6$, which is consistent with the cosmic SFRD ever known. Our samples of emission-line galaxies covering a sufficiently large survey volume are useful to investigate the evolution of star-forming galaxies since the cosmic noon in a wide range of environments including galaxy clusters, filaments, and voids.

preprint2020arXiv

ALMA uncovers the [CII] emission and warm dust continuum in a z = 8.31 Lyman break galaxy

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

preprint2020arXiv

AstroSat detection of Lyman continuum emission from a z=1.42 galaxy

One of the outstanding problems of current observational cosmology is to understand the nature of sources that produced the bulk of the ionizing radiation after the Cosmic Dark Age. Direct detection of these reionization sources is practically infeasible at high redshift due to the steep decline of intergalactic medium transmission. However, a number of low-redshift analogs emitting Lyman continuum at~900 Angstrom ~restframe are now detected at $z< 0.4$ and there are detections in the range $ 2.5< z< 3.5$ also. Here, we report the detection of Lyman continuum emission with a high escape fraction (>20%) from a low-mass clumpy galaxy at z=1.42, in the middle of the redshift range where no detection has been made before and near the peak of the Cosmic Star-formation history. The observation was made in the Hubble Extreme Deep field by the wide-field Ultra-Violet-Imaging Telescope onboard AstroSat. This is the first detection of Extreme Ultraviolet radiation from a distant galaxy at a rest-frame wavelength of 600 Angstrom, and it opens up a new window to constrain the shape of the ionization spectrum. Further observations with AstroSat should significantly increase the sample of Lyman continuum leaking galaxies at Cosmic Noon.

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

Radiative equilibrium estimates of dust temperature and mass in high-redshift galaxies

Estimating the temperature and mass of dust in high-$z$ galaxies is essential for discussions of the origin of dust in the early Universe. However, this suffers from limited sampling of the infrared spectral-energy distribution. Here we present an algorithm for deriving the temperature and mass of dust in a galaxy, assuming dust to be in radiative equilibrium. We formulate the algorithm for three geometries: a thin spherical shell, a homogeneous sphere, and a clumpy sphere. We also discuss effects of the mass absorption coefficients of dust at ultraviolet and infrared wavelengths, $κ_{\rm UV}$ and $κ_{\rm IR}$, respectively. As an example, we apply the algorithm to a normal, dusty star-forming galaxy at $z=7.5$, A1689zD1, for which three data points in the dust continuum are available. Using $κ_{\rm UV}=5.0\times10^4$ cm$^2$ g$^{-1}$ and $κ_{\rm IR}=30(λ/100μm)^{-β}$ cm$^2$ g$^{-1}$ with $β=2.0$, we obtain dust temperatures of 38--70~K and masses of $10^{6.5-7.3}$ M$_\odot$ for the three geometries considered. We obtain similar temperatures and masses from just a single data point in the dust continuum, suggesting the usefulness of the algorithm for high-$z$ galaxies with limited infrared observations. In the clumpy-sphere case, the temperature becomes equal to that of the usual modified black-body fit, because an additional parameter describing the clumpiness works as an adjuster. The best-fit clumpiness parameter is $ξ_{\rm cl}=0.1$, corresponding to $\sim10$\% of the volume filling factor of the clumps in this high-$z$ galaxy if the clump size is $\sim10$ pc, similar to that of giant molecular clouds in the local Universe.

preprint2020arXiv

The Lyman Continuum Escape Fraction of Galaxies and AGN in the GOODS Fields

We present our analysis of the LyC emission and escape fraction of 111 spectroscopically verified galaxies with and without AGN from $2.26<z<4.3$. We extended our ERS sample from Smith et al. (2018; arXiv:1602.01555) with 64 galaxies in the GOODS North and South fields using WFC3/UVIS F225W, F275W, and F336W mosaics we independently drizzled using the HDUV, CANDELS, and UVUDF data. Among the 17 AGN from the 111 galaxies, one provided a LyC detection in F275W at $m_{AB}=23.19$ mag (S/N $\simeq$ 133) and $GALEX$ NUV at $m_{AB}=23.77$ mag (S/N $\simeq$ 13). We simultaneously fit $SDSS$ and $Chandra$ spectra of this AGN to an accretion disk and Comptonization model and find $f_{esc}$ values of $f_{esc}^{F275W}\simeq 28^{+20}_{-4}$% and $f_{esc}^{NUV}\simeq 30^{+22}_{-5}$%. For the remaining 110 galaxies, we stack image cutouts that capture their LyC emission using the F225W, F275W, and F336W data of the GOODS and ERS samples, and both combined, as well as subsamples of galaxies with and without AGN, and $all$ galaxies. We find the stack of 17 AGN dominate the LyC production from $\langle z\rangle\simeq 2.3-4.3$ by a factor of $\sim$10 compared to all 94 galaxies without AGN. While the IGM of the early universe may have been reionized mostly by massive stars, there is evidence that a significant portion of the ionizing energy came from AGN.

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.

preprint2016arXiv

Detection of an oxygen emission line from a high redshift galaxy in the reionization epoch

The physical properties and elemental abundances of the interstellar medium in galaxies during cosmic reionization are important for understanding the role of galaxies in this process. We report the Atacama Large Millimeter/submillimeter Array detection of an oxygen emission line at a wavelength of 88 micrometers from a galaxy at an epoch about 700 million years after the Big Bang. The oxygen abundance of this galaxy is estimated at about one-tenth that of the Sun. The non-detection of far-infrared continuum emission indicates a deficiency of interstellar dust in the galaxy. A carbon emission line at a wavelength of 158 micrometers is also not detected, implying an unusually small amount of neutral gas. These properties might allow ionizing photons to escape into the intergalactic medium.

preprint2016arXiv

Lyman continuum leaking AGN in the SSA22 field

Subaru/SuprimeCam narrowband photometry of the SSA22 field reveals the presence of four Lyman continuum (LyC) candidates among a sample of 14 AGN. Two show offsets and likely have stellar LyC in nature or are foreground contaminants. The remaining two LyC candidates are Type I AGN. We argue that the average LyC escape fraction of high redshift low luminosity AGN is not likely to be unity, as often assumed in the literature. From direct measurement we obtain the average LyC-to-UV flux density ratio and ionizing emissivity for a number of AGN classes and find it at least a factor of two lower than values obtained assuming f_esc = 1. Comparing to recent Lyα forest measurements, AGNs at redshift z~3 make up at most ~12% and as little as ~5% of the total ionizing budget. Our results suggest that AGNs are unlikely to dominate the ionization budget of the Universe at high redshifts.

preprint2016arXiv

Nebular line emission from z > 7 galaxies in cosmological simulations: rest-frame UV to Optical lines

We have performed very large and high resolution cosmological hydrodynamic simulations in order to investigate detectability of nebular lines in the rest-frame ultraviolet (UV) to optical wavelength range from galaxies at $z > 7$. We find that the expected line fluxes are very well correlated with apparent UV magnitudes. The C IV $1549 \rm Å$ and the C III] $1909 \rm Å$ lines of galaxies brighter than $26~ \rm AB$ are detectable with current facilities such as the VLT/XShooter and the Keck/MOSFIRE. Metal lines such as C IV $1549 \rm Å$, C III] $1909 \rm Å$, [O II] $3727 \rm Å$ and [O III] $4959/5007 \rm Å$ can be good targets for the spectroscopic observation with Thirty Meter Telescope (TMT), European Extremely Large Telescope (E-ELT), Giant Magellan Telescope (GMT) and James Webb Space Telescope (JWST). H$α$ and H$β$ lines are also expected to be detectable with these telescopes. We also predict detectability of the nebular lines for $z > 10$ galaxies that will be found with JWST, Wide-Field Infrared Survey Telescope (WFIRST) and First Light And Reionization Explorer (FLARE) $(11 \leq z \leq 15$). We conclude that the C IV $1549 \rm Å$, C III] $1909 \rm Å$, [O III] $4959/5007 \rm Å$ and H$β$ lines even from $z \sim 15$ galaxies can be strong targets for TMT, ELT and JWST. We also find that the magnification by gravitational lensing is of great help to detect such high-$z$ galaxies. According to our model, the C III] $1909 \rm Å$ line in $z > 9$ galaxy candidates is detectable even using the current facilities.

preprint2015arXiv

Primordial star clusters at extreme magnification

Gravitationally lensed galaxies with magnification ~10-100 are routinely detected at high redshifts, but magnifications significantly higher than this are hampered by a combination of low probability and large source sizes. Magnifications of ~1000 may nonetheless be relevant in the case of intrinsically small, high-redshift objects with very high number densities. Here, we explore the prospects of detecting compact (< 10 pc), high-redshift (z > 7) Population III star clusters at such extreme magnifications in large-area surveys with planned telescopes like Euclid, WFIRST and WISH. We find that the planned WISH 100 sq. deg ultradeep survey may be able to detect a small number of such objects, provided that the total stellar mass of these star clusters is > 10000 solar masses. If candidates for such lensed Population III star clusters are found, follow-up spectroscopy of the surrounding nebula with the James Webb Space Telescope or groundbased Extremely Large Telescopes should be able to confirm the Population III nature of these objects. Multiband photometry of these objects with the James Webb Space Telescope also has the potential to confirm that the stellar initial mass function in these Population III star clusters is top-heavy, as supported by current simulations.

preprint2014arXiv

An updated analytic model for the attenuation by the intergalactic medium

We present an updated version of the so-called Madau model for the attenuation by the intergalactic neutral hydrogen against the radiation from distant objects. First, we derive a distribution function of the intergalactic absorbers from the latest observational statistics of the Ly$α$ forest, Lyman limit systems, and damped Ly$α$ systems. The distribution function excellently reproduces the observed redshift evolutions of the Ly$α$ depression and of the mean-free-path of the Lyman continuum simultaneously. Then, we derive a set of the analytic functions which describe the mean intergalactic attenuation curve for objects at $z>0.5$. The new model predicts less (or more) Ly$α$ attenuation for $z\simeq3$--5 ($z>6$) sources through usual broad-band filters relative to the original Madau model. This may cause a systematic difference in the photometric redshift estimates, which is, however, still as small as about 0.05. Finally, we find a more than 0.5 mag overestimation of the Lyman continuum attenuation in the original Madau model at $z>3$, which causes a significant overcorrection against direct observations of the Lyman continuum of galaxies.

preprint2014arXiv

Physical Properties of UDF12 Galaxies in Cosmological simulations

We have performed a large cosmological hydrodynamics simulation tailored to the deep survey with the Hubble Space Telescope made in 2012, the so-called UDF12 campaign. After making a light-cone output, we have applied the same color selection criteria as the UDF12 campaign to select galaxies from our simulation, and then, have examined the physical properties of them as a proxy of the real observed UDF12 galaxies at $z > 7$. As a result, we find that the halo mass is almost linearly proportional to the observed ultraviolet (UV) luminosity ($4 \times 10^{11}~{\rm M_{\odot}}$ at $M_{\rm UV} = -21$). The dust attenuation and UV slope $β$ well correlates with the observed UV luminosity, which is consistent with observations quantitatively. The star formation rate (SFR) is also linearly proportional to the stellar mass and the specific SFR shows only a weak dependency on the mass. We also find an increasing star formation history with a time-scale of $\sim100$ Myr in the high-$z$ galaxies. An average metallicity weighted by the Lyman continuum luminosity reaches up to $>0.1$ Solar even at $z \sim 10$, suggesting a rapid metal enrichment. We also expect $\geq 0.1$ mJy at 350 GHz of the dust thermal emission from the galaxies with $H_{160} \leq 27$, which can be detectable with the Atacama Large Milimetre-submilimetre Array. The galaxies selected by the UDF12 survey contribute to only $52--12\%$ of the cosmic SFR density from $z \sim 7$ to $z \sim 10$, respectively. The James Webb Space Telescope will push the detection fraction up to $77--72\%$.

preprint2014arXiv

The critical radiation intensity for direct collapse black hole formation: dependence on the radiation spectral shape

It has been proposed that supermassive black holes (SMBHs) are originated from direct-collapse black holes (DCBHs) that are formed at z gtrsim 10 in the primordial gas in the case that H2 cooling is suppressed by strong external radiation. In this work, we study the critical specific intensity J^crit required for DCBH formation for various radiation spectral shapes by a series of one-zone calculations of a collapsing primordial- gas cloud. We calculate the critical specific intensity at the Lyman-Werner (LW) bands J^crit_LW,21 (in units of 10^-21 erg s^-1 Hz^-1 sr^-1 cm^-2) for realistic spectra of metal-poor galaxies. We find J^crit is not sensitive to the age or metallicity for the constant star formation galaxies with J^crit_LW,21 = 1300-1400, while J^crit decreases as galaxies become older or more metal-enriched for the instantaneous starburst galaxies. However, such dependence for the instantaneous starburst galaxies is weak for the young or extremely metal-poor galaxies: J^crit_LW,21 = 1000-1400 for the young galaxies and J^crit_LW,21 approx 1400 for the extremely metal-poor galaxies. The typical value of J^crit for the realistic spectra is higher than those expected in the literature, which affects the estimated DCBH number density n_DCBH. By extrapolating the result of Dijkstra, Ferrara and Mesinger, we obtain n_DCBH sim 10^-10 cMpc^-3 at z = 10, although there is still large uncertainty in this estimation.

preprint2014arXiv

The environments of Ly$α$ blobs I: Wide-field Ly$α$ imaging of TN J1338-1942, a powerful radio galaxy at $z\simeq 4.1$ associated with a giant Ly$α$ nebula

We exploit wide-field Ly$alpha$ imaging with Subaru to probe the environment around TN J1338-1942, a powerful radio galaxy with a >100 kpc Ly$α$ halo at z=4.11. We used a sample of Ly$α$ emitters (LAEs) down to $\log(L_{\rm Lyα} [erg\, s^{-1}])\sim 42.8$ to measure the galaxy density around TNJ1338, compared to a control sample from a blank field taken with the same instrument. We found that TNJ1338 resides in a region with a peak overdensity of $δ_{\rm LAE}=2.8\pm 0.5$ on scales of $8\, h^{-1}\rm Mpc$ (on the sky) and $112\, h^{-1}\rm Mpc$ (line of sight) in comoving coordinates. Adjacent to this overdensity, we found a strong underdensity where virtually no LAEs are detected. We used a semi-analytical model of LAEs derived from the Millennium Simulation to compare our results with theoretical predictions. While the theoretical density distribution is consistent with the blank field, overdense regions such as that around TNJ1338 are very rare, with a number density of $6.4\times 10^{-8}\rm Mpc^{-3}$ (comoving), corresponding to the densest < 0.4 percentile at $z\simeq 4.1$. We also found that the Ly$α$ luminosity function in the TNJ1338 field differs from that in the blank field: the number of bright LAEs ($\log(L_{\rm Lyα}[erg\,s^{-1}]) \gtrsim 43.3$) is enhanced, while the number of fainter LAEs is relatively suppressed. These results suggest that some powerful radio galaxies associated with Ly$α$ nebulae reside in extreme overdensities on $\sim 3$--$6\, \rm Mpc$ scales, where star-formation and AGN activity may be enhanced via frequent galaxy mergers or high rates of gas accretion from the surroundings.

preprint2014arXiv

Variability of Disk Emission in Pre-Main Sequence and Related Stars. III. Exploring Structural Changes in the Pre-transitional Disk in HD 169142

We present near-IR and far-UV observations of the pre-transitional (gapped) disk in HD 169142 using NASA's Infrared Telescope Facility and Hubble Space Telescope. The combination of our data along with existing data sets into the broadband spectral energy distribution reveals variability of up to 45% between ~1.5-10 μm over a maximum timescale of 10 years. All observations known to us separate into two distinct states corresponding to a high near-IR state in the pre-2000 epoch and a low state in the post-2000 epoch, indicating activity within the <1 AU region of the disk. Through analysis of the Pa β and Br γ lines in our data we derive a mass accretion rate in May 2013 of (1.5 - 2.7) x 10^-9 Msun/yr. We present a theoretical modeling analysis of the disk in HD 169142 using Monte-Carlo radiative transfer simulation software to explore the conditions and perhaps signs of planetary formation in our collection of 24 years of observations. We find that shifting the outer edge (r = 0.3 AU) of the inner disk by 0.05 AU toward the star (in simulation of accretion and/or sculpting by forming planets) successfully reproduces the shift in NIR flux. We establish that the ~40-70 AU dark ring imaged in the NIR by Quanz et al. (2013) and Momose et al. (2013) and at 7 mm by Osorio et al. (2014) may be reproduced with a 30% scaled density profile throughout the region, strengthening the link to this structure being dynamically cleared by one or more planetary mass bodies.

preprint2013arXiv

ALMA will determine the spectroscopic redshift z>8 with FIR [O III] emission lines

We investigate the potential use of nebular emission lines in the rest-frame far-infrared (FIR) for determining spectroscopic redshift of z>8 galaxies with the Atacama Large Millimeter/sub-millimeter Array (ALMA). After making a line emissivity model as a function of metallicity, especially for the [O III] 88 micron line which is likely to be the strongest FIR line from H II regions, we predict the line fluxes from high-z galaxies based on a cosmological hydrodynamics simulation of galaxy formation. Since the metallicity of galaxies reaches at ~0.2 Zsun even at z>8 in our simulation, we expect the [O III] 88 micron line as strong as 1.3 mJy for 27 AB objects, which is detectable at a high significance by <1 hour integration with ALMA. Therefore, the [O III] 88 micron line would be the best tool to confirm the spectroscopic redshifts beyond z=8.

preprint2013arXiv

Revisiting the Cosmic Star Formation History: Caution on the Uncertainties in Dust Correction and Star Formation Rate Conversion

The cosmic star formation rate density (CSFRD) has been observationally investigated out to redshift z~10. However, most of theoretical models for galaxy formation underpredict the CSFRD at z>1. Since the theoretical models reproduce the observed luminosity functions (LFs), luminosity densities (LDs), and stellar mass density at each redshift, this inconsistency does not simply imply that theoretical models should incorporate some missing unknown physical processes in galaxy formation. Here, we examine the cause of this inconsistency in UV wavelengths by using a mock catalog of galaxies generated by a semi-analytic model of galaxy formation. We find that this inconsistency is due to two observational uncertainties: dust obscuration correction and conversion from UV luminosity to star formation rate (SFR). The methods for correction of obscuration and SFR conversion used in observational studies result in the overestimation of CSFRD by ~ 0.1-0.3 dex and ~ 0.1-0.2 dex, respectively, compared to the results obtained directly from our mock catalog. We present new empirical calibrations for dust attenuation and conversion from observed UV LFs and LDs into CSFRD.

preprint2013arXiv

The Effect of the Remnant Mass in Estimating Stellar Mass of Galaxies

The definition of the galactic stellar mass estimated from the spectral energy distribution is ambiguous in the literature; whether the stellar mass includes the mass of the stellar remnants, i.e. white dwarfs, neutron stars, and black holes, is not well described. The remnant mass fraction in the total (living+remnant) stellar mass of a simple stellar population monotonically increases with the age of the population, and the initial mass function and metallicity affect the increasing rate. Since galaxies are composed of a number of stellar populations, the remnant mass fraction may depend on the total stellar mass of galaxies in a complex way. As a result, the shape of the stellar mass function of galaxies may change, depending on the definition of the stellar mass. In order to explore this issue, we ran a cosmological hydrodynamical simulation, and then, we have found that the remnant mass fraction indeed correlates with the total stellar mass of galaxies. However, this correlation is weak and the remnant fraction can be regarded as a constant which depends only on the redshift. Therefore, the shape of the stellar mass function is almost unchanged, but it simply shifts horizontally if the remnant mass includes or not. The shift is larger at lower redshift and it reaches 0.2-dex at z=0 for a Chabrier IMF. Since this causes a systematic difference, we should take care of the definition of the `stellar' mass, when comparing one's result with others.

preprint2012arXiv

A possible signature of primordial stellar populations in $z=3$ Lyman $α$ emitters

Observations with Subaru telescope have detected surprisingly strong Lyman continuum (LyC; $\sim900$ Å in the rest-frame) from some Lyman $α$ emitters (LAEs) at $z=3.1$. We have examined the stellar population which simultaneously accounts for the strength of the LyC and the spectral slope of non-ionizing ultraviolet of the LAEs. As a result, we have found that stellar populations with metallicity $Z\geq1/50 Z_\odot$ can explain the observed LyC strength only with a very top-heavy initial mass function (IMF; $< m > \sim 50 M_\odot$). However, the critical metallicity for such an IMF is expected to be much lower. A very young ($\sim1$ Myr) and massive ($\sim100$ $M_\odot$) extremely metal-poor ($Z\leq5\times10^{-4}Z_\odot$) or metal-free (so-called Population III) stellar population can also reproduce the observed LyC strength if the mass fraction of such 'primordial' stellar population is $\sim1$% in total stellar mass of the LAEs.

preprint2012arXiv

Dust formation history of galaxies: a critical role of metallicity for the dust mass growth by accreting materials in the interstellar medium

This paper investigate what is the main driver of the dust mass growth in the interstellar medium (ISM) by using a chemical evolution model of galaxy with metals (elements heavier than helium) in dust phase in addition to the total amount of metals. We consider asymptotic giant branch (AGB) stars, type II supernovae (SNe II) and the dust mass growth in the ISM as the sources of dust, and SN shocks as the destruction mechanism of dust. Further, to describe the dust evolution precisely, our model takes into account the age and metallicity (the ratio of metal mass to ISM mass) dependence of the sources of dust. We particularly focused on the dust mass growth, and found that the dust mass growth in the ISM is regulated by the metallicity. To quantify this aspect, we introduce a "critical metallicity", which is a metallicity at which the contribution of stars (AGB stars and SNe II) equals that of the dust mass growth in the ISM. If the star formation timescale is shorter, the value of the critical metallicity is higher, but the galactic age at which the metallicity reaches the critical metallicity is shorter. From observations, it was expected that the dust mass growth was the dominant source of dust in the Milky Way and dusty QSOs at high redshifts. By introducing the critical metallicity, it is clearly shown that the dust mass growth is the main source of dust in such galaxies with various star formation timescales and ages. The dust mass growth in the ISM is regulated by metallicity, and we stress that the critical metallicity works as an indicator to judge whether the grain growth in the ISM is dominant source of dust in a galaxy, especially because of the strong and nonlinear dependence on the metallicity.

preprint2012arXiv

Effect of Photodesorption on Snow Line at the Surface of Optically Thick Circumstellar Disks around Herbig Ae/Be Stars

We investigate the effect of photodesorption on the snow line position at the surface of a protoplanetary disk around a Herbig Ae/Be star, motivated by the detection of water ice particles at the surface of the disk around HD142527 by Honda et al. For this aim, we obtain the density and temperature structure in the disk with a 1+1D radiative transfer and determine the distribution of water ice particles in the disk by the balance between condensation, sublimation, and photodesorption. We find that photodesorption induced by the far-ultraviolet radiation from the central star depresses the ice-condensation front toward the mid-plane and pushes the surface snow line outward significantly when the stellar effective temperature exceeds a certain critical value. This critical effective temperature depends on the stellar luminosity and mass, the water abundance in the disk, and the yield of photodesorption. We present an approximate analytic formula for the critical temperature. We separate Herbig Ae/Be stars into two groups on the HR diagram according to the critical temperature; one is the disks where photodesorption is effective and from which we may not find ice particles at the surface, and the other is the disks where photodesorption is not effective. We estimate the snow line position at the surface of the disk around HD142527 to be 100--300 AU which is consistent with the water ice detection at >140 AU in the disk. All results depend on the dust grain size by a complex way and this point requires more work in future.

preprint2012arXiv

Reexamination of the Infrared Excess-Ultraviolet Slope Relation of Local Galaxies

The relation between the ratio of infrared (IR) and ultraviolet (UV) flux densities (the infrared excess: IRX) and the slope of the UV spectrum (β) of galaxies plays a fundamental role in the evaluation of the dust attenuation of star forming galaxies especially at high redshifts. Many authors, however, pointed out that there is a significant dispersion and/or deviation from the originally proposed IRX-βrelation depending on sample selection. We reexamined the IRX-βrelation by measuring the far- and near-UV flux densities of the original sample galaxies with GALEX and AKARI imaging data, and constructed a revised formula. We found that the newly obtained IRX values were lower than the original relation because of the significant underestimation of the UV flux densities of the galaxies, caused by the small aperture of IUE, Further, since the original relation was based on IRAS data which covered a wavelength range of λ= 42--122\mum, using the data from AKARI which has wider wavelength coverage toward longer wavelengths, we obtained an appropriate IRX-βrelation with total dust emission (TIR): \log(L_{\rm TIR}/L_{\rm FUV}) = \log [10^{0.4(3.06+1.58β)}-1] +0.22. This new relation is consistent with most of the preceding results for samples selected at optical and UV, though there is a significant scatter around it. We also found that even the quiescent class of IR galaxies follows this new relation, though luminous and ultraluminous IR galaxies distribute completely differently as well known before.

preprint2012arXiv

The origin of dust in galaxies revisited: the mechanism determining dust content

The origin of cosmic dust is a fundamental issue in planetary science. This paper revisits the origin of dust in galaxies, in particular, in the Milky Way, by using a chemical evolution model of a galaxy composed of stars, interstellar medium, metals (elements heavier than helium), and dust. We start from a review of time-evolutionary equations of the four components, and then, we present simple recipes for the stellar remnant mass and yields of metal and dust based on models of stellar nucleosynthesis and dust formation. After calibrating some model parameters with the data from the solar neighborhood, we have confirmed a shortage of the stellar dust production rate relative to the dust destruction rate by supernovae if the destruction efficiency suggested by theoretical works is correct. If the dust mass growth by material accretion in molecular clouds is active, the observed dust amount in the solar neighborhood is reproduced. We present a clear analytic explanation of the mechanism for determining dust content in galaxies after the activation of accretion growth: a balance between accretion growth and supernova destruction. Thus, the dust content is independent of the uncertainty of the stellar dust yield after the growth activation. The timing of the activation is determined by a critical metal mass fraction which depends on the growth and destruction efficiencies. The solar system formation seems to have occurred well after the activation and plenty of dust would have existed in the proto-solar nebula.

preprint2011arXiv

Rest-frame ultraviolet-to-optical spectral characteristics of extremely metal-poor and metal-free galaxies

Finding the first generation of galaxies in the early Universe is the greatest step forward for understanding galaxy formation and evolution. For strategic survey of such galaxies and interpretation of the obtained data, this paper presents an ultraviolet-to-optical spectral model of galaxies with a great care of the nebular emission. In particular, we present a machine-readable table of intensities of 119 nebular emission lines from Ly$α$ to the rest-frame 1 $μ$m as a function of metallicity from zero to the Solar one. Based on the spectral model, we present criteria of equivalent widths of Ly$α$, He {\sc ii} $\lambda1640$, H$α$, H$β$, [O {\sc iii}] $\lambda5007$ to select extremely metal-poor and metal-free galaxies although these criteria have uncertainty caused by the Lyman continuum escape fraction and the star formation duration. We also present criteria of broad-band colours which will be useful to select candidates for spectroscopic follow-up from drop-out galaxies. We propose the line intensity ratio of [O {\sc iii}] $\lambda5007$ to H$β$ $<0.1$ as the most robust criterion for $<1/1000$ of the Solar metallicity. This ratio of a galaxy with a few $M_\odot$ yr$^{-1}$ at $z\sim8$ is detectable by spectroscopy with the James Webb Space Telescope within a reasonable exposure time.

preprint2010arXiv

Lyman `bump' galaxies - II. A possible signature of massive extremely metal-poor or metal-free stars in z=3.1 Ly-alpha emitters

(Abridged) Deep NB359 imaging with Subaru by Iwata et al. have detected surprisingly strong Lyman continuum (LyC; ~900A in the rest-frame) from some LAEs at z=3.1. However, the redshifts might be misidentified due to a narrow wavelength coverage in previous spectroscopy. We here present new deep spectroscopy covering the observed 4,000-7,000A with VLT/VIMOS and Subaru/FOCAS of 8 LAEs detected in NB359. All the 8 objects have only one detectable emission line around 4,970A which is most likely to be Ly-A at z=3.1, and thus, the objects are certainly LAEs at the redshift. However, 5 of them show a ~0.''8 spatial offset between the Ly-A emission and the source detected in NB359. No indications of the redshifts of the NB359 sources are found although it is statistically difficult that all the 5 LAEs have a foreground object accounting for the NB359 flux. The rest 3 LAEs show no significant offset from the NB359 position. Therefore, they are truly LyC emitting LAEs at z=3.1. We also examine the stellar population which simultaneously accounts for the strength of the LyC and the spectral slope of non-ionizing ultraviolet of the LAEs. We consider the latest statistics of Lyman limit systems to estimate the LyC optical depth in the IGM and an additional contribution of the bound-free LyC from photo-ionized nebulae to the LyC emissivity. As a result, we find that stellar populations with metallicity Z>=1/50Z_sun can explain the observed LyC strength only with a very top-heavy initial mass function (IMF; <m>~50 M_sun). However, the critical metallicity for such an IMF is expected to be much lower. A very young (~1 Myr) and massive (~100 M_sun) extremely metal-poor (Z<=5e-4Z_sun) or metal-free (so-called Population III) stellar population can reproduce the observed LyC strength. The required mass fraction of such `primordial' stellar population is ~1--10% in total stellar mass of the LAEs.

preprint2009arXiv

Lyman "bump" galaxies - I. Spectral energy distribution of galaxies with an escape of nebular Lyman continuum

It is essential to know galactic emissivity and spectrum of Lyman continuum (LyC) in order to understand the cosmic reionization. Here we consider an escape of nebular LyC from galaxies and examine the consequent spectral energy distribution. It is usually assumed that hydrogen nebular LyC mostly produced by bound-free transitions is consumed within photo-ionized nebulae (so-called "on-the-spot" approximation). However, an escape of the continuum should be taken into account if stellar LyC escapes from galaxies through "matter-bounded" nebulae. We show that the escaping hydrogen bound-free LyC makes a strong bump just below the Lyman limit. Such a galaxy would be observed as a Lyman "bump" galaxy. This bump results from the radiation energy re-distribution of stellar LyC by nebulae. The strength of the bump depends on electron temperature in nebulae, escape fraction of stellar and nebular LyC, hardness of stellar LyC (i.e. metallicity, initial mass function, age and star formation history), and IGM attenuation. We can use the bump to find very young ($\sim1$ Myr), massive ($\sim100$ $M_\odot$), and extremely metal-poor (or metal-free) stellar populations at $z<4$. Because of the bump, 900 Å-to-1500 Åluminosity density ratio (per Hz) becomes maximum (2--3 times larger than the stellar intrinsic ratio) when about 40% of the stellar LyC is absorbed by nebulae. The total number of escaping LyC photons increases due to the escape of nebular LyC but does not exceed the stellar intrinsic one. The radiation energy re-distribution by nebulae reduces the mean energy of escaping LyC only by $\approx10$% relative to that of stellar LyC. Therefore, the effect of the escape of nebular LyC on the reionization process may be small.

preprint2008arXiv

Observational Possibility of the "Snow Line" on the Surface of Circumstellar Disks with the Scattered Light

We discuss how we obtain the spatial distribution of ice on the surface of the circumstellar disk around young stars. Ice in the disks plays a very important role in various issues, for instance, on the disk structure, on the planet formation, on the isotopic anomaly in meteorites, and on the origin of the sea on the Earth. Therefore, the spatially resolved observation of the condensation/sublimation front of ice, so-called ``snow line'' is strongly required. Here, we propose a new method for obtaining the spatially resolved ``snow line'' on the circumstellar disks by observing 3 \micron H$_2$O ice feature in the scattered light. Based on radiative transfer considerations, we show that the feature is clearly imprinted in the spectrum of the scattered light from both optically thick and thin circumstellar disks. We also show that the scattered light and the H$_2$O ice feature from protoplanetary disks are detectable and spatially resolvable with the current instruments through a $H_2O$ narrowband filter around 3 \micron. Finally, we present a diagnostics of disk dust properties on the $K-H_2O$ and $K-L'$ two color diagram.

preprint2003arXiv

Evolution of Dust-to-Metal Ratio in Galaxies

This paper investigates the evolution of the dust-to-metal ratio in galaxies based on a simple evolution model for the amount of metal and dust with infall. We take into account grain formation in stellar mass-loss gas, grain growth by the accretion of metallic atoms in a cold dense cloud, and grain destruction by SNe shocks. Especially, we propose that the accretion efficiency is independent of the star-formation history. This predicts various evolutionary tracks in the metallicity ($Z$)--dust-to-gas ratio ($\cal D$) plane depending on the star-formation history. In this framework, the observed linear $Z$--$\cal D$ relation of nearby spiral galaxies can be interpreted as a sequence of a constant galactic age. We emphasize that an observational study of the $Z$--$\cal D$ relation of galaxies at $z\sim 1$ is very useful to constrain the efficiencies of dust growth and destruction. We also suggest that the Lyman break galaxies at $z\sim 3$ have a very low dust-to-metal ratio, typically $\ltsim 0.1$. Although the effect of infall on the evolutionary tracks in the $Z$--$\cal D$ plane is quite small, the dispersion of the infall rate can disturb the $Z$--$\cal D$ relation with a constant galactic age.

preprint2002arXiv

Star Formation Rate from Dust Infrared Emission

We examine what types of galaxies the conversion formula from dust infrared (IR) luminosity into the star formation rate (SFR) derived by Kennicutt (1998) is applicable to. The ratio of the observed IR luminosity, $L_{\rm IR}$, to the intrinsic bolometric luminosity of the newly ($\la$ 10 Myr) formed stars, $L_{\rm SF}$, of a galaxy can be determined by a mean dust opacity in the interstellar medium and the activity of the current star formation. We find that these parameters area being $0.5 \le L_{\rm IR}/L_{\rm SF} \le 2.0$ is very large, and many nearby normal and active star-forming galaxies really fall in this area. It results from offsetting two effects of a small dust opacity and a large cirrus contribution of normal galaxies relative to starburst galaxies on the conversion of the stellar emission into the dust IR emission. In conclusion, the SFR determined from the IR luminosity under the assumption of $L_{\rm IR}=L_{\rm SF}$ like Kennicutt (1998) is reliable within a factor of 2 for all galaxies except for dust rich but quiescent galaxies and extremely dust poor galaxies.