We present optical and near-infrared observations of SN~Ib~2019ehk. We show that it evolved to a Ca-rich transient according to its spectral properties and evolution in late phases. It, however, shows a few distinguishable properties from the canonical Ca-rich transients: a short-duration first peak in the light curve, high peak luminosity, and association with a star-forming environment. Indeed, some of these features are shared with iPTF14gqr and iPTF16hgs, which are candidates for a special class of core-collapse SNe (CCSNe): the so-called ultra-stripped envelope SNe, i.e., a relatively low-mass He (or C+O) star explosion in a binary as a precursor of double neutron star binaries. The estimated ejecta mass ($0.43 M_\odot$) and explosion energy ($1.7 \times 10^{50} $~erg) are consistent with this scenario. The analysis of the first peak suggests existence of dense circumstellar material in the vicinity of the progenitor, implying a CCSN origin. Based on these analyses, we suggest SN 2019ehk is another candidate for an ultra-stripped envelope SN. These ultra-stripped envelope SN candidates seem to form a subpopulation among Ca-rich transients, associated with young population. We propose that the key to distinguishing this population is the early first peak in their light curves.

We present optical observations of type Ia supernova (SN) 2019ein, starting at 2 days after the estimated explosion date. The spectra and the light curves show that SN 2019ein belongs to the High-Velocity (HV) and Bload Line groups with relatively rapid decline in the light curves (Delta m15(B) = 1.36 +- 0.02 mag) and the short rise time (15.37 +- 0.55 days). The Si II 6355 velocity, associated with a photospheric component but not with a detached high-velocity feature, reached ~ 20,000 km s-1 at 12 days before the B-band maximum. The line velocity however decreased very rapidly and smoothly toward the maximum light, where it was ~ 13,000 km s-1 as relatively low among HV SNe. This indicates that the speed of the spectral evolution of HV SNe Ia is correlated not only to the velocity at the maximum light, but also to the light curve decline rate like the case for Normal-Velocity (NV) SNe Ia. Spectral synthesis modeling shows that the outermost layer at > 17,000 km s-1 is well described by the O-Ne-C burning layer extending to at least 25,000 km s-1, and there is no unburnt carbon below 30,000 km s-1; these properties are largely consistent with the delayed detonation scenario, and are shared with the prototypical HV SN 2002bo despite the large difference in Delta m15(B). This structure is strikingly different from that derived for the well-studied NV SN 2011fe. We suggest that the relation between the mass of 56Ni (or Delta m15) and the extent of the O-Ne-C burning layer provides an important constraint on the explosion mechanism(s) of HV and NV SNe.

PTFO 8-8695b (CVSO 30b) is a young planet candidate whose host star is a $\sim$ 2.6 Myr-old T-Tauri star, and there have been continuous discussions about the nature of this system. To unveil the mystery of this system, we observed PTFO8-8695 for around five years at optical and infrared bands simultaneously using Kanata telescope at the Higashi-Hiroshima Observatory. Through our observations, we found that the reported fading event split into two: deeper but phase-shifted "dip-A" and shallower but equiphase "dip-B". These dips disappeared at different epochs, and then, dip-B reappeared. Based on the observed wavelength dependence of dip depths, a dust clump and a precessing planet are likely origins of dip-A and B, respectively. Here we propose "a precessing planet associated with a dust cloud" scenario for this system. This scenario is consistent with the reported change in the depth of fading events, and even with the reported results, which were thought to be negative evidence to the planetary hypothesis, such as the past non-detection of the Rossiter-McLaughlin effect. If this scenario is correct, this is the third case of a young (<3 Myr) planet around a pre-main sequence star. This finding implies that a planet can be formed within a few Myr.

