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 high-resolution ($R \simeq 20,000$) near-infrared (9,100-13,500 $\AA$) long-slit spectrum of P Cygni obtained with the newly commissioned WINERED spectrograph in Japan. In the obtained spectrum, we have found that the velocity profiles of the [Fe II] emission lines are resolved into two peaks at a velocity of $\simeq220$ km s$^{-1}$ with a moderate dip in between and with additional sub-peaks at $\simeq\pm100$~km~s$^{-1}$. The sub-peak component is confirmed with the long-slit echellogram to originate in the known shell with a radius of $\simeq10^{\prime\prime}$, which was originally created by the outburst in 1600 AD. On the other hand, the $\simeq220$ km s$^{-1}$ component, which dominates the [Fe II] flux from P Cygni, is found to be concentrated closer to the central star with an apparent spatial extent of $\simeq3^{\prime\prime}$. The extent is much larger than the compact ($<0^{\prime\prime}.1$) regions traced with hydrogen, helium, and metal permitted lines. The velocity, estimated mass, and dynamical time of the extended emission-line region suggest that the region is an outer part of the stellar wind region. We suggest that the newly-identified emission-line region may trace a reverse shock due to the stellar wind overtaking the outburst shell.

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 an optical and near-infrared (hereafter NIR) polarimetric study of a comet C/2013 US10 (Catalina) observed on UT 2015 December 17-18 at phase angles of $\alpha$=52.1 deg - 53.1 deg. Additionally, we obtained an optical spectrum and multi-band images to examine the influence of gas emission. We find that the observed optical signals are significantly influenced by gas emission, that is, the gas-to-total intensity ratio varies from 5 to 30 % in the $R_{\rm C}$ and 3 to 18 % in the $I_{\rm C}$ bands, depending on the position in the coma. We derive the `gas-free dust polarization degrees' of 13.8$\pm$1.0 % in the $R_{\rm C}$ and 12.5$\pm$1.1 % in the $I_{\rm C}$ bands and a gray polarimetric color, i.e., -8.7$\pm$9.9 % $\mu \mathrm{m}$$^{-1}$ in optical and 1.6$\pm$0.9 % $\mu \mathrm{m}$$^{-1}$ in NIR. The increments of polarization obtained from the gas correction show that the polarimetric properties of the dust in this low-polarization comet are not different from those in high-polarization comets. In this process, the cometocentric distance dependence of polarization has disappeared. We also find that the $R_{\rm C}$-band polarization degree of the southeast dust tail, which consists of large dust particles (100 $\mu \mathrm{m}$ - 1 mm), is similar to that in the outer coma where small and large ones are mixed. Our study confirms that the dichotomy of cometary polarization does not result from the difference of dust properties, but from depolarizing gas contamination. This conclusion can provide a strong support for similarity in origin of comets.

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

Multiple outbursts of a Jupiter-family comet, 15P/Finlay, occurred from late 2014 to early 2015. We conducted an observation of the comet after the first outburst and subsequently witnessed another outburst on 2015 January 15.6-15.7. The gas, consisting mostly of C2 and CN, and dust particles expanded at speeds of 1,110 +/- 180 m/s and 570 +/- 40 m/s at a heliocentric distance of 1.0 AU. We estimated the maximum ratio of solar radiation pressure with respect to the solar gravity beta_max = 1.6 +/- 0.2, which is consistent with porous dust particles composed of silicates and organics. We found that 10^8-10^9 kg of dust particles (assumed to be 0.3 micron - 1 mm) were ejected through each outburst. Although the total mass is three orders of magnitude smaller than that of the 17P/Holmes event observed in 2007, the kinetic energy per unit mass (104 J/kg) is equivalent to the estimated values of 17P/Holmes and 332P/2010 V1 (Ikeya-Murakami), suggesting that the outbursts were caused by a similar physical mechanism. From a survey of cometary outbursts on the basis of voluntary reports, we conjecture that 15P/Finlay-class outbursts occur >1.5 times annually and inject dust particles from Jupiter-family comets and Encke-type comets into interplanetary space at a rate of ~10 kg/s or more.

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 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 obtained the near-infrared (NIR) high-resolution ($R\equiv\lambda/\Delta\lambda\sim20,000$) spectra of the seven brightest early-type stars in the Cygnus OB2 association for investigating the environmental dependence of diffuse interstellar bands (DIBs). The WINERED spectrograph mounted on the Araki 1.3m telescope in Japan was used to collect data. All 20 of the known DIBs within the wavelength coverage of WINERED ($0.91<\lambda<1.36\mu$m) were clearly detected along all lines of sight because of their high flux density in the NIR wavelength range and the large extinction. The equivalent widths (EWs) of DIBs were not correlated with the column densities of C$_2$ molecules, which trace the patchy dense component, suggesting that the NIR DIB carriers are distributed mainly in the diffuse component. On the basis of the correlations among the NIR DIBs both for stars in Cyg OB2 and stars observed previously, $\lambda\lambda$10780, 10792, 11797, 12623, and 13175 are found to constitute a "family", in which the DIBs are correlated well over the wide EW range. In contrast, the EW of $\lambda$10504 is found to remain almost constant over the stars in Cyg OB2. The extinction estimated from the average EW of $\lambda$10504 ($A_V\sim3.6$mag) roughly corresponds to the lower limit of the extinction distribution of OB stars in Cyg OB2. This suggests that $\lambda$10504 is absorbed only by the foreground clouds, implying that the carrier of $\lambda$10504 is completely destroyed in Cyg OB2, probably by the strong UV radiation field. The different behaviors of the DIBs may be caused by different properties of the DIB carriers.

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 present a comprehensive survey of diffuse interstellar bands (DIBs) in $0.91-1.32\mu$m with a newly-developed near-infrared (NIR) spectrograph, WINERED, mounted on the Araki 1.3 m Telescope in Japan. We obtained high-resolution ($R=28,300$) spectra of 25 early-type stars with color excesses of $0.07<E(B-V)<3.4$. In addition to the five DIBs previously detected in this wavelength range, we identified 15 new DIBs, 7 of which were reported as DIB "candidates" by Cox. We analyze the correlations among NIR DIBs, strong optical DIBs, and the reddening of the stars. Consequently, we found that all NIR DIBs show weaker correlations with the reddening rather than the optical DIBs, suggesting that the equivalent widths of NIR DIBs depend on some physical conditions of the interstellar clouds, such as UV flux. Three NIR DIBs, $\lambda\lambda$10780, 10792, and 11797, are found to be classifiable as a "family," in which the DIBs are well correlated with each other, suggesting that the carriers of these DIBs are connected with some chemical reactions and/or have similar physical properties such as ionization potential. We also found that three strongest NIR DIBs $\lambda\lambda10780, 11797, 13175$ are well correlated with the optical DIB $\lambda 5780.5$, whose carrier is proposed to be a cation molecule with high ionization potential, indicating that the carriers of the NIR DIBs could be cation molecules.

WINERED is a newly built high-efficiency (throughput$ > 25-30\%$) and high-resolution spectrograph customized for short NIR bands at 0.9-1.35 ${\rm \mu}$m. WINERED is equipped with ambient temperature optics and a cryogenic camera using a 1.7 ${\rm \mu}$m cut-off HgCdTe HAWAII-2RG array detector. WINERED has two grating modes: one with a conventional reflective echelle grating (R$\sim$28,300), which covers 0.9-1.35 $\mu$m simultaneously, the other with ZnSe or ZnS immersion grating (R$\sim$100,000). We have completed the development of WINERED except for the immersion grating, and started engineering and science observations at the Nasmyth platform of the 1.3 m Araki Telescope at Koyama Astronomical Observatory of Kyoto-Sangyo University in Japan. We confirmed that the spectral resolution ($R\sim$ 28,300) and the throughput ($>$ 40\% w/o telescope/atmosphere/array QE) meet our specifications. We measured ambient thermal backgrounds (e.g., 0.06 ${\rm [e^{-}/sec/pixel]}$ at 287 K), which are roughly consistent with that we expected. WINERED is a portable instrument that can be installed at any telescope with Nasmyth focus as a PI-type instrument. If WINERED is installed on a 10 meter telescope, the limiting magnitude is expected to be J=18-19, which can provide high-resolution spectra with high quality even for faint distant objects.

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.