We present optical multi-colour photometry of V404 Cyg during the outburst from December, 2015 to January, 2016 together with the simultaneous X-ray data. This outburst occurred less than 6 months after the previous outburst in June-July, 2015. These two outbursts in 2015 were of a slow rise and rapid decay-type and showed large-amplitude ($\sim$2 mag) and short-term ($\sim$10 min-3 hours) optical variations even at low luminosity (0.01-0.1$L_{\rm Edd}$). We found correlated optical and X-ray variations in two $\sim$1 hour time intervals and performed Bayesian time delay estimations between them. In the previous version, the observation times of X-ray light curves were measured at the satellite and their system of times was Terrestrial Time (TT), while those of optical light curves were measured at the Earth and their system of times was Coordinated Universal Time (UTC). In this version, we have corrected the observation times and obtained a Bayesian estimate of an optical delay against the X-ray emission, which is $\sim$30 s, during those two intervals. In addition, the relationship between the optical and X-ray luminosity was $L_{\rm opt} \propto L_{\rm X}^{0.25-0.29}$ at that time. These features can be naturally explained by disc reprocessing.

How black holes accrete surrounding matter is a fundamental, yet unsolved question in astrophysics. It is generally believed that matter is absorbed into black holes via accretion disks, the state of which depends primarily on the mass-accretion rate. When this rate approaches the critical rate (the Eddington limit), thermal instability is supposed to occur in the inner disc, causing repetitive patterns of large-amplitude X-ray variability (oscillations) on timescales of minutes to hours. In fact, such oscillations have been observed only in sources with a high mass accretion rate, such as GRS 1915+105. These large-amplitude, relatively slow timescale, phenomena are thought to have physical origins distinct from X-ray or optical variations with small amplitudes and fast ($\lesssim$10 sec) timescales often observed in other black hole binaries (e.g., XTE J1118+480 and GX 339-4). Here we report an extensive multi-colour optical photometric data set of V404 Cygni, an X-ray transient source containing a black hole of nine solar masses (and a conpanion star) at a distance of 2.4 kiloparsecs. Our data show that optical oscillations on timescales of 100 seconds to 2.5 hours can occur at mass-accretion rates more than ten times lower than previously thought. This suggests that the accretion rate is not the critical parameter for inducing inner-disc instabilities. Instead, we propose that a long orbital period is a key condition for these large-amplitude oscillations, because the outer part of the large disc in binaries with long orbital periods will have surface densities too low to maintain sustained mass accretion to the inner part of the disc. The lack of sustained accretion -- not the actual rate -- would then be the critical factor causing large-amplitude oscillations in long-period systems.

We carried out the photometric observations of the SU UMa-type dwarf nova ER UMa during 2011 and 2012, which showed the existence of persistent negative superhumps even during the superoutburst. We performed two-dimensional period analysis of its light curves by using a method called "least absolute shrinkage and selection operator" (Lasso) and "phase dispersion minimization" (PDM) analysis, and we found that the period of negative superhumps systematically changed between a superoutburst and the next superoutburst. The trend of the period change can beinterpreted as reflecting the change of the disk radius. This change of the disk radius is in good agreement with the predicted change of the disk radius by the thermal-tidal instability (TTI) model. The normal outbursts within a supercycle showed a general trend that the rising rate to maximum becomes slower as the next superoutburst approaches. The change can be interpreted as the consequence of the increased gas-stream flow onto the inner region of the disk as the result of the tilted disk. Some of the superoutbursts were found to be triggered by a precursor normal outburst when the positive superhumps appeared to develop. The positive and negative superhumps co-existed during the superoutburst. The positive superhumps were prominent only during four or five days after the supermaximum, while the signal of the negative superhumps became strong after the middle phase of the superoutburst plateau. A simple combination of the positive and negative superhumps was found to be insufficient in reproducing the complex profile variation. We were able to detect the developing phase of positive superhumps (stage A superhumps) for the first time in ER UMa-type dwarf novae. Using the period of stage A superhumps, we obtained a mass ratio of 0.100(15), which indicates that ER UMa is on the ordinary evolutional track of CVs.

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

Continuing the project described by Kato et al. (2009, arXiv:0905.1757), we studied 86 SU UMa-type dwarf novae. We confirmed the general trends such as the relation between period derivatives and orbital periods. There are some systems showing positive period derivatives despite the long orbital periods. We observed the 2011 outburst of the WZ Sge-type dwarf nova BW Scl, and recorded an O-C diagram similar to those of previously known WZ Sge-type dwarf novae. The WZ Sge-type dwarf nova OT J184228.1+483742 showed an unusual pattern of double outbursts composed of an outburst with early superhumps and one with ordinary superhumps. We propose an interpretation that a very small growth rate of the 3:1 resonance due to an extremely low mass-ratio led to a quenching of the superoutburst before the ordinary superhumps appeared. We studied ER UMa-type dwarf novae and found that V1159 Ori showed positive superhumps similar to ER UMa in the 1990s. The recently recognized ER UMa-type object BK Lyn dominantly showed negative superhumps, and its behavior was very similar to the present-day state of ER UMa. The pattern of period variations in AM CVn-type objects were very similar to short-period hydrogen-rich SU UMa-type dwarf novae, making them helium analogue of hydrogen-rich SU UMa-type dwarf novae. SBS 1108+574, a peculiar hydrogen-rich dwarf nova below the period minimum, showed a very similar pattern of period variations to those of short-period SU UMa-type dwarf novae. The mass-ratio derived from the detected orbital period suggests that this secondary is a somewhat evolved star whose hydrogen envelope was mostly stripped during the mass-exchange. CC Scl, MASTER OT J072948.66+593824.4 and OT J173516.9+154708 showed only low-amplitude superhumps with complex profiles. These superhumps are likely a combination of closely separated two periods.

