We use MMT spectroscopy and deep Subaru Hyper Suprime-Cam (HSC) imaging to compare the spectroscopic central stellar velocity dispersion of quiescent galaxies with the effective dispersion of the dark matter halo derived from the stacked lensing signal. The spectroscopic survey (the Smithsonian Hectospec Lensing Survey) provides a sample of 4585 quiescent galaxy lenses with measured line-of-sight central stellar velocity dispersion ($\sigma_{\rm SHELS}$) that is more than 85% complete for $R < 20.6$, $D_{n}4000> 1.5$ and $M_* > 10^{9.5}{\rm M}_{\odot}$. The median redshift of the sample of lenses is 0.32. We measure the stacked lensing signal from the HSC deep imaging. The central stellar velocity dispersion is directly proportional to the velocity dispersion derived from the lensing $\sigma_{\rm Lens}$, $\sigma_{\rm Lens} = (0.98\pm 0.14)\sigma_{\rm SHELS}+(19.89\pm31.89)$. The independent spectroscopic and weak lensing velocity dispersions probe different scales, $\sim3$kpc and $\gtrsim$ 100 kpc, respectively, and strongly indicate that the observable central stellar velocity dispersion for quiescent galaxies is a good proxy for the velocity dispersion of the dark matter halo. We thus demonstrate the power of combining high-quality imaging and spectroscopy to shed light on the connection between galaxies and their dark matter halos.

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.

Hyper Suprime-Cam (HSC) is a wide-field imaging camera on the prime focus of the 8.2m Subaru telescope on the summit of Maunakea in Hawaii. A team of scientists from Japan, Taiwan and Princeton University is using HSC to carry out a 300-night multi-band imaging survey of the high-latitude sky. The survey includes three layers: the Wide layer will cover 1400 deg$^2$ in five broad bands ($grizy$), with a $5\,\sigma$ point-source depth of $r \approx 26$. The Deep layer covers a total of 26~deg$^2$ in four fields, going roughly a magnitude fainter, while the UltraDeep layer goes almost a magnitude fainter still in two pointings of HSC (a total of 3.5 deg$^2$). Here we describe the instrument, the science goals of the survey, and the survey strategy and data processing. This paper serves as an introduction to a special issue of the Publications of the Astronomical Society of Japan, which includes a large number of technical and scientific papers describing results from the early phases of this survey.

The first detected gravitational wave from a neutron star merger was GW170817. In this study, we present J-GEM follow-up observations of SSS17a, an electromagnetic counterpart of GW170817. SSS17a shows a 2.5-mag decline in the $z$-band from 1.7 days to 7.7 days after the merger. Such a rapid decline is not comparable with supernovae light curves at any epoch. The color of SSS17a also evolves rapidly and becomes redder for later epochs; the $z-H$ color changed by approximately 2.5 mag in the period of 0.7 days to 7.7 days. The rapid evolution of both the optical brightness and the color are consistent with the expected properties of a kilonova that is powered by the radioactive decay of newly synthesized $r$-process nuclei. Kilonova models with Lanthanide elements can reproduce the aforementioned observed properties well, which suggests that $r$-process nucleosynthesis beyond the second peak takes place in SSS17a. However, the absolute magnitude of SSS17a is brighter than the expected brightness of the kilonova models with the ejecta mass of 0.01 $\Msun$, which suggests a more intense mass ejection ($\sim 0.03 \Msun$) or possibly an additional energy source.

The Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) is a three-layered imaging survey aimed at addressing some of the most outstanding questions in astronomy today, including the nature of dark matter and dark energy. The survey has been awarded 300 nights of observing time at the Subaru Telescope and it started in March 2014. This paper presents the first public data release of HSC-SSP. This release includes data taken in the first 1.7 years of observations (61.5 nights) and each of the Wide, Deep, and UltraDeep layers covers about 108, 26, and 4 square degrees down to depths of i~26.4, ~26.5, and ~27.0 mag, respectively (5sigma for point sources). All the layers are observed in five broad bands (grizy), and the Deep and UltraDeep layers are observed in narrow bands as well. We achieve an impressive image quality of 0.6 arcsec in the i-band in the Wide layer. We show that we achieve 1-2 per cent PSF photometry (rms) both internally and externally (against Pan-STARRS1), and ~10 mas and 40 mas internal and external astrometric accuracy, respectively. Both the calibrated images and catalogs are made available to the community through dedicated user interfaces and database servers. In addition to the pipeline products, we also provide value-added products such as photometric redshifts and a collection of public spectroscopic redshifts. Detailed descriptions of all the data can be found online. The data release website is https://hsc-release.mtk.nao.ac.jp/.

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.

