We present visible and near-infrared photometric and spectroscopic observations of interstellar object 2I/Borisov taken from 2019 September 10 to 2019 November 29 using the GROWTH, the APO ARC 3.5 m and the NASA/IRTF 3.0 m combined with post and pre-discovery observations of 2I obtained by ZTF from 2019 March 17 to 2019 May 5. Comparison with imaging of distant Solar System comets shows an object very similar to mildly active Solar System comets with an out-gassing rate of $\sim$10$^{27}$ mol/sec. The photometry, taken in filters spanning the visible and NIR range shows a gradual brightening trend of $\sim0.03$ mags/day since 2019 September 10 UTC for a reddish object becoming neutral in the NIR. The lightcurve from recent and pre-discovery data reveals a brightness trend suggesting the recent onset of significant H$_2$O sublimation with the comet being active with super volatiles such as CO at heliocentric distances $>$6 au consistent with its extended morphology. Using the advanced capability to significantly reduce the scattered light from the coma enabled by high-resolution NIR images from Keck adaptive optics taken on 2019 October 04, we estimate a diameter of 2I's nucleus of $\lesssim$1.4 km. We use the size estimates of 1I/'Oumuamua and 2I/Borisov to roughly estimate the slope of the ISO size-distribution resulting in a slope of $\sim$3.4$\pm$1.2, similar to Solar System comets and bodies produced from collisional equilibrium.

We report on the variations of the physical parameters of the jet observed in the blazar Mrk 421, and discuss the origin of X-ray flares in the jet, based on the analysis of the several spectral energy distributions (SEDs). The SEDs are modeled using the one-zone synchrotron self-Compton (SSC) model and its parameters determined using a Markov chain Monte Carlo method. The lack of data at TeV energies means many of the parameters cannot be uniquely determined and are correlated. These are studied in detail. We found that the optimal solution can be uniquely determined only when we apply a constraint to one of four parameters: the magnetic field (B), Doppler factor, size of the emitting region, and normalization factor of the electron energy distribution. We used 31 sets of SED from 2009 to 2014 with optical-UV data observed with UVOT/Swift and the Kanata telescope, X-ray data with XRT/Swift, and gamma-ray data with the Fermi Large Area Telescope (LAT). The result of our SED analysis suggests that, in the X-ray faint state, the emission occurs in a relatively small area (~ 10^16 cm) with relatively strong magnetic field (B~10^-1 G). The X-ray bright state shows a tendency opposite to that of the faint state, that is, a large emitting area (~10^18 cm), probably in the downstream of the jet and weak magnetic field (B~10^-3 G). The high X-ray flux was due to an increase in the maximum energy of electrons. On the other hand, the presence of two kinds of emitting areas implies that the one-zone model is unsuitable to reproduce, at least a part of the observed SEDs.

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.

Gaia16aye was a binary microlensing event discovered in the direction towards the northern Galactic disc and was one of the first microlensing events detected and alerted to by the Gaia space mission. Its light curve exhibited five distinct brightening episodes, reaching up to I=12 mag, and it was covered in great detail with almost 25,000 data points gathered by a network of telescopes. We present the photometric and spectroscopic follow-up covering 500 days of the event evolution. We employed a full Keplerian binary orbit microlensing model combined with the motion of Earth and Gaia around the Sun to reproduce the complex light curve. The photometric data allowed us to solve the microlensing event entirely and to derive the complete and unique set of orbital parameters of the binary lensing system. We also report on the detection of the first-ever microlensing space-parallax between the Earth and Gaia located at L2. The properties of the binary system were derived from microlensing parameters, and we found that the system is composed of two main-sequence stars with masses 0.57$\pm$0.05 $M_\odot$ and 0.36$\pm$0.03 $M_\odot$ at 780 pc, with an orbital period of 2.88 years and an eccentricity of 0.30. We also predict the astrometric microlensing signal for this binary lens as it will be seen by Gaia as well as the radial velocity curve for the binary system. Events such as Gaia16aye indicate the potential for the microlensing method of probing the mass function of dark objects, including black holes, in directions other than that of the Galactic bulge. This case also emphasises the importance of long-term time-domain coordinated observations that can be made with a network of heterogeneous telescopes.

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.

We report results from X-ray and optical observations of the Galactic black hole candidate MAXI J1828-249, performed with Suzaku and the Kanata telescope around the X-ray flux peak in the 2013 outburst. The time-averaged X-ray spectrum covering 0.6--168 keV was approximately characterized by a strong multi-color disk blackbody component with an inner disk temperature of ~0.6 keV, and a power-law tail with a photon index of ~2.0. We detected an additional structure at 5-10 keV, which can be modelled neither with X-ray reflection on the disk, nor relativistic broadening of the disk emission. Instead, it was successfully reproduced with a Comptonization of disk photons by thermal electrons with a relatively low temperature (<~10 keV). We infer that the source was in the intermediate state, considering its long-term trend in the hardness intensity diagram, the strength of the spectral power-law tail, and its variability properties. The low-temperature Comptonization component could be produced in a boundary region between the truncated standard disk and the hot inner flow, or a Comptonizing region that somehow developed above the disk surface. The multi-wavelength spectral energy distribution suggests that the optical and UV fluxes were dominated by irradiated outer disk emission.

We performed deep spectropolarimetric observations of a prototypical starburst galaxy M82 with the Subaru Telescope in order to study the kinematics of the dust outflow. We obtained optical polarized emission-line spectra up to $\sim$4~kpc away from the nucleus of the galaxy along three position angles within the conical outflowing wind (superwind). The H$\alpha$ emission line in the superwind is strongly polarized and the polarization pattern shows dust scattering of central light sources, being consistent with the previous works. The intensity weighted polarization degree of the H$\alpha$ line reaches $\sim$30\%\ at maximum. There are at least two light sources at the central region of the galaxy; one of which is located at the near-infrared nucleus and the other resides at one of the peaks of the 3 mm radio and molecular gas emission. The outer (> 1 kpc) dust is illuminated by the former, whereas the inner dust is scattering the light from the latter. We investigated as well the dust motion from the velocity field of the polarized H$\alpha$ emission line. The dust is accelerated outward on the northwest side of the nucleus. A simple bi-conical dust outflow model shows that the outflow velocity of the dust reaches $\gtrsim 300-450$ km/s at $\sim 4$ kpc from the nucleus, suggesting that some portion of the dust escapes from the gravitational potential of M82 into the intergalactic space. At some regions on the southeast side, in particular along the position angle of 138$^\circ$, the dust has radial velocity slower than the systemic velocity of the galaxy, apparently suggesting inflowing motion toward the nucleus. These components are spatially consistent with a part of the molecular gas stream, which is kinematically independent of the outflow gas, thus the apparent inflow motion of the dust reflects the streaming motion associated with the molecular gas stream.

We report the Subaru Hyper Suprime-Cam (HSC) discovery of two Ly$\alpha$ blobs (LABs), dubbed z70-1 and z49-1 at $z=6.965$ and $z=4.888$ respectively, that are Ly$\alpha$ emitters with a bright ($\log L_{\rm Ly\alpha}/{\rm [erg\ s^{-1}]}>43.4$) and spatially-extended Ly$\alpha$ emission, and present the photometric and spectroscopic properties of a total of seven LABs; the two new LABs and five previously-known LABs at $z=5.7-6.6$. The z70-1 LAB shows the extended Ly$\alpha$ emission with a scale length of $1.4\pm 0.2$ kpc, about three times larger than the UV continuum emission, making z70-1 the most distant LAB identified to date. All of the 7 LABs, except z49-1, exhibit no AGN signatures such as X-ray emission, {\sc Nv}$\lambda$1240 emission, or Ly$\alpha$ line broadening, while z49-1 has a strong {\sc Civ}$\lambda$1548 emission line indicating an AGN on the basis of the UV-line ratio diagnostics. We carefully model the point-spread functions of the HSC images, and conduct two-component exponential profile fitting to the extended Ly$\alpha$ emission of the LABs. The Ly$\alpha$ scale lengths of the core (star-forming region) and the halo components are $r_{\rm c}=0.6-1.2$ kpc and $r_{\rm h}=2.0-13.8$ kpc, respectively. The average $r_{\rm h}$ of the LABs falls on the extrapolation of the $r_{\rm h}$-Ly$\alpha$ luminosity relation of the Ly$\alpha$ halos around VLT/MUSE star-forming galaxies at the similar redshifts, suggesting that typical LABs at $z\gtrsim5$ are not special objects, but star-forming galaxies at the bright end.

We describe a dynamic science portal called the GROWTH Marshal that allows time-domain astronomers to define science programs, program filters to save sources from different discovery streams, co-ordinate follow-up with various robotic or classical telescopes, analyze the panchromatic follow-up data and generate summary tables for publication. The GROWTH marshal currently serves 137 scientists, 38 science programs and 67 telescopes. Every night, in real-time, several science programs apply various customized filters to the 10^5 nightly alerts from the Zwicky Transient Facility. Here, we describe the schematic and explain the functionality of the various components of this international collaborative platform.

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.

