We present the stellar population and ionized-gas outflow properties of ultra-luminous IR galaxies (ULIRGs) at $z=$ 0.1-1.0, which are selected from AKARI FIR all-sky survey. We construct a catalog of 1077 ULIRGs to examine feedback effect after major mergers. 202 out of the 1077 ULIRGs are spectroscopically identified by SDSS and Subaru/FOCAS observations. Thanks to deeper depth and higher resolution of AKARI compared to the previous IRAS survey, and reliable identification from WISE MIR pointing, the sample is unique in identifying optically-faint (i$\sim$20) IR-bright galaxies, which could be missed in previous surveys. A self-consistent spectrum-SED decomposition method, which constrains stellar population properties in SED modeling based on spectral fitting results, has been employed for 149 ULIRGs whose optical continua are dominated by host galaxies. They are massive galaxies ($M_{\rm star}\sim10^{11}$-$10^{12}$ M$_{\odot}$), associated with intense star formation activities (SFR $\sim$ 200-2000 M$_{\odot}$ yr$^{-1}$). The sample covers a range of AGN bolometric luminosity of $10^{10}$-$10^{13}$ L$_{\odot}$, and the outflow velocity measured from [OIII] 5007A line shows a correlation with AGN luminosity. Eight galaxies show extremely fast outflows with velocity up to 1500-2000 km s$^{-1}$. However, the co-existence of vigorous starbursts and strong outflows suggests the star formation has not been quenched during the ULIRG phase. By deriving stellar mass and mass fraction of young stellar population, we find no significant discrepancies between stellar properties of ULIRGs with weak and powerful AGNs. The results are not consistent with the merger-induced evolutionary scenario, which predicts that SF-dominated ULIRGs show smaller stellar mass and younger stellar populations compared to AGN-dominated ULIRGs.

The dust reverberation mapping is one of powerful methods to investigate the structure of the dusty tori in AGNs, and it has been performed on more than a hundred type 1 AGNs. However, no clear results have been reported on type 2 AGNs because their strong optical-UV extinction completely hides their accretion disc emission. Here we focus on an X-ray-bright type 2 AGN, NGC 2110, and utilize 2-20 keV X-ray variation monitored by MAXI to trace disc emission, instead of optical-UV variation. Comparing it with light curves in the WISE infrared (IR) W1 band ($\lambda = 3.4$ $\mu$m) and W2 band ($\lambda = 4.6$ $\mu$m) with cross-correlation analyses, we found candidates of the dust reverberation time lag at $\sim60$ days, $\sim130$ days, and $\sim1250$ days between the X-ray flux variation and those of the IR bands. By examining the best-fitting X-ray and IR light curves with the derived time lags, we found that the time lag of $\sim130$ days is most favoured. With this time lag, the relation between the time lag and luminosity of NGC 2110 is consistent with those in type 1 AGNs, suggesting that the dust reverberation in NGC 2110 mainly originates in hot dust in the torus innermost region, the same as in type 1 AGNs. As demonstrated by the present study, X-ray and IR simultaneous monitoring can be a promising tool to perform the dust reverberation mapping on type 2 AGNs.

We study variability of active galactic nuclei (AGNs) by using the deep optical multiband photometry data obtained from the Hyper Suprime-Cam Subaru Strategic Program (HSC SSP) survey in the COSMOS field. The images analyzed here were taken with 8, 10, 13, and 15 epochs over three years in the $g$, $r$, $i$, and $z$ bands, respectively. We identified 491 robust variable AGN candidates, down to $i=25$ mag and with redshift up to $4.26$. Ninety percent of the variability-selected AGNs are individually identified with the X-ray sources detected in the Chandra COSMOS Legacy survey. We investigate their properties in variability by using structure function analysis and find that the structure function for low-luminosity AGNs ($L_{\mathrm{bol}}\lesssim10^{45}$ erg s$^{-1}$) shows a positive correlation with luminosity, which is the opposite trend for the luminous quasars. This trend is likely to be caused by larger contribution of the host galaxy light for lower-luminosity AGNs. Using the model templates of galaxy spectra, we evaluate the amount of host galaxy contribution to the structure function analysis and find that dominance of the young stellar population is needed to explain the observed luminosity dependence. This suggests that low-luminosity AGNs at $0.8\lesssim z\lesssim1.8$ are predominantly hosted in star-forming galaxies. The X-ray stacking analysis reveals the significant emission from the individually X-ray undetected AGNs in our variability-selected sample. The stacked samples show very large hardness ratios in their stacked X-ray spectrum, which suggests that these optically variable sources have large soft X-ray absorption by dust-free gas.