We present the results of flux density, spectral index, and polarization intra-night monitoring studies of a sample of eight optically bright blazars, carried out by employing several small to moderate aperture (0.4\,m to 1.5\,m diameter) telescopes fitted with CCDs and polarimeters located in Europe, India, and Japan. The duty cycle of flux variability for the targets is found to be $\sim 45$ percent, similar to that reported in earlier studies. The computed two-point spectral indices are found to be between 0.65 to 1.87 for our sample, comprised of low- and intermediate frequency peaked blazars, with one exception; they are also found to be statistically variable for about half the instances where `confirmed' variability is detected in flux density. In the analysis of the spectral evolution of the targets on hourly timescale, a counter-clockwise loop (soft-lagging) is noted in the flux-spectral index plane on two occasions, and in one case a clear spectral flattening with the decreasing flux is observed. In our data set, we also observe a variety of flux-polarization degree variability patterns, including instances with a relatively straightforward anti-correlation, correlation, or counter-clockwise looping. These changes are typically reflected in the flux-polarization angle plane: the anti-correlation between the flux and polarization degree is accompanied by an anti-correlation between the polarization angle and flux, while the counter-clockwise flux-PD looping behaviour is accompanied by a clockwise looping in the flux-polarization angle representation. We discuss our findings in the framework of the internal shock scenario for blazar sources.

We followed-up the massive young stellar object (MYSO) S255-NIRS3 (=S255-IRS1b) during its recent accretion outburst event in the Ks band with Kanata/HONIR for four years after its burst and obtained a long-term light curve. This is the most complete NIR light-curve of the S255-NIRS3 burst event that has ever been presented. The light curve showed a steep increase reaching a peak flux that was 3.4 mag brighter than the quiescent phase and then a relatively moderate year-scale fading until the last observation, similar to that of the accretion burst events such as EXors found in lower-mass young stellar objects. The behavior of the Ks band light curve is similar to that observed in 6.7 GHz class II methanol maser emission, with a sudden increase followed by moderate year-scale fading. However, the maser emission peaks appear 30-50 days earlier than that of the Ks band emission. The similarities confirmed that the origins of the maser emission and the Ks band continuum emission is common as previously shown from another infrared and radio observations by Stecklum et al. (2016); Caratti o Garatti et al. (2017a); Moscadelli et al. (2017). However, the differences in energy transfer paths, such as the exciting/emitting/scattering structures, may cause the delay in the flux-peak dates.

High cadence ultraviolet, optical and near-infrared photometric and low-resolution spectroscopic observations of the peculiar Type II supernova (SN) 2018hna are presented. The early phase multiband light curves exhibit the adiabatic cooling envelope emission following the shock breakout up to ~14 days from the explosion. SN~2018hna has a rise time of $\sim$\,88 days in the V-band, similar to SN 1987A. A $\rm^{56}Ni$ mass of ~0.087$\pm$0.004 $\rm M_{\odot}$ is inferred for SN 2018hna from its bolometric light curve. Hydrodynamical modelling of the cooling phase suggests a progenitor with a radius ~50 $\rm R_{\odot}$, a mass of ~14-20 $\rm M_{\odot}$ and explosion energy of ~1.7-2.9$\rm \times$ $\rm 10^{51}\ erg$. The smaller inferred radius of the progenitor than a standard red supergiant is indicative of a blue supergiant progenitor of SN 2018hna. A sub-solar metallicity (~0.3 $\rm Z_{\odot}$) is inferred for the host galaxy UGC 07534, concurrent with the low-metallicity environments of 1987A-like events.

We present the results based on photometric ($Swift$ UVOT), broad-band polarimetric ($V$ and $R$-band) and optical spectroscopic observations of the Type IIn supernova (SN) 2017hcc. Our study is supplemented with spectropolarimetric data available in literature for this event. The post-peak light curve evolution is slow ($\sim$0.2 mag 100 d$^{-1}$ in $b$-band). The spectrum of $\sim$+27 d shows a blue continuum with narrow emission lines, typical of a Type IIn SN. Archival polarization data along with the $Gaia$ DR2 distances have been utilized to evaluate the interstellar polarization (ISP) towards the SN direction which is found to be $P_{ISP}$ = 0.17 $\pm$ 0.02 per cent and $\theta_{ISP}$ = 140$^{\circ}$ $\pm$ 3$^{\circ}$. To extract the intrinsic polarization of SN 2017hcc, both the observed and the literature polarization measurements were corrected for ISP. We noticed a significant decline of $\sim$3.5 per cent ($V$-band) in the intrinsic level of polarization spanning a period of $\sim$2 months. In contrast, the intrinsic polarization angles remain nearly constant at all epochs. Our study indicates a substantial variation in the degree of asymmetry in either the ejecta and/or the surrounding medium of SN 2017hcc. We also estimate a mass-loss rate of $\dot M$ = 0.12 M$_{\odot}$ yr$^{-1}$ (for $v_w$ = 20 km s$^{-1}$) which suggests that the progenitor of SN 2017hcc is most likely a Luminous Blue Variable.

We present optical and near-infrared observations of the rapidly evolving supernova (SN) 2017czd that shows hydrogen features. The optical light curves exhibit a short plateau phase ($\sim 13$ days in the $R$-band) followed by a rapid decline by $4.5$ mag in $\sim 20 \mathrm{days}$ after the plateau. The decline rate is larger than those of any standard SNe, and close to those of rapidly evolving transients. The peak absolute magnitude is $-16.8$ mag in the $V$-band, which is within the observed range for SNe IIP and rapidly evolving transients. The spectra of SN 2017czd clearly show the hydrogen features and resemble those of SNe IIP at first. The H$\alpha$ line, however, does not evolve much with time and it becomes similar to those in SNe IIb at decline phase. We calculate the synthetic light curves using a SN IIb progenitor which has 16 M$_{\odot}$ at the zero-age main sequence and evolves in a binary system. The model with a low explosion energy ($5\times 10^{50}$ erg) and a low ${}^{56}$Ni mass ($0.003 \mathrm{M}_{\odot}$) can reproduce the short plateau phase as well as the sudden drop of the light curve as observed in SN 2017czd. We conclude that SN 2017czd might be the first identified weak explosion from a SN IIb progenitor. We suggest that some rapidly evolving transients can be explained by such a weak explosion of the progenitors with little hydrogen-rich envelope.

We report on a recent multi-band optical photometric and polarimetric observational campaign of the blazar OJ 287 which was carried out during September 2016 -- December 2017. We employed nine telescopes in Bulgaria, China, Georgia, Japan, Serbia, Spain and the United States. We collected over 1800 photometric image frames in BVRI bands and over 100 polarimetric measurements over ~175 nights. In 11 nights with many quasi-simultaneous multi-band (V, R, I) observations, we did not detect any genuine intraday variability in flux or color. On longer timescales, multiple flaring events were seen. Large changes in color with respect to time and in a color--magnitude diagram were seen, and while only a weak systematic variability trend was noticed in color with respect to time, the color--magnitude diagram shows a bluer-when-brighter trend. Large changes in the degree of polarization, and substantial swings in the polarization angle were detected. The fractional Stokes parameters of the polarization showed a systematic trend with time in the beginning of these observations, followed by chaotic changes and then an apparently systematic variation at the end. These polarization changes coincide with the detection and duration of the source at very high energies as seen by VERITAS. The spectral index shows a systematic variation with time and V-band magnitude. We briefly discuss possible physical mechanisms that could explain the observed flux, color, polarization, and spectral variability.

We report X-ray, optical, and near-infrared monitoring of the new X-ray transient MAXI J1820$+$070 discovered with MAXI on 2018 March 11. Its X-ray intensity reached $\sim 2$ Crab in 2--20 keV at the end of March, and then gradually decreased until the middle of June. In this period, the X-ray spectrum was described by Comptonization of the disk emission, with a photon index of $\sim$1.5 and an electron temperature of $\sim$50 keV, which is consistent with a black hole X-ray binary in the low/hard state. The electron temperature and the photon index were slightly decreased and increased with increasing flux, respectively. The source showed significant X-ray flux variation on timescales of seconds. This short-term variation was found to be associated with changes in the spectral shape, and the photon index became slightly harder at higher fluxes. This suggests that the variation was produced by a change in the properties of the hot electron cloud responsible for the strong Comptonization. Modeling a multi-wavelength SED around the X-ray flux peak at the end of March, covering the near-infrared to X-ray bands, we found that the optical and near-infrared fluxes were likely contributed substantially by the jet emission. Before this outburst, the source was never detected in the X-ray band with MAXI (with a 3$\sigma$ upper limit of $\sim$0.2 mCrab in 4--10 keV, obtained from the 7-year data in 2009--2016), whereas weak optical and infrared activity was found at their flux levels $\sim$3 orders of magnitude lower than the peak fluxes in the outburst.

Spectropolarimetry is a powerful technique for investigating the physical properties of gas and solid materials in cometary comae without mutual contamination, but there have been few spectropolarimetric studies to extract each component. We attempt to derive the continuum polarization degree of comet 2P/Encke, free from influence of molecular emissions. The target is unique in that it has an orbit dynamically decoupled from Jupiter like main-belt asteroids, while ejecting gas and dust like ordinary comets. We observed the comet using the Higashi-Hiroshima Optical and Near-Infrared Camera attached to the Cassegrain focus of the 150-cm Kanata telescope on UT 2017 February 21 when the comet was at the solar phase angle of 75.7 deg. We find that the continuum polarization degree with respect to the scattering plane is 33.8+/-2.7 % at the effective wavelength of 0.815 um, which is significantly higher than those of cometary dust in a high-Pmax group at similar phase angles. Assuming that an ensemble polarimetric response of 2P/Encke's dust as a function of phase angle is morphologically similar with those of other comets, its maximum polarization degree is estimated to > 40 % at the phase angle of ~100 deg. In addition, we obtain the polarization degrees of the C2 swan bands (0.51-0.56 um), the NH2 alpha bands (0.62-0.69 um) and the CN-red system (0.78-0.94 um) in a range of 3-19 %, which depend on the molecular species and rotational quantum numbers of each branch. The polarization vector aligns nearly perpendicularly to the scattering plane with the average of 0.4 deg over a wavelength range of 0.50-0.97 um. From the observational evidence, we conjecture that the large polarization degree of 2P/Encke would be attributable to a dominance of large dust particles around the nucleus, which have remained after frequent perihelion passages near the Sun.

We present optical and near-infrared observations of the nearby Type Iax supernova (SN) 2014dt from 14 to 410 days after the maximum light. The velocities of the iron absorption lines in the early phase indicated that SN 2014dt showed slower expansion than the well-observed Type Iax SNe 2002cx, 2005hk and 2012Z. In the late phase, the evolution of the light curve and that of the spectra were considerably slower. The spectral energy distribution kept roughly the same shape after ~100 days, and the bolometric light curve flattened during the same period. These observations suggest the existence of an optically thick component that almost fully trapped the {\gamma}-ray energy from 56 Co decay. These findings are consistent with the predictions of the weak deflagration model, leaving a bound white dwarf remnant after the explosion.

We present a multi-wavelength spectral and temporal investigation of OJ 287 emission during its strong optical-to-X-ray activity between July 2016 - July 2017. The daily $\gamma$-ray fluxes from \emph{Fermi}-LAT are consistent with no variability. The strong optical-to-X-ray variability is accompanied by a change in power-law spectral index of the X-ray spectrum from $< 2$ to $>2$, with variations often associated with changes in optical polarization properties. Cross-correlations between optical-to-X-ray emission during four continuous segments show simultaneous optical-ultraviolet (UV) variations while the X-ray and UV/optical are simultaneous only during the middle two segments. In the first segment, the results suggest X-rays lag the optical/UV, while in the last segment X-rays lead by $\sim$ 5-6 days. The last segment also shows a systematic trend with variations appearing first at higher energies followed by lower energy ones. The LAT spectrum before the VHE activity is similar to preceding quiescent state spectrum while it hardens during VHE activity period and is consistent with the extrapolated VHE spectrum during the latter. Overall, the broadband spectral energy distributions (SEDs) during high activity periods are a combination of a typical OJ 287 SED and an HBL SED, and can be explained in a two-zone leptonic model, with the second zone located at parsec scales, beyond the broad line region, being responsible for the HBL-like spectrum. The change of polarization properties from systematic to chaotic and back to systematic, before, during and after the VHE activity, suggest dynamic roles for magnetic fields and turbulence.

We present optical and near-infrared observations of a low-luminosity Type IIP supernova (SN) 2016bkv from the initial rising phase to the plateau phase. Our observations show that the end of the plateau is extended to $\gtrsim 140$ days since the explosion, indicating that this SN takes one of the longest time to finish the plateau phase. among Type IIP SNe (SNe IIP), including low-luminosity (LL) SNe IIP. The line velocities of various ions at the middle of the plateau phase are as low as 1,000--1,500 km s$^{-1}$, which is the lowest even among LL SNe IIP. These measurements imply that the ejecta mass in SN 2016bkv is larger than that of the well-studied LL IIP SN 2003Z. In the early phase, SN 2016bkv shows a strong bump in the light curve. In addition, the optical spectra in this bump phase exhibit a blue continuum accompanied with a narrow H$\alpha$ emission line. These features indicate an interaction between the SN ejecta and the circumstellar matter (CSM) as in SNe IIn. Assuming the ejecta-CSM interaction scenario, the mass loss rate is estimated to be $\sim 1.7 \times 10^{-2} M_{\odot}$ yr$^{-1}$ within a few years before the SN explosion. This is comparable to or even larger than the largest mass loss rate observed for the Galactic red supergiants ($\sim 10^{-3} M_{\odot}$ yr$^{-1}$ for VY CMa). We suggest that the progenitor star of SN 2016bkv experienced a violent mass loss just before the SN explosion.

Recent detection of gravitational waves from a neutron star (NS) merger event GW170817 and identification of an electromagnetic counterpart provide a unique opportunity to study the physical processes in NS mergers. To derive properties of ejected material from the NS merger, we perform radiative transfer simulations of kilonova, optical and near-infrared emissions powered by radioactive decays of r-process nuclei synthesized in the merger. We find that the observed near-infrared emission lasting for > 10 days is explained by 0.03 Msun of ejecta containing lanthanide elements. However, the blue optical component observed at the initial phases requires an ejecta component with a relatively high electron fraction (Ye). We show that both optical and near-infrared emissions are simultaneously reproduced by the ejecta with a medium Ye of ~ 0.25. We suggest that a dominant component powering the emission is post-merger ejecta, which exhibits that mass ejection after the first dynamical ejection is quite efficient. Our results indicate that NS mergers synthesize a wide range of r-process elements and strengthen the hypothesis that NS mergers are the origin of r-process elements in the Universe.

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.

We present a multi-wavelength spectral and temporal analysis of the blazar OJ 287 during its recent activity between December 2015 -- May 2016, showing strong variability in the near-infrared (NIR) to X-ray energies with detection at $\gamma$-ray energies as well. Most of the optical flux variations exhibit strong changes in polarization angle and degree. All the inter-band time lags are consistent with simultaneous emissions. Interestingly, on days with excellent data coverage in the NIR--UV bands, the spectral energy distributions (SEDs) show signatures of bumps in the visible--UV bands, never seen before in this source. The optical bump can be explained as accretion-disk emission associated with the primary black hole of mass $\sim \rm 1.8 \times10^{10} M_{\odot}$ while the little bump feature in the optical-UV appears consistent with line emission. Further, the broadband SEDs extracted during the first flare and during a quiescent period during this span show very different $\gamma$-ray spectra compared to previously observed flare or quiescent spectra. The probable thermal bump in the visible seems to have been clearly present since May 2013, as found by examining all available NIR-optical observations, and favors the binary super-massive black hole model. The simultaneous multi-wavelength variability and relatively weak $\gamma$-ray emission that shows a shift in the SED peak is consistent with $\gamma$-ray emission originating from inverse Compton scattering of photons from the line emission that apparently contributes to the little blue bump.

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

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

V404 Cygni is a well-known black hole binary candidate thought to have relativistic jets. It showed extreme outbursts on June 2015, characterized by a large amplitude and short time variation of flux in the radio, optical, and X-ray bands. Not only disk emission, but also synchrotron radiation from the relativistic jets were suggested by radio observations. However, it is difficult to measure the accurate spectral shape in the optical/near infrared band because there are uncertainties of interstellar extinction. To estimate the extinction value for V404 Cygni, we performed photopolarimetric and spectroscopic observations of V404 Cygni and nearby field stars. Here, we estimate the Galactic extinction using interstellar polarization based on the observation that the origin of the optical polarization is the interstellar medium, and investigate the properties of interstellar polarization around V404~Cygni. We found a good correlation between the color excess and polarization degree in the field stars. We also confirmed that the wavelength dependence of the polarization degree in the highly polarized field stars was similar to that of V404~Cygni. Using the highly polarized field stars, we estimated the color excess and the extinction, $E(B-V)=1.2 \pm 0.2$ and $3.0 < A(V) < 3.6$, respectively. A tendency for a bluer peak of polarization ($\lambda_{\rm max}<5500$ \AA) was commonly seen in the highly polarized field stars, suggesting that the dust grains toward this region are generally smaller than the Galactic average. The corrected spectral energy distribution of V404~Cygni in the near infrared (NIR) and optical bands in our results indicated a spectral break between $2.5 \times 10^{14}$ Hz and $3.7 \times 10^{14}$ Hz, which might be originated in the synchrotron self absorption.

We present recent optical photometric observations of the blazar OJ 287 taken during September 2015 -- May 2016. Our intense observations of the blazar started in November 2015 and continued until May 2016 and included detection of the large optical outburst in December 2016 that was predicted using the binary black hole model for OJ 287. For our observing campaign, we used a total of 9 ground based optical telescopes of which one is in Japan, one is in India, three are in Bulgaria, one is in Serbia, one is in Georgia, and two are in the USA. These observations were carried out in 102 nights with a total of ~ 1000 image frames in BVRI bands, though the majority were in the R band. We detected a second comparably strong flare in March 2016. In addition, we investigated multi-band flux variations, colour variations, and spectral changes in the blazar on diverse timescales as they are useful in understanding the emission mechanisms. We briefly discuss the possible physical mechanisms most likely responsible for the observed flux, colour and spectral variability.

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.

Blazars are highly variable active galactic nuclei which emit radiation at all wavelengths from radio to gamma-rays. Polarized radiation from blazars is one key piece of evidence for synchrotron radiation at low energies and it also varies dramatically. The polarization of blazars is of interest for understanding the origin, confinement, and propagation of jets. However, even though numerous measurements have been performed, the mechanisms behind jet creation, composition and variability are still debated. We performed simultaneous gamma-ray and optical photopolarimetry observations of 45 blazars between Jul. 2008 and Dec. 2014 to investigate the mechanisms of variability and search for a basic relation between the several subclasses of blazars. We identify a correlation between the maximum degree of optical linear polarization and the gamma-ray luminosity or the ratio of gamma-ray to optical fluxes. Since the maximum polarization degree depends on the condition of the magnetic field (chaotic or ordered), this result implies a systematic difference in the intrinsic alignment of magnetic fields in pc-scale relativistic jets between different types blazars (FSRQs vs. BL Lacs), and consequently between different types of radio galaxies (FR Is vs. FR IIs).

We observed RZ LMi, which is renowned for the extremely (~19d) short supercycle and is a member of a small, unusual class of cataclysmic variables called ER UMa-type dwarf novae, in 2013 and 2016. In 2016, the supercycles of this object substantially lengthened in comparison to the previous measurements to 35, 32, 60d for three consecutive superoutbursts. We consider that the object virtually experienced a transition to the novalike state (permanent superhumper). This observed behavior extremely well reproduced the prediction of the thermal-tidal instability model. We detected a precursor in the 2016 superoutburst and detected growing (stage A) superhumps with a mean period of 0.0602(1)d in 2016 and in 2013. Combined with the period of superhumps immediately after the superoutburst, the mass ratio is not as small as in WZ Sge-type dwarf novae, having orbital periods similar to RZ LMi. By using least absolute shrinkage and selection operator (Lasso) two-dimensional power spectra, we detected possible negative superhumps with a period of 0.05710(1)d. We estimated the orbital period of 0.05792d, which suggests a mass ratio of 0.105(5). This relatively large mass ratio is even above ordinary SU UMa-type dwarf novae, and it is also possible that the exceptionally high mass-transfer rate in RZ LMi may be a result of a stripped core evolved secondary which are evolving toward an AM CVn-type object.

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.

We report on two superoutbursts of the AM CVn-type object CR Boo in 2014 April--March and 2015 May--June. A precursor outburst acompanied both of these superoutbursts. During the rising branch of the main superoutburst in 2014, we detected growing superhumps (stage A superhumps) whose period was $0.017669(24)$ d. Assuming that this period reflects the dynamical precession rate at the radius of the 3:1 resonance, we could estimate the mass ratio ($q=M_2/M_1$) of 0.101(4) by using the stage A superhump period and the orbital one of 0.0170290(6) d. This mass ratio is consistent with that expected by the theoretical evolutionary model of AM CVn-type objects. The detection of precursor outbursts and stage A superhumps is the second case in AM CVn-type objects. There are two interpretations of the outbursts of AM CVn-type objects. One is a dwarf nova (DN) outbursts analogy, which is caused by thermal and tidal instabilities. Another is the VY Scl-type variation, which is caused by the variation of the mass-transfer rate of the secondary. This detection of the superhump variations strongly suggests the former interpretation.

We present extensively dense observations of the super-Chandrasekhar supernova (SC SN) candidate SN 2012dn from $-11$ to $+140$ days after the date of its $B$-band maximum in the optical and near-infrared (NIR) wavelengths conducted through the OISTER ToO program. The NIR light curves and color evolutions up to 35 days after the $B$-band maximum provided an excellent match with those of another SC SN 2009dc, providing a further support to the nature of SN 2012dn as a SC SN. We found that SN 2012dn exhibited strong excesses in the NIR wavelengths from $30$ days after the $B$-band maximum. The $H$ and $K_{s}$-band light curves exhibited much later maximum dates at $40$ and $70$ days after the $B$-band maximum, respectively, compared with those of normal SNe Ia. The $H$ and $K_{s}$-band light curves subtracted by those of SN 2009dc displayed plateaued evolutions, indicating a NIR echo from the surrounding dust. The distance to the inner boundary of the dust shell is limited to be $4.8 - 6.4\times10^{-2}$ pc. No emission lines were found in its early phase spectrum, suggesting that the ejecta-CSM interaction could not occur. On the other hand, we found no signature that strongly supports the scenario of dust formation. The mass loss rate of the pre-explosion system is estimated to be $10^{-6}-10^{-5}$ M$_{\odot}$ yr$^{-1}$, assuming that the wind velocity of the system is $10-100$ km~s$^{-1}$, which suggests that the progenitor of SN 2012dn could be a recurrent nova system. We conclude that the progenitor of this SC SN could be explained by the single-degenerate scenario.

We present simultaneous optical and near-infrared (NIR) polarimetric results for the black hole binary V404 Cygni spanning the duration of its 7-day long optically-brightest phase of its 2015 June outburst. The simultaneous R and Ks-band light curves showed almost the same temporal variation except for the isolated (~30 min duration) orphan Ks-band flare observed at MJD 57193.54. We did not find any significant temporal variation of polarization degree (PD) and position angle (PA) in both R and Ks bands throughout our observations, including the duration of the orphan NIR flare. We show that the observed PD and PA are predominantly interstellar in origin by comparing the V404 Cyg polarimetric results with those of the surrounding sources within the 7'x7' field-of-view. The low intrinsic PD (less than a few percent) implies that the optical and NIR emissions are dominated by either disk or optically-thick synchrotron emission, or both. We also present the broadband spectra of V404 Cyg during the orphan NIR flare and a relatively faint and steady state by including quasi-simultaneous Swift/XRT and INTEGRAL fluxes. By adopting a single-zone synchrotron plus inverse-Compton model as widely used in modeling of blazars, we constrained the parameters of a putative jet. Because the jet synchrotron component cannot exceed the Swift/XRT disk/corona flux, the cutoff Lorentz factor in the electron energy distribution is constrained to be <10^2, suggesting particle acceleration is less efficient in this microquasar jet outburst compared to AGN jets. We also suggest that the loading of the baryon component inside the jet is inevitable based on energetic arguments.

We present broadband photometric and polarimetric observations of two type II supernovae (SNe) 2013hj and 2014G. SN 2014G is a spectroscopically classified type IIL event, which we also confirm photometrically as its light curve show characteristic features (plateau slope of 2.55 mag (100 d)$ ^{-1} $ in V-band and duration of $ \sim77 $d) of a generic IIL SN. On the other hand SN 2013hj also shows high plateau decline rate of $ 1.5 $ mag (100 d)$ ^{-1} $ in V-band, similar to SNe IIL, but marginally lower than SNe IIL template light curves. Our high cadence photometric observations of SNe 2013hj and 2014G enables us to cover all characteristic phases up to radioactive tail of optical light curves. Broadband polarimetric observations reveal some polarization in SN 2013hj with subtle enhancement as SN evolves towards plateau end, however the polarization angle remains constant throughout the evolution. This characteristic is consistent with the idea that the evolving SN with recombining hydrogen envelope is slowly revealing more asymmetric central region of explosion. Modelling of bolometric light curve yields a progenitor mass of $ \sim11 $M$_{\odot}$ with a radius of $ \sim700 $R$_{\odot}$ for SN 2013hj, while for SN 2014G model estimated progenitor mass is $ \sim9 $M$_{\odot}$ with a radius of $ \sim630 $R$_{\odot}$, both having a typical energy budget of $ \sim2\times10^{51} $ erg.

Continuing the project described by Kato et al. (2009, arXiv:0905.1757), we collected times of superhump maxima for 102 SU UMa-type dwarf novae observed mainly during the 2014-2015 season and characterized these objects. Our project has greatly improved the statistics of the distribution of orbital periods, which is a good approximation of the distribution of cataclysmic variables at the terminal evolutionary stage, and confirmed the presence of a period minimum at a period of 0.053 d and a period spike just above this period. The number density monotonically decreased toward the longer period and there was no strong indication of a period gap. We detected possible negative superhumps in Z Cha. It is possible that normal outbursts are also suppressed by the presence of a disk tilt in this system. There was no indication of enhanced orbital humps just preceding the superoutburst, and this result favors the thermal-tidal disk instability as the origin of superoutbursts. We detected superhumps in three AM CVn-type dwarf novae. Our observations and recent other detections suggest that 8% of objects showing dwarf nova-type outbursts are AM CVn-type objects. AM CVn-type objects and EI Psc-type object may be more abundant than previously recognized. OT J213806, a WZ Sge-type object, exhibited a remarkably different feature between the 2010 and 2014 superoutbursts. Although the 2014 superoutburst was much fainter the plateau phase was shorter than the 2010 one, the course of the rebrightening phase was similar. This object indicates that the O-C diagrams of superhumps can be indeed variable at least in WZ Sge-type objects. Four deeply eclipsing SU UMa-type dwarf novae (ASASSN-13cx, ASASSN-14ag, ASASSN-15bu, NSV 4618) were identified. We studied long-term trends in supercycles in MM Hya and CY UMa and found systematic variations of supercycles of ~20%.

We report the late-time evolution of Type IIb Supernova (SN IIb) 2013df. SN 2013df showed a dramatic change in its spectral features at ~1 year after the explosion. Early on it showed typical characteristics shared by SNe IIb/Ib/Ic dominated by metal emission lines, while later on it was dominated by broad and flat-topped Halpha and He I emissions. The late-time spectra are strikingly similar to SN IIb 1993J, which is the only previous example clearly showing the same transition. This late-time evolution is fully explained by a change in the energy input from the $^{56}$Co decay to the interaction between the SN ejecta and dense circumstellar matter (CSM). The mass loss rate is derived to be ~(5.4 +- 3.2) x 10^{-5} Msun/yr (for the wind velocity of ~20 km/s), similar to SN 1993J but larger than SN IIb 2011dh by an order of magnitude. The striking similarity between SNe 2013df and 1993J in the (candidate) progenitors and the CSM environments, and the contrast in these natures to SN 2011dh, infer that there is a link between the natures of the progenitor and the mass loss: SNe IIb with a more extended progenitor have experienced a much stronger mass loss in the final centuries toward the explosion. It might indicate that SNe IIb from a more extended progenitor are the explosions during a strong binary interaction phase, while those from a less extended progenitor have a delay between the strong binary interaction and the explosion.

We presented optical and near-infrared multi-band linear polarimetry of the highly reddened Type Ia SN~2014J appeared in M82. SN~2014J exhibits large polarization at shorter wavelengths, e.g., $4.8$\% in $B$ band, and the polarization decreases rapidly at longer wavelengths, with the position angle of the polarization remaining at approximately $40^{\circ}$ over the observed wavelength range. These polarimetric properties suggest that the observed polarization is likely to be caused predominantly by the interstellar dust within M82. Further analysis shows that the polarization peaks at a wavelengths much shorter than those obtained for the Galactic dust. The wavelength dependence of the polarization can be better described by an inverse power law rather than by Serkowski law for Galactic interstellar polarization. These suggests that the nature of the dust in M82 may be different from that in our Galaxy, with polarizing dust grains having a mean radius of $<0.1\ \mu$m.

We report on a superoutburst of the AM CVn-type object SDSS J090221.35+381941.9 [J0902; orbital period 0.03355(6) d] in 2014 March-April. The entire outburst consisted of a precursor outburst and the main superoutburst, followed by a short rebrightening. During the rising branch of the main superoutburst, we detected growing superhumps (stage A superhumps) with a period of 0.03409(1) d. During the plateau phase of the superoutburst, superhumps with a shorter period (stage B superhumps) were observed. Using the orbital period and the period of the stage A superhumps, we were able to measure the dynamical precession rate of the accretion disk at the 3:1 resonance, and obtained a mass ratio (q) of 0.041(7). This is the first successful measurement of the mass ratio in an AM CVn-type object using the recently developed stage A superhump method. The value is generally in good agreement with the theoretical evolutionary model. The orbital period of J0902 is the longest among the outbursting AM CVn-type objects, and the borderline between the outbursting systems and systems with stable cool disks appears to be longer than had been supposed.

Continuing the project described by Kato et al. (2009, PASJ, 61, S395, arXiv:0905.1757), we collected times of superhump maxima for 56 SU UMa-type dwarf novae mainly observed during the 2013-2014 season and characterized these objects. We detected negative superhumps in VW Hyi and indicated that the low number of normal outbursts in some supercycle can be interpreted as a result of the disk tilt. This finding, combined with the Kepler observation of V1504 Cyg and V344 Lyr, suggests that the disk tilt is responsible for modulating the outburst pattern in SU UMa-type dwarf novae. We also studied the deeply eclipsing WZ Sge-type dwarf nova MASTER OT J005740.99+443101.5 and found evidence of a sharp eclipse during the phase of early superhumps. The profile can be reproduced by a combination of the eclipse of the axisymmetric disk and the uneclipsed light source of early superhumps. This finding confirms the lack of evince of a greatly enhanced hot spot during the early stage of WZ Sge-type outburst. We detected growing (stage A) superhumps in MN Dra and give a suggestion that some of SU UMa-type dwarf novae situated near the critical condition of tidal instability may show long-lasting stage A superhumps. The large negative period derivatives reported in such systems can be understood a result of the combination of stage A and B superhumps. The WZ Sge-type dwarf novae AL Com and ASASSN-13ck showed a long-lasting (plateau-type) rebrightening. In the early phase of the rebrightening, both objects showed a precursor-like outburst, suggesting that the long-lasting rebrightening is triggered by a precursor outburst.

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 report photometric and spectroscopic observations of the nearby Type Ia Supernova (SN Ia) 2012ht from $-15.8$ days to $+49.1$ days after $B$-band maximum. The decline rate of the light curve is $\Delta m_{15}$($B$)$=1.39~\pm~0.05$ mag, which is intermediate between normal and subluminous SNe Ia, and similar to that of the `transitional' Type Ia SN 2004eo. The spectral line profiles also closely resemble those of SN 2004eo. We were able to observe SN 2012ht at very early phase, when it was still rising and was about three magnitudes fainter than at the peak. The rise time to the $B$-band maximum is estimated to be $17.6 \pm 0.5$ days and the time of the explosion is MJD $56277.98 \pm 0.13$. SN 2012ht is the first transitional SN Ia whose rise time is directly measured without using light curve templates, and the fifth SN Ia overall. This rise time is consistent with those of the other four SNe within the measurement error, even including the extremely early detection of SN 2013dy. The rising part of the light curve can be fitted by a quadratic function, and shows no sign of a shock-heating component due to the interaction of the ejecta with a companion star. The rise time is significantly longer than that inferred for subluminous SNe such as SN 1991bg, which suggests that a progenitor and/or explosion mechanism of transitional SNe Ia are more similar to normal SNe Ia rather than subluminous SNe Ia.