We report the results of optical--infrared follow-up observations of the gravitational wave (GW) event GW151226 detected by the Advanced LIGO in the framework of J-GEM (Japanese collaboration for Gravitational wave ElectroMagnetic follow-up). We performed wide-field optical imaging surveys with Kiso Wide Field Camera (KWFC), Hyper Suprime-Cam (HSC), and MOA-cam3. The KWFC survey started at 2.26 days after the GW event and covered 778 deg$^2$ centered at the high Galactic region of the skymap of GW151226. We started the HSC follow-up observations from 12 days after the event and covered an area of 63.5 deg$^2$ of the highest probability region of the northern sky with the limiting magnitudes of 24.6 and 23.8 for i band and z band, respectively. MOA-cam3 covered 145 deg$^2$ of the skymap with MOA-red filter 2.5 months after the GW alert. Total area covered by the wide-field surveys was 986.5 deg$^2$. The integrated detection probability of all the observed area was $\sim$29%. We also performed galaxy-targeted observations with six optical and near-infrared telescopes from 1.61 days after the event. Total of 238 nearby (<100 Mpc) galaxies were observed with the typical I band limiting magnitude of $\sim$19.5. We detected 13 supernova candidates with the KWFC survey, and 60 extragalactic transients with the HSC survey. Two third of the HSC transients were likely supernovae and the remaining one third were possible active galactic nuclei. With our observational campaign, we found no transients that are likely to be associated with GW151226.

We describe a weak lensing view of the downsizing of star forming galaxies based on cross correlating a weak lensing ($\kappa$) map with a predicted map constructed from a redshift survey. Moderately deep and high resolution images with Subaru/Hyper Suprime-Cam covering the 4 deg^2 DLS F2 field provide a $\kappa$ map with 1 arcmin resolution. A dense complete redshift survey of the F2 field including 12,705 galaxies with $R\leq20.6$ is the basis for construction of the predicted map. The zero-lag cross-correlation between the \kappa and predicted maps is significant at the $30\sigma$ level. The width of the cross-correlation peak is comparable with the angular scale of rich cluster at $z\sim0.3$, the median depth of the redshift survey. Slices of the predicted map in $\delta{z} = 0.05$ redshift bins enable exploration of the impact of structure as a function of redshift. The zero-lag normalised cross-correlation has significant local maxima at redshifts coinciding with known massive X-ray clusters. Even in slices where there are no known massive clusters, there is significant signal in the cross-correlation originating from lower mass groups that trace the large-scale of the universe. Spectroscopic $D_n4000$ measurements enable division of the sample into star-forming and quiescent populations. The significance of the cross-correlation with structure containing star-forming galaxies increases with redshift from $5\sigma$ at $z = 0.3$ to $7 \sigma$ at $z = 0.5$. The weak lensing results are consistent with the downsizing view of galaxy evolution established on the basis of many other independent studies.

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.

A large reverberation mapping study of the Seyfert 1 galaxy NGC 7469 has yielded emission-line lags for Hbeta 4861 and He II 4686 and a central black hole mass measurement of about 10 million solar masses, consistent with previous measurements. A very low level of variability during the monitoring campaign precluded meeting our original goal of recovering velocity-delay maps from the data, but with the new Hbeta measurement, NGC 7469 is no longer an outlier in the relationship between the size of the Hbeta-emitting broad-line region and the AGN luminosity. It was necessary to detrend the continuum and Hbeta and He II 4686 line light curves and those from archival UV data for different time-series analysis methods to yield consistent results.

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 have pioneered a new method for the measurement of extragalactic distances. This method uses the time-lag between variations in the short wavelength and long wavelength light from an active galactic nucleus (AGN), based on a quantitative physical model of dust reverberation that relates the time-lag to the absolute luminosity of the AGN. We use the large homogeneous data set from intensive monitoring observations in optical and near-infrared wavelength bands with the dedicated 2-m MAGNUM telescope to obtain the distances to 17 AGNs in the redshift range z=0.0024 to z=0.0353. These distance measurements are compared with distances measured using Cepheid variable stars, and are used to infer that H_0= 73 +- 3 (random) km/s/Mpc. The systematic error in H_0 is examined, and the uncertainty in the size distribution of dust grains is the largest source of the systematic error, which is much reduced for a sample of AGNs for which their parameter values in the model of dust reverberation are individually measured. This AGN time-lag method can be used beyond 30 Mpc, the farthest distance reached by extragalactic Cepheids, and can be extended to high-redshift quasi-stellar objects.

We investigate the UV-optical (longward of Ly$\alpha$ 1216\AA) spectral variability of nearly 9000 quasars ($0<z<4$) using multi-epoch photometric data within the SDSS Stripe 82 region. The regression slope in the flux-flux space of a quasar light curve directly measures the color of the flux difference spectrum, then the spectral shape of the flux difference spectra can be derived by taking a careful look at the redshift dependence of the regression slopes. First, we confirm that the observed quasar spectrum becomes bluer when the quasar becomes brighter. We infer the spectral index of the composite difference spectrum as $\alpha_{\nu}^{\text{dif}}\sim +1/3$ (in the form of $f_{\nu}\propto \nu^{\alpha_{\nu}}$), which is significantly bluer than that of the composite spectrum $\alpha_{\nu}^{\text{com}}\sim -0.5$. We also show that the continuum variability cannot be explained by the accretion disk models with varying mass accretion rate. Second, we examine the effects of broad emission line variability on the color-redshift space. The variability of the "Small Blue Bump" is extensively discussed. We show that the low-ionization lines of MgII and FeII are less variable compared to Balmer emission lines and high-ionization lines, and the Balmer continuum is the dominant variable source around $\sim 3000$\AA. These results are compared with previous studies, and the physical mechanisms of the variability of the continuum and emission lines are discussed.

