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.

We argue that the global signal of neutral hydrogen 21cm line can be a powerful probe of primordial power spectrum on small scales. Since the amplitude of small scale primordial fluctuations is important to determine the early structure formation and the timing when the sources of Lyman ${\alpha}$ photons are produced, they in turn affect the neutral hydrogen 21cm line signal. We show that the information of the position of the absorption trough can severely constrain the small scale amplitude of primordial fluctuations once astrophysical parameters relevant to the 21cm line signal are fixed. We also discuss how the uncertainties of astrophysical parameters affect the constraints.

Pulsar search with time-domain observation is very computationally expensive and data volume will be enormous with the next generation telescopes such as the Square Kilometre Array. We apply artificial neural networks (ANNs), a machine learning method, for efficient selection of pulsar candidates from radio continuum surveys, which are much cheaper than time-domain observation. With observed quantities such as radio fluxes, sky position and compactness as inputs, our ANNs output the "score" that indicates the degree of likeliness of an object to be a pulsar. We demonstrate ANNs based on existing survey data by the TIFR GMRT Sky Survey (TGSS) and the NRAO VLA Sky Survey (NVSS) and test their performance. Precision, which is the ratio of the number of pulsars classified correctly as pulsars to that of any objects classified as pulsars, is about 96$\%$. Finally, we apply the trained ANNs to unidentified radio sources and our fiducial ANN with five inputs (the galactic longitude and latitude, the TGSS and NVSS fluxes and compactness) generates 2,436 pulsar candidates from 456,866 unidentified radio sources. These candidates need to be confirmed if they are truly pulsars by time-domain observations. More information such as polarization will narrow the candidates down further.

We compute the spherically-averaged power spectrum from four seasons of data obtained for the Epoch of Reionisation (EoR) project observed with the Murchison Widefield Array (MWA). We measure the EoR power spectrum over $k= 0.07-3.0~h$Mpc$^{-1}$ at redshifts $z=6.5-8.7$. The largest aggregation of 110 hours on EoR0 high-band (3,340 observations), yields a lowest measurement of (43~mK)$^2$ = 1.8$\times$10$^3$ mK$^2$ at $k$=0.14~$h$Mpc$^{-1}$ and $z=6.5$ (2$\sigma$ thermal noise plus sample variance). Using the Real-Time System to calibrate and the CHIPS pipeline to estimate power spectra, we select the best observations from the central five pointings within the 2013--2016 observing seasons, observing three independent fields and in two frequency bands. This yields 13,591 2-minute snapshots (453 hours), based on a quality assurance metric that measures ionospheric activity. We perform another cut to remove poorly-calibrated data, based on power in the foreground-dominated and EoR-dominated regions of the two-dimensional power spectrum, reducing the set to 12,569 observations (419 hours). These data are processed in groups of 20 observations, to retain the capacity to identify poor data, and used to analyse the evolution and structure of the data over field, frequency, and data quality. We subsequently choose the cleanest 8,935 observations (298 hours of data) to form integrated power spectra over the different fields, pointings and redshift ranges.

The compact configuration of Phase II of the Murchison Widefield Array (MWA) consists of both a redundant subarray and pseudo-random baselines, offering unique opportunities to perform sky-model and redundant interferometric calibration. The highly redundant hexagonal cores give improved power spectrum sensitivity. In this paper, we present the analysis of nearly 40 hours of data targeting one of the MWA's EoR fields observed in 2016. We use both improved analysis techniques presented in Barry et al. (2019) as well as several additional techniques developed for this work, including data quality control methods and interferometric calibration approaches. We show the EoR power spectrum limits at redshift 6.5, 6.8 and 7.1 based on our deep analysis on this 40-hour data set. These limits span a range in $k$ space of $0.18$ $h$ $\mathrm{Mpc^{-1}}$ $<k<1.6$ $h$ $\mathrm{Mpc^{-1}}$, with a lowest measurement of $\Delta^2\leqslant2.39\times 10^3$ $\mathrm{mK}^2$ at $k=0.59$ $h$ $\mathrm{Mpc^{-1}}$ and $z=6.5$.

Measurements of 21 cm Epoch of Reionization (EoR) structure are subject to systematics originating from both the analysis and the observation conditions. Using 2013 data from the Murchison Widefield Array (MWA), we show the importance of mitigating both sources of contamination. A direct comparison between results from Beardsley et al. 2016 and our updated analysis demonstrates new precision techniques, lowering analysis systematics by a factor of 2.8 in power. We then further lower systematics by excising observations contaminated by ultra-faint RFI, reducing by an additional factor of 3.8 in power for the zenith pointing. With this enhanced analysis precision and newly developed RFI mitigation, we calculate a noise-dominated upper limit on the EoR structure of $\Delta^2 \leq 3.9 \times 10^3$ mK$^2$ at $k=0.20$ $\textit{h}$ Mpc$^{-1}$ and $z=7$ using 21 hr of data, improving previous MWA limits by almost an order of magnitude.

