Waseda University Nasu telescope array is a spatial fast Fourier transform (FFT) interferometer consisting of eight linearly aligned antennas with 20 m spherical dishes. This type of interferometer was developed to survey transient radio sources with an angular resolution as high as that of a 160 m dish with a field of view as wide as that of a 20 m dish. We have been performing drift-scan-mode observations, in which the telescope scans the sky around a selected declination as the earth rotates. The black hole X-ray binary V404 Cygni underwent a new outburst in 2015 June after a quiescent period of 26 years. Because of the interest in black hole binaries, a considerable amount of data on this outburst at all wavelengths was accumulated. Using the above telescope, we had been monitoring V404 Cygni daily from one month before the X-ray outburst, and two radio flares at 1.4 GHz were detected on June 21.73 and June 26.71. The flux density and time-scale of each flare were 313+/-30 mJy and 1.50+/-0.49 days, 364+/-30 mJy and 1.70+/-0.16 days, respectively. We have also confirmed the extreme variation of radio spectra within a short period by collecting other radio data observed with several radio telescopes. Such spectral behaviors are considered to reflect the change in the opacity of the ejected blobs associated with these extreme activities in radio and X-ray. Our 1.4 GHz radio data are expected to be helpful for studying the physics of the accretion and ejection phenomena around black holes.

We present the results of a dust-reverberation survey of quasars at redshifts z<0.6. We found a delayed response of the K-band flux variation after the optical flux variation in 25 out of 31 targets, and obtained the lag time between them for 22 targets. Combined with the results for nearby Seyfert galaxies, we provide the largest homogeneous collection of K-band dust-reverberation data for 36 type 1 active galactic nuclei (AGNs). This doubles the sample and includes the most distant AGN and the largest lag so far measured. We estimated the optical luminosity of the AGN component of each target using three different methods: spectral decomposition, the flux-variation-gradient method, and image decomposition. We found a strong correlation between the reverberation radius for the innermost dust torus and the optical luminosity over a range of approximately four orders of magnitude in luminosity, as is already known for Seyfert galaxies. We estimated the luminosity distances of the AGNs based on their dust-reverberation lags, and found that the data in the redshift-distance diagram are consistent with the current standard estimates of the cosmological parameters. We also present the radius-luminosity relations for isotropic luminosity indicators such as the hard X-ray (14--195 keV), [OIV] 25.89 um, and mid-infrared (12 um) continuum luminosities, which are applicable to obscured AGNs.

Using a prototype of the Tomo-e Gozen wide-field CMOS mosaic camera, we acquire wide-field optical images at a cadence of 2 Hz and search them for transient sources of duration 1.5 to 11.5 seconds. Over the course of eight nights, our survey encompasses the equivalent of roughly two days on one square degree, to a fluence equivalent to a limiting magnitude about $V = 15.6$ in a 1-second exposure. After examining by eye the candidates identified by a software pipeline, we find no sources which meet all our criteria. We compute upper limits to the rate of optical transients consistent with our survey, and compare those to the rates expected and observed for representative sources of ephemeral optical light.

We report observations of a stellar occultation by the classical Kuiper belt object (50000) Quaoar occurred on 28 June 2019. A single-chord high-cadence (2 Hz) photometry dataset was obtained with the Tomo-e Gozen CMOS camera mounted on the 1.05 m Schmidt telescope at Kiso Observatory. The obtained ingress and egress data do not show any indication of atmospheric refraction and allow to set new $1\sigma$ and $3\sigma$ upper limits of 6 and 16 nbar, respectively, for the surface pressure of a pure methane atmosphere. These upper limits are lower than the saturation vapor pressure of methane at Quaoar's expected mean surface temperature ($T \sim 44$ K) and imply the absence of a $\sim$10 nbar-level global atmosphere formed by methane ice on Quaoar's surface.

An observational result of a radio continuum and H92$\alpha$ radio recombination line of the Galactic Center Lobe (GCL), using the Yamaguchi 32 m radio telescope, is reported. The obtained spatial intensity distribution of the radio recombination line shows two distinctive ridge-like structures extending from the galactic plane vertically to the north at the eastern and western sides of the galactic center, which are connected to each other at a latitude of $1.2^{\circ}$ to form a loop-like structure as a whole. This suggests that most of the radio continuum emission of the GCL is free-free emission, and that the GCL is filled with thermal plasma. The east ridge of the GCL observed with the radio recombination line separates 30 pc from the radio arc, which has been considered as a part of the GCL, but coincides with a ridge of the radio continuum at a galactic longitude of $0^{\circ}$. The radial velocity of the radio recombination line is found to be between $-4$ and $+10$ km s$^{-1}$ across the GCL. This velocity is much smaller than the one expected from the galactic rotation, and hence indicates that the GCL exists apart from the galactic center. These characteristics of the GCL suggest that the long-standing hypothesis that the GCL was created by an explosive activity in the galactic center is unlikely, but favor that the GCL is a giant HII region.

