We report photometric and spectroscopic observations and analysis of the 2019 superoutburst of TCP J21040470+4631129. This object showed a 9-mag superoutburst with early superhumps and ordinary superhumps, which are the features of WZ Sge-type dwarf novae. Five rebrightenings were observed after the main superoutburst. The spectra during the post-superoutburst stage showed the Balmer, He I and possible sodium doublet features. The mass ratio is derived as 0.0880(9) from the period of the superhump. During the third and fifth rebrightenings, growing superhumps and superoutbursts were observed, which have never been detected during a rebrightening phase among WZ Sge-type dwarf novae with multiple rebrightenings. To induce a superoutburst during the brightening phase, the accretion disk was needed to expand beyond the 3:1 resonance radius of the system again after the main superoutburst. These peculiar phenomena can be explained by the enhanced viscosity and large radius of the disk suggested by the higher luminosity and the presence of late-stage superhumps during the post-superoutburst stage, plus by more mass supply from the cool mass reservoir and/or from the secondary because of the enhanced mass transfer than those of other WZ Sge-type dwarf novae.

The star V426 Sge (HBHA~1704-05), originally classified as an emission-line object and a semi-regular variable, brightened at the beginning of August 2018, showing signatures of a symbiotic star outburst. We aim to confirm the nature of V426 Sge as a classical symbiotic star, determine the photometric ephemeris of the light minima, and suggest the path from its 1968 symbiotic nova outburst to the following 2018 Z And-type outburst. We re-constructed an historical light curve (LC) of V426 Sge from approximately the year 1900, and used original low- and high-resolution spectroscopy complemented with Swift-XRT and UVOT, optical UBVRcIc and near-infrared JHKL photometry obtained during the 2018 outburst and the following quiescence. The historical LC reveals no symbiotic-like activity from 1900 to 1967. In 1968, V426 Sge experienced a symbiotic nova outburst that ceased around 1990. From approximately 1972, a wave-like orbitally related variation with a period of $493.4\pm 0.7$ days developed in the LC. This was interrupted by a Z And-type outburst from the beginning of August 2018 to the middle of February 2019. At the maximum of the 2018 outburst, the burning white dwarf (WD) increased its temperature to $>2\times 10^5$ K, generated a luminosity of $\sim 7\times 10^{37}(d/3.3kpc)^2$ erg/s, and blew a wind at the rate of $\sim 3\times 10^{-6}$ M$_{\odot}$/yr. The donor is a normal M4-5 III giant and the accretor is a low-mass $\sim$0.5 M$_{\odot}$ WD. During the transition from the symbiotic nova outburst to the quiescent phase, a pronounced sinusoidal variation along the orbit develops in the LC of most symbiotic novae. The following eventual outburst is of Z And-type, when the accretion by the WD temporarily exceeds the upper limit of the stable burning. At this point the system becomes a classical symbiotic star.

Continuing the project described by Kato et al. (2009, PASJ, 61, S395, arXiv/0905.1757), we collected times of superhump maxima for 102 SU UMa-type dwarf novae observed mainly during the 2017 season and characterized these objects. WZ Sge-type stars identified in this study are PT And, ASASSN-17ei, ASASSN-17el, ASASSN-17es, ASASSN-17fn, ASASSN-17fz, ASASSN-17hw, ASASSN-17kd, ASASSN-17la, PNV J20205397+2508145 and TCP J00332502-3518565. We obtained new mass ratios for 7 objects using growing superhumps (stage A). ASASSN-17gf is an EI Psc-type object below the period minimum. CRTS J080941.3+171528 and DDE 51 are objects in the period gap and both showed long-lasting phase of stage A superhumps. We also summarized the recent advances in understanding of SU UMa-type and WZ Sge-type dwarf novae.

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.

Using 172 plates taken with the 40-cm astrograph of the Sternberg Astronomical Institute (Lomonosov Moscow University) in 1976-1994 and digitized with the resolution of 2400 dpi, we discovered and studied 275 new variable stars. We present the list of our new variables with all necessary information concerning their brightness variations. As in our earlier studies, the new discoveries show a rather large number of high-amplitude Delta Scuti variables, predicting that many stars of this type remain not detected in the whole sky. We also performed automated classification of the newly discovered variable stars based on the Random Forest algorithm. The results of the automated classification were compared to traditional classification and showed that automated classification was possible even with noisy photographic data. However, further improvement of automated techniques is needed, which is especially important having in mind the very large numbers of new discoveries expected from all-sky surveys.

Photographic plate archives contain a wealth of information about positions and brightness celestial objects had decades ago. Plate digitization is necessary to make this information accessible, but extracting it is a technical challenge. We develop algorithms used to extract photometry with the accuracy of better than ~0.1m in the magnitude range 13<B<17 from photographic images obtained in 1948-1996 with the 40cm Sternberg institute's astrograph (30x30cm plate size, 10x10deg field of view) and digitized using a flatbed scanner. The extracted photographic lightcurves are used to identify thousands of new high-amplitude (>0.2m) variable stars. The algorithms are implemented in the free software VaST available at http://scan.sai.msu.ru/vast/

