Detecting $\rm H_I$ 21cm line in the intergalactic medium (IGM) during the Epoch of Reionization (EoR) suffers from foreground contamination such as Galactic synchrotron and extragalactic radio sources. Cross-correlation between the 21cm line and Lyman-$\alpha$ emitter (LAE) galaxies is a powerful tool to identify the 21cm signal since the 21cm line emission has correlation with LAEs while the LAEs are statistically independent of the foregrounds. So far, the detectability of 21cm-LAE cross-power spectrum has been investigated with simple LAE models where the observed Ly$\alpha$ luminosity is proportional to the dark matter halo mass. However, the previous models were inconsistent with the latest observational data of LAEs obtained with Subaru/Hyper Suprime-Cam (HSC). Here, we revisit the detectability of 21cm-LAE cross-power spectrum adopting a state-of-the-art LAE model consistent with all Subaru/HSC observations such as the Ly$\alpha$ luminosity function, LAE angular auto-correlation, and the LAE fractions in the continuum selected galaxies. We find that resultant cross-power spectrum with the updated LAE model is reduced at small scales ($k\sim 1\ \rm Mpc^{-1}$) compared to the simple models, while the amplitudes at large scales ($k \lesssim 0.2 \ \rm Mpc^{-1}$) are not affected so much. We conclude that the large-scale signal would be detectable with Square Kilometre Array (SKA) and HSC LAE cross-correlation but detecting the small scale signal would require an extended HSC LAE survey with an area of $\sim 75\ \rm deg^2$ or 3000 hrs observation time of 21cm line with SKA.

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$.

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.

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}}$.

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.

We report results of the 2006 April multi-wavelengths campaign of SS 433, focusing on X-ray data observed with Suzaku at two orbital phases (in- and out-of- eclipse) and simultaneous optical spectroscopic observations. By analyzing the Fe25 K_alpha lines originating from the jets, we detect rapid variability of the Doppler shifts, dz/dt ~ 0.019/0.33 day^-1, which is larger than those expected from the precession and/or nodding motion. This phenomenon probably corresponding to "jitter" motions observed for the first time in X-rays, for which significant variability both in the jet angle and intrinsic speed is required. From the time lag of optical Doppler curves from those of X-rays, we estimate the distance of the optical jets from the base to be ~(3-4) \times 10^14 cm. Based on the radiatively cooling jet model, we determine the innermost temperature of the jets to be T_0 = 13 +/- 2 keV and 16 +/- 3 keV (the average of the blue and red jets) for the out-of-eclipse and in-eclipse phase, respectively, from the line intensity ratio of Fe25 K_alpha and Fe26 K_alpha. While the broad band continuum spectra over the 5--40 keV band in eclipse is consistent with a multi-temperature bremsstrahlung emission expected from the jets, and its reflection component from cold matter, the out-of-eclipse spectrum is harder than the jet emission with the base temperature determined above, implying the presence of an additional hard component.

We report the observation of a new dwarf nova, OT J055717+683226, during its first-ever recorded superoutburst in December 2006. Our observation shows that this object is an SU UMa-type dwarf nova having a very short superhump period of 76.67+/- 0.03 min (0.05324+/-0.00002 d). The next superoutburst was observed in March 2008. The recurrence time of superoutbursts (supercycle) is, hence, estimated to be ~480 d. The supercycle is much shorter than those of WZ Sge-type dwarf novae having supercycles of >~ 10 yr, which are a major population of dwarf novae in the shortest orbital period regime (<~85 min). Using a hierarchical cluster analysis, we identified seven groups of dwarf novae in the shortest orbital period regime. We identified a small group of objects that have short supercycles, small outburst amplitudes, and large superhump period excesses, compared with those of WZ Sge stars. OT J055717+683226 probably belongs to this group.

We present results of optical spectroscopic observations of the mass donor star in SS 433 with Subaru and Gemini, with an aim to best constrain the mass of the compact object. Subaru/FOCAS observations were performed on October 6-8 and 10, 2007, covering the orbital phase of phi=0.96-0.26. We first calculate cross correlation function of these spectra with that of the reference star HD 9233 in the wavelength range of 4740-4840 Angstrom. This region is selected to avoid 'strong' absorption lines accompanied with contaminating emission components. The same analysis is applied to archive data of Gemini/GMOS taken at phi=0.84-0.30 by Hillwig & Gies (2008). From the Subaru and Gemini CCF results, the amplitude of radial velocity curve of the donor star is determined to be 58.3+/-3.8 km s-1 with a systemic velocity of 59.2+/-2.5 km s-1. Together with the radial velocity curve of the compact object, we derive the mass of the donor star and compact object to be M_O=12.4+/-1.9 M_sun and M_X=4.3+/-0.6 M_sun, respectively. We conclude, however, that these values should be taken as upper limits. From the analysis of the averaged absorption line profiles of strong lines and weak lines observed with Subaru, we find evidence for heating effects from the compact object. Using a simple model, we find that the true radial velocity amplitude of the donor star could be as low as 40+/-5 km s-1 in order to produce the observed absorption-line profiles. Taking into account the heating of the donor star may lower the derived masses to M_O=10.4 +2.3/-1.9 M_sun and M_X=2.5 +0.7/-0.6 M_sun. Our final constraint, 1.9 M_sun< M_X <4.9 M_sun, indicates that the compact object in SS 433 is most likely a low mass black hole, although the possibility of a massive neutron star cannot be firmly excluded.