Rapidly evolving transients form a new class of transients which show shorter timescales of the light curves than those of typical core-collapse and thermonuclear supernovae. We performed a systematic search for rapidly evolving transients using the deep data taken with the Hyper Suprime-Cam Subaru Strategic Program Transient Survey. By measuring the timescales of the light curves of 1824 transients, we identified 5 rapidly evolving transients. Our samples are found in a wide range of redshifts (0.3 $\le$ z $\le$ 1.5) and peak absolute magnitudes ($-$17 $\ge$ $M_i$ $\ge$ $-$20). The properties of the light curves are similar to those of the previously discovered rapidly evolving transients. They show a relatively blue spectral energy distribution, with the best-fit blackbody of 8,000 - 18,000 K. We show that some of the transients require power sources other than the radioactive decays of $^{56}$Ni because of their high peak luminosities and short timescales. The host galaxies of all the samples are star-forming galaxies, suggesting a massive star origin for the rapidly evolving transients. The event rate is roughly estimated to be $\sim$4,000 events yr$^{-1}$ Gpc$^{-3}$, which is about 1 $\%$ of core-collapse supernovae.

Observed extinction curves of active galactic nuclei (AGNs) are significantly different from those observed in the Milky Way. The observations require preferential removal of small grains at the AGN environment; however, the physics for this remains unclear. In this paper, we propose that dust destruction by charging, or Coulomb explosion, may be responsible for AGN extinction curves. Harsh AGN radiation makes a dust grain highly charged through photoelectric emission, and grain fission via the Coulomb explosion occurs when the electrostatic tensile stress of a charge grain exceeds its tensile strength. We show that the Coulomb explosion can preferentially remove both small silicate and graphite grains and successfully reproduce both flat extinction curves and the absence of 2175\AA~bump.

We examine the long-term optical/near-infrared (NIR) flux variability of a "changing-look" active galactic nucleus (AGN) Mrk 590 between 1998 and 2007. Multi-band multi-epoch optical/NIR photometry data from the SDSS Stripe 82 database and the Multicolor Active Galactic Nuclei Monitoring (MAGNUM) project reveal that Mrk 590 experienced a sudden luminosity decrease during the period from 2000 to 2001. Detection of dust reverberation lag signals between $V$- and $K$-band light curves obtained by the MAGNUM project during the faint state in $2003-2007$ suggests that the dust torus innermost radius $R_\text{dust}$ of Mrk 590 had become very small [$R_\text{dust} \simeq 32$~ light-days (lt-days)] by the year 2004 according to the aforementioned significant decrease in AGN luminosity. The $R_\text{dust}$ in the faint state is comparable to the H$\beta$ broad line region (BLR) radius of $R_{\text{H}\beta, \text{BLR}} \simeq 26$ lt-days measured by previous reverberation mapping observations during the bright state of Mrk 590 in $1990-1996$. These observations indicate that the innermost radius of the dust torus in Mrk 590 decreased rapidly after the AGN ultraviolet-optical luminosity drop, and that the replenishment time scale of the innermost dust distribution is less than 4 years, which is much shorter than the free fall time scale of BLR gas or dust clouds. We suggest that rapid replenishment of the innermost dust distribution can be accomplished either by new dust formation in radiatively-cooled BLR gas clouds or by new dust formation in the disk atmosphere and subsequent vertical wind from the dusty disk as a result of radiation pressure.

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.