Faraday tomography is thought to be a powerful tool to explore cosmic magnetic field. Broadband radio polarimetric data is essential to ensure the quality of Faraday tomography, but such data is not easy to obtain because of radio frequency interferences (RFIs). In this paper, we investigate optimum frequency coverage of Faraday tomography so as to explore Faraday rotation measure (RM) due to the intergalactic magnetic field (IGMF) in filaments of galaxies. We adopt a simple model of the IGMF and estimate confidence intervals of the model parameters using the Fisher information matrix. We find that meaningful constraints for RM due to the IGMF are available with data at multiple narrowbands which are scattered over the ultra-high frequency (UHF, 300 MHz - 3000 MHz). The optimum frequency depends on the Faraday thickness of the Milky Way foreground. These results are obtained for a wide brightness range of the background source including fast radio bursts (FRBs). We discuss the relation between the polarized-intensity spectrum and the optimum frequency.

Imaging observations of faint meteors were carried out on April 11 and 14, 2016 with a wide-field CMOS mosaic camera, Tomo-e PM, mounted on the 105-cm Schmidt telescope at Kiso Observatory, the University of Tokyo. Tomo-e PM, which is a prototype model of Tomo-e Gozen, can monitor a sky of ${\sim}1.98\,\mathrm{deg^2}$ at 2\,Hz. The numbers of detected meteors are 1514 and 706 on April 11 and 14, respectively. The detected meteors are attributed to sporadic meteors. Their absolute magnitudes range from $+4$ to $+10\,\mathrm{mag}$ in the $V$-band, corresponding to about $8.3{\times}10^{-2}$ to $3.3{\times}10^{-4}\,\mathrm{g}$ in mass. The present magnitude distributions we obtained are well explained by a single power-law luminosity function with a slope parameter $r = 3.1{\pm}0.4$ and a meteor rate $\log_{10}N_0 = -5.5{\pm}0.5$. The results demonstrate a high performance of telescopic observations with a wide-field video camera to constrain the luminosity function of faint meteors. The performance of Tomo-e Gozen is about two times higher than that of Tomo-e PM. A survey with Tomo-e Gozen will provide a more robust measurement of the luminosity function.

To search for giant X-ray pulses correlated with the giant radio pulses (GRPs) from the Crab pulsar, we performed a simultaneous observation of the Crab pulsar with the X-ray satellite Hitomi in the 2 -- 300 keV band and the Kashima NICT radio observatory in the 1.4 -- 1.7 GHz band with a net exposure of about 2 ks on 25 March 2016, just before the loss of the Hitomi mission.The timing performance of the Hitomi instruments was confirmed to meet the timing requirement and about 1,000 and 100 GRPs were simultaneously observed at the main and inter-pulse phases, respectively, and we found no apparent correlation between the giant radio pulses and the X-ray emission in either the main or inter-pulse phases.All variations are within the 2 sigma fluctuations of the X-ray fluxes at the pulse peaks, and the 3 sigma upper limits of variations of main- or inter- pulse GRPs are 22\% or 80\% of the peak flux in a 0.20 phase width, respectively, in the 2 -- 300 keV band.The values become 25\% or 110\% for main or inter-pulse GRPs, respectively, when the phase width is restricted into the 0.03 phase.Among the upper limits from the Hitomi satellite, those in the 4.5-10 keV and the 70-300 keV are obtained for the first time, and those in other bands are consistent with previous reports.Numerically, the upper limits of main- and inter-pulse GRPs in the 0.20 phase width are about (2.4 and 9.3) $\times 10^{-11}$ erg cm$^{-2}$, respectively. No significant variability in pulse profiles implies that the GRPs originated from a local place within the magnetosphere and the number of photon-emitting particles temporally increases.However, the results do not statistically rule out variations correlated with the GRPs, because the possible X-ray enhancement may appear due to a $>0.02$\% brightening of the pulse-peak flux under such conditions.

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.

In Yoshii et al. (2014), we described a new method for measuring extragalactic distances based on dust reverberation in active galactic nuclei (AGNs), and we validated our new method with Cepheid variable stars. In this paper, we validate our new method with Type Ia supernovae (SNe Ia) which occurred in two of the AGN host galaxies during our AGN monitoring program: SN 2004bd in NGC 3786 and SN 2008ec in NGC 7469. Their multicolor light curves were observed and analyzed using two widely accepted methods for measuring SN distances, and the distance moduli derived are $\mu=33.47\pm 0.15$ for SN 2004bd and $33.83\pm 0.07$ for SN 2008ec. These results are used to obtain independently the distance measurement calibration factor, $g$. The $g$ value obtained from the SN Ia discussed in this paper is $g_{\rm SN} = 10.61\pm 0.50$ which matches, within the range of 1$\sigma$ uncertainty, $g_{\rm DUST} = 10.60$, previously calculated ab initio in Yoshii et al. (2014). Having validated our new method for measuring extragalactic distances, we use our new method to calibrate reverberation distances derived from variations of H$\beta$ emission in the AGN broad line region (BLR), extending the Hubble diagram to $z\approx 0.3$ where distinguishing between cosmologies is becoming possible.

