Context. It has become evident that one-zone synchrotron self-Compton models are not always adequate for very-high-energy (VHE) gamma-ray emitting blazars. While two-component models are performing better, they are difficult to constrain due to the large number of free parameters. Aims. In this work, we make a first attempt to take into account the observational constraints from Very Long Baseline Interferometry (VLBI) data, long-term light curves (radio, optical, and X-rays) and optical polarisation to limit the parameter space for a two-component model and test if it can still reproduce the observed spectral energy distribution (SED) of the blazars. Methods. We selected five TeV BL Lac objects based on the availability of VHE gamma-ray and optical polarisation data. We collected constraints for the jet parameters from VLBI observations. We evaluated the contributions of the two components to the optical flux by means of decomposition of long-term radio and optical light curves as well as modeling of the optical polarisation variability of the objects. We selected eight epochs for these five objects, based on the variability observed at VHE gamma rays, for which we constructed the SEDs that we then modeled with a two-component model. Results. We found parameter sets which can reproduce the broadband SED of the sources in the framework of two-component models considering all available observational constraints from VLBI observations. Moreover, the constraints obtained from the long-term behavior of the sources in the lower energy bands could be used to determine the region where the emission in each band originates. Finally, we attempted to use optical polarisation data to shed new light on the behavior of the two components in the optical band. Our observationally constrained two zone model allows explanation of the entire SED from radio to VHE with two co-located emission regions.

We report on the results of a multiwavelength monitoring campaign of the bright, nearby Seyfert galaxy, Ark 120 using a ~50-day observing programme with Swift and a ~4-month co-ordinated ground-based observing campaign, predominantly using the Skynet Robotic Telescope Network. We find Ark 120 to be variable at all optical, UV, and X-ray wavelengths, with the variability observed to be well-correlated between wavelength bands on short timescales. We perform cross-correlation analysis across all available wavelength bands, detecting time delays between emission in the X-ray band and the Swift V, B and UVW1 bands. In each case, we find that the longer-wavelength emission is delayed with respect to the shorter-wavelength emission. Within our measurement uncertainties, the time delays are consistent with the \tau ~ \lambda^{4/3} relation, as predicted by a disc reprocessing scenario. The measured lag centroids are \tau_{cent} = 11.90 +/- 7.33, 10.80 +/- 4.08, and 10.60 +/- 2.87 days between the X-ray and V, B, and UVW1 bands, respectively. These time delays are longer than those expected from standard accretion theory and, as such, Ark 120 may be another example of an active galaxy whose accretion disc appears to exist on a larger scale than predicted by the standard thin-disc model. Additionally, we detect further inter-band time delays: most notably between the ground-based I and B bands (\tau_{cent} = 3.46 +/- 0.86 days), and between both the Swift XRT and UVW1 bands and the I band (\tau_{cent} = 12.34 +/- 4.83 and 2.69 +/- 2.05 days, respectively), highlighting the importance of co-ordinated ground-based optical observations.

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 decade-long (2008-2018) $\gamma$-ray to 1 GHz radio monitoring of the blazar 3C 279, including GASP/WEBT, $\it{Fermi}$ and $\it{Swift}$ data, as well as polarimetric and spectroscopic data. The X-ray and $\gamma$-ray light curves correlate well, with no delay > 3 hours, implying general co-spatiality of the emission regions. The $\gamma$-ray-optical flux-flux relation changes with activity state, ranging from a linear to a more complex dependence. The behaviour of the Stokes parameters at optical and radio wavelengths, including 43 GHz VLBA images, supports either a predominantly helical magnetic field or motion of the radiating plasma along a spiral path. Apparent speeds of emission knots range from 10 to 37c, with the highest values requiring bulk Lorentz factors close to those needed to explain $\gamma$-ray variability on very short time scales. The Mg II emission line flux in the `blue' and `red' wings correlates with the optical synchrotron continuum flux density, possibly providing a variable source of seed photons for inverse Compton scattering. In the radio bands we find progressive delays of the most prominent light curve maxima with decreasing frequency, as expected from the frequency dependence of the $\tau=1$ surface of synchrotron self-absorption. The global maximum in the 86 GHz light curve becomes less prominent at lower frequencies, while a local maximum, appearing in 2014, strengthens toward decreasing frequencies, becoming pronounced at $\sim5$ GHz. These tendencies suggest different Doppler boosting of stratified radio-emitting zones in the jet.

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.

We present a multiwavelength study of the flat-spectrum radio quasar CTA 102 during 2013-2017. We use radio-to-optical data obtained by the Whole Earth Blazar Telescope, 15 GHz data from the Owens Valley Radio Observatory, 91 and 103 GHz data from the Atacama Large Millimeter Array, near-infrared data from the Rapid Eye Monitor telescope, as well as data from the Swift (optical-UV and X-rays) and Fermi ($\gamma$ rays) satellites to study flux and spectral variability and the correlation between flux changes at different wavelengths. Unprecedented $\gamma$-ray flaring activity was observed during 2016 November-2017 February, with four major outbursts. A peak flux of (2158 $\pm$ 63)$\times$10$^{-8}$ ph cm$^{-2}$ s$^{-1}$, corresponding to a luminosity of (2.2 $\pm$ 0.1)$\times$10$^{50}$ erg s$^{-1}$, was reached on 2016 December 28. These four $\gamma$-ray outbursts have corresponding events in the near-infrared, optical, and UV bands, with the peaks observed at the same time. A general agreement between X-ray and $\gamma$-ray activity is found. The $\gamma$-ray flux variations show a general, strong correlation with the optical ones with no time lag between the two bands and a comparable variability amplitude. This $\gamma$-ray/optical relationship is in agreement with the geometrical model that has successfully explained the low-energy flux and spectral behaviour, suggesting that the long-term flux variations are mainly due to changes in the Doppler factor produced by variations of the viewing angle of the emitting regions. The difference in behaviour between radio and higher energy emission would be ascribed to different viewing angles of the jet regions producing their emission.

