We performed high-resolution optical spectroscopy and X-ray observations of the recently identified Mira-type symbiotic star EF Aql. Based on high-resolution optical spectroscopy obtained with SALT, we determine the temperature ($\sim $55 000 K) and the luminosity ($\sim$ 5.3 $L_\odot$) of the hot component in the system. The heliocentric radial velocities of the emission lines in the spectra reveal possible stratification of the chemical elements. We also estimate the mass-loss rate of the Mira donor star. Our Swift observation did not detect EF Aql in X-rays. The upper limit of the X-ray observations is 10$^{-12}$ erg cm$^{-2}$ s$^{-1}$, which means that EF Aql is consistent with the faintest X-ray systems detected so far. Otherwise we detected it with the UVOT instrument with an average UVM2 magnitude of 14.05. During the exposure, EF Aql became approximately 0.2 UVM2 magnitudes fainter. The periodogram analysis of the V-band data reveals an improved period of 320.4$\pm$0.3 d caused by the pulsations of the Mira-type donor star. The spectra are available upon request from the authors.

We analyse optical photometric data of short term variability (flickering) of the accreting white dwarf in the jet-ejecting symbiotic star MWC560. The observations are obtained in 17 nights during the period November 2011 - October 2019. The colour-magnitude diagram shows that the hot component of the system becomes redder as it gets brighter. For the flickering source we find that it has colour 0.14 < B-V < 0.40, temperature in the range 6300 < T_fl < 11000 K, and radius 1.2 < R_fl < 18 Rsun. We find a strong correlation (correlation coefficient 0.76, significance < 0.001) between B band magnitude and the average radius of the flickering source - as the brightness of the system increases the size of the flickering source also increases. The estimated temperature is similar to that of the bright spot of cataclysmic variables. In 2019 the flickering is missing, and the B-V colour of the hot component becomes bluer.

We study the photometric behavior of the recurrent nova RS Oph by 58 monitoring light curves (LCs), taken by 5 telescopes. All LCs show repeating time structures with some quasi-periods (QPs) in time scales from minutes to hours. In our previous work 97 QPs were detected in the LCs by local minimums of structure functions and local maximums of auto-correlation functions. The distribution of the QPs shows modes at 8, 13, 21, 30, 48 and 73 min, where the mode at 8 min is poorly unveiled. These modes follow a power function with base $1.55\approx 3/2$ with standard deviation 4.7%. This function predicts modes also at 5.3 and 3.5 min, which are not detected in the full MLCs. In the present work we analyze simple small parts from high resolution LCs. We confirm the QPs modes at 8.0, 5.3 and 3.5 min. Generally, we found 8 QP modes with regular logarithmic distribution in the time interval 3.5-73 min. We also show typical intra-night evolutions of QP modes in the minute scale -- sharp or gradual transitions from one QP mode to other. In the end we find that the parts of the LCs carry out the properties of the whole LCs at short time scale. This lead to two well pronounces dependences - between the range deviation and standard deviation of the LC, as well as between the quasi-period and the relevant level of the density function of the LC.

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 report photometric observations in Johnson UBV bands of the short term variability of Mira. The amplitude detected is 0.16 mag in B band. Adopting interstellar extinction E(B-V)=0, we find for the flickering source colour B-V=1.3, temperature T=3700 K, and radius R=0.77 Rsun. The colour of the flickering source is considerably redder than the average B-V colour of the cataclysmic variables.

We analyzed 29 pairs of time series in B and V bands of the recurrent nova RS Oph. The observations were carried out in 2008-2017 with duration 0.6 - 3.6 hours, with time resolution 0.5 - 3.3 min. We scanned digitally each series by data windows with various sizes Theta and derived two of the simplest fractal parameters for every Theta - standard deviation D and structural deviation S. Using the local minima of the structural function log S = f_S(log Theta) we unveiled 80 time structures, 42 in B band and 38 in V band, with time sizes 10-120 min. About 3/4 of the time sizes belong to the interval 10-40 min and about 1/4 lie in the interval 60-120 min. The respective cycles per day are 144-36 c/d and 24-15 c/d. On logarithmic scale, the distribution of the time sizes shows maximums at about 10, 21, 36 and 74 min. The 10 min flickering is poorly detectable in our series and we found the most widespread time structures (in about 1/5 of the cases) have time sizes about 21 min (about 69 c/d). Using the deviation function log D = f_D(log Theta) we estimated the relative cumulative energy (including the energy of the shorter structures in it), associated with the detected structure sizes, to be in the interval of the relative fluxes 2-11 %. The energies correlate weakly with the logarithms of the structure sizes, with correlation coefficients 0.60 and 0.57, under slope coefficients 0.04 and 0.03 in B and V band, respectively. The distributions of the energies occur bimodal, with maximums about 4% and 6% in B band, as well as about 3% and 5% in V band. The left and right modes of the distributions may be associated with the structure sizes 10 - 21 min and 37 - 74 min, respectively.