The Square Kilometre Array will revolutionize pulsar studies with its wide field-of-view, wide-band observation and high sensitivity, increasing the number of observable pulsars by more than an order of magnitude. Pulsars are of interest not only for the study of neutron stars themselves but for their usage as tools for probing fundamental physics such as general relativity, gravitational waves and nuclear interaction. In this article, we summarize the activity and interests of SKA-Japan Pulsar Science Working Group, focusing on an investigation of modified gravity theory with the supermassive black hole in the Galactic Centre, gravitational-wave detection from cosmic strings and binary supermassive black holes, a study of the physical state of plasma close to pulsars using giant radio pulses and determination of magnetic field structure of Galaxy with pulsar pairs.

We have developed a near-infrared camera called ANIR (Atacama Near-InfraRed camera) for the University of Tokyo Atacama Observatory 1.0m telescope (miniTAO) installed at the summit of Cerro Chajnantor (5640 m above sea level) in northern Chile. The camera provides a field of view of 5'.1 $\times$ 5'.1 with a spatial resolution of 0".298 /pixel in the wavelength range of 0.95 to 2.4 $\mu$m. Taking advantage of the dry site, the camera is capable of hydrogen Paschen-$\alpha$ (Pa$\alpha$, $\lambda=$1.8751 $\mu$m in air) narrow-band imaging observations, at which wavelength ground-based observations have been quite difficult due to deep atmospheric absorption mainly from water vapor. We have been successfully obtaining Pa$\alpha$ images of Galactic objects and nearby galaxies since the first-light observation in 2009 with ANIR. The throughputs at the narrow-band filters ($N1875$, $N191$) including the atmospheric absorption show larger dispersion (~10%) than those at broad-band filters (a few %), indicating that they are affected by temporal fluctuations in Precipitable Water Vapor (PWV) above the site. We evaluate the PWV content via the atmospheric transmittance at the narrow-band filters, and derive that the median and the dispersion of the distribution of the PWV are 0.40+/-0.30 mm for $N1875$ and 0.37+/-0.21 mm for $N191$, which are remarkably smaller (49+/-38% for $N1875$ and 59+/-26% for $N191$) than radiometry measurements at the base of Cerro Chajnantor (5100 m alt.). The decrease in PWV can be explained by the altitude of the site when we assume that the vertical distribution of the water vapor is approximated at an exponential profile with scale heights within 0.3-1.9 km (previously observed values at night). We thus conclude that miniTAO/ANIR at the summit of Cerro Chajnantor indeed provides us an excellent capability for a "ground-based" Pa$\alpha$ observation.

Luminous infrared galaxies (LIRGs) are enshrouded by a large amount of dust, produced by their active star formation, and it is difficult to measure their activity in the optical wavelength. We have carried out Pa$\alpha$ narrow-band imaging observations of 38 nearby star-forming galaxies including 33 LIRGs listed in $IRAS$ RBGS catalog with the Atacama Near InfraRed camera (ANIR) on the University of Tokyo Atacama Observatory (TAO) 1.0 m telescope (miniTAO). Star formation rates (SFRs) estimated from the Pa$\alpha$ fluxes, corrected for dust extinction using the Balmer Decrement Method (typically $A_V$ $\sim$ 4.3 mag), show a good correlation with those from the bolometric infrared luminosity of $IRAS$ data within a scatter of 0.27 dex. This suggests that the correction of dust extinction for Pa$\alpha$ flux is sufficient in our sample. We measure the physical sizes and the surface density of infrared luminosities ($\Sigma_{L(\mathrm{IR})}$) and $SFR$ ($\Sigma_{SFR}$) of star-forming region for individual galaxies, and find that most of the galaxies follow a sequence of local ultra luminous or luminous infrared galaxies (U/LIRGs) on the $L(\mathrm{IR})$-$\Sigma_{L(\mathrm{IR})}$ and $SFR$-$\Sigma_{SFR}$ plane. We confirm that a transition of the sequence from normal galaxies to U/LIRGs is seen at $L(\mathrm{IR})=8\times10^{10}$ $L_{\odot}$. Also, we find that there is a large scatter in physical size, different from those of normal galaxies or ULIRGs. Considering the fact that most of U/LIRGs are merging or interacting galaxies, this scatter may be caused by strong external factors or differences of their merging stage.

We present our discovery of dramatic variability in SDSS J1100+4421 by the high-cadence transient survey Kiso Supernova Survey (KISS). The source brightened in the optical by at least a factor of three within about half a day. Spectroscopic observations suggest that this object is likely a narrow-line Seyfert 1 galaxy (NLS1) at z=0.840, however with unusually strong narrow emission lines. The estimated black hole mass of ~ 10^7 Msun implies bolometric nuclear luminosity close to the Eddington limit. SDSS J1100+4421 is also extremely radio-loud, with a radio loudness parameter of R ~ 4 x 10^2 - 3 x 10^3, which implies the presence of relativistic jets. Rapid and large-amplitude optical variability of the target, reminiscent of that found in a few radio- and gamma-ray loud NLS1s, is therefore produced most likely in a blazar-like core. The 1.4 GHz radio image of the source shows an extended structure with a linear size of about 100 kpc. If SDSS J1100+4421 is a genuine NLS1, as suggested here, this radio structure would then be the largest ever discovered in this type of active galaxies.

