ASASSN-14dx showed an extraordinary outburst whose features are the small outburst amplitude (~ 2.3 mag) and long duration (> 4 years). Because we found a long observational gap of 123 d before the outburst detection, we propose that the main outburst plateau was missed and that this outburst is just a "fading tail" often seen after the WZ Sge-type superoutbursts. In order to distinguish between WZ Sge and SU UMa-type dwarf novae (DNe), we investigated Gaia DR2 statistically. We applied a logistic regression model and succeeded in classifying by using absolute Gaia magnitudes $M_{G}$ and Gaia colors $G_{\rm BP}-G_{\rm RP}$. Our new classifier also suggests that ASASSN-14dx is the best candidate of a WZ Sge-type DN. We estimated distances from the earth of known WZ Sge stars by using Gaia DR2 parallaxes. The result indicates that ASASSN-14dx is the third nearest WZ Sge star (next to WZ Sge and V455 And), and hence the object can show the third brightest WZ Sge-type superoutburst whose maximum is $V$ = 8-9 mag.

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.

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.

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

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

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

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

We report on photometric observations of two dwarf novae, OT J075418.7+381225 and OT J230425.8+062546, which underwent superoutburst in 2013 (OT J075418) and in 2011 (OT J230425). Their mean period of the superhump was 0.0722403(26) d (OT J074518) and 0.067317(35) d (OT J230425). These objects showed a very long growing stage of the superhump (stage A) and a large period decrease in stage A-B transition. The long stage A suggests slow evolution of the superhump due to very small mass ratios of these objects. The decline rates during the plateau phase in the superoutburst of these objects were lower than those of SU UMa-type DNe with a similar superhump period. These properties were similar to those of SSS J122221.7-311523, the best candidate for the period bouncer. Therefore, these two DNe are regarded as good candidates for the period bouncer. We estimated the number density of period bouncers roughly from our observations in the recent five years. There is a possibility that these WZ Sge-type dwarf novae with unusual outburst properties can account for the missing population of the period bouncer expected from the evolutionary scenario.

The analysis of 14 periodograms of EZ Lyn for the data spaced over 565 d in 2012--2014 (2-3.5 yr after 2010 outburst) yielded the existence of the stable signals around 100 c/d and three signals around 310 c/d, 338 c/d and 368 c/d (the corresponding periods are 864 s, 279 s, 256 s and 235 s). We interpret them as independent non-radial pulsations of the white dwarf in EZ Lyn, but a possibility that a linear combination of frequency at 100 c/d and harmonic of orbital period could produce the frequency at 368 c/d also cannot be excluded. The signal at 100 c/d was detected during the first stay in the instability strip as a transient one. The period at 338 c/d, is a known non-radial pulsation EZ Lyn entered the instability strip after the 2010 outburst. We detected the signals around 310 c/d and 368 c/d for the first time. We applied the two-dimensional least absolute shrinkage and selection operator (Lasso) analysis for the first time to explore the behavior of these signals on the scale of hours for nightly runs of observations having duration of 6-12 hr. The Lasso analysis revealed the simultaneous existence of all three frequencies (310 c/d, 338 c/d and 368 c/d) for majority of nights of observations, but with variable amplitudes and variable drifts of frequencies by 2-6 percents on a time scale of ~5-7 hr. The largest drift we detected corresponded to 17.5 s in period in ~5 hours.

We present photometric study of NY Ser, an in-the-gap SU UMa-type nova, in 2002 and 2013. We determined the duration of the superoutburst and the mean superhump period to be 18 d and 0.10458 d, respectively. We detected in 2013 that NY Ser showed two distinct states separated by the superoutburst. A state of rather infrequent normal outbursts lasted at least 44 d before the superoutburst and a state of frequent outbursts started immediately after the superoutburst and lasted at least for 34 d. Unlike a typical SU UMa star with bimodal distribution of the outbursts duration, NY Ser displayed a diversity of normal outbursts. In the state of infrequent outbursts, we detected a wide ~12 d outburst accompanied by 0.098 d orbital modulation but without superhumps ever established in NY Ser. We classified this as the "wide normal outburst". The orbital period dominated both in quiescence and during normal outbursts in this state. In the state of the most frequent normal outbursts, the 0.10465 d positive superhumps dominated and co-existed with the orbital modulation. In 2002 we detected the normal outburst of "intermediate" 5-6 d duration that was also accompanied by orbital modulations.

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.

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

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

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

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

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

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.

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.