The Murchison Widefield Array (MWA) is an open access telescope dedicated to studying the low frequency (80$-$300 MHz) southern sky. Since beginning operations in mid 2013, the MWA has opened a new observational window in the southern hemisphere enabling many science areas. The driving science objectives of the original design were to observe 21\,cm radiation from the Epoch of Reionisation (EoR), explore the radio time domain, perform Galactic and extragalactic surveys, and monitor solar, heliospheric, and ionospheric phenomena. All together 60$+$ programs recorded 20,000 hours producing 146 papers to date. In 2016 the telescope underwent a major upgrade resulting in alternating compact and extended configurations. Other upgrades, including digital back-ends and a rapid-response triggering system, have been developed since the original array was commissioned. In this paper we review the major results from the prior operation of the MWA, and then discuss the new science paths enabled by the improved capabilities. We group these science opportunities by the four original science themes, but also include ideas for directions outside these categories.

We apply two methods to estimate the 21~cm bispectrum from data taken within the Epoch of Reionisation (EoR) project of the Murchison Widefield Array (MWA). Using data acquired with the Phase II compact array allows a direct bispectrum estimate to be undertaken on the multiple redundantly-spaced triangles of antenna tiles, as well as an estimate based on data gridded to the $uv$-plane. The direct and gridded bispectrum estimators are applied to 21 hours of high-band (167--197~MHz; $z$=6.2--7.5) data from the 2016 and 2017 observing seasons. Analytic predictions for the bispectrum bias and variance for point source foregrounds are derived. We compare the output of these approaches, the foreground contribution to the signal, and future prospects for measuring the bispectra with redundant and non-redundant arrays. We find that some triangle configurations yield bispectrum estimates that are consistent with the expected noise level after 10 hours, while equilateral configurations are strongly foreground-dominated. Careful choice of triangle configurations may be made to reduce foreground bias that hinders power spectrum estimators, and the 21~cm bispectrum may be accessible in less time than the 21~cm power spectrum for some wave modes, with detections in hundreds of hours.

We explore methods for robust estimation of the 21 cm signal from the Epoch of Reionisation (EoR). A Kernel Density Estimator (KDE) is introduced for measuring the spatial temperature fluctuation power spectrum from the EoR. The KDE estimates the underlying probability distribution function of fluctuations as a function of spatial scale, and contains different systematic biases and errors to the typical approach to estimating the fluctuation power spectrum. Extraction of histograms of visibilities allows moments analysis to be used to discriminate foregrounds from 21 cm signal and thermal noise. We use the information available in the histograms, along with the statistical dis-similarity of foregrounds from two independent observing fields, to robustly separate foregrounds from cosmological signal, while making no assumptions about the Gaussianity of the signal. Using two independent observing fields to robustly discriminate signal from foregrounds is crucial for the analysis presented in this paper. We apply the techniques to 13 hours of Murchison Widefield Array (MWA) EoR data over two observing fields. We compare the output to that obtained with a comparative power spectrum estimation method, and demonstrate the reduced foreground contamination using this approach. Using the second moment obtained directly from the KDE distribution functions yields a factor of 2-3 improvement in power for k < 0.3hMpc^{-1} compared with a matched delay space power estimator, while weighting data by additional statistics does not offer significant improvement beyond that available for thermal noise-only weights.

Neutral hydrogen (HI) intensity mapping is capable of measuring redshift evolution of HI density parameter $\Omega_{\rm HI}$, which is an important parameter to understand structure formation in the post-reionization epoch. Future HI observation with Square Kilometre Array (SKA) can significantly improve constraints on the parameter. However, the observation of HI suffers from the contamination from extremely bright foreground emissions, and it is necessary to consider a signal validation method complementary to the measurement of HI auto power spectrum. In this work, we propose to take a cross correlation between 21cm line intensity map and convergence map reconstructed from Planck observation and forecast a constraint on the $\Omega_{\rm HI}$. We find that the SKA-Mid operated as single-dish mode has a sufficient capability to detect the cross correlation and constrain $\Omega_{\rm HI}$ with $10-20\%$ precision at wide range of redshifts $0.5<z<3$, when combined with the Planck constraints on the cosmological parameters.

Observation of the 21cm line signal from neutral hydrogen during the Epoch of Reionization is challenging due to extremely bright Galactic and extragalactic foregrounds and complicated instrumental calibration. A reasonable approach for mitigating these problems is the cross correlation with other observables. In this work, we present the first results of the cross power spectrum (CPS) between radio images observed by the Murchison Widefield Array and the cosmic microwave background (CMB), measured by the Planck experiment. We study the systematics due to the ionospheric activity, the dependence of CPS on group of pointings, and frequency. The resulting CPS is consistent with zero because the error is dominated by the foregrounds in the 21cm observation. Additionally, the variance of the signal indicates the presence of unexpected systematics error at small scales. Furthermore, we reduce the error by one order of magnitude with application of a foreground removal using a polynomial fitting method. Based on the results, we find that the detection of the 21cm-CMB CPS with the MWA Phase I requires more than 99.95% of the foreground signal removed, 2000 hours of deep observation and 50% of the sky fraction coverage.