Continuing the project described by Kato et al. (2009a, arXiv:0905.1757), we collected times of superhump maxima for SU UMa-type dwarf novae mainly observed during the 2012-2013 season. We found three objects (V444 Peg, CSS J203937 and MASTER J212624) having strongly positive period derivatives despite the long orbital period (Porb). By using the period of growing stage (stage A) superhumps, we obtained mass ratios for six objects. We characterized nine new WZ Sge-type dwarf novae. We made a pilot survey of the decline rate of slowly fading part of SU UMa-type and WZ Sge-type outbursts. The decline time scale was found to generally follow the expected Porb^(1/4) dependence and WZ Sge-type outbursts also generally follow this trend. There are some objects which show slower decline rates, and we consider these objects good candidates for period bouncers. We also studied unusual behavior in some objects, including BK Lyn which made a transition from an ER UMa-type state to the novalike (standstill) state in 2013 and unusually frequent occurrence of superoutbursts in NY Ser and CR Boo. We applied least absolute shrinkage and selection operator (Lasso) power spectral analysis, which has been proven to be very effective in analyzing the Kepler data, to ground-based photometry of BK Lyn and detected the dramatic disappearance of the signal of negative superhumps in 2013. We suggested that the mass-transfer rates did not vary strongly between the ER UMa-type state and novalike state in BK Lyn, and this transition was less likely caused by a systematic variation of the mass-transfer rate.