Wide-field optical surveys have begun to uncover large samples of fast (t_rise < 5d), luminous (M_peak < -18), blue transients. While commonly attributed to the breakout of a supernova shock into a dense wind, the great distances to the transients of this class found so far have hampered detailed investigation of their properties. We present photometry and spectroscopy from a comprehensive worldwide campaign to observe AT2018cow (ATLAS18qqn), the first fast-luminous optical transient to be found in real time at low redshift. Our first spectra (<2 days after discovery) are entirely featureless. A very broad absorption feature suggestive of near-relativistic velocities develops between 3-8 days, then disappears. Broad emission features of H and He develop after >10 days. The spectrum remains extremely hot throughout its evolution, and the photospheric radius contracts with time (receding below R<10^14 cm after 1 month). This behaviour does not match that of any known supernova, although a relativistic jet within a fallback supernova could explain some of the observed features. Alternatively, the transient could originate from the disruption of a star by an intermediate-mass black hole, although this would require long-lasting emission of highly super-Eddington thermal radiation. In either case, AT2018cow suggests that the population of fast luminous transients represents a new class of astrophysical event. Intensive follow-up of this event in its late phases, and of any future events found at comparable distance, will be essential to better constrain their origins.

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.

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 the Ly$\alpha$ luminosity function (LF) derived from 34 Ly$\alpha$ emitters (LAEs) at $z=7.0$ on the sky of $3.1$ deg$^2$, the largest sample compared to those in the literature obtained at a redshift $z\gtrsim7$. The LAE sample is made by deep large-area Subaru narrowband observations conducted by the Cosmic HydrOgen Reionization Unveiled with Subaru (CHORUS) project. The $z=7.0$ Ly$\alpha$ LF of our project is consistent with those of the previous DECam and Subaru studies at the bright and faint ends, respectively, while our $z=7.0$ Ly$\alpha$ LF has uncertainties significantly smaller than those of the previous study results. Exploiting the small errors of our measurements, we investigate the shape of the faint to bright-end Ly$\alpha$ LF. We find that the $z=7.0$ Ly$\alpha$ LF shape can be explained by the steep slope of $\alpha \simeq -2.5$ suggested at $z=6.6$, and that there is no clear signature of a bright-end excess at $z\simeq 7$ claimed by the previous work, which was thought to be made by the ionized bubbles around bright LAEs whose Ly$\alpha$ photons could easily escape from the partly neutral IGM at $z \simeq 7$. We estimate the Ly$\alpha$ luminosity densities (LDs) with Ly$\alpha$ LFs at $z\simeq 6-8$ given by our and the previous studies, and compare the evolution of the UV-continuum LD estimated with dropouts. The Ly$\alpha$ LD monotonically decreases from $z\sim 6$ to $8$, and evolves stronger than the UV-continuum LD, indicative of the Ly$\alpha$ damping wing absorption of the IGM towards the heart of the reionization epoch.

We analyze the linear polarization of the relativistic jet in BL Lacertae object OJ~287 as revealed by multi-epoch Very Long Baseline Array (VLBA) images at 43 GHz and monitoring observations at optical bands. The electric-vector position angle (EVPA) of the optical polarization matches that at 43 GHz at locations that are often in the compact millimeter-wave "core" or, at other epochs, coincident with a bright, quasi-stationary emission feature $\sim0.2$~milliarcsec ($\sim$0.9~pc projected on the sky) downstream from the core. This implies that electrons with high enough energies to emit optical synchrotron and $\gamma$-ray inverse Compton radiation are accelerated both in the core and at the downstream feature, the latter of which lies $\geq10$~pc from the central engine. The polarization vector in the stationary feature is nearly parallel to the jet axis, as expected for a conical standing shock capable of accelerating electrons to GeV energies.

We have monitored the R-band optical linear polarisation of ten jetted NLSy1 galaxies with the aim to quantify their variability and search for candidate long rotation of the polarisation plane. In all cases for which adequate datasets are available we observe significant variability of both the polarisation fraction and angle. In the best-sampled cases we identify candidate long rotations of the polarisation plane. We present an approach that assesses the probability that the observed phenomenology is the result of pure noise. We conclude that although this possibility cannot be excluded it is much more likely that the EVPA undergoes an intrinsic evolution. We compute the most probable parameters of the intrinsic event which forecasts events consistent with the observations. In one case we find that the EVPA shows a preferred direction which, however, does not imply any dominance of a toroidal or poloidal component of the magnetic field at those scales.

The long gamma-ray burst GRB 161017A was detected by Fermi and Swift, and its afterglow was observed by the MITSuME 50-cm optical telescope promptly about 50 s after the burst. Early optical observations revealed that the optical lightcurve exhibits a plateau and rebrightening in the early afterglow phase about 500 and 5000 s after the trigger, respectively. By investigating the behavior of the spectral and temporal flux variation, it was found that the plateau and re-brightening cannot be explained in the context of the simple standard afterglow model. These observational features can be explained with two independent refreshed shocks, which indicate the long-acting central engine. We evaluated the physical parameters of the subsequent shells, and we then determined the kinetic energy ratio of the two colliding shells to the leading shell to be roughly 1 and 8, respectively. In addition, two prominent X-ray flares about 200 s after the trigger may be signatures of delayed ejections of the energetic jets responsible for the refreshed shocks. Such late activity of the central engine and X-ray flares play a crucial role in understanding the mechanisms for jet formation and photon emission.

Narrow line Seyfert 1 galaxies (NLSy1s) constitute the AGN subclass associated with systematically smaller black hole masses. A few radio loud ones have been detected in MeV -- GeV energy bands by Fermi and evidence for the presence of blazar-like jets has been accumulated. In this study we wish to quantify the temporal behaviour of the optical polarisation, fraction and angle, for a selected sample of radio loud NLSy1s. We also search for rotations of the polarisation plane similar to those commonly observed in blazars. We have conducted R-band optical polarisation monitoring of a sample of 10 RL NLSy1s 5 of which have been previously detected by Fermi. The dataset includes observations with the RoboPol, KANATA, Perkins and Steward polarimeters. In the cases where evidences for long rotations of the polarisation plane are found, we carry out numerical simulations to assess the probability that they are caused by intrinsically evolving EVPAs instead of observational noise. Even our moderately sampled sources show indications of variability, both in polarisation fraction and angle. For the four best sampled objects in our sample we find multiple periods of significant polarisation angle variability. In the two best sampled cases, namely J1505+0326 and J0324+3410, we find indications for three long rotations. We show that although noise can induce the observed behaviour, it is much more likely that the apparent rotation is caused by intrinsic evolution of the EVPA. To our knowledge this is the very first detection of such events in this class of sources. In the case of the largest dataset (J0324+3410) we find that the EVPA concentrates around a direction which is at 49.3\degr to the 15-GHz radio jet implying a projected magnetic field at an angle of 40.7\degr to that axis.

The BL Lac object S5~0716+714, a highly variable blazar, underwent an impressive outburst in January 2015 (Phase A), followed by minor activity in February (Phase B). The MAGIC observations were triggered by the optical flux observed in Phase A, corresponding to the brightest ever reported state of the source in the R-band. The comprehensive dataset collected is investigated in order to shed light on the mechanism of the broadband emission. Multi-wavelength light curves have been studied together with the broadband Spectral Energy Distributions (SEDs). The data set collected spans from radio, optical photometry and polarimetry, X-ray, high-energy (HE, 0.1 GeV < E < 100 GeV) with \textit{Fermi}-LAT to the very-high-energy (VHE, E>100 GeV) with MAGIC. The flaring state of Phase A was detected in all the energy bands, providing for the first time a multi-wavelength sample of simultaneous data from the radio band to the VHE. In the constructed SED the \textit{Swift}-XRT+\textit{NuSTAR} data constrain the transition between the synchrotron and inverse Compton components very accurately, while the second peak is constrained from 0.1~GeV to 600~GeV by \textit{Fermi}+MAGIC data. The broadband SED cannot be described with a one-zone synchrotron self-Compton model as it severely underestimates the optical flux in order to reproduce the X-ray to $\gamma$-ray data. Instead we use a two-zone model. The EVPA shows an unprecedented fast rotation. An estimation of the redshift of the source by combined HE and VHE data provides a value of $z = 0.31 \pm 0.02_{stats} \pm 0.05_{sys}$, confirming the literature value. The data show the VHE emission originating in the entrance and exit of a superluminal knot in and out a recollimation shock in the inner jet. A shock-shock interaction in the jet seems responsible for the observed flares and EVPA swing. This scenario is also consistent with the SED modelling.

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 investigate Lya, [OIII]5007, Ha, and [CII]158um emission from 1124 galaxies at z=4.9-7.0. Our sample is composed of 1092 Lya emitters (LAEs) at z=4.9, 5.7, 6.6, and 7.0 identified by Subaru/Hyper Suprime-Cam (HSC) narrowband surveys covered by Spitzer large area survey with Subaru/HSC (SPLASH) and 34 galaxies at z=5.148-7.508 with deep ALMA [CII]158um data in the literature. Fluxes of strong rest-frame optical lines of [OIII] and Ha (Hb) are constrained by significant excesses found in the SPLASH 3.6 and 4.5um photometry. At z=4.9, we find that the rest-frame Ha equivalent width and the Lya escape fraction f_Lya positively correlate with the rest-frame Lya equivalent width EW^0_Lya. The f_Lya-EW^0_Lya correlation is similarly found at z~0-2, suggesting no evolution of the correlation over z~0-5. The typical ionizing photon production efficiency of LAEs is logxi_ion/[Hz erg^-1]~25.5 significantly (60-100%) higher than those of LBGs at a given UV magnitude. At z=5.7-7.0, there exists an interesting turn-over trend that the [OIII]/Ha flux ratio increases in EW^0_Lya~0-30 A, and then decreases out to EW^0_Lya~130 A. We also identify an anti-correlation between a [CII] luminosity to star-formation rate ratio (L_[CII]/SFR) and EW^0_Lya at the >99% confidence level. We carefully investigate physical origins of the correlations with stellar-synthesis and photoionization models, and find that a simple anti-correlation between EW_Lya^0 and metallicity explains self-consistently all of the correlations of Lya, Ha, [OIII]/Ha, and [CII] identified in our study, indicating detections of metal-poor (~0.03 Zo) galaxies with EW^0_Lya~200 A.

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.