Photometric measurements are prone to systematic errors presenting a challenge to low-amplitude variability detection. In search for a general-purpose variability detection technique able to recover a broad range of variability types including currently unknown ones, we test 18 statistical characteristics quantifying scatter and/or correlation between brightness measurements. We compare their performance in identifying variable objects in seven time series data sets obtained with telescopes ranging in size from a telephoto lens to 1m-class and probing variability on time-scales from minutes to decades. The test data sets together include lightcurves of 127539 objects, among them 1251 variable stars of various types and represent a range of observing conditions often found in ground-based variability surveys. The real data are complemented by simulations. We propose a combination of two indices that together recover a broad range of variability types from photometric data characterized by a wide variety of sampling patterns, photometric accuracies, and percentages of outlier measurements. The first index is the interquartile range (IQR) of magnitude measurements, sensitive to variability irrespective of a time-scale and resistant to outliers. It can be complemented by the ratio of the lightcurve variance to the mean square successive difference, 1/h, which is efficient in detecting variability on time-scales longer than the typical time interval between observations. Variable objects have larger 1/h and/or IQR values than non-variable objects of similar brightness. Another approach to variability detection is to combine many variability indices using principal component analysis. We present 124 previously unknown variable stars found in the test data.

Continuing the project described by Kato et al. (2009, arXiv:0905.1757), we collected times of superhump maxima for 128 SU UMa-type dwarf novae observed mainly during the 2015-2016 season and characterized these objects. The data have improved the distribution of orbital periods, the relation between the orbital period and the variation of superhumps, the relation between period variations and the rebrightening type in WZ Sge-type objects. Coupled with new measurements of mass ratios using growing stages of superhumps, we now have a clearer and statistically greatly improved evolutionary path near the terminal stage of evolution of cataclysmic variables. Three objects (V452 Cas, KK Tel, ASASSN-15cl) appear to have slowly growing superhumps, which is proposed to reflect the slow growth of the 3:1 resonance near the stability border. ASASSN-15sl, ASASSN-15ux, SDSS J074859.55+312512.6 and CRTS J200331.3-284941 are newly identified eclipsing SU UMa-type (or WZ Sge-type) dwarf novae. ASASSN-15cy has a short (~0.050 d) superhump period and appears to belong to EI Psc-type objects with compact secondaries having an evolved core. ASASSN-15gn, ASASSN-15hn, ASASSN-15kh and ASASSN-16bu are candidate period bouncers with superhump periods longer than 0.06 d. We have newly obtained superhump periods for 79 objects and 13 orbital periods, including periods from early superhumps. In order that the future observations will be more astrophysically beneficial and rewarding to observers, we propose guidelines how to organize observations of various superoutbursts.

We identify 339 known and 316 new variable stars of various types among 250000 lightcurves obtained by digitizing 167 30x30cm photographic plates of the Moscow collection. We use these data to conduct a comprehensive test of 18 statistical characteristics (variability indices) in search for the best general-purpose variability detection statistic. We find that the highest peak on the DFT periodogram, interquartile range, median absolute deviation, and Stetson's L index are the most efficient in recovering variable objects from the set of photographic lightcurves used in our test.

By digitizing astronomical photographic plates one may extract full information stored on them, something that could not be practically achieved with classical analogue methods. We are developing algorithms for variable objects search using digitized photographic images and apply them to 30cm (10x10 deg. field of view) plates obtained with the 40cm astrograph in 1940-90s and digitized with a flatbed scanner. Having more than 100 such plates per field, we conduct a census of high-amplitude (>0.3m) variable stars changing their brightness in the range 13<m<17 on timescales from hours to years in selected sky regions. This effort led to discovery of ~1000 new variable stars. We estimate that 1.2 +/- 0.1% of all stars show easily-detectable light variations; 0.7 +/- 0.1% of the stars are eclipsing binaries (64 +/- 4% of them are EA type, 22 +/- 2% are EW type and 14 +/- 2% are EB type); 0.3 +/- 0.1% of the stars are red variable giants and supergiants of M, SR and L types.

We present low-resolution phase-resolved spectra of the polar MT Dra during its high states in 2006 and 2009. Balmer series, He I and He II 4686 \AA\ emission lines have complex shapes and show similar profile variations over an orbital cycle. The radial velocities vary with orbital period and display motions of the gas falling close to the magnetic poles of the white dwarf.

The multi-site photometric observations of MN Dra were made over 77 nights in August-November, 2009. The total exposure was 433 hours. During this time the binary underwent two superoutbursts and five normal outbursts. During the course of first superoutburst period of positive superhumps decreased with extremely large $\dot P = -1.5 \times 1.0^{-4}$ for SU UMa-like dwarf novae, confirming known behavior of MN Dra [1]. Between the superoutbursts MN Dra displayed negative superhumps. Their period changed cyclically around 0.096-day value.

We present the results of photometric observations of SDSS J080434.20+510349.2 in its low state and during an outburst and spectroscopy during the outburst. We found such peculiarities as a long-term outburst with amplitude probably not less than 6m, rarity of the outbursts, eleven rebrightenings, and a short (0.059713(7) d) superhump period. We conclude that this star belongs to the WZ Sge-type subclass of cataclysmic variables. The spectrum shows both emission and absorption lines of H and He superimposed on the blue continuum. We also found 8 -- 9 min. brightness variations during the end of the super-outburst plateau that could be related to pulsations of the white dwarf in the binary.