In order to construct a sample of ultra-luminous infrared galaxies (ULIRGs, with infrared luminosity, $L_{\rm IR} > 10^{12}$ L$_{\odot}$) at 0.5 < z < 1, we are conducting an optical follow-up program for bright 90-$\mu$m FIR sources with a faint optical (i < 20) counterpart selected in the AKARI Far-Infrared Surveyor (FIS) Bright Source catalog (Ver.2). AKARI-FIS-V2 J0916248+073034, identified as a ULIRG at z = 0.49 in the spectroscopic follow-up observation, indicates signatures of an extremely strong outflow in its emission line profiles. Its [OIII] 5007\AA\ emission line shows FWHM of 1830 km s$^{-1}$ and velocity shift of -770 km s$^{-1}$ in relative to the stellar absorption lines. Furthermore, low-ionization [OII] 3726\AA\ 3729\AA\ doublet also shows large FWHM of 910 km s$^{-1}$ and velocity shift of -380 km s$^{-1}$. After the removal of an unresolved nuclear component, the long-slit spectroscopy 2D image possibly shows that the outflow extends to radius of 4 kpc. The mass outflow and energy ejection rates are estimated to be 500 M$_{\odot}$ yr$^{-1}$ and $4\times10^{44}$ erg s$^{-1}$, respectively, which imply that the outflow is among the most powerful ones observed in ULIRGs and QSOs at 0.3 < z < 1.6. The co-existence of the strong outflow and intense star formation (star formation rate of 990 M$_{\odot}$ yr$^{-1}$) indicates that the feedback of the strong outflow has not severely affect the star-forming region of the galaxy.

We report the discovery of an SN1988Z-like type IIn supernova KISS15s found in a low-mass star-forming galaxy at redshift z=0.038 during the course of the Kiso Supernova Survey (KISS). KISS15s shows long-duration optical continuum and emission line light curves, indicating that KISS15s is powered by a continuous interaction between the expanding ejecta and dense circumstellar medium (CSM). The H$\alpha$ emission line profile can be decomposed into four Gaussians of narrow, intermediate, blue-shifted intermediate, and broad velocity width components, with a full width at half maximum of $\lesssim 100$, $\sim 2,000$, and $\sim 14,000$ km s${}^{-1}$ for the narrow, intermediate, and broad components, respectively. The presence of the blue-shifted intermediate component, of which the line-of-sight velocity relative to the systemic velocity is about $-5,000$ km s${}^{-1}$, suggests that the ejecta-CSM interaction region has an inhomogeneous morphology and anisotropic expansion velocity. We found that KISS15s shows increasing infrared continuum emission, which can be interpreted as hot dust thermal emission of $T \sim 1,200$ K from newly formed dust in a cool, dense shell in the ejecta-CSM interaction region. The progenitor mass-loss rate, inferred from bolometric luminosity, is $\dot{M} \sim 0.4 M_{\odot} \text{yr}^{-1} (v_{w}/40 \text{km}~\text{s}^{-1})$, where $v_{w}$ is the progenitor's stellar wind velocity. This implies that the progenitor of KISS15s was a red supergiant star or a luminous blue variable that had experienced a large mass-loss in the centuries before the explosion.

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.