We performed photometry with a 1 minute time resolution of the symbiotic stars EF Aquilae, AG Pegasi and SU Lyncis in Johnson B and V band. Our observations of the symbiotic Mira-type star EF Aql demonstrate the presence of stochastic light variations with an amplitude of about 0.25 magnitudes on a time scale of 5 minutes. The observations prove the white dwarf nature of the hot component in the binary system. It is the 11th symbiotic star (among more than 200 symbiotic stars known in our Galaxy) which displays optical flickering. For SU Lyn we do not detect flickering with an amplitude above 0.03 mag in B band. For AG Peg, the amplitude of variability in B and V band is smaller than 0.05 mag and 0.04 mag respectively.

We performed 48.6 hours (in 28 nights) of simultaneous B and V band observations of the flickering variability of the recurrent nova RS Oph in quiescence. During the time of our observations the brightness of the system varied between 13.2 > B > 11.1 and the colour in the range 0.86 < B-V < 1.33 . We find that RS~Oph becomes more blue, as it becomes brighter, however the hot component becomes more red as it becomes brighter (assuming that the red giant is non-variable). During all the runs RS Oph exhibits flickering with amplitude 0.16 - 0.59 mag in B band. For the flickering source we find that it has colour -0.14 < B-V < 0.40, temperature in the range 7200 < T_fl < 18900, and average radius 1.1 < R_fl < 6.7 R_sun. We do not find a correlation between the temperature of the flickering and the brightness. However, we do find a strong correlation (correlation coefficient 0.81, significance 1.1x10^{-7} ) between B band magnitude and the average radius of the flickering source - as the brightness of the system increases the size of the flickering source also increases. The estimated temperature is similar to that of the bright spot of cataclysmic variables. The persistent presence of flickering indicates that the white dwarf is actively accreting material for the next outburst.

We obtained and analyzed more than 100 hours of multicolour optical time series of the blazar CTA 102 during its 2012 and 2016 outbursts. The object reached almost 11-th mag at the end of 2016, which is perhaps the brightest blazar state ever observed! During both outbursts, CTA 102 showed significant and rapid variability on intra-night time scales, reaching up to 0.2 mag for 30 min on some occasions. The "rms-flux" relation, built for all datasets, shows a large scatter and no apparent saturation on the magnitude scale. The ensemble structure function of the light curves can be fitted well with a straight line of a slope of $\sim$0.4. The time lags between the different optical bands appear to be consistent with zero, taking into account our time resolution. We discuss different variability scenarios and favor the changing Doppler factor of the emitting blobs as the most plausible one to account for the observed intra-night variability.

We report optical CCD photometry of the recently identified symbiotic star EF Aql. Our observations in Johnson V and B bands clearly show the presence of stochastic light variations with an amplitude of about 0.2 mag on a time scale of minutes. The observations point toward a white dwarf (WD) as the hot component in the system. It is the 11-th object among more than 200 symbiotic stars known with detected optical flickering. Estimates of the mass accretion rate onto the WD and the mass loss rate in the wind of the Mira secondary star lead to the conclusion that less than 1 per cent of the wind is captured by the WD. Eight further candidates for the detection of flickering in similar systems are suggested.

We report simultaneous observations of the flaring behaviour of the cataclysmic variable star AE Aqr. The observations are in Johnson B and V bands. The colour-magnitude diagrams (B-V versus V and B-V vs. B) show that the star becomes blues as it becomes brighter. In our model AE Aqr behaviour can be explained with flares (fireballs) with 0.03 < B-V < 0.30 and temperature in the interval 8000 < T < 12000.