The Kiso Supernova Survey (KISS) is a high-cadence optical wide-field supernova (SN) survey. The primary goal of the survey is to catch the very early light of a SN, during the shock breakout phase. Detection of SN shock breakouts combined with multi-band photometry obtained with other facilities would provide detailed physical information on the progenitor stars of SNe. The survey is performed using a 2.2x2.2 deg field-of-view instrument on the 1.05-m Kiso Schmidt telescope, the Kiso Wide Field Camera (KWFC). We take a three-minute exposure in g-band once every hour in our survey, reaching magnitude g~20-21. About 100 nights of telescope time per year have been spent on the survey since April 2012. The number of the shock breakout detections is estimated to be of order of 1 during our 3-year project. This paper summarizes the KISS project including the KWFC observing setup, the survey strategy, the data reduction system, and CBET-reported SNe discovered so far by KISS.

We present the results of a dust reverberation survey for 17 nearby Seyfert 1 galaxies, which provides the largest homogeneous data collection for the radius of the innermost dust torus. A delayed response of the K-band light curve after the V-band light curve was found for all targets, and 49 measurements of lag times between the flux variation of the dust emission in the K band and that of the optical continuum emission in the V band were obtained. The lag times strongly correlated with the optical luminosity in the luminosity range of M_V=-16 to -22 mag, and the regression analysis was performed to obtain the correlation log $\Delta t$ (days) = -2.11 -0.2 M_V assuming $\Delta t \propto L^{0.5}$, which was theoretically expected. We discuss the possible origins of the intrinsic scatter of the dust lag-luminosity correlation, which was estimated to be about 0.13 dex, and we find that the difference of internal extinction and delayed response of changes in lag times to the flux variations could have partly contributed to intrinsic scatter. However, we could not detect any systematic change of the correlation with the subclass of the Seyfert type or the Eddington ratio. Finally, we compare the dust reverberation radius with the near-infrared interferometric radius of the dust torus and the reverberation radius of broad Balmer emission lines. The interferometric radius in the K band was found to be systematically larger than the dust reverberation radius in the same band by about a factor of two, which could be interpreted by the difference between the flux-weighted radius and the response-weighted radius of the innermost dust torus. The reverberation radius of the broad Balmer emission lines was found to be systematically smaller than the dust reverberation radius by about a factor of 4-5, which strongly supports the unified scheme of the Seyfert type of active galactic nuclei. (Abridged)

We carried out the photometric observations of the SU UMa-type dwarf nova ER UMa during 2011 and 2012, which showed the existence of persistent negative superhumps even during the superoutburst. We performed two-dimensional period analysis of its light curves by using a method called "least absolute shrinkage and selection operator" (Lasso) and "phase dispersion minimization" (PDM) analysis, and we found that the period of negative superhumps systematically changed between a superoutburst and the next superoutburst. The trend of the period change can beinterpreted as reflecting the change of the disk radius. This change of the disk radius is in good agreement with the predicted change of the disk radius by the thermal-tidal instability (TTI) model. The normal outbursts within a supercycle showed a general trend that the rising rate to maximum becomes slower as the next superoutburst approaches. The change can be interpreted as the consequence of the increased gas-stream flow onto the inner region of the disk as the result of the tilted disk. Some of the superoutbursts were found to be triggered by a precursor normal outburst when the positive superhumps appeared to develop. The positive and negative superhumps co-existed during the superoutburst. The positive superhumps were prominent only during four or five days after the supermaximum, while the signal of the negative superhumps became strong after the middle phase of the superoutburst plateau. A simple combination of the positive and negative superhumps was found to be insufficient in reproducing the complex profile variation. We were able to detect the developing phase of positive superhumps (stage A superhumps) for the first time in ER UMa-type dwarf novae. Using the period of stage A superhumps, we obtained a mass ratio of 0.100(15), which indicates that ER UMa is on the ordinary evolutional track of CVs.

We present new data for five under-luminous type II-plateau supernovae (SNe IIP), namely SN 1999gn, SN 2002gd, SN 2003Z, SN 2004eg and SN 2006ov. This new sample of low-luminosity SNe IIP (LL SNe IIP) is analyzed together with similar objects studied in the past. All of them show a flat light curve plateau lasting about 100 days, an under luminous late-time exponential tail, intrinsic colours that are unusually red, and spectra showing prominent and narrow P-Cygni lines. A velocity of the ejected material below 10^3 km/s is inferred from measurements at the end of the plateau. The 56Ni masses ejected in the explosion are very small (less than 10^-2 solar masses). We investigate the correlations among 56Ni mass, expansion velocity of the ejecta and absolute magnitude in the middle of the plateau, confirming the main findings of Hamuy (2003), according to which events showing brighter plateau and larger expansion velocities are expected to produce more 56Ni. We propose that these faint objects represent the low luminosity tail of a continuous distribution in parameters space of SNe IIP. The physical properties of the progenitors at the explosion are estimated through the hydrodynamical modeling of the observables for two representative events of this class, namely SN 2005cs and SN 2008in. We find that the majority of LL SNe IIP, and quite possibly all, originate in the core-collapse of intermediate mass stars, in the mass range 10-15 solar masses.

