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.

Line intensity mapping (LIM) is an emerging observational method to study the large-scale structure of the Universe and its evolution. LIM does not resolve individual sources but probes the fluctuations of integrated line emissions. A serious limitation with LIM is that contributions of different emission lines from sources at different redshifts are all confused at an observed wavelength. We propose a deep learning application to solve this problem. We use conditional generative adversarial networks to extract designated information from LIM. We consider a simple case with two populations of emission line galaxies; H$\rm\alpha$ emitting galaxies at $z = 1.3$ are confused with [OIII] emitters at $z = 2.0$ in a single observed waveband at 1.5 $\rm\mu$m. Our networks trained with 30,000 mock observation maps are able to extract the total intensity and the spatial distribution of H$\rm\alpha$ emitting galaxies at $z = 1.3$. The intensity peaks are successfully located with 74% precision. The precision increases to 91% when we combine the results of 5 networks. The mean intensity and the power spectrum are reconstructed with an accuracy of $\sim$10%. The extracted galaxy distributions at a wider range of redshift can be used for studies on cosmology and on galaxy formation and evolution.

We study statistics of the 21-cm signal from the epoch of reionization. We propose to use [OIII] line emitting galaxies to cross-correlate with the 21-cm signal from $z = 7 - 10$. To this aim, we employ simulations of reionization obtained post-processing the high-resolution cosmological hydrodynamic simulation Massive Black-II with the three-dimensional (3D) radiative transfer code CRASH to follow the propagation of ionizing photons from a variety of sources. We show that, during the early phases of reionization, the 21-cm signal is positively correlated with the spatial distribution of the [OIII] emitters on large scales $(k < 1h~\rm cMpc^{-1})$. This positive correlation is generated by the temperature-galaxy correlation and it is a few times larger than when we assume that the heating is saturated. As the reionized regions expand, the correlation changes its sign to negative from $z = 10$ to 8. The signals at this epoch can be detected by combining the Square Kilometre Array (SKA) and a wide-field [OIII] emitter survey. We also calculate the cross-power spectrum with a 3D [OIII] intensity field, aiming at exploiting future intensity mapping observations. We conclude that high-redshift [OIII] line emitters can be used to probe the reionization process when the intergalactic medium is largely neutral.

Recent observations with the Atacama Large Millimeter/submillimeter Array (ALMA) detected far-infrared emission lines such as the [OIII] 88 \mu m line from galaxies at $z \sim 7 - 9$. Far-infrared lines can be used to probe the structure and kinematics of such high-redshift galaxies as well as to accurately determine their spectroscopic redshifts. We use a cosmological simulation of galaxy formation to study the physical properties of [OIII] 88 \mu m emitters. In a comoving volume of 50 $h^{-1}$ Mpc on a side, we locate 34 galaxies with stellar masses greater than $10^8\ {\rm M_\odot}$ at $z = 9$, and more than 270 such galaxies at $z = 7$. We calculate the [OIII] 88 \mu m luminosities ($L_{\rm OIII}$) by combining a physical model of HII regions with emission line calculations using the photoionization code CLOUDY. We show that the resulting $L_{\rm OIII}$, for a given star formation rate, is slightly higher than predicted from the empirical relation for local galaxies, and is consistent with recent observations of galaxies at redshifts 7 - 9. Bright [OIII] emitters with $L_{\rm OIII} > 10^8 {\rm L_\odot}$ have stellar masses greater than $10^9\ {\rm M_\odot}$, star formation rates higher than $3\ {\rm M_\odot\ yr}^{-1}$, and the typical metallicity is $\sim 0.1\ {\rm Z_\odot}$. The galaxies are hosted by dark matter halos with masses greater than $10^{10.5}\ {\rm M_\odot}$. Massive galaxies show characteristic structure where the [OIII] emitting gas largely overlaps with young stars, but the emission peak is separated from the main stellar population, suggesting the stochastic and localized nature of star formation in the first galaxies. We propose to use the [OIII] 5007 \AA\ line, to be detected by James Webb Space Telescope (JWST), to study the properties of galaxies whose [OIII] 88 \mu m line emission has been already detected with ALMA.

We have analyzed the OCS, H$_2$CS, CH$_3$OH, and HCOOCH$_3$ data observed toward the low-mass protostar IRAS 16293--2422 Source B at a sub-arcsecond resolution with ALMA. A clear chemical differentiation is seen in their distributions; OCS and H$_2$CS are extended with a slight rotation signature, while CH$_3$OH and HCOOCH$_3$ are concentrated near the protostar. Such a chemical change in the vicinity of the protostar is similar to the companion (Source A) case. The extended component is interpreted by the infalling-rotating envelope model with a nearly face-on configuration. The radius of the centrifugal barrier of the infalling-rotating envelope is roughly evaluated to be ($30-50$) au. The observed lines show the inverse P-Cygni profile, indicating the infall motion with in a few 10 au from the protostar. The nearly pole-on geometry of the outflow lobes is inferred from the SiO distribution, and thus, the infalling and outflowing motions should coexist along the line-of-sight to the protostar. This implies that the infalling gas is localized near the protostar and the current launching points of the outflow have an offset from the protostar. A possible mechanism for this configuration is discussed.