We analyze optical photometric data of short term variability (flickering) of accreting white dwarfs in cataclysmic variables (KR Aur, MV Lyr, V794 Aql, TT Ari, V425 Cas), recurrent novae (RS Oph and T CrB) and jet-ejecting symbiotic stars (CH Cyg and MWC 560). We find that the amplitude-flux relationship is visible over four orders of magnitude, in the range of fluxes from $10^{29}$ to $10^{33}$ erg s$^{-1}$ \AA$^{-1}$, as a "statistically perfect" correlation with correlation coefficient 0.96 and p-value $ \sim 10^{-28}$. In the above range, the amplitude of variability for any of our 9 objects is proportional to the flux level with (almost) one and the same factor of proportionality for all 9 accreting white dwarfs with $\Delta F = 0.36 (\pm 0.05) F_{av}$, $\sigma_{rms} = 0.086(\pm 0.011) F_{av}$, and $\sigma_{rms} / \Delta F = 0.24 \pm 0.02$. Over all, our results indicate that the viscosity in the accretion discs is practically the same for all 9 objects in our sample, in the mass accretion rate range $2 \times 10^{-11} - 2\times10^{-7}$ $M_\odot$ yr$^{-1}$.

We analyzed the multi-band optical behaviour of the BL Lacertae object, S5 0716+714, during its outburst state from 2014 November - 2015 March. We took data on 23 nights at three observatories, one in India and two in Bulgaria, making quasi-simultaneous observations in B, V, R, and I bands. We measured multi-band optical fluxes, colour and spectral variations for this blazar on intraday and short timescales. The source was in a flaring state during the period analyzed and displayed intense variability in all wavelengths. R band magnitude of 11.6 was attained by the target on 18 Jan 2015, which is the brightest value ever recorded for S5 0716+714. The discrete correlation function method yielded good correlation between the bands with no measurable time lags, implying that radiation in these bands originate from the same region and by the same mechanism. We also used the structure function technique to look for characteristic timescales in the light curves. During the times of rapid variability, no evidence for the source to display spectral changes with magnitude was found on either of the timescales. The amplitude of variations tends to increase with increasing frequency with a maximum of $\sim$ 22% seen during flaring states in B band. A mild trend of larger variability amplitude as the source brightens was also found. We found the duty cycle of our source during the analyzed period to be $\sim$ 90%. We also investigated the optical spectral energy distribution of S5 0716+714 using B, V, R, and I data points for 21 nights. We briefly discuss physical mechanisms most likely responsible for its flux and spectral variations.

To search for optical variability on a wide range of timescales, we have carried out photometric monitoring of two flat spectrum radio quasars, 3C 454.3 and 3C 279, plus one BL Lac, S5 0716+714, all of which have been exhibiting remarkably high activity and pronounced variability at all wavelengths. CCD magnitudes in B, V, R and I pass-bands were determined for $\sim$ 7000 new optical observations from 114 nights made during 2011 - 2014, with an average length of $\sim$ 4 h each, at seven optical telescopes: four in Bulgaria, one in Greece, and two in India. We measured multiband optical flux and colour variations on diverse timescales. Discrete correlation functions were computed among B, V, R, and I observations, to search for any time delays. We found weak correlations in some cases with no significant time lags. The structure function method was used to estimate any characteristic time-scales of variability. We also investigated the spectral energy distribution of the three blazars using B, V, R, I, J and K pass-band data. We found that the sources almost always follows a bluer-when-brighter trend. We discuss possible physical causes of the observed spectral variability.

We report observations of the flickering variability of the symbiotic recurrent nova RS~Oph at quiescence in five bands ($UBVRI$). We find evidence of a correlation between the peak-to-peak flickering amplitude ($\Delta F$) and the average flux of the hot component ($F_{\rm av}$). The correlation is highly significant, with a correlation coefficient of 0.85 and a $p$-value of~$\sim 10^{-20}$. Combining the data from all wavebands, we find a dependence of the type $\Delta F \propto F^k_{\rm av}$, with power-law index $k = 1.02 \pm 0.04$ for the $UBVRI$ flickering of RS~Oph. Thus, the relationship between the amplitude of variability and the average flux of the hot component is consistent with linearity. The rms amplitude of flickering is on average 8 per cent ($\pm2$ per cent) of $F_{\rm av}$. The detected correlation is similar to that found in accreting black holes/neutron stars and cataclysmic variables. The possible reasons are briefly discussed. The data are available upon request from the authors.

