Researcher profile

Takuya Hashimoto

Takuya Hashimoto contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

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

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

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

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

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

Building this graph slice

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

Published work

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

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.