We report on photometric observations of WZ Sge-type dwarf novae, MASTER OT J211258.65+242145.4 and MASTER OT J203749.39+552210.3 which underwent outbursts in 2012. Early superhumps were recorded in both systems. During superoutburst plateau, ordinary superhumps with a period of 0.060291(4) d (MASTER J211258) and of 0.061307(9) d (MASTER J203749) in average were observed. MASTER J211258 and MASTER J203749 exhibited eight and more than four post-superoutburst rebrightenings, respectively. In the final part of the superoutburst, an increase in the superhump periods was seen in both systems. We have made a survey of WZ Sge-type dwarf novae with multiple rebrightenings, and confirmed that the superhump periods of WZ Sge-type dwarf novae with multiple rebrightenings were longer than those of WZ Sge-type dwarf novae without a rebrightening. Although WZ Sge-type dwarf novae with multiple rebrightenings have been thought to be the good candidates for period bouncers based on their low mass ratio (q) from inferred from the period of fully grown (stage B) superhumps, our new method using the period of growing superhumps (stage A superhumps), however, implies higher q than those expected from stage B superhumps. These q values appear to be consistent with the duration of the stage A superoutbursts, which likely reflects the growth time of the 3:1 resonance. We present a working hypothesis that the small fractional superhump excesses for stage B superhumps in these systems may be explained as a result that a higher gas pressure effect works in these systems than in ordinary SU UMa-type dwarf novae. This result leads to a new picture that WZ Sge-type dwarf novae with multiple rebrightenings and SU UMa-type dwarf novae without a rebrightening (they are not period bouncers) are located in the same place on the evolutionary track.