Continuing the project described by Kato et al. (2009, arXiv:0905.1757), we studied 86 SU UMa-type dwarf novae. We confirmed the general trends such as the relation between period derivatives and orbital periods. There are some systems showing positive period derivatives despite the long orbital periods. We observed the 2011 outburst of the WZ Sge-type dwarf nova BW Scl, and recorded an O-C diagram similar to those of previously known WZ Sge-type dwarf novae. The WZ Sge-type dwarf nova OT J184228.1+483742 showed an unusual pattern of double outbursts composed of an outburst with early superhumps and one with ordinary superhumps. We propose an interpretation that a very small growth rate of the 3:1 resonance due to an extremely low mass-ratio led to a quenching of the superoutburst before the ordinary superhumps appeared. We studied ER UMa-type dwarf novae and found that V1159 Ori showed positive superhumps similar to ER UMa in the 1990s. The recently recognized ER UMa-type object BK Lyn dominantly showed negative superhumps, and its behavior was very similar to the present-day state of ER UMa. The pattern of period variations in AM CVn-type objects were very similar to short-period hydrogen-rich SU UMa-type dwarf novae, making them helium analogue of hydrogen-rich SU UMa-type dwarf novae. SBS 1108+574, a peculiar hydrogen-rich dwarf nova below the period minimum, showed a very similar pattern of period variations to those of short-period SU UMa-type dwarf novae. The mass-ratio derived from the detected orbital period suggests that this secondary is a somewhat evolved star whose hydrogen envelope was mostly stripped during the mass-exchange. CC Scl, MASTER OT J072948.66+593824.4 and OT J173516.9+154708 showed only low-amplitude superhumps with complex profiles. These superhumps are likely a combination of closely separated two periods.

Continuing the project described by Kato et al. (2009, PASJ 61, S395, arXiv:0905.1757), we collected times of superhump maxima for 51 SU UMa-type dwarf novae mainly observed during the 2010-2011 season. Although most of the new data for systems with short superhump periods basically confirmed the findings by Kato et al. (2009) and Kato et al. (2010, PASJ 62, 1525, arXiv:1009.5444), the long-period system GX Cas showed an exceptionally large positive period derivative. An analysis of public Kepler data of V344 Lyr and V1504 Cyg yielded less striking stage transitions. In V344 Lyr, there was prominent secondary component growing during the late stage of superoutbursts, and the component persisted at least for two more cycles of successive normal outbursts. We also investigated the superoutbursts of two conspicuous eclipsing objects: HT Cas and the WZ Sge-type object SDSS J080434.20+510349.2. Strong beat phenomena were detected in both objects, and late-stage superhumps in the latter object had an almost constant luminosity during the repeated rebrightenings. The WZ Sge-type object SDSS J133941.11+484727.5 showed a phase reversal around the rapid fading from the superoutburst. The object showed a prominent beat phenomenon even after the end of the superoutburst. A pilot study of superhump amplitudes indicated that the amplitudes of superhumps are strongly correlated with orbital periods, and the dependence on the inclination is weak in systems with inclinations smaller than 80 deg.

We report the first observational search for tau neutrinos from gamma ray bursts (GRBs) using one of the Ashra light collectors. The earth-skimming tau-neutrino technique of imaging Cherenkov tau showers was applied as a detection method. We set stringent upper limits on the tau-neutrino fluence in PeV-EeV region for 3780 s (between 2.83 and 1.78 hours before) and another 3780 s (between 21.2 and 22.2 hours after) surrounding GRB081203A triggered by the Swift satellite. This first search for PeV-EeV tau neutrino complements other experiments in energy range and methodology, and suggests the prologue of "multi-particle astronomy" with a precise determination of time and location.

We examine whether the spectral energy distribution of optical continuum emission of active galactic nuclei (AGNs) changes during flux variation, based on accurate and frequent monitoring observations of 11 nearby Seyfert galaxies and QSOs carried out in the B, V, and I bands for seven years by the MAGNUM telescope. The multi-epoch flux data in any two different bands obtained on the same night show a very tight linear flux to flux relationship for all target AGNs. The flux of the host galaxy within the photometric aperture is carefully estimated by surface brightness fitting to available high-resolution HST images and MAGNUM images. The flux of narrow emission lines in the photometric bands is also estimated from available spectroscopic data. We find that the non-variable component of the host galaxy plus narrow emission lines for all target AGNs is located on the fainter extension of the linear regression line of multi-epoch flux data in the flux to flux diagram. This result strongly indicates that the spectral shape of AGN continuum emission in the optical region does not systematically change during flux variation. The trend of spectral hardening that optical continuum emission becomes bluer as it becomes brighter, which has been reported by many studies, is therefore interpreted as the domination of the variable component of the nearly constant spectral shape of an AGN as it brightens over the non-variable component of the host galaxy plus narrow lines, which is usually redder than AGN continuum emission.