We performed simultaneous observations in 3 bands (UBV) of the flickering variability of the recurrent novae RS Oph and T CrB at quiescence. Using new and published data, we compare the colours of the flickering in cataclysmic variables and symbiotic recurrent novae. We find a difference between the colours of the flickering source in these two types of accreting white dwarfs. The detected difference is highly significant with $ p-value \approx 2 \times 10^{-6}$ for the distributions of $(U-B)_0$ colour and $p \approx 3 \times 10^{-5}$ on (U-B) versus (B-V) diagram. The possible physical reasons are briefly discussed. The data are available upon request from the authors.

By using standard broad-band VRI photometry we were able to discriminate the variations of the broad hydrogen alpha line from the continuum variations for the active galaxy Mkn 279. Cross-correlating both light curves enabled us to determine the time lag of the broad line variations behind the continuum and thus to determine the BLR size (about 8 light days). Our preliminary results are rather consistent with the spectroscopic reverberation mapping results (about 12/17 days). This study is a part of an ambitious program to perform photometric reverberation mapping and determine BLR sizes (respectively - the central black hole masses) for more that 100 nearby AGN.

We report simultaneous multicolour observations in 5 bands (UBVRI) of the flickering variability of the cataclysmic variable AE Aqr. Our aim is to estimate the parameters (colours, temperature, size) of the fireballs that produce the optical flares. The observed rise time of the optical flares is in the interval 220 - 440 sec. We estimate the dereddened colours of the fireballs: (U-B)_0 in the range 0.8-1.4, (B-V)_0 ~ 0.03-0.24, (V-I)_0 ~ 0.26-0.78. We find for the fireballs a temperature in the range 10000 - 25000 K, mass (7-90).10^{19} g, size (3-7).10^9 cm (using a distance of d=86 pc). These values refer to the peak of the flares observed in UBVRI bands. The data are available upon request from the authors.

We report the results of optical monitoring for a sample of 11 blazars including 10 BL Lacs and 1 Flat Spectrum Radio Quasar (FSRQ). We have measured the multiband optical flux and colour variations in these blazars on intra-day and short-term timescales of months and have limited data for 2 more blazars. These photometric observations were made during 2009 to 2011, using six optical telescopes, four in Bulgaria, one in Greece and one in India. On short-term timescales we found significant flux variations in 9 of the sources and colour variations in 3 of them. Intra-day variability was detected on 6 nights for 2 sources out of the 18 nights and 4 sources for which we collected such data. These new optical observations of these blazars plus data from our previous published papers (for 3 more blazars) were used to analyze their spectral flux distributions in the optical frequency range. Our full sample for this purpose includes 6 high-synchrotron-frequency-peaked BL Lacs (HSPs), 3 intermediate-synchrotron-frequency-peaked BL Lacs (ISPs) and 6 low-synchrotron-frequency-peaked BL Lacs (LSPs; including both BL Lacs and FSRQs). We also investigated the spectral slope variability and found that the average spectral slopes of LSPs show a good accordance with the Synchrotron Self-Compton (SSC) loss dominated model. Our analysis supports previous studies that found that the spectra of the HSPs and FSRQs have significant additional emission components. The spectra of all these HSPs and LSPs get flatter when they become brighter, while for FSRQs the opposite appears to hold. This supports the hypothesis that there is a significant thermal contribution to the optical spectrum for FSRQs.

In this paper we present results of a short-term optical monitoring of 13 blazars. The objects were monitored mostly in the R-band for a total of ~ 160 hours between 2006 and 2011. We study the nature of the short-term variations and show that most of them could be described as slow, smooth, and (almost) linear changes of up to ~ 0.1 mag/hour, but many objects show no short-term variations at all. In fact, we found only ~ 2 per cent chance to observe variability of more than 0.1 mag/hour for the sample we observed. Hints for quasi-periodic oscillations at very low amplitude levels are also found for some objects. We briefly discuss some of the possible mechanisms to generate the intra-night variability and the quasi-periodic oscillations.

