We study the future prospects of the 21cm forest observations on the axion-like dark matter when the spontaneous breaking of the global Peccei-Quinn (PQ) symmetry occurs after the inflation. The large isocurvature perturbations of order unity sourced from axion-like particles can result in the enhancement of minihalo formation, and the subsequent hierarchical structure formation can affect the minihalo abundance whose masses can exceed ${\cal O}(10^4) M_{\odot}$ relevant for the 21cm forest observations. We show that the 21cm forest observations are capable of probing the axion-like particle mass in the range $10^{-18}\lesssim m_a \lesssim 10^{-12}$ eV for the temperature independent axion mass. For the temperature dependent axion mass, the zero temperature axion mass scale for which the 21cm forest measurements can be affected is extended further to as big as of order $10^{-6}$ eV.

The radio observation of 21cm-line signal from the Epoch of Reionization (EoR) enables us to explore the evolution of galaxies and intergalactic medium in the early universe. However, the detection and imaging of 21cm-line signal are tough due to the foreground and instrumental systematics. In order to overcome these obstacles, as a new approach, we propose to take a cross correlation between observed 21cm-line data and 21cm-line images generated from the distribution of the Lyman-$\alpha$ emitters (LAEs) through machine learning. In order to create 21cm-line maps from LAE distribution, we apply conditional Generative Adversarial Network (cGAN) trained with the results of our numerical simulations. We find that 21cm-line brightness temperature maps and the neutral fraction maps can be well reproduced at large scales. Furthermore, we show that the cross correlation is detectable at $k < 0.2~{\rm Mpc}^{-1}$ by combing 400 hours of MWA Phase II observation and LAE deep survey of the Subaru Hyper Suprime Camera. Our new approach of cross correlation with image construction using the cGAN can not only boost the detectability of EoR 21cm-line signal but also allow us to estimate the 21cm-line auto-power spectrum.

The bubble size distribution of ionized hydrogen regions, which could be derived from the tomographic imaging data of the redshifted 21~cm signal from the epoch of reionization, probes the information about the morphology of \HII\ bubbles during the reionization. However, 21~cm imaging is observationally very challenging even for the upcoming large radio interferometers. Given that these interferometers promise to measure the 21~cm power spectrum accurately, we propose a new method, which is based on the artificial neural networks (ANN), to reconstruct the \HII\ bubble size distribution from the 21~cm power spectrum. We demonstrate that the reconstruction from the 21~cm power spectrum can be almost as accurate as directly measured from the imaging data with the fractional error $\lesssim 10\%$, even with thermal noise at the sensitivity level of the Square Kilometre Array. Nevertheless, systematic errors might arise from approximations made in reionization simulations used for training the ANN. This paper, as the first in a series, exemplifies the possibility of recovering from the 21~cm power spectrum with ANN additional statistics of cosmic reionization that could not otherwise be inferred from the power spectrum analysis directly in the conventional methods.

The ultra-light scalar fields can arise ubiquitously, for instance, as a result of the spontaneous breaking of an approximate symmetry such as the axion and more generally the axion-like particles. In addition to the particle physics motivations, these particles can also play a major role in cosmology by contributing to dark matter abundance and affecting the structure formation at sub-Mpc scales. In this paper, we propose to use the 21cm forest observations to probe the nature of ultra-light dark matter. The 21cm forest can probe much smaller scales than the Lyman-$\alpha$ forest, that is, $k\gtrsim 10\mathrm{Mpc}^{-1}$. We explore the range of the ultra-light dark matter mass $m_{u}$ and $f_u$, the fraction of ultra-light dark matter with respect to the total matter, which can be probed by the 21cm forest. We find that 21cm forest can potentially put the dark matter mass lower bound $m_u \gtrsim 10^{-18}$ eV for $f_u=1$, which is 3 orders of magnitude bigger mass scale than those probed by the current Lyman-$\alpha$ forest observations.While the effects of the ultra-light particles on the structure formation become smaller when the dominant component of dark matter is composed of the conventional cold dark matter, we find that the 21cm forest is still powerful enough to probe the sub-component ultra-light dark matter mass up to the order of $10^{-19}$ eV. The Fisher matrix analysis shows that $(m_u,f_u)\sim (10^{-20}\mathrm{eV}, 0.3)$ is the most optimal parameter set which the 21cm forest can probe with the minimal errors for a sub-component ultra-light dark matter scenario.