We systematically surveyed period variations of superhumps in SU UMa-type dwarf novae based on newly obtained data and past publications. In many systems, the evolution of superhump period are found to be composed of three distinct stages: early evolutionary stage with a longer superhump period, middle stage with systematically varying periods, final stage with a shorter, stable superhump period. During the middle stage, many systems with superhump periods less than 0.08 d show positive period derivatives. Contrary to the earlier claim, we found no clear evidence for variation of period derivatives between superoutburst of the same object. We present an interpretation that the lengthening of the superhump period is a result of outward propagation of the eccentricity wave and is limited by the radius near the tidal truncation. We interpret that late stage superhumps are rejuvenized excitation of 3:1 resonance when the superhumps in the outer disk is effectively quenched. Many of WZ Sge-type dwarf novae showed long-enduring superhumps during the post-superoutburst stage having periods longer than those during the main superoutburst. The period derivatives in WZ Sge-type dwarf novae are found to be strongly correlated with the fractional superhump excess, or consequently, mass ratio. WZ Sge-type dwarf novae with a long-lasting rebrightening or with multiple rebrightenings tend to have smaller period derivatives and are excellent candidate for the systems around or after the period minimum of evolution of cataclysmic variables (abridged).

We report time-resolved CCD photometry of the cataclysmic variable EG Aquarii during the 2006 November outburst During the outburst, superhumps were unambiguously detected with a mean period of 0.078828(6) days, firstly classifying the object as an SU UMa-type dwarf nova. It also turned out that the outburst contained a precursor. At the end of the precursor, immature profiles of humps were observed. By a phase analysis of these humps, we interpreted the features as superhumps. This is the second example that the superhumps were shown during a precursor. Near the maximum stage of the outburst, we discovered an abrupt shift of the superhump period by ${\sim}$ 0.002 days. After the supermaximum, the superhump period decreased at the rate of $\dot{P}/P$=$-8.2{\times}10^{-5}$, which is typical for SU UMa-type dwarf novae. Although the outburst light curve was characteristic of SU UMa-type dwarf novae, long-term monitoring of the variable shows no outbursts over the past decade. We note on the basic properties of long period and inactive SU UMa-type dwarf novae.

We observed an X-ray afterglow of GRB 060904A with the Swift and Suzaku satellites. We found rapid spectral softening during both the prompt tail phase and the decline phase of an X-ray flare in the BAT and XRT data. The observed spectra were fit by power-law photon indices which rapidly changed from $\Gamma = 1.51^{+0.04}_{-0.03}$ to $\Gamma = 5.30^{+0.69}_{-0.59}$ within a few hundred seconds in the prompt tail. This is one of the steepest X-ray spectra ever observed, making it quite difficult to explain by simple electron acceleration and synchrotron radiation. Then, we applied an alternative spectral fitting using a broken power-law with exponential cutoff (BPEC) model. It is valid to consider the situation that the cutoff energy is equivalent to the synchrotron frequency of the maximum energy electrons in their energy distribution. Since the spectral cutoff appears in the soft X-ray band, we conclude the electron acceleration has been inefficient in the internal shocks of GRB 060904A. These cutoff spectra suddenly disappeared at the transition time from the prompt tail phase to the shallow decay one. After that, typical afterglow spectra with the photon indices of 2.0 are continuously and preciously monitored by both XRT and Suzaku/XIS up to 1 day since the burst trigger time. We could successfully trace the temporal history of two characteristic break energies (peak energy and cutoff energy) and they show the time dependence of $\propto t^{-3} \sim t^{-4}$ while the following afterglow spectra are quite stable. This fact indicates that the emitting material of prompt tail is due to completely different dynamics from the shallow decay component. Therefore we conclude the emission sites of two distinct phenomena obviously differ from each other.