We propose a new method to reconstruct the structure of accretion disks in dwarf novae using multi-band light curves of early superhumps. Our model assumes that early superhumps are caused by the rotation effect of non-axisymmetrically flaring disks. We have developed a Bayesian model for this reconstruction, in which a smoother disk-structure tends to have a higher prior probability. We analyzed simultaneous optical and near-infrared photometric data of early superhumps of the dwarf nova, V455 And using this technique. The reconstructed disk has two flaring parts in the outermost region of the disk. These parts are responsible for the primary and secondary maxima of the light curves. The height-to-radius ratio is h/r=0.20-0.25 in the outermost region. In addition to the outermost flaring structures, flaring arm-like patterns can be seen in an inner region of the reconstructed disk. The overall profile of the reconstructed disk is reminiscent of the disk structure that is deformed by the tidal effect. However, an inner arm-like pattern, which is responsible for the secondary minimum in the light curve, cannot be reproduced only by the tidal effect. It implies the presence of another mechanism that deforms the disk structure. Alternatively, the temperature distribution of the disk could be non-axisymmetric. We demonstrate that the disk structure with weaker arm-like patterns is optimal in the model including the irradiation effect. However, the strongly irradiated disk gives quite blue colors, which may conflict with the observation. Our results suggest that the amplitude of early superhumps depends mainly on the height of the outermost flaring regions of the disk. We predict that early superhumps can be detected with an amplitude of >0.02 mag in about 90% of WZ Sge stars.

We report on a discovery of "negative" superhumps during the 2011 January superoutburst of ER UMa. During the superoutburst which started on 2011 January 16, we detected negative superhumps having a period of 0.062242(9) d, shorter than the orbital period by 2.2%. No evidence of positive superhumps was detected during this observation. This finding indicates that the disk exhibited retrograde precession during this superoutburst, contrary to all other known cases of superoutbursts. The duration of this superoutburst was shorter than those of ordinary superoutbursts and the intervals of normal outbursts were longer than ordinary ones. We suggest a possibility that such unusual outburst properties are likely a result of the disk tilt, which is supposed to be a cause of negative superhumps: the tilted disk could prevent the disk from being filled with materials in the outmost region which is supposed to be responsible for long-duration superoutbursts in ER UMa-type dwarf novae. The discovery signifies the importance of the classical prograde precession in sustaining long-duration superoutbursts. Furthermore, the presence of pronounced negative superhumps in this system with a high mass-transfer rate favors the hypothesis that hydrodynamical lift is the cause of the disk tilt.

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.

We systematically surveyed period variations of superhumps in SU UMa-type dwarf novae based on newly obtained data and past publications. In many systems, the evolution of superhump period are found to be composed of three distinct stages: early evolutionary stage with a longer superhump period, middle stage with systematically varying periods, final stage with a shorter, stable superhump period. During the middle stage, many systems with superhump periods less than 0.08 d show positive period derivatives. Contrary to the earlier claim, we found no clear evidence for variation of period derivatives between superoutburst of the same object. We present an interpretation that the lengthening of the superhump period is a result of outward propagation of the eccentricity wave and is limited by the radius near the tidal truncation. We interpret that late stage superhumps are rejuvenized excitation of 3:1 resonance when the superhumps in the outer disk is effectively quenched. Many of WZ Sge-type dwarf novae showed long-enduring superhumps during the post-superoutburst stage having periods longer than those during the main superoutburst. The period derivatives in WZ Sge-type dwarf novae are found to be strongly correlated with the fractional superhump excess, or consequently, mass ratio. WZ Sge-type dwarf novae with a long-lasting rebrightening or with multiple rebrightenings tend to have smaller period derivatives and are excellent candidate for the systems around or after the period minimum of evolution of cataclysmic variables (abridged).