We present the results of photometric observations of three TeV blazars, 3C 66A, S5 0954+658 and BL Lacertae, during the period 2013--2017. Our extensive observations were performed in a total of 360 nights which produced $\sim$6820 image frames in BVRI bands. We study flux and spectral variability of these blazars on these lengthy timescales. We also examine the optical Spectral Energy Distributions of these blazars, which are crucial in understanding the emission mechanism of long-term variability in blazars. All three TeV blazars exhibited strong flux variability during our observations. The colour variations are mildly chromatic on long timescales for two of them. The nature of the long-term variability of 3C 66A and S5 0954+658 is consistent with a model of a non-thermal variable component that has a continuous injection of relativistic electrons with power law distributions around 4.3 and 4.6, respectively. However, the long-term flux and colour variability of BL Lac suggests that these can arise from modest changes in velocities or viewing angle toward the emission region, leading to variations in the Doppler boosting of the radiation by a factor ~1.2 over the period of these observations.

We report on a recent multi-band optical photometric and polarimetric observational campaign of the blazar OJ 287 which was carried out during September 2016 -- December 2017. We employed nine telescopes in Bulgaria, China, Georgia, Japan, Serbia, Spain and the United States. We collected over 1800 photometric image frames in BVRI bands and over 100 polarimetric measurements over ~175 nights. In 11 nights with many quasi-simultaneous multi-band (V, R, I) observations, we did not detect any genuine intraday variability in flux or color. On longer timescales, multiple flaring events were seen. Large changes in color with respect to time and in a color--magnitude diagram were seen, and while only a weak systematic variability trend was noticed in color with respect to time, the color--magnitude diagram shows a bluer-when-brighter trend. Large changes in the degree of polarization, and substantial swings in the polarization angle were detected. The fractional Stokes parameters of the polarization showed a systematic trend with time in the beginning of these observations, followed by chaotic changes and then an apparently systematic variation at the end. These polarization changes coincide with the detection and duration of the source at very high energies as seen by VERITAS. The spectral index shows a systematic variation with time and V-band magnitude. We briefly discuss possible physical mechanisms that could explain the observed flux, color, polarization, and spectral variability.

The flat-spectrum radio quasar 4C $+$71.07 is a high-redshift ($z=2.172$), $\gamma$-loud blazar whose optical emission is dominated by the thermal radiation from accretion disc. 4C $+$71.07 has been detected in outburst twice by the AGILE $\gamma$-ray satellite during the period end of October - mid November 2015, when it reached a $\gamma$-ray flux of the order of $F_{\rm E>100\,MeV} = (1.2 \pm 0.3)\times 10^{-6}$ photons cm$^{-2}$ s$^{-1}$ and $F_{\rm E>100\,MeV} = (3.1 \pm 0.6)\times 10^{-6}$ photons cm$^{-2}$ s$^{-1}$, respectively, allowing us to investigate the properties of the jet and of the emission region. We investigated its spectral energy distribution by means of almost simultaneous observations covering the cm, mm, near-infrared, optical, ultra-violet, X-ray and $\gamma$-ray energy bands obtained by the GASP-WEBT Consortium, the Swift and the AGILE and Fermi satellites. The spectral energy distribution of the second $\gamma$-ray flare (the one whose energy coverage is more dense) can be modelled by means of a one-zone leptonic model, yielding a total jet power of about $4\times10^{47}$ erg s$^{-1}$. During the most prominent $\gamma$-ray flaring period our model is consistent with a dissipation region within the broad-line region. Moreover, this class of high-redshift, large-mass black-hole flat-spectrum radio quasars might be good targets for future $\gamma$-ray satellites such as e-ASTROGAM.

We present a multi-wavelength spectral and temporal investigation of OJ 287 emission during its strong optical-to-X-ray activity between July 2016 - July 2017. The daily $\gamma$-ray fluxes from \emph{Fermi}-LAT are consistent with no variability. The strong optical-to-X-ray variability is accompanied by a change in power-law spectral index of the X-ray spectrum from $< 2$ to $>2$, with variations often associated with changes in optical polarization properties. Cross-correlations between optical-to-X-ray emission during four continuous segments show simultaneous optical-ultraviolet (UV) variations while the X-ray and UV/optical are simultaneous only during the middle two segments. In the first segment, the results suggest X-rays lag the optical/UV, while in the last segment X-rays lead by $\sim$ 5-6 days. The last segment also shows a systematic trend with variations appearing first at higher energies followed by lower energy ones. The LAT spectrum before the VHE activity is similar to preceding quiescent state spectrum while it hardens during VHE activity period and is consistent with the extrapolated VHE spectrum during the latter. Overall, the broadband spectral energy distributions (SEDs) during high activity periods are a combination of a typical OJ 287 SED and an HBL SED, and can be explained in a two-zone leptonic model, with the second zone located at parsec scales, beyond the broad line region, being responsible for the HBL-like spectrum. The change of polarization properties from systematic to chaotic and back to systematic, before, during and after the VHE activity, suggest dynamic roles for magnetic fields and turbulence.

We report the AGILE detection and the results of the multifrequency follow-up observations of a bright $\gamma$-ray flare of the blazar 3C 279 in June 2015. We use AGILE-GRID and Fermi-LAT $\gamma$-ray data, together with Swift-XRT, Swift-UVOT, and ground-based GASP-WEBT optical observations, including polarization information, to study the source variability and the overall spectral energy distribution during the $\gamma$-ray flare. The $\gamma$-ray flaring data, compared with as yet unpublished simultaneous optical data which allow to set constraints on the big blue bump disk luminosity, show very high Compton dominance values of $\sim 100$, with a ratio of $\gamma$-ray to optical emission rising by a factor of three in a few hours. The multi-wavelength behavior of the source during the flare challenges one-zone leptonic theoretical models. The new observations during the June 2015 flare are also compared with already published data and non-simultaneous historical 3C 279 archival data.