Continuing the project described by Kato et al. (2009, PASJ 61, S395, arXiv:0905.1757), we collected times of superhump maxima for 51 SU UMa-type dwarf novae mainly observed during the 2010-2011 season. Although most of the new data for systems with short superhump periods basically confirmed the findings by Kato et al. (2009) and Kato et al. (2010, PASJ 62, 1525, arXiv:1009.5444), the long-period system GX Cas showed an exceptionally large positive period derivative. An analysis of public Kepler data of V344 Lyr and V1504 Cyg yielded less striking stage transitions. In V344 Lyr, there was prominent secondary component growing during the late stage of superoutbursts, and the component persisted at least for two more cycles of successive normal outbursts. We also investigated the superoutbursts of two conspicuous eclipsing objects: HT Cas and the WZ Sge-type object SDSS J080434.20+510349.2. Strong beat phenomena were detected in both objects, and late-stage superhumps in the latter object had an almost constant luminosity during the repeated rebrightenings. The WZ Sge-type object SDSS J133941.11+484727.5 showed a phase reversal around the rapid fading from the superoutburst. The object showed a prominent beat phenomenon even after the end of the superoutburst. A pilot study of superhump amplitudes indicated that the amplitudes of superhumps are strongly correlated with orbital periods, and the dependence on the inclination is weak in systems with inclinations smaller than 80 deg.

As an extension of the project in Kato et al. (2009, arXiv:0905.1757), we collected times of superhump maxima for 61 SU UMa-type dwarf novae mainly observed during the 2009-2010 season. The newly obtained data confirmed the basic findings reported in Kato et al. (2009): the presence of stages A-C, as well as the predominance of positive period derivatives during stage B in systems with superhump periods shorter than 0.07 d. There was a systematic difference in period derivatives for systems with superhump periods longer than 0.075 d between this study and Kato et al. (2009). We suggest that this difference is possibly caused by the relative lack of frequently outbursting SU UMa-type dwarf novae in this period regime in the present study. We recorded a strong beat phenomenon during the 2009 superoutburst of IY UMa. The close correlation between the beat period and superhump period suggests that the changing angular velocity of the apsidal motion of the elliptical disk is responsible for the variation of superhump periods. We also described three new WZ Sge-type objects with established early superhumps and one with likely early superhumps. We also suggest that two systems, VX For and EL UMa, are WZ Sge-type dwarf novae with multiple rebrightenings. The O-C variation in OT J213806.6+261957 suggests that the frequent absence of rebrightenings in very short-Porb objects can be a result of sustained superoutburst plateau at the epoch when usual SU UMa-type dwarf novae return to quiescence preceding a rebrightening. We also present a formulation for a variety of Bayesian extension to traditional period analyses.

During the period when the Sun was intensely active on October-November 2003, two remarkable solar neutron events were observed by the ground-based neutron monitors. On October 28, 2003, in association with an X17.2 large flare, solar neutrons were detected with high statistical significance (6.4 sigma) by the neutron monitor at Tsumeb, Namibia. On November 4, 2003, in association with an X28 class flare, relativistic solar neutrons were observed by the neutron monitors at Haleakala in Hawaii and Mexico City, and by the solar neutron telescope at Mauna Kea in Hawaii simultaneously. Clear excesses were observed at the same time by these detectors, with the significance calculated as 7.5 sigma for Haleakala, and 5.2 sigma for Mexico City. The detector onboard the INTEGRAL satellite observed a high flux of hard X-rays and gamma-rays at the same time in these events. By using the time profiles of the gamma-ray lines, we can explain the time profile of the neutron monitor. It appears that neutrons were produced at the same time as the gamma-ray emission.

The black hole X-ray binary V4641 Sgr experienced an outburst in 2002 May which was detected at X-ray, optical, and radio wavelengths. The outburst lasted for only 6 days, but the object remained active for the next several months. Here we report on the detailed properties of light curves during the outburst and the post-outburst active phase. We reveal that rapid optical variations of ~100 s became more prominent when a thermal flare weakened and the optical spectrum flattened in the Ic, Rc, and V-band region. In conjunction with the flat spectrum in the radio range, this strongly indicates that the origin of rapid variations is not thermal emission, but synchrotron emission. Just after the outburst, we detected repeated flares at optical and X-ray wavelengths. The optical and X-ray light curves exhibited a strong correlation, with the X-rays, lagging by about 7 min. The X-ray lag can be understood in terms of a hot region propagating into the inner region of the accretion flow. The short X-ray lag, however, requires modifications of this simple scenario to account for the short propagation time. We also detected rapid optical variations with surprisingly high amplitude 50 days after the outburst, which we call optical flashes. During the most prominent optical flash, the object brightened by 1.2 mag only within 30 s. The released energy indicates that the emission source should be at the innermost region of the accretion flow.