The long-term optical and near infrared monitoring observations for a type 1 act ive galactic nucleus NGC 4151 were carried out for six years from 2001 to 2006 b y using the MAGNUM telescope, and delayed response of flux variations in the $K(2.2\mu m)$ band to those in the $V(0.55\mu m)$ band was clearly detected. Based on cross correlation analysis, we precisely measured a lag time $\Delta t$ for eight separate periods, and we found that $\Delta t$ is not constant changing be tween 30 and 70 days during the monitoring period. Since $\Delta t$ is the ligh t travel time from the central energy source out to the surrounding dust torus, this is the first convincing evidence that the inner radius of dust torus did ch ange in an individual AGN. In order to relate such a change of $\Delta t$ with a change of AGN luminosity $L$, we presented a method of taking an average of th e observed $V$-band fluxes that corresponds to the measured value of $\Delta t$, and we found that the time-changing track of NGC 4151 in the $\Delta t$ versus $L$ diagram during the monitoring period deviates from the relation of $\Delta t \propto L^{0.5}$ expected from dust reverberation. This result, combined with t he elapsed time from period to period for which $\Delta t$ was measured, indicat es that the timescale of dust formation is about one year, which should be taken into account as a new constraint in future studies of dust evolution in AGNs.

Using homogeneous CCD photometric data from the 105-cm Kiso Schmidt telescope covering a 50' x 50' field, we study the mass functions (MFs) of nine open clusters. The ages and Galactocentric distances of the target clusters vary from 16 - 2000 Myr and 9-10.8 kpc, respectively. The values of MF slopes vary from -1.1 to -2.1. The classical value derived by Salpeter (1955) for the slope of the IMF is \Gamma = -1.35. The MFs in the outer regions of the clusters are found to be steeper than in the inner regions, indicating the presence of mass segregation in the clusters.The MF slopes (in the outer region as well as the whole cluster) undergo an exponential decay with the evolutionary parameter \tau (= age/ relaxation time). It seems that the evaporation of low-mass members from outer regions of the clusters is not significant at larger Galactocentric distances. It is concluded that the initial mass function (IMF) in the anticentre direction of the Galaxy might have been steeper than the IMF in the opposite direction. A comparison of the observed CMDs of the clusters with synthetic CMDs gives a photometric binary content of ~40%.

We study the ``normal'' decay phase of the X-ray afterglows of gamma-ray bursts (GRBs), which follows the shallow decay phase, using the events simultaneously observed in the R-band. The classical external shock model -- in which neither the delayed energy injection nor time-dependency of shock micro-physics is considered -- shows that the decay indices of the X-ray and R-band light curves, $\alpha_{\rm X}$ and $\alpha_{\rm O}$, obey a certain relation, and that in particular, $\alpha_{\rm O}-\alpha_{\rm X}$ should be larger than -1/4 unless the ambient density increases with the distance from the central engine. For our selected 14 samples, we have found that 4 events violate the limit at more than the 3$\sigma$ level, so that a fraction of events are outliers of the classical external shock model at the ``normal'' decay phase.

We present the most successful infrared cloud monitor for a robotic telescope. This system was originally developed for the MAGNUM 2-m telescope, which has been achieving unmanned and automated monitoring observation of active galactic nuclei at Haleakala on the Hawaiian island of Maui since 2001. Using a thermal imager and two aspherical mirrors, it at once sees almost the whole sky at a wavelength of $\lambda\sim 10\mu{\rm m}$. Its outdoor part is weather-proof and is totally maintenance-free. The images obtained every one or two minutes are analysed immediately into several ranks of weather condition, from which our automated observing system not only decides to open or close the dome, but also selects what types of observations should be done. The whole-sky data accumulated over four years show that 50$-$60 % of all nights are photometric, and about 75 % are observable with respect to cloud condition at Haleakala. Many copies of this system are now used all over the world such as Mauna Kea in Hawaii, Atacama in Chile, and Okayama and Kiso in Japan.

Observations of the optical afterglow of GRB 041006 with the Kiso Observatory 1.05 m Schmidt telescope, the Lulin Observatory 1.0 m telescope and the Xinglong Observatory 0.6 m telescope. Three-bands (B, V and R) of photometric data points were obtained on 2004 October 6, 0.025-0.329 days after the burst. These very early multi band light curves imply the existence of a color dependent plateau phase. The B-band light curve shows a clear plateau at around 0.03 days after the burst. The R band light curve shows the hint of a plateau, or a possible slope change, at around 0.1 days after the burst. The overall behavior of these multi-band light curves may be interpreted in terms of the sum of two separate components, one showing a monotonic decay the other exhibiting a rising and a falling phase, as described by the standard afterglow model.

Multi-wavelength B, V, R, I observations of the optical afterglow of GRB 050319 were performed by the 1.05-m telescope at Kiso Observatory and the 1.0-m telescope at Lulin Observatory from 1.31 hours to 9.92 hours after the burst. Our R band lightcurves, combined with other published data, can be described by the smooth broken power-law function, with $\alpha_1$ = -0.84 $\pm$0.02 to $\alpha_2$ = -0.48$\pm$0.03, 0.04 days after the GRB. The optical lightcurves are characterized by shallow decays-- as was also observed in the X-rays-- which may have a similar origin, related to energy injection. However, our observations indicate that there is still a puzzle concerning the chromatic breaks in the R band lightcurve (at 0.04 days) and the X-ray lightcurve (at 0.004 days) that remains to be solved.

