Source author record

Takuya Hashimoto

Takuya Hashimoto 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

22works
2topics
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

22 published item(s)

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

ALMA reveals extended cool gas and hot ionized outflows in a typical star-forming galaxy at $z=7.13$

We present spatially-resolved morphological properties of [CII] 158 $μ$m, [OIII] 88 $μ$m, dust, and rest-frame ultraviolet (UV) continuum emission for A1689-zD1, a strongly lensed, sub-L* galaxy at $z=7.13$, by utilizing deep Atacama Large Millimeter/submillimeter Array (ALMA) and Hubble Space Telescope (HST) observations. While the [OIII] line and UV continuum are compact, the [CII] line is extended up to a radius of $r \sim 12$ kpc. Using multi-band rest-frame far-infrared (FIR) continuum data ranging from 52-400 $μ$m, we find an average dust temperature and emissivity index of $T_{\rm dust} = 41^{+17}_{-14}$ K and $β= 1.7^{+1.1}_{-0.7}$, respectively, across the galaxy. We find slight differences in the dust continuum profiles at different wavelengths, which may indicate that the dust temperature decreases with distance. We map the star-formation rate (SFR) via IR and UV luminosities and determine a total SFR of $37\pm 1~M_\odot~{\rm yr}^{-1}$ with an obscured fraction of $87\%$. While the [OIII] line is a good tracer of the SFR, the [CII] line shows deviation from the local $L_{\rm [CII]}$-SFR relations in the outskirts of the galaxy. Finally, we observe a clear difference in the line profile between [CII] and [OIII], with significant residuals ($\sim 5σ$) in the [OIII] line spectrum after subtracting a single Gaussian model. This suggests a possible origin of the extended [CII] structure from the cooling of hot ionized outflows. The extended [CII] and high-velocity [OIII] emission may both contribute in part to the high $L_{\rm [OIII]}$/$L_{\rm [CII]}$ ratios recently reported in $z>6$ galaxies.

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

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

preprint2021arXiv

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

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

preprint2021arXiv

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

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

Elevated ionizing photon production efficiency in faint high-equivalent-width Lyman-alpha emitters

While low-luminosity galaxies dominate number counts at all redshifts, their contribution to cosmic Reionization is poorly understood due to a lack of knowledge of their physical properties. We isolate a sample of 35 z~4-5 continuum-faint Lyman-alpha emitters from deep VLT/MUSE spectroscopy and directly measure their Halpha emission using stacked Spitzer/IRAC Ch. 1 photometry. Based on Hubble Space Telescope imaging, we determine that the average UV continuum magnitude is fainter than -16 (~0.01 L_star), implying a median Lyman-alpha equivalent width of 249 Angstroms. By combining the Halpha measurement with the UV magnitude we determine the ionizing photon production efficiency, xi_ion, a first for such faint galaxies. The measurement of log (xi_ion [Hz/erg]) = 26.28 (+0.28; -0.40) is in excess of literature measurements of both continuum- and emission line-selected samples, implying a more efficient production of ionizing photons in these lower-luminosity, Lyman-alpha-selected systems. We conclude that this elevated efficiency can be explained by stellar populations with metallicities between 4e-4 and 0.008, with light-weighted ages less than 3 Myr.

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 MUSE Hubble Ultra Deep Field Survey XIII. Spatially resolved spectral properties of Lyman alpha haloes around star-forming galaxies at z > 3

