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.

We report on a multi-band variability and correlation study of the TeV blazar Mrk 421 during an exceptional flaring activity observed from 2013 April 11 to 2013 April 19. The study uses, among others, data from GASP-WEBT, Swift, NuSTAR, Fermi-LAT, VERITAS, and MAGIC. The large blazar activity, and the 43 hours of simultaneous NuSTAR and MAGIC/VERITAS observations, permitted variability studies on 15 minute time bins, and over three X-ray bands (3-7 keV, 7-30 keV and 30-80 keV) and three very-high-energy (>0.1 TeV, hereafter VHE) gamma-ray bands (0.2-0.4 TeV, 0.4-0.8 TeV and >0.8 TeV). We detected substantial flux variations on multi-hour and sub-hour timescales in all the X-ray and VHE gamma-ray bands. The characteristics of the sub-hour flux variations are essentially energy-independent, while the multi-hour flux variations can have a strong dependence on the energy of the X-ray and the VHE gamma rays. The three VHE bands and the three X-ray bands are positively correlated with no time-lag, but the strength and the characteristics of the correlation changes substantially over time and across energy bands. Our findings favour multi-zone scenarios for explaining the achromatic/chromatic variability of the fast/slow components of the light curves, as well as the changes in the flux-flux correlation on day-long timescales. We interpret these results within a magnetic reconnection scenario, where the multi-hour flux variations are dominated by the combined emission from various plasmoids of different sizes and velocities, while the sub-hour flux variations are dominated by the emission from a single small plasmoid moving across the magnetic reconnection layer.

We report the results of our optical (VRI) photometric observations of the TeV blazar 1ES 0806$+$524 on 153 nights during 2011-2019 using seven optical telescopes in Europe and Asia. We investigated the variability of the blazar on intraday as well as on long-term timescales. We examined eighteen intraday light curves for flux and color variations using the most reliable power-enhanced F-test and the nested ANOVA test. Only on one night was a small, but significant, variation found, in both $V$ band and $R$ band light curves. The $V-R$ color index was constant on every one of those nights. Flux density changes of around 80 % were seen over the course of these eight years in multiple bands. We found a weighted mean optical spectral index of 0.639$\pm$0.002 during our monitoring period by fitting a power law ($F_{\nu} \propto \nu^{-\alpha}$) in 23 optical ($VRI$) spectral energy distributions of 1ES 0806$+$524. We discuss different possible mechanisms responsible for blazar variability on diverse timescales.

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.

Here we present quasi-simultaneous observations of the flickering of the symbiotic binary star CH Cyg in U, B and V bands. We calculate the flickering source parameters and discuss the possible reason for the flickering cessation in the period 2010-2013.

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 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 present results of optical observations of the extremely violently variable blazar S4 0954+65 on intra-night time scales. The object showed flux changes of up to 100% within a few hours. Time delays between optical bands, color changes and "rms-flux" relations are investigated and the results are discussed in terms of existing models of blazar variability.