We propose a new method of analysing a variable component for type 1 active galactic nuclei (AGNs) in the near-infrared wavelength region. This analysis uses a multiple regression technique and divides the variable component into two components originating in the accretion disk at the center of AGNs and from the dust torus that far surrounds the disk. Applying this analysis to the long-term $VHK$ monitoring data of MCG+08-11-011 that were obtained by the MAGNUM project, we found that the $(H-K)$-color temperature of the dust component is $T = 1635$K $\pm20$K, which agrees with the sublimation temperature of dust grains, and that the time delay of $K$ to $H$ variations is $\Delta t\approx 6$ days, which indicates the existence of a radial temperature gradient in the dust torus. As for the disk component, we found that the power-law spectrum of $f_\nu \propto \nu^\alpha$ in the $V$ to near-infrared $HK$ bands varies with a fixed index of $\alpha\approx -0.1$ -- +0.4, which is broadly consistent with the irradiated standard disk model. The outer part of the disk therefore extends out to a radial distance where the temperature decreases to radiate the light in the near-infrared.

We carried out a one-night optical V and near-infrared JHK monitoring observation of the least luminous Seyfert 1 galaxy, NGC4395, on 2004 May 1, and detected for the first time the intraday flux variations in the J and H bands, while such variation was not clearly seen for the K band. The detected J and H variations are synchronized with the flux variation in the V band, which indicates that the intraday-variable component of near-infrared continuum emission of the NGC4395 nucleus is an extension of power-law continuum emission to the near-infrared and originates in an outer region of the central accretion disk. On the other hand, from our regular program of long-term optical BVI and near-infrared JHK monitoring observation of NGC4395 from 2004 February 12 until 2005 January 22, we found large flux variations in all the bands on time scales of days to months. The optical BVI variations are almost synchronized with each other, but not completely with the near-infrared JHK variations. The color temperature of the near-infrared variable component is estimated to be T=1320-1710 K, in agreement with thermal emission from hot dust tori in active galactic nuclei (AGNs). We therefore conclude that the near-infrared variation consists of two components having different time scales, so that a small K-flux variation on a time scale of a few hours would possibly be veiled by large variation of thermal dust emission on a time scale of days.

Late-time BVRIJHK photometry of the peculiar Type Ic SN 2002ap, taken between 2002 June 12 and 2003 August 29 with the MAGNUM telescope, is presented. The light curve decline rate is derived in each band and the color evolution is studied through comparison with nebular spectra and with SN 1998bw. Using the photometry, the OIR bolometric light curve is built, extending from before light maximum to day 580 after explosion. The light curve has a late-time shape strikingly similar to that of the hypernova SN 1998bw. The decline rate changes from 0.018 mag/day between day 130 and 230 to 0.014 mag/day between day 270 and 580. To reproduce the late-time light curve, a dense core must be added to the 1-D hypernova model that best fits the early-time observations, bringing the ejecta mass from 2.5 Msun to 3 Msun without much change in the kinetic energy, which is 4 times 10^51 ergs. This is similar to the case of other hypernovae and suggests asymmetry. A large H-band bump developed in the spectral energy distribution after about day 300, probably caused by strong [Si I] 1.646 micron and 1.608 micron emissions. The near-infrared flux contribution increased simultaneously from <30% to >50% at day 580. The near-infrared light curves were compared with those of other Type Ib/c supernovae, among which SN 1983I seems similar to SN 2002ap both in the near-infrared and in the optical.

The most intense monitoring observations yet made in the optical (UBV) and near-infrared (JHK) wave bands were carried out for nearby Seyfert1 galaxies of NGC 5548, NGC 4051, NGC 3227, and NGC 7469. Over three years of observations with MAGNUM telescope since early 2001, clear time-delayed response of the K-band flux variations to the V-band flux variations was detected for all of these galaxies. Their H-K color temperature was estimated to be 1500-1800 K from the observed flux variation gradients, which supports a view that the bulk of the K flux should originate in the thermal radiation of hot dust that surrounds the central engine. Cross-correlation analysis was performed to quantify the lag time corresponding to the light-travel distance of the hot dust region from the central engine. The measured lag time is 47-53 days for NGC 5548, 11-18 days for NGC 4051, about 20 days for NGC 3227, and 65-87 days for NGC 7469. We found that the lag time is tightly correlated with the optical luminosity as expected from dust reverberation ($\Delta t \propto L^{0.5}$), while only weakly with the central virial mass, which suggests that an inner radius of the dust torus around the active nucleus has a one-to-one correspondence to central luminosity. In the lag time versus central luminosity diagram, the K-band lag times place an upper boundary on the similar lag times of broad-emission lines in the literature. This not only supports the unified scheme of AGNs, but also implies a physical transition from the BLR out to the dust torus that encircles the BLR. Furthermore, our V-band flux variations of NGC 5548 on timescales of up to 10 days are found to correlate with X-ray variations and delay behind them by one or two days, indicating the thermal reprocessing of X-ray emission by the central accretion flow.

We have established a GRB follow-up observation system at Kiso observatory (Japan) in 2001. Since the east Asian area had been blank for the GRB follow-up observational network, this observational system is very important in studying the temporal and spectral evolution of early afterglows. Using this system, we have performed quick observations for optical afterglows from early phase based on HETE-2 and INTEGRAL alerts. Thanks to the quick follow-up observation system, we have been able to use the Kiso observatory in 20 events, and conduct their follow-up observations in optical and near infrared wavelengths.