Merging neutron stars offer an exquisite laboratory for simultaneously studying strong-field gravity and matter in extreme environments. We establish the physical association of an electromagnetic counterpart EM170817 to gravitational waves (GW170817) detected from merging neutron stars. By synthesizing a panchromatic dataset, we demonstrate that merging neutron stars are a long-sought production site forging heavy elements by r-process nucleosynthesis. The weak gamma-rays seen in EM170817 are dissimilar to classical short gamma-ray bursts with ultra-relativistic jets. Instead, we suggest that breakout of a wide-angle, mildly-relativistic cocoon engulfing the jet elegantly explains the low-luminosity gamma-rays, the high-luminosity ultraviolet-optical-infrared and the delayed radio/X-ray emission. We posit that all merging neutron stars may lead to a wide-angle cocoon breakout; sometimes accompanied by a successful jet and sometimes a choked jet.

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 report on the variation in the optical polarization of the blazar PKS 1749+096 observed in 2008--2015. The degree of polarization (PD) tends to increase in short flares having a time-scale of a few days. The object favors a polarization angle (PA) of $40^\circ$--$50^\circ$ at the flare maxima, which is close to the position angle of the jet ($20^\circ$--$40^\circ$). Three clear polarization rotations were detected in the negative PA direction associated with flares. In addition, a rapid and large decrease in the PA was observed in the other two flares, while another two flares showed no large PA variation. The light curve maxima of the flares possibly tend to lag behind the PD maxima and color-index minima. The PA became $-50^\circ$ to $-20^\circ$ in the decay phase of active states, which is almost perpendicular to the jet position angle. We propose a scenario to explain these observational features, where transverse shocks propagate along curved trajectories. The favored PA at the flare maxima suggests that the observed variations were governed by the variations in the Doppler factor, $\delta$. Based on this scenario, the minimum viewing angle of the source, $\theta_\mathrm{min}=4.8^\circ$--$6.6^\circ$, and the location of the source, $\Delta r\gtrsim 0.1$pc, from the central black hole were estimated. In addition, the acceleration of electrons by the shock and synchrotron cooling would have a time-scale similar to that of the change in $\delta$. The combined effect of the variation in $\delta$ and acceleration/cooling of electrons is probably responsible for the observed diversity of the polarization variations in the flares.

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 multi-wavelength follow-up campaigns by the AstroSat-CZTI and GROWTH collaborations to search for an electromagnetic counterpart to the gravitational wave event GW170104. At the time of the GW170104 trigger, the AstroSat CZTI field-of-view covered 50.3\% of the sky localization. We do not detect any hard X-ray (>100 keV) signal at this time, and place an upper limit of $\approx 4.5 \times 10^{-7}~{\rm erg~cm}^{-2}{\rm~s}^{-1}$ for a 1\,s timescale. Separately, the ATLAS survey reported a rapidly fading optical source dubbed ATLAS17aeu in the error circle of GW170104. Our panchromatic investigation of ATLAS17aeu shows that it is the afterglow of an unrelated long, soft GRB~170105A, with only a fortuitous spatial coincidence with GW170104. We then discuss the properties of this transient in the context of standard long GRB afterglow models.

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 explore physical properties of the shocked external medium (i.e., a shell) in 3C 84 associated with the recurrent radio lobe born around 1960. In the previous work of Ito et al., we investigated a dynamical and radiative evolution of such a shell after the central engine stops the jet launching and we found that a fossil shell emission overwhelms that of the rapidly fading radio lobe. We apply this model to 3C 84 and find the followings: (i) The fossil shell made of shocked diffuse ambient matter with the number density of 0.3 cm$^{-3}$ radiates bright Inverse-Compton (IC) emission with the seed photons of the radio emission from the central compact region and the IC emission is above the sensitivity threshold of the Cherenkov Telescope Array (CTA). (ii) When the fossil shell is produced in a geometrically thick ionized plasma with the number density of $10^{3}$ cm$^{-3}$ and the field strength in the shell may reach about 17 mG in the presence of magnetic fields amplification and the radio emission becomes comparable to the sensitivity of deep imaging VLBI observations. A possible production of ultra high energy cosmic-rays (UHECRs) in the dense shocked plasma is also argued.