We present spatially resolved maps of six individually-detected Lyman alpha haloes (LAHs) as well as a first statistical analysis of the Lyman alpha (Lya) spectral signature in the circum-galactic medium of high-redshift star-forming galaxies using MUSE. Our resolved spectroscopic analysis of the LAHs reveals significant intrahalo variations of the Lya line profile. Using a three-dimensional two-component model for the Lya emission, we measure the full width at half maximum (FWHM), the peak velocity shift and the asymmetry of the Lya line in the core and in the halo of 19 galaxies. We find that the Lya line shape is statistically different in the halo compared to the core for ~40% of our galaxies. Similarly to object-by-object based studies and a recent resolved study using lensing, we find a correlation between the peak velocity shift and the width of the Lya line both at the interstellar and circum-galactic scales. While there is a lack of correlation between the spectral properties and the spatial scale lengths of our LAHs, we find a correlation between the width of the line in the LAH and the halo flux fraction. Interestingly, UV bright galaxies show broader, more redshifted and less asymmetric Lya lines in their haloes. The most significant correlation found is for the FWHM of the line and the UV continuum slope of the galaxy, suggesting that the redder galaxies have broader Lya lines. The generally broad and red line shapes found in the halo component suggests that the Lya haloes are powered either by scattering processes through an outflowing medium, fluorescent emission from outflowing cold clumps of gas, or a mix of both. Considering the large diversity of the Lya line profiles observed in our sample and the lack of strong correlation, the interpretation of our results is still broadly open and underlines the need for realistic spatially resolved models of the LAHs.

preprint2020arXiv

The MUSE Hubble Ultra Deep Field Survey XIV. The evolution of the Lya emitter fraction from z=3 to z=6

The Lya emitter (LAE) fraction, X_LAE, is a potentially powerful probe of the evolution of the intergalactic neutral hydrogen gas fraction. However, uncertainties in the measurement of X_LAE are still debated. Thanks to deep data obtained with MUSE, we can measure the evolution of X_LAE homogeneously over a wide redshift range of z~3-6 for UV-faint galaxies (down to M_1500~-17.75). This is significantly fainter than in former studies, and allows us to probe the bulk of the population of high-z star-forming galaxies. We construct a UV-complete photo-redshift sample following UV luminosity functions and measure the Lya emission with MUSE using the second data release from the MUSE HUDF Survey. We derive the redshift evolution of X_LAE for M_1500 in [-21.75;-17.75] for the first time with a equivalent width range EW(Lya)>=65 A and find low values of X_ LAE<~30% at z<~6. For M_1500 in [-20.25;-18.75] and EW(Lya)<~25 A, our X_LAE values are consistent with those in the literature within 1sigma at z<~5, but our median values are systematically lower than reported values over the whole redshift range. In addition, we do not find a significant dependence of X_LAE on M_1500 for EW(Lya)>~50 A at z~3-4, in contrast with previous work. The differences in X_LAE mainly arise from selection biases for Lyman Break Galaxies (LBGs) in the literature: UV-faint LBGs are more easily selected if they have strong Lya emission, hence X_LAE is biased towards higher values. Our results suggest either a lower increase of X_LAE towards z~6 than previously suggested, or even a turnover of X_LAE at z~5.5, which may be the signature of a late or patchy reionization process. We compared our results with predictions from a cosmological galaxy evolution model. We find that a model with a bursty star formation (SF) can reproduce our observed X_LAE much better than models where SF is a smooth function of time.

preprint2020arXiv

The MUSE Hubble Ultra Deep Field Survey. XV. The mean rest-UV spectra of Ly-alpha emitters at z>3