We report optical spectroscopic observations of the Be/gamma-ray binaries LSI+61303, MWC 148 and MWC 656. The peak separation and equivalent widths of prominent emission lines (H-alpha, H-beta, H-gamma, HeI, and FeII) are measured. We estimated the circumstellar disc size, compared it with separation between the components, and discussed the disc truncation. We find that in LSI+61303 the compact object comes into contact with the outer parts of the circumstellar disc at periastron, in MWC 148 the compact object goes deeply into the disc during the periastron passage, and in MWC 656 the black hole is accreting from the outer parts of the circumstellar disc along the entire orbit. The interstellar extinction was estimated using interstellar lines. The rotation of the mass donors appears to be similar to the rotation of the mass donors in Be/X-ray binaries. We suggest that X-ray/optical periodicity of about 1 day deserves to be searched for.

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 present coordinated multiwavelength observations of the bright, nearby BL Lac object Mrk 421 taken in 2013 January-March, involving GASP-WEBT, Swift, NuSTAR, Fermi-LAT, MAGIC, VERITAS, and other collaborations and instruments, providing data from radio to very-high-energy (VHE) gamma-ray bands. NuSTAR yielded previously unattainable sensitivity in the 3-79 keV range, revealing that the spectrum softens when the source is dimmer until the X-ray spectral shape saturates into a steep power law with a photon index of approximately 3, with no evidence for an exponential cutoff or additional hard components up to about 80 keV. For the first time, we observed both the synchrotron and the inverse-Compton peaks of the spectral energy distribution (SED) simultaneously shifted to frequencies below the typical quiescent state by an order of magnitude. The fractional variability as a function of photon energy shows a double-bump structure which relates to the two bumps of the broadband SED. In each bump, the variability increases with energy which, in the framework of the synchrotron self-Compton model, implies that the electrons with higher energies are more variable. The measured multi-band variability, the significant X-ray-to-VHE correlation down to some of the lowest fluxes ever observed in both bands, the lack of correlation between optical/UV and X-ray flux, the low degree of polarization and its significant (random) variations, the short estimated electron cooling time, and the significantly longer variability timescale observed in the NuSTAR light curves point toward in-situ electron acceleration, and suggest that there are multiple compact regions contributing to the broadband emission of Mrk 421 during low-activity states.

We report on simultaneous broadband observations of the TeV-emitting blazar Markarian 501 between 1 April and 10 August 2013, including the first detailed characterization of the synchrotron peak with Swift and NuSTAR. During the campaign, the nearby BL Lac object was observed in both a quiescent and an elevated state. The broadband campaign includes observations with NuSTAR, MAGIC, VERITAS, the Fermi Large Area Telescope (LAT), Swift X-ray Telescope and UV Optical Telescope, various ground-based optical instruments, including the GASP-WEBT program, as well as radio observations by OVRO, Mets\"ahovi and the F-Gamma consortium. Some of the MAGIC observations were affected by a sand layer from the Saharan desert, and had to be corrected using event-by-event corrections derived with a LIDAR (LIght Detection And Ranging) facility. This is the first time that LIDAR information is used to produce a physics result with Cherenkov Telescope data taken during adverse atmospheric conditions, and hence sets a precedent for the current and future ground-based gamma-ray instruments. The NuSTAR instrument provides unprecedented sensitivity in hard X-rays, showing the source to display a spectral energy distribution between 3 and 79 keV consistent with a log-parabolic spectrum and hard X-ray variability on hour timescales. None (of the four extended NuSTAR observations) shows evidence of the onset of inverse-Compton emission at hard X-ray energies. We apply a single-zone equilibrium synchrotron self-Compton model to five simultaneous broadband spectral energy distributions. We find that the synchrotron self-Compton model can reproduce the observed broadband states through a decrease in the magnetic field strength coinciding with an increase in the luminosity and hardness of the relativistic leptons responsible for the high-energy emission.

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.

We present new radial velocities of the high-mass X-ray binary star 4U 2206+54 based on optical spectra obtained with the Coude spectrograph at the 2m RCC telescope at the Rozhen National Astronomical Observatory, Bulgaria in the period November 2011 -- July 2013. The radial velocity curve of the HeI $\lambda$6678 \AA\ line is modeled with an orbital period P$_{orb}$ = 9.568~d and an eccentricity of $e$ = 0.3. These new measurements of the radial velocity resolve the disagreements of the orbital period discussions.

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 138 spectral observations of the H-alpha emission line of the radio- and gamma-ray-emitting Be/X-ray binary LSI+61303 obtained during the period of September 1998 -- January 2013. From measuring various H-alpha parameters, we found that the orbital modulation of the H-alpha is best visible in the equivalent width ratio EW(B)/EW(R), the equivalent width of the blue hump, and in the radial velocity of the central dip. The periodogram analysis confirmed that the H-alpha emission is modulated with the orbital and superorbital periods. For the past 20 years the radius of the circumstellar disk is similar to the Roche lobe size at the periastron. It is probably truncated by a 6:1 resonance. The orbital maximum of the equivalent width of H-alpha emission peaks after the periastron and coincides on average with the X-ray and gamma-ray maxima. All the spectra are available upon request from the authors and through the CDS.

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.