S5 0716+714 is a well known BL Lac object, one of the brightest and most active blazars. The discovery in the Very High Energy band (VHE, E > 100 GeV) by MAGIC happened in 2008. In January 2015 the source went through the brightest optical state ever observed, triggering MAGIC follow-up and a VHE detection with 13{\sigma} significance (ATel 6999). Rich multiwavelength coverage of the flare allowed us to construct the broad-band spectral energy distribution of S5 0716+714 during its brightest outburst. In this work we will present the preliminary analysis of MAGIC and Fermi-LAT data of the flaring activity in January and February 2015 for the HE (0.1 < HE < 300 GeV) and VHE band, together with radio (Mets\"ahovi, OVRO, VLBA, Effelsberg), sub-millimeter (SMA), optical (Tuorla, Perkins, Steward, AZT-8+ST7, LX-200, Kanata), X-ray and UV (Swift-XRT and UVOT), in the same time-window and discuss the time variability of the multiwavelength light curves during this impressive outburst.

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.

Here we report on the results of the WEBT photo-polarimetric campaign targeting the blazar S5~0716+71, organized in March 2014 to monitor the source simultaneously in BVRI and near IR filters. The campaign resulted in an unprecedented dataset spanning $\sim 110$\,h of nearly continuous, multi-band observations, including two sets of densely sampled polarimetric data mainly in R filter. During the campaign, the source displayed pronounced variability with peak-to-peak variations of about $30\%$ and "bluer-when-brighter" spectral evolution, consisting of a day-timescale modulation with superimposed hourlong microflares characterized by $\sim 0.1$\,mag flux changes. We performed an in-depth search for quasi-periodicities in the source light curve; hints for the presence of oscillations on timescales of $\sim 3$\,h and $\sim 5$\,h do not represent highly significant departures from a pure red-noise power spectrum. We observed that, at a certain configuration of the optical polarization angle relative to the positional angle of the innermost radio jet in the source, changes in the polarization degree led the total flux variability by about 2\,h; meanwhile, when the relative configuration of the polarization and jet angles altered, no such lag could be noted. The microflaring events, when analyzed as separate pulse emission components, were found to be characterized by a very high polarization degree ($> 30\%$) and polarization angles which differed substantially from the polarization angle of the underlying background component, or from the radio jet positional angle. We discuss the results in the general context of blazar emission and energy dissipation models.

We analyzed Suzaku/XIS data of 2006--2015 observations of a gamma-ray emitting radio galaxy NGC 1275, and brightening of the nucleus in the X-ray band was found in 2013--2015, correlating with GeV Gamma-ray brightening. This is the first evidence of variability with correlation between GeV gamma-ray and X-ray for NGC 1275. We also analyzed Swift/XRT data of NGC 1275, and found that X-ray was flaring by a factor of $\sim$5 in several days in 2006, 2010, and 2013. The X-ray spectrum during the flare was featureless and somewhat steeper with a photon index of $\sim$2 against $\sim$1.7 in the normal state, indicating that a synchrotron component became brighter. A large Xray to GeV gamma-ray flux ratio in the flare could be explained by the shock-in-jet scenario. On the other hand, a long-term gradual brightening of radio, X-ray, and GeV gamma-ray with a larger gamma-ray amplitude could be origin of other than internal shocks, and then we discuss some possibilities.

Over the past few years, several occasions of large, continuous rotations of the electric vector position angle (EVPA) of linearly polarized optical emission from blazars have been reported. These events are often coincident with high energy gamma-ray flares and they have attracted considerable attention, as they could allow one to probe the magnetic field structure in the gamma-ray emitting region of the jet. The flat-spectrum radio quasar 3C279 is one of the most prominent examples showing this behaviour. Our goal is to study the observed EVPA rotations and to distinguish between a stochastic and a deterministic origin of the polarization variability. We have combined multiple data sets of R-band photometry and optical polarimetry measurements of 3C279, yielding exceptionally well-sampled flux density and polarization curves that cover a period of 2008-2012. Several large EVPA rotations are identified in the data. We introduce a quantitative measure for the EVPA curve smoothness, which is then used to test a set of simple random walk polarization variability models against the data. 3C279 shows different polarization variation characteristics during an optical low-flux state and a flaring state. The polarization variation during the flaring state, especially the smooth approx. 360 degrees rotation of the EVPA in mid-2011, is not consistent with the tested stochastic processes. We conclude that during the two different optical flux states, two different processes govern the polarization variation, possibly a stochastic process during the low-brightness state and a deterministic process during the flaring activity.

After several years of quiescence, the blazar CTA 102 underwent an exceptional outburst in 2012 September-October. The flare was tracked from gamma-ray to near-infrared frequencies, including Fermi and Swift data as well as photometric and polarimetric data from several observatories. An intensive GASP-WEBT collaboration campaign in optical and NIR bands, with an addition of previously unpublished archival data and extension through fall 2015, allows comparison of this outburst with the previous activity period of this blazar in 2004-2005. We find remarkable similarity between the optical and gamma-ray behaviour of CTA 102 during the outburst, with a time lag between the two light curves of ~1 hour, indicative of co-spatiality of the optical and gamma-ray emission regions. The relation between the gamma-ray and optical fluxes is consistent with the SSC mechanism, with a quadratic dependence of the SSC gamma-ray flux on the synchrotron optical flux evident in the post-outburst stage. However, the gamma-ray/optical relationship is linear during the outburst; we attribute this to changes in the Doppler factor. A strong harder-when-brighter spectral dependence is seen both the in gamma-ray and optical non-thermal emission. This hardening can be explained by convexity of the UV-NIR spectrum that moves to higher frequencies owing to an increased Doppler shift as the viewing angle decreases during the outburst stage. The overall pattern of Stokes parameter variations agrees with a model of a radiating blob or shock wave that moves along a helical path down the jet.

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 our optical follow-up observations to search for an electromagnetic counterpart of the first gravitational wave source GW150914 in the framework of the Japanese collaboration for Gravitational wave ElectroMagnetic follow-up (J-GEM), which is an observing group utilizing optical and radio telescopes in Japan, as well as those in New Zealand, China, South Africa, Chile, and Hawaii. We carried out a wide-field imaging survey with Kiso Wide Field Camera (KWFC) on the 1.05-m Kiso Schmidt telescope in Japan and a galaxy-targeted survey with Tripole5 on the B&C 61-cm telescope in New Zealand. Approximately 24 deg2 regions in total were surveyed in i-band with KWFC and 18 nearby galaxies were observed with Tripole5 in g-, r-, and i-bands 4-12 days after the gravitational wave detection. Median 5-sigma depths are i~18.9 mag for the KWFC data and g~18.9 mag, r~18.7 mag, and i~18.3 mag for the Tripole5 data. Probability for a counterpart to be in the observed area is 1.2% in the initial skymap and 0.1% in the final skymap. We do not find any transient source associated to an external galaxy with spatial offset from its center, which is consistent with the local supernova rate. We summarize future prospects and ongoing efforts to pin down electromagnetic counterparts of binary black hole mergers as well as neutron star mergers.

We report on Fermi Large Area Telescope (LAT) and multi-wavelength results on the recently-discovered very-high-energy (VHE, $E>$ 100 GeV) blazar S4 0954+65 ($z=0.368$) during an exceptionally bright optical flare in 2015 February. During the time period (2015 February, 13/14, or MJD 57067) when the MAGIC telescope detected VHE $\gamma$-ray emission from the source, the Fermi-LAT data indicated a significant spectral hardening at GeV energies, with a power-law photon index of $1.8 \pm 0.1$---compared with the 3FGL value (averaged over four years of observation) of $2.34 \pm 0.04$. In contrast, Swift/XRT data showed a softening of the X-ray spectrum, with a photon index of $1.72 \pm 0.08$ (compared with $1.38 \pm 0.03$ averaged during the flare from MJD 57066 to 57077), possibly indicating a modest contribution of synchrotron photons by the highest-energy electrons superposed on the inverse Compton component. Fitting of the quasi-simultaneous ($<1$ day) broadband spectrum with a one-zone synchrotron plus inverse-Compton model revealed that GeV/TeV emission could be produced by inverse-Compton scattering of external photons from the dust torus. We emphasize that a flaring blazar showing high flux of $\gtrsim 1.0 \times 10^{-6}$ photons cm$^{-2}$ s$^{-1}$ ($E>$ 100 MeV) and a hard spectral index of $\Gamma_{\rm GeV} < 2.0$ detected by Fermi-LAT on daily time scales is a promising target for TeV follow-up by ground-based Cherenkov telescopes to discover high-redshift blazars, investigate their temporal variability and spectral features in the VHE band, and also constrain the intensity of the extragalactic background light.

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.

Detailed investigation of broadband flux variability in the blazar 3C 273 allows us to probe the location and size of emission regions and their physical conditions. We report the results on correlation studies of the flaring activity observed between 2008 and 2012. The observed broadband variations were investigated using the structure function and the discrete correlation function, and power spectral density analysis (PSD) methods. The PSD analysis showed that the optical/IR light curve slopes are consistent with the slope of white noise processes, while, the PSD slopes at radio, X-ray and gamma-ray energies are consistent with red-noise processes. The flux variations at gamma-ray and mm-radio bands are found to be significantly correlated. Using the estimated time lag of (110\pm27) days between gamma-ray and radio light curves, we constrained the location of the gamma-ray emission region at a de-projected distance of 1.2\pm0.9 pc from the jet apex. Flux variations at X-ray bands were found to have a significant correlation with variations at both radio and \gamma-rays energies. The correlation between X-rays and gamma-rays light curves suggests presence of two components responsible for the X-ray emission. A negative time lag of -(50\pm20) days, where the X-rays are leading the emission, suggests X-rays are emitted closer to the jet apex from a compact region at a distance of ~(0.5\pm0.4) pc from the jet apex. A positive time lag of (110\pm20) days suggests jet-base origin of the other X-ray component at ~(4--5)~pc from the jet apex. The flux variations at radio frequencies were found to be well correlated with each other such that the variations at higher frequencies are leading the lower frequencies, which could be expected in the standard shock-in-jet model.

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 investigate shell emission associated with dying radio loud AGNs. First, based on our recent work by Ito et al. (2015), we describe the dynamical and spectral evolutions of shells after stopping the jet energy injection. We find that the shell emission overwhelms that of the radio lobes soon after stopping the jet energy injection because fresh electrons are continuously supplied into the shell via the forward shock while the radio lobes rapidly fade out without jet energy injection. We find that such fossil shells can be a new class of target sources for SKA telescope. Next, we apply the model to the nearby radio source 3C84. Then, we find that the fossil shell emission in 3C84 is less luminous in radio band while it is bright in TeV gamma-ray band and it can be detectable by CTA.

The magnetar 4U~0142+61 has been well studied at optical and infrared wavelengths and is known to have a complicated broad-band spectrum over the wavelength range. Here we report the result from our linear imaging polarimetry of the magnetar at optical $I$-band. From the polarimetric observation carried out with the 8.2-m Subaru telescope, we determine the degree of linear polarization $P=1.0\pm$3.4\%, or $P\leq$5.6\% (90\% confidence level). Considering models suggested for optical emission from magnetars, we discuss the implications of our result. The upper limit measurement indicates that different from radio pulsars, magnetars probably would not have strongly polarized optical emission if the emission arises from their magnetosphere as suggested.

We conducted an optical and near-infrared polarimetric observation of the highly dormant Jupiter-Family Comet, 209P/LINEAR. Because of its low activity, we were able to determine the linear polarization degrees of the coma dust particles and nucleus independently, that is $P_n$=30.3$^{+1.3}_{-0.9}$% at $\alpha$=92.2$^\circ$ and $P_n$=31.0$^{+1.0}_{-0.7}$% at $\alpha$=99.5$^\circ$ for the nucleus, and $P_c$=28.8$^{+0.4}_{-0.4}$% at $\alpha$=92.2$^\circ$ and 29.6$^{+0.3}_{-0.3}$% at $\alpha$=99.5$^\circ$ for the coma. We detected no significant variation in $P$ at the phase angle coverage of 92.2$^\circ$-99.5$^\circ$, which may imply that the obtained polarization degrees are nearly at maximum in the phase-polarization curves. By fitting with an empirical function, we obtained the maximum values of linear polarization degrees $P_\mathrm{max}$=30.8% for the nucleus and $P_\mathrm{max}$=29.6% for the dust coma. The $P_\mathrm{max}$ of the dust coma is consistent with those of dust-rich comets. The low geometric albedo of $P_v$=0.05 was derived from the slope-albedo relationship and was associated with high $P_\mathrm{max}$. We examined $P_\mathrm{max}$-albedo relations between asteroids and 209P, and found that the so-called Umov law seems to be applicable on this cometary surface.

We perform a multi-wavelength polarimetric study of the quasar CTA 102 during an extraordinarily bright $\gamma$-ray outburst detected by the {\it Fermi} Large Area Telescope in September-October 2012 when the source reached a flux of F$_{>100~\mathrm{MeV}} =5.2\pm0.4\times10^{-6}$ photons cm$^{-2}$ s$^{-1}$. At the same time the source displayed an unprecedented optical and NIR outburst. We study the evolution of the parsec scale jet with ultra-high angular resolution through a sequence of 80 total and polarized intensity Very Long Baseline Array images at 43 GHz, covering the observing period from June 2007 to June 2014. We find that the $\gamma$-ray outburst is coincident with flares at all the other frequencies and is related to the passage of a new superluminal knot through the radio core. The powerful $\gamma$-ray emission is associated with a change in direction of the jet, which became oriented more closely to our line of sight ($\theta\sim$1.2$^{\circ}$) during the ejection of the knot and the $\gamma$-ray outburst. During the flare, the optical polarized emission displays intra-day variability and a clear clockwise rotation of EVPAs, which we associate with the path followed by the knot as it moves along helical magnetic field lines, although a random walk of the EVPA caused by a turbulent magnetic field cannot be ruled out. We locate the $\gamma$-ray outburst a short distance downstream of the radio core, parsecs from the black hole. This suggests that synchrotron self-Compton scattering of near-infrared to ultraviolet photons is the probable mechanism for the $\gamma$-ray production.

The occurrence of low-amplitude flux variations in blazars on hourly timescales, commonly known as microvariability, is still a widely debated subject in high-energy astrophysics. Several competing scenarios have been proposed to explain such occurrences, including various jet plasma instabilities leading to the formation of shocks, magnetic reconnection sites, and turbulence. In this letter we present the results of our detailed investigation of a prominent, five-hour-long optical microflare detected during recent WEBT campaign in 2014, March 2-6 targeting the blazar 0716+714. After separating the flaring component from the underlying base emission continuum of the blazar, we find that the microflare is highly polarized, with the polarization degree $\sim (40-60)\%$$\pm (2-10)\%$, and the electric vector position angle $\sim (10 - 20)$deg$\pm (1-8)$deg slightly misaligned with respect to the position angle of the radio jet. The microflare evolution in the $(Q,\,U)$ Stokes parameter space exhibits a looping behavior with a counter-clockwise rotation, meaning polarization degree decreasing with the flux (but higher in the flux decaying phase), and approximately stable polarization angle. The overall very high polarization degree of the flare, its symmetric flux rise and decay profiles, and also its structured evolution in the $Q-U$ plane, all imply that the observed flux variation corresponds to a single emission region characterized by a highly ordered magnetic field. As discussed in the paper, a small-scale but strong shock propagating within the outflow, and compressing a disordered magnetic field component, provides a natural, though not unique, interpretation of our findings.

We report observations of the Type Iax supernova (SN Iax) 2012Z at optical and near-infrared wavelengths from immediately after the explosion until $\sim$ $260$ days after the maximum luminosity using the Optical and Infrared Synergetic Telescopes for Education and Research (OISTER) Target-of-Opportunity (ToO) program and the Subaru telescope. We found that the near-infrared (NIR) light curve evolutions and color evolutions are similar to those of SNe Iax 2005hk and 2008ha. The NIR absolute magnitudes ($M_{J}\sim-18.1$ mag and $M_{H}\sim-18.3$ mag) and the rate of decline of the light curve ($\Delta$ $m_{15}$($B$)$=1.6 \pm 0.1$ mag) are very similar to those of SN 2005hk ($M_{J}\sim-17.7$ mag, $M_{H}\sim$$-18.0$ mag, and $\Delta$ $m_{15}$($B$)$\sim1.6$ mag), yet differ significantly from SNe 2008ha and 2010ae ($M_{J}\sim-14 - -15$ mag and $\Delta$ $m_{15}$($B$)$\sim2.4-2.7$ mag). The estimated rise time is $12.0 \pm 3.0$ days, which is significantly shorter than that of SN 2005hk or any other Ia SNe. The rapid rise indicates that the $^{56}$Ni distribution may extend into the outer layer or that the effective opacity may be lower than that in normal SNe Ia. The late-phase spectrum exhibits broader emission lines than those of SN 2005hk by a factor of 6--8. Such high velocities of the emission lines indicate that the density profile of the inner ejecta extends more than that of SN 2005hk. We argue that the most favored explosion scenario is a `failed deflagration' model, although the pulsational delayed detonations is not excluded.

We report on monitoring observations of the TeV gamma-ray binary HESS J0632+057, which were carried out to constrain the interaction between the Be circumstellar disk and the compact object of unknown nature, and provide for the first time high-dispersion (R > 50000) optical spectra in the second half of the orbital cycle, from apastron through periastron. The Halpha, Hbeta, and Hgamma line profiles are found to exhibit remarkable short-term variability for ~1 month after the apastron (phase 0.6--0.7), whereas they show little variation near the periastron. These emission lines show "S-shaped" variations with timescale of ~150 days, which is about twice that reported previously. In contrast to the Balmer lines, no profile variability is seen in any FeII emission line. We estimate the radii of emitting regions of the Halpha, Hbeta, Hgamma, and FeII emission lines to be ~30, 11, 7, and 2 stellar radii (R_*), respectively. The amplitudes of the line profile variations in different lines indicate that the interaction with the compact object affects the Be disk down to, at least, the radius of 7 R_* after the apastron. This fact, together with little profile variability near the periastron, rules out the tidal force as the major cause of disk variability. Although this leaves the pulsar wind as the most likely candidate mechanism for disk variations, understanding the details of the interaction, particularly the mechanism for causing a large-scale disk disturbance after the apastron, remains an open question.

We performed a systematic X-ray study of eight nearby $\gamma$-ray bright radio galaxies with {\em Suzaku} for understanding the origin of their X-ray emissions. The {\em Suzaku} spectra for five of those have been presented previously, while the remaining three (M\,87, PKS\,0625$-$354, and 3C\,78) are presented here for the first time. Based on the Fe-K line strength, X-ray variability, and X-ray power-law photon indices, and using additional information on the [O III] line emission, we argue for a jet origin of the observed X-ray emission in these three sources. We also analyzed five years of {\em Fermi} Large Area Telescope (LAT) GeV gamma-ray data on PKS\,0625$-$354 and 3C\,78 to understand these sources within the blazar picture. We found significant $\gamma$-ray variability in the former object. Overall, we note that the {\em Suzaku} spectra for both PKS\,0625$-$354 and 3C\,78 are rather soft, while the LAT spectra are unusually hard when compared with other $\gamma$-ray detected low-power (FR\,I) radio galaxies. We demonstrate that the constructed broad-band spectral energy distributions of PKS\,0625$-$354 and 3C\,78 are well described by a one-zone synchrotron/synchrotron self-Compton model. The results of the modeling indicate lower bulk Lorentz factors compared to those typically found in other BL Lac objects, but consistent with the values inferred from modeling other LAT-detected FR\,I radio galaxies. Interestingly, the modeling also implies very high peak ($\sim 10^{16}$\,Hz) synchrotron frequencies in the two analyzed sources, contrary to previously-suggested scenarios for FR I/BL Lac unification. We discuss the implications of our findings in the context of the FR\,I/BL Lac unification schemes.

We report the results of a multi-band observing campaign on the famous blazar 3C 279 conducted during a phase of increased activity from 2013 December to 2014 April, including first observations of it with NuSTAR. The $\gamma$-ray emission of the source measured by Fermi-LAT showed multiple distinct flares reaching the highest flux level measured in this object since the beginning of the Fermi mission, with $F(E > 100\,{\rm MeV})$ of $10^{-5}$ photons cm$^{-2}$ s$^{-1}$, and with a flux doubling time scale as short as 2 hours. The $\gamma$-ray spectrum during one of the flares was very hard, with an index of $\Gamma_\gamma = 1.7 \pm 0.1$, which is rarely seen in flat spectrum radio quasars. The lack of concurrent optical variability implies a very high Compton dominance parameter $L_\gamma/L_{\rm syn} > 300$. Two 1-day NuSTAR observations with accompanying Swift pointings were separated by 2 weeks, probing different levels of source activity. While the 0.5$-$70 keV X-ray spectrum obtained during the first pointing, and fitted jointly with Swift-XRT is well-described by a simple power law, the second joint observation showed an unusual spectral structure: the spectrum softens by $\Delta\Gamma_{\rm X} \simeq 0.4$ at $\sim$4 keV. Modeling the broad-band SED during this flare with the standard synchrotron plus inverse Compton model requires: (1) the location of the $\gamma$-ray emitting region is comparable with the broad line region radius, (2) a very hard electron energy distribution index $p \simeq 1$, (3) total jet power significantly exceeding the accretion disk luminosity $L_{\rm j}/L_{\rm d} \gtrsim 10$, and (4) extremely low jet magnetization with $L_{\rm B}/L_{\rm j} \lesssim 10^{-4}$. We also find that single-zone models that match the observed $\gamma$-ray and optical spectra cannot satisfactorily explain the production of X-ray emission.

We present six-year multi-wavelength monitoring result for broad-line radio galaxy 3C 120. The source was sporadically detected by Fermi-LAT and after the MeV/GeV gamma-ray detection the 43 GHz radio core brightened and a knot ejected from an unresolved core, implying that the radio-gamma phenomena are physically connected. We show that the gamma-ray emission region is located at sub-pc distance from the central black hole, and MeV/GeV gamma-ray emission mechanism is inverse-Compton scattering of synchrotron photons. We also discuss future perspective revealed by next-generation X-ray satellite Astro-H.

Broadband spectrum of AGN consists of multiple components such as jet emission and accretion disk emission. Temporal correlation study is useful to understand emission components and their physical origins. We have performed optical monitoring using Kanata telescope for 4 radio galaxies and 6 radio-loud Narrow-Line Seyfert 1 (RL-NLSy1): 2 gamma-ray-loud RL-NLSy1s, 1H 0323+342 and PMN J0948+0022, and 4 gamma-ray-quiet RL-NLSy1s. From these results, it is suggested that RL-NLSy1s show a disk-dominant phase and a jet-dominant phase in the optical band, but it is not well correlated with brightness.

1H0323+342 is one of narrow-line radio-loud Seyfert 1 galaxies (RL-NLS1), which is a new class of gamma-ray emitting AGNs. Narrow-line Seyfert 1 galaxies (NLS1) have a small-mass black hole, but its mass accretion rate is almost as high as Eddington limit. Therefore, by observing NLS1s, we can know the evolution of supermassive black holes at the center of galaxies. Some of NLS1s are radio-loud and we call them RL-NLS1. From past observations, multi-wavelength spectrum of RL-NLS1s is similar to that of typical blazars; the synchrotron emission in the lower energy band up to the optical band, and inverse Compton scattering of low energy photons from disk, torus, and broad line region. X-ray band is a transittion region between the synchrotron and inverse Compton, and also there is a possible disk/corona emission. Therefore, we studied the energy-dependence of time variability of the X-ray emission of 1H0323+342, which have been observed by Suzaku in 2009 and 2013, in order to constrain the emission mechanism. We found that the lower energy below 1 keV and the higher energy above 7 keV show a different variability from the middle energy band, indicating at least two emission components in the X-ray band. X-ray spectrum is not a simple power-law, but requires an additional features; a broken power-law plus flat hard component, or a power-law plus a relativistic reflection component. Each spectral component seems to vary independently.

NGC 1275 is a gamma-ray-emitting radio galaxy at the center of the Perseus cluster. Its multi-wavelength spectrum is similar to that of blazers, and thus a jet-origin of gamma-ray emissions is believed. In the optical and X-ray region, NGC 1275 also shows a bright core, but their origin has not been understood, since a disk emission is not ruled out. In fact, NGC 1275 exhibits optical broad emission lines and a X-ray Fe-K line, which are typical for Seyfert galaxies. In our precious studies of NGC 1275 with Suzaku/XIS, no X-ray time variability was found from 2006 to 2011, regardless of moderate gamma-ray variability observed by {it Fermi}-LAT~\cite{Yamazaki}. We have continued monitoring observations of NGC 1275 with Suzaku/XIS. In 2013-2014, MeV/GeV gams-ray flux of NGC 1275 gradually increased and reached the maximum at the beginning of 2014. Correlated with this recent gamma-ray activity, we found that X-ray flux also increased, and this is the first evidence of X-ray variability of NGC 1275. Following these results, we discuss the emission component during the time variability, but we cannot decide the origin of X-ray variability correlating with gamma-ray. Therefore, for future observation, it is important to observe NGC 1275 by using Fermi gamma-ray, XMM-Newton, NuStar, ASTRO-H X-ray, CTA TeV gamma-ray and Kanata optical telescope.

We report on long-term multi-wavelength monitoring of blazar Mrk~421 from 2010 to 2011. The source exhibited extreme X-ray flares in 2010. Our research group performed optical photopolarimetric follow-up observations using the Kanata telescope. In 2010, the variability in the X-ray band was significant, while the optical and ultraviolet (UV) flux decreased gradually. Polarization properties also exhibited unique variability in 2010, suggesting the presence of systematic component of polarization and magnetic field alignment for the emergence of a new polarized emission region. In contrast, in 2011 the variability in the X-ray band was smaller, and the variability in the optical and UV bands was larger, than in 2010. To explore the reasons for these differences, spectral fitting analysis was performed via simple synchrotron-self Compton modelling; the results revealed different behaviors in terms of spectral evolution between these periods, suggesting different variability mechanisms between 2010 and 2011. In 2010, the radiation was likely the result of energy injection into the emitting regions with an aligned magnetic field. In contrast, in 2011 the superposition of different emission regions may have contributed to the low degree of observed polarization. It also implies that high-energy electron which were not accelerated to ultra-relativistic velocities were injected in 2011.

Blazars are thought to possess a relativistic jet that is pointing toward the direction of the Earth and the elect of relativistic beaming enhances its apparent brightness. They radiate in all wavebands from the radio to the gamma-ray bands via the synchrotron and the inverse Compton scattering process. Numerous observations are performed but the mechanism of variability, creation and composition of jets are still controversial. We performed multi-wavelength monitoring with optical polarization for 3C 66A, Mrk 421, CTA 102 and PMN J0948+0022 to investigate the mechanisms of variability and research the emission region in the relativistic jets. Consequently, an emergence of new emission component in flaring state is suggested in each object. The most significant aspect of these results is its wide range of sizes of emission regions from $10^{14}-10^{16}$ cm, which implies the model with a number of independent emission regions with variety sizes and randomly orientation. The "shock-in-jet" scenario can explain high PD and direction of PA in each objects. It might reflect the common mechanism of flares in the relativistic jets.

We present multi-wavelength monitoring results for the broad-line radio galaxy 3C 120 in the MeV/GeV, sub-millimeter, and 43 GHz bands over six years. Over the past two years, Fermi-LAT sporadically detected 3C 120 with high significance and the 230 GHz data also suggest an enhanced activity of the source. After the MeV/GeV detection from 3C 120 in MJD 56240-56300, 43 GHz VLBA monitoring revealed a brightening of the radio core, followed by the ejection of a superluminal knot. Since we observed the gamma-ray and VLBA phenomena in temporal proximity to each other, it is naturally assumed that they are physically connected. This assumption was further supported by the subsequent observation that the 43 GHz core brightened again after a gamma-ray flare occurred around MJD 56560. We can then infer that the MeV/GeV emission took place inside an unresolved 43 GHz core of 3C 120 and that the jet dissipation occurred at sub-parsec distances from the central black hole, if we take the distance of the 43 GHz core from the central black hole as ~ 0.5 pc, as previously estimated from the time lag between X-ray dips and knot ejections (Marscher et al. 2002; Chatterjee et al. 2009). Based on our constraints on the relative locations of the emission regions and energetic arguments, we conclude that the gamma rays are more favorably produced via the synchrotron self-Compton process, rather than inverse Compton scattering of external photons coming from the broad line region or hot dusty torus. We also derived the electron distribution and magnetic field by modeling the simultaneous broadband spectrum.

We performed a systematic X-ray study of eight nearby gamma-ray bright radio galaxies with Suzaku for understanding the origin of their X-ray emissions. The Suzaku spectra for five of those have been presented previously, while the remaining three (M87, PKS0625-354, and 3C78) are presented here for the first time. Based on the Fe-K line strength, X-ray variability, and X-ray power-law photon indices, and using additional information on the [O III] line emission, we argue for a jet origin of the observed X-ray emission in these three sources. We also analyzed five years of Fermi Large Area Telescope (LAT) GeV gamma-ray data on PKS0625-354 and 3C78 to understand these sources within the blazar picture. We found significant gamma-ray variability in the former object. Overall, we note that the Suzaku spectra for both PKS0625-354 and 3C78 are rather soft, while the LAT spectra are unusually hard when compared with other gamma-ray detected low-power (FR I) radio galaxies. We demonstrate that the constructed broad-band spectral energy distributions of PKS0625-354 and 3C78 are well described by a one-zone synchrotron/synchrotron self-Compton model. The results of the modeling indicate lower bulk Lorentz factors compared to those typically found in other BL Lac objects, but consistent with the values inferred from modeling other LAT-detected FR\,I radio galaxies. Interestingly, the modeling also implies very high peak ($\sim 10^{16}$ Hz) synchrotron frequencies in the two analyzed sources, contrary to previously-suggested scenarios for FR I/BL Lac unification. We discuss the implications of our findings in the context of the FR I/BL Lac unification schemes.

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.

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

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

We present results of optical polarimetric and multi-band photometric observations for gamma-ray loud narrow-line Seyfert 1 galaxy 1H 0323+342. This object has been monitored by 1.5 m Kanata telescope since 2012 September but following a gamma-ray flux enhancement detected by Fermi-LAT on MJD 56483 (2013 July 10) dense follow-up was performed by ten 0.5-2.0 m telescopes in Japan over one week. The 2-year R_C-band light curve showed clear brightening corresponding to the gamma-ray flux increase and then decayed gradually. The high state as a whole lasted for ~20 days, during which we clearly detected optical polarization from this object. The polarization degree (PD) of the source increased from 0-1% in quiescence to ~3% at maximum and then declined to the quiescent level, with the duration of the enhancement of less than 10 days. The moderate PD around the peak allowed us to precisely measure the daily polarization angle (PA). As a result, we found that the daily PAs were almost constant and aligned to the jet axis, suggesting that the magnetic field direction at the emission region is transverse to the jet. This implies either a presence of helical/toroidal magnetic field or transverse magnetic field compressed by shock(s). We also found small-amplitude intra-night variability during the 2-hour continuous exposure on a single night. We discuss these findings based on the turbulent multi-zone model recently advocated by Marscher (2014). Optical to ultraviolet (UV) spectrum showed a rising shape in the higher frequency and the UV magnitude measured by Swift/UVOT was steady even during the flaring state, suggesting that thermal emission from accretion disk is dominant in that band.

Previously unremarkable, the extragalactic radio source GB 1310+487 showed a gamma-ray flare on 2009 November 18, reaching a daily flux of ~10^-6 photons/cm^2/s at energies E>100 MeV and becoming one of the brightest GeV sources for about two weeks. Its optical spectrum is not typical for a blazar, instead, it resembles those of narrow emission-line galaxies. We investigate changes of the object's radio-to-GeV spectral energy distribution (SED) during and after the prominent GeV flare with the aim to determine the nature of the object and constrain the origin of the variable high-energy emission. The data collected by the Fermi and AGILE satellites at gamma-ray energies, Swift at X-ray and ultraviolet, Kanata, NOT, and Keck telescopes at optical, OAGH and WISE at infrared, and IRAM 30m, OVRO 40m, Effelsberg 100m, RATAN-600, and VLBA at radio, are analysed together to trace the SED evolution on timescales of months. The gamma-ray/radio-loud narrow-line active galactic nucleus (AGN) is located at redshift z=0.638. It is shining through an unrelated foreground galaxy at z=0.500. The AGN light is likely amplified by a factor of a few because of gravitational lensing. The AGN SED shows a two-humped structure typical of blazars and gamma-ray-loud NLSy1 galaxies, with the high-energy (inverse-Compton) emission dominating by more than an order of magnitude over the low-energy (synchrotron) emission during gamma-ray flares. The difference between the two SED humps is smaller during the low-activity state. Fermi observations reveal a strong correlation between the gamma-ray flux and spectral index, with the hardest spectrum observed during the brightest gamma-ray state. If the gamma-ray flux is a mixture of synchrotron self-Compton (SSC) and external Compton (EC) emission, the observed GeV spectral variability may result from varying relative contributions of these two emission components.

The gamma-ray-detected blazar 3C 454.3 exhibits dramatic flux and polarization variations in the optical and near-infrared bands. In December 2010, the object emitted a very bright outburst. We monitored it for approximately four years (including the 2010 outburst) by optical and near-infrared photopolarimetry. During the 2010 outburst, the object emitted two rapid, redder brightenings, at which the polarization degrees (PDs) in both bands increased significantly and the bands exhibited a frequency-dependent polarization. The observed frequency-dependent polarization leads us to propose that the polarization vector is composed of two vectors. Therefore, we separate the observed polarization vectors into short and long-term components that we attribute to the emissions of the rapid brightenings and the outburst that varied the timescale of days and months, respectively. The estimated PD of the short-term component is greater than the maximum observed PD and is close to the theoretical maximum PD. We constrain the bulk Lorentz factors and inclination angles between the jet axis and the line of sight from the estimated PDs. In this case, the inclination angle of the emitting region of short-term component from the first rapid brightening should be equal to 90$^{\circ}$, because the estimated PD of the short-term component was approximately equal to the theoretical maximum PD. Thus, the Doppler factor at the emitting region of the first rapid brightening should be equal to the bulk Lorentz factor.

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.

We present the results of photometric (V band) and polarimetric observations of the blazar BL Lac during 2008--2010 using TRISPEC attached to the KANATA 1.5-m telescope in Japan. The data reveal a great deal of variability ranging from days to months with detection of strong variations in fractional polarization. The V band flux strongly anti-correlates with the degree of polarization during the first of two observing seasons but not during the second. The direction of the electric vector, however, remained roughly constant during all our observations. These results are consistent with a model with at least two emission regions being present, with the more variable component having a polarization direction nearly perpendicular to that of the relatively quiescent region so that a rising flux can produce a decline in degree of polarization. We also computed models involving helical jet structures and single transverse shocks in jets and show that they might also be able to agree with the anti-correlations between flux and fractional polarization.

Quasar 3C 279 is known to exhibit episodes of optical polarization angle rotation. We present new, well-sampled optical polarization data for 3C 279 and introduce a method to distinguish between random and deterministic electric vector position angle (EVPA) variations. We observe EVPA rotations in both directions with different amplitudes and find that the EVPA variation shows characteristics of both random and deterministic cases. Our analysis indicates that the EVPA variation is likely dominated by a random process in the low brightness state of the jet and by a deterministic process in the flaring state.

We report on optical photopolarimetric results of the radio-loud narrow line Seyfert 1 (RL-NLSy1) galaxy PMN J0948+0022 on 2012 December to 2013 February triggered by flux enhancements in near infrared and gamma-ray bands. Thanks to one-shot polarimetry of the HOWPol installed to the Kanata telescope, we have detected very rapid variability in the polarized-flux light curve on MJD 56281 (2012 December 20). The rise and decay times were about 140 sec and 180 sec, respectively. The polarization degree (PD) reached 36 +/- 3% at the peak of the short-duration pulse, while polarization angle (PA) remained almost constant. In addition, temporal profiles of the total flux and PD showed highly variable but well correlated behavior and discrete correlation function analysis revealed that no significant time lag of more than 10 min was present. The high PD and minute-scale variability in polarized flux provides a clear evidence of synchrotron radiation from a very compact emission region of 10^14 cm size with highly ordered magnetic field. Such micro variability of polarization are also observed in several blazar jets, but its complex relation between total flux and PD are explained by multi-zone model in several blazars. The implied single emission region in PMN J0948+0022 might be reflecting a difference of jets between RL-NLSy1s and blazars.

We present a set of photometric and spectroscopic observations of a bright Type Ib supernova SN 2012au from -6d until ~+150d after maximum. The shape of its early R-band light curve is similar to that of an average Type Ib/c supernova. The peak absolute magnitude is M_R=-18.7+-0.2 mag, which suggests that this supernova belongs to a very luminous group among Type Ib supernovae. The line velocity of He I {\lambda}5876 is about 15,000 km/s around maximum, which is much faster than that in a typical Type Ib supernova. From the quasi-bolometric peak luminosity of (6.7+-1.3)x10^(42) erg/s, we estimate the \Ni mass produced during the explosion as ~0.30 Msun. We also give a rough constraint to the ejecta mass 5-7 Msun and the kinetic energy (7-18)x10^(51) erg. We find a weak correlation between the peak absolute magnitude and He I velocity among Type Ib SNe. The similarities to SN 1998bw in the density structure inferred from the light curve model as well as the large peak bolometric luminosity suggest that SN 2012au had properties similar to energetic Type Ic supernovae.

CTA 102, classified as a flat spectrum radio quasar at z=1.037, produced exceptionally bright optical flare in 2012 September. Following Fermi-LAT detection of enhanced gamma-ray activity, we densely monitored this source in the optical and near-infrared bands for the subsequent ten nights using twelve telescopes in Japan and South-Africa. On MJD 56197 (2012 September 27, 4-5 days after the peak of bright gamma-ray flare), polarized flux showed a transient increase, while total flux and polarization angle remained almost constant during the "orphan polarized-flux flare". We also detected an intra-night and prominent flare on MJD 56202. The total and polarized fluxes showed quite similar temporal variations, but PA again remained constant during the flare. Interestingly, the polarization angles during the two flares were significantly different from the jet direction. Emergence of a new emission component with high polarization degree (PD) up to 40% would be responsible for the observed two flares, and such a high PD indicates a presence of highly ordered magnetic field at the emission site. We discuss that the well-ordered magnetic field and even the observed directions of polarization angle which is grossly perpendicular to the jet are reasonably accounted for by transverse shock(s) propagating down the jet.

We present a study of SN 2009js in NGC 918. Multi-band Kanata optical photometry covering the first ~120 days show the source to be a Type IIP SN. Reddening is dominated by that due to our Galaxy. One-year-post-explosion photometry with the NTT, and a Subaru optical spectrum 16 days post-discovery, both imply a good match with the well-studied subluminous SN 2005cs. The plateau phase luminosity of SN 2009js and its plateau duration are more similar to the intermediate luminosity IIP SN 2008in. Thus, SN 2009js shares characteristics with both subluminous and intermediate luminosity SNe. Its radioactive tail luminosity lies between SN 2005cs and SN 2008in, whereas its quasi-bolometric luminosity decline from peak to plateau (quantified by a newly-defined parameter Delta[logL] measuring adiabatic cooling following shock breakout) is much smaller than both the others. We estimate the ejected mass of 56Ni to be low (~0.007 Msun). The SN explosion energy appears to have been small, similar to SN 2005cs. SN 2009js is the first subluminous SN IIP to be studied in the mid-infrared. It was serendipitously caught by Spitzer at very early times. In addition, it was detected by WISE 105 days later with a significant 4.6 micron flux excess above the photosphere. The infrared excess luminosity relative to the photosphere is clearly smaller than that of SN 2004dj extensively studied in the mid-infrared. The excess may be tentatively assigned to heated dust with mass ~3e-5 Msun, or to CO fundamental emission as a precursor to dust formation.

We present the results of a series of radio, optical, X-ray and gamma-ray observations of the BL Lac object S50716+714 carried out between April 2007 and January 2011. The multi-frequency observations were obtained using several ground and space based facilities. The intense optical monitoring of the source reveals faster repetitive variations superimposed on a long-term variability trend at a time scale of ~350 days. Episodes of fast variability recur on time scales of ~ 60-70 days. The intense and simultaneous activity at optical and gamma-ray frequencies favors the SSC mechanism for the production of the high-energy emission. Two major low-peaking radio flares were observed during this high optical/gamma-ray activity period. The radio flares are characterized by a rising and a decaying stage and are in agreement with the formation of a shock and its evolution. We found that the evolution of the radio flares requires a geometrical variation in addition to intrinsic variations of the source. Different estimates yield a robust and self-consistent lower limits of \delta > 20 and equipartition magnetic field B_eq > 0.36 G. Causality arguments constrain the size of emission region \theta < 0.004 mas. We found a significant correlation between flux variations at radio frequencies with those at optical and gamma-rays. The optical/GeV flux variations lead the radio variability by ~65 days. The longer time delays between low-peaking radio outbursts and optical flares imply that optical flares are the precursors of radio ones. An orphan X-ray flare challenges the simple, one-zone emission models, rendering them too simple. Here we also describe the spectral energy distribution modeling of the source from simultaneous data taken through different activity periods.

3C 66A is an intermediate-frequency-peaked BL Lac object detected by the Large Area Telescope onboard the Fermi Gamma-ray Space Telescope. We present a study of the long-term variations of this blazar seen over 2 years at GeV energies with Fermi and in the optical (flux and polarization) and near infrared with the Kanata telescope. In 2008, the first year of the study, we find a correlation between the gamma-ray flux and the measurements taken with the Kanata telescope. This is in contrast to the later measurements performed during 2009--2010 which show only a weak correlation along with a gradual increase of the optical flux. We calculate an external seed photon energy density assuming that the gamma-ray emission is due to external Compton scattering. The energy density of the external photons is found to be higher by a factor of two in 2008 compared to 2009--2010. We conclude that the different behaviors observed between the first year and the later years might be explained by postulating two different emission components.

The blazar AO 0235+164 (z = 0.94) has been one of the most active objects observed by Fermi Large Area Telescope (LAT) since its launch in Summer 2008. In addition to the continuous coverage by Fermi, contemporaneous observations were carried out from the radio to {\gamma} -ray bands between 2008 September and 2009 February. In this paper, we summarize the rich multi-wavelength data collected during the campaign (including F-GAMMA, GASP- WEBT, Kanata, OVRO, RXTE, SMARTS, Swift, and other instruments), examine the cross-correlation between the light curves measured in the different energy bands, and interpret the resulting spectral energy distributions in the context of well-known blazar emission models. We find that the {\gamma} -ray activity is well correlated with a series of near-IR/optical flares, accompanied by an increase in the optical polarization degree. On the other hand, the X-ray light curve shows a distinct 20 day high state of unusually soft spectrum, which does not match the extrapolation of the optical/UV synchrotron spectrum. We tentatively interpret this feature as the bulk Compton emission by cold electrons contained in the jet, which requires an accretion disk corona with an effective covering factor of 19% at a distance of 100 Rg . We model the broadband spectra with a leptonic model with external radiation dominated by the infrared emission from the dusty torus.

We present time-resolved broad-band observations of the quasar 3C 279 obtained from multi-wavelength campaigns conducted during the first two years of the Fermi Gamma-ray Space Telescope mission. While investigating the previously reported gamma-ray/optical flare accompanied by a change in optical polarization, we found that the optical emission appears delayed with respect to the gamma-ray emission by about 10 days. X-ray observations reveal a pair of `isolated' flares separated by ~90 days, with only weak gamma-ray/optical counterparts. The spectral structure measured by Spitzer reveals a synchrotron component peaking in the mid-infrared band with a sharp break at the far-infrared band during the gamma-ray flare, while the peak appears in the mm/sub-mm band in the low state. Selected spectral energy distributions are fitted with leptonic models including Comptonization of external radiation produced in a dusty torus or the broad-line region. Adopting the interpretation of the polarization swing involving propagation of the emitting region along a curved trajectory, we can explain the evolution of the broad-band spectra during the gamma-ray flaring event by a shift of its location from ~ 1 pc to ~ 4 pc from the central black hole. On the other hand, if the gamma-ray flare is generated instead at sub-pc distance from the central black hole, the far-infrared break can be explained by synchrotron self-absorption. We also model the low spectral state, dominated by the mm/sub-mm peaking synchrotron component, and suggest that the corresponding inverse-Compton component explains the steady X-ray emission.

We report that the optical polarization in the afterglow of GRB 091208B is measured at t = 149 - 706 s after the burst trigger, and the polarization degree is P = 10.4% +/- 2.5%. The optical light curve at this time shows a power-law decay with index -0.75 +/- 0.02, which is interpreted as the forward shock synchrotron emission, and thus this is the first detection of the early-time optical polarization in the forward shock (rather than that in the reverse shock reported by Steele et al. (2009). This detection disfavors the afterglow model in which the magnetic fields in the emission region are random on the plasma skin depth scales, such as amplified by the plasma instabilities, e.g., Weibel instability. We suggest that the fields are amplified by the magnetohydrodynamic instabilities, which would be tested by future observations of the temporal changes of the polarization degrees and angles for other bursts.

In December 2009, the bright blazar, 3C 454.3 exhibited a strong outburst in the optical, X-ray and gamma-ray regions. We performed photometric and polarimetric monitoring of this outburst in the optical and near-infrared bands with TRISPEC and HOWPol attached to the Kanata telescope. We also observed this outburst in the infrared band with AKARI, and the radio band with the 32-m radio telescope of Yamaguchi University. The object was in an active state from JD 2455055 to 2455159. It was 1.3 mag brighter than its quiescent state before JD 2455055 in the optical band. After the end of the active state in JD 2455159, a prominent outburst was observed in all wavelengths. The outburst continued for two months. Our optical and nearinfrared polarimetric observations revealed that the position angle of the polarization (PA) apparently rotated clockwise by 240 degrees within 11 d in the active state (JD 2455063-2455074), and after this rotation, PA remained almost constant during our monitoring. In the outburst state, PA smoothly rotated counterclockwise by 350 degrees within 35 d (JD 2455157-2455192). Thus, we detected two distinct rotation events of polarization vector in opposite directions. We discuss these two events compared with the past rotation events observed in 2005, 2007 and 2008.

We report on multi-band photometric and polarimetric observations of the blazars AO 0235+164 and PKS 1510-089. These two blazars were active in 2008 and 2009, respectively. In these active states, prominent short flares were observed in both objects, having amplitudes of >1 mag within 10 d. The $V-J$ color became bluer when the objects were brighter in these flares. On the other hand, the color of PKS 1510-089 exhibited a trend that it became redder when it was brighter, except for its prominent flare. This redder-when-brighter trend can be explained by the strong contribution of thermal emission from an accretion disk. The polarization degree increased at the flares, and reached >25 % at the maxima. We compare these flares in AO 0235+164 and PKS 1510-089 with other short flares which were detected by our monitoring of 41 blazars. Those two flares had one of the largest variation amplitudes in both flux and polarization degree. Furthermore, we found a significant positive correlation between the amplitudes of the flux and polarization degree in the short flares. It indicates that the short flares originate from the region where the magnetic field is aligned.

We present early spectroscopy of the recurrent nova U~Sco during the outburst in 2010. We successfully obtained time-series spectra at $t_{\rm d}=$0.37--0.44~d, where $t_{\rm d}$ denotes the time from the discovery of the present outburst. This is the first time-resolved spectroscopy on the first night of U Sco outbursts. At $t_{\rm d}\sim 0.4$~d the H$\alpha$ line consists of a blue-shifted ($-5000$ km s$^{-1}$) narrow absorption component and a wide emission component having triple peaks, a blue ($\sim -3000$ km s$^{-1}$), a central ($\sim 0$ km s$^{-1}$) and a red ($\sim +3000$ km s$^{-1}$) ones. The blue and red peaks developed more rapidly than the central one during the first night. This rapid variation would be caused by the growth of aspherical wind produced during the earliest stage of the outburst. At $t_{\rm d}=1.4$~d the H$\alpha$ line has a nearly flat-topped profile with weak blue and red peaks at $\sim \pm 3000$ km s$^{-1}$. This profile can be attributed to a nearly spherical shell, while the asphericity growing on the first night still remains. The wind asphericity is less significant after $t_{\rm d}=9$ d.

During the month of December, 2009 the blazar 3C 454.3 became the brightest gamma-ray source in the sky, reaching a peak flux F ~2000E-8 ph/cm2/s for E > 100 MeV. Starting in November, 2009 intensive multifrequency campaigns monitored the 3C 454 gamma-ray outburst. Here we report the results of a 2-month campaign involving AGILE, INTEGRAL, Swift/XRT, Swift/BAT, RossiXTE for the high-energy observations, and Swift/UVOT, KANATA, GRT, REM for the near-IR/optical/UV data. The GASP/WEBT provided radio and additional optical data. We detected a long-term active emission phase lasting ~1 month at all wavelengths: in the gamma-ray band, peak emission was reached on December 2-3, 2009. Remarkably, this gamma-ray super-flare was not accompanied by correspondingly intense emission in the optical/UV band that reached a level substantially lower than the previous observations in 2007-2008. The lack of strong simultaneous optical brightening during the super-flare and the determination of the broad-band spectral evolution severely constrain the theoretical modelling. We find that the pre- and post-flare broad-band behavior can be explained by a one-zone model involving SSC plus external Compton emission from an accretion disk and a broad-line region. However, the spectra of the Dec. 2-3, 2009 super-flare and of the secondary peak emission on Dec. 9, 2009 cannot be satisfactorily modelled by a simple one-zone model. An additional particle component is most likely active during these states.

(Abridged) We have conducted a detailed investigation of the broad-band spectral properties of the \gamma-ray selected blazars of the Fermi LAT Bright AGN Sample (LBAS). By combining our accurately estimated Fermi gamma-ray spectra with Swift, radio, infra-red, optical and other hard X-ray/gamma-ray data, collected within three months of the LBAS data taking period, we were able to assemble high-quality and quasi-simultaneous Spectral Energy Distributions (SED) for 48 LBAS blazars.The SED of these gamma-ray sources is similar to that of blazars discovered at other wavelengths, clearly showing, in the usual Log $\nu $ - Log $\nu$ F$_\nu$ representation, the typical broad-band spectral signatures normally attributed to a combination of low-energy synchrotron radiation followed by inverse Compton emission of one or more components. We have used these SEDs to characterize the peak intensity of both the low and the high-energy components. The results have been used to derive empirical relationships that estimate the position of the two peaks from the broad-band colors (i.e. the radio to optical and optical to X-ray spectral slopes) and from the gamma-ray spectral index. Our data show that the synchrotron peak frequency $\nu_p^S$ is positioned between 10$^{12.5}$ and 10$^{14.5}$ Hz in broad-lined FSRQs and between $10^{13}$ and $10^{17}$ Hz in featureless BL Lacertae objects.We find that the gamma-ray spectral slope is strongly correlated with the synchrotron peak energy and with the X-ray spectral index, as expected at first order in synchrotron - inverse Compton scenarios. However, simple homogeneous, one-zone, Synchrotron Self Compton (SSC) models cannot explain most of our SEDs, especially in the case of FSRQs and low energy peaked (LBL) BL Lacs. (...)

We developed a method to separate a long-term trend from observed temporal variations of polarization in blazars using a Bayesian approach. The temporal variation of the polarization vector is apparently erratic in most blazars, while several objects occasionally exhibited systematic variations, for example, an increase of the polarization degree associated with a flare of the total flux. We assume that the observed polarization vector is a superposition of distinct two components, a long-term trend and a short-term variation component responsible for short flares. Our Bayesian model estimates the long-term trend which satisfies the condition that the total flux correlates with the polarized flux of the short-term component. We demonstrate that assumed long-term polarization components are successfully separated by the Bayesian model for artificial data. We applied this method to photopolarimetric data of OJ 287, S5 0716+714, and S2 0109+224. Simple and systematic long-term trends were obtained in OJ 287 and S2 0109+224, while no such a trend was identified in S5 0716+714. We propose that the apparently erratic variations of polarization in OJ 287 and S2 0109+224 are due to the presence of the long-term polarization component. The behavior of polarization in S5 0716+714 during our observation period implies the presence of a number of polarization components having a quite short time-scale of variations.

We present early phase observations in optical and near-infrared wavelengths for the extremely luminous Type Ia supernova (SN Ia) 2009dc. The decline rate of the light curve is $\Delta m_{15}(B)=0.65\pm 0.03$, which is one of the slowest among SNe Ia. The peak $V$-band absolute magnitude is $M_{V}=-19.90\pm 0.15$ mag even if the host extinction is $A_{V}=0$ mag. It reaches $M_{V}=-20.19\pm 0.19$ mag for the host extinction of $A_{V}=0.29$ mag as inferred from the observed Na {\sc i} D line absorption in the host. Our $JHK_{s}$-band photometry shows that the SN is one of the most luminous SNe Ia also in near-infrared wavelengths. These results indicate that SN 2009dc belongs to the most luminous class of SNe Ia, like SN 2003fg and SN 2006gz. We estimate the ejected $^{56}$Ni mass of $1.2\pm 0.3$ $\Msun$ for no host extinction case (or 1.6$\pm$ 0.4 M$_{\odot}$ for the host extinction of $A_{V}=0.29$ mag). The C {\sc ii} $\lambda$6580 absorption line keeps visible until a week after maximum, which diminished in SN 2006gz before its maximum brightness. The line velocity of Si {\sc ii} $\lambda$6355 is about 8000 km s$^{-1}$ around the maximum, being considerably slower than that of SN 2006gz, while comparable to that of SN 2003fg. The velocity of the C {\sc ii} line is almost comparable to that of the Si {\sc ii}. The presence of the carbon line suggests that thick unburned C+O layers remain after the explosion. SN 2009dc is a plausible candidate of the super-Chandrasekhar mass SNe Ia.