(abridged) The scale of the acoustic oscillation of baryons at the baryon-photon decoupling is imprinted on the spatial distribution of galaxies in the Universe, known as the baryon acoustic oscillation (BAO). The correlation functions and power spectrum are used as a central tool for the studies on the BAO analysis. In this work, we analyzed the spatial distribution of galaxies with a method from the topological data analysis (TDA), in order to detect and examine the BAO signal in the galaxy distribution. The TDA provides a method to treat various types of "holes" in point set data, by constructing the persistent homology (PH) group from the geometric structure of data points and handling the topological information of the dataset. We can obtain the information on the size, position, and statistical significance of the holes in the data. A particularly strong point of the persistent homology is that it can classify the holes by their spatial dimension, i.e., a 0-dim separation, 1-dim loop, 2-dim shell, etc. We first analyzed the simulation datasets with and without the baryon physics to examine the performance of the PH method. We found that the PH is indeed able to detect the BAO signal: simulation data with baryon physics present a prominent signal from the BAO, while data without baryon physics does not show this signal. Then, we applied the PH to a quasar sample at $z <1.0$ from extended Baryon Oscillation Spectroscopic Survey in Sloan Digital Sky Survey Data Release 14. We discovered a characteristic hole (a hollow shell) at a scaler $\sim150\ [{\rm Mpc}]$. This exactly corresponds to the BAO signature imprinted in the galaxy/quasar distribution. We performed this analysis on a small subsample of 2000 quasars. This clearly demonstrates that the PH analysis is very efficient in finding this type of topological structures even if the sampling is very sparse.

Many studies have shown that there are clear sequences in the period-luminosity relationship (PLR) for Mira variables and semiregular variables (SRVs) in the Large Magellanic Cloud (LMC). To investigate the PLR for SRVs in our galaxy, we examined the annual parallax measurement and conducted K'-band photometric monitoring of an SRV star SV Pegasus (SV Peg). We measured the position change of the associating H$_2$O maser spots by phase-referencing VLBI observations with VERA at 22 GHz, spanning approximately 3 yr, and detected an annual parallax of $\pi = 3.00 \pm 0.06$ mas, corresponding to a distance of $D=333 \pm 7$ pc. This result is in good agreement with the Hipparcos parallax and improves the accuracy of the distance from 35 $\%$ to 2 $\%$. However, the GAIA DR2 catalog gave a parallax of $\pi=1.12\pm0.28$ mas for SV Peg. This indicates that the GAIA result might be blurred by the effect of the stellar size because the estimated stellar radius was $\sim 5$ mas, which is comparable to the parallax. We obtained a K'-band mean magnitude of $m_{K'} = -0.48$ mag and a period of $P=177$ days from our photometric monitoring with a 1-m telescope. Using the trigonometric distance, we derived an absolute magnitude of $M_{K'}=-8.09 \pm 0.05$ mag. This result shows that the position of SV Peg in the PLR falls on the C' sequence found in the PLR in the LMC, which is similar to other SRVs in our galaxy.

We present monitoring campaign observations at optical and near-infrared (NIR) wavelengths for a radio-loud active galactic nucleus (AGN) at z=0.840, SDSS~J110006.07+442144.3 (hereafter, J1100+4421), which was identified during a flare phase in late February, 2014. The campaigns consist of three intensive observing runs from the discovery to March, 2015, mostly within the scheme of the OISTER collaboration. Optical-NIR light curves and simultaneous spectral energy distributions (SEDs) are obtained. Our measurements show the strongest brightening in March, 2015. We found that the optical-NIR SEDs of J1100+4421 show an almost steady shape despite the large and rapid intranight variability. This constant SED shape is confirmed to extend to $\sim5~\mu$m in the observed frame using the archival WISE data. Given the lack of absorption lines and the steep power-law spectrum of $\alpha_{\nu}\sim-1.4$, where $f_{\nu}\propto\nu^{\alpha_{\nu}}$, synchrotron radiation by a relativistic jet with no or small contributions from the host galaxy and the accretion disk seems most plausible as an optical-NIR emission mechanism. The steep optical-NIR spectral shape and the large amplitude of variability are consistent with this object being a low $\nu_{\rm{peak}}$ jet-dominated AGN. In addition, sub-arcsec resolution optical imaging data taken with Subaru Hyper Suprime-Cam does not show a clear extended component and the spatial scales are significantly smaller than the large extensions detected at radio wavelengths. The optical spectrum of a possible faint companion galaxy does not show any emission lines at the same redshift and hence a merging hypothesis for this AGN-related activity is not supported by our observations.