Understanding the detailed process of cosmic reionization is one of the remaining problems in astrophysics and cosmology. Here we construct a model of cosmic reionization that includes contribution from high-$z$ galaxies and active galactic nuclei (AGNs), and calculate the reionization and thermal histories with the model. To keep the model to be general and realistic, we vary the escape fraction of ionizing photons, $f_{\rm esc}$, and the faint-end slope of the AGN luminosity function at high redshifts, $\alpha_{\rm hz}$, within the constraints from the observed cosmic star formation history and the observed bright-end UV luminosity functions at $z \leq 6$. Additionally, we model the spectral energy distribution (SED) of AGNs which depends on the Eddington ratio and the black-hole mass. By comparing the computed reionization histories with the observed \ion{H}{i} fractions and the optical depth for Thomson scattering from Planck, we find that $\alpha_{\rm hz} > -1.5$ and $f_{\rm esc} < 0.15$ are favored when we employ the bright-end luminosity function obtained by Giallongo et al. Our result suggests that an AGN-dominated model with the abundance of faint-AGNs being as large as the estimate by Giallongo et al. is allowed only if the contribution from high-z galaxies is almost negligible, while a galaxy-dominant model is also allowed. We also find that the shape of SED has a significant impact on the thermal history. Therefore it is expected that measurements of the thermal state of the IGM provide useful information on properties of ionizing sources.

The 21 cm signal from the Epoch of Reionization should be observed within the next decade. While a simple statistical detection is expected with SKA pathfinders, the SKA will hopefully produce a full 3D mapping of the signal. To extract from the observed data constraints on the parameters describing the underlying astrophysical processes, inversion methods must be developed. For example, the Markov Chain Monte Carlo method has been successfully applied. Here we test another possible inversion method: artificial neural networks (ANN). We produce a training set which consists of 70 individual sample. Each sample is made of the 21 cm power spectrum at different redshifts produced with the 21cmFast code plus the value of three parameters used in the semi-numerical simulations that describe astrophysical processes. Using this set we train the network to minimize the error between the parameter values it produces as an output and the true values. We explore the impact of the architecture of the network on the quality of the training. Then we test the trained network on the new set of 54 test samples with different values of the parameters. We find that the quality of the parameter reconstruction depends on the sensitivity of the power spectrum to the different parameters at a given redshift, that including thermal noise and sample variance decreases the quality of the reconstruction and that using the power spectrum at several redshifts as an input to the ANN improves the quality of the reconstruction. We conclude that ANNs are a viable inversion method whose main strength is that they require a sparse exploration of the parameter space and thus should be usable with full numerical simulations.

We perform a Fisher analysis to estimate the expected constraints on the Epoch of Reionization (EoR) model parameters (i.e., minimum virial temperature, the ionizing efficiency and the mean free path of ionizing photons) taking into account the thermal noise of the exisitting telescopes, MWA and LOFAR. We consider how the inclusion of the 21cm bispectrum improves the constraints compared to using the power spectrum alone. Assuming that we perfectly remove the foreground, we found that the bispectrum, which is calculated by the 21cmFAST code, can constrain the EoR model parameters more tightly than the power spectrum since the bispectrum is more sensitive to the EoR model parameters than the power spectrum. We also found that degeneracy between the EoR model parameters can be broken by combining the bispectrum with the power spectrum.

In the past several decades, the standard cosmological model has been established and its parameters have been measured to a high precision, while there are still many of the fundamental questions in cosmology; such as the physics in the very early Universe, the origin of the cosmic acceleration and the nature of the dark matter. The future world's largest radio telescope, Square Kilometre Array (SKA), will be able to open the new frontier of cosmology and will be one of the most powerful tools for cosmology in the next decade. The cosmological surveys conducted by the SKA would have the potential not only to answer these fundamental questions but also deliver the precision cosmology. In this article we briefly review the role of the SKA from the view point of the modern cosmology. The cosmology science led by the SKA-Japan Consortium (SKA-JP) Cosmology Science Working Group is also discussed.

