Recent observations have successfully detected [O III] $88.3\,{\rm \mu m}$ and [C II] $157.6\,{\rm \mu m}$ lines from galaxies in the early Universe with the Atacama Large Millimeter Array (ALMA). Combining cosmological hydrodynamic simulations and radiative transfer calculations, we present relations between the metal line emission and galaxy evolution at $z=6-15$. We find that galaxies during their starburst phases have high [O III] luminosity of $\sim 10^{42}~\rm erg~s^{-1}$. Once supernova feedback quenches star formation, [O III] luminosities rapidly decrease and continue to be zero for $\sim 100\,{\rm Myr}$. The slope of the relation between $\log{(\rm SFR/M_{\odot}~yr^{-1})}$ and $\log{(L_{\rm [O_{III}]}/{\rm L_{\odot}})}$ at $z=6-9$ is 1.03, and 1.43 for $\log{(L_{\rm [C_{II}]}/{\rm L_{\odot}})}$. As gas metallicity increases from sub-solar to solar metallicity by metal enrichment from star formation and feedback, the line luminosity ratio $L_{\rm [O_{III}]} / L_{\rm [C_{II}]}$ decreases from $\sim 10$ to $\sim 1$ because the O/C abundance ratio decreases due to carbon-rich winds from AGB stars. Therefore, we suggest that the combination of [O III] and [C II] lines is a good probe to investigate the relative distribution of ionized and neutral gas in high-$z$ galaxies. In addition, we show that deep [C II] observations with a sensitivity of $\sim 10^{-2}~{\rm mJy~arcsec^{-2}}$ can probe the extended neutral gas disks of high-$z$ galaxies.

We report the discovery of 10-kpc scale [CII] 158um halos surrounding star-forming galaxies in the early Universe. We choose deep ALMA data of 18 galaxies each with a star-formation rate of ~ 10-70 Msun with no signature of AGN whose [CII] lines are individually detected at z=5.153-7.142, and conduct stacking of the [CII] lines and dust-continuum in the uv-visibility plane. The radial profiles of the surface brightnesses show a 10-kpc scale [CII] halo at the 9.2sigma level significantly extended more than the HST stellar continuum data by a factor of ~5 on the exponential-profile basis, as well as the dust continuum. We also compare the radial profiles of [CII] and Lya halos universally found in star-forming galaxies at this epoch, and find that the scale lengths agree within the 1sigma level. While two independent hydrodynamical zoom-in simulations match the dust and stellar continuum properties, the simulations cannot reproduce the extended [CII] line emission. The existence of the extended [CII] halo is the evidence of outflow remnants in the early galaxies and suggest that the outflows may be dominated by cold-mode outflows expelling the neutral gas.

Recent observations have successfully detected UV or infrared flux from galaxies at the epoch of reionization. However, the origin of their radiative properties has not been fully understood yet. Combining cosmological hydrodynamic simulations and radiative transfer calculations, we present theoretical predictions of multi-wavelength radiative properties of the first galaxies at z=6-15. We find that most of the gas and dust are ejected from star-forming regions due to supernova (SN) feedback, which allows UV photons to escape. We show that the peak of SED rapidly shifts between UV and infrared wavelengths on a timescale of 100 Myr due to intermittent star formation and feedback. When dusty gas covers the star-forming regions, the galaxies become bright in the observed-frame sub-millimeter wavelengths. In addition, we find that the escape fraction of ionizing photons also changes between 1-40% at z>10. The mass fraction of HII region changes with the star formation history, resulting in the fluctuations of metal lines and Lyman-alpha line luminosities. In the starbursting phase of galaxies with the halo mass $\sim 10^{11}\,{\rm M_{\odot}}$ ($10^{12}\,{\rm M_{\odot}}$), the simulated galaxy has $L_{\rm [O\,III]} \sim 10^{42}\, (10^{43})\, {\rm erg\,s^{-1}}$, which is consistent with the observed star-forming galaxies at z>7. Our simulations suggest that deep [C II] observation with ALMA can trace the distribution of neutral gas extending over $\sim 20$ physical kpc. We also find that the luminosity ratio $L_{\rm [O\,III]}/L_{\rm [C\,II]}$ decreases with bolometric luminosity due to metal enrichment. Our simulations show that the combination of multi-wavelength observations by ALMA and JWST will be able to reveal the multi-phase ISM structure and the transition from starbursting to outflowing phases of high-z galaxies.

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

Recent observations have successfully detected UV-bright and infrared-bright galaxies in the epoch of reionization. However, the origin of their radiative properties has not been understood yet. Combining cosmological hydrodynamic simulations and radiative transfer calculations, we present predictions of multi-wavelength radiative properties of the first galaxies at $z\sim 6-15$. Using zoom-in initial conditions, we investigate three massive galaxies and their satellites in different environment and halo masses: $M_{\rm h} = 2.4\times 10^{10}{\rm M_{\odot}}$ (Halo-10), $1.6\times 10^{11}{\rm M_{\odot}}$ (Halo-11) and $0.7\times 10^{12}{\rm M_{\odot}}$ (Halo-12) at $z=6$. We find that most of gas and dust are ejected from star-forming regions by supernova feedback, which allows UV photons to escape. We show that the peak of the spectral energy distribution (SED) rapidly changes between UV and infrared wavelengths on a time-scale of $\sim$\,100\,Myrs due to intermittent star formation and feedback, and the escape fraction of UV photons fluctuates in the range of $0.2-0.8$ at $z<10$ with a time-averaged value of 0.3. When dusty gas covers the star-forming regions, the galaxies become bright in the observed-frame sub-millimeter wavelengths. We predict the detectability of high-$z$ galaxies with the Atacama Large Millimeter Array (ALMA). For a sensitivity limit of $0.1\,{\rm mJy}$ at $850\,{\rm \mu m}$, the detection probability of galaxies in halos $M_{\rm h} \gtrsim 10^{11}\,{\rm M_{\odot}}$ at $z\lesssim 7$ exceeds fifty per cent. We argue that supernova feedback can produce the observed diversity of SEDs for high-$z$ galaxies.

Star formation in high-redshift dwarf galaxies is a key to understand early galaxy evolution in the early Universe. Using the three-dimensional hydrodynamics code GIZMO, we study the formation mechanism of cold, high-density gas clouds in interacting dwarf galaxies with halo masses of $\sim 3 \times 10^{7}~M_{\odot}$, which are likely to be the formation sites of early star clusters. Our simulations can resolve both the structure of interstellar medium on small scales of $\lesssim 0.1$ pc and the galactic disk simultaneously. We find that the cold gas clouds form in the post-shock region via thermal instability due to metal-line cooling, when the cooling time is shorter than the galactic dynamical time. The mass function of cold clouds shows almost a power-law initially with an upper limit of thermally unstable scale. We find that some clouds merge into more massive ones with $\gtrsim 10^{4}~M_{\odot}$ within $\sim 2~{\rm Myr}$. Only the massive cold clouds with $\gtrsim 10^{3}~M_{\odot}$ can keep collapsing due to gravitational instability, resulting in the formation of star clusters. In addition, we investigate the dependence of cloud mass function on metallicity and ${\rm H_{2}}$ abundance, and show that the cases with low metallicities ($\lesssim 10^{-2}~Z_{\odot}$) or high ${\rm H_{2}}$ abundance ($\gtrsim 10^{-3}$) cannot form massive cold clouds with $\gtrsim 10^{3}~M_{\odot}$.