Two years' worth of u-, g-, r-, i-, and z-band optical light curves were obtained for the massive type 1 quasar PG 2308+098 at z=0.433 using the 1.05-m Kiso Schmidt telescope/Kiso Wide Field Camera, and inter-band time lags of the light curves were measured. Wavelength-dependent continuum reverberation lag signals of several tens of days relative to the u-band were detected at g-, r-, i-, and z-bands, where the longer wavelength bands showed larger lags. From the wavelength-dependent lags, and assuming the standard disk temperature radial profile $T \propto R_{\rm disk}^{-3/4}$ and an X-ray/far-ultraviolet reprocessing picture, a constraint on the radius of the accretion disk responsible for the rest-frame 2500 \AA\ disk continuum emission was derived as $R_{\rm disk} = 9.46^{+0.29}_{-3.12}$ light-days. The derived disk size is slightly (1.2-1.8 times) larger than the theoretical disk size of $R_{\rm disk} = 5.46$ light-days predicted from the black hole mass ($M_{\rm BH}$) and Eddington ratio estimates of PG 2308+098. This result is roughly in accordance with previous studies of lower mass active galactic nuclei (AGNs), where measured disk sizes have been found to be larger than the standard disk model predictions by a factor of $\sim 3$; however, the disk size discrepancy is more modest in PG 2308+098. By compiling literature values of the disk size constraints from continuum reverberation and gravitational microlensing observations for AGNs/quasars, we show that the $M_{\rm BH}$ dependence of $R_{\rm disk}$ is weaker than that expected from the standard disk model. These observations suggest that the standard Shakura-Sunyaev accretion disk theory has limitations in describing AGN/quasar accretion disks.

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 H$\alpha$ intensity map of the host galaxy of the repeating fast radio burst FRB 121102 at a redshift of z=0.193 obtained with the AO-assisted Kyoto 3DII optical integral-field unit mounted on the 8.2-m Subaru Telescope. We detected a compact H$\alpha$-emitting (i.e., star-forming) region in the galaxy, which has a much smaller angular size [$< 0".57$ (1.9 kpc) at full width at half maximum (FWHM)] than the extended stellar continuum emission region determined by the Gemini/GMOS z'-band image [$\simeq 1".4$ (4.6 kpc) at FWHM with ellipticity b/a=0.45]. The spatial offset between the centroid of the H$\alpha$ emission region and the position of the radio bursts is $0".08 \pm 0".02$ ($0.26 \pm 0.07$ kpc), indicating that FRB 121102 is located within the star-forming region. This close spatial association of FRB 121102 with the star-forming region is consistent with expectations from young pulsar/magnetar models for FRB 121102, and it also suggests that the observed H$\alpha$ emission region can make a major dispersion measure (DM) contribution to the host galaxy DM component of FRB 121102. Nevertheless, the largest possible value of the DM contribution from the H$\alpha$ emission region inferred from our observations still requires a significant amount of ionized baryons in intergalactic medium (the so-called `missing' baryons) as the DM source of FRB 121102, and we obtain a 90\% confidence level lower limit on the cosmic baryon density in the intergalactic medium in the low-redshift universe as $\Omega_{IGM} > 0.012$.

We examine the optical photometric and polarimetric variability of the luminous type 1 non-blazar quasar 3C 323.1 (PG 1545+210). Two optical spectro-polarimetric measurements taken during the periods 1996$-$98 and 2003 combined with a $V$-band imaging polarimetric measurement taken in 2002 reveal that (1) as noted in the literature, the polarization of 3C 323.1 is confined only to the continuum emission, that is, the emission from the broad line region is unpolarized; (2) the polarized flux spectra show evidence of a time-variable broad absorption feature in the wavelength range of the Balmer continuum and other recombination lines; (3) weak variability in the polarization position angle ($PA$) of $\sim$ 4 deg over a time-scale of 4$-$6 years is observed; and (4) the V-band total flux and the polarized flux show highly correlated variability over a time-scale of one year. Taking the above-mentioned photometric and polarimetric variability properties and the results from previous studies into consideration, we propose a geometrical model for the polarization source in 3C 323.1, in which an equatorial absorbing region and an axi-asymmetric equatorial electron-scattering region are assumed to be located between the accretion disc and the broad line region. The scattering/absorbing regions can perhaps be attributed to the accretion disc wind or flared disc surface, but further polarimetric monitoring observations for 3C~323.1 and other quasars with continuum-confined polarization are needed to probe the true physical origins of these regions.

From 2013 April to 2014 April, we performed an X-ray and optical simultaneous monitoring of the type 1.5 Seyfert galaxy NGC 3516. It employed Suzaku, and 5 Japanese ground-based telescopes, the Pirka, Kiso Schmidt, Nayuta, MITSuME, and the Kanata telescopes. The Suzaku observations were conducted seven times with various intervals ranging from days, weeks, to months, with an exposure of $\sim50$ ksec each. The optical $B$-band observations not only covered those of Suzaku almost simultaneously, but also followed the source as frequently as possible. As a result, NGC 3516 was found in its faint phase with the 2-10 keV flux of $0.21-2.70 \times 10^{-11}$ erg s$^{-1}$ cm$^{-2}$. The 2-45 keV X-ray spectra were composed of a dominant variable hard power-law continuum with a photon index of $\sim1.7$, and a non-relativistic reflection component with a prominent Fe-K$\alpha$ emission line. Producing the $B$-band light curve by differential image photometry, we found that the $B$-band flux changed by $\sim2.7 \times 10^{-11}$ erg s$^{-1}$ cm$^{-2}$, which is comparable to the X-ray variation, and detected a significant flux correlation between the hard power-law component in X-rays and the $B$-band radiation, for the first time in NGC 3516. By examining their correlation, we found that the X-ray flux preceded that of $B$ band by $2.0^{+0.7}_{-0.6}$ days ($1\sigma$ error). Although this result supports the X-ray reprocessing model, the derived lag is too large to be explained by the standard view which assumes a "lamppost"-type X-ray illuminator located near a standard accretion disk. Our results are better explained by assuming a hot accretion flow and a truncated disk.

Optical spectropolarimetry carried out by Kishimoto et al. (2004) has shown that several luminous type 1 quasars show a strong decrease of the polarized continuum flux in the rest-frame near-UV wavelengths of $\lambda<4000$\AA. In the literature, this spectral feature is interpreted as evidence of the broadened hydrogen Balmer absorption edge imprinted in the accretion disk thermal emission due to the disk atmospheric opacity effect. On the other hand, the quasar flux variability studies have shown that the variable continuum component in UV-optical spectra of quasars, which is considered to be a good indicator of the intrinsic spectral shape of the accretion disk emission, generally have significantly flat spectral shape throughout the near-UV to optical spectral range. To examine whether the disk continuum spectral shapes revealed as the polarized flux and as the variable component spectra are consistent with each other, we carry out multi-band photometric monitoring observations for a sample of four polarization-decreasing quasars of Kishimoto et al. (2004) (4C09.72, 3C323.1, Ton 202, and B2 1208+32) to derive the variable component spectra and compare the spectral shape of them with that of the polarized flux spectra. Contrary to expectation, we confirm that the two spectral components of these quasars have totally different spectral shape in that the variable component spectra are significantly bluer compared to the polarized flux spectra. This discrepancy in the spectral shape may imply either (1) the decrease of polarization degree in the rest-frame UV wavelengths is not indicating the Balmer absorption edge feature but is induced by some unknown (de)polarization mechanisms, or (2) the UV-optical flux variability is occurring preferentially at the hot inner radii of the accretion disk and thus the variable component spectra do not reflect the whole accretion disk emission.

We present a new method to effectively select objects which may be low-mass active black holes (BHs) at galaxy centers using high-cadence optical imaging data, and our first spectroscopic identification of an active 2.7x10^6 Msun BH at z=0.164. This active BH was originally selected due to its rapid optical variability, from a few hours to a day, based on Subaru Hyper Suprime-Cam~(HSC) g-band imaging data taken with 1-hour cadence. Broad and narrow H-alpha and many other emission lines are detected in our optical spectra taken with Subaru FOCAS, and the BH mass is measured via the broad H-alpha emission line width (1,880 km s^{-1}) and luminosity (4.2x10^{40} erg s^{-1}) after careful correction for the atmospheric absorption around 7,580-7,720A. We measure the Eddington ratio to be as low as 0.05, considerably smaller than those in a previous SDSS sample with similar BH mass and redshift, which indicates one of the strong potentials of our Subaru survey. The g-r color and morphology of the extended component indicate that the host galaxy is a star-forming galaxy. We also show effectiveness of our variability selection for low-mass active BHs.

The physical mechanisms of the quasar ultraviolet (UV)-optical variability are not well understood despite the long history of observations. Recently, Dexter & Agol presented a model of quasar UV-optical variability, which assumes large local temperature fluctuations in the quasar accretion discs. This inhomogeneous accretion disc model is claimed to describe not only the single-band variability amplitude, but also microlensing size constraints and the quasar composite spectral shape. In this work, we examine the validity of the inhomogeneous accretion disc model in the light of quasar UV-optical spectral variability by using five-band multi-epoch light curves for nearly 9 000 quasars in the Sloan Digital Sky Survey (SDSS) Stripe 82 region. By comparing the values of the intrinsic scatter $\sigma_{\text{int}}$ of the two-band magnitude-magnitude plots for the SDSS quasar light curves and for the simulated light curves, we show that Dexter & Agol's inhomogeneous accretion disc model cannot explain the tight inter-band correlation often observed in the SDSS quasar light curves. This result leads us to conclude that the local temperature fluctuations in the accretion discs are not the main driver of the several years' UV-optical variability of quasars, and consequently, that the assumption that the quasar accretion discs have large localized temperature fluctuations is not preferred from the viewpoint of the UV-optical spectral variability.

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