We monitored BL Lacertae frequently during 2014 - 2016 when it was generally in a high state. We searched for intra-day variability for 43 nights using quasi-simultaneous measurements in the B, V, R, and I bands (totaling 143 light curves); the typical sampling interval was about eight minutes. On hour-like timescales, BL Lac exhibited significant variations during 13 nights in various optical bands. Significant spectral variations are seen during most of these nights such that the optical spectrum becomes bluer when brighter. The amplitude of variability is usually greater for longer observations but is lower when BL Lac is brighter. No evidence for periodicities or characteristic variability time-scales in the light curves was found. The color variations are mildly chromatic on long timescales.

We present a multi-wavelength spectral and temporal analysis of the blazar OJ 287 during its recent activity between December 2015 -- May 2016, showing strong variability in the near-infrared (NIR) to X-ray energies with detection at $\gamma$-ray energies as well. Most of the optical flux variations exhibit strong changes in polarization angle and degree. All the inter-band time lags are consistent with simultaneous emissions. Interestingly, on days with excellent data coverage in the NIR--UV bands, the spectral energy distributions (SEDs) show signatures of bumps in the visible--UV bands, never seen before in this source. The optical bump can be explained as accretion-disk emission associated with the primary black hole of mass $\sim \rm 1.8 \times10^{10} M_{\odot}$ while the little bump feature in the optical-UV appears consistent with line emission. Further, the broadband SEDs extracted during the first flare and during a quiescent period during this span show very different $\gamma$-ray spectra compared to previously observed flare or quiescent spectra. The probable thermal bump in the visible seems to have been clearly present since May 2013, as found by examining all available NIR-optical observations, and favors the binary super-massive black hole model. The simultaneous multi-wavelength variability and relatively weak $\gamma$-ray emission that shows a shift in the SED peak is consistent with $\gamma$-ray emission originating from inverse Compton scattering of photons from the line emission that apparently contributes to the little blue bump.

The blazar 3C454.3 exhibited a strong flare seen in gamma-rays, X-rays, and optical/NIR bands during 3--12 December 2009. Emission in the V and J bands rose more gradually than did the gamma-rays and soft X-rays, though all peaked at nearly the same time. Optical polarization measurements showed dramatic changes during the flare, with a strong anti-correlation between optical flux and degree of polarization (which rose from ~ 3% to ~ 20%) during the declining phase of the flare. The flare was accompanied by large rapid swings in polarization angle of ~ 170 degree. This combination of behaviors appear to be unique. We have cm-band radio data during the same period but they show no correlation with variations at higher frequencies. Such peculiar behavior may be explained using jet models incorporating fully relativistic effects with a dominant source region moving along a helical path or by a shock-in-jet model incorporating three-dimensional radiation transfer if there is a dominant helical magnetic field. We find that spectral energy distributions at different times during the flare can be fit using modified one-zone models where only the magnetic field strength and particle break frequencies and normalizations need change. An optical spectrum taken at nearly the same time provides an estimate for the central black hole mass of ~ 2.3 * 10^9 M_sun. We also consider two weaker flares seen during the $\sim 200$ d span over which multi-band data are available. In one of them, the V and J bands appear to lead the $\gamma$-ray and X-ray bands by a few days; in the other, all variations are simultaneous.

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 examine the 2008-2016 $\gamma$-ray and optical light curves of three bright BL Lac objects, 0716+714, MRK 421, BL Lac, which exhibit large structured variability. We searched for periodicities by using a fully Bayesian approach. For two out of three sources investigated no significant periodic variability was found. In the case of BL Lac we detected a periodicity of ~ 680 days. Although the signal related to this is modest, the coincidence of the periods in both gamma and optical bands is indicative of a physical relevance. Considering previous literature results, possibly related $\gamma$-ray and optical periodicities of about one year time scale are proposed in 4 bright $\gamma$-ray blazars out of the 10 examined in detail. Comparing with results from periodicity search of optical archives of quasars, the presence of quasi-periodicities in blazars might be more frequent by a large factor. This suggests the intriguing possibility that the basic conditions for their observability are related to the relativistic jet in the observer direction, but the overall picture remains uncertain.

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.

Here we report on the results of the WEBT photo-polarimetric campaign targeting the blazar S5~0716+71, organized in March 2014 to monitor the source simultaneously in BVRI and near IR filters. The campaign resulted in an unprecedented dataset spanning $\sim 110$\,h of nearly continuous, multi-band observations, including two sets of densely sampled polarimetric data mainly in R filter. During the campaign, the source displayed pronounced variability with peak-to-peak variations of about $30\%$ and "bluer-when-brighter" spectral evolution, consisting of a day-timescale modulation with superimposed hourlong microflares characterized by $\sim 0.1$\,mag flux changes. We performed an in-depth search for quasi-periodicities in the source light curve; hints for the presence of oscillations on timescales of $\sim 3$\,h and $\sim 5$\,h do not represent highly significant departures from a pure red-noise power spectrum. We observed that, at a certain configuration of the optical polarization angle relative to the positional angle of the innermost radio jet in the source, changes in the polarization degree led the total flux variability by about 2\,h; meanwhile, when the relative configuration of the polarization and jet angles altered, no such lag could be noted. The microflaring events, when analyzed as separate pulse emission components, were found to be characterized by a very high polarization degree ($> 30\%$) and polarization angles which differed substantially from the polarization angle of the underlying background component, or from the radio jet positional angle. We discuss the results in the general context of blazar emission and energy dissipation models.