The brightness temperature of the redshifted 21cm line brings rich information on the Inter Galactic Medium (IGM) from the Cosmic Dawn and Epoch of Reionization (EoR). While the power spectrum is a useful tool to statistically investigate the 21cm signal, the 21cm brightness temperature field is highly non-Gaussian, and the power spectrum is inadequate to characterize the non-Gaussianity. The Minkowski Functionals (MFs) are promising tools to extract non-gaussian features of the 21cm signal and give topological information such as morphology of ionized bubbles. In this work, we study the 21cm line signal in detail with MFs. To promote understanding of basic features of the 21cm signal, we calculate the MFs of not only the hydrogen neutral fraction but the matter density and spin temperature, which contribute to the brightness temperature fluctuations. We find that the structure of the brightness temperature depends mainly on the ionized fraction and the spin temperature at late and early stages of the EoR, respectively. Further, we investigate the redshift evolution of the MFs at $7 < z < 20$. Then, we find that, after the onset of reionization, the MFs reflect mainly the ionized bubble property. In addition, the MFs are sensitive to model parameters which are related to the topology of ionized bubbles and we consider the possibility of constraining the parameters by the future 21cm signal observations.

Redshifted 21cm-line signal from neutral hydrogens in the intergalactic medium (IGM) gives a direct probe of the epoch of reionization (EoR). In this paper, we investigate the potential of the variance and skewness of the probability distribution function of the 21cm brightness temperature for constraining EoR models. These statistical quantities are simple, easy to calculate from the observed visibility and thus suitable for the early exploration of the EoR with ongoing telescopes such as the Murchison Widefield Array (MWA) and LOw Frequency ARray (LOFAR). We show, by performing Fisher analysis, that the variance and skewness at $z=7-9$ are complementary to each other to constrain the EoR model parameters such as the minimum virial temperature of halos which host luminous objects, ionizing efficiency and mean free path of ionizing photons in the IGM. Quantitatively, the constraining power highly depends on the quality of the foreground subtraction and calibration. We give a best case estimate of the constraints on the parameters, neglecting the systematics other than the thermal noise.

Redshifted 21cm signal is a promising tool to investigate the state of intergalactic medium (IGM) in the Cosmic Dawn (CD) and Epoch of Reionization(EoR). In our previous work (Shimabukuro et al 2015), we studied the variance and skewness to give a clear interpretation of 21cm power spectrum and found that skewness is a good indicator of the epoch when X-ray heating becomes effective. Thus, the non-Gaussian feature of the spatial distribution of the 21cm signal is expected to be useful to investigate the astrophysical effects in the CD and EoR. In this paper, in order to investigate such a non-Gaussian feature in more detail, we focus on the bispectrum of the 21cm signal. It is expected that the 21cm brightness temperature bispectrum is produced by non-gaussianity due to the various astrophysical effects such as Wouthysen-Field (WF) effect, X-ray heating and reionization. We study the various properties of 21cm bispectrum such as scale dependence, shape dependence and redshift evolution. And also we study the contribution from each component of 21cm bispectrum. We find that the contribution from each component has characteristic scale-dependent feature, and it is expected that we could obtain more detailed information on the IGM in the CD and EoR by using the 21cm bispectrum in the future experiments, combined with the power spectrum and skewness.