We investigate the impact of the recent observation of the 21cm global signal by EDGES on primordial power spectrum, particularly focusing on the running parameters $\alpha_s$ and $ \beta_s$ which characterize the detailed scale dependence of the primordial spectrum. When primordial power spectrum is enhanced/suppressed on small scales, the structure formation proceeds faster/slower and changes the abundance of small size halos, which affects the sources of Lyman $\alpha$ radiation at high redshifts to alter the position of the absorption line. Recent observation of EDGES detected the 21cm absorption line at $z=17.2$ and this result also indicates that the brightness temperature is consistent with zero for $z \lesssim 14$ and $z \gtrsim 22$, which can exclude a scenario giving the absorption line at such redshifts. We argue that the bound on the running parameters can be obtained by requiring that the absorption line should not exceed observational bounds at such redshift ranges and found that some parameter space of $\alpha_s$ and $\beta_s$ allowed by Planck may be disfavored for some values of astrophysical parameters.

Cross-correlation between the redshifted 21 cm signal and Lyman-{\alpha} emitters (LAEs) is powerful tool to probe the Epoch of Reionization (EoR). Although the cross-power spectrum (PS) has an advantage of not correlating with foregrounds much brighter than the 21 cm signal, the galactic and extra-galactic foregrounds prevent detection since they contribute to the variance of the cross PS. Therefore, strategies for mitigating foregrounds are required. In this work, we study the impact of foreground avoidance on the measurement of the 21 cm-LAE cross-correlation. We then simulate the 21 cm observation as observed by the Murchison Widefield Array (MWA). The point source foreground is modelled from the GaLactic and Extragalactic All-sky Murchison Widefield Array (GLEAM) survey catalogue, and the diffuse foreground is evaluated using a parametric model. For LAE observations, we assume a large survey of the Subaru Hyper Supreme-Cam (HSC), with spectroscopic observations of the Prime Focus Spectrograph (PFS). To predict the 21 cm signal, we employ a numerical simulation combining post processed radiative transfer and radiation hydrodynamics. Using these models, the signal-to-noise ratio of 2D PS shows the foreground contamination dominates the error of cross-PS even in the so-called `EoR window'. We find that at least 99% of the point source foreground and 80% of the galactic diffuse foreground must be removed to measure the EoR signal at large scales $k<0.5 h \rm Mpc^{-1}$. Additionally, a sensitivity 80 times larger than that of the MWA operating with 128 tiles and 99% of the point source foreground removal are required for a detection at small scales.

Detection of the redshifted 21cm-line signal from neutral hydrogen in the intergalactic medium (IGM) during the Epoch of Reionization (EoR) is complicated by intense foregrounds such as galactic synchrotron and extragalactic radio galaxies. The 21cm-Lyman-$\alpha$ emitter(LAE) cross-correlation is one of the tools available to reduce the foreground effects because the foreground emission from such radio sources is statistically independent of LAE distribution. LAE surveys during the EoR at redshifts $z=6.6$ and $7.3$ are ongoing by the Subaru Hyper Suprime-Cam (HSC). Additionally, Prime Focus Spectrograph (PFS) will provide precise redshift information of the LAEs discovered by the HSC survey. In this paper, we investigate the detectability of the 21cm signal with the 21cm-LAE cross-correlation by using our improved reionization simulations. We also focus on the error budget and evaluate it quantitatively in order to consider a strategy to improve the signal-to-noise ratio. In addition, we explore an expansion of the LAE survey to suggest optimal survey parameters and show a potential to measure a characteristic size of ionized bubbles via the turnover scale of the cross-power spectrum. As a result, we find that the Murchison Widefield Array (MWA) has ability to detect the cross-power spectrum signal on large scales by combining LAE Deep field survey of HSC. We also show that the sensitivity is improved dramatically at small scales by adding redshift information from the PFS measurements. The Square Kilometre Array (SKA) has a potential to measure the turnover scale with an accuracy of $6\times10^{-3}~{\rm Mpc^{-1}}$.

The dispersion measure (DM) of high-redshift $(z \gtrsim 6)$ transient objects such as Fast Radio Bursts can be a powerful tool to probe the intergalactic medium during the Epoch of Reionization. In this paper, we study the variance of the DMs of objects with the same redshift as a potential probe of the size distribution of ionized bubbles. We calculate the DM variance with a simple model with randomly-distributed spherical bubbles. It is found that the DM variance reflects the characteristics of the probability distribution of the bubble size. We find the variance can be measured precisely enough to obtain the information on the typical size with a few hundred sources at a single redshift.

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.

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.

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$