We present the early-phase spectra and the light curves of the broad-lined supernova (SN) 2016coi from $t=7$ to $67$ days after the estimated explosion date. This SN was initially reported as a broad-lined Type SN Ic (SN Ic-BL). However, we found that spectra up to $t=12$ days exhibited the He~{\sc i} $\lambda$5876, $\lambda$6678, and $\lambda$7065 absorption lines. We show that the smoothed and blueshifted spectra of normal SNe Ib are remarkably similar to the observed spectrum of SN 2016coi. The line velocities of SN 2016coi were similar to those of SNe Ic-BL and significantly faster than those of SNe Ib. Analyses of the line velocity and light curve suggest that the kinetic energy and the total ejecta mass of SN 2016coi are similar to those of SNe Ic-BL. Together with broad-lined SNe 2009bb and 2012ap for which the detection of He~{\sc i} were also reported, these SNe could be transitional objects between SNe Ic-BL and SNe Ib, and be classified as broad-lined Type `Ib' SNe (SNe `Ib'-BL). Our work demonstrates the diversity of the outermost layer in broad-lined SNe, which should be related to the variety of the evolutionary paths.

Solar neutrons have been detected using the neutron monitor located at Mt. Chacaltaya, Bolivia, in association with a large solar flare on November 24, 2000. This is the first detection of solar neutrons by the neutron monitor that have been reported so far in solar cycle 23. The statistical significance of the detection is 5.5 sigma. In this flare, the intense emission of hard X-rays and gamma-rays has been observed by the Yohkoh Hard X-ray Telescope (HXT) and Gamma Ray Spectrometer (GRS), respectively. The production time of solar neutrons is better correlated with those of hard X-rays and gamma-rays than with the production time of soft X-rays. The observations of the solar neutrons on the ground have been limited to solar flares with soft X-ray class greater than X8 in former solar cycles. In this cycle, however, neutrons were detected associated with an X2.3 solar flare on November 24, 2000. This is the first report of the detection of solar neutrons on the ground associated with a solar flare with its X-ray class smaller than X8.

With the Akeno Giant Air Shower Array (AGASA), 581 cosmic rays above 10^19eV, 47 above 4 x 10^19eV, and 7 above 10^20eV are observed until August 1998. Arrival direction distribution of these extremely high energy cosmic rays has been studied. While no significant large-scale anisotropy is found on the celestial sphere, some interesting clusters of cosmic rays are observed. Above 4 x 10^19eV, there are one triplet and three doublets within separation angle of 2.5^o and the probability of observing these clusters by a chance coincidence under an isotropic distribution is smaller than 1 %. Especially the triplet is observed against expected 0.05 events. The cos(\theta_GC) distribution expected from the Dark Matter Halo model fits the data as well as an isotropic distribution above 2 x 10^19eV and 4 x 10^19eV, but is a poorer fit than isotropy above 10^19eV. Arrival direction distribution of seven 10^20eV cosmic rays is consistent with that of lower energy cosmic rays and is uniform. Three of seven are members of doublets above about 4 x 10^19eV.

Anisotropy in the arrival directions of cosmic rays around 10^{18}eV is studied using data from the Akeno 20 km^2 array and the Akeno Giant Air Shower Array (AGASA), using a total of about 216,000 showers observed over 15 years above 10^{17}eV. In the first harmonic analysis, we have found significant anisotropy of $\sim$ 4 % around 10^{18}eV, corresponding to a chance probability of $\sim 10^{-5}$ after taking the number of independent trials into account. With two dimensional analysis in right ascension and declination, this anisotropy is interpreted as an excess of showers near the directions of the Galactic Center and the Cygnus region. This is a clear evidence for the existence of the galactic cosmic ray up to the energy of 10^{18}eV. Primary particle which contribute this anisotropy may be proton or neutron.

Anisotropy in the arrival directions of cosmic rays with energies above 10$^{17}$eV is studied using data from the Akeno 20 km$^2$ array and the Akeno Giant Air Shower Array (AGASA), using a total of about 117,000 showers observed during 11 years. In the first harmonic analysis, we have found strong anisotropy of $\sim$ 4% around 10$^{18}$eV, corresponding to a chance probability of 0.2%. With two dimensional analysis in right ascension and declination, this anisotropy is interpreted as an excess of showers near the directions of the Galactic Center and the Cygnus region.

The cosmic-ray energy spectrum above 10^{18.5} eV is reported using the updated data set of the Akeno Giant Air Shower Array (AGASA) from February 1990 to October 1997. The energy spectrum extends beyond 10^{20} eV and the energy gap between the highest energy event and the others is being filled up with recently observed events. The spectral shape suggests the absence of the 2.7 K cutoff in the energy spectrum or a possible presence of a new component beyond the 2.7 K cutoff.