We observed the optical afterglow of the bright gamma-ray burst GRB030329 on the nights of 2003 March 29, using the Kiso observatory (the University of Tokyo) 1.05 m Schmidt telescope. Data were taken from March 29 13:21:26 UT to 17:43:16 (0.072 to 0.253 days after the burst), using an $Rc$-band filter. The obtained $Rc$-band light curve has been fitted successfully by a single power law function with decay index of $0.891\pm0.004$. These results remain unchanged when incorporating two early photometric data points at 0.065 and 0.073 days, reported by Price et al.(2003) using the SSO 40 inch telescope, and further including RTT150 data (Burenin et al. 2003) covering at about 0.3 days. Over the period of 0.065-0.285 days after the burst, any deviation from the power-law decay is smaller than $\pm$0.007 mag. The temporal structure reported by Uemura et al. (2003) does not show up in our $R$-band light curve.

The most intense monitoring observations yet made were carried out on the Seyfert 1 galaxy NGC 4151 in the optical and near-infrared wave-bands. A lag from the optical light curve to the near-infrared light curve was measured. The lag-time between the V and K light curves at the flux minimum in 2001 was precisely 48+2-3 days, as determined by a cross-correlation analysis. The correlation between the optical luminosity of an active galactic nucleus (AGN) and the lag-time between the UV/optical and the near-infrared light curves is presented for NGC 4151 in combination with previous lag-time measurements of NGC 4151 and other AGNs in the literature. This correlation is interpreted as thermal dust reverberation in an AGN, where the near-infrared emission from an AGN is expected to be the thermal re-radiation from hot dust surrounding the central engine at a radius where the temperature equals to that of the dust sublimation temperature. We find that the inner radius of the dust torus in NGC 4151 is $\sim $ 0.04 pc corresponding to the measured lag-time, well outside the broad line region (BLR) determined by other reverberation studies of the emission lines.

Observations were made of the optical afterglow of GRB020813 (Fox, Blake & Price, 2002) with the KISO observatory 1.05 m Schmidt telescope and the Bisei astronomical observatory 1.01 m telescope. Four-band ($B, V, R$, and $I$) photometric data points were obtained from 2002, August 13 10:52 to 16:46 UT, or 0.346$-$0.516 days after the burst. In order to investigate the early-time ($<$1 day) evolution of the afterglow, four-band light curves were produced by analyzing the data taken at these two astronomical observatories, as well as publicly released data taken by the Magellan Baade telescope (Gladders and Hall, 2002c). The light curves can be approximated by a broken power law, of which the indices are approximately 0.46 and 1.33 before and after a break at $\sim$0.2 days, respectively. The optical spectral index stayed approximately constant at $\sim$0.9 over 0.17 $\sim$ 4.07 days after the burst. Since the temporal decay index after the break and the spectral index measured at that time are both consistent with those predicted by a spherical expansion model, the early break is unlikely to be a jet break, but likely to represent the end of an early bump in the light curve as was observed in the optical afterglow of GRB021004.

Supernova (SN) 2002ap in M74 was observed in the $UBVRIJHK$ bands for the first 40 days following its discovery (2002 January 29) until it disappeared because of solar conjunction, and then in June after it reappeared. The magnitudes and dates of peak brightness in each band were determined. While the rate of increase of the brightness before the peak is almost independent of wavelength, the subsequent rate of decrease becomes smaller with wavelength from the $U$ to the $R$ band, and is constant at wavelengths beyond $I$. The photometric evolution is faster than in the well-known ``hypernovae'' SNe~1998bw and 1997ef, indicating that SN 2002ap ejected less mass. The bolometric light curve of SN 2002ap for the full period of observations was constructed. The absolute magnitude is found to be much fainter than that of SN 1998bw, but is similar to that of SN 1997ef, which lies at the faint end of the hypernova population. The bolometric light curve at the early epochs was best reproduced with the explosion of a C+O star that ejects $2.5~M_\sun$ with kinetic energy $E_{\rm K}=4\times 10^{51}~{\rm ergs}$. A comparison of the predicted brightness of SN 2002ap with that observed after solar conjunction may imply that $\gamma$-ray deposition at the later epochs was more efficient than in the model. This may be due to an asymmetric explosion.

By image processing and color excess analyses of B, V, R, and I-band CCD images of the central region of M31 taken with the Kiso 105-cm Schmidt telescope, we found a bar of 200 pc length. We also found a more extended ``face-on" spiral feature of dark clouds, which appear t o be connected to this bar and is probably an out-of-plane structure. The \co-line emission was detected from the dark clouds using the Nobeyama 45m telescope, and the virial mass of individual clouds was estimated to be about 10E6 solar masses. The CO intensity-to-virial mass ratio is anomalously small, and so is the color excess-to-mass ratio, which indicates a low temperature and a small amount dust per unit gas mass in the cloud, respectively. This implies that the molecular clouds in the central few hundred pc of M31 are more primeval, suggesting a low star-formation and low dust supply. We propose a possible scenario for the origin of the face-on spiral containing such gas on the basis of a ram-pressure accretion model of stripped gas clouds from the companion galaxies. Key words: Dark clouds -- CO line -- Galaxies; individual (M31) -- Interstellar dust -- Nuclei of galaxies.