After several years of quiescence, the blazar CTA 102 underwent an exceptional outburst in 2012 September-October. The flare was tracked from gamma-ray to near-infrared frequencies, including Fermi and Swift data as well as photometric and polarimetric data from several observatories. An intensive GASP-WEBT collaboration campaign in optical and NIR bands, with an addition of previously unpublished archival data and extension through fall 2015, allows comparison of this outburst with the previous activity period of this blazar in 2004-2005. We find remarkable similarity between the optical and gamma-ray behaviour of CTA 102 during the outburst, with a time lag between the two light curves of ~1 hour, indicative of co-spatiality of the optical and gamma-ray emission regions. The relation between the gamma-ray and optical fluxes is consistent with the SSC mechanism, with a quadratic dependence of the SSC gamma-ray flux on the synchrotron optical flux evident in the post-outburst stage. However, the gamma-ray/optical relationship is linear during the outburst; we attribute this to changes in the Doppler factor. A strong harder-when-brighter spectral dependence is seen both the in gamma-ray and optical non-thermal emission. This hardening can be explained by convexity of the UV-NIR spectrum that moves to higher frequencies owing to an increased Doppler shift as the viewing angle decreases during the outburst stage. The overall pattern of Stokes parameter variations agrees with a model of a radiating blob or shock wave that moves along a helical path down the jet.

We search for low-dimensional chaotic signatures in the optical lightcurve of the $Kepler$ field blazar W2R 1926+42. The frequently used correlation integral method is employed in our analysis. We find no apparent evidence for the presence of low-dimensional chaos in the lightcurve. If further confirmed, these results could be of importance for modeling the blazar emission mechanisms.

We present our optical photometric observations of three TeV blazars, PKS 1510-089, PG 1553+113 and Mrk 501 taken using two telescopes in India, one in Bulgaria, one in Greece and one in Serbia during 2012 - 2014. These observations covered a total of 95 nights with a total of 202 B filter frames, 247 images in V band, 817 in R band while 229 images were taken in the I filter. This work is focused on multi-band flux and colour variability studies of these blazars on diverse timescales which are useful in understanding the emission mechanisms. We studied the variability characteristics of above three blazars and found all to be active over our entire observational campaigns. We also searched for any correlation between the brightness of the sources and their colour indices. During the times of variability, no significant evidence for the sources to display spectral changes correlated with magnitude was found on timescales of a few months. We briefly discuss the possible physical mechanisms most likely responsible for the observed flux variability.

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 study the long term Kepler light curve of the blazar W2R 1926+42 ($\sim$ 1.6 years) which indicates a variety of variability properties during different intervals of observation. The normalized excess variance, $F_{\rm var}$ ranges from 1.8 % in the quiescent phase and 43.3 % in the outburst phase. We find no significant deviation from linearity in the $F_{\rm var}$-flux relation. Time series analysis is conducted using the Fourier power spectrum and the wavelet analysis methods to study the power spectral density (PSD) shape, infer characteristic timescales and statistically significant quasi-periodic oscillations (QPOs). A bending power law with an associated timescale of $T_B = 6.2^{+6.4}_{-3.1}$ hours is inferred in the PSD analysis. We obtain a black hole mass of $M_\bullet = (1.5 - 5.9) \times 10^7 M_\odot$ for the first time using $F_{\rm var}$ and the bend timescale for this source. From a mean outburst lifetime of days, we infer a distance from the jet base $r \leq 1.75$ pc indicating that the outburst originates due to a shock. A possible QPO peaked at 9.1 days and lasting 3.4 cycles is inferred from the wavelet analysis. Assuming that the QPO is a true feature, $r = (152 - 378)~ G M_\bullet/c^2$ and supported by the other timing analysis products such as a weighted mean PSD slope of $-1.5 \pm 0.2$ from the PSD analysis, we argue that the observed variability and the weak and short duration QPO could be due to jet based processes including orbital features in a relativistic helical jet and others such as shocks and turbulence.

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.

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 observed the prototype blazar, BL Lacertae, extensively in optical and radio bands during an active phase in the period 2010--2013 when the source showed several prominent outbursts. We searched for possible correlations and time lags between the optical and radio band flux variations using multifrequency data to learn about the mechanisms producing variability. During an active phase of BL Lacertae, we searched for possible correlations and time lags between multifrequency light curves of several optical and radio bands. We tried to estimate any possible variability timescales and inter-band lags in these bands. We performed optical observations in B, V, R and I bands from seven telescopes in Bulgaria, Georgia, Greece and India and obtained radio data at 36.8, 22.2, 14.5, 8 and 4.8 GHz frequencies from three telescopes in Ukraine, Finland and USA. Significant cross-correlations between optical and radio bands are found in our observations with a delay of cm-fluxes with respect to optical ones of ~250 days. The optical and radio light curves do not show any significant timescales of variability. BL Lacertae showed many optical 'mini-flares' on short time-scales. Variations on longer term timescales are mildly chromatic with superposition of many strong optical outbursts. In radio bands, the amplitude of variability is frequency dependent. Flux variations at higher radio frequencies lead the lower frequencies by days or weeks. The optical variations are consistent with being dominated by a geometric scenario where a region of emitting plasma moves along a helical path in a relativistic jet. The frequency dependence of the variability amplitude supports an origin of the observed variations intrinsic to the source.

We present the results of extensive multi-band intra-night optical monitoring of BL Lacertae during 2010--2012. BL Lacertae was very active in this period and showed intense variability in almost all wavelengths. We extensively observed it for a total for 38 nights; on 26 of them observations were done quasi-simultaneously in B, V, R and I bands (totaling 113 light curves), with an average sampling interval of around 8 minutes. BL Lacertae showed significant variations on hour-like timescales in a total of 19 nights in different optical bands. We did not find any evidence for periodicities or characteristic variability time-scales in the light curves. The intranight variability amplitude is generally greater at higher frequencies and decreases as the source flux increases. We found spectral variations in BL Lacertae in the sense that the optical spectrum becomes flatter as the flux increases but in several flaring states deviates from the linear trend suggesting different jet components contributing to the emission at different times.