We investigate the ultraviolet (UV) spectral properties of faint Lyman-$α$ emitters (LAEs) in the redshift range 2.9<z<4.6 and provide material to prepare future observations of the faint Universe. We use data from the MUSE Hubble Ultra Deep Survey to construct mean rest-frame spectra of continuum-faint (median M$_{UV}$ of -18 and down to M$_{UV}$ of -16), low stellar mass (median value of $10^{8.4}$ and down to $10^{7}M_{\odot}$) LAEs at redshift z>3. We compute various averaged spectra of LAEs sub-sampled on the basis of their observational (e.g., Ly$α$ strength, UV magnitude and spectral slope) and physical (e.g., stellar mass and star-formation rate) properties. We search for UV spectral features other than Ly$α$, such as higher-ionization nebular emission lines and absorption features. We successfully observe the OIII]1666 and CIII]909 collisionally excited emission lines and the HeII1640 recombination feature, as well as the resonant CIV1550 doublet either in emission or P-Cygni. We compare the observed spectral properties of the different mean spectra and find the emission lines to vary with the observational and physical properties of the LAEs. In particular, the mean spectra of LAEs with larger Ly$α$ equivalent widths, fainter UV magnitudes, bluer UV spectral slopes and lower stellar masses show the strongest nebular emission. The line ratios of these lines are similar to those measured in the spectra of local metal-poor galaxies, while their equivalent widths are weaker compared to the handful of extreme values detected in individual spectra of z>2 galaxies. This suggests that weak UV features are likely ubiquitous in high z, low-mass and faint LAEs. We publicly release the stacked spectra as they can serve as empirical templates for the design of future observations, such as those with the James Webb Space Telescope and the Extremely Large Telescope.

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

preprint2018arXiv

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

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

preprint2016arXiv

A Close Comparison between Observed and Modeled Lyα Lines for z ~ 2.2 Lyman Alpha Emitters

We present the results of a Lya profile analysis of 12 Lya emitters (LAEs) at z = 2.2 with high-resolution Lya spectra. We find that all 12 objects have a Lya profile with the main peak redward of the systemic redshift defined by nebular lines, and five have a weak, secondary peak blueward of the systemic redshift (blue bump). The average velocity offset of the red main peak (the blue bump, if any) with respect to the systemic redshift is Delta_v_Lya,r = 174+/- 19 km s-1 (Delta_v_Lya,b = -316+/-45 km s-1), which is smaller than (comparable to) that of Lyman-break galaxies (LBGs). The outflow velocities inferred from metal absorption lines in three individual and one stacked spectra are comparable to those of LBGs. The expanding shell model constructed by Verhamme et al. (2006) reproduces not only the Lya profiles but also other observed quantities including the outflow velocity and the FWHM of nebular lines for the non-blue bump objects. On the other hand, the model predicts too high FWHMs of nebular lines for the blue bump objects, although this discrepancy may disappear if we introduce additional Lya photons produced by gravitational cooling. We show that the small Delta_v_Lya,r values of our sample can be explained by low neutral-hydrogen column densities of log(NHI) = 18.9 cm-2 on average. This value is more than one order of magnitude lower than those of LBGs but is consistent with recent findings that LAEs have high ionization parameters and low Hi gas masses. This result suggests that low NHI values, giving reduced numbers of resonant scattering of Lya photons, are the key to the strong Lya emission of LAEs.

preprint2016arXiv

Ly$α$ Emitters with Very Large Ly$α$ Equivalent Widths, EW$_{\rm 0}$(Ly$α$) $\simeq 200-400$ Å, at $z\sim 2$

We present physical properties of spectroscopically confirmed Ly$α$ emitters (LAEs) with very large rest-frame Ly$α$ equivalent widths EW$_{\rm 0}$(Ly$α$). Although the definition of large EW$_{\rm 0}$(Ly$α$) LAEs is usually difficult due to limited statistical and systematic uncertainties, we identify six LAEs selected from $\sim 3000$ LAEs at $z\sim 2$ with reliable measurements of EW$_{\rm 0}$ (Ly$α$) $\simeq 200-400$ Å given by careful continuum determinations with our deep photometric and spectroscopic data. These large EW$_{\rm 0}$(Ly$α$) LAEs do not have signatures of AGN, but notably small stellar masses of $M_{\rm *} = 10^{7-8}$ $M_{\rm \odot}$ and high specific star-formation rates (star formation rate per unit galaxy stellar mass) of $\sim 100$ Gyr$^{-1}$. These LAEs are characterized by the median values of $L({\rm Lyα})=3.7\times 10^{42}$ erg s$^{-1}$ and $M_{\rm UV}=-18.0$ as well as the blue UV continuum slope of $β= -2.5\pm0.2$ and the low dust extinction $E(B-V)_{\rm *} = 0.02^{+0.04}_{-0.02}$, which indicate a high median Ly$α$ escape fraction of $f_{\rm esc}^{\rm Lyα}=0.68\pm0.30$. This large $f_{\rm esc}^{\rm Lyα}$ value is explained by the low {\sc Hi} column density in the ISM that is consistent with FWHM of the Ly$α$ line, ${\rm FWHM (Lyα)}=212\pm32$ km s$^{-1}$, significantly narrower than those of small EW$_{\rm 0}$(Ly$α$) LAEs. Based on the stellar evolution models, our observational constraints of the large EW$_{\rm 0}$ (Ly$α$), the small $β$, and the rest-frame He{\sc ii} equivalent width imply that at least a half of our large EW$_{\rm 0}$(Ly$α$) LAEs would have young stellar ages of $\lesssim 20$ Myr and very low metallicities of $Z<0.02 Z_\odot$ regardless of the star-formation history.