We report on time-resolved CCD photometry of four outbursts of a short-period SU UMa-type dwarf nova, V844 Herculis. We successfully determined the mean superhump periods to be 0.05584(64) days, and 0.055883(3) for the 2002 May superoutburst, and the 2006 April-May superoutburst, respectively. During the 2002 October observations, we confirmed that the outburst is a normal outburst, which is the first recorded normal outburst in V844 Her. We also examined superhump period changes during 2002 May and 2006 April-May superoutbursts, both of which showed increasing superhump period over the course of the plateau stage. In order to examine the long-term behavior of V844 Her, we analyzed archival data over the past ten years since the discovery of this binary. Although photometry is not satisfactory in some superoutbursts, we found that V844 Her showed no precursors and rebrightenings. Based on the long-term light curve, we further confirmed V844 Her has shown almost no normal outbursts despite the fact that the supercycle of the system is estimated to be about 300 days. In order to explain the long-term light curves of V844 Her, evaporation in the accretion disk may play a role in the avoidance of several normal outbursts, which does not contradict with the relatively large X-ray luminosity of V844 Her.

Development of MITSuME is reported. Two 50-cm optical telescopes have been built at Akeno in Yamanashi prefecture and at Okayama Astrophysical Observatory (OAO) in Okayama prefecture. Three CCD cameras for simultaneous g'RcIc photometry are to be mounted on each focal plane, covering a wide FOV of about 30" x 30". The limiting magnitude at V is fainter than 18. In addition to these two optical telescopes, a 91-cm IR telescope with a 1 deg x 1 deg field of view is being built at OAO, which performs photometry in YJHK bands. These robotic telescopes can start the observation of counterparts of a GRB within a minute from an alert. We aim to obtain photometric redshifts exceeding 10 with these telescopes. The performance and the current construction status of the telescopes are presented.

A radio-IR-optical-X-ray observation campaign of SS 433 has been performed in April 2006, when the jet axis is almost perpendicular to the line of sight. Five flares have been detected during the campaign by radio monitoring observation with RATAN-600. The X-ray astronomical satellite Suzaku observed the source in and out of eclipse. In the X-ray data out of eclipse, the flux shows a significant variation with a time scale of hours. The source seems to be in the active state during the campaign. The observation logs and preliminary results are presented.

Results are presented of early X-ray afterglow observations of GRB 060105 by Swift and Suzaku. The bright, long gamma-ray burst GRB 060105 triggered the Swift Burst Alert Telescope (BAT) at 06:49:28 on 5 January 2006. The Suzaku team commenced a pre-planned target of opportunity observation at 19 ks (5.3 hr) after the Swift trigger. Following the prompt emission and successive very steep decay, a shallow decay was observed from T_0+187 s to T_0+1287 s. After an observation gap during T_0 +(1.5-3) ks, an extremely early steep decay was observed in T_0+(4-30) ks. The lightcurve flattened again at T_0+30 ks, and another steep decay followed from T_0+50 ks to the end of observations. Both steep decays exhibited decay indices of 2.3 - 2.4. This very early break, if it is a jet break, is the earliest case among X-ray afterglow observations, suggesting a very narrow jet whose opening angle is well below 1 degree. The unique Suzaku/XIS data allow us to set very tight upper limits on line emission or absorption in this GRB. For the reported pseudo-redshift of z=4.0+/-1.3 the upper limit on the iron line equivalent width is 50 eV.

We report on the time-resolved CCD photometry of a newly discovered variable star, TSS J022216.4+412259.9 during the outburst in 2005 November-December brightening. The obtained light curves unambiguously showed 0.2-0.3 mag modulations, which we confirmed to be the superhump observed among SU UMa-type dwarf novae. We also performed a period search for the data obtained during the outburst plateau phase, and revealed the existence of the two periodicities: 0.054868(98) days for the first two nights and 0.055544(26) days for the following plateau phase. This bi-periodicity is hardly observed in usual SU UMa-type dwarf novae, but characteristic of WZ Sge-type stars. We undoubtedly detected a rebrightening in the post-outburst stage, which is typical of short-period SU UMa-type dwarf novae including WZ Sge-type stars. These observations suggests that TSS J022216.4+412259.9 may be a new WZ Sge stars breaking the shortest superhump period of 0.05648 days for V592 Her among this class with a known superhump period so far.

We report on time-resolved photometry of the superoutburst of an SU UMa-type dwarf nova, SDSS J013701.06-091234.9 in 2003 December-2004 January. The obtained light curves definitely show superhumps with a period of 0.056686 (12) d, which is one of the shortest superhump periods among those of SU UMa-type dwarf novae ever observed. Considering quiescent photometric studies, we estimated the fractional superhump excess to be 0.024. Spectroscopic observations by Szkody et al. (2003) provided evidence for TiO bands despite the short orbital period, implying that the system has a luminous secondary star. We draw a color-color diagram of SU UMa-type dwarf novae in quiescence using 2MASS archives, revealing that the location of this star in the color-color diagram is deviated from the general trend. The distance to the system was roughly estimated to be 300${\pm}$80 pc, using the empirical period-absolute magnitude relation and based on the proper motion.