We report the results of quasi-simultaneous two filter optical monitoring of two high-energy peaked blazars, 1ES 1959+650 and 1ES 2344+514, to search for microvariability and short-term variability (STV). We carried out optical photometric monitoring of these sources in an alternating sequence of B and R pass-bands, and have 24 and 19 nights of new data for these two sources, respectively. No genuine microvariability (intra-night variability) was detected in either of these sources. This non-detection of intra-night variations is in agreement with the conclusions of previous studies that high-energy peaked BL Lacs are intrinsically less variable than low-energy peaked BL Lacs in the optical bands. We also report the results of STV studies for these two sources between July 2009 and August 2010. Genuine STV is found for the source 1ES 1959+650 but not for 1ES 2344+514. We briefly discuss possible reasons for the difference between the intra-night variability behaviour of high- and low-energy peaked blazars.

We present the results of our study of spectral energy distributions (SEDs) of a sample of ten low- to intermediate-synchrotron-peaked blazars. We investigate some of the physical parameters most likely responsible for the observed short-term variations in blazars. To do so, we focus on the study of changes in the SEDs of blazars corresponding to changes in their respective optical fluxes. We model the observed spectra of blazars from radio to optical frequencies using a synchrotron model that entails a log-parabolic distribution of electron energies. A significant correlation among the two fitted spectral parameters ($a$, $b$) of log-parabolic curves and a negative trend among the peak frequency and spectral curvature parameter, $b$, emphasize that the SEDs of blazars are fitted well by log-parabolic curves. On considering each model parameter that could be responsible for changes in the observed SEDs of these blazars, we find that changes in the jet Doppler factors are most important.

We have measured multi-band optical flux and colour variations for a sample of 12 low energy peaked blazars on short, day-to-month, timescales. Our sample contains six BL Lacertae objects and six flat spectrum radio quasars. These photometric observations, made during September 2008 to June 2009, used five optical telescopes, one in India and four in Bulgaria. We detected short term flux variations in eleven of these blazars and colour variability in eight of them. Our data indicate that six blazars (3C 66A, AO 0235+164, S5 0716+714, PKS 0735+178, OJ 287 and 3C 454.3) were observed in pre- or post-outburst states, that five (PKS 0420-014, 4C 29.45, 3C 279, PKS 1510-089 and BL Lac) were in a low state, while one (3C 273) was in an essentially steady state. The duty cycles for flux and colour variations on short timescales in these low energy peaked blazars are ~ 92 percent and ~ 33 percent, respectively. The colour vs magnitude correlations seen here support the hypothesis that BL Lac objects tend to become bluer with increase in brightness; however, flat spectrum radio quasars may show the opposite trend, and there are exceptions to these trends in both categories of blazar. We briefly discuss emission models for active galactic nuclei that might explain our results.

We report observations of the flickering variability of the recurrent nova RS Oph at quiescence on the basis of simultaneous observations in 5 bands (UBVRI). RS Oph has flickering source with (U-B)_0=-0.62 \pm 0.07, (B-V)_0=0.15 \pm 0.10, (V-R)_0=0.25 \pm 0.05. We find for the flickering source a temperature T_fl = 9500 \pm 500 K, and luminosity L_fl = 50 - 150 L_sun (using a distance of d=1.6kpc). We also find that on a (U-B) vs (B-V) diagram the flickering of the symbiotic stars differs from that of the cataclysmic variables. The possible source of the flickering is discussed. The data are available upon request from the authors and on the web www.astro.bas.bg/~rz/RSOph.UBVRI.2010.MNRAS.tar.gz.

We present the results from a monitoring campaign of the Narrow-Line Seyfert~1 galaxy PG 1211+143. The object was monitored with ground-based facilities (UBVRI photometry; from February to July, 2007) and with Swift (X-ray photometry/spectroscopy and UV/Optical photometry; between March and May, 2007). We found PG 1211+143 in a historical low X-ray flux state at the beginning of the Swift monitoring campaign in March 2007. It is seen from the light curves that while violently variable in X-rays, the quasar shows little variations in optical/UV bands. The X-ray spectrum in the low state is similar to other Narrow-Line Seyfert 1 galaxies during their low-states and can be explained by a strong partial covering absorber or by X-ray reflection onto the disk. With the current data set, however, it is not possible to distinguish between both scenarios. The interband cross-correlation functions indicate a possible reprocessing of the X-rays into the longer wavelengths, consistent with the idea of a thin accretion disk, powering the quasar. The time lags between the X-ray and the optical/UV light curves, ranging from ~2 to ~18 days for the different wavebands, scale approximately as ~lambda^(4/3), but appear to be somewhat larger than expected for this object, taking into account its accretion disk parameters. Possible implications for the location of the X-ray irradiating source are discussed.