preprint2014arXiv

What is the Physical Origin of Strong Lya Emission? I. Demographics of Lya Emitter Structures

We present the results of structure analyses for a large sample of 426 Lya emitters (LAEs) at z~2.2 that are observed with HST/ACS and WFC3-IR by deep extra-galactic legacy surveys. We confirm that the merger fraction and the average ellipticity of LAE's stellar component are 10-30 % and 0.4-0.6, respectively, that are comparable with previous study results. We successfully identify that some LAEs have a spatial offset between Lya and stellar-continuum emission peaks, d_Lya, by ~2.5-4 kpc beyond our statistical errors. To uncover the physical origin of strong Lya emission found in LAEs, we investigate Lya equivalent width (EW) dependences of these three structural parameters, merger fraction, d_Lya, and ellipticity of stellar distribution in the range of EW(Lya)=20-250A. Contrary to expectations, we find that merger fraction does not significantly increase with Lya EW. We reveal an anti-correlation between d_Lya and EW(Lya) by Kolmogorov-Smirnov (KS) test. There is a trend that the LAEs with a large Lya EW have a small ellipticity. This is consistent with the recent theoretical claims that Lya photons can more easily escape from face-on disks having a small ellipticity, due to less inter-stellar gas along the line of sight, although our KS test indicates that this trend is not statistically significant. Our results of Lya-EW dependence generally support the idea that an HI column density is a key quantity determining Lya emissivity.

preprint2014arXiv

What is the physical origin of strong Lya emission? II. Gas Kinematics and Distribution of Lya Emitters

We present a statistical study of velocities of Lya, interstellar (IS) absorption, and nebular lines and gas covering fraction for Lya emitters (LAEs) at z~2. We make a sample of 22 LAEs with a large Lya equivalent width (EW) of > 50A based on our deep Keck/LRIS observations, in conjunction with spectroscopic data from the Subaru/FMOS program and the literature. We estimate the average velocity offset of Lya from a systemic redshift determined with nebular lines to be dv_Lya=234+-9 km s-1. Using a Kolmogorv-Smirnov test, we confirm the previous claim of Hashimoto et al. (2013) that the average dv_Lya of LAEs is smaller than that of LBGs. Our LRIS data successfully identify blue-shifted multiple IS absorption lines in the UV continua of four LAEs on an individual basis. The average velocity offset of IS absorption lines from a systemic redshift is dv_IS=204+-27 km s-1, indicating LAE's gas outflow with a velocity comparable to typical LBGs. Thus, the ratio, R^Lya_ IS = dv_Lya/dv_IS of LAEs, is around unity, suggestive of low impacts on Lya transmission by resonant scattering of neutral hydrogen in the IS medium. We find an anti-correlation between Lya EW and the covering fraction, f_c, estimated from the depth of absorption lines, where f_c is an indicator of average neutral hydrogen column density, N_HI. The results of our study support the idea that N_HI is a key quantity determining Lya emissivity.

