Blazars can be divided into two subtypes, flat spectrum radio quasars (FSRQs) and BL Lac objects, which have been distinguished phenomenologically by the strength of their optical emission lines, while their physical nature and relationship are still not fully understood. In this paper, we focus on the differences in their variability. We characterize the blazar variability using the Ornstein-Uhlenbeck (OU) process, and investigate the features that are discriminative for the two subtypes. We used optical photometric and polarimetric data obtained with the 1.5-m Kanata telescope for 2008-2014. We found that four features, namely the variation amplitude, characteristic timescale, and non-stationarity of the variability obtained from the light curves and the median of the degree of polarization (PD), are essential for distinguishing between FSRQs and BL Lac objects. FSRQs are characterized by rare and large flares, while the variability of BL Lac objects can be reproduced with a stationary OU process with relatively small amplitudes. The characteristics of the variability are governed not by the differences in the jet structure between the subtypes, but by the peak frequency of the synchrotron emission. This implies that the nature of the variation in the jets is common in FSRQs and BL Lac objects. We found that BL Lac objects tend to have high PD medians, which suggests that they have a stable polarization component. FSRQs have no such component, possibly because of a strong Compton cooling effect in sub-pc scale jets.

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

Japanese Venus Climate Orbiter/AKATSUKI was proposed in 2001 with strong support by international Venus science community and approved as an ISAS (The Institute of Space and Astronautical Science) mission soon after the proposal. The mission life we expected was more than two Earth years in Venus orbit. AKATSUKI was successfully launched at 06:58:22JST on May 21, 2010, by H-IIA F17. After the separation from H-IIA, the telemetry from AKATSUKI was normally detected by DSN Goldstone station (10:00JST) and the solar cell paddles' deployment was confirmed. After a successful cruise, the malfunction happened on the propulsion system during the Venus orbit insertion (VOI) on Dec 7, 2010. The engine shut down before the planned reduction in speed to achieve. The spacecraft did not enter the Venus orbit, but entered an orbit around the Sun with a period of 203 days. Most of the fuel still had remained, but the orbital maneuvering engine was found to be broken and unusable. However, we have found an alternate way of achieving orbit by using only the reaction control system (RSC). We had adopted the alternate way for orbital maneuver and three minor maneuvers in Nov 2011 were successfully done so that AKATSUKI would meet Venus in 2015. We are considering several scenarios for VOI using only RCS.

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

We present 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 have observed the lightcurves of 13 V-type asteroids ((1933) Tinchen, (2011) Veteraniya, (2508) Alupka, (3657) Ermolova, (3900) Knezevic, (4005) Dyagilev, (4383) Suruga, (4434) Nikulin, (4796) Lewis, (6331) 1992 $\mathrm{FZ_{1}}$, (8645) 1998 TN, (10285) Renemichelsen, and (10320) Reiland). Using these observations we determined the rotational rates of the asteroids, with the exception of Nikulin and Renemichelsen. The distribution of rotational rates of 59 V-type asteroids in the inner main belt, including 29 members of the Vesta family that are regarded as ejecta from the asteroid (4) Vesta, is inconsistent with the best-fit Maxwellian distribution. This inconsistency may be due to the effect of thermal radiation Yarkovsky--O'Keefe--Radzievskii--Paddack (YORP) torques, and implies that the collision event that formed V-type asteroids is sub-billion to several billion years in age.

We propose the Hyper-Kamiokande (Hyper-K) detector as a next generation underground water Cherenkov detector. It will serve as a far detector of a long baseline neutrino oscillation experiment envisioned for the upgraded J-PARC, and as a detector capable of observing -- far beyond the sensitivity of the Super-Kamiokande (Super-K) detector -- proton decays, atmospheric neutrinos, and neutrinos from astronomical origins. The baseline design of Hyper-K is based on the highly successful Super-K, taking full advantage of a well-proven technology. (to be continued)

We present the yield and spectra of stable secondary particles (gamma, e^+/-, nu_e, nubar_e, nu_mu, and nubar_mu) of p-p interaction in parameterized formulae to facilitate calculations involving them in astronomical environments. The formulae are derived on the up-to-date p-p interaction model by [Kamae05] which incorporates the logarithmically rising inelastic cross-section, the diffraction dissociation process, and the Feynman scaling violation. To improve fidelity to experimental data in lower energies, two baryon resonance contributions have been added: one representing Delta(1232) and the other multiple resonances around 1600MeV/c^2. The parametrized formulae predict that all secondary particle spectra be harder by about 0.05 in power-law indices than that of the incident proton and their inclusive cross-sections be larger than those predicted by p-p interaction models based on the Feynman scaling.

We present here a new calculation of the gamma-ray spectrum from pp->pi^0 in the Galactic ridge environment. The calculation includes the diffractive pp interaction and incorporates the Feynman scaling violation for the first time. Galactic diffuse gamma-rays come, predominantly, from pi^0->gamma gamma in the sub-GeV to multi-GeV range. Hunter et al. found, however, an excess in the GeV range ("GeV Excess") in the EGRET Galactic diffuse spectrum above the prediction based on experimental pp->pi^0 cross-sections and the Feynman scaling hypothesis. We show, in this work, that the diffractive process makes the gamma-ray spectrum harder than the incident proton spectrum by ~0.05 in power-law index, and, that the scaling violation produces 30-80% more pi^0 than the scaling model for incident proton energies above 100GeV. Combination of the two can explain about a half of the "GeV Excess" with the local cosmic proton (power-law index ~2.7). The excess can be fully explained if the proton spectral index in the Galactic ridge is a little harder (~0.2 in power-law index) than the local spectrum. Given also in the paper is that the diffractive process enhances e^+ over e^- and the scaling violation gives 50-100% higher p-bar yield than without the violation, both in the multi-GeV range.