Cosmic reionization is known to be a major phase transition of the gas in the Universe. Since astronomical objects formed in the early Universe, such as the first stars, galaxies and black holes, are expected to have caused cosmic reionization, the formation history and properties of such objects are closely related to the reionization process. In spite of the importance of exploring reionization, our understandings regarding reionization is not sufficient yet. Square Kilometre Array (SKA) is a next-generation large telescope that will be operated in the next decade. Although several programs of next-generation telescopes are currently scheduled, the SKA will be the unique telescope with a potential to directly observe neutral hydrogen up to z~30, and provide us with valuable information on the Cosmic Dawn (CD) and the Epoch of Reionization (EoR). The early science with the SKA will start in a few years; it is thus the time for us to elaborate a strategy for CD/EoR Science with the SKA. The purpose of this document is to introduce Japanese scientific interests in the SKA project and to report results of our investigation.

The redshifted 21cm line signal from neutral hydrogens is a promising tool to probe the cosmic dawn and the epoch of reionization (EoR). Ongoing and future low-frequency radio experiments are expected to detect its fluctuations, especially through the power spectrum. In this paper, we give a physical interpretation of the time evolution of the power spectrum of the 21cm brightness temperature fluctuations, which can be decomposed into dark matter density, spin temperature and neutral fraction of hydrogen fluctuations. From the one-point statistics of the fluctuations, such as variance and skewness, we find that the peaks and dips in the time evolution are deeply related to X-ray heating of the intergalactic gas, which controls the spin temperature. We suggest the skewness of the brightness temperature distribution is a key observable to identify the onset of X-ray heating.

The 21cm line brightness temperature brings rich information about Epoch of Reionizaton (EoR) and high-$z$ universe (Cosmic Dawn and Dark Age). While the power spectrum is a useful tool to investigate the EoR signal statistically, higher-order statistics such as bispectrum are also valuable because the EoR signal is expected to be highly non-Gaussian. In this paper, we develop a formalism to calculate the bispectrum contributed from the thermal noise taking array configularion of telescopes into account, by extending a formalism for the power spectrum \cite{2006ApJ...653..815M}. We apply our formalism to the ongoing and future telescopes such as expanded Murchison Widefield Array (MWA), LOw Frequency ARray (LOFAR) and Square Kilometre Array (SKA). We find that expanded MWA does not have enough sensitivity to detect the bispectrum signal. On the other hand, LOFAR has better sensitivity and will be able to detect the peaks of the bispectrum as a function of redshift at large scales with comoving wavenumber $k \lesssim 0.03 {\rm Mpc}^{-1}$. The SKA has enough sensitivity to detect the bispectrum for much smaller scales $k \lesssim 0.3 {\rm Mpc}^{-1}$ and redshift $z \lesssim 20$

Although the cosmological paradigm based on cold dark matter and adiabatic, nearly scale-invariant primordial fluctuations is consistent with a wide variety of existing observations, it has yet to be sufficiently tested on scales smaller than those of massive galaxies, and various alternatives have been proposed that differ significantly in the consequent small-scale power spectrum (SSPS) of large-scale structure. Here we show that a powerful probe of the SSPS at $k\gtrsim 10$ Mpc$^{-1}$ can be provided by the 21 cm forest, that is, systems of narrow absorption lines due to intervening, cold neutral hydrogen in the spectra of high-redshift background radio sources in the cosmic reionization epoch. Such features are expected to be caused predominantly by collapsed gas in starless minihalos, whose mass function can be very sensitive to the SSPS. As specific examples, we consider the effects of neutrino mass, running spectral index (RSI) and warm dark matter (WDM) on the SSPS, and evaluate the expected distribution in optical depth of 21 cm absorbers out to different redshifts. Within the current constraints on quantities such as the sum of neutrino masses $\sum m_\nu$, running of the primordial spectral index $d n_s/d \ln k$ and WDM particle mass $m_{\rm WDM}$, the statistics of the 21 cm forest manifest observationally significant differences that become larger at higher redshifts. In particular, it may be possible to probe the range of $m_{\rm WDM} \gtrsim 10$ keV that may otherwise be inaccessible. Future observations of the 21 cm forest by the Square Kilometer Array may offer a unique and valuable probe of the SSPS, as long as radio sources such as quasars or Population III gamma-ray bursts with sufficient brightness and number exist at redshifts of $z \gtrsim$ 10 - 20, and the astrophysical effects of reionization and heating can be discriminated.