The occurrence of low-amplitude flux variations in blazars on hourly timescales, commonly known as microvariability, is still a widely debated subject in high-energy astrophysics. Several competing scenarios have been proposed to explain such occurrences, including various jet plasma instabilities leading to the formation of shocks, magnetic reconnection sites, and turbulence. In this letter we present the results of our detailed investigation of a prominent, five-hour-long optical microflare detected during recent WEBT campaign in 2014, March 2-6 targeting the blazar 0716+714. After separating the flaring component from the underlying base emission continuum of the blazar, we find that the microflare is highly polarized, with the polarization degree $\sim (40-60)\%$$\pm (2-10)\%$, and the electric vector position angle $\sim (10 - 20)$deg$\pm (1-8)$deg slightly misaligned with respect to the position angle of the radio jet. The microflare evolution in the $(Q,\,U)$ Stokes parameter space exhibits a looping behavior with a counter-clockwise rotation, meaning polarization degree decreasing with the flux (but higher in the flux decaying phase), and approximately stable polarization angle. The overall very high polarization degree of the flare, its symmetric flux rise and decay profiles, and also its structured evolution in the $Q-U$ plane, all imply that the observed flux variation corresponds to a single emission region characterized by a highly ordered magnetic field. As discussed in the paper, a small-scale but strong shock propagating within the outflow, and compressing a disordered magnetic field component, provides a natural, though not unique, interpretation of our findings.

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.

In this letter we present results from intra-night monitoring in three colors (BRI) of the blazar S4 0954+65 during its recent (2015 February) unprecedented high state. We find violent variations on very short time scales, reaching magnitude change levels of 0.1-0.2 mag/h. On some occasions, changes of ~0.1 mag are observed even within ~10 min. During the night of 14.02.2015 an exponential drop of ~0.7 magnitudes is detected for about 5 hours. Cross-correlation between the light curves does not reveal any detectable wavelength-dependent time delays, larger than ~5 min. Color changes "bluer-when-brighter" are observed on longer time scales. Possible variability mechanisms to explain the observations are discussed and a preference to the geometrical one is given.

In this work we search for signatures of low-dimensional chaos in the temporal behavior of the Kepler-field blazar W2R 1946+42. We use a publicly available, ~160000 points long and mostly equally spaced, light curve of W2R 1946+42. We apply the correlation integral method to both -- real datasets and phase randomized "surrogates". We are not able to confirm the presence of low-dimensional chaos in the light curve. This result, however, still leads to some important implications for blazar emission mechanisms, which are discussed.

The UV spectral range (1100 - 3000 A) contains the strongest resonance lines observed in active galactic nuclei (AGN). Analysis of UV line intensity ratios and profile shapes in quasar spectra provide diagnostics of physical and dynamical conditions in the broad line emitting region. This paper discusses properties of UV lines in type-1 AGN spectra, and how they lead an estimate of ionizing photon flux, chemical abundances, radius of the broad line emitting region and central black hole mass. These estimates are meaningfully contextualised through the 4D "eigenvector-1" (4DE1) formalism.

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.

Context. The international whole earth blazar telescope (WEBT) consortium planned and carried out three days of intensive micro-variability observations of S5 0716+714 from February 22, 2009 to February 25, 2009. This object was chosen due to its bright apparent magnitude range, its high declination, and its very large duty cycle for micro-variations. Aims. We report here on the long continuous optical micro-variability light curve of 0716+714 obtained during the multi-site observing campaign during which the Blazar showed almost constant variability over a 0.5 magnitude range. The resulting light curve is presented here for the first time. Observations from participating observatories were corrected for instrumental differences and combined to construct the overall smoothed light curve. Methods. Thirty-six observatories in sixteen countries participated in this continuous monitoring program and twenty of them submitted data for compilation into a continuous light curve. The light curve was analyzed using several techniques including Fourier transform, Wavelet and noise analysis techniques. Those results led us to model the light curve by attributing the variations to a series of synchrotron pulses. Results. We have interpreted the observed microvariations in this extended light curve in terms of a new model consisting of individual stochastic pulses due to cells in a turbulent jet which are energized by a passing shock and cool by means of synchrotron emission. We obtained an excellent fit to the 72-hour light curve with the synchrotron pulse model.

We report results from a 1 week multi-wavelength campaign to monitor the BL Lac object S5 0716+714 (on December 9-16, 2009). In the radio bands the source shows rapid (~ (0.5-1.5) day) intra-day variability with peak amplitudes of up to ~ 10 %. The variability at 2.8 cm leads by about 1 day the variability at 6 cm and 11 cm. This time lag and more rapid variations suggests an intrinsic contribution to the source's intraday variability at 2.8 cm, while at 6 cm and 11 cm interstellar scintillation (ISS) seems to predominate. Large and quasi-sinusoidal variations of ~ 0.8 mag were detected in the V, R and I-bands. The X-ray data (0.2-10 keV) do not reveal significant variability on a 4 day time scale, favoring reprocessed inverse-Compton over synchrotron radiation in this band. The characteristic variability time scales in radio and optical bands are similar. A quasi-periodic variation (QPO) of 0.9 - 1.1 days in the optical data may be present, but if so it is marginal and limited to 2.2 cycles. Cross-correlations between radio and optical are discussed. The lack of a strong radio-optical correlation indicates different physical causes of variability (ISS at long radio wavelengths, source intrinsic origin in the optical), and is consistent with a high jet opacity and a compact synchrotron component peaking at ~= 100 GHz in an ongoing very prominent flux density outburst. For the campaign period, we construct a quasi-simultaneous spectral energy distribution (SED), including gamma-ray data from the FERMI satellite. We obtain lower limits for the relativistic Doppler-boosting of delta >= 12-26, which for a BL\,Lac type object, is remarkably high.

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

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.