In order to test a recent hypothesis that the dispersion in the Schmidt-Kennicutt law arises from variations in the evolutionary stage of star forming molecular clouds, we compared molecular gas and recent star formation in an early-phase merger galaxy pair, Taffy I (UGC\ 12915/UGC\ 12914, VV\ 254) which went through a direct collision 20 Myr ago and whose star forming regions are expected to have similar ages. Narrow-band Pa-alpha image is obtained using the ANIR near-infrared camera on the mini-TAO 1m telescope. The image enables us to derive accurate star formation rates within the galaxy directly. The total star formation rate, 22.2 M_sun/yr, was found to be much higher than previous estimates. Ages of individual star forming blobs estimated from equivalent widths indicate that most star forming regions are ~7 Myr old, except for a giant HII region at the bridge which is much younger. Comparison between star formation rates and molecular gas masses for the regions with the same age exhibits a surprisingly tight correlation, a slope of unity, and star formation efficiencies comparable to those of starburst galaxies. These results suggest that Taffy I has just evolved into a starburst system after the collision, and the star forming sites are at a similar stage in their evolution from natal molecular clouds except for the bridge region. The tight Schmidt-Kennicutt law supports the scenario that dispersion in the star formation law is in large part due to differences in evolutionary stage of star forming regions.

We present the results from a detailed analysis of photometric and spectrophotometric data on five Seyfert 1 galaxies observed as a part of a recent reverberation mapping program. The data were collected at several observatories over a 140-day span beginning in 2010 August and ending in 2011 January. We obtained high sampling-rate light curves for Mrk 335, Mrk 1501, 3C120, Mrk 6, and PG2130+099, from which we have measured the time lag between variations in the 5100 Angstrom continuum and the H-beta broad emission line. We then used these measurements to calculate the mass of the supermassive black hole at the center of each of these galaxies. Our new measurements substantially improve previous measurements of MBH and the size of the broad line-emitting region for four sources and add a measurement for one new object. Our new measurements are consistent with photoionization physics regulating the location of the broad line region in active galactic nuclei.

We examine whether the spectral energy distribution of UV continuum emission of active galactic nuclei changes during flux variation. We used multi-epoch photometric data of QSOs in the Stripe 82 observed by the SDSS Legacy Survey and selected 10 bright QSOs observed with high photometric accuracies, in the redshift range of z = 1.0-2.4 where strong broad emission lines such as Ly\alpha and CIV do not contaminate SDSS filters, to examine spectral variation of the UV continuum emission with broadband photometries. All target QSOs showed clear flux variations during the monitoring period 1998-2007, and the multi-epoch flux data in two different bands obtained on the same night showed a linear flux-to-flux relationship for all target QSOs. Assigning the flux in the longer wavelength to the x-axis in the flux-to-flux diagram, the x-intercept of the best-fit linear regression line was positive for most targets, which means that their colors in the observing bands become bluer as they become brighter. Then, the host-galaxy flux was estimated on the basis of the correlation between the stellar mass of the bulge of the host galaxy and the central black hole mass. We found that the longer-wavelength flux of the host galaxy was systematically smaller than that of the fainter extension of the best-fit regression line at the same shorter-wavelength flux for most targets. This result strongly indicates that the spectral shape of the continuum emission of QSOs in the UV region usually becomes bluer as it becomes brighter. We found that the multi-epoch flux-to-flux plots could be fitted well with the standard accretion disk model changing the mass accretion rate with a constant black hole mass for most targets. This finding strongly supports the standard accretion disk model for UV continuum emission of QSOs.

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.

We present the results of concurrent X-ray and optical monitoring of the Seyfert 1 galaxy Mrk 79 over a period of more than five years. We find that on short to medium time-scales (days to a few tens of days) the 2-10 keV X-ray and optical u and V band fluxes are significantly correlated, with a delay between the bands consistent with zero days. We show that most of these variations may be well reproduced by a model where the short-term optical variations originate from reprocessing of X-rays by an optically thick accretion disc. The optical light curves, however, also display long time-scale variations over thousands of days, which are not present in the X-ray light curve. These optical variations must originate from an independent variability mechanism and we show that they can be produced by variations in the (geometrically) thin disc accretion rate as well as by varying reprocessed fractions through changes in the location of the X-ray corona.

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.

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.

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.