preprint2013arXiv

First Spectroscopic Evidence for High Ionization State and Low Oxygen Abundance in Lya Emitters

We present results from Keck/NIRSPEC and Magellan/MMIRS follow-up spectroscopy of Lya emitters (LAEs) at z=2.2 identified in our Subaru narrowband survey. We successfully detect Ha emission from seven LAEs, and perform a detailed analysis of six LAEs free from AGN activity, two out of which, CDFS-3865 and COSMOS-30679, have [OII] and [OIII] line detections. They are the first [OII]-detected LAEs at high-z, and their [OIII]/[OII] ratios and R23-indices provide the first simultaneous determinations of ionization parameter and oxygen abundance for LAEs. CDFS-3865 has a very high ionization parameter (q_{ion}=2.5^{+1.7}_{-0.8}x10^8 cm s^{-1}) and a low oxygen abundance (12+log(O/H)=7.84^{+0.24}_{-0.25}) in contrast with moderate values of other high-z galaxies such as Lyman-break galaxies (LBGs). COSMOS-30679 also possesses a relatively high ionization parameter (q_{ion}=8^{+10}_{-4}x10^7 cm s^{-1}) and a low oxygen abundance (12+log(O/H)=8.18^{+0.28}_{-0.28}). Both LAEs appear to fall below the mass-metallicity relation of z~2 LBGs. Similarly, a low metallicity of 12+log(O/H)<8.4 is independently indicated for typical LAEs from a composite spectrum and the [NII]/Ha index. Such high ionization parameters and low oxygen abundances can be found in local star-forming galaxies, but this extreme local population occupies only ~0.06% of the SDSS spectroscopic galaxy sample with a number density ~100 times smaller than that of LAEs. With their high ionization parameters and low oxygen abundances, LAEs would represent an early stage of galaxy formation dominated by massive stars in compact star-forming regions. High-q_{ion} galaxies like LAEs would produce ionizing photons efficiently with a high escape fraction achieved by density-bounded HII regions, which would significantly contribute to cosmic reionization at z>6.

preprint2013arXiv

Gas Motion Study of Lya Emitters at z~2 Using UV and Optical Spectral Lines

We present the results of Magellan/MMIRS and Keck/NIRSPEC spectroscopy for five Lya emitters (LAEs) at z=2.2 for which high-resolution FUV spectra from Magellan/MagE are available. We detect nebular emission lines including Ha on the individual basis and low-ionization interstellar (LIS) absorption lines in a stacked FUV spectrum, and measure average offset velocities of the Lya line, Delta_v_Lya, and LIS absorption lines, Delta_v_abs, with respect to the systemic velocity defined by the nebular lines. For a sample of eight z~2-3 LAEs without AGN from our study and the literature, we obtain Delta_v_Lya = 175+/-35 km s^{-1}, which is significantly smaller than that of Lyman-break Galaxies (LBGs), Delta_v_Lya=400 km s^{-1}. The stacked FUV spectrum gives Delta_v_abs = -179 +/- 73 km s^{-1}, comparable to that of LBGs. These positive Delta_v_Lya and negative Delta_v_abs suggest that LAEs also have outflows. In contrast to LBGs, however, the LAEs' Delta_v_Lya is as small as |Delta_v_abs|, suggesting low neutral hydrogen column densities. Such a low column density with a small number of resonant scattering may cause the observed strong Lya emission of LAEs. We find an anti-correlation between Lya equivalent width (EW) and Delta_v_Lya in a compilation of LAE and LBG samples. Although its physical origin is not clear, this anti-correlation result appears to challenge the hypothesis that a strong outflow, by means of a reduced number of resonant scattering, produces a large EW. If LAEs at z>6 have similarly small Delta_v_Lya values, constraints on the reionization history derived from the Lya transmissivity may need to be revised.