Context. 3C 454.3 is a very active flat spectrum radio quasar (blazar) that has undergone a recent outburst in all observed bands, including the optical. Aims. In this work we explore the short-term optical variability of 3C 454.3 during its outburst by searching for time delays between different optical bands. Finding one would be important for understanding the evolution of the spectrum of the relativistic electrons, which generate the synchrotron jet emission. Methods. We performed photometric monitoring of the object by repeating exposures in different optical bands (BVRI). Occasionally, different telescopes were used to monitor the object in the same band to verify the reliability of the smallest variations we observed. Results. Except on one occasion, where we found indications of a lag of the blue wavelengths behind the red ones, the results are inconclusive for most of the other cases. There were either no structures in the light curves to be able to search for patterns, or else different approaches led to different conclusions.

[Abridged] We compare broad emission line profiles and estimate line ratios for all major emission lines between Ly-alpha and H-beta in a sample of six quasars. The sources were chosen with two criteria in mind: the existence of high quality optical and UV spectra as well as the possibility to sample the spectroscopic diversity in the 4D Eigenvector 1 context . In the latter sense each source occupies a region (bin) in the FWHM(H-beta) vs. optical FeII strength plane that is significantly different from the others. High S/N H-beta emission line profiles are used as templates for modeling the other lines (Ly-alpha, CIV 1549, HeII 1640, Al III 1860, Si III] 1892, and Mg II 2800). We can adequately model all broad lines assuming the existence of three components distinguished by blueshifted, unshifted and redshifted centroids (indicated as blue, broad and very broad component respectively). The broad component (high electron density, low ionization parameter; high column density) is present in almost all type-1 quasars and therefore corresponds most closely to the classical broad line emitting region (the reverberating component). The blue component emission (lower electron density; high ionization; low column density) arises in less optically thick gas; it is often thought to arise in an accretion disk wind. The least understood component involves the very broad component (high ionization and large column density). It is perhaps the most distinguishing characteristic of quasars with FWHM H-beta > 4000 km/s that belong to the so-called Population B of our 4DE1 space. Population A quasars (FWHM H-beta < 4000 km/s) are dominated by broad component emission in H-beta and blue component emission in CIV 1549 and other high ionization lines. 4DE1 appears to be the most useful current context for revealing and unifying spectral diversity in type-1 quasars.

We report on 18 months of multiwavelength observations of the blazar 3C 454.3 (Crazy Diamond) carried out in July 2007-January 2009. We show the results of the AGILE campaigns which took place on May-June 2008, July-August 2008, and October 2008-January 2009. During the May 2008-January 2009 period, the source average flux was highly variable, from an average gamma-ray flux F(E>100MeV) > 200E-8 ph/cm2/s in May-June 2008, to F(E>100MeV)~80E-8 ph/cm2/s in October 2008-January 2009. The average gamma-ray spectrum between 100 MeV and 1 GeV can be fit by a simple power law (Gamma_GRID ~ 2.0 to 2.2). Only 3-sigma upper limits can be derived in the 20-60 keV energy band with Super-AGILE. During July-August 2007 and May-June 2008, RXTE measured a flux of F(3-20 keV)= 8.4E-11 erg/cm2/s, and F(3-20 keV)=4.5E-11 erg/cm2/s, respectively and a constant photon index Gamma_PCA=1.65. Swift/XRT observations were carried out during all AGILE campaigns, obtaining a F(2-10 keV)=(0.9-7.5)E-11 erg/cm2/s and a photon index Gamma_XRT=1.33-2.04. BAT measured an average flux of ~5 mCrab. GASP-WEBT monitored 3C 454.3 during the whole 2007-2008 period from the radio to the optical. A correlation analysis between the optical and the gamma-ray fluxes shows a time lag of tau=-0.4 days. An analysis of 15 GHz and 43 GHz VLBI core radio flux observations shows an increasing trend of the core radio flux, anti- correlated with the higher frequency data. The modeling SEDs, and the behavior of the long-term light curves in different energy bands, allow us to compare the jet properties during different emission states, and to study the geometrical properties of the jet on a time-span longer than one year.

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.

Aims. The main purpose of this paper is to study time delays between the light variations in different wavebands for a sample of quasars. Measuring a reliable time delay for a large number of quasars may help constraint the models of their central engines. The standard accretion disk irradiation model predicts a delay of the longer wavelengths behind the shorter ones, a delay that depends on the fundamental quasar parameters. Since the black hole masses and the accretion rates are approximately known for the sample we use, one can compare the observed time delays with the expected ones. Methods. We applied the interpolation cross-correlation function (ICCF) method to the Giveon et al. sample of 42 quasars, monitored in two (B and R) colors, to find the time lags represented by the ICCF peaks. Different tests were performed to assess the influence of photometric errors, sampling, etc., on the final result. Results. We found that most of the objects show a delay in the red light curve behind the blue one (a positive lag), which on average for the sample is about +4 days (+3 for the median), although the scatter is significant. These results are broadly consistent with the reprocessing model, especially for the well-sampled objects. The normalized time-lag deviations do not seem to correlate significantly with other quasar properties, including optical, radio, or X-ray measurables. On the other hand, many objects show a clear negative lag, which, if real, may have important consequences for the variability models.

The main purpose of this work is to improve the existing knowledge about the most powerful engines in the Universe - quasars. Although a lot is already known, we still have only a vague idea how these engines work exactly, why they behave as they do, and what the relation is between their evolution and the evolution of their harboring galaxy. Methods we used are based on optical spectroscopy of visually bright quasars, many of which have recently been discovered as X-ray sources, but eventually missed in color-selected surveys. The spectra typically cover the 4200-7000 AA region, allowing measurements of the characteristics of the hydrogen lines, the FeII contribution, and other lines of interest. We present accurate redshift estimates and Seyfert type classification of the objects. We also show that the contribution of the host galaxy to the optical continuum is non-negligible in many cases, as is the intrinsic AGN absorption. Consequences of not correcting for those factors when estimating different quasar parameters are discussed. We also find some evidence of a non-unity slope in the relation between the internal extinction based on the Balmer decrement and the one on the optical continuum slope, implying, if further confirmed, the intriguing possibility that some absorbing material might actually be located between the continuum source and the broad-line region.

We have measured the projected rotational velocities (vsini) in a number of symbiotic stars and M giants using high resolution spectroscopic observations. On the basis of our measurements and data from the literature, we compare the rotation of mass-donors in symbiotics with vsini of field giants and find that: (1) the K giants in S-type symbiotics rotate at vsini>4.5 km/s, which is 2-4 times faster than the field K giants; (2) the M giants in S-type symbiotics rotate on average 1.5 times faster than the field M giants. Statistical tests show that these differences are highly significant: p-value < 0.001 in the spectral type bins K2III-K5III, M0III-M6III, and M2III-M5III; (3) our new observations of D'-type symbiotics also confirm that they are fast rotators. As a result of the rapid rotation, the cool giants in symbiotics should have 3-30 times larger mass loss rates. Our results suggest also that bipolar ejections in symbiotics seem to happen in objects where the mass donors rotate faster than the orbital period. All spectra used in our series of papers can be obtained upon request from the authors.

We present the results of the blazar 3C 345 monitoring in Johnson-Cousins BVRI bands for the period 1996 - 2006. We have collected 29 V and 43 R data points for this period; the BI light curves contain a few measurements only. The accuracy of our photometry is not better than 0.03 mag in the VR bands. The total amplitude of the variability obtained from our data is 2.06 mag in the V band and 2.25 mag in the R one. 3C 345 showed periods of flaring activity during 1998/99 and 2001: a maximum of the blazar brightness was detected in 2001 February - 15.345 mag in the V band and 14.944 mag in the R one. We confirm that during brighter stages 3C 345 becomes redder; for higher fluxes the colour index seems to be less dependent on the magnitude. The intra-night monitoring of 3C 345 in three consecutive nights in 2001 August revealed no significant intra-night variability; 3C 345 did not show evident flux changes over timescales of weeks around the period of the intra-night monitoring. This result supports the existing facts that intra-night variability is correlated with rapid flux changes rather than with specific flux levels.

We have measured the projected rotational velocities (v sin i) of the mass donors for 29 S-type symbiotic stars using high resolution spectroscopic observations and the cross-correlation function (CCF) method. The results of the CCF have been controlled with synthetic spectra. The typical rotational velocity of the K and M giants in S-type symbiotics appeared to be 4.5 < vsini < 11.7 km/s. In a sub-sample of 16 S--type symbiotic stars (with known orbital periods and well measured vsini) 15 have deviations from synchronization less then the 3-sigma level. This means that we did not find evidence for a statistically significant deviation from the synchronization for any of these 15 objects. The deviation from synchronization is statistically significant (at confidence level >99%) only for the recurrent nova RS Oph. For 22 S-type symbiotics we give clues as to what their orbital periods could be.

[Abridged] We have been exploring a spectroscopic unification for all known types of broad line emitting AGN. The 4D Eigenvector 1 (4DE1) parameter space shows promise as a unification capable of organizing quasar diversity on a sequence primarily governed by Eddington ratio. This paper considers the role of CIV 1549 measures with special emphasis on the CIV 1549 line shift as a principal 4DE1 diagnostic. We use HST archival spectra for 130 sources with S/N high enough to permit reliable CIV 1549 broad component (BC) measures. We find a CIV 1549 BC profile blueshift that is strongly concentrated among (largely radio-quiet: RQ) sources with FWHM(H beta BC) < 4000 km/s (which we call Population A). Narrow line Seyfert 1 (NLSy1, with FWHM H beta < 2000 km/s) sources belong to this population but do not emerge as a distinct class. The systematic blueshift, widely interpreted as arising in a disk wind/outflow, is not observed in broader lined AGN which we call Population B. We find new correlations between FWHM(CIV 1549 BC) and CIV 1549 line shift as well as the equivalent width of CIV 1549. They are seen only in Pop. A sources. CIV 1549 measures enhance the apparent dichotomy at FWHM(Hbeta BC) approx. 4000 \kms\ (Sulentic et al. 2000) suggesting that it has more significance in the context of Broad Line Region structure than the more commonly discussed RL vs. RQ dichotomy. Black hole masses computed from FWHM CIV 1549 BC for about 80 AGN indicate that the CIV 1549 width is a poor virial estimator. Comparison of mass estimates derived from Hbeta BC and CIV 1549 reveals that the latter show different and nonlinear offsets for population A and B sources. A significant number of sources also show narrow line CIV 1549 emission. We present a recipe for CIV 1549 narrow component extraction.

In this work we connect some measurable properties of the CIV1549 emission line with the quasar accretion rates (Eddington ratios). A tight correlation is found for a sample of more than a hundred nearby objects, suggesting a possible method for a relatively accurate estimate of the Eddington ratio of high-redshift quasars, at least for the radio-quiet ones. This paper further confirms the existing notion that the CIV changes (shifts) are mostly driven by the accretion rate.

This paper presents results of a search for short-term variability in the optical band of selected high-luminosity, high-redshift radio-quiet quasars. Each quasar has been monitored typically for 2 - 4 hours with a time resolution of 2 - 5 minutes and a photometric accuracy of about 0.01 - 0.02 mag. Due to the significant redshift (z>2), the covered wavelength range falls into the UV region (typically 1500 - 2500A). We found no statistical evidence for any continuum variations larger than 0.01 - 0.02 for any of the monitored objects. Our results suggest that the presence of a short-term variability in radio-quiet quasars is unlikely even in the UV region, contrary to reports by other authors. This conclusion holds true at least for high-luminosity (large black hole mass and accretion rate?) objects. The results are consistent with the idea that significant short-term (less than 1 hour) variations in AGN, where observed, should be attributed primarily to processes in a relativistic jet.

We present composite UV spectra for low redshift Type 1 AGN binned to exploit the information content of the Eigenvector 1 (E1) parameter space. Composite spectra allow a decomposition of the CIV1549 line profile - one of the strongest high-ionization lines. The simplest CIV decomposition into narrow (NLR), broad (BLR) and very broad (VBLR) components suggests that different components have an analog in Hb with two major exceptions. VBLR emission is seen only in population B (FWHM(Hb)>4000 km/s) sources. A blue shifted/asymmetric BLR component is seen only in pop. A (FWHM(Hb)<4000 km/s) HIL such as CIV. The blueshifted component is thought to arise in a wind or outflow. Our analysis suggests that such a wind can only be produced in pop. A (almost all radio-quiet) sources where the accretion rate is relatively high. Comparison between broad UV lines in radio-loud (RL) and radio-quiet (RQ) sources shows few significant differences. Clear evidence is found for a narrow CIV component in most radio-loud sources. We find also some indirect indications that the black hole (BH) spin, rather than BH mass or accretion rate is a key trigger in determining whether an object will be RL or RQ. We find a ten-fold decrease in EW CIV with Eddington ratio (decreasing from ~1 to \~0.01) while NV shows no change. These trends suggest a luminosity-independent "Baldwin effect" where the physical driver may be the Eddington ratio.

We present results from a recent broad-band monitoring in optics of the Seyfert 1 type galaxy Mrk 279. We build and analyse the BVRI light curve covering a period of seven years (1995 - 2002). We also show some evidence for the existence of two different states in brightness and suggest, based on a modelling of the optical continuum, that these states may result from transition between a thin disk and an ADAF accretion modes. We assume that the short-term variability is due to a reprocessing of a variable X-ray emission from an inner ADAF part of the flow, while the long-term one may be a result from a change of the transition radius. Our tests show a good match with the observations for a reasonable set of accretion parameters, rather close to the ones, expected for Mrk 279.

We consider the properties of radio-loud (RL) AGN in the context of the Eigenvector 1 (E1) parameter space. RL sources show a restricted E1 parameter space occupation relative to the radio-quiet (RQ) majority. The Fanaroff-Riley II ``parent population'' of relatively un-boosted RL sources (median radio/optical flux ratio ~490) shows the most restricted occupation. RL sources have different broad line properties (and inferred black hole masses and Eddington ratios). FWHM H_beta for the broad line component in RL sources are at least twice as large as the RQ majority. The average broad FeII emission line strength is also about half that for RQ sources. Our sample suggests that the RL cutoff occurs near R_k=70 or logP(6cm)=32.0 ergs/s/Hz. Sources below this cutoff are RQ although we cannot rule out the existence of a distinct intermediate population. We show that the Doppler boosted core-dominated RL sources (median flux ratio ~1000) lie towards smaller FWHM(H_beta_bc) and stronger FeII in E1 as expected if the lines arise in an accretion disk. Our subsample of superluminal sources, with orientation inferred from the synchrotron self Compton model, reinforce this general E1 trend and allow us to estimate the role of source orientation in driving E1 domain occupation.

We investigate the radial velocity difference between the [OIII]5007,4959 and H-beta lines for a sample of 200 low redshift AGN. We identify seven objects showing an [OIII]5007 blueshift relative to H-beta with amplitude larger than 250 km/s (blue "outliers"). These line shifts are found in sources where the broad high ionization lines (e.g. CIV1549) also show a large systematic blueshift. Such blueshifts occur only in the population A region of the Eigenvector 1 parameter domain (that also contains NLSy1 sources). We suggest that [OIII] blueshifts are also associated with the high ionization outflow originating in these sources. This is a direct kinematic linkage between narrow and broad line emitting gas.

[LaTeX removed] Recent work has shown that it is possible to systematize quasar spectral diversity in a parameter space called ``Eigenvector 1'' (E1). We present median AGN spectra for fixed regions of the E1 (optical) parameter space (FWHM(H-beta) vs. equivalent width ratio RFE=W(FeII4570)/W(H-beta). Comparison of the median spectra for different regions show considerable differences. We suggest that an E1-driven approach to median/average spectra emphasizes significant differences between AGN, and offers more insights into AGN physics and dynamics than a single population median/average derived from a large and heterogeneous sample of sources. We find that the H-beta broad component line profile changes in shape along the E1 sequence both in average centroid shift and asymmetry. While objects with FWHM(H-beta)< 4000 km/s are well fitted by a Lorentz function, AGN with FWHM(H-beta)> 4000 km/s are well fitted if two broad line components are used: a broad (the "classical" broad line component) and a very broad/redshifted component.

The broad line profiles from a nonplanar (warped) accretion disk as a result from reprocessing of central hard radiation are studied in this paper. The disk is warped because of the Bardeen-Petterson effect. Obtained profiles depend significantly on the line of sight - they are asymmetric as well as frequency-shifted in most cases. Since such frequency shifted asymmetric profiles are generally not observed in AGN, the most plausible interpretation is that either the major part of the broad emission